Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

A stochastic-dynamic model for global atmospheric mass-field statistics

Global atmospheric mass field error correlations based on satellite observations and on numerical forecasts show strong and systematic latitude dependence. A model for the latitude dependent spatial correlation structure of mass field forecast errors is derived from dynamical considerations. Three methods of solution were tested. In the first method, the equation was solved by expansion in spherical harmonics, and the correlation function was computed analytically using the expansion coefficients. In the second method, the finite difference equivalent of the equation was solved using a fast poisson solver. The correlation function was computed using stratified sampling of the individual realizations. In the third method, a higher order equation was derived, and solved directly in finite differences by two successive applications of the fast poisson solver. The three methods were compared for accuracy and efficiency, and the third method was chosen as clearly superior.

Ghil, M.↗
DOE OSTI2012

Integrable string models of WZNW model type with constant SU(2), SO(3), SP(2) and SU(3) torsion and hydrodynamic chains

The integrability of string model of WZNW model type with constant SU (2), SO (3), SP (2) torsion is investigated. The closed boson string model in the background gravity and antisymmetric B-field is considered as integrable system in terms of initial chiral currents. The model is considered under assumption that internal torsion related with metric of Riemann-Cartan space and external torsion related with antisymmetric B-field are (anti)coincide. New equation of motion and exact solution this equation was obtained for string model with constant SU (2), SO (3), SP (2) torsion. New equations of motion and new Poisson brackets (PB) for infinite dimensional hydrodynamic chains was obtained for string model with constant SU (n ), SO (n ), SP (n ) torsion for n → ∞.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1974

Computation of solar perturbations with Poisson series

Description of a project for computing first-order perturbations of natural or artificial satellites by integrating the equations of motion on a computer with automatic Poisson series expansions. A basic feature of the method of solution is that the classical variation-of-parameters formulation is used rather than rectangular coordinates. However, the variation-of-parameters formulation uses the three rectangular components of the disturbing force rather than the classical disturbing function, so that there is no problem in expanding the disturbing function in series. Another characteristic of the variation-of-parameters formulation employed is that six rather unusual variables are used in order to avoid singularities at the zero eccentricity and zero (or 90 deg) inclination. The integration process starts by assuming that all the orbit elements present on the right-hand sides of the equations of motion are constants. These right-hand sides are then simple Poisson series which can be obtained with the use of the Bessel expansions of the two-body problem in conjunction with certain interation methods. These Poisson series can then be integrated term by term, and a first-order solution is obtained.

Broucke, R.↗
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 OSTI2013

Conservation of energy and momentum in nonrelativistic plasmas

Conservation laws of energy and momentum for nonrelativistic plasmas are derived from applying Noether's theorem to the action integral for the Vlasov-Poisson-Ampere system [Sugama, Phys. Plasmas 7, 466 (2000)]. The symmetric pressure tensor is obtained from modifying the asymmetric canonical pressure tensor with using the rotational symmetry of the action integral. Differences between the resultant conservation laws and those for the Vlasov-Maxwell system including the Maxwell displacement current are clarified. These results provide a useful basis for gyrokinetic conservation laws because gyrokinetic equations are derived as an approximation of the Vlasov-Poisson-Ampere system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Fluctuating hydrodynamics of reactive liquid mixtures

Fluctuating hydrodynamics (FHD) provides a framework for modeling microscopic fluctuations in a manner consistent with statistical mechanics and nonequilibrium thermodynamics. This paper presents an FHD formulation for isothermal reactive incompressible liquid mixtures with stochastic chemistry. Fluctuating multispecies mass diffusion is formulated using a Maxwell-Stefan description without assuming a dilute solution, and momentum dynamics is described by a stochastic Navier-Stokes equation for the fluid velocity. We consider a thermodynamically consistent generalization for the law of mass action for non-dilute mixtures and use it in the chemical master equation (CME) to model reactions as a Poisson process. The FHD approach provides remarkable computational efficiency over traditional reaction-diffusion master equation methods when the number of reactive molecules is large, while also retaining accuracy even when there are as few as ten reactive molecules per hydrodynamic cell. We present a numerical algorithm to solve the coupled FHD and CME equations and validate it on both equilibrium and nonequilibrium problems. We simulate a diffusively driven gravitational instability in the presence of an acid-base neutralization reaction, starting from a perfectly flat interface. We demonstrate that the coupling between velocity and concentration fluctuations dominates the initial growth of the instability.

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

Filtering with Marked Point Process Observations via Poisson Chaos Expansion

We study a general filtering problem with marked point process observations. The motivation comes from modeling financial ultra-high frequency data. First, we rigorously derive the unnormalized filtering equation with marked point process observations under mild assumptions, especially relaxing the bounded condition of stochastic intensity. Then, we derive the Poisson chaos expansion for the unnormalized filter. Based on the chaos expansion, we establish the uniqueness of solutions of the unnormalized filtering equation. Moreover, we derive the Poisson chaos expansion for the unnormalized filter density under additional conditions. To explore the computational advantage, we further construct a new consistent recursive numerical scheme based on the truncation of the chaos density expansion for a simple case. The new algorithm divides the computations into those containing solely system coefficients and those including the observations, and assign the former off-line.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Homogeneous systems with quadratic integrals, Lie-Poisson quasibrackets, and Kovalevskaya's method

We consider differential equations with quadratic right-hand sides that admit two quadratic first integrals, one of which is a positive-definite quadratic form. We indicate conditions of general nature under which a linear change of variables reduces this system to a certain 'canonical' form. Under these conditions, the system turns out to be divergenceless and can be reduced to a Hamiltonian form, but the corresponding linear Lie-Poisson bracket does not always satisfy the Jacobi identity. In the three-dimensional case, the equations can be reduced to the classical equations of the Euler top, and in four-dimensional space, the system turns out to be superintegrable and coincides with the Euler-Poincaré equations on some Lie algebra. In the five-dimensional case we find a reducing multiplier after multiplying by which the Poisson bracket satisfies the Jacobi identity. In the general case for n>5 we prove the absence of a reducing multiplier. As an example we consider a system of Lotka-Volterra type with quadratic right-hand sides that was studied by Kovalevskaya from the viewpoint of conditions of uniqueness of its solutions as functions of complex time. Bibliography: 38 titles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Correlation functions from a unified variational principle: Trial Lie groups

Time-dependent expectation values and correlation functions for many-body quantum systems are evaluated by means of a unified variational principle. It optimizes a generating functional depending on sources associated with the observables of interest. It is built by imposing through Lagrange multipliers constraints that account for the initial state (at equilibrium or off equilibrium) and for the backward Heisenberg evolution of the observables. The trial objects are respectively akin to a density operator and to an operator involving the observables of interest and the sources. We work out here the case where trial spaces constitute Lie groups. This choice reduces the original degrees of freedom to those of the underlying Lie algebra, consisting of simple observables; the resulting objects are labeled by the indices of a basis of this algebra. Explicit results are obtained by expanding in powers of the sources. Zeroth and first orders provide thermodynamic quantities and expectation values in the form of mean-field approximations, with dynamical equations having a classical Lie–Poisson structure. At second order, the variational expression for two-time correlation functions separates–as does its exact counterpart–the approximate dynamics of the observables from the approximate correlations in the initial state. Two building blocks are involved: (i) a commutation matrix which stems from the structure constants of the Lie algebra; and (ii) the second-derivative matrix of a free-energy function. The diagonalization of both matrices, required for practical calculations, is worked out, in a way analogous to the standard RPA. The ensuing structure of the variational formulae is the same as for a system of non-interacting bosons (or of harmonic oscillators) plus, at non-zero temperature, classical Gaussian variables. This property is explained by mapping the original Lie algebra onto a simpler Lie algebra. The results, valid for any trial Lie group, fulfill consistency properties and encompass several special cases: linear responses, static and time-dependent fluctuations, zero- and high-temperature limits, static and dynamic stability of small deviations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2021

Higher Ion Selectivity with Lower Energy Usage Promoted by Electro-osmotic Flow in the Transport through Conical Nanopores

The tradeoff between selectivity and throughput presents fundamental challenges to improve desalination and charge storage, salinity gradient-based energy harvesting, memory device/circuit development, and so forth. The well-known ion-current rectification and the recently resolved time-dependent transport hysteresis in conical nanopores or asymmetric nanointerfaces offer new opportunities for the selective transport of matter. This report shows that electro-osmotic flow (EOF) is an overlooked factor that increases ion selectivity while maintaining enhanced transport throughput in rectified nanoscale electrokinetic transport. The increased selectivity originates primarily from the suppression of anion flux by the fluid flow in the opposite direction under the applied electrical field. By solving the Poisson and Nernst–Planck (PNP) equations without and with coupled Navier–Stokes (PNP–NS) equations, the EOF effects on cation and anion transport are unequivocally revealed in asymmetric nanopipettes. The flux of cations and anions as well as the transference number and flow velocity are elucidated using the models and boundary conditions validated by previous experiments. A dimensionless parameter, radius over the Debye length, reveals optimal ion selectivity and energy cost at intermediate ion concentrations and nanopore sizes, up to hundreds of millimolars and tens of nanometers under time-dependent potential stimulus. Furthermore, the fundamental insights into EOF at nanointerfaces suggest new routes/strategies for better separation, analysis, and energy applications.

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

Hydrogen and Electric Power Generation from Liquid Microjets: Design Principles for Optimizing Conversion Efficiency

Historically, Liquid water microjets have been successfully employed for both electrical power generation and gaseous hydrogen production, but the demonstrated efficiencies have been low. Here, we employ a combination of a modified Poisson-Boltzmann description, continuum hydrodynamic equations, and microjet electrokinetic experiments to gain detailed insight into the origin of the streaming currents produced in pure water. We identify the contributions to the streaming current from specific ion adsorption at the solid/liquid interface and from long-ranged electrostatic interactions, finding that the portion originating from the latter dominate at charged surfaces. The detailed understanding afforded by theory and the close agreement with experimental results elucidates design principles for optimizing hydrogen production and power generation. Changing the sign of the surface charge density through targeted use of surface coatings via silanization switches the primary charge carrier between hydronium and hydroxide and therefore switches the corresponding production of molecular hydrogen to oxygen at the target electrode. Moreover, hydrophobic surface coatings reduce dissipation due to fluid/solid friction, thereby increasing the conversion efficiency.

08 HYDROGEN↗
DOE OSTI2019

Why the particle-in-cell method captures instability enhanced collisions

While the particle-in-cell (PIC) method has been the subject of years of theoretical study, the common misconception that PIC approximates the collisionless Vlasov-Maxwell system of equations, or the Vlasov-Poisson system for the electrostatic case, is widely cited. Here, a PIC-relevant generalization of the instability enhanced Lenard-Balescu collision operator is derived. The analysis of this collision operator demonstrates that while coulomb collisions are truncated by the electric field grid resolution, longer range collisions that arise from imperfect shielding in the presence of instabilities persist. These instability enhanced collisions are completely captured by PIC as long as there is sufficient grid resolution to resolve the wavenumber of the unstable modes of the instability. The inclusion of this behavior is akin to particle wave interactions in the Vlasov based quasilinear theory but with one important difference. Quasilinear theory requires an initial spectral energy density which cannot be supplied self-consistently from within the theory because its initial value is determined by non-Vlasov collisional effects. With the proper electric field grid resolution, the initial spectral energy density is included self-consistently, along with the generation of plasma waves originating from discrete particle motion. Finally, predictions of the grid resolution effect are found to be in agreement with PIC simulations at varying grid resolutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

A complete quasiclassical map for the dynamics of interacting fermions

We present a strategy for mapping the dynamics of a fermionic quantum system to a set of classical dynamical variables. The approach is based on imposing the correspondence relation between the commutator and the Poisson bracket, preserving Heisenberg's equation of motion for one-body operators. In order to accommodate the effect of two-body terms, we further impose quantization on the spin-dependent occupation numbers in the classical equations of motion, with a parameter that is determined self-consistently. Expectation values for observables are taken with respect to an initial quasiclassical distribution that respects the original quantization of the occupation numbers. The proposed classical map becomes complete under the evolution of quadratic Hamiltonians and is extended for all even order observables. We show that the map provides an accurate description of the dynamics for an interacting quantum impurity model in the coulomb blockade regime, at both low and high temperatures. In conclusion, the numerical results are aided by a novel importance sampling scheme that employs a reference system to reduce significantly the sampling effort required to converge the classical calculations.

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↗
NASA NTRS1991

Self-consistent analysis of lattice-matched and pseudomorphic quantum-well emission transistors

A self-consistent analysis of the quantum-well emission transistor (QWET) is presented allowing an exact calculation of the device quantum properties. Poisson's and Schroedinger's equation are solved numerically using a finite-difference method on a self-consistent basis. Pseudomorphic AlGaAs/InGaAs designs with 15-20 percent excess In are suggested for improving the device performance. Design with doping in various parts of the QWET are also studied. This analysis reveals that the device performance is less optimistic than previously predicted by analytic approaches. By introducing the pseudomorphic channel principle, while maintaining a reasonably low Al content for the gate and collector layers, it is, however, possible to obtain satisfactory performance. Optimum pseudomorphic designs showed high current driving capability (200,000 A/sq cm), high transconductance (3S/mm), and small intrinsic delay time (2 ps).

Hong, Kyushik↗