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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Predictive turbulence-driven flux model of scrape-off layer widths across confinement regimes in tokamaks

Reliable scrape-off layer (SOL) profile decay lengths predictions are needed to design and operate future tokamaks. The present manuscript describes a new model based on turbulent transport that is able to predict SOL widths for both L-mode and H-mode plasmas. The model is based upon the sheared-spectral filament paradigm (Peret et al (WEST Team) 2022 Phys. Plasmas 29 072306), however, incorporating the effects of thermal transport in order to calculate the parallel heat fluxes. The effects of magnetic shear and ExB shear on the cross-field transport are crucial to explain the shorter SOL decay lengths found in H-mode. The model is validated against a database of thousands of DIII-D L-mode and H-mode SOL profiles. We also calculate SOL decay length predictions in terms of plasma and engineer control parameters, which are in agreement with the multi-machine empirical H-mode scaling (Eich et al (ASDEX Upgrade Team and JET EFDA Contributors) 2013 Nucl. Fusion 53 093031), however, with an additional device geometry dependence. ITER SOL width predictions by the model are 3 times higher than the empirical scaling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Design and Modeling of an 80 bar Oxy-Combustor for Direct Fired Supercritical CO2 Applications

As supercritical CO2 power cycles for fossil energy power generation continue to generate interest, there are significant issues and unanswered questions regarding injector design, flame stabilization, wall heat transfer, CO emissions, combustion dynamics and other combustion phenomenon. For natural gas, direct-fired cycles with carbon capture it is believed that Computational Fluid Dynamics (CFD) modeling will play an essential role in the combustor design process. To accurately model turbulent reacting flows at these unique conditions, experimental data is needed to validate CFD codes and sub-models and is currently lacking at conditions relevant for these cycles. This paper presents the conceptual design and CFD simulations of an experimental 80 bar oxy-combustion facility and test article currently under construction at NETL. The facility is targeted towards the testing of a single injector, direct-fired sCO2 combustor at the 100 kW thermal output level. While these conditions do not reflect the actual Allam cycle operating conditions (300 bar) they are viewed as a stepping stone in the model validation process at supercritical conditions. Reynolds Averaged Navier-Stokes as well as Large Eddy Simulation are used to model the turbulent combustion process and aid in the design of the combustor and injector. Process parameters including oxidizer preheat temperature and combustor flowrate are investigated.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1995

Computational analysis of forebody tangential slot blowing on the high alpha research vehicle

A numerical analysis of forebody tangential slot blowing as a means of generating side force and yawing moment is conducted using an aircraft geometry. The Reynolds-averaged, thin-layer, Navier-Stokes equations are solved using a partially flux-split, approximately-factored algorithm. An algebraic turbulence model is used to determine the turbulent eddy viscosity values. Solutions are obtained using both patched and overset grid systems. In the patched grid model, and actuator plane is used to introduce jet variables into the flow field. The overset grid model is used to model the physical slot geometry and facilitate modeling of the full aircraft configuration. A slot optimization study indicates that a short slot located close to the nose of the aircraft provided the most side force and yawing moment per unit blowing coefficient. Comparison of computed surface pressure with that obtained in full-scale wind tunnel tests produce good agreement, indicating the numerical method and grid system used in the study are valid. Full aircraft computations resolve the changes in vortex burst point due to blowing. A time-accurate full-aircraft solution shows the effect of blowing on the changes in the frequency of the aerodynamic loads over the vertical tails. A study of the effects of freestream Mach number and various jet parameters indicates blowing remains effective through the transonic Mach range. An investigation of the force onset time lag associated with forebody blowing shows the lag to be minimal. The knowledge obtained in this study may be applied to the design of a forebody tangential slot blowing system for use on flight aircraft.

Gee, Ken↗
NASA NTRS1990

A spectrum model for weakly anisotropic turbulence

A simple model, parameterized by the Reynolds stress anisotropy, is proposed for the spectrum of weakly anisotropic turbulence. It contains a model constant that affects its region of realizability. This spectrum model is used to derive a one-point closure to the rapid pressure-strain term. The derived pressure-strain model is linear in the Reynolds stress anisotropy and is of the same form as the closure model of Launder et al. (1975). The spectrum model becomes unrealizable in some regions of wave space for sufficiently large anisotropy of the Reynolds stress, and this is used to infer the region of validity of the linear closure model. It is found that the extent of the valid region is very small when the model constant is set to match rapid distortion theory, and largest for a model-constant set close to the value suggested by LRR. However, even the largest valid domain does not extend very far from isotropy, suggesting inherent weakness in the linear pressure-strain models.

Shih, Tsan-Hsing↗
NASA NTRS1989

Distortion of homogeneous turbulence by axisymmetric strain and dilatation

Rapid distortion theory is applied to study distortion of homogeneous turbulence subject to two different axisymmetric strain modes: the axisymmetric contraction (AC, nozzle-type flow), and the axisymmetric expansion (AE, diffuser-type flow). The paper explores the differences in effects of the two axisymmetric strain modes on the anisotropy of correlations and structures of turbulence; examines the effect of dilatation on the distortion of turbulence; and provides a theoretical background for turbulence model development. It is found that velocity and vorticity fluctuations are enhanced more efficiently by contraction than by expansion; contraction produces much higher anisotropy in velocity and vorticity than expansion; root-mean-square pressure is slightly reduced during contraction, whereas it increases rapidly during expansion; and vortical structures of rodlike shape develop in a contraction flow, while disklike structures develop in an expansion flow. A simple model that reflects the dependence of turbulence evolution on structural parameters such as the Reynolds-stress anisotropy and total strain is proposed, and is shown to outperform all other models for all cases examined, regardless of the mean strain rate.

Lee, Moon Joo↗
NASA NTRS1990

On the Subgrid-Scale Modeling of Compressible Turbulence

A new sub-grid scale model is presented for the large-eddy simulation of compressible turbulence. In the proposed model, compressibility contributions have been incorporated in the sub-grid scale eddy viscosity which, in the incompressible limit, reduce to a form originally proposed by Smagorinsky (1963). The model has been tested against a simple extension of the traditional Smagorinsky eddy viscosity model using simulations of decaying, compressible homogeneous turbulence. Simulation results show that the proposed model provides greater dissipation of the compressive modes of the resolved-scale velocity field than does the Smagorinsky eddy viscosity model. For an initial r.m.s. turbulence Mach number of 1.0, simulations performed using the Smagorinsky model become physically unrealizable (i.e., negative energies) because of the inability of the model to sufficiently dissipate fluctuations due to resolved scale velocity dilations. The proposed model is able to provide the necessary dissipation of this energy and maintain the realizability of the flow. Following Zeman (1990), turbulent shocklets are considered to dissipate energy independent of the Kolmogorov energy cascade. A possible parameterization of dissipation by turbulent shocklets for Large-Eddy Simulation is also presented.

Squires, Kyle↗
DOE OSTI2020

Three-dimensional signatures of self-similarity in a high-energy-density plasma shear-driven mixing layer

A hydrodynamic shear mixing layer experiment at the National Ignition Facility had previously demonstrated Eulerian scaling of integrated, late-time quantities, including turbulent kinetic energy. In this manuscript, the experiment is repeated with new materials. In this work, using the new dataset, we demonstrate that Euler-number scalings hold not just for late time, but dynamically throughout the experiment, for measurements in all three spatial dimensions. Incorporating the dynamic scaling leads to an enhanced calculation that the heavier of the two scaled experiments has approached three generations of mergers of its primary instability's structures and a consistent observation of such a merger in action in the lighter of the two scaled experiments. Furthermore, the improved scrutiny of the time evolution of instability structures leads to sharper estimates of turbulent kinetic energy, including a demonstration of different behaviors correlating with surface roughness (quantitatively consistent with transitions between laminar and turbulent initial states), as predicted by a Reynolds-averaged turbulent model, which evidently correctly handles the differing shock-roughness interactions to drive its internal state of the model into different regimes. Altogether, a picture arises of the analytical improvements in treating these variations (of times, densities, and roughnesses) as a unified whole and of multiple ways by which deviations from the scaling could indicate an onset of non-hydrodynamic behavior. Such deviations were not expected for these experiments (which models correctly indicated would remain hydrodynamic) but could be introduced by, for example, imposing external fields or increasing drive energy to test conditions relevant to inertial confinement fusion or other high-energy-density experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

Simulation of Rotor Blade Element Turbulence

A turbulence model has been developed for blade-element helicopter simulation. This model, called Simulation of Rotor Blade Element Turbulence (SORBET), uses an innovative temporal and geometrical distribution algorithm that preserves the statistical characteristics of the turbulence spectra over the rotor disc, while providing velocity components in real time to each of five blade-element stations along each of four blades. An initial investigation of SORBET has been performed using a piloted, motion-based simulation of the Sikorsky UH60A Black Hawk. Although only the vertical component of stochastic turbulence was used in this investigation, vertical turbulence components induce vehicle responses in all translational and rotational degrees of freedom of the helicopter. The single-degree-of-freedom configuration of SORBET was compared to a conventional full 6-degrees-of-freedom baseline configuration, where translational velocity inputs are superimposed at the vehicle center of gravity, and rotational velocity inputs are created from filters that approximate the immersion rate into the turbulent field. For high-speed flight the vehicle responses were satisfactory for both models. Test pilots could not distinguish differences between the baseline configuration and SORBET. In low-speed flight the baseline configuration received criticism for its high frequency content, whereas the SORBET model elicited favorable pilot opinion. For this helicopter, which has fully articulated blades, results from SORBET show that vehicle responses to turbulent blade-station disturbances are severely attenuated. This is corroborated by in-flight observation of the rotor tip path plane as compared to vehicle responses.

McFarland, R. E.↗
NASA NTRSDate not supplied

CFD Model Development of a Cryogenic Storage Tank Self-Pressurization in Normal Gravity and Validation against SHIIVER Experiment

Two-Phase flow and heat transfer simulations with interfacial phase change of the Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER) self-pressurization experiment were conducted using storage tank CFD model in the framework of the ANSYS Fluent CFD code. The simulations were performed for the 70% fill level case with MLI on domes and no vapor cooling. All the phase change calculations in these simulations were generated by in-house Schrage-based evaporation-condensation model. The calculations were performed using both the VOF and Sharp Interface multiphase 2Daxisymmetric models. A number of parametric and sensitivity studies were performed to check the various aspects of the CFD model. These studies helped to understand the effects of varying several parameters on the tank pressure and temperature during self-pressurization. These parameters include: the effect of turbulence modeling and turbulence damping at the interface; the effect of constant vs. temperature dependent fluid properties; the effect of accommodation coefficient; the effect of modeling phase change at the interface and the effect of the fill level.

Computational Fluid Dynamics↗
DOE OSTI2022

Toward the core-edge coupling of delta-f and total-f gyrokinetic models

The coupling of core reduced delta-f and edge total-f gyrokinetic models could enable a significant speed-up of the whole device modeling simulations. It is thus interesting to compare these models in the core where they will be coupled. The main difference between these models is the presence of the zeroth order term on the right-hand side of the total-f gyrokinetic equation. This term is responsible for driving a radial electric field that is also associated with Geodesic acoustic mode-like oscillations. To investigate the coupling between these models, the subsequent large-scale gyrokinetic 3D turbulence simulations have been performed starting from an already saturated gyrokinetic axisymmetric equilibria (collisions are not included). This work has also been extended to couple different axisymmetric and turbulent models available in XGC. It is shown that the radial electric field and its drive have to be taken care of consistently while coupling different models together.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

Structure of turbulence in three-dimensional boundary layers

This report provides an overview of the three dimensional turbulent boundary layer concepts and of the currently available experimental information for their turbulence modeling. It is found that more reliable turbulence data, especially of the Reynolds stress transport terms, is needed to improve the existing modeling capabilities. An experiment is proposed to study the three dimensional boundary layer formed by a 'sink flow' in a fully developed two dimensional turbulent boundary layer. Also, the mean and turbulence field measurement procedure using a three component laser Doppler velocimeter is described.

Subramanian, Chelakara S.↗
NASA NTRS1989

Turbulent eddy viscosity modeling in transonic shock/boundary-layer interactions

The treatment of turbulence effects on transonic shock/turbulent boundary layer interaction is addressed within the context of a triple deck approach valid for arbitrary practical Reynolds numbers between 1000 and 10 billion. The modeling of the eddy viscosity and basic turbulent boundary profile effects in each deck is examined in detail using Law-of-the-Wall/Law-of-the-Wake concepts as the foundation. Results of parametric studies showing how each of these turbulence model aspects influences typical interaction zone property distributions (wall pressure, displacement thickness and local skin friction) are presented and discussed.

Inger, G. R.↗
NASA NTRS1981

Mixing of a passive scalar in isotropic and sheared homogeneous turbulence

In order to calculate the velocity and scalar fields, the three dimensional, time-dependent equations of motion and the diffusion equation were solved numerically. The following cases were treated: isotropic, homogeneous turbulence with decay of a passive scalar; and homogeneous turbulent shear flow with a passive scalar whose mean varies linearly in the spanwise direction. The solutions were obtained at relatively low Reynolds numbers so that all of the turbulent scales could be resolved without modeling. Turbulent statistics such as integral length scales, Taylor microscales, Kolmogorov length scale, one- and two-point correlations of velocity-velocity and velocity-scalar, turbulent Prandtl/Schmidt number, r.m.s. values of velocities, the scalar quantity and pressure, skewness, decay rates, and decay exponents were calculated. The results are compared with the available expermental results, and good agreement is obtained.

Shirani, E.↗
NASA NTRS1986

Computational simulation of liquid rocket injector anomalies

A computer model has been developed to analyze the three-dimensional two-phase reactive flows in liquid fueled rocket combustors. The model is designed to study the influence of liquid propellant injection nonuniformities on the flow pattern, combustion and heat transfer within the combustor. The Eulerian-Lagrangian approach for simulating polidisperse spray flow, evaporation and combustion has been used. Full coupling between the phases is accounted for. A nonorthogonal, body fitted coordinate system along with a conservative control volume formulation is employed. The physical models built into the model include a kappa-epsilon turbulence model, a two-step chemical reaction, and the six-flux radiation model. Semiempirical models are used to describe all interphase coupling terms as well as chemical reaction rates. The purpose of this study was to demonstrate an analytical capability to predict the effects of reactant injection nonuniformities (injection anomalies) on combustion and heat transfer within the rocket combustion chamber. The results show promising application of the model to comprehensive modeling of liquid propellant rocket engines.

Przekwas, A. J.↗
NASA NTRS1988

Numerical study of turbulent internal shear layer flow in an axisymmetric U-duct

This paper presents a simple algebraic turbulence model developed for internal flows which is based on the vorticity length and magnitude. The length scale is defined using the distance between the maximum and the minimum points of the absolute vorticity. This model is numerically tested in a turbulent internal layer flow through an axisymmetric U-duct with very sharp curvature; the ratio of the boundary layer thickness to the radius of curvature of the duct is of order 1. In this U-duct flow, strong adverse and favorable pressure gradients coexist and interact with each other. Satisfactoy agreement with experimental results is obtained.

Chang, James L. C.↗
NASA NTRS1990

Fully elliptic incompressible flow calculations on regular grid by a new pressure substitution method

A new method is presented for the solution of incompressible flow in generalized coordinates. This method is based on the substitution of the pressure weighted form of the momentum equations into the continuity equation. The algorithm is rigorously derived and a Fourier analysis is used to assess its suitability to act as an error smoother. Linear stability analysis results indicate that the performance of the new pressure substitution method (PSM) and the pressure correction method (PCM) is about the same at low Reynolds numbers, with no significant pressure gradient. At high Reynolds numbers the PSM shows much faster convergence. Likewise prediction of various flows indicate that the PSM has better accuracy for high Reynolds number flows with significant pressure gradients. Since most practical aerodynamic flows have significant pressure gradients, the PSM seems to be attractive for such flows. Solutions for both laminar and turbulent flow are compared with the experimental data. A two-equation low Reynolds number turbulence model is used to resolve the turbulent flowfield.

Hobson, G. V.↗
NASA NTRS1991

Numerical analysis of ignition transients in solid rocket motors

A model to analyze the unstready multidimensional turbulent flow in the head-end, star slot section of a solid rocket motor during the ignition transient is developed. The present paper examines the fluid dynamic aspects of the starting transient. A finite-difference solution of the unsteady, compressible, full Navier-Stokes equations, together with a two-equation k-epsilon turbulence model, is obtained utilizing the MacCormack explicit, predictor-corrector technique. Computed results for the flow in the slot are compared with experimental data obtained in cold flow tests of a scaled model of the Space Shuttle SRM head-end, star grain section. The agreement between the calculated flowfield and existing experimental data is very good.

Ciucci, A.↗