Partial Differential Equations - Modelling and ... - ResearchGate
Partial Differential Equations - Modelling and ... - ResearchGate Partial Differential Equations - Modelling and ... - ResearchGate
A Numerical Method for Fluid Flows with Complex Free Surfaces Andrea Bonito 1∗ and Alexandre Caboussat 2 , Marco Picasso 3 , and Jacques Rappaz 3 1 Department of Mathematics, University of Maryland, College Park, MD 20742-4015, USA andrea.bonito@a3.epfl.ch 2 Department of Mathematics, University of Houston, 77204-3008, Houston, TX, USA caboussat@math.uh.edu 3 Institute of Analysis and Scientific Computing, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland {marco.picasso;jacques.rappaz}@epfl.ch Summary. A numerical method for the simulation of fluid flows with complex free surfaces is presented. The liquid is assumed to be a Newtonian or a viscoelastic fluid. The compressible effects of the surrounding gas are taken into account, as well as surface tension forces. An Eulerian approach based on the volume-of-fluid formulation is chosen. A time splitting algorithm, together with a two-grids method, allows the various physical phenomena to be decoupled. A chronological approach is adopted to highlight the successive improvements of the model and the wide range of applications. Numerical results show the potentialities of the method. 1 Introduction Complex free surface phenomena involving Newtonian and/or non-Newtonian flows are nowadays a topic of active research in many fields of physics, engineering or bioengineering. The literature contains numerous models for complex liquid-gas free surfaces problems, see, e.g., [FCD + 06, SZ99]. For instance, when considering the injection of a liquid in a complex cavity initially filled with gas, an Eulerian approach is generally adopted in order to catch the topology changes of the liquid region. Such two-phases flows are computationally expensive in three space dimensions since (at least) both the velocity and pressure must be computed at each grid point of the whole liquid-gas domain. The purpose of this article is to review a numerical model in order to compute complex free surface flows in three space dimensions. The features ∗ Partially supported by the Swiss National Science Foundation Fellowship PBEL2– 114311
188 A. Bonito et al. of the model are the following. A volume-of-fluid method is used to track the liquid domain, which can exhibit complex topology changes. The velocity field is computed only in the liquid region. The incompressible liquid can be modeled either as a Newtonian or as a viscoelastic fluid. The ideal gas law is used to compute the external pressure in the surrounding gas and the resulting force is added on the liquid-gas free surface. Surface tension effects can also be taken into account on the liquid-gas free surface. The complete description of the model can be found in [MPR99, MPR03, CPR05, Cab06, BPL06]. The numerical model is based on a time-splitting approach [Glo03] and a two-grids method. This allows advection, diffusion and viscoelastic phenomena to be decoupled, as well as the treatment of the liquid and gas phases. Finite element techniques [FF92] are used to solve the diffusion phenomena using an unstructured mesh of the cavity containing the liquid. A forward characteristic method [Pir89] on a structured grid allows advection phenomena to be solved efficiently. The article is structured as follows. In Section 2, the simplest model is presented: the liquid is an incompressible Newtonian fluid, the effects of the surrounding gas and surface tension are neglected. The effects of the surrounding gas are described in Section 3, those of the surface tension in Section 4. Finally, the case of a viscoelastic liquid is considered in Section 5. Numerical results are presented throughout the text and illustrate the capabilities and improvements of the model. 2 Modeling of an Incompressible Newtonian Fluid with a Free Surface 2.1 Governing Equations The model presented in this section has already been published in [MPR99, MPR03]. Let Λ, with a boundary ∂Λ, be a cavity of R 3 in which a liquid must be confined, and let T>0 be the final time of simulation. For any given time t ∈ (0,T), let Ω t , with a boundary ∂Ω t , be the domain occupied by the liquid, let Γ t = ∂Ω t \ ∂Λ be the free surface between the liquid and the surrounding gas and let Q T be the space-time domain containing the liquid, i.e. Q T = {(x, t) :x ∈ Ω t , 0
- Page 137 and 138: 132 S. Lapin et al. Ω R γ Ω 2 Γ
- Page 139 and 140: 134 S. Lapin et al. ∫ ∂ 2 ∫
- Page 141 and 142: 136 S. Lapin et al. Ω R γ Fig. 3.
- Page 143 and 144: 138 S. Lapin et al. 4 Energy Inequa
- Page 145 and 146: 140 S. Lapin et al. 5 Numerical Exp
- Page 147 and 148: 142 S. Lapin et al. Fig. 6. Contour
- Page 149 and 150: 144 S. Lapin et al. Fig. 9. Obstacl
- Page 151 and 152: Domain Decomposition and Electronic
- Page 153 and 154: Domain Decomposition Approach for C
- Page 155 and 156: Domain Decomposition Approach for C
- Page 157 and 158: Domain Decomposition Approach for C
- Page 159 and 160: Domain Decomposition Approach for C
- Page 161 and 162: Domain Decomposition Approach for C
- Page 163 and 164: Domain Decomposition Approach for C
- Page 165 and 166: Domain Decomposition Approach for C
- Page 167 and 168: Domain Decomposition Approach for C
- Page 169 and 170: Numerical Analysis of a Finite Elem
- Page 171 and 172: Numerical Analysis of a Finite Elem
- Page 173 and 174: Numerical Analysis of a Finite Elem
- Page 175 and 176: Numerical Analysis of a Finite Elem
- Page 177 and 178: Numerical Analysis of a Finite Elem
- Page 179 and 180: Numerical Analysis of a Finite Elem
- Page 181 and 182: so that |u| 2 1,ω h ≤ 2 Numerica
- Page 183 and 184: Numerical Analysis of a Finite Elem
- Page 185 and 186: Numerical Analysis of a Finite Elem
- Page 187: Numerical Analysis of a Finite Elem
- Page 191 and 192: 190 A. Bonito et al. Fig. 1. The sp
- Page 193 and 194: 192 A. Bonito et al. 3 1 16 4 1 1 1
- Page 195 and 196: 194 A. Bonito et al. ∫ v n+1 h
- Page 197 and 198: 196 A. Bonito et al. −pn +2µD(v)
- Page 199 and 200: 198 A. Bonito et al. each of its pa
- Page 201 and 202: 200 A. Bonito et al. The normal vec
- Page 203 and 204: 202 A. Bonito et al. with initial c
- Page 205 and 206: 204 A. Bonito et al. Fig. 8. Jet bu
- Page 207 and 208: 206 A. Bonito et al. References [AM
- Page 209 and 210: 208 A. Bonito et al. [Set96] J. A.
- Page 211 and 212: 210 J. Hao et al. due to shear flow
- Page 213 and 214: 212 J. Hao et al. The backward reac
- Page 215 and 216: 214 J. Hao et al. where u and p den
- Page 217 and 218: 216 J. Hao et al. * * * * * * * * *
- Page 219 and 220: 218 J. Hao et al. and solve for V n
- Page 221 and 222: 220 J. Hao et al. Table 2. The calc
- Page 223 and 224: 222 J. Hao et al. References [ASS80
- Page 225 and 226: Computing the Eigenvalues of the La
- Page 227 and 228: Eigenvalues of the Laplace-Beltrami
- Page 229 and 230: Eigenvalues of the Laplace-Beltrami
- Page 231 and 232: Eigenvalues of the Laplace-Beltrami
- Page 233 and 234: A Fixed Domain Approach in Shape Op
- Page 235 and 236: Shape Optimization Problems with Ne
- Page 237 and 238: Shape Optimization Problems with Ne
188 A. Bonito et al.<br />
of the model are the following. A volume-of-fluid method is used to track<br />
the liquid domain, which can exhibit complex topology changes. The velocity<br />
field is computed only in the liquid region. The incompressible liquid can be<br />
modeled either as a Newtonian or as a viscoelastic fluid. The ideal gas law is<br />
used to compute the external pressure in the surrounding gas <strong>and</strong> the resulting<br />
force is added on the liquid-gas free surface. Surface tension effects can also<br />
be taken into account on the liquid-gas free surface. The complete description<br />
of the model can be found in [MPR99, MPR03, CPR05, Cab06, BPL06].<br />
The numerical model is based on a time-splitting approach [Glo03] <strong>and</strong> a<br />
two-grids method. This allows advection, diffusion <strong>and</strong> viscoelastic phenomena<br />
to be decoupled, as well as the treatment of the liquid <strong>and</strong> gas phases. Finite<br />
element techniques [FF92] are used to solve the diffusion phenomena using an<br />
unstructured mesh of the cavity containing the liquid. A forward characteristic<br />
method [Pir89] on a structured grid allows advection phenomena to be solved<br />
efficiently.<br />
The article is structured as follows. In Section 2, the simplest model is<br />
presented: the liquid is an incompressible Newtonian fluid, the effects of the<br />
surrounding gas <strong>and</strong> surface tension are neglected. The effects of the surrounding<br />
gas are described in Section 3, those of the surface tension in Section 4.<br />
Finally, the case of a viscoelastic liquid is considered in Section 5. Numerical<br />
results are presented throughout the text <strong>and</strong> illustrate the capabilities <strong>and</strong><br />
improvements of the model.<br />
2 Modeling of an Incompressible Newtonian Fluid<br />
with a Free Surface<br />
2.1 Governing <strong>Equations</strong><br />
The model presented in this section has already been published in [MPR99,<br />
MPR03]. Let Λ, with a boundary ∂Λ, be a cavity of R 3 in which a liquid<br />
must be confined, <strong>and</strong> let T>0 be the final time of simulation. For any given<br />
time t ∈ (0,T), let Ω t , with a boundary ∂Ω t , be the domain occupied by<br />
the liquid, let Γ t = ∂Ω t \ ∂Λ be the free surface between the liquid <strong>and</strong> the<br />
surrounding gas <strong>and</strong> let Q T be the space-time domain containing the liquid,<br />
i.e. Q T = {(x, t) :x ∈ Ω t , 0