Engineering Papers⌕ Search

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.

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

At least 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Digital Lunar Exploration Sites (DLES)

It has been almost 50 years since humans last set foot on the Moon. With NASA’s Artemis program, the United States and its international and commercial partners are embarking on a new endeavor to explore the lunar surface. Before we return, we will have simulated every aspect of these future missions. Many of these models and simulations (M&S) will rely on well-known and commonly-used technologies, some of which trace their origins back to the Apollo program. However, M&S has advanced significantly, as have the underlying computational capabilities. As a result, we are able to model many more aspects of the Artemis vehicles and support systems with significantly improved detail and confidence. Accurately and efficiently modeling the lunar environment will be critical to simulating the Artemis elements and mission activities. This includes characterizing and modeling lunar topography, smaller craters, exposed surface rocks, lunar regolith, surface lighting, and ambient thermal environment. These are all necessary for understanding fundamental behaviors and performance of vehicles and support systems in the lunar environment and are often determining factors in the selection of exploration sites and defining mission profiles. The Astromaterials Research and Exploration Sciences (ARES) team in the Exploration Integration and Science Directorate at NASA’s Johnson Space Center (JSC) and the NASA Exploration Systems Simulations (NExSyS) team in the Simulation and Graphics Branch at JSC are developing and maintaining the Digital Lunar Exploration Sites (DLES) data, documentation, and software packages. This data is being fed directly into a diverse collection of graphics and simulation environments, where it is used to construct the closest known truth for numerous potential Lunar South Pole landing sites. DLES is available to projects across NASA and particularly the Artemis program to support coordinated digital representations of the lunar environment. This paper describes the fundamental need for DLES, the science data sets that are going into DLES, some of the processes used to integrate this data into DLES products, the basic products that constitute DLES, and some examples of DLES in use.

Lunar↗
NASA NTRS2023

Implementing a VLBI Time Delay Model for Earth-orbiting Satellites: Partial derivatives and Verification

This document describes the partial derivatives for commonly estimated parameters including antenna positions, satellite position, and satellite velocity from a near-field VLBI delay model for Earth-orbiting satellites. This model was presented in the Journal of Geodesy by Jaron & Nothnagel (2019). In this context, we present a streamlined version of the near-field VLBI delay. From this simplified model, we deduce a delay rate expression and calculate partial derivatives, maintaining only the terms with significant impact on the computed derivatives’ magnitude. To verify the accuracy of the simplified model and the partial derivatives computed from it, we have created a simple simulation in Matlab of an Earth-orbiting satellite at the altitude of a typical Global Navigation Satellite Systems (GNSS) satellite. From this simulation, we compare the simplified and original delays, and we verify the magnitude and direction of the partial derivatives against the numerically computed derivatives from the original delay model. The partial derivatives and simplified VLBI delay model detailed here are implemented in Fortran in the open-source library VTD.

Geodesy↗
NASA NTRS1981

The potential impact of scatterometry on oceanography - A wave forecasting case

A series of observing system simulation experiments have been performed in order to assess the potential impact of marine surface wind data on numerical weather prediction. In addition to conventional data, the experiments simulated the time-continuous assimilation of remotely sensed marine surface wind or temperature sounding data. The wind data were fabricated directly for model grid points intercepted by a Seasat-1 scatterometer swath and were assimilated into the lowest active level (945 mb) of the model using a localized successive correction method. It is shown that Seasat wind data can greatly improve numerical weather forecasts due to better definition of specific features. The case of the QE II storm is examined.

Cane, M. A.↗
NASA NTRS1990

Surface entropy of liquids via a direct Monte Carlo approach - Application to liquid Si

Two methods are presented for a direct Monte Carlo evaluation of the surface entropy S(s) of a liquid interacting by specified, volume-independent potentials. The first method is based on an application of the approach of Ferrenberg and Swendsen (1988, 1989) to Monte Carlo simulations at two different temperatures; it gives much more reliable results for S(s) in liquid Si than previous calculations based on numerical differentiation. The second method expresses the surface entropy directly as a canonical average at fixed temperature.

Wang, Z. Q.↗
NASA NTRS2006

Opposed Jet Burner Approach for Characterizing Flameholding Potentials of Hydrocarbon Scramjet Fuels

Opposed Jet Burner (OJB) tools have been used extensively by the authors to measure Flame Strength (FS) extinction limits of laminar H2/N2 air and (recently) hydrocarbon (HC) air Counterflow Diffusion Flames (CFDFs) at one atm. This paper details normalization of FSs of N2- diluted H2 and HC systems to account for effects of fuel composition, temperature, pressure, jet diameter, inflow Reynolds number, and inflow velocity profile (plug, contoured nozzle; and parabolic, straight tube). Normalized results exemplify a sensitive accurate means of validating, globally, reduced chemical kinetic models at approx. 1 atm and the relatively low temperatures approximating the loss of non-premixed idealized flameholding, e.g., in scramjet combustors. Laminar FS is defined locally as maximum air input velocity, U(sub air), that sustains combustion of a counter-jet of g-fuel at extinction. It uniquely characterizes a fuel. And global axial strain rate at extinction (U(sub air) normalized by nozzle or tube diameter, D(sub n or (sub t)) can be compared directly with computed extinction limits, determined using either a 1-D Navier Stokes stream-function solution, using detailed transport and finite rate chemistry, or (better yet) a detailed 2-D Navier Stokes numerical simulation. The experimental results define an idealized flameholding reactivity scale that shows wide ranging (50 x) normalized FS s for various vaporized-liquid and gaseous HCs, including, in ascending order: JP-10, methane, JP-7, n-heptane, n-butane, propane, ethane, and ethylene. Results from H2 air produce a unique and exceptionally strong flame that agree within approx. 1% of a recent 2-D numerically simulated FS for a 3 mm tube-OJB. Thus we suggest that experimental FS s and/or FS ratios, for various neat and blended HCs w/ and w/o additives, offer accurate global tests of chemical kinetic models at the Ts and Ps of extinction. In conclusion, we argue the FS approach is more direct and fundamental, for assessing, e.g., idealized scramjet flameholding potentials, than measurements of laminar burning velocity or blowout in a Perfectly Stirred Reactor, because the latter characterize premixed combustion in the absence of aerodynamic strain. And FS directly measures a chemical kinetic characteristic of non-premixed combustion at typical flameholding temperatures. It mimics conditions where gfuels are typically injected into a subsonic flameholding recirculation zone that captures air, where the effects of aerodynamic strain and associated multi-component diffusion become important.

Pellett, Gerald L.↗
DOE OSTI2024

"Hall" transport of liquid crystal solitons in Couette flow

Topology establishes a unifying framework for a diverse range of scientific areas including particle physics, cosmology, and condensed matter physics. One of the most fascinating manifestations of topology in the context of condensed matter is the topological Hall effect, and its relative: the Skyrmion Hall effect. Skyrmions are stable vortex-like spin configurations in certain chiral magnets, and when subject to external electric currents can drift in the transverse direction to the current. These quasi-particles are characterised by a conserved topological charge which in the Skyrmion Hall effect plays the role of electric charges in the ordinary Hall effect. Recently, it has been shown that liquid crystals endowed with chiral properties serve as an ideal testbed for the fundamental investigation of topological solitons, including their two- and three-dimensional realisations. Here, we show experimentally and numerically that three-dimensional solitons aka "torons" exhibit a Hall-like effect when driven by shear flows: the torons are deflected in the direction perpendicular to the shear plane. The experimental results are rationalised in terms of the dynamic Ericksen-Leslie equations, which predict the emergence of the transverse component of the net mass flow, the magnitude of which scales as the 3rd power of the shear rate. The perturbation analysis highlights an interplay of the viscous and chiral elastic torques as the mechanism for the emergence of net transverse currents. Numerical simulations demonstrate, however, that torons are not merely dragged by the flow but move with their own transverse speed, much larger than the average flow velocity in the transverse direction. Our findings may enable responsive microfluidic applications relying on soft topological solitons.

FOS: Physical sciences↗
DOE OSTI2022

Drag reduction and the Vogel exponent of a flexible beam in transient shear flows

Interactions between a flexible beam and a fluid in a channel are of great relevance to biological hairy surfaces, aquatic vegetation, marine life (e.g., fish gills), and many industrial systems alike. While steady state response of a beam to such flows is fairly well-explored, their behavior in the transient regime is not fully understood. Here, a series of numerical simulations are performed to study the laminar Couette flow of an incompressible viscous fluid past an elastic beam in a two-dimensional channel. The flexible beam is perpendicular to the direction of flow, and its base is fixed to the stationary bottom of the channel. We measure the evolution of the Vogel exponent, drag reduction, and reconfiguration number during the transient and steady-state response of the fluid–structure system for different geometrical and physical properties. Our benchmark shows a good agreement between numerical and experimental observations. Our results show that the system's steady-state response at different bulk-fluid velocities can be reproduced by investigating the shear flow response during the transient regime. We define a new variable that characterizes the evolution of the local velocity profile in the proximity of the free end of the beam and use that to characterize the transient-regime response. The analysis yields insight into the competing effects of elasticity of the beam and non-linear flow response.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Validation of SAS4A/SASSYS-1 for Steady-State Single-Phase Natural Circulation

The development of numerical software for engineering applications requires the validation of the code against experimental benchmark datasets. To support the validation of SAS4A/SASSYS-1 in simulating the single-phase natural circulation within nuclear reactor systems, high-precision steady-state experiments are performed to capture single-phase natural circulation phenomena on an existing scaled facility with comprehensive instrumentation. Forced convection tests are designed precisely capturing the facility’s critical thermal-hydraulic parameters required for one-dimensional modeling, and a comprehensive single-phase natural circulation dataset is obtained with well-documented experimental uncertainty and facility description. Analysis and discussion based on the finalized dataset confirm the dataset’s ability in capturing dominant phenomena and address important physical interpretation of parameters. The dataset reported in this project provides a valuable benchmark resource for the validation of system analysis codes under single-phase natural circulation inside nuclear reactors. Validation of SAS4A/SASSYS-1 is then performed against the obtained benchmark dataset to confirm the capability of its physics model in simulating the steady-state single-phase natural circulation. The experimental facility is modeled in the candidate code. Solution verification is performed using Richardson-extrapolation-based estimators to quantify and restrict numerical errors from discretization. Numerical uncertainty originating from finite maximum pseudo-transient time is also quantified and restricted. Input uncertainty provided by the benchmark dataset is forward propagated through the candidate code, directly quantifying the simulation output in a Monte Carlo approach. The composition of the output uncertainty is further determined by the estimators of Sobol’ indices through a variance-based sensitivity analysis. With uncertainty quantified for each individual condition, detailed comparison between the simulation results and experimental data is performed covering the whole dataset, which shows consistent agreement for all important quantities. The validation activity reported in this project demonstrates that SAS4A/SASSYS-1 can predict the primary parameters with satisfactorily accuracy under steady-state single-phase natural circulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Laboratory measurements of the physics of auroral electron acceleration by Alfvén waves

While the aurora has attracted attention for millennia, important questions remain unanswered. Foremost is how auroral electrons are accelerated before colliding with the ionosphere and producing auroral light. Powerful Alfvén waves are often found traveling Earthward above auroras with sufficient energy to generate auroras, but there has been no direct measurement of the processes by which Alfvén waves transfer their energy to auroral electrons. Here, we show laboratory measurements of the resonant transfer of energy from Alfvén waves to electrons under conditions relevant to the auroral zone. Experiments are performed by launching Alfvén waves and simultaneously recording the electron velocity distribution. Numerical simulations and analytical theory support that the measured energy transfer process produces accelerated electrons capable of reaching auroral energies. The experiments, theory, and simulations demonstrate a clear causal relationship between Alfvén waves and accelerated electrons that directly cause auroras.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Synchrotron intensity gradient revealing magnetic fields in galaxy clusters

Magnetic fields and their dynamical interplay with matter in galaxy clusters contribute to the physical properties and evolution of the intracluster medium. However, the current understanding of the origin and properties of cluster magnetic fields is still limited by observational challenges. In this article, we map the magnetic fields at hundreds-kpc scales of five clusters RXC J1314.4-2515, Abell 2345, Abell 3376, MCXC J0352.4-7401, and El Gordo using the synchrotron intensity gradient technique in conjunction with high-resolution radio observations from the Jansky Very Large Array (JVLA) and the Karoo Array Telescope (MeerKAT). We demonstrate that the magnetic field orientation of radio relics derived from synchrotron intensity gradient is in agreement with that obtained with synchrotron polarization. Most importantly, the synchrotron intensity gradient is not limited by Faraday depolarization in the cluster central regions and allows us to map magnetic fields in the radio halos of RXC J1314.4-2515 and El Gordo. We find that magnetic fields in radio halos exhibit a preferential direction along the major merger axis and show turbulent structures at higher angular resolution. The results are consistent with expectations from numerical simulations, which predict turbulent magnetic fields in cluster mergers that are stirred and amplified by matter motions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Scaling patch analysis of turbulent planar plume

Proper scaling in turbulent planar plumes is investigated in this study using a scaling patch approach. Based on the scaled boundary conditions, a proper velocity scale for the mean axial flow is the plume centerline velocity U ref = U ctr , and a proper temperature scale for the temperature excess is Θ ref = T ctr – T ∞ , where T ctr is the plume centerline temperature and T ∞ is the ambient fluid temperature. By seeking an admissible scaling, a key concept in the scaling patch approach, for the mean continuity, mean momentum, and mean energy equations, respectively, the following is found: (1) a proper scale for the mean transverse flow is V ref = ( d δ / d x ) U ctr , where d δ / d x is the growth rate of the plume width. (2) A proper scale for the Reynolds shear stress is R vu , ref = U ctr V ref = ( d δ / d x ) U ctr 2 , a mix of the scales for the mean axial and transverse flows. (3) A proper scale for the turbulent heat flux is R v θ , ref = V ref Θ ctr , a mix of the scales for the mean transverse flow and mean temperature excess. The mean transverse flow thus plays a critical role in the scaling of turbulent planar plumes. Approximate functions are developed for the scaled mean transverse flow, Reynolds shear stress, and turbulent temperature flux, and are found to agree favorably with experimental and numerical simulation data. The integral analysis of the mean momentum equation yields a Richardson number Ri, which remains invariant in the axial direction. The Richardson number is defined as Ri = def g β Θ ctr δ t / ( U ctr V ref ) ≈ 1 / 2 , where g is the gravitational acceleration, β is the thermal expansion coefficient, and δ t is the plume half-width based on the mean temperature profile. This Richardson number arises directly from the scaling patch analysis of the mean momentum equation, including both the streamwise and transverse velocity scales.

42 ENGINEERING↗
DOE OSTI2022

On the physical size of the Milky Way globular cluster NGC 7089 (M2)

ABSTRACT We study the outer regions of the Milky Way globular cluster NGC 7089 based on new Dark Energy Camera observations. The resulting background-cleaned stellar density profile reveals the existence of an extended envelope. We confirm previous results that cluster stars are found out up to ∼1° from the cluster’s centre, which is nearly three times the value of the most robust tidal radii estimations. We also used results from direct N-body simulations in order to compare with the observations. We found a fairly good agreement between the observed and numerically generated stellar density profiles. Because of the existence of gaps and substructures along globular cluster tidal tails, we closely examined the structure of the outer cluster region beyond the Jacobi radius. We extended the analysis to a sample of 35 globular clusters, 20 of them with observed tidal tails. We found that if the stellar density profile follows a power law ∝ r−α, the α slope correlates with the globular cluster present mass, in the sense that, the more massive the globular cluster, the smaller the α value. This trend is not found in globular clusters without observed tidal tails. The origin of such a phenomenon could be related, among other reasons, to the proposed so-called potential escapers or to the formation of globular clusters within dark matter minihaloes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Physics-constrained, low-dimensional models for magnetohydrodynamics: First-principles and data-driven approaches

Plasmas are highly nonlinear and multiscale, motivating a hierarchy of models to understand and describe their behavior. However, there is a scarcity of plasma models of lower fidelity than magnetohydrodynamics (MHD), although these reduced models hold promise for understanding key physical mechanisms, efficient computation, and real-time optimization and control. Galerkin models, obtained by projection of the MHD equations onto a truncated modal basis, and data-driven models, obtained by modern machine learning and system identification, can furnish this gap in the lower levels of the model hierarchy. This work develops a reduced-order modeling framework for compressible plasmas, leveraging decades of progress in projection-based and data-driven modeling of fluids. We begin by formalizing projection-based model reduction for nonlinear MHD systems. To avoid separate modal decompositions for the magnetic, velocity, and pressure fields, we introduce an energy inner product to synthesize all of the fields into a dimensionally consistent, reduced-order basis. Next, we obtain an analytic model by Galerkin projection of the Hall-MHD equations onto these modes. We illustrate how global conservation laws constrain the model parameters, revealing symmetries that can be enforced in data-driven models, directly connecting these models to the underlying physics. We demonstrate the effectiveness of this approach on data from high-fidelity numerical simulations of a three-dimensional spheromak experiment. Finally, this manuscript builds a bridge to the extensive Galerkin literature in fluid mechanics and facilitates future principled development of projection-based and data-driven models for plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Enabling Predictive Scale-Bridging Simulations through Active Learning (Institutional Computing Annual Report (Project w20_alscalebridging)) [Slides]

The goal of this project was to develop, demonstrate, and provide a new capability to achieve greater physical fidelity in large-scale simulations, rather than the usual brute-force increases in the number of mesh elements or particles. This was done by using machine learning (ML) techniques to develop emulators for subscale physics that can be used in coarse-scale continuum simulations, trained on fine-scale molecular dynamics (MD) simulations that are launched on-the-fly via active learning. raditional inertial confinement fusion (ICF) simulations rely on numerical diffusion to simulate molecular effects such as non-local transport and mixing without truly accounting for molecular interactions; our approach directly accounts for this physics.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1987

Numerical simulation of the flow field and fuel sprays in an IC engine

A two-dimensional model for axisymmetric piston-cylinder configurations is developed to study the flow field in two-stroke direct-injection Diesel engines under motored conditions. The model accounts for turbulence by a two-equation model for the turbulence kinetic energy and its rate of dissipation. A discrete droplet model is used to simulate the fuel spray, and the effects of the gas phase turbulence on the droplets is considered. It is shown that a fluctuating velocity can be added to the mean droplet velocity every time step if the step is small enough. Good agreement with experimental data is found for a range of ambient pressures in Diesel engine-type microenvironments. The effects of the intake swirl angle in the spray penetration, vaporization, and mixing in a uniflow-scavenged two-stroke Diesel engine are analyzed. It is found that the swirl increases the gas phase turbulence levels and the rates of vaporization.

Nguyen, H. L.↗
NASA NTRS1989

Transonic Symposium: Theory, Application, and Experiment, volume 1, part 2

In order to assess the state of the art in transonic flow disciplines and to glimpse at future directions, NASA-Langley held a Transonic Symposium. Emphasis was placed on steady, three dimensional external, transonic flow and its simulation, both numerically and experimentally. The symposium included technical sessions on wind tunnel and flight experiments; computational fluid dynamic applications; inviscid methods and grid generation; viscous methods and boundary layer stability; and wind tunnel techniques and wall interference. This, being volume 1, is unclassified.

Foughner, Jerome T., Jr.↗
NASA NTRS1988

Lateral charge transport from heavy-ion tracks in integrated circuit chips

A 256K DRAM has been used to study the lateral transport of charge (electron-hole pairs) induced by direct ionization from heavy-ion tracks in an IC. The qualitative charge transport has been simulated using a two-dimensional numerical code in cylindrical coordinates. The experimental bit-map data clearly show the manifestation of lateral charge transport in the creation of adjacent multiple-bit errors from a single heavy-ion track. The heavy-ion data further demonstrate the occurrence of multiple-bit errors from single ion tracks with sufficient stopping power. The qualitative numerical simulation results suggest that electric-field-funnel-aided (drift) collection accounts for single error generated by an ion passing through a charge-collecting junction, while multiple errors from a single ion track are due to lateral diffusion of ion-generated charge.

Zoutendyk, J. A.↗
NASA NTRS1990

Numerical experiments in supersonic boundary-layer stability

The three-dimensional (3-D) time-dependent compressible Navier-Stokes equations are numerically solved by a Fourier-Chebyshev collocation method to study the stability of supersonic flows over a flat plate. Several direct simulations carried out in this study suggest the existence of a secondary instability that might provide a route to transition. The interaction of the modes involved in the secondary instability is possibly amenable to a Floquet-type analysis. Pertinent differences between this instability and the analogous incompressible K-type secondary instability are pointed out. Some preliminary results of a 2-D direct simulation of the nonlinear evolution of a second mode perturbation wave are also discussed.

Erlebacher, Gordon↗