Engineering Papers⌕ Search

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.

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

At least 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Groundwater Circulation Within the Mountain Block: Combining Flow and Transport Models With Magnetotelluric Observations to Untangle Its Nested Nature

Abstract Mountains are vital water sources for humans and ecosystems, continuously replenishing lowland aquifers through surface runoff and mountain recharge. Quantifying these fluxes and their relative importance is essential for sustainable water resource management. However, our mechanistic understanding of the flow and transport processes determining the connection between the mountain block and the basin aquifer remains limited. Traditional conceptualizations assume groundwater circulation within the mountain block is predominantly shallow. This view neglects the role of deep groundwater flowpaths significantly contributing to the water, solute, and energy budgets. Overcoming these limitations requires a holistic characterization of the multiscale nature of groundwater flow along the mountain‐to‐valley continuum. As a proof‐of‐concept, we use a coupled groundwater flow and transport model to design a series of numerical experiments that explore the role of geology, topography, and weathering rates in groundwater circulation and their resulting resistivity patterns. Our results show that accumulating solutes near stagnation zones create contrasting electrical resistivity patterns that separate local, intermediate, and regional flow cells, presenting a target for magnetotelluric observations. To demonstrate the sensitivity of magnetotelluric data to features in our resistivity models, we use the MARE2DEM electromagnetic modeling code to perform forward and inverse simulations. This study highlights the potential of magnetotelluric surveys to image the resistivity structure resulting from multiscale groundwater circulation through relatively impervious crystalline basement rocks in mountainous terrains. This capability could change our understanding of the critical zone, offering a holistic perspective that includes deep groundwater circulation and its role in conveying solutes and energy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Smoothed aggregation for difficult stretched mesh and coefficient variation problems

Abstract Four adaptations of the smoothed aggregation algebraic multigrid (SA‐AMG) method are proposed with an eye toward improving the convergence and robustness of the solver in situations when the discretization matrix contains many weak connections. These weak connections can cause higher than expected levels of fill‐in within the coarse discretization matrices and can also give rise to suboptimal smoothing within the prolongator smoothing phase. These smoothing drawbacks are due to the relatively small size of some diagonal entries within the filtered matrix that one obtains after dropping the weak connections. The new algorithms consider modifications to the Jacobi‐like step that defines the prolongator smoother, modifications to the filtered matrix, and also direct modifications to the resulting grid transfer operators. Numerical results are given illustrating the potential benefits of the proposed adaptations.

Hu, Jonathan J.↗
DOE OSTI2024

Preserving Tracer Correlations in Moment-Based Atmospheric Transport Models

A linear non-diffusive algorithm for advective transport is developed that greatly improves the detail at which aerosols and clouds can be represented in atmospheric models. Linear advection schemes preserve tracer correlations but the most basic linear scheme is rarely used by atmospheric modelers on account of its excessive numerical diffusion. Higher-order schemes are in widespread use, but these present new problems as nonlinear adjustments are required to avoid occurrences of negative concentrations, spurious oscillations, and other non-physical effects. Generally successful at reducing numerical diffusion during the advection of individual tracers, for example, particle number or mass, the higher-order schemes fail to preserve even the simplest of correlations between interrelated tracers. As a result, important attributes of aerosol and cloud populations including radial moments of particle size distributions, molecular precursors related through chemical equilibria, aerosol mixing state, and distribution of cloud phase are poorly represented. We introduce a new transport scheme, minVAR, that is both non-diffusive and preservative of tracer correlations, thereby combining the best features of the basic and higher-order schemes while enabling new features such as the tracking of sub-grid information at arbitrarily fine scales with high computational efficiency.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Thermodynamic approach for enhancing superconducting critical current performance

Abstract The addition of artificial pinning centers has led to an impressive increase in the critical current density ( J c ) of superconductors, enabling record-breaking all-superconducting magnets and other applications. The J c of superconductors has reached ~0.2–0.3 J d , where J d is the depairing current density, and the numerical factor depends on the pinning optimization. By modifying λ and/or ξ, the penetration depth and coherence length, respectively, we can increase J d . For (Y 0.77 Gd 0.23 )Ba 2 Cu 3 O y ((Y,Gd)123), we can achieve this by controlling the carrier density, which is related to λ and ξ. We can also tune λ and ξ by controlling the chemical pressure in Fe-based superconductors, i.e., BaFe 2 (As 1− x P x ) 2 films. The variation in λ and ξ leads to an intrinsic improvement in J c via J d , allowing extremely high values of J c of 130 MA/cm 2 and 8.0 MA/cm 2 at 4.2 K, consistent with an enhancement in J d of a factor of 2 for both incoherent nanoparticle-doped (Y,Gd)123 coated conductors (CCs) and BaFe 2 (As 1− x P x ) 2 films, showing that this new material design is useful for achieving high critical current densities in a wide array of superconductors. The remarkably high vortex-pinning force in combination with this thermodynamic and pinning optimization route for the (Y,Gd)123 CCs reached ~3.17 TN/m 3 at 4.2 K and 18 T ( H || c ), the highest values ever reported for any superconductor.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Challenges in the detection of gases in exoplanet atmospheres

Claims of detections of gases in exoplanet atmospheres often rely on comparisons between models including and excluding specific chemical species. However, the space of molecular combinations available for model construction is vast and highly degenerate. Only a limited subset of these combinations is typically explored for any given detection. As a result, apparent detections of trace gases risk being artefacts of incomplete modelling rather than robust identification of atmospheric constituents, especially in the low-signal-to-noise regime. Here, using the sub-Neptune K2-18 b as a case study, we show that recent biosignature claims vanish when the model space is expanded, with numerous alternatives providing equally good or better fits. We demonstrate that the significance of a claimed detection relies on the choice of models being compared, and that model preference does not in itself imply the presence of a specific gas. In conclusion, we recommend treating model comparisons instead as relative adequacy tests, which should be supported by theoretical predictions and complementary metrics of statistical significance to attribute a signal to a particular gas.

Astronomy and AstroPhysics↗
DOE OSTI2021

Analytic and numerical vibronic spectra from quasi-classical trajectory ensembles

The truncated Wigner approximation to quantum dynamics in phase space is explored in the context of computing vibronic line shapes for monomer linear optical spectra. We consider multiple model potential forms including a shifted harmonic oscillator with both equal and unequal frequencies on the ground and excited state potentials as well as a shifted Morse potential model. For the equal-frequency shifted harmonic oscillator model, we derive an analytic expression for the exact vibronic line shape that emphasizes the importance of using a quantum mechanical distribution of phase space initial conditions. For the unequal-frequency shifted harmonic oscillator model, we are no longer able to obtain an exact expression for the vibronic line shape in terms of independent deterministic classical trajectories. We show how one can rigorously account for corrections to the truncated Wigner approximation through nonlinear responses of the line shape function to momentum fluctuations along a classical trajectory and demonstrate the qualitative improvement in the resulting spectrum when the leading-order quantum correction is included. Finally, we numerically simulate absorption spectra of a highly anharmonic shifted Morse potential model. We find that, while finite quantization and the dissociation limit are captured with reasonable accuracy, there is a qualitative breakdown of the quasi-classical trajectory ensemble’s ability to describe the vibronic line shape when the relative shift in Morse potentials becomes large. The work presented here provides clarity on the origin of unphysical negative features known to contaminate absorption spectra computed with quasi-classical trajectory ensembles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

On the discretization error of the discrete generalized quantum master equation

The transfer tensor method (TTM) [Cerrillo and Cao, Phys. Rev. Lett. 112 , 110401 (2014)] can be considered a discrete-time formulation of the Nakajima–Zwanzig quantum master equation (NZ-QME) for modeling non-Markovian quantum dynamics. A recent paper [Makri, J. Chem. Theory Comput. 21 , 5037 (2025)] raised concerns regarding the consistency of the TTM discretization, particularly a spurious term at the initial time t = 0. Here, this work presents a detailed analysis of the discretization structure of the TTM, clarifying the origin of the initial-time correction and establishing a consistent relationship between the TTM discrete-time memory kernel K N and the continuous-time NZ-QME kernel $\mathscr{K}$( N Δ t ). This relationship is validated numerically using the spin-boson model, demonstrating convergence of reconstructed memory kernels and accurate dynamical evolution as Δ t → 0. While the TTM provides a consistent discretization, we note that alternative schemes are also viable, such as the midpoint derivative/midpoint integral scheme proposed in Makri’s work. The relative performance of various schemes for either computing accurate $\mathscr{K}$( N Δ t ) from exact dynamics or obtaining accurate dynamics from exact $\mathscr{K}$( N Δ t ) warrants further investigation.

Density-matrix↗
DOE OSTI2022

Insulating antiferromagnetism in VTe

Here, we report a detailed theoretical and experimental study on the vanadium monotelluride VTe, which crystallizes in the NiAs hexagonal structure. First-principles calculations reveal a complex hierarchy of magnetic interactions and energy scales, with the ground state theoretically determined as an ($\frac{1}{2}$, 0, $\frac{1}{2}$) antiferromagnetic ordering with insulating character and a band gap of 0.5 eV. Experimental synthesis and characterization efforts find a substantially off-stoichiometric orthorhombic structure (a defect NiAs structure) with composition V 0.85 Te, and an apparent Néel point of some 45 K. First-principles calculations find good agreement with the observed Néel point. We also give an extended examination of the effects of off-stoichiometry on the calculated energetics, finding significant volume-related effects. Our first-principles calculations find the stoichiometric phase VTe to have a negative vanadium defect formation energy of over 1 eV, thus explaining the formation of the off-stoichiometric phase. Finally, we provide a structural explanation for the formation of defect structures in this and numerous other NiAs-structure materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Minimizing Ground Risk in Cellular-Connected Drone Corridors With mmWave Links

Unmanned Aircraft Systems (UASs) have been receiving significant interest and support from academia, industry, and regulatory bodies over the past decade due to their various use cases. To safely integrate UAS operations into the national airspace, particularly overpopulated regions, the risk posed to ground users, buildings, and vehicles due to unmanned aerial vehicle (UAV) flight should be minimized. This risk can be represented by a numerical metric, which we refer to in this article as the “ground risk.” Many UAS applications also depend on the presence of a reliable wireless communication link between the UAV and a control station for the transmission of UAV position, surveillance video, UAV payload commands, and other mission-related data. Such wireless communication requirements also need to be considered in the design of UAS operations. In this article, we consider both these aspects and study the design of nonintersecting trajectories for UAS operations to minimize ground risk, subject to constraints on the wireless signal strength and geometry of the trajectory, specified in terms of: 1) an enclosing cylinder within which the trajectory must lie and 2) an integrated angular change along the UAV's trajectory. The performance of a computationally expensive optimal algorithm is compared with that of a computationally faster heuristic approach within the dense urban environment of Manhattan, NY, USA. Performance evaluation using ray-tracing simulations shows that the heuristic approach performs close to the optimal algorithm at a reduced computation cost. In conclusion, this research can be utilized to make UAS operations safe and reliable and accelerate their adoption.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

A Novel Design for Improving the Control on the Stainless-Steel Vessel Welding Process for Superconducting Magnets

Stainless steel vessels see widespread use in superconducting magnets for particle accelerator applications. Their function varies in different magnet designs: they always provide the necessary liquid helium containment, but in some cases are also used to provide azimuthal prestress and can also be welded to the magnet end plate to provide additional longitudinal stiffness. A magnet designed with the bladder and key technology does not rely on the structural role of the vessel. They are structurally supported using azimuthally prestressed aluminum shells, and the longitudinal constraint by rods. In this case, the magnet designer would generally like to minimize the interaction between the magnet and the stainless-steel vessel and to minimize the coil stress variation due to the vessel. The stress state in the vessel and in the coil is a function of the circumferential interference, defined as the vessel azimuthal length minus the magnet circumference. The vessel and the magnet azimuthal length machining tolerances are relatively large resulting in significant stress variations in the superconducting coils. Here in this paper we introduce an interference-control shim, which can significantly limit the stress variation of the coils for a given variation of the interference. The effectiveness of the interference-control shim is evaluated numerically on the MQXF, the low-β quadrupole for the High Luminosity LHC.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Scaled Vecchia Approximation for Fast Computer-Model Emulation

Many scientific phenomena are studied using computer experiments consisting of multiple runs of a computer model while varying the input settings. Gaussian processes (GPs) are a popular tool for the analysis of computer experiments, enabling interpolation between input settings, but direct GP inference is computationally infeasible for large datasets. We adapt and extend a powerful class of GP methods from spatial statistics to enable the scalable analysis and emulation of large computer experiments. Specifically, we apply Vecchia’s ordered conditional approximation in a transformed input space, with each input scaled according to how strongly it relates to the computer-model response. The scaling is learned from the data by estimating parameters in the GP covariance function using Fisher scoring. Our methods are highly scalable, enabling estimation, joint prediction, and simulation in near-linear time in the number of model runs. In several numerical examples, our approach substantially outperformed existing methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Climate Influences on Capacity Expansion Planning with Application to the Western U.S

Electric power system planners utilize a variety of planning tools to inform decisions concerning generation and transmission additions to the electric grid, the need for operational changes, and to evaluate potential stressors on the system. Numerous factors contribute to the planning process including projected fuel and technology costs, policy and load profiles. There is also a growing recognition of the interdependency of the electric grid with other natural and engineered systems. Here we explore how future climate change and hydropower operability might influence decisions related to electricity capacity expansion planning and operations. To do so we assemble a multi-model framework. Specifically, water resource modeling is used to simulate climate impacts on future water supply for thermoelectric and hydropower generation. Separately, temperature impacts on electricity load are evaluated. Together, these climate factors spatially constrain a capacity expansion model that projects generation and transmission additions to the grid. The projected new capacity-builds are then evaluated on their operations, reliability, and cost under average and extreme climate conditions using production cost modeling. This coupled framework is demonstrated on the electric grid in the Western U.S., supporting capacity expansion planning by WECC, the North American Electric Reliability Corporation (NERC) regional entity responsible for reliability assurance of the Western Interconnection. This region was selected in part because the West is unique in that it has high potential for renewable penetrations and is experiencing large retirements/displacements of baseload resources, primarily coal, leading to possible operational challenges in terms of changing resource mix and the need for resource flexibility. Toward this challenge, planning scenarios encompass a range of alternative energy, climate and drought futures. In this context we explore answers to two strategic questions: 1) How does changing climate influence electricity expansion planning (generation and transmission) and future operations, including type and capacity of new builds, system reliability, cost and environmental impacts? 2) How does the representation of hydropower in the modeling framework influence the evaluation of bulk power system operations? Results indicate that climate has a measurable influence on recommendations concerning the capacity, type and location of new generation and transmission additions, with up to 17 GW additional capacity needed by 2038 to meet peak loads (~6.6% increase over capacity-builds based on historical climate). The extent of additional infrastructure needs is strongly influenced by future water availability for hydropower and the potential deployment of demand response technologies. Systems designed for future climate conditions were found to maintain high system reliability under a range of electricity and water availability scenarios (including significant drought), with minimal system curtailments. Additional capacity needs due to higher load tend to increase cumulative 20-year investment and operating costs by $\$$5-$\$$17 billion and generation costs increase by 9 to 19%. Finally, changing the representation of hydropower flexibility has a relatively small influence on capacity expansion in the Western Interconnection through 2038, but hydropower flexibility impacts generation costs to a similar extent as climate.

13 HYDRO ENERGY↗
DOE OSTI2025

System Analysis of an Internal Combustion Engine (ICE) - Solid Oxide Fuel Cell (SOFC) Hybrid Cycle

The variability of renewable energy sources poses challenges for reliable grid operation. Conventional thermal power sources, though reliable, often lack operational flexibility. Hybrid energy systems that integrate Solid Oxide Fuel cells (SOFC) with Internal Combustion Engines (ICE) offer a promising solution by achieving relatively higher efficiency and grid-following capability. This study investigates performance of a 100-kW pressurized SOFC-ICE hybrid cycle. In this configuration, unutilized SOFC fuel is used to drive the engine, with a turbocharger providing air supply and an external reformer generating syngas. System components were numerically modeled using MATLAB/SIMULINK for the SOFC and reformer, and EBSILON® for the ICE and balance of plant. Parametric studies varied fuel utilization (70-90%), reformer temperature (600 – 1000K), anode off-gas recirculation (0 – 70%), and current density (0.2 – 0.55 A/cm2). Results show that the SOFC and ICE operate as thermally independent topping and bottoming cycles, achieving peak efficiency of 62% under optimized conditions.

hybrid↗
DOE OSTI2025

2024 Incubator Final Report: Active Hybrid Mooring

This Incubator project has defined and developed a representation of an active hybrid mooring (AHM) system in Simulink to assess whether the speed and accuracy of AHM is sufficient to warrant further investigation in a larger project. We first built a representative mooring system for the VolturnUS-S 15MW floating offshore wind turbine (FOWT) in 1,000-meter water depth to provide baseline loads and displacements needed in the AHM system, with the results showing off-the-shelf actuators would be practical. We then defined model configurations for the numerical and physical substructures, using standalone MoorDyn v2 with a novel input strategy and chained hydraulic actuators represented in Simulink Simscape, respectively. The remainder of the project has focused on verifying the existing models for speed and accuracy against OpenFAST simulations, increasing their complexity to enable six degree of freedom (6DOF) operation, and connecting the numerical and physical substructures into a single holistic model. The verification phase of the project has shown promising results, though further refinement is needed to resolve accuracy discrepancies and for robust performance across a wider range of metocean conditions. The numerical model we developed is computationally stable and operates faster than real time using a time step of 0.01 seconds or shorter, indicating the model is likely fast enough for our AHM design to operate successfully. Additionally, the Simscape model of a hydraulic actuator operates with less than 0.3% relative error when actuating mooring loads from an OpenFAST simulation of the VolturnUS-S 1,000-meter mooring system. The aforementined refinements will be performed in the next year as a Continuing Incubator, with key remaining tasks include programming the predictor-corrector loop to minimize movement error, testing against more complex nonlinear wave cases, and expanding to enable loads from connected shared mooring lines.

17 WIND ENERGY↗
DOE OSTI2026

Juvenile Salmon and Their Habitats in the Columbia River Estuary: A Review and Synthesis of Knowledge Development 2000–2025

[This is a 90% discussion draft.] This is the third Synthesis Memorandum funded by the U.S. Army Corps of Engineers and developed for the Columbia Estuary Ecosystem Restoration Program (CEERP) on the topic of habitat restoration in the Columbia River Estuary (CRE) from Bonneville Dam to the river mouth. While the first two were developed by PNNL and NOAA without the benefit of stakeholder participation, for the current memo, two key activities were initiated: (1) review, by the Expert Regional Technical Group (ERTG), of status and trends monitoring and action effectiveness monitoring funded by CEERP, and (2) a workshop including representatives of the Bonneville Power Administration and the U.S. Army Corps of Engineers (the action agencies [AAs]), the National Oceanic and Atmospheric Administration (NOAA), major research agencies contributing to CEERP, and sponsors who implement CEERP restoration actions. A systematic literature review was conducted using ClarivateTM Web of ScienceTM database. The topics of interest for CRE relevant research included salmon ecology, physical processes, and wetland habitats, and therefore required the use of broad search terms. Our final search criteria included a combination of Boolean operators and an approach to combine different sets of search terms. The final search result yielded 669 records. The records were classified by groups and assigned to the relevant disciplinary expert for review. The review identified substantive advances in understanding the provision of salmon habitat functions through spatiotemporally dynamic physical and ecological processes, and the use of CRE habitats by numerous stocks of juvenile salmon. It also uncovered heretofore unincorporated historical documentation of riparian habitats across the CRE. The characterization of the structural components of floodplain habitat including plant associations and channel networks has advanced considerably, together with the understanding of seasonal changes and long-term trends. The relative influence of salmon-habitat location in the CRE as compared with temporal factors, mainly season, has been well described, which affects the prioritization of restoration. Stressors on the ecosystem and fish, and the drivers of these stressors, have been more carefully elucidated and predictive models are in various stages of development. The vision, aims, and design of restoration projects have advanced together with methods of data collection, analysis, and modeling that have seen substantial improvements. Experiments intended to inform the design of restoration projects are underway or have been completed. An important outstanding area of research that has lagged behind the advances in fundamental understanding of the ecosystem and salmon habitat functions remains the peer-reviewed documentation of the outcomes of restoration for both habitats and fish functions.

estuary↗
DOE OSTI2020

Influence of soluble oligomeric aluminum on precipitation in the Al-KOH-H2O system

The role of oligomeric aluminate (Al(OH)4-) species in the precipitation of aluminum phases such as gibbsite (a-Al(OH)3) from aqueous hydroxide solutions remains unclear and difficult to probe directly, despite its importance to developing accurate predictions of Al solubility in highly alkaline systems. Precipitation in this system entails a transition from predominantly tetrahedrally coordinated Al species in solution to octahedrally coordinated Al in gibbsite. Here we report a quantitative study of dissolved Al in the Al-KOH-H2O system using a combination of molecular spectroscopies. We established a relationship between changes in 27Al NMR chemical shifts and the relative intensity of Raman vibrational bands, indicative of variations in the ensemble speciation of Al in solution, and the formation of unique contact ion pair interactions with the aluminate dimer, Al2O(OH)62 . A strong correlation between the extent of Al oligomerization and the amount of solvated Al was demonstrated by systematically varying the KOH:Al molar ratio. The concentration of dissolved oligomeric Al in solution also directly impacted the particle size and morphology of gibbsite; high concentrations of dimeric Al2O(OH)62-, yielded smaller and more numerous anhedral to subhedral gibbsite particles, while low concentrations yielded fewer and larger euhedral gibbsite platelets. The collective observations suggest a key role of the Al2O(OH)62- dimer in promoting gibbsite precipitation from solution, where the potassium ion-paired dimer possibly catalyzes a more rapid transformation of Al coordination from tetrahedral in solution to octahedral in gibbsite.

Dembowski, Mateusz↗
DOE OSTI2024

Thermodynamics of the near-extremal Kerr spacetime

We examine the thermodynamics of a near-extremal Kerr black hole, and demonstrate that the geometry behaves as an ordinary quantum system with a vanishingly small degeneracy at low temperatures. This is in contrast with the classical analysis, which instead predicts a macroscopic entropy for the extremal Kerr black hole. Our results follow from a careful analysis of the gravitational path integral. Specifically, the low temperature canonical partition function behaves as Z ~ T$\frac{3}{2}$ e S 0 +c log S 0 , with S 0 the classical degeneracy and c a numerical coefficient we compute. This is in line with the general expectations for non-supersymmetric near-extremal black hole thermodynamics, as has been clarified in the recent past, although cases without spherical symmetry have not yet been fully analyzed until now. We also point out some curious features relating to the rotational zero modes of the near-extremal Kerr black hole background that affects the coefficient c. This raises a puzzle when considering similar black holes in string theory. Our results generalize to other rotating black holes, as we briefly exemplify.

79 ASTRONOMY AND ASTROPHYSICS↗