Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulent Simulation”

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 235 records · Page 13

Direct simulations of turbulent flow using finite-difference schemes

A high-order accurate finite-difference approach is presented for calculating incompressible turbulent flow. The methods used include a kinetic energy conserving central difference scheme and an upwind difference scheme. The methods are evaluated in test cases for the evolution of small-amplitude disturbances and fully developed turbulent channel flow. It is suggested that the finite-difference approach can be applied to complex geometries more easilty than highly accurate spectral methods. It is concluded that the upwind scheme is a good candidate for direct simulations of turbulent flows over complex geometries.

Rai, Man Mohan↗

Verification of Eulerian-Eulerian and Eulerian-Lagrangian simulations for turbulent fluid-particle flows

Here, we present a verification study of three simulation techniques for fluid–particle flows, including an Euler–Lagrange approach (EL) inspired by Jackson's seminal work on fluidized particles, a quadrature–based moment method based on the anisotropic Gaussian closure (AG), and the traditional two-fluid model. We perform simulations of two problems: particles in frozen homogeneous isotropic turbulence (HIT) and cluster-induced turbulence (CIT). For verification, we evaluate various techniques for extracting statistics from EL and study the convergence properties of the three methods under grid refinement. The convergence is found to depend on the simulation method and on the problem, with CIT simulations posing fewer difficulties than HIT. Specifically, EL converges under refinement for both HIT and CIT, but statistics exhibit dependence on the postprocessing parameters. For CIT, AG produces similar results to EL. For HIT, converging both TFM and AG poses challenges. Overall, extracting converged, parameter-independent Eulerian statistics remains a challenge for all methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Large Eddy Simulations of Turbulent and Buoyant Flows in Urban and Complex Terrain Areas Using the Aeolus Model

Fast and accurate predictions of the flow and transport of materials in urban and complex terrain areas are challenging because of the heterogeneity of buildings and land features of different shapes and sizes connected by canyons and channels, which results in complex patterns of turbulence that can enhance material concentrations in certain regions. To address this challenge, we have developed an efficient three-dimensional computational fluid dynamics (CFD) code called Aeolus that is based on first principles for predicting transport and dispersion of materials in complex terrain and urban areas. The model can be run in a very efficient Reynolds average Navier–Stokes (RANS) mode or a detailed large eddy simulation (LES) mode. The RANS version of Aeolus was previously validated against field data for tracer gas and radiological dispersal releases. As a part of this work, we have validated the Aeolus model in LES mode against two different sets of data: (1) turbulence quantities measured in complex terrain at Askervein Hill; and (2) wind and tracer data from the Joint Urban 2003 field campaign for urban topography. As a third set-up, we have applied Aeolus to simulate cloud rise dynamics for buoyant plumes from high-temperature explosions. For all three cases, Aeolus LES predictions compare well to observations and other models. These results indicate that Aeolus LES can be used to accurately simulate turbulent flow and transport for a wide range of applications and scales.

42 ENGINEERING↗

Simulation of turbulent transonic separated flow over an airfoil

A code developed for simulating high Reynolds number transonic flow fields of arbitrary configuration is described. This code, in conjunction with laboratory experiments, is used to devise and test turbulence transport models which may be suitable in the prediction of such flow fields, with particular emphasis on regions of flow separation. The solutions describe the flow field, including both the shock-induced and trailing-edge separation regions, in sufficient detail to provide the profile and friction drag.

Deiwert, G. S.↗

Validation of new and existing methods for time-domain simulations of turbulence and loads

We seek to obtain a second-by-second match between the simulated and measured structural loads of a utility-scale wind turbine. To obtain the one-to-one load simulations, we start with the furthest upstream component of the modeling chain: the turbulent inflow. We consider new and existing methods to generate constrained-turbulence flow fields. The new method is based on large-eddy simulations (LES) and machine learning (ML). The existing methods include Kaimal-based TurbSim and the superstatistical wind field model. The inflow measurements used to constrain these simulations are obtained with a nacelle-mounted scanning lidar. We compare the flow fields for the different inflow simulation approaches and validate their associated load predictions against measurements collected in the Rotor Aero-dynamics, Aeroelastics, and Wake (RAAW) field campaign. We find that the rotor-position control developed for this study is key in enabling the time match between measurements and simulations. When this control approach is used, the load simulation performance tracks with the inflow simulation fidelity, with LES+ML yielding errors ≤ 4% for the damage-equivalent loads of flapwise bending moment, and tower fore-aft bending moments.

17 WIND ENERGY↗

Direct simulation of turbulent combustion

Understanding and modeling of turbulent combustion are key-problems in the computation of numerous practical systems. Because of the lack of analytical theories in this field and of the difficulty of performing precise experiments, direct simulation appears to be one of the most attractive tools to use in addressing this problem. The present work can be split into two parts: (1) Development and validation of a direct simulation method for turbulent combustion; (2) Applications of the method to premixed turbulent combustion problems. The goal of part 1 is to define and to test a numerical method for direct simulation of reacting flows. A high level of confidence should be attached to direct simulation results, and this can only be achieved through extensive validation tests. In part 2, direct simulation is used to address some of the many critical problems related to turbulent combustion. At the present time, I have limited this work to premixed combustion and considered only four basic issues: (1) The effect of pressure waves on flame propagation; (2) The interaction between flame fronts and vortices; (3) The influence of curvature on premixed flame fronts; and (4) The validation of flamelet models for premixed turbulent combustion.

Poinsot, T. J.↗

Numerical simulation of turbulent flows around airfoil and wing

During the last years the simulation of compressible viscous flows has received much attention. While the numerical methods were improved drastically, a satisfactory modeling of the Reynolds stresses is still missing. In this paper, after a short description of the numerical procedure used for solving the Reynolds equations, experiments with a promising simple turbulence model are discussed.

Marx, Yves P.↗

Numerical simulations of turbulent molecular clouds regulated by radiation feedback forces. I. Star formation rate and efficiency

Radiation feedback from stellar clusters is expected to play a key role in setting the rate and efficiency of star formation in giant molecular clouds. To investigate how radiation forces influence realistic turbulent systems, we have conducted a series of numerical simulations employing the Hyperion radiation hydrodynamics solver, considering the regime that is optically thick to ultraviolet and optically thin to infrared radiation. Our model clouds cover initial surface densities between Σ{sub cl,0}∼10--300 M{sub ⊙} pc{sup −2}, with varying initial turbulence. We follow them through turbulent, self-gravitating collapse, star cluster formation, and cloud dispersal by stellar radiation. All our models display a log-normal distribution of gas surface density Σ; for an initial virial parameter α{sub vir,0}=2, the log-normal standard deviation is σ{sub lnΣ}=1--1.5 and the star formation rate coefficient ε{sub ff, ρ-bar }=0.3--0.5, both of which are sensitive to turbulence but not radiation feedback. The net star formation efficiency (SFE) ε{sub final} increases with Σ{sub cl,0} and decreases with α{sub vir,0}. We interpret these results via a simple conceptual framework, whereby steady star formation increases the radiation force, such that local gas patches at successively higher Σ become unbound. Based on this formalism (with fixed σ{sub lnΣ}), we provide an analytic upper bound on ε{sub final}, which is in good agreement with our numerical results. The final SFE depends on the distribution of Eddington ratios in the cloud and is strongly increased by the turbulent compression of gas.

79 ASTRONOMY AND ASTROPHYSICS↗

Direct numerical simulation of turbulent flow in a rotating square duct

A fully developed turbulent flow in a rotating straight square duct is simulated by direct numerical simulations at Re{sub τ} = 300 and 0 ≤ Ro{sub τ} ≤ 40. The rotating axis is parallel to two opposite walls of the duct and normal to the main flow. Variations of the turbulence statistics with the rotation rate are presented, and a comparison with the rotating turbulent channel flow is discussed. Rich secondary flow patterns in the cross section are observed by varying the rotation rate. The appearance of a pair of additional vortices above the pressure wall is carefully examined, and the underlying mechanism is explained according to the budget analysis of the mean momentum equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Reduced-model scrape-off layer turbulence (nSOLT) simulations comparing three fueling scenarios

The 2D scrape-off-layer turbulence code (nSOLT) includes 1D Boltzmann neutral–plasma interactions, a model of divertor recycling (introduced here), and a fixed source of plasma concentrated at the core-side boundary. Three fueling methods are considered herein: (1) neutral injection in the far-SOL is accomplished by specifying the density of Franck–Condon distributed neutrals streaming in from the boundary. (2) Divertor recycling is modeled by injecting a fraction of the particle parallel flux in the scrape-off layer (SOL) back into the edge region as a source of plasma. (3) A constant source fuels the edge plasma from the core-side boundary to model pellet injection. For machine parameters (B, R m , and L // ) illustrative of the MAST-U device, and for a deuterium plasma, turbulent equilibria are obtained that share the same plasma fueling rate for each of the three fueling methods, with only one of the sources on in each case. In the presence of self-consistent turbulence, quasi-steady plasma and neutral (deuterium) profiles, fueling efficiencies, SOL transparencies, and heat flux widths are compared. Characteristics of the turbulent fluctuations, including skewness, cross-phases, and power spectra, are described. The calculated fueling efficiencies, SOL transparencies to neutral penetration, and many of the turbulent properties are remarkably similar for all three fueling methods despite significant differences in the plasma profiles. The nonlinear states of the three cases are dominated by separatrix-spanning vortex cells that control particle and heat losses into the SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Implementation of turbulence models into simulators

Simulation of turbulence as it relates to the flight training environment is discussed in general terms. Simulators that put random white noise into the system and simulators that put random motion into the equations of motion are discussed, as are simulators that incorporate pitch and roll moments into vertical turbulence. Wind shear models and simplified models of microburst phenomena are covered.

Ireland, Robert L.↗

Steady-State Gyrokinetics Transport Code (SSGKT), A Scientific Application Partnership with the Framework Application for Core-Edge Transport Simulations, Final Report

This project initiated the development of TGYRO - a steady-state Gyrokinetic transport code (SSGKT) that integrates micro-scale GYRO turbulence simulations into a framework for practical multi-scale simulation of conventional tokamaks as well as future reactors. Using a lightweight master transport code, multiple independent (each massively parallel) gyrokinetic simulations are coordinated. The capability to evolve profiles using the TGLF model was also added to TGYRO and represents a more typical use-case for TGYRO. The goal of the project was to develop a steady-state Gyrokinetic transport code (SSGKT) that integrates micro-scale gyrokinetic turbulence simulations into a framework for practical multi-scale simulation of a burning plasma core ? the International Thermonuclear Experimental Reactor (ITER) in particular. This multi-scale simulation capability will be used to predict the performance (the fusion energy gain, Q) given the H-mode pedestal temperature and density. At present, projections of this type rely on transport models like GLF23, which are based on rather approximate fits to the results of linear and nonlinear simulations. Our goal is to make these performance projections with precise nonlinear gyrokinetic simulations. The method of approach is to use a lightweight master transport code to coordinate multiple independent (each massively parallel) gyrokinetic simulations using the GYRO code. This project targets the practical multi-scale simulation of a reactor core plasma in order to predict the core temperature and density profiles given the H-mode pedestal temperature and density. A master transport code will provide feedback to O(16) independent gyrokinetic simulations (each massively parallel). A successful feedback scheme offers a novel approach to predictive modeling of an important national and international problem. Success in this area of fusion simulations will allow US scientists to direct the research path of ITER over the next two decades. The design of an efficient feedback algorithm is a serious numerical challenge. Although the power source and transport balance coding in the master are standard, it is nontrivial to design a feedback loop that can cope with outputs that are both intermittent and extremely expensive. A prototypical feedback scheme has already been successfully demonstrated for a single global GYRO simulation, although the robustness and efficiency are likely far from optimal. Once the transport feedback scheme is perfected, it could, in principle, be embedded into any of the more elaborate transport codes (ONETWO, TRANSP, and CORSICA), or adopted by other FSP-related multi-scale projects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulations of Turbulent Momentum and Scalar Transport in Confined Swirling Coaxial Jets

This paper presents the numerical simulations of confined three dimensional coaxial water jets. The objectives are to validate the newly proposed nonlinear turbulence models of momentum and scalar transport, and to evaluate the newly introduced scalar APDF and DWFDF equation along with its Eulerian implementation in the National Combustion Code (NCC). Simulations conducted include the steady RANS, the unsteady RANS (URANS), and the time-filtered Navier-Stokes (TFNS) with and without invoking the APDF or DWFDF equation. When the APDF or DWFDF equation is invoked, the simulations are of a hybrid nature, i.e., the transport equations of energy and species are replaced by the APDF or DWFDF equation. Results of simulations are compared with the available experimental data. Some positive impacts of the nonlinear turbulence models and the Eulerian scalar APDF and DWFDF approach are observed.

Shih, Tsan-Hsing↗

Efficient Parallel Algorithm For Direct Numerical Simulation of Turbulent Flows

A distributed algorithm for a high-order-accurate finite-difference approach to the direct numerical simulation (DNS) of transition and turbulence in compressible flows is described. This work has two major objectives. The first objective is to demonstrate that parallel and distributed-memory machines can be successfully and efficiently used to solve computationally intensive and input/output intensive algorithms of the DNS class. The second objective is to show that the computational complexity involved in solving the tridiagonal systems inherent in the DNS algorithm can be reduced by algorithm innovations that obviate the need to use a parallelized tridiagonal solver.

Moitra, Stuti↗

Numerical Simulation of Turbulent Jets with Rectangular Cross-Section

Three-dimensional turbulent jets with rectangular cross-section are simulated with a finite-difference numerical method. The full Navier-Stokes equations are solved at low Reynolds numbers, whereas at the high Reynolds numbers filtered forms of the equations are solved along with a sub-grid scale model to approximate effects of the unresolved scales. A 2-N storage, third-order Runge-Kutta scheme is used for temporal discretization and a fourth-order compact scheme is used for spatial discretization. Computations are performed for different inlet conditions which represent different types of jet forcing. The phenomenon of axis-switching is observed, and it is confirmed that this is based on self-induction of the vorticity field. Budgets of the mean streamwise velocity show that convection is balanced by gradients of the Reynolds stresses and the pressure.

Wilson, Robert V.↗

Numerical simulation of turbulent flows

Computational models of turbulence in incompressible Newtonian fluids governed by the Navier-Stokes equations are reviewed. The governing equations are presented, and both direct and large-eddy-simulations are examined. Resolution requirements and numerical techniques of spatial representation, definition of initial and boundary conditions, and time advancement are considered. Results of simulations of homogeneous turbulence in uniform shear, the evolution of a turbulent mixing layer, and turbulent channel flow are presented graphically and discussed.

Rogallo, R. S.↗

Efficient implicit LES method for the simulation of turbulent cavitating flows

We present a numerical method for efficient large-eddy simulation of compressible liquid flows with cavitation based on an implicit subgrid-scale model. Phase change and subgrid-scale interface structures are modeled by a homogeneous mixture model that assumes local thermodynamic equilibrium. Unlike previous approaches, emphasis is placed on operating on a small stencil (at most four cells). The truncation error of the discretization is designed to function as a physically consistent subgrid-scale model for turbulence. We formulate a sensor functional that detects shock waves or pseudo-phase boundaries within the homogeneous mixture model for localizing numerical dissipation. In smooth regions of the flow field, a formally non-dissipative central discretization scheme is used in combination with a regularization term to model the effect of unresolved subgrid scales. The new method is validated by computing standard single- and two-phase test-cases. Comparison of results for a turbulent cavitating mixing layer obtained with the new method demonstrates its suitability for the target applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗