Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale 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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

AMORE-Isoprene v1.0: A New Reduced Mechanism for Gas-Phase Isoprene Oxidation

Gas-phase oxidation of isoprene by ozone (O 3 ) and the hydroxyl (OH) and nitrate (NO 3 ) radicals significantly impacts tropospheric oxidant levels and secondary organic aerosol formation. The most comprehensive and up-to-date chemical mechanism for isoprene oxidation consists of several hundred species and over 800 reactions. Therefore, the computational expense of including the entire mechanism in large-scale atmospheric chemical transport models is usually prohibitive, and most models employ reduced isoprene mechanisms ranging in size from ∼ 10 to ∼ 200 species. We have developed a new reduced isoprene oxidation mechanism using a directed-graph path-based automated model reduction approach, with minimal manual adjustment of the output mechanism. The approach takes as inputs a full isoprene oxidation mechanism, the environmental parameter space, and a list of priority species which are protected from elimination during the reduction process. Our reduced mechanism, AMORE-Isoprene (where AMORE stands for Automated Model Reduction), consists of 12 species which are unique to the isoprene mechanism as well as 22 reactions. We demonstrate its performance in a box model in comparison with experimental data from the literature and other current isoprene oxidation mechanisms. AMORE-Isoprene's performance with respect to predicting the time evolution of isoprene oxidation products, including isoprene epoxydiols (IEPOX) and formaldehyde, is favorable compared with other similarly sized mechanisms. When AMORE-Isoprene is included in the Community Regional Atmospheric Chemistry Multiphase Mechanism 1.0 (CRACMM1AMORE) in the Community Multiscale Air Quality Model (CMAQ, v5.3.3), O 3 and formaldehyde agreement with Environmental Protection Agency (EPA) Air Quality System observations is improved. O 3 bias is reduced by 3.4 ppb under daytime conditions for O 3 concentrations over 50 ppb. Formaldehyde bias is reduced by 0.26 ppb on average for all formaldehyde measurements compared with the base CRACMM1. There was no significant change in computation time between CRACMM1AMORE and the base CRACMM. AMORE-Isoprene shows a 35 % improvement in agreement between simulated IEPOX concentrations and chamber data over the base CRACMM1 mechanism when compared in the Framework for 0-D Atmospheric Modeling (F0AM) box model framework. This work demonstrates a new highly reduced isoprene mechanism and shows the potential value of automated model reduction for complex reaction systems.

Automated Model Reduction↗
NASA NTRS2014

On the Evaluation of Global Sea-Salt Aerosol Models at Coastal/Orographic Sites

Sea-salt aerosol global models are typically evaluated against concentration observations at coastal stations that are unaffected by local surf conditions and thus considered representative of open ocean conditions. Despite recent improvements in sea-salt source functions, studies still show significant model errors in specific regions. Using a multiscale model, we investigated the effect of high model resolution (0.1deg x 0.1deg vs. 1deg x 1.4deg) upon sea-salt patterns in four stations from the University of Miami Network: Baring Head, Chatam Island, and Invercargill in New Zealand, and Marion Island in the sub-antarctic Indian Ocean. Normalized biases improved from +63.7% to +3.3% and correlation increased from 0.52 to 0.84. The representation of sea/land interfaces, mesoscale circulations, and precipitation with the higher resolution model played a major role in the simulation of annual concentration trends. Our results recommend caution when comparing or constraining global models using surface concentration observations from coastal stations.

sea-salt aerosol↗
NASA NTRSDate not supplied

Advancement of Entry System Modeling to Support Exploration of Giant Planets

This paper describes NASA’s efforts to advance entry system modeling and simulation capabilities to support future exploration of Giant planets. The Giant planets are key destinations of interest to the planetary science community for their potential to provide insight into the formation and evolution of our Solar System, as well as extrasolar planetary systems. To date, the Galileo atmospheric probe is the only purpose-built entry probe to a Giant planet. Post-flight analysis of Galileo’s performance showed that there was significant recession of the thermal protection system (TPS), well beyond what was anticipated on the flank, and this was due in part to insufficiently accurate capability for estimating the flight environment and TPS response. While Galileo ultimately survived its flight, the example serves to highlight the great challenge of designing successful missions for environments that are poorly understood or where models have not yet been validated. An important means to reduce mission risks is the incorporation of physics-based modeling with well-quantified uncertainties. The emphasis on physics-based modeling – in contrast to empirically-driven models – is motivated by the fact that it is impossible to completely replicate entry environments through ground tests and, therefore, extrapolation to the flight environment is required. Basing analysis in fundamental physics removes the bias of ground test limitations, though one must then be careful to properly characterize model inputs, simplifying assumptions, and the limits wherein the model is valid. NASA’s Entry Systems Modeling (ESM) Project is tasked with investigating such considerations for planetary science missions across the Solar System, and in recent years has begun to do so for Giant planets. The most distinctive features of the Giant planets, from an entry system perspective, are the atmospheres composed primarily of hydrogen and helium. The entry velocities of proposed missions are generally very large and can therefore be expected to result in significant convective and radiative heating generated by the vehicle’s shock layer. Yet thermochemical behavior of the hydrogen-helium system is not well understood under such conditions. The ESM project is leading efforts to develop accurate thermochemical databases based on state-of-the-art measurements in the Electric Arc Shock Tube and detailed computational chemistry. The large heat fluxes anticipated by missions has driven interest in new TPS materials, in particular woven materials, which may be enabling but have never been flown before. Consequently, multiscale models are in development to describe properties and performance of the materials from micro- to system-scale. The goal is to not only provide accurate thermal response but also to inform thermostructural reliability predictions for extreme entries. Additionally, new computational models have been developed to evaluate performance of non-destructive evaluation techniques which are vital to establishing acceptance of systems to be free of manufacturing faults like material cracking, voids, and debonding. Finally, in the area of guidance and control, aerocapture has been shown conceptually to provide a number of mission benefits, including reducing transit time and increasing payload fraction. The ESM project is building a launch-to-landing trajectory simulation capability to enable detailed studies of aerocapture maneuvers in the context of Giant planets missions. The final presentation and paper will describe each of these topics in detail, including discussion of specific gaps and the technical approach to solving them. In addition, the final paper will briefly discuss ongoing coordination between ESM project work and an ESA-funded technology development activity comprised of validation testing in the Oxford T6, IRS PWK and IST ESTHER tunnels, as well as state-to-state modeling of the shock layer to better represent non-Boltzmann energy distributions leading to non-equilibrium radiation.

Entry systems↗
NASA NTRS2016

Nitrogen Dioxide Observations from the Geostationary Trace Gas and Aerosol Sensor Optimization (GeoTaso) Airborne Instrument: Retrieval Algorithm and Measurements During DISCOVER-AQ Texas 2013

The Geostationary Trace gas and Aerosol Sensor Optimization (GeoTASO) airborne instrument is a test bed for upcoming air quality satellite instruments that will measure backscattered ultraviolet, visible and near-infrared light from geostationary orbit. GeoTASO flew on the NASA Falcon aircraft in its first intensive field measurement campaign during the Deriving Information on Surface Conditions from Column and Vertically Resolved Observations Relevant to Air Quality (DISCOVER-AQ) Earth Venture Mission over Houston, Texas, in September 2013. Measurements of backscattered solar radiation between 420 and 465 nm collected on 4 days during the campaign are used to determine slant column amounts of NO2 at 250 m x 250 m spatial resolution with a fitting precision of 2.2 x 10(exp 15) molecules/sq cm. These slant columns are converted to tropospheric NO2 vertical columns using a radiative transfer model and trace gas profiles from the Community Multiscale Air Quality (CMAQ) model. Total column NO2 from GeoTASO is well correlated with ground-based Pandora observations (r = 0.90 on the most polluted and cloud-free day of measurements and r = 0.74 overall), with GeoTASO NO2 slightly higher for the most polluted observations. Surface NO2 mixing ratios inferred from GeoTASO using the CMAQ model show good correlation with NO2 measured in situ at the surface during the campaign (r = 0.85). NO2 slant columns from GeoTASO also agree well with preliminary retrievals from the GEO-CAPE Airborne Simulator (GCAS) which flew on the NASA King Air B200 (r = 0.81, slope = 0.91). Enhanced NO2 is resolvable over areas of traffic NOx emissions and near individual petrochemical facilities.

Aerosol Sensor Optimization↗
NASA NTRS2019

Electron Scattering by Low-frequency Whistler Waves at Earth’s Bow Shock

Electrons are accelerated to nonthermal energies at shocks in space and astrophysical environments. While shock drift acceleration (SDA) has been considered a key process of electron acceleration at Earth’s bow shock, it has also been recognized that SDA needs to be combined with an additional stochastic process to explain the observed power-law energy spectra. Here, we show mildly energetic (∼0.5 keV) electrons are locally scattered (and accelerated while being confined) by magnetosonic-whistler waves within the shock transition layer, especially when the shock angle is large (θ(sub Bn) approximately equal or greater than 70°). When measured by the Magnetospheric Multiscale mission at a high cadence, ∼0.5 keV electron flux increased exponentially in the shock transition layer. However, the flux profile was not entirely smooth and the fluctuation showed temporal/spectral association with large-amplitude (δB/B ~ 0.3), low-frequency (approximately equal or less than 0.1 Ω(sub ce) where Ω(sub ce) is the cyclotron frequency), obliquely propagating (θ(sub kB) ~ 30°–60°, where θ(sub kB) is the angle between the wave vector and background magnetic field) whistler waves, indicating that the particles were interacting with the waves. Particle simulations demonstrate that, although linear cyclotron resonances with ∼0.5 keV electrons are unlikely due to the obliquity and low frequencies of the waves, the electrons are still scattered beyond 90° pitch angle by (1) resonant mirroring (transit-time damping), (2) non-resonant mirroring, and (3) subharmonic cyclotron resonances. Such coupled nonlinear scattering processes are likely to provide the stochasticity needed to explain the power-law formation.

Oka, M.↗
NASA NTRS2013

The Effect of Scale Dependent Discretization on the Progressive Failure of Composite Materials Using Multiscale Analyses

A multiscale modeling methodology, which incorporates a statistical distribution of fiber strengths into coupled micromechanics/ finite element analyses, is applied to unidirectional polymer matrix composites (PMCs) to analyze the effect of mesh discretization both at the micro- and macroscales on the predicted ultimate tensile (UTS) strength and failure behavior. The NASA code FEAMAC and the ABAQUS finite element solver were used to analyze the progressive failure of a PMC tensile specimen that initiates at the repeating unit cell (RUC) level. Three different finite element mesh densities were employed and each coupled with an appropriate RUC. Multiple simulations were performed in order to assess the effect of a statistical distribution of fiber strengths on the bulk composite failure and predicted strength. The coupled effects of both the micro- and macroscale discretizations were found to have a noticeable effect on the predicted UTS and computational efficiency of the simulations.

Ricks, Trenton M.↗
NASA NTRS2016

Multispacecraft Observations and Modeling of the 22/23 June 2015 Geomagnetic Storm

The magnetic storm of 22-23 June 2015 was one of the largest in the current solar cycle. We present in situ observations from the Magnetospheric Multiscale Mission (MMS) and the Van Allen Probes (VAP) in the magnetotail, field-aligned currents from AMPERE (Active Magnetosphere and Planetary Electrodynamics Response), and ionospheric flow data from Defense Meteorological Satellite Program (DMSP). Our real-time space weather alert system sent out a "red alert," correctly predicting Kp indices greater than 8. We show strong outflow of ionospheric oxygen, dipolarizations in the MMS magnetometer data, and dropouts in the particle fluxes seen by the MMS Fast Plasma Instrument suite. At ionospheric altitudes, the AMPERE data show highly variable currents exceeding 20 MA. We present numerical simulations with the Block Adaptive Tree-Solarwind - Roe - Upwind Scheme (BATS-R-US) global magnetohydrodynamic model linked with the Rice Convection Model. The model predicted the magnitude of the dipolarizations, and varying polar cap convection patterns, which were confirmed by DMSP measurements.

Reiff, P. H.↗
NASA NTRS2022

Superclustering with the Atacama Cosmology Telescope and Dark Energy Survey. I. Evidence for Thermal Energy Anisotropy Using Oriented Stacking

The cosmic web contains filamentary structure on a wide range of scales. On the largest scales, superclustering aligns multiple galaxy clusters along intercluster bridges, visible through their thermal Sunyaev–Zel'dovich signal in the cosmic microwave background. We demonstrate a new, flexible method to analyze the hot gas signal from multiscale extended structures. We use a Compton y-map from the Atacama Cosmology Telescope (ACT) stacked on redMaPPer cluster positions from the optical Dark Energy Survey (DES). Cutout images from the y-map are oriented with large-scale structure information from DES galaxy data such that the superclustering signal is aligned before being overlaid. We find evidence of an extended quadrupole moment of the stacked y signal at the 3.5σ level, demonstrating that the large-scale thermal energy surrounding galaxy clusters is anisotropically distributed. We compare our ACT × DES results with the Buzzard simulations, finding broad agreement. Using simulations, we highlight the promise of this novel technique for constraining the evolution of anisotropic, non-Gaussian structure using future combinations of microwave and optical surveys.

M. Lokken↗
NASA NTRS2012

Expected Navigation Flight Performance for the Magnetospheric Multiscale (MMS) Mission

The Magnetospheric Multiscale (MMS) mission consists of four formation-flying spacecraft placed in highly eccentric elliptical orbits about the Earth. The primary scientific mission objective is to study magnetic reconnection within the Earth s magnetosphere. The baseline navigation concept is the independent estimation of each spacecraft state using GPS pseudorange measurements (referenced to an onboard Ultra Stable Oscillator) and accelerometer measurements during maneuvers. State estimation for the MMS spacecraft is performed onboard each vehicle using the Goddard Enhanced Onboard Navigation System, which is embedded in the Navigator GPS receiver. This paper describes the latest efforts to characterize expected navigation flight performance using upgraded simulation models derived from recent analyses.

Olson, Corwin↗
NASA NTRS2021

Mineral Dust Cycle in the Multiscale Online Nonhydrostatic AtmospheRe CHemistry Model (Monarch) Version 2.0

We present the dust module in the Multiscale Online Non-hydrostatic AtmospheRe CHemistry model (MONARCH) version 2.0, a chemical weather prediction system that can be used for regional and global modeling at a range of resolutions. The representations of dust processes in MONARCH were upgraded with a focus on dust emission (emission parameterizations, entrainment thresholds, considerations of soil moisture and surface cover), lower boundary conditions (roughness, potential dust sources), and dust–radiation interactions. MONARCH now allows modeling of global and regional mineral dust cycles using fundamentally different paradigms, ranging from strongly simplified to physics-based parameterizations. We present a detailed description of these updates along with four global benchmark simulations, which use conceptually different dust emission parameterizations, and we evaluate the simulations against observations of dust optical depth. We determine key dust parameters, such as global annual emission/deposition flux, dust loading, dust optical depth, mass-extinction efficiency, single-scattering albedo, and direct radiative effects. For dust-particle diameters up to 20 µm, the total annual dust emission and deposition fluxes obtained with our four experiments range between about 3500 and 6000 Tg, which largely depend upon differences in the emitted size distribution. Considering ellipsoidal particle shapes and dust refractive indices that account for size-resolved mineralogy, we estimate the global total (longwave and shortwave) dust direct radiative effect (DRE) at the surface to range between about −0.90 and −0.63 W m−2 and at the top of the atmosphere between −0.20 and −0.28 W m−2. Our evaluation demonstrates that MONARCH is able to reproduce key features of the spatiotemporal variability of the global dust cycle with important and insightful differences between the different configurations.

Dust cycle↗
NASA NTRS2016

Ion-Scale Secondary Flux Ropes Generated by Magnetopause Reconnection as Resolved by MMS

New Magnetospheric Multiscale (MMS) observations of small-scale (approx. 7 ion inertial length radius) flux transfer events (FTEs) at the dayside magnetopause are reported. The 1O km MMS tetrahedron size enables their structure and properties to be calculated using a variety of multispacecraft techniques, allowing them to be identified as flux ropes, whose flux content is small (approx. 22 kWb).The current density, calculated using plasma and magnetic field measurements independently, is found to be filamentary. lntercomparison of the plasma moments with electric and magnetic field measurements reveals structured non-frozen-in ion behavior. The data are further compared with a particle-in-cell simulation. It is concluded that these small-scale flux ropes, which are not seen to be growing, represent a distinct class of FTE which is generated on the magnetopause by secondary reconnection.

Eastwood, J. P.↗
NASA NTRS2018

MMS Extended Mission Design: Evaluation of a Lunar Gravity Assist Option

This paper will describe a study that was carried out on the design of a set of maneuvers that were considered for the later stages of an extended mission of the Magnetospheric Multiscale (MMS) mission. The goal of these maneuvers was to put MMS into a significantly different orbit from those flown heretofore, so allowing science collection in a different region of the magnetosphere. This study was made feasible by the fact that the rate at which fuel is being consumed to maintain small formations on the MMS high-apogee orbit is less than expected pre-flight: the current consumption rate is only about 2 kg/yr/spacecraft. In addition, the spacecraft finished the prime mission with a significant amount of fuel remaining: this was about 1-sigma above the mean when compared with pre-launch Monte Carlo simulations. The resulting situation is similar to that of a libration orbit mission, where station-keeping requires so little fuel that any margin at all will lead to an extensive mission lifetime. In the case of MMS, the spacecraft could, if desired, perform formation flying in the current orbit for several decades. Alternatively, the spacecraft could use a significant fraction of the remaining fuel to perform major orbit modifications, while still leaving enough to conduct formation flying for on the order of a decade. The extended mission maneuvers studied here are further apogee-raises, with the goal of setting up one or more lunar gravity assists. Geometry dictates that a lunar encounter is only achievable when the MMS apogee vector lies approximately in the lunar orbit plane: this limits the possible dates to mid-2021 or early 2027.

Design↗
NASA NTRS2019

Mechanism of Reconnection on Kinetic Scales Based on Magnetospheric Multiscale Mission Observations

We examine the role that ions and electrons play in reconnection using observations from the Magnetospheric Multiscale (MMS) mission on kinetic ion and electron scales, which are much shorter than magnetohydrodynamic scales. This study reports observations with unprecedented high resolution that MMS provides for magnetic field (7.8 ms) and plasma (30 ms for electrons and 150 ms for ions). We analyze and compare approaches to the magnetopause in 2016 November, to the electron diffusion region in the magnetotail in 2017 July followed by a current sheet crossing in 2018 July. Besides magnetic field reversals, changes in the direction of the flow velocity, and ion and electron heating, MMS observed large fluctuations in the electron flow speeds in the magnetotail. As expected from numerical simulations, we have verified that when the field lines and plasma become decoupled a large reconnecting electric field related to the Hall current (1–10 mV/m) is responsible for fast reconnection in the ion diffusion region. Although inertial accelerating forces remain moderate (1–2 mV/m), the electric fields resulting from the divergence of the full electron pressure tensor provide the main contribution to the generalized Ohm’s law at the neutral sheet (as large as 200 mV/m). In our view, this illustrates that when ions decouple electron physics dominates. The results obtained on kinetic scales may be useful for better understanding the physical mechanisms governing reconnection processes in various magnetized laboratory and space plasmas.

W. M. Macek↗
NASA NTRS2008

New Developments in the Embedded Statistical Coupling Method: Atomistic/Continuum Crack Propagation

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain has been enhanced. The concurrent MD-FEM coupling methodology uses statistical averaging of the deformation of the atomistic MD domain to provide interface displacement boundary conditions to the surrounding continuum FEM region, which, in turn, generates interface reaction forces that are applied as piecewise constant traction boundary conditions to the MD domain. The enhancement is based on the addition of molecular dynamics-based cohesive zone model (CZM) elements near the MD-FEM interface. The CZM elements are a continuum interpretation of the traction-displacement relationships taken from MD simulations using Cohesive Zone Volume Elements (CZVE). The addition of CZM elements to the concurrent MD-FEM analysis provides a consistent set of atomistically-based cohesive properties within the finite element region near the growing crack. Another set of CZVEs are then used to extract revised CZM relationships from the enhanced embedded statistical coupling method (ESCM) simulation of an edge crack under uniaxial loading.

Saether, E.↗
NASA NTRSDate not supplied

Hierarchical Multiscale Process Modeling of A Textile Composite Y-Joint for the Aurora D8 Aircraft

Woven polymer matrix composites (PMCs) are widely used in aerospace applications for their favorable specific properties. However, the behavior of the inter- and intra-tow matrix during manufacturing of woven PMCs is difficult to isolate experimentally due to geometric and processing complexities. Classical closed-form homogenization approaches allow for composite thermos-chemo-mechanical properties to be characterized at any fiber volume fraction. In this work, virtual curing is applied to the inter- and intra-tow matrix materials of a plain weave repeating unit cell (RUC). Thermo-chemo-mechanical property evolution within the tows is calculated using the Concentric Cylinder Model (CCM) and the Rule of Mixtures (ROM). Curing is implemented through user-written subroutines within the commercial finite element method software Abaqus. Virtual curing results are compared with those from a linear-elastic thermal cool down simulation in Abaqus. Results indicate that tow property evolution during curing has a non-negligible effect on the final stress state.

curing↗
NASA NTRSDate not supplied

Porosity Model Parameter Effects on Metal Additive Manufacturing Simulations

This work investigates the effects of model parameters on thermomechanical simulations of a metal additive manufacturing process. The first half of this work explores the sensitivity and effects of the two porosity parameters used in the evolving constitutive model of the simulation. Varying the two parameters of the model show that for the simulations conducted here, the linear coefficient drives the predicted residual stress results while the constant term has little impact. The second half of the work considers the size, order, and the number of element layers in the finite element mesh. The stress distribution plots show that second-order elements along with a dense mesh marginally improve mesh convergence at a significant cost.

Additive Manufacturing↗
NASA NTRSDate not supplied

NavCube3-mini Lunar GNSS Receiver

This paper describes development and testing of NASA Goddard Space Flight Center’s new NavCube3-mini (NC3m) spaceborne, weak-signal GNSS receiver, which targets all Earth orbit regimes with special focus on lunar applications. NC3m derives from the ground-breaking Magnetospheric Multiscale (MMS) mission Navigator GPS receiver. The MMS-Navigator (launched 2015) holds a Guinness World Record for highest altitude GPS fix and is currently in a highly elliptic orbit with a 29 Earth radii apogee, nearly half lunar distance. NC3m has reduced size, weight, and power compared to the MMS-Navigator, making it suitable for smallsat applications, and adds multi-frequency and multi-GNSS capabilities, among other improvements. A NC3m engineering test unit was subjected to and successfully completed a comprehensive Technology Readiness Level 6 (system/subsystem model or prototype demonstration in a relevant environment) testing campaign in the second half of 2022. This involved high-fidelity simulations in LEO, GEO, and lunar regimes and full environmental testing, including vibration, thermal vacuum, and electromagnetic compatibility. An overview of the NC3m receiver, the test setup, and a sample of results is presented. The results include predicted performance at high-altitude and in the lunar regime.

NavCube3↗
NASA NTRS2010

Multiscale Modeling of Solar Coronal Magnetic Reconnection

Magnetic reconnection is widely believed to be the primary process by which the magnetic field releases energy to plasma in the Sun's corona. For example, in the breakout model for the initiation of coronal mass ejections/eruptive flares, reconnection is responsible for the catastrophic destabilizing of magnetic force balance in the corona, leading to explosive energy release. A critical requirement for the reconnection is that it have a "switch-on' nature in that the reconnection stays off until a large store of magnetic free energy has built up, and then it turn on abruptly and stay on until most of this free energy has been released. We discuss the implications of this requirement for reconnection in the context of the breakout model for CMEs/flares. We argue that it imposes stringent constraints on the properties of the flux breaking mechanism, which is expected to operate in the corona on kinetic scales. We present numerical simulations demonstrating how the reconnection and the eruption depend on the effective resistivity, i.e., the effective Lundquist number, and propose a model for incorporating kinetic flux-breaking mechanisms into MHO calculation of CMEs/flares.

Antiochos, Spiro K.↗