Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase field method”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Out of time order correlation of the Hubbard model with random local disorder

The out-of-time-order correlator (OTOC) serves as a powerful tool for investigating quantum information spreading and chaos in complex systems. We present a method employing non-equilibrium dynamical mean-field theory and coherent potential approximation combined with diagrammatic perturbation on the Schwinger–Keldysh contour to calculate the OTOC for correlated fermionic systems subjected to both random disorder and electron interaction. Furthermore, our key finding is that random disorder enhances the OTOC decay in the Hubbard model for the metallic phase in the weakly interacting limit. However, the current limitation of our perturbative solver restricts the applicability to weak interaction regimes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Phase Field Modeling of Corrosion for Design of Next-Generation Magnesium-Aluminum Vehicle Joints

In this study, WPI and its partners sought to understand galvanic corrosion in Mg-Al friction stir weld (FSW) joints. Unlike fusion welds, FSW is done in the solid state. This limits (but does not eliminate) the formation of brittle Al-Mg intermetallics, resulting in better joint strength and corrosion resistance. This study aimed to understand corrosion behavior of Mg-Al FSW joints by using a Cahn-Hilliard phase field electrochemistry model based on the work of Pongsaksawad et al. Its context challenge problem was an ultra-light door designed by Magna with 6061-Al sheet outer and ZEK100-Mg sheet inner panels, joined by an adhesive and FSW. Unlike most FSW joints, which weld through the softer material into the harder one, the Al sheet’s hem joint around Mg required the FSW tool to through the harder Al into the softer Mg, which was itself a challenge. Despite this challenge, PNNL succeeded in consistently making very strong FSW lap joints, by using a triflute tool and power control. Characterization showed mostly separate Al and Mg regions in the weld, usually with a hook protruding from Al into Mg, which likely contributed to high joint strength. Nanohardness mapping showed higher hardness in thin Mg swirls into the Al nugget, likely due to the intermetallics. The study used the industry standard SAE J2334 Cyclic Corrosion Test (CCT) to simulate corrosion conditions over the life of a vehicle, as well as linear polarization testing and ASTM G71 pitting corrosion testing. Corrosion reactions were: H 2 O + ½O 2 + 2e - → 2OH - at the Al cathode, and either Mg → Mg 2+ + 2e - or Mg + 2OH - → Mg(OH) 2 + 2e - at the Mg anode, with electron transfer through the joint. Mass loss increased with Al section length, supporting the hypothesis that galvanic corrosion in Al-Mg joints is limited by the cathodic reaction. Both pitting corrosion (first anode reaction) and hydroxide film formation (second reaction) were observed. The J2334 CCT test method was slightly modified, adding a chromium solution cleaning step after each week of corrosion testing, in order to remove Mg oxides and hydroxides and accurately measure mass loss over time. That said, FSW joints treated in this way did not did not exhibit significant reduction in lap shear strength of the joints vs. newly welded samples which had not undergone corrosion. Corrosion appeared to consist mainly of pitting in Mg in a way which did not directly affect the bond between the materials where the joints failed. A Cahn-Hilliard phase field model coupled with electrical potential described the electronically mediated galvanic corrosion electrochemical reactions described above in the four-component Al-Mg-H 2 O-(H 2 O+½O 2 ) system. This 2-D model successfully predicted Mg(OH)₂ formation at the anode in some circumstances, and Mg pitting corrosion in others, and showed the correct electric field directions in both cases. The model predicted Mg pit corrosion depth within a factor of two of measured pits.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Bulk-like first-order magnetoelastic transition in FeRh particles

Near-equiatomic, chemically-ordered iron-rhodium (FeRh) alloy is a fundamentally interesting material that may become useful in niche applications making use of its unique magneto functional phenomena, for example, the giant inverse magnetocaloric effect near room temperature that is associated with a sharp first-order magnetic phase transition. The nearly discontinuous antiferromagnetic-ferromagnetic phase transformation in bulk FeRh is well-known; however, the transition broadens considerably in fine particles and films with thickness less than 50 nm, precluding their potential applications. Here, we report an abrupt, bulk-like first-order magnetoelastic transformation in powders consisting of sub-micron particles of nearly equiatomic FeRh compound synthesized via solid-state mechanochemical co-reduction of FeF 2 and RhCl3 and subsequent heat treatments. In this work, we demonstrate that annealing at temperatures ranging from 600 °C to 800 °C enables tailoring phase content, particle size, and magnetic properties of the powders. A maximum magnetic-field-induced entropy change of ~10 J/kg K at μ0ΔH = 1 T has been achieved in powders annealed at 800 °C. The retention of extraordinary responsiveness in sub-micron particles of FeRh is likely to open doors for system component fabrication using additive manufacturing methods, along with new opportunities to employ FeRh in theranostics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Metastable Materials Discovery in the Age of Large-Scale Computation

Computational materials discovery has been successful in predicting novel, technologically relevant materials. However, it has remained focused almost exclusively on finding ground-state structures. Now that the lower-hanging fruit has been found in many fields of application, materials exploration is moving toward metastable materials: higher energy phases that are stable at practical time scales. Because of the challenges associated with predicting which phases are realistic, this class of materials has remained relatively unexplored, despite numerous examples of metastable structures with unmatched properties (e.g., diamond). This article highlights recent advances in developing computational and theoretical methods for predicting useful and realizable metastable materials. Topics discussed cover (1) the latest strategies for identifying potential metastable phases, (2) methodologies for assessing which phases can be realized experimentally, and (3) current approaches to estimate the lifetime of metastable materials.

amorphous materials↗
DOE OSTI2021

Eulerian simulation of complex suspensions and biolocomotion in three dimensions

Significance Fluid–structure interactions are ubiquitous in many natural and man-made environments. They are difficult to study analytically, and therefore accurate and flexible computational methods are an indispensable tool in the field. Typically, fluids are simulated with a fixed background computational mesh, whereas a solid is simulated with a mesh that moves with it, making it challenging to couple the two. Here we develop a three-dimensional computational method where both fluid and solid can be represented on a fixed computational grid, which simplifies the coupling between the two phases considerably. Our method can simulate scenarios that are difficult to do with other methods, such as complex suspensions containing soft particles that are both heavier and lighter than the fluid.

3D fluid-structure interaction↗
DOE OSTI2020

Magnetotransport and electronic structure of the antiferromagnetic semimetal YbAs

A number of rare-earth monopnictides have topologically nontrivial band structures together with magnetism and strong electronic correlations. In order to examine whether the antiferromagnetic (AFM) semimetal YbAs (TN=0.5 K) exhibits such a scenario, we have grown high-quality single crystals using a flux method, and characterized the magnetic properties and electronic structure using specific heat, magnetotransport, and angle-resolved photoemission spectroscopy (ARPES) measurements, together with density functional theory (DFT) calculations. Both ARPES and DFT calculations find no evidence for band inversions in YbAs, indicating a topologically trivial electronic structure. From low-temperature magnetotransport measurements, we map the field-temperature phase diagram, where we find the presence of a field stabilized phase distinct from the AFM phase at low temperatures. Additionally, an extremely large magnetoresistance (XMR) for both YbAs and the nonmagnetic counterpart LuAs is also observed, which can consistently be accounted for by the presence of electron-hole compensation. Moreover, an angle-dependent study of the Shubnikov-de Haas effect oscillations reveals very similar Fermi surfaces between YbAs and LuAs, with light effective masses down to at least 0.5 K, indicating that the Yb-4f electrons are well localized, and do not contribute to the Fermi surface. However, the influence of the localized Yb-4f electrons on the magnetotransport of YbAs can be discerned from the distinct temperature dependence of the XMR compared to that of LuAs, which we attribute to the influence of short-ranged spin correlations that appear well above TN.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

YbV 3 ⁢Sb 4 and EuV 3 ⁢Sb 4 vanadium-based kagome metals with Yb 2+ and Eu 2+ zigzag chains

Here, in this work, we present Yb⁢V 3 ⁢Sb 4 and Eu⁢V 3 ⁢Sb 4 , two compounds exhibiting slightly distorted vanadium-based kagome nets interleaved with zigzag chains of divalent Yb 2+ and Eu 2+ ions. Single crystal growth methods are reported alongside magnetic, electronic, and heat capacity measurements. Yb⁢V 3⁢ Sb 4 is a nonmagnetic metal with no collective phase transitions observed between 60 mK and 300 K. Conversely, Eu⁢V 3 ⁢Sb 4 is a magnetic kagome metal exhibiting easy-plane ferromagneticlike order below T C =32 K with hints of modulated spin texture under low field. With Yb⁢V 3 ⁢Sb 4 and Eu⁢V 3 ⁢Sb 4 we demonstrate another direction for the discovery and development of vanadium-based kagome metals while incorporating the chemical and magnetic degrees of freedom offered by a rare-earth sublattice.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Photovoltaic Analysis and Response Support (PARS) Platform for Solar Situational Awareness and Resiliency Services

The project's primary objective is to develop a digital-twin based Photovoltaic (PV) Analysis and Response Support (PARS) platform, which aims to provide real-time situational awareness and optimal response plans. This platform is designed to enhance the performance of hybrid PV systems, making them competitive with or even superior to conventional generation resources. The PARS platform enabled the project team to develop and evaluate an extensive suite of grid support functionalities for the hybrid PV systems to enhance grid performance, across key areas including visibility, dispatchability, security, resilience, and reliability. Given the global push toward achieving 100% clean energy by 2035, there is a significant increase in the integration of inverter-based resources (IBRs) throughout the energy grid. Effectively managing the inherent variability and uncertainty associated with IBRs is crucial for ensuring cost-effectiveness, reliability, and security in both the main grid and islanded microgrids. Constrained to a limited array of IEEE test systems or standard feeder models, traditional IBR modeling struggles to assimilate new field data, accurately reflect system dynamics, and adapt to the evolving energy landscape. In our project, we embraced a Digital Twin (DT) strategy for crafting the PARS platform. A digital twin acts as a precise virtual counterpart of a physical system, built on historical data and continuously honed with real-time insights. This enables the high-fidelity DT to accurately mirror current system operations and forecast future scenarios. Consequently, the PARS platform becomes an ideal environment for testing and refining monitoring, control, power, and energy management algorithms designed to boost hybrid PV system performance. The defining feature of the PARS platform, distinguishing it from other advanced simulation tools, is its exceptional adaptability. This is achieved by employing actual network topologies and utilizing real-time field data for fine-tuning and calibration, ensuring a close emulation of real-world conditions. The project deliverables include: 1) High-fidelity IBR models and tools for real-time parameterization, utilizing real-time field measurements to refine IBR models for enhanced accuracy and performance; 2) Grid-forming and Grid-following capabilities to deliver resilience services, including blackstart, voltage and frequency support, cold-load pick-up, power reserves, and three-phase load balancing across grid-connected and microgrid settings; 3) Machine learning-based forecasting tools and methods for generating synthetic data and topologies, creating diverse and realistic simulation environments for evaluating varied operational scenarios; 4) Advanced microgrid power and energy management algorithms for optimizing the integration and operation of PV, storage, and demand response resources within both feeder and community scales. The power grid data sets are provided by four utility companies in North Carolina and the New York Power Administration. Acting as industry advisors, our industry partners communicated stakeholder needs and regulatory standards to the research teams, aiding technology transfer by incorporating the developed methodologies into their daily operations. This collaboration ensures that the PARS platform, functioning as a power system digital twin, enhances our understanding of IBR dynamic behaviors and enables the development and evaluation of IBR control functions that match or exceed the capabilities of conventional synchronous generators.

14 SOLAR ENERGY↗
DOE OSTI2024

GENERATING CORCOS-COHERENT SIGNAL SERIES USING PHASE PERTURBATION​

Phase perturbations are used to simulate the time-domain behavior of pressure loads from turbulent boundary layers while faithfully representing the coherence decay behaviors of turbulence-generated pressure fields. Turbulence "evolves" as it proceeds downstream, and coherence decays as the decaying exponential of the product of a flow-related constant with frequency and distance. Time domain methods for a line, a plane, a closed cyclic surface, and most of a cone were given.

Swift, Stephen Hales↗
DOE OSTI2021

Fabrication of Low-Cost Large-Volume Ceramic A 2 HfX 6 (A= Cs or Tl, X = Cl, Br, or I) Scintillators for Gamma Ray Detection (SBIR Phase I Final Technical Report)

Scintillator crystals play an important role in the radiation detection field. Widespread use of scintillators as gamma-ray detectors is largely generated by their extensive availability and tunable properties, such high light output, high stopping power (Z eff ), fast decay time, and good proportionality. Additionally, the cost for manufacturing a scintillation detector like NaI:Tl is usually considerably lower than the cost for manufacturing a semiconductor detector like CdZnTe. Because there is no such thing as an ideal scintillation material, an application requiring certain detection characteristics may incorporate a scintillator tailored to its specific properties. The vast variety of applications and requirements necessitates more research into new scintillation materials and/or better methods of producing existing materials.The goal of this project was to grow low cost and environmentally stable inorganic transparent ceramic scintillators with excellent gamma ray resolution, excellent energy proportionality, excellent detection efficiency due to high density (>5 g/cm 3 ) and very high Z eff (55-80), and good light yields (>40,000 ph/MeV). In Phase I Xtallized Intelligence, Inc. (XI, Inc) developed a novel ceramic fabrication technique to produce low cost and environmentally stable highly efficient inorganic transparent ceramic scintillators of various dimensions. XI, Inc., collaborating with Fisk University (Fisk), investigated the scintillation properties of these new ceramic scintillators and compared them to in their single crystal counterparts. The results of this Phase I project show that successful production of high-quality inorganic halide ceramic scintillators Cs 2 HfCl 6 (CHC) and Tl 2 HfCl 6 (THC). Both ceramic CHC and THC scintillators have achieved good performance close to the performance of their single crystal counterparts. Fabricating these inorganic ceramic scintillators mitigate many issues encountered during conventional bulk crystal growth by melt methods. Additional benefits of the ceramic fabrication technique include high production yield, low production cost, fast production time, and no material waste Inorganic transparent ceramic scintillators produced in this project will enhance cost effectiveness at the instrument level based on low projected cost of the proposed compounds, as much smaller crystal sizes would be required to achieve similar efficiency as current radioisotope identification devices (RIID’s) used in homeland security applications as well as spectrometers in high energy physics applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Machine-learning force-field models for dynamical simulations of metallic magnets

We review recent advances in machine-learning (ML) force-field methods for Landau–Lifshitz–Gilbert simulations of itinerant electron magnets, focusing on their scalability and transferability. Built on the principle of locality, a deep neural-network model is developed to efficiently and accurately predict electron-mediated forces governing spin dynamics. Symmetry-aware descriptors constructed through a group-theoretical approach ensure rigorous incorporation of both lattice and spin-rotation symmetries. The framework is demonstrated using the prototypical s-d exchange model widely employed in spintronics. ML-enabled large-scale simulations reveal novel nonequilibrium phenomena, including anomalous coarsening of tetrahedral spin order on the triangular lattice and the freezing of phase-separation dynamics in lightly hole-doped, strong-coupling square-lattice systems. These results establish ML force-field frameworks as scalable, accurate, and versatile tools for modeling nonequilibrium spin dynamics in itinerant magnets.

Artificial neural networks↗
DOE OSTI2020

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗
DOE OSTI2024

Extraction of Drell-Yan Angular Parameters in $pp$ Collisions with a 120 GeV Beam Energy Using a Deep-Learning Unfolding Algorithm

Dilepton production in pp collisions through the Drell-Yan process provides a crucial tool for studying the internal quark-gluon structure of the nucleon. By precisely measuring the $\cos2\phi$ asymmetry, where $\phi$ represents the azimuthal angle of the $l^{+}l^{-}$ pair in the Collins-Soper frame, we can gain valuable insights into the proton’s structure and the transverse momentum ($q_{T}$) dependence of the $\cos2\phi$ asymmetry. SeaQuest, a fixed-target Drell-Yan experiment at Fermilab, involved an unpolarized proton beam colliding with unpolarized LH$_{2}$ and LD$_{2}$ targets. Measurements obtained from experiments typically require corrections for detector inefficiencies, smearing, and acceptance. Traditionally, these corrections involve “unfolding” the detector-level measurements through matrix operations. However, in higher-dimensional phase space, these conventional methods fail to scale effectively. To overcome these limitations, we employ an unbinned unfolding method that utilizes deep neural networks for unfolding higher-dimensional phase space. In this presentation, we will explain the design of the neural network architecture, our training strategies, and outline our plans to achieve conclusive results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Predictive Chemical Kinetic Modeling: Where We Succeed, Where We Struggle, and What Comes Next

Chemical kinetic modeling plays a foundational role in fields ranging from energy to environmental science, pharmaceuticals, and advanced materials. The past two decades have seen remarkable progress, particularly in modeling gas-phase reactions for thermochemical processes, leading to impactful industrial applications such as steam cracking and air quality management. However, new challenges are emerging. The successful development of systematic methodologies for the description of gas-phase kinetics opens the possibility to apply the same approach to the study of more challenging systems. Here, we review recent advances, including ab initio transition state theory-based master equation estimation of elementary rates, automated mechanism generation, machine-learning-assisted kinetics, and uncertainty quantification, and discuss the advances needed to apply the same methodological approach in areas such as heterogeneous catalysis, electrochemistry, liquid-phase and solid-state reactivity, and multiscale model integration. We advocate for the development of targeted tools, especially methods that go beyond empirical tuning toward first-principles-based predictions. We highlight the need for accessible software and AIaugmented workflows to democratize modeling for industry and academia alike. In this perspective, we call attention to not only what has worked but also what remains unsolved, advocating to avoid overemphasizing successes in scientific works at the expense of realism. The next decade should focus on predictive capability, physical accuracy, and community infrastructure (e.g., databases and services) to enable innovation across diverse fields. We argue that kinetic modeling, properly equipped, can accelerate discovery far beyond its traditional domains.

ab initio calculations↗
DOE OSTI2025

Dynamic response of 17-4 stainless steel as a function of manufacturing method and heat treatment

We present a series of plate-impact experiments on 17-4 stainless steel to study the effect of manufacturing method and heat treatment on the Hugoniot elastic limit (HEL), Hugoniot, phase transformation stress, and spallation strength. Two traditional manufacturing methods were considered, wrought processing and casting, as well as two additive manufacturing methods, laser powder-bed fusion (LPBF) and wire-fed electron beam (EBAM). For both LPBF and EBAM 17-4 stainless steel variants, two billets were printed, enabling the application of two unique heat treatments. The HEL stress depended heavily on the thermal history, with the HEL increasing after the formation of Cu-rich precipitates via heat treatment. The Hugoniot response both below and above the phase transition was unaffected by the manufacturing method or heat treatment. The phase transition stress depended heavily on the thermal history, with its variation being attributed to the presence of various microstructural features. This is supported by a marked increase in the phase transition stress after precipitation hardening. These results suggest that the notion of the phase transition stress being dictated by bulk composition is an oversimplification and the stress fields generated by the meso-scale structure are a dominant force. The spallation strength was lower in the cast material compared to all other 17-4 stainless steel variants due to the presence of brittle δ-ferrite inclusions. Additionally, a drop in the tensile strain-rate was observed in the spallation response above the phase transition stress, which was hypothesized to stem from the kinetics of the reversion to the low-pressure phase during spall.

Compressive stress↗
DOE OSTI2023

Spin Squeezing by Rydberg Dressing in an Array of Atomic Ensembles

Here, we report on the creation of an array of spin-squeezed ensembles of cesium atoms via Rydberg dressing, a technique that offers optical control over local interactions between neutral atoms. We optimize the coherence of the interactions by a stroboscopic dressing sequence that suppresses super-Poissonian loss. We thereby prepare squeezed states of N = 200 atoms with a metrological squeezing parameter ξ 2 = 0.77⁢(9) quantifying the reduction in phase variance below the standard quantum limit. We realize metrological gain across three spatially separated ensembles in parallel, with the strength of squeezing controlled by the local intensity of the dressing light. Our method can be applied to enhance the precision of tests of fundamental physics based on arrays of atomic clocks and to enable quantum-enhanced imaging of electromagnetic fields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Validation and Verification for INL Modelica-based TEDS models Via Experimental Results

This report provides an overview on the verification and validation (V&V) of the Thermal Energy Distribution System (TEDS) model developed in the Modelica process modeling ecosystem using experimental data. Model development has led to the creation of a dynamic process model of the experimental TEDS facility housed within the Energy Systems Laboratory (ESL) at Idaho National Laboratory (INL). The model was then used during the preconstruction phase of the experimental effort to inform experimental design (e.g., insulation requirements, bypass line placement, expected performance of components) and to test innovative control schemes prior to the initial operation. The TEDS model developed in Modelica includes the primary components of the TEDS experimental unit: a 200kW Chromalox heater; a single-tank packed-bed thermal energy storage system filled with 0.125-inch alumina (Al2O3) beads; an ethylene-glycol-to-Therminol-66 heat exchanger; system piping; five control valves; and all associated temperature, pressure, and volumetric flow sensors. Using the Institute of Electrical and Electronics Engineers (IEEE) V&V methodologies, considered the gold standard in the engineering field, the model was verified using a combination of static analysis, spatial convergence, and regression tests. Then using dynamic time warping (DTW) initial runs to validate and tune the TEDS model versus the experiment were conducted. This tuning method was accomplished using the INL Risk Analysis Virtual ENvironment (RAVEN) software package. Tuning is required to account for physical phenomena that are less understood within the empirical heat transfer correlations. Through the commencement of this work, a systems-level model of TEDS with associated control systems, sensors, piping diameters, and component capabilities has been created. This model was utilized in the pre-experimental phase to inform system design, insulation thicknesses, and potential control schemes to operate the system effectively and safely. Then, initial experimental startup and operational data were used to demonstrate the validation and tuning methodology. This process demonstrates the classical two-step approach of a model informing experimental design followed by the experiment validation and tuning the model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗