Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation models”

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,711 records · Page 95

Particle-in-cell simulations of magnetically driven reconnection using laser-powered capacitor coils

In this study, we propose an experimental scheme to fulfill magnetically driven reconnections. Here, two laser beams are focused on a capacitor-coil target and then strong currents are wired in two parallel circular coils. Magnetic reconnection occurs between the two magnetic bubbles created by the currents in the two parallel circular coils. A two-dimensional particle-in-cell simulation model in the cylindrical coordinate is used to investigate such a process, and the simulations are performed in the ( r , z ) plane. The results show that with the increase of the currents in the two coils, the associated magnetic bubbles expand and a current sheet is formed between the two bubbles. Magnetic reconnection occurs when the current sheet is sufficiently thin. A quadrupole structure of the magnetic field in the θ direction ( B θ ) is generated in the diffusion region and a strong electron current along the r direction ( J e r ) is also formed due to the existence of the high-speed electron flow away from the X line in the center of the outflow region. Finally, because the X line is a circle along the θ direction, the convergence of the plasma flow around r = 0 will lead to the asymmetry of J e r and B θ between the two outflow regions of magnetic reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamical regimes of vortex flow in type-II superconductors with parallel twin boundaries

We explore the dynamics of driven magnetic flux lines in disordered type-II superconductors in the presence of twin boundaries oriented parallel to the direction of the applied magnetic field, using a three-dimensional elastic line model simulated with Langevin molecular dynamics. The lines are driven perpendicular to the planes to model the effect of an electric current applied parallel to the planes and perpendicular to the magnetic field. A study of the long-time non-equilibrium steady states for varying sample orientation and thickness reveals a rich collection of dynamical regimes spanning the depinning crossover region that separates the pinned and moving-lattice states of vortex matter. We observe the emergence of a preferred direction for the ordering of the Abrikosov lattice in the free- owing vortex regime due to asymmetric pinning by the planar defects. Furthermore, we have performed novel direct measurements of flux line excitations such as half-loops and double kinks to aid the characterization of the topologically rich flux flow profile.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Development of a Two-Fluid Molten Salt Reactor Safeguards Model

This work outlines the development of a two-fluid molten salt reactor process and safeguards model. The model is split into two parts consisting of a process model and a safeguards model. The process model is based on a design by Flibe Energy, the Liquid-Fluoride Thorium Reactor, which is a two-fluid molten salt reactor that performs full salt processing on-site to remove fission products and re-fuel the reactor. The model simulates feed and consumption rates of the reactor fuel and blanket salts. The process model includes the reactor core and salt processing loops. The reactor core model has a robust architecture that allows for integration with other tools and data sets as they become available. A majority of the effort to date has been focused on the process model, and the safeguards model will be developed in detail in future work. A preliminary safeguards analysis was performed based on actinide inventories, and a preliminary materials accountancy approach was initialized. The results of this analysis are presented along with a description of the model development.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Test Information Management System (Report FY 2020 and FY 2021)

Delivery Environments supported work in Fiscal Years 2020 and 2021 to develop the Test Information Management System (TIMS). This effort involved creating a suite of tools consisting of a database and the uploading and downloading scripts. It has been demonstrated that the TIMS database is suitable for archiving raw measurements performed on various measured test articles. The TIMS database has been populated with data from the TRUST projects in FY 20 and FY 21. This included the functional data in the Sensors records as well as metadata stored in records in several other tables. The performed work demonstrated that a database archiving validation test data and integrating it with engineering analysis simulations has a great potential to increase the efficiency, responsiveness, and confidence in the modeling, simulation, and validation efforts for future (and current) weapon systems and assemblies in normal and abnormal environments. In this report, we also discuss how the records belonging to the same project should be arranged, what are the minimum requirements for linking the records using the tabular links, and we provide several recommendations to the projects engineers and the database administrator to improve the TIMS database.

97 MATHEMATICS AND COMPUTING↗

ORNL Design Engineering Library: Validated Thermoelectric Digital Twins

Radioisotope thermoelectric generators (RTGs) serve a crucial role in supplying thermal and electrical energy for reliable and long-duration power in remote and extreme environments. The thermal energy is provided by the decay of radioisotopes and is converted into electrical energy using the steady-state thermoelectric process by exploiting the Seebeck effect. Maintaining a thermal gradient through the material is necessary to drive current production. Modern predictive multi-physics tools are able to effectively describe the complex and strongly coupled phenomena needed for accurate model-based system engineering efforts needed to develop higher performing RTGs. Experimental validation of model simulations is essential to confirm the accuracy, reliability, and applicability of predictive multiphysics tools. Validation fosters confidence among users, helps meet regulatory requirements, and contributes to the ongoing improvement of simulation techniques, ultimately leading to better, safer, and more efficient designs and processes. Validated models can then be used as digital twins and can be interrogated to understand and quantify the performance gaps between theoretical and physical systems. These insights can be used to build higher performing systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

FAST Low-Energy Beamline Studies: Toward High-Peak 5D Brightness Beams for FAST-GREENS

The FAST beamline is the injector for the planned Gamma-Ray Electron ENhanced Source (GREENS) program, which aims to achieve the demonstration and first application of a high-efficiency, high-average-power free-electron laser at 515 nm. FAST-GREENS requires high 5D peak brightness; transverse normalized projected emittances of 3 mm-mrad and a peak current of 600 A are the minimum beam requirements for the FEL to reach the 10% efficiency goal. In this work, studies of the low-energy section of the FAST beamline are presented toward these ends, including preliminary measurements of beam compression and beam emittance. An effort toward developing a high-fidelity simulation model that could be later optimized for FAST-GREENS is presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The Structure of Liquid and Amorphous Hafnia

Understanding the atomic structure of amorphous solids is important in predicting and tuning their macroscopic behavior. Here, we use a combination of high-energy X-ray diffraction, neutron diffraction, and molecular dynamics simulations to benchmark the atomic interactions in the high temperature stable liquid and low-density amorphous solid states of hafnia. The diffraction results reveal an average Hf–O coordination number of ~7 exists in both the liquid and amorphous nanoparticle forms studied. The measured pair distribution functions are compared to those generated from several simulation models in the literature. We have also performed ab initio and classical molecular dynamics simulations that show density has a strong effect on the polyhedral connectivity. The liquid shows a broad distribution of Hf–Hf interactions, while the formation of low-density amorphous nanoclusters can reproduce the sharp split peak in the Hf–Hf partial pair distribution function observed in experiment. The agglomeration of amorphous nanoparticles condensed from the gas phase is associated with the formation of both edge-sharing and corner-sharing HfO6,7 polyhedra resembling that observed in the monoclinic phase.

36 MATERIALS SCIENCE↗

Clay Generic Disposal System Model - Sensitivity Analysis for 32 PWR Assembly Canisters (+2 associated model files).

The Used Fuel Disposition Campaign (UFDC), as part of the DOE Office of Nuclear Energy’s (DOE-NE) Fuel Cycle Technology program (FCT) is investigating the disposal of high level radioactive waste (HLW) and spent nuclear fuela (SNF) in a variety of geologic media. The feasibility of disposing SNF and HLW in clay media has been investigated and has been shown to be promising [Ref. 1]. In addition the disposal of these wastes in clay media is being investigated in Belgium, France, and Switzerland. Thus, Argillaceous media is one of the environments being considered by UFDC. As identified by researchers at Sandia National Laboratory, potentially suitable formations that may exist in the U.S. include mudstone, clay, shale, and argillite formations [Ref. 1]. These formations encompass a broad range of material properties. In this report, reference to clay media is intended to cover the full range of material properties. This report presents the status of the development of a simulation model for evaluating the performance of generic clay media. The clay Generic Disposal System Model (GDSM) repository performance simulation tool has been developed with the flexibility to evaluate not only different properties, but different waste streams/forms and different repository designs and engineered barrier configurations/ materials that could be used to dispose of these wastes.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Atomistically resolved hot exciton relaxation dynamics in CdSe quantum dots: Experiment and theory

Semiconductor quantum dots (QDs) are well known to give rise to a quantum confined structure of excitons. Because of this quantum confinement, new physics of hot exciton relaxation dynamics arises. Decades of work using transient absorption (TA) spectroscopy have yielded initial simple observations, such as estimates of the cooling rate from single pump photon energy experiments. More detailed TA experiments employed variable pump photon energies to measure excitonic state-resolved transition rates. These TA measurements, usually the simplest form, have been employed to characterize QDs and their relaxation dynamics to this day. Yet, these TA measurements are fundamentally lacking in their ability to measure energy-resolved hot exciton cooling, which requires observation of the full cooling history through the real excitonic manifold. Here, we employ coherent multi-dimensional spectroscopy (CMDS) to perform an atomistically directed study of hot exciton cooling in CdSe QDs, revealing energy resolved relaxation dynamics. CMDS experiments are compared with simulations and prior TA measurements and simpler theories. Our findings reveal a hot exciton relaxation dynamics landscape. This relaxation dynamics landscape is a linear or sub-linear function of excess energy for different structures of QDs, with a strong size dependence. Our model simulations parameterized by the empirical pseudopotential model reproduces the experimental functional form and the dependence upon QD diameter and shell.

Atomistic simulations↗

Mathematical model of a personalized neoantigen cancer vaccine and the human immune system

Cancer vaccines are an important component of the cancer immunotherapy toolkit enhancing immune response to malignant cells by activating CD4 + and CD8 + T cells. Multiple successful clinical applications of cancer vaccines have shown good safety and efficacy. Despite the notable progress, significant challenges remain in obtaining consistent immune responses across heterogeneous patient populations, as well as various cancers. We present a mechanistic mathematical model describing key interactions of a personalized neoantigen cancer vaccine with an individual patient’s immune system. Specifically, the model considers the vaccine concentration of tumor-specific antigen peptides and adjuvant, the patient’s major histocompatibility complexes I and II copy numbers, tumor size, T cells, and antigen presenting cells. We parametrized the model using patient-specific data from a clinical study in which individualized cancer vaccines were used to treat six melanoma patients. Model simulations predicted both immune responses, represented by T cell counts, to the vaccine as well as clinical outcome (determined as change of tumor size). This model, although complex, can be used to describe, simulate, and predict the behavior of the human immune system to a personalized cancer vaccine.

60 APPLIED LIFE SCIENCES↗

A New Signal Model for Axion Cavity Searches from N -body Simulations

Signal estimates for direct axion dark matter (DM) searches have used the isothermal sphere halo model for the last several decades. While insightful, the isothermal model does not capture effects from a halo’s infall history nor the influence of baryonic matter, which has been shown to significantly influence a halo’s inner structure. The high resolution of cavity axion detectors can make use of modern cosmological structure-formation simulations, which begin from realistic initial conditions, incorporate a wide range of baryonic physics, and are capable of resolving detailed structure. This work uses a state-of-the-art cosmological N -body+Smoothed-Particle Hydrodynamics simulation to develop an improved signal model for axion cavity searches. Signal shapes from a class of galaxies encompassing the Milky Way are found to depart significantly from the isothermal sphere. A new signal model for axion detectors is proposed and projected sensitivity bounds on the Axion DM eXperiment (ADMX) data are presented.

79 ASTRONOMY AND ASTROPHYSICS↗

An innovative method to build a comprehensive kinetic model for air injection using TGA/DSC experiments

In order to study the AIP (air injection process), a practical method was developed to obtain the kinetic data for an AIP by matching the simulation model results with the TGA (thermogravimetry analysis) experiments in our previous work. However, one shortcoming for that method is the non-uniqueness of the kinetic data obtained. The other shortcoming is that TGA cannot measure the exothermic characteristics of a crude oil. The enthalpy values of the reactions are crucial for developing a comprehensive kinetic model as the enthalpy represents the heat generation of the chemical reaction and dominants the heat effect during the AIP. Here in this study, an innovative method is proposed to build a comprehensive kinetic model by combining both TGA/DSC (differential scanning calorimetry) experiments and simulation. In the low temperature interval, the distillation dominant stage is differentiated from the LTO (low temperature oxidation) stage by comparing the nitrogen purging TGA experiments and the air purging TGA experiments. Then, the LTO stage is further divided into several sub-divided reaction stages with corresponding temperature intervals by analyzing the TGA/DSC data. The kinetic data are estimated by applying the Arrhenius method to sub-divided reaction stages of the TGA experiments. In order to obtain the enthalpy values for oxidation reactions, the negative heat flow of the distillation process is eliminated by subtracting the nitrogen purging DSC data from the air purging DSC data, and the enthalpy values of corresponding reaction stages are obtained from the subtracted heat flow curve whose area of the exothermic peak represents the enthalpy value of a stage. Moreover, the pseudo components are defined based on the nitrogen purging TGA experiments which are associated with the sub-divided reaction stages defined previously to complete the kinetic model. Finally, the kinetic model is validated with the air purging TGA experiments. This study shows a detailed workflow to build a comprehensive kinetic model for AIP by TGA/DSC experiments. An application to a crude oil is also presented to demonstrate its practicability.

42 ENGINEERING↗

Sensitivity of Near-Surface Variables in the RUC Land Surface Model in the Weather Research and Forecasting Model

In this study, we investigate the parametric sensitivity of near-surface variables, such as sensible heat flux, latent heat flux, ground heat flux, hub-height wind speed and land surface temperature, to the parameters used in the Rapid Update Cycle (RUC) land surface model (LSM) during a wintertime and summertime period. The model simulations are compared with observations collected from the second Wind Forecast Improvement Project (WFIP2) field campaign. The results suggest that parameters related to snow/ice and thermal processes can have significant impact on the simulated near-surface variables. Out of the 11 examined parameters, only 6 of them have considerable influences on the model behaviors and explain about 60 ~ 80 % of the estimated total variance of the simulated variables. In addition, the magnitude of the parametric sensitivity varies with season. For instance, parameters associated with snow/ice processes are dominant during the wintertime whereas those associated with thermal processes are more important during the summertime. Furthermore, the impact of the identified parameters on the simulated variables is highly related to the topography. There is a high degree of sensitivity to the parameter values over the slope region. This points out the importance of collecting field observations over steep areas to better quantity the appropriate values of key parameters. Overall, our findings provide a better understanding of the RUC LSM behavior associated with parameter uncertainties and can be used to improve the forecasting skill of land surface processes via calibration of the most uncertain model parameters.

17 WIND ENERGY↗

An information theoretic approach to use high-fidelity codes to calibrate low-fidelity codes

For many simulation models, it can be prohibitively expensive or physically infeasible to obtain a complete set of experimental data to calibrate model parameters. In such cases, one can alternatively employ validated higher-fidelity codes to generate simulated data, which can be used to calibrate the lower-fidelity code. In this paper, we employ an information-theoretic framework to determine the reduction in parameter uncertainty that is obtained by evaluating the high-fidelity code at a specific set of design conditions. These conditions are chosen sequentially, based on the amount of information that they contribute to the low-fidelity model parameters. The goal is to employ Bayesian experimental design techniques to minimize the number of high-fidelity code evaluations required to accurately calibrate the low-fidelity model. We illustrate the performance of this framework using heat and diffusion examples, a 1-D kinetic neutron diffusion equation, and a particle transport model, and include initial results from the integration of the high-fidelity thermal-hydraulics code Hydra-TH with a low-fidelity exponential model for the friction correlation factor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Experimental validation of a co-simulation architecture for modeling whole-building and detailed electrical distribution performance

This article presents an experimental validation of a co-simulation architecture for simultaneously modeling whole-building energy performance and detailed building electrical distribution system performance. The co-simulation architecture consists of a whole-building energy model (EnergyPlus®) embedded within a Modelica-based building electrical distribution system library called the Building Electrical Efficiency Analysis Model (BEEAM) using the Functional Mock-up Interface standard. We validate the model using experimental data collected at a full-scale test cell within Lawrence Berkeley National Laboratory’s FLEXLAB® facility. In conclusion, we show that the co-simulation model accurately predicts the electrical, mechanical, and thermal performance of the test cell for typical loads with both an AC and a DC electrical distribution topology.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Uncertain Spatial Pattern of Future Land Use and Land Cover Change and Its Impacts on Terrestrial Carbon Cycle Over the Arctic–Boreal Region of North America

Land use and land cover change (LULCC) represents a key process of human-Earth system interaction and has profound impacts on terrestrial ecosystem carbon cycling. As a key input for ecosystem models, future gridded LULCC data is typically spatially downscaled from regional LULCC projections by integrated assessment models, such as the Global Change Analysis Model (GCAM). The uncertainty associated with the different spatial downscaling methods and its impacts on the subsequent model projections have been historically ignored and rarely examined. This study investigated this problem using two representative spatial downscaling methods and focused on their impacts on the carbon cycle over the Arctic-Boreal Vulnerability Experiment (ABoVE) domain, where extensive LULCC is expected. Specifically, we used the Future Land Use Simulation model (FLUS) and the Demeter model to generate 0.25° gridded LULCC data (i.e., LULCC FLUS and LULCC Demeter , respectively) with the same input of regional LULCC projections from GCAM, under both the low (i.e., SSP126) and high (i.e., SSP585) greenhouse gas emission scenarios. The two sets of downscaled LULCC were used to drive the Community Land Model version 5 and prognostically simulate the terrestrial carbon cycle dynamics over the 21st century. The results suggest large spatial-temporal differences between LULCC FLUS and LULCC Demeter , and the spatial distributions of the needleleaf evergreen boreal tree, broadleaf deciduous boreal tree, broadleaf deciduous boreal shrub, and C3 arctic grass are particularly different under both SSP126 and SSP585. Additionally, the spatiotemporal differences are larger under SSP126 than SSP585, due to more intensive LULCC under SSP126 than SSP585 from GCAM projection. The differences in LULCC further lead to large discrepancies in the spatial patterns of projected gross primary productivity, ecosystem respiration, and net ecosystem exchange, which represent more than 79% of the contributions of future LULCC in 2100. Additionally, the difference in carbon flux under SSP126 is generally larger than those under SSP585. This study highlights the importance of considering the uncertainties induced by the spatial downscaling process in future LULCC projections and carbon cycle simulations.

54 ENVIRONMENTAL SCIENCES↗

A sharp interface model for deterministic simulation of dendrite growth

A high order level set model is developed for deterministic simulation of dendritic growth in unstable solidifying systems. The model captures motion of the front implicitly on a structured finite difference grid, enables calculation of its geometric properties and also applies boundary conditions on the immersed interface. Interfacial capillary effect is incorporated in the model through the Gibbs-Thomson condition. Canonical problems for evaluating grid convergence of the numerical method and validation tests for stability of a growing nucleus in the presence of isotropic surface tension are presented. The growth morphology of solidifying nuclei in undercooled metallic melts is quantitatively analyzed. Effects of crystal anisotropy and melt undercooling on the front geometry, propagation speed and formation of branched dendritic structures are examined. In conclusion, the complex morphological changes, such as remelting of secondary perturbations under specific conditions of undercooling, are also captured during the quantitative analysis.

36 MATERIALS SCIENCE↗

Biomass sorghum production in Iowa, USA

This dataset includes field measurement and model simulation data from Josh Bendorf's (former MS student with Drs. VanLoocke & Heaton, ISU Dept. of Agronomy) thesis project focusing on the simulation of biomass sorghum ( Sorghum bicolor L. Moench) in the state of Iowa. Measured field data are provided in CSV file format for 2019 & 2020. Model output files are provided as NetCDF file format for the state of Iowa (1998-2007) & the larger US region (1969-2018).

Agricultural production systems simulation↗