You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Copy file name to clipboardExpand all lines: doc/source/examples/multiphysics/tracer-in-static-mixer/tracer-in-static-mixer.rst
+15-14
Original file line number
Diff line number
Diff line change
@@ -45,6 +45,8 @@ In this example, we use the auxiliary tracer physics to quantify the mixing of t
45
45
#. Simulating the transient flow with a small ``time step`` until the flow has reached the steady state;
46
46
#. Simulate the flow of the tracer with a larger ``time step`` while reusing the computed velocity field from the previous simulation.
47
47
48
+
This strategy reduces the computational cost of the simulation as the fluid flow solution is reused in the second step.
49
+
48
50
--------------
49
51
Parameter File
50
52
--------------
@@ -70,7 +72,7 @@ Simulation control
70
72
set output frequency = 5
71
73
end
72
74
73
-
#. We use a ``bdf1`` time integration scheme with a short ``time step = 1e-4`` seconds. After ``time end = 40e-4`` seconds, we consider that the velocity field has reached a steady state.
75
+
We use a ``bdf1`` time integration scheme with a short ``time step = 1e-4`` :math:`\text{s}`. After ``time end = 40e-4`` :math:`\text{s}`, we consider that the velocity field has reached a steady state. According to :doc:`../../sharp-immersed-boundary/3d-rbf-static-mixer/3d-rbf-static-mixer`, the pressure drop varies by less than :math:`5` % in the range [30, 40] :math:`\text{s}`.
74
76
75
77
Restart
76
78
******************
@@ -86,13 +88,12 @@ We use Lethe's checkpoint/restart mechanism to feed the steady-state flow inform
86
88
set restart = false
87
89
end
88
90
89
-
#. We ``checkpoint`` the simulation at every ``5`` time steps, both for safety (each time step takes a long time to complete) and for reuse in the second simulation.
90
-
#. We do not ``restart`` the simulation unless stopped before it reaches ``time end``.
91
+
We ``checkpoint`` the simulation at every ``5`` time steps, both for safety (each time step takes a long time to complete) and for reuse in the second simulation.
91
92
92
93
Multiphysics
93
94
******************
94
95
95
-
Both ``fluid dynamics`` and ``tracer`` are enabled for the first simulation. However, the tracer injection only begins at the second step.
96
+
Both ``fluid dynamics`` and ``tracer`` are enabled for the first simulation to ensure initialization of the ``tracer`` field. However, the tracer injection only begins at the second step.
96
97
97
98
.. code-block:: text
98
99
@@ -104,7 +105,7 @@ Both ``fluid dynamics`` and ``tracer`` are enabled for the first simulation. How
104
105
Physical Properties
105
106
*******************
106
107
107
-
In this case we consider that we have a passive tracer in water. The units used for the physical properties are centimeters and seconds.
108
+
In this case we consider that we have a passive tracer in water. The units used for the physical properties are :math:`\text{cm}` and :math:`\text{s}`.
108
109
109
110
.. code-block:: text
110
111
@@ -121,9 +122,9 @@ In this case we consider that we have a passive tracer in water. The units used
121
122
end
122
123
123
124
#. The ``tracer diffusivity model`` is ``immersed solid tanh``. This model is used in ``lethe-fluid-sharp`` for tracer flow percolating immersed solids.
124
-
#. The ``tracer diffusivity outside`` is ``1e-5``, as this is a typical value for a passive tracer in a liquid.
125
-
#. The ``tracer diffusivity inside`` is set to ``1e-10``. The low value prevents diffusivity inside the solid while providing numerical stability.
126
-
#. The ``thickness`` is ``5e-1``. At the scale of the problem, this provides a smooth transition without generating oscillations between liquid and solid phases.
125
+
#. The ``tracer diffusivity outside`` is ``1e-5``:math:`\text{cm²/s}`, as this is a typical value for a passive tracer in a liquid.
126
+
#. The ``tracer diffusivity inside`` is set to ``1e-10``:math:`\text{cm²/s}`. The low value prevents diffusivity inside the solid while providing numerical stability (:math:`> 0`).
127
+
#. The ``thickness`` is ``5e-1``:math:`\text{cm}`. At the scale of the problem, this provides a smooth transition without generating oscillations between liquid and solid phases. The thickness is of the order of magnitude of the smallest cell length to restrict the transition to one cell thickness.
127
128
128
129
Tracer Boundary Conditions
129
130
***************************
@@ -163,7 +164,7 @@ Tracer Boundary Conditions
163
164
end
164
165
165
166
#. We use ``time dependent`` boundary conditions, defined at :math:`y<0` and :math:`11 < t < 61` to inject a pulse on the lower half of the inlet.
166
-
#. All other boundary conditions are ``outlet``. This condition is natural for the outlet of the problem. For lateral walls, this condition represents an impermeable wall.
167
+
#. All other boundary conditions are ``outlet``. This condition is natural for the outlet of the problem. For lateral walls, this condition represents an impermeable wall since velocity is perpendicular.
167
168
168
169
Post-processing
169
170
*******************
@@ -181,7 +182,7 @@ Post-processing
181
182
set tracer flow rate name = tracer_flow_rate
182
183
end
183
184
184
-
#. ``calculate average velocities`` is enabled, beginning at ``initial time = 30e-4``. This means that the last :math:`25\%` of the first simulation will be used to calculate the time-averaged fluid velocity profile.
185
+
#. ``calculate average velocities`` is enabled, beginning at ``initial time = 30e-4``:math:`\text{s}`. This means that the last :math:`25\%` of the first simulation will be used to calculate the time-averaged fluid velocity profile.
185
186
#. ``calculate tracer flow rate`` is enabled to provide data for mixing quantification.
186
187
187
188
Simulation 2: Tracer transport
@@ -201,8 +202,8 @@ Simulation control
201
202
...
202
203
end
203
204
204
-
#. We use a longer ``time step`` of ``10`` seconds.
205
-
#. We simulate until ``time end = 500`` seconds to allow the tracer to flow through the entire length of the domain.
205
+
#. We use a longer ``time step`` of ``10`` :math:`\text{s}`.
206
+
#. We simulate until ``time end = 500`` :math:`\text{s}` to allow the tracer to flow through the entire length of the domain.
206
207
207
208
Restart
208
209
******************
@@ -262,7 +263,7 @@ The simulation can be launched on multiple cores using ``mpirun`` and the ``leth
262
263
Results
263
264
-------
264
265
265
-
The following movie shows the evolution of tracer through the static mixer, both as a colored slice and colored streamlines:
266
+
The following movie shows the flow of the tracer through the static mixer, both as a colored slice and colored streamlines:
266
267
267
268
.. raw:: html
268
269
@@ -285,7 +286,7 @@ The tracer evolution through the inlet and outlet can be monitored by plotting t
As the Figure shows, the concentration of the tracer flattens as it flows. The gap between between the inlet and outlet peaks is of :math:`170` seconds. When compared to the theoretical time of :math:`150` seconds (:math:`d_x/u_x`, with :math:`d_x` the domain length and :math:`u_x` the inlet velocity), this difference can be explained by retention effects and the tortuous paths that the tracer travels through.
289
+
As the Figure shows, the concentration of the tracer flattens as it flows. The gap between between the inlet and outlet peaks is of :math:`170` :math:`\text{s}`. When compared to the theoretical time of :math:`150` :math:`\text{s}` (:math:`d_x/u_x`, with :math:`d_x` the domain length and :math:`u_x` the inlet velocity), this difference can be explained by retention effects and the tortuous paths that the tracer travels through.
0 commit comments