This paper presents composable behavioral models and a schematic-based simulation methodology to enable topdown design of electrokinetic (EK) lab-on-a-chip (LoC). Complex EK LoCs are shown to be decomposable into a system of elements with simple geometry and specific function. Parameterized and analytical models are developed to describe the electric and biofluidic behavior within each element. Electric and biofluidic pins at element terminals support the communication between adjacent elements in a simulation schematic. An analog hardware description language implementation of the models is used to simulate LoC subsystems for micromixing and electrophoretic separation. Both direct current (dc) and transient analysis can be performed to capture the influence of system topology, element sizes, material properties, and operational parameters on LoC system performance. Accuracy (relative error generally less than 5%) and speedup (> 100×) of the schematic-based simulation methodology are demonstrated by comparison to experimental measurements and continuum numerical simulation. Index Terms-Behavioral model, dispersion, lab-on-a-chip, micromixing, schematic-based simulation, top-down design.NOMENCLATURE Variables A Maximum concentration of species band. c Concentration of species band. c m Channel width-averaged concentration of species band. c p pth moment of the concentration in a longitudinal filament of the species band. C e Electric conductivity of buffer, S/m. d n Mixing concentration coefficients. d m Mixing concentration coefficients at inlet of converging/diverging mixing intersection. D Molecular diffusivity of species, m 2 /s. E Electric field strength, V/cm. h Depth of microchannels, µm. I Electric current through buffer, A. L Length of microchannels, µm. L det Detector path length, µm.