FEATool Multiphysics
v1.17.1
Finite Element Analysis Toolbox
|
This multiphysics example examines microfluidic flow and coupled mass transport due to electroosmosis in a channel with a T-shaped junction. With application of an electric field in a micro channel a flow effect is induced along the walls due to chemical reactions between the liquid and the wall material. This effect is called electro-osmotic flow (EOF), and can for example be used in chemical analysis to separate reactants and chemical species. The flow effect is here modeled with the Helmholtz-Smoluchowski boundary condition [1,2] which prescribes a tangential slip velocity as
where μEOF = εε0ζ/μ is the electroosmotic mobility. In this model, the walls are aligned with the coordinate axes, and the boundary condition can therefore be simplified to -μEOFVx and -μEOFVy for the horizontal and vertical boundaries, respectively (where Vx and Vy are the gradients of the electric potential V). A passive scalar representing a chemical analyte is injected at the inlet for a prescribed duration (concentration 1 for t <= 30 and linearly decreasing for the next 10 µs) whereby the transport, diffusion, and path of the species is tracked through the channel.
This model is available as an automated tutorial by selecting Model Examples and Tutorials... > Multiphysics > Electro-Osmotic Flow from the File menu. Or alternatively, follow the step-by-step instructions below.
In the New Model dialog box, select 2D for the Space Dimensions, and select Conductive Media DC from the Select Physics drop-down menu (the physics modes for fluid flow and mass transport will be added later). Leave the space dimension and dependent variable names to their default values and press OK to finish the physics mode selection.
The geometry consists of a 45 by 55 micrometer T-shaped junction. One can either join two rectangles or use the polygon tool to create the shape of the geometry.
Enter the following data into the Point coordinates table.
x | y | |
---|---|---|
1 | 0 | 2e-05 |
2 | 5e-05 | 2e-05 |
3 | 5e-05 | 0 |
4 | 5.5e-05 | 0 |
5 | 5.5e-05 | 4.5e-05 |
6 | 5e-05 | 4.5e-05 |
7 | 5e-05 | 2.5e-05 |
8 | 0 | 2.5e-05 |
Press OK to finish and close the dialog box.
Press the Grid mode button in the Mode Toolbar to change mode and generate a default grid. The initial grid is usually quite coarse and will not give a very accurate solution, so click on the Refine button three or more times to create a finer and more accurate grid.
1e3
into the Density edit field.1e-3
into the Viscosity edit field.The species diffusion is assumed to be constant and small, while convection is due to the osmotic-velocity given by the Navier-Stokes physics mode (coupling u and v for the x and y-velocities, respectively).
1e-10
into the Diffusion coefficient edit field.u
into the Convection velocity in x-direction edit field.Enter v
into the Convection velocity in y-direction edit field.
For convection dominated transport problems numerical diffusion can help to ensure a more physical and non-oscillatory solution. In this case isotropic diffusion is added (most diffusive), but more accurate streamline diffusion (less diffusive), or high order shock-capturing artificial diffusion can also be used, but are more costly and may take longer to compute a solution.
Select the Check to enable isotropic artificial diffusion check box.
Boundary conditions are defined in Boundary Mode and describes how the model interacts with the surrounding environment.
The electric potential should be set to 100 V and the left side inlet, and 30 V and 0 at the top and bottom boundaries, respectively. Neutral current flow conditions are used for the walls.
100
into the Electric potential edit field.30
into the Electric potential edit field.0
into the Electric potential edit field.There should not be any forced flows or velocities prescribed, so first select Neutral conditions for all the inlets and outlets.
Helmholtz-Smoluchowski boundary condition couples the electric field and gradient of the potential to the velocity at the walls and is driving the flow here. Since the walls are parallel to the coordinate axis the tangential component boundary conditions can be reduced to -5e-8*Vy and -5e-8*Vx in the x and y-directions.
Enter -5e-8*Vx
into the Velocity in x-direction edit field.
-5e-8*Vy
into the Velocity in y-direction edit field.No species transport should be allowed through the walls and insulation conditions are therefore appropriate there.
The inlet concentration should be 1 for t < 3e-5 and a linear decreasing function between 3e-5 < t < 4e-5. This can be represented as the switch condition (t<=3e-5)+(4e-5-t)/1e-5*(t>3e-5&t<4e-5) (where the logical conditions <, <=, and > either evaluate to zero or one.)
Enter (t<=3e-5)+(4e-5-t)/1e-5*(t>3e-5&t<4e-5)
into the Concentration edit field.
Now that the problem is fully specified, press the Solve Mode Toolbar button to switch to solve mode. Open the Solver Settings dialog box and select the time dependent solver, time step 1e-5 and simulation time 1.2e-3 seconds.
1e-5
into the Time step size edit field.Enter 1.2e-3
into the Duration of time-dependent simulation (maximum time) edit field.
Press OK to apply the selected solver settings, and then press the = Toolbar button to call the solver.
After the problem has been solved FEATool will automatically switch to postprocessing mode and show the electric potential at the final time step.
Also plot the velocity field to see how the majority of the flow takes the lower path.
30
into the Arrow spacing in y-direction (integer or coordinate vector) edit field.30
into the Arrow spacing in x-direction (integer or coordinate vector) edit field.Finally plot the concentration at different times to see how the species is transported from the inlet to the outlet.
Press OK to plot and visualize the selected postprocessing options.
To create an animation of the whole process one can use the Animate/Playback Solution... option from the Post menu.
Also note that in order to save computational resources it is possible to first compute stationary solutions for the electric potential and velocity field (as they will not change with time), and keep them constant (de-activated subdomains) while computing the mass transport problem (thus reducing the memory and time required for the time-dependent solver). See the Multi-Simulation Heat Exchanger for an example of this solution approach.
The electro-osmotic flow multiphysics model has now been completed and can be saved as a binary (.fea) model file, or exported as a programmable MATLAB m-script text file, or GUI script (.fes) file.
[1] G.E. Karniadakis and A. Beskok, Microflows: fundamentals and simulation. Springer-Verlag, New York, Berlin, Heidelberg, 2001.
[2] R.J. Yang, L.M. Fu, Y.C. Lin, Electroosmotic Flow in Microchannels. J. Colloid and Interface Science, 239:98-105, November 2001.