Search NASA⌕ Search

SEARCH · Search NASA

Results for “box modeling”

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 19 records

The Messy Aerosol Submodel MADE3 (v2.0b): Description and a Box Model Test

We introduce MADE3 (Modal Aerosol Dynamics model for Europe, adapted for global applications, 3rd generation), an aerosol dynamics submodel for application within the MESSy framework (Modular Earth Submodel System). MADE3 builds on the predecessor aerosol submodels MADE and MADE-in. Its main new features are the explicit representation of coarse particle interactions both with other particles and with condensable gases, and the inclusion of hydrochloric acid (HCl)chloride (Cl) partitioning between the gas and condensed phases. The aerosol size distribution is represented in the new submodel as a superposition of nine lognormal modes: one for fully soluble particles, one for insoluble particles, and one for mixed particles in each of three size ranges (Aitken, accumulation, and coarse mode size ranges). In order to assess the performance of MADE3 we compare it to its predecessor MADE and to the much more detailed particle-resolved aerosol model PartMC-MOSAIC in a box model simulation of an idealized marine boundary layer test case. MADE3 and MADE results are very similar, except in the coarse mode, where the aerosol is dominated by sea spray particles. Cl is reduced in MADE3 with respect to MADE due to the HClCl partitioning that leads to Cl removal from the sea spray aerosol in our test case. Additionally, aerosol nitrate concentration is higher in MADE3 due to the condensation of nitric acid on coarse particles. MADE3 and PartMC- MOSAIC show substantial differences in the fine particle size distributions (sizes about 2 micrometers) that could be relevant when simulating climate effects on a global scale. Nevertheless, the agreement between MADE3 and PartMC-MOSAIC is very good when it comes to coarse particle size distribution, and also in terms of aerosol composition. Considering these results and the well-established ability of MADE in reproducing observed aerosol loadings and composition, MADE3 seems suitable for application within a global model.

Mixing state↗

Box models for the evolution of atmospheric oxygen: an update

A simple 3-box model of the atmosphere/ocean system is used to describe the various stages in the evolution of atmospheric oxygen. In Stage I, which probably lasted until redbeds began to form about 2.0 Ga ago, the Earth's surface environment was generally devoid of free O2, except possibly in localized regions of high productivity in the surface ocean. In Stage II, which may have lasted for less than 150 Ma, the atmosphere and surface ocean were oxidizing, while the deep ocean remained anoxic. In Stage III, which commenced with the disappearance of banded iron formations around 1.85 Ga ago and has lasted until the present, all three surface reservoirs contained appreciable amounts of free O2. Recent and not-so-recent controversies regarding the abundance of oxygen in the Archean atmosphere are identified and discussed. The rate of O2 increase during the Middle and Late Proterozoic is identified as another outstanding question.

NASA Program Exobiology↗

Electron propagation in the leaky box model with a truncated pathlength distribution

A study of electron propagation using the leaky box model is discussed. It is shown that a truncated pathlength distribution due to a lack of nearby sources is responsible for the steepening of the electron spectrum. The electron spectrum is broken into three regions: a low-energy region where electron storage is dominated by escape, a medium energy region in which the electron energy loss lifetime is sufficiently short to dominate propagation, but not so short as to prevent electrons from propagating throughout the storage region, and a third energy region in which the energy loss lifetime is shorter than the time it takes for cosmic rays to diffuse to earth from the nearest source. The asymptotic slope of the electron spectrum is shown to be steeper than that of the injection spectrum by more than one power of energy.

Mauger, B. G.↗

Global Sensitivity Analyses for Test Planning with Black-Box Models for Mars Sample Return

This work describes sensitivity analyses performed on complex black-box models used to support experimental test planning under limited resources in the context of the Mars Sample Return program, which aims at bringing to Earth rock and atmospheric samples from Mars. We develop a systematic workflow that allows the analysts to simultaneously obtain quantitative insights on key drivers of uncertainty, on the direction of impact, and the presence of interactions. We apply novel optimal transport-based global sensitivity measures to tackle the multivariate nature of the output. On the modeling side, we apply multi-fidelity techniques that leverage low-fidelity models to speed up the calculations and make up for the limited amount of high-fidelity samples, while keeping these in the loop for accuracy guarantees. The sensitivity analysis reveals insights useful for the analysts to understand the model's behavior and identify the factors to focus on during testing in order to maximize the value of information extracted from them to ensure mission success when limited resources are available.

Giuseppe Cataldo↗

Comparison of distributed reacceleration and leaky-box models of cosmic-ray abundances (Z = 3-28)

A large collection of elemental and isotopic cosmic-ray data has been analyzed using the leaky-box transport model with and without reacceleration in the interstellar medium. Abundances of isotopes and elements with charges Z = 3-28 and energies E = 10 MeV/nucleon-1 TeV/nucleon were explored. Our results demonstrate that reacceleration models make detailed and accurate predictions with the same number of parameters or fewer as standard leaky-box models. Ad hoc fitting parameters in the standard model are replaced by astrophysically significant reacceleration parameters. Distributed reacceleration models explain the peak in secondary-to-primary ratios around 1 GeV/nucleon. They diminish the discrepancy between rigidity-dependent leakage and energy-independent anisotropy. They also offer the possibility of understanding isotopic anomalies at low energy.

Letaw, John R.↗

Coupled Particle Transport and Pattern Formation in a Nonlinear Leaky-Box Model

Effects of particle-particle coupling on particle characteristics in nonlinear leaky-box type descriptions of the acceleration and transport of energetic particles in space plasmas are examined in the framework of a simple two-particle model based on the Fokker-Planck equation in momentum space. In this model, the two particles are assumed coupled via a common nonlinear source term. In analogy with a prototypical mathematical system of diffusion-driven instability, this work demonstrates that steady-state patterns with strong dependence on the magnetic turbulence but a rather weak one on the coupled particles attributes can emerge in solutions of a nonlinearly coupled leaky-box model. The insight gained from this simple model may be of wider use and significance to nonlinearly coupled leaky-box type descriptions in general.

Barghouty, A. F.↗

Why do leaky-box models work so fine?

By introducing the concept of the age distribution of cosmic rays it is possible to decouple spatial from momentum transport, and simple leaky-box type equations result. The influence of spatial inhomogeneities, geometries, and source distributions enters the spatially homogeneous, infinite (i.e., leaky box) problem through appropriate mean lifetimes. A precise prescription of how to obtain these mean lifetimes, i.c., for comparison with data measured in the vicinity of the solar system they have to be calculated from the age distribution at the spatial position of the observer.

Schlickeiser, R.↗

Stochastic acceleration of cometary pickup ions - The classic leaky box model

The acceleration of cometary pickup ions by magnetohydrodynamic waves at P/Giacobini-Zinner is examined in a model where acceleration predominantly occurs downstream of the bow shock throughout the cometosheath where intense magnetic turbulence exists. The mean free path for scattering by the magnetic fluctuations in this region is less than the characteristic dimension of the cometosheath so that pickup ions are rendered isotropic and energized by a modest amount in the process. This principal loss mechanism for the ions is spatial diffusion out of the acceleration region moderated by the self-same accelerating waves. This particular feature constrains the model in a way that the predicted ion spectrum is uniquely determined by the power spectrum of the magnetic turbulence. At both P/Giacobini-Zinner and P/Halley, the turbulence is non-Kolmogoroff with a spectral index of 2 resulting in an ion spectral behavior that is approximately an exponential in ion speed, consistent with a recent analysis of the Giacobini-Zinner data (see Richardson et al.).

Barbosa, D. D.↗

Analytic solution of the box diffusion model for a global ocean

The paper presents an analytic solution for the temperature of the mixed layer for a box diffusion model of a global ocean. It is shown that the exact solution approaches the appropriate limits of a mixed layer ocean model and of a diffusive ocean without a mixed layer on top. It is also shown that, over a large range of parameters of the box diffusion model, the difference between the temperature calculated using a box diffusion model and the temperature calculated using a purely diffusive ocean is buried in the noise.

Lebedeff, Sergej A.↗

A Fast and Efficient Version of the TwO-Moment Aerosol Sectional (TOMAS) Global Aerosol Microphysics Model

This study develops more computationally efficient versions of the TwO-Moment Aerosol Sectional (TOMAS) microphysics algorithms, collectively called Fast TOMAS. Several methods for speeding up the algorithm were attempted, but only reducing the number of size sections was adopted. Fast TOMAS models, coupled to the GISS GCM II-prime, require a new coagulation algorithm with less restrictive size resolution assumptions but only minor changes in other processes. Fast TOMAS models have been evaluated in a box model against analytical solutions of coagulation and condensation and in a 3-D model against the original TOMAS (TOMAS-30) model. Condensation and coagulation in the Fast TOMAS models agree well with the analytical solution but show slightly more bias than the TOMAS-30 box model. In the 3-D model, errors resulting from decreased size resolution in each process (i.e., emissions, cloud processing wet deposition, microphysics) are quantified in a series of model sensitivity simulations. Errors resulting from lower size resolution in condensation and coagulation, defined as the microphysics error, affect number and mass concentrations by only a few percent. The microphysics error in CN70CN100 (number concentrations of particles larger than 70100 nm diameter), proxies for cloud condensation nuclei, range from 5 to 5 in most regions. The largest errors are associated with decreasing the size resolution in the cloud processing wet deposition calculations, defined as cloud-processing error, and range from 20 to 15 in most regions for CN70CN100 concentrations. Overall, the Fast TOMAS models increase the computational speed by 2 to 3 times with only small numerical errors stemming from condensation and coagulation calculations when compared to TOMAS-30. The faster versions of the TOMAS model allow for the longer, multi-year simulations required to assess aerosol effects on cloud lifetime and precipitation.

Lee, Yunha↗

Spatial Pattern of Cell Damage in Tissue from Heavy Ions

A new Monte Carlo algorithm was developed that can model passage of heavy ions in a tissue, and their action on the cellular matrix for 2- or 3-dimensional cases. The build-up of secondaries such as projectile fragments, target fragments, other light fragments, and delta-rays was simulated. Cells were modeled as a cell culture monolayer in one example, where the data were taken directly from microscopy (2-d cell matrix). A simple model of tissue was given as abstract spheres with close approximation to real cell geometries (3-d cell matrix), as well as a realistic model of tissue was proposed based on microscopy images. Image segmentation was used to identify cells in an irradiated cell culture monolayer, or slices of tissue. The cells were then inserted into the model box pixel by pixel. In the case of cell monolayers (2-d), the image size may exceed the modeled box size. Such image was is moved with respect to the box in order to sample as many cells as possible. In the case of the simple tissue (3-d), the tissue box is modeled with periodic boundary conditions, which extrapolate the technique to macroscopic volumes of tissue. For real tissue, specific spatial patterns for cell apoptosis and necrosis are expected. The cell patterns were modeled based on action cross sections for apoptosis and necrosis estimated based on BNL data, and other experimental data.

Ponomarev, Artem L.↗

Intercomparison of the representations of the atmospheric chemistry of pre-industrial methane and ozone in earth system and other global chemistry-transport models

An intercomparison has been set up to study the representation of the atmospheric chemistry of the pre-industrial troposphere in earth system and other global tropospheric chemistry-transport models. The intercomparison employed a constrained box model and utilised tropospheric trace gas composition data for the pre-industrial times at ninety mid-latitude surface locations. Incremental additions of four organic compounds: methane, ethane, acetone and propane, were used to perturb the constrained box model and generate responses in hydroxyl radicals and tropospheric ozone at each location and with each chemical mechanism. Although the responses agreed well across the chemical mechanisms from the selected earth system and other global tropospheric chemistry-transport models, there were differences in the detailed responses between the chemical mechanisms that could be tracked down by sensitivity analysis to differences in the representation of C1–C3 chemistry. Inter-mechanism ranges in NOx compensation points were about 0.17 ± 0.12 when expressed relative to the inter-mechanism average. Monte Carlo uncertainty analysis carried out with a single chemical mechanism put the intra-mechanism range a factor of three higher at 0.50 ± 0.12. Similar differences between inter-mechanism and intra-mechanism ranges were found for hydroxyl radical depletion but were up to a factor of six wider for ozone formation from incremental additions of organic compounds. The cause of the discrepancies between the inter- and intra-mechanism ranges was found to be the large uncertainties that are present in the laboratory determinations of the rate coefficients and product channel branching ratios of some key chemical reactions involving organic peroxy radicals and hydroperoxides. Whilst these large uncertainties are present in the laboratory determinations, there will be irreducible uncertainties in the predictions from the earth system and other chemistry-transport models of methane and tropospheric ozone trends since pre-industrial times and hence their contributions to the radiative forcing of climate change. Further definitive laboratory studies of the reaction rates and product yields of the reactions of the simple organic peroxy radicals and hydroperoxides are required to resolve and reduce current uncertainties in earth system and chemistry-transport model predictions.

atmospheric chemistry↗

Thermal Analysis of a High-Speed Aircraft Wing Using p-Version Finite Elements

This paper presents the results of conceptual level thermal analyses of a High Speed Civil Transport (HSCT) wing using p-version finite elements. The work was motivated by a thermal analysis of a HSCT wing structure which showed the importance of radiation heat transfer throughout the structure. The analysis also showed that refining a traditional finite element mesh to accurately capture the temperature distribution on the internal structure led to very large meshes with unacceptably long execution times. Further study indicated using p-version finite elements might improve computation performance for this class of problem. Methods for determining internal radiation heat transfer were then developed and demonstrated on test problems representative of the geometry found in an aircraft wing structure. This paper presents the results of the application of these new methods to the analysis of a high speed aircraft wing. Results for both a wing box model as well as a full wing model are presented. 'Me reduced wing box model allows for a comparison of the traditional finite element method with mesh refinement (h-refinement) to the new p-version finite elements while the full wing model demonstrates the applicability and efficiency of p-version finite elements for large models.

Gould, Dana C.↗

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↗

Methane and Environmental Change during the Paleocene-Eocene Thermal Maximum (PETM): Modeling the PETM Onset as a Two-stage Event

An atmospheric CH4 box model coupled to a global carbon cycle box model is used to constrain the carbon emission associated with the PETM and assess the role of CH4 during this event. A range of atmospheric and oceanic emission scenarios representing different amounts, rates, and isotopic signatures of emitted carbon are used to model the PETM onset. The first 3 kyr of the onset, a pre-isotope excursion stage, is simulated by the atmospheric release of 900 to 1100 Pg C CH4 with a delta C-13 of -22 to - 30 %. For a global average warming of 3 deg C, a release of CO2 to the ocean and CH4 to the atmosphere totalling 900 to 1400 Pg C, with a delta C-13 of -50 to -60%, simulates the subsequent 1 -kyr isotope excursion stage. To explain the observations, the carbon must have been released over at most 500 years. The first stage results cannot be associated with any known PETM hypothesis. However, the second stage results are consistent with a methane hydrate source. More than a single source of carbon is required to explain the PETM onset.

Carozza, David A.↗

Cosmic ray antiprotons in closed galaxy model

The flux of secondary antiprotons expected for the leaky-box model was calculated as well as that for the closed galaxy model of Peters and Westergard (1977). The antiproton/proton ratio observed at several GeV is a factor of 4 higher than the prediction for the leaky-box model but is consistent with that predicted for the closed galaxy model. New low energy data is not consistent with either model. The possibility of a primary antiproton component is discussed.

Protheroe, R.↗

Internal Structural Design of the Common Research Model Wing Box for Aeroelastic Tailoring

This work explores the use of alternative internal structural designs within a full-scale wing box structure for aeroelastic tailoring, with a focus on curvilinear spars, ribs, and stringers. The baseline wing model is a fully-populated, cantilevered wing box structure of the Common Research Model (CRM). Metrics of interest include the wing weight, the onset of dynamic flutter, and the static aeroelastic stresses. Twelve parametric studies alter the number of internal structural members along with their location, orientation, and curvature. Additional evaluation metrics are considered to identify design trends that lead to lighter-weight, aeroelastically stable wing designs. The best designs of the individual studies are compared and discussed, with a focus on weight reduction and flutter resistance. The largest weight reductions were obtained by removing the inner spar, and performance was maintained by shifting stringers forward and/or using curvilinear ribs: 5.6% weight reduction, a 13.9% improvement in flutter speed, but a 3.0% increase in stress levels. Flutter resistance was also maintained using straight-rotated ribs although the design had a 4.2% lower flutter speed than the curved ribs of similar weight and stress levels were higher. For some configurations, the differences between curved and straight ribs were smaller, which provides motivation for future optimization-based studies to fully exploit the trade-offs.

Jutte, Christine V.↗