Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model based development”

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 973 records · Page 54

EMIS-CompetitiveEquilibrium.jl [SWR-19-56]

The Electricity Markets Investment Suite - Competitive Equilibrium (EMIS-CE) model is an optimization-based power system investment model developed at NREL. It attempts to approximate a market-driven investment equilibrium outcome under the assumption of idealized or perfect competition (e.g. no participants are able to exert market power or leverage exclusive cost advantages) and uncertain future conditions.

Stephen, Gord [National Renewable Energy Lab. (NRE↗

“Gray Zone” Simulations Using a Three-Dimensional Planetary Boundary Layer Parameterization in the Weather Research and Forecasting Model

Abstract Generating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal gridcell spacing (Δ x ) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δ x ≈ 100–1000 m). At these gridcell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional mesoscale simulations) nor fully resolved [as in traditional large-eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting Model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation. Significance Statement The spatial resolution of weather models continues to increase at a rapid rate in accordance with the enhancement of computing power. As a result, smaller-scale atmospheric features become more explicitly resolved. However, most numerical models still ignore the impact of horizontal weather variations on boundary layer flows, which becomes more important at these smaller spatial scales. To address this issue, we have implemented a new modeling approach, using fundamental principles, which accounts for horizontal variability. Our results show that including three-dimensional effects of turbulence is necessary to achieve realistic boundary layer characteristics. This novel technique may be useful for many applications including complex terrain flows, pollutant dispersion, and surface–atmosphere interaction studies.

Juliano, Timothy W.↗

Chemical Blast Standard (1 kg)

Chemical explosions create blast waves with large overpressure disturbances. It is important to develop a standard blast model based on data to accurately predict acoustic blast-wave amplitudes near detonations and invert for explosion energy from distant observations of blast-wave signals. However, open data from large, controlled chemical explosions with reliable ground truth can be challenging to find. The lack of access to such data could limit the number of contributions to related research and potentially stifle the rate of discoveries or validation of existing models. Here, to address these data scarcity problem, we have curated and compiled a standardized set of 817 blast-wave waveforms from 19 distinct high-explosive events. The blast-wave waveforms are standardized to a 1 kg trinitrotoluene explosion using scaling laws and corrections for location effects. A brief overview of the dataset is presented along with explosion feature models as well as recommendations for extracting explosion features. The resulting dataset is distributed to an open repository in both Seismic Analysis Code and pandas DataFrame formats containing the waveforms, the scaled distances, and the sample rates.

58 GEOSCIENCES↗

Study of Conduction-Cooled Superconducting Quadrupole Magnets Combined with Dipole Correctors for the ILC Main Linac

A superconducting rf (SRF) cryomodule for International Linear Collider(ILC) Main Linac equips focus/steering magnets. The magnets are "superferric" magnets with four superconducting (SC) race track coils conductively cooled from the cryomodule LHe supply pipe. The quadrupole field gradient and dipole field are 40 T/m and 0.1 T, respectively. The magnet length and iron-pole radius are 1 m and 0.045 m, respectively. It is known that dark current is generated at SRF cavities and accelerated through the following linac string. The dark current reaches and heats the SC magnets. It is estimated that the power deposition in the magnet may reach more than a few watts and temperature of the SC coils may locally reach to critical temperature of NbTi. It is important to make the magnet not reach quench with sufficient conduction cooling. We aim to realize the SC magnet which can stably operate under such condition. We plan to develop test coils made of three types of SC materials, NbTi, Nb₃Sn, and MgB₂ and study thermal characteristics and stability . We will develop a short model magnet, based on the test coil results. Here, we will present the magnet design study and the R&D plan.

43 PARTICLE ACCELERATORS↗

Extending and Automating a Systems-Theoretic Hazard Analysis for Requirements Generation and Analysis

Systems Theoretic Process Analysis (STPA) is a powerful new hazard analysis method designed to go beyond traditional safety techniques - such as Fault Tree Analysis (FTA) - that overlook important causes of accidents like flawed requirements, dysfunctional component interactions, and software errors. While proving to be very effective on real systems, no formal structure has been defined for STPA and its application has been ad-hoc with no rigorous procedures or model-based design tools. This report defines a formal mathematical structure underlying STPA and describes a procedure for systematically performing an STPA analysis based on that structure. A method for using the results of the hazard analysis to generate formal safety-critical, model-based system and software requirements is also presented. Techniques to automate both the analysis and the requirements generation are introduced, as well as a method to detect conflicts between the safety and other functional model-based requirements during early development of the system.

99 GENERAL AND MISCELLANEOUS↗

Review of Data from Shock Physics Experiments SPE-1, SPE-2 and SPE-3: Data Corrections and Phenomenology Analysis

The NNSA is conducting a series of high explosive (HE) tests – the Source Physics Experiments – at the National Nuclear Security Site (NNSS). These experiments are intended to investigate the generation of geophysical signals from controlled sources. These signals will aid in the development of physics-based models to aid in the discrimination of signals emitted from these and other possible sources, such as earthquakes and clandestine nuclear events. The SPE events include an array of accelerometers near to the source intended to provide measurements of the strong ground motion, or “near field,” regime. LANL performed an extensive review of the near field measurements from the first two events (SPE-1 and SPE-2) in this series. This review revealed that the near field gauges likely rotated during installation, causing measurement uncertainty. 2 This conclusion led to the installation of new gauges for the third experiment, SPE-3, utilizing new physical design and improved installation practices to achieve better quality control on their orientation. A subsequent analysis of SPE-3 data confirmed that rotations occurred and provided the basis for corrections to those data where the amount of rotation could be credibly determined. The corrected data supports a more confident analysis of the entire near field data set. This report summarizes the data errors and provides results of the analysis of the corrected dataset.

42 ENGINEERING↗

Sorption Modeling and Verification for Off-Gas Treatment

The project has made progress toward developing a comprehensive modeling capability for the capture of target species in off gas evolved during the reprocessing of nuclear fuel. The effort has integrated experimentation, model development, and computer code development for adsorption and absorption processes. For adsorption, a modeling library has been initiated to include (a) equilibrium models for uptake of off-gas components by adsorbents, (b) mass transfer models to describe mass transfer to a particle, diffusion through the pores of the particle and adsorption on the active sites of the particle, and (c) interconnection of these models to fixed bed adsorption modeling which includes advection through the bed. For single-component equilibria, a Generalized Statistical Thermodynamic Adsorption (GSTA) code was developed to represent experimental data from a broad range of isotherm types; this is equivalent to a Langmuir isotherm in the two-parameter case, and was demonstrated for Kr on INL-engineered sorbent HZ PAN, water sorption on molecular sieve A sorbent material (MS3A), and Kr and Xe capture on metal-organic framework (MOF) materials. The GSTA isotherm was extended to multicomponent systems through application of a modified spreading pressure surface activity model and generalized predictive adsorbed solution theory; the result is the capability to estimate multicomponent adsorption equilibria from single-component isotherms. This advance, which enhances the capability to simulate systems related to off-gas treatment, has been demonstrated for a range of real-gas systems in the literature and is ready for testing with data currently being collected for multicomponent systems of interest, including iodine and water on MS3A. A diffusion kinetic model for sorbent pellets involving pore and surface diffusion as well as external mass transfer has been established, and a methodology was developed for determining unknown diffusivity parameters from transient uptake data. Two parallel approaches have been explored for integrating the kernels described above into a mass-transport model for adsorption in fixed beds. In one, the GSTA isotherm kernel has been incorporated into the MOOSE framework; in the other approach, a focused finite-difference framework and PDE kernels have been developed. Issues, including oscillatory behavior in MOOSE solutions to advection-diffusion problems, and opportunities have been identified for each approach, and a path forward has been identified toward developing a stronger modeling platform. Experimental systems were established for collection of microscopic kinetics and equilibria data for single and multicomponent uptake of gaseous species on solid sorbents. The systems, which can operate at ambient temperature to 250°C and dew points from -69 to 17°C, are useful for collecting data needed for modeling performance of sorbents of interest. Experiments were conducted to determine applicable models and parameters for isotherms and mass transfer for water and/or iodine adsorption on MS3A. Validation experiments were also conducted for water adsorption on fixed beds of MS3A. For absorption, work involved modeling with supportive experimentation. A dynamic model was developed to simulate CO 2 absorption with chemical reaction using high alkaline content water solutions. A computer code was developed to implement the model based upon transient mass and energy balances. Experiments were conducted in a laboratory-scale column to determine model parameters. The influence of geometric parameters and operating variables on CO 2 absorption was studied over a wide range of conditions. This project has resulted in 7 publications, with 3 manuscripts in preparation. Also, 15 presentations were given at national meetings of ANS and AIChE and at Material Recovery and Waste Forms Campaign Working Group meetings.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Planning for Pre-Exascale Platform Environment (Fiscal Year 2015 Level 2 Milestone 5216)

This Plan for ASC Pre-Exascale Platform Environments document constitutes the deliverable for the fiscal year 2015 (FY15) Advanced Simulation and Computing (ASC) Program Level 2 milestone Planning for Pre-Exascale Platform Environment. It acknowledges and quantifies challenges and recognized gaps for moving the ASC Program towards effective use of exascale platforms and recommends strategies to address these gaps. This document also presents an update to the concerns, strategies, and plans presented in the FY08 predecessor document that dealt with the upcoming (at the time) petascale high performance computing (HPC) platforms. With the looming push towards exascale systems, a review of the earlier document was appropriate in light of the myriad architectural choices currently under consideration. The ASC Program believes the platforms to be fielded in the 2020s will be fundamentally different systems that stress ASC’s ability to modify codes to take full advantage of new or unique features. In addition, the scale of components will increase the difficulty of maintaining an errorfree system, thus driving new approaches to resilience and error detection/correction. The code revamps of the past, from serial- to vector-centric code to distributed memory to threaded implementations, will be revisited as codes adapt to a new message passing interface (MPI) plus “x” or more advanced and dynamic programming models based on architectural specifics. Development efforts are already underway in some cases, and more difficult or uncertain aspects of the new architectures will require research and analysis that may inform future directions for program choices. In addition, the potential diversity of system architectures may require parallel if not duplicative efforts to analyze and modify environments, codes, subsystems, libraries, debugging tools, and performance analysis techniques as well as exploring new monitoring methodologies. It is difficult if not impossible to selectively eliminate some of these activities until more information is available through simulations of potential architectures, analysis of systems designs, and informed study of commodity technologies that will be the constituent parts of future platforms.

97 MATHEMATICS AND COMPUTING↗

MODAL CLASSIFICATION: Processing and Anomaly Detection in a Noisy Vibrational Environment

n this report we summarize the development of a model-based detection scheme capable of the on-line detection of structural anomalies applying both system identification methods to extract a modal model and state estimation/detection methods to both detect an anomaly and locate all affected sensor positions specified by measurement location. Background information evolving from a simple single input/single output (SISO) mechanical system introducing the state-space formulation which leads directly to a multiple input/multiple output (MIMO) state-space formulation is discussed.From that representation the modal state-space model is introduced through an eigenvalue-eigenvector formulation leading to a similarity transformation and complex modes. Next the sampled-data state-space system is introduced leading to discrete state-space systems which provide the basis for the system identification methods to follow.

42 ENGINEERING↗

Dynamic Mode Decomposition of Solids

Dynamic mode decomposition (DMD) is a method that has gained notoriety in the field of turbulent fluid flows as a method for decomposing the flow field into modes that could be further deconstructed to understand their influence on the overall dynamics of the system. Forays into solids and non-linear systems were considered, but not applied. In this work, DMD was applied for the first time to the heat-diffusion and reactive heat diffusion equations on a random particle pack of uniform solid spheres. A verification of a linear heat-diffusion test problem was successful, showing equality between the normal modes and Koopman modes obtained from DMD. Further application to a non-linear reactive system revealed stability limits of the underlying modes which are dependent on microstructure and chemical kinetics. This work will enable the development of reactive material models based on further analysis using DMD to quantify the statistical dependencies of transient response on microstructural characteristics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Irradiation enhanced diffusion and diffusional creep in U₃Si₂

U₃Si₂ is a candidate as an advanced fuel for light water reactors (LWRs). Compared to UO₂, there are significant data gaps for the thermophysical and thermomechanical properties of U₃Si₂. Atomic scale simulations can help fill in those gaps, while also providing a more mechanistic understanding of the material properties. Creep is an important property for fuel performance as it influences the pellet-cladding gap, which in turn affects the temperature of the fuel (a key parameter in many other aspects of fuel behavior). In this work, we employed a DFT-informed cluster dynamics framework to predict point defect concentrations in U₃Si₂. The concentrations, diffusivities, and volumes of these defects were used to develop a diffusional creep model based on bulk (Nabarro-Herring) and grain boundary (Coble) diffusion processes. It was found that irradiation enhanced Nabarro-Herring creep is dominant at low temperatures, while Coble creep dominates for high temperatures. The model compares well against available experimental data and has been implemented in BISON. A demonstration case using simple power profiles has been carried out, showing that negligible creep occurs due to the low temperatures experienced by U₃Si₂, a consequence of its high thermal conductivity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A spatially explicit individual-based modeling approach to evaluate the cumulative effects of wind energy development on the greater sage-grouse. Final Report

In order to achieve the national vision of 20% of energy being produced by wind by 2030, issues associated with wind energy development’s potential impacts on ecological resources must be addressed. In particular, wind energy and associated transmission system development may affect wildlife populations through direct mortality, habitat fragmentation, habitat degradation, and disturbance. Although individual projects may have little effect on the sustainability of wildlife populations, cumulative impacts could be significant. The wind energy industry and regulators have developed processes for assessing the impacts of individual projects and identifying mitigation requirements, but no one has developed a methodology for considering cumulative impacts in a systematic and holistic way.

17 WIND ENERGY↗

NEAMS Problem Draft (FY2020 Assessment)

This report documents the physical configuration required to model the assessment problem posed by the Reactor Physics focus area within the Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. This model is based on the Empire model developed by staff at Los Alamos National Laboratory under the MEITNER (Modeling-Enhanced Innovations Trailblazing Nuclear Energy Reinvigoration) program. The reference report does not provide dimensions for development of computational models; those data are provided elsewhere in published. Thus, this report has been prepared to provide a single reference that can be used to create consistent models by assessment participants.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Prediction of the influence of strain magnitude and strain path on the resulting microstructure and post-forming mechanical behavior of RA containing high formability steels

In third generation advanced high strength steels (AHSS) there is considerable potential to alter the transformation characteristics through variation of process parameters and through changes in chemical composition. This proposal aims to utilize this known phenomenon to develop microstructure-based transformation models to predict austenite stability as a function of chemical composition, strain level, deformation mode. The ultimate intent is to predict microstructure and performance during room temperature forming to produce automotive structural components with superior strength and failure resistance.

36 MATERIALS SCIENCE↗

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↗

Conversion Kinetics and Ionic Conductivity in Na-β”-Alumina + YSZ (Naβ”AY) Sodium Solid Electrolyte via Vapor Phase Conversion Process

Sodium ion batteries have been receiving increasing attention and may see potential revival in the near future, particularly in large-scale grid energy storage coupling with wind and solar power generation, due to the abundant sodium resources, low cost, and sufficiently high energy density. Among the known sodium ion conductors, the Na-β”-alumina electrolyte remains highly attractive because of its high ionic conductivity. This study focuses on the vapor phase synthesis of a Na-β”-Alumina + YSZ (Naβ”AY) composite sodium electrolyte, which has higher mechanical strength and stability than conventional single phase β”-Alumina. The objectives are the measurement of conversion kinetics through a newly developed weight-gain based model and the determination of sodium ionic conductivity in the composite electrolyte. Starting samples contained ~70 vol% α-Alumina and ~30 vol% YSZ (3 mol% Y 2 O 3 stabilized Zirconia) with and without a thin alumina surface layer made by sintering in air at 1600 °C. The sintered samples were placed in a powder of Na-β”-alumina and heat-treated at 1250 °C for various periods. Sample dimensions and weight were measured as a function of heat treatment time. The conversion of α-Alumina in the α-Alumina + YSZ composite into Naβ”AY occurred by coupled diffusion of sodium ions through Na-β”-alumina and of oxygen ions through YSZ, effectively diffusing Na2O. From the analysis of the time dependence of sample mass and dimensions, the effective diffusion coefficient of Na 2 O through the sample, D eff , was estimated to be 1.74 x 10 -7 cm2 s -1 , and the effective interface transfer parameter, k eff , was estimated as 2.33 x 10 -6 cm s -1 . By depositing a thin alumina coating layer on top of the bulk composite, the chemical diffusion coefficient of oxygen through single phase Na-β”-alumina was estimated as 4.35 x 10 -10 cm 2 s -1 . An AC impedance measurement was performed on a fully converted Naβ”AY composite, and the conductivity of the composite electrolyte was 1.3 x 10 -1 S cm -1 at 300 °C and 1.6 x 10 -3 S cm -1 at 25 °C, indicating promising applications in solid state or molten salt batteries at low to intermediate temperatures.

36 MATERIALS SCIENCE↗

A DYNAMIC DENSITY FUNCTIONAL THEORY APPROACH TO DIFFUSION IN WHITE DWARFS AND NEUTRON STAR ENVELOPES

We develop a multicomponent hydrodynamic model based on moments of the Born–Bogolyubov–Green–Kirkwood–Yvon hierarchy equations for physical conditions relevant to astrophysical plasmas. These equations incorporate strong correlations through a density functional theory closure, while transport enters through a relaxation approximation. This approach enables the introduction of Coulomb coupling correction terms into the standard Burgers equations. The diffusive currents for these strongly coupled plasmas is self-consistently derived. The settling of impurities and its impact on cooling can be greatly affected by strong Coulomb coupling, which we show can be quantified using the direct correlation function.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimal Voltage Regulation for Unbalanced Distribution Networks Considering Distributed Energy Resources

With increasing penetration of distributed generation in the distribution networks (DN), the secure and optimal operation of DN has become an important concern. In this paper, an iterative quadratic constrained quadratic programming model to minimize voltage deviations and maximize distributed energy resource (DER) active power output in a three phase unbalanced distribution system is developed. The optimization model is based on the linearized sensitivity coefficients between controlled variables (e.g., node voltages) and control variables (e.g., real and reactive power injections of DERs). To avoid the oscillation of solution when it is close to the optimum, a golden search method is introduced to control the step size. Numerical simulations on modified IEEE 13 nodes test feeders show the efficiency of the proposed model. Compared to the results solved by heuristic search (harmony algorithm), the proposed model converges quickly to the global optimum.

distributed energy resources↗