Time domain modeling of constant power loads for electromagnetic transient simulations
Not Available
Alboaouh, Kamel↗
SEARCH · Engineering Papers
Results for “time dependent behavior”
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.
Spatially-resolved lithiation dynamics from operando X-ray diffraction and electrochemical modeling of lithium-ion cells
In this study, energy dispersive X-ray diffraction is used to profile the time evolution of ordered Li x C 6 phases in solid electrodes of lithium-ion cells charged at rates between 0.2 and 4.7C (where 1C corresponds to full discharge in 1 h). The methods for quantifying lithium concentration in these phases from the acquired diffraction patterns are described. Compact expressions for time-dependent concentration gradients in the solid electrodes using orthogonal polynomial expansions are presented. Experimentally, these gradients persisted in lithiated graphite electrodes even after the cells rested at open-circuit for over 9 h. A multiphase electrochemical model of graphite intercalation captured many of the observed behaviors, including the progression of phase transitions and the persistent gradients at zero current. However, the magnitude of concentration gradients in both the oxide cathode and graphite anode is underestimated by the model, even at moderate currents.
25 ENERGY STORAGE↗
Effects of anisotropic shale creep on the stress and permeability evolution of a geological nuclear waste repository
To ensure long-term safety and performance, geological nuclear waste repositories require low-permeability barriers such as bentonite buffers and/or shale host rock. Shale is not only known for its low permeability but also for its trend to undergo time-dependent deformation (i.e., creep), which could heal damage, but the effects of shale creep on the long-term performance of nuclear waste repositories have not been clearly understood. In particular, the anisotropic nature of shale (i.e., bedding) could have a significant effect on its creep behavior, and consequently, on the long-term performance of nuclear waste repositories. In this research, numerical simulations were carried out with the objective of showing the effects of anisotropic shale creep on the stress and permeability evolution of a generic geological nuclear waste repository in shale. The TOUGH-FLAC simulator, a thermo-hydromechanically (THM) coupled numerical code, was used for the simulations. To achieve the objective, comparisons were performed between the results of anisotropic shale creep simulations and those of different simulation cases, namely, no creep (i.e., elastic), isotropic creep, and long-term creep shale cases. Results of the comparisons show that the elastic and isotropic creep shale cases respectively led to the overestimation and underestimation of stress and permeability in the repository, whereas the long-term creep shale case, which accumulated greater creep in later periods than in earlier periods, helped to keep large shear and tensile stresses from developing while maintaining compressive spherical stress, resulting in consistently low permeability levels. These results indicate that performance assessments with elastic and isotropic creep formation models will provide the upper and lower bound estimates of stress and permeability, while more reasonable estimates will be given by an anisotropic creep formation model, and that shale with long-term creep characteristics will be beneficial in many aspects of the safety and performance of nuclear waste repositories.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Predicting Quantum Criticality in Single-Crystalline Ce 2 Ru 3 Ge 5
Strongly correlated f-electron systems are known to host exotic quantum states, such as quantum criticality, complex order parameters, and unconventional superconductivity. However, the appearance of these exotic states is difficult to predict, making the study of quantum critical behavior challenging, especially in ferromagnetic materials. Herein, we report a structure–property map for Ce 2 M 3 X 5 (M = transition metal; X = main group element) that aids in the targeted design of materials likely to exhibit quantum criticality. Guided by this map, we report on the synthesis of single-crystalline Ce 2 Ru 3 Ge 5 and provide, for the first time, magnetic susceptibility, heat capacity, resistivity, and magnetoresistance measurements on single crystals. Here, we observe a weak ferromagnetic-like response at 7.5 K, which is contrasted with the bulk ferromagnetic ordering that appears in polycrystalline samples. Non-Fermi liquid behavior is seen in the temperature dependent electrical resistivity and heat capacity of the single crystals, suggesting proximity to a ferromagnetic quantum critical point without chemical or physical pressure. Given the contrast with previous reports of polycrystals, these results lead us to propose that single crystalline Ce 2 Ru 3 Ge 5 is intrinsically tuned into the vicinity of a ferromagnetic quantum critical point.
36 MATERIALS SCIENCE↗
Viscoelastic Material Characterization and Modeling of Photovoltaic Module Packaging Materials for Direct Finite-Element Method Input
Numerical tools, such as the finite-element method, are increasingly used to design and evaluate the photovoltaic (PV) modules, providing for the reduction of development time and improved performance and reliability. However, high-fidelity material models are necessary to accurately model the complex structural behavior of the involved packaging materials. A common simplification used in recent years is to model the polymer materials (i.e., encapsulant and backsheet) as linear elastic, which will lead to inaccurate results. Therefore, in this work, we present a thorough characterization of the time- and temperature-dependent mechanical response of predominant PV module encapsulant and backsheet materials. Based on this material characterization, we developed and experimentally validated generalized Maxwell models to describe each material's viscoelastic response. In addition, we included measurements of the coefficient of thermal expansion and presented all material models in such a fashion for direct input into commercial finite-element method modeling software.
FEM↗
Transient Analysis of a Micro-reactor using the DireWolf Code Suite
Transient analyses of heat pipe micro-reactors are necessary to ensure that hypothetical accident scenarios do not comprise reactor safety. Due to its small size and reliance on heat-pipes for cooling, the micro-reactor design introduced in this paper is a tightly coupled system which requires multi-physics tools to accurately model transient events. Idaho National Laboratory’s DireWolf code suite based on the MOOSE framework is tailor-built to model heat-pipe reactors. This paper demonstrates DireWolf’s ability to simulate the coupled thermal-neutronics transient behavior of a heat-pipe micro-reactor. The transient events presented here include an inadvertent rotation of all control drums simultaneously and a sudden complete ejection of a single control drum. A detailed description of each event is provided along with simulation results, including time dependent power and temperature distributions, and discussion. This is a Westinghouse Electric Company (WEC) led publication.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Transient analysis of a micro-reactor using the DireWolf code suite
Transient analyses of heat pipe micro-reactors are necessary to ensure that hypothetical accident scenarios do not comprise reactor safety. Due to its small size and reliance on heat-pipes for cooling, the micro-reactor design introduced in this paper is a tightly coupled system which requires multi-physics tools to accurately model transient events. Idaho National Laboratory's DireWolf code suite based on the MOOSE framework is tailor-built to model heat-pipe reactors. This paper demonstrates DireWolf's ability to simulate the coupled thermal-neutronics transient behavior of a heat-pipe micro-reactor. The transient events presented here include an inadvertent rotation of all control drums simultaneously and a sudden complete rotation of a single control drum. A detailed description of each event is provided along with simulation results, including time dependent power and temperature distributions, and discussion. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Data-driven Mori–Zwanzig modeling of Lagrangian particle dynamics in turbulent flows
The dynamics of Lagrangian particles in turbulence play a crucial role in mixing, transport, and dispersion in complex flows. Their trajectories exhibit highly nontrivial statistical behavior, motivating the development of surrogate models that can reproduce these trajectories without incurring the high computational cost of direct numerical simulations of the full Eulerian field. This task is particularly challenging because reduced-order models typically lack access to the full set of interactions with the underlying turbulent field. Novel data-driven machine learning techniques can be powerful in capturing and reproducing complex statistics of the reduced-order/surrogate dynamics. In this work, we show how one can learn a surrogate dynamical system that is able to evolve a turbulent Lagrangian trajectory in a way that is point-wise accurate for short-time predictions (with respect to Kolmogorov time) and stable and statistically accurate at long times. This approach is based on the Mori–Zwanzig formalism, which prescribes a mathematical decomposition of the full dynamical system into resolved dynamics that depend on the current state and the past history of a reduced set of observables, and the unresolved orthogonal dynamics due to unresolved degrees of freedom of the initial state. We show how by training this reduced order model on a point-wise error metric on short time-prediction, we are able to correctly learn the dynamics of Lagrangian turbulence, such that also the long-time statistical behavior is stably recovered at test time. This opens up a range of applications, for example, for the control of active Lagrangian agents in turbulence.
97 MATHEMATICS AND COMPUTING↗
Electronic mobility, doping, and defects in epitaxial BaZrS 3 chalcogenide perovskite thin films
We present the electronic transport properties of BaZrS 3 thin films grown epitaxially by gas-source molecular beam epitaxy. We observe n-type behavior in all samples, with carrier concentration ranging from 4 × 10 18 to 4 × 10 20 cm −3 at room temperature (RT). We observe a champion RT Hall mobility of 11.1 cm 2 V −1 s −1 , which is competitive with established thin-film photovoltaic absorbers. Temperature-dependent Hall mobility data show that phonon scattering dominates at room temperature, in agreement with computational predictions. X-ray diffraction data illustrate a correlation between mobility and antiphase boundary concentration, illustrating how microstructure can affect transport. Despite the well-established environmental stability of chalcogenide perovskites, we observe significant changes to electronic properties as a function of storage time in ambient conditions. With the help of secondary ion mass spectrometry measurements, we propose and support a defect mechanism that explains this behavior: as-grown films have a high concentration of sulfur vacancies that are shallow donors (V ⋅ S or V ⋅⋅ S ), which are converted into neutral oxygen defects (O$^{×}_{S}$) upon air exposure. We discuss the relevance of this defect mechanism within the larger context of chalcogenide perovskite research, and we identify means to stabilize the electronic properties.
Van Sambeek, Jack [Massachusetts Institute of Tech↗
Insights into mixing of non-isothermal multi-polymer melts for complex plastics recycling
Catalytic recycling or upcycling of plastics is often limited not by catalyst performance, but by transport, arising from highly viscous, non-Newtonian polymer melts. In this work, we develop a reactor-scale framework that integrates rheological measurements, constitutive modeling, computational fluid dynamics (CFD), and experiments to quantify mixing, heat transfer, and dispersion in surrogate hydrocarbon melts representing mixed plastics systems. Temperature- and shear rate-dependent viscosity of low-density polyethylene (LDPE) and high-density polyethylene (HDPE) is measured to create two surrogate polymers (PLD and PHD) that capture the dominant shear-thinning flow behavior while neglecting strong elastic effects, enabling tractable simulation of non-isothermal, polymer-melt mixing using a Carreau-Arrhenius generalized Newtonian framework. Three-dimensional CFD simulations are employed to evaluate impeller performance in PLD using mixing time, cavern volume, thermal uniformity, and interfacial area for regimes in which viscoelastic effects are not dominant. We show that magnetic stir bars commonly used in lab-scale studies produce large thermal gradients (~60 °C) and poor mixing, even under idealized power delivery and polymer flow conditions. In contrast, close-clearance anchor impellers achieve near-isothermal operation, reduce mixing times by up to 5×, and provide >90% active circulation volume. We further demonstrate that, at low pseudo-Deborah number (De*), motor power requirements can be predicted directly from shear rate-dependent rheology using the Carreau-Arrhenius framework, enabling rational selection of operating conditions. Extension to surrogate immiscible multi-polymer systems based on PLD and PHD shows that interfacial area is highly sensitive to operating conditions and impeller design, with coaxial anchor-turbine configurations enhancing dispersion by up to 4 × .
Close-clearance impellers↗
Dynamics of argon metastables in Ar–CH 4 radio frequency capacitively-coupled plasma: real-time monitoring with neural network-augmented broadband optical emission spectroscopy
In moderate-pressure radio frequency (RF) capacitively coupled plasmas generated in argon–methane mixtures, the density of argon metastable atoms (Ar 1 s 5 ) exhibits a non-monotonic dependence on methane (CH 4 ) concentration. Laser-induced fluorescence (LIF) was used to measure and compare local Ar 1 s 5 densities in Ar and Ar–CH 4 plasmas at 2.6 Torr and RF powers of 17–117 W. The addition of 1% CH 4 increases the metastable density, and 2% CH 4 triggers a strong depletion by an order of magnitude, compared to 1% CH 4 case. This non-monotonic behavior demonstrates the sensitivity of metastable populations to small gas admixtures, which is critical for processes where metastables drive precursor dissociation. For real-time monitoring of metastable population, broadband optical emission spectroscopy (OES) is augmented with a feedforward neural network (NN) to predict Ar s1 5 densities from spectral features. When trained on LIF data, the NN replicates the absolute densities and the dynamic trends of Ar 1 s 5 density variation. The NN-augmented broadband OES approach can be used as a simple and cost-effective tool for tracking Ar metastables in Ar-rich plasmas, facilitating industrial-scale optimization.
Yatom, Shurik [Princeton Plasma Physics Laboratory↗
Influence of Kosmotrope and Chaotrope Salts on Water Structural Relaxation
The structural relaxation in water solutions of kosmotrope (structure maker) and chaotrope (structure breaker) salts, namely sodium chloride, potassium chloride, and cesium chloride, were studied through quasielastic neutron scattering measurements. We found that the collective dynamics relaxation time at the structure factor peak obtained using heavy water solutions shows a distinctively different behavior in the kosmotrope as opposed to the chaotrope solutions, increasing with the salt concentration in the former and decreasing in the latter. In both cases the trends are proportional to the concentration dependence of the relative viscosity of the solutions. Furthermore, these results indicate that kosmotropes and chaotropes influence the solution’s viscosity by impacting in opposite ways the hydrogen bond network of water, strengthening it in one case and softening it in the other.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Quantitative kinetic rules for plastic strain-induced α - ω phase transformation in Zr under high pressure
Plastic strain-induced phase transformations (PTs) and chemical reactions under high pressure are broadly spread in modern technologies, friction and wear, geophysics, and astrogeology. However, because of very heterogeneous fields of plastic strain $E$ p and stress σ tensors and volume fraction c of phases in a sample compressed in a diamond anvil cell (DAC) and impossibility of measurements of σ and $E$ p , there are no strict kinetic equations for them. Here, we develop a kinetic model, finite element method (FEM) approach, and combined FEM-experimental approaches to determine all fields in strongly plastically predeformed Zr compressed in DAC, and specific kinetic equation for α-ω PT consistent with experimental data for the entire sample. Since all fields in the sample are very heterogeneous, data are obtained for numerous complex 7D paths in the space of 3 components of the plastic strain tensor and 4 components of the stress tensor. Kinetic equation depends on accumulated plastic strain (instead of time) and pressure and is independent of plastic strain and deviatoric stress tensors, i.e., it can be applied for various above processes. Our results initiate kinetic studies of strain-induced PTs and provide efforts toward more comprehensive understanding of material behavior in extreme conditions.
36 MATERIALS SCIENCE↗
Breakdown front dynamics of streamer-like discharge in electron-charged polymethyl methacrylate
We present the results of an experimental study of fast breakdown in electron-charged polymethyl methacrylate. We irradiate bulk polymethyl methacrylate disks with diameters up to one meter at different implanted charge densities and measure the discharge current during the forced electrical breakdown of the material. We infer the breakdown dynamics from these current waveforms, including the velocity and time dependence of the electric field driving the breakdown, and compare these results with the physical electrical tree patterns left behind in the material. We find that the dynamics and physical characteristics of the breakdown channels in electron-irradiated solids depart from typical expectations of electrical treeing behavior in solid materials. Thus, we interpret these results as an expression of streamer discharges in dense gases motivated by the existence of trapped gases in the solid due to radiation damage. We show that the dynamics of the breakdown channels in the electron-charged solid dielectric material is well described by applying standard streamer physics to this physical system. The results show that the dynamics of breakdown channels in solid dielectric material is a promising avenue to further understand streamer discharges in different media under extreme conditions.
36 MATERIALS SCIENCE↗
Plasmonic Hot-Carrier Generation and Catalysis in Ti 3 C 2 O 2 from Real-Time TDDFT Simulations
Photoinduced hot electrons are central to plasmon-driven catalysis. Atomically thin Ti 3 C 2 O 2 , with high carrier density and broad optical absorption, offers a promising platform for plasmon-driven reactions. However, comprehensive investigations of its plasmon resonance, hot-carrier generation, and plasmonic catalytic performance remain limited. In this work, real-time time-dependent density functional theory (rt-TDDFT) was employed to study Ti 3 C 2 O 2 ’s plasmon excitation and hot-carrier generation from nonradiative plasmon damping. The temporal evolution of the dipole moment reveals plasmon resonance in Ti 3 C 2 O 2 , followed by strong plasmon damping that redistributes the stored energy to generate hot carriers. Ti 3 C 2 O 2 with low oxygen vacancy concentration (O v -Ti 3 C 2 O 2 ) exhibits plasmonic behavior resembling the pristine surface, and the plasmon-generated hot electrons can markedly reduce the dissociation barrier of CO 2 at the oxygen vacancy. These findings provide fundamental insights into the plasmonic properties of Ti 3 C 2 O 2 and how they drive its catalytic performance in surface reactions, which is valuable for advancing plasmon-driven catalysis.
CO2 reduction↗
Alluvial fan morphology: A self-similar free boundary problem description
In this work, we examine approximate geometrically self-similar solutions to a parabolic free boundary value problem applied to alluvial fan surface morphology and growth. Alluvial fans are fan- or cone-shaped sedimentary deposits caused by the rapid deposition of sediment from a canyon discharging onto a flatter plain. Longitudinal, topographic profiles of fans can be readily described by a seemingly time independent dimensionless profile (DeChant et al., 1999). However, because an alluvial fan can be expected to grow over time, it is clear that this “steady” profile is certainly time dependent and can be described using a space-time self-similar solution. In an experimental and theory-based study, Guerit et al. (2014) developed a self-similar (or as they describe it a self-affine) linear solution based upon an approximate first order small parameter expansion solution for a 1-d homogeneous nonlinear diffusion equation. Direct substitution of this result into a linear diffusion equation suggests that this first order expression may not fully satisfy the associated governing equation. In contrast, we develop a more complete solution based upon a modeled approximation for the axi-symmetric formulation such that the associated temporal behavior is consistent with a 1/3 time power-law as described by Reitz and Jerolmack (2014). The resulting expression is an exact solution to a linear heat equation. Furthermore, we emphasize that a small parameter is not inherent to the resulting profile result and is not included in our model development. Though developed using rather different approaches, the formal solution developed here is in good agreement with the simple polynomial described by DeChant et al. (1999) suggesting that this self-similar solution is a suitable time dependent representation of alluvial fan longitudinal profile form and improves on earlier work.
58 GEOSCIENCES↗
Behavior of Tidally Triggered Earthquakes Depends on Fluid Conditions
Seismicity at The Geysers geothermal field and in Oklahoma is heavily influenced by industrial activities related to energy production, though the mechanism in which earthquakes are induced or triggered is different. At The Geysers, much of the seismicity is linked to thermoelastic stresses caused by injecting cold water into hot rocks, whereas in Oklahoma, the seismicity is linked to a reduction of confining stress on faults due to increasing pore pressure resulting from wastewater injections. Here, in this study, we show that these contrasting conditions are also evident in tidally triggered earthquakes. At The Geysers, earthquakes preferentially occur during maximum extensional stress, which does not occur at the same time as maximum shear stress on optimally oriented faults in the regional stress field. In Oklahoma, earthquakes preferentially occur during maximum shear stress on optimally oriented faults rather than maximum extensional stress. The magnitude of tidal extensional stress is naturally much greater than tidal shear stress. However, in a fluid-saturated environment, pore pressure responds to changes in volume, which can counteract or reduce the effect of the applied stress. The difference in behavior at these two sites is indicative of the level of coupling between applied stress and pore pressure, corresponding to unsaturated conditions at The Geysers and high pore pressure in Oklahoma.
58 GEOSCIENCES↗
Charge order driven multiferroic behavior in Sr 4 Fe 6 O 12 : An ab-initio study
Here, in this paper, we report the structural, electronic, and ferroelectric properties of the layered mixed-valent transition-metal compound Sr 4 Fe 6 O 12 (SFO). We demonstrate how SFO undergoes a phase transition from a high-temperature (𝑇) centrosymmetric tetragonal phase (𝑃4 2 /𝑚𝑛𝑚) to a low-𝑇 polar orthorhombic phase (𝑃𝑚𝑛2 1 ). The transition is primarily driven by charge ordering at tetrahedral Fe-layer creating Fe 3+ and Fe 2+ cations between two edge-sharing tetrahedra. This charge ordering induces electronic polarization, which is remarkably larger (3.5 times) in magnitude than ionic polarization and oppositely directed, giving a net polarization of 0.05 C/m 2 which is comparable to the state-of-the-art rare-earth nickelets and manganite perovskites. The direction of structural distortion, governed by the polar mode irrep $Γ^−_5$, depends sensitively on the type of magnetic ordering in the Fe-octahedral layer. Consequently, both the ionic and electronic polarization directions are influenced by magnetic ordering, suggesting the potential for multiferroic behavior with strong magnetoelectric coupling in this material.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗