Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation modelling”

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,153 records · Page 64

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↗

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↗

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↗

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.↗

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↗

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 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.↗

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↗

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↗

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.↗

Passive and Active Microwave Remote Sensing of Precipitation and Latent Heating Distributions in the Tropics from TRMM

Passive and active microwave remote sensing data are analyzed to identify signatures of precipitation and vertical motion in tropical convection. A database of cloud/radiative model simulations is used to quantify surface rain rates and latent heating profiles that are consistent with these signatures. At satellite footprint-scale (approximately 10 km), rain rate and latent heating estimates are subject to significant random errors, but by averaging the estimates in space and time, random errors are substantially reduced, Bias errors have been minimized by improving the microphysics in the supporting cloud/radiative model simulations, and by imposing a consistent definition of remotely-sensed and model-simulated convective/stratiform rain coverage. Remotely-sensed precipitation and latent heating distributions in the tropics are derived from Tropical Rainfall Measuring Mission (TRMM) and Special Sensor Microwave/ Imager (SSM/ I) sensor data. The prototype Version 6 TRMM passive microwave algorithm typically yields average heating profiles with maxima between 6 and 7 km altitude for organized mesoscale convective systems. Retrieved heating profiles for individual convective systems are compared to coincident estimates based upon a combination of dual-Doppler radar and rawinsonde data. Also, large-scale latent heating distributions are compared to estimates derived from a simpler technique that utilizes observations of surface rain rate and stratiform rain proportion to infer vertical heating structure. Results of these tests will be presented at the conference.

Olson, William S.↗

PARAGRAPH: A graphics tool for performance and reliability analysis

PARAGRAPH is an animated graphics display package. It consists of two parts: an interface to CSIM (a process based simulation language) and a graphic display system. When a simulation model is executed on CSIM the interface collects pertinent performance analysis data and writes them to a file. This file is then fed to the graphic display system which depicts the execution of the simulation model visually. This report focuses on the graphical display system. Specifically it describes the user interface and features of the display system. It also explains how Interviews (which is based on X-windows) is used as the basis for the design of PARAGRAPH. The last section contains an example in which a basic load balancing simulation model is used to demonstrate the features and the capability of PARAGRAPH.

Lee, Kevin Douglas↗

Volcanic Contribution to Decadal Changes in Tropospheric Temperature

Despite continued growth in atmospheric levels of greenhouse gases, global mean surface and tropospheric temperatures have shown slower warming since 1998 than previously. Possible explanations for the slow-down include internal climate variability, external cooling influences and observational errors. Several recent modelling studies have examined the contribution of early twenty-first-century volcanic eruptions to the muted surface warming. Here we present a detailed analysis of the impact of recent volcanic forcing on tropospheric temperature, based on observations as well as climate model simulations. We identify statistically significant correlations between observations of stratospheric aerosol optical depth and satellite-based estimates of both tropospheric temperature and short-wave fluxes at the top of the atmosphere. We show that climate model simulations without the effects of early twenty-first-century volcanic eruptions overestimate the tropospheric warming observed since 1998. In two simulations with more realistic volcanic influences following the 1991 Pinatubo eruption, differences between simulated and observed tropospheric temperature trends over the period 1998 to 2012 are up to 15% smaller, with large uncertainties in the magnitude of the effect. To reduce these uncertainties, better observations of eruption-specific properties of volcanic aerosols are needed, as well as improved representation of these eruption-specific properties in climate model simulations.

Atmospheric temperature↗

Constructing Neural Network Based Models for Simulating Dynamical Systems

Dynamical systems see widespread use in natural sciences like physics, biology, and chemistry, as well as engineering disciplines such as circuit analysis, computational fluid dynamics, and control. For simple systems, the differential equations governing the dynamics can be derived by applying fundamental physical laws. However, for more complex systems, this approach becomes exceedingly difficult. Data-driven modeling is an alternative paradigm that seeks to learn an approximation of the dynamics of a system using observations of the true system. In recent years, there has been an increased interest in applying data-driven modeling techniques to solve a wide range of problems in physics and engineering. Here this article provides a survey of the different ways to construct models of dynamical systems using neural networks. In addition to the basic overview, we review the related literature and outline the most significant challenges from numerical simulations that this modeling paradigm must overcome. Based on the reviewed literature and identified challenges, we provide a discussion on promising research areas.

97 MATHEMATICS AND COMPUTING↗

Transport and chemistry of isoprene and its oxidation products in deep convective clouds

Deep convective clouds can transport trace gases from the planetary boundary layer into the upper troposphere where subsequent chemistry may impact aerosol particle formation and growth. In this modelling study, we investigate processes that affect isoprene and its oxidation products injected into the upper troposphere by an isolated deep convective cloud in the Amazon. We run a photochemical box model with coupled cloud microphysics along hundreds of individual air parcel trajectories sampled from a cloud-resolving model simulation of a convective event. The box model simulates gas-phase chemical reactions, gas scavenging by liquid and ice hydrometeors, and turbulent dilution inside a deep convective cloud. The results illustrate the potential importance of gas uptake to anvil ice in regulating the intensity of the isoprene oxidation and associated low volatility organic vapour concentrations in the outflow. Isoprene transport and fate also depends on the abundance of lightning-generated nitrogen oxide radicals (NO x = NO + NO 2 ). If gas uptake on ice is efficient and lightning activity is low, around 30% of the boundary layer isoprene will survive to the cloud outflow after approximately one hour of transport, while all the low volatile oxidation products will be scavenged by the cloud hydrometeors. If lightning NO x is abundant and gas uptake by ice is inefficient, then all isoprene will be oxidised during transport or in the immediate outflow region, while several low volatility isoprene oxidation products will have elevated concentrations in the cloud outflow. Reducing uncertainties associated with the uptake of vapours on ice hydrometeors, especially HO 2 and oxygenated organics, is essential to improve predictions of isoprene and its oxidation products in deep convective outflows and their potential contribution to new particle formation and growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental investigation of abnormal transverse flow enhancement of α particles in heavy-ion collisions

The mass dependence of the transverse flow for Z = 1–5 fragments from the collisions of 40 Ar + 27 Al, 40 Ar + 48 Ti, and 40 Ar + 58 Ni at 47 MeV/nucleon is investigated experimentally in this article. The transverse flow values are determined using the in-plane components of the fragment transverse momenta, where three conventional methods, i.e., the kinetic flow tensor method, the transverse momentum analysis method, and the azimuthal correlation method, are applied to reconstruct the reaction plane in an event-by-event basis. It is demonstrated from the comparison of the present experimental mass dependent flow measurements and the model simulations using an improved antisymmetrized molecular dynamics model that the experimentally observed abnormal α transverse flow enhancement is closely related to the reaction plane reconstruction procedure in the flow extraction. We further investigate the physical existence of the abnormal α flow behavior using a two-particle azimuthal correlation method, which allows us to provide the relative flow magnitude information with an identification of fragment charge number without the knowledge of the reaction plane differing from the three conventional methods. It is found that the relative flow magnitudes deduced from the two-particle azimuthal correlation functions with an identification of Z, with the correction for the recoil effect imposed by the momentum conservation, show a monotonically increasing trend as a function of fragment charge number, with no exception of the α flow enhancement. These results, in addition to those from the improved antisymmetrized molecular dynamics model simulations, definitely provide experimental evidences for the inexistence of the abnormal α flow behavior in the heavy-ion collisions at the present incident energy region in nature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Biomedical systems analysis program

Biomedical monitoring programs which were developed to provide a system analysis context for a unified hypothesis for adaptation to space flight are presented and discussed. A real-time system of data analysis and decision making to assure the greatest possible crew safety and mission success is described. Information about man's abilities, limitations, and characteristic reactions to weightless space flight was analyzed and simulation models were developed. The predictive capabilities of simulation models for fluid-electrolyte regulation, erythropoiesis regulation, and calcium regulation are discussed.

Source record↗