Search NASA⌕ Search

SEARCH · Search NASA

Results for “Modeling & Simulation (M&S)”

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 271 records · Page 15

Factors Controlling the Properties of Multi-Phase Arctic Stratocumulus Clouds

The 2004 Multi-Phase Arctic Cloud Experiment (M-PACE) IOP at the ARM NSA site focused on measuring the properties of autumn transition-season arctic stratus and the environmental conditions controlling them, including concentrations of heterogeneous ice nuclei. Our work aims to use a large-eddy simulation (LES) code with embedded size-resolved aerosol and cloud microphysics to identify factors controlling multi-phase arctic stratus. Our preliminary simulations of autumn transition-season clouds observed during the 1994 Beaufort and Arctic Seas Experiment (BASE) indicated that low concentrations of ice nuclei, which were not measured, may have significantly lowered liquid water content and thereby stabilized cloud evolution. However, cloud drop concentrations appeared to be virtually immune to changes in liquid water content, indicating an active Bergeron process with little effect of collection on drop number concentration. We will compare these results with preliminary simulations from October 8-13 during MPACE. The sensitivity of cloud properties to uncertainty in other factors, such as large-scale forcings and aerosol profiles, will also be investigated. Based on the LES simulations with M-PACE data, preliminary results from the NASA GlSS single-column model (SCM) will be used to examine the sensitivity of predicted cloud properties to changing cloud drop number concentrations for multi-phase arctic clouds. Present parametrizations assumed fixed cloud droplet number concentrations and these will be modified using M-PACE data.

Fridlind, Ann↗

Design of a Mars Airplane Propulsion System for the Aerial Regional-Scale Environmental Survey (ARES) Mission Concept

The Aerial Regional-Scale Environmental Survey (ARES) is a Mars exploration mission concept that utilizes a rocket propelled airplane to take scientific measurements of atmospheric, surface, and subsurface phenomena. The liquid rocket propulsion system design has matured through several design cycles and trade studies since the inception of the ARES concept in 2002. This paper describes the process of selecting a bipropellant system over other propulsion system options, and provides details on the rocket system design, thrusters, propellant tank and PMD design, propellant isolation, and flow control hardware. The paper also summarizes computer model results of thruster plume interactions and simulated flight performance. The airplane has a 6.25 m wingspan with a total wet mass of 185 kg and has to ability to fly over 600 km through the atmosphere of Mars with 45 kg of MMH / MON3 propellant.

Kuhl, Christopher A.↗

Development of a “Nature Run” for Observing System Simulation Experiments (OSSEs) for Snow Mission Development

Snow is a fundamental component of global and regional water budgets, particularly in mountainous areas and regions downstream that rely on snowmelt for water resources. Land surface models (LSMs) are commonly used to develop spatially distributed estimates of snow water equivalent (SWE) and runoff. However, LSMs are limited by uncertainties in model physics and parameters, among other factors. In this study, we describe the use of model calibration tools to improve snow simulations within the Noah-MP LSM as the first step in an observing system simulation experiment (OSSE). Noah-MP is calibrated against the University of Arizona (UA) SWE product over a western Colorado domain. With spatially varying calibrated parameters, we run calibrated and default Noah-MP simulations for water years 2010–20. By evaluating both simulations against the UA dataset, we show that calibration decreases domain averaged temporal RMSE and bias for snow depth from 0.15 to 0.13 m and from −0.036 to −0.0023 m, respectively, and improves the timing of snow ablation. Increased snow simulation performance also improves estimates of model-simulated runoff in four of six study basins, though only one has statistically significant improvement. Spatially distributed Noah-MP snow parameters perform better than default uniform values. We demonstrate that calibrating variables related to snow albedo calculations and rain–snow partitioning, among other processes, is a necessary step for creating a nature run that reasonably approximates true snow conditions for the OSSEs. Additionally, the inclusion of a snowfall scaling term can address biases in precipitation from meteorological forcing datasets, further improving the utility of LSMs for generating reliable spatiotemporal estimates of snow.

Melissa L. Wrzesien↗

The extended atmosphere and circumstellar environment of the cool evolved star VX Sagittarii as seen by MATISSE

Context. VX Sgr is a cool, evolved, and luminous red star whose stellar parameters are difficult to determine, which affects its classification. Aims. We aim to spatially resolve the photospheric extent as well as the circumstellar environment. Methods. We used interferometric observations obtained with the MATISSE instrument in the L (3–4 μm), M (4.5–5 μm), and N (8–13 μm) bands. We reconstructed monochromatic images using the MIRA software. We used 3D radiation-hydrodynamics simulations carried out with CO5BOLD and a uniform disc model to estimate the apparent diameter and interpret the stellar surface structures. Moreover, we employed the radiative transfer codes OPTIM3D and RADMC3D to compute the spectral energy distribution for the L, M, and N bands, respectively. Results. MATISSE observations unveil, for the first time, the morphology of VX Sgr across the L, M, and N bands. The reconstructed images show a complex morphology with brighter areas whose characteristics depend on the wavelength probed. We measured the angular diameter as a function of the wavelength and showed that the photospheric extent in the L and M bands depends on the opacity through the atmosphere. In addition to this, we also concluded that the observed photospheric inhomogeneities can be interpreted as convection-related surface structures. The comparison in the N band yielded a qualitative agreement between the N-band spectrum and simple dust radiative transfer simulations. However, it is not possible to firmly conclude on the interpretation of the current data because of the difficulty in constraing the model parameters using the limited accuracy of our absolute flux calibration. Conclusions. MATISSE observations and the derived reconstructed images unveil the appearance of VX Sgr’s stellar surface and circumstellar environment across a very large spectral domain for the first time.

stars: atmospheres↗

Application and Evaluation of a Snowmelt Runoff Model in the Tamor River Basin, Eastern Himalaya Using a Markov Chain Monte Carlo (MCMC) Data Assimilation Approach

Previous studies have drawn attention to substantial hydrological changes taking place in mountainous watersheds where hydrology is dominated by cryospheric processes. Modelling is an important tool for understanding these changes but is particularly challenging in mountainous terrain owing to scarcity of ground observations and uncertainty of model parameters across space and time. This study utilizes a Markov Chain Monte Carlo data assimilation approach to examine and evaluate the performance of a conceptual, degree-day snowmelt runoff model applied in the Tamor River basin in the eastern Nepalese Himalaya. The snowmelt runoff model is calibrated using daily streamflow from 2002 to 2006 with fairly high accuracy (average Nash-Sutcliffe metric approx. 0.84, annual volume bias <3%). The Markov Chain Monte Carlo approach constrains the parameters to which the model is most sensitive (e.g. lapse rate and recession coefficient) and maximizes model fit and performance. Model simulated streamflow using an interpolated precipitation data set decreases the fractional contribution from rainfall compared with simulations using observed station precipitation. The average snowmelt contribution to total runoff in the Tamor River basin for the 2002-2006 period is estimated to be 29.7+/-2.9% (which includes 4.2+/-0.9% from snowfall that promptly melts), whereas 70.3+/-2.6% is attributed to contributions from rainfall. On average, the elevation zone in the 4000-5500m range contributes the most to basin runoff, averaging 56.9+/-3.6% of all snowmelt input and 28.9+/-1.1% of all rainfall input to runoff. Model simulated streamflow using an interpolated precipitation data set decreases the fractional contribution from rainfall versus snowmelt compared with simulations using observed station precipitation. Model experiments indicate that the hydrograph itself does not constrain estimates of snowmelt versus rainfall contributions to total outflow but that this derives from the degree-day melting model. Lastly, we demonstrate that the data assimilation approach is useful for quantifying and reducing uncertainty related to model parameters and thus provides uncertainty bounds on snowmelt and rainfall contributions in such mountainous watersheds.

runoff↗

Exploitation of a Validation Hierarchy for Modeling and Simulation

Across engineering there is an evolving need to increase reliance on physics-based simulation to develop, design and optimize engineering systems. This increased reliance on modeling and simulation has highlighted a growing need to transform the confidence that modeling and simulation analysts have in their results into credibility for systems engineers to design and field systems more quickly and with less physical testing. For isolated components of a complex system, where a single discipline may drive product design, this is less of a concern as the relationship is often straightforward and easy to explain. However, when these isolated components are integrated, and are expected to operate in a multi-disciplinary context in which safety critical systems are involved, new concepts and model assurance standards are required. In this paper we address this challenge by showing how a model validation hierarchy can be exploited to identify those model validation experiments that will contribute most to increasing confidence and credibility of modeling and simulation predictions. The approach that is adopted contains four main steps. The first step is the construction of a model validation hierarchy that links subsystems, assemblies, and components to a hierarchy of physical experiments that can be used support model validation. This hierarchy connects the concerns of systems engineers to those of the modeling and simulation analyst in a clear and logical way. The structure and content of this hierarchy is then used in a second step to establish which physical phenomena have the greatest impact on overall system performance metrics. A gap analysis technique, based upon modeling and simulation concerns, is then used to prioritize the important physical phenomenon. Unfortunately, a common outcome of such gap analyses is the identification of many important gaps and so, in the final step of our process, we advocate the use of a global sensitivity analysis as a means to complete the prioritization.

Verification and Validation↗

Calculation of Dose Deposition in 3D Voxels by Heavy Ions and Simulation of gamma-H2AX Experiments

The biological response to high-LET radiation is different from low-LET radiation due to several factors, notably difference in energy deposition and formation of radiolytic species. Of particular importance in radiobiology is the formation of double-strand breaks (DSB), which can be detected by -H2AX foci experiments. These experiments has revealed important differences in the spatial distribution of DSB induced by low- and high-LET radiations [1,2]. To simulate -H2AX experiments, models based on amorphous track with radial dose are often combined with random walk chromosome models [3,4]. In this work, a new approach using the Monte-Carlo track structure code RITRACKS [5] and chromosome models have been used to simulate DSB formation. At first, RITRACKS have been used to simulate the irradiation of a cubic volume of 5 m by 1) 450 1H+ ions of 300 MeV (LET 0.3 keV/ m) and 2) by 1 56Fe26+ ion of 1 GeV/amu (LET 150 keV/ m). All energy deposition events are recorded to calculate dose in voxels of 20 m. The dose voxels are distributed randomly and scattered uniformly within the volume irradiated by low-LET radiation. Many differences are found in the spatial distribution of dose voxels for the 56Fe26+ ion. The track structure can be distinguished, and voxels with very high dose are found in the region corresponding to the track "core". These high-dose voxels are not found in the low-LET irradiation simulation and indicate clustered energy deposition, which may be responsible for complex DSB. In the second step, assuming that DSB will be found only in voxels where energy is deposited by the radiation, the intersection points between voxels with dose > 0 and simulated chromosomes were obtained. The spatial distribution of the intersection points is similar to -H2AX foci experiments. These preliminary results suggest that combining stochastic track structure and chromosome models could be a good approach to understand radiation-induced DSB and chromosome aberrations.

Plante, I.↗

Soil Moisture Data Assimilation in the NASA Land Information System for Local Modeling Applications and Improved Situational Awareness

As part of the NASA Soil Moisture Active Passive (SMAP) Early Adopter (EA) program, the NASA Shortterm Prediction Research and Transition (SPoRT) Center has implemented a data assimilation (DA) routine into the NASA Land Information System (LIS) for soil moisture retrievals from the European Space Agency's Soil Moisture Ocean Salinity (SMOS) satellite. The SMAP EA program promotes application‐driven research to provide a fundamental understanding of how SMAP data products will be used to improve decision‐making at operational agencies. SPoRT has partnered with select NOAA/NWS Weather Forecast Offices (WFOs) that use output from a real‐time regional configuration of LIS, without soil moisture DA, to initialize local numerical weather prediction (NWP) models and enhance situational awareness. Improvements to local NWP with the current LIS have been demonstrated; however, a better representation of the land surface through assimilation of SMOS (and eventually SMAP) retrievals is expected to lead to further model improvement, particularly during warm‐season months. SPoRT will collaborate with select WFOs to assess the impact of soil moisture DA on operational forecast situations. Assimilation of the legacy SMOS instrument data provides an opportunity to develop expertise in preparation for using SMAP data products shortly after the scheduled launch on 5 November 2014. SMOS contains a passive L‐band radiometer that is used to retrieve surface soil moisture at 35‐km resolution with an accuracy of 0.04 cu cm cm (exp -3). SMAP will feature a comparable passive L‐band instrument in conjunction with a 3‐km resolution active radar component of slightly degraded accuracy. A combined radar‐radiometer product will offer unprecedented global coverage of soil moisture at high spatial resolution (9 km) for hydrometeorological applications, balancing the resolution and accuracy of the active and passive instruments, respectively. The LIS software framework manages land surface model (LSM) simulations and includes an Ensemble Kalman Filter for conducting land surface DA. SPoRT has added a module to read, quality‐control and bias‐correct swaths of Level II SMOS soil moisture retrievals prior to assimilation within LIS. The impact of SMOS DA is being tested using the Noah LSM. Experiments are being conducted to examine the impacts of SMOS soil moisture DA on the resulting LISNoah fields and subsequent NWP simulations using the Weather Research and Forecasting (WRF) model initialized with LIS‐Noah output. LIS‐Noah soil moisture will be validated against in situ observations from Texas A&M's North American Soil Moisture Database to reveal the impact and possible improvement in soil moisture trends through DA. WRF model NWP case studies will test the impacts of DA on the simulated near‐surface and boundary‐layer environments, and precipitation during both quiescent and disturbed weather scenarios. Emphasis will be placed on cases with large analysis increments, especially due to contributions from regional irrigation patterns that are not represented by precipitation input in the baseline LIS‐Noah run. This poster presentation will describe the soil moisture DA methodology and highlight LIS‐Noah and WRF simulation results with and without assimilation.

Case, Jonathan L.↗

Computer simulation of on-orbit manned maneuvering unit operations

Simulation of spacecraft on-orbit operations is discussed in reference to Martin Marietta's Space Operations Simulation laboratory's use of computer software models to drive a six-degree-of-freedom moving base carriage and two target gimbal systems. In particular, key simulation issues and related computer software models associated with providing real-time, man-in-the-loop simulations of the Manned Maneuvering Unit (MMU) are addressed with special attention given to how effectively these models and motion systems simulate the MMU's actual on-orbit operations. The weightless effects of the space environment require the development of entirely new devices for locomotion. Since the access to space is very limited, it is necessary to design, build, and test these new devices within the physical constraints of earth using simulators. The simulation method that is discussed here is the technique of using computer software models to drive a Moving Base Carriage (MBC) that is capable of providing simultaneous six-degree-of-freedom motions. This method, utilized at Martin Marietta's Space Operations Simulation (SOS) laboratory, provides the ability to simulate the operation of manned spacecraft, provides the pilot with proper three-dimensional visual cues, and allows training of on-orbit operations. The purpose here is to discuss significant MMU simulation issues, the related models that were developed in response to these issues and how effectively these models simulate the MMU's actual on-orbiter operations.

Stuart, G. M.↗

Utilization of UARS Data in Validation of Photochemical and Dynamical Mechanisms in Stratospheric Models

The proposed work utilized Upper Atmosphere Research Satellite (UARS) measurements of short-lived and long-lived species, in conjunction with existing photochemical "box" models, trajectory models, and two-dimensional global models, to elucidate outstanding questions in our understanding of photochemical and dynamical mechanisms in the stratosphere. Particular emphasis was given to arriving at the best possible understanding of the chemical and dynamical contribution to the stratospheric ozone budget. Such understanding will increase confidence in the simulations carried out by assessment models.

Ko, Malcolm K. W.↗

3D Modeling of Ultrasonic Wave Interaction with Disbonds and Weak Bonds

Ultrasonic techniques, such as the use of guided waves, can be ideal for finding damage in the plate and pipe-like structures used in aerospace applications. However, the interaction of waves with real flaw types and geometries can lead to experimental signals that are difficult to interpret. 3-dimensional (3D) elastic wave simulations can be a powerful tool in understanding the complicated wave scattering involved in flaw detection and for optimizing experimental techniques. We have developed and implemented parallel 3D elastodynamic finite integration technique (3D EFIT) code to investigate Lamb wave scattering from realistic flaws. This paper discusses simulation results for an aluminum-aluminum diffusion disbond and an aluminum-epoxy disbond and compares results from the disbond case to the common artificial flaw type of a flat-bottom hole. The paper also discusses the potential for extending the 3D EFIT equations to incorporate physics-based weak bond models for simulating wave scattering from weak adhesive bonds.

Leckey, C.↗

Simulation of Sweep-Jet Flow Control, Single Jet and Full Vertical Tail

This work is a simulation technology demonstrator, of sweep jet flow control used to suppress boundary layer separation and increase the maximum achievable load coefficients. A sweep jet is a discrete Coanda jet that oscillates in the plane parallel to an aerodynamic surface. It injects mass and momentum in the approximate streamwise direction. It also generates turbulent eddies at the oscillation frequency, which are typically large relative to the scales of boundary layer turbulence, and which augment mixing across the boundary layer to attack flow separation. Simulations of a fluidic oscillator, the sweep jet emerging from a nozzle downstream of the oscillator, and an array of sweep jets which suppresses boundary layer separation are performed. Simulation results are compared to data from a dedicated validation experiment of a single oscillator and its sweep jet, and from a wind tunnel test of a full-scale Boeing 757 vertical tail augmented with an array of sweep jets. A critical step in the work is the development of realistic time-dependent sweep jet inflow boundary conditions, derived from the results of the single-oscillator simulations, which create the sweep jets in the full-tail simulations. Simulations were performed using the computational fluid dynamics (CFD) solver Overow, with high-order spatial discretization and a range of turbulence modeling. Good results were obtained for all flows simulated, when suitable turbulence modeling was used.

flow control↗

Multiscale Simulation of Deployable Composite Structures

In this paper, a multiscale simulation method for analyzing deployable composite structures is presented. Effective shell properties of the composites are obtained based on Mechanics of Structure Genome (MSG) homogenization, and then implemented into a user-subroutine UGENS for structural simulation with shell elements in Abaqus. The column bending test(CBT) of a flat thin flexure and lenticular composite boom in a simplified deployer structure are studied for demonstration. A viscoelastic material model with direct integration is adopted in this paper. The CBT simulation shows good agreement with experiments during relaxation, while errors are observed when comparing residual deformation. It is shown that this CBT model can be calibrated to CBT test results. For the lenticular boom analysis, the complete process of flattening, coiling, stowage, deployment and recovery is simulated with the viscoelastic shell model.

Finite element analysis↗

Magnetic Shielding of the Acceleration Channel Walls in a Long-Life Hall Thruster

In a Qualification Life Test (QLT) of the BPT-4000 Hall thruster that recently accumulated greater than 10,000 h it was found that the erosion of the acceleration channel practically stopped after approximately 5,600 h. Numerical simulations of this thruster using a 2-D axisymmetric, magnetic field-aligned-mesh (MFAM) plasma solver reveal that the process that led to this significant reduction of the erosion was multifaceted. It is found that when the channel receded from its early-in-life geometry to its steady-state configuration several changes in the near-wall plasma and sheath were induced by the magnetic field that, collectively, constituted an effective shielding of the walls from any significant ion bombardment. Because all such changes in the behavior of the ionized gas near the eroding surfaces were caused by the topology of the magnetic field there, we term this process "magnetic shielding."

life modeling↗

X-Ray Morphology of Cluster-Mass Haloes in Self-Interacting Dark Matter

We perform cosmological zoom-in simulations of 19 relaxed cluster-mass haloes with the inclusion of adiabatic gas in the cold dark matter (CDM) and self-interacting dark matter (SIDM) models. These clusters are selected as dynamically relaxed clusters from a parent simulation with M200(1-3) × 10 15 M. Both the dark matter and the intracluster gas distributions in SIDM appear more spherical than their CDM counterparts. Mock X-ray images are generated based on the simulations and are compared to the real X-ray images of 84 relaxed clusters selected from the Chandra and ROSAT archives. We perform ellipse fitting for the isophotes of mock and real X-ray images and obtain the ellipticities at cluster-centric radii of r0.1–0.2R 200 .The X-ray isophotes in SIDM models with increasing cross-sections are rounder than their CDM counterparts, which manifests as asystematic shift in the distribution function of ellipticities. Unexpectedly, the X-ray morphology of the observed non-cool-coreclusters agrees better with SIDM model swith cross-section (σ/m)=0.5–1cm2g−1 than CDM and SIDM with (σ/m)=0.1cm2g−1. Our statistical analysis indicates that the latter two models are disfavoured at the 68 percent confidence level (as conservative estimates). This conclusion is not altered by shifting the radial range of measurements or applying a temperature selection criterion. However, the primary uncertainty originates from the lack of baryonic physics in the adiabatic model, such as cooling, star formation and feedback effects, which still have the potential to reconcile CDM simulations with observations.

dark matter – cosmology↗