Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model initialization”

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 919 records · Page 51

SRF Cavity Fault Classification and Prediction at Jefferson Lab

Over the last few years several machine learning projects at Jefferson Lab have had a common focus to optimize operation of superconducting RF (SRF) cavities in the Continuous Electron Beam Accelerator Facility (CEBAF). In this talk we highlight work to identify and classify types of faults from C100-type cavities and then to extend those capabilities to provide real-time fault prediction. Early prediction may enable mitigation strategies to prevent some types of faults. In our approach we apply a two-step fault prediction pipeline. In the first step, a model distinguishes between faulty and normal signals. In the second step, signals flagged as faulty by the first model are classified into one of seven fault types based on learned signatures in the data. Initial results show that our model can successfully predict most fault types 200 ms before onset. In additional to model performance, we also highlight challenges in working with real-world data and challenges for deploying models.

Tennant, C.↗

Modeling of the Jupiter Lase Facility Using Virtual Beamline ++ (Vbl++) Code

The Jupiter Laser Facility (JLF) at Lawrence Livermore National Laboratory (LLNL) is a mid-scale laser facility, with three different laser platforms: Janus (Target Area 1), Titan, and COMET. To help the facility and its users better understand the laser performance, a model of Janus and Titan has been developed using the Virtual Beamline ++ code (VBL++), which was initially developed by LLNL to model the National Ignition Facility beamlines and predict their performance.

43 PARTICLE ACCELERATORS↗

Modeling of the TRIGA IPR-R1 research reactor with the Serpent2/RINNOVO Nodal core analysis package

The Serpent2/RINNOVO nodal core analysis code system, a dedicated tool for modeling research reactors, has the capability to accurately predict important core physics parameters involved in the safety of reactor operation, such as various reactivity coefficients, control rod and bank worths including the shutdown margin, power distributions as well as local neutron flux predictions at various core locations of high importance, e.g., at irradiation rigs. In this work, a hexagonal model of the unrodded initial core of the Brazilian IPR-R1 Mark I type TRIGA nuclear reactor has been created using the Serpent2/RINNOVO code system. The choice of employing a hexagonal core geometry representation has mainly been made to facilitate subsequent fuel shuffling operations and core follow calculations of this reactor. Numerical results in terms of the core eigenvalue and the assembly power distribution have then been compared against corresponding full core Serpent2 results to prove feasibility of using the Serpent2/RINNOVO code system for modeling small and highly heterogeneous TRIGA reactors. Overall, RINNOVO predicts the power distribution very accurately but the eigenvalue error still remains quite large. In light of being the very first evaluation of a TRIGA reactor with the RINNOVO nodal core simulator, these results are also considered to be very preliminary. Improved accuracy is expected by incorporating a proper methodology in the current code system for computing discontinuity factors for multi-assembly configurations. Furthermore, improved results are expected by increasing the number of energy groups used by RINNOVO in these core calculations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Causal mechanism of injection-induced earthquakes through the M w 5.5 Pohang earthquake case study

Causal mechanisms for fluid injection-induced earthquakes remain a challenge to identify. Past studies largely established spatiotemporal correlations. Here, we propose a multi-process causal mechanism for injection-induced earthquakes through a case study of the 2017 M w 5.5 induced earthquake near Pohang Enhanced Geothermal System, Korea, where detailed hydraulic stimulation and on-site seismicity monitoring data provide an unprecedented opportunity. Pore pressure modeling reveals that pore pressure changes initiate seismicity on critically stressed faults and Coulomb static stress transfer modeling reveals that earthquake interactions promote continued seismicity, leading to larger events. On the basis of these results, we propose the following causal mechanism for induced seismicity: pore pressure increase and earthquake interactions lead to fault weakening and ultimately triggering larger earthquakes later in the process. We suggest that it is prudent that pore pressure change, initial seismicity locations, and Coulomb static stress transfer from seismicity earlier in the sequence are assessed in real-time.

58 GEOSCIENCES↗

Update of Residential Load Profile for WECC Load Composition Model Using Cross-Correlation Method

The Western Electricity Coordinating Council (WECC) initiated the development of a composite load model for transmission-level power system simulators. To more precisely utilize the model to represent the power consumption pattern for various locations of the WECC network, an up-to-date load composition dataset is required. The study in this paper aims to reconstruct the residential sector of the load composition model (LCM) dataset using latest survey or measurement data. The LCM residential dataset contains load profiles for twelve climate zones and five season types. The northwest energy efficiency alliance (NEEA) launched a study to generate the load profiles for four climate zones in the Pacific Northwest. Based on the NEEA database, a component-wise cross-correlation-based approach has been developed in this paper to update the obsolete residential load profiles of the four pacific northwest climate zones and populate the new load shape data into the other eight climate zones.

Liu, Yuan↗

Report on the sensitivity kernel construction and updating the SPiRaL model with waveform data

The WAVEFORMS Initiative in the Ground-based Nuclear Detonation Detection (GNDD) program includes research leading towards the prediction of entire seismic and acoustic waveforms produced by natural and manmade events, including explosions. A key aspect of this research is the development of Earth (seismic) models at multiple scales including crustal, regional, and global. The work described here pertains to the global-scale seismic tomography effort led by LLNL.

58 GEOSCIENCES↗

An Assessment of the Error Due to Computing Waste Isolation Pilot Plant Porosity Using the Porosity Response Surface Approach

The Waste Isolation Pilot Plant Performance Assessment (WIPP PA) must predict the likelihood that radionuclides will escape into the biosphere via mechanisms that depend on geohydraulic flow. Ideally, one would predict the geohydraulic flow using coupled geohydraulic and geomechanical simulations, but such coupled simulations are not computationally tractable. Instead, Sandia has historically used a look-up table of porosities for a given fluid pressure and time, called the porosity response surface, but this approach can introduce porosity errors because it largely ignores the porosity’s dependence on the past fluid pressure history. This report discusses efforts to quantify these porosity errors for both the legacy and new porosity response surfaces. Six hundred different fluid pressure histories were fed through the legacy/new geomechanical model and the legacy/new porosity response surface to generate six hundred porosity error histories. The error associated with the legacy porosity surface was substantial, while the error associated with the new porosity surface was typically small, except when fluid pressures exceeded the lithostatic pressure at the repository. In response to the errors at high pressures, a preliminary study of the WIPP PA’s sensitivity to these porosity errors was conducted. The study found that reducing the porosity errors at high pressures negligibly affected predictions of radionuclide releases. Finally, an initial machine-learned model for porosity was developed. This ML model significantly reduced the porosity error at high pressures, but sizable errors remained, so more development is necessary before coupling an ML model to the geohydraulic model.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

High Yield, Economical and Environmentally Benign Production of Rare Earth Elements from Coal Ash (Phase II Final Summary Report)

Fly ash stored in landfills and ponds across the United States is an attractive, abundant domestic resource for the cost-effective recovery of rare earth elements (REE) and other critical minerals (CM). Physical Sciences Inc. (PSI) and its team members, Winner Water Services (WWS) and University of Kentucky/Center for Applied Energy Research (UK/CAER) successfully executed a multiphase program that developed technologies and their implementation in a pilot plant. We demonstrated plant operations for cost-effective and environmentally-friendly production of rare earth element oxide (REO) concentrates, and the critical minerals scandium and aluminum (in the forms of salts or oxide products), from coal ash. We also constructed and demonstrated a research-scale (0.5 kg/day) micropilot facility to validate the key physical and chemical processing operations, predict yields, and troubleshoot process bottlenecks. The project team then designed, constructed and operated two decoupled pilot plants: (1) an operational pilot plant for physical separation processes with capacity of 0.4 metric tons per day (tpd), where we optimized processes to produce selected ash fractions as the feedstock for chemical processing and as valuable byproducts such as cenospheres, magnetic ash, and secondary fuel carbon, and (2) an operational pilot plant for chemical ash processing with a capacity of 0.5 tpd that developed optimized processes for the production of: (a) REO concentrates, (b) critical minerals (Sc, Al), and (c) beneficiated ash as a valuable byproduct suitable for cement applications. In Phase I, the project team (with Equinox Chemicals in place of WWS) developed and demonstrated the feasibility of the physical and chemical separation processes, developed the design of a pilot plant, and began the development of a preliminary techno-economic model. In the baseline (initial) Phase II program, the project team developed and demonstrated the above pilot scale plant, producing salable REE concentrates, including Y and Sc (REYSc), plus commercially viable byproducts, using environmentally safe and high-yield physical and chemical enrichment processes. The team successfully demonstrated chemical pilot design, construction, shakedown, and operations of the plant. We produced the Phase II deliverable REYSc concentrate ((50 g of >60 wt.% purity REYSc salts on elemental basis), generated the feed for the Phase II follow-on program, identified processing challenges for future optimizations, and refined the techno-economic model. In the Phase II follow-on program, the project team: (1) developed and demonstrated processes to increase the REE amount by 3X (content basis) and convert the Phase II REE salt mixture to an oxide mixture, (2) produced/delivered >38 g of REO mixture with >85 wt.% purity (elemental basis); (3) developed processes to recover critical minerals scandium and aluminum from intermediate streams; (4) produced/delivered > 1 g of scandium salt mixture with >85 wt. % purity (elemental basis); (5) produced/delivered > 100 g of aluminum oxide type material with >70% wt. purity (elemental basis); and (6) updated the techno-economic model from the baseline Phase II program to assess CAPEX and OPEX of a commercial operation. This program has developed extensive databases on process chemistry, unit operations, plant engineering, and techno-economics that will enable further scale-up toward commercial plant design. Specific future developments will be focused on achieving dramatic savings in energy, reagent usage, and operating costs. The combined results will contribute significantly for maturing the technologies of REE recovery from coal byproducts and promote the establishment of domestic REE and CM supply chains.

01 COAL, LIGNITE, AND PEAT↗

Nonequilibrium effects and transverse spherocity in ultrarelativistic proton-nucleus collisions

We investigate the effects of nonequilibrium dynamics in small colliding systems by comparing a nonequilibrium transport approach, Parton-Hadron-String Dynamics (PHSD), with a (2+1)-dimensional viscous hydrodynamic model, VISHNew. Focusing on p + Pb collisions at energies available at the Large Hadron Collider, we extract the initial conditions for the hydrodynamic model from PHSD, in order to reduce the impact of the early out-of-equilibrium dynamics and focus on the traces of nonequilibrium in the ensuing medium evolution. We find that in the transport approach quantities like energy density and bulk viscous pressure are highly inhomogeneous on the transverse plane during the whole evolution, whereas the hydrodynamic simulations dissolve more efficiently the initial spatial irregularities, still keeping a high degree of inhomogeneity due to the smaller size and lifetime of the medium produced in p + Pb collisions with respect to heavy-ion reactions. As a first step that will help to identify the impact of these nonequilibrium effects on final observables in proton-nucleus collisions, we perform an analysis of the transverse spherocity, an event-shape observable able to distinguish between jetty and isotropic configurations of transverse momenta. We found that the spherocity distribution in PHSD is slightly shifted towards the isotropic limit with respect to the hydrodynamic result. Even though this dissimilarity is partially due to the difference in the final charged particle production, it mainly comes from the different description within the two frameworks of the medium produced in small colliding systems. This finding supports the idea that multidifferential measurements, such as those based on event categorization according to multiplicity and spherocity, are useful to study final-state observables in ultrarelativistic proton-nucleus collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

FY18 J23X PSDT Final Report: Shock Initiation of PBX 9502 Lot HOL88H891-008

This report details a series of plate impact experiments on the insensitive explosive PBX 9502 (95% TATB, 5% Kel-F 800) with lot number HOL88H891-008. This is a ’virgin’ lot (no recycled material used), for which many detonation performance measurements are available; the addition of these 1-D shock initiation experiments will aid model calibration efforts for this material. Samples were machined from isostatically-pressed billets with anominal density of 1.89 g/cc. The evolution of the shock to a detonation was characterized by 9 embedded magnetic particle velocity gauges, and 116 embedded magnetic time of arrival gauges (AKA shock trackers). Pho-ton Doppler velocimetry (PDV) with a LiF window was used at the back surface of the samples to capture the wave profile exiting the sample. Six experiments were performed, covering a pressure range of 10.8–18.6 GPa and run distances to detonation of 2.7–12.2+ mm. The shock initiation of this lot is found to be consistent with Pop plots reported for other lots of PBX 9502. All data have been uploaded to the LANL Small Scale Database for use by the research community.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Surrogate Modelling of 3rd Integer Resonant Extraction at Fermilab Delivery Ring

We present an ongoing work in which a surrogate model is being developed to reproduce the response dynamics of the third-integer resonant extraction process in the Delivery Ring (DR) at Fermilab. This effort is in pursuit of smoothly extracting circulating beam to the Mu2e Experiment s production target, wherein the goal is to extract a uniform slice of the circulating $1e12$ protons in the DR over 25,000 turns (43~ms). The DR contains 3 harmonic sextupoles which excite a third-integer resonance as well as three fast, tune-ramping quadrupole magnets which drive the horizontal tune towards the $29/3$ resonance. In our initial work the surrogate model trains on a semi-analytical simulation provided in the same format as live data. Using Reinforcement Learning (and other potential ML methods), the trained surrogate acts as the environment in which a simple ML control agent could learn to dynamically adjust the quadrupole ramp at 430 break points within the 43 microsecond spill window. The control agent will be hosted on a dedicated Arria 10 FPGA, introducing its own requirements on control agent architecture. In this work we report the accuracy and fidelity of surrogate models in comparison to the response dynamics of the physics simulator.

Narayanan, Aakaash [Fermilab] (ORCID:0000000157944↗

Cross-scale observations of the 2015 St. Patrick's day storm: THEMIS, Van Allen Probes, and TWINS

In this paper, we present cross-scale magnetospheric observations of the 17 March 2015 (St. Patrick's Day) storm, by Time History of Events and Macroscale Interactions during Substorms (THEMIS), Van Allen Probes (Radiation Belt Storm Probes), and Two Wide-angle Imaging Neutral-atom Spectrometers (TWINS), plus upstream ACE/Wind solar wind data. THEMIS crossed the bow shock or magnetopause 22 times and observed the magnetospheric compression that initiated the storm. Empirical models reproduce these boundary locations within 0.7 R E . Van Allen Probes crossed the plasmapause 13 times; test particle simulations reproduce these encounters within 0.5 R E . Before the storm, Van Allen Probes measured quiet double-nose proton spectra in the region of corotating cold plasma. About 15 min after a 0605 UT dayside southward turning, Van Allen Probes captured the onset of inner magnetospheric convection, as a density decrease at the moving corotation-convection boundary (CCB) and a steep increase in ring current (RC) proton flux. During the first several hours of the storm, Van Allen Probes measured highly dynamic ion signatures (numerous injections and multiple spectral peaks). Sustained convection after ~1200 UT initiated a major buildup of the midnight-sector ring current (measured by RBSP A), with much weaker duskside fluxes (measured by RBSP B, THEMIS a and THEMIS d). A close conjunction of THEMIS d, RBSP A, and TWINS 1 at 1631 UT shows good three-way agreement in the shapes of two-peak spectra from the center of the partial RC. A midstorm injection, observed by Van Allen Probes and TWINS at 1740 UT, brought in fresh ions with lower average energies (leading to globally less energetic spectra in precipitating ions) but increased the total pressure. Finally, the cross-scale measurements of 17 March 2015 contain significant spatial, spectral, and temporal structure.

79 ASTRONOMY AND ASTROPHYSICS↗

A large-scale, agent-based simulation of metropolitan freight movements with passenger and freight market interactions

This study summarizes the first stage in the implementation of an agent-based freight modeling system that has a global representation of agents and detailed modeling of a large-scale transportation network. The model is used to evaluate the transportation and energy impacts of goods movement across urban and national scales. The framework is implemented within POLARIS, a C++-based Planning and Operations Language for Agent-based Regional Integrated Simulation, which consists of an activity-based modeling (ABM) and dynamic traffic assignment (DTA) system that has robust features for passenger travel. This platform provides a tool to model interactions among consumers, producers, and the transportation system. The main objective of this initial implementation is to implement a freight model within POLARIS following an agent-based paradigm with behavioral and simulation methods. This paper presents the initial framework and illustrates the application of the model. Building upon earlier works, a parcel location assignment algorithm for business establishments in the population is documented, along with a method for estimating establishment production and consumption volumes. In addition to population generation, other features of the model include push-pull supply chains, multimodal path choice, choice of transportation logistics node, and dynamic traffic assignment. A module with e-commerce supply and demand was also developed to analyze the effects of e-commerce delivery on last-mile energy use and congestion.

agent-based freight↗

Simulations and models for the Richtmyer–Meshkov instability with broadband perturbations

The Richtmyer–Meshkov instability (RMI) is shock driven and affects many phenomena from inertial fusion to supernova explosions. The behavior of single-modes in the RMI has been studied extensively but less is known with the broadband perturbations that occur in applications. Here, we describe extensive numerical simulations and modeling of the RMI with broadband perturbations with an initial power spectrum of the form P ∝ km, where k is the wavenumber. The hydrodynamic condition is the same as in the θ-Group Collaboration [Thornber et al., “Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The θ-group collaboration,” Phys. Fluids 29, 105107 (2017)] with a Mach ∼1.86 shock and gamma-law = 5/3 fluids with Atwood number A = −0.5. The bubble amplitude hB is found to grow in two stages. Initially, hB undergoes a phase-inversion and grows linearly in time (t) at a rate consistent with a linear theory. Asymptotically, hB grows as a power law ∼tθ when k|hB| > O(1), where k is the average wavenumber for the initial spectrum. The RMI behavior in simulations and modeling agree over a wide range of exponent m, spectral width, initial amplitude, and time. The agreement is quantified objectively using statistical analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulations and models for the Richtmyer– Meshkov instability with broadband perturbations

The Richtmyer–Meshkov instability (RMI) is shock driven and affects many phenomena from inertial fusion to supernova explosions. The behavior of single-modes in the RMI has been studied extensively but less is known with the broadband perturbations that occur in applications. Here, we describe extensive numerical simulations and modeling of the RMI with broadband perturbations with an initial power spectrum of the form P∝ k m , where k is the wavenumber. The hydrodynamic condition is the same as in the $\theta$-Group Collaboration [Thornber et al., “Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The $\theta$-group collaboration,” Phys. Fluids 29, 105107 (2017)] with a Mach ~1.86 shock and gamma-law = 5/3 fluids with Atwood number A = 0.5. The bubble amplitude h B is found to grow in two stages. Initially, h B undergoes a phase-inversion and grows linearly in time (t) at a rate consistent with a linear theory. Asymptotically, h B grows as a power law $~t^{\theta}$ when $\langle{k}\rangle{|h_B|}$ > O(1), where $\langle{k}\rangle$ is the average wavenumber for the initial spectrum. The RMI behavior in simulations and modeling agree over a wide range of exponent m, spectral width, initial amplitude, and time. The agreement is quantified objectively using statistical analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

How to sustain a CO 2 -thermosiphon in a partially saturated geothermal reservoir: Lessons learned from field experiment and numerical modeling

CO 2 has been proposed as a working fluid for geothermal energy production because of its ability to establish a self-sustaining CO 2 thermosiphon, taking advantage of the strong temperature dependence of CO 2 density. To test the concept of CO 2 heat extraction, in January 2015 a CO 2 thermosiphon was operated at the SECARB Cranfield Site, Cranfield, Mississippi, where a brine-saturated sand at a depth of 3.2 km has been under near continuous CO 2 flood since December 2009 as part of a U.S. Department of Energy demonstration of CO 2 sequestration, resulting in a partially saturated reservoir surrounding a well pair. The lateral distance between the producer and injector was 112 m at reservoir depth, a distance considered pre-commercial in scale, but great enough that thermal breakthrough was still not significant after several years of injection. Instead of producing power with a turbine, heat was extracted heat from recirculated fluid using a heat exchanger and portable chiller. The well field and surface equipment were instrumented to compare field observations with predicted responses from numerical models. Thermosiphon flow could be initiated by venting, but thereafter flow rate steadily declined, indicating that the thermosiphon was not sustainable. To model the system, the capability of T2Well, a fully coupled wellbore/reservoir numerical simulator, was expanded to enable simulation of the entire loop of fluid circulation in the fully-coupled system consisting of the injection/production wells, the reservoir, and the surface devices (heat exchanger, flow-rate regulator etc.). Combined with the newly developed TOUGH2 equation of state module called EOS7CMA, the enhanced T2 Well was used prior to the field experiment to simulate the circulation of a CO 2 -H 2 O-CH 4 mixture in a model geothermal system patterned after the Cranfield demonstration test. The model predicted that a sustainable thermosiphon could be achieved. After the field thermosiphon did not achieve the pre-test prediction of flow rates and thermosiphon sustainability, the numerical model was modified to improve realism and calibrate certain processes; it was then able to reproduce the major phenomena observed in the field. In a series of sensitivity studies, many factors were found that could potentially contribute to the failing of a sustainable thermosiphon. These factors could be categorized as two types: factors that increase the resistance to flow and factors that increase heat loss of the working fluid. In conclusion, the lessons learned can be applied to both future modeling and to achieving CO 2 -based geothermal reservoir exploitation.

15 GEOTHERMAL ENERGY↗

Skyrmion soliton motion on periodic substrates by atomistic and particle-based simulations

Abstract We compare the dynamical behavior of magnetic skyrmions interacting with square and triangular defect arrays just above commensuration using both atomistic and particle-based models. Under applied drives, the initial motion is a kink traveling through the pinned skyrmion lattice. For the square defect array, both models agree well and show a regime in which the soliton motion is locked along 45°. The atomistic model also produces locking of a soliton along 30°, which is absent in the particle-based model. For the triangular defect array, the atomistic model exhibits 30° soliton motion over a wide region of external current values. In contrast, the particle-based model gives 45° soliton motion over a small range of external driving force values. The difference arises because the nondeforming particle model facilitates meandering skyrmion orbits while the deformable atomistic model enables stronger skyrmion-skyrmion interactions that reduce the meandering. Our results indicate that soliton motion through pinned skyrmion lattices on a periodic substrate is a robust effect.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Orbit modeling of fast particle redistribution induced by sawtooth instability

Initial tests on the National Spherical Torus Experiment Upgrade (NSTX-U) (Menard et al 2012 Nucl. Fusion 083015) device suggest that introducing energy selectivity for sawtooth induced fast ion redistribution is required to improve the agreement between experimental and simulated quantities such as neutron rate and Fast-Ion D-Alpha profiles. The aim of this work is to assess the requirements to properly describe the behaviour of fast ions during a sawtooth crash for predictive sawtooth simulations. As the first step, in this work, we use the particle-following Orbit code to characterize the redistribution of fast particles. In order for a sawtooth crash to be simulated, a spatial and temporal displacement is implemented into the Orbit code. The perturbation amplitude is determined by comparison with experimental measurement of the neutron rate drop. The characteristics of fast ions with different orbit types are investigated in phase and real space. Due to a sawtooth crash, fast ion energy and angular momentum are modified resulting in the redistribution in phase space and orbit type change. The redistribution of fast ions in real space shows that the sawtooth instability brings different effect on fast particles with different orbit types as observed in experiments. Finally, the initial interpretative Transp simulation using the so-called kick model based on the Orbit modeling result shows an improvement of fast ion redistribution before and after a sawtooth crash but the neutron rate still has discrepancy compared to the experimental measurement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗