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.

SourceResearch connectionFrom yearThrough year

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

At least 1,657 records · Page 92

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Analysis of Fire-Induced Circulations during the FireFlux2 Experiment

Despite recent advances in both coupled fire modeling and measurement techniques to sample the fire environment, the fire–atmosphere coupling mechanisms that lead to fast propagating wildfires remain poorly understood. This knowledge gap adversely affects fire management when wildland fires propagate unexpectedly rapidly and shift direction due to the fire impacts on local wind conditions. In this work, we utilized observational data from the FireFlux2 prescribed burn and numerical simulations performed with a coupled fire–atmosphere model WRF-SFIRE to assess the small-scale impacts of fire on local micrometeorology under moderate wind conditions (10–12 m/s). The FireFlux2 prescribed burn provided a comprehensive observational dataset with in situ meteorological observations as well as IR measurements of fire progression. To directly quantify the effects of fire–atmosphere interactions, two WRF-SFIRE simulations were executed. One simulation was run in a two-way coupled mode in which the heat and moisture fluxes emitted from the fire were injected into the atmosphere, and the other simulation was performed in a one-way coupled mode for which the atmosphere was not affected by the fire. The difference between these two simulations was used to analyze and quantify the fire impacts on the atmospheric circulation at different sections of the fire front. The fire-released heat fluxes resulted in vertical velocities as high as 10.8 m/s at the highest measurement level (20 m above ground level) gradually diminishing with height and dropping to 7.9 m/s at 5.77 m. The fire-induced horizontal winds indicated the strongest fire-induced flow at the lowest measurement levels (as high as 3.3 m/s) gradually decreasing to less than 1 m/s at 20 m above ground level. The analysis of the simulated flow indicates significant differences between the fire-induced circulation at the fire head and on the flanks. The fire-induced circulation was much stronger near the fire head than at the flanks, where the fire did not produce particularly strong cross-fire flow and did not significantly change the lateral fire progression. However, at the head of the fire the fire-induced winds blowing across the front were the strongest and significantly accelerated fire progression. The two-way coupled simulation including the fire-induced winds produced 36.2% faster fire propagation than the one-way coupled run, and more realistically represented the fire progression.

Benik, Jeremy T.↗

DOE OSTI2017

Comparison of Coarse-Grained Approaches in Predicting Polymer Nanocomposite Phase Behavior

Because of the considerable parameter space, efficient theoretical and simulation methods are required to predict the morphology and guide experiments in polymer nanocomposites (PNCs). Unfortunately, theoretical and simulation methods are restricted in their ability to accurately map to experiments based on necessary approximations and numerical limitations. In this study, we provide direct comparisons of two recently developed coarse-grained approaches for modeling polymer nanocomposites (PNCs): polymer nanocomposite field theory (PNC-FT) and dynamic mean-field theory (DMFT). These methods are uniquely suited to efficiently capture mesoscale phase behavior of PNCs in comparison to other theoretical and simulation frameworks. We demonstrate the ability of both methods to capture macrophase separation and describe the thermodynamics of PNCs. We systematically test how the nanoparticle morphology in PNCs is affected by a uniform probability distribution of grafting sites, common in field-based methods, versus random discrete grafting sites on the nanoparticle surface. We also analyze the accuracy of the mean-field approximation in capturing the phase behavior of PNCs. Moreover, the DMFT method introduces the ability to describe nonequilibrium phase behavior while the PNC-FT method is strictly an equilibrium method. With the DMFT method we are able to show the evolution of nonequilibrium states toward their equilibrium state and a qualitative assessment of the dynamics in these systems. These simulations are compared to experiments consisting of polystyrene grafted gold nanorods in a poly(methyl methacrylate) matrix to ensure the model gives results that qualitatively agree with the experiments. This study reveals that nanoparticles in a relatively high matrix molecular weight are trapped in a nonequilibrium state and demonstrates the utility of the DMFT framework in capturing nonequilibrium phase behavior of PNCs. In conclusion, both the PNC-FT and DMFT framework are promising methods to describe the thermodynamic and nonequilibrium phase behavior of PNCs.

36 MATERIALS SCIENCE↗

NASA NTRS1997

Global MHD Magnetospheric Simulation of January 10, 1997 Encounter with Magnetic Cloud

We present a global MHD magnetospheric simulation of the encounter of the magnetic cloud with the terrestrial magnetosphere on January 10-11, 1997. The MHD simulation is driven by IMF and solar wind plasma measurements provided by Wind situated about 100 $R_E$ upstream from Earth. Field-aligned currents produced in the model are mapped down to the ionosphere and are directly compared to auroral images from the Polar UVI. Particular attention will be paid to the interval following the the initial shock wave arrival (around 0100 UT on the 1Oth) but preceding the passage of the magnetic cloud proper (commencing at about 0430 UT) for which there was continuous viewing of the entire auroral oval from apogee by Polar. This turbulent period is characterized by numerous dynamic pressure rises and dips and several northward and southward turnings of the IMF, all of which generate dynamic activity in the simulation that is reflected in the mapped field-aligned current patterns. As discussed by Brittnacher et al. in this session, the auroral morphology imaged by UVI during this period includes: shock wave-induced brightening of the oval followed by a pseudo-onset near midnight, several sun-aligned and curled arcs within the polar cap, an expanding polar cap cleared of arcs, and finally substorm onset at 0337 UT. These features will be directly compared to region 1 and 2 current systems as well as cusp currents in the simulation. Spacecraft magnetopause crossings are also predicted by the global simulation and will be compared to observed crossings, including numerous dayside crossings by Geotail on the 1Oth, and geosynchronous crossings early on January 11.

Elsen, R.↗

DOE OSTI2017

GeneratorSE: A Sizing Tool for Variable-Speed Wind Turbine Generators

This report documents a set of analytical models employed by the optimization algorithms within the GeneratorSE framework. The initial values and boundary conditions employed for the generation of the various designs and initial estimates for basic design dimensions, masses, and efficiency for the four different models of generators are presented and compared with empirical data collected from previous studies and some existing commercial turbines. These models include designs applicable for variable-speed, high-torque application featuring direct-drive synchronous generators and low-torque application featuring induction generators. In all of the four models presented, the main focus of optimization is electromagnetic design with the exception of permanent-magnet and wire-wound synchronous generators, wherein the structural design is also optimized. Thermal design is accommodated in GeneratorSE as a secondary attribute by limiting the winding current densities to acceptable limits. A preliminary validation of electromagnetic design was carried out by comparing the optimized magnetic loading against those predicted by numerical simulation in FEMM4.2, a finite-element software for analyzing electromagnetic and thermal physics problems for electrical machines. For direct-drive synchronous generators, the analytical models for the structural design are validated by static structural analysis in ANSYS.

17 WIND ENERGY↗

NASA NTRS1990

Numerical simulation of rarefied gas flow through a slit

Two different approaches, the finite-difference method coupled with the discrete-ordinate method (FDDO), and the direct-simulation Monte Carlo (DSMC) method, are used in the analysis of the flow of a rarefied gas from one reservoir to another through a two-dimensional slit. The cases considered are for hard vacuum downstream pressure, finite pressure ratios, and isobaric pressure with thermal diffusion, which are not well established in spite of the simplicity of the flow field. In the FDDO analysis, by employing the discrete-ordinate method, the Boltzmann equation simplified by a model collision integral is transformed to a set of partial differential equations which are continuous in physical space but are point functions in molecular velocity space. The set of partial differential equations are solved by means of a finite-difference approximation. In the DSMC analysis, three kinds of collision sampling techniques, the time counter (TC) method, the null collision (NC) method, and the no time counter (NTC) method, are used.

Keith, Theo G., Jr.↗

NASA NTRS2003

Prospect of Using Numerical Dynamo Model for Prediction of Geomagnetic Secular Variation

Modeling of the Earth's core has reached a level of maturity to where the incorporation of observations into the simulations through data assimilation has become feasible. Data assimilation is a method by which observations of a system are combined with a model output (or forecast) to obtain a best guess of the state of the system, called the analysis. The analysis is then used as an initial condition for the next forecast. By doing assimilation, not only we shall be able to predict partially secular variation of the core field, we could also use observations to further our understanding of dynamical states in the Earth's core. One of the first steps in the development of an assimilation system is a comparison between the observations and the model solution. The highly turbulent nature of core dynamics, along with the absence of any regular external forcing and constraint (which occurs in atmospheric dynamics, for example) means that short time comparisons (approx. 1000 years) cannot be made between model and observations. In order to make sensible comparisons, a direct insertion assimilation method has been implemented. In this approach, magnetic field observations at the Earth's surface have been substituted into the numerical model, such that the ratio of the multiple components and the dipole component from observation is adjusted at the core-mantle boundary and extended to the interior of the core, while the total magnetic energy remains unchanged. This adjusted magnetic field is then used as the initial field for a new simulation. In this way, a time tugged simulation is created which can then be compared directly with observations. We present numerical solutions with and without data insertion and discuss their implications for the development of a more rigorous assimilation system.

Kuang, Weijia↗

DOE OSTI2017

Hydrophilic domain structure in polymer exchange membranes: Simulations of NMR spin diffusion experiments to address ability for model discrimination

Here, we detail the development of a flexible simulation program (NMR_DIFFSIM) that solves the nuclear magnetic resonance (NMR) spin diffusion equation for arbitrary polymer architectures. The program was used to explore the proton ( 1 H) NMR spin diffusion behavior predicted for a range of geometrical models describing polymer exchange membranes. These results were also directly compared with the NMR spin diffusion behavior predicted for more complex domain structures obtained from molecular dynamics (MD) simulations. The numerical implementation and capabilities of NMR_DIFFSIM were demonstrated by evaluating the experimental NMR spin diffusion behavior for the hydrophilic domain structure in sulfonated Diels-Alder Poly(Phenylene) (SDAPP) polymer membranes. The impact of morphology variations as a function of sulfonation and hydration level on the resulting NMR spin diffusion behavior were determined. These simulations allowed us to critically address the ability of NMR spin diffusion to discriminate between different structural models, and to highlight the extremely high fidelity experimental data required to accomplish this. A direct comparison of experimental double-quantum-filtered 1 H NMR spin diffusion in SDAPP membranes to the spin diffusion behavior predicted for MD-proposed morphologies revealed excellent agreement, providing experimental support for the MD structures at low to moderate hydration levels.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Photon phase-space dynamics in a plasma wakefield accelerator

Frequency up-shifting of laser light in a beam-driven plasma wakefield has the potential to provide high-intensity sources of short wavelength radiation. Simulations have demonstrated that a laser pulse can undergo large frequency shifts, limited only by the drive beam energy, when the plasma density is tailored to match the accelerating phase of the wake to the group velocity of the pulse. Here, we study the dynamical evolution of photons in the 1D1P phase-space vicinity of the plasma wake-phase matching condition. Numerical calculations using a photon kinetic model are validated by direct comparison with 1D and quasi-3D particle-in-cell simulations. These calculations form the basis of a linear theory of the photon dynamics which reveals several important results, including scalings for the properties of the witness pulse and a self-similar solution for the photon phase-space dynamics. One prediction of the analytic theory is that the pulse can be compressed indefinitely, though the amount of compression would ultimately be limited by practical constraints. These results suggest that photon acceleration can provide a novel source of sub-femtosecond, short wavelength radiation.

XUV generation↗

NASA NTRS2019

Numerical investigations of background pressure effects and channel erosion in the SPT-140 Hall thruster for the Psyche mission

Four SPT-140 Hall thrusters are the main components of the electric propulsion system for the Psyche mission. We employ the 2-D (r-z) code Hall2De to conduct numerical simulations of the SPT-140 with the goal of estimating erosion of the acceleration channel rings for the Psyche mission profile. The simulations are informed by direct measurements of the plasma conditions in the acceleration channel that were obtained using laser-induced fluorescence. We validate our simulation results with additional plasma measurements, wear test erosion rates, and performance data as a function of background pressure. We then make use of a series of numerical simulations to compute the progressive erosion of the channel walls, using as input parameters a profile for the Psyche mission that includes operation at 4.5, 2.5 and 1.0 kW with duration of all segments increased by a 50% margin. We find that the channel walls are not completely eroded at the end of the mission (plus 50% margin). We also observe a progressive decrease in thrust as the channel walls erode that is consistent with the wear test and flight data, with a rate of ~4.5 mN/kh for the first 2,500 h and constant thrust after that. The reasons for these thrust trends are explained.

Lenguito, Giovanni↗

DOE OSTI2020

Experimental and numerical studies of seismic fluid-structure interaction in a base-supported cylindrical vessel

Seismic design and qualification of advanced reactors will rely heavily on the use of verified and validated numerical models capable of capturing the interaction of the vessel, its contained fluid, and the internal equipment: fluid-structure interaction (FSI) analysis. Analytical solutions can be used for preliminary sizing and design of such vessels but their application is limited to simple geometries and boundary conditions, and small amplitude, translational (and rotational) inputs. To validate numerical models for seismic FSI analysis in finite element codes, a comprehensive set of experiments was performed on a liquid-filled cylindrical vessel, using a 6 degree-of-freedom earthquake simulator. Results in terms of sloshing frequency, damping ratio in sloshing modes, and hydrodynamic responses (wave height, hydrodynamic pressure, base shear, and base moment) for multi-directional earthquake simulator inputs are reported and compared with analytical solutions for liquid-filled vessels. The impact of seismic (base) isolation on hydrodynamic responses was studied using earthquake simulator inputs generated using a virtual isolation system. Data from the experiments are used to validate a numerical model of the fluid-filled vessel using the Arbitrary Lagrangian Eulerian (ALE) solver in the commercial finite element program LS-DYNA. Validation studies are presented for multi-directional seismic inputs, including rocking motions. Lagrangian modeling approaches using an elastic material formulation for the fluid are also investigated and their limitations and possible applications are identified. Here, the results are broadly applicable to the seismic response of base supported, liquid-filled vessels.

58 GEOSCIENCES↗

DOE OSTI2017

Influence of capillary end effects on steady-state relative permeability estimates from direct pore-scale simulations

In this study, we investigate and characterize the influence of capillary end effects on steady-state relative permeabilities obtained in pore-scale numerical simulations of two-phase flows. Our study is motivated by the observation that capillary end effects documented in two-phase laboratory-scale experiments can significantly influence permeability estimates. While numerical simulations of two-phase flows in reconstructed pore-spaces are increasingly employed to characterize relative permeabilities, a phenomenon which is akin to capillary end effects can also arise in such analyses due to the constraints applied at the boundaries of the computational domain. We profile the relative strength of these capillary end effects on the calculation of steady-state relative permeabilities obtained within randomly generated porous micro-structures using a finite volume-based two-phase flow solver. Finally, we suggest a procedure to estimate the extent of the regions influenced by these capillary end effects, which in turn allows for the alleviation of bias in the estimation of relative permeabilities.

42 ENGINEERING↗

NASA NTRS1992

Experimental and numerical investigations of low-density nozzle and plume flows of nitrogen

New experimental data are used to show that the direct simulation Monte Carlo (DSMC) method provides an accurate description of a low-density flow in the nozzle and the near-field expansion of a small rocket for two slightly different experimental configurations. These results verify the DSMC method in an expansion flow for the first time. From a number of different gas/surface interaction models, it is found that fully diffuse reflection gives the best agreement with experiment.

Boyd, Iain D.↗

DOE OSTI2020

Fully compressible simulations of waves and core convection in main-sequence stars

Context: Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and requires rather high artificial viscosity and thermal diffusivity for numerical stability. A direct comparison with observations is thus limited. Aims: In this study, we explore the capabilities of our compressible multidimensional Seven-League Hydro (SLH) code to simulate stellar oscillations. Methods: We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic 3 M ⊙ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results: Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with the observed spectra of massive stars. Conclusion: The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows for the coupling of gravity and pressure waves in the outer convective envelopes of evolved stars to be studied in the future.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2001

Formation of Electrostatic Potential Drops in the Auroral Zone

In order to examine the self-consistent formation of large-scale quasi-static parallel electric fields in the auroral zone on a micro/meso scale, a particle in cell simulation has been developed. The code resolves electron Debye length scales so that electron micro-processes are included and a variable grid scheme is used such that the overall length scale of the simulation is of the order of an Earth radii along the magnetic field. The simulation is electrostatic and includes the magnetic mirror force, as well as two types of plasmas, a cold dense ionospheric plasma and a warm tenuous magnetospheric plasma. In order to study the formation of parallel electric fields in the auroral zone, different magnetospheric ion and electron inflow boundary conditions are used to drive the system. It has been found that for conditions in the primary (upward) current region an upward directed quasi-static electric field can form across the system due to magnetic mirroring of the magnetospheric ions and electrons at different altitudes. For conditions in the return (downward) current region it is shown that a quasi-static parallel electric field in the opposite sense of that in the primary current region is formed, i.e., the parallel electric field is directed earthward. The conditions for how these different electric fields can be formed are discussed using satellite observations and numerical simulations.

Schriver, D.↗

DOE OSTI2018

Effects of fuel viscosity on the primary breakup dynamics of a high-speed liquid jet with comparison to X-ray radiography

One of the major concerns in combustion engines is the sensitivity of engine performance to fuel prop- erties. Recent works have shown that even slight differences in fuel properties can cause significant changes in performance and emission of an engine. In order to design the combustion engines with multi-fuel flexi- bilities, the precise assessment of fuel sensitivity on liquid jet atomization process is a prerequisite since the resulting fuel/air mixture is critical to the subsequent combustion process. The present study is focusing on the effect of physical fuel properties, mostly viscosity difference, on the breakup process of the liquid jet injected into still air. Two different jet fuels, CAT-A2 and CAT-C3, are considered here as surrogates for a fossil-based fuel and a bio-derived high-viscosity alternative fuel. The simulations are performed using the volume-of-fluid (VoF) interface tracking method coupled to Lagrangian particle method in order to capture the breakup instabilities of jets and the resulting droplets. The investigations take the actual geometry of the injector into account to resolve the unsteady flow phenomena inside the nozzle that impact the turbu- lence transition and atomization. The simulation results are compared to the experimental measurement using X-ray radiography. Both simulation and X-ray measurements consistently describe the effects of different fuels on the fundamental properties of atomization including the breakup length, transverse liquid volume fraction and the droplet sauter-mean-diameter. The application of a Detailed Numerical Simulation approach complemented by unique X-ray diagnostics is novel and providing new understanding and research directions in engine spray dynamics

Breakup instability↗

DOE OSTI2021

First observation of increased DT yield over prediction due to addition of hydrogen

Here, in a number of reported instances, implosions utilizing fuel mixtures have resulted in anomalously low fusion yields below those predicted by radiation-hydrodynamics simulations. Inter-species ion diffusion has been suggested as a possible cause of the observed yield degradation in fuel mixture implosions. An experimental platform utilizing hydro-equivalent deuterium–tritium (DT), deuterium–tritium–hydrogen (DTH), and deuterium-tritium-helium3 (DT 3 He) capsule implosions was developed to determine whether the inter-species ion diffusion theory may describe the resulting fuel mixture implosion behavior. The implosion experiments were performed at the Omega laser facility. X-ray images and shell areal density diagnostics results show that the hydro-equivalent three capsules (DT, DTH, and DT 3 He) have similar compression behavior. However, nuclear yield deviation was observed from the scaling determined using a fusion yield formula. In the DT 3 He mixture, a reduced yield of a factor of 0.65 ± 0.13 was observed, which is similar to a yield reduction observed in D 3 He mixture by Rygg et al. (i.e., Rygg effect). In contrast, in the DTH mixture, a factor of 1.17 ±0.15 yield increase was observed, which we named the inverse Rygg effect. The yield increase observed in the DTH mixture is consistent with the inter-species ion diffusion theory where lighter H diffuses away from the core and concentrated DT in the core produces higher yield. An inter-species ion diffusion model, the Zimmerman–Paquette–Kagan–Zhdanov model, implemented in a Lagrangian radiation-hydrodynamics fluid code, was also used to analyze the present data, without the need to assume hydrodynamic equivalence of the capsules, but it does not completely explain the DTH or DT 3 He capsules although its effects are in the correct direction. Simulation-based Bayesian inference was used in the latter analysis to quantify the uncertainty in the numerical simulations. The simulation-based analysis resulted in an inferred Rygg-effect yield decrease factor of 0.91 ± 0.02 for the DT 3 He mixture, and an inferred inverse-Rygg yield increase factor of 1.21 ± 0.04 for the DTH mixture, based on simulations ignoring ion diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Integrating micromagnets and hybrid nanowires for topological quantum computing

Majorana zero modes are expected to arise in semiconductor-superconductor hybrid systems, with potential topological quantum computing applications. One limitation of this approach is the need for a relatively high external magnetic field that should also change direction at the nanoscale. This proposal considers devices that incorporate micromagnets to address this challenge. We perform numerical simulations of stray magnetic fields from different micromagnet configurations, which are then used to solve for Majorana wavefunctions. Several devices are proposed, starting with the basic four-magnet design to align magnetic field with the nanowire and scaling up to nanowire T-junctions. The feasibility of the approach is assessed by performing magnetic imaging of prototype patterns.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1979

A two dimensional model of viscous slot flow

The effects of viscosity on flows through slots are investigated. Using Briley's formulation of the Alternating Direction Implicit method together with the Thompson transformation, numerical solution of viscous incompressible flow through a slot is simulated, for Reynolds number 1000 and openness ratio one tenth. The problem is formulated in terms of stream function and vorticity. A wall iteration of Israeli type is employed to iteratively induce the no slip condition on velocity.

Cooke, C. H.↗