Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulence modeling”

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 793 records · Page 44

Modeling of Turbulence Effects on Liquid Jet Atomization and Breakup

Recent experimental investigations and physical modeling studies have indicated that turbulence behaviors within a liquid jet have considerable effects on the atomization process. This study aims to model the turbulence effect in the atomization process of a cylindrical liquid jet. Two widely used models, the Kelvin-Helmholtz (KH) instability of Reitz (blob model) and the Taylor-Analogy-Breakup (TAB) secondary droplet breakup by O Rourke et al, are further extended to include turbulence effects. In the primary breakup model, the level of the turbulence effect on the liquid breakup depends on the characteristic scales and the initial flow conditions. For the secondary breakup, an additional turbulence force acted on parent drops is modeled and integrated into the TAB governing equation. The drop size formed from this breakup regime is estimated based on the energy balance before and after the breakup occurrence. This paper describes theoretical development of the current models, called "T-blob" and "T-TAB", for primary and secondary breakup respectivety. Several assessment studies are also presented in this paper.

Trinh, Huu P.↗

Numerical Modeling of Turbulent Combustion

The work in numerical modeling is focused on the use of the random vortex method to treat turbulent flow fields associated with combustion while flame fronts are considered as interfaces between reactants and products, propagating with the flow and at the same time advancing in the direction normal to themselves at a prescribed burning speed. The latter is associated with the generation of specific volume (the flame front acting, in effect, as the locus of volumetric sources) to account for the expansion of the flow field due to the exothermicity of the combustion process. The model was applied to the flow in a channel equipped with a rearward facing step. The results obtained revealed the mechanism of the formation of large scale turbulent structure in the wake of the step, while it showed the flame to stabilize on the outer edges of these eddies.

Ghoneim, A. F.↗

A review of Reynolds stress models for turbulent shear flows

A detailed review of recent developments in Reynolds stress modeling for incompressible turbulent shear flows is provided. The mathematical foundations of both two-equation models and full second-order closures are explored in depth. It is shown how these models can be systematically derived for two-dimensional mean turbulent flows that are close to equilibrium. A variety of examples are provided to demonstrate how well properly calibrated versions of these models perform for such flows. However, substantial problems remain for the description of more complex turbulent flows where there are large departures from equilibrium. Recent efforts to extend Reynolds stress models to nonequilibrium turbulent flows are discussed briefly along with the major modeling issues relevant to practical naval hydrodynamics applications.

Speziale, Charles G.↗

The effects of streamline curvature and swirl on turbulent flows in curved ducts

A technique for improving the numerical predictions of turbulent flows with the effect of streamline curvature is developed. Separated flows, the flow in a curved duct, and swirling flows are examples of flow fields where streamline curvature plays a dominant role. A comprehensive literature review on the effect of streamline curvature was conducted. New algebraic formulations for the eddy viscosity incorporating the kappa-epsilon turbulence model are proposed to account for various effects of streamline curvature. The loci of flow reversal of the separated flows over various backward-facing steps are employed to test the capability of the proposed turbulence model in capturing the effect of local curvature. The inclusion of the effect of longitudinal curvature in the proposed turbulence model is validated by predicting the distributions of the static pressure coefficients in an S-bend duct and in 180 degree turn-around ducts. The proposed turbulence model embedded with transverse curvature modification is substantiated by predicting the decay of the axial velocities in the confined swirling flows. The numerical predictions of different curvature effects by the proposed turbulence models are also reported.

Cheng, Chih-Hsiung↗

Pressure Loss Predictions of the Reactor Simulator Subsystem at NASA Glenn Research Center

Testing of the Fission Power System (FPS) Technology Demonstration Unit (TDU) is being conducted at NASA Glenn Research Center. The TDU consists of three subsystems: the reactor simulator (RxSim), the Stirling Power Conversion Unit (PCU), and the heat exchanger manifold (HXM). An annular linear induction pump (ALIP) is used to drive the working fluid. A preliminary version of the TDU system (which excludes the PCU for now) is referred to as the "RxSim subsystem" and was used to conduct flow tests in Vacuum Facility 6 (VF 6). In parallel, a computational model of the RxSim subsystem was created based on the computer-aided-design (CAD) model and was used to predict loop pressure losses over a range of mass flows. This was done to assess the ability of the pump to meet the design intent mass flow demand. Measured data indicates that the pump can produce 2.333 kg/sec of flow, which is enough to supply the RxSim subsystem with a nominal flow of 1.75 kg/sec. Computational predictions indicated that the pump could provide 2.157 kg/sec (using the Spalart-Allmaras (S‒A) turbulence model) and 2.223 kg/sec (using the k- turbulence model). The computational error of the predictions for the available mass flow is ‒0.176 kg/sec (with the S-A turbulence model) and -0.110 kg/sec (with the k- turbulence model) when compared to measured data.

Mass Flow↗

Numerical computation of shock wave-turbulent boundary layer interaction in transonic flow over an axisymmetric curved hill

A control-volume based finite difference computation of a turbulent transonic flow over an axisymmetric curved hill is presented. The numerical method is based on the SIMPLE algorithm, and hence the conservation of mass equation is replaced by a pressure correction equation for compressible flows. The turbulence is described by a k-epsilon turbulence model supplemented by a near-wall turbulence model. In the method, the dissipation rate in the region very close to the wall is obtained from an algebraic equation and that for the rest of the flow domain is obtained by solving a partial differential equation for the dissipation rate. The other flow equations are integrated up to the wall. It is shown that the present turbulence model yields the correct location of the compression shock. The other computational results are also in good agreement with experimental data.

Kim, S.-W.↗

Pressure Loss Predictions of the Reactor Simulator Subsystem at NASA GRC

Testing of the Fission Power System (FPS) Technology Demonstration Unit (TDU) is being conducted at NASA GRC. The TDU consists of three subsystems: the Reactor Simulator (RxSim), the Stirling Power Conversion Unit (PCU), and the Heat Exchanger Manifold (HXM). An Annular Linear Induction Pump (ALIP) is used to drive the working fluid. A preliminary version of the TDU system (which excludes the PCU for now), is referred to as the RxSim subsystem and was used to conduct flow tests in Vacuum Facility 6 (VF 6). In parallel, a computational model of the RxSim subsystem was created based on the CAD model and was used to predict loop pressure losses over a range of mass flows. This was done to assess the ability of the pump to meet the design intent mass flow demand. Measured data indicates that the pump can produce 2.333 kg/sec of flow, which is enough to supply the RxSim subsystem with a nominal flow of 1.75 kg/sec. Computational predictions indicated that the pump could provide 2.157 kg/sec (using the Spalart-Allmaras turbulence model), and 2.223 kg/sec (using the k- turbulence model). The computational error of the predictions for the available mass flow is -0.176 kg/sec (with the S-A turbulence model) and -0.110 kg/sec (with the k-epsilon turbulence model) when compared to measured data.

Turbulence Models↗

Numerical study of boundary layer interaction with shocks: Method and code validation

A major problem in modeling of turbulent supersonic flows is the correct assessment of viscous-inviscid interaction problems. Of particular interest is the interaction of boundary layers with shocks. Present turbulence models give in most cases unsatisfactory results in the region of rapid distortion and in the separation region (if one is present) in particular with regard to mean flow profiles and turbulence quantities. The objective of the present work is the direct numerical simulation of shock boundary layer interaction. This report summarizes the first phase during which a numerical method suitable for this problem has been developed and a computer code has been written and tested.

Adams, Nikolaus A.↗

Supersonic Combusting Jet Experiments for Code Development and Validation

Computational fluid dynamics (CFD) methods based on the Reynolds averaged Navier-Stokes (RANS) equations are extensively employed in the design of hypersonic airbreathing engines. A fundamental weakness in these methods is the accurate mathematical modeling of turbulence and turbulence-combustion interactions. This paper is one of a companion pair of papers which describe the use of the combined dual pump CARS and newly developed interferometric Rayleigh scattering optical system to acquire time-accurate turbulence and mixing data with which to aid calibration and verification of the current turbulence models. This test program consisted of a pilot laboratory-scale experiment used to develop the optical system and data acquisition techniques, and a significantly larger experiment to acquire the necessary turbulence data. The flow is an axially-symmetric, supersonic, combusting, free jet that provides good optical access, consisting of a central jet of hot "vitiated air" and a coflow jet of hydrogen or ethylene fuel. In this paper the development of both experiments is described. Facility and flow visualization data are presented for various types of flames, including the flames selected for detailed study with the CARS/Rayleigh optical techniques.

Cutler, Andrew D.↗

Computational Fluid Dynamics Combustion Modeling for Rotating Detonation Engines

This paper focuses on the development and validation of a combustion model for Computational Fluid Dynamics (CFD) modeling of Rotating Detonation Engines. A zero-dimensional Partially Stirred Reactor (PaSR) with a detailed chemical kinetic mechanism for hydrogen and air is used to model turbulent combustion. The model is computationally efficient and is based on the notion of partial mixing at the sub-grid level with turbulent exchange between mixed and unmixed regions. The ability of the PaSR model to accurately represent both detonative and deflagrative combustion is assessed by validating the results against experimental data. The effects of mesh resolution on the solution are also studied in order to determine if a mesh independent solution is obtainable with the Large Eddy Simulation (LES) approach to modeling turbulence. A comparison is made between the PaSR model and simply ignoring turbulence chemistry interactions which assumes that all species are perfectly mixed at the sub-grid level.

Strakey, Peter↗

Navier-Stokes analysis of turbomachinery blade external heat transfer

The two-dimensional, compressible, thin-layer Navier-Stokes and energy equations were solved numerically to obtain heat transfer rates on turbomachinery blades. The Baldwin-Lomax algebraic model and the q-omega low Reynolds number two-equation model were used for modeling of turbulence. For the numerical solution of the governing equations a four-stage Runge-Kutta solver was employed. The turbulence model equations were solved using an implicit scheme. Numerical solutions are presented for two-dimensional flow within two vane cascades. The heat transfer results and the pressure distributions were compared with published experimental data. The agreement between the numerical calculations and the experimental values were found to be generally favorable. The position of transition from laminar to turbulent flow was also predicted accurately.

Gorla, Rama S. R.↗

Navier-Stokes analysis of turbomachinery blade external heat transfer

The two-dimensional, compressible, thin-layer Navier-Stokes and energy equations were solved numerically to obtain heat transfer rates on turbomachinery blades. The Baldwin-Lomax algebraic model and the q - omega low Reynolds number, two-equation model were used for modeling of turbulence. For the numerical solution of the governing equations a four-stage Runge-Kutta solver was employed. The turbulence model equations were solved using an implicit scheme. Numerical solutions are presented for two-dimensional flow within two vane cascades. The heat transfer results and the pressure distributions were compared with published experimental data. The agreement between the numerical calculations and the experimental values were found to be generally favorable. The position of transition from laminar to turbulent flow was also predicted accurately.

Ameri, A. A.↗

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗

Military specifications

The current situation relative to the military specification is that there is not one specific model of turbulence which people are using. Particular disagreement exists on how turbulence levels will vary with qualitative analysis. It does not tie one down to specifics. When it comes to flying quality specifications, many feel that one should stay with the definitions of the Cooper-Harper rating scale but allow the levels to shift depending on the level of turbulence. There is a ride quality specification in the MIL-SPEC having to do with flight control systems design that is related to a turbulence model. This spec (MIL-F8785C) and others are discussed.

Reynolds, Philip↗

An eddy viscosity for three-dimensional boundary-layer flows

This paper proposes an isotropic eddy viscosity model for three-dimensional boundary-layer flows which is an improved version of the Cebeci-Smith model. Both turbulence models are tested by comparison with two experiments. The proposed model is found to perform much better than the Cebeci-Smith model. This improvement is due to a reduction of the eddy viscosity coefficient.

Abid, Ridha↗

A cell-vertex multigrid method for the Navier-Stokes equations

A cell-vertex scheme for the Navier-Stokes equations, which is based on central difference approximations and Runge-Kutta time stepping, is described. Using local time stepping, implicit residual smoothing, a multigrid method, and carefully controlled artificial dissipative terms, very good convergence rates are obtained for a wide range of two- and three-dimensional flows over airfoils and wings. The accuracy of the code is examined by grid refinement studies and comparison with experimental data. For an accurate prediction of turbulent flows with strong separations, a modified version of the nonequilibrium turbulence model of Johnson and King is introduced, which is well suited for an implementation into three-dimensional Navier-Stokes codes. It is shown that the solutions for three-dimensional flows with strong separations can be dramatically improved, when a nonequilibrium model of turbulence is used.

Radespiel, R.↗

Numerical aspects of computing high Reynolds number flows on unstructured meshes

An edge-data structure describing a mesh edge-wise given the vertices of each edge and neighboring cell information is used developing algorithms for the Navier-Stokes equations on triangular meshes. Edge formulas for the Galerkin and finite-element discretization of gradient, divergence, Hessian, and Laplacian operators are derived. A simple edge formula is derived for the discretization of the Laplacian operator, where precise theoretical conditions for a discrete maximum principle can be ascertained. Practical issues associated with solving the Navier-Stokes equations on unstructured meshes are addressed, along with issues concerning the generation of highly stretched triangular meshes and the modeling of turbulence on unstructured meshes. A turbulence modeling strategy is proposed, and numerical results for a high-Reynolds-number flow about single- and multielement airfoils are discussed.

Barth, Timothy J.↗