Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 1,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Preliminary Structural Design Considerations and Mass Efficiencies for Lunar Surface Manipulator Concepts

The mass and sizing characteristics of manipulators for Lunar and Mars planetary surface applications are investigated by analyzing three structural configurations: a simple cantilevered boom with a square tubular cross-section; a hybrid cable/boom configuration with a square tubular cross-section support structure; and a hybrid cable/boom configuration with a square truss cross-section support structure. Design procedures are developed for the three configurations and numerical examples are given. A new set of performance parameters are developed that relate the mass of manipulators and cranes to a loading parameter. These parameters enable the masses of different manipulator configurations to be compared over a wide range of design loads and reach envelopes (radii). The use of these parameters is demonstrated in the form of a structural efficiency chart using the newly considered manipulator configurations. To understand the performance of Lunar and Mars manipulators, the design procedures were exercised on the three manipulator configurations assuming graphite/epoxy materials for the tubes and trusses. It is also assumed that the actuators are electric motor, gear reduction systems. Numerical results for manipulator masses and sizes are presented for a variety of manipulator reach and payload mass capabilities. Results are presented that demonstrate the sensitivity of manipulator mass to operational radius, tip force, and actuator efficiency. The effect of the value of gravitational force on the ratio of manipulator-mass to payload-mass is also shown. Finally, results are presented to demonstrate the relative mass reduction for the use of graphite/epoxy compared to aluminum for the support structure.

Dorsey, John T.↗
DOE OSTI2018

A note on the S-matrix bootstrap for the 2d O(N) bosonic model

In this work we apply the S-matrix bootstrap maximization program to the 2d bosonic O ( N ) integrable model which has N species of scalar particles of mass m and no bound states. Since in previous studies theories were defined by maximizing the coupling between particles and their bound states, the main problem appears to be to find what other functional can be used to define this model. Instead, we argue that the defining property of this integrable model is that it resides at a vertex of the convex space determined by the unitarity and crossing constraints. Thus, the integrable model can be found by maximizing any linear functional whose gradient points in the general direction of the vertex, namely within a cone determined by the normals to the faces intersecting at the vertex. This is a standard problem in applied mathematics, related to semi-definite programming and solvable by fast available numerical algorithms. The information provided by the numerical solution is enough to reproduce the known analytical solution without using integrability, namely the Yang-Baxter equation. This situation seems quite generic so we expect that other theories without continuous parameters can also be found by maximizing linear functionals in the convex space of allowed S-matrices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

A semi-discrete scheme for solving fourth-order partial integro-differential equation with a weakly singular kernel using Legendre wavelets method

In this paper, a semi-discrete scheme is presented for solving fourth-order partial integro-differential equations with a weakly singular kernel. The second-order backward difference formula is used to discretize the temporal derivatives. After discretizing the temporal derivatives, the considered problems are converted into a set of ordinary differential equations which are solved by using Legendre wavelets collocation method. The stability and convergence properties related to the time discretization are discussed and theoretically proven. Several numerical examples are included to demonstrate the accuracy and efficiency of the proposed method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Phase controlled metal–insulator transition in multi-leg quasiperiodic optical lattices

A tight-binding model of a multi-leg ladder network with a continuous quasiperiodic modulation in both the site potential and the inter-arm hopping integral is considered. The model mimics optical lattices where ultra-cold fermionic or bosonic atoms are trapped in double well potentials. It is observed that, the relative phase difference between the on-site potential and the inter-arm hopping integral, which can be controlled by the tuning of the interfering laser beams trapping the cold atoms, can result in a mixed spectrum of one or more absolutely continuous subband(s) and point like spectral measures. This opens up the possibility of a re-entrant metal–insulator transition. The subtle role played by the relative phase difference mentioned above is revealed, and we corroborate it numerically by working out the multi-channel electronic transmission for finite two-, and three-leg ladder networks. The extension of the calculation beyond the two-leg case is trivial, and is discussed in the work. - Graphical abstract: ▪ - Highlights: • Phase controlled metal–insulator transition is discussed. • An analytical prescription is given to understand MI transition. • Our work provides a way of designing experiments involving laser beams.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2018

Direct numerical simulation of reactor two-phase flows enabled by high-performance computing

Nuclear reactor two-phase flows remain a great engineering challenge, where the high-resolution two-phase flow database which can inform practical model development is still sparse due to the extreme reactor operation conditions and measurement difficulties. Owing to the rapid growth of computing power, the direct numerical simulation (DNS) is enjoying a renewed interest in investigating the related flow problems. A combination between DNS and an interface tracking method can provide a unique opportunity to study two-phase flows based on first principles calculations. More importantly, state-of-the-art high-performance computing (HPC) facilities are helping unlock this great potential. This paper reviews the recent research progress of two-phase flow DNS related to reactor applications. The progress in large-scale bubbly flow DNS has been focused not only on the sheer size of those simulations in terms of resolved Reynolds number, but also on the associated advanced modeling and analysis techniques. Specifically, the current areas of active research include modeling of sub-cooled boiling, bubble coalescence, as well as the advanced post-processing toolkit for bubbly flow simulations in reactor geometries. A novel bubble tracking method has been developed to track the evolution of bubbles in two-phase bubbly flow. Also, spectral analysis of DNS database in different geometries has been performed to investigate the modulation of the energy spectrum slope due to bubble-induced turbulence. In addition, the single-and two-phase analysis results are presented for turbulent flows within the pressurized water reactor (PWR) core geometries. The related simulations are possible to carry out only with the world leading HPC platforms. Furthermore, these simulations are allowing more complex turbulence model development and validation for use in 3D multiphase computational fluid dynamics (M-CFD) codes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Validation of plasma shape reconstruction by Cauchy condition surface method in KSTAR

Cauchy Condition Surface (CCS) method is a numerical approach to reconstruct the plasma boundary and calculate the quantities related to plasma shape using the magnetic diagnostics in real time. It has been applied to the KSTAR plasma in order to establish the plasma shape reconstruction with the high elongation of plasma shape and the large effect of eddy currents flowing in the tokamak structures for the first time. For applying the CCS calculation to the KSTAR plasma, the effects by the eddy currents and the ferromagnetic materials on the plasma shape reconstruction are studied. The CCS calculation includes the effect of eddy currents and excludes the magnetic diagnostics, which is expected to be influenced largely by ferromagnetic materials. Calculations have been performed to validate the plasma shape reconstruction in 2012 KSTAR experimental campaign. Comparison between the CCS calculation and non-magnetic measurements revealed that the CCS calculation can reconstruct the accurate plasma shape even with a small I{sub P}.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Novel signal inversion of laser beam induced current for femtosecond-laser-drilling-induced junction on vacancy-doped p-type HgCdTe

In this paper, experimental results of temperature-dependent signal inversion of laser beam induced current (LBIC) for femtosecond-laser-drilling-induced junction on vacancy-doped p-type HgCdTe are reported. LBIC characterization shows that the traps induced by femtosecond laser drilling are sensitive to temperature. Theoretical models for trap-related p-n junction transformation are proposed and demonstrated using numerical simulations. The simulations are in good agreement with the experimental results. The effects of traps and mixed conduction are possibly the main reasons that result in the novel signal inversion of LBIC microscope at room temperature. The research results provide a theoretical guide for practical applications of large-scale array HgCdTe infrared photovoltaic detectors formed by femtosecond laser drilling, which may act as a potential new method for fabricating HgCdTe photodiodes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Accurate modeling of plasma acceleration with arbitrary order pseudo-spectral particle-in-cell methods

Particle in Cell (PIC) simulations are a widely used tool for the investigation of both laser- and beam-driven plasma acceleration. It is a known issue that the beam quality can be artificially degraded by numerical Cherenkov radiation (NCR) resulting primarily from an incorrectly modeled dispersion relation. Pseudo-spectral solvers featuring infinite order stencils can strongly reduce NCR - or even suppress it - and are therefore well suited to correctly model the beam properties. For efficient parallelization of the PIC algorithm, however, localized solvers are inevitable. Arbitrary order pseudo-spectral methods provide this needed locality. Yet, these methods can again be prone to NCR. Here in this paper, we show that acceptably low solver orders are sufficient to correctly model the physics of interest, while allowing for parallel computation by domain decomposition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

On numerical study of coupled plasma-wall instability

In this work, a rather simple approach to the numerical study of instabilities, associated with coupled plasma and wall-related transport processes, is suggested. It allows coupling of different plasma and wall-related transport models without performing additional simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2015

Climate change impacts on US agriculture and forestry: benefits of global climate stabilization

Increasing atmospheric carbon dioxide levels, higher temperatures, altered precipitation patterns, and other climate change impacts have already begun to affect US agriculture and forestry, with impacts expected to become more substantial in the future. There have been numerous studies of climate change impacts on agriculture or forestry, but relatively little research examining the long-term net impacts of a stabilization scenario relative to a case with unabated climate change. We provide an analysis of the potential benefits of global climate change mitigation for US agriculture and forestry through 2100, accounting for landowner decisions regarding land use, crop mix, and management practices. The analytic approach involves a combination of climate models, a crop process model (EPIC), a dynamic vegetation model used for forests (MC1), and an economic model of the US forestry and agricultural sector (FASOM-GHG). We find substantial impacts on productivity, commodity markets, and consumer and producer welfare for the stabilization scenario relative to unabated climate change, though the magnitude and direction of impacts vary across regions and commodities. Although there is variability in welfare impacts across climate simulations, we find positive net benefits from stabilization in all cases, with cumulative impacts ranging from $32.7 billion to $54.5 billion over the period 2015-2100. Our estimates contribute to the literature on potential benefits of GHG mitigation and can help inform policy decisions weighing alternative mitigation and adaptation actions.

agricultural modeling↗
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 OSTI2018

ExaTN

A challenging scientific problem that will benefit from exascale computing is predictive simulation of two- and three-dimensional quantum many-body Hamiltonians in condensed matter physics, together with the related strong entanglement problems in quantum chemistry and quantum computing. Efficient numerical simulations will require new insights into wavefunction compression algorithms that scale on extremely heterogeneous HPC architectures. To address this need, we propose to develop a scalable math library of parallel multi-linear algebraic numerical primitives for hierarchical tensor compression. Our library will subsequently be used for building an efficient, performance-portable HPC framework for simulating strongly correlated quantum systems. In particular, our will will be the first to deliver a massively parallel implementation of the multiscale entanglement renormalization ansatz (MERA), a hierarchical tensor decomposition quantum scheme capable of expressing local expectation values of strongly entangled higher-dimensional quantum systems efficiently. To address the extreme complexity and high computational cost of such schemes, we will automate their parallel implementation by delivering an integrated software framework, ExaTN (Exascale Tensor Networks), that will process hierarchical tensor representations on exascale HPC systems via an asymchronous task-based programming model. As a specific application, our framework will provide scientists with a capability of simulating strongly entangled quantum systems, thus positioning ORNL as a strategic leader in modeling of strongly correlated materials.

McCaskey, AlexanderJ↗