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.

267 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Investigation of Lunar-Inspired Geopolymer Concrete Formulations Mixed and Cured in Microgravity on the International Space Station (ISS)

The research outlined in this presentation investigates the use of various lunar regolith simulants in geopolymer lunar concrete mixes mixed and cured on the International Space Station (ISS). The motivation for this work is to study the effects of gravity on the microstructure of alkali-activated materials cured with heat, and to develop materials for the construction of long-term infrastructure on the lunar surface with in-situ resource utilization (ISRU). ISRU for construction materials reduces the cost and mass of payloads related to lunar construction. The advantage of geopolymer concrete as opposed to traditional portland cement concrete is that water acts as a medium for the polymerization reaction and leaves the system throughout the process, reducing its demand. Twelve samples of lunar regolith simulant and a solution composed of sodium hydroxide and sodium silicate were sent to the ISS. The three simulants were OPRH2N, OPRL2N, and JSC-1AF, using only particles less than 53 µm in diameter to increase reactivity of the simulant. Simulant to solution ratios were determined by workability while mixing. The simulant and solution were sealed Burst Pouches® along with 2 other sealed bags to prevent material from leaking. Crew member F-14 conducted testing on the ISS by introducing the solution to the simulant in the Burst Pouch®, mixing the sample with a spatula, and then clamping the specimen in the fresh state to prevent flow inside the Burst Pouch®. These specimens were then put in a thermos heated to 80C via sealed drinking water bags to cure for 24 hours with a temperature logger. The cured specimens remained in microgravity for at least 28 days and were returned from the ISS in February 2025. The specimens were then brought to the NASA Marshall Space Flight Center (MSFC) to analyze. Material characterization consisted of conducting Micro-CT tests of entire samples in their sealed apparatus to a resolution of 25µm. 2D image slices were saved in each orthogonal direction of each specimen at a 0.03 mm step size from the 3D model to conduct analytical porosity calculations. Representative samples from each specimen were sampled to perform helium gas pycnometery and were then mounted in resin for SEM imaging, EDS, and nanoindentation. Porosity was analyzed analytically using micromechanics modelling with the assistance of the NASA Multiscale Analysis Tool (NASMAT), as well as the NASA Advanced Supercomputing (NAS) servers (V. Saseendran & N. Yamamoto, 2024). Density was measured using helium gas pycnometery and was then compared to the theoretical density for experimental porosity calculation. Due to the samples’ non-uniform shape being cured in a pouch, traditional compression and tensile strength testing could not be performed. Nanoindentation was conducted at Clarkson University to determine the microhardness and reduced modulus of elasticity. Results from flight samples can be compared to ground samples currently in DLR’s possession to determine the effect on microstructure from being mixed and cured in microgravity. This study gives further insight and understanding of geopolymer lunar concrete and its viability as a lunar construction material with ISRU.

Adam Johnson↗
NASA NTRS2005

The Effects of Thrust Uncertainty and Attitude Knowledge Errors on the MMS Formation Maintenance Maneuver

The Magnetospheric Multiscale (MMS) Mission will use a formation of four spinning spacecraft to study the Earth s magnetosphere. The science objectives of MMS require a near-regular tetrahedron formation to be maintained with side lengths ranging from ten kilometers to several thousand kilometers at orbit apogee. To reduce spacecraft complexity and cost, the current mission concept assumes MMS can achieve its formation goals through open-loop orbit control from the ground, rather than in-flight, closed-loop formation control that has been the subject of recent study. Significant analysis has been performed to provide optimal reference orbit and relative orbit designs. However, the feasibility of achieving these orbits, and maintaining them for an extended period of time in the presence of real world errors and perturbations has not been investigated. In particular, attitude knowledge and control errors, which may have a negligible effect on orbit control for conventional missions with spinning spacecraft, can contribute significant errors to the MMS orbits. In this work, a 6 degree-of-freedom (DOF) simulation is developed and used to analyze the effects of realistic errors on formation maintenance maneuver accuracy. Several realistic considerations including a finite-burn model, attitude perturbations, and thrust uncertainty are studied. The primary objective is to quantify the effects of realistic attitude and orbit control, knowledge, and actuator errors on the formation geometry by observing representative maneuver errors of a single spacecraft.

Tsai, Dean C.↗
NASA NTRS2018

Investigating Extreme Events in the NASA GEOS Multiscale Modeling System

Advances in computing capabilities and scientific development have come together to evolve general circulation models into multi-scale Earth system modeling tools. The Goddard Earth Observing System (GEOS) model is one such example of this evolution. The GEOS model is driven by the finite-volume cubed-sphere (FV3) non-hydrostatic dynamical core. Surrounding FV3 is a scale-aware physics package and data assimilation capability permitting multi-scale application of GEOS for sub-seasonal to seasonal climate prediction, medium range weather prediction, and global mesoscale modeling at convection allowing resolutions. GEOS also includes a comprehensive chemistry package representing a range of capabilities from basic chemistry and interactve aerosols and gaseous species, to carbon emissions and uptake, and complex ozone photochemistry. In this study, we apply the GEOS model to study the fidelity of these processes with increasing horizontal resolution and scale-aware processes in GEOS on extreme events. The GEOS model is run for 40-days beginning in August 2016 at three uniform global resolutions of 13-km (c768), 6-km (c1536) and 3-km (c3072) with 72 vertical levels up to 0.01mb. The model physics use the Grell-Freitas scale-aware convection scheme to dynamically reduce the role of parameterized deep convection as resolved scale processes in the model take over at higher resolutions. We will include high-resolution global emissions and fluxes of aerosols and carbon downscaled from recent satellite observations. We will compare these simulations with reanalyses and observations, focusing on rainfall, clouds and radiative forcing at hourly to monthly timescales. We will closely examine the probability distribution of precipitation intensities and radiative properties of clouds on a daily time scale. In addition, we will focus on extreme events, in particular the diurnal cycle of convection over the US and the frequency and physical nature of organized convection and heavy rain events across the globe.

Putman, William↗
NASA NTRS2008

Calculations of High-Temperature Jet Flow Using Hybrid Reynolds-Average Navier-Stokes Formulations

Two multiscale-type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds-averaged Navier Stokes equations. The first scheme is a hybrid Reynolds-averaged- Navier Stokes/large-eddy-simulation model using the two-equation k(epsilon) model with a Reynolds-averaged-Navier Stokes/large-eddy-simulation transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier Stokes model in which the unresolved kinetic energy parameter f(sub k) is allowed to vary as a function of grid spacing and the turbulence length scale. This parameter is estimated based on a novel two-stage procedure to efficiently estimate the level of scale resolution possible for a given flow on a given grid for partially averaged Navier Stokes. It has been found that the prescribed scale resolution can play a major role in obtaining accurate flow solutions. The parameter f(sub k) varies between zero and one and is equal to one in the viscous sublayer and when the Reynolds-averaged Navier Stokes turbulent viscosity becomes smaller than the large-eddy-simulation viscosity. The formulation, usage methodology, and validation examples are presented to demonstrate the enhancement of PAB3D's time-accurate turbulence modeling capabilities. The accurate simulations of flow and turbulent quantities will provide a valuable tool for accurate jet noise predictions. Solutions from these models are compared with Reynolds-averaged Navier Stokes results and experimental data for high-temperature jet flows. The current results show promise for the capability of hybrid Reynolds-averaged Navier Stokes and large eddy simulation and partially averaged Navier Stokes in simulating such flow phenomena.

Abdol-Hamid, Khaled S.↗
NASA NTRS2009

Magnetospheric Multiscale Mission (MMS) Phase 2B Navigation Performance

The Magnetospheric Multiscale (MMS) formation flying mission, which consists of four spacecraft flying in a tetrahedral formation, has challenging navigation requirements associated with determining and maintaining the relative separations required to meet the science requirements. The baseline navigation concept for MMS is for each spacecraft to independently estimate its position, velocity and clock states using GPS pseudorange data provided by the Goddard Space Flight Center-developed Navigator receiver and maneuver acceleration measurements provided by the spacecraft's attitude control subsystem. State estimation is performed onboard in real-time using the Goddard Enhanced Onboard Navigation System flight software, which is embedded in the Navigator receiver. The current concept of operations for formation maintenance consists of a sequence of two maintenance maneuvers that is performed every 2 weeks. Phase 2b of the MMS mission, in which the spacecraft are in 1.2 x 25 Earth radii orbits with nominal separations at apogee ranging from 30 km to 400 km, has the most challenging navigation requirements because, during this phase, GPS signal acquisition is restricted to less than one day of the 2.8-day orbit. This paper summarizes the results from high-fidelity simulations to determine if the MMS navigation requirements can be met between and immediately following the maintenance maneuver sequence in Phase 2b.

Scaperoth, Paige Thomas↗
NASA NTRS2013

Orbital Manuvering System Design and Performance For the Magnetosperic Multiscale Constellation

The Magnetospheric Multiscale (MMS) mission, launched on March 13, 2015, is the fourth mission of NASA's Solar Terrestrial Probe program. The MMS mission consists of four identically instrumented observatories that function as a constellation to provide the first definitive study of magnetic reconnection in space. Since it is frequently desirable to isolate electric and magnetic field sensors from stray effects caused by the spacecraft's core-body, the suite of instruments on MMS includes six radial and two axial instrument-booms with deployed lengths ranging from 5-60 meters (see Figure 1). The observatory is spin-stabilized about its positive z-axis with a nominal rate slightly above 3 rev/min (RPM). The spin is also used to maintain tension in the four radial wire-booms. Each observatory's Attitude Control System (ACS) consists of digital sun sensors, star cameras, accelerometers, and mono-propellant hydrazine thrusters-responsible for orbital adjustments, attitude control, and spin adjustments. The sections that follow describe performance requirements, the hardware and algorithms used for 6-DOF estimation, and then similarly for 6-DOF control. The paper concludes with maneuver performance based on both simulated and on-orbit telem.

Determination and Control↗
NASA NTRS2011

Magnetic Reconnection: Theoretical and Observational Perspectives: Preface

Magnetic reconnection is a fundamental plasma-physical process by which energy stored in a magnetic field is converted, often explosively, into heat and the kinetic energy of the charged particles that constitute the plasma. It occurs in a variety of astrophysical settings, ranging from the solar corona to pulsar magnetospheres and winds, as well as in laboratory fusion experiments, where it is responsible for sawtooth crashes. First proposed by R.G. Giovanelli in the late I 940s as the mechanism responsible for solar flares, magnetic reconnection was invoked at the beginning of the space age to explain not just solar flares but also the transfer of energy, mass, and momentum from the solar wind to Earth's magnetosphere and the subsequent storage and release of the transferred energy in the magnetotai\. During the half century or so that has followed the seminal theoretical works by J.W. Dungey, P.A. Sweet, E.N. Parker, and H.E. Petschek, in-situ measurements by Earth-orbiting satellites and remote-sensing observations of the solar corona have provided a growing body of evidence for the occurrence of reconnection at the Sun, in the solar wind, and in the near-Earth space environment. The last thirty years have also seen the development of laboratory reconnection experiments at a number of institutions. In parallel with the efforts of experimentalists in both space and laboratory plasma physics, theorists have investigated, analytically and with the help of increasingly powerful MHD, hybrid, and kinetic numerical simulations, the structure of the diffusion region, the factors controlling the rate, onset, and cessation of reconnection, and the detailed physics that enables the demagnetization of the ions and electrons and the topological reconfiguration of the magnetic field. Moreover, the scope of theoretical reconnection studies has been extended well beyond solar system and laboratory plasmas to include more exotic astrophysical plasma systems whose strong (10(exp 14)-10(exp 15) G) magnetic fields require that models of reconnection in these systems incorporate quantum electrodynamical, special relativistic, and radiative effects. The papers collected in this topical issue of Space Science Reviews cover different aspects of recent theoretical and observational work on magnetic reconnection in solar and space physics, astrophysics, and laboratory plasma physics. They derive from presentations given at a workshop on magnetic reconnection held in the Yosemite National Park, February 8-12,2010. The intent of the workshop was to stimulate, through a combination of tutorial talks, shorter focused talks, and extensive informal discussions, an interdisciplinary dialogue among members of the different research communities working on the problem of magnetic reconnection. One of the motivating considerations for holding the workshop was its relevance to NASA's Magnetospheric Multiscale (MMS) mission, scheduled for launch in 2014. The four identically instrumented MMS spacecraft are designed to study reconnect ion in Earth's magnetosphere and, specifically, to probe the electron diffusion region in order to determine the microphysical processes that enable the change in the topology of the magnetic field. Building on the achievements of the multi spacecraft Cluster and THEMIS missions, MMS will use the magnetosphere as an astrophysical plasma laboratory in which to test, through in-situ measurement of the plasma, energetic particles, and electric and magnetic fields, various models and theories that have emerged during the past twenty years, a period of extraordinarily productive theoretical and observational work.

Lewis, W. S.↗
NASA NTRS2021

An Encounter with the Ion and Electron Diffusion Regions at a Flapping and Twisted Tail Current Sheet

We analyze data returned by the Magnetospheric Multiscale mission (MMS) constellation during a rapid (~1.5 s) traversal of a flapping and reconnecting current sheet (CS) in the near-Earth magnetotail (X ~ -20 RE). The CS was highly tilted, with its normal pointing strongly duskward. Its extreme thinness was confirmed by a curvature analysis of the magnetic field lines. The event was associated with a guide field of 8% of the reconnecting components. From the pitch angle distributions of low-energy electrons we infer that the crossing occurred earthward of the X-line. Traveling practically normal to the CS, MMS encountered an ion diffusion region (IDR) in which was embedded an electron diffusion region (EDR). IDR signatures included breaking of the ion frozen-in condition in the presence of Hall B and E fields. EDR signatures included a strong out of-plane current associated with a super-Alfvenic electron jet, evidence of positive energy transfer, and a temperature anisotropy (TeII > Te1) which disappeared at the field reversal. Derived scale sizes normal to the CS are: ~6.9 de (EDR) and ~0.4 di (IDR) (40 and 100 km). We estimate the average dimensionless reconnection rate as 0.077 + 0.050. The observations and inferences are supported by PIC numerical simulations. We find very good agreement in the reconnection rates. We also discuss the effects of asymmetries in the density, temperature and magnetic field strength on the Hall fields and length of the outflow jets. The event is associated with a substorm onset which began 7 min after the MMS observations.

C J Farrugia↗
NASA NTRS2004

Numerical Study of Flow Past a Circular Cylinder Using RANS, Hybrid RANS/LES and PANS Formulations

Two multiscale type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds Averaged Navier-Stokes (RANS) equations. The first scheme is a hybrid RANS/LES model utilizing the two-equation (k(sub epsilon)) model with a RANS/LES transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier-Stokes (PANS) model, where the unresolved kinetic energy parameter (f(sub k)) is allowed to vary as a function of grid spacing and the turbulence length scale. Solutions from these models are compared to RANS results and experimental data for a stationary and rotating cylinder. The parameter f(sub k) varies between zero and one and has the characteristic to be equal to one in the viscous sub layer, and when the RANS turbulent viscosity becomes smaller than the LES viscosity. The formulation, usage methodology, and validation example are presented to demonstrate the enhancement of PAB3D's time-accurate and turbulence modeling capabilities. The models are compared to RANS results and experimental data for turbulent separated flows (TS) and laminar separated flows (LS) around stationary and rotating cylinders. For a stationary cylinder, the TS case is accurately simulated using the general two-equation k(sub epsilon) turbulence model (eddy-viscosity model). PAB3D accurately predicts the drag coefficient (CD), lift coefficient (CL) and the Strouhal number (St). The LS case was a challenge for the RANS computation with an eddy-viscosity turbulence model. The RANS/LES and PANS performed well and showed marked improvements over the RANS solution. The modified PANS model was the most accurate. For the rotating cylinder, the spin ratio varied from zero to one, and the PANS results were in good agreement with published experimental data. RANS/LES and PANS capture both temporal and spatial fluctuations and produce large-scale structures that do not occur in the RANS simulation. The current results show promise for the capability of PANS in simulating unsteady and complex flow phenomena.

Elmiligui, Alaa A.↗
NASA NTRS2013

Higher Surface Ozone Concentrations Over the Chesapeake Bay than Over the Adjacent Land: Observations and Models from the DISCOVER-AQ and CBODAQ Campaigns

Air quality models, such as the Community Multiscale Air Quality (CMAQ) model, indicate decidedly higher ozone near the surface of large interior water bodies, such as the Great Lakes and Chesapeake Bay. In order to test the validity of the model output, we performed surface measurements of ozone (O3) and total reactive nitrogen (NOy) on the 26-m Delaware II NOAA Small Research Vessel experimental (SRVx), deployed in the Chesapeake Bay for 10 daytime cruises in July 2011 as part of NASA's GEO-CAPE CBODAQ oceanographic field campaign in conjunction with NASA's DISCOVER-AQ air quality field campaign. During this 10-day period, the EPA O3 regulatory standard of 75 ppbv averaged over an 8-h period was exceeded four times over water while ground stations in the area only exceeded the standard at most twice. This suggests that on days when the Baltimore/Washington region is in compliance with the EPA standard, air quality over the Chesapeake Bay might exceed the EPA standard. Ozone observations over the bay during the afternoon were consistently 10-20% higher than the closest upwind ground sites during the 10-day campaign; this pattern persisted during good and poor air quality days. A lower boundary layer, reduced cloud cover, slower dry deposition rates, and other lesser mechanisms, contribute to the local maximum of ozone over the Chesapeake Bay. Observations from this campaign were compared to a CMAQ simulation at 1.33 km resolution. The model is able to predict the regional maximum of ozone over the Chesapeake Bay accurately, but NOy concentrations are significantly overestimated. Explanations for the overestimation of NOy in the model simulations are also explored

boundary layer↗
NASA NTRS2023

Comparison of Dust Optical Depth From Multi-Sensor Products and MONARCH (Multiscale Online Non-hydrostatic AtmospheRe CHemistry) Dust Reanalysis Over North Africa, the Middle East, and Europe

Aerosol reanalysis datasets are model-based, observationally constrained, continuous 3D aerosol fields with a relatively high temporal frequency that can be used to assess aerosol variations and trends, climate effects, and impacts on socioeconomic sectors, such as health. Here we compare and assess the recently published MONARCH (Multiscale Online Non-hydrostatic AtmospheRe CHemistry) high-resolution regional desert dust reanalysis over northern Africa, the Middle East, and Europe (NAMEE) with a combination of ground-based observations and space-based dust retrievals and products. In particular, we compare the total and coarse dust optical depth (DOD) from the new reanalysis with DOD products derived from MODIS (MODerate resolution Imaging Spectroradiometer), MISR (Multi-angle Imaging SpectroRadiometer), and IASI (Infrared Atmospheric Sounding Interferometer) spaceborne instruments. Despite the larger uncertainties, satellite-based datasets provide a better geographical coverage than ground-based observations, and the use of different retrievals and products allows at least partially overcoming some single-product weaknesses in the comparison. Nevertheless, limitations and uncertainties due to the type of sensor, its operating principle, its sensitivity, its temporal and spatial resolution, and the methodology for retrieving or further deriving dust products are factors that bias the reanalysis assessment. We, therefore, also use ground-based DOD observations provided by 238 stations of the AERONET (AErosol RObotic NETwork) located within the NAMEE region as a reference evaluation dataset. In particular, prior to the reanalysis assessment, the satellite datasets were evaluated against AERONET, showing moderate underestimations in the vicinities of dust sources and downwind regions, whereas small or significant overestimations, depending on the dataset, can be found in the remote regions. Taking these results into consideration, the MONARCH reanalysis assessment shows that total and coarse-DOD simulations are consistent with satellite- and ground-based data, qualitatively capturing the major dust sources in the area in addition to the dust transport patterns. Moreover, the MONARCH reanalysis reproduces the seasonal dust cycle, identifying the increased dust activity that occurred in the NAMEE region during spring and summer. The quantitative comparison between the MONARCH reanalysis DOD and satellite multi-sensor products shows that the reanalysis tends to slightly overestimate the desert dust that is emitted from the source regions and underestimate the transported dust over the outflow regions, implying that the model's removal of dust particles from the atmosphere, through deposition processes, is too effective. More specifically, small positive biases are found over the Sahara desert (0.04) and negative biases over the Atlantic Ocean and the Arabian Sea (−0.04), which constitute the main pathways of the long-range dust transport. Considering the DOD values recorded on average there, such discrepancies can be considered low, as the low relative bias in the Sahara desert (< 50 %) and over the adjacent maritime regions (< 100 %) certifies. Similarly, over areas with intense dust activity, the linear correlation coefficient between the MONARCH reanalysis simulations and the ensemble of the satellite products is significantly high for both total and coarse DOD, reaching 0.8 over the Middle East, the Atlantic Ocean, and the Arabian Sea and exceeding it over the African continent. Moreover, the low relative biases and high correlations are associated with regions for which large numbers of observations are available, thus allowing for robust reanalysis assessment.

Michail Mytilinaios↗
NASA NTRS2019

Revealing the Multiscale Nature of Turbulence in Space Plasmas with an Innovative Swarm of Spacecraft

Abstract: Turbulence is fundamentally a multiscale phenomena, with energy, mass, and momentum being transported across both spatial and temporal scales. In a magnetized collisionless plasma, there are unresolved questions about the structure of the turbulent transfer of energy, as well as how the energy is extracted from the cascade and dissipated as heat in the constituent charged particles. The plasma in the solar wind acts as an accessible natural laboratory to study these processes, and much progress has indeed been made since the dawn of the space age in understanding turbulence via in situ observations of turbulent plasmas near Earth. However, to date, these observations have been limited to a single, or at best, a tight cluster of points, leading to ambiguities in at what scales energy is contained, how it is transported, and by what mechanism it is dissipated. In this presentation, we describe a heliophysics mission concept aimed at understanding turbulence that is enabled by a swarm of small satellites. The proposed "HelioSwarm" mission will measure turbulent fields and flows and charged particles simultaneously at many points spanning size and time scales from the fluid to sub-ion regime. In doing so, we will be able to disentangle how the turbulence depends on time and space, directly observe the change in internal energy in the plasma, and definitively capture the dynamic relation between turbulence and structures. While the processes under examination are universal, arising throughout our solar system and universe, they are difficult to reproduce in either terrestrial laboratories or numerical simulations, meaning that the only a multipoint observatory in the near-Earth heliosphere will be able to study them in sufficient detail to discern the underlying physics. In this talk, we highlight the enabling role that small satellites play in providing closure on these long-standing but critically important science questions.

Spence, H. E.↗
NASA NTRS2026

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗
NASA NTRS2021

Harmonized Emissions Component (HEMCO) 3.0 as a Versatile Emissions Component for Atmospheric Models: Application in the GEOS-Chem, NASA GEOS, WRF-GC, CESM2, NOAA GEFS-Aerosol, and NOAA UFS Models

Emissions are a central component of atmospheric chemistry models. The Harmonized Emissions Component (HEMCO) is a software component for computing emissions from a user-selected ensemble of emission inventories and algorithms. It allows users to re-grid, combine, overwrite, subset, and scale emissions from different inventories through a configuration file and with no change to the model source code. The configuration file also maps emissions to model species with appropriate units. HEMCO can operate in offline stand-alone mode, but more importantly it provides an online facility for models to compute emissions at runtime. HEMCO complies with the Earth System Modeling Framework (ESMF) for portability across models. We present a new version here, HEMCO 3.0, that features an improved three-layer architecture to facilitate implementation into any atmospheric model and improved capability for calculating emissions at any model resolution including multiscale and unstructured grids. The three-layer architecture of HEMCO 3.0 includes (1) the Data Input Layer that reads the configuration file and accesses the HEMCO library of emission inventories and other environmental data, (2) the HEMCO Core that computes emissions on the user-selected HEMCO grid, and (3) the Model Interface Layer that re-grids (if needed) and serves the data to the atmospheric model and also serves model data to the HEMCO Core for computing emissions dependent on model state (such as from dust or vegetation). The HEMCO Core is common to the implementation in all models, while the Data Input Layer and the Model Interface Layer are adaptable to the model environment. Default versions of the Data Input Layer and Model Interface Layer enable straightforward implementation of HEMCO in any simple model architecture, and options are available to disable features such as re-gridding that may be done by independent couplers in more complex architectures. The HEMCO library of emission inventories and algorithms is continuously enriched through user contributions so that new inventories can be immediately shared across models. HEMCO can also serve as a general data broker for models to process input data not only for emissions but for any gridded environmental datasets. We describe existing implementations of HEMCO 3.0 in (1) the GEOS-Chem “Classic” chemical transport model with shared-memory infrastructure, (2) the high-performance GEOS-Chem (GCHP) model with distributed-memory architecture, (3) the NASA GEOS Earth System Model (GEOS ESM), (4) the Weather Research and Forecasting model with GEOS-Chem (WRF-GC), (5) the Community Earth System Model Version 2 (CESM2), and (6) the NOAA Global Ensemble Forecast System – Aerosols (GEFS-Aerosols), as well as the planned implementation in the NOAA Unified Forecast System (UFS). Implementation of HEMCO in CESM2 contributes to the Multi-Scale Infrastructure for Chemistry and Aerosols (MUSICA) by providing a common emissions infrastructure to support different simulations of atmospheric chemistry across scales.

Haipeng Lin↗
NASA NTRS2016

Tool for Generation of MAC/GMC Representative Unit Cell for CMC/PMC Analysis

This document describes a recently developed analysis tool that enhances the resident capabilities of the Micromechanics Analysis Code with the Generalized Method of Cells (MAC/GMC) 4.0. This tool is especially useful in analyzing ceramic matrix composites (CMCs), where higher fidelity with improved accuracy of local response is needed. The tool, however, can be used for analyzing polymer matrix composites (PMCs) as well. MAC/GMC 4.0 is a composite material and laminate analysis software developed at NASA Glenn Research Center. The software package has been built around the concept of the generalized method of cells (GMC). The computer code is developed with a user friendly framework, along with a library of local inelastic, damage, and failure models. Further, application of simulated thermomechanical loading, generation of output results, and selection of architectures to represent the composite material have been automated to increase the user friendliness, as well as to make it more robust in terms of input preparation and code execution. Finally, classical lamination theory has been implemented within the software, wherein GMC is used to model the composite material response of each ply. Thus, the full range of GMC composite material capabilities is available for analysis of arbitrary laminate configurations as well. The primary focus of the current effort is to provide a graphical user interface (GUI) capability that generates a number of different user-defined repeating unit cells (RUCs). In addition, the code has provisions for generation of a MAC/GMC-compatible input text file that can be merged with any MAC/GMC input file tailored to analyze composite materials. Although the primary intention was to address the three different constituents and phases that are usually present in CMCs-namely, fibers, matrix, and interphase-it can be easily modified to address two-phase polymer matrix composite (PMC) materials where an interphase is absent. Currently, the tool capability includes generation of RUCs for square packing, hexagonal packing, and random fiber packing as well as RUCs based on actual composite micrographs. All these options have the fibers modeled as having a circular cross-sectional area. In addition, a simplified version of RUC is provided where the fibers are treated as having a square cross section and are distributed randomly. This RUC facilitates a speedy analysis using the higher fidelity version of GMC known as HFGMC. The first four mentioned options above support uniform subcell discretization. The last one has variable subcell sizes due to the primary intention of keeping the RUC size to a minimum to gain the speed ups using the higher fidelity version of MAC. The code is implemented within the MATLAB (The Mathworks, Inc., Natick, MA) developmental framework; however, a standalone application that does not need a priori MATLAB installation is also created with the aid of the MATLAB compiler.

Materials Engineering↗