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 325 records · Page 18

RWPV bioreactor mass transport: earth-based and in microgravity

Mass transport and mixing of perfused scalar quantities in the NASA Rotating Wall Perfused Vessel bioreactor are studied using numerical models of the flow field and scalar concentration field. Operating conditions typical of both microgravity and ground-based cell cultures are studied to determine the expected vessel performance for both flight and ground-based control experiments. Results are presented for the transport of oxygen with cell densities and consumption rates typical of colon cancer cells cultured in the RWPV. The transport and mixing characteristics are first investigated with a step change in the perfusion inlet concentration by computing the time histories of the time to exceed 10% inlet concentration. The effects of a uniform cell utilization rate are then investigated with time histories of the outlet concentration, volume average concentration, and volume fraction starved. It is found that the operating conditions used in microgravity produce results that are quite different then those for ground-based conditions. Mixing times for microgravity conditions are significantly shorter than those for ground-based operation. Increasing the differential rotation rates (microgravity) increases the mixing and transport, while increasing the mean rotation rate (ground-based) suppresses both. Increasing perfusion rates enhances mass transport for both microgravity and ground-based cases, however, for the present range of operating conditions, above 5-10 cc/min there are diminishing returns as much of the inlet fluid is transported directly to the perfusion exit. The results show that exit concentration is not a good indicator of the concentration distributions in the vessel. In microgravity conditions, the NASA RWPV bioreactor with the viscous pump has been shown to provide an environment that is well mixed. Even when operated near the theoretical minimum perfusion rates, only a small fraction of the volume provides less than the required oxygen levels. 2002 Wiley Periodicals, Inc.

Models, Biological↗

ISS Destiny Laboratory Smoke Detection Model

Smoke transport and detection were modeled numerically in the ISS Destiny module using the NIST, Fire Dynamics Simulator code. The airflows in Destiny were modeled using the existing flow conditions and the module geometry included obstructions that simulate the currently installed hardware on orbit. The smoke source was modeled as a 0.152 by 0.152 m region that emitted smoke particulate ranging from 1.46 to 8.47 mg/s. In the module domain, the smoke source was placed in the center of each Destiny rack location and the model was run to determine the time required for the two smoke detectors to alarm. Overall the detection times were dominated by the circumferential flow, the axial flow from the intermodule ventilation and the smoke source strength.

Brooker, John E.↗

Multiscale Modeling of Damage Processes in Aluminum Alloys: Grain-Scale Mechanisms

This paper has two goals related to the development of a physically-grounded methodology for modeling the initial stages of fatigue crack growth in an aluminum alloy. The aluminum alloy, AA 7075-T651, is susceptible to fatigue cracking that nucleates from cracked second phase iron-bearing particles. Thus, the first goal of the paper is to validate an existing framework for the prediction of the conditions under which the particles crack. The observed statistics of particle cracking (defined as incubation for this alloy) must be accurately predicted to simulate the stochastic nature of microstructurally small fatigue crack (MSFC) formation. Also, only by simulating incubation of damage in a statistically accurate manner can subsequent stages of crack growth be accurately predicted. To maintain fidelity and computational efficiency, a filtering procedure was developed to eliminate particles that were unlikely to crack. The particle filter considers the distributions of particle sizes and shapes, grain texture, and the configuration of the surrounding grains. This filter helps substantially reduce the number of particles that need to be included in the microstructural models and forms the basis of the future work on the subsequent stages of MSFC, crack nucleation and microstructurally small crack propagation. A physics-based approach to simulating fracture should ultimately begin at nanometer length scale, in which atomistic simulation is used to predict the fundamental damage mechanisms of MSFC. These mechanisms include dislocation formation and interaction, interstitial void formation, and atomic diffusion. However, atomistic simulations quickly become computationally intractable as the system size increases, especially when directly linking to the already large microstructural models. Therefore, the second goal of this paper is to propose a method that will incorporate atomistic simulation and small-scale experimental characterization into the existing multiscale framework. At the microscale, the nanoscale mechanics are represented within cohesive zones where appropriate, i.e. where the mechanics observed at the nanoscale can be represented as occurring on a plane such as at grain boundaries or slip planes at a crack front. Important advancements that are yet to be made include: 1. an increased fidelity in cohesive zone modeling; 2. a means to understand how atomistic simulation scales with time; 3. a new experimental methodology for generating empirical models for CZMs and emerging materials; and 4. a validation of simulations of the damage processes at the nano-micro scale. With ever-increasing computer power, the long-term ability to employ atomistic simulation for the prognosis of structural components will not be limited by computation power, but by our lack of knowledge in incorporating atomistic models into simulations of MSFC into a multiscale framework.

Hochhalter, J. D.↗

Base Flow and Heat Transfer Characteristics of a Four-Nozzle Clustered Rocket Engine: Effect of Nozzle Pressure Ratio

The base pressure and heating characteristics of a four-nozzle clustered rocket configuration is studied numerically with the aid of OVERFLOW Navier-Stokes code. A pressure ratio (chamber pressure to freestream static pressure) range of 990 to 5,920 and a freestream Mach number range of 2.5 to 3.5 are studied. The qualitative trends of decreasing base pressure with increasing pressure ratio and increasing base heat flux with increasing pressure ratio are correctly predicted. However, the predictions for base pressure and base heat flux show deviations from the wind tunnel data. The differences in absolute values between the computation and the data are attributed to factors such as perfect gas (thermally and calorically perfect) assumption, turbulence model inaccuracies in the simulation, and lack of grid adaptation.

Nallasamy, R.↗

Changes in Extratropical Storm Track Cloudiness 1983-2008: Observational Support for a Poleward Shift

Climate model simulations suggest that the extratropical storm tracks will shift poleward as a consequence of global warming. In this study the northern and southern hemisphere storm tracks over the Pacific and Atlantic ocean basins are studied using observational data, primarily from the International Satellite Cloud Climatology Project, ISCCP. Potential shifts in the storm tracks are examined using the observed cloud structures as proxies for cyclone activity. Different data analysis methods are employed, with the objective to address difficulties and uncertainties in using ISCCP data for regional trend analysis. In particular, three data filtering techniques are explored; excluding specific problematic regions from the analysis, regressing out a spurious viewing geometry effect, and excluding specific cloud types from the analysis. These adjustments all, to varying degree, moderate the cloud trends in the original data but leave the qualitative aspects of those trends largely unaffected. Therefore, our analysis suggests that ISCCP data can be used to interpret regional trends in cloudiness, provided that data and instrumental artefacts are recognized and accounted for. The variation in magnitude between trends emerging from application of different data correction methods, allows us to estimate possible ranges for the observational changes. It is found that the storm tracks, here represented by the extent of the midlatitude-centered band of maximum cloud cover over the studied ocean basins, experience a poleward shift as well as a narrowing over the 25 year period covered by ISCCP. The observed magnitudes of these effects are larger than in current generation climate models (CMIP3). The magnitude of the shift is particularly large in the northern hemisphere Atlantic. This is also the one of the four regions in which imperfect data primarily prevents us from drawing firm conclusions. The shifted path and reduced extent of the storm track cloudiness is accompanied by a regional reduction in total cloud cover. This decrease in cloudiness can primarily be ascribed to low level clouds, whereas the upper level cloud fraction actually increases, according to ISCCP. Independent satellite observations of radiative fluxes at the top of the atmosphere are consistent with the changes in total cloud cover. The shift in cloudiness is also supported by a shift in central position of the mid-troposphere meridional temperature gradient. We do not find support for aerosols playing a significant role in the satellite observed changes in cloudiness. The observed changes in storm track cloudiness can be related to local cloud-induced changes in radiative forcing, using ERBE and CERES radiative fluxes. The shortwave and the longwave components are found to act together, leading to a positive (warming) net radiative effect in response to the cloud changes in the storm track regions, indicative of positive cloud feedback. Among the CMIP3 models that simulate poleward shifts in all four storm track areas, all but one show decreasing cloud amount on a global mean scale in response to increased CO2 forcing, further consistent with positive cloud feedback. Models with low equilibrium climate sensitivity to a lesser extent than higher-sensitivity models simulate a poleward shift of the storm tracks.

Bender, Frida A-M.↗

Simulated Phase-dependent Spectra of Terrestrial Aquaplanets in M Dwarf Systems

Orbital phase-dependent variations in thermal emission and reflected stellar energy spectra can provide meaningful constraints on the climate states of terrestrial extrasolar planets orbiting M dwarf stars. Spatial distributions of water vapor, clouds, and surface ice are controlled by climate. In turn, water, in each of its thermodynamic phases, imposes significant modulations to thermal and reflected planetary spectra. Here we explore these characteristic spectral signals, based on 3D climate simulations of Earth-sized aquaplanets orbiting M dwarf stars near the habitable zone. By using 3D models, we can self-consistently predict surface temperatures and the location ofwater vapor, clouds, and surface ice in the climate system. Habitable zone planets in M dwarf systems are expected to be in synchronous rotation with their host star and thus present distinct differences in emitted and reflected energy fluxes depending on the observed hemisphere. Here we illustrate that icy, temperate, and incipient runaway greenhouse climate states exhibit phase-dependent spectral signals that enable their characterization.

Wolf, E. T.↗

Wall-Modeled Large-Eddy Simulation of Supersonic Parachute Inflation

Supersonic parachutes have been used in nearly every robotic mission to another planetary body with an atmosphere because they are one of the most mass-efficient ways to decelerate a payload to land on the surface. Short of performing flight tests in the upper Earth atmosphere, we currently cannot reliably predict a novel parachute system’s performance or potential failure modes, or even confidently explain it after the fact, as was the case with the Low Density Supersonic Decelerator flight tests and subsequent investigations. Fluid-structure interaction (FSI) simulations have the potential to bridge this gap. To this end, we present several improvements to the state-of-the-art for simulating supersonic parachutes using FSI: a higher effective resolution convective flux, an immersed boundary turbulent wall layer modeling approach to capture viscous effects, and a novel method to obtain a more realistic initial parachute shape. A recent supersonic parachute flight test is simulated and compared to measurements for the purposes of model validation.

ESM↗

Electrification Within Wintertime Stratiform Regions Sampled During the 2020/2022 NASA IMPACTS Field Campaign

Two nor'easter events—sampled during the NASA Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) field campaign—were examined to characterize the microphysics in relation to the underlying electrification processes within wintertime stratiform regions. A theoretical model was developed to determine whether accretion or diffusion growth regimes were preferential during periods of greatest electrification. Model simulation with electrification parameterization was used to provide supplemental context to the physical processes of in-cloud microphysics and electrification. The strongest electric fields (i.e., ∼80 V m −1 at 20 km) during the 2020 NASA IMPACTS deployment was associated with large non-rimed ice crystals colliding with each other. During the 29–30 January 2022 science flight, the NASA P-3 microphysical probe data demonstrated that non-inductive charging was possible off the coastline of Cape Cod, Massachusetts. Later in the science flight, when the NASA P-3 and ER-2 were coordinating with each other, measured electric fields consistently were less than 8 V m −1 and electrification was subdued owing to reduced concentrations of graupel and large ice hydrometeors. Altogether, the in-situ observations provide evidence for the non-riming collisional charging mechanism and demonstrates that graupel and supercooled liquid water may not be necessary for weak electrification within wintertime stratiform regions. Model output from simulation of both events suggested that the main synoptic snowbands were associated with elevated hydrometeor snow charge density and electric fields.

snowfall↗

Computer simulation of fibrillation threshold measurements and electrophysiologic testing procedures

A finite element model of cardiac conduction was used to simulate two experimental protocols: 1) fibrillation threshold measurements and 2) clinical electrophysiologic (EP) testing procedures. The model consisted of a cylindrical lattice whose properties were determined by four parameters: element length, conduction velocity, mean refractory period, and standard deviation of refractory periods. Different stimulation patterns were applied to the lattice under a given set of lattice parameter values and the response of the model was observed through a simulated electrocardiogram. The studies confirm that the model can account for observations made in experimental fibrillation threshold measurements and in clinical EP testing protocols.

NASA Discipline Cardiopulmonary↗

Biomass Estimation from Simulated GEDI, ICESat-2 and NISAR Across Environmental Gradients in Sonoma County, California

Estimates of the magnitude and distribution of aboveground carbon in Earth’s forests remain uncertain, yet knowledge of forest carbon content at a global scale is critical for forest management in support of climate mitigation. In light of this knowledge gap, several upcoming spaceborne missions aim to map forest aboveground biomass, and many new biomass products are expected from these datasets. As these new missions host different technologies, each with relative strengths and weaknesses for biomass retrieval, as well as different spatial resolutions, consistently comparing or combining biomass estimates from these new datasets will be challenging. This paper presents a demonstration of an inter-comparison of biomass estimates from simulations of three NASA missions (GEDI, ICESat-2 and NISAR) over Sonoma county in California, USA. We use a high resolution, locally calibrated airborne lidar map as our reference dataset, and emphasize the importance of considering uncertainties in both reference maps and spaceborne estimates when conducting biomass product validation. GEDI and ICESat-2 were simulated from airborne lidar point clouds, while UAVSAR’s L-band backscatter was used as a proxy for NISAR. To estimate biomass for the lidar missions we used GEDI’s footprint-level biomass algorithms, and also adapted these for application to ICESat-2. For UAVSAR, we developed a locally trained biomass model, calibrated against the ALS reference map. Each mission simulation was evaluated in comparison to the local reference map at its native product resolution (25 m, 100m transect, and 1 ha) yielding RMSEs of 57%, 75%, and 89% for GEDI, NISAR, and ICESat-2 respectively. RMSE values increased for GEDI’s power beam during simulated daytime conditions (64%), coverage beam during nighttime conditions (72%), and coverage beam daytime conditions (87%). We also test the application of GEDI’s biomass modeling framework for estimation of biomass from ICESat-2, and fine that ICESat-2 yields reasonable biomass estimates, particularly in relatively short, open canopies. Results suggest that while all three missions will produce datasets useful for biomass mapping, tall, dense canopies such as those found in Sonoma County present the greatest challenges for all three missions, while steep slopes also prove challenging for single-date SAR based biomass retrieval. Our methods provide guidance for the inter-comparison and validation of spaceborne biomass estimates through the use of airborne lidar reference maps, and could be repeated with on-orbit estimates in any area with high quality field plot and ALS data. These methods allow for regional interpretations and filtering of multi-mission biomass estimates toward improved wall-to-wall biomass maps through data fusion.

GEDI↗

Changes in Anthropogenic Precursor Emissions Drive Shifts in the Ozone Seasonal Cycle Throughout the Northern Midlatitude Troposphere

Simulations by six Coupled Model Intercomparison Project Phase 6 (CMIP6) Earth system models indicate that the seasonal cycle of baseline tropospheric ozone at northern midlatitudes has been shifting since the mid-20th century. Beginning in ∼ 1940, the magnitude of the seasonal cycle increased by ∼10 ppb (measured from seasonal minimum to maximum), and the seasonal maximum shifted to later in the year by about 3 weeks. This shift maximized in the mid-1980s, followed by a reversal – the seasonal cycle decreased in amplitude and the maximum shifted back to earlier in the year. Similar changes are seen in measurements collected from the 1970s to the present. The timing of the seasonal cycle changes is generally concurrent with the rise and fall of anthropogenic emissions that followed industrialization and the subsequent implementation of air quality emission controls. A quantitative comparison of the temporal changes in the ozone seasonal cycle at sites in both Europe and North America with the temporal changes in ozone precursor emissions across the northern midlatitudes found a high degree of similarity between these two temporal patterns. We hypothesize that changing precursor emissions are responsible for the shift in the ozone seasonal cycle; this is supported by the absence of such seasonal shifts in southern midlatitudes where anthropogenic emissions are much smaller. We also suggest a mechanism by which changing emissions drive the changing seasonal cycle: increasing emissions of NOx allow summertime photochemical production of ozone to become more important than ozone transported from the stratosphere, and increasing volatile organic compounds (VOCs) lead to progressively greater photochemical ozone production in the summer months, thereby increasing the amplitude of the seasonal ozone cycle. Decreasing emissions of both precursor classes then reverse these changes. The quantitative parameter values that characterize the seasonal shifts provide useful benchmarks for evaluating model simulations, both against observations and between models.

Ozone↗

The Challenge of Configuring Model-Based Space Mission Planners

Mission planning is central to space mission operations, and has benefited from advances in model-based planning software. Constraints arise from many sources, including simulators and engineering specification documents, and ensuring that constraints are correctly represented in the planner is a challenge. As mission constraints evolve, planning domain modelers need help with modeling constraints efficiently using the available source data, catching errors quickly, and correcting the model. This paper describes the current state of the practice in designing model-based mission planning tools, the challenges facing model developers, and a proposed Interactive Model Development Environment (IMDE) to configure mission planning systems. We describe current and future technology developments that can be integrated into an IMDE.

model-based systems↗

Advanced Supersonic Parachute Inflation Research and Experiment-2 (ASPIRE2) Flight Mechanics Modeling and Simulation

Introduction: The Advanced Supersonic Para-chute Inflation Research Experiment-2 (ASPIRE2) program is a sounding rocket flight test to be conducted at Wallops Island, VA in early 2025. This program is a risk mitigation exercise for the Mars Sample Retrieval Lander (MSRL) program and builds upon the success of the 2017 & 2018 ASPIRE program, [1,2]. While the ASPIRE program certified the strengthened 21.5 m diameter disk-gap-band (DGB) deployed at Mach 1.7, ASPIRE2 will certify a 24 m diameter DGB deployed at Mach 2.1. The need to certify this increased parachute performance is driven the increased lander mass for MSRL; over 50% increase in comparison to M2020. Modeling: Like its predecessor, ASPIRE2 is developing a multi-body flight dynamics simulation to predict parachute dynamics and aide in designing the flight test that will target Mars-relevant flight conditions, as shown in Figure 1. This work de-scribes the parachute modeling, flight mechanics simulations (from payload separation to splash-down), and design trades used to prepare for the 2025 ASPIRE2 flight. Discussed herein are comparisons between ASPIRE and ASPIRE2, noting key differences in the parachute modeling and vehicle configuration. A study on the attitude control system performance impacts with respect to payload section design will also be presented. This work will quantify the pre-flight parachute performance in the presence of uncertainties, such as those associated with the separation from the sounding rocket, atmosphere, the parachute system, and vehicle mass. The pre-flight predictions will include Monte Carlo analyses, powered by the flight mechanics simulations to show the ASPIRE2 vehicle performance in meeting program requirements on parachute deployment conditions (Mach & dynamic pressure), parachute loads, vehicle attitude at key milestones, and the vehicle splashdown conditions.

Entry Descent and Landing↗

Effect of gaseous and solid simulated jet plumes on an 040A space shuttle launch configuration at m=1.6 to 2.2

The effect of plume-induced flow separation and aspiration effects due to operation of both orbiter and the solid rocket motors on a 0.019-scale model of the launch configuration of the Space Shuttle Vehicle is determined. Longitudinal and lateral-directional stability data were obtained at Mach numbers of 1.6, 2.0, and 2.2 with and without the engines operating. The plumes exiting from the engines were simulated by a cold-gas jet supplied by an auxiliary 200-atm air supply system and solid-body plume simulators. The aerodynamic effects produced by these two simulation procedures are compared. The parameters most significantly affected by the jet plumes are pitching moment, elevon control effectiveness, axial force, and orbiter wing loads. The solid rocket motor (SRM) plumes have the largest effect on the aerodynamic characteristics. The effect of the orbiter plumes in combination with the SRM plumes is also significant. Variations in the nozzle design parameters and configuration changes can reduce the jet plume-induced aerodynamic effects.

Dods, J. B., Jr.↗

The PMIP4 Last Glacial Maximum experiments: preliminary results and comparison with the PMIP3 simulations

The Last Glacial Maximum (LGM, ∼ 21 000 years ago) has been a major focus for evaluating how well state-of-the-art climate models simulate climate changes as large as those expected in the future using paleoclimate reconstructions. A new generation of climate models has been used to generate LGM simulations as part of the Paleoclimate Modelling Intercomparison Project (PMIP) contribution to the Coupled Model Intercomparison Project (CMIP). Here, we provide a preliminary analysis and evaluation of the results of these LGM experiments (PMIP4, most of which are PMIP4-CMIP6) and compare them with the previous generation of simulations (PMIP3, most of which are PMIP3-CMIP5). We show that the global averages of the PMIP4 simulations span a larger range in terms of mean annual surface air temperature and mean annual precipitation compared to the PMIP3-CMIP5 simulations, with some PMIP4 simulations reaching a globally colder and drier state. However, the multi-model global cooling average is similar for the PMIP4 and PMIP3 ensembles, while the multi-model PMIP4 mean annual precipitation average is drier than the PMIP3 one. There are important differences in both atmospheric and oceanic circulations between the two sets of experiments, with the northern and southern jet streams being more poleward and the changes in the Atlantic Meridional Overturning Circulation being less pronounced in the PMIP4-CMIP6 simulations than in the PMIP3-CMIP5 simulations. Changes in simulated precipitation patterns are influenced by both temperature and circulation changes. Differences in simulated climate between individual models remain large. Therefore, although there are differences in the average behaviour across the two ensembles, the new simulation results are not fundamentally different from the PMIP3-CMIP5 results. Evaluation of large-scale climate features, such as land–sea contrast and polar amplification, confirms that the models capture these well and within the uncertainty of the paleoclimate reconstructions. Nevertheless, regional climate changes are less well simulated: the models underestimate extratropical cooling, particularly in winter, and precipitation changes. These results point to the utility of using paleoclimate simulations to understand the mechanisms of climate change and evaluate model performance.

Last Glacial Maximum↗

Synthesis of a control model for a liquid nitrogen cooled, closed circuit, cryogenic nitrogen wind tunnel and its validation

The details of the efforts to synthesize a control-compatible multivariable model of a liquid nitrogen cooled, gaseous nitrogen operated, closed circuit, cryogenic pressure tunnel are presented. The synthesized model was transformed into a real-time cryogenic tunnel simulator, and this model is validated by comparing the model responses to the actual tunnel responses of the 0.3 m transonic cryogenic tunnel, using the quasi-steady-state and the transient responses of the model and the tunnel. The global nature of the simple, explicit, lumped multivariable model of a closed circuit cryogenic tunnel is demonstrated.

Balakrishna, S.↗

Transonic flutter and gust-response tests and analyses of a wind-tunnel model of a torsion free wing airplane

An exploratory study of a 1/5.5 size, complete airplane version of a torsion free wing (TFW) fighter aircraft was conducted. The TFW consisted of a wing/boom/canard assembly on each fuselage side that was interconnected by a common pivot shaft so that the TFW could rotate freely in pitch. The effect of the TFW was evaluated by comparing data obtained with the TFW free and the TFW locked to the fuselage. With the model mounted on cables to simulate an airplane free flying condition, flutter boundaries were measured at Mach number (M) from 0.85 to 1.0 and gust responses at M = 0.65 and 0.90. The critical flutter mode for the TFW free configuration was found experimentally to occur at M = 0.95 and had the rigid TFW pitch mode as its apparent aerodynamic driver.

Ruhlin, C. L.↗