Search NASA⌕ Search

SEARCH · Search NASA

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 289 records · Page 16

Direct Numerical Simulation of Receptivity to Roughness in a Swept-Wing Boundary Layer at High Reynolds Numbers

Direct numerical simulations (DNS) are performed to examine the receptivity to roughness in a spatially developing three-dimensional boundary layer over an in finite-swept natural-laminar-flow wing at a free stream Mach number of 0:75 and a chord Reynolds number of approximately 25 million based on the long, swept chord. Stationary cross ow disturbances are excited by applying either critically spaced discrete cylinders of micron size or naturally occurring distributed roughness in the leading-edge region. The DNS data show that the spanwise spectral content of the excited cross ow disturbances is highly dependent upon the shape of roughness elements, and the initial growth of the cross ow structures is a nonlinear function of the element height. The linear growth rate of the excited cross ow disturbances predicted by DNS shows good agreement with linear parabolized stability equations. The receptivity study lays the foundation for investigating the stabilization of the naturally most unstable steady cross ow mode by using spanwise periodic DREs.

Nicholson, Gary↗

Large eddy numerical simulations of turbulent flows

Large eddy simulations are a numerical technique in which large scale turbulent structures are explicitly computed and the small structures are modelled. Arguments for believing this method to be superior to more conventional approaches are given, the basis of the method is given, and some typical results displayed. The results show that the method does have enormous promise, but much further development is required.

Ferziger, J. H.↗

Numerical simulation in three space dimensions of time-dependent thermal convection in a rotating fluid

Three-dimensional time-dependent convection in a plane layer of fluid, uniformly heated from below and subject to vertical shear and to rotation about an axis tilted from the vertical, was simulated by the numerical solution of the Boussinesq equations, including all Coriolis terms. Rotation about a vertical axis produces smaller convection cells with diminished heat fluxes and considerable vorticity. When the rotation axis is tilted from the vertical to represent tropical latitudes, the convection cells become elongated in a N-S direction. Imposed flows with constant vertical shear produce convective rolls aligned with the mean flow. When the rotation vector is tilted from the vertical, the competing effects due to rotation and shear can stabilize the convective motions.

Hathaway, D. H.↗

Numerical simulation methods of incompressible flows and an application to the Space Shuttle main engine

This paper discusses incompressible Navier-Stokes solution methods with an emphasis on the pseudocompressibility method. A steady-state flow solver based on the pseudocompressibility approach is then described. This flow solver code has been used to analyze the internal flow in the Space Shuttle main engine hot-gas manifold. Salient features associated with this three-dimensional realistic flow simulation are discussed. Numerical solutions relevant to the current engine analysis and the redesign effort are discussed along with experimental results. This example demonstrates the potential of computational fluid dynamics as a design tool for aerospace applications.

Chang, J. L. C.↗

Numerical simulation methods of incompressible flows and an application to the space shuttle main engine

Incompressible Navier-Stokes solution methods are discussed with an emphasis on the pseudocompressibility method. A steady-state flow solver based on the pseudocompressibility approach is then described. This flow-solver code was used to analyze the internal flow in the Space Shuttle main engine hot-gas manifold. Salient features associated with this three-dimensional realistic flow simulation are discussed. Numerical solutions relevant to the current engine analysis and the redesign effort are discussed along with experimental results. This example demonstrates the potential of computational fluid dynamics as a design tool for aerospace applications.

Chang, J. L. C.↗

Studying Turbulence Using Numerical Simulation Databases - IX: Proceedings of the 2002 Summer Program

The ninth Summer Program of the Center for Turbulence Research was held during the period July 29th - August 23rd, 2002. The increase in number of participants, noted in the Preface to the Proceedings of the 2000 Program, continues: this year there were 50 participants from ten countries, and 30 hosts from Stanford and NASA-Ames. This Proceedings volume contains 32 papers that span a wide range of topics and an enormous range of physical scales. The papers have been divided into seven groups: Acoustics, RANS modeling, Combustion, Large-eddy simulation (LES), LES Numerics, Stratified Flows, and Fundamentals, In several cases, a paper could have fitted in more than one group so the classification is somewhat arbitrary.

Bradshaw, Peter↗

Numerical simulation of transonic propeller flow using a three-dimensional small disturbance code employing novel helical coordinates

The numerical simulation of three-dimensional transonic flow about propeller blades is discussed. The equations for the unsteady potential flow about propellers is given for an arbitrary coordinate system. From this the small disturbance form of the equation is derived for a new helical coordinate system. The new coordinate system is suited to propeller flow and allows cascade boundary conditions to be applied straightforward. A numerical scheme is employed which solves the steady flow as an asymptotic limit of unsteady flow. Solutions are presented for subsonic and transonic flow about a 5 percent thick bicircular arc blade of an eight bladed cascade. Both high and low advance ratio cases are given which include a lifting case as well as nonlifting cases. The nonlifting cases are compared to solutions from a Euler code.

Snyder, Aaron↗

The void spectrum in two-dimensional numerical simulations of gravitational clustering

An algorithm for deriving a spectrum of void sizes from two-dimensional high-resolution numerical simulations of gravitational clustering is tested, and it is verified that it produces the correct results where those results can be anticipated. The method is used to study the growth of voids as clustering proceeds. It is found that the most stable indicator of the characteristic void 'size' in the simulations is the mean fractional area covered by voids of diameter d, in a density field smoothed at its correlation length. Very accurate scaling behavior is found in power-law numerical models as they evolve. Eventually, this scaling breaks down as the nonlinearity reaches larger scales. It is shown that this breakdown is a manifestation of the undesirable effect of boundary conditions on simulations, even with the very large dynamic range possible here. A simple criterion is suggested for deciding when simulations with modest large-scale power may systematically underestimate the frequency of larger voids.

Kauffmann, Guinevere↗

Numerical Simulation of Lean Blowout of Alternative Fuels in 7-element Lean Direct Injector

This research presents the result of numerically simulating 7-element swirl-venturi Lean Direct Injector (SV-LDI) lean blowout (LBO) experiments conducted at NASA Glenn Research Center in May of 2019. After simulating a cold flow case to confirm the pressure drop agrees well with the experiment, additional cases with two different fuels (an average jet fuel and a Gevo alcohol-to-jet fuel) from the National Jet Fuels Combustion Program (NJFCP) were computed to numerically determine the LBO condition. The procedure to approach the LBO follows the method used in the experiment where the air mass flow rate is gradually increased while the fuel supply is maintained. Transient history of global heat release rate as a function of air flow rate is presented, as well as temperature contours at different conditions to give a visual representation of the flame state. The Open National Combustion Code (OpenNCC) used in this research adopted reduced HyChem (Hybrid Chemistry) models along with k-LES turbulence model and a Lagrangian spray model that takes into account droplet internal temperature distribution affected by the shear force on the droplet surface. The transport equations and chemical reaction terms are integrated together to enhance conservation of chemical species that are especially important in the near LBO conditions. After showing the computed range of LBO agrees well with the experimental measurements, time averaged solutions of both fuels at their initial condition and at their limiting condition just before LBO are compared in detail to facilitate the understanding of LBO mechanism.

combustion↗

Numerical simulations using a dynamic solution-adaptive grid algorithm, with applications to unsteady internal flows

An investigation into the numerical simulation of unsteady flows is undertaken using a two-stage Runge-Kutta scheme coupled with the dynamic solution-adaptive grid algorithm developed by the authors. The inviscid fluxes are described by a modified Advective Upwind Split Method to eliminate the need for artificial dissipation. A well-documented numerical example containing moving discontinuities is presented that demonstrates the ability of the coupled grid/solver scheme to accurately capture unsteady flowfield phenomena. Applications are to a typical inlet diffuser configuration at Mach 3.0 with excessive back pressure inducing inlet unstart.

Benson, Rusty A.↗

Numerical simulation of transition in a decelerating boundary layer

Transition in a decelerating flat-plate boundary layer is numerically simulated up to the beginning of three-dimensional breakdown, and the results are compared with an experiment. The adverse pressure gradient induced by deceleration increases the growth rate of disturbances and allows transition at lower Reynolds numbers. The primary instability is characterized by a wave packet, which undergoes three-dimensional distortion. Lambda vortices are locally observed, but they are not aligned with respect to the flow direction.

Yang, Kyung Soo↗

Numerical simulations - Some results for the 2- and 3-D Hubbard models and a 2-D electron phonon model

Numerical simulations on the half-filled three-dimensional Hubbard model clearly show the onset of Neel order. Simulations of the two-dimensional electron-phonon Holstein model show the competition between the formation of a Peierls-CDW state and a superconducting state. However, the behavior of the partly filled two-dimensional Hubbard model is more difficult to determine. At half-filling, the antiferromagnetic correlations grow as T is reduced. Doping away from half-filling suppresses these correlations, and it is found that there is a weak attractive pairing interaction in the d-wave channel. However, the strength of the pair field susceptibility is weak at the temperatures and lattice sizes that have been simulated, and the nature of the low-temperature state of the nearly half-filled Hubbard model remains open.

Scalapino, D. J.↗

High-Resolution Numerical Simulation Of Shock Waves

NASA technical memorandum compares results of upwind and symmetric shock-capturing methods in numerical simulation of gas-dynamic flows. Both methods find shocks as sharp variations in fluid properties over few grid points. Methods differ in type of artificial viscosity introduced to stabilize computations. Symmetric method shown as accurate as upwind method, but with fewer and simpler time steps in transient case or fewer iteration steps in steady-state case.

Yee, H. C.↗

Multigrid methods for numerical simulation of laminar diffusion flames

This paper documents the result of a computational study of multigrid methods for numerical simulation of 2D diffusion flames. The focus is on a simplified combustion model, which is assumed to be a single step, infinitely fast and irreversible chemical reaction with five species (C3H8, O2, N2, CO2 and H2O). A fully-implicit second-order hybrid scheme is developed on a staggered grid, which is stretched in the streamwise coordinate direction. A full approximation multigrid scheme (FAS) based on line distributive relaxation is developed as a fast solver for the algebraic equations arising at each time step. Convergence of the process for the simplified model problem is more than two-orders of magnitude faster than other iterative methods, and the computational results show good grid convergence, with second-order accuracy, as well as qualitatively agreement with the results of other researchers.

Liu, C.↗

Comparison of numerically-simulated microphysical characteristics of convective storms with multiparameter radar observations during COHMEX

The microphysical/dynamical nature and evolution of numerically simulated convective storms observed on July 11, 1986 during the Cooperative Huntsville Meteorological Experiments are investigated with the Goddard Cumulus Ensemble model. The model-simulated hydrometeor and radar reflectivity structure is verified with results derived from the multiparameter radar observations, which includes a description of the reflectivity, differential reflectivity, linear depolarization ratio, and hail signal for two thunderstorm cases analyzed by Fulton and Heymsfield (1990). Radar analysis shows the storm to split with the newer more vigorous cell in the SE quadrant. Echo tops are near 15 km while peak reflectivities exceed 50 dBZ. Preliminary model simulations likewise show a splitting storm with regeneration occurring on the SSE flank. Echo tops are around 14 km with peak reflectivities over 45 dBZ.

Tao, W.-K.↗

Numerical Simulation of Rocket Exhaust Interaction with Lunar Soil

This technology development originated from the need to assess the debris threat resulting from soil material erosion induced by landing spacecraft rocket plume impingement on extraterrestrial planetary surfaces. The impact of soil debris was observed to be highly detrimental during NASA s Apollo lunar missions and will pose a threat for any future landings on the Moon, Mars, and other exploration targets. The innovation developed under this program provides a simulation tool that combines modeling of the diverse disciplines of rocket plume impingement gas dynamics, granular soil material liberation, and soil debris particle kinetics into one unified simulation system. The Unified Flow Solver (UFS) developed by CFDRC enabled the efficient, seamless simulation of mixed continuum and rarefied rocket plume flow utilizing a novel direct numerical simulation technique of the Boltzmann gas dynamics equation. The characteristics of the soil granular material response and modeling of the erosion and liberation processes were enabled through novel first principle-based granular mechanics models developed by the University of Florida specifically for the highly irregularly shaped and cohesive lunar regolith material. These tools were integrated into a unique simulation system that accounts for all relevant physics aspects: (1) Modeling of spacecraft rocket plume impingement flow under lunar vacuum environment resulting in a mixed continuum and rarefied flow; (2) Modeling of lunar soil characteristics to capture soil-specific effects of particle size and shape composition, soil layer cohesion and granular flow physics; and (3) Accurate tracking of soil-borne debris particles beginning with aerodynamically driven motion inside the plume to purely ballistic motion in lunar far field conditions. In the earlier project phase of this innovation, the capabilities of the UFS for mixed continuum and rarefied flow situations were validated and demonstrated for lunar lander rocket plume flow impingement under lunar vacuum conditions. Applications and improvements to the granular flow simulation tools contributed by the University of Florida were tested against Earth environment experimental results. Requirements for developing, validating, and demonstrating this solution environment were clearly identified, and an effective second phase execution plan was devised. In this phase, the physics models were refined and fully integrated into a production-oriented simulation tool set. Three-dimensional simulations of Apollo Lunar Excursion Module (LEM) and Altair landers (including full-scale lander geometry) established the practical applicability of the UFS simulation approach and its advanced performance level for large-scale realistic problems.

Liever, Peter↗

Direct numerical simulations of three-component Rayleigh–Taylor mixing and an improved model for multicomponent reacting mixtures

We present direct numerical simulations of a three-layer Rayleigh–Taylor instability (RTI) problem with a configuration based on the experiments of Suchandra & Ranjan ( J. Fluid Mech. , vol. 974, 2023, A35) and Jacobs & Dalziel ( J. Fluid Mech. , vol. 542, 2005, pp. 251–279). The problem consists of a layer of light fluid between two layers of heavy fluid with an Atwood number of 0.3. These simulations are first validated through comparison with available experimental data. The validated simulations are then utilized to analyse statistics in this three-component flow. First, length scales are examined utilizing spectra and two-point spatial correlations of velocity and species concentration fluctuations. Next, joint probability density functions (p.d.f.s) of species concentration are compared against several model p.d.f.s representing generalizations of the bivariate beta distribution. Notably, the joint p.d.f.s do not appear to be accurately described by a Dirichlet distribution, indicating the marginal distributions do not conform to a beta distribution. Finally, similarity of the present configuration to three-component mixing found in inertial confinement fusion (ICF) applications is exploited to develop and validate an improved model for the impact of multicomponent mixing on thermonuclear (TN) reaction rates. A single time instant from the present simulations is chosen for a TN burn calculation under the hypothetical assumption of ICF materials and temperatures. Total TN output from this second calculation is then compared against the prediction of the improved model. The new model is found to accurately predict TN reaction rates in both premixed and non-premixed configurations.

42 ENGINEERING↗