Search NASA⌕ Search

SEARCH · Search NASA

Results for “Poisson equation”

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

Ion Transport in Charged Membranes: Linking Electric-Field-Driven Mechanisms to Pore Size via Perturbation Analysis

Ion-exchange membranes are a critical component in electrochemical systems. Nevertheless, the understanding and modeling of ion transport within these porous structures have been limited by particular complexity reductions, either ignoring the dimensionality of their porous network architecture or imposing geometric assumptions (i.e., overlapping double layers). Before addressing this morphology-transport gap, a framework that relates the driving forces of transport to the geometry of a single pore is required. In this work, our modeling domain consists of a two-dimensional single pore with charged walls, connecting two identical electrolyte reservoirs. Using the Poisson-Nernst-Planck equations and regular perturbation theory, we decouple the electric fields and analyze the driving forces of ion transport, specifically electromigration and induced electroosmosis within the pore. These processes are described as analytical functions of the interaction aspect ratio, ?, defined as the ratio of the pore radius to the Debye length. Using this parameter, our study (i) describes the interplay between electromigrative and electroosmotic mechanisms that set ionic conductivity, (ii) identifies a dimensionless group of intrinsic electrolyte properties that indicates the predominant driving force, and (iii) provides a qualitative, confinement-dependent perspective on selectivity in ion-conducting membranes.

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

Large deviations of ionic currents in dilute electrolytes

Here, we evaluate the exponentially rare fluctuations of the ionic current for a dilute electrolyte by means of macroscopic fluctuation theory. We consider the fluctuating hydrodynamics of a fluid electrolyte described by a stochastic Poisson–Nernst–Planck equation. We derive the Euler–Lagrange equations that dictate the optimal concentration profiles of ions conditioned on exhibiting a given current, whose form determines the likelihood of that current in the long-time limit. For a symmetric electrolyte under small applied voltages, number density fluctuations are small, and ionic current fluctuations are Gaussian with a variance determined by the Nernst–Einstein conductivity. Under large applied potentials, the ionic current distribution is generically non-Gaussian. Its structure is constrained thermodynamically by Gallavotti–Cohen symmetry and the thermodynamic uncertainty principle.

Farhadi, Jafar [University of California, Berkeley↗
DOE OSTI2025

Cosmological perturbation theory for large scale structure in phase space

We develop a framework for Large Scale Structure (LSS) perturbation theory, that solves the Vlasov-Poisson system of equations for the distribution function in full phase space. This approach relaxes the usual apriori assumption of negligible velocity dispersion underlying the Standard Perturbation Theory (SPT). We apply the new method to rederive the usual SPT kernels up to third order in the perturbative expansion. We also show that a counterterm, identical to the one introduced by standard Effective Field Theory (EFT) methods, naturally arises within our framework. We finish by making a precise connection to EFT techniques, which reveals the necessity of the EFTofLSS to self-consistently model the long-wavelength fluid, and illustrates the importance of having theoretical control over short distance fluctuations.

Cosmological perturbation theory in GR and beyond↗
DOE OSTI2025

Unraveling membrane electrode assembly design for electrochemical conversion of carbon dioxide to formate/formic acid

This work presents a one-dimensional continuum modeling approach to investigate various cell architectures used for electrochemical conversion of CO 2 to formate/formic acid. Ion transport is simulated by a system of generalized modified Poisson–Nernst–Planck (GMPNP) equations that reflect the reactive transport phenomena including steric effects as the electrolyte solutions become concentrated. In the cathode catalyst layer, ionic current contributions from both the supporting electrolyte and solid-state ionomer are considered. Voltage and CO 2 utilization breakdowns are utilized to deconvolute the impacts of the cell architecture. The origins of (bi)carbonate formation in the cathode are explored, as the subsequent decrease in CO 2 availability is a key reason for low faradaic efficiencies to formate/formic acid. In addition, the role of a supporting electrolyte (KOH) is investigated to understand its tradeoffs: while the K + ions can improve both conductivity and electrochemically active surface area in the cathode, the presence of OH − ions raises the pH and leads to deleterious formation of (bi)carbonates. To this end, we also present parametric studies on the concentration and flow rate of supplied KOH to the cell, to establish a path towards eliminating the need for a supporting electrolyte.

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

Proton diffusion and hydrogen/deuterium exchange in amorphous solid water at temperatures from 114 to 134 K

The reaction coefficient for hydrogen/deuterium (H/D) exchange and the diffusion of hydrated excess protons within amorphous solid water (ASW) are characterized as a function of temperature. For these experiments, water films are deposited on a Pt(111) substrate at 108 K, and reactions with pre-adsorbed hydrogen atoms produce hydrated protons. Upon heating, protons diffuse within the water, and H/D exchange occurs when they encounter D2O probe molecules deposited in the films. The time-dependent concentration of D2O is monitored with infrared spectroscopy, and it indicates the protons diffusion from the substrate and establish an equilibrium distribution prior to significant H/D exchange for temperatures 114 K ≤T≤ 134 K. By controlling the distance between the D2O molecules and the substrate, we probe the distribution of protons within the film. It decays as x−2 for the examined range of x (12–52 nm) due to the electric field that develops between the diffusing protons and their image charges in the metal substrate. This agrees with the theoretical distance scaling for the equilibrated proton concentration in a dielectric near a metal boundary. From the proton concentration and the measured D2O decay rate, a lower bound for the proton diffusion coefficient ranging from 10−20 m2/s at 114 K to 10−18 m2/s at 134 K is estimated. The diffusion coefficient has an activation energy of 0.40 eV, which is comparable to energies reported for molecular translations and rotations of H2O, suggesting they may play a critical role in the proton diffusion mechanism within ASW.

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

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

Kinetic plasma simulations solve the Vlasov-Poisson or Vlasov-Maxwell equations to evolve scalar-variable distribution functions in position-velocity phase space and vector-variable electromagnetic fields in configuration space. The immense computational cost of evolving high-dimensional variables, and their large number of degrees of freedom, often limits the utility of continuum kinetic simulations and presents a challenge when it comes to accurately simulating real-world physical phenomena. To address this challenge, we present techniques that accelerate and minimize the computational work required for a scalable Vlasov-Poisson solver. We show theoretical hardware compute and communication bounds for solving a fourth-order finite-volume Vlasov-Poisson system. These bounds are then used to inform and evaluate the design of performance portable algorithms for a multiple graphics processing unit (GPU) accelerated version of the Vlasov-Poisson solver VCK-CPU [1]. We demonstrate that the multi-GPU Vlasov solver implementation, VCK-GPU, simultaneously minimizes required inter-process data transfer while also being bounded by the machine network performance limits. This results in an overall strong scaling speedup per timestep of up to 40x in three-dimensional phase space (one position, two velocity coordinates) and 54x in four dimensional phase space (two position, two velocity coordinates) and a 341x increase in simulation throughput of the GPU accelerated code over the existing CPU code. The GPU code is also able to weak scale up to 256 compute nodes and 1024 GPUs. In conclusion, we demonstrate that the improved compute performance enables exploring configurations which were previously computationally infeasible, including resolving fine-scale distribution function filamentation and multi-species dynamics with realistic electron-proton mass ratios.

Continuum kinetics↗
DOE OSTI2026

Vidyut3d: A GPU accelerated fluid solver for non-equilibrium plasmas on adaptive grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure three-electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate ~ 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Vidyut3d: A Gpu Accelerated Fluid Solver for Non-Equilibrium Plasmas on Adaptive Grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure twin electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate approximately 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

Sitaraman, Hariswaran↗
DOE OSTI2024

Development of a Performance Portable Non-Equilibrium Plasma Fluid Solver on Adaptive Grids

This presentation will describe the numerical techniques, programming paradigms, verification, and performance of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures. Our plasma fluid model solves the conservation equations for self-consistent electrostatic Poisson, electron and heavy species transport, and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive mesh management library, AMReX (Zhang et al., JOSS, 4 (37) 1370, 2019), and can be built and run on widely available vendor specific GPU architectures (NVIDIA/AMD/Intel). We utilize a non-subcycled second order semi-implicit time-stepping method where all adaptive mesh refinement (AMR) levels are advanced with the same time step. The composite multi-level multigrid solver from within AMReX is used for each of the governing equations that are cast into a Helmholtz equation form. We have also developed a python based chemical mechanism parser framework that uses a similar format as CANTERA (Goodwin et al., Zenodo, 2018) yaml files as input. Our custom parser reads the yaml file and provides C++ files with transport and production rate functions that can be executed on both host (CPU) and device (GPU). We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on low-pressure capacitive and high-pressure streamer discharges. Our initial performance studies indicate 10X speed-up using 20 NVIDIA GPUs versus 200 CPUs for an atmospheric streamer discharge problem solved on a 512 x 1024 x 512 grid.

graphics processing units↗
DOE OSTI2025

CLPNets: Coupled Lie–Poisson neural networks for multi-part Hamiltonian systems with symmetries

To accurately compute data-based prediction of Hamiltonian systems, it is essential to utilize methods that preserve the structure of the equations over time. We consider a particularly challenging case of systems with interacting parts that do not reduce to pure momentum evolution. Such systems are essential in scientific computations, such as discretization of a continuum elastic rod, which can be viewed as the group of rotations and translations $SE(3)$. The evolution involves not only the momenta but also the relative positions and orientations of the particles. The presence of Lie group-valued elements, such as relative positions and orientations, poses a problem for applying previously derived methods for data-based computing. We develop a novel method of data-based computation and complete phase space learning of such systems. We follow the original framework of SympNets (Jin et al., 2020) and LPNets (Eldred et al., 2024), building the neural network from phase space mappings that preserve the Lie–Poisson structure. We derive a novel system of mappings that are built into neural networks describing the evolution of such systems. We call such networks Coupled Lie–Poisson Neural Networks, or CLPNets. We consider increasingly complex examples for the applications of CLPNets, starting with the rotation of two rigid bodies about a common axis, progressing to the free rotation of two rigid bodies, and finally to the evolution of two connected and interacting $SE(3)$ components, describing the discretization of an elastic rod into two elements. Our method preserves all Casimir invariants to machine precision, preserves energy to high accuracy, and shows good resistance to the curse of dimensionality, requiring only a few thousand data points for all cases studied (three to eighteen dimensions). Additionally, the method is highly economical in memory requirements, requiring only about 200 parameters for the most complex case considered.

Data-based modeling↗
DOE OSTI2026

A Statistician’s Overview of Physics-Informed Neural Networks for Spatio-Temporal Data

The recent success of deep neural network models with physical constraints (so-called, Physics-Informed Neural Networks, PINNs) has led to renewed interest in the incorporation of mechanistic information in predictive models. Statisticians and others have long been interested in this problem, which has led to several practical and innovative solutions dating back decades. In this overview, we focus on the problem of data-driven prediction and inference of dynamic spatio-temporal processes that include mechanistic information, such as would be available from partial differential equations, with a strong focus on the quantification of uncertainty associated with data, process, and parameters. Here, we give a brief review of several paradigms and focus our attention on Bayesian implementations given they naturally accommodate uncertainty quantification. We then show that it is straight-forward to include the Bayesian PINN (B-PINN) within the Bayesian hierarchical model (BHM) framework that has long been considered for modeling dynamic spatio-temporal processes. Such a BHM-PINN is illustrated via a simulation study in which a latent nonlinear Burgers’ equation PDE governs the dynamics of Poisson distributed spatio-temporal data. Supplementary materials for this article are available online, including a standardized description of the materials available for reproducing the work.

Bayesian↗
DOE OSTI2025

Integrable higher-spin deformations of sigma models from auxiliary fields

We construct a new infinite family of integrable deformations of the principal chiral model (PCM) parametrized by an interaction function of several variables, which extends the formalism of [C. Ferko and L. Smith, An infinite family of integrable sigma models using auxiliary fields, .] and includes deformations of the PCM by functions of both the stress tensor and higher-spin conserved currents. We show in detail that every model in this class admits a Lax representation for its equations of motion, and that the Poisson bracket of the Lax connection takes the Maillet form, establishing the existence of an infinite set of Poisson-commuting conserved charges. We argue that the non-Abelian T-dual of any model in this family is classically integrable, and that T-duality “commutes” with a general deformation in this class, in a sense which we make precise. Finally, we demonstrate that these higher-spin auxiliary field deformations can be extended to accommodate the addition of a Wess-Zumino term, and we exhibit the Lax connection in this case. Published by the American Physical Society 2025

Bielli, Daniele (ORCID:0009000640034768)↗
DOE OSTI2025

High-pressure elasticity and equation of state of the fluoroelastomer Viton® A-500

Viton® A is a semi-crystalline copolymer of polyvinylidene fluoride and hexafluoropropylene used in various engineering applications due to its mechanical properties and chemical inertness. In situ ultrasonic spectroscopy and x-ray radiography measurements were performed in a Paris–Edinburgh press to measure the pressure dependence of the transverse and longitudinal acoustic velocity of the fluoroelastomer A-500 from 2.7 to 5.7 GPa at 296 K. In addition, we performed high-pressure Brillouin scattering measurements to obtain acoustic velocities from ambient pressure to 5.7 GPa to supplement the ultrasonic measurements, especially at low pressures. The acoustic velocities were then used to calculate a pressure–volume (P–V) equation of state, the bulk and shear moduli, and the Poisson's ratio. These quantities are compared with the reported pressure-dependent properties of related polymers over this range of pressures.

Acoustical properties↗
DOE OSTI2024

Physics-Based Machine Learning Methods for U-235 Forensics Signatures

Signatures of low-intensity U-235 sources have been recently studied by utilizing a variety of machine learning (ML) classifiers using features derived from gamma spectral measurements collectedunder structured campaigns. Several ML classifiers, such as ensemble of tress and classification trees, revealed misleadingly-optimistic training error due to over-fitting, and furthermore,their performance is not directly relatable to the physical properties due to their data-driven, opaque designs. We present a regression-based ML method that first estimates the inverse distanceto the source and then utilizes a threshold to infer its presence, by representing the background as a source located at an infinite distance. For the inverse distance estimation, we study the ensembleof trees and Gaussian process regression methods, and a hyper parameter auto-tuning and selection method that employs five regression estimators. These methods avoid the over-fittingobserved in several ML classifiers, while providing the classification error nearly comparable to them based on independent test data. Their error is directly related to estimates of the inversephysical distance to source, and the precision of error determines the seperability property that determines the false alarm and missed detection rates. The property of monotonic decrease of thesource strength with increasing detector distance combined with Poisson distribution of measurements is utilized to analytically validate these methods by deriving the generalization equations ofunderlying regression methods.

Rao, Nageswara↗
DOE OSTI2024

Physics-Based Machine Learning Methods for U-235 Forensics Signatures

Signatures of low-intensity U-235 sources have been recently studied by utilizing a variety of machine learning (ML) classifiers using features derived from gamma spectral measurements collected under structured campaigns. Several ML classifiers, such as ensemble of tress and classification trees, revealed misleadingly-optimistic training error due to over-fitting, and furthermore, their performance is not directly relatable to the physical properties due to their data-driven, opaque designs. We present a regression-based ML method that first estimates the inverse distance to the source and then utilizes a threshold to infer its presence, by representing the background as a source located at an infinite distance. For the inverse distance estimation, we study the ensemble of trees and Gaussian process regression methods, and a hyper parameter auto-tuning and selection method that employs five regression estimators. These methods avoid the over-fitting observed in several ML classifiers, while providing the classification error nearly comparable to them based on independent test data. Their error is directly related to estimates of the inverse physical distance to source, and the precision of error determines the seperability property that determines the false alarm and missed detection rates. The property of monotonic decrease of the source strength with increasing detector distance combined with Poisson distribution of measurements is utilized to analytically validate these methods by deriving the generalization equations of underlying regression methods.

Rao, Nageswara↗
DOE OSTI2025

Sound Velocities in Vanadium Reveal Complex Elastic Behavior at High Pressures

Compressional (VP) and shear (VS) wave velocities of polycrystalline vanadium were measured simultaneously up to 11.5 GPa at room temperature using ultrasonic interferometry in a multi-anvil press. Complex softening behavior in VS and resulting shear moduli are discovered, possibly revealing a precursor to the reported phase transition within 30–60 GPa. The current data enables a comprehensive assessment of the elastic and mechanical properties of vanadium at high pressures, including bulk and shear moduli, Young’s modulus, Poisson’s ratio, and Pugh’s ratio. Through fitting to the 3rd-order finite strain equations, the elastic moduli and their pressure derivatives were determined to be K S0 = 151 (2) GPa, G 0 = 46.9 (8) GPa, K$^{’}_{S0}$ = 3.47 (5), and G$^{’}_{0}$ = 0.62 (1). These experimental results allow us to compare with and benchmark the existing Steinberg–Guinan models for extrapolations to extreme pressure and temperature conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Analysis of the weighted shifted boundary method for the Poisson and Stokes problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. Accuracy is maintained by properly shifting the location and values of the boundary conditions. This avoids integration over cut cells and the associated implementation issues. Recently, the Weighted SBM (WSBM) was proposed for the Navier-Stokes equations with free surfaces and the Stokes flow with moving boundaries. The attribute “weighted” in the name WSBM stems from the fact that its variational form is weighted with the elemental volume fraction of active fluid. The motivation for the development of the WSBM was the preservation of the volume of active fluid to a higher degree of accuracy, which in turn resulted in improved stability and robustness characteristics in moving-boundary, time-dependent simulations. In this article, we present the numerical analysis of the WSBM formulations for the Poisson and Stokes problems. We give mathematical conditions under which the bilinear forms defining the discrete variational formulations are uniformly coercive (Poisson problem) or inf-sup stable (Stokes problem). By these results, stability and optimal convergence is proven in the natural norm; L2-error estimates can also be derived.

Approximate domain boundaries↗
DOE OSTI2025

Non-linear dynamics of jet quenching

We develop a comprehensive analytic framework for jet quenching in QCD media, based on a medium-induced parton cascade sourced by collinear virtual splittings. We show that the energy flow out of the jet cone, driven by turbulent gluon cascades, is governed by a non-linear rate equation that resums gluon splittings at arbitrary angles and is enhanced by the medium length, L. The solution of this equation sets the initial condition for a non-linear DGLAP-like evolution equation, which describes the collinear early vacuum cascade resolved by the medium at angles exceeding the medium resolution angle, θ c . For asymptotic jet energies, the medium-induced cascade displays an exponential behavior that generalizes the Poisson-like distribution of parton energy loss. This formulation enables the resummation of leading contributions in α s ln(1/R), and α s ln(R/θ c ), and powers of α s L. We briefly explore the limit of strong quenching, where analytic treatments are feasible, offering insights into the impact of parton cascades on jet quenching. These results provide guidance for future numerical simulations and analytical investigations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗