SEARCH · Engineering Papers

Results for “Direct 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,279 records · Page 71

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Comparing Growth Rates of Simulated Moist and Dry Convective Thermals

The spreading rates of convective thermals are linked to their net entrainment, and previous literature has suggested differences in spreading rates between moist and dry thermals. Here, growth rates of idealized numerically simulated axisymmetric dry and moist convective thermals are directly compared. In an environment with dry-neutral stratification, the increase of thermal radius with height dR / dz is a larger by a factor of 1.7 for dry thermals relative to moist thermals. The fractional change in thermal volume ε is also greater for dry thermals within a distance of ~4 radii from the initial thermal height. Values of dR / dz are nearly constant with height for both moist and dry thermals consistent with classical theory based on dimensional analysis. The simulations are also consistent with theory relating impulse, circulation, and spreading rate for dry thermals proposed in previous papers and extended here to moist thermals, predicting they will spread less than dry thermals. Tests adding heating to dry thermals, either spread uniformly across the thermal volume or concentrated in the inner core, indicate that dR / dz and ε are smaller for moist thermals because latent heating is confined mostly to their cores. Additional axisymmetric moist simulations using modified lapse rates and large-eddy simulations support this analysis. Overall, these results indicate that slow spreading rates are a fundamental feature of moist thermals caused by latent heating that alters the spatial distribution of buoyancy within them relative to dry thermals.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1979

Numerical modelling of the formation process of planets from protoplanetary cloud

Evolution of the plane protoplanetary cloud, consisting of a great number of gravitationally interacting and uniting under collision bodies (protoplanets) moving in the central field of a large mass (the Sun or a planet), is considered. It is shown that in the course of protoplanetary cloud evolution the ring zones of matter expansion and compression occur with the subsequent development leading to formation of planets, rotating about their axes mainly directly. The principal numerical results were obtained through digital simulation of planetary accumulation.

Kozlov, N. N.↗

NASA NTRS1988

Crater size estimates for large-body terrestrial impact

Calculating the effects of impacts leading to global catastrophes requires knowledge of the impact process at very large size scales. This information cannot be obtained directly but must be inferred from subscale physical simulations, numerical simulations, and scaling laws. Schmidt and Holsapple presented scaling laws based upon laboratory-scale impact experiments performed on a centrifuge (Schmidt, 1980 and Schmidt and Holsapple, 1980). These experiments were used to develop scaling laws which were among the first to include gravity dependence associated with increasing event size. At that time using the results of experiments in dry sand and in water to provide bounds on crater size, they recognized that more precise bounds on large-body impact crater formation could be obtained with additional centrifuge experiments conducted in other geological media. In that previous work, simple power-law formulae were developed to relate final crater diameter to impactor size and velocity. In addition, Schmidt (1980) and Holsapple and Schmidt (1982) recognized that the energy scaling exponent is not a universal constant but depends upon the target media. Recently, Holsapple and Schmidt (1987) includes results for non-porous materials and provides a basis for estimating crater formation kinematics and final crater size. A revised set of scaling relationships for all crater parameters of interest are presented. These include results for various target media and include the kinematics of formation. Particular attention is given to possible limits brought about by very large impactors.

Schmidt, Robert M.↗

NASA NTRS1989

A review of quasi-coherent structures in a numerically simulated turbulent boundary layer

Preliminary results of a comprehensive study of the structural aspects of a numerically simulated number turbulent boundary layer are presented. A direct Navier-Stokes simulation of a flat-plate, zero pressure gradient boundary layer at Re0 = 670 was used. Most of the known nonrandom, coherent features of turbulent boundary layers are confirmed in the simulation, and several new aspects of their spatial character are reported. The spatial relationships between many of the various structures are described, forming the basis for a more complete kinematical picture of boundary layer physics than has been previously known. In particular, the importance of vortex structures of various forms to the generation of Reynolds shear stress is investigated.

Robinson, S. K.↗

NASA NTRS1992

Interferometric reconstruction of three-dimensional high-speed aerodynamic flows

Holographic interferometry can be a very useful diagnostic tool in high-speed aerodynamic testing. During this summer research period, various possible approaches for accurately reconstructing three-dimensional flows from limited data were examined. The approach based on the combination of the following three techniques appears to be promising: (1) Continuous Local Basis Function Method - this computational tomographic method has a power to accurately reconstruct continuous regions and is appropriate from well-conditioned to moderately limited data; (2) Variable Basis Method - this computational tomographic method provides accuracy near discontinuities, i.e., shock regions, and is appropriate from moderately-limited to severely-limited data; and (3) Complementary Field Method - this is a general iterative reconstructor that can be coupled with any computational tomographic techniques. Mathematically, it can be shown that this method can provide better accuracy than the direct reconstruction as in a conventional approach. Our numerical simulation of experiments demonstrated improved reconstruction results even when these techniques were individually tested.

Cha, Soyoung Stephen↗

NASA NTRS1993

Surrogates for numerical simulations; optimization of eddy-promoter heat exchangers

Although the advent of fast and inexpensive parallel computers has rendered numerous previously intractable calculations feasible, many numerical simulations remain too resource-intensive to be directly inserted in engineering optimization efforts. An attractive alternative to direct insertion considers models for computational systems: the expensive simulation is evoked only to construct and validate a simplified, input-output model; this simplified input-output model then serves as a simulation surrogate in subsequent engineering optimization studies. A simple 'Bayesian-validated' statistical framework for the construction, validation, and purposive application of static computer simulation surrogates is presented. As an example, dissipation-transport optimization of laminar-flow eddy-promoter heat exchangers are considered: parallel spectral element Navier-Stokes calculations serve to construct and validate surrogates for the flowrate and Nusselt number; these surrogates then represent the originating Navier-Stokes equations in the ensuing design process.

Patera, Anthony T.↗

NASA NTRS2001

Textbook Multigrid Efficiency for Computational Fluid Dynamics Simulations

Considerable progress over the past thirty years has been made in the development of large-scale computational fluid dynamics (CFD) solvers for the Euler and Navier-Stokes equations. Computations are used routinely to design the cruise shapes of transport aircraft through complex-geometry simulations involving the solution of 25-100 million equations; in this arena the number of wind-tunnel tests for a new design has been substantially reduced. However, simulations of the entire flight envelope of the vehicle, including maximum lift, buffet onset, flutter, and control effectiveness have not been as successful in eliminating the reliance on wind-tunnel testing. These simulations involve unsteady flows with more separation and stronger shock waves than at cruise. The main reasons limiting further inroads of CFD into the design process are: (1) the reliability of turbulence models; and (2) the time and expense of the numerical simulation. Because of the prohibitive resolution requirements of direct simulations at high Reynolds numbers, transition and turbulence modeling is expected to remain an issue for the near term. The focus of this paper addresses the latter problem by attempting to attain optimal efficiencies in solving the governing equations. Typically current CFD codes based on the use of multigrid acceleration techniques and multistage Runge-Kutta time-stepping schemes are able to converge lift and drag values for cruise configurations within approximately 1000 residual evaluations. An optimally convergent method is defined as having textbook multigrid efficiency (TME), meaning the solutions to the governing system of equations are attained in a computational work which is a small (less than 10) multiple of the operation count in the discretized system of equations (residual equations). In this paper, a distributed relaxation approach to achieving TME for Reynolds-averaged Navier-Stokes (RNAS) equations are discussed along with the foundations that form the basis of this approach. Because the governing equations are a set of coupled nonlinear conservation equations with discontinuities (shocks, slip lines, etc.) and singularities (flow- or grid-induced), the difficulties are many. This paper summarizes recent progress towards the attainment of TME in basic CFD simulations.

Brandt, Achi↗

NASA NTRS2005

GRAVITATIONAL WAVE EXTRACTION FROM AN INSPIRALING CONFIGURATION OF MERGING BLACK HOLES

We present new techniques for evolving binary black hole systems which allow the accurate determination of gravitational waveforms directly from the wave zone region of the numerical simulations. Rather than excising the black hole interiors, our approach follows the "puncture" treatment of black holes, but utilizing a new gauge condition which allows the black holes to move successfully through the computational domain. We apply these techniques to an inspiraling binary, modeling the radiation generated during the final plunge and ringdown. We demonstrate convergence of the waveforms and and good conservation of mass-energy, with just over 3% of the system s mass converted to gravitational radiation.

Baker, John G.↗

DOE OSTI2022

Hybrid simulation of Alfvén wave parametric decay instability in a laboratory relevant plasma

Large-amplitude Alfvén waves are subject to parametric decays, which can have important consequences in space, astrophysical, and fusion plasmas. Although the Alfvén wave parametric decay instability (PDI) was predicted decades ago, observational evidence is limited, stimulating considerable interest in laboratory demonstration of the instability and associated numerical modeling. Here, we report an important step toward direct hybrid simulation of the Alfvén wave PDI in a laboratory plasma, using antenna-like wave injection of a circularly polarized wave and realistic wave-plasma parameters. Considering collisionless damping, we identify the threshold Alfvén wave amplitudes and frequencies required for triggering the instability. These threshold behaviors are corroborated by simple theoretical analysis, considering the balance between PDI growth and ion Landau damping and PDI development in a bounded plasma. Other effects not included in the present model such as finite transverse wave scale and ion–neutral collisions are briefly discussed. These hybrid simulations promise to be a useful tool for investigating laboratory Alfvén wave dynamics and may provide guidance for future laboratory demonstration of the PDI process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1982

Numerical simulation of wall-bounded turbulent shear flows

Developments in three dimensional, time dependent numerical simulation of turbulent flows bounded by a wall are reviewed. Both direct and large eddy simulation techniques are considered within the same computational framework. The computational spatial grid requirements as dictated by the known structure of turbulent boundary layers are presented. The numerical methods currently in use are reviewed and some of the features of these algorithms, including spatial differencing and accuracy, time advancement, and data management are discussed. A selection of the results of the recent calculations of turbulent channel flow, including the effects of system rotation and transpiration on the flow are included.

Moin, P.↗

NASA NTRS1982

Numerical simulation of wall-bounded turbulent shear flows

Moin, P.↗

NASA NTRS1988

Numerical investigations in three-dimensional internal flows

Renewed interest in simulation of high enthalpy flows was brought about by interest in the design and testing of hypersonic propulsion components. This renewed interest has led to refurbishment of high speed ground based wind tunnel facilities. Simulation of the high enthalpy flows occurring in the NASP propulsion system is of particular interest. Studying the combustor presentes some of the more challenging problems in the ground based simulation of the fluid flow. A difficulty in simulating the required conditions at the entrance to the combustor is related to achieveing a sufficiently high static pressure such that combustion would be possible at a high Mach number. An analysis of the flow within the original nozzle was not performed previously using numerical simulations, although it was operated experimentally. A Navier-Stokes code was chosen in the present study to examine analytically the flow within representative selected nozzles for use in combustion experiments. The numerical simulation of these fluid flows is popular and may be able to shed light on the direction that redesign of the shock tube nozzles should take. The numerical codes were modified to account for the axially symmetric flow field, for the inclusion of equilibrium real gas effects and for ideal gases whose thermodynamic properties are frozen at the stagnation conditions. Numerical simulation of flows within the original and two proposed nozzles, including the energy lost to the relatively cold walls, are investigated.

Rose, William C.↗

NASA NTRS1990

Simulations and observations of cloudtop processes

Turbulent entrainment at zero mean shear stratified interfaces has been studied extensively in the laboratory and theoretically for the classical situation in which density is a passive tracer of the mixing and the turbulent motions producing the entrainment are directed toward the interface. It is the purpose of the numerical simulations and data analysis to investigate these processes and, specifically, to focus on the following questions: (1) Can local cooling below cloudtop play an important role in setting up convective circulations within the cloud, and bringing about entrainment; (2) Can Cloudtop Entrainment Instability (CEI) alone lead to runaway entrainment under geophysically realistic conditions; and (3) What are the important mechanisms of entrainment at cloudtop under zero or low mean shear conditions.

Siems, S. T.↗

DOE OSTI2023

Enhanced cutoff energies for direct and rescattered strong-field photoelectron emission of plasmonic nanoparticles

Abstract The efficient generation, accurate detection, and detailed physical tracking of energetic electrons are of applied interest for high harmonics generation, electron-impact spectroscopy, and femtosecond time-resolved scanning tunneling microscopy. We here investigate the generation of photoelectrons (PEs) by exposing plasmonic nanostructures to intense laser pulses in the infrared (IR) spectral regime and analyze the sensitivity of PE spectra to competing elementary interactions for direct and rescattered photoemission pathways. Specifically, we measured and numerically simulated emitted PE momentum distributions from prototypical spherical gold nanoparticles (NPs) with diameters between 5 and 70 nm generated by short laser pulses with peak intensities of 8.0 × 10 12 and 1.2 × 10 13 W/cm 2 , demonstrating the shaping of PE spectra by the Coulomb repulsion between PEs, accumulating residual charges on the NP, and induced plasmonic electric fields. Compared to well-understood rescattering PE cutoff energies for strong-field photoemission from gaseous atomic targets (10× the ponderomotive energy), our measured and simulated PE spectra reveal a dramatic cutoff-energy increase of two orders of magnitude with a significantly higher contribution from direct photoemission. Our findings indicate that direct PEs reach up to 93 % of the rescattered electron cutoff energy, in contrast to 20 % for gaseous atoms, suggesting a novel scheme for the development of compact tunable tabletop electron sources.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1990

Global simulation of the solar wind-magnetosphere interaction - The importance of its numerical validity

This paper examines the global interaction of the solar wind with the geomagnetic field, using a newly developed three-dimensional high-precision MHD simulation code which employs the fourth-order Runge-Kutta-Gill scheme in time and the direct finite-difference method in space, and which allows much less numerical error than the conventional code. The simulation reveals in detail the formation processes of the magnetosphere and the bow shock, as well as the plasma sheet formation process. It is shown that, contrary to the conventional understanding, no Kelvin-Helmholtz instability occurs along the magnetopause in an ideal MHD.

Watanabe, K.↗

NASA NTRS1986

Numerical simulations of rarefied reentry flows

A review of recent calculations obtained with Bird's direct simulation Monte Carlo (DSMC) method is presented for the transitional flowfield encountered at reentry conditions. Consequently, the emphasis is the real-gas effects resulting from a highly energetic nonequilibrium flow. The DSMC calculations for both wide-angle and slender bodies are compared with continuum calculations. The wide-angle body calculations simulate the windward nose of the Shuttle Orbiter and projected aeroassisted orbital transfer vehicles. For the slender, blunted configurations, both two-dimensional and axisymmetric results are presented. Comparisons between the DSMC and continuum calculations show the altitude range where differences in flowfield structure and surface quantities become significant.

Moss, James N.↗

NASA NTRS1974

Laser Doppler velocimeter system simulation for sensing aircraft wake vortices

A hydrodynamic model of aircraft vortex wakes in an irregular wind shear field near the ground is developed and used as a basis for modeling the characteristics of a laser Doppler detection and vortex location system. The trailing vortex sheet and the wind shear are represented by discrete free vortices distributed over a two-dimensional grid. The time dependent hydrodynamic equations are solved by direct numerical integration in the Boussinesq approximation. The ground boundary is simulated by images, and fast Fourier Transform techniques are used to evaluate the vorticity stream function. The atmospheric turbulence was simulated by constructing specific realizations at time equal to zero, assuming that Kolmogoroff's law applies, and that the dissipation rate is constant throughout the flow field. The response of a simulated laser Doppler velocimeter is analyzed by simulating the signal return from the flow field as sensed by a simulation of the optical/electronic system.

Thomson, J. A. L.↗

DOE OSTI2021

Modelling the complex evaporated gas flow and its impact on particle spattering during laser powder bed fusion

The additive manufacturing (AM) of metals is becoming an increasingly important production process with the potential to replace traditional techniques such as casting. Laser Powder Bed Fusion (LPBF) is used in many applications to print metal parts from powder. The metal powder is heated locally with sufficient laser radiation that the liquid melt easily reaches its boiling temperature, which leads to a metallic vapour jet that can entrain both powder bed particles and molten droplets. The small size of laser-matter interaction site makes a detailed experimental analysis of the process challenging. Synchrotron X-ray imaging experiments are one of the few methods which can capture the dynamic melting and solidification processes. Comparing such experiments with computer simulations of the process is an important approach in order to better understand the manufacturing process and to analyse the influence of process parameters on the evaporated gas jet and the subsequent impact on particle ejection, leading to potentially reduced AM component quality. The melting and solidification of the metal powder is simulated using an Eulerian multiphase approach based on a control volume discretization of powder bed and substrate and a volume of liquid separation from melt and gas phase. The gas phase modelled as an ideal gas reaches velocities up to 100 m/s. Lagrangian particle tracking in the simulation demonstrates that the velocity fields calculated by the Eulerian multi-phase approach in combination with a standard drag-force model lead to particle accelerations in good agreement with those measured experimentally. In order to avoid numerical laborious Lagrangian calculations, a direct method to compare an Eulerian multiphase simulation with synchrotron X-ray experiments was introduced and validated. This approach is used to analyse the influence of process parameters including laser power and laser speed on the maximal acceleration of particles from the melt pool area. While the particle acceleration increases linearly with line energy in the conduction mode, a linear decrease of the acceleration with increasing line energy can be found in the transition mode before the accel-eration increases again with line energy in the keyhole mode.

36 MATERIALS SCIENCE↗