Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulator Study”

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 703 records · Page 39

Ground state property calculations of LiH n complexes using IBM Qiskit’s quantum simulator

In this study, the variational quantum eigensolver (VQE) on a quantum simulator is used in calculating ground state electronic structure properties of the LiH n , n = 1–3, complexes including their singly charged ions. Results calculated using classical electronic structure algorithms are also included. We investigate the use of the unitary coupled cluster with singles and doubles (UCCSD) Ansatz using VQE within Qiskit and compare results to full configuration interaction (FCI) calculations. Computed ground state energies, electron affinities, ionization potentials, and dipole moments are considered. We report the first-of-its-kind simulated quantum computing results of selected LiH n species and use the parity orbital to qubit mapping scheme. We find that VQE/UCCSD results are comparable to classical coupled clusters with singles and doubles for all considered systems with respect to FCI. A VQE calculation cost evaluation is included in which we evaluate performance using both Jordan–Wigner and parity orbital to qubit mapping schemes. We also discuss some of the current limitations of utilizing VQE for the study of chemical systems.

97 MATHEMATICS AND COMPUTING↗

A Review of Tank 48H Treatment of Tetraphenylborate with Permanganate

Tank 48H contains roughly 270,000 gallons of radioactive waste material. The waste stored in this tank was to be processed in multiple stages utilizing facilities at the Savannah River Site (SRS) almost 30 years ago. The In-Tank Precipitation Process (ITP) was initiated in Tank 48H, which precipitated highly radioactive cesium-137 using sodium tetraphenylborate (NaTPB). While the process succeeded in precipitating the Cs, Sr, and other actinides, it also generated an unexpectedly large amount of benzene. The evolved benzene created a safety concern and made the waste incompatible with further downstream processing. This halted the ITP, and a new method of treatment was required to continue processing the legacy waste contained in Tank 48H. A myriad of treatment options have been proposed with multiple teams of researchers assembled to work on this highly complex issue over the last few decades. This review investigated one potentially viable treatment option, in-tank oxidation with sodium permanganate. Permanganate has been well known in literature as a strong oxidizing agent for organic compounds, as well as being utilized at the Savannah River Site (SRS) in other processes. A small number of studies have been conducted utilizing waste simulants to evaluate the use of permanganate as an oxidant for tetraphenylborate (TPB). A search of the literature and data from these studies indicates that permanganate could be a viable treatment for destruction of TPB in Tank 48H. While the scoping studies had a small number of individual experiments and nonideal conditions, the permanganate decomposed up to 90% of the TPB. A free hydroxide concentration above 1.0 M is required for tank corrosion control. Simulant studies indicate no decrease in TPB decomposition by permanganate until pH 14. The simulant studies show an increase in TPB decomposition as the temperature of the solution is increased to 40 °C. The post-reaction analysis of previous simulant tests did not look at all of the organic degradation products. Investigations into what these organic products are and in what quantity will help guide determinations as to whether the downstream processing facilities are able to handle the material that will be generated. Study on the time frame for the reaction between permanganate and TPB should be investigated as the literature reports only extend out to two weeks reaction time. The permanganate treatment conditions indicate no corrosion control concerns and a longer timescale reaction may be needed for in-tank treatment. In addition, further study would be useful to identify a lower boundary condition for the ratio of TPB and oxidant. The simulant tests applied large excesses of permanganate, and this may be unnecessary. The size constraint of the tank means that there will be practical limitations on the amount of sodium permanganate that can be added to the tank. The amount of permanganate should be minimized as much as possible while still ensuring decomposition of the TPB to minimize the amount of manganese dioxide solids generated.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Dark Energy Studies with LSST using the Photon Simulator (Final Report)

This grant funded the development and dissemination of the Photon Simulator (PhoSim) for the purpose of studying dark energy at high precision with the upcoming Large Synoptic Survey Telescope (LSST) astronomical survey. The work was in collaboration with the Dark Energy Science Collaboration (DESC). Several detailed physics improvements were made in the optics, atmosphere and sensor, a number of validation studies were performed, and a significant number of usability features were implemented. Future work in DESC will use PhoSim as the image simulation tool for data challenges used by the analysis groups.

79 ASTRONOMY AND ASTROPHYSICS↗

Effect of surface properties and polymer chain length on polymer adsorption in solution

In polymer nanoparticle composites (PNCs) with attractive interactions between nanoparticles (NPs) and polymers, a bound layer of the polymer forms on the NP surface, with significant effects on the macroscopic properties of the PNCs. The adsorption and wetting behaviors of polymer solutions in the presence of a solid surface are critical to the fabrication process of PNCs. In this study, we use both classical density functional theory (cDFT) and molecular dynamics (MD) simulations to study dilute and semi-dilute solutions of short polymer chains near a solid surface. Using cDFT, we calculate the equilibrium properties of polymer solutions near a flat surface while varying the solvent quality, surface–fluid interactions, and the polymer chain lengths to investigate their effects on the polymer adsorption and wetting transitions. Using MD simulations, we simulate polymer solutions near solid surfaces with three different curvatures (a flat surface and NPs with two radii) to study the static conformation of the polymer bound layer near the surface and the dynamic chain adsorption process. Additionally, we find that the bulk polymer concentration at which the wetting transition in the poor solvent system occurs is not affected by the difference in surface–fluid interactions; however, a threshold value of surface–fluid interaction is needed to observe the wetting transition. We also find that with good solvent, increasing the chain length or the difference in the surface–polymer interaction relative to the surface–solvent interaction increases the surface coverage of polymer segments and independent chains for all surface curvatures. Finally, we demonstrate that the polymer segmental adsorption times are heavily influenced only by the surface–fluid interactions, although polymers desorb more quickly from highly curved surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Characterization of partial wetting by CMAS droplets using multiphase many-body dissipative particle dynamics and data-driven discovery based on PINNs

The molten sand that is a mixture of calcia, magnesia, alumina and silicate, known as CMAS, is characterized by its high viscosity, density and surface tension. The unique properties of CMAS make it a challenging material to deal with in high-temperature applications, requiring innovative solutions and materials to prevent its buildup and damage to critical equipment. Here, we use multiphase many-body dissipative particle dynamics simulations to study the wetting dynamics of highly viscous molten CMAS droplets. The simulations are performed in three dimensions, with varying initial droplet sizes and equilibrium contact angles. We propose a parametric ordinary differential equation (ODE) that captures the spreading radius behaviour of the CMAS droplets. The ODE parameters are then identified based on the physics-informed neural network (PINN) framework. Subsequently, the closed-form dependency of parameter values found by the PINN on the initial radii and contact angles are given using symbolic regression. Finally, we employ Bayesian PINNs (B-PINNs) to assess and quantify the uncertainty associated with the discovered parameters. In brief, this study provides insight into spreading dynamics of CMAS droplets by fusing simple parametric ODE modelling and state-of-the-art machine-learning techniques.

Mechanics↗

Identification of Collectible Items in the Rancor Microworld Simulator Compared to Full-scope Studies

Most studies collecting reliability data in human reliability analysis (HRA) have concentrated on research using full-scope simulators and actual operators. However, several researchers have found that there are challenges to collecting the variety of items as well as the amount of data needed to support HRA. As an opposite and complementing approach to the full-scope study, Idaho National Laboratory (INL) has developed a simplified simulator, i.e., Rancor Microworld, for generating HRA data based on student subjects. This paper aims to identify collectible items in Rancor Microworld versus full-scope simulators. In this paper, collectible items and their different analysis levels in Rancor Microworld are identified in comparison with full-scope simulators. Then, how to treat the items is suggested with experiment directions that will be treated in the future.

99 GENERAL AND MISCELLANEOUS↗

Response of extreme precipitation to uniform surface warming in quasi-global aquaplanet simulations at high resolution

Projections of precipitation extremes in simulations with global climate models are very uncertain in the tropics, in part because of the use of parameterizations of deep convection and model deficiencies in simulating convective organization. Here, we analyse precipitation extremes in high-resolution simulations that are run without a convective parameterization on a quasi-global aquaplanet. The frequency distributions of precipitation rates and precipitation cluster sizes in the tropics of a control simulation are similar to the observed distributions. In response to climate warming, 3 h precipitation extremes increase at rates of up to 9 % K - 1 in the tropics because of a combination of positive thermodynamic and dynamic contributions. The dynamic contribution at different latitudes is connected to the vertical structure of warming using a moist static stability. When the precipitation rates are first averaged to a daily timescale and coarse-grained to a typical global climate-model resolution prior to calculating the precipitation extremes, the response of the precipitation extremes to warming becomes more similar to what was found previously in coarse-resolution aquaplanet studies. However, the simulations studied here do not exhibit the high rates of increase of tropical precipitation extremes found in projections with some global climate models. This article is part of a discussion meeting issue ‘Intensification of short-duration rainfall extremes and implications for flash flood risks’.

Science & Technology - Other Topics↗

A Simulation Framework for Analyzing the Impact of Stochastic Occupant Behaviors on Demand Flexibility in Typical Commercial Buildings

As one of the primary users of the electric grid, buildings and building equipment, including heating, ventilation, and air conditioning (HVAC) systems, can be leveraged to provide the flexible demand needed to balance the grid. Typical strategies to achieve demand flexibility are to reduce electricity use during peak or critical periods by shutting down equipment or relaxing system setpoints, which will inevitably impact the occupants’ comfort. When occupants feel uncomfortable, they may take actions to regain their comfort, and some of those actions (such as turning on a personal fan) may have a negative impact on meeting the demand response goal. Therefore, it is important to incorporate occupant behaviors into the assessment ofthe building demand flexibility potential. In this study, a simulation framework that includes simulation of zone thermal loads, an HVAC system, and occupant behaviors, was developed to investigate the impact of occupant behaviors on demand flexibility. A case study was conducted using a small office model from the U.S. Department of Energy (DOE) Commercial Prototype Building Models to simulate the building envelope and zone loads. An agent-based occupant thermal behavior model was adapted to forecast occupants’ thermal comfort and their resulting thermal behaviors. An artificial neural network (ANN) based airflow model trained from a computational fluid dynamics (CFD) model of the zone was adopted to better predict the ambient environment of each occupant. An air-source heat pump simulation model that was calibrated from a real two-stage air-source heat pump system was used as the HVAC system. A typical load shedding event during peak hours was studied. Repeated simulations were conducted to capture the stochastic effects of occupant behaviors. The interplay between the demand flexibility, occupant comfort and behavior were analyzed by evaluating key performance indicators, including the energy use, occupant discomfort duration, and occupant behavior duration during the peak period. The results suggest that this framework can be used to analyze typical commercial buildings and their HVAC systems in terms ofdemand flexibility potential under the impact of occupant behaviors.

Chen, Zhelun↗

Evaluating High-Halide Waste Form Options for Salt-Based Nuclear Waste Simulants

In this study, waste forms were being explored for an electrochemical salt simulant, , referred to as ERV3, which is a high-LiCl/KCl salt containing simulated fission products (i.e., Sr, Cs, Nd) and Na to represent bond sodium from Experimental Breeder Reactor-II metallic fast reactor fuel. The goal was to find glassy systems that could be used to immobilize the salt in a single-step process. The envisionment of this process would be find a frit glass that could be added to the salt waste, heat treated, poured into waste canisters, and then stored for disposal. To perform this study, a literature review was conducted, the most promising seven systems were fabricated without the salt, and then mixed with salt simulant and heat treated under different processes. The criteria that were used to screen potential compositions included demonstrated alkali incorporation, could be melted at reasonably low temperature ( T ≤ 1000°C), and if the compositions had some demonstrated data for waste-form-related properties, that was a benefit. High marks were given for compositions that showed amorphous nature after slow cooling of samples containing ERV3.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Damage to Polystyrene Polymer Film by Shock Wave Induced Bubble Collapse

Metallic surfaces that are in contact with solutions are commonly used in numerous applications where these surfaces can be damaged by shock wave induced bubble collapse. Use of polymer films that coat such surfaces to prevent them from damage requires a better understanding of how much harm collapsing bubbles produce in the films. In this study, we report the results from coarse-grained molecular dynamics simulations to study the damage to polystyrene (PS) films coating a hard surface. The damage was caused by a collapsing nanobubble located in the proximity of the film and interacting with an impinging shock wave. This collapse produces a high-speed water jet that impacts the PS film with a greater force than the shock front and creates cavities/pits in the PS film. We observed that polymer molecules located in the jet vicinity undergo conformational extension in the direction perpendicular to the jet motion, while chain molecules in the rest of the film undergo compression. Furthermore, we also observed that damage to the film is sensitive to the strength of the shock wave.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulation of femtosecond interferometry when studying the effect of intense laser irradiation of a transparent medium

A method for comparing experimental and simulation results obtained when transmitting a femtosecond laser pulse through a material is proposed. This method is based on computer reproduction of photographs recorded using femtosecond interferometry. A numerical analysis of the spatial and temporal evolution of an intense (∼10{sup 13} W cm{sup −2}) pulse propagating in the bulk of yttrium aluminium garnet is presented. (paper)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A simplified CFD approach for modeling mass transport in catalytic open-cell foams

A simplified macroscopic CFD approach is presented to model mass transport including chemical reactions in washcoated open-cell foams. The foam is treated as a porous medium. Species conversion during chemical reactions is modeled using appropriate source terms based on reaction rate expressions and modified to account for the mass transport resistances occurring at the fluid-washcoat interfaces and within the washcoat layers. As example, the catalytic CO oxidation over platinum is studied. The simulation results show good agreement with experimental data from literature. A parametric study on washcoat parameters, such as thickness, tortuosity, porosity, and size, is carried out. Increasing the washcoat thickness from 5 to 100 µm or decreasing the tortuosity to porosity ratio from 5 to 20 decreases the CO conversion by 10 %. The proposed model is found to be reliable and has the advantage of lower computational cost, making it a suitable tool for foam-based catalytic reactor design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluation of Acoustic Channel Capacity for Complex Piping Topology

Transmission of information using elastic waves on existing metallic pipes provides an alternative communication option for a nuclear facility. The advantages of this approach consist of transmitting information through physical barriers, such as the containment building wall. In addition, this approach provides physical cybersecurity since the pipe becomes an analog communication channel which is difficult to eavesdrop on and difficult to sever. In prior work, a laboratory bench-scale system consisting of a six foot-long straight section of nuclear grade stainless steel pipe and high-temperature ultrasonic LiNbO 3 transducers was assembled. The lab systems was used for communication protocol development and preliminary data transmission demonstration. Next phase of project development involves deployment of the ultrasonic communication system in a representative environment. Mechanisms Engineering Test Loop (METL) facility at ANL was chosen as a viable candidate for communication system demonstration. Analysis of piping manifolds at METL has shown that in a typical scenario, ultrasonic signals have to be transmitted over straight and bent piping sections. This report describes preliminary studies on transmission of information with ultrasonic shear waves across bent piping sections. A test article consisting of a stainless steel pipe bent at 90° was developed for laboratory analysis of ultrasonic signal transmission. The bent piping test article developed by welding two straight pipes to an elbow. For consistency, the diameter and wall thickness of the bent pipe is the same as that of the straight pipe utilized in prior studies. COMSOL computer simulations were performed to study ultrasonic refracted shear wave coupling and transmission across bent piping. Preliminary evaluations of ultrasonic shear wave transmission were conducted with piezoelectric (PZT) and LiNbO 3 transducers, with the signals recorded with previously developed LabView interface. Preliminary results have shown significant dispersion of transmitted signals, most likely due to reflections and scattering by piping welds. Signal distortion compensation algorithm based on time reversal modulation (TRM) was initially investigated for transmission through solids, and demonstrated for signals transmitted over the bent pipe. Using TRM, it was shown that communication bitrates of 100Kbps were achievable, in principle.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Final Report for Transmission of Information by Acoustic Communication along Metal Pathways in Nuclear Facilities

Transmission of information using elastic waves on existing metallic pipes provides an alternative communication option for a nuclear facility. The advantages of this approach consist of transmitting information through barriers, such as the containment building wall. A viable candidate for acoustic communication channel is a chemical volume control system (CVCS) stainless steel pipe, which penetrates through the containment building wall. A laboratory bench-scale system consisting of a nuclear grade CVCS-like pipe and ultrasonic transducers was assembled for a preliminary communication system analysis. Because of low bandwidth and spectral dispersion of ultrasonic transducers, on off keying (OOK) protocol was chosen for data communication. Operation of the acoustic communication during normal and post-accident conditions involves data transmission on the CVCS pipe, which according to literature guidelines, will be at 50°C and 150°C, respectively. Preliminary communication system development was performed using paintbrush PZT transducers capable of continuous operation at temperatures up to 80°C. Next, high-temperature LiNbO 3 acoustic transducers, originally developed for flow metering applications of EBRII reactor at Argonne, were utilized to develop a high-temperature acoustic communication system. A Gaussian pulse shape filter was introduced to suppress ringing and thus reduce inter-symbol interference. This resulted in enhancement of data transmission bitrate. Heating tapes, temperature controllers, and thermal insulation were installed on the laboratory pipe to study acoustic communication at elevated temperature. This included development of custom insertion heating elements for controlling temperature at pipe ends. Demonstrations of communication at high temperature included transmission of images and text files. Main achievements thus far include demonstration of transmission of 90KB image at the bitrate of 10Kbps with bit error rate smaller than 10 -3 across the pipe heated to 50°C and 150°C. Final phase of project development involved investigating strategies for deployment of the ultrasonic communication system in a representative environment. Mechanisms Engineering Test Loop (METL) facility at ANL was chosen as a viable candidate for communication system demonstration. Analysis of piping manifolds at METL has shown that in a typical scenario, ultrasonic signals have to be transmitted over straight and bent piping sections. This report describes preliminary studies on transmission of information with ultrasonic shear waves across bent piping sections. A test article consisting of a stainless steel pipe bent at 90° was developed for laboratory analysis of ultrasonic signal transmission. The bent piping test article developed by welding two straight pipes to an elbow. For consistency, the diameter and wall thickness of the bent pipe is the same as that of the straight pipe utilized in prior studies. COMSOL computer simulations were performed to study ultrasonic refracted shear wave coupling and transmission across bent piping. Preliminary evaluations of ultrasonic shear wave transmission were conducted with piezoelectric (PZT) and LiNbO 3 transducers, with the signals recorded with previously developed LabView interface. Preliminary results have shown significant dispersion of transmitted signals, most likely due to reflections and scattering by piping welds. Signal distortion compensation algorithm based on time reversal modulation (TRM) was initially investigated for transmission through solids, and demonstrated for signals transmitted over the bent pipe. Using TRM, it was shown that communication bitrates of 100Kbps were achievable, in principle.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Simulating the CMS High Granularity Calorimeter with ML

Detector simulation is a key component of physics analysis and related activities in CMS. In the upcoming High Luminosity LHC era, simulation will be required to use a smaller fraction of computing in order to satisfy resource constraints. At the same time, CMS will be upgraded with the new High Granularity Calorimeter (HGCal), which requires significantly more resources to simulate than the existing CMS calorimeters. This computing challenge motivates the use of generative machine learning models as surrogates to replace full physics-based simulation. We study the application of state-of-the-art diffusion models to simulate particle showers in the CMS HGCal. We will discuss methods to overcome the challenges posed by the high-dimensional, irregular geometry of the HGCal. The quality of the showers produced by the diffusion model will be assessed by comparison to the full GEANT4-based simulation. The increase in simulation throughput will be quantified and methods to accelerate the diffusion model inference will also be discussed.

Amram, Oz↗

Multi-Timescale Integrated Dynamics and Scheduling for Solar (MIDAS-Solar) (Final Technical Report)

Solar photovoltaic (PV) installations have experienced unprecedented growth in the United States. PV will become not only an energy producer but also a necessary provider of ancillary services at multiple timescales. Conventional methods to simulate power systems operations - such as long-term production simulation (which typically considers schedules from hours to minutes by using an optimization framework) and short-term transient studies (which simulate dynamics from seconds to sub-seconds using state variables and differential equations) - are not sufficient for studying the multiple-timescale variation of solar generation and its impact on system reliability. Long-term system economics and short-term system dynamics are highly coupled, particularly when the penetration level of renewable generation is extremely high, because the uncertainty and variability of solar generation will impact both power system steady-state and dynamic performance. This project helps meet and exceed the U.S. Department of Energy Office of Energy Efficiency and Renewable Energy Solar Energy Technologies Office goal of systems integration by directly addressing this stability and reliability challenge for power grid planning and operation. We have developed a temporally comprehensive, closed-loop simulation model, named Multi-timescale Integrated Dynamics and Scheduling (MIDAS), that seamlessly simulates power system operations from economic scheduling (day-ahead to hours) to dynamic response analysis (seconds to sub-seconds). For schedules with very high levels of inverter-based resources (IBRs), up to and including 100%, the stability of grid controls has been evaluated through electromagnetic transient (EMT) simulations and power-hardware-in-the-loop (PHIL) simulations of key transient events at key schedule points. Specifically, MIDAS provides: 1) a closed-loop simulation framework for simulating timescales from economic scheduling to dynamic stability analysis; 2) machine learning-based stability assessment; 3) EMT modeling and analysis for large-scale power systems; 4) MIDAS PHIL test bed. We worked with Hawaii Electric Companies to apply the MIDAS study framework to a Maui grid study. The entire island's transmission system was modeled in detail - from a yearly scheduling model, to a second-level frequency dynamic model, down to a sub-second-scale EMT model to address critical stability issues. The project demonstrated how MIDAS can help system planners and operators assess system reliability and stability while the power grid is marching toward a high-renewable, high-IBR future. In this Maui grid study, we found that 100% instantaneous IBR operation is achievable in EMT simulation and PHIL testing, and grid planners and operators might need new analysis/simulation tools to assess grid reliability and stability in the scheduling stage. MIDAS will bring Maui and other systems closer to 100% clean and stable energy futures. (In this study, we examined transient stability. Other topics necessary for 100% IBR operation, such as protection and resource adequacy, were not examined.)

100% Renewables↗

Progress on Understanding Rayleigh–Taylor Flow and Mixing Using Synergy Between Simulation, Modeling, and Experiment

Simultaneous advances in numerical methods and computing, theoretical techniques, and experimental diagnostics have all led independently to better understanding of Rayleigh–Taylor (RT) instability, turbulence, and mixing. In particular, experiments have provided significant motivation for many simulation and modeling studies, as well as validation data. Numerical simulations have also provided data that is not currently measurable or very difficult to measure accurately in RT unstable flows. Thus, simulations have also motivated new measurements in this class of buoyancy-driven flows. This overview discusses simulation and modeling studies synergistic with experiments and examples of how experiments have motivated simulations and models of RT instability, flow, and mixing. First, a brief summary of measured experimental and calculated simulation quantities, of experimental approaches, and of issues and challenges in the simulation and modeling of RT experiments is presented. Implicit large-eddy, direct numerical, and large-eddy simulations validated using RT experimental data are then discussed. This is followed by a discussion of modeling using analytical, modal, buoyancy–drag, and turbulent transport models of RT mixing experiments. The discussion will focus on three-dimensional RT mixing arising from multimode perturbations. Finally, this focused review concludes with a perspective on future simulation, modeling, and experimental directions for further research. Research in simulation and modeling of RT unstable flows, coupled with experiments, has made significant progress over the past several decades. This overview serves as an opportunity to both discuss progress and to stimulate future research on simulation and modeling of this unique class of hydrodynamically unstable turbulent flows.

42 ENGINEERING↗