Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Integrated 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 289 records · Page 16

Final Report on the NA-22 sponsored project from FY17-19: Next-Generation Correlated Fission Event Simulations

Through this multi-organization correlated fission physics research and development project, the CGMF and FREYA fission event generators have been implemented into the MCNP6.2 code and released to the public while a variety of dedicated differential experiments on correlated fission signatures have been performed to help validate the complex fission theory and physics implemented in these codes. In addition to enabling these new capabilities, a multitude of efforts in modeling and transport simulation has been directed toward a variety of applications including general nuclear nonproliferation, subcritical neutron multiplication, and nuclear criticality safety applications. While there is still much more to be done in the area of correlated fission physics, including the study of b-delayed fission observables, more robust validation and an extension into more application spaces, one of the key successes of this project can be summarized as the delivery of a newly demonstrated and consistently integrated simulation capability spanning several complex nuclear physics codes. From the special nuclear materials source term to the radiation transport to the time-coincident detector response modeling, new and existing users around the world, especially throughout the nuclear nonproliferation communities, can now begin to model these kinds of highly complex special nuclear material applications. Another successful aspect of this project included fruitful collaborations, mentoring opportunities and community outreach. One extended collaboration workshop was held in Santa Fe, NM in March 2017. Two other limited collaboration meetings took place including a brief gathering during the mid-project Schubert review at LANL and a single all-day special session on “Correlated Fission Physics, Transport and Applications” held during the 2019 Annual Meeting of the Institute of Nuclear Materials Management (INMM) in Palm Desert, CA. Finally, as part of this project, we have mentored 14 undergraduate students, PhD candidates, and postdoctoral researchers listed in Appendix A. The remainder of this report contains detailed information followed by highlights of the goals and accomplishments in Appendix B broken down by task within this project.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

GISS Model E2.2: A Climate Model Optimized for the Middle Atmosphere. Part 2: Validation of Large-Scale Transport and Evaluation of Climate Response

Here we examine the large‐scale transport characteristics of the new ``Middle Atmosphere" NASA Goddard Institute for Space Studies climate model (E2.2). First we evaluate the stratospheric transport circulation in historical atmosphere‐only simulations integrated with interactive trace gas and aerosol chemistry. Compared to lower vertical resolution model versions, E2.2 exhibits improved tropical ascent and older stratospheric mean ages that are more consistent with observed values. In the troposphere, poleward transport to the Arctic and interhemispheric mean ages in E2.2 are comparable to models participating in the Chemistry Climate Modeling Initiative. In addition to validating E2.2 we also assess its ``transport sensitivity" using the coupled atmosphere‐ocean abrupt 4xCO2 and transient 1%CO2 simulations submitted to the Coupled Model Intercomparison Project, Phase 6, along with a 2xCO2 simulation used to evaluate the linearity of the transport circulation's response to increased CO2. We show that decreases (increases) in a stratospheric mean age (idealized surface loss) tracer scale linearly with increased lower stratospheric upwelling, which also increases linearly with warming tropical sea surface temperatures (SSTs). Abrupt 2xCO2 and 4xCO2 experiments constrained with (fixed) pre‐industrial SSTs are also used to quantify the relative importance of rapid adjustments versus SST feedbacks to the transport circulation responses in the model. Finally, sensitivity experiments are presented to illustrate the impact of changes in the convective parameterization on stratospheric transport.

GISS Model E2.2↗

Integration of Advanced Probabilistic Analysis Techniques with Multi-Physics Models

An integrated simulation platform that couples probabilistic analysis-based tools with model-based simulation tools can provide valuable insights for reactive and proactive responses to plant operating conditions. The objective of this work is to demonstrate the benefits of a partial implementation of the Small Modular Reactor (SMR) Probabilistic Risk Assessment (PRA) Detailed Framework Specification through the coupling of advanced PRA capabilities and accurate multi-physics plant models. Coupling a probabilistic model with a multi-physics model will aid in design, operations, and safety by providing a more accurate understanding of plant behavior. This represents the first attempt at actually integrating these two types of analyses for a control system used for operations, on a faster than real-time basis. This report documents the development of the basic communication capability to exchange data with the probabilistic model using Reliability Workbench (RWB) and the multi-physics model using Dymola. The communication pathways from injecting a fault (i.e., failing a component) to the probabilistic and multi-physics models were successfully completed. This first version was tested with prototypic models represented in both RWB and Modelica. First, a simple event tree/fault tree (ET/FT) model was created to develop the software code to implement the communication capabilities between the dynamic-link library (dll) and RWB. A program, written in C#, successfully communicates faults to the probabilistic model through the dll. A systems model of the Advanced Liquid-Metal Reactor–Power Reactor Inherently Safe Module (ALMR-PRISM) design developed under another DOE project was upgraded using Dymola to include proper interfaces to allow data exchange with the control application (ConApp). A program, written in C+, successfully communicates faults to the multi-physics model. The results of the example simulation were successfully plotted.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Determining Optimal Asset Location for Rapid and Efficient Wildfire Suppression: A Simulation-Based Approach

The impact of wildfire incidents has been growing in recent years, posing a serious threat to communities at the urban-wildland interface. To address this problem, there have been growing calls to use UAVs to increase the capacity of responsible agencies to quickly and effectively suppress fires and to reduce risks associated with firefighting. One of the opportunities associated with UAVs is the ability to rapidly and autonomously operate from limited-access air bases where fires are expected to burn. This study provides an approach to determine where these air bases should be placed in order to most rapidly extinguish fires, given provided fuel distributions. This approach uses an integrated simulation of fire propagation and UAV-based suppression actions to determine how much of a given environment was burned over a range of scenarios. It then uses an optimization method to explore the space and determine the location with the least burned area. Results show the approach to efficiently and effectively provide optimal bases for single-base placements over a range of scenarios, though future work is required to adequately calibrate the model and study how it can be used in multiple-base placement problems.

Daniel Hulse↗

Determining Optimal Asset Location for Rapid and Efficient Wildfire Suppression: A Simulation-Based Approach

The impact of wildfire incidents has been growing in recent years, posing a serious threat to communities at the urban-wildland interface. To address this problem, there have been growing calls to use UAVs to increase the capacity of responsible agencies to quickly and effectively suppress fires and to reduce risks associated with firefighting. One of the opportunities associated with UAVs is the ability to rapidly and autonomously operate from limited-access air bases where fires are expected to burn. This study provides an approach to determine where these air bases should be placed in order to most rapidly extinguish fires, given provided fuel distributions. This approach uses an integrated simulation of fire propagation and UAV-based suppression actions to determine how much of a given environment was burned over a range of scenarios. It then uses an optimization method to explore the space and determine the location with the least burned area. Results show the approach to efficiently and effectively provide optimal bases for single-base placements over a range of scenarios, though future work is required to adequately calibrate the model and study how it can be used in multiple-base placement problems.

Daniel Hulse↗

Measuring the shock impedance mismatch between high-density carbon and deuterium at the National Ignition Facility

Fine-grained diamond, or high-density carbon (HDC), is being used as an ablator for inertial confinement fusion (ICF) research at the National Ignition Facility (NIF). Accurate equation of state (EOS) knowledge over a wide range of phase space is critical in the design and analysis of integrated ICF experiments. Here in this paper, we report shock and release measurements of the shock impedance mismatch between HDC and liquid deuterium conducted during shock-timing experiments having a first shock in the ablator ranging between 8 and 14 Mbar. Using ultrafast Doppler imaging velocimetry to track the leading shock front, we characterize the shock velocity discontinuity upon the arrival of the shock at the HDC/liquid deuterium interface. Comparing the experimental data with tabular EOS models used to simulate integrated ICF experiments indicates the need for an improved multiphase EOS model for HDC in order to achieve a significant increase in neutron yield in indirect-driven ICF implosions with HDC ablators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Program Development Tools and Infrastructures: Activities by ASC CSSE/CCE and Partners

Exascale class machines will exhibit a new level of complexity: they will feature an unprecedented number of cores and threads, will most likely be heterogeneous and deeply hierarchical, and offer a range of new hardware techniques (such as speculative threading, transactional memory, programmable prefetching, and programmable accelerators), which all have to be utilized for an application to realize the full potential of the machine. Additionally, users will be faced with less memory per core, fixed total power budgets, and sharply reduced MTBFs. At the same time, it is expected that the complexity of applications will rise sharply for exascale systems, both to implement new science possible at exascale and to exploit the new hardware features necessary to achieve exascale performance. This is particularly true for many of the NNSA codes, which are large and often highly complex integrated simulation codes that push the limits of everything in the system including language features. To overcome these limitations and to enable users to reach exascale performance, users will expect a new generation of tools that address the bottlenecks of exascale machines, that work seamlessly with the (set of) programming models on the target machines, that scale with the machine, that provide automatic analysis capabilities, and that are flexible and modular enough to overcome the complexities and changing demands of the exascale architectures. Further, any tool must be robust enough to handle the complexity of large integrated codes while keeping the user's learning curve low. With the ASC program, in particular the CSSE (Computational Systems and Software Engineering) and CCE (Common Compute Environment) projects, we are working towards a new generation of tools that fulfill these requirements and that provide our users as well as the larger HPC community with the necessary tools, techniques, and methodologies required to make exascale performance a reality.

43 PARTICLE ACCELERATORS↗

Energetic particle transport in 3D equilibrium. Final report

In this project, we have further developed the state-of-the-art fusion simulation code, gyrokinetic toroidal code (GTC), for studying the energetic particle (EP) transport in the 3D equilibrium. We have applied the new capabilities to simulate the EP transport by the Alfven eigenmodes in stellarators and by 3D magnetic islands in tokamaks. As part of V&V efforts for 3D capability, the 3D capabilities of GTC have been applied to study thermal transport in stellarators and tokamak with 3D magnetic fields by microturbulence, which could cause additional EP loss. Effects magnetic islands on neoclassical bootstrap current have also been studied. The computational capability developed and physics learned in this project are in support of fusion SciDAC Center for Integrated Simulation of Energetic Particles in Burning Plasmas (ISEP). During this project, the PI has supervised 6 postdoctoral researchers and 4 PhD students at the University of California, Irvine. A complete list of 49 publications and 23 invited talks partially supported by this grant since 7/2015 appears at the end of this report. Five key publications closely related to this project are highlighted in this report.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The development of a high-resolution Eulerian radiation-hydrodynamics simulation capability for laser-driven Hohlraums

Hohlraums are hollow cylindrical cavities with high-Z material walls used to convert laser energy into uniform x-ray radiation drives for inertial confinement fusion capsule implosions and high energy density physics experiments. Credible computational modeling of hohlraums requires detailed modeling and coupling of laser physics, hydrodynamics, radiation transport, heat transport, and atomic physics. We report on improvements to Los Alamos National Laboratory's xRAGE radiation-hydrodynamics code in order to enable hohlraum modeling. xRAGE's Eulerian hydrodynamics and adaptive mesh refinement make it uniquely well suited to study the impacts of multiscale features in hohlraums. In order to provide confidence in this new modeling capability, we demonstrate xRAGE's ability to produce reasonable agreement with data from several benchmark hohlraum experiments. We also use xRAGE to perform integrated simulations of a recent layered high density carbon capsule implosion on the National Ignition Facility in order to evaluate the potential impacts of the capsule support tent, mixed cell conductivity methodologies, plasma transport, and cross-beam energy transfer (XBT). We find that XBT, seeded by plasma flows in the laser entrance hole (LEH), causes a slight decrease in energy coupling to the capsule and that all of these impact the symmetry of the x-ray drive such that they have an appreciable impact on the capsule implosion shape.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Integrated energy balance analysis for Space Station Freedom

An integrated simulation model is described which characterizes the dynamic interaction of the energy transport subsystems of Space Station Freedom for given orbital conditions and for a given set of power and thermal loads. Subsystems included in the model are the Electric Power System (EPS), the Internal Thermal Control System (ITCS), the External Thermal Control System (ETCS), and the cabin Temperature and Humidity Control System (THC) (which includes the avionics air cooling, cabin air cooling, and intermodule ventilation systems). Models of the subsystems were developed in a number of system-specific modeling tools and validated. The subsystem models are then combined into integrated models to address a number of integrated performance issues involving the ability of the integrated energy transport system of Space Station Freedom to provide power, controlled cabin temperature and humidity, and equipment thermal control to support operations.

Tandler, John↗

Materials Design using an Active Subspace-based Batch Bayesian Optimization Approach

Integrated computational materials engineering (ICME) calls for integrating simulation tools and/or experiments to develop new materials and materials systems. However, implementation of ICME approaches is challenging mainly due to the considerable computational expense of such frameworks and large dimensionality of the design space. Addressing these challenges is thus critical to the success of ICME initiatives. We present here a specific Bayesian optimization framework designed to address these two challenges. In particular, we propose an active subspace batch Bayesian optimization framework. The framework makes use of dimension reduction via the active subspace method and makes use of the ability to query in parallel via the batch Bayesian optimization approach. Here, the integration of these techniques leads to significant efficiency improvements while maintaining accuracy.

36 MATERIALS SCIENCE↗

Saturn Radiation Dose Environment Characterization

To understand the environment where a time-resolved hard x-ray spectrometer (AXIOM) might be fielded, experiments and simulations were performed to analyze the radiation dose environment underneath the Saturn vacuum dome. Knowledge of this environment is critical to the design and placement of the spectrometer. Experiments demonstrated that the machine performance, at least in terms of on-axis dose, has not significantly changed over the decades. Simulations of the off-axis dose were performed to identify possible spectrometer locations of interest. The effects from the source and dome hardware as well as source distributions and angles of incidence on the radiation environment were also investigated. Finally, a unified radiation transport model was developed for two widely used radiation transport codes to investigate the off-axis dose profiles and the time-dependent x-ray energy spectrum. The demonstrated equivalence of the unified radiation transport model between the radiation transport codes allows the team to tie future time-dependent x-ray environment calculations to previous integral simulations for the Saturn facility.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Transient modeling/analysis of hyperbolic heat conduction problems employing mixed implicit-explicit alpha method

This paper describes the evaluation of mixed implicit-explicit finite element formulations for hyperbolic heat conduction problems involving non-Fourier effects. In particular, mixed implicit-explicit formulations employing the alpha method proposed by Hughes et al. (1987, 1990) are described for the numerical simulation of hyperbolic heat conduction models, which involves time-dependent relaxation effects. Existing analytical approaches for modeling/analysis of such models involve complex mathematical formulations for obtaining closed-form solutions, while in certain numerical formulations the difficulties include severe oscillatory solution behavior (which often disguises the true response) in the vicinity of the thermal disturbances, which propagate with finite velocities. In view of these factors, the alpha method is evaluated to assess the control of the amount of numerical dissipation for predicting the transient propagating thermal disturbances. Numerical test models are presented, and pertinent conclusions are drawn for the mixed-time integration simulation of hyperbolic heat conduction models involving non-Fourier effects.

Tamma, Kumar K.↗

Department of Energy's Biological and Environmental Research Strategic Data Roadmap for Earth System Science

Rapid advances in experimental, sensor, and computational technologies and techniques are driving exponential growth in the volume, acquisition rate, variety, and complexity of scientific data. This wealth of scientifically meaningful data has tremendous potential to lead to scientific discovery. However, to achieve scientific breakthroughs, these data must be exploitable—they must be analyzed effectively and efficiently and the results shared and communicated easily within the wider Department of Energy’s (DOE’s) Biological and Environmental Research (BER) Climate and Environmental Sciences Division (CESD) community. The explosion in data complexity and scale makes these tasks exceedingly difficult to achieve, particularly given that an increasing number of disciplines are working across techniques, integrating simulation and experimental or observational results (see Table 5 in Appendix 2). Consequently, we need new approaches to data management, analysis, and visualization that provide research teams with easy-to-use and scalable end-to-end solutions. These solutions must facilitate (and where feasible, automate and capture) every stage in the data lifecycle (shown in Figure 1), from collection to management, annotation, sharing, discovery, analysis, and visualization. In addition, the core functionalities are the same across climate science communities, but they require customization to adapt to specific needs and fit into research and analysis workflows. To this end, the mission of CESD’s Data and Informatics Program is to integrate all existing and future distributed CESD data holdings into a seamless and unified environment for the acceleration of Earth system science.

54 ENVIRONMENTAL SCIENCES↗

Modelling Curvilinear Beamline Effects on Beam- and Spin-Dynamics in the Fermilab Muon g-2 Storage Ring.

A model of the Fermilab Muon g-2 Experiment storage ring is developed with the BMAD simulation library for beam- and spin-dynamics studies. Particular focus is given to the electrostatic quadrupoles, the plates of which follow the curvature of the particle design trajectory of the magnetically uniform storage ring. To model these, standard electric multipoles and a novel method based on cylindrical harmonics are investigated and found to accurately match the electric field maps produced for the bulk and fringe regions of the curved quadrupoles respectively, despite these methods being specifically for straight beamlines. The toroidal multipole description is then developed, to include curvature and longitudinal field variations (to model the fringe field regions). Various combinations of the different descriptions (including use of the field maps) are used for the bulk and fringe fields of the quadrupoles, and the toroidal multipoles are found to provide the best agreement with the measured tunes of the storage ring. Two integrators are also developed, based on the method of Wu, Forest and Robin. The first of these describes motion in a uniform magnetic dipole field with arbitrary electric field, and produces results comparable to those of the adaptive-step fourth-order Runge–Kutta integrator implemented in BMAD, suitable for tracking over thousands of turns. The second is a generalisation of the first to arbitrary magnetic vector potentials superimposed on that of the uniform dipole, with arbitrary electric fields, and was also found to agree with the Runge–Kutta integrator. Simulations of antimuon storage for 2000 turns of the storage ring were performed, where losses were studied along with the evolution of beam polarisation. The storage losses, important for reducing the systematic error in the measurement of the muon g-2, were found to be largely due to collimator shape. The correlation between momentum and spin precession frequency reflect the effects expected from the electric field correction in the T-BMT equation, on the order of 10 -6 . Particle injection was also studied, for various injection angles, quadrupole strengths and strengths of injection kick, and in the cases where experimental data was provided, the model was found to be consistent with the experiment.

43 PARTICLE ACCELERATORS↗

Mixing effects in the crystallization of supercooled quantum binary liquids

By means of Raman spectroscopy of liquid microjets, we have investigated the crystallization process of supercooled quantum liquid mixtures composed of parahydrogen (pH{sub 2}) or orthodeuterium (oD{sub 2}) diluted with small amounts of neon. We show that the introduction of the Ne impurities affects the crystallization kinetics in terms of a significant reduction of the measured pH{sub 2} and oD{sub 2} crystal growth rates, similarly to what found in our previous work on supercooled pH{sub 2}-oD{sub 2} liquid mixtures [Kühnel et al., Phys. Rev. B 89, 180201(R) (2014)]. Our experimental results, in combination with path-integral simulations of the supercooled liquid mixtures, suggest in particular a correlation between the measured growth rates and the ratio of the effective particle sizes originating from quantum delocalization effects. We further show that the crystalline structure of the mixtures is also affected to a large extent by the presence of the Ne impurities, which likely initiate the freezing process through the formation of Ne-rich crystallites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗