Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Numerical and boundary condition effects on the prediction of detonation engine behavior using detailed numerical simulations

High-fidelity numerical simulations of an experimental rotating detonation engine with discrete fuel/air injection were conducted. A series of configurations with different feed-plenum pressures but with constant equivalence ratio were studied. Detailed chemical kinetics for the hydrogen/air system is used. A resolution study for the full rotating detonation engine (RDE) system simulation is also conducted. Two kinds of boundary conditions, a total pressure boundary and a constant mass flow rate boundary, are used to assess the effects of the inlet boundary. As mass flow rate is increased, the total pressure boundary causes more error in the axial pressure distribution while the constant mass flow rate gives a better solution for all cases ran. The simulations confirm experimental findings, and reproduce qualitative as well as some of the quantitative trends. These results demonstrate that a) fuel-air mixing is highly non-uniform within the detonation chamber, leading to variations in local equivalence ratio, b) the fuel and oxidizer injectors experience significant backflow as the detonation wave passes over, but recover at different rates which further augments the inefficiencies in mixing, and c) parasitic combustion in the mixing region makes the detonation wave weak by extending the reaction zone across the wave.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Conservative numerical schemes with optimal dispersive wave relations: Part I. Derivation and analysis

An energy-conserving and an energy-and-enstrophy conserving numerical schemes are derived by approximating the Hamiltonian formulation of the inviscid shallow water flows based on the vorticity-divergence variables. These schemes also conserve the first-order moments such as mass and vorticity, as usual. The conservative properties of the schemes stem from the skew-symmetry and singularities of the Poisson brackets, which are carefully retained in the discrete approximations. Here, the schemes operate on unstructured orthogonal dual meshes, over bounded or unbounded domains, and they are also shown to possess the same optimal dispersive wave relations as those of the Z-grid scheme, which is a consequence of the use of the vorticity and divergence variables.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Bridging the Transition from Mesoscale to Microscale Turbulence in Numerical Weather Prediction Models

With a focus towards developing multiscale capabilities in numerical weather prediction models, the specific problem of the transition from the mesoscale to the microscale is investigated. For that purpose, idealized one-way nested mesoscale to large-eddy simulation (LES) experiments were carried out using the Weather Research and Forecasting model framework. It is demonstrated that switching from one-dimensional turbulent diffusion in the mesoscale model to three-dimensional LES mixing does not necessarily result in an instantaneous development of turbulence in the LES domain. However, very large fetches are needed for the natural transition to turbulence to occur. The computational burden imposed by these long fetches necessitates the development of methods to accelerate the generation of turbulence on a nested LES domain forced by a smooth mesoscale inflow. To that end, four new methods based upon finite amplitude perturbations of the potential temperature field along the LES inflow boundaries are developed, and investigated under convective conditions. Each method accelerated the development of turbulence within the LES domain, with two of the methods resulting in a rapid generation of production and inertial range energy content associated to microscales that is consistent with non-nested simulations using periodic boundary conditions. The cell perturbation approach, the simplest and most efficient of the best performing methods, was investigated further under neutral and stable conditions. Successful results were obtained in all the regimes, where satisfactory agreement of mean velocity, variances and turbulent fluxes, as well as velocity and temperature spectra, was achieved with reference non-nested simulations. In contrast, the non-perturbed LES solution exhibited important energy deficits associated to a delayed establishment of fully-developed turbulence. The cell perturbation method has negligible computational cost, significantly accelerates the generation of realistic turbulence, and requires minimal parameter tuning, with the necessary information relatable to mean inflow conditions provided by the mesoscale solution.

17 WIND ENERGY↗
DOE OSTI2016

Suppressing the numerical Cherenkov radiation in the Yee numerical scheme

The next generation of laser facilities will routinely produce relativistic particle beams from the interaction of intense laser pulses with solids and/or gases. Their modeling with Particle-In-Cell (PIC) codes needs dispersion-free Maxwell solvers in order to properly describe the interaction of electromagnetic waves with relativistic particles. A particular attention is devoted to the suppression of the numerical Cherenkov instability, responsible for the noise generation. It occurs when the electromagnetic wave is artificially slowed down because of the finite mesh size, thus allowing for the high energy particles to propagate with super-luminous velocities. In the present paper, we show how a slight increase of the light velocity in the Maxwell's equations enables to suppress this instability while keeping a good overall precision of calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Numerical Simulation and Experimental Validation of Joint Performance in Aluminum-Steel Lap Welds Formed by Friction Stir Dovetailing

Friction stir dovetailing (FSD) is a new dissimilar material joining process that simultaneously forms a mechanical interlock and metallurgical bond at the dissimilar material interface. This work presents development of a modeling and simulation approach to predict mechanical performance of FSD thick section aluminum to steel joints. A finite element analysis (FEA) was carried out in order to predict the load-carrying capacity and failure location for aluminum thicknesses between 12.7 mm and 50.8 mm with different numbers of dovetails. The numerical results and corresponding experimental investigation are in agreement and suggest that the developed methods can be used effectively for design and analysis of FSD joints using standard FEA tools.

Friction stir dovetailing, Friction stir welding, ↗
DOE OSTI2019

Numerical study and test of the APS linac transverse deflecting cavity

The transverse deflecting cavity can be used to transform particle distributions in the 6D phase space, which makes it a promising component in phase space beam diagnostics and beam manipulations. In the Advanced Photon Source (APS) Linac, a LOLA-type traveling wave deflecting cavity was installed for the diagnostics of beam characteristics, such as the bunch length, bunch temporal profile, time-dependent energy spread and slice (time-correlated) transverse emittance. In this paper, this deflecting cavity is modeled and analyzed with the numerical method. The effects of the center cell and coupler cell dimensions on the performance of the whole structure are studied, which shows the coupler cell radius has a dominant effect over the coupler cell length and slot width on the global reflection coefficient and field flatness. Important RF parameters, such as S-11, field flatness and phase advance etc., of the LOLA-II cavity in tuning are calculated and discussed. After optimization, the field flatness of the cavity is 2.1%, the phase advance is 119.66 degrees degrees with a standard deviation smaller than 0.5 degrees and the bandwidth is 15.738 MHz when VSWR < 1.1 (i.e. the S-11, < -26 dB). Bead-pull measurements and RF conditioning have been also performed and discussed. Finally, the initial cavity commissioning with electron beam of 375 MeV is reported, and compared with simulation.

coupler↗
DOE OSTI2017

Direct numerical simulations of a high Karlovitz number laboratory premixed jet flame – an analysis of flame stretch and flame thickening [Direct numerical simulations of a high Ka laboratory premixed jet flame - an analysis of flame stretch and flame thickening]

This article reports an analysis of the first detailed chemistry direct numerical simulation (DNS) of a high Karlovitz number laboratory premixed flame. The DNS results are first compared with those from laser-based diagnostics with good agreement. The subsequent analysis focuses on a detailed investigation of the flame area, its local thickness and their rates of change in isosurface following reference frames, quantities that are intimately connected. The net flame stretch is demonstrated to be a small residual of large competing terms: the positive tangential strain term and the negative curvature stretch term. The latter is found to be driven by flame speed–curvature correlations and dominated in net by low probability highly curved regions. Flame thickening is demonstrated to be substantial on average, while local regions of flame thinning are also observed. The rate of change of the flame thickness (as measured by the scalar gradient magnitude) is demonstrated, analogously to flame stretch, to be a competition between straining tending to increase gradients and flame speed variations in the normal direction tending to decrease them. The flame stretch and flame thickness analyses are connected by the observation that high positive tangential strain rate regions generally correspond with low curvature regions; these regions tend to be positively stretched in net and are relatively thinner compared with other regions. Finally, high curvature magnitude regions (both positive and negative) generally correspond with lower tangential strain; these regions are in net negatively stretched and thickened substantially.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Simple and General Unitarity Conserving Numerical Real-Time Propagators of the Time-Dependent Schrödinger Equation Based on Magnus Expansion

Magnus expansion (ME) provides a general way to expand the real-time propagator of a time-dependent Hamiltonian within the exponential such that the unitarity is satisfied at any order. Here, we use this property and explicit integration of Lagrange interpolation formulas for the time-dependent Hamiltonian within each time interval and derive approximations that preserve unitarity for the differential time evolution operators of general time-dependent Hamiltonians. The resulting second-order approximation is the same as using the average of Hamiltonians for two end points of time. We identify three fourth-order approximations involving commutators of Hamiltonians at different times and also derive a sixth-order expression. A test of these approximations along with other available expressions for a two-state time-dependent Hamiltonian with sinusoidal time dependences provides information on the relative performance of these approximations and suggests that the derived expressions can serve as useful numerical tools for time evolution in time-resolved spectroscopy, quantum control, quantum sensing, real-time ab initio quantum dynamics, and open system quantum dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Numerical Simulation of Vibrational Sum Frequency Generation Intensity for Non-Centrosymmetric Domains Interspersed in an Amorphous Matrix: A Case Study for Cellulose in Plant Cell Wall

Vibrational sum frequency generation (SFG) spectroscopy can specifically probe molecular species non-centrosymmetrically arranged in a centrosymmetric or isotropic medium. This capability has been extensively utilized to detect and study molecular species present at the two-dimensional (2D) interface at which the centrosymmetry or isotropy of bulk phases is naturally broken. The same principle has been demonstrated to be very effective for the selective detection of non-centrosymmetric crystalline nanodomains interspersed in three-dimensional (3D) amorphous phases. However, the full spectral interpretation of SFG features has been difficult due to the complexity associated with the theoretical calculation of SFG responses of such 3D systems. This paper describes a numerical method to predict the relative SFG intensities of non-centrosymmetric nanodomains in 3D systems as functions of their size and concentration as well as their assembly patterns, i.e., the distributions of tilt, azimuth, and rotation angles with respect to the lab coordinate. We applied the developed method to predict changes in the CH and OH stretch modes characteristic to crystalline cellulose microfibrils distributed with various orders, which are relevant to plant cell wall structures. As a result, the same algorithm can also be applied to any SFG-active nanodomains interspersed in 3D amorphous matrices.

36 MATERIALS SCIENCE↗
DOE OSTI2012

An investigation of the effect of pore scale flow on average geochemical reaction rates using direct numerical simulation

The scale-dependence of geochemical reaction rates hinders their use in continuum scale models intended for the interpretation and prediction of chemical fate and transport in subsurface environments such as those considered for geologic sequestration of CO 2 . Processes that take place at the pore scale, especially those involving mass transport limitations to reactive surfaces, may contribute to the discrepancy commonly observed between laboratory-determined and continuum-scale or field rates. In this study we investigate the dependence of mineral dissolution rates on the pore structure of the porous media by means of pore scale modeling of flow and multicomponent reactive transport. The pore scale model is composed of high-performance simulation tools and algorithms for incompressible flow and conservative transport combined with a general-purpose multicomponent geochemical reaction code. The model performs direct numerical simulation of reactive transport based on an operator-splitting approach to coupling transport and reactions. The approach is validated with a Poiseuille flow single-pore experiment and verified with an equivalent 1-D continuum-scale model of a capillary tube packed with calcite spheres. Using the case of calcite dissolution as an example, the high-resolution model is used to demonstrate that nonuniformity in the flow field at the pore scale has the effect of decreasing the overall reactivity of the system, even when systems with identical reactive surface area are considered. In conclusion, the effect becomes more pronounced as the heterogeneity of the reactive grain packing increases, particularly where the flow slows sufficiently such that the solution approaches equilibrium locally and the average rate becomes transport-limited.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Convergence of Cloud Droplet Spectral Relative Dispersion During Entrainment‐Mixing Based on Particle‐Resolved Direct Numerical Simulations

Entrainment-mixing processes critically impact cloud microphysical properties, but their effects on the relative dispersion (d) of cloud droplet size distributions (CDSDs) remain elusive. A direct numerical simulation model is initialized with different CDSDs to fill the gap. These results show that d decreases for broad CDSDs and increases for narrow ones, ultimately converging to approximately 0.5 regardless of initial CDSDs during the evaporation-dominated mixing stage. The supersaturation fluctuation and the shape of CDSDs jointly influence the convergence behavior of d. Further sensitivity tests show that the initial microphysical/dynamical/thermodynamical conditions exert negligible effects on the final converged value of d but affect the convergence rate (k). The k generally increases with increasing droplet number concentration and dissipation rate, and increases with decreasing liquid water content, relative humidity of entrained air, and mixing fraction of cloudy air. A conceptual model with two timescales is proposed; k and the timescales are negatively correlated, meaning that slow mixing and/or evaporation process results in slow convergence of d. In conclusion, this finding provides an important reference for improving understanding and parameterization of d during the entrainment-mixing processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

Numerical analysis of three complex three-dimensional shock wave / turbulent boundary layer interaction flows

A set of Reynolds-averaged Navier-Stokes (RANS) computations with turbulence closure provided by the Spalart-Allmaras (SA) model have been carried out for prediction of shock-induced three-dimensional (3D) turbulent separated flows. The experimental data for different aerodynamic test configurations have been used for assessing credibility of the numerical method employed. Particularly, shock wave / turbulent boundary layer interactions (SWTBLI) in the vicinity of an asymmetric sharp double-fin (DF) with different (7° and 11°) deflection angles mounted on a flat plate and two conically sharp cylindrical bodies at varying interbody distances and nose cone angles 60° and 40° mounted over a flat plate at freestream Mach number 4, as well as a transonic fan stage operating in the near-stall regime are considered. The gas dynamic structure and topology of 3D separated flows, surface flow patterns, and pressure distributions as well as body aerodynamic force prediction are analyzed. A transonic fan stage operating in the near-stall regime and a possibility of applying flow control is investigated. High Performance Computing was employed to make high resolution computations of these flows possible, and advanced 3D visualization techniques were employed in order to improve understanding of the separating flow phenomena.

42 ENGINEERING↗
DOE OSTI2024

How numerical simulations helped to achieve breakeven on the NIF

The inertial confinement fusion program relies upon detailed simulations with inertial confinement fusion (ICF) codes to design targets and to interpret the experimental results. These simulations treat as much physics from essential principles as is practical, including laser deposition, cross beam energy transfer, x-ray production and transport, nonlocal thermal equilibrium kinetics, thermal transport, hydrodynamic instabilities, thermonuclear burn, and transport of reaction products. Improvements in radiation hydrodynamic code capabilities and vast increases in computing power have enabled more realistic, accurate 3D simulations that treat all known asymmetry sources. We describe how numerical simulations helped to guide the program, assess the impediments to breakeven, and optimize every aspect of target design. A preshot simulation of the first National Ignition Facility experiment that surpassed breakeven predicted an increased yield that matches the experimental result, within the preshot predicted uncertainty, with a target gain of 1.5. We will cover the key developments in Lawrence Livermore National Laboratory ICF codes that enabled these simulations and give specific examples of how they helped to guide the program.

Marinak, M. M. (ORCID:0009000331127696)↗
DOE OSTI2016

Numerical solution of the Sinanoǧlu equation using a multicentre radial-angular grid

The pair functions that minimize the correlation energy of second-order many-body perturbation (MP2) theory (the Hylleraas functional) are obtained as solutions to the corresponding Sinanoglu equation by expanding them on a six-dimensional, multicenter, radial-angular grid of two electrons. Cusps in the pair functions at the nuclei are described numerically accurately by the multicenter grid. A cusp in each singlet pair function at the coalescence of the two electrons is taken into account analytically by a correlation factor. With a grid of approximately 10,000 points per atom, the MP2 correlation energies for atoms and polyatomic molecules are obtained usually within 0.1 mE h of the complete-basis-set results. The correlation factor, auxiliary basis functions, and a judicious choice of integration algorithms are all necessary to stabilize the grid-based MP2 and underlying Hartree–Fock (HF) calculations. The auxiliary basis set, in particular, largely restores the Hermiticity and diagonal dominance of the Fock matrix as well as furnishes virtual orbitals used in a resolution-of-the-identity approximation to lower the dimension of some integrals. The results of the grid-based HF and MP2 calculations without a correlation factor are found to suffer from large, nonsystematic errors frequently.

Many-body perturbation theory↗
DOE OSTI2016

Comprehensive numerical methodology for direct numerical simulations of compressible Rayleigh-Taylor instability

A comprehensive numerical methodology has been developed that handles the challenges introduced by considering the compressive nature of Rayleigh-Taylor instability (RTI) systems, which include sharp interfacial density gradients on strongly stratified background states, acoustic wave generation and removal at computational boundaries, and stratification-dependent vorticity production. The computational framework is used to simulate two-dimensional single-mode RTI to extreme late-times for a wide range of flow compressibility and variable density effects. The results show that flow compressibility acts to reduce the growth of RTI for low Atwood numbers, as predicted from linear stability analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗