Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuous”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Scattering in black hole backgrounds and higher-spin amplitudes. Part II

We continue to investigate correspondences between, on the one hand, scattering amplitudes for massive higher-spin particles and gravitons in appropriate quantum-to-classical limits, and on the other hand, classical gravitational interactions of spinning black holes according to general relativity. We first construct an ansatz for a gravitational Compton amplitude, at tree level, constrained only by locality, crossing symmetry, unitarity and consistency with the linearized-Kerr 3-point amplitude, to all orders in the black hole’s spin. We then explore the extent to which a unique classical Compton amplitude can be identified by comparing with the results of the classical process of scattering long-wavelength gravitational waves off an exact Kerr black hole, determined by appropriate solutions of the Teukolsky equation. Up to fourth order in spin, we find complete agreement with a previously conjectured exponential form of the tree-level Compton amplitude. At higher orders, we extract tree-level contributions from the Teukolsky amplitude by an analytic continuation from a physical (a/GM < 1) to a particle-like (a/GM > 1) regime. Up to the sixth order in spin, we identify a unique conservative part of the amplitude which is insensitive both to the choice of boundary conditions at the black hole horizon and to branch choices in the analytic continuation. The remainder of the amplitude is determined modulo an overall sign from a branch choice, with the sign flipping under exchanging purely ingoing and purely outgoing boundary conditions at the horizon. Along the way, we make contact with novel applications of massive spinor-helicity variables pertaining to their relation to EFT operators and (spinning) partial amplitudes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Low-dimensional de Sitter quantum gravity

We study aspects of Jackiw-Teitelboim (JT) quantum gravity in two-dimensional nearly de Sitter (dS) spacetime, as well as pure de Sitter quantum gravity in three dimensions. These are each theories of boundary modes, which include a reparameterization field on each connected component of the boundary as well as topological degrees of freedom. In two dimensions, the boundary theory is closely related to the Schwarzian path integral, and in three dimensions to the quantization of coadjoint orbits of the Virasoro group. Using these boundary theories we compute loop corrections to the wavefunction of the universe, and investigate gravitational contributions to scattering. Along the way, we show that JT gravity in dS2 is an analytic continuation of JT gravity in Euclidean AdS2, and that pure gravity in dS3 is a continuation of pure gravity in Euclidean AdS3. We define a genus expansion for de Sitter JT gravity by summing over higher genus generalizations of surfaces used in the Hartle-Hawking construction. Assuming a conjecture regarding the volumes of moduli spaces of such surfaces, we find that the de Sitter genus expansion is the continuation of the recently discovered AdS genus expansion. Then both may be understood as coming from the genus expansion of the same double-scaled matrix model, which would provide a non-perturbative completion of de Sitter JT gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

(2, 2) Scattering and the celestial torus

Analytic continuation from Minkowski space to (2, 2) split signature spacetime has proven to be a powerful tool for the study of scattering amplitudes. Here we show that, under this continuation, null infinity becomes the product of a null interval with a celestial torus (replacing the celestial sphere) and has only one connected component. Spacelike and timelike infinity are time-periodic quotients of AdS 3 . These three components of infinity combine to an S 3 represented as a toric fibration over the interval. Privileged scattering states of scalars organize into SL(2, $\mathbb{R}$) L ×SL(2, $\mathbb{R}$) R conformal primary wave functions and their descendants with real integral or half-integral conformal weights, giving the normally continuous scattering problem a discrete character.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Black holes in Klein space

The analytic continuation of the general signature (1, 3) Lorentzian Kerr-Taub-NUT black holes to signature (2, 2) Kleinian black holes is studied. Their global structure is characterized by a toric Penrose diagram resembling their Lorentzian counterparts. Kleinian black holes are found to be self-dual when their mass and NUT charge are equal for any value of the Kerr rotation parameter a. Remarkably, it is shown that the rotation a can be eliminated by a large diffeomorphism; this result also holds in Euclidean signature. The continuation from Lorentzian to Kleinian signature is naturally induced by the analytic continuation of the S-matrix. Indeed, we show that the geometry of linearized black holes, including Kerr-Taub-NUT, is captured by (2, 2) three-point scattering amplitudes of a graviton and a massive spinning particle. This stands in sharp contrast to their Lorentzian counterparts for which the latter vanishes kinematically and enables a direct link to the S-matrix.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Volcano infrasound: progress and future directions

Over the past two decades (2000–2020), volcano infrasound (acoustic waves with frequencies less than 20 Hz propagating in the atmosphere) has evolved from an area of academic research to a useful monitoring tool. As a result, infrasound is routinely used by volcano observatories around the world to detect, locate, and characterize volcanic activity. It is particularly useful in confirming subaerial activity and monitoring remote eruptions, and it has shown promise in forecasting paroxysmal activity at open-vent systems. Fundamental research on volcano infrasound is providing substantial new insights on eruption dynamics and volcanic processes and will continue to do so over the next decade. The increased availability of infrasound sensors will expand observations of varied eruption styles, and the associated increase in data volume will make machine learning workflows more feasible. More sophisticated modeling will be applied to examine infrasound source and propagation effects from local to global distances, leading to improved infrasound-derived estimates of eruption properties. Future work will use infrasound to detect, locate, and characterize moving flows, such as pyroclastic density currents, lahars, rockfalls, lava flows, and avalanches. Infrasound observations will be further integrated with other data streams, such as seismic, ground- and satellite-based thermal and visual imagery, geodetic, lightning, and gas data. The volcano infrasound community should continue efforts to make data and codes accessible and to improve diversity, equity, and inclusion in the field. In summary, the next decade of volcano infrasound research will continue to advance our understanding of complex volcano processes through increased data availability, sensor technologies, enhanced modeling capabilities, and novel data analysis methods that will improve hazard detection and mitigation.

58 GEOSCIENCES↗
DOE OSTI2021

Convergence of sum-up rounding schemes for cloaking problems governed by the Helmholtz equation

In this work, we consider the problem of designing a cloak for waves described by the Helmholtz equation from an integer programming point of view. The problem can be modeled as a PDE-constrained optimization problem with integer-valued control inputs that are distributed in the computational domain. A first-discretize-then-optimize approach results in a large-scale mixed-integer nonlinear program that is in general intractable because of the large number of integer variables that arise from the discretization of the domain. Instead, we propose an efficient algorithm that is able to approximate the local infima of the underlying nonconvex infinite-dimensional problem arbitrarily close without the need to solve the discretized finite-dimensional integer programs to optimality. We optimize only the continuous relaxations of the approximations for local minima and then apply the sum-up rounding methodology to obtain integer-valued controls. If the solutions of the discretized continuous relaxations converge to a local minimizer of the continuous relaxation, then the resulting discrete-valued control sequence converges weakly \(^*\) in \(L^\infty\) to the same local minimizer. These approximation properties follow under suitable refinements of the involved discretization grids. Our results use familiar concepts arising from the analytical properties of the underlying PDE and complement previous results, derived from a topology optimization point of view.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Phase Inversion in a Two-phase, BCC+B2, Refractory High Entropy Alloy

A phenomenon of "phase inversion", presumably the first ever experimental evidence in metallic alloys, is shown in a refractory high entropy alloy (RHEA), Al0.5NbTa0.8Ti1.5V0.2Zr. Phase inversion in crystalline solid systems is driven by the differences in elastic modulus of the two phases. Quenching from a high-temperature single phase field, the RHEA exhibits a co-continuous mixture of a disordered BCC and an ordered B2 phase, that upon isothermal annealing at 600°C develops via spinodal decomposition into a continuous B2 matrix with discrete cuboidal BCC precipitates aligned along the <001> directions. Longer term annealing at 600°C results in the development of necking constrictions along the B2 channels, eventually pinching-off these channels and making the BCC phase continuous with discrete B2 precipitates. This inversion the process can be related to the simultaneous operation of two processes:(i) spheroidization of the initially discrete cuboidal BCC precipitates driven by a reduction in the total interface energy and (ii) an increase in the stiffness of the B2 phase, relative to the BCC phase, due to chemical composition changes during annealing, forcing the B2 regions to become discrete driven by the reduction in the total elastic strain energy.

Soni, Vishal↗
DOE OSTI2022

Scaling up carboxylic acid production from cheese whey and brewery wastewater via methane-arrested anaerobic digestion

In a circular economy, organic waste streams are valuable resources for sustainably producing chemicals and fuels. This study investigates a new methane-arrested anaerobic digestion (MAAD) process that converts high-strength cheese whey and brewery wastewater into carboxylic acids. The process was developed and optimized under various bench-scale semi-continuous (fed-batch) operating conditions (e.g., retention time, organic loading rate, pH, and feed/harvest frequency). The MAAD responses to various control-failure scenarios were also systematically investigated. The highest total acid productivity was 26g/(L liq ·d) with a substrate conversion of 0.79g COD digested /g COD fed at a hydraulic retention time (HRT) of approximately 2 d in a 14-L digester. The most stable conditions for digester operation (HRT 3 d at pH 6.0 and 40°C) were selected for process scale-up to 100gal (~380 L). Semi-continuous, pilot-scale MAAD successfully produced a total acid concentration of 40.6+/-1.1 g/L with 8.1% acetic acid, 45.1% butyric acid, and 44.5% lactic acid. The links between wastewater characteristics, operation mode, digester scale, and microbial community structure were statistically analyzed. The results show genera Sporolactobacillus and Clostridium positively correlate with butyric acid production (Pearson correlation coefficient>0.5). Moreover, four kinetic models were developed and fit to batch MAAD experimental datasets (R 2 >95%) and were successfully applied to predict the total acid production in both bench-and pilot-scale semi-continuous MAADs. In conclusion, this study shows MAAD has the potential for industrial-scale applications and is a robust platform to valorize low- or negative-value waste streams into high-value bioproducts.

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

Isogeometric Bézier dual mortaring: The Kirchhoff–Love shell problem

Here in this paper we develop an isogeometric Bézier dual mortar method for coupling multi-patch Kirchhoff–Love shell structures. The proposed approach weakly enforces the continuity of the solution at patch interfaces through a dual mortar method and can be applied to both conforming and non-conforming discretizations. As the employed dual basis functions have local supports and satisfy the biorthogonality property, the resulting stiffness matrix is sparse. In addition, the coupling accuracy is optimal because the dual basis possesses the polynomial reproduction property. We also formulate the continuity constraints through the Rodrigues’ rotation operator which gives a unified framework for coupling patches that are intersected with G 1 continuity as well as patches that meet at a kink. Several linear and nonlinear examples demonstrated the performance and robustness of the proposed coupling techniques.

42 ENGINEERING↗
DOE OSTI2025

Sparse chronology strategy for integrating seasonal energy storage in capacity expansion models

Here, this study develops the sparse chronology method to enhance the representative period framework in capacity expansion models, enabling the effective integration of long-duration energy storage modeling. Traditional representative period methods cannot capture the state of charge of seasonal energy storage systems because they do not establish effective inter-day linkages to connect the state of charge between periods. The sparse chronology approach addresses this limitation by establishing inter-day linkages that allow state of charge to shift inter-seasonally. At the same time, it groups identical representative days into partitions, applying constraints sparsely and implicitly to reduce computational load further. Validation results demonstrate that this method successfully simulates long-duration energy storage patterns, achieving close alignment with a continuous yearly benchmark model, with seasonal trends and state of charge cycles clearly represented. The computational load analysis reveals that the sparse chronology method efficiently applies constraints on maximum and minimum state of charge limits within the representative day framework, eliminating the need for detailed constraints on each individual day. By partitioning representative days and constraining only the start and end of each partition, the method significantly decreases computational requirements. Simulation results show that sparse chronology closely approximates the continuous yearly method's accuracy, even with as few as 20 representative days, achieving correlation values with the benchmark of nearly 0.9 in state of charge plots. Furthermore, it maintains computational efficiency, requiring only 4 % of the solver time compared to the continuous yearly method with 20 representative days. This approach allows capacity expansion models to incorporate long-duration energy storage with high temporal, spatial, and technological resolution, enabling more detailed modeling for large-scale power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

A Physics-Constrained Deep Learning Model for Simulating Multiphase Flow in 3D Heterogeneous Porous Media

Physics-based simulators for multiphase flow in porous media emulate nonlinear processes with coupled physics, and usually require extensive computational resources for software development, maintenance and simulation execution. As a result, a huge demand exists for fast modeling of coupled processes in a wide range of subsurface applications including geological sequestration, hydrocarbon recovery and geothermal energy extraction. In this work, an efficient physics-constrained deep learning model is developed for solving multiphase flow in 3-Dimensional (3D) heterogeneous porous media. The model fully leverages the spatial topology predictive capability of convolutional neural networks, specifically U-Net with successive contracting and expansive steps, and is coupled with an efficient continuity-based smoother to predict flow responses that need spatial continuity. Furthermore, the transient regions are penalized to steer the training process such that the model can accurately capture flow in these regions. The model takes inputs including properties of porous media, fluid properties and well controls, and predicts the temporal-spatial evolution of the state variables (pressure and saturation). While maintaining the continuity of fluid flow, the 3D spatial domain is decomposed into 2D images for reducing training cost, and the decomposition results in an increased number of training data samples and better training efficiency. Additionally, a surrogate model is separately constructed as a postprocessor to calculate well flow rate based on the predictions of state variables from the deep learning model. We use the example of CO 2 injection into saline aquifers, and apply the physics-constrained deep learning model that is trained from physics-based simulation data and emulates the physics process. The model performs prediction with a speedup of ~ 1400 times compared to physics-based simulations, and the average temporal errors of predicted pressure and saturation plumes are 0.27% and 0.099% respectively. Furthermore, water production rate is efficiently predicted by a surrogate model for well flow rate, with a mean error less than 5%. Therefore, with its unique scheme to cope with the fidelity in fluid flow in porous media, the physics-constrained deep learning model can become an efficient predictive model for computationally demanding inverse problems or other coupled processes.

58 GEOSCIENCES↗
DOE OSTI2026

Reservoir-scale model of geologic hydrogen production from serpentinization: Cyclic injection in a dual-permeability fracture-matrix system

Geologic hydrogen (GeoH 2 ) from serpentinization is a promising low-carbon resource, but its reservoir-scale behavior remains poorly understood. We develop a dual-permeability reactive transport model for an injector–producer well pair in ultramafic rock that couples multiphase flow, heat transfer, geochemistry, and porosity–permeability evolution. Here, the model is calibrated to olivine flow-through experiments and upscaled to two-year long simulations with continuous injection and cyclic injection with shut-in-to-injection ratios (SIR = 1, 0.5, 0.1). Olivine reacts along high-flux pathways to form lizardite and magnetite, increasing pH; H 2 (aq) and H 2 (g) peak early and then decline as exsolution and advective export outpace local generation. Continuous injection yields the highest cumulative H 2 but the lowest water-use efficiency (1.138 x 10 –6 mol/kgw). Cyclic injection increases this ratio to 1.35 x 10 –6 , 1.377 x 10 –6 , and 1.182 x 10 –6 mol/kgw for SIR = 1, 0.5, and 0.1, respectively; SIR = 0.5 provides a ~20% improvement over continuous injection and the best compromise for pilot design.

08 HYDROGEN↗
DOE OSTI2021

Modeling of fission and activation products in molten salt reactors and their potential impact on the radionuclide monitoring stations of the International Monitoring System

Molten Salt Reactors (MSRs) are one of six Generation IV reactor designs currently under development around the world. Because of the unique operating conditions of MSRs, which include molten fuel and the continuous removal of gaseous fission products during operation, work was performed to analyze the potential impact of emissions on the International Monitoring System (IMS) of the Comprehensive Nuclear-Test-Ban Treaty (CTBT). Simulations were performed to predict the production of IMS-relevant radionuclides in four MSR designs operating under two scenarios: (1) a sealed reactor with releases only during operational shutdown, and (2) continuous reprocessing or sparging of the fuel salt. From these production estimates the radioxenon and radioiodine signatures were extracted and compared to three current reactor designs (Pressurized Water Reactor (PWR), BWR, RBMK). In cases where continuous reprocessing of the fuel salt occurred, both the radioxenon and radioiodine signatures were nearly indistinguishable from a nuclear explosion. Estimates were also made of the potential emission rate of radioxenon for three reactor designs and it was found that MSRs have the potential to emit radioxenon isotopes at a rate of 10^15-8×10^16 Bq/d for 133Xe if no abatement is used. An assessment was made of activation products using a candidate fuel salt (FLiBe) mixed with corrosion products for the Thorium Molten Salt Reactor (TMSR-LF1).

Johnson, Christine M.↗
DOE OSTI2024

Accelerating Multivariate Functional Approximation Computation with Domain Decomposition Techniques⋆

Modeling large datasets through Multivariate Functional Approximations (MFA) provide an elegant way to handle many visualization and scientific analysis workflows. The process necessitates scalable data partitioning methods to compute MFA representations efficiently without compromising the accuracy or continuity of the reconstructed solution. We propose a domain -decomposed method for computing the MFA with B -spline bases, which reduces the total work per task and uses a restricted Additive Schwarz (RAS) method to converge the control point data degrees -of -freedom along subdomain boundaries. We provide an in-depth analysis of the parallel approach with domain decomposition solvers, aiming to minimize local subdomain error residuals and recover high -order continuity at subdomain interfaces with appropriate choices of knot overlaps. The communication cost, determined by the overlap regions in the RAS implementation, is optimized to recover the numerical error profile of the single subdomain case. Our proposed method stands in contrast to previous methods, which typically only recover either C 0 or at best C 1 continuity for arbitrary B -spline degree expansions, or those that require post -processing to blend discontinuities in the reconstructed data. We demonstrate the effectiveness of our approach using analytical and real -world datasets in 1D, 2D, and 3D through both strong and weak scaling studies. The performance results indicate that the overall cost of computing the approximation is directly proportional to the underlying nearest -neighbor communication implementation, and is only weakly dependent on the overlap region size that determines the size of the messages. This finding underscores the efficiency and scalability of our proposed method, making it a promising solution for handling large datasets in scientific workflows.

additive Schwarz solvers↗
DOE OSTI2022

Discharge intermittency considerably changes ZnO spatial distribution in porous Zn anodes

Porous Zn anodes are ubiquitous in primary batteries and are under development for low-cost rechargeable batteries. Spatial distribution of ZnO discharge product is a critical factor in these, because it can passivate the active material. In rechargeable cells this is related to a major failure mechanism called shape change, in which ZnO is relocated to inactive locations. In this work we demonstrate rest steps during discharge of primary Zn anodes dramatically alter the placement of ZnO in the anode. In alkaline electrolyte, ZnO discharge product is typically modeled as precipitating close to Zn particles, forming a porous ZnO shell around the Zn core. Further, anodes discharged continuously at low-rate are compared to anodes similarly discharged intermittently, using in situ computed tomography from a synchrotron source. Zn–ZnO core-shell structures are produced during continuous discharge but are not found in cells discharged intermittently. Continuously discharged cells showed that ZnO was formed most strongly near the separator, in agreement with Zn anode battery models. In pulse-discharged cells, ZnO was more radially distributed and was in large formations not physically connected to Zn particles. Thus, discharge intermittency changes spatial distribution of ZnO in ways that are unpredicted by Zn anode models.

25 ENERGY STORAGE↗
DOE OSTI2026

Evaluation of high-temperature properties of binary, eutectic Ni-X alloys (X = B, Ca, Ce, La, Y, Zr)

We investigated the high-temperature behavior of binary eutectic nickel-base alloys with low-solubility elements, towards developing alloys with new combinations of structural and functional (e.g. thermal or corrosion) properties. Based on thermodynamic screening, we cast binary near-eutectic Ni-B, Ni-Ca, Ni-Ce, Ni-La, Ni-Y, and Ni-Zr alloys and characterized their microstructure and high-temperature properties. All alloys (except Ni-Ca) contain a fine eutectic microstructure, containing submicron alternating lamellae of γ-Ni and continuous intermetallic phases. The eutectic microstructure did not significantly coarsen at 700°C (except for Ni-B), but measurably coarsened in all alloys at 900°C, most quickly in Ni-B, followed by Ni-Ce, Ni-La, Ni-Y, and Ni-Zr. Coarsening in these alloys occurs by fault migration; diffusion within the intermetallic, and its interfacial energy, likely influence coarsening kinetics. At room-temperature, all alloys possess negligible ductility due to the continuous, brittle intermetallic. At 700 and 900°C, the Ni-Zr alloy showed the best combination of strength and ductility. Oxidation experiments at 900°C on Ni-Zr and Ni-Ce revealed heavy internal oxidation of the continuous intermetallic phase. Overall, the Ni-Zr alloy showed good high-temperature coarsening resistance and strength, though ductility and oxidation resistance require improvement for use in structural applications. Possible routes for improving alloy properties via additive manufacturing are discussed.

Ekaputra, Clement [ORNL] (ORCID:0000000307166479)↗
DOE OSTI2023

Moisture-driven CO 2 pump for direct air capture

Direct Air Capture, the chemical separation of carbon dioxide (CO 2 ) from air, is considered a necessary approach to generate negative carbon emissions and in turn limit global warming, though energetic and monetary costs must continue to decline to meet the needed gigaton scale. The promise of energy-efficient and continuous carbon dioxide membrane separation has only recently been investigated with the development of facilitated chemical transport and novel electrochemically driven methods. This work demonstrates the potential for continuous pumping of CO 2 from air using commercial anion exchange membranes driven solely by a water vapor gradient. A reactive transport model that couples ionic transport with moisture sensitive bicarbonate chemistry in charged polymers, also termed moisture swing sorption, was developed alongside experimentally measured equilibrium and transport properties of CO 2 and water vapor in Fumasep FAA-3 to determine the mass transport limits of such a system. The results show that the commercial membrane is kinetically limited by the moisture sensitive chemistry, enabling a CO 2 current of 1.1 μmol/m 2 • s. In diffusion limited materials, a CO 2 pumping flux can reach 11 μmol/m 2 • s in and peak performance is expected with a dry air relative humidity around 50% pumping into a saturated water vapor permeate. In all configurations, water loss becomes the major cost in such a separation and must be managed through optimal polymer design or water recovery mechanisms. The work demonstrates the applicability of a new driving force, the evaporation of water, to meet the energetic demand of CO 2 separation from air.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Necessary and sufficient conditions for quadratic stabilizability of switched systems on non-uniform time domains

In this paper, we consider the quadratic stabilizability via state feedback for a particular class of switched systems that evolve on a non-uniform time domain by introducing time scales theory. The system considered switches between a continuous-time subsystem with variable lengths and a discrete-time subsystem with variable discrete step sizes. Necessary and sufficient conditions are derived to guarantee the quadratic stability of this class of switched systems via a switching state feedback law based on the existence of a common positive definite matrix satisfying the quadratic stabilizability condition by considering that the two subsystems are unstable. By state feedback, we mean that the switching among subsystems depends on the system states. Current results for this kind of state switching feedback control are derived only for switched systems evolving on a continuous time domain or a discrete time domain with fixed step’s size. These results are not applicable for the particular class of switched systems where there is a mixing between the continuous and discrete dynamics. This motivates the derivation of a new and more general state feedback control law for switched systems in this work. Here, a numerical example illustrating the results is presented.

42 ENGINEERING↗