Search NASA⌕ Search

SEARCH · Search NASA

Results for “analytical and numerical solutions”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Verification of numerical models for seismic fluid-structure interaction analysis of internal components in liquid-filled advanced reactors

Earthquake shaking of a liquid-filled advanced reactor induces fluid-structure interaction (FSI) between the reactor vessel (tank), its internal components, and the contained liquid. Seismic design, qualification, and risk assessment of such reactors must consider fluid-structure responses, their geometries and support conditions, and three-directional seismic inputs, all of which require numerical simulations. Numerical models used for analysis of safety-related nuclear equipment must be verified and validated. Here, this paper verifies numerical models of submerged components using analytical solutions. Seismic FSI analysis of the numerical models is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Prior analytical solutions are reworked in this paper, and calculation errors are identified and corrected. The solutions address frequencies of two concentric cylindrical pipes filled with liquid: the inner pipe is a submerged component, and the outer pipe contains the liquid. A unitless frequency coefficient is defined and presented for pipes of different materials, filled with different liquids, and with a range of dimensions for application to advanced reactors. The numerical models are verified here by comparing the lateral frequencies of submerged components with those calculated using the corrected analytical solutions. Recommendations for verification of numerical models of internal components in advanced reactors for seismic FSI analysis are provided. Although the unitless frequency coefficients and verification procedures are developed for application to advanced nuclear reactors, they are broadly applicable to FSI analysis of submerged components in liquid-filled vessels such as storage tanks, boilers, and steam generators.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Relativistic levels of mesic atoms

We revisit the derivation of the covariant two-body scalar-fermion equation with a Coulomb interaction, presented in a previous paper. We show that it can be given the formal aspect of a Dirac equation, but for the fact that the eigenvalue is also contained in one of the coefficients and thus it is not linearly included. The discussion of the boundary value problem is therefore different, although some properties of the Dirac equation can be recovered in an approximation bringing back to the concept of reduced mass. We discuss a mixed analytic-numerical method of solution which allows to obtain very accurate results and we calculate the lowest levels, states and QED corrections for pionic and kaonic atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

HARD: A performance portable radiation hydrodynamics code based on FleCSI framework

Hydrodynamics And Radiation Diffusion (HARD) is an open-source application for high-performance simulations of compressible hydrodynamics with radiation-diffusion coupling. Built on the FleCSI (Bergen et al., 2021 [1]) (Flexible Computational Science Infrastructure) framework, HARD expresses its computational units as tasks whose execution can be orchestrated by multiple back-end runtimes, including Legion (Bauer et al., 2012 [2]), MPI (Forum, 1994 [3]), and HPX (Kaiser et al., 2020 [4]). Node-level parallelism is handled through Kokkos (Edwards et al., 2014 [5]), providing a single-source, portable code base that runs efficiently on laptops, small homogeneous clusters, and the largest heterogeneous supercomputers currently available. To ensure scientific reliability, HARD includes a regression test suite that automatically reproduces canonical verification problems such as the Sod and LeBlanc shock tubes, and the Sedov blast wave, comparing numerical solutions against known analytical results. The project is distributed under an OSI-approved license, hosted on GitHub, and accompanied by reproducible build scripts and continuous integration workflows. This combination of performance portability, verification infrastructure, and community-focused development makes HARD a sustainable platform for advancing radiation hydrodynamics research across multiple domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Sachdev-Ye-Kitaev models and beyond: Window into non-Fermi liquids

In this work is a review of the Sachdev-Ye-Kitaev (SYK) model of compressible quantum many-body systems without quasiparticle excitations, and its connections to various theoretical studies of non-Fermi liquids in condensed matter physics. The review is placed in the context of numerous experimental observations on correlated electron materials. Strong correlations in metals are often associated with their proximity to a Mott transition to an insulator created by the local Coulomb repulsion between the electrons. The phase diagrams of a number of models of such a local electronic correlation are explored, employing a dynamical mean-field theory in the presence of random spin exchange interactions. Numerical analyses and analytical solutions, using renormalization group methods and expansions in large spin degeneracy, lead to critical regions that display SYK physics. The models studied include the single-band Hubbard model, the t–J model, and the two-band Kondo-Heisenberg model in the presence of random spin exchange interactions. Also examined are non-Fermi liquids obtained by considering each SYK model with random four-fermion interactions to be a multiorbital atom, with the SYK atoms arranged in an infinite lattice. Connections are made to theories of sharp Fermi surfaces without any low-energy quasiparticles in the absence of spatial disorder, obtained by coupling a Fermi liquid to a gapless boson; a systematic large-N theory of such a critical Fermi surface, with SYK characteristics, is obtained by averaging over an ensemble of theories with random boson-fermion couplings. Finally, an overview of the links between the SYK model and quantum gravity is presented, and the review ends with an outlook on open questions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Heart Shape to Fracture Distance: Characterizing Hydraulic Fracture Propagation before Hits

Estimating the distance from the hydraulic fracture tip to the monitor well can be useful for fracture characterization, well spacing optimization, and preventing parent-child well interference. A heart-shaped signal is referred to as the extensional precursor of a fracture hit recorded by crosswell strain measurements and can serve as a vital tool for such estimation. This study incorporates the 3D displacement discontinuity method (DDM) to understand the impact of fracture geometry and monitor well offset on the heart-shaped signal’s characteristics. Results from numerical simulation and analytical solutions reveal a strong linear correlation between the spatial extent of the heart-shaped signal and the fracture tip distance. This relationship was further developed to predict tip distance using field data from the Hydraulic Fracture Test Site 2 (HFTS2). A reasonable approximation result from field data further validates the methodology. In addition, it is worth noting that the estimation accuracy depends on the ratio between fracture dimension and tip distance. The findings of this study offer a novel approach for real-time monitoring and characterizing hydraulic fracture propagation, which can be further used for well spacing optimization in unconventional and enhanced geothermal system reservoir development, as well as caprock integrity monitoring for carbon sequestration projects.

58 GEOSCIENCES↗
DOE OSTI2024

Characterizing Hydraulic Fracture Propagation Before Fracture Hits

Estimating the distance from hydraulic fracture tip to monitor well can be very useful for fracture characterization, well spacing optimization, and preventing parent-child well interference. Heart-shape signal is referred to as the extensional precursor of fracture hit recorded by cross-well strain measurements and can be served as a vital tool to make such estimation. This study incorporates the 3D Displacement Discontinuity Method to understand the impact of fracture geometry and monitor well offset on the heart-shape signal’s characteristics. Results from numerical simulation and analytical solutions reveal a strong linear correlation between the spatial extent of heart shape signal and fracture tip distance. This relationship was further developed to predict tip distance using field data from the Hydraulic Fracture Test Site 2. A reasonable approximation result from field data further validates the methodology. In addition, it is worth noting that the estimation accuracy depends on the ratio between fracture dimension and tip distance. The findings of this study offer a novel approach for real-time monitoring and characterizing hydraulic fracture propagation, which can be further used for well spacing optimization in unconventional and Enhanced Geothermal System reservoir development, as well as cap rock integrity monitoring for carbon sequestration projects.

Jin, Ge↗
DOE OSTI2025

Testing horizontal momentum boundary conditions in MPAS-Ocean

A rising concern within E3SM’s MPAS-Ocean is the implementation of boundary conditions. Currently, MPAS-Ocean utilizes a no-slip boundary condition, meaning that the velocity of fluid at the boundaries is zero; however, the model should should satisfy both no-slip and free-slip boundary conditions given the correct parameters. This study focuses on testing these boundary conditions in depth, determining the accuracy of the current implementation. To accomplish this, we focus on two test cases: the barotropic gyre and barotropic channel. By manipulating parameters and comparing the numerical to corresponding analytical solutions (when available), we show that MPAS-Ocean’s ability to satisfy both boundary conditions is promising.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Theoretical modeling of a bottom-raised oscillating surge wave energy converter structural loadings and power performances

Here, this study presents theoretical formulations to evaluate the fundamental parameters and performance characteristics of a bottom-raised oscillating surge wave energy converter (OSWEC) device. Employing a flat plate assumption and potential flow formulation in elliptical coordinates, closed-form equations for the added mass, radiation damping, and excitation forces/torques in the relevant pitch-pitch and surge-pitch directions of motion are developed and used to calculate the system's response amplitude operator and the forces and moments acting on the foundation. The model is benchmarked against numerical simulations using WAMIT and WEC-Sim, showcasing excellent agreement. The sensitivity of plate thickness on the analytical hydrodynamic solutions is investigated over several thickness-to-width ratios ranging from 1:80 to 1:10. The results show that as the thickness of the benchmark OSWEC increases, the deviation of the analytical hydrodynamic coefficients from the numerical solutions grows from 3% to 25%. Differences in the excitation forces and torques, however, are contained within 12%. While the flat plate assumption is a limitation of the proposed analytical model, the error is within a reasonable margin for use in the design space exploration phase before a higher-fidelity (and thus more computationally expensive) model is employed. A parametric study demonstrates the ability of the analytical model to quickly sweep over a domain of OSWEC dimensions, illustrating the analytical model's utility in the early phases of design.

13 HYDRO ENERGY↗
DOE OSTI2025

Acausality-driven instabilities in relativistic viscous hydrodynamics

We investigate non-linear instabilities stemming from superluminal propagation of information in Israel-Stewart-like models of relativistic viscous fluid dynamics. In relativity, the characteristic speed of propagation of information, $w$, and the speed of the fluid, $v$, allow us to differentiate between regimes of the hydrodynamic equations that are acausal but stable ($w>1$), unstable ($v^{2} w^{2} \geq 1$), and covariantly ill-posed ($w^{2} \leq 0$). As an analytical benchmark, we present a new solution that illustrates these distinct regimes. We compare this analytical solution to the result of a numerical relativistic viscous fluid dynamics solver, and confirm that the analytical result can be recovered numerically in the stable regime, whether causal or acausal. The onset of numerical instabilities is further found to occur in the regime predicted by the analytical solution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Should We Conserve Entropy or Energy when Computing CAPE with Mixed-Phase Precipitation Physics?

Abstract The rapidly increasing resolution of global atmospheric reanalysis and climate model datasets necessitates finding methods for computing convective available potential energy (CAPE) both efficiently and accurately. To this end, this article compares two common methods for computing CAPE which conserve either energy or entropy. Inaccuracies in these computations arise from both physical and numerical errors. For instance, computing CAPE with entropy conserved results in physical errors from nonequilibrium phase transitions but minimizes numerical errors because solutions are analytic at each height. In contrast, computing CAPE with energy conserved avoids these physical errors, but accumulates numerical errors that are grid-resolution-dependent because the numerical integration of a differential equation is required. Analysis of CAPE computed with large databases of soundings from the tropical Amazon and midlatitude storm environments shows that physical errors from the entropy method are typically 1%–3% as large as CAPE, which is comparable to the numerical errors from conserving energy with grid spacing of 25 and 250 m using explicit first-order and second-order integration schemes, respectively. Errors in entropy-based CAPE calculations are also insensitive to vertical grid spacing, in contrast to energy-based calculations whose error strongly scales with the grid spacing. It is shown that entropy-based methods are advantageous when intercomparing datasets with differing vertical resolution because they produce accurate and reasonably fast results that are insensitive to grid resolution, whereas a second-order energy-based method is advantageous when analyzing data with a consistent vertical resolution because of its superior computational efficiency. Significance Statement Convective available potential energy (CAPE) is a measure of instability in the atmosphere that helps forecasters and researchers understand when and where thunderstorms will form. The purpose of this article is to identify the most efficient and accurate methods for computing CAPE. Two methods are considered here, one that relates to the entropy (a measure of thermodynamic disorder) of an air parcel and one that relates to the energy of an air parcel. Results indicate that the entropy method is most accurate and insensitive to the resolution of the data used for the calculation (which can vary considerably), whereas the energy method uses the least computation time.

Peters, John M.↗
DOE OSTI2021

Analytical solutions of the Arrhenius-Semenov problem for constant volume burn

Analytical solutions to the Semenov thermal ignition problem for constant volume burn governed by Arrhenius reaction kinetics are derived. Specifically, an approximate analytical solution technique for the Arrhenius-Semenov differential equation is derived for reaction orders n ϵ R> 0 and exact solutions are also constructed for reaction orders n ϵ N : n ≤ 3. The approximation technique relies on expansion of the respective nondominant terms in the differential equation at the lower and upper bounds of the reaction progress variable in order to create a pair of integrable series. The two integrated series are then connected to create a single continuous analytical solution. Excellent agreement is observed between the analytical approximation and solutions obtained numerically. The presented approximation constitutes a simple and robust strategy for solving the Arrhenius-Semenov problem analytically.

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

Interfacial electroneutrality controls transport of asymmetric salts through charge-patterned mosaic membranes

Membranes that selectively enhance target solute permeation while rejecting competing species are essential for precision separations. This study introduces charge-patterned mosaic membranes (CMMs) that selectively transport divalent asymmetric salts by leveraging a net-neutral membrane–solution interface. This mechanism, dictated by the charge ratio of positive and negative domains on the membrane surface and the balance of cations and anions in the salt, is supported by analytical, numerical, and experimental results. Analytical solutions identified cationic domain coverages (f + ) of 33%, 50%, and 66% as optimal for the selective transport of +2:−1 salts, +1:−1 salts, and +1:−2 salts, respectively, under conditions where the pattern size (L) is significantly larger than the Debye length. Numerical simulations and experiments using CMMs with alternating charged-stripes inkjet-printed onto nanostructure copolymer substrates confirmed these findings. By varying stripe widths to control f + , pressure-driven filtration experiments demonstrated selective enrichment of MgCl 2 and K 2 SO 4 at the predicted f + values, with deviations from these values leading to salt rejection. These results highlight the pivotal role of a net-neutral interface in enabling asymmetric salt enrichment. This study positions CMMs as a versatile platform for tuning ion selectivity, addressing challenges in resource recovery, water treatment, and precision separations.

additive manufacturing↗
DOE OSTI2024

Electric fields near undulating dielectric membranes

Dielectric interfaces are crucial to the behavior of charged membranes, from graphene to synthetic and biological lipid bilayers. Understanding electrolyte behavior near these interfaces remains a challenge, especially in the case of rough dielectric surfaces. A lack of analytical solutions consigns this problem to numerical treatments. Here, we report an analytic method for determining electrostatic potentials near curved dielectric membranes in a two-dimensional periodic “slab” geometry using a periodic summation of Green’s functions. This method is amenable to simulating arbitrary groups of charges near surfaces with two-dimensional deformations. We concentrate on one-dimensional undulations. We show that increasing membrane undulation increases the asymmetry of interfacial charge distributions due to preferential ionic repulsion from troughs. In the limit of thick membranes, we recover results mimicking those for electrolytes near a single interface. Our work demonstrates that rough surfaces generate charge patterns in electrolytes of charged molecules or mixed-valence ions.

Chemistry↗
DOE OSTI2025

Consistent solutions of the radiation diffusion equation in spherical and cylindrical geometries

We have extended the radiation diffusion model of Hammer and Rosen [Phys. Plasmas 10, 1829 (2003)] to diverging spherical and cylindrical geometries. The effect of curvilinear geometry on the supersonic, expanding wavefront increases as the internal radius of a spherical or cylindrical shell approaches zero. Small spherical geometries are important for modeling systems at the size scale of inertial confinement fusion capsules, at these scales existing quasi-analytic models for planar geometry significantly disagree with the results of simulation. With this method, the benefits of rapid iteration can be applied to common spherical systems at much smaller length scales. We present comparisons between numerical diffusion solutions and the analytic model to give ranges of applicability for the model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Verification and Demonstration of One-Dimensional Freezing Model in SAM for Salt-Cooled Reactor Analysis Applications

This work presented the development and implementation of the one-dimensional freezing model in system analysis code, SAM, as well as code verification, and code demonstration during a postulated overcooling transient, for fluoride salt-cooled high-temperature reactor (FHR) system and safety analysis applications. The paper at first summarized the freezing model, finite element numerical method, and special numerical treatment for handling phase appearance/disappearance. Analytical solutions were derived for two cases (with and without solid walls) for code verifications purpose. As expected, numerical results predicted by the SAM code agreed very well with the analytical solution. A code demonstration was then performed on a postulated protected overcooling event transient of a generic reference PB-FHR design. The code was found to successfully predict salt freezing during such a postulated event. However, due to lack of salt freezing testing data, code validation has not been performed in this work, which will be pursued in later studies when such data becomes available.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A correction procedure for self-induced velocity of a finite-sized particle in two-way coupled Euler–Lagrange simulations

The importance of incorporating a correction to undo the self-induced perturbation velocity of a particle, when its size becomes comparable to the Eulerian grid, in a two-way coupled Euler–Lagrange (EL) simulation is now well appreciated. The present work improves upon the prior correction procedures in a few important ways. First, the past correction procedures have been scalar-based with the assumption that the quasi-steady force is the source of self-induced velocity perturbation. Here we generalize to a vector correction procedure and thereby the directions of feedback force and relative velocity can be different. This allows the correction procedure to be used even in the presence of added-mass, history, and lift forces. Second, the effect of a nearby wall has been systematically included in the correction procedure. The correction procedure depends on fundamental Oseen solutions of streamwise and transverse regularized feedback forces. We present a Fourier transform-based analytical approach to obtaining these regularized Oseen solutions. We also present a step-by-step numerical procedure for obtaining the Oseen solutions in any EL code. With the analytical or numerical Oseen functions, the correction procedure can be easily implemented in any EL code. Iterations are required in solving the implicit correction equations and it is demonstrated that the correction procedure converges rapidly within three or four iterations. In conclusion, a simple empirical approach is also presented to account for unsteady effects in the correction procedure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Viscous Gubser flow with conserved charges to benchmark fluid simulations

We present semi-analytical solutions for the evolution of both the temperature and chemical potentials for viscous Gubser flow with conserved charges. Such a solution can be especially useful in testing numerical codes intended to simulate relativistic fluids with large chemical potentials. The freeze-out hypersurface profiles for constant energy density are calculated, along with the corresponding normal vectors, and presented as a new unit test for numerical codes. We also compare the influence of the equation of state on the semi-analytical solutions. We benchmark the newly developed smoothed particle hydrodynamics code ccake that includes both shear viscosity and three conserved charges. Here, the numerical solutions are in excellent agreement with the semi-analytical solution and are also able to accurately reproduce the hypersurface at freeze-out.

Hydrodynamic models↗