Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulent flows”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

Tetrahedral-Mesh Simulation of Turbulent Flows with the Space-Time Conservative Schemes

Direct numerical simulations of turbulent flows are predominantly carried out using structured, hexahedral meshes despite decades of development in unstructured mesh methods. Tetrahedral meshes offer ease of mesh generation around complex geometries and the potential of an orientation free grid that would provide un-biased small-scale dissipation and more accurate intermediate scale solutions. However, due to the lack of consistent multi-dimensional numerical formulations in conventional schemes for triangular and tetrahedral meshes at the cell interfaces, numerical issues exist when flow discontinuities or stagnation regions are present. The space-time conservative conservation element solution element (CESE) method - due to its Riemann-solver-free shock capturing capabilities, non-dissipative baseline schemes, and flux conservation in time as well as space - has the potential to more accurately simulate turbulent flows using unstructured tetrahedral meshes. To pave the way towards accurate simulation of shock/turbulent boundary-layer interaction, a series of wave and shock interaction benchmark problems that increase in complexity, are computed in this paper with triangular/tetrahedral meshes. Preliminary computations for the normal shock/turbulence interactions are carried out with a relatively coarse mesh, by direct numerical simulations standards, in order to assess other effects such as boundary conditions and the necessity of a buffer domain. The results indicate that qualitative agreement with previous studies can be obtained for flows where, strong shocks co-exist along with unsteady waves that display a broad range of scales, with a relatively compact computational domain and less stringent requirements for grid clustering near the shock. With the space-time conservation properties, stable solutions without any spurious wave reflections can be obtained without a need for buffer domains near the outflow/farfield boundaries. Computational results for the isotropic turbulent flow decay, at a relatively high turbulent Mach number, show a nicely behaved spectral decay rate for medium to high wave numbers. The high-order CESE schemes offer very robust solutions even with the presence of strong shocks or widespread shocklets. The explicit formulation in conjunction with a close to unity theoretical upper Courant number bound has the potential to offer an efficient numerical framework for general compressible turbulent flow simulations with unstructured meshes.

Chang, Chau-Lyan↗

High Order Numerical Methods for LES of Turbulent Flows with Shocks

Simulation of turbulent flows with shocks employing explicit subgrid-scale (SGS) filtering may encounter a loss of accuracy in the vicinity of a shock. In this work we perform a comparative study of different approaches to reduce this loss of accuracy within the framework of the dynamic Germano SGS model. One of the possible approaches is to apply Harten's subcell resolution procedure to locate and sharpen the shock, and to use a one-sided test filter at the grid points adjacent to the exact shock location. The other considered approach is local disabling of the SGS terms in the vicinity of the shock location. In this study we use a canonical shock-turbulence interaction problem for comparison of the considered modifications of the SGS filtering procedure. For the considered test case both approaches show a similar improvement in the accuracy near the shock.

High Order↗

High Order Numerical Methods for LES of Turbulent Flows with Shocks

Simulation of turbulent flows with shocks employing explicit subgrid-scale (SGS) filtering may encounter a loss of accuracy in the vicinity of a shock. In this work we perform a comparative study of different approaches to reduce this loss of accuracy within the framework of the dynamic Germano SGS model. One of the possible approaches is to apply Harten’s subcell resolution procedure to locate and sharpen the shock, and to use a one-sided test filter at the grid points adjacent to the exact shock location. The other considered approach is local disabling of the SGS terms in the vicinity of the shock location. In this study we use a canonical shock-turbulence interaction problem for comparison of the considered modifications of the SGS filtering procedure. For the considered test case both approaches show a similar improvement in the accuracy near the shock.

LES↗

A three-dimensional incompressible flow simulation method and its application to the Space Shuttle main engine. II Turbulent flow

An implicit finite difference code cast in general curvilinear coordinates is further developed for three-dimensional incompressible turbulent flows. The code is based on the method of pseudocompressibility and utilizes the Beam and Warming implicit approximate factorization algorithm to achieve computational efficiency. A multiple-zone method is further extended to include composite-grids to overcome the excessive computer memory required for solving turbulent flows in complex three-dimensional geometries. A simple turbulence model is proposed for internal flows. The code is being used for the Space Shuttle Main Engine (SSME) internal flow analyses.

Chang, J. L. C.↗

Second-order closure models for supersonic turbulent flows

Recent work on the development of a second-order closure model for high-speed compressible flows is reviewed. This turbulent closure is based on the solution of modeled transport equations for the Favre-averaged Reynolds stress tensor and the solenoidal part of the turbulent dissipation rate. A new model for the compressible dissipation is used along with traditional gradient transport models for the Reynolds heat flux and mass flux terms. Consistent with simple asymptotic analyses, the deviatoric part of the remaining higher-order correlations in the Reynolds stress transport equations are modeled by a variable density extension of the newest incompressible models. The resulting second-order closure model is tested in a variety of compressible turbulent flows which include the decay of isotropic turbulence, homogeneous shear flow, the supersonic mixing layer, and the supersonic flat-plate turbulent boundary layer. Comparisons between the model predictions and the results of physical and numerical experiments are quite encouraging.

Speziale, Charles G.↗

Second-order closure models for supersonic turbulent flows

Recent work by the authors on the development of a second-order closure model for high-speed compressible flows is reviewed. This turbulence closure is based on the solution of modeled transport equations for the Favre-averaged Reynolds stress tensor and the solenoidal part of the turbulent dissipation rate. A new model for the compressible dissipation is used along with traditional gradient transport models for the Reynolds heat flux and mass flux terms. Consistent with simple asymptotic analyses, the deviatoric part of the remaining higher-order correlations in the Reynolds stress transport equation are modeled by a variable density extension of the newest incompressible models. The resulting second-order closure model is tested in a variety of compressible turbulent flows which include the decay of isotropic turbulence, homogeneous shear flow, the supersonic mixing layer, and the supersonic flat-plate turbulent boundary layer. Comparisons between the model predictions and the results of physical and numerical experiments are quite encouraging.

Speziale, Charles G.↗

An approximate inertial manifold (AIM) based closure for turbulent flows

A closure model for turbulent flows is developed based on a dynamical system theory. An appropriately discretized formulation of the governing equations is considered for this process. The key ingredient is an approximation of the system’s attractor, where all the trajectories in phase space are confined. This approximate inertial manifold based approach provides a path to track trajectories of the system in a lower-dimensional subspace. Unlike conventional coarse-graining approaches, the turbulent field is decomposed into resolved and unresolved dynamics using the properties of the governing equations. The novelty of the approach relies on the reconstruction of the unresolved field constrained by the governing equations. A posteriori tests for homogeneous isotropic turbulence and the Kuramoto–Sivashinsky equation show promising results for considerable dimension reduction with strong convergence properties. The proposed model outperforms the dynamic Smagorinsky model, and the computational overhead is competitive with similar approaches.

97 MATHEMATICS AND COMPUTING↗

Studies in Transition and Time Varying Turbulent Flows

The research focused on two areas: (a) the dynamics of forced turbulent flows and (b) time filtered Large Eddy Simulations (TLES). The dynamics of turbulent flows arising from external forcing of the turbulence are poorly understood. In particular, here are many unanswered questions relating the basic dynamical balances and the existence or nonexistence of statistical equilibrium of forced turbulent flows. The research used direct numerical simulations to explore these questions. The properties of the temporally filtered Navier-Stokes equations were also studied.

Grosch, Chester E.↗

Studies in Forced and Time Varying Turbulent Flows

The reasearch focused on two areas; (a) the dynamics of forced turbulent flows and (b) time filtered Large Eddy Simulations (TLES). The dynamics of turbulent flows arising from external forcing of the turbulence are poorly understood. In particular, here are many unanswered questions relating the basic dynamical balances and the existence or nonexistence of statistical equilibrium of forced turbulent flows. This research used rapid distortion theory and direct numerical simulations to explore these questions. The properties of the temporally filtered Navier-Stokes equations were also studied.

Grosch, Chester E.↗

Numerical investigation of an internal layer in turbulent flow over a curved hill

A numerical investigation of incompressible and compressible turbulent flows over strongly curved surfaces is presented. The turbulent flow equations are solved by a pressure based Navier-Stokes equations solver. In the method, the conservation of mass equation is replaced by a pressure correction equation applicable for both compressible and incompressible flows. The turbulence is described by a multiple time scale turbulence model supplemented with a near-wall turbulence model. The numerical results show that the internal layer is a strong turbulence field which is developed beneath the external boundary layer and is located very close to the wall. The development of the internal layer is attributed to the enormous mean flow strain rate caused by the streamline curvature. The external boundary layer flow responds rather slowly to the streamline curvature. Thus, the turbulence field of the forward corner of the curved hill is characterized by two turbulence fields interacting with each other. The turbulence intensity of the internal layer is much stronger than that of the external boundary layer, so that the development of a new boundary layer in the downstream region of the curved hill depends mostly on the internal layer. These numerical results are in good agreement with the measured data, and show that the turbulence model can resolve the turbulence field subjected to the strong streamline curvature.

Kim, S.-W.↗

Anisotropic Turbulent Flows Observed in Above-the-loop-top Regions during Solar Flares

Abstract Solar flare above-the-loop-top (ALT) regions are vital for understanding solar eruptions and fundamental processes in plasma physics. Recent advances in three-dimensional (3D) magnetohydrodynamic (MHD) simulations have revealed unprecedented details on turbulent flows and MHD instabilities in flare ALT regions. Here, for the first time, we examine the observable anisotropic properties of turbulent flows in ALT by applying a flow-tracking algorithm on narrow-band extreme-ultraviolet images that are observed from the face-on viewing perspective. First, the results quantitatively confirm the previous observation that vertical motions dominate and that the anisotropic flows are widely distributed in the entire ALT region with the contribution from both upflows and downflows. Second, the anisotropy shows height-dependent features, with the most substantial anisotropy appearing at a certain middle height in ALT, which agrees well with the MHD modeling results where turbulent flows are caused by Rayleigh–Taylor-type instabilities in the ALT region. Finally, our finding suggests that supra-arcade downflows (SADs), the most prominently visible dynamical structures in ALT regions, are only one aspect of turbulent flows. Among these turbulent flows, we also report the antisunward-moving underdense flows that might develop due to MHD instabilities, as suggested by previous 3D flare models. Our results indicate that the entire flare fan displays group behavior of turbulent flows where the observational bright spikes and relatively dark SADs exhibit similar anisotropic characteristics.

Xie, Xiaoyan (ORCID:0009000705827807)↗

Self similarity of two point correlations in wall bounded turbulent flows

The structure of turbulence at a height y from a wall is affected by the local mean shear at y, by the direct effect of the wall on the eddies, and by the action of other eddies close to or far from the wall. Some researchers believe that a single one of these mechanisms is dominant, while others believe that these effects have to be considered together. It is important to understand the relative importance of these effects in order to develop closure models, for example for the dissipation or for the Reynolds stress equation, and to understand the eddy structure of cross correlation functions and other measures. The specific objective was to examine the two point correlation, R sub vv, of the normal velocity component v near the wall in a turbulent channel flow and in a turbulent boundary layer. The preliminary results show that even in the inhomogeneous turbulent boundary layer, the two-point correlation function may have self similar forms. The results also show that the effects of shear and of blocking are equally important in the form of correlation functions for spacing normal to the wall. But for spanwise spacing, it was found that the eddy structure is quire different in these near flows. So any theory for turbulent structure must take both these effects into account.

Hunt, J. C. R.↗

Numerical study of laminar and turbulent flows inside a turnaround duct with and without guide vanes

The purpose of this study is to examine in detail incompressible laminar and turbulent flows inside a turnaround duct with and without guide vanes and to investigate the effects of vanes on the flow characteristics. To perform this study, an implicit finite difference code cast in general curvilinear coordinates is further developed. The code is based on the method of pseudo-compressibility and utilize ADI or implicit approximate factorization algorithm to achieve computational efficiency. Method of segmental sweeping is developed to overcome the multiple-zone problem due to imposition of guide vanes. In the present study, several test cases have been computed. These test cases include laminar and turbulent flows inside a turnaround duct without and with two or three guide vanes. The study reveals that: (1) there exists large recirculation zones inside the duct if no vanes are present; (2) properly shaped and positional guide vanes are effective in eliminating flow separation; and (3) laminar and turbulent flows have similar flow features. But turbulent flow has less total pressure drop.

Lin, S.-J.↗

The Turbulent Flow in Diffusers of Small Divergence Angle

The turbulent flow in a conical diffuser represents the type of turbulent boundary layer with positive longitudinal pressure gradient. In contrast to the boundary layer problem, however, it is not necessary that the pressure distribution along the limits of the boundary layer(along the axis of the diffuser) be given, since this distribution can be obtained from the computation. This circumstance, together with the greater simplicity of the problem as a whole, provides a useful basis for the study of the extension of the results of semiempirical theories to the case of motion with a positive pressure gradient. In the first part of the paper,formulas are derived for the computation of the velocity and.pressure distributions in the turbulent flow along, and at right angles to, the axis of a diffuser of small cone angle. The problem is solved.

Gourzhienko, G. A.↗

CFD modeling of turbulent flows around the SSME main injector assembly using porosity formulation

Hot gas turbulent flow distribution around the main injector assembly of the Space Shuttle Main Engine (SSME) and LOX flow distribution through the LOX posts have a great effect on the combustion phenomena inside the main combustion chamber. In order to design a CFD model to be an effective engineering analysis tool with good computational turn-around time (especially for 3-D flow problems) and still maintain good accuracy in describing the flow features, the concept of porosity was employed to describe the effects of blockage and drag force due to the presence of the LOX posts in the turbulent flow field around the main injector assembly of the SSME. Two-dimensional numerical studies were conducted to identify the drag coefficients of the flows, both through tube banks and round the shielded posts, over a wide range of Reynolds numbers. Empirical, analytical expressions of the drag coefficients as a function of local flow Reynolds number were then deduced. The porosity model was applied to the turbulent flow around the main injector assembly of the SSME, and analyses were performed. The 3-D CFD analysis was divided into three parts: LOX dome, hot gas injector assembly, and hydrogen cavity. The numerical results indicate that the mixture ratio at the downstream of injector face was close to stoichiometric around baffle elements.

Cheng, Gary C.↗