Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time domain modeling”

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 397 records · Page 22

NRC Multiphysics Analysis Capability Deployment (FY2020--Part 3)

This report details progress and activities of Idaho National Laboratory (INL) on the NRC project ”Development and Modeling Support for Advanced Non-Light Water Reactors.” The tasks completed for this report are as follows: First, Task 2c: Explicit modeling of pebble transient temperature response. In this simulation, the PBMR-400 reactor experiences a 20 second power ramp from 100% to 150% power. This is followed by a similar reduction in the power back to 100%. Several multiscale pebble coupling approaches are tested with one pebble per mesh element in the active core region. The results show good conservation behavior and the stability of the coupling. Next, Extended scope part 1: An assessment of the computational efficiency of the DFEM heat transfer solver shows good scalability. The DFEM solver is a factor of 4 more expensive in solution time than the FEM solver for heat transfer problems due to the increased number of degrees of freedom. Nonetheless, the DFEM approach provides the user with the flexibility to model gap heat transfer problems. Then Extended scope part 2: the GapHeatTransferInterfaceMaterial was improved to give the user increased flexibility with the modeling of heat transfer through gaps with the DFEM solver. A number of gap parameters can now be coupled both through functions and variables. Finally, Extended scope part 3: demonstration of how the gap width between hexagonal fuel cells can be calculated during a heat up transient and used in the GapHeat- TransferInterface model. A full domain DFEM model with gap expansion is coupled to a SubApp that models the thermal expansion of the base plate. The results show the expected physical behavior, although have not been fully bench- marked at this point in time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Quantifying the Role of Ocean Dynamics in SST Variability across GCMs and Observations

Abstract Midlatitude SSTs forced by mesoscale oceanic processes can affect the large-scale atmosphere, pointing to the ocean’s crucial role outside the tropics. Previous studies have shown oceanic mesoscale processes’ effect on global and regional climate variability. This study quantifies the local contribution of ocean dynamics to mixed-layer temperature across the globe by directly estimating the ocean heat flux divergence resolved by state-of-the-art ocean reanalysis, eddy-resolving, and eddy-parameterized versions of two U.S. national climate models and indirectly from air–sea flux satellite-based estimates. Our results show that the eddy-resolving climate simulations resolve mixed-layer temperature variances that are larger and closer to those inferred from observations than both their eddy-parameterized counterparts and ECCO over much of the extratropics. The observations and the eddy-resolving models indicate a more significant role of ocean dynamics in the mixed-layer temperature variability than the surface fluxes over most extratropics compared to their eddy-parameterized versions. A frequency domain analysis shows that the better-resolved ocean mesoscale and thermal gradients enhance the variance over a time scale from 2 months to 30 years. Results show agreement in the ocean’s contribution among satellite-based estimates, ocean reanalysis products, and ocean eddy-resolving simulations. At the same time, differences emerge for ECCO and the eddy-parameterized models, suggesting that surface fluxes account for a larger fraction of the mixed-layer temperature variability in most of the extratropics.

Siqueira, L.↗

JACC: Leveraging HPC Meta-Programming and Performance Portability with the Just-in-Time and LLVM-based Julia Language

We present JACC (Julia for Accelerators), the first high-level, and performance-portable model for the just-in-time and LLVM-based Julia language. JACC provides a unified and lightweight front end across different back ends available in Julia, enabling the same Julia code to run efficiently on many HPC CPU and GPU targets. We evaluated the performance of JACC for common HPC kernels as well as for the most computationally demanding kernels used in applications, HPCCG, a supercomputing benchmark test for sparse domains, and HARVEY, a blood flow simulator to assist in the diagnosis and treatment of patients suffering from vascular diseases. We carried out the performance analysis on the most advanced US DOE supercomputers: Aurora, Frontier, and Perlmutter. Overall, we show that JACC has a negligible overhead versus vendor-specific solutions, reporting GPU speedups with no extra cost to programmability.

Valero-Lara, Pedro↗

Thermal coupling-decoupling mechanism of heat transfer across van der Waals interfaces in n-eicosane

Progress toward improving thermal transport properties of n-eicosane, a promising phase change material for thermal energy storage near human body temperature, is hindered by lack of a general theoretical framework to model thermal interfacial conductance (TIC) of the van der Waals (vdW) force bonded molecular interfaces, which have little compositional/structural mismatch. Combining molecular dynamics simulations at temperatures up to the melting point and a time-domain analysis of interfacial heat currents in individual molecules, we unveil that (1) the heat flux across a molecular vdW interface has an alternating coupling-decoupling pattern at the single molecule level due to thermal fluctuations, and (2) the net heat transfer during the thermally coupled periods is orders of magnitude larger than that during the thermally decoupled periods. Furthermore, we demonstrate that a vdW interface's TIC correlates almost linearly with two atomistic-scale coupling-decoupling parameters: αc and σc, that represent the duration and strength of the thermal coupling, respectively. This thermal coupling-decoupling mechanism proposes a new paradigm that explicitly accounts for the effects of structural flexibility and dynamical fluctuation on heat transfer across molecular interfaces, and will likely provide new insights into modeling thermal transport properties in a wide-ranging soft materials, including self-assembled monolayers, paraffins or lipids.

42 ENGINEERING↗

Electrical Submersible Pumps: A System Modeling Approach for Power Quality Analysis With Variable Frequency Drives

This article proposes a high-level engineering guide to develop an integrated system model for power quality analysis in electrical submersible pump (ESP) applications with variable frequency drives (VFDs). Such analyses are troublesome in the industry due to the complexity of these systems. In this article, simple steps to perform system integration analyses of such arrangements, including torsional analysis, are developed. A simplified VFD-ESP model suitable for coupled electrical and mechanical analysis in steady state is proposed in all their configurations of practical interest. Such a model can be easily implemented in common simulation software, significantly reducing engineering efforts for implementation and analysis. The focus of the proposed model is the prediction of ESP failures that might result from a poor power quality caused by VFDs. Analytical expressions of different types of harmonics in these systems, as well as their accurate locations in the frequency domain, including their interharmonics and common-mode harmonics, are derived for this purpose. Further, the effectiveness of the proposed model is verified through offline and real-time hybrid simulation results. Finally, a comparison between simulation results obtained using the proposed model and measurements collected on a down-scale laboratory prototype is carried out to demonstrate the accuracy of the suggested modeling approach.

42 ENGINEERING↗

Hybrid quantum-classical approach for coupled-cluster Green's function theory

The three key elements of a quantum simulation are state preparation, time evolution, and measurement. While the complexity scaling of time evolution and measurements are well known, many state preparation methods are strongly system-dependent and require prior knowledge of the system's eigenvalue spectrum. Here, we report on a quantum-classical implementation of the coupled-cluster Green's function (CCGF) method, which replaces explicit ground state preparation with the task of applying unitary operators to a simple product state. While our approach is broadly applicable to many models, we demonstrate it here for the Anderson impurity model (AIM). The method requires a number of T gates that grows as O ( N 5 ) per time step to calculate the impurity Green's function in the time domain, where N is the total number of energy levels in the AIM. Since the number of T gates is analogous to the computational time complexity of a classical simulation, we achieve an order of magnitude improvement over a classical CCGF calculation of the same order, which requires O ( N 6 ) computational resources per time step.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Model-Free Dynamic Voltage Control of Distributed Energy Resource (DER)-Based Microgrids

In this paper, we present a new control technique for sustaining dynamic voltage stability by effective reactive power control and coordination of distributed energy resources (DERs) in microgrids. The proposed control technique is based on model-free control (MFC), which has shown successful operation and improved performance in different domains and applications. This paper presents its first use in the voltage stability of a microgrid setting employing multiple synchronous generator (SG)-based and power electronic (PE)-based DERs. MFC is a computationally efficient, data-driven control technique that does not require modelling of the different components and disturbances in the power system. It is utilized as an online controller to achieve the dynamic voltage stability of a microgrid system under different disturbances and fault conditions. A 21-bus microgrid system fed by multiple DERs is considered as a case study and the overall dynamic voltage stability is investigated using time-domain dynamic simulations. Numerical results show that the proposed MFC provides improvements on the dynamic load bus voltage profiles and requires less computational time as compared to the traditional enhanced microgrid voltage stabilizer (EMGVS) scheme. Due to its simplicity and low computational requirement, MFC can be easily implemented in resource-constrained computing devices such as smart inverters.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Advances in Machine and Deep Learning for Modeling and Real-time Detection of Multi-Messenger Sources

We live in momentous times. The science community is empowered with an arsenal of cosmic messengers to study the universe in unprecedented detail. Gravitational waves, electromagnetic waves, neutrinos, and cosmic rays cover a wide range of wavelengths and timescales. Combining and processing these datasets that vary in volume, speed, and dimensionality requires new modes of instrument coordination, funding, and international collaboration with a specialized human and technological infrastructure. In tandem with the advent of large-scale scientific facilities, the last decade has experienced an unprecedented transformation in computing and signal-processing algorithms. The combination of graphics processing units, deep learning, and the availability of open source, high-quality datasets has powered the rise of artificial intelligence. This digital revolution now powers a multibillion dollar industry, with far-reaching implications in technology and society. In this chapter, we describe pioneering efforts to adapt artificial intelligence algorithms to address computational grand challenges in multi-messenger astrophysics. We review the rapid evolution of these disruptive algorithms, from the first class of algorithms introduced in early 2017 to the sophisticated algorithms that now incorporate domain expertise in their architectural design and optimization schemes. We discuss the importance of scientific visualization and extreme-scale computing in reducing time-to-insight and obtaining new knowledge from the interplay between models and data.

Artificial Intelligence↗

Improved Modeling of Waste Isolation Pilot Plant Disposal Room Porosity

Disposal rooms at the Waste Isolation Pilot Plant (WIPP) contain waste and gas, and their porosity evolves over time. This report presents several improvements to the disposal room porosity model and presents new porosity predictions for use in future WIPP Performance Assessment activities. The improvements pertain to three sub-models: the geomechanical model, the waste compaction model, and the gas pressurization model. The impacts of each major improvement were quantified and the new porosity predictions were shown to be both mesh and domain size converged. Also, sensitivity studies on the disposal room horizon, clay seam friction coefficients, and homogenized waste representation were performed to support assumptions in the disposal room porosity model. To compare the legacy and new porosity predictions, the simulation results were plotted as a response surface, where gas pressure and time are inputs and porosity is the output. The new porosity response surface is insensitive to pressures beneath lithostatic pressure and highly sensitive to pressures above lithostatic pressure. The legacy porosity response surface, on the other hand, has moderate porosity gradients over all pressures. The new porosity response surface has a stronger scientific foundation than the legacy surface and may now be used for Compliance Decision Analyses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Active Learning for Anomaly Detection in Environmental data

Due to the growing amount of data from in-situ sensors in environmental monitoring, it becomes necessary to automatically detect anomalous data points. Nowadays, this is mainly performed using supervised machine learning models, which need a fully labelled data set for their training process. However, the process of labelling data is typically cumbersome and, as a result, a hindrance to the adoption of machine learning methods for automated anomaly detection. In this work, we propose to address this challenge by means of active learning. This method consists of querying the domain expert for the labels of only a selected subset of the full data set. We show that this reduces the time and costs associated to labelling while delivering the same or similar anomaly detection performances. Finally, we also show that machine learning models providing a nonlinear classification boundary are to be recommended for anomaly detection in complex environmental data sets.

54 ENVIRONMENTAL SCIENCES↗

Improved Weld Residual Stress Modeling System in BlackBear

This report presents enhancements to the MOOSE-based BlackBear application aimed at improving its capability to simulate welding and other thermo-mechanical manufacturing processes. Two primary avenues of improvement are pursued. First, to enhance user accessibility, we introduce a centralized default block restriction mechanism that ensures coverage checks are performed within user-specified default blocks. This default setting is applied consistently to all block-describable objects, such as variables, kernels, and more. In addition, we develop a modular action for moving heat source simulations, which integrates path file parsing, subdomain modification, and heat source kernel enforcement into a single, streamlined configuration. Second, to improve solver robustness, we implement an alternative method for assigning initial conditions to the updated active domain during the simulation, thereby enhancing convergence behavior. To validate the framework, we design and conduct several benchmark simulations, including heat conduction with progressive material addition, linear elasticity with time-dependent material deposition, and viscoplasticity model with isotropic hardening under similar conditions. Finally, we demonstrate the effectiveness of the proposed framework through large-scale thermo-mechanical welding simulations in both two and three dimensions.

42 ENGINEERING↗

NRC Multiphysics Analysis Capability Deployment FY2020: Part 3

This report details progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Sup- port for Advanced Non-Light Water Reactors.” The tasks completed for this report are: Task2c: Explicit modeling of pebble transient temperature response. In this simulation, the 400 MWth Pebble-Bed Modular Reactor (PBMR) design, PBMR- 400, experiences a 20-second power ramp from 100% to 150% power. This is followed by a similar reduction in the power back to 100%. Several multiscale pebble coupling approaches are tested with one pebble per mesh element in the active core region. The results show good conservation behavior and the stability of the coupling.; Extended scope part 1: An assessment of the computational efficiency of the Discontinuous Finite Element Method (DFEM) heat transfer solver shows good scalability. The DFEM solver is a factor of 4 more expensive in solution time than the Finite Element Method (FEM) solver for heat transfer problems due to the increased number of degrees of freedom. Nonetheless, the DFEM approach provides the user with the flexibility to model gap heat transfer problems.; Extended scope part 2: The GapHeatTransferInterfaceMaterial was improved to give the user increased flexibility with the modeling of heat transfer through gaps with the DFEM solver. A number of gap parameters can now be coupled both through functions and variables.; Extended scope part 3: Demonstration of how the gap width between hexagonal fuel cells can be calculated during a heat-up transient and used in the GapHeatTransferInterface model. A full-domain DFEM model with gap expansion is coupled to a SubApp that models the thermal expansion of the base plate. The results show the expected physical behavior, although have not been fully bench-marked at this point in time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

2021 Summer SPE Project

SPE is a project in order to develop new, more physics-based, seismic models of explosions (see Nelson et al. 2013). One key component of this effort is numerical modeling enabled by modern state-of-art code/software. Accurate modeling of the shape and amplitude of seismic waves from their generation to their propagation to remote monitoring seismic stations is important to our ability to determine the origin and strength of the source from remote recording. In this project, the modeling is performed by coupling two codes HOSS and SPECFEM3D. HOSS models the dynamic nonlinear processes happening near the explosion. SPECFEM3D computes the propagation of seismic waves as they travel through 3D complex Earth models. SPECFEM3D was modified in order to be driven by a set of time-series calculated by HOSS in lieu of a seismic source. The goal of this summer project is the investigation of several questions pertinent to the establishment of a full end-to-end modeling capability from the high-rate strain regime area to remote distances where seismic station record seismic waves generated by explosions. The investigated questions are: (1) How to perform proper filtering? Direct-solution modeling only sustains a limited range of frequency depending on the grid size. As we go from one modeling domain to the other via coupling, the mesh size is getting coarser to allow modeling at large scale but also to account for the fact that high-frequency waves do not physically travel to large distances. So filtering of the time-series generated by the near-field hydrodynamic modeling is a current practice often employed but its effects on the modeled waveforms has to be investigated. (2) Quantitative assessment of the efficiency of attenuation to remove high-frequency content of the wavefield. This assessment will allow to create meshes with a grid size appropriate to the actual physics of wave propagation for a given explosion. (3) Checking that the coupling process and the two codes respect the cylindrical symmetry that is expected in the case of a pure explosion in a half-space. (4) The effect of the state-of-stress in the near-source area on the modeled seismic waveforms. These questions will be investigated through the modeling of SPE-4P, the fourth explosion of this series because it was designed to have little interaction with the geologic surrounding and the free surface of the Earth so that it is the most explosion-like experiment, with the most symmetries to be verified.

58 GEOSCIENCES↗

Measurements of nonequilibrium interatomic forces using time-domain x-ray scattering

Here, we demonstrate an experimental approach to determining the excited-state interatomic forces using femtosecond x-ray pulses from an x-ray free-electron laser. We determine experimentally the excited-state interatomic forces that connect photoexcited carriers to the nonequilibrium lattice dynamics in the prototypical Peierls-distorted material, bismuth. The forces are obtained by a constrained least-squares fit of a pairwise interatomic force model to the excited-state phonon dispersion relation as measured by the time- and momentum-resolved x-ray diffuse scattering. We find that photoexcited carriers weaken predominantly the nearest-neighbor forces, which drives the measured softening of the transverse acoustic modes throughout the Brillouin zone as well as the zone-center A 1g optical mode. This demonstrates a bond-selective approach to measuring electron-phonon coupling relevant to a broad range of photoinduced phase transitions and transient light-driven states in quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Octofitter: Fast, Flexible, and Accurate Orbit Modeling to Detect Exoplanets

As next-generation imaging instruments and interferometers search for planets closer to their stars, they must contend with increasing orbital motion and longer integration times. These compounding effects make it difficult to detect faint planets but also present an opportunity. Increased orbital motion makes it possible to move the search for planets into the orbital domain, where direct images can be freely combined with the radial velocity and proper motion anomaly, even without a confirmed detection in any single epoch. In this paper, we present a fast and differentiable multimethod orbit-modeling and planet detection code called Octofitter. This code is designed to be highly modular and allows users to easily adjust priors, change parameterizations, and specify arbitrary function relations between the parameters of one or more planets. Octofitter further supplies tools for examining model outputs including prior and posterior predictive checks and simulation-based calibration. We demonstrate the capabilities of Octofitter on real and simulated data from different instruments and methods, including HD 91312, simulated JWST/NIRISS aperture masking interferometry observations, radial velocity curves, and grids of images from the Gemini Planet Imager. We show that Octofitter can reliably recover faint planets in long sequences of images with arbitrary orbital motion. This publicly available tool will enable the broad application of multiepoch and multimethod exoplanet detection, which could improve how future targeted ground- and space-based surveys are performed. Finally, its rapid convergence makes it a useful addition to the existing ecosystem of tools for modeling the orbits of directly imaged planets.

79 ASTRONOMY AND ASTROPHYSICS↗

Emergence of robust anthropogenic increase of heat stress-related variables projected from CORDEX-CORE climate simulations

The information of when and where region-specific patterns in both mean and extreme temperatures leading to heat stress will emerge from the present-day climate variability is important to plan adaptation options, but to date studies on this issue still remain limited and fragmented. Here, we estimate the time of emergence (TOE) of temperature and wet-bulb temperature (Tw), a better indication of heat stress, using fine-scale, long-term regional climate model projections under the RCP2.6 and RCP8.5 scenarios across six different domains. Differently from previous studies, the TOE is determined using three methods applied on impact-relevant variables: two different signal-to-noise frameworks based on summer mean temperature and Tw and a statistical test to identify significant differences in daily extreme distributions. The TOE response to RCP2.6 and RCP8.5 with respect to the end of 20th century variability differs significantly regardless of which TOE metric is applied. For summer mean temperature, the land fraction reaching TOE is expected to exceed 90% by the 2050s under the RCP8.5, whereas the increase rate of land exposure to TOE tends to stagnate over time under the RCP2.6 so that more than 40% of land will not experience TOE by the end of the 21st century. Compared to temperature, the TOE of Tw is reached earlier in most of the wet tropics but is delayed in hot and dry regions because of the nonlinear response of Tw to humidity. For both temperature and Tw, the TOE appears earlier in regions with low baseline variability, such as in the tropics. Despite the uncertainties arising from the choice of TOE metrics, the vast majority of regions in Africa and southeast Asia experience TOE in the early 21st century under both the RCP2.6 and RCP8.5 scenarios, which stresses the urgent need for developing adequate adaptation strategies in these regions.

54 ENVIRONMENTAL SCIENCES↗

Picosecond volume expansion drives a later-time insulator–metal transition in a nano-textured Mott insulator

There is significant technological interest in developing ever faster switching between different electronic and magnetic states of matter. Manipulating properties at terahertz rates requires accessing the intrinsic timescales of both electrons and associated phonons, which is possible with short-pulse photoexcitation. However, in many Mott insulators, the electronic transition is accompanied by the nucleation and growth of percolating domains of the changed lattice structure, leading to empirical timescales dominated by slowly coarsening dynamics. Here, in this study, we use time-resolved X-ray diffraction and reflectivity measurements to show that the photoinduced insulator-to-metal transition in an epitaxially strained Mott insulating thin film occurs without observable domain formation and coarsening effects, allowing the study of the intrinsic electronic and lattice dynamics. Above a fluence threshold, the initial electronic excitation drives a fast lattice rearrangement, which is followed by a slower electronic evolution into a metastable nonequilibrium state. Microscopic model calculations based on time-dependent dynamical mean-field theory and semiclassical lattice dynamics explain the threshold behaviour and elucidate the delayed onset of the electronic phase transition. This work highlights the importance of combined electronic and structural studies in unravelling the physics of dynamic transitions and the timescales of photoinduced processes. During a photoinduced phase transition, electronic rearrangements are usually faster than lattice ones. Time-resolved measurements now show that the insulator-to-metal transition in a thin-film Mott insulator is preceded by lattice reconfiguration.

36 MATERIALS SCIENCE↗

Assessing free tropospheric quasi-equilibrium for different GCM resolutions using a cloud-resolving model simulation of tropical convection

Abstract This study examines the free-tropospheric quasi-equilibrium at different global climate model (GCM) resolutions using the simulation of tropical convection by a cloud-resolving model during the Tropical Western Pacific International Cloud Experiment. The simulated dynamic and thermodynamic fields within the model domain are averaged over subdomains of different sizes equivalent to different GCM resolutions. These coarse-grained fields are then used to compute CAPE and its change with time, and their relationships with simulated convection. Results show that CAPE change with time is controlled predominantly by variations of thermodynamic properties in the planetary boundary layer for all subdomain sizes ranging from 64 to 4 km. Lag correlation analysis shows that CAPE generation by the free-tropospheric dynamical advection (dCAPE ls ) leads convective precipitation but is in phase with convective mass flux at 600 mb and 500 mb vertical velocity for all subdomain sizes. However, the correlation coefficients and regression slopes decrease as the subdomain size decreases for subdomain sizes smaller than 16 km. This is probably due to increased randomness of convection and more scale-dependence of the relationships when the subdomain size reaches the grey zone. By examining the sensitivity of the relationships of convection with dCAPE ls to temporal scales in different subdomain size, it shows that the quasi-equilibrium between dCAPE ls and convection holds well for timescales of 30 min or longer at all subdomain sizes. These results suggest that the free tropospheric quasi-equilibrium assumption may still be useable even for GCM resolutions in the grey zone.

54 ENVIRONMENTAL SCIENCES↗