SEARCH · Search NASA
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.
The objective of these notes is to demonstrate that the formulation of the MHD Boussinesq convection equations that follows the conventions of the turbulence community used in our previous work is consistent and compatible with the formulation from the convection community used in Busse and Pesch.
We developed a biomechanics-informed online learning framework to learn the dynamics with ground truth generated with multiscale modeling simulation. It was built on Summit-like supercomputers, which were also used to benchmark and validate our framework on one physiologically significant modeling of deformable biological cells. We generalized the century-old equation of Jeffery orbits to a new equation of motion with additional parameters to account for the flow conditions and the cell deformability. Using simulation data at particle-based resolutions for flowing cells and the learned parameters from our framework, we validated the new equation by the motions, mostly rotations, of a human platelet in shear blood flow at various shear stresses and platelet deformability. Our online framework, which surrogates redundant computations in the conventional multiscale modeling by solutions of our learned equation, accelerates the conventional modeling by three orders of magnitude without visible loss of accuracy.
The deep operator network (DeepONet) has shown remarkable potential in solving partial differential equations (PDEs) by mapping between infinite-dimensional function spaces using labeled datasets. However, in scenarios lacking labeled data, the physics-informed DeepONet (PI-DeepONet) approach, which utilizes the residual loss of the governing PDE to optimize the network parameters, faces significant computational challenges, particularly due to the curse of dimensionality. This limitation has hindered its application to high-dimensional problems, making even standard 3D spatial with 1D temporal problems computationally prohibitive. Additionally, the computational requirement increases exponentially with the discretization density of the domain. Here, to address these challenges and enhance scalability for high-dimensional PDEs, we introduce the Separable physics-informed DeepONet (Sep-PI-DeepONet). This framework employs a factorization technique, utilizing sub-networks for individual one-dimensional coordinates, thereby reducing the number of forward passes and the size of the Jacobian matrix required for gradient computations. By incorporating forward-mode automatic differentiation (AD), we further optimize computational efficiency, achieving linear scaling of computational cost with discretization density and dimensionality, making our approach highly suitable for high-dimensional PDEs. We demonstrate the effectiveness of Sep-PI-DeepONet through three benchmark PDE models: the viscous Burgers’ equation, Biot’s consolidation theory, and a parameterized heat equation. Our framework maintains accuracy comparable to the conventional PI-DeepONet while reducing training time by two orders of magnitude. Notably, for the heat equation solved as a 4D problem, the conventional PI-DeepONet was computationally infeasible (estimated 289.35 h), while the Sep-PI-DeepONet completed training in just 2.5 h. These results underscore the potential of Sep-PI-DeepONet in efficiently solving complex, high-dimensional PDEs, marking a significant advancement in physics-informed machine learning.
Dust poses a serious threat to tokamak operation and safety. It is important to study the behaviour of dust grains under tokamak's discharge conditions, which depends heavily on their size and charge. Existing simulations mainly address issues on dust grains with radii larger than 1 μm, in which case, the drift effect due to electromagnetic fields can be safely ignored. For nanometer scale dust grains, however, the drift effect becomes significant and a new model based on guiding-centre system needs to be established. In this work, the NDS has been done under BOUT++ framework. The simulation contains two parts. Part one, NDS evaluates the charging and ablation processes of the dust grains. In the second part, the guiding-centre orbits of dust particles are tracked in tokamak plasmas, whose parameters are obtained from BOUT++, a highly desirable C++ code package for performing parallel plasma fluid simulations with an arbitrary number of equations in 3D curvilinear coordinates. The orbit of nanodust dynamics is described by guiding centre equations for simplicity, and these equations are numerically solved by conventional fourth-order Runge Kutta method. Simulations provide results such as trajectories and evolutions of dust particles with different sizes and velocities for different tokamak geometries. Results show tungsten dust grains with a radius of a few nanometers launched from outer midplane will oscillate before totally ablated in C-Mod. The oscillation in this case is driven by the ion drag force. Larger Nanodust with a radius of 100 nm, on the contrary, cannot be completely constrained by the electromagnetic field. In conclusion, the high plasma temperature and density in the seperatrix region causes severe dust ablation, resulting in total ablation within several ms.
Full-waveform inversion (FWI) is an accurate imaging approach for modeling the velocity structure by minimizing the misfit between recorded and predicted seismic waveforms. However, the strong nonlinearity of FWI resulting from fitting oscillatory waveforms can trap the optimization in local minima. We have adopted a neural-network-based full-waveform inversion (NNFWI) method that integrates deep neural networks with FWI by representing the velocity model with a generative neural network. Neural networks can naturally introduce spatial correlations as regularization to the generated velocity model, which suppresses noise in the gradients and mitigates local minima. Furthermore, the velocity model generated by neural networks is input to the same partial differential equation (PDE) solvers used in conventional FWI. The gradients of the neural networks and PDEs are calculated using automatic differentiation, which back propagates gradients through the acoustic PDEs and neural network layers to update the weights of the generative neural network. Experiments on 1D velocity models, the Marmousi model, and the 2004 BP model determine that NNFWI can mitigate local minima, especially for imaging high-contrast features such as salt bodies, and it significantly improves the inversion in the presence of noise. Adding dropout layers to the neural network model also allows analyzing the uncertainty of the inversion results through Monte Carlo dropout. NNFWI opens a new pathway to combine deep learning and FWI for exploiting the characteristics of deep neural networks and the high accuracy of PDE solvers. Because NNFWI does not require extra training data and optimization loops, it provides an attractive and straightforward alternative to conventional FWI.
A technique originating in quantum dynamics is used to derive a trajectory-based, semi-analytical solution for two-phase flow. The partial differential equation governing the evolution of the aqueous phase is equivalent to a family of ordinary differential equations defined along a path through the porous medium. The trajectories may be found by solving the differential equations directly or by post-processing the output of a numerical solution to the full set of governing equations. The trajectories, which differ from conventional streamlines, are found to bend downward in response to gravitational forces. The curvature is more pronounced as the dip of the porous layer containing the flow increases. Subtle changes in the relative permeability curve can lead to significant variations in the trajectories. The ordinary differential equation for the trajectory provides an expression for the travel time along the path. The expression produces a semi-analytical approximation to the model parameter sensitivities, the partial derivatives of the travel times with respect to changes in the permeability model. The semi-analytical trajectory-based sensitivities generally agree with those computed using a numerical reservoir simulator and a perturbation approach. The sensitivities are useful in tomographic imaging algorithms designed to estimate the spatial variation in permeability within a porous medium using multiphase observations.
While recent advances in deep learning have shown promising efficiency gains in solving time-dependent partial differential equations (PDEs), matching the accuracy of conventional numerical solvers still remains a challenge. One strategy to improve the accuracy of deep learning-based solutions for time-dependent PDEs is to use the learned model as the coarse propagator in the Parareal method and a traditional numerical method as the fine solver. However, successful integration of deep learning into the Parareal method requires consistency between the coarse and fine solvers, particularly for PDEs exhibiting rapid changes such as sharp transitions. Here, to ensure this consistency, we propose using convolutional neural networks (CNNs) to learn the fully discrete time-stepping operator defined by the same numerical scheme employed as the fine solver. We demonstrate the effectiveness of the proposed method in solving the classical and mass-conservative Allen–Cahn (AC) equations. Through iterative updates in the Parareal algorithm, our approach achieves a significant computational speedup compared to traditional fine solvers while converging to high-accuracy solutions. Our results highlight that the proposed hybrid Parareal algorithm effectively accelerates simulations, particularly when implemented on multiple GPUs, and converges to the desired accuracy in only a few iterations. Another advantage of our method is that the CNN model is trained on trajectory-based data generated from random initial conditions, such that the trained model can be used to solve the AC equations with various initial conditions without retraining. This work demonstrates the potential of integrating neural network methods into parallel-in-time frameworks for efficient and accurate simulations of time-dependent PDEs.
This short paper is organized around a few discrete questions. How have U.S. perspectives on conventional and nuclear interactions evolved in recent years? How have these changing perspectives shaped U.S. thinking on perceived nuclear escalation risks with Russia, both long standing and newly emerged? And finally, what do these U.S. perceptions of interactions and risks mean for future potential arms control either in the conventional or nuclear arena? The answers to these questions appear to show some promise for future work in this area. Conventional and nuclear interactions are increasing, and thus the two domains cannot be kept as deliberately or artificially separated as in the past. There is a growing mutual recognition regarding the interplay between conventional and nuclear capabilities in the strategic stability equation. There is also an increasingly shared recognition in the United States and Russia that nuclear escalation risks are likely to originate in conventional crises or conflict, placing an imperative on some form of conflict prevention procedures at the lower end of the spectrum. This requirement suggests some potential areas for risk reduction measures outside of those traditionally considered in conventional and nuclear arms control over the past several decades.
The behavior of suspended particles in turbulent flows is a recalcitrant problem spanning wide-ranging fields including geomorphology, hydrology, and dispersion of particulate matter in the atmosphere. One key mechanism underlying particle suspension is the difference between particle settling velocity (w s ) in turbulence and its still water counterpart (w so ). This difference is explored here for a range of particle-to-fluid densities (1 – 10) and particle diameter to Kolmogorov micro-eddy sizes (0.1 – 10). Conventional models of particle fluxes that equate w s to w so result in eddy diffusivities and turbulent Schmidt numbers contradictory to laboratory experiments. Incorporating virtual mass and Basset history forces resolves these inconsistencies, providing clarity as to why w s /w so is sub-unity for the aforementioned conditions. The proposed formulation can be imminently used to model particle settling in turbulence, especially when sediment distribution outcomes over extended time scales far surpassing turbulence time scales are sought.
With only less than 10% recovery, the primary production of hydrocarbon from shale reservoirs has redefined the energy equation in the world. Similar to conventional reservoirs, Enhanced Oil Recovery (EOR) techniques could be devised to enhance the current recovery factors. However, shale reservoirs possess unique characteristics that significantly affect the fluid properties. Therefore, we are adopting a molecular simulation approach that is well-suited to account for these effects to evaluate the performance of three different gases, methane, carbon dioxide and nitrogen, to recover the hydrocarbons from rough pore surfaces. Our hydrocarbon systems consists of either a single component (decane) or more than one component (decane and pentane). We simulated cases where concurrent and countercurrent displacement is studied. For concurrent displacement (injected fluids displace hydrocarbons towards the production region), we found that nitrogen and methane yielded similar recovery; however nitrogen exhibited a faster breakthrough. On the other hand, carbon dioxide was more effective in extracting the hydrocarbons when sufficient pressure was maintained. For countercurrent displacement (gases are injected and hydrocarbons are produced from the same direction), methane was found to be more effective, followed by carbon dioxide and nitrogen. In all cases, confinement reduced the recovery factor of all gases. This work provides insights to devise strategies to improve the current recovery factors observed in shale reservoirs.
A method for machine learning and serving of discrete field theories in physics is developed. The learning algorithm trains a discrete field theory from a set of observational data on a spacetime lattice, and the serving algorithm uses the learned discrete field theory to predict new observations of the field for new boundary and initial conditions. The approach of learning discrete field theories overcomes the difficulties associated with learning continuous theories by artificial intelligence. The serving algorithm of discrete field theories belongs to the family of structure-preserving geometric algorithms, which have been proven to be superior to the conventional algorithms based on discretization of differential equations. The effectiveness of the method and algorithms developed is demonstrated using the examples of nonlinear oscillations and the Kepler problem. In particular, the learning algorithm learns a discrete field theory from a set of data of planetary orbits similar to what Kepler inherited from Tycho Brahe in 1601, and the serving algorithm correctly predicts other planetary orbits, including parabolic and hyperbolic escaping orbits, of the solar system without learning or knowing Newton’s laws of motion and universal gravitation. The proposed algorithms are expected to be applicable when the effects of special relativity and general relativity are important.
Here, we implement Crossing Statistics to reconstruct in a model-agnostic manner the expansion history of the universe and properties of dark energy, using DESI Data Release 1 (DR1) BAO data in combination with one of three different supernova compilations (PantheonPlus, Union3, and DES-SN5YR) and Planck CMB observations. Our results hint towards an evolving and emergent dark energy behaviour, with negligible presence of dark energy at z ≳ 1, at varying significance depending on data sets combined. In all these reconstructions, the cosmological constant lies outside the 95% confidence intervals for some redshift ranges. This dark energy behaviour, reconstructed using Crossing Statistics, is in agreement with results from the conventional w 0 –w a dark energy equation of state parametrization reported in the DESI Key cosmology paper. Our results add an extensive class of model-agnostic reconstructions with acceptable fits to the data, including models where cosmic acceleration slows down at low redshifts. We also report constraints on H 0 r d from our model-agnostic analysis, independent of the pre-recombination physics.
Supercritical carbon dioxide (sCO 2 ) power cycles show promising potential of higher plant efficiencies and power densities for a wide range of power generation applications such as fossil fuel power plants, nuclear power production, solar power, and geothermal power generation. sCO 2 leakage through the turbomachinery has been one of the main concerns in such applications. To offer a potential solution, we propose an elastohydrodynamic (EHD) seal that can work at elevated pressures and temperatures with low leakage and minimal wear. The EHD seal has a very simple, sleeve-like structure, wrapping on the rotor with minimal initial clearance at micron levels. Here, in this work, a proof-of-concept study for the proposed EHD seal was presented by using the simplified Reynolds equation and Lame’s formula for the fluid flow in the clearance and for seal deformation, respectively. The set of nonlinear equations was solved by using both the conventional Prediction–Correction (PC) method and modern Physics-Informed Neural Network (PINN). It was shown that the physics-informed deep learning method provided good computational efficiency in resolving the steep pressure gradient in the clearance with good accuracy. The results showed that the leakage rates increased quadratically with working pressures and reached a steady-state at high-pressure values of 15~20 MPa, where Q = 300 g/s at 20 MPa for an initial seal clearance of 255 μm. This indicates that the EHD seal could be tailored to become a potential solution to minimize the sCO 2 discharge in power plants.
It has been suggested that Fe abundances of K dwarfs derived from Fe i and Fe ii lines show considerable discrepancies, and oxygen abundances determined from high-excitation O i 7771–5 triplet lines are appreciably overestimated (the problem becoming more serious toward lower T {sub eff}), which, however, has not yet been widely confirmed. With the aim of clarifying this issue, we spectroscopically determined the atmospheric parameters of 148 G–K dwarfs (Hyades cluster stars and field stars) by assuming the classical Fe i/Fe ii ionization equilibrium as usual, and determined their oxygen abundances by applying the non-local thermal equilibrium spectrum fitting analysis to O i 7771–5 lines. It turned out that the resulting parameters did not show any significant inconsistency with those determined by other methods (for example, the mean differences in T {sub eff} and logg from the well-determined solutions of Hyades dwarfs are mostly ≲100 K and ≲0.1 dex). Likewise, the oxygen abundances of Hyades stars are around [O/H] ∼ +0.2 dex (consistent with the metallicity of this cluster) without exhibiting any systematic T {sub eff}-dependence. Accordingly, we conclude that parameters can be spectroscopically evaluated to a sufficient precision in the conventional manner (based on the Saha–Boltzmann equation for Fe i/Fe ii) and oxygen abundances can be reliably determined from the O i 7771–5 triplet for K dwarfs as far as stars of T {sub eff} ≳ 4500 K are concerned. We suspect that previously reported strongly T {sub eff}-dependent discrepancies may have stemmed mainly from overestimation of weak-line strengths and/or improper T {sub eff} scale.
We recast superfluid hydrodynamics as the hydrodynamic theory of a system with an emergent anomalous higher-form symmetry. The higher-form charge counts the winding planes of the superfluid - its constitutive relation replaces the Josephson relation of conventional superfluid hydrodynamics. This formulation puts all hydrodynamic equations on equal footing. The anomalous Ward identity can be used as an alternative starting point to prove the existence of a Goldstone boson, without reference to spontaneous symmetry breaking. This provides an alternative characterization of Landau phase transitions in terms of higher-form symmetries and their anomalies instead of how the symmetries are realized. This treatment is more general and, in particular, includes the case of BKT transitions. As an application of this formalism we construct the hydrodynamic theories of conventional (0-form) and 1-form superfluids.
The gyrokinetic (GK) field equation is a three-dimensional (3D) elliptic equation, but it is often simplified to a set of two-dimensional (2D) equations by assuming that the field does not vary along a specific direction. However, this simplification can introduce inevitable 0th-order numerical errors, as nonlinear mode coupling in toroidal geometry can produce undesirable harmonic modes that violate the assumption. In this work, we propose a novel directional finite difference method (FDM) with a local coordinate transformation to better resolve the target field of interest. The directional FDM can accurately solve 3D GK field equations without simplifications, which can overcome the limitations of conventional methods. The accuracy and efficiency of different FDMs are analyzed in great detail for a variety of geometries, from simple 2D Cartesian coordinates to realistic 3D curvilinear coordinates. The 0th-order numerical errors of simplified 2D GK equations were found to be more problematic for low-harmonic modes and low aspect ratio geometries such as spherical tokamaks. On the other hand, the directional 3D FDM can accurately resolve a much wider range of harmonic modes aligned to the direction of interest, including the low-harmonic modes. In conclusion, we demonstrate that the directional 3D FDM is a highly effective algorithm for solving the 3D GK field equations, achieving accuracy improvements of 10 to 100 times or more, particularly for low-harmonic modes in spherical tokamaks.
Here, we present a novel finite element analysis of inelastic structures containing Shape Memory Alloys (SMAs). Phenomenological constitutive models for SMAs lead to material nonlinearities, that require substantial computational effort to resolve. Finite element analysis methods, which rely on Gauss quadrature integration schemes, must solve two sets of coupled differential equations: one at the global level and the other at the local, i.e. Gauss point level. In contrast to the conventional return mapping algorithm, which solves these two sets of coupled differential equations separately using a nested Newton procedure, we propose a scheme to solve the local and global differential equations simultaneously. In the process we also derive closed-form expressions used to update the internal/constitutive state variables, and unify the popular closest-point and cutting plane methods with our formulas. Numerical testing indicates that our method allows for larger thermomechanical loading steps and provides increased computational efficiency, over the standard return mapping algorithm.