Abstract
We introduce a new regularization model for incompressible fluid flow, which is a regularization of the EMAC (energy, momentum, and angular momentum conserving) formulation of the Navier–Stokes equations (NSE) that we call EMAC-Reg. The EMAC formulation has proved to be a useful formulation because it conserves energy, momentum, and angular momentum even when the divergence constraint is only weakly enforced. However, it is still a NSE formulation and so cannot resolve higher Reynolds number flows without very fine meshes. By carefully introducing regularization into the EMAC formulation, we create a model more suitable for coarser mesh computations but that still conserves the same quantities as EMAC, i.e., energy, momentum, and angular momentum. We show that EMAC-Reg, when semi-discretized with a finite element spatial discretization is well-posed and optimally accurate. Numerical results are provided that show EMAC-Reg is a robust coarse mesh model.
Funding statement: The work was partially supported by NSF grant DMS 2011490.
Acknowledgment
ClemsonUniversity is acknowledged for allowing a generous amount of computation time on the Palmetto cluster.
References
[1] R. V. Abramov and A. J.Majda, Discrete approximations with additional conserved quantities: deterministic and statistical behavior, Methods Appl. Anal. 10 (2003), No. 2, 151–190.10.4310/MAA.2003.v10.n2.a1Search in Google Scholar
[2] N. Adams and S. Stolz, A subgrid-scale deconvolution approach for shock capturing, J. Comput. Phys. (2002).10.1006/jcph.2002.7034Search in Google Scholar
[3] N. A. Adams and S. Stolz, On the approximate deconvolution procedure for LES, Physics of Fluids 2 (1999), 1699–1701.10.1063/1.869867Search in Google Scholar
[4] A. Arakawa, Computational design for long-term numerical integration of the equations of fluid motion: Two dimensional incompressible flow, Part I, J. Comput. Phys. 1 (1966), 119–143.10.1016/0021-9991(66)90015-5Search in Google Scholar
[5] D. Arndt, W. Bangerth, B. Blais, M. Fehling, R. Gassmöller, T. Heister, et al., The deal. II library, Version 9.3, J. Numer. Math. 29 (2021), No. 3, 171–186.10.1515/jnma-2021-0081Search in Google Scholar
[6] D. Arnold and J. Qin, Quadratic velocity/linear pressure Stokes elements, In: Advances in Computer Methods for Partial Differential Equations VII (Eds. R. Vichnevetsky, D. Knight, and G. Richter), IMACS, 1992, pp. 28–34.Search in Google Scholar
[7] L. C. Berselli, T. Iliescu, and W. J. Layton, Mathematics of Large Eddy Simulation of Turbulent Flows, Scientific Computation, Springer-Verlag, Berlin, 2006.Search in Google Scholar
[8] A. Bowers, T. Kim, M. Neda, L. Rebholz, and E. Fried, The Leray-αβ-deconvolution model: energy analysis and numerical algorithms, Appl. Math. Modelling 37 (2013), No. 3, 1225–1241.10.1016/j.apm.2012.03.040Search in Google Scholar
[9] S. Breckling, M. Neda, and F. Pahlevani, A sensitivity study of the Navier–Stokes-α model, Computers & Math. Appl. 75 (2018), No. 2, 666–689.10.1016/j.camwa.2017.09.036Search in Google Scholar
[10] T. Chacon-Rebello and R. Lewandowski, Mathematical and Numerical Foundations of Turbulence Models and Applications, Springer, New York, 2014.10.1007/978-1-4939-0455-6Search in Google Scholar
[11] S. Charnyi, T. Heister, M.Olshanskii, and L. Rebholz, On conservation laws of Navier–Stokes Galerkin discretizations, J. Comput. Phys. 337 (2017), 289–308.10.1016/j.jcp.2017.02.039Search in Google Scholar
[12] S. Charnyi, T. Heister, M.Olshanskii, and L. Rebholz, Eflcient discretizations for the EMAC formulation of the incompressible Navier–Stokes equations, Appl. Numer. Math. 141 (2019), 220–233.10.1016/j.apnum.2018.11.013Search in Google Scholar
[13] S. Chen, C. Foias, D. D. Holm, E.Olson, E. S. Titi, and S.Wynne, The Camassa–Holm equations as a closure model for turbulent channel and pipe flow, Phys. Review Letters 81 (1998), 5338–5341.10.1103/PhysRevLett.81.5338Search in Google Scholar
[14] S. Chen, C. Foias, E.Olson, E. S. Titi, and W.Wynne, The Camassa–Holm equations and turbulence, Physica D 133 (1999), 49–65.10.1016/S0167-2789(99)00098-6Search in Google Scholar
[15] A. Cheskidov, D. Holm, E.Olson, and E. S. Titi, On a Leray-α model of turbulence, Proceedings of The Royal Society A 461 (2005), 629–649.10.1098/rspa.2004.1373Search in Google Scholar
[16] A. J. Chorin, Numerical solution for the Navier–Stokes equations, Math. Comp. 22 (1968), 745–762.10.1090/S0025-5718-1968-0242392-2Search in Google Scholar
[17] J. Connors, Convergence analysis and computational testing of the finite element discretization of the Navier–Stokes α model, Numer. Meth. Partial Differ. Equ. 26 (2010), No. 6, 1328–1350.10.1002/num.20493Search in Google Scholar
[18] V. Cuff, A. Dunca, C.Manica, and L. Rebholz, The reduced order NS-α model for incompressible flow: theory, numerical analysis and benchmark testing, ESAIM: Mathematical Modelling and Numerical Analysis 49 (2015), No. 3, 641–662.10.1051/m2an/2014053Search in Google Scholar
[19] A. Dunca and Y. Epshteyn, On the Stolz–Adams deconvolution model for the large-eddy simulation of turbulent flows, SIAM J. Math. Anal. 37 (2005), No. 6, 1890–1902.10.1137/S0036141003436302Search in Google Scholar
[20] J. A. Evans and T. J. R. Hughes, Isogeometric divergence-conforming B-splines for the steady Navier–Stokes equations, Math. Models Methods Appl. Sci. 23 (2013), No. 08, 1421–1478.10.1142/S0218202513500139Search in Google Scholar
[21] J. A. Evans and T. J. R. Hughes, Isogeometric divergence-conforming B-splines for the unsteady Navier–Stokes equations, J. Comput. Phys. 241 (2013), 141–167.10.21236/ADA560939Search in Google Scholar
[22] O. Lehmkuhl F. Sacco, B. Paun, T. Iles, P. Iaizzo, G. Houzeaux, M. Vzaquex, C. Butakoff, and J. Aguado-Sierra, Evaluating the roles of detailed endocardial structures on right ventricular haemodynamics by means of CFD simulations, Int. J. Numer. Methods Biomed. Engrg. 34 (2018), 1–14.10.1002/cnm.3115Search in Google Scholar
[23] O. Lehmkuhl F. Sacco, B. Paun, T. Iles, P. Iaizzo, G. Houzeaux, M. Vzaquex, C. Butakoff, and J. Aguado-Sierra, Left ventricular trabeculations decrease the wall shear stress and increase the intra-ventricular pressure drop in CFD simulations, Frontiers Physiol. 9 (2018), 1–15.10.3389/fphys.2018.00458Search in Google Scholar
[24] G. Fix, Finite element models for Ocean circulation problems, SIAM J. Appl. Math. 29 (1975), No. 3, 371–387.10.1137/0129031Search in Google Scholar
[25] C. Foias, D. D. Holm, and E. S. Titi, The Navier–Stokes-α model of fluid turbulence, Physica D 152 (2001), 505–519.10.1016/S0167-2789(01)00191-9Search in Google Scholar
[26] C. Foias, D. D. Holm, and E. S. Titi, The three dimensional viscous Camassa–Holm equations, and their relation to the Navier–Stokes equations and turbulence theory, J. Dynam. Differ. Equ. 14 (2002), 1–35.10.1023/A:1012984210582Search in Google Scholar
[27] U. Frisch, Turbulence, Cambridge University Press, 1995.10.1017/CBO9781139170666Search in Google Scholar
[28] B. J. Geurts and D. Holm, Regularization modeling for large eddy simulation, Phys. Fluids 15 (2003), L13.10.1063/1.1529180Search in Google Scholar
[29] B. J. Geurts and D. Holm, Leray and LANS-α modeling of turbulent mixing, J. Turbulence 7 (2006), 1–33.10.1080/14685240500307389Search in Google Scholar
[30] J. Guzman and M. Neilan, Conforming and divergence-free Stokes elements in three dimensions, IMA J. Numer. Anal. 34 (2014), No. 4, 1489–1508.10.1093/imanum/drt053Search in Google Scholar
[31] J. Guzman and M. Neilan, Conforming and divergence-free Stokes elements on general triangular meshes, Math. Comput. 83 (2014), 15–36.10.1090/S0025-5718-2013-02753-6Search in Google Scholar
[32] D. Holm and B. T. Nadiga, Modeling mesoscale turbulence in the barotropic double-gyre circulation, J. Phys. Oceanography 33 (2003), 2355–2365.10.1175/1520-0485(2003)033<2355:MMTITB>2.0.CO;2Search in Google Scholar
[33] A. N. Kolmogorov, Dissipation of energy in isotropic turbulence, Dokl. Akad. Nauk. SSSR 32 (1941), 19–21.Search in Google Scholar
[34] A. N. Kolmogorov, The local structure of turbulence in incompressible viscous fluids at very large Reynolds numbers, Dokl. Akad. Nauk. SSSR 30 (1941), 299–303.Search in Google Scholar
[35] A. N. Kolmogorov, On the degeneration of isotropic turbulence in an incompressible viscous fluids, Dokl. Akad. Nauk. SSSR 31 (1941), 538–541.Search in Google Scholar
[36] W. Layton, Introduction to the Numerical Analysis of Incompressible Viscous Flows, SIAM, 2008.10.1137/1.9780898718904Search in Google Scholar
[37] W. Layton, C.Manica, M. Neda, and L. Rebholz, Numerical analysis and computational testing of a high accuracy Leray-deconvolution model of turbulence, Numer. Meth. Partial Differ. Equ. 24 (2008), No. 2, 555–582.10.1002/num.20281Search in Google Scholar
[38] W. Layton and L. Rebholz, Approximate Deconvolution Models of Turbulence: Analysis, Phenomenology and Numerical Analysis, Springer-Verlag, 2012.10.1007/978-3-642-24409-4Search in Google Scholar
[39] O. Lehmkuhl, U. Piomelli, and G. Houzeaux, On the extension of the integral length-scale approximation model to complex geometries, Int. J. Heat Fluid Flow 78 (2019), No. 108422, 1–12.10.1016/j.ijheatfluidflow.2019.108422Search in Google Scholar
[40] J. Liu and W.Wang, Energy and helicity preserving schemes for hydro and magnetohydro-dynamics flows with symmetry, J. Comput. Phys. 200 (2004), 8–33.10.1016/j.jcp.2004.03.005Search in Google Scholar
[41] R.Martin, M. Soria, O. Lehmkuhl, A. Gorobets, and A. Duben, Noise radiated by an open cavity at low Mach number: Effect of the cavity oscillation mode, Int. J. Aeroacoustics 18 (2019), No. 6-7, 647–668.10.1177/1475472X19871534Search in Google Scholar
[42] M.Olshanskii and L. Rebholz, Longer time accuracy for incompressible Navier–Stokes simulations with the EMAC formulation, Comp. Methods Appl. Mech. Engrg. 372 (2020), 113369.10.1016/j.cma.2020.113369Search in Google Scholar
[43] H. Owen, G. Chrysokentis, M. Avila, D. Mira, G. Houzeaux, R. Borrell, J. C. Cajas, and O. Lehmkuhl, Wall-modeled large-eddy simulation in a finite element framework, Int. J. Numer. Methods Fluids 92 (2020), No. 1, 20–37.10.1002/fld.4770Search in Google Scholar
[44] A. Palha and M. Gerritsma, A mass, energy, enstrophy and vorticity conserving (MEEVC) mimetic spectral element discretization for the 2D incompressible Navier–Stokes equations, J. Comput. Phys. 328 (2017), 200–220.10.1016/j.jcp.2016.10.009Search in Google Scholar
[45] D. Pastrana, J. C. Cajas, O. Lehmkuhl, I. Rodríguez, and G. Houzeaux, Large-eddy simulations of the vortex-induced vibration of a low mass ratio two-degree-of-freedom circular cylinder at subcritical Reynolds numbers, Comput. Fluids 173 (2018), 118–132.10.1016/j.compfluid.2018.03.016Search in Google Scholar
[46] S. Pope, Turbulent Flows, Cambridge University Press, 2000.10.1017/CBO9780511840531Search in Google Scholar
[47] L. Rebholz, An energy and helicity conserving finite element scheme for the Navier–Stokes equations, SIAM J. Numer. Anal. 45 (2007), No. 4, 1622–1638.10.1137/060651227Search in Google Scholar
[48] L. Rebholz, Conservation laws of turbulence models, J. Math. Anal. Appl. 326 (2007), No. 1, 33–45.10.1016/j.jmaa.2006.02.026Search in Google Scholar
[49] L. Rebholz, T.-Y. Kim, and Y.-L. Byon, On an accurate α model for coarse mesh turbulent channel flow simulation, Appl. Math. Modelling 43 (2017), 139–154.10.1016/j.apm.2016.10.059Search in Google Scholar
[50] T. C. Rebollo and R. Lewandowski, Mathematical and Numerical Foundations of Turbulence Models and Applications, Springer, 2014.Search in Google Scholar
[51] R. Salmon and L. D. Talley, Generalizations of Arakawa’s Jacobian, J. Comput. Phys. 83 (1989), 247–259.10.1016/0021-9991(89)90118-6Search in Google Scholar
[52] P. Schroeder, V. John, P. Lederer, C. Lehrenfeld, G. Lube, and J. Schoberl, On reference solutions and the sensitivity of the 2D Kelvin–Helmholtz instability problem, Comput. Math. Appl. 77 (2019), No. 4, 1010–1028.10.1016/j.camwa.2018.10.030Search in Google Scholar
[53] P. Schroeder and G. Lube, Pressure-robust analysis of divergence-free and conforming FEM for evolutionary incompressible Navier–Stokes flows, J. Numer. Math. 25 (2017), No. 4, 249–276.10.1515/jnma-2016-1101Search in Google Scholar
[54] J. S.Smagorinsky, General circulation experiments with the primitive equations, Mon. Weather Rev. 91 (1963), 99–164.10.1175/1520-0493(1963)091<0099:GCEWTP>2.3.CO;2Search in Google Scholar
[55] C. Sorgentone, S. La Cognata, and J. Nordstrom, A new high order energy and enstrophy conserving Arakawa-like Jacobian differential operator, J. Comput. Physics 301 (2015), 167–177.10.1016/j.jcp.2015.08.028Search in Google Scholar
[56] S. Stolz, N. Adams, and L. Kleiser, The approximate deconvolution model for large-eddy simulations of compressible flows and its application to shock-turbulent-boundary-layer interaction, Phys. Fluids 13 (2001), No. 10, 2985–3001.10.1063/1.1397277Search in Google Scholar
[57] X. Xie, D. Wells, Z.Wang, and T. Iliescu, Numerical analysis of the Leray reduced order model, J. Comput. Appl. Math. 328 (2018), 12–29.10.1016/j.cam.2017.06.026Search in Google Scholar
[58] S. Zhang, A new family of stable mixed finite elements for the 3D Stokes equations, Math. Comput. 74 (2005), 543–554.10.1090/S0025-5718-04-01711-9Search in Google Scholar
A Momentum/angular momentum conservation of NS-α and Leray-α formulations
Here we show that the NS-α and Leray-α formulations do not conserve momentum or angular momentum if div u ≠ 0 and div w ≠ 0 where w represents the filtered velocity ū. Further note u = Fw where F = −α2ΔI + I.
A.1 NS-α
Recall the nonlinear term of the NS-α formulation is
Test (A.1) with ei for i = 1, 2, 3. After applying the space–time divergence theorem and rearranging some we get
Assuming ν = f = 0, (A.2) simplifies into
If the nonlinear term is equal to zero, then we will have momentum conservation. We now check this:
where the above two equalities come from vector identities. Also note that because ei is a vector of scalars, ((∇ ⋅ ei)w, u) = (w ⋅ ∇ei , u) = 0. This leaves us with
which we cannot conclude is zero, hence we cannot say that the NS-α formulation preserves momentum. For angular momentum, we test (A.1) with φi and the algebra works out similar to momentum,
Since ∇⋅ φi = 0 for i = 1, 2, 3, we have
Also recall using (3.15) in Theorem 3.3, we have
This gives us
Much like with momentum, we cannot conclude that this quantity is zero, and we expect it is not zero.
A.2 Leray-α
Recall the nonlinear term of the Leray-α formulation is
We test (A.3) with ei for i = 1, 2, 3 and integrate. Similar to (A.2)
Assuming ν = f = 0, (A.4) simplifies to
If the nonlinear term is equal to zero, thenwe will have momentum conservation. Using (2.6) on the nonlinear term we get
which is not zero when ∇⋅ w ≠ 0. Hence momentum is not necessarily conserved.
For angular momentum we test (A.3) with φi for i = 1, 2, 3 and it simplifies to
Now similarly to the momentum proof, we have for the nonlinear term
where the first term disappears by applying (3.15) similarly to Theorem 3.3 (and the angular momentum proof for NS-α in Appendix A.1). Thus the nonlinear term does not vanish, so angular momentum is not conserved.
© 2022 Walter de Gruyter GmbH, Berlin/Boston