Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “modeling simulation”

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 1,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Molten salt reactor neutronics and fuel cycle modeling and simulation with SCALE

Current interest in advanced nuclear energy and molten salt reactor (MSR) concepts has enhanced interest in building the tools necessary to analyze these systems. A Python script known as ChemTriton has been developed to simulate equilibrium MSR fuel cycle performance by modeling the changing isotopic composition of an irradiated fuel salt using SCALE for neutron transport and depletion calculations. Some capabilities in ChemTriton that have improved, include a generic geometry capable of modeling multi-zone and multi-fluid systems, enhanced time-dependent feed and separations, and a critical concentration search. Although more generally applicable, the capabilities developed to date are illustrated in this paper in three applied problems: (1) simulating the startup of a thorium-based MSR fuel cycle (a likely scenario requires the first of these MSRs to be started without available 233U); (2) determining the effect of the removal of different fission products on MSR operations; and (3) obtaining the equilibrium concentration of a mixed-oxide light-water reactor fuel in a two-stage fuel cycle with a sodium fast reactor. Moreover, the third problem is chosen to demonstrate versatility in an application to analyze the fuel cycle of a non-MSR system. During the first application, the initial fuel salt compositions fueled with different sources of fissile material are made feasible after (1) removing the associated nonfissile actinides after much of the initial fissile isotopes have burned and (2) optimizing the thorium concentration to maintain a critical configuration without significantly reducing breeding capability. In the second application, noble metal, volatile gas, and rare earth element fission products are shown to have a strong negative effect on criticality in a uranium-fueled thermal-spectrum MSR; their removal significantly increases core lifetime (by 30%) and fuel utilization. In the third application, the fuel of a mixed-oxide light-water reactor approaches an equilibrium composition after 20 depletion steps, demonstrating the potential for the longer time scales required to achieve equilibrium for solid-fueled systems over liquid fuel systems. This time to equilibrium can be reduced by starting with an initial fuel composition closer to that of the equilibrium fuel, reducing the need to handle time-dependent fuel compositions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Radar imaging of shallow water bathymetry: A case study in the Yangtze Estuary

This study focuses on 2-dimensional (2-D) radar imaging of bathymetric features in the shallow water of the Yangtze Estuary using synthetic aperture radar (SAR) observations and model simulations. A validated 2-D shallow water numerical model simulates the barotropic current velocity, and the simulated current fields together with the relevant parameters of radar observations are then invoked in the radar imaging model as the input. The results show that variations in the simulated image intensity are mainly dominated by distinct radar backscatter anomalies caused by wave-current interactions in the vicinity of rapidly changing underwater topographies. The comparison between the simulated and observed SAR images shows a reasonable agreement, demonstrating that our approach may be implemented to monitor changes in the shallow water bathymetry of the Yangtze Estuary in the future.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2011

Engine Icing Modeling and Simulation (Part 2): Performance Simulation of Engine Rollback Phenomena

Ice buildup in the compressor section of a commercial aircraft gas turbine engine can cause a number of engine failures. One of these failure modes is known as engine rollback: an uncommanded decrease in thrust accompanied by a decrease in fan speed and an increase in turbine temperature. This paper describes the development of a model which simulates the system level impact of engine icing using the Commercial Modular Aero-Propulsion System Simulation 40k (C-MAPSS40k). When an ice blockage is added to C-MAPSS40k, the control system responds in a manner similar to that of an actual engine, and, in cases with severe blockage, an engine rollback is observed. Using this capability to simulate engine rollback, a proof-of-concept detection scheme is developed and tested using only typical engine sensors. This paper concludes that the engine control system s limit protection is the proximate cause of iced engine rollback and that the controller can detect the buildup of ice particles in the compressor section. This work serves as a feasibility study for continued research into the detection and mitigation of engine rollback using the propulsion control system.

May, Ryan D.↗
NASA NTRS2014

How Do Various Maize Crop Models Vary in Their Responses to Climate Change Factors?

Potential consequences of climate change on crop production can be studied using mechanistic crop simulation models. While a broad variety of maize simulation models exist, it is not known whether different models diverge on grain yield responses to changes in climatic factors, or whether they agree in their general trends related to phenology, growth, and yield. With the goal of analyzing the sensitivity of simulated yields to changes in temperature and atmospheric carbon dioxide concentrations [CO2], we present the largest maize crop model intercomparison to date, including 23 different models. These models were evaluated for four locations representing a wide range of maize production conditions in the world: Lusignan (France), Ames (USA), Rio Verde (Brazil) and Morogoro (Tanzania). While individual models differed considerably in absolute yield simulation at the four sites, an ensemble of a minimum number of models was able to simulate absolute yields accurately at the four sites even with low data for calibration, thus suggesting that using an ensemble of models has merit. Temperature increase had strong negative influence on modeled yield response of roughly -0.5 Mg ha(sup 1) per degC. Doubling [CO2] from 360 to 720 lmol mol 1 increased grain yield by 7.5% on average across models and the sites. That would therefore make temperature the main factor altering maize yields at the end of this century. Furthermore, there was a large uncertainty in the yield response to [CO2] among models. Model responses to temperature and [CO2] did not differ whether models were simulated with low calibration information or, simulated with high level of calibration information.

carbon dioxide↗
DOE OSTI2018

Soil moisture drives microbial controls on carbon decomposition in two subtropical forests

Knowledge of microbial mechanisms is critical to understand Earth's biogeochemical cycle under climate and environmental changes. However, large uncertainties remain in model simulations and predictions due to the lack of explicit parameterization of microbial data and few applications beyond the laboratory. In addition, most experimental and modeling studies of warming-induced changes in soil carbon (C) focus on temperature sensitivity, neglecting concomitant effects of changes in soil moisture. Soil microbes are sensitive to moisture, and their responses can dramatically impact soil biogeochemical cycles. Here we represent microbial and enzymatic functions in response to changes in moisture in the Microbial-ENzyme Decomposition (MEND) model. Through modeling with long-term field observations from subtropical forests, we demonstrate that parameterization with microbial data in addition to respiration fluxes greatly increases confidence in model simulations. We further employ the calibrated model to simulate the responses of soil organic C (SOC) under multiple environmental change scenarios. The model shows significant increases in SOC in response to decreasing soil moisture and only minor changes in SOC in response to increasing soil temperature. Increasing litter inputs also cause a significant increase in SOC in the pine forest, whereas an insignificant negative effect is simulated in the broadleaf forest. We also demonstrate the co-metabolism mechanism for the priming effects, i.e., more labile inputs to soil could stimulate microbial and enzymatic growth and activity. Furthermore our study provides strong evidence of microbial control over soil C decomposition and suggests the future trajectory of soil C may be more responsive to changes in soil moisture than temperature, particularly in tropical and subtropical environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Towards a theory of economic value for modeling and simulation: incremental cost of parallel simulation (wip)

A useful perspective on the economic value of simulation in engineering comes from considering the role simulation plays in constructing a system. We propose a conceptual model that combines an embryonic concept of the value of an information bit with cost estimation tools used in software engineering. The model provides an economic justification for the common conception that models should be built to address specific questions rather than to comprehensively represent a referent system. Extending these ideas, we offer a model for the incremental value of parallel execution to explain why simulation applications that take advantage of sophisticated parallel discrete event simulation tools form a relatively narrow niche despite the ubiquity of parallel computing hardware. The proposed model enables us to bring value of information criteria into model development and to make predictions about the future of parallel simulation applications.

Nutaro, James J↗
DOE OSTI2020

Hybrid predictive modeling approach simulated evapotranspiration and ecosystem respiration data

This dataset includes measured data used for developing hybrid-predictive-modeling (HPM) approach and simulated evapotranspiration and ecosystem respiration data across several Fluxnet sites, SNOTEL sites and East River locations (Chen et al., 2020 in review). Fluxnet sites considered in this study include: CA-OAS, CA-OBS, US-NR1, US-TON, US-VAR, US-SRM, US-WHS, US-WKG. Snotel sites Butte (ER-BT), Porphyry Creek (ER-PK) and Schofield Pass (ER-SP) are included as well 16 east river locations with different vegetation types denoted by evergreen forests (EF), deciduous forests (DF), riparian shrublands (RS) and meadow grassland (MS). This data package includes: 1) A summary spread sheet of the sites and acronyms for data variables. 2) Processed Fluxnet datasets that are obtained from Fluxnet.fluxdata.org with HPM estimated evapotranspiration (ET) and ecosystem respiration (RECO). Trained models and metrics for parameters for the Fluxnet sites were added in an update on November 2, 2020; 3) Processed SNOTEL datasets obtained from www.nrcs.usda.gov/snow/ and HPM estimated ET and RECO; 4) Weather data at east river locations obtained from DAYMET (daymet.ornl.gov) as well as HPM estimated ET and RECO. NDVI datasets are obtained from Landsat 5, 7 and 8 on Google Earth Engine platform. Detailed QA/QC and sampling methods of measured data are available at the original data source.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1992

Preliminary subsonic aerodynamic model for simulation studies of the HL-20 lifting body

A nonlinear, six-degree-of-freedom aerodynamic model for an early version of the HL-20 lifting body is described and compared with wind tunnel data upon which it is based. Polynomial functions describing most of the aerodynamic parameters are given and tables of these functions are presented. Techniques used to arrive at these functions are described. Basic aerodynamic coefficients were modeled as functions of angles of attack and sideslip. Vehicle lateral symmetry was assumed. Compressibility (Mach) effects were ignored. Control-surface effectiveness was assumed to vary linearly with angle of deflection and was assumed to be invariant with the angle of sideslip. Dynamic derivatives were obtained from predictive aerodynamic codes. Landing-gear and ground effects were scaled from Space Shuttle data. The model described is provided to support pilot-in-the-loop simulation studies of the HL-20. By providing the data in tabular format, the model is suitable for the data interpolation architecture of many existing engineering simulation facilities. Because of the preliminary nature of the data, however, this model is not recommended for study of the absolute performance of the HL-20.

Jackson, E. Bruce↗
DOE OSTI2019

A simulation-based model for understanding the time dependent x-ray drive asymmetries and error bars in indirectly driven implosions on the National Ignition Facility

Time-dependent low-mode asymmetries are believed to play a leading role in limiting the performance of current inertial confinement fusion implosions on the National Ignition Facility (NIF) [E. I. Moses et al., Phys. Plasmas 16, 041006 (2009)]. These long wavelength modes are initiated and driven by asymmetries in the x-ray flux from the hohlraum; however, the underlying hydrodynamics of the implosion also act to modify and amplify these asymmetries. Here we present here a simulation-based model connecting the time-dependent drive asymmetry seen by the capsule to the measured inflight and hot spot symmetries. This approach is based on a Green's function analysis for which we evaluate the response of the capsule to impulses of drive asymmetry at a series of times. Our model sheds new light on the sensitivity to the drive asymmetry of an imploded capsule, giving a new tool for design. Inverting the problem and finding the drive asymmetry needed to match the experimental data allow us to tightly constrain the drive asymmetry seen by the capsule, providing an error estimate on the result. Doing so, we are able to point out when and how the complex hohlraum simulations start to deviate from what they should obtain to match the experimental data. Ultimately, we project to use this model to make some experimental recommendations to fix the time-dependent low-mode asymmetry of indirectly driven implosions and identify additional measurements to further constrain the asymmetries with a view to improving target design on the NIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Head-Down Tilt with Balanced Traction as a Model for Simulating Spinal Acclimation to Microgravity

Astronauts experience total body height increases of 4 to 7 cm in microgravity. Thus, stretching of the spinal cord, nerve roots, and muscular and ligamentous tissues may be responsible for the hyperreflexia, back pain, and muscular atrophy associated with exposure to microgravity. Axial compression of the spine makes 6 deg. head-down tilt (HDT) an unsuitable model for spinal acclimation to microgravity. However, this axial compression may be counteracted by balanced traction consisting of 10% body weight (sin 6 deg. = 0.1) applied to the legs. Six healthy male subjects underwent 3 days each of 60 HDT with balanced traction and horizontal bed rest (HBR), with a 2 week recovery period between treatments. Total body and spine length, lumbar disc height, back pain, erector spinae intramuscular pressure, and ankle joint torque were measured before, during and after each treatment. Total body and spine (processes of L5 - C7) lengths increased significantly more during HDT with balanced traction (22 +/- 8 mm and 25 +/- 8 mm, respectively) than during HBR (16 +/- 4 mm and 14 +/- 9 mm, respectively). Back and leg pain were significantly greater during HDT with balanced traction than during HBR. The distance between the lower end plate of L4 and the upper endplate of S1, as measured by sonography, increased significantly in both treatments to the same degree (2.9 +/- 1.9 mm, HDT with balanced traction; 3.3 +/- 1.5 mm, HBR). Intramuscular pressure of the erector spinae muscles and maximal ankle joint torque were unaltered with both models. While neither model increased height to the magnitude observed in microgravity, HDT with balanced traction may be a better model for simulating the body lengthening and back pain experienced in microgravity.

Ballard, R. E.↗
DOE OSTI2022

The contribution of Humboldt Glacier, northern Greenland, to sea-level rise through 2100 constrained by recent observations of speedup and retreat

Abstract. Humboldt Glacier, northern Greenland, has retreated and accelerated through the 21st century, raising concerns that it could be a significant contributor to future sea-level rise. We use a data-constrained ensemble of three-dimensional higher-order ice sheet model simulations to estimate the likely range of sea-level rise from the continued retreat of Humboldt Glacier. We first solve for basal traction using observed ice thickness, bed topography, and ice surface velocity from the year 2007 in a PDE-constrained (partial differential equation) optimization. Next, we impose calving rates to match mean observed retreat rates from winter 2007–2008 to winter 2017–2018 in a transient calibration of the exponent in the power-law basal friction relationship. We find that power-law exponents in the range of 1/7–1/5 – rather than the commonly used 1/3–1 – are necessary to reproduce the observed speedup over this period. We then tune an iceberg calving parameterization based on the von Mises stress yield criterion in another transient-calibration step to approximate both observed ice velocities and terminus position in 2017–2018. Finally, we use the range of basal friction relationship exponents and calving parameter values to generate the ensemble of model simulations from 2007–2100 under three climate forcing scenarios from CMIP5 (two RCP8.5 forcings, Representative Concentration Pathway) and CMIP6 (one SSP5-8.5 forcing, Shared Socioeconomic Pathway). Our simulations predict 5.2–8.7 mm of sea-level rise from Humboldt Glacier, significantly higher than a previous estimate (∼ 3.5 mm) and equivalent to a substantial fraction of the 40–140 mm predicted by ISMIP6 from the whole Greenland Ice Sheet. Our larger future sea-level rise prediction results from the transient calibration of our basal friction law to match the observed speedup, which requires a semi-plastic bed rheology. In many simulations, our model predicts the growth of a sizable ice shelf in the middle of the 21st century. Thus, atmospheric warming could lead to more retreat than predicted here if increased surface melt promotes hydrofracture of the ice shelf. Our data-constrained simulations of Humboldt Glacier underscore the sensitivity of model predictions of Greenland outlet glacier response to warming to choices of basal shear stress and iceberg calving parameterizations. Further, transient calibration of these parameterizations, which has not typically been performed, is necessary to reproduce observed behavior. Current estimates of future sea-level rise from the Greenland Ice Sheet could, therefore, contain significant biases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

MBMS1.0: An Open-Source Code for Modeling and Simulation of Membrane-Based Dehumidification and Energy Recovery

Air dehumidification processes are important and significantly impact both the performance of building ventilation and air conditioning systems, and the health of occupants. In a membrane-based dehumidifier, a permselective membrane is adopted to separate the two flow streams at feed and permeate sides. The feed side is continuously supplied with humid air at atmospheric pressure, while the permeate side is connected to a vacuum pump to produce vacuum or a low permeate vapor pressure as a driving force for water vapor permeation. The software MBMS 1.0 enables modeling and simulation of the design and configuration optimization of a membrane-based air-to-air (A2A) or air-to-vapor (A2V) device. A first-order physical-constant model was developed to reasonably account for the fundamental heat and mass transfer processes involved.

Gao, Zhiming↗
DOE OSTI2017

Franck–Condon Models for Simulating the Band Shape of Electronic Absorption Spectra

Band shape is an essential ingredient in the simulation of electronic absorption spectra. The excitation of multiple series of vibrational levels during an electronic excitation is a main contributor to band shapes. Here we present two simple models based on the Franck–Condon displaced-harmonic-oscillator model. The models are both derived from the time-dependent formulation of electronic spectroscopy. They assume that the transition dipoles do not depend on geometry and that the potential energy surfaces are locally quadratic; one model is second order in time and is called LQ2, and the other is third order in time and is called LQ3. These models are suitable for simulating the unresolved vibronic band shapes of electronic spectra that involve many vibrational modes. The models are straightforward and can be easily applied to simulate absorption spectra that are composed of many electronic transitions. As compared to carrying out molecular dynamics simulations, they require relatively few electronic structure calculations, and the additional cost for constructing the spectra is negligible. Therefore, the models are suitable for simulating the spectra of complex systems such as transition-metal complexes.

UV/Vis spectroscopy↗
DOE OSTI2020

Functional Requirements for the Modeling and Simulation of Advanced (Non-LWR) Reactor Mechanistic Source Term

A number of industry vendors are leading a resurgence in advanced reactor development. As the assurance of public safety from the accidental release of radionuclides to the environment is central to regulatory licensing and the reactor design process, the development of a mechanistic source term (MST) assessment has been identified by the advanced reactor industry as a foremost priority. In contrast to historic light-water reactor (LWR) source term analysis efforts, an MST attempts to realistically model the release and transport of radionuclides from the source to the environment for specific scenarios, while accounting for retention or transmutation phenomena and associated uncertainties. This work seeks to aid in the development of MST modeling and simulation tools by performing a preliminary identification of functional requirements for the major advanced reactor types.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2009

NASA Standard for Models and Simulations (M and S): Development Process and Rationale

After the Columbia Accident Investigation Board (CAIB) report. the NASA Administrator at that time chartered an executive team (known as the Diaz Team) to identify the CAIB report elements with Agency-wide applicability, and to develop corrective measures to address each element. This report documents the chronological development and release of an Agency-wide Standard for Models and Simulations (M&S) (NASA Standard 7009) in response to Action #4 from the report, "A Renewed Commitment to Excellence: An Assessment of the NASA Agency-wide Applicability of the Columbia Accident Investigation Board Report, January 30, 2004".

Zang, Thomas A.↗
NASA NTRS2014

Reply to Comment by Cole-dai Et Al. on "Climatic Impact of the Long-lasting Laki Eruption: Inapplicability of Mass-independent Sulfur Isotope Composition Measurements"

Here we respond to the comments by Cole-Dai et al. [2014] on our article Schmidt et al. [2012]. Specifically, in response to section 2 of their reply, we argued in Schmidt et al. [2012] that based on previously published estimates of the volatile release height during the 1783-1784 C.E. Laki eruption, the lack of a sulfur massindependent fractionation (MIF) anomaly is expected. In other words, no previous study on Laki ever argued that this eruption emitted SO2 into altitudes >13-15 km. In section 2.3, Cole-Dai et al. [2014] argue that the nonzero Δ33S value of their Laki sample 1 may be explained by a short-lived explosive phase at Laki during which volatiles reached the stratosphere. In Schmidt et al. [2012] in section 2, we argued in agreement with Cole-Dai et al. [2014] (section 3.1) that for a MIF anomaly to be preserved, the Laki volatiles would have had to be emitted in >20 km altitude. Our main point is that eruption column heights >20 km are unlikely based on the historical accounts and plume-rise modeling for the Laki eruption [Stothers et al., 1986; Woods, 1993; Thordarson and Self, 2003]. In Schmidt et al. [2012], we argued that to deduce a short-lived climatic impact of the Laki eruption based on the lack of a MIF anomaly and the length of the sulfate deposition in Greenland ice cores may be misleading because the climatic impact will outlast the radiative forcing of the Laki aerosol cloud. Cole-Dai et al. [2014] acknowledge the latter in their reply in section 4.2. We agreewith Cole-Dai et al. [2014] in that themagnitude and length of the climatic impact during the winter of 1783-1784 depends on the altitude of the volatile release during the eruption (sections 2.3 and 4.2). However, even if we assumed that during Laki all sulfur dioxide (SO2) would have been released in the troposphere, then the aerosol cloud would still be present in the upper troposphere during March 1784, as is evident from independent model simulations of this "tropospheric-only" scenario [Stevenson et al., 2003]. We acknowledge that there is uncertainty on the volatile release height for Laki; however, it is worth considering that those climate model simulations that used an injection altitude between 9 km and 13 km for the Laki SO2 [Highwood and Stevenson, 2003; Oman et al., 2006a, 2006b; Schmidt et al., 2012] best match the observed temperature changes during summer of 1783 [Angell and Korshover, 1985; Brázdil et al., 2010; Briffa et al., 1998; D'Arrigo and Jacoby, 1999; Jacoby et al., 1999; Kington, 1988; Manley, 1974; Parker et al., 1992; Thordarson and Self, 2003]. Based on these model simulations, a climatic impact during the winter of 1783-1784, albeit weaker than during the climactic phases of Laki, is expected (and our argument here does not exclude the role of natural variability in contributing to the cold winter of 1783-1784 as discussed in Schmidt et al. [2012]). Therefore, we continue to argue that for high-latitude eruptions such as Laki, the applicability of sulfur isotopic measurements to interpret the climatic relevance has yet to be demonstrated. Itmay transpire that the interpretation of MIF signals for the climate-relevance of an eruption is valid and unambiguous only for short-lived explosive eruptions in the tropics. In terms of the processes producing a MIF anomaly (section 3.3 in Cole-Dai et al. [2014]), the works by Hattori et al. [2013] and Ono et al. [2013] suggest that there are remaining issues not discussed by Cole-Dai et al. [2014], for instance, self-shielding of SO2 due to high column densities typical for eruptions of Pinatubo-scale and greater, and the preservation of the MIF signature in general.

general circulation model↗
DOE OSTI2026

Characterizing Turbulence at a Forest Edge: Comparing Sub-Filter Scale Turbulence Models in Simulations of Flow over a Canopy

In wildfires, atmospheric turbulence plays a major role in the transfer of turbulent kinetic energy. Understanding how turbulence feeds back into a dynamical system is important, down to the varying small scales of fuel structures (i.e. pine needles, grass). Large eddy simulations (LES) are a common way of numerically representing turbulence. The Smagorinsky model (1963) serves as one of the most studied sub-grid scale representations in LES. In this investigation, the Smagorinsky model was implemented in HIGRAD/FIRETEC, LANL’s coupled fire-atmosphere model. This study was motivated by the need to quantitatively investigate the vorticity budget equation in HIGRAD/FIRETEC. The Smagorinsky turbulent kinetic energy (TKE) was compared to FIRETEC’s 1.5-order TKE eddy-viscosity subgrid-scale model, known as the Linn turbulence model. This was done in simulations of flow over flat terrain with a homogeneous, cuboidal canopy in the center of the domain. Examinations of the modeled vertical TKE profile and turbulent statistics at the leading edge, and throughout the canopy, show that the Smagorinsky model provides comparable results to that of the original closure model posed in FIRETEC.

58 GEOSCIENCES↗
NASA NTRS2011

Standards in Modeling and Simulation: The Next Ten Years MODSIM World Paper 2010

The world has moved on since the introduction of the Distributed Interactive Simulation (DIS) standard in the early 1980s. The cold-war maybe over but there is still a requirement to train for and analyze the next generation of threats that face the free world. With the emergence of new and more powerful computer technology and techniques means that modeling and simulation (M&S) has become an important and growing, part in satisfying this requirement. As an industry grows, the benefits from standardization within that industry grow with it. For example, it is difficult to imagine what the USA would be like without the 110 volts standard for domestic electricity supply. This paper contains an overview of the outcomes from a recent workshop to investigate the possible future of M&S standards within the federal government.

Collins, Andrew J.↗