SEARCH · Engineering Papers

Results for “numerical relativity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Risk-Averse Bi-Level Stochastic Network Interdiction Model for Cyber-Security Risk Management

Security of cyber networks is crucial; recent severe cyber-attacks have had a devastating effect on many large organizations. The attack graph, which maps the potential attack paths of a cyber network, is a popular tool for analyzing cyber system vulnerability. In this study, we propose a bi-level stochastic network interdiction model on an attack graph to enable a risk-averse, resource constrained cyber network defender to optimally deploy security countermeasures that protect against attackers with an uncertain budget. This risk- averse conditional-value-at-risk (CVaR) model minimizes a weighted sum of the expected maximum loss over all scenarios and the expected maximum loss from the most damaging attack scenarios. We develop a customized constraint and column generation algorithm to solve our model as well as several acceleration techniques to improve the computational efficiency. Numerical experiments demonstrate that the acceleration techniques enable the solution of relatively large problems within a reasonable amount of time: applying all the acceleration techniques also reduces the average computation time of the basic algorithm by 71% for 100-node graphs. Using metrics called mean-risk value of stochastic solution and value of risk-aversion, computational results suggest that our stochastic risk-averse model significantly outperforms deterministic and risk-neutral models when 1) the distribution of attacker budget is heavy-right-tailed and 2) the defender is highly risk-averse.

Bhuiyan, Tanveer H.↗

DOE OSTI2021

A screening model for predicting injection well pressure buildup and plume extent in CO 2 geologic storage projects

We present a new screening model for predicting injection well pressure buildup and CO 2 plume migration for CO 2 geologic sequestration projects. The model requires only limited information and is quite accurate when compared to detailed simulation results. Such models can assist project developers during the early days of project planning (e.g., for 45Q related projects), and also help regulators perform simple checks against detailed numerical models. The screening model consists of two correlations: one for injectivity index in terms of the slope of the CO 2 fractional flow curve, and the second for total storage efficiency within the footprint of plume as a function of gravity number and slope of the CO 2 fractional flow curve. Using these two correlations, as well as a knowledge of some basic reservoir characteristics and estimates of fluid properties from standard correlations, the injection-well pressure buildup and CO 2 plume extent in the formation can be readily estimated. Finally, the new correlations show a good match with the results of the underlying simulations, and also provide good agreement with independent calculations for two example problems.

42 ENGINEERING↗

DOE OSTI2023

Assessing the source of error in the Thomas–Fermi–von Weizsäcker density functional

We investigate the source of error in the Thomas–Fermi–von Weizsäcker (TFW) density functional relative to Kohn–Sham density functional theory (DFT). In particular, through numerical studies on a range of materials, for a variety of crystal structures subject to strain and atomic displacements, we find that while the ground state electron density in TFW orbital-free DFT is close to the Kohn–Sham density, the corresponding energy deviates significantly from the Kohn–Sham value. Here, we show that these differences are a consequence of the poor representation of the linear response within the TFW approximation for the electronic kinetic energy, confirming conjectures in the literature. In so doing, we find that the energy computed from a non-self-consistent Kohn–Sham calculation using the TFW electronic ground state density is in very good agreement with that obtained from the fully self-consistent Kohn–Sham solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Colliding Pulse Mode Locking and Nonequilibrium Carrier Dynamics within Vertical External-Cavity Surface-Emitting Lasers

The complex nonlinear dynamics of colliding counterpropagating ultrashort laser pulses circulating in a vertical external-cavity surface-emitting ring laser cavity are captured by a first-principles quantum-mechanical microscopic model coupled to an intracavity beam propagator. Furthermore, our numerical simulations establish that overall mode locking stability is sensitive to relative phases and delays of the pulses impinging on the semiconductor saturable absorber mirror. Delay offsets of the counterpropagating pulses lead to asymmetric carrier recovery and instability, often leading in turn to single unidirectional circulating pulses. Although computationally challenging, we also include transverse effects and show that the influence of interference patterns from spatial gratings written by the counterpropagating pulses is significant in stabilizing the modelocked pulses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Soft pions and transport near the chiral critical point

During the expansion of a heavy ion collision, the system passes close to the $\textit{O}$(4) critical point of QCD, and thus the fluctuations of the order parameter $(σ, \overrightarrow{π}$) are expected to be enhanced. In this work, our goal is to compute how these enhanced fluctuations modify the transport coefficients of QCD near the pseudocritical point. We also make a phenomenological estimate for how chiral fluctuations could effect the momentum spectrum of soft pions. We first formulate the appropriate stochastic hydrodynamic equations close to the $\textit{O}$(4) critical point. Then, working in mean field, we determine the correlation functions of the stress tensor and the currents which result from this stochastic real-time theory, and use these correlation functions to determine the scaling behavior of the transport coefficients. The hydrodynamic theory also describes the propagation of pion waves, fixing the scaling behavior of the dispersion curve of soft pions. We present scaling functions for the shear viscosity and the charge conductivities near the pseudocritical point and estimate the absolute magnitude of the critical fluctuations to these parameters and the bulk viscosity. Using the calculated pion dispersion curve, we estimate the expected critical enhancement of soft pion yields, and this estimate provides a plausible explanation for the excess seen in experiment relative to ordinary hydrodynamic computations. Our results motivate further phenomenological and numerical work on the implications of chiral symmetry on real-time properties of thermal QCD near the pseudocritical point.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

From Failure to Insight: Analyzing Disk Breakdowns in Large-Scale HPC Environments

Disk failure data provides valuable insights for preventing failures, enhancing storage robustness, guiding system design and deployment, and ensuring reliable operations at data centers. This paper introduces two disk failure datasets collected from large-scale HPC production environments over the past five years, comprising over 5,000 failure records from more than 40,000 disks. We analyzed these datasets across multiple dimensions, including temporal, spatial, and relational trends, and performed a comprehensive reliability assessment. Our analysis yielded numerous observations and insights that influence various operational aspects of HPC storage systems. We believe this study offers a holistic understanding of disk failure trends likely to interest the HPC storage community.

George, Anjus↗

DOE OSTI2020

Revealing Three-Dimensional Morphology in Nanoporous Gold Using Three-Dimensional X-Ray Fresnel Coherent Diffractive Imaging Tomography

Nanoporous metals fabricated by dealloying have a unique bi-continuous, sponge-like porous structure with ultra-high surface area. The unique properties of these materials, especially nanoporous gold, have numerous potential applications in sensors and actuators and in energy-related applications such as catalytic materials, super-capacitors, and battery supports. The degree of porosity and size of the metal ligaments are critical parameters that determine many properties and thus govern the functionalities of nanoporous metals in many applications including energy storage and conversion. Additionally, we used Fresnel coherent diffractive imaging combined with tomographic reconstruction to quantify the nanoscale three-dimensional spatial distribution and homogeneity of the porosity and ligament size within a bulk sample of nanoporous gold. The average porosity and its standard deviation along the axial direction through the sample were determined, as well as the characteristic feature size and its standard deviation. The result shows that free corrosion is an effective way to create homogeneous nanoporous metals with sample sizes on the order of 1 µm.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Temporal reflection from short pump pulses inside a dispersive nonlinear medium: the impact of pump parameters

We have studied, through a series of experiments and numerical simulations, how temporal reflection from an intense pump pulse inside a photonic crystal fiber is affected by parameters of the pump pulse used to form a moving high-index boundary. We used femtosecond pump pulses, which slow down inside the fiber as their spectrum red-shifts because of intrapulse Raman scattering. Temporal reflection of probe pulses occurs from such decelerating pump pulses. We changed the width and chirp of our pump pulses with a 4f pulse shaper capable of providing both spectral filtering and frequency chirping. We found that temporal refection exhibited novel features, to our knowledge, when pump pulses were made wider or chirped. In both cases, two or more reflected pulses were produced at different wavelengths in a specific range of the initial pump-probe delays. Furthermore, numerical simulations reveal that the origin of such novel features is related to the complex nonlinear evolution of pump pulses inside optical fibers.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

Studies for the selection of fully contained muon neutrino CC events in the ICARUS T600 detector at Fermilab

The ICARUS T600 LAr-TPC detector has started its new physics runs in June 2022 at Fermilab within the SBN program. This detector is recording neutrino interactions from the Booster BNB neutrino beam in order to definitively clarify the open questions related to the possible existence of sterile neutrinos, as suggested by numerous observed experimental anomalies. In addition the neutrinos events recorded from the NuMI off-axis beam are under study in order to perform neutrino-Argon cross section measurements. At Fermilab, ICARUS is facing a challenging experimental condition: the detector is presently installed essentially at Earth surface where cosmic ray particles can become a serious source of background for the neutrino event search. This condition makes it necessary to deploy suitable automatic tools for the identification, selection, and measurement of the neutrino events among the millions of events triggered by cosmics. In this contribution, two automatic selection procedures devoted to the identification of fully contained muon neutrino interactions in the BNB neutrino beam will be presented, together with the first results of their application on simulated and on recently recorded events in the T600 detector.

Farnese, Christian [INFN, Padua] (ORCID:0000000163↗

DOE OSTI2025

Studies for the Selection of Fully Contained Muon Neutrino CC Events in the ICARUS T600 Detector at Fermilab

Farnese, Christian [INFN, Padua]↗

DOE OSTI2024

Research Reactors Division Infrastructure Investment Plan for the High Flux Isotope Reactor

The High Flux Isotope Reactor (HFIR) is a unique national asset. Operational for nearly 60 years, continued investment into the aging infrastructure is necessary to ensure operation for another 6 decades. Additionally, growing missions require HFIR as well as important upgrades. Consequently, carefully integrated planning is required to ensure that infrastructure investments are timely executed to ensure long-term, reliable operation of HFIR. Concerns about challenges to the operational reliability of HFIR resulted in a recommendation from the 2023 Operations Review by the US Department of Energy (DOE) Office of Basic Energy Sciences that a HFIR management strategy be developed to address the infrastructure needs. This report defines the investment needs, which are evolving as new upgrade efforts are better defined. HFIR is part of the three-source strategy within the Neutron Sciences Directorate (NScD) and contributes to the five strategic science areas outlined in the NScD 10 Year Strategic Science Plan: quantum materials, soft matter, materials and engineering, chemistry, and biosciences. Fundamental to this strategy are three core values: operational excellence, responsible stewardship, and servant leadership. These values guide our mission of safe and reliable operation of the reactor and require a strong and just nuclear safety culture, a solemn respect for responsible care of the facility, good workforce development, robust procedures and processes, an effective communication strategy, world-class asset management, a determined customer focus, and a commitment to protecting the environment, the safety and health of the public and our people, and the quality of work performed within our facility. These principles are all essential to operate HFIR at a world-class level. The Research Reactors Division (RRD) will lead a new era of neutron science and isotope production at HFIR through responsible and purposeful leadership and unwavering support of the science community. The approach outlined in this plan highlights the direction leadership is taking to ensure that HFIR is ready to support the science challenges and national needs of the future and that the United States maintains world leadership in neutron sciences. The plan is in alignment with the DOE’s desire to continue operating HFIR and with the NScD strategic science goals for the future. HFIR is an aging facility with numerous infrastructure challenges and needs. It has an aging workforce in relation to the general population of Oak Ridge National Laboratory (ORNL), with many expected retirements over the next 5–10 years. With an increase in work scope caused by changing national priorities and science goals, several critical hires have been identified. To manage HFIR’s infrastructure needs, a prioritized list of equipment upgrades has been identified along with an analysis of future staffing requirements. A desire to operate HFIR at eight cycles per year will necessarily require some significant changes to procedures and processes currently in place as well as targeted staffing additions. Many of the equipment upgrades identified in this plan will significantly increase the reliability of the plant, thus contributing to the effort to reach the goal of safely operating eight cycles per year. A plan to attain eight-cycle operation is being prepared in parallel with the activities identified in this plan, although the actions identified to satisfy both plans will overlap. This plan identifies new infrastructure needs—for both plant equipment and staffing—thus necessitating formulation of future budget requests to fund the increased work scope and improvement activities. Some activities are currently being scheduled with the expectation that funding will be received. Any delays to funding or reductions of funding from the identified cost estimations will directly and negatively affect the plan’s implementation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2021

A Lamination Model for Pressure-Assisted Sintering of Multilayered Porous Structures

This work describes a lamination model for pressure-assisted sintering of thin, multilayered, and porous structures based on the linear viscous constitutive theory of sintering and the classical laminated plate theory of continuum mechanics. A constant out-of-plane normal stress is assumed in the constitutive relation. The lamination relations between the force/moment resultants and the strain/curvature rates are presented. Numerical simulations were performed for a symmetric tri-layer laminate consisting of a 10% gadolinia doped ceria (Ce 0.9 Gd 0.1 O 1.95-δ ) composite structure, where porous layers were adhered to the top and bottom of a denser layer under uniaxially-applied pressures and the sinter forging conditions. The numerical results show that, compared with free sintering, the applied pressure can significantly reduce the sintering time required to achieve given layer thicknesses and porosities. Unlike free sintering, which results in a monotonic decrease of the laminate in-plane dimension, pressure-assisted sintering may produce an in-plane dimension increase or decrease, depending on the applied pressure and sintering time. Finally, the individual layers in the laminate exhibit different stress characteristics under pressure-assisted sintering.

42 ENGINEERING↗

DOE OSTI2022

Persistent mysteries of jet engines, formation, propagation, and particle acceleration: Have they been addressed experimentally?

Here, the physics of astrophysical jets can be divided into three regimes: (i) engine and launch (ii) propagation and collimation, (iii) dissipation and particle acceleration. Since astrophysical jets comprise a huge range of scales and phenomena, practicality dictates that most studies of jets intentionally or inadvertently focus on one of these regimes, and even therein, one body of work may be simply boundary condition for another. We first discuss long standing persistent mysteries that pertain the physics of each of these regimes, independent of the method used to study them. This discussion makes contact with frontiers of plasma astrophysics more generally. While observations theory, and simulations, and have long been the main tools of the trade, what about laboratory experiments? Jet related experiments have offered controlled studies of specific principles, physical processes, and benchmarks for numerical and theoretical calculations. We discuss what has been accomplished on these fronts. Although experiments have indeed helped us to understand certain processes, proof of principle concepts, and benchmarked codes, they have yet to solved an astrophysical jet mystery on their own. A challenge is that experimental tools used for jet-related experiments so far, are typically not machines originally designed for that purpose, or designed with specific astrophysical mysteries in mind. This presents an opportunity for a different way of thinking about the development of future platforms: start with the astrophysical mystery and build an experiment to address it.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

General approaches for shear-correcting coordinate transformations in Bragg coherent diffraction imaging. Part II

X-ray Bragg coherent diffraction imaging (BCDI) has been demonstrated as a powerful 3D microscopy approach for the investigation of sub-micrometre-scale crystalline particles. The approach is based on the measurement of a series of coherent Bragg diffraction intensity patterns that are numerically inverted to retrieve an image of the spatial distribution of the relative phase and amplitude of the Bragg structure factor of the diffracting sample. This 3D information, which is collected through an angular rotation of the sample, is necessarily obtained in a non-orthogonal frame in Fourier space that must be eventually reconciled. To deal with this, the approach currently favored by practitioners (detailed in Part I) is to perform the entire inversion in conjugate non-orthogonal real- and Fourier-space frames, and to transform the 3D sample image into an orthogonal frame as a post-processing step for result analysis. In this article, which is a direct follow-up of Part I, two different transformation strategies are demonstrated, which enable the entire inversion procedure of the measured data set to be performed in an orthogonal frame. The new approaches described here build mathematical and numerical frameworks that apply to the cases of evenly and non-evenly sampled data along the direction of sample rotation ( i.e. the rocking curve). The value of these methods is that they rely on the experimental geometry, and they incorporate significantly more information about that geometry into the design of the phase-retrieval Fourier transformation than the strategy presented in Part I. Finally, two important outcomes are (1) that the resulting sample image is correctly interpreted in a shear-free frame and (2) physically realistic constraints of BCDI phase retrieval that are difficult to implement with current methods are easily incorporated. Computing scripts are also given to aid readers in the implementation of the proposed formalisms.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Water transport properties of virtual fractal porous media: Implications for the unsaturated transport properties of cement-based materials

Cement-based materials are more and more recognized as fractal materials. From a practical point of view, fractality means that the pore size distribution can be described using fractal scaling law (i.e. power functions). Here, the simplest scaling law making use of a single fractal dimension was used to generate simple and virtual fractal porous media (based on bundles of parallel cylindrical pores). The capillary curve and permeability (relative and intrinsic) of the virtual porous media were then estimated through theoretical and numerical approaches. The results show that: (1) although bereft of any physical basis, van Genuchten equation bears some fractal information; (2) although highly flexible, van Genuchten equation fails to perfectly fit the capillary curves and can thus generate error in inverse analysis (for intrinsic permeability evaluation) and (3) there seems to be a relation between the intrinsic permeability and van Genuchten pressure parameter.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Integral X-ray diffuse scattering for studying irradiation-induced dislocation loops in single crystals

Numerical diffuse scattering cross-section calculations are used to establish a rigorous basis for determining the concentration and size distribution of dislocation loops in irradiated single crystals from integral X-ray diffuse scattering (XRDS) measurements. Differential XRDS intensities for prismatic {111} type dislocation loops are numerically calculated as a function of loop radius R and wavevectors q relative to Bragg reflections in tungsten. The results show the well known 1/q 2 Huang scattering form at small q that transitions to a ∼1/q 4 dependence associated with the Stokes–Wilson approximation for q ≳ 1/R. More importantly, they show further that the 1/q 4 falloff is not the asymptotic large-q form of the diffuse scattering for small loops (R < 200 Å) as has often been assumed. Rather, for loop sizes as small as R ≃ 5 Å with strong curvature, the calculations show definitively that the scattering transitions to a robust 1/q 5 falloff at larger q that arises due to the local strains near the dislocation core defining the circumference of the dislocation loops. The presence of this 1/q 5 asymptotic form for both small and large loops is experimentally confirmed using an integral XRDS measurement around the 110 reflection on self-ion-irradiated tungsten combined with numerically calculated integral XRDS cross-sections. Accordingly, the historical two-region theoretical treatment of the cross-sections for integral XRDS is extended to a three-region model that has direct sensitivity to the (first-moment) dislocation line lengths of dislocation loops. In conclusion, these developments enable the use of both numerical and analytically modeled cross-sections to make accurate integral XRDS determinations of dislocation loop sizes and concentrations using modest-intensity laboratory X-ray sources.

42 ENGINEERING↗

DOE OSTI2024

Weak baselines and reporting biases lead to overoptimism in machine learning for fluid-related partial differential equations

One of the most promising applications of machine learning in computational physics is to accelerate the solution of partial differential equations (PDEs). The key objective of machine-learning-based PDE solvers is to output a sufficiently accurate solution faster than standard numerical methods, which are used as a baseline comparison. Here, we first perform a systematic review of the ML-for-PDE-solving literature. Out of all of the articles that report using ML to solve a fluid-related PDE and claim to outperform a standard numerical method, we determine that 79% (60/76) make a comparison with a weak baseline. Second, we find evidence that reporting biases are widespread, especially outcome reporting and publication biases. We conclude that ML-for-PDE-solving research is overoptimistic: weak baselines lead to overly positive results, while reporting biases lead to under-reporting of negative results. To a large extent, these issues seem to be caused by factors similar to those of past reproducibility crises: researcher degrees of freedom and a bias towards positive results. We call for bottom-up cultural changes to minimize biased reporting as well as top-down structural reforms to reduce perverse incentives for doing so.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Differentiable programming for online training of a neural artificial viscosity function within a staggered grid Lagrangian hydrodynamics scheme

Lagrangian methods to solve the inviscid Euler equations produce numerical oscillations near shock waves. A common approach to reducing these oscillations is to add artificial viscosity (AV) to the discrete equations. The AV term acts as a dissipative mechanism that attenuates oscillations by smearing the shock across a finite number of computational cells. However, AV introduces several control parameters that are not determined by the underlying physical model, and hence, in practice are tuned to the characteristics of a given problem. We seek to improve the standard quadratic-linear AV form by replacing it with a learned neural function that reduces oscillations relative to exact solutions of the Euler equations, resulting in a hybrid numerical-neural hydrodynamic solver. Because AV is an artificial construct that exists solely to improve the numerical properties of a hydrodynamic code, there is no offline ‘viscosity data’ against which a neural network can be trained before inserting into a numerical simulation, thus requiring online training. We achieve this via differentiable programming, i.e. end-to-end backpropagation or adjoint solution through both the neural and differential equation code, using automatic differentiation of the hybrid code in the Julia programming language to calculate the necessary loss function gradients. A novel offline pre-training step accelerates training by initializing the neural network to the default numerical AV scheme, which can be learned rapidly by space-filling sampling over the AV input space. We find that online training over early time steps of simulation is sufficient to learn a neural AV function that reduces numerical oscillations in long-term hydrodynamic shock simulations. These results offer an early proof-of-principle that online differentiable training of hybrid numerical schemes with novel neural network components can improve certain performance aspects existing in purely numerical schemes.

97 MATHEMATICS AND COMPUTING↗