SEARCH · Engineering Papers

Results for “Direct Numerical Simulations”

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,801 records · Page 100

All sources NASA NTRS DOE OSTI

NASA NTRS2000

Rainfall Morphology in Semi-Tropical Convergence Zones

Central Florida is the ideal test laboratory for studying convergence zone-induced convection. The region regularly experiences sea breeze fronts and rainfall-induced outflow boundaries. The focus of this study is the common yet poorly-studied convergence zone established by the interaction of the sea breeze front and an outflow boundary. Previous studies have investigated mechanisms primarily affecting storm initiation by such convergence zones. Few have focused on rainfall morphology yet these storms contribute a significant amount precipitation to the annual rainfall budget. Low-level convergence and mid-tropospheric moisture have both been shown to correlate with rainfall amounts in Florida. Using 2D and 3D numerical simulations, the roles of low-level convergence and mid-tropospheric moisture in rainfall evolution are examined. The results indicate that time-averaged, vertical moisture flux (VMF) at the sea breeze front/outflow convergence zone is directly and linearly proportional to initial condensation rates. This proportionality establishes a similar relationship between VMF and initial rainfall. Vertical moisture flux, which encompasses depth and magnitude of convergence, is better correlated to initial rainfall production than surface moisture convergence. This extends early observational studies which linked rainfall in Florida to surface moisture convergence. The amount and distribution of mid-tropospheric moisture determines how rainfall associated with secondary cells develop. Rainfall amount and efficiency varied significantly over an observable range of relative humidities in the 850- 500 mb layer even though rainfall evolution was similar during the initial or "first-cell" period. Rainfall variability was attributed to drier mid-tropospheric environments inhibiting secondary cell development through entrainment effects. Observationally, 850-500 mb moisture structure exhibits wider variability than lower level moisture, which is virtually always present in Florida. A likely consequence of the variability in 850-500 moisture is a stronger statistical correlation to rainfall, which observational studies have noted. The study indicates that vertical moisture flux forcing at convergence zones is critical in determining rainfall in the initial stage of development but plays a decreasing role in rainfall evolution as the system matures. The mid-tropospheric moisture (e.g. environment) plays an increasing role in rainfall evolution as the system matures. This suggests the need to improve measurements of magnitude/depth of convergence and mid-tropospheric moisture distribution. It also highlights the need for better parameterization of entrainment and vertical moisture distribution in larger-scale models.

Shepherd, J. Marshall↗

NASA NTRS2003

High-Resolution Simulation of Hurricane Bonnie (1998): The Organization of Vertical Motion - Part 1

Hurricanes are well known for their strong winds and heavy rainfall, particularly in the intense rainband (eyewall) surrounding the calmer eye of the storm. In some hurricanes, the rainfall is distributed evenly around the eye so that it has a donut shape on radar images. In other cases, the rainfall is concentrated on one side of the eyewall and nearly absent on the other side and is said to be asymmetric. This study examines how the vertical air motions that produce the rainfall are distributed within the eyewall of an asymmetric hurricane and the factors that cause this pattern of rainfall. We use a sophisticated numerical forecast model to simulate Hurricane Bonnie, which occurred in late August of 1998 during a special NASA field experiment designed to study hurricanes. The simulation results suggest that vertical wind shear (a rapid change in wind speed or direction with height) caused the asymmetric rainfall and vertical air motion patterns by tilting the hurricane vortex and favoring upward air motions in the direction of tilt. Although the rainfall in the hurricane eyewall may surround more than half of the eye, the updrafts that produce the rainfall are concentrated in very small-scale, intense updraft cores that occupy only about 10% of the eyewall area. The model simulation suggests that the timing and location of individual updraft cores are controlled by intense, small-scale vortices (regions of rapidly swirling flow) in the eyewall and that the updrafts form when the vortices encounter low-level air moving into the eyewall.

Braun, Scott A.↗

NASA NTRS2004

Over-Luminous Elliptical Galaxies

The first paper from our work has been completed and accepted for publication. Another paper presents a study of the ESO 30601 70 galaxy group, combining Chandra, XMM-Newton, and optical observations. We find that the system is a true fossil galaxy group - a group whose optical light is dominated by a single galaxy. The group X-ray emission is composed of a central, dense, cool core (10 kpc in radius) and an isothermal medium beyond the central 10 kpc. The region between 10 and 50 kpc (the cooling radius) has the same temperature as the gas from 50 to 400 kpc, although the gas cooling time between 10 and 50 kpc (2-6 Gyr) is shorter than the Hubble time. Thus, the ESO 3060170 group does not have a group-sized cooling core. We suggest that the group cooling core may have been heated by a central active galactic nucleus (AGN) outburst in the past and that the small, dense, cool core is the truncated relic of a previous cooling core. The Chandra observations also reveal a variety of X-ray features in the central region, including a finger, an edge-like feature, and a small tail, all aligned along a north-south axis, as are the galaxy light and group galaxy distribution. The proposed AGN outburst may cause gas to slosh around the center and produce these asymmetric features. The observed flat temperature profile to 1/3rvir is not consistent with the predicted temperature profile in recent numerical simulations. We compare the entropy profile of the ESO 3060170 group with those of three other groups and find a flatter relation than that predicted by simulations involving only shock heating, S approximately r approximately 0.85. This is direct evidence of the importance of non-gravitational processes in group centers. We derive the mass profiles within 1/3rvir and find that the ESO 3060170 group is the most massive fossil group known.

Forman, William↗

NASA NTRS2018

Multi-Scale Modeling of Low-Density Carbon-Phenolic Ablators

Protecting a spacecraft during atmospheric entry is one of highest risk factors that needs to be mitigated during design of a space exploration mission. At entry speeds from space, air turns into high-temperature plasma, and spacecraft Thermal Protection Systems (TPS) are needed to protect the vehicle payload. Modern successful material architectures of spacecraft shields use a porous carbon fiber substrate impregnated with phenolic as an ablator material. In the lecture, efforts to build a Predictive Material Modeling framework for porous ablators from micro-scale to macro-scale will be presented. Several numerical methods and techniques will be summarized that use voxelized images to compute geometrical properties of the porous substrate. These computed properties include porosity, specific surface area and tortuosity that are otherwise indirectly measured through experimental techniques. Direct simulation Monte Carlo (DSMC), a particle-based method for approximating the Boltzmann equation, is used to compute the permeability coefficient of the porous substrate based on its digitized representation. The method computes the flow within the microstructure, where the size of the pores may approach the mean-free-path of the flow. Finally, a high-fidelity model implemented in PATO (Porous-material Analysis Toolbox) is discussed, and some examples of ablative material response are presented, including for the first time 3D simulations of the full tiled heat shield for the Mars Science Laboratory (MSL) capsule.

Ablators↗

DOE OSTI2024

Rethinking Coherent Diffraction Instruments for High Brightness Sources

Over the past two decades, exceptional progress has been made providing coherent x-ray beams at both high-brightness synchrotron sources and x-ray free electron lasers (XFEL). The availability of these coherent x-rays has led to a surge in instruments that exploit x-ray coherence for either x-ray photon correlation spectroscopy (XPCS) or coherent diffraction imaging (CDI). A key to these techniques has been to resolve, or at least nearly resolve, the speckles associated with the scattering from structural disorder. Given practical limits on detector pixel size, this requires long sample-to-detector distances. This paper presents numerical simulations of an instrument that can efMiciently provide a long sample-detector distance while maintaining the ability to rapidly set momentum transfer by inserting crystals with different order in the incident beam, thus changing in incident direction instead of moving the detector.

Fuoss, Paul H↗

DOE OSTI2023

Accurate simulation of direct laser acceleration in a laser wakefield accelerator

In a laser wakefield accelerator (LWFA), an intense laser pulse excites a plasma wave that traps and accelerates electrons to relativistic energies. When the pulse overlaps the accelerated electrons, it can enhance the energy gain through direct laser acceleration (DLA) by resonantly driving the betatron oscillations of the electrons in the plasma wave. The traditional particle-in-cell (PIC) algorithm, although often the tool of choice to study DLA, contains inherent errors due to numerical dispersion and the time staggering of the electric and magnetic fields. Furthermore, conventional PIC implementations cannot reliably disentangle the fields of the plasma wave and laser pulse, which obscures interpretation of the dominant acceleration mechanism. Here, a customized field solver that reduces errors from both numerical dispersion and time staggering is used in conjunction with a field decomposition into azimuthal modes to perform PIC simulations of DLA in an LWFA. Comparisons with traditional PIC methods, model equations, and experimental data show improved accuracy with the customized solver and convergence with an order-of-magnitude fewer cells. Furthermore, the azimuthal-mode decomposition reveals that the most energetic electrons receive comparable energy from DLA and LWFA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

NUMERICAL INVESTIGATION OF FUEL PROPERTY EFFECTS ON MIXED-MODE COMBUSTION IN A SPARK-IGNITION ENGINE

In the present study, mixed-mode combustion of an E30 fuel in a direct-injection spark-ignition engine is numerically investigated at a fuel-lean operating condition using multidimensional computational fluid dynamics (CFD). A fuel surrogate matching Research Octane Number (RON) and Motor Octane Number (MON) of E30 is first developed using neural network based non-linear regression model. To enable efficient 3D engine simulations, a 164-species skeletal reaction mechanism incorporating NO, chemistry is reducedfrom a detailed chemical kinetic model. A hybrid approach that incorporates the G-equation model for tracking turbulentflame front, and the multi-zone well-stirred reactor model for predicting auto-ignition in the end gas, is employed to account for turbulent combustion interactions in the engine cylinder. Predicted in-cylinder pressure and heat release rate traces agree well with experimental measurements. The proposed modelling approach also captures moderated cyclic variability. Two different types of combustion cycles, corresponding to purely deflagrative and mixed-mode combustion, are observed. In contrast to the purely deflagrative cycles, mixed-mode combustion cycles feature early flame propagation followed by end-gas auto-ignition, leading to two distinctive peaks in heat release rate traces. The positive correlation between mixed-mode combustion cycles and early flame propagation is well captured by simulations. With the validated numerical setup, effects of NOx chemistry on mixed-mode combustion predictions are investigated. NOx chemistry is found to promote auto-ignition through residual gas recirculation, while the deflagrative flame propagation phase remains largely unaffected. Local sensitivity analysis is then performed to understand effects of physical and chemical properties of the fuel, i.e., heat of evaporation (HoV) and laminar flame speed (S-L). An increased HoV tends to suppress end-gas auto-ignition due to increased vaporization cooling, while the impact of HoV on flame propagation is insignificant. In contrast, an increased SL is found to significantly promote both flame propagation and auto-ignition. The promoting effect of SL on auto-ignition is not a direct chemical effect; it is rather caused by an advancement of the combustion phasing, which increases compression heating of the end gas.

detailed chemical kinetics↗

NASA NTRS2001

Analysis of High Contrast Imaging Techniques for Space Based Direct Planetary Imaging

We report on our ongoing investigations of a number of techniques for direct detection and imaging of Earth-like planets around nearby stellar sources. Herein, we give a quantitative analysis of these techniques and compare and contrast them via computer simulations. The techniques we will be reporting on are Bracewell Interferometry, Nisenson Apodized Square Aperture, and Coronagraphic masking techniques. We parameterize our results with respect to wavelength, aperture size, effects of mirror speckle, both mid- and high-spatial frequency, detector and photon noise as well pointing error. The recent numerous detections of Jupiter and Saturn like planets has driven a resurgence in research of space based high contrast imaging techniques for direct planetary imaging. Work is currently ongoing for concepts for NASA's Terrestrial Planet Finder mission and a number of study teams have been funded. The authors are members of one team.

Lyon, Richard G.↗

NASA NTRS2003

Numerical simulation of the geographical sources of water for Continental Scale Experiments (CSEs) Precipitation

There are several important research questions that the Global Energy and Water Cycle Experiment (GEWEX) is actively pursuing, namely: What is the intensity of the water cycle and how does it change? And what is the sustainability of water resources? Much of the research to address these questions is directed at understanding the atmospheric water cycle. In this paper, we have used a new diagnostic tool, called Water Vapor Tracers (WVTs), to quantify the how much precipitation originated as continental or oceanic evaporation. This shows how long water can remain in the atmosphere and how far it can travel. The model-simulated data are analyzed over regions of interest to the GEWEX community, specifically, their Continental Scale Experiments (CSEs) that are in place in the United States, Europe, Asia, Brazil, Africa and Canada. The paper presents quantitative data on how much each continent and ocean on Earth supplies water for each CSE. Furthermore, the analysis also shows the seasonal variation of the water sources. For example, in the United States, summertime precipitation is dominated by continental (land surface) sources of water, while wintertime precipitation is dominated by the Pacific Ocean sources of water. We also analyze the residence time of water in the atmosphere. The new diagnostic shows a longer residence time for water (9.2 days) than more traditional estimates (7.5 days). We emphasize that the results are based on model simulations and they depend on the model s veracity. However, there are many potential uses for the new diagnostic tool in understanding weather processes and large and small scales.

Bosilovich, Michael G.↗

DOE OSTI2025

Innovative Symmetrical Electrolyte Architecture Enables Ultra-Flat and Thermal Resilient Protonic Ceramic Electrochemical Cells

Protonic ceramic electrochemical cells (PCECs) represent a promising class of solid‐state energy conversion devices capable of high‐efficiency hydrogen production and power generation. However, the practical deployment of planar PCECs is fundamentally constrained by severe structural deformation and mechanical failure during fabrication, stemming from asymmetric shrinkage between the thin electrolyte and the thick NiO‐based support layer. Here, in this work, a functionally integrated, symmetry‐engineered double‐sided electrolyte (DE) design is unveiled, which not only suppresses thermally induced curvature but also unlocks significant gains in electrochemical performance and stability. This architecture intrinsically balances shrinkage dynamics across the cell bilaterally, enabling the fabrication of ultra‐flat 5 × 5 cm 2 cells with sub‐100 µm thickness variation. A numerical solid mechanics simulation is introduced to investigate and interpret this achievement. Beyond structural advantages, the DE configuration enhances the cell operational stability, delivering a low open‐circuit voltage degradation of 9.5 mV/100 h across 80 thermal cycles. This work establishes a compelling paradigm wherein architectural symmetry directly translates to both mechanical fidelity and functional enhancement, offering a promising route toward PCECs scale‐up.

08 - HYDROGEN↗

NASA NTRS2011

A Fully Nonlinear, Dynamically Consistent Numerical Model for Solid-Body Ship Motion. I. Ship Motion with Fixed Heading

In this paper, we describe the details of our numerical model for simulating ship solidbody motion in a given environment. In this model, the fully nonlinear dynamical equations governing the time-varying solid-body ship motion under the forces arising from ship wave interactions are solved with given initial conditions. The net force and moment (torque) on the ship body are directly calculated via integration of the hydrodynamic pressure over the wetted surface and the buoyancy effect from the underwater volume of the actual ship hull with a hybrid finite-difference/finite-element method. Neither empirical nor free parametrization is introduced in this model, i.e. no a priori experimental data are needed for modelling. This model is benchmarked with many experiments of various ship hulls for heave, roll and pitch motion. In addition to the benchmark cases, numerical experiments are also carried out for strongly nonlinear ship motion with a fixed heading. These new cases demonstrate clearly the importance of nonlinearities in ship motion modelling.

Lin, Ray-Quing↗

DOE OSTI2021

Scoping study of detecting high harmonic fast waves in NSTX-U hot core plasma directly using beam emission spectroscopy

In this paper, we used radio frequency (RF) full wave simulation to explore the feasibility of generating a density beat oscillation formed by high harmonic fast waves (HHFWs) in the National Spherical Tokamak Experiment-Upgrade (NSTX-U) and of measuring it numerically with 2D beam emission spectroscopy (BES). This oscillation is of interest as it can be used to determine the injected RF fields. We began by computing HHFW fields in NSTX-U under a variety of experimental conditions and antenna configurations with Petra-M, a 3D RF wavefield solver. These fields were then used to evaluate the amplitude of δne,beat/ne0 in the volume of measurement, which we predict to be between 5 × 10−10 and 10−9, although it is highly dependent on plasma parameters. Trends in magnetic field, core temperature, core density, and antenna phasing are presented, and the parameters most suitable to this measurement are identified. Finally, we developed a synthetic BES diagnostic to compare the experimental BES signal against.

Rutherford, Grant (ORCID:0000000163620037)↗

DOE OSTI2024

In-cylinder spray evolution in a motored central-injection gasoline engine: Imaging and simulating the effects of flash-boiling and intake crossflow

Accurate predictions of fuel spray behavior and mixture formation in simulations of direct-injection spark-ignition (DISI) engines are fundamental to ensure proper description of all subsequent processes including ignition, combustion, and emissions. In this work, the spray evolution in a single-cylinder optical DISI engine was studied experimentally and numerically with the goal of enabling predictive computational fluid dynamics (CFD) modeling of in-cylinder sprays. The authors explored a wide range of operating conditions characterized by several fuel injection temperatures and engine speeds, using a well-characterized nine-component gasoline surrogate known as PACE-20. The effect of flash boiling and intake crossflow on the spray is discussed, with a focus on evaluating the ability of the spray models to capture highly transient spray behavior. In the experiments, the fuel temperature was varied between 20°C and 80°C, allowing for non-flash- to flash-boiling transition to emerge with enhanced flashing intensity at the highest temperatures. Spray collapse resulted in vapor-rich regions, owing to the locally lower inertia of the fluid. Varying the engine speed from 650 to 1950 rpm promoted increasingly more turbulent in-cylinder crossflow which interacted with the spray during the injection event and resulted in enhanced spray dispersion. The CFD model was able to capture the spray morphology transition at different fuel temperatures and engine speeds adequately. Further, it is shown that the spray breakup model could capture the transitional spray behavior induced by flash boiling atomization and intake flow via proper initialization of the spray cone angle and calibration of the spray models’ constants.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2025

A high-order computational framework for particle-resolved simulations of disperse multiphase flows

This work presents a high-order numerical approach for particle-resolved simulations of disperse multiphase flows, where the Navier-Stokes equations for fluid flow are solved using a high-order spectral element method in the Eulerian framework, and the particle phase is directly simulated with a discrete element method. The coupling between particles and fluids is explicitly handled using an adapted direct-forcing immersed boundary method. Unlike the conventional schemes, a high-order barycentric Lagrange interpolation method and a Gaussian projection kernel are used to ensure accurate momentum exchange between local boundary points and surrounding fluid nodes in the framework of high-order fluid solver. Benchmark tests of increasing complexity are conducted to demonstrate the accuracy and efficiency of our method. Here, it is found that our approach exhibits an excellent convergence performance, as the fluid element/grid is refined and the number of boundary points increases. Compared to conventional low-order methods, the proposed high-order framework enables the use of substantially larger fluid elements while maintaining high accuracy in modeling fluid-particle interactions, owing to the enhanced resolution of high-order basis functions. Moreover, since the primary unknowns are stored at element or grid nodes, the high-order approach offers improved efficiency in both CPU memory usage and total computational cost.

42 ENGINEERING↗

DOE OSTI2022

SNICAR-ADv4: a physically based radiative transfer model to represent the spectral albedo of glacier ice

Abstract. Accurate modeling of cryospheric surface albedo is essential for our understanding of climate change as snow and ice surfaces regulate the global radiative budget and sea-level through their albedo and mass balance. Although significant progress has been made using physical principles to represent the dynamic albedo of snow, models of glacier ice albedo tend to be heavily parameterized and not explicitly connected with physical properties that govern albedo, such as the number and size of air bubbles, specific surface area (SSA), presence of abiotic and biotic light absorbing constituents (LACs), and characteristics of any overlying snow. Here, we introduce SNICAR-ADv4, an extension of the multi-layer two-stream delta-Eddington radiative transfer model with the adding–doubling solver that has been previously applied to represent snow and sea-ice spectral albedo. SNICAR-ADv4 treats spectrally resolved Fresnel reflectance and transmittance between overlying snow and higher-density glacier ice, scattering by air bubbles of varying sizes, and numerous types of LACs. SNICAR-ADv4 simulates a wide range of clean snow and ice broadband albedo (BBA), ranging from 0.88 for (30 µm) fine-grain snow to 0.03 for bare and bubble-free ice under direct light. Our results indicate that representing ice with a density of 650 kg m−3 as snow with no refractive Fresnel layer, as done previously, generally overestimates the BBA by an average of 0.058. However, because most naturally occurring ice surfaces are roughened “white ice”, we recommend modeling a thin snow layer over bare ice simulations. We find optimal agreement with measurements by representing cryospheric media with densities less than 650 kg m−3 as snow and larger-density media as bubbly ice with a Fresnel layer. SNICAR-ADv4 also simulates the non-linear albedo impacts from LACs with changing ice SSA, with peak impact per unit mass of LACs near SSAs of 0.1–0.01 m2 kg−1. For bare, bubble-free ice, LACs actually increase the albedo. SNICAR-ADv4 represents smooth transitions between snow, firn, and ice surfaces and accurately reproduces measured spectral albedos of a variety of glacier surfaces. This work paves the way for adapting SNICAR-ADv4 to be used in land ice model components of Earth system models.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1981

Simulation of three-dimensional incompressible flows with a vortex-in-cell method

A new method for the numerical simulation of three-dimensional incompressible flows is described. The vortex-in-cell (VIC) method presented traces the motion of the vortex filaments in the velocity field which these filaments create. The velocity field is not calculated directly by the Biot-Savart law of interaction but by creating a mesh-record of the vorticity field, then integrating a Poisson's equation via the fast Fourier transform to generate a mesh-record of the velocity field. The computed scales of motion are assumed to be essentially inviscid. Viscous of subgrid-scale effects are incorporated into a filtering procedure in wave vector space. Results of tracing a periodic array of single vortex rings are compared with a Green's function calculation. The agreement is very good.

Couet, B.↗

NASA NTRS1993

Rotational effects on convection simulated at different latitudes

We simulate numerically convection inside the solar convection zone under the influence of rotation at different latitudes. The computational domain is a small rectangular box with stress-free upper and lower boundaries, and with periodicity assumed in the lateral directions. We study the transport of angular momentum, which is important for the generation of differential rotation. The sign and the latitudinal dependence of the horizontal Reynolds stress component turn out to be in good agreement with correlation measurements of sunspot proper motions and with predictions from the theory of the Lambda effect. We also investigate the other components of the Reynolds stress as well as the eddy heat flux tensor, both of which are needed in mean field models of differential rotation.

Pulkkinen, Pentti↗

DOE OSTI2023

Examination of a Methane/Diesel RCCI Engine Using Pele: Preprint

Multi-fuel, advanced injection strategies have become increasingly promising as a strategy to mitigate the emissions generated from internal combustion engines. By carefully controlling the combustion phasing in-cylinder, these new multi-pulse, multi-fuel injection strategies are able to burn in the low-temperature combustion regime where both NOx and soot are not readily produced, reducing the need for extensive exhaust gas recirculation systems. In this study, we examine a reactivity-controlled compression ignition (RCCI) strategy that uses an early pre-filled methane-air mixture with low turbulence background as the low-reactivity fuel and a direct injection of four discrete dodecane jets as a surrogate for the high-reactivity diesel fuel. We use the Pele software suite, a highly optimized, exascale-ready, adaptive mesh refinement codebase to perform high-resolution numerical simulations of a scaled down, single cylinder from the RCCI engine. Here, we resolve the ignition kernels down to micrometer scales and present several statistical quantities evaluating the development of the flow and detailing the onset of ignition and subsequent flame development. Particular attention is paid to the conditions surrounding the onset of the first ignition kernels and discussing what led to the development of those conditions.

CFD↗

Refine your search to explore more results.