Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical 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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Direct numerical simulations of fluid flow, heat transfer and phase changes

Direct numerical simulations of fluid flow, heat transfer, and phase changes are presented. The simulations are made possible by a recently developed finite difference/front tracking method based on the one-field formulation of the governing equations where a single set of conservation equations is written for all the phases involved. The conservation equations are solved on a fixed rectangular grid, but the phase boundaries are kept sharp by tracking them explicitly by a moving grid of lower dimension. The method is discussed and applications to boiling heat transfer and the solidification of drops colliding with a wall are shown.

Phase Transformations↗
NASA NTRS2015

The Importance of the Cathode Plume and Its Interactions with the Ion Beam in Numerical Simulations of Hall Thrusters

Hall2De is a first-principles, 2-D axisymmetric code that solves the equations of motion for ions, electrons, and neutrals on a magnetic-field-aligned grid. The computational domain downstream of the acceleration channel exit plane is large enough to include self-consistently the cathode boundary. In this paper, we present results from numerical simulations of the H6 laboratory thruster with an internally mounted cathode, with the aim of highlighting the importance of properly accounting for the interactions between the ion beam and cathode plume. The anomalous transport of electrons across magnetic field lines in Hall2De is modelled using an anomalous collision frequency, νanom, yielding νanom approximately equal to omega ce (i.e., the electron cyclotron frequency) in the plume. We first show that restricting the anomalous collision frequency to only regions where the current density of ions is large does not alter the plasma discharge in the Hall thruster as long as the interaction between the ion beam and the cathode plume is captured properly in the computational domain. This implies that the boundary conditions must be placed sufficiently far as to not interfere with the electron transport in this region. These simulation results suggest that electron transport across magnetic field lines occurs largely inside the beam and may be driven by the interactions between beam ions and electrons. A second finding that puts in relevance the importance of including the cathode plume in numerical simulations is on the significance of accounting for the ion acoustic turbulence (IAT), now known to occur in the vicinity of the cathode exit. We have included in the Hall2De simulations a model of the IAT-driven anomalous collision frequency based on Sagdeev's model for saturation of the ion-acoustic instability. This implementation has allowed us to achieve excellent agreement with experimental measurements in the near plume obtained during the operation of the H6 thruster at nominal conditions (300V, 20A) and chamber background pressure of approximately 1.5 x 10(exp -5) Torr. In addition, the numerical results obtained with the latter approach exhibit less sensitivity to background pressure than previous attempts at explaining the features of the plasma properties in the near plume.

Lopez Ortega, Alejandro↗
DOE OSTI2018

Investigation of Turbulent Entrainment-Mixing Processes With a New Particle-Resolved Direct Numerical Simulation Model

Here, a new particle-resolved three dimensional direct numerical simulation (DNS) model is developed that combines Lagrangian droplet tracking with the Eulerian field representation of turbulence near the Kolmogorov microscale. Six numerical experiments are performed to investigate the processes of entrainment of clear air and subsequent mixing with cloudy air and their interactions with cloud microphysics. The experiments are designed to represent different combinations of three configurations of initial cloudy area and two turbulence modes (decaying and forced turbulence). Five existing measures of microphysical homogeneous mixing degree are examined, modified, and compared in terms of their ability as a unifying measure to represent the effect of various entrainment-mixing mechanisms on cloud microphysics. Also examined and compared are the conventional Damköhler number and transition scale number as a dynamical measure of different mixing mechanisms. Relationships between the various microphysical measures and dynamical measures are investigated in search for a unified parameterization of entrainment-mixing processes. The results show that even with the same cloud water fraction, the thermodynamic and microphysical properties are different, especially for the decaying cases. Further analysis confirms that despite the detailed differences in cloud properties among the six simulation scenarios, the variety of turbulent entrainment-mixing mechanisms can be reasonably represented with power-law relationships between the microphysical homogeneous mixing degrees and the dynamical measures.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Direct Numerical Simulation of Particle Behaviour in a Gas-Solid Three Dimensional Plane Jet

In this paper, direct numerical simulations (DNS) of a three-dimensional (3D), non-reacting, temporally evolving planar jet laden with mono-dispersed solid particles in the two-way coupling (TWC) regime are performed. Three different particles Stokes numbers (St = 0.1, 1, 10) have been considered. This has been achieved by varying the particle diameter while keeping the particle mass loading (fm = 1) and the jet Reynolds number (Rejet = 2000) unchanged. The objective is to study the effect of the particle Stokes number TWC regime on the temporal development of the planar jet. Two-way coupled momentum and heat transfer has been studied by investigating mean relative velocity and temperature. Results indicate that the relative parameters are more pronounced on the edges of the jet and decrease in time in general. At the center of the jet however, the mean value first increases and then decreases again. Additionally, lighter particles spread farther than heavier particles from the center of the jet. Furthermore, the heavier particles delay the development of the jet due to TWC effects.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Critical comparison of second-order closures with direct numerical simulations of homogeneous turbulence

Recently, several models have been proposed for closing the second-moment equations, in which the velocity-pressure gradient tensor and the dissipation rate tensor are two of the most important terms. In the literature, these correlation tensors are usually decomposed into a so-called rapid term and a return-to-isotropy term. Models of these terms have been used in global flow calculations together with other modeled terms. However, their individual behaviors in different flows have not been fully examined because they are unmeasurable in the laboratory. Recently, the development of direct numerical simulation (DNS) of turbulence has given us the possibility to do this kind of study. With direct numerical simulation, we may use the solution to calculate exactly the values of these correlation terms and then directly compare them with the values from their modeled formulations. In this paper, we make direct comparisons of five representative rapid models and eight return-to-isotropy models using the DNS data of 45 homogeneous flows, which were done by Rogers et al. (1986) and Lee and Reynolds (1985). The purpose of these direct comparisons is to explore the performance of these models in different flows and identify the ones that give the best performance. The paper also describes the modeling procedure, model constraints, and the various evaluated models. The detailed results of the direct comparisons are discussed, and a few concluding remarks on turbulence models are given.

Shih, Tsan-Hsing↗
NASA NTRS1986

A numerical simulation of the inviscid flow through a counterrotating propeller

The results of a numerical simulation of the time-averaged inviscid flow field through the blade rows of a multiblade row turboprop configuration are presented. The governing equations are outlined along with a discussion of the solution procedure and coding strategy. Numerical results obtained from a simulation of the flow field through a modern high-speed turboprop will be shown.

Celestina, M. L.↗
NASA NTRS1983

A numerical simulation of transition in plane channel flow

This paper involves a numerical simulation of the final stages of transition to turbulence in plane channel flow at a Reynolds number of 7500. Three-dimensional, incompressible Navier-Stokes equations are numerically integrated to obtain the time-evolution of two- and three-dimensional finite-amplitude disturbances. Computations are performed on the CYBER-203 vector processor for a 32 x 33 x 32 grid. Solutions indicate the existence of structures similar to those observed in the laboratory and which are characteristic of various stages of transition that lead to final breakdown. Details of the resulting flow field after breakdown indicate the evolution of streak-like formations found in turbulent flows. Although the flow field does approach a steady-state (turbulent channel flow), implementation of subgrid-scale terms are necesary to obtain proper turbulent statistics.

Biringen, S.↗
NASA NTRS1991

A critical comparison of second order closures with direct numerical simulation of homogeneous turbulence

Recently, several second order closure models have been proposed for closing the second moment equations, in which the velocity-pressure gradient (and scalar-pressure gradient) tensor and the dissipation rate tensor are two of the most important terms. In the literature, these correlation tensors are usually decomposed into a so called rapid term and a return-to-isotropy term. Models of these terms have been used in global flow calculations together with other modeled terms. However, their individual behavior in different flows have not been fully examined because they are un-measurable in the laboratory. Recently, the development of direct numerical simulation (DNS) of turbulence has given us the opportunity to do this kind of study. With the direct numerical simulation, we may use the solution to exactly calculate the values of these correlation terms and then directly compare them with the values from their modeled formulations (models). Here, we make direct comparisons of five representative rapid models and eight return-to-isotropy models using the DNS data of forty five homogeneous flows which were done by Rogers et al. (1986) and Lee et al. (1985). The purpose of these direct comparisons is to explore the performance of these models in different flows and identify the ones which give the best performance. The modeling procedure, model constraints, and the various evaluated models are described. The detailed results of the direct comparisons are discussed, and a few concluding remarks on turbulence models are given.

Shih, Tsan-Hsing↗
DOE OSTI2023

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part I: Methodology and Challenges

Abstract High-fidelity simulation of transitional and turbulent flows over multi-scale surface roughness presents several challenges. For instance, the complex and irregular geometrical nature of surface roughness makes it impractical to employ conforming structured grids, commonly adopted in large-scale numerical simulations due to their high computational efficiency. One possible solution to overcome this problem is offered by immersed boundary methods, which allow wall boundary conditions to be enforced on grids that do not conform to the geometry of the solid boundary. To this end, a three-dimensional, second-order accurate boundary data immersion method (BDIM) is adopted. A novel mapping algorithm that can be applied to general three-dimensional surfaces is presented, together with a newly developed data-capturing methodology to extract and analyze on-surface flow quantities of interest. A rigorous procedure to compute gradient quantities such as the wall shear stress and the heat flux on complex non-conforming geometries is also introduced. The new framework is validated by performing a direct numerical simulation (DNS) of fully developed turbulent channel flow over sinusoidal egg-carton roughness in a minimal-span domain. For this canonical case, the averaged streamwise velocity profiles are compared against results from the literature obtained with a body-fitted grid. General guidelines on the BDIM resolution requirements for multi-scale roughness simulation are given. Momentum and energy balance methods are used to validate the calculation of the overall skin friction and heat transfer at the wall. The BDIM is then employed to investigate the effect of irregular homogeneous surface roughness on the performance of an LS-89 high-pressure turbine blade at engine-relevant conditions using DNS. This is the first application of the BDIM to realize multi-scale roughness for transitional flow in transonic conditions in the context of high-pressure turbines. The methodology adopted to generate the desired roughness distribution and to apply it to the reference blade geometry is introduced. The results are compared to the case of an equivalent smooth blade.

Engineering↗
NASA NTRS2014

Numerical Simulation of a High-Lift Configuration with Embedded Fluidic Actuators

Numerical simulations have been performed for a vertical tail configuration with deflected rudder. The suction surface of the main element of this configuration is embedded with an array of 32 fluidic actuators that produce oscillating sweeping jets. Such oscillating jets have been found to be very effective for flow control applications in the past. In the current paper, a high-fidelity computational fluid dynamics (CFD) code known as the PowerFLOW(Registered TradeMark) code is used to simulate the entire flow field associated with this configuration, including the flow inside the actuators. The computed results for the surface pressure and integrated forces compare favorably with measured data. In addition, numerical solutions predict the correct trends in forces with active flow control compared to the no control case. Effect of varying yaw and rudder deflection angles are also presented. In addition, computations have been performed at a higher Reynolds number to assess the performance of fluidic actuators at flight conditions.

Vatsa, Veer N.↗
DOE OSTI2018

Direct numerical simulation of flame stabilization assisted by autoignition in a reheat gas turbine combustor

A three-dimensional direct numerical simulation (DNS) is performed for a turbulent hydrogen-air flame, represented with detailed chemistry, stabilized in a model gas-turbine combustor. The combustor geometry consists of a mixing duct followed by a sudden expansion and a combustion chamber, which represents a geometrically simplified version of Ansaldo Energia’s GT26/GT36 sequential combustor design. In this configuration, a very lean blend of hydrogen and vitiated air is prepared in the mixing duct and convected into the combustion chamber, where the residence time from the inlet of the mixing duct to the combustion chamber is designed to coincide with the ignition delay time of the mixture. The results show that when the flame is stabilized at its design position, combustion occurs due to both autoignition and flame propagation (deflagration) modes at different locations within the combustion chamber. A chemical explosive mode analysis (CEMA) reveals that most of the fuel is consumed due to autoignition in the bulk-flow along the centerline of the combustor, and lower amounts of fuel are consumed by flame propagation near the corners of the sudden expansion, where the unburnt temperature is reduced by the thermal wall boundary layers. An unstable operating condition is also identified, wherein periodic auto-ignition events occur within the mixing duct. These events appear upstream of the intended stabilization position, due to positive temperature fluctuations induced by pressure waves originating from within the combustion chamber. This DNS investigation represents the initial step of a comprehensive research effort aimed at gaining detailed physical insight into the rate-limiting processes that govern the sequential combustor behavior and avoid the insurgence of the off-design auto-ignition events.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2015

Numerical Simulations of Subscale Wind Turbine Rotor Inboard Airfoils at Low Reynolds Number

New blade designs are planned to support future research campaigns at the SWiFT facility in Lubbock, Texas. The sub-scale blades will reproduce specific aerodynamic characteristics of utility-scale rotors. Reynolds numbers for megawatt-, utility-scale rotors are generally above 2-8 million. The thickness of inboard airfoils for these large rotors are typically as high as 35-40%. The thickness and the proximity to three-dimensional flow of these airfoils present design and analysis challenges, even at the full scale. However, more than a decade of experience with the airfoils in numerical simulation, in the wind tunnel, and in the field has generated confidence in their performance. Reynolds number regimes for the sub-scale rotor are significantly lower for the inboard blade, ranging from 0.7 to 1 million. Performance of the thick airfoils in this regime is uncertain because of the lack of wind tunnel data and the inherent challenge associated with numerical simulations. This report documents efforts to determine the most capable analysis tools to support these simulations in an effort to improve understanding of the aerodynamic properties of thick airfoils in this Reynolds number regime. Numerical results from various codes of four airfoils are verified against previously published wind tunnel results where data at those Reynolds numbers are available. Results are then computed for other Reynolds numbers of interest.

17 WIND ENERGY↗
NASA NTRS2002

Numerical Simulation of One- And Two-Phase Flows In Propulsion Systems

In this report, we present some results of problems investigated during joint research between the Hampton University (HU) Fluid Mechanics and Acoustics Laboratory (FM&AL), NASA Glenn Research Center (GRC) and the Hyper-X Program of the NASA Langley Research Center (LaRC). This work is supported by joint research between the NASA GRC/HU FM&AL and the Institute of Mechanics at Moscow State University (IM/MSU) in Russia under a Civilian Research and Development Foundation (CRDF) grant, #RE1-2068. The main areas of current scientific interest of the FM&AL include an investigation of the proposed and patented advanced methods for aircraft engine thrust and noise benefits. These methods are based on nontraditional 3D (three dimensional) corrugated and composite nozzle, inlet, propeller and screw designs such as the Bluebell and Telescope nozzles, Mobius-shaped screws, etc. These are the main subject of our other projects, of which one is the NASA MURED's (Minority University Research and Education Division) FAR (Faculty Awards for Research) Award, #NAG-3-2249. Working jointly with this project team, our team also analyzes additional methods for exhaust jet noise reduction. These methods are without essential thrust loss and even with thrust augmentation. The research is focused on a wide regime of problems in the propulsion field as well as in experimental testing and theoretical and numerical simulation analyses for advanced aircraft and rocket engines. The FM&AL Team uses analytical methods, numerical simulations and experimental tests at the Hampton University campus, NASA and IM/MSU. The main results obtained by FM&AL team were published in the papers and patents.

Gilinsky, Mikhail M.↗
NASA NTRS1998

Numerical Simulation of the Generation of Axisymmetric Mode Jet Screech Tones

An imperfectly expanded supersonic jet, invariably, radiates both broadband noise and discrete frequency sound called screech tones. Screech tones are known to be generated by a feedback loop driven by the large scale instability waves of the jet flow. Inside the jet plume is a quasi-periodic shock cell structure. The interaction of the instability waves and the shock cell structure, as the former propagates through the latter, is responsible for the generation of the tones. Presently, there are formulas that can predict the tone frequency fairly accurately. However, there is no known way to predict the screech tone intensity. In this work, the screech phenomenon of an axisymmetric jet at low supersonic Mach number is reproduced by numerical simulation. The computed mean velocity profiles and the shock cell pressure distribution of the jet are found to be in good agreement with experimental measurements. The same is true with the simulated screech frequency. Calculated screech tone intensity and directivity at selected jet Mach number are reported in this paper. The present results demonstrate that numerical simulation using computational aeroacoustics methods offers not only a reliable way to determine the screech tone intensity and directivity but also an opportunity to study the physics and detailed mechanisms of the phenomenon by an entirely new approach.

Shen, Hao↗
DOE OSTI2024

Towards accelerating particle-resolved direct numerical simulation with neural operators

In this paper, we present our ongoing work aimed at accelerating a particle-resolved direct numerical simulation model designed to study aerosol–cloud–turbulence interactions. The dynamical model consists of two main components—a set of fluid dynamics equations for air velocity, temperature, and humidity, coupled with a set of equations for particle (i.e., cloud droplet) tracing. Rather than attempting to replace the original numerical solution method in its entirety with a machine learning (ML) method, we consider developing a hybrid approach. We exploit the potential of neural operator learning to yield fast and accurate surrogate models and, in this study, develop such surrogates for the velocity and vorticity fields. We discuss results from numerical experiments designed to assess the performance of ML architectures under consideration as well as their suitability for capturing the behavior of relevant dynamical systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Numerical simulation of evaporating wavy falling liquid films in laminar gas streams

Numerical simulations are performed to investigate the interfacial heat and mass transfer from evaporating wavy falling liquid films in interaction with laminar gas streams. The OpenFOAM solver has been used to conduct the simulations where the liquid-gas interface is resolved using the Volume of Fluid method of solving for three phases, i.e., liquid, vapor, and air. The configuration considered is a falling liquid (water) film on a heated vertical plate with a confined laminar moist air (gas) flow that is either (a) co-current or (b) counter-current to the downward liquid flow. The evaporation at the liquid-gas interface is driven by the interfacial gradient of the vapor mass fraction. Interfacial waves are triggered using a monochromatic forcing disturbance that leads to sinusoidal or solitary waves forming at the liquid-gas interface under respective forcing frequencies. Further, the numerical model is validated with the available experimental data. The results show nearly a 15% enhancement in time-averaged Sherwood number (Sh) due to film waviness (sinusoidal or solitary) at the lower volumetric gas flow rate, Q g = +50 (co-current) and Q g = -50 (counter-current). This enhancement in the Sh for both the waves further increases by 11% with Q g = +800 and 196% with Q g = -800. A closer examination of the mass transfer process over a wave demonstrates that with Q g = +50, the concentration of the gas side streamlines at the trough locations of the wave leads to higher values of Sh at these locations. However, with Q g = +800, although the overall Sh increases, vortices appear at the wave trough locations, leading to a corresponding decrease in the local Sh values. Correlations are proposed for predicting Sh under co-current and counter-current gas flow effects.

36 MATERIALS SCIENCE↗
NASA NTRS1980

Numerical simulation of the final stages of terrestrial planet formation

Three representative numerical simulations of the growth of the terrestrial planets by accretion of large protoplanets are considered. The mass and relative-velocity distributions of the bodies are free to evolve simultaneously in response to close gravitational encounters and occasional collisions between bodies. The collisions between bodies arise therefore in a natural way and the assumption of expressions for the relative-velocity distribution and the gravitational collision cross section is unnecessary. These simulations indicate that the growth of bodies with final masses approaching those of Venus and earth is possible, at least for the case of a two-dimensional system

Cox, L. P.↗