Search NASA⌕ Search

SEARCH · Search NASA

Results for “BOUNDARY VALUE”

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 OSTI2024

Fast permeability measurement for tight reservoir cores using only initial data of the one chamber pressure pulse decay test

Here, in this study, a mathematical model for fast determination of the permeabilities of tight rocks using measurements taken from the initial period of the One Chamber Pressure Pulse Decay (OC-PPD) test is presented. The model applies to measurements taken both before and after the pressure pulse front has reached the downstream end of the specimen. The analytical solutions for the pressure decay in the upstream chamber are derived based on a parabolic arc approximation of pore pressure distribution along the test specimen. This approximation allows converting the initial–boundary value problem of fluid diffusion in the specimen, governed by partial differential equations, to a system of ordinary differential equations that can be easily solved by explicit formulae. Thus, an explicit formula for the pressure decay rate is obtained, which enables inverse analysis of the initial experimental data to estimate the rock permeability. The proposed method expedites the pulse decay test as it does not require the system to reach equilibrium. The method is validated with three sets of experimental data of the OC-PPD test using helium as the diffusing fluid, for which the relative error of the permeability is found to be less than 6%. This method is particularly useful if the equilibrium time of the pulse decay test for rock specimens with permeabilities in the range of nano-Darcy takes hours or days.

early-time solution↗
DOE OSTI2024

Encoding of linear kinetic plasma problems in quantum circuits via data compression

We propose an algorithm for encoding linear kinetic plasma problems in quantum circuits. The focus is on modelling electrostatic linear waves in a one-dimensional Maxwellian electron plasma. The waves are described by the linearized Vlasov–Ampère system with a spatially localized external current that drives plasma oscillations. This system is formulated as a boundary-value problem and cast in the form of a linear vector equation Aψ = b to be solved by using the quantum signal processing algorithm. The latter requires encoding of matrix A in a quantum circuit as a sub-block of a unitary matrix. We propose how to encode A in a circuit in a compressed form and discuss how the resulting circuit scales with the problem size and the desired precision.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Constructing field-aligned coordinate systems for gyrokinetic simulations of tokamaks in X-point geometries

Structures in tokamak plasmas are elongated along the direction of the magnetic field and short in the directions perpendicular to the magnetic field. Many tokamak simulation codes take advantage of this by using a field-aligned coordinate system. However, field-aligned coordinate systems have a coordinate singularity at magnetic X-points where the poloidal magnetic field vanishes, which makes it difficult to use field-aligned coordinate systems when simulating the core and scrape-off layer simultaneously. Here, we present an algorithm for grid generation and computing geometric quantities in a standard field-aligned coordinate system that avoids the singularity and allows one to conduct two-dimensional gyrokinetic axisymmetric simulations in X-point geometries. Convergence tests of advection, boundary value problems and geometric quantities all show greater than first-order convergence even in the vicinity of the X-point. We also demonstrate the geometric consistency of our algorithm with an example simulation of the spherical tokamak for energy production, which shows machine-precision particle conservation.

fusion plasma↗
DOE OSTI2025

On the motion of compact objects in relativistic viscous fluids

We present a world-line effective field theory of compact objects moving relativistically through a viscous fluid. The theory is valid when velocity gradients are small compared to the inverse size of the object. Working within the EFT eliminates the need to solve a boundary value problem by turning all interactions between the fluid and the object into a source term in the action. We use the EFT to derive the relativistic equations of motion for a compact object immersed in a viscous fluid in a curved background, when the relative velocity of the object and the fluid is small compared to the speed of light.

astrophysical black holes↗
DOE OSTI2025

Extended Galerkin Neural Network Approximation of Singular Variational Problems with Error Control

We present extended Galerkin neural networks, a variational framework for approximating general boundary value problems (BVPs) with error control. The main contributions of this work are (1) a rigorous theory guiding the construction of new weighted least squares variational formulations suitable for use in neural network approximation of general BVPs, and (2) an “extended” feedforward network architecture which incorporates and is even capable of learning singular solution structures, thus greatly improving approximability of singular solutions. Furthermore, numerical results are presented for several problems, including steady Stokes flow around reentrant corners and in convex corners with Moffatt eddies in order to demonstrate efficacy of the method.

a posteriori error estimate↗
DOE OSTI2024

LLNL FESP Theory Highlights: October 2024

I. Novikau, I. Y. Dodin, E. A. Startsev, I. Joseph, Quantum algorithms for simulating dissipative linear and nonlinear dynamics of plasmas. Invited talk at the 66th Annual Meeting of the APS Division of Plasma Physics, Atlanta, Georgia. Novikau I., Dodin I.Y., Startsev E.A., Encoding of linear kinetic plasma problems in quantum circuits via data compression, Journal of Plasma Physics. 2024;90(4):805900401, doi:10.1017/S0022377824000795. We propose an algorithm for encoding linear kinetic plasma problems in quantum circuits. The focus is on modelling electrostatic linear waves in a one-dimensional Maxwellian electron plasma. The waves are described by the linearized Vlasov–Ampère system with a spatially localized external current that drives plasma oscillations. This system is formulated as a boundary-value problem and cast in the form of a linear vector equation to be solved by using the quantum signal processing algorithm. The latter requires encoding of a matrix in a quantum circuit as a sub-block of a unitary matrix. We propose how to encode in a circuit in a compressed form and discuss how the resulting circuit scales with the problem size and the desired precision.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Model Development and Analysis of a High-Fidelity Neutron Transport Sensor: The Quadrupole Detector Concept for Measurement of the Neutron Flux Gradient

Accurate reconstruction of the neutron flux distribution within a reactor core is essential for safe and efficient reactor operation. Traditional power shape synthesis in Light Water Reactors relies on hundreds of in-core detectors. However, this approach becomes impractical for Advanced Reactors and Microreactors due to limited space and harsh environments. To address this challenge, we propose a data-driven methodology that combines high-fidelity modeling with real-time ex-core sensor measurements, enabling the reconstruction of core power distribution while minimizing the reliance on intrusive in-core instrumentation. This project began in FY24 and achieved two initial milestones: (1) the definition of a three-year development plan for a Digital Twin framework and (2) the development of high-fidelity neutronics models of the Purdue University Reactor One (PUR-1) using both MCNP6 and OpenMC. The PUR-1 reactor, a zero-power facility, was selected due to its suitability for neutronics-focused modeling and the availability of experimental data for validation. Both models were benchmarked using neutron flux measurements obtained from irradiated gold foils, which were strategically placed within the core during a dedicated campaign in July 2024. This report marks the continuation and completion of those foundational tasks. The OpenMC model has been refined (improved geometric accuracy, expanded cross-section libraries, and refined sampling) and validated using additional experimental data. An updated sensor design—based on quadrupole configuration—was designed to measure both ex-core flux and its spatial gradient. These measurements will serve as inputs to a neural network-based reconstruction algorithm. Finally, the methodology was demonstrated on a two-dimensional test case representative of the heterogeneous material composition of the PUR-1 reactor core. A neural network implementation of the Kirchhoff-Helmholtz integral equation was employed to solve the boundary value problem using peripheral sensor measurements. The preliminary results confirm the strong potential of the proposed approach for accurate and minimally invasive neutron flux reconstruction.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Charge And Dynamic Current On Tubular Antennas For Various Drive Conditions

The mixed boundary value problem of a tubular conductor is solved using an approximate representation of its Fourier coefficients. A two term solution is derived, which represents the solution over an extremely broad range of aspect ratios. This representation is used to find the electrostatic solution and capacitance of a charged tube as well as the solution of a tube in a uniform field and its dipole moment. This second case is directly useful as a model for a monopole electric field probe. This approximation is a special case of a representation using a combination of Chebyshev and Legendre polynomials. Combining the charged tube and tube in a uniform field allows the solution of voltage driven tubular antennas. Comparisons are made with numerical solutions using piecewise sinusoidal representations of the current. The results are also generalized to the dynamic case and up to and beyond the first resonance. Simple corrections for finite gap and delta gap drives to magnetic frill drives are examined using infinite tube integral transform representations. Corrections between magnetic frill drives and coaxial drives are also given. Approximate drive corrections using conformal mapping and an effective radius are also discussed. Finally, this efficient current representation is applied to the magnetic problem involving simple tubular solenoids.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Planetary Boundary-Layer Height (PBLHT) Value-Added Product: Remote-Sensing Retrievals

The planetary boundary layer (PBL) is fundamental to numerous atmospheric processes, including aerosol mixing and transport, cloud evolution, and precipitation formation. A critical parameter in these studies is the PBL height (PBLHT). This vertical depth is essential for characterizing PBL structures in numerical simulations and serves as a primary metric for estimating flux exchanges between the Earth’s surface and the atmosphere. Radiosonde (SONDE) observations provide high-vertical-resolution measurements of temperature and moisture profiles and are widely used to estimate PBLHT (Liu and Liang 2010, Seidel et al. 2010). The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility’s PBLHT value-added product (VAP) for radiosonde measurements, known as PBLHTSONDE, applies three commonly used methods—the Heffter (1980) method, the Liu and Liang (2010) method, and the bulk Richardson number approach (Seibert et al. 2000)—to derive PBLHT. The PBLHTSONDE VAP operates routinely at ARM observatories and mobile facilities, with data available from the ARM Data Center shortly after sounding observations are collected (Sivaraman et al. 2013). However, radiosonde observations are limited by their low temporal resolution. Most stations launch soundings only twice daily, which constrains the ability to investigate and characterize the temporal evolution of the PBL using radiosonde data alone. The use of continuous remote-sensing observations provides high temporal resolution of PBLHT estimates. These observations include aerosol lidars (Dang et al. 2019, Su et al. 2020), Doppler lidar (DL; Tucker et al. 2009, Krishnamurthy et al. 2021), and water vapor and/or temperature lidars and radiometers (Turner et al. 2014). These observations provide valuable data on the PBL’s thermodynamic properties (e.g., water vapor and/or temperature lidars and radiometers), dynamic properties (e.g., DL), and distribution of tracer substances (e.g., aerosol lidars), all of which can be used to estimate PBLHT. ARM developed PBLHT estimates from the micropulse lidar (MPL; PBLHTMPL), Doppler lidar (PBLHTDL), and combined Raman lidar (RL)/atmospheric emitted radiance interferometer (AERI) thermodynamic profiles (PBLHTTHERMO). Each estimate captures different physical characteristics of the boundary layer—aerosol tracers, vertical velocity turbulence, and thermodynamic structure—and exhibits distinct strengths and limitations depending on the PBL regime and time of day. In addition, the ARM ceilometer (CEIL) provides three potential PBLHT candidates derived from the vendor's built-in algorithm. Building on these individual retrievals, ARM developed the PBLHTBEML VAP, which combines the four remote-sensing-based estimates with ancillary meteorological variables using the machine learning approach of Zhang et al. (2025) to produce a best-estimate PBLHT at 10-minute resolution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Analysis of Control Behavior in Eco-Driving Speed Optimization Using Pontryagin’s Minimum Principle

The energy efficiency of autonomous vehicles can be improved by selecting an optimized speed profile. Energy savings can be maximized by performing control optimization with knowledge of the powertrain characteristics and future driving conditions. Previous studies have shown that Pontryagin’s minimum principle (PMP) performs well in vehicle speed optimization problems. Building on the methods proposed in previous studies, the contribution of this study is to derive meaningful observations from the concepts and results of PMP to enhance the understanding of the control problem. In particular, the switching behavior of the control mode is analyzed with supportive variables, such as ξ and mv, which dictates the changes in the control modes. Additionally, the existence of the singular control is analyzed, which helps in understanding the cruise driving in the control problem. Finally, we obtain several solutions that satisfy various boundary conditions along with a map of the reachable states, and discuss the impact of cruise driving. This is helpful for designing practical control concepts for real-world applications based on this map. Previous studies have contributed significantly to this control problem; however, this study provides a better understanding of the issue and offers guidance and inspiration for future real-world applications based on these meaningful observations.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Temperature Measurements in Hypersonic Wind Tunnels via Femtosecond Coherent Anti-Stokes Raman Scattering

A femtosecond coherent anti-Stokes Raman scattering (fs CARS) instrument is developed to perform gas-phase thermometry in cold-flow hypersonic wind tunnels. Measurements are reported for Mach 8 and 14 pure-nitrogen flows. The fs CARS instrument includes a 100 fs pump/Stokes pulse and a spectrally narrow probe pulse from a second harmonic bandwidth compressor. Important experimental considerations such as limits on the pump/Stokes pulse energy are discussed. The fs CARS focusing and collimating optics are mounted on a two-axis translation stage system to scan the measurement location during a 30 second wind tunnel run. Single-laser-shot rotational CARS spectra are recorded at the laser repetition rate of 1 kHz in the wind tunnel freestream and near simple cone models. Spectral fitting is used to determine quantitative gas temperatures. Freestream temperatures at Mach 8 and 14 spanned ranges of 40–75 and 35–50 K, respectively, depending on tunnel operating conditions. Temperature variations across the central 100 mm span of the wind tunnel were quantified. Measured temperature jumps across conical bow shocks from various models varied by less than 1% from predicted values. Hypersonic boundary layer measurements were demonstrated. In conclusion, these measurements illustrate the utility and robustness of this instrument for the study of complex fluid flow phenomena in challenging ground test facilities.

Aerodynamics↗
DOE OSTI2024

Spatially distributed atmospheric boundary layer properties in Houston – A value-added observational dataset

Abstract In 2022, Houston, TX became a nexus for field campaigns aiming to further our understanding of the feedbacks between convective clouds, aerosols and atmospheric boundary layer (ABL) properties. Houston’s proximity to the Gulf of Mexico and Galveston Bay motivated the collection of spatially distributed observations to disentangle coastal and urban processes. This paper presents a value-added ABL dataset derived from observations collected by eight research teams over 46 days between 2 June - 18 September 2022. The dataset spans 14 sites distributed within a ~80-km radius around Houston. Measurements from three types of instruments are analyzed to objectively provide estimates of nine ABL parameters, both thermodynamic (potential temperature, and relative humidity profiles and thermodynamic ABL depth) and dynamic (horizontal wind speed and direction, mean vertical velocity, updraft and downdraft speed profiles, and dynamical ABL depth). Contextual information about cloud occurrence is also provided. The dataset is prepared on a uniform time-height grid of 1 h and 30 m resolution to facilitate its use as a benchmark for forthcoming numerical simulations and the fundamental study of atmospheric processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Revealing the complex chemistry of grain boundaries in K-doped BaFe 2 As 2 with atom probe tomography

Iron-based superconductors have attractive properties for high-field applications, but there is a lack of understanding of the effect of grain boundary chemistry on the in-field performance. The near atomic-scale resolution, ppm sensitivity and 3D analysis offered by atom probe tomography make it a powerful tool to investigate the nanoscale structure and chemistry of these defects in fine-grained K-doped BaFe 2 As 2 samples. A computational method to systematically extract and compare the Gibbsian interfacial excess of chemical species across grain boundaries has been explored in this work. The robustness of the method has been tested by evaluating the effects of selected variables on simulated APT datasets. The accuracy and precision of the calculated Gibbsian interfacial excess were found to be stable over a range of analysis conditions: varying grain boundary widths and detection efficiencies, spatial precisions below 1.5 nm, and bin widths between 1.2 and 1.6 nm. For the K-doped BaFe 2 As 2 samples studied, segregation of As, Ba, K and impurities of O, Na, and Sb were found at grain boundaries. The Gibbsian excess values were found to vary widely between different boundaries, showing the complexity of the grain boundary chemistry in this material. Possible links between the observed critical current density (Jc) of these samples and their nano- and micro-structure have also been investigated and discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Absolute entropy and the observer’s no-boundary state

We investigate the no-boundary proposal for closed universes with an observer. We argue that the observer’s no-boundary state is the identity operator on the physical Hilbert space, i.e., the maximum entropy state and show this explicitly in Jackiw-Teitelboim gravity. Geometrically, the no-boundary state is a bra-ket wormhole. Expectation values in the no-boundary state provide a trace for the observer’s algebra, which allows one to define von Neumann entropy for observers in different universes as the relative entropy with respect to the no-boundary state. This result is consistent with all previously discussed cases of traces for invariantly defined regions.

Cosmological models↗
DOE OSTI2024

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗
DOE OSTI2025

Investigation of tearing mode stability near ideal stability boundaries via asymptotic matching techniques

A number of improvements to the TJ toroidal tearing mode code [Fitzpatrick, Phys. Plasmas 31, 102507 (2024)] are documented. The TJ code is also successfully benchmarked against the STRIDE toroidal tearing mode code [Glasser and Koleman, Phys. Plasmas 25, 082502 (2018)]. Finally, the new capabilities of the TJ code are used to investigate the stability of tearing modes in tokamak plasmas as an ideal stability boundary, associated with either an external-kink or an internal-kink mode, is approached. All elements of the tearing stability matrix are found to tend to infinity as an ideal stability boundary is approached. Furthermore, as the stability boundary is approached, the eigenfunctions of the various tearing modes in the plasma, which are decoupled by sheared plasma rotation, are all found to morph into that of the marginally stable ideal mode. However, the growth-rates and real frequencies of the various “ideal-tearing-modes” are different from one another. Moreover, the growth-rate of the ideal-tearing-mode that reconnects magnetic flux at the rational surface that lies closest to the edge of the plasma is the one that tends to a very large value as the stability boundary is approached. A relatively simple test for ideal stability that is capable of detecting stability boundaries for external-kink and internal-kink modes, even in the presence of a very close-fitting ideal wall, is described and verified.

FOS: Physical sciences↗
DOE OSTI2025

Resistive hose modes in tokamak runaway electron beams II

Resistive hose instabilities of runaway electron (RE) beams immersed in resistive background plasmas are examined with analytic and numerical calculations. The RE beam-plasma equilibria considered are characteristic of the situation observed post-thermal quench in a tokamak disruption. An analytic linear dispersion relation is presented for the case of a uniform RE current density profile with a sharp boundary in cylindrical geometry. Initial value linear calculations for a more general profile in toroidal geometry find that reducing aspect ratio increases the resistive hose mode growth rates with fixed safety factor profile. Nonlinear calculations in cylindrical and toroidal geometry find that the resistive hose instability-driven fluctuations relax the gradient of the current density profile. In toroidal geometry, changes to the magnetic topology are observed as a result of the resistive hose activity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗