Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “current conservation”

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 289 records · Page 16

Numerical challenges for energy conservation in N -body simulations of collapsing self-interacting dark matter halos

Dark matter (DM) halos can be subject to gravothermal collapse if the DM is not collisionless, but engaged in strong self-interactions instead. When the scattering is able to efficiently transfer heat from the centre to the outskirts, the central region of the halo collapses and reaches densities much higher than those for collisionless DM. This phenomenon is potentially observable in studies of strong lensing. Current theoretical efforts are motivated by observations of surprisingly dense substructures. However, a comparison with observations requires accurate predictions. One method to obtain such predictions is to use N-body simulations. Collapsed halos are extreme systems that pose severe challenges when applying state-of-the-art codes to model self-interacting dark matter (SIDM). In this work, we investigate the root of such problems, with a focus on energy non-conservation. Moreover, we discuss possible strategies to avoid them. We ran N-body simulations, both with and without SIDM, of an isolated DM-only halo and we adjusted the numerical parameters to check the accuracy of the simulation. We find that not only the numerical scheme for SIDM can lead to energy non-conservation, but also the modelling of gravitational interaction and the time integration are problematic. The main issues we find are: (a) particles changing their time step in a non-time-reversible manner; (b) the asymmetry in the tree-based gravitational force evaluation; and (c) SIDM velocity kicks breaking the time symmetry. Tuning the parameters of the simulation to achieve a high level of accuracy allows us to conserve energy not only at early stages of the evolution, but also later on. However, the cost of the simulations becomes prohibitively large as a result. Some of the problems that make the simulations of the gravothermal collapse phase inaccurate can be overcome by choosing appropriate numerical schemes. However, other issues still pose a challenge. Our findings motivate further works on addressing the challenges in simulating strong DM self-interactions.

dark matter↗

Combined Preoperative Hypofractionated Radiotherapy With Doxorubicin-Ifosfamide Chemotherapy in Marginally Resectable Soft Tissue Sarcomas: Results of a Phase 2 Clinical Trial

There is no standard treatment for marginally resectable soft tissue sarcomas (STSs) of the extremities and trunk wall, and current approaches produce unsatisfactory results. We hypothesized that the combination of doxorubicin-ifosfamide (AI) chemotherapy and 5 × 5 Gy hypofractionated radiotherapy can generate a higher ratio of limb-sparing or conservative surgeries with negative microscopic margins (R0) and acceptable treatment toxicity.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

A conservative approach to scaling magneto-inertial fusion concepts to larger pulsed-power drivers

The Magnetized Liner Inertial Fusion (MagLIF) experimental platform [M. R. Gomez et al., Phys. Rev. Lett. 113, 155003 (2014)] represents the most successful demonstration of magneto-inertial fusion (MIF) techniques to date in pursuit of ignition and significant fusion yields. The pressing question remains regarding how to scale MIF concepts like MagLIF to more powerful pulsed-power drivers while avoiding significant changes in physical regimes that could adversely impact performance. In this work, we propose a conservative approach for scaling general MIF implosions, including MagLIF. Underpinning our scaling approach is a theoretical framework describing the evolution of the trajectory and thickness of a thin-walled, cylindrical, current-driven shell imploding on preheated, adiabatic fuel. By imposing that scaled implosions remain self-similar, we obtain a set of scaling rules expressing key target design parameters and performance metrics as functions of the maximum driver current I max . We identify several scaling paths offering unique, complementary benefits and trade-offs in terms of physics risks and driver requirements. Remarkably, when scaling present-day experiments to higher coupled energies, these paths are predicted to preserve or reduce the majority of known performance-degrading effects, including hydrodynamic instabilities, impurity mix, fuel energy losses, and laser-plasma interactions, with notable exceptions clearly delineated. In the absence of α heating, our scaling paths exhibit neutron yield per-unit-length scaling as $\tilde{Y}$∝[I$3\atop{max}$,I${4.14}\atop{max}$] and ignition parameter scaling as χ∝[I max ,I${2.14}\atop{max}$]. By considering the specific physics risks unique to each scaling path, we provide a roadmap for future investigations to evaluate different scaling options through detailed numerical studies and scaling-focused experiments on present-day facilities. Overall, these results highlight the potential of MIF as a key component of the national ignition effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Using action space clustering to constrain the recent accretion history of Milky Way-like galaxies

ABSTRACT In the currently favoured cosmological paradigm galaxies form hierarchically through the accretion of satellites. Since a satellite is less massive than the host, its stars occupy a smaller volume in action space. Actions are conserved when the potential of the host halo changes adiabatically, so stars from an accreted satellite would remain clustered in action space as the host evolves. In this paper, we identify recently disrupted accreted satellites in three Milky Way-like disc galaxies from the cosmological baryonic FIRE-2 simulations by tracking satellites through simulation snapshots. We try to recover these satellites by applying the cluster analysis algorithm Enlink to the orbital actions of accreted star particles in the z = 0 snapshot. Even with completely error-free mock data we find that only 35 per cent (14/39) satellites are well recovered while the rest (25/39) are poorly recovered (i.e. either contaminated or split up). Most (10/14 ∼70 per cent) of the well-recovered satellites have infall times <7.1 Gyr ago and total mass >4 × 108M⊙ (stellar mass more than 1.2 × 106 M⊙, although our upper mass limit is likely to be resolution dependent). Since cosmological simulations predict that stellar haloes include a population of in situ stars, we test our ability to recover satellites when the data include 10–50 per cent in situ contamination. We find that most previously well-recovered satellites stay well recovered even with 50 per cent contamination. With the wealth of 6D phase space data becoming available we expect that cluster analysis in action space will be useful in identifying the majority of recently accreted and moderately massive satellites in the Milky Way.

79 ASTRONOMY AND ASTROPHYSICS↗

Piezosuperconductivity: Novel effects in noncentrosymmetric superconductors

Here we study effects in noncentrosymmetric superconductors arising from their unique coupling of Cooper-pair condensate and elasticity. We show that although the much discussed Lifshitz coupling is not observable in a uniform bulk state, it strikingly endows dislocations with a fractional magnetic flux. We also predict a generation of voltage-free strain by a DC current in a P- and T-breaking Josephson junction. Viewing superconductors through the lens of higher-form symmetries we identify the Lifshitz coupling as a chemical potential for the approximately conserved winding number, drawing an analogy with pyroelectric insulators.

36 MATERIALS SCIENCE↗

Field‐scale analysis of miscanthus production indicates climate change may increase the opportunity for water quality improvement in a key Iowa watershed

Abstract The Raccoon River Basin is the primary source for drinking water in Iowa's largest city and plays a major role in the Mississippi River Basin's high nutrient exports. Future climate change may have major impacts on the biological, physiological, and agronomic processes imposing a threat to ecosystem services. Efforts to reduce nitrogen (N) loads within this basin have included local litigation and the implementation of the Iowa Nutrient Reduction Strategy, which suggest incorporating bioenergy crops (i.e., miscanthus) within the current corn–soybean landscape to reach a 41% reduction in nitrate loads. This study focuses on simulating N export for historical and future land use scenarios by using an agroecosystem model (Agro‐IBIS) and a hydrology model (THMB) at the 500‐m resolution, similar to the scale of agricultural fields. Model simulations are driven by CMIP5 climate data for historical, mid‐century, and late‐century under the RCP 4.5 and 8.5 warming projections. Using recent crop profit analyses for the state of Iowa, profitability maps were generated and nitrogen leaching thresholds were used to determine where miscanthus should replace corn–soybean area to maximize reductions in N pollution. Our results show that miscanthus inclusion on low profit and high N leaching areas can result in a 4% reduction of N loss under current climate conditions and may reduce N loss by 21%–26% under future climate conditions, implying that water quality has the potential continue to improve under future climate conditions when strategically implemented conservation practices are included in future farm management plans.

54 ENVIRONMENTAL SCIENCES↗

Volume Change from Solidification Correlation for Hanford Tank Waste.

The U.S. Department of Energy (DOE), Hanford Field Office’s primary mission is to safely and effectively treat Hanford’s tank waste and deliver environmental remediation. Mixed radioactive waste is stored in the underground tanks at the Hanford Site. It was recently estimated that retrieval of the waste in the 200 West Area underground tanks in the SY, S, SX, and U Tank Farms will result in about 41 million gallons of mixed low-level waste (MLLW) (RPP-RPT-65147, Rev. 1). The current plan is to retrieve at least 22 S, SX, and U Farm tanks and pretreat1 to produce pretreated tank waste (PTW) which will be further treated (including solidification /immobilization) for Resource Conservation and Recovery Act (RCRA) Land Disposal Restriction (LDR) organics and inorganics before being transferred to an offsite out-of-state facility for disposal (RPP-PLAN-66135, Rev. 2). In addition, solidified PTW from the 200 East Area may also be transferred to an offsite out-of-state facility for disposal to ensure availability of critical Double-Shell Tank space, meet retrieval obligations, and optimize 200 East Area operations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

An Innovative Circuit Breaker for Renewable Power Grids

A key step towards a carbon-neutral future could be reached through dispersed power grids, featuring networks of local-scale renewable energy and battery storage plants. To prevent these power grids from damaging themselves and their surroundings when electrical faults arise, they must be integrated with ‘circuit breakers’, which temporarily interrupt the current flowing through them. However, currently available circuit breakers cannot handle the medium-voltage direct currents best suited for these grids. Through an innovative new circuit breaker design, Steve Schmalz and his colleagues at Eaton Corporation, the Illinois Institute of Technology (IIT), Virginia Tech and the National Energy Technology Laboratory (NETL), hope that this challenge will soon be overcome.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Spark Channel Dynamics of Electrostatic Discharges

When two differently-charged objects are brought in close proximity to each other, the resulting high electric fields can cause electron avalanche breakdown of the air gap separating the objects, a process known as electrostatic discharge (ESD). If enough initial charge is stored on the objects, the electrical breakdown can proceed to ionize the air to such a degree that a highly conductive filament of plasma forms in the gap, known as a spark channel. The spark electrically bridges the air gap, resulting in a rapid pulse of current that neutralizes the charge difference. The current pulse produces significant heating of the gas in the spark, resulting in dissociation, ionization, thermal radiation, and hydrodynamic expansion. ESD presents a hazard to electrically-sensitive devices, with consequences such as economic losses (e.g. damaged electronics) or unsafe response (e.g. unintended ignition of flammable gas mixtures, initiation of detonators, etc.). For this thesis, the ESD spark is taken to occur between two conducting electrodes, with the spark channel being axisymmetric in a cylindrical coordinate system centered on the channel. An RLC-type circuit is used for the discharge model of the ESD event. The spark is treated as a time-dependent resistance that is in series with a capacitance, an inductance, and (optionally) a load resistance representing a “victim” component under threat from the ESD event. The primary motivation of this work is to use a numerical hydrodynamic model to understand the energy dissipation and transport processes in the spark. The model consists of the compressible Euler equations of mass, momentum, and energy conservation together with an Eddington/P1 approximation for thermal radiation transport. To close the hydrodynamic system, an equation of state (EOS) was fitted from tabular data for air that accounts for the dissociation and ionization of air species. The hydrodynamic equations are solved using a conservative Lagrangian finite volume method. These partial differential equations are coupled to the circuit equations by calculation of the spark resistance via numerical integration of the electrical conductivity of the channel. Computational results are compared against experimental measurements of discharge current and radial density of the spark channel.

42 ENGINEERING↗

Overview of the SPARC tokamak

The SPARC tokamak is a critical next step towards commercial fusion energy. SPARC is designed as a high-field ( $B_0 = 12.2$ T), compact ( $R_0 = 1.85$ m, $a = 0.57$ m), superconducting, D-T tokamak with the goal of producing fusion gain $Q>2$ from a magnetically confined fusion plasma for the first time. Currently under design, SPARC will continue the high-field path of the Alcator series of tokamaks, utilizing new magnets based on rare earth barium copper oxide high-temperature superconductors to achieve high performance in a compact device. The goal of $Q>2$ is achievable with conservative physics assumptions ( $H_{98,y2} = 0.7$ ) and, with the nominal assumption of $H_{98,y2} = 1$ , SPARC is projected to attain $Q \approx 11$ and $P_{\textrm {fusion}} \approx 140$ MW. SPARC will therefore constitute a unique platform for burning plasma physics research with high density ( $\langle n_{e} \rangle \approx 3 \times 10^{20}\ \textrm {m}^{-3}$ ), high temperature ( $\langle T_e \rangle \approx 7$ keV) and high power density ( $P_{\textrm {fusion}}/V_{\textrm {plasma}} \approx 7\ \textrm {MW}\,\textrm {m}^{-3}$ ) relevant to fusion power plants. SPARC's place in the path to commercial fusion energy, its parameters and the current status of SPARC design work are presented. This work also describes the basis for global performance projections and summarizes some of the physics analysis that is presented in greater detail in the companion articles of this collection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Accurately simulating core-collapse self-interacting dark matter halos

The properties of satellite halos provide a promising probe for dark matter (DM) physics. Observations have motivated current efforts to explain surprisingly compact DM halos. If DM is not collisionless, but has strong self-interactions, halos can undergo gravothermal collapse, leading to higher densities in the central region of the halo. However, it is challenging to model this collapse phase from first principles. To improve on this, we sought to better understand the numerical challenges and convergence properties of self-interacting dark matter (SIDM) N-body simulations in the collapse phase. Especially, our aim was to better understand the evolution of satellite halos. To do so, we ran SIDM N-body simulations of a low-mass halo in isolation and within an external gravitational potential. The simulation set-up was motivated by the perturber of the stellar stream GD-1. We find that the halo evolution is very sensitive to energy conservation errors, and a SIDM kernel size that is too large can artificially speed up the collapse. Moreover, we demonstrate that the King model can describe the density profile at small radii for the late stages that we have simulated. Furthermore, for our most highly resolved simulation (N = 5 × 10 7 ) we have made the data public. It can serve as a benchmark. Overall, we find that the current numerical methods do not suffer from convergence problems in the late collapse phase and provide guidance on how to choose numerical parameters, for example that the energy conservation error is better kept well below 1%. This allows simulations to be run of halos that become concentrated enough to explain observations of GD-1-like stellar streams or strong gravitational lensing systems.

dark matter↗

Whole-genome demography of COVID-19 virus during its pandemic period and on “panvalent” vaccine design

With over 16 million submitted genomic sequences, the SARS-CoV-2 (SC2) virus, the cause of the most recent worldwide COVID-19 pandemic, has become the most sequenced genome of all known viruses, revealing, for example, a vast number of expanding viral lineages. Since the pandemic phase appears to be over, we performed a retrospective re-examination of the demographic grouping pattern and their genomic characteristics during the entire pandemic period up to the peak of the last pandemic wave. For our study, we extracted from the NCBI only unique viral sequences and converted each sequence data to a relational vector, indicating the presence/absence of each variational event compared to a “reference” sequence. Our study revealed several genomic features that are unexpected or different from those of previous studies. For example, approximately 44,000 variants with unique sequences emerged during the pandemic period; they group into only four major viral-genomic groups and each has a set of mostly unique highly-conserved variant-genotypes (HCVGs); and a small set from the first (“ancestral”) group was inherited by the three (“descendant”) groups, suggesting that HCVGs in the next group may be predictable from the current group(s). Such a concept may be potentially important in designing “panvalent” vaccines against the current and future waves of viral infections.

60 APPLIED LIFE SCIENCES↗

Xyce(™) Parallel Electronic Simulator v.7.5

The Xyce Parallel Electronic Simulator simulates electronic circuit behavior in DC, AC, HB, MPDE and transient mode using standard analog (DAE) and/or device (PDE) device models including several age and radiation aware devices. It supports a variety of computing platforms (both serial and parallel) computers. Lastly, it uses a variety of modern solution algorithms dynamic parallel load-balancing and iterative solvers.! ! Xyce is primarily used to simulate the voltage and current behavior of a circuit network (a network of electronic devices connected via a conductive network). As a tool, it is mainly used for the design and analysis of electronic circuits.! ! Kirchoff's conservation laws are enforced over a network using modified nodal analysis. This results in a set of differential algebraic equations (DAEs). The resulting nonlinear problem is solved iteratively using a fully coupled Newton method, which in turn results in a linear system that is solved by either a standard sparse-direct solver or iteratively using Trilinos linear solver packages, also developed at Sandia National Laboratories.

Source record↗

Dialogue Between Lighting and HVAC Systems: Improving building system integration

Lighting systems have long been capable of sensing when someone enters or exits a room and using that knowledge to turn lights on or off. More recently, connected lighting systems with sensors integrated into every luminaire have become broadly available, facilitating highly granular occupancy detection. Similarly, HVAC systems have long been able to use an understanding of building occupancy to adjust temperature setpoints and reduce energy use without significant impacts to occupant comfort. Energy codes (e.g., ASHRAE/IES Standard 90.1, IECC, Title 24) now require “occupied standby HVAC control,” whereby systems adjust both temperature and ventilation setpoints in zones that are determined to be unoccupied during normal occupancy hours. Here, this article discusses current issues that stymie the integration of Lighting and HVAC systems, and DOE activities focused on addressing them.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Reactor Pressure Vessel Fracture Mechanics Development and Concrete Application Testing for Grizzly

The Grizzly code is being developed to address degradation issues in nuclear reactor structures and components. For light-water reactors, Grizzly currently has capabilities to simulate degradation processes and their effects on structural integrity in two key areas: reactor pressure vessels (RPVs) and reinforced concrete structures. This report documents improvements made to Grizzly’s ability to address both of these structural systems. For RPVs, the reduced-order models (ROMs) used in fracture mechanics calculations have been expanded to allow their application over a broader range of the parameter space than was permitted by the previous models. The ROMs currently used in Grizzly for the evaluation of flaws that are fully embedded within the RPV (as opposed to surface-breaking flaws) are based on a model that is known to be conservative, indicating higher stress intensity factors than would be obtained from direct simulations. A more accurate model that eliminates these excess conservatisms has been recently included in the American Society of Mechanical Engineers Boiler and Pressure Vessel Code but was not applicable for flaws near the RPV surface, which is where the most critical flaws are usually located. That model has recently been extended for increased applicability in this near-surface region. The ROMs for embedded flaws in the Grizzly code have been expanded to include these recent extensions, which permit their use in a much broader set of cases than previously possible. Direct 3D simulations have been used to check these ROMs and have shown good agreement in most cases, although there are still some cases that need further investigation. There are considerable benefits to using these these more accurate and less conservative ROMs for embedded flaws. On a benchmark probabilistic fracture mechanics problem tested here, the conditional probability of fracture initiation computed for a population of flaws in a single plate in an RPV decreased by over a factor of 3. To address aging in reinforced concrete structures, a capability to simulate multiple degradation mechanisms, including alkali-silica reaction and radiation-induced volumetric expansion has been developed in Grizzly over the past several years. This had previously been demonstrated on laboratory-scale specimens but not on full-scale nuclear concrete structures with reinforcement. To demonstrate the applicability of Grizzly to the analysis of large-scale structures of interest, a full 3D model of a representative reinforced concrete structure, including a complex arrangement of reinforcing bars, was developed and demonstrated in Grizzly.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Computer Science Research Needs for Parallel Discrete Event Simulation (PDES)

Historically, scientific computing efforts have demonstrated the clear need for, and effective use of, supercomputing with traditional time-stepped simulations. Nevertheless, there are several areas in the mission spaces of the U.S. Department of Energy and other agencies waiting to tap advanced computing research using a different, discrete event style of modeling, simulation, and analysis. These span a wide spectrum of applications including energy grid resilience, urban planning and policy, transportation science, building technologies, emergency response and planning, environmental impact analysis, computational epidemiology, Internet communications, cyber security, and cyber-physical systems, to name only a few. Even within traditional scientific applications, the role of discrete event modes of execution is increasing in the form of new event-based mathematical solvers such as quantized state integration methods and discrete-continuous hybrid system solvers. Co-design of advanced supercomputing hardware systems is another area that exploits discrete event simulation at its core for effective analyses. Complex systems, entity behaviors and interconnections play a significant role in all these applications, which are mapped to large-scale models with discrete event formulations. To make advancements in all the aforementioned scientific areas, many technical aspects need to be more thoroughly studied and deeply understood in parallel discrete event simulation (PDES). The unique dynamics inherent in a discrete event modeling approach, by their very nature, intersect and influence the entire stack of the computing system, including (a) the unique nature of the instruction sets exercised in PDES workloads without a predominance of high-precision floating point operations, (b) virtual time-constrained multi-threaded execution of many logical processes per processor, (c) extremely variable and difficult to predict network traffic characteristics, (d) interfaces and inter-dependencies with machine learning and artificial intelligence codes at higher software layers, and (e) highly challenging load balancing needs, especially in effectively accounting for accelerated/extremely heterogeneous computing in current and future high-performance computing systems. Efficient and accurate parallel execution of PDES workloads is also dominated by challenges in dealing with their asynchronous concurrency fundamentally present at the model level. Conservative synchronization, optimistic/speculative synchronization, and their hybrid schemes open new questions in fundamental computer science with respect to reversibility of computation and prediction (lookahead) of behaviors inherent within model codes. On the implementation front, there are relatively few scalable, general-purpose parallel discrete event simulators in the world, and even fewer have been studied on emerging hardware platforms. To enable scientific advances using PDES, the research needs in computer science must also be pursued and met in the intersection of the algorithmic and hardware-aware aspects of scalable PDES engines. This report is aimed at capturing a computer science-oriented view of this important area of research in PDES, presenting a sample of important applications with their inherent discrete event technology elements. Needs are outlined in core areas of parallel discrete event research as well as cross-cutting directions in computer science research that positively impact scientific advancements across several important application areas. A selection of priority research opportunities in advanced computing for PDES is identified to serve as reference for key research topics and their order of importance for scientific advancements.

97 MATHEMATICS AND COMPUTING↗

Effects of Borate and Organics on U(VI) Solubility in WIPP Brine

The solubility of uranium (VI) in Waste Isolation Pilot Plant (WIPP)-relevant brine was determined to support ongoing WIPP recertification activities. This research was performed by the Los Alamos National Laboratory Carlsbad Operations (LANL-CO) Actinide Chemistry and Repository Science Program (ACRSP). The WIPP Actinide Source Term Program (ASTP) did not develop a model for the solubility of actinides in the VI oxidation state. The solubility of UO 2 2+ , in the absence of WIPP specific data, is presently set to be equal to a conservatively high 1 mM within the WIPP Performance Assessment (PA) for all expected WIPP conditions (SOTERM, 2019) as selected at the recommendation of the Environment Protection Agency (EPA) (EPA, 2005). According to the current WIPP chemistry model assumptions and conditions, the expected pC H+ is about 9.5 and controlled by MgO buffering CO 3 2- . The goal of this study is to perform screening experiments that account for the contributions of organics and borate on uranium solubility. In this report, the solubility of U(VI) was determined at pC H+ 9 WIPP brine in the absence or presence of borate and organics at under-saturation approach. Experiments were equilibrated for about 135 days. Organic compounds present in WIPP waste can form strong complexes with actinides and can affect the oxidation states of actinides. The organic compounds addressed in WIPP performance assessment include EDTA (Ethylenediaminetetraacetic Acid), oxalate, citrate, and acetate (SOTERM, 2019). These data quantify the effects of WIPP-relevant concentrations of borate and organics effects on the solubility of U(VI) to challenge the predictions of the WIPP actinide model and inform decisions and recommendations made in the upcoming recertification of the WIPP (CRA-2024). The experiments performed were done according to the U.S. Department of Energy (DOE) approved Test Plan entitled “Experimental Strategy to Challenge Actinide Solubility Predictions” and designated LCO-ACP-26. All data reported were obtained under the LANL-CO Quality Assurance Program, which is compliant with the DOE Carlsbad Field Office, Quality Assurance Program Document (CBFO/QAPD) (QAPD, 2017).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗