Next Article in Journal
Quantized Graph Neural Networks for Image Classification
Previous Article in Journal
Robust and Sparse Portfolio: Optimization Models and Algorithms
Previous Article in Special Issue
Exploring the Influence of Induced Magnetic Fields and Double-Diffusive Convection on Carreau Nanofluid Flow through Diverse Geometries: A Comparative Study Using Numerical and ANN Approaches
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Mathematically Exact and Well-Determined System of Equations to Close Reynolds-Averaged Navier–Stokes Equations

Department of Mechanical Engineering, Incheon National University, Academy-ro 119, Incheon 22012, Republic of Korea
Mathematics 2023, 11(24), 4926; https://doi.org/10.3390/math11244926
Submission received: 12 November 2023 / Revised: 7 December 2023 / Accepted: 8 December 2023 / Published: 11 December 2023

Abstract

:
Since Sir Osborne Reynolds presented the Reynolds-averaged Navier–Stokes (RANS) equations in 1895, the construction of complete closure for RANS equations has been regarded as extremely challenging. Taking into account that the Navier–Stokes equations are not coherent for instantaneous and mean flows, a body of knowledge outside the scope of classical mechanics may be amenable to the closure problem. In this regard, the methodology of physics-to-geometry transformation, which is coherent for both flows, is applied to RANS equations to construct six additional equations. The proposed equations stand out from existing RANS closure models and turbulence quantity transport equations in two respects: they are mathematically exact and well-determined.

1. Introduction

Since the Reynolds-averaged Navier–Stokes (RANS) equations were first introduced [1], constructing a complete closure for them has remained an open problem. Causal relations over time between mean momenta and stresses in turbulent flow are not generally describable in terms of Navier–Stokes equations, so these are not coherent for instantaneous and mean flows:
N . N . ,
where N is the Navier–Stokes operator and . denotes a time- or ensemble-averaging operator. This suggests that mean flow is fundamentally different from instantaneous flow in the scope of classical mechanics, which is a plausible reason why the physical concepts and laws established for instantaneous flows cannot directly apply to mean flows. The notorious difficulty of the closure problem may originate from this difference.
Mathematical approaches make it possible to understand mean flows via instantaneous flow physics, despite the above-described incoherence. For linear flow quantities and operators, the numerical equality derived by the use of the Reynolds decomposition and the fundamental properties of Reynolds-averaging allows for mathematical mimics; examples of this include the mean rate of strain and rotation tensors. For nonlinear ones, such a mimic is not possible since additional terms that require physical analysis arise from the process of Reynolds averaging. In particular, for mathematical approaches to construct a RANS closure, such an appearance of additional unclosed terms is definitely an obstacle even though numerical equality is obtainable. On the other hand, it is also not easy to mathematically derive an exact relation between the Reynolds stress tensor and mean flow quantities without additional unclosed terms. This is a common dilemma encountered whilst constructing a RANS closure.
Modeling has been widely adopted as an alternative methodology since it allows for an explicit formulation of Reynolds stress tensor and is therefore tractable for numerical simulation. However, it resorts to approaches in which mathematical exactness is not ensured a priori, such as physical intuition and empiricism, which inevitably lead to the issue of universality. Examples include the linear stress–strain relation with the eddy-viscosity, which was generalized from the Boussinesq hypothesis [2], and the mixing length theory [3]. In particular, for the former, the isotropy eddy-viscosity assumption has represented a key concept for constructing RANS closure models, but its scope of application is limited; it is strictly correct under the condition of τ S , where τ is the turbulent time scale and S is the mean shear time scale [4], while otherwise the linear relation between the Reynolds stress tensor and the mean rate-of-strain tensor is violated in anisotropic turbulent flows (see, e.g., Ref. [5]). Correspondingly, the following efforts were focused on the development of a RANS closure model to obtain high predictive capability for anisotropic turbulent flows. Actually, advanced RANS closure models were devised by incorporating anisotropy into an eddy-viscosity form or considering higher-order stress–strain relations; examples of this include quadratic [6] and cubic k ϵ models [7], V2F model [8], and quadratic constitutive relation-based models [9,10]. Recently, it was demonstrated that the above-described shortcoming of linear eddy-viscosity models can be remedied by adopting data processing techniques into existing RANS closure models (see, e.g., Refs. [11,12,13]). For a comprehensive review of the related literature on data-driven RANS models, see the review by Duraisamy et al. [14]. Despite the continuous development of RANS closure models to date, none of them ensure mathematical exactness a priori. Therefore, their scope of application is not deterministic, but inevitably open to an empirical study.
Mathematical exactness is a necessary condition for a complete closure of RANS equations, but it is definitely not sufficient. In practice, the mathematically exact relations associated with the Reynolds stress tensor and turbulence quantity transport equations, which have been reported to date, suffer due to the appearance of additional unclosed terms. The needs for modeling were thus inevitably evoked to resolve this under-determinedness. A representative example is the Reynolds stress transport equations, which have unclosed source terms with the exception of the production term. The under-determinedness of these equations seems to have motivated the development of algebraic stress models. Rodi [15] constructed an implicit algebraic stress model by applying an equilibrium hypothesis to the modeled Reynolds stress transport equation of Launder et al. [16]. Since this implicit model requires numerical iterations, it could not be computationally robust [17]. Pope [18] derived an explicit relation for the Reynolds stress tensor by combining the tensorial polynomial expansion thereof with the modeled Reynolds stress transport equation of Rodi [15], but the proposed model is limited to two-dimensional turbulent flow. Gatski and Speziale [19] extended the approach of Pope [18] to three-dimensional turbulent flow by fully considering the ten integrity bases. Further developments were made by introducing near-wall treatments [20] and considering streamline curvature effects [21]. In light of the fundamental reason for the emergence of these algebraic stress models, one may know that the non-appearance of additional unclosed terms is also a necessary condition for the construction of a complete RANS closure.
Most recently, Ryu [22] constructed three identities between Reynolds stresses via spatial mapping and application of the differential Gauss–Bonnet formula [23]. Their methodology is noticeable in two respects: mathematical exactness is ensured by the mathematically proven numerical equality, and the appearance of additional unclosed terms can be prevented by a selective choice of the unit conversion function. However, above all, it is only concerned to the analyticity of a physical scalar field in R 3 , which suggests that the methodology can be coherent for both instantaneous and mean flow fields.
In the present paper, the objective is to construct a mathematically exact and well-determined RANS closure. For this task, the incompressible RANS equations are developed into six additional ones by means of the mathematical framework of Ryu [22].

2. Physics-to-Geometry Transformation

We elaborate on the construction process of the mathematical framework of Ryu [22] to help the reader better understand our application of it to RANS equations. Their mathematical framework is constructed in two steps: (i) a physical scalar field of interest is mapped into R 3 and (ii) the differential version of the Gauss–Bonnet formula is applied to the coordinate plane subfields of the mapped field.
In regard to the former, it is the key to convert the unit of the associated physical quantity into that of length. Hereafter, let G ( t , x ) denote a time-dependent, physical scalar field and D R 3 be a spatial domain for G ( t , x ) , where t is time and x = ( x 1 , x 2 , x 3 ) are the Cartesian coordinates. The unit of G ( t , x ) may be converted into that of length, as follows:
G ˜ ( t , x ) = G ( t , x ) σ ( t , x ) ,
where σ is the unit-conversion function. Hereafter, we use the overtilde sign ( . ) ˜ to indicate such a converted field. Then, the respective three subfields of G ˜ ( t , x ) are defined on the three coordinate planes intersecting at each point of D. For example, the converted field G ˜ has three coordinate plane subfields at time t o and p o ( x o , y o , z o ) D : G ˜ ( t o , x , y , z o ) , G ˜ ( t o , x , y o , z ) , G ˜ ( t o , x o , y , z ) .
The differential version of the Gauss–Bonnet formula, which has been mathematically proven in Ref. [23], yields three differential equations at p o for the above-defined three graphs, respectively. The geometric entities of each of those equations are then expanded as introduced in Appendix A. Consequently, this affords three identities in terms only of G ˜ at p o :
G ˜ x x G ˜ y y G ˜ x y 2 1 + G ˜ x 2 + G ˜ y 2 3 / 2 + x G ˜ y y G ˜ x 1 + G ˜ y 2 1 + G ˜ x 2 + G ˜ y 2 y G ˜ x x G ˜ y 1 + G ˜ x 2 1 + G ˜ x 2 + G ˜ y 2                                                                                                                                                                       + 2 x y arccos G ˜ x G ˜ y 1 + G ˜ x 2 1 + G ˜ y 2 = 0 ,
G ˜ z z G ˜ x x G ˜ z x 2 1 + G ˜ z 2 + G ˜ x 2 3 / 2 + z G ˜ x x G ˜ z 1 + G ˜ x 2 1 + G ˜ z 2 + G ˜ x 2 x G ˜ z z G ˜ x 1 + G ˜ z 2 1 + G ˜ z 2 + G ˜ x 2                                                                                                                                                                       + 2 z x arccos G ˜ z G ˜ x 1 + G ˜ z 2 1 + G ˜ x 2 = 0 ,
G ˜ y y G ˜ z z G ˜ y z 2 1 + G ˜ y 2 + G ˜ z 2 3 / 2 + y G ˜ z z G ˜ y 1 + G ˜ z 2 1 + G ˜ y 2 + G ˜ z 2 z G ˜ y y G ˜ z 1 + G ˜ y 2 1 + G ˜ y 2 + G ˜ z 2                                                                                                                                                                       + 2 y z arccos G ˜ y G ˜ z 1 + G ˜ y 2 1 + G ˜ z 2 = 0 ,
where the subscripts denote the derivatives of G ˜ with respect to x, y, and z, respectively. Note that application of the differential Gauss–Bonnet formula is first considered only at a point of D. For convenience, the left-hand sides of these identities are briefly expressed as the sum of the three budgets:
G x , y ; G ˜ t o , x , y , z o + Π x , y ; G ˜ t o , x , y , z o + Φ x , y ; G ˜ t o , x , y , z o = 0 ,
G z , x ; G ˜ t o , x , y o , z + Π z , x ; G ˜ t o , x , y o , z + Φ z , x ; G ˜ t o , x , y o , z = 0 ,
G y , z ; G ˜ t o , x o , y , z + Π y , z ; G ˜ t o , x o , y , z + Φ y , z ; G ˜ t o , x o , y , z = 0 ,
where G , Π , and Φ are the respective differential operators in order. By adding up the operators of the above equations, respectively, we obtain
G s x ; G ˜ | ( t o , p o ) + Π s x ; G ˜ | ( t o , p o ) + Φ s x ; G ˜ | ( t o , p o ) = 0 .
Extending the above-described application of the differential Gauss–Bonnet formula to every point of D makes this combined identity applicable as a whole on D:
G s x ; G ˜ ( t o , x ) + Π s x ; G ˜ ( t o , x ) + Φ s x ; G ˜ ( t o , x ) = 0 , x D .

3. Invariance Properties

The above-introduced identity of Equation (10) originates from the differential version of the Gauss–Bonnet formula that concerns only the geometry of an orientable smooth surface in R 3 . This implies intriguing invariance properties, which are not dealt with in Ryu [22]. Hereafter, let R s denote the sum of G s , Π s , and Φ s .
Lemma 1. 
If the coordinate plane subfields of G ˜ are orientable and smooth in R 3 at every instant in time, then the combined identity of Equation (10) is invariant under a shift in time:
R s ( x ; G ˜ ( t , x ) ) = R s ( x ; G ˜ ( t + Δ t , x ) ) .
Lemma 2. 
Let G * ( x ) be a statistic of G ( t , x ) over time and let G * ˜ ( x ) be the converted field of G * ( x ) by a time-independent unit conversion function σ ( x ) . If the coordinate plane subfields of G * ˜ are orientable and smooth in R 3 , then the combined identity of Equation (10) is indifferent to taking a statistic over time:
R s ( x ; G ˜ ( t , x ) ) = R s ( x ; G * ˜ ( x ) ) .
Proof. 
These two lemmas are proven in Theorem A1. □
Lemma 3. 
The combined identity of Equation (10) is invariant under Galilean coordinate transformation.
Proof. 
Let x = ( x , y , z ) and t be spatial coordinates and time in the moving frame of reference S with constant velocity V along the x-direction. Then, the coordinates of x and time t are written in terms of the coordinates of x = ( x , y , z ) and time t in stationary frame of reference S, as follows:
x = x V t , y = y , z = z , t = t .
These relations give the respective expressions for x, y, and z. In substituting those into Equation (10) and manipulating the derivatives therein, the differential operator R s is form-invariant. □

4. Mathematical Formulation

4.1. Generic Formulation

Prior to the construction of a system of equations as a RANS closure, we devise a generic method to derive six additional equations from an established one. To elaborate on this method, first let A = B be an equation with rank r 3 tensor fields A ( t , x ) and B ( t , x ) in R 3 .
First, two scalar quantities A and B are defined consistently from A and B , respectively. For example, in the case of a vector, its magnitude may be chosen as a scalar quantity, and, in the case of a rank-2 tensor, its invariants may afford the scalar quantities. The original equation is then developed into A = B .
Second, A and B are converted into a quantity of length by the same unit conversion:
A ˜ = A ( t , x ) σ r ( t , x ) ,
B ˜ = B ( t , x ) σ r ( t , x ) .
Note that the numerical equality of the developed equation that A = B is retained due to the same unit conversion function:
A ˜ ( t , x ) = B ˜ ( t , x ) .
Since there are many possible options about the determination of a unit conversion function, a selective choice is possible for a specific purpose.
The mathematical framework of Equation (10) gives two identities in terms only of A ˜ and B ˜ , respectively:
G s x ; A ˜ + Π s x ; A ˜ + Φ s x ; A ˜ = 0 ,
G s x ; B ˜ + Π s x ; B ˜ + Φ s x ; B ˜ = 0 .
To derive additional relations between A ˜ and B ˜ , it may be a way to selectively substitute A ˜ for B ˜ and vice versa in Equations (17) and (18), respectively. Note that the numerical equality in Equations (17) and (18) is not affected by these substitutions, owing to A ˜ = B ˜ . In other words, the numerical equality of A ˜ = B ˜ is a necessary condition for the selectively substituted equations to hold true. The number of all possible cases for such a combinatorial substitution is six. The resulting equations are written as follows:
G s x ; B ˜ + Π s x ; A ˜ + Φ s x ; A ˜ = 0 ,
G s x ; A ˜ + Π s x ; B ˜ + Φ s x ; A ˜ = 0 ,
G s x ; A ˜ + Π s x ; A ˜ + Φ s x ; B ˜ = 0 ,
G s x ; A ˜ + Π s x ; B ˜ + Φ s x ; B ˜ = 0 ,
G s x ; B ˜ + Π s x ; A ˜ + Φ s x ; B ˜ = 0 ,
G s x ; B ˜ + Π s x ; B ˜ + Φ s x ; A ˜ = 0 .
The form of these six equations depends on a chosen unit conversion function. Nonetheless, if it consists only of known variables or simply a constant, the appearance of additional unclosed terms can be avoided. In addition, the mathematical exactness of these equations is ensured by the mathematically proven numerical equality of the differential Gauss–Bonnet formula [23], if the coordinate plane subfields of A ˜ and B ˜ are orientable and smooth in R 3 .

4.2. Application to RANS Equations

Taking into consideration that the above-described method gives a total of six additional equations, it is most instructive to apply it to the governing equations for RANS or large-eddy simulation. For the former, the incompressible RANS equations are considered as a target of application in the Cartesian coordinate system:
U i t + U j U i x j = P x i + ν 2 U i x j x j τ i j x j ,
where U i is the averaged velocity component along the x i direction; P is the averaged pressure divided by the fluid density ρ ; ν is the kinematic viscosity; and τ i j is the Reynolds stress tensor. The incompressible RANS equations are then developed into a single one by taking divergence on both sides of Equation (25):
2 τ i j x i x j = 2 P x i 2 U j x i U i x j ,
The left- and right-hand sides of this reduced equation are briefly written as follows:
M ( t , x ) = F ( t , x ) .
Then, M and F are mapped into R 3 by means of the same unit conversion function:
M ˜ ( t , x ) = M ( t , x ) σ r ( t , x ) ,
F ˜ ( t , x ) = F ( t , x ) σ r ( t , x ) .
We choose σ r = 1 to avoid appearances of additional unclosed terms and to make the two conversions simplest. Accordingly,
M ˜ ( t , x ) = F ˜ ( t , x ) .
Then, this equation is developed into six additional ones through the above-described combinatorial substitutions.
The two constraints for the differential Gauss–Bonnet formula, orientability and smoothness, are not strong to over-constrain the system in question. Since a smooth real scalar function has continuous derivatives of all orders [24] and a surface given by the graph of a differentiable real function is orientable [25], any coordinate subfields of a smooth scalar field in R 3 must be orientable. Accordingly, smoothness is the key necessary condition for the mathematical framework of Equations (3)–(5) to hold true. Taking into consideration that smooth flow fields are common in turbulence, it is difficult to see that the physical quantities described by the proposed equations are significantly distorted by the smoothness constraint.
The numerical equality of the incompressible RANS equations is a necessary condition for the proposed six equations to hold true, as mentioned earlier. Further, the mean continuity equation is applied whilst deriving the mean Poisson equation from the incompressible RANS equations. Therefore, mass and momentum conservation laws are involved in the derivation process of the presented six equations. Moreover, those equations are supposed to be coupled with the mean continuity equation and the incompressible RANS equations in order to obtain ten unknown quantities: one averaged pressure, three averaged velocities, and six Reynolds stress components. Therefore, the proposed six equations are not isolated from mass and momentum conservation laws in solving the ten simultaneous equations.

5. Summary and Concluding Remarks

Six equations have been constructed as a closure for the Reynolds stress tensor by utilizing the mathematical framework of Ryu [22]. The target of application in the present study is distinct from their one; the objects of interest are not Reynolds stresses but the system of RANS equations. The proposed RANS closure is mathematically exact and well-determined, in which it stands out from existing RANS closure models and the Reynolds stress transport equations.
It is extremely challenging to analytically solve the incompressible RANS equations in combination with the presented six equations and the mean continuity equation, owing to their implicitness and double-embedded nonlinearity. By contrast, it seems more feasible to numerically solve these simultaneous equations. Regarding the latter, it would be a first step to set physical boundary conditions for Reynolds stresses in order to close the system of equations. Then, we could aim to find a possibility for obtaining numerically stable and accurate solutions for a benchmark turbulent flow. Despite such a barrier to entry, this is a challenge worth taking on as it could be a pioneering attempt to essentially overcome the inherent limitation of modeling.

Funding

This research received no external funding.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. Expansion of Geometric Entity

This Appendix describes the expansions of the geometric entities of the differential Gauss–Bonnet formula for an orientable smooth surface given by parametrization r = ( u , v , f ( u , v ) ) , ( u , v ) R 2 . First, we recall the differential Gauss–Bonnet theorem from Ref. [23].
Theorem A1. 
Let S be an orientable smooth surface in R 3 , and let r : U S be a parametrization of S in an open set U R 2 . Then for each ( u , v ) U
K | N | + F b u F a v 2 ϕ u v = 0 ,
where K is the Gaussian curvature over S, N is the normal to S, F a , and F b are the products of the geodesic curvatures of the coordinate curves v = const and u = const and the speeds of those curves, respectively, and ϕ is the positively oriented angle of intersection from the coordinate curve v = const to u = const on S.
The differential Gauss–Bonnet formula consists of seven geometric entities: K, | N | , κ g , u = const , κ g , v = const , | r u | , | r v | , and ϕ , where κ g is the geodesic curvature over a coordinate curve on the map of r and the subscripts u and v denote the first-order derivatives of r with respect to u and v, respectively. In particular, for a surface given by parametrization ( u , v , f ( u , v ) ) , these are expanded in terms of the derivatives of the two-variable function f ( u , v ) , as follows [25]:
K = f u u f v v f u v 2 1 + f u 2 + f v 2 2 ,
| N | = 1 + f u 2 + f v 2 ,
| r u | = 1 + f u 2 ,
| r v | = 1 + f v 2 ,
κ g , v = const = r u u , n r u | r u | 3 = f u u f v ( 1 + f u 2 ) 3 / 2 1 + f u 2 + f v 2 ,
κ g , u = const = r v v , n r v | r v | 3 = f v v f u ( 1 + f v 2 ) 3 / 2 1 + f u 2 + f v 2 ,
ϕ = arccos r u , r v | r u | | r v | = arccos f u f v 1 + f u 2 1 + f v 2 ,
where n is the unit normal to S. The geometric entities of the differential Gauss–Bonnet formula are similarly expandable for the surface given by parametrization ( g ( u , v ) , u , v ) or ( v , h ( u , v ) , u ) .

References

  1. Reynolds, O. On the dynamical theory of incompressible viscous fluids and the determination of the criterion. Philos. Trans. R. Soc. Lond. Ser. A 1895, 186, 123–164. [Google Scholar]
  2. Boussinesq, J.V. Essai sur la théorie des eaux courantes. In Mémoires Présentés par Divers Savants à I’Académie des Sciences; Imprimerie Nationale: Paris, France, 1877; Volume 23, pp. 1–680. [Google Scholar]
  3. Prandtl, L. 7. Bericht über Untersuchungen zur ausgebildeten Turbulenz. Z. Angew. Math. Mech. 1925, 5, 136–139. [Google Scholar] [CrossRef]
  4. Pope, S.B. Turbulent Flows; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  5. Mompean, G.; Gavrilakis, S.; Machiels, L.; Deville, M.O. On predicting the turbulence-induced secondary flows using nonlinear k-ϵ models. Phys. Fluid 1996, 8, 1856–1868. [Google Scholar] [CrossRef]
  6. Speziale, C.G. On nonlinear k-l and k-ϵ models of turbulence. J. Fluid Mech. 1987, 178, 459–475. [Google Scholar] [CrossRef]
  7. Craft, T.J.; Launder, B.E.; Suga, K. Development and application of a cubic eddy-viscosity model of turbulence. Int. J. Heat Fluid Flow 1996, 17, 108–115. [Google Scholar] [CrossRef]
  8. Durbin, P. Separated flow computations with the kϵv2 model. AIAA J. 1995, 33, 659–664. [Google Scholar] [CrossRef]
  9. Spalart, P.R. Strategies for turbulence modelling and simulations. Int. J. Heat Fluid Flow 2000, 21, 252–263. [Google Scholar] [CrossRef]
  10. Rumsey, C.L.; Carlson, J.R.; Pulliam, T.H.; Spalart, P.R. Improvements to the quadratic constitutive relation based on NASA juncture flow data. AIAA J. 2020, 58, 4374–4384. [Google Scholar] [CrossRef]
  11. Ling, J.; Kurzawski, A.; Templeton, J. Reynolds averaged turbulence modelling using deep neural networks with embedded invariance. J. Fluid Mech. 2016, 807, 155–166. [Google Scholar] [CrossRef]
  12. Xiao, H.; Wu, J.L.; Wang, J.X.; Sun, R.; Roy, C.J. Quantifying and reducing model-form uncertainties in Reynolds-averaged Navier–Stokes simulations: A data-driven, physics-informed Bayesian approach. J. Comput. Phys. 2016, 324, 115–136. [Google Scholar] [CrossRef]
  13. Singh, A.P.; Medida, S.; Duraisamy, K. Machine-learning-augmented predictive modeling of turbulent separated flows over airfoils. AIAA J. 2017, 55, 2215–2227. [Google Scholar] [CrossRef]
  14. Duraisamy, K.; Iaccarino, G.; Xiao, H. Turbulence modeling in the age of data. Annu. Rev. Fluid Mech. 2019, 51, 357–377. [Google Scholar] [CrossRef]
  15. Rodi, W. The Prediction of Free Turbulent Boundary Layers by Use of a Two-Equation Model of Turbulence. Ph.D. Thesis, University of London, London, UK, 1972. [Google Scholar]
  16. Launder, B.E.; Reece, G.J.; Rodi, W. Progress in the development of a Reynolds-stress turbulence closure. J. Fluid Mech. 1975, 68, 537–566. [Google Scholar] [CrossRef]
  17. Gatski, T.B.; Jongen, T. Nonlinear eddy viscosity and algebraic stress models for solving complex turbulent flows. Prog. Aerosp. Sci. 2000, 36, 655–682. [Google Scholar] [CrossRef]
  18. Pope, S.B. A more general effective-viscosity hypothesis. J. Fluid Mech. 1975, 72, 331–340. [Google Scholar] [CrossRef]
  19. Gatski, T.B.; Speziale, C.G. On explicit algebraic stress models for complex turbulent flows. J. Fluid Mech. 1993, 254, 59–78. [Google Scholar] [CrossRef]
  20. Wallin, S.; Johansson, A. An explicit algebraic Reynolds stress model for incompressible and compressible turbulent flows. J. Fluid Mech. 2000, 403, 89–132. [Google Scholar] [CrossRef]
  21. Wallin, S.; Johansson, A. Modeling streamline curvature effects in explicit algebraic Reynolds stress turbulence models. Int. J. Heat Fluid Flow 2002, 23, 721–730. [Google Scholar] [CrossRef]
  22. Ryu, S. Physics-to-Geometry Transformation to Construct Identities between Reynolds Stresses. Mathematics 2023, 11, 3698. [Google Scholar] [CrossRef]
  23. Ryu, S. A differential relation of metric properties for orientable smooth surfaces in ℝ3. Mathematics 2023, 11, 2337. [Google Scholar] [CrossRef]
  24. Krantz, S.G.; Parks, H.R. A Primer of Real Analytic Functions, 2nd ed.; Academic Press: Boston, MA, USA, 2002. [Google Scholar]
  25. Do Carmo, M.P. Differential Geometry of Curves and Surfaces, 2nd ed.; Dover: New York, NY, USA, 2016. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ryu, S. A Mathematically Exact and Well-Determined System of Equations to Close Reynolds-Averaged Navier–Stokes Equations. Mathematics 2023, 11, 4926. https://doi.org/10.3390/math11244926

AMA Style

Ryu S. A Mathematically Exact and Well-Determined System of Equations to Close Reynolds-Averaged Navier–Stokes Equations. Mathematics. 2023; 11(24):4926. https://doi.org/10.3390/math11244926

Chicago/Turabian Style

Ryu, Sungmin. 2023. "A Mathematically Exact and Well-Determined System of Equations to Close Reynolds-Averaged Navier–Stokes Equations" Mathematics 11, no. 24: 4926. https://doi.org/10.3390/math11244926

APA Style

Ryu, S. (2023). A Mathematically Exact and Well-Determined System of Equations to Close Reynolds-Averaged Navier–Stokes Equations. Mathematics, 11(24), 4926. https://doi.org/10.3390/math11244926

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop