Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “modelling and 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,495 records · Page 83

Pitfalls and guidelines for the numerical evaluation of moderate-order system frequency response

The design and evaluation of a feedback control system via frequency response methods relies heavily upon numerical methods. In application, one can usually develop low order simulation models which for the most part are devoid of numerical problems. However, when complex feedback interactions, for example, between instrument control systems and their flexible mounting structure, must be evaluated, simulation models become moderate to large order and numerical problems become common. A large body of relevant numerical error analysis literature is summarized in a large language understandable to nonspecialists. The intent is to provide engineers using simulation models with an engineering feel for potential numerical problems without getting intertwined in the complexities of the associated mathematical theory. Guidelines are also provided by suggesting alternate state of the art methods which have good numerical evaluation characteristics.

Frisch, H. P.↗

Agent-based modeling and simulation for the circular economy: Lessons learned and path forward

Circular economy aims at decoupling human activities from resource use and creating wealth. However, many have questioned the link between increased circularity and sustainability, resulting in several methodological approaches being developed to answer that question. This article analyzes and discusses the insights gained from applying agent-based modeling and simulation to study the techno-economic and social conditions promoting circularity and sustainability. This article analyzes the benefits and limitations of this technology and discusses future methodology developments within the circular economy context. Moreover, six limits of the circular economy concept are used to interpret insights from the literature: thermodynamic limits, system boundary limits, limits posed by the physical scale of the economy, limits posed by path dependencies and lock-in, limits of governance and management, and limits of social and cultural definitions. Promising research avenues are to use this methodology with machine learning, industrial ecology methods, and detailed geographic information.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

WSMR FBR Modeling

A dynamic system simulation model of the WSMR Fast Burst Reactor (FBR) has been developed to simulate the time-dependent behavior of prompt burst excursions. The model combines neutron kinetics with fuel material hydrodynamics to track both fission power and the thermo-physical response of the fuel. This model is used to characterize prompt burst behavior, such as spike yield and pulse FWHM, and to serve as a time-dependent source for predicting associated FBR phenomena, such as electric field generation and atmospheric fluorescence.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Developments in Atmosphere Revitalization Modeling and Simulation

"NASA's Advanced Exploration Systems (AES) program is pioneering new approaches for rapidly developing prototype systems, demonstrating key capabilities, and validating operational concepts for future human missions beyond Earth orbit" (NASA 2012). These forays beyond the confines of earth's gravity will place unprecedented demands on launch systems. They must not only blast out of earth's gravity well as during the Apollo moon missions, but also launch the supplies needed to sustain a crew over longer periods for exploration missions beyond earth's moon. Thus all spacecraft systems, including those for the separation of metabolic carbon dioxide and water from a crewed vehicle, must be minimized with respect to mass, power, and volume. Emphasis is also placed on system robustness both to minimize replacement parts and ensure crew safety when a quick return to earth is not possible. Current efforts are focused on improving the current state-of-the-art systems utilizing fixed beds of sorbent pellets by evaluating structured sorbents, seeking more robust pelletized sorbents, and examining alternate bed configurations to improve system efficiency and reliability. These development efforts combine testing of sub-scale systems and multi-physics computer simulations to evaluate candidate approaches, select the best performing options, and optimize the configuration of the selected approach, which is then implemented in a full-scale integrated atmosphere revitalization test. This paper describes the development of atmosphere revitalization models and simulations. A companion paper discusses the hardware design and sorbent screening and characterization effort in support of the Atmosphere Revitalization Recovery and Environmental Monitoring (ARREM) project within the AES program.

Knox, James C.↗

Tropospheric Jet Response to Antarctic Ozone Depletion: An Update with Chemistry-Climate Model Initiative (CCMI) Models

The Southern Hemisphere (SH) zonal-mean circulation change in response to Antarctic ozone depletion is re-visited by examining a set of the latest model simulations archived for the Chemistry-Climate Model Initiative (CCMI) project. All models reasonably well reproduce Antarctic ozone depletion in the late 20th century. The related SH-summer circulation changes, such as a poleward intensification of westerly jet and a poleward expansion of the Hadley cell, are also well captured. All experiments exhibit quantitatively the same multi-model mean trend, irrespective of whether the ocean is coupled or prescribed. Results are also quantitatively similar to those derived from the Coupled Model Intercomparison Project phase 5 (CMIP5) high-top model simulations in which the stratospheric ozone is mostly prescribed with monthly- and zonally-averaged values. These results suggest that the ozone-hole-induced SH-summer circulation changes are robust across the models irrespective of the specific chemistry-atmosphere-ocean coupling.

chemistry-climate models↗

Material Protection, Accounting, and Control Technologies (MPACT): Modeling and Simulation Roadmap

The development of sustainable advanced nuclear fuel cycles is a long-term goal of the Office of Nuclear Energy’s (DOE-NE) Fuel Cycle Technologies program. The Material Protection, Accounting, and Control Technologies (MPACT) campaign is supporting research and development (R&D) of advanced instrumentation, analysis tools, and integration methodologies to meet this goal. This advanced R&D is intended to facilitate safeguards and security by design of fuel cycle facilities. The lab-scale demonstration of a virtual facility, distributed test bed, that connects the individual tools being developed at National Laboratories and university research establishments, is a key program milestone for 2020. These tools will consist of instrumentation and devices as well as computer software for modeling. To aid in framing its long-term goal, during FY16, a modeling and simulation roadmap is being developed for three major areas of investigation: (1) radiation transport and sensors, (2) process and chemical models, and (3) shock physics and assessments. For each area, current modeling approaches are described, and gaps and needs are identified.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Turbulent, extreme multi-zone model for simulating flux and polarization variability in blazars

The author presents a model for variability of the flux and polarization of blazars in which turbulent plasma flowing at a relativistic speed down a jet crosses a standing conical shock. The shock compresses the plasma and accelerates electrons to energies up to γ{sub max} ≳ 10{sup 4} times their rest-mass energy, with the value of γ{sub max} determined by the direction of the magnetic field relative to the shock front. The turbulence is approximated in a computer code as many cells, each with a uniform magnetic field whose direction is selected randomly. The density of high-energy electrons in the plasma changes randomly with time in a manner consistent with the power spectral density of flux variations derived from observations of blazars. The variations in flux and polarization are therefore caused by continuous noise processes rather than by singular events such as explosive injection of energy at the base of the jet. Sample simulations illustrate the behavior of flux and linear polarization versus time that such a model produces. The variations in γ-ray flux generated by the code are often, but not always, correlated with those at lower frequencies, and many of the flares are sharply peaked. The mean degree of polarization of synchrotron radiation is higher and its timescale of variability shorter toward higher frequencies, while the polarization electric vector sometimes randomly executes apparent rotations. The slope of the spectral energy distribution exhibits sharper breaks than can arise solely from energy losses. All of these results correspond to properties observed in blazars.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling and Simulation of Variable Mass, Flexible Structures

The advent of the new Ares I launch vehicle has highlighted the need for advanced dynamic analysis tools for variable mass, flexible structures. This system is composed of interconnected flexible stages or components undergoing rapid mass depletion through the consumption of solid or liquid propellant. In addition to large rigid body configuration changes, the system simultaneously experiences elastic deformations. In most applications, the elastic deformations are compatible with linear strain-displacement relationships and are typically modeled using the assumed modes technique. The deformation of the system is approximated through the linear combination of the products of spatial shape functions and generalized time coordinates. Spatial shape functions are traditionally composed of normal mode shapes of the system or even constraint modes and static deformations derived from finite element models of the system. Equations of motion for systems undergoing coupled large rigid body motion and elastic deformation have previously been derived through a number of techniques [1]. However, in these derivations, the mode shapes or spatial shape functions of the system components were considered constant. But with the Ares I vehicle, the structural characteristics of the system are changing with the mass of the system. Previous approaches to solving this problem involve periodic updates to the spatial shape functions or interpolation between shape functions based on system mass or elapsed mission time. These solutions often introduce misleading or even unstable numerical transients into the system. Plus, interpolation on a shape function is not intuitive. This paper presents an approach in which the shape functions are held constant and operate on the changing mass and stiffness matrices of the vehicle components. Each vehicle stage or component finite element model is broken into dry structure and propellant models. A library of propellant models is used to describe the distribution of mass in the fuel tank or Solid Rocket Booster (SRB) case for various propellant levels. Based on the mass consumed by the liquid engine or SRB, the appropriate propellant model is coupled with the dry structure model for the stage. Then using vehicle configuration data, the integrated vehicle model is assembled and operated on by the constant system shape functions. The system mode shapes and frequencies can then be computed from the resulting generalized mass and stiffness matrices for that mass configuration. The rigid body mass properties of the vehicle are derived from the integrated vehicle model. The coupling terms between the vehicle rigid body motion and elastic deformation are also updated from the constant system shape functions and the integrated vehicle model. This approach was first used to analyze variable mass spinning beams and then prototyped into a generic dynamics simulation engine. The resulting code was tested against Crew Launch Vehicle (CLV-)class problems worked in the TREETOPS simulation package and by Wilson [2]. The Ares I System Integration Laboratory (SIL) is currently being developed at the Marshall Space Flight Center (MSFC) to test vehicle avionics hardware and software in a hardware-in-the-loop (HWIL) environment and certify that the integrated system is prepared for flight. The Ares I SIL utilizes the Ares Real-Time Environment for Modeling, Integration, and Simulation (ARTEMIS) tool to simulate the launch vehicle and stimulate avionics hardware. Due to the presence of vehicle control system filters and the thrust oscillation suppression system, which are tuned to the structural characteristics of the vehicle, ARTEMIS must incorporate accurate structural models of the Ares I launch vehicle. The ARTEMIS core dynamics simulation models the highly coupled nature of the vehicle flexible body dynamics, propellant slosh, and vehicle nozzle inertia effects combined with mass and flexible body properties that vary significant with time during the flight. All forces that act on the vehicle during flight must be simulated, including deflected engine thrust force, spatially distributed aerodynamic forces, gravity, and reaction control jet thrust forces. These forces are used to excite an integrated flexible vehicle, slosh, and nozzle dynamics model for the vehicle stack that simulates large rigid body translations and rotations along with small elastic deformations. Highly effective matrix math operations on a distributed, threaded high-performance simulation node allow ARTEMIS to retain up to 30 modes of flex for real-time simulation. Stage elements that separate from the stack during flight are propagated as independent rigid six degrees of freedom (6DOF) bodies. This paper will present the formulation of the resulting equations of motion, solutions to example problems, and describe the resulting dynamics simulation engine within ARTEMIS.

Tobbe, Patrick A.↗

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↗

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↗

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

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↗

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↗

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↗

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↗

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↗

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