Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Effects of operator splitting and low Mach-number correction in turbulent mixing transition simulations

In this paper, transition and turbulence decay with the Taylor–Green vortex have been effectively used to demonstrate emulation of high Reynolds-number ( R e ) physical dissipation through numerical convective effects of various non-oscillatory finite-volume algorithms for implicit large eddy simulation (ILES), e.g. using the Godunov-based Eulerian adaptive mesh refinement code xRAGE. The inverse-chevron shock tube experiment simulations have been also used to assess xRAGE based ILES for shock driven turbulent mixing, compared with available simulation and laboratory data. The previous assessments are extended to evaluate new directionally-unsplit high-order algorithms in xRAGE, including a correction to address the well-known issue of excessive numerical diffusion of shock-capturing (e.g., Godunov-type) schemes for low Mach numbers. The unsplit options for hydrodynamics in xRAGE are discussed in detail, followed by fundamental tests with representative shock problems. Basic issues of transition to turbulence and turbulent mixing are discussed, and results of simulations of high- R e turbulent flow and mixing in canonical test cases are reported. Finally, compared to the directional-split cases, and for each grid resolution considered, unsplit results exhibit transition to turbulence with much higher effective R e —and significantly more so with the low Mach number correction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Mesoscale shock structure in particulate composites

Multiscale experiments in heterogeneous materials and the knowledge of their physics under shock compression are limited. Here, this study examines the multiscale shock response of particulate composites comprised of soda-lime glass particles in a PMMA matrix using full-field high speed digital image correlation (DIC) for the first time. Normal plate impact experiments, and complementary numerical simulations, are conducted at stresses ranging from 1.1 - 3.1 GPa to elucidate the mesoscale mechanisms responsible for the distinct shock structure observed in particulate composites. The particle velocity from the macroscopic measurement at continuum scale shows a relatively smooth velocity profile, with shock thickness decreasing with an increase in shock stress, and the composite exhibits strain rate scaling as the second power of the shock stress. In contrast, the mesoscopic response was highly heterogeneous, which led to a rough shock front and the formation of a train of weak shocks traveling at different velocities. Additionally, the normal shock was seen to diffuse the momentum in the transverse direction, affecting the shock rise and the rounding-off observed at the continuum scale measurements. The numerical simulations indicate that the reflections at the interfaces, wave scattering, and interference of these reflected waves are the primary mechanisms for the observed rough shock fronts.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Direct numerical study of leading-edge contamination

Instability, turbulence, and relaminarization in the attachment-line region of swept and unswept cylindrical bodies are studied by numerical solution of the full Navier-Stokes equations. The flow is simulated over a strip containing the attachment-line and treated as homogeneous in the spanwise direction; the disturbances decay exponentially upstream. Transpiration through the wall may be prescribed. The new method, which admits completely general disturbance, agrees with published linear-stability results, which were limited to an apparently restrictive form of disturbance. Fully developed turbulent solutions with sweep are generated and compare well with the experiment. The turbulence is subcritical (except for blowing), resulting in large hysteresis loops. By lowering the sweep Reynolds number, or increasing the suction, the turbulent flow is made to relaminarize. The relaminarization Reynolds number is much less sensitive to suction than the linear-stability Reynolds number. Extensive attempts to detect the postulated nonlinear instability of the unswept flow failed, suggesting that this flow is linearly and nonlinearly stable.

Spalart, P. R.↗
NASA NTRS1988

Mesoscale spiral vortex embedded within a Lake Michigan snow squall band - High resolution satellite observations and numerical model simulations

It is known that Great Lakes snow squall convection occurs in a variety of different modes depending on various factors such as air-water temperature contrast, boundary-layer wind shear, and geostrophic wind direction. An exceptional and often neglected source of data for mesoscale cloud studies is the ultrahigh resolution multispectral data produced by Landsat satellites. On October 19, 1972, a clearly defined spiral vortex was noted in a Landsat-1 image near the southern end of Lake Michigan during an exceptionally early cold air outbreak over a still very warm lake. In a numerical simulation using a three-dimensional Eulerian hydrostatic primitive equation mesoscale model with an initially uniform wind field, a definite analog to the observed vortex was generated. This suggests that intense surface heating can be a principal cause in the development of a low-level mesoscale vortex.

Lyons, Walter A.↗
NASA NTRS2006

Atmospheric Teleconnection over Eurasia Induced by Aerosol Radiative Forcing during Boreal Spring

The direct effects of aerosols on global and regional climate during boreal spring are investigated based on numerical simulations with the NASA Global Modeling and Assimilation Office finite-volume general circulation model (fvGCM) with Microphyics of Clouds with the Relaxed Arakawa Schubert Scheme (McRAS), using aerosol forcing functions derived from the Goddard Ozone Chemistry Aerosol Radiation and Transport model (GOCART). The authors find that anomalous atmospheric heat sources induced by absorbing aerosols (dust and black carbon) excite a planetary-scale teleconnection pattern in sea level pressure, temperature, and geopotential height spanning North Africa through Eurasia to the North Pacific. Surface cooling due to direct effects of aerosols is found in the vicinity and downstream of the aerosol source regions, that is, South Asia, East Asia, and northern and western Africa. Significant atmospheric heating is found in regions with large loading of dust (over northern Africa and the Middle East) and black carbon (over Southeast Asia). Paradoxically, the most pronounced feature in aerosol-induced surface temperature is an east west dipole anomaly with strong cooling over the Caspian Sea and warming over central and northeastern Asia, where aerosol concentrations are low. Analyses of circulation anomalies show that the dipole anomaly is a part of an atmospheric teleconnection pattern driven by atmospheric heating anomalies induced by absorbing aerosols in the source regions, but the influence was conveyed globally through barotropic energy dispersion and sustained by feedback processes associated with the regional circulations. The surface temperature signature associated with the aerosol-induced teleconnection bears striking resemblance to the spatial pattern of observed long-term trend in surface temperature over Eurasia. Additionally, the boreal spring wave train pattern is similar to that reported by Fukutomi et al. associated with the boreal summer precipitation seesaw between eastern and western Siberia. The results of this study raise the possibility that global aerosol forcing during boreal spring may play an important role in spawning atmospheric teleconnections that affect regional and global climates.

Kim, Maeng-Ki↗
NASA NTRS1986

Simulation of the turbulent Rayleigh-Benard problem using a spectral/finite difference technique

The three-dimensional, incompressible Navier-Stokes and energy equations with the Bousinesq assumption have been directly simulated at a Rayleigh number of 3.8 x 10 to the 5th power and a Prandtl number of 0.76. In the vertical direction, wall boundaries were used and in the horizontal, periodic boundary conditions were used. A spectral/finite difference numerical method was used to simulate the flow. The flow at these conditions is turbulent and a sufficiently fine mesh was used to capture all relevant flow scales. The results of the simulation are compared to experimental data to justify the conclusion that the small scale motion is adequately resolved.

Eidson, T. M.↗
NASA NTRS1979

Large eddy simulation of turbulent channel flow: ILLIAC 4 calculation

The three-dimensional time dependent equations of motion were numerically integrated for fully-developed turbulent channel flow. A large scale flow field was obtained directly from the solution of these equations, and small scale field motions were simulated through an eddy viscosity model. The calculations were carried out on the ILLIAC 4 computer. The computed flow patterns show that the wall layer consists of coherent structures of low speed and high speed streaks alternating in the spanwise direction. These structures were absent in the regions away from the wall. Hot spots, small localized regions of very large turbulent shear stress, were frequently observed. The profiles of the pressure velocity-gradient correlations show a significant transfer of energy from the normal to the spanwise component of turbulent kinetic energy in the immediate neighborhood of the wall ('the splatting effect').

Kim, J.↗
NASA NTRS1990

Staggered solution procedures for multibody dynamics simulation

A computational procedure is presented for the direct integration of the multibody dynamical (MBD) equations with constraints. Numerical experiments conducted with this MBD procedure indicate that it yields robust solutions in cases where the step size furnishes more than 20 samples for the period of highest apparent response frequency in a given multibody system. The MBD solution procedure is implemented in two separate modules: the generalized coordinated solver, 'CINT', and the constraint Lagrange multiplier solver, 'LINT'.

Park, K. C.↗
NASA NTRS2000

Numerical Study of Rarefied Hypersonic Flow Interacting with a Continuum Jet

An uncoupled Computational Fluid Dynamics-Direct Simulation Monte Carlo (CFD-DSMC) technique is developed and applied to provide solutions for continuum jets interacting with rarefied external flows. The technique is based on a correlation of the appropriate Bird breakdown parameter for a transitional-rarefied condition that defines a surface within which the continuum solution is unaffected by the external flow-jet interaction. The method is applied to two problems to assess and demonstrate its validity; one of a jet interaction in the transitional-rarefied flow regime and the other in the moderately rarefied regime. Results show that the appropriate Bird breakdown surface for uncoupling the continuum and non-continuum solutions is a function of a non-dimensional parameter relating the momentum flux and collisionality between the two interacting flows. The correlation is exploited for the simulation of a jet interaction modeled for an experimental condition in the transitional-rarefied flow regime and the validity of the correlation is demonstrated. The uncoupled technique is also applied to an aerobraking flight condition for the Mars Global Surveyor spacecraft with attitude control system jet interaction. Aerodynamic yawing moment coefficients for cases without and with jet interaction at various angles-of-attack were predicted, and results from the present method compare well with values published previously. The flow field and surface properties are analyzed in some detail to describe the mechanism by which the jet interaction affects the aerodynamics.

Glass, Christopher E.↗
DOE OSTI2026

Prediction of dynamic behavior of a system undergoing an exchange of components using a transmission simulator method

In systems with modular design, it is of interest to perform component exchanges to upgrade components and explore new system configurations. If the dynamics of the components are known, then dynamic substructuring can be used to obtain a prediction of system performance for the new configuration. The transmission simulator method is a substructuring technique that utilizes an additional fixture to apply an elastic boundary condition to a test article such that an appropriate dynamic response can be obtained for a given substructuring application, and the technique has been used widely for performing coupling and decoupling operations. Here, in this work, a modified transmission simulator method is proposed that allows a direct exchange of components to be performed. In a numerical demonstration, the proposed method is shown to accurately predict the dynamics of an assembly that has undergone a component exchange in which the exchanged components have different material properties and different geometric properties. This approach can be applied to dynamic qualification frameworks in which structures undergo modifications or upgrades rather than complete redesigns.

Component exchange↗
NASA NTRS2000

DSMC Simulation of Separated Flows About Flared Bodies at Hypersonic Conditions

This paper describes the results of a numerical study of interacting hypersonic flows at conditions that can be produced in ground-based test facilities. The computations are made with the direct simulation Monte Carlo (DSMC) method of Bird. The focus is on Mach 10 flows about flared axisymmetric configurations, both hollow cylinder flares and double cones. The flow conditions are those for which experiments have been or will be performed in the ONERA R5Ch low-density wind tunnel and the Calspan-University of Buffalo Research Center (CUBRC) Large Energy National Shock (LENS) tunnel. The range of flow conditions, model configurations, and model sizes provides a significant range of shock/shock and shock/boundary layer interactions at low Reynolds number conditions. Results presented will highlight the sensitivity of the calculations to grid resolution, contrast the differences in flow structure for hypersonic cold flows and those of more energetic but still low enthalpy flows, and compare the present results with experimental measurements for surface heating, pressure, and extent of separation.

Moss, James N.↗
DOE OSTI2021

Multi-Area Model-Free State Estimation via Distributed Tensor Decomposition

This paper proposes a model-free method for distribution system state estimation based on tensor completion using canonical polyadic decomposition. In particular, we consider a setting where the network is divided into multiple areas. The measured physical quantities at buses located in the same area are processed by an area controller. A third-order tensor is constructed to collect these measured quantities. The measurements are analyzed locally to recover the full state information of the network. A closed-form iterative algorithm based on the alternating direction method of multipliers is developed to obtain the low-rank factors of the whole network state tensor where information exchange happens only between neighboring areas. To demonstrate the efficacy of the developed algorithm, numerical simulations are carried out using an IEEE test system.

alternating direction method of multipliers↗
DOE OSTI2025

Experimental and computational evaluation of alpha particle production from laser-driven proton–boron nuclear reaction in hole-boring scheme

The majority of studies on laser-driven proton–boron nuclear reaction is based on the measurement of α-particles with solid-state nuclear tracks detector (Cr39). However, Cr39's interpretation is difficult due to the presence of several other accelerated particles which can bias the analysis. Furthermore, in some laser irradiation geometries, cross-checking measurements are almost impossible. In this case, numerical simulations can play a very important role in supporting the experimental analysis. In our work, we exploited different laser irradiation schemes (pitcher–catcher and direct irradiation) during the same experimental campaign, and we performed numerical analysis, allowing to obtain conclusive results on laser-driven proton–boron reactions. A direct comparison of the two laser irradiation schemes, using the same laser parameters is presented.

Huault, M. (ORCID:0000000186999166)↗
NASA NTRS2009

Chemical Source Inversion using Assimilated Constituent Observations in an Idealized Two-dimensional System

We present a source inversion technique for chemical constituents that uses assimilated constituent observations rather than directly using the observations. The method is tested with a simple model problem, which is a two-dimensional Fourier-Galerkin transport model combined with a Kalman filter for data assimilation. Inversion is carried out using a Green's function method and observations are simulated from a true state with added Gaussian noise. The forecast state uses the same spectral spectral model, but differs by an unbiased Gaussian model error, and emissions models with constant errors. The numerical experiments employ both simulated in situ and satellite observation networks. Source inversion was carried out by either direct use of synthetically generated observations with added noise, or by first assimilating the observations and using the analyses to extract observations. We have conducted 20 identical twin experiments for each set of source and observation configurations, and find that in the limiting cases of a very few localized observations, or an extremely large observation network there is little advantage to carrying out assimilation first. However, in intermediate observation densities, there decreases in source inversion error standard deviation using the Kalman filter algorithm followed by Green's function inversion by 50% to 95%.

Tangborn, Andrew↗
DOE OSTI2018

Hopf bifurcation of a delayed chemostat model with general monotone response functions

In this paper, a chemostat model with general monotone response functions and two discrete time delays is proposed to describe the dynamical behavior about predator–prey system. First, by analyzing the characteristic equation associated with the model, we obtain the conditions of the existence and stability of extinction equilibrium and positive equilibrium. Choosing delays as bifurcation parameters, the existence of Hopf bifurcations is investigated in detail. Second, by virtue of the Poincaré normal form method and center manifold theorem, explicit formulas are derived to determine the direction of Hopf bifurcation and the stability of the bifurcating periodic solutions. Finally, some numerical simulations are carried out to illustrate the theoretical results and the biological significance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Ptychographic phase retrieval by proximal algorithms

We derive a set of ptychography phase-retrieval iterative engines based on proximal algorithms originally developed in convex optimization theory, and discuss their connections with existing ones. The use of proximal operator creates a simple frame work that allows us to incorporate the effect of noise from a maximum-likelihood (ML) principle. We focus on three particular algorithms, namely proximal minimization, alternating direction method of multiplier and accelerated proximal gradient (APG). We benchmark their performance with numerical simulations, and discuss their optimal conditions for convergence and accuracy. An experimental dataset is used to demonstrate their effectiveness as well, in which case an array of cubic Au nanoparticles with a size of 50 nm is imaged. We show that with the presence of Poisson noise, a dataset with photon counts up to 10 4 at one detector pixel already requires ML-based methods to achieve a stable convergence. Among the three algorithms derived in this work, APG method is reported first time for its application in ptychographic reconstruction and shows superior performance in terms of both accuracy and convergence rate with a noisy dataset.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Turbulence-level dependence of cosmic ray parallel diffusion

ABSTRACT Understanding the transport of energetic cosmic rays belongs to the most challenging topics in astrophysics. Diffusion due to scattering by electromagnetic fluctuations is a key process in cosmic ray transport. The transition from a ballistic to a diffusive-propagation regime is presented in direct numerical calculations of diffusion coefficients for homogeneous magnetic field lines subject to turbulent perturbations. Simulation results are compared with theoretical derivations of the parallel diffusion coefficient’s dependences on the energy and the fluctuation amplitudes in the limit of weak turbulence. The present study shows that the widely used extrapolation of the energy scaling for the parallel diffusion coefficient to high turbulence levels predicted by quasi-linear theory does not provide a universally accurate description in the resonant-scattering regime. It is highlighted here that the numerically calculated diffusion coefficients can be polluted for low energies due to missing resonant interaction possibilities of the particles with the turbulence. Five reduced-rigidity regimes are established, which are separated by analytical boundaries derived in this work. Consequently, a proper description of cosmic ray propagation can only be achieved by using a turbulence-level-dependent diffusion coefficient and can contribute to solving the Galactic cosmic ray gradient problem.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Propagation of Plasma Bunches through a Transverse Magnetic Barrier

The injection of a plasma bunch into a multipolar trap can be applied to fill the trap with a plasma. The injection of the bunch into a tokamak-like trap can be considered an additional means for controlling the processes of plasma heating and fuel delivery to the central zone of a thermonuclear reactor. In both cases, the bunch is injected normally to the magnetic field of the trap. It has been shown theoretically, experimentally, and by numerical simulation that the depth of plasma bunch penetration into the magnetic field varies in direct proportion to the bunch energy and in inverse proportion to the magnetic pressure and the cross-sectional area of the plasma bunch. The data of this work allow researchers to estimate the values of plasma bunch parameters at which the bunch will be trapped. As a result, the process of plasma bunch trapping has been optimized.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗