Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transient sources”

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 361 records · Page 20

An Optically Discovered Outburst from XTE J1859+226

Abstract Using the Zwicky Transient Facility, in 2021 February we identified the first known outburst of the black hole X-ray transient XTE J1859+226 since its discovery in 1999. The outburst was visible at X-ray, UV, and optical wavelengths for less than 20 days, substantially shorter than its full outburst of 320 days in 1999, and the observed peak luminosity was 2 orders of magnitude lower. Its peak bolometric luminosity was only 2 × 10 35 erg s −1 , implying an Eddington fraction of about 3 × 10 −4 . The source remained in the hard spectral state throughout the outburst. From optical spectroscopy measurements we estimate an outer disk radius of 10 11 cm. The low observed X-ray luminosity is not sufficient to irradiate the entire disk, but we observe a surprising exponential decline in the X-ray light curve. These observations highlight the potential of optical and infrared synoptic surveys to discover low-luminosity activity from X-ray transients.

Astronomy & Astrophysics↗

Tree hydrodynamic modelling of the soil–plant–atmosphere continuum using FETCH3

Modelling the water transport along the soil–plant–atmosphere continuum is fundamental to estimating and predicting transpiration fluxes. A Finite-difference Ecosystem-scale Tree Crown Hydrodynamics model (FETCH3) for the water fluxes across the soil–plant–atmosphere continuum is presented here. The model combines the water transport pathways into one vertical dimension, and assumes that the water flow through the soil, roots, and above-ground xylem can be approximated as flow in porous media. This results in a system of three partial differential equations, resembling the Richardson–Richards equation, describing the transport of water through the plant system and with additional terms representing sinks and sources for the transfer of water from the soil to the roots and from the leaves to the atmosphere. The numerical scheme, developed in Python 3, was tested against exact analytical solutions for steady state and transient conditions using simplified but realistic model parameterizations. The model was also used to simulate a previously published case study, where observed transpiration rates were available, to evaluate model performance. With the same model setup as the published case study, FETCH3 results were in agreement with observations. Through a rigorous coupling of soil, root xylem, and stem xylem, FETCH3 can account for variable water capacitance, while conserving mass and the continuity of the water potential between these three layers. FETCH3 provides a ready-to-use open access numerical model for the simulation of water fluxes across the soil–plant–atmosphere continuum.

54 ENVIRONMENTAL SCIENCES↗

Enhanced Large-Signal Stability Method for Grid-Forming Inverters During Current Limiting: Preprint

Grid-forming (GFM) inverters are a promising technology for the widespread integration of renewable energy sources in future power systems. As a key element of GFM inverter control, the primary controller governs the internal reference voltage and angle. During contingencies in the grid---such as faults, voltage drops, or frequency and phase jumps---an inverter can be forced into a current-limiting mode of operation modulating inverter dynamics, and, as a result, it is prone to losing synchronism with the grid. In this paper, we propose a novel GFM primary control method with an additional synchronization term that naturally activates during contingencies to improve the dynamic response. The method allows the inverter to remain synchronized with the grid, which improves the inverter's dynamic behavior both during and after current-limiting grid conditions and enhances grid support, including voltage support using full current capacity. The method is demonstrated for voltage, frequency, and phase jumps both in a single-machine-to-infinite-bus and a network-wide electromagnetic transient simulation of the IEEE 14-bus system with 5 GFM inverters. The simulations provide insights into the proposed synchronization method and confirm the high potential of the method, which robustly secures synchronism under severe contingencies.

current limiting↗

Hardware Demonstration of Weak Grid Oscillations in Grid-Following Converters

The increasing penetration of inverter-based resources (IBRs) introduces some unexpected dynamic issues, including low-frequency oscillations. In this paper, a laboratory-scale grid-following voltage-source converter (VSC) system is implemented to demonstrate weak grid oscillations. Grid-following control is applied in VSC to provide active power, reactive power or ac voltage control. The test bed is also replicated in electromagnetic transient (EMT) simulation environment (MAT-LAB/SimPowerSystems) for benchmark purpose. Case studies are carried out to demonstrate low-frequency oscillations under real power/ac voltage or real power/reactive power control. The excellent agreement between the simulation and experimental results ensures that the two test beds are exchangeable and can be used for testing oscillation mitigation strategies in future research.

Bao, Li↗

Modeling a ring magnet in ALEGRA

We show here that Sandia's ALEGRA software can be used to model a permanent magnet in 2D and 3D, with accuracy matching that of the open-source commercial software FEMM. This is done by conducting simulations and experimental measurements for a commercial-grade N42 neodymium alloy ring magnet with a measured magnetic field strength of approximately 0.4 T in its immediate vicinity. Transient simulations using ALEGRA and static simulations using FEMM are conducted. Comparisons are made between simulations and measurements, and amongst the simulations, for sample locations in the steady-state magnetic field. The comparisons show that all models capture the data to within 7%. The FEMM and ALEGRA results agree to within approximately 2%. The most accurate solutions in ALEGRA are obtained using quadrilateral or hexahedral elements. In the case where iron shielding disks are included in the magnetized space, ALEGRA simulations are considerably more expensive because of the increased magnetic diffusion time, but FEMM and ALEGRA results are still in agreement. The magnetic field data are portable to other software interfaces using the Exodus file format.

36 MATERIALS SCIENCE↗

Recovered supernova Ia rate from simulated LSST images

Aims.TheVera C. RubinObservatory’s Legacy Survey of Space and Time (LSST) will revolutionize time-domain astronomy by detecting millions of different transients. In particular, it is expected to increase the number of known type Ia supernovae (SN Ia) by a factor of 100 compared to existing samples up to redshift ∼1.2. Such a high number of events will dramatically reduce statistical uncertainties in the analysis of the properties and rates of these objects. However, the impact of all other sources of uncertainty on the measurement of the SN Ia rate must still be evaluated. The comprehension and reduction of such uncertainties will be fundamental both for cosmology and stellar evolution studies, as measuring the SN Ia rate can put constraints on the evolutionary scenarios of different SN Ia progenitors. Methods.We used simulated data from the Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2) and LSST Data Preview 0 to measure the SN Ia rate on a 15 deg 2 region of the “wide-fast-deep” area. We selected a sample of SN candidates detected in difference images, associated them to the host galaxy with a specially developed algorithm, and retrieved their photometric redshifts. We then tested different light-curve classification methods, with and without redshift priors (albeit ignoring contamination from other transients, as DC2 contains only SN Ia). We discuss how the distribution in redshift measured for the SN candidates changes according to the selected host galaxy and redshift estimate. Results.We measured the SN Ia rate, analyzing the impact of uncertainties due to photometric redshift, host-galaxy association and classification on the distribution in redshift of the starting sample. We find that we are missing 17% of the SN Ia, on average, with respect to the simulated sample. As 10% of the mismatch is due to the uncertainty on the photometric redshift alone (which also affects classification when used as a prior), we conclude that this parameter is the major source of uncertainty. We discuss possible reduction of the errors in the measurement of the SN Ia rate, including synergies with other surveys, which may help us to use the rate to discriminate different progenitor models.

Astronomy & Astrophysics↗

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)↗

Improved PV System Control Strategies to Reduce Power Management Costs in Nanogrids

An improved PV system control method is proposed to reduce nanogrid operation costs in this paper. A model including it various components such as photovoltaic (PV) systems, energy storage systems (ESSs), gateways, and household loads is considered with the constraints of the ESS, PV irradiance from real data, and household load. We designed an optimal economic dispatch strategy with an improved PV system control method. Combining the proposed optimal economic dispatch and PV system control strategies, it can improve the control performance for both transient and steady-state responses thereby enabling the maximum power to be extracted from the PV. Consequently, the PV power is maximized, which allows the ESS to use less power and sell the surplus to external power sources, which means the proposed method decreases the nanogrid operation costs. Furthermore, this performance is verified via nanogrid simulations and PV experimental kit.

PV system control↗

Efficient SABRE-SHEATH Hyperpolarization of Potent Branched-Chain-Amino-Acid Metabolic Probe [1-13C]ketoisocaproate

Efficient 13C hyperpolarization of ketoisocaproate is demonstrated in natural isotopic abundance and [1-13C]enriched forms via SABRE-SHEATH (Signal Amplification By Reversible Exchange in SHield Enables Alignment Transfer to Heteronuclei). Parahydrogen, as the source of nuclear spin order, and ketoisocaproate undergo simultaneous chemical exchange with an Ir-IMes-based hexacoordinate complex in CD3OD. SABRE-SHEATH enables spontaneous polarization transfer from parahydrogen-derived hydrides to the 13C nucleus of transiently bound ketoisocaproate. 13C polarization values of up to 18% are achieved at the 1-13C site in 1 min in the liquid state at 30 mM substrate concentration. The efficient polarization build-up becomes possible due to favorable relaxation dynamics. Specifically, the exponential build-up time constant (14.3 ± 0.6 s) is substantially lower than the corresponding polarization decay time constant (22.8 ± 1.2 s) at the optimum polarization transfer field (0.4 microtesla) and temperature (10 °C). The experiments with natural abundance ketoisocaproate revealed polarization level on the 13C-2 site of less than 1%—i.e., one order of magnitude lower than that of the 1-13C site—which is only partially due to more-efficient relaxation dynamics in sub-microtesla fields. We rationalize the overall much lower 13C-2 polarization efficiency in part by less favorable catalyst-binding dynamics of the C-2 site. Pilot SABRE experiments at pH 4.0 (acidified sample) versus pH 6.1 (unaltered sodium [1-13C]ketoisocaproate) reveal substantial modulation of SABRE-SHEATH processes by pH, warranting future systematic pH titration studies of ketoisocaproate, as well as other structurally similar ketocarboxylate motifs including pyruvate and alpha-ketoglutarate, with the overarching goal of maximizing 13C polarization levels in these potent molecular probes. Finally, we also report on the pilot post-mortem use of HP [1-13C]ketoisocaproate in a euthanized mouse, demonstrating that SABRE-hyperpolarized 13C contrast agents hold promise for future metabolic studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Enabling Load Following Capability in the Transatomic Power MSR

This project is dedicated towards designing a fuel processing system that enables liquid-fueled molten salt reactors (MSR) to load follow by removing the dissolved xenon in the fuel salt. As one of the Gen-IV nuclear reactor concepts, the molten salt reactor receives increasing development interests in the recent years. One distinguishing feature of the liquid-fueled molten salt reactor is its improved ability to operate in a load-following mode by including the unique online fission product removal system. Load-following means that the reactor changes its power output based on the demand on the grid. Most of the current operating nuclear reactors have limited load-following ability and operate as the base load on the grid. Due to the rapid increase of solar energy, the requirement on load-following capacity is significantly increased because of the varying power output of the solar panels, yet the traditional load-following capacity is expected to decrease as the decarbonization of the grid continues. Therefore, the ability to perform load-following operation for the nuclear reactors will greatly enhance the resilience of the grid and make nuclear energy more economically competitive. This load-following feature is included in many commercial molten salt reactor designs, such as the designs by Transatomic Power, Terrestrial Energy, and Flibe Energy. Unfortunately, detailed analysis of the fuel processing system for commercial scale MSRs is still lacking, as well as how the fuel processing quantitively impacts the load-following operation. Moreover, experimental data for many of the underlying physics of fuel processing is limited. This project aims to pave the way for the fuel processing technology to advance to the commercial stage by performing combined experimental and simulation research. During the project period, four interconnected aspects of the development of the fuel processing system in liquid-fueled molten salt reactors are investigated. These aspects are the simulation and analysis of the fission product removal system, the fuel cycle simulation, the coupled reactor neutronics and thermal hydraulics transient simulation, and the gaseous fission product removal experiment. Multiphase CFD simulations are performed for components of the processing systems, and simplified air-water experiments are carried out to provide validation data. It is concluded that the CFD simulation can satisfactorily predict the system level performance of the components, and engineering models are constructed based on this success. Fuel cycle analysis is performed for two representative MSR design, the MSBR and the Transatomic Power MSR. Open-source code SaltProc is developed to incorporate the unique fuel processing system of the MSRs. It is concluded that the removal of xenon is essential for load-following operation in thermal spectrum MSR and Molten Salt Breeder Reactor. For the Transatomic Power MSR, the xenon poisoning effect is negligible due to its relatively fast neutron spectrum, though the overall fuel cycle economics still benefits from the removal of xenon. Coupled reactor neutronics and thermal hydraulics transient simulation is performed specifically for the Transatomic Power MSR. It is concluded that the reactor core design could perform power ramping fast enough to satisfy load-following operation. Combining the findings from each aspect, it is concluded that the load-following operation of a thermal neutron MSR is dependent upon the removal of xenon, which could be achieved for a commercial sized reactor using continuous inert gas sparging in a separate system with reasonable dimensions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Reflections on bubble walls

We discuss the dynamics of expanding bubble walls in the presence of massive dark photons whose mass changes as they cross the wall. For sufficiently thin walls, we show that there exists a transient kinematic regime characterized by a constant reflection probability of longitudinal — but not transverse — modes. This effect can have important implications for the dynamics of expanding vacuum bubbles in the early Universe. Most notably, it leads to a new source of pressure on the expanding interface, featuring a non-monotonic dependence on the γ-factor of the bubble walls and reaching a peak at intermediate γ-factors that we dub Maximum Dynamic Pressure. When this pressure is large enough to halt the acceleration of the bubble walls, the difference in vacuum energy densities goes into making a fraction of the dark photons relativistic, turning them into dark radiation. If the dark radiation remains relativistic until late times, an observable contribution to ΔNeff is possible for phase transitions with strength α ~ 10 -2 - 10 -1 .

79 ASTRONOMY AND ASTROPHYSICS↗

Proton beam power limits for stationary water-cooled tungsten target with different cladding materials

The proton beam power limit for a solid-tungsten spallation target is largely determined by beam induced thermomechanical structural loads and decay heat power deposition, while its lifetime is limited by radiation damage and fatigue life of the target materials. In this paper, we studied the power limits of a stationary water-cooled solid tungsten target concept. Tantalum clad tungsten was considered as a reference case. Being a low activation material, zircaloy 2 cladding option was studied and its decay heat driven power limit was compared with the reference case. Zirconium alloys have proven operations records in spallation target and nuclear fission environments, supported by materials data obtained from post irradiation examinations. Recent study also demonstrated feasibility of diffusion bonding zirconium to tungsten using vanadium foil inter layer. Particle transport simulations code FLUKA was used to calculate energy deposition and decay heat power deposition in the target, based on the beam parameters technically feasible at the Second Target Station of the Spallation Neutron Source at Oak Ridge National Laboratory. The energy deposition data were used for flow, thermal, and structural analyses to determine the beam intensity limit on the target concept studied. The decay heat deposition data were used to calculate the transient temperature evolution in the tungsten volumes in a loss of coolant accident (LOCA) scenario to determine its beam power limit. For a 1.3 GeV proton beam, the power limit on a stationary target was 400 kW for a tantalum clad target model and 800 kW for a zircaloy 2 clad target model.

Lee, Yong Joong↗

Performance Evaluation of Drain Water Heat Recovery Exchangers for Heat Pump Water Heaters

Water heating constitutes approximately 18% of energy consumption and is the second largest energy expenses in United States homes. Heat pump water heaters (HPWHs) are energy efficient technologies with lower carbon footprints as compared to conventional water heating technologies, such as gas and electrical resistance heaters. The performance of water-source HPWHs can be improved by recovering heat from blackwater using drain heat recovery heat exchangers. Depending on water draw patterns of single family and multifamily residences, the operation and heat transfer performance of these heat exchangers are highly transient. This paper examines the thermo-hydraulic performance of a drain heat recovery heat exchanger during its transient and steady-state operations. The heat recovery performance of the exchanger was evaluated at different inlet temperatures and flow rates ranging from 9°C to 36°C and 0.03 kg/s to 0.28 kg/s, respectively. The obtained effectiveness varies from 45%–85% depending on the operating conditions. The test facility, steady-state and transient experimental procedures are reported in detail. Further, the effect of operating conditions on the effectiveness is discussed. The experimental approach and results of the study will provide insights into the transient operation of drain recovery heat exchangers as well as, guide sizing for various steady state conditions.

Krishnan, Easwaran↗

Early Research in Load-Following Management for HPC-Nuclear Integration

With the rising demand for high performance computing (HPC) and artificial intelligence (AI) systems, maintaining a stable and efficient power supply is increasingly critical. The HPC team at Idaho National Laboratory is spearheading efforts to seamlessly integrate HPC systems with nuclear reactors. This lightning talk explores one early strategy for managing power fluctuations using software-defined controls. To effectively harness nuclear reactors for power generation, control mechanisms are essential to address the slow load-following capabilities of reactors, which are typically around 5% per minute. While this rate is sufficient for many uses, large HPC systems can experience rapid power consumption changes by tens of megawatts when jobs start or stop running. A reactor could overproduce power and match the peak power rating for the HPC system, however when the system is not running a job or a job unexpectedly stops, the load-following of the system would be affected leading to power being wasted and the likelihood of power transient occurrences increases. Controlling the increase or decrease of power consumption on these systems at the same rate as the load-following of reactors is one piece of the puzzle to properly utilizing nuclear reactors as a power source for HPC systems.

97 - MATHEMATICS AND COMPUTING↗

Validation of URANS and STRUCT-ε turbulence models for stratified sodium flow

Simulations of a transient stratified sodium experiment are carried out using a classic unsteady RANS model and the second-generation URANS model, STRUCT-ε. Turbulence modeling challenges and their implications to stratified flow prediction are discussed in the context of other sources of error. Input errors are discussed and addressed; discretization error is calculated to be less than 5% of the inlet velocity, for 80% of the domain; and remaining errors in temperature distributions are attributed to the turbulence model. Qualitative flow features from the simulations are presented and discussed. Compared to the experiment, the STRUCT-ε turbulence model provides a more physically accurate prediction of temperature and momentum mixing in key regions of the domain. Quantitative measures such as the L 2 norm of the temperature discrepancy demonstrate the improved performance of the STRUCT-ε approach. In conclusion, the magnitude of the temperature fluctuations is very well-predicted by the STRUCT-ε, while URANS overpredicts them by approximately 50%.

42 ENGINEERING↗

Sub-Diurnal Methane Variations on Mars Driven by Barometric Pumping and Planetary Boundary Layer Evolution

In recent years, the Tunable Laser Spectrometer within the Sample Analysis at Mars (TLS-SAM) instrument on board the Mars Science Laboratory (MSL) Curiosity rover has detected methane variations in the atmosphere at Gale crater. Methane concentrations appear to fluctuate seasonally as well as sub-diurnally, which is difficult to reconcile with an as-yet-unknown transport mechanism delivering the gas from underground to the atmosphere. To potentially explain the fluctuations, we consider barometrically induced transport of methane from an underground source to the surface, modulated by temperature-dependent adsorption. The subsurface fractured-rock seepage model is coupled to a simplified 1-D atmospheric mixing model to provide insights on the pattern of atmospheric methane concentrations in response to transient surface methane emissions, as well as to predict sub-diurnal variation in methane abundance for the northern summer period, which is a candidate time frame for a MSL Curiosity sampling campaign. Our analysis suggests that there is a lower limit to the subsurface fracture density that can produce the observed methane patterns, below which the atmospheric methane variations would be out of phase with the observations. The best-performing model scenarios indicate a significant, short-lived methane pulse just prior to sunrise, the detection of which by TLS-SAM would be a potential indicator of the contribution of barometric pumping to Mars' atmospheric methane variations.

54 ENVIRONMENTAL SCIENCES↗

Operability of a Power System with Synchronous Condensers and Grid-Following Inverters

Growing shares of inverter-based resources gener-ally correlate with a reduction in inertia as synchronous genera-tors are displaced. Along the path to high shares of inverter-based resources, in particular with only conventional grid-following inverter controls, a proffered solution is the use of synchronous condensers as a source of inertia to maintain the frequency-power balance relationship. An outstanding question is whether only grid-following inverters and synchronous condensers yield a viable power system; i.e., all frequency response is derived from inverters. A validated electromagnetic transient model of the Maui system, with many nonlinear elements such as load-shedding, line tripping, and inverter ride through criteria disabled, is used to investigate the stability of such a system. Two types of perturbations were applied, a 15% generation loss and a fault event, and it was found that the system remains stable.

grid-following inverters↗

Electromagnetic Analysis of ITER Electron Cyclotron Emission Model

The ITER electron cyclotron emission (ECE) diagnostic system is located at diagnostics shielding module 2 (DSM2), equatorial port 9 (EP9), to measure electron temperature profile and electron temperature fluctuations, and also assess nonthermal electron distributions via the oblique view. Therefore, ECE has both radial and oblique views with two hot sources for calibration and two couples of mirrors for two different optical views. It shall receive that the large electromagnetic (EM) loads up to 100-MN/m3 force density due to the eddy currents on the copper mirrors generated by the short transient plasma disruption. The simplified EM analysis model using magnetic field ( B ) data and flux variations (dB/dt) method based on the worst case of plasma disruption, MD_DW_EXP16MS_CATIII, has been developed to calculate the force and moments for each ECE component and support structure for EP9 DSM2 Bay2 and Bay3. The B and dB/dt methods use the constant B and dB/dt assumed for each local bay space, but the real B and dB/dt are basically varied along the radial direction of each bay. The complex global model with the local ECE components and support structure using Maxwell transient has been developed to do EM analysis to compare the results of the B and dB/dt models. The volumetric EM force density of the whole ECE components and support structure can be used for the subsequent structural analysis. In conclusion, combined with the thermal and nuclear loads and seismic and initial loads, the ECE integration analysis can be finalized for the operation case.

Fang, J. [Princeton Plasma Physics Laboratory (PPP↗