Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transients (Dynamics)”

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

Development and Porting of Nuclear Reactor Computational Models for the NRIC Virtual Test Bed in FY23

The U.S. DOE’s National Reactor Innovation Center (NRIC) hosts an open-source website and associated GitHub repository called the Virtual Test Bed (VTB) on which computational models for advanced reactors are documented and shared with the reactor community. In FY23 under the NRIC program, computational models for nuclear reactor analysis were developed and contributed to the VTB Github repository during FY23 by contributors at Argonne National Laboratory. Other models developed outside of NRIC were also ported to the VTB as part of this work, serving a broader mission to make computational reactor analysis models more widely available to the reactor community. The model development activities focused on demonstration of gas-cooled microreactor models including computational fluid dynamics simulations of flow through an industry-inspired air jacket design, and simulation of multiphysics transients for a gas-cooled microreactor assembly. Additionally, development of a molten salt reactor (led at Idaho National Laboratory) was supported through Argonne’s expertise in multigroup cross sections generation. Models for these two reactor types were targeted due to their relevancy to NRIC’s Demonstration of Microreactor Experiments (DOME) and Laboratory for Operation and Testing in the U.S. (LOTUS) physical test beds, which are slated to host microreactor and molten salt reactor experiments. The model porting activities consisted of developing detailed documentation for several physics models which originate from the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. We include only contributions from Argonne National Laboratory in this report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

PV-Fed Smart Inverters for Mitigation of Voltage and Frequency Fluctuations in Islanded Microgrids

The presence of low-inertia distributed generation (DG) units makes islanded microgrids vulnerable to voltage and frequency variations. Besides, the considerable difference between the inertia of synchronous-based and inverter-based DGs results in a power imbalance between supply and demand during abnormal conditions. As a result, both voltage and frequency of microgrid ac-bus start oscillating which might have severe consequences such as blackouts. This paper deploys the traditional controller of photovoltaic (PV) units to improve the dynamic behavior of islanded microgrids by suppressing the voltage and frequency fluctuations. To this end, an adaptive piecewise droop (APD) characteristic is proposed and employed in PV units to attain a faster balance between generation and consumption during transients, leading to an enhanced frequency response. Besides, the reactive-power control loop is equipped with a droop characteristic which enables the PV units to inject/absorb reactive power during transients and participate in voltage-profile enhancement of the system. Case study results are presented using PSCAD/EMTDC to confirm the validity of proposed method in improving the dynamics of islanded microgrids.

Pilehvar, Mohsen S.↗

Performance Analysis and Simulaion of the Hydraulic Scram System in TREAT Reactor

The Transient Reactor Test Facility (TREAT) at Idaho National Laboratory (INL) serves a vital role in nuclear fuel safety research, enabling transient experiments that simulate reactivity excursions and accident scenarios. Central to these operations is the transient control rod drive system (TCRDS), which drives rapid motion of the transient control rods such that TREAT can simulate rapid power changes typical of reactor accidents. The reliability and performance of this system are critical for protecting both fuel specimens and reactor infrastructure. This study presents the initial phase of a two-year investigation into the dynamics and reliability of the TREAT hydraulic TCRDS. Conducted in collaboration with INL, the research employs a combined computational and experimental approach to analyze the system's response time, pressure transients, and potential failure modes. Emphasis is placed on understanding how fluid characteristics influence the TCRDS’s ability to achieve both rapid power changes and mechanical stability. The TRDS and the skid that powers it will be analyzed throughout this investigation. Computational modeling using computational fluid dynamics (CFD) will simulate the hydraulic response under varying conditions. In parallel, experimental testing planned at INL will validate these models and capture key performance metrics. This paper outlines the system design, analytical framework, and modeling strategies that form the foundation for later testing. Ultimately, this work aims to support improvements to the TCRDS’s design and reliability, contributing to the broader goal of enhancing nuclear fuel safety and sustaining TREAT’s mission as a premier nuclear fuel test facility.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Electronic Tuning of Photoexcited Dynamics in Heteroleptic Cu(I) Complex Photosensitizers

Photoexcited dynamics of heteroleptic Cu(I) complexes as noble-metal-free photosensitizers are closely intertwined with the nature of their ligands. By utilizing ultrafast optical and X-ray transient absorption spectroscopies, we characterized a new set of heteroleptic Cu(I) complexes [Cu(PPh 3 )( 2 )(BPy R )] + (R = CH 3 , H, Br to COOCH 3 ), with an increase in the electron-withdrawing ability of the functional group (R). Here we found that after the transient photooxidation of Cu(I) to Cu(II), the increasing electron-withdrawing ability of R barely affects the internal conversion (IC) (e.g., Jahn-taller (JT) distortion) between singlet MLCT states. However, it does accelerate the dynamics of intersystem crossing (ISC) between singlet and triplet MLCT states and the subsequent decay from the triplet MLCT state to the ground state. The associated lifetime constants are reduced by up to 300%. Our understanding of the photoexcited dynamics in heteroleptic Cu(I) complexes through ligand electronic tuning provides valuable insight into the rational design of efficient Cu(I) complex photosensitizers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine Learning Assisted Safety Modeling and Analysis of Advanced Reactors

With the advances in computational power and numerical methods, analysts can now rely on first-principle simulations to predict ultra-fine details in a variety of applications. Advances in machine learning (ML) have produced algorithms that can now learn high-level abstractions via hierarchical models. This project aims to leverage advances in ML techniques and the available high-resolution simulation data to develop a novel modeling and simulation (M\&S) methodology for reactor safety analysis. While application-agnostic ML techniques are available, complex physics constraints need to be incorporated into ML techniques to build ML-based closures for computationally efficient predictive simulations. This project intends to develop a physics-guided data-driven multi-scale methodology for M\&S of advanced reactors. The project focuses on thermal fluid (T/F) phenomena, which play major roles in advanced reactor safety. Specifically, we propose a data-driven coarse-mesh turbulence model based on local flow features for the transient analysis of thermal mixing and stratification in a sodium-cooled fast reactor (SFR). The model has a coarse-mesh setup to ensure computational efficiency, while it is trained by fine-mesh computational fluid dynamics (CFD) data with Reynolds-averaged Navier-Stokes (RANS) turbulence model to ensure accuracy. Three different neural networks are developed and tested for loss-of-flow transients in the hot pool of SFR, i.e. the densely connected convolutional neural network (DCNN), long-short-term-memory network based on proper orthogonal decomposition (POD-LSTM), and the DCNN informed by LSTM (DCNN-LSTM). The performances of these three neural networks are evaluated based on baseline models. The DCNN-LSTM model has been chosen for further hyperparameter optimization. Furthermore, based on a simplified two-dimensional case, uncertainty quantification (UQ) of the developed ML-based closure are investigated with three methods, i.e. Monte Carlo dropout, deep ensemble, and Bayesian neural network. The developed ML-based turbulent viscosity closure relation based on deep ensemble is then integrated into the system analysis module SAM and serves as a term in the conservation equations. Such a SAM-ML based procedure guarantees that the obtained results are consistent with the physical constraints of the thermal-fluid system. The SAM-ML simulation on the same loss-of-flow transient showed comparable accuracy with the CFD simulation but with a much coarser mesh setup. Last but not least, the ML-based closure improvement with the support of higher-fidelity data from large eddy simulation (LES) is discussed. As a first step towards this direction, a baseline LES simulation is performed to obtain comparable data with RANS results. Based on the early results, future investigation on further improving the ML-based closure is discussed. We believe the developed approach that combines scientific machine learning with nuclear system analysis code can benefit the advanced reactor community as more accurate safety analyses will better characterize reactor safety margins and reduce licensing efforts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Virtual Self-Excited Induction Generator-Based Grid-Forming Inverter Control for Robust Voltage Regulation Under Nonideal Loading

This paper presents a generator-inspired control methodology for grid-forming (GFM) inverters that deliberately emulates a self-excited induction generator so that the inverter can hold its voltage and frequency under difficult loading and severe terminal disturbances across wide voltage and frequency ranges. The design integrates a Lyapunov energy function-based inner loop to provide high bandwidth and strong disturbance rejection, and it complements this with a passivity-based argument that furnishes a coherent large-signal stability guarantee beyond small-signal limits. Analytical insights are developed via the Krylov-Bogoliubov-Mitropolsky averaging method, which reveals an intrinsic resistive droop characteristic; these closed-form relations both explain the observed dynamics and yield simple, decentralized tuning rules. The methodology is validated on a controller-hardware-in-the-loop platform and exercised in real time across balanced, unbalanced, and nonlinear loads, as well as during parallel operation. Across these scenarios, the inverter maintains balanced three-phase voltages, limits harmonic content, settles quickly with well-damped transients, and remains resilient when multiple units operate in parallel. The contributions are a self-excited-machine-inspired GFM controller with enhanced dynamic performance and robustness, a single stability rationale grounded in passivity, closed-form expressions that guide tuning, and comprehensive hardware-in-the-loop validations demonstrating effectiveness and superiority under challenging operating conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Terahertz-Driven Local Dipolar Correlation in a Quantum Paraelectric

Light-induced ferroelectricity in quantum paraelectrics is a new avenue of achieving dynamic stabilization of hidden orders in quantum materials. In this Letter, we explore the possibility of driving a transient ferroelectric phase in the quantum paraelectric KTaO 3 via intense terahertz excitation of the soft mode. We observe a long-lived relaxation in the terahertz-driven second harmonic generation (SHG) signal that lasts up to 20 ps at 10 K, which may be attributed to light-induced ferroelectricity. Through analyzing the terahertz-induced coherent soft-mode oscillation and finding its hardening with fluence well described by a single-well potential, we demonstrate that intense terahertz pulses up to 500 kV=cm cannot drive a global ferroelectric phase in KTaO 3 . Instead, we find the unusual long-lived relaxation of the SHG signal comes from a terahertz-driven moderate dipolar correlation between the defect-induced local polar structures. We discuss the impact of our findings on current investigations of the terahertz-induced ferroelectric phase in quantum paraelectrics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Permutationally Invariant Polynomial Expansions with Unrestricted Complexity

A general strategy is presented for constructing and validating permutationally invariant polynomial (PIP) expansions for chemical systems of any stoichiometry. Demonstrations are made for three categories of gas-phase dynamics and kinetics: collisional energy-transfer trajectories for predicting pressure-dependent kinetics, three-body collisions for describing transient van der Waals adducts relevant to atmospheric chemistry, and nonthermal reactivity via quasiclassical trajectories. In total, 30 systems are considered with up to 15 atoms and 39 degrees of freedom. Permutational invariance is enforced in PIP expansions with as many as 13 million terms and 13 permutationally distinct atom types by taking advantage of petascale computational resources. The quality of the PIP expansions is demonstrated through the systematic convergence of in-sample and out-of-sample errors with respect to both the number of training data and the order of the expansion, and these errors are shown to predict errors in the dynamics for both reactive and nonreactive applications. Here, the parallelized code distributed as part of this work enables the automation of PIP generation for complex systems with multiple channels and flexible user-defined symmetry constraints and for automatically removing unphysical unconnected terms from the basis set expansions, all of which are required for simulating complex reactive systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Barrier-Certificated Reinforcement Learning Approach for Enhancing Power System Transient Stability

Increasing integration of renewable resources brings more flexibility and poses new challenges to modern power systems, leading to highly nonlinear and complex dynamics. Here, this paper aims to provide a general solution framework to traditional control problems, such as frequency control and voltage control, which attempt to maintain the stability of either synchronous generators-governed or inverter-governed systems when subjected to a disturbance and simultaneously guarantee operational constraints, providing a complete complement to existing works on control design. Building on reinforcement learning (RL) and control barrier functions, the framework includes two subsystems, i.e., a model-free controller and a barrier-certification system, which discover RL-based control actions and sequentially filter them using a barrier certificate to satisfy operational constraints. Calculating a barrier function is generally challenging for a complex power system. This is addressed by representing the barrier function using neural networks (NNs) and data-based approaches. An adaptive method is introduced to certify the neural barrier function that perseveres barrier conditions, which is more compatible with online implementation. The proposed framework synthesizes a stabilizing controller that satisfies predefined safety regions. The effectiveness of the proposed framework is demonstrated via several comparative case studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Forward Feedback Control Scheme for a Solar Thermochemical Moving Bed Counter-Current Flow Reactor

Abstract Pelletized thermochemical energy storage media has a potential for long-duration energy storage. Production of solid-state energy storage media can be done within a cavity chemical reactor that captures concentrated solar radiation from a solar thermal field. The temperature stability of a solar reactor is directly influenced by the solar flux intercepted. This paper presents a low-order physical model to simulate the dynamic response of temperature inside a tubular plug-flow reactor prototype. Solid granular particles are fed to the reactor from the top whereas a counter-current flowing gas enters the reactor from the bottom. An in-house code was developed to model transient heat transfer of the reactor wall, gas, and moving particles. The model was preliminarily validated with packed beds for different temperature ranges and two gas flowrates. Dynamic response of the reactor temperature is simulated for different input power and gas/particle flowrates. The results show that the system response can be controlled efficiently by utilizing input power (solar flux) as a control parameter. A conventional proportional integral (PI) controller is designed to control the temperature inside the reactor and to maintain it during the solar flux intermittency. The controller parameters are tuned using the Ziegler–Nichols method to ensure optimal system response. The results show that the feedback control model is successful in tracking different reference reactor temperatures within a reasonable settling time of 30 min and eliminated overshoot. This study can be extended to include a hybrid reactor with a multi-input, multi-output variable system.

Energy & Fuels↗

Origin of layered perovskite device efficiencies revealed by multidimensional time-of-flight spectroscopy

Mixtures of layered perovskite quantum wells with different sizes form prototypical light-harvesting antenna structures in solution-processed films. Gradients in the bandgaps and energy levels are established by concentrating the smallest and largest quantum wells near opposing electrodes in photovoltaic devices. Whereas short-range energy and charge carrier funneling behaviors have been observed in layered perovskites, our recent work suggests that such light-harvesting processes do not assist long-range charge transport due to carrier trapping at interfaces between quantum wells and interstitial organic spacer molecules. Here, we apply a two-pulse time-of-flight technique to a family of layered perovskite systems to explore the effects that interstitial organic molecules have on charge carrier dynamics. In these experiments, the first laser pulse initiates carrier drift within the active layer of a photovoltaic device, whereas the second pulse probes the transient concentrations of photoexcited carriers as they approach the electrodes. The instantaneous drift velocities determined with this method suggest that the rates of trap-induced carrier deceleration increase with the concentrations of organic spacer cations. Overall, our experimental results and model calculations suggest that the layered perovskite device efficiencies primarily reflect the dynamics of carrier trapping at interfaces between quantum wells and interstitial organic phases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulation of the supercritical CO 2 recompression Brayton power cycle with a high-temperature regenerator

The supercritical carbon dioxide (sCO 2 ) recompression Brayton cycle promises higher efficiency and lower capital cost than traditional steam Rankine power cycles. However, achieving high efficiency requires large, highly effective recuperators. Regenerators may be a low-cost alternative to printed circuit and micro-tube heat exchangers for recuperation in sCO 2 power cycles. Regenerators are a periodic heat exchanger in which thermal energy is extracted from the hot stream, stored in solid media, and then released to the cold stream at a later time. Fixed bed regenerators with valves to direct fluid are the preferred method for implementing regenerators in power cycles, but the inherently transient nature of these systems has not been characterized for this application. This study presents the simulation of a high-temperature regenerator within a 10 MWe sCO 2 recompression Brayton cycle. A transient, one-dimensional regenerator model presented in a previous study is used to simulate the regenerator. Dynamic heat exchanger models are also developed for the precooler, low-temperature recuperator, and primary heat exchanger, and the compressors and the turbine are modeled with off-design performance maps. We assess two regenerator-valve design options; one for fast switching, and one for reduced flow rate fluctuations. System simulation finds that both designs see significant fluctuations in turbomachinery flow rate, turbomachinery and system power, and regenerator discharge process outlet temperature. While designing the regenerator-valve subsystem for lower fluctuations is possible, the regenerator cold discharge temperature and net power still fluctuate by ±77.6°C and 6%, respectively. Increasing buffer volume is not effective at sufficiently reducing these fluctuations, but adding a packed bed in between the regenerator and the primary heat exchanger can reduce regenerator discharge process outlet temperature fluctuations to 6.4°C. Further reductions could be possible by increasing the size of this packed bed.

42 ENGINEERING↗

River Dynamics Control Transit Time Distributions and Biogeochemical Reactions in a Dam-Regulated River Corridor

Residence Time Distributions (RTDs) exert an important control on biogeochemical translation in the hydrological systems. Previous tracer studies have revealed that RTDs often followed time-invariant exponential, lognormal, or heavy-tailed RTDs that have power-law behavior for long tails in headwater or low-order streams. However, there is a recognition that RTDs can be more complicated and time-variable in response to dynamic hydrological forcing. Here, we use particle tracking to estimate RTDs along the Hanford Reach of the Columbia River and to quantify the influences of river stage fluctuations. Particle tracking is conducted using the velocity field from high-resolution 3D groundwater flow simulations. The effects of hydrological forcing on the residence time distribution were evaluated by varying river flow boundary conditions and releasing particles in different time windows. Our results revealed that dynamic stage fluctuations created rapidly changing losing-gaining conditions in the river and led to highly transient RTDs, which contributes to multiple modes of RTDs. Dam-induced high-frequency (sub-daily) flow variation contributes to the short-time (sub-daily) component of the RTDs. Deviation of the reactant consumption under the single-mode assumption compared to the multimodal RTDs is relatively small (~5%) and appears when the Damköhler number was close to one.More specifically, high-frequency (daily) dam-induced stage variations have higher impacts on biogeochemical reactions with faster reaction rate, and can potentially increase the reactant consumption by 26.68%. Based on these findings, we suggest that current river basin models could be improved by including bank storage and more complex RTDs influenced by both short and long term river stage fluctuations.

54 ENVIRONMENTAL SCIENCES↗

Scaling Behavior of Thermally Driven Fractures in Deep Low‐Permeability Formations: A Plane Strain Model With 1‐D Heat Conduction

Abstract Injection of cold fluids through/into deep formations may cause significant cooling, thermal stress, and possible thermal fracturing. In this study, the thermal fracturing of low‐permeability formations under one‐dimensional heat conduction was investigated using a plane strain model. Dimensionless governing equations, with dimensionless fracture length , aperture , spacing , time , and effective confining stress , were derived. Solution of single thermal fracture was derived analytically, while solution of multiple fractures with constant (or dynamic) spacing were obtained using the displacement discontinuity method (and stability analysis). For single fracture, increases nonlinearly with and then transitions to scaling law , indicating that late‐time fracture length increases linearly with the square root of cooling time. For constantly spaced fractures, deviates from the single‐fracture solution at a later for a larger , showing slower propagation under inter‐fracture stress interaction . For dynamically spaced fractures, fracture arrest induced by stress interaction was determined by the stability analysis; the fully transient solution provides evolution of dimensionless fracture length, spacing, aperture, and pattern; a similar scaling law, with , obtained shows the effect of both stress interaction and fracture arrest. The solution and scaling law provide fast predictions for all reservoir and cooling conditions using (single) model parameter . Application to a geothermal site with demonstrates that thermal fractures reach 0.67, 6.25, and 78.00 m in length, 0.49, 2.30, and 13.00 m in spacing, and 0.43, 2.09, and 12.19 mm in aperture at 1, 100, and 10,000 days.

Chen, Bin↗

Large Divergence of Projected High Latitude Vegetation Composition and Productivity Due To Functional Trait Uncertainty

Abstract Vegetation distribution and composition are expected to change in northern high latitudes under rapid warming, which regulates ecosystem functions but remains challenging to predict. Vegetation change arises from the interplay of chronic climate trends such as warming and transient demographic processes of recruitment, growth, competition, and mortality. Most predictive models overlooked the role of demographic dynamics controlled by plant traits. Here, we simulate vegetation dynamics at the Kougarok Hillslope site in Alaska under historical and future climates using the E3SM Land Model coupled to the Functionally Assembled Terrestrial Simulator (ELM‐FATES). To evaluate the roles of plant traits, we parameterize the model with 5,265 trait configurations representing diverse physiological and demographic strategies. Results show current modeled biomass, composition, and productivity are most sensitive to traits controlling photosynthetic capacity, carbon allocation, allometry, and phenology. Among all trait configurations, ∼5% reproduce in situ biomass and plant functional type (PFT) composition measured in 2016, that are indistinguishable from these two observed ecosystem states. Notably, these same trait configurations produce diverging biomass, composition, and productivity under future climate, where the uncertainty attributable to traits is twice the change attributable to climate change. The variation of projected productivity arises from emerging PFT composition under novel climate regimes, primarily explained by traits controlling cold‐induced mortality, recruitment, and allometry. Our findings highlight the importance and uncertainty of demographic dynamics and its interaction with climate change in shaping Arctic vegetation change. Improved model predictions will likely benefit from explicit consideration of vegetation demography and better constraints of critical traits.

54 ENVIRONMENTAL SCIENCES↗

Polarization-dependent photoemission electron microscopy for domain imaging of inorganic and molecular materials

Polarization-dependent photoemission electron microscopy (PD-PEEM) exploits spatial variation in the optical selection rules of materials to image domain formation and material organization on the nanoscale. In this Perspective, we discuss the mechanism of PD-PEEM that results in the observed image contrast in experiments and provide examples of a wide range of material domain structures that PD-PEEM has been able to elucidate, including molecular and polymer domains, local electronic structure and defect symmetry, (anti)ferroelectricity, and ferromagnetism. As a result, we discuss challenges and new directions that are possible with this tool for probing domain structure in materials, including investigating the formation of transient ordered states, multiferroics, and the influence of molecular and polymer order and disorder on excited state dynamics and charge transport.

36 MATERIALS SCIENCE↗

Controller Hardware-in-the-Loop Evaluation of a Microgrid Controller for a Microgrid System with Multiple Grid-Forming Inverters

This paper presents the laboratory evaluation of a commercial Microgrid Management System (MGMS) implemented in the real-world Bronzeville Microgrid which features a futuristic scenario with high renewable energy integration and the use of multiple Grid-Forming (GFM) inverters. The primary objective of the performance evaluation for the MGMS is to assess the MGMS's capability to dispatch GFM units, including a GFM PV unit and two GFM battery units, to maintain the system stability and ensure economic operation, thus guaranteeing the microgrid's resilience during prolonged outages and dynamic events. The laboratory controller hardware-in-the-loop provides realistic testing environment through detailed electromagnetic transient modeling of the microgrid system, hardware MGMS, and standard communication protocols (DNP3). The CHIL evaluation shows how the MGMS effectively manages the GFM inverters, highlighting its performance in maintaining stability, reliability, and survivability in a microgrid environment with a high penetration of renewable energy sources.

controller hardware-in-the-loop↗

Halted-pendulum Relaxation: Application to White Dwarf Binary Initial Data

Studying compact-star binaries and their mergers is integral to determining progenitors for observable transients. Today, compact-star mergers are typically studied via state-of-the-art computational fluid dynamics codes. One such numerical technique, smoothed particle hydrodynamics (SPH), is frequently chosen for its excellent mass, energy, and momentum conservation. The natural treatment of vacuum and the ability to represent highly irregular morphologies make SPH an excellent tool for the study of compact-star binaries and mergers. For many scenarios, including binary systems, the outcome of simulations is only as accurate as the initial conditions. For SPH, it is essential to ensure that the particles are distributed regularly, representing the initial density profile but without long-range correlations. Particle noise in the form of high-frequency local motion and low-frequency global dynamics must be damped out. Damping the latter can be as computationally intensive as the actual simulation. We discuss a new and straightforward relaxation method, halted-pendulum relaxation (HPR), to remove global oscillation modes of SPH particle configurations. In combination with effective external potentials representing gravitational and orbital forces, we show that HPR has an excellent performance in efficiently relaxing SPH particles to the desired density distribution and removing global oscillation modes. We compare the method to frequently used relaxation approaches and test it on a white dwarf binary model at its Roche-lobe overflow limit. We highlight the importance of our method in achieving accurate initial conditions and its effect on achieving circular orbits and realistic accretion rates when compared with other general relaxation methods.

79 ASTRONOMY AND ASTROPHYSICS↗