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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

NASA NTRS1982

A numerical method for solving the Vlasov equation

A numerical procedure is derived for the solution of the Vlasov-Poisson system of equations in two phase-space variables. Derivatives with respect to the phase-space variables are approximated by a weighted sum of the values of the distribution function at property chosen neighboring points. The resulting set of ordinary differential equations is then solved by using an appropriate time intergration scheme. The accuracy of the proposed method is tested with some simple model problems. The results for the free streaming case, linear Landau damping, and nonlinear Landau damping are investigated and compared with those of the splitting scheme. The proposed method is found to be very accurate and efficient.

Satofuka, N.↗

DOE OSTI2015

Theory of Current Transients in Planar Semiconductor Devices: Insights and Applications to Organic Solar Cells

Time-domain current measurements are widely used to characterize semiconductor material properties, such as carrier mobility, doping concentration, carrier lifetime, and the static dielectric constant. It is therefore critical that these measurements be theoretically understood if they are to be successfully applied to assess the properties of materials and devices. In this paper, we derive generalized relations for describing current-density transients in planar semiconductor devices at uniform temperature. By spatially averaging the charge densities inside the semiconductor, we are able to provide a rigorous, straightforward, and experimentally relevant way to interpret these measurements. The formalism details several subtle aspects of current transients, including how the electrode charge relates to applied bias and internal space charge, how the displacement current can alter the apparent free-carrier current, and how to understand the integral of a charge-extraction transient. We also demonstrate how the formalism can be employed to derive the current transients arising from simple physical models, like those used to describe charge extraction by linearly increasing voltage (CELIV) and time-of-flight experiments. In doing so, we find that there is a nonintuitive factor-of-2 reduction in the apparent free-carrier concentration that can be easily missed, for example, in the application of charge-extraction models. Finally, to validate our theory and better understand the different current contributions, we perform a full time-domain drift-diffusion simulation of a CELIV trace and compare the results to our formalism. As expected, our analytic equations match precisely with the numerical solutions to the drift-diffusion, Poisson, and continuity equations. Furthermore, overall, our formalism provides a straightforward and general way to think about how the internal space-charge distribution, the electrode charge, and the externally applied bias translate into a measured current transient in a planar semiconductor device.

36 MATERIALS SCIENCE↗

NASA NTRS1988

Consistent boundary conditions for Reduced Navier-Stokes (RNS) scheme applied to 3-dimensional internal viscous flows

A consistent and efficient set of boundary conditions is developed for the multi-sweep space-marching pressure-elliptic Reduced Navier-Stokes (RNS) scheme as applied for 3-D internal viscous flow problems. No-slip boundary conditions are directly imposed on the solid walls. There is no iteration procedure required in the cross plane to ensure mass conservation across each marching plane. The finite difference equations forming the coefficient matrix are ordered such that the surface normal velocity is specified on all the solid walls; unlike external flows, a pressure boundary condition in the cross plane is not required. Since continuity is directly satisfied at all points in the flow domain, the first order momentum equations can be solved directly for the pressure without the need for a Poisson pressure correction equation. The procedure developed herein can also be applied with periodic boundary conditions. The analysis is given for general compressible flows. Incompressible flow solutions are obtained, for straight and curved ducts of square cross section, to validate the procedure. These solutions are used to demonstrate the applicability of the RNS scheme, with the improved boundary conditions for internal flows with strong interaction, as would be encountered in ducts and turbomachinery geometries.

Reddy, D. R.↗

DOE OSTI2016

cDPD: A new dissipative particle dynamics method for modeling electrokinetic phenomena at the mesoscale

In this work, we develop a “charged” dissipative particle dynamics (cDPD) model for simulating mesoscopic electrokinetic phenomena governed by the stochastic Poisson-Nernst-Planck and the Navier-Stokes equations. In particular, the transport equations of ionic species are incorporated into the DPD framework by introducing extra degrees of freedom and corresponding evolution equations associated with each DPD particle. Diffusion of ionic species driven by the ionic concentration gradient, electrostatic potential gradient, and thermal fluctuations is captured accurately via pairwise fluxes between DPD particles. The electrostatic potential is obtained by solving the Poisson equation on the moving DPD particles iteratively at each time step. For charged surfaces in bounded systems, an effective boundary treatment methodology is developed for imposing both the correct hydrodynamic and electrokinetics boundary conditions in cDPD simulations. To validate the proposed cDPD model and the corresponding boundary conditions, we first study the electrostatic structure in the vicinity of a charged solid surface, i.e., we perform cDPD simulations of the electrostatic double layer and show that our results are in good agreement with the well-known mean-field theoretical solutions. Furthermore, we simulate the electrostatic structure and capacity densities between charged parallel plates in salt solutions with different salt concentrations. Moreover, we employ the proposed methodology to study the electro-osmotic and electro-osmotic/pressure-driven flows in a micro-channel. In the latter case, we simulate the dilute poly-electrolyte solution drifting by electro-osmotic flow in a micro-channel, hence demonstrating the flexibility and capability of this method in studying complex fluids with electrostatic interactions at the micro- and nano-scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Train small, model big: Scalable physics simulators via reduced order modeling and domain decomposition

Numerous cutting-edge scientific technologies originate at the laboratory scale, but transitioning them to practical industry applications is a formidable challenge. Traditional pilot projects at intermediate scales are costly and time-consuming. An alternative, the pilot-scale model, relies on high-fidelity numerical simulations, but even these simulations can be computationally prohibitive at larger scales. To overcome these limitations, we propose a scalable, physics-constrained reduced order model (ROM) method. The ROM identifies critical physics modes from small-scale unit components, projecting governing equations onto these modes to create a reduced model that retains essential physics details. We also employ Discontinuous Galerkin Domain Decomposition (DG-DD) to apply ROM to unit components and interfaces, enabling the construction of large-scale global systems without data at such large scales. Here this method is demonstrated on the Poisson and Stokes flow equations, showing that it can solve equations about 15–40 times faster with only ~1% relative error. Furthermore, ROM takes one order of magnitude less memory than the full order model, enabling larger scale predictions at a given memory limitation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2013

Theory of spatially non-symmetric kinetic equilibria for collisionless plasmas

The problem posed by the possible existence/non-existence of spatially non-symmetric kinetic equilibria has remained unsolved in plasma theory. For collisionless magnetized plasmas, this involves the construction of stationary solutions of the Vlasov-Maxwell equations. In this paper, the issue is addressed for non-relativistic plasmas both in astrophysical and laboratory contexts. The treatment is based on a Lagrangian variational description of single-particle dynamics. Starting point is a non-perturbative formulation of gyrokinetic theory, which allows one to construct 'a posteriori' with prescribed order of accuracy an asymptotic representation for the magnetic moment. In terms of the relevant particle adiabatic invariants generalized bi-Maxwellian equilibria are proved to exist. These are shown to recover, under suitable assumptions, a Chapman-Enskog form which permits an analytical treatment of the corresponding fluid moments. In particular, the constrained posed by the Poisson and the Ampere equations are analyzed, both for quasi-neutral and non-neutral plasmas. The conditions of existence of the corresponding non-symmetric kinetic equilibria are investigated. As a notable feature, both astrophysical and laboratory plasmas are shown to exhibit, under suitable conditions, a kinetic dynamo, whereby the equilibrium magnetic field can be self-generated by the equilibrium plasma currents.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1990

SINDA '85/FLUINT, Version 2.2

SINDA, Systems Improved Numerical Differencing Analyzer, software system developed for solving physical problems governed by diffusion-type equations modeled by lumped-parameter representations. Used as general thermal analyzer with resistor-and-capacitor-network representations, adapted to wide range of problems represented by such differential equations as Fourier, Poisson, or LaPlace. Software system updated. FLUINT, FLUid INTegrator, is advanced, one-dimensional fluid-analysis program solving equations of arbitrary fluid-flow networks. Written in FORTRAN and ASSEMBLER.

Cullimore, Brent↗

NASA NTRS1993

Incompressible spectral-element method: Derivation of equations

A fractional-step splitting scheme breaks the full Navier-Stokes equations into explicit and implicit portions amenable to the calculus of variations. Beginning with the functional forms of the Poisson and Helmholtz equations, we substitute finite expansion series for the dependent variables and derive the matrix equations for the unknown expansion coefficients. This method employs a new splitting scheme which differs from conventional three-step (nonlinear, pressure, viscous) schemes. The nonlinear step appears in the conventional, explicit manner, the difference occurs in the pressure step. Instead of solving for the pressure gradient using the nonlinear velocity, we add the viscous portion of the Navier-Stokes equation from the previous time step to the velocity before solving for the pressure gradient. By combining this 'predicted' pressure gradient with the nonlinear velocity in an explicit term, and the Crank-Nicholson method for the viscous terms, we develop a Helmholtz equation for the final velocity.

Deanna, Russell G.↗

DOE OSTI2026

A tensor train-based isogeometric solver for large-scale 3D poisson problems

We introduce a three-dimensional (3D), fully tensor train (TT) assembled isogeometric analysis (IGA) framework, TT-IGA, for solving partial differential equations (PDEs). Our method reformulates IGA discrete operators into TT format, enabling efficient compression and computation. Geometry evaluations use the original NURBS description at sampling points and TT approximation is applied to geometry-derived coefficient fields and discrete operators. We demonstrate the effectiveness of the proposed TT-IGA framework on the three-dimensional Poisson equation, achieving substantial reductions in memory and computational cost without compromising solution quality.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2014

Higher-order Hamiltonian fluid reduction of Vlasov equation

From the Hamiltonian structure of the Vlasov equation, we build a Hamiltonian model for the first three moments of the Vlasov distribution function, namely, the density, the momentum density and the specific internal energy. We derive the Poisson bracket of this model from the Poisson bracket of the Vlasov equation, and we discuss the associated Casimir invariants.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2013

A mixed quantum-classical Liouville study of the population dynamics in a model photo-induced condensed phase electron transfer reaction

We apply two approximate solutions of the quantum-classical Liouville equation (QCLE) in the mapping representation to the simulation of the laser-induced response of a quantum subsystem coupled to a classical environment. These solutions, known as the Poisson Bracket Mapping Equation (PBME) and the Forward-Backward (FB) trajectory solutions, involve simple algorithms in which the dynamics of both the quantum and classical degrees of freedom are described in terms of continuous variables, as opposed to standard surface-hopping solutions in which the classical degrees of freedom hop between potential energy surfaces dictated by the discrete adiabatic state of the quantum subsystem. The validity of these QCLE-based solutions is tested on a non-trivial electron transfer model involving more than two quantum states, a time-dependent Hamiltonian, strong subsystem-bath coupling, and an initial energy shift between the donor and acceptor states that depends on the strength of the subsystem-bath coupling. In particular, we calculate the time-dependent population of the photoexcited donor state in response to an ultrafast, on-resonance pump pulse in a three-state model of an electron transfer complex that is coupled asymmetrically to a bath of harmonic oscillators through the optically dark acceptor state. Within this approach, the three-state electron transfer complex is treated quantum mechanically, while the bath oscillators are treated classically. When compared to the more accurate QCLE-based surface-hopping solution and to the numerically exact quantum results, we find that the PBME solution is not capable of qualitatively capturing the population dynamics, whereas the FB solution is. However, when the subsystem-bath coupling is decreased (which also decreases the initial energy shift between the donor and acceptor states) or the initial shift is removed altogether, both the PBME and FB results agree better with the QCLE-based surface-hopping results. These findings highlight the challenges posed by various conditions such as a time-dependent external field, the strength of the subsystem-bath coupling, and the degree of asymmetry on the accuracy of the PBME and FB algorithms.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1974

Computer solutions of Wu's equations for compressible flow through turbomachines

Two computers programs, known as Matrix Through-Flow and Matrix Blade-To-Blade, for analyzing the meridional and blade-to-blade flow patterns are described. The numerical solutions are obtained by finite difference approximations to the governing Poisson-type differential equations for the stream function. Solutions for several turbomachines, giving flow patterns and velocity distributions, are included.

Smith, D. J. L.↗

NASA NTRS1975

Cost performance satellite design using queueing theory

A modified Poisson arrival, infinite server queuing model is used to determine the effects of limiting the number of broadcast channels (C) of a direct broadcast satellite used for public service purposes (remote health care, education, etc.). The model is based on the reproductive property of the Poisson distribution. A difference equation has been developed to describe the change in the Poisson parameter. When all initially delayed arrivals reenter the system a (C plus 1) order polynomial must be solved to determine the effective value of the Poisson parameter. When less than 100% of the arrivals reenter the system the effective value must be determined by solving a transcendental equation. The model was used to determine the minimum number of channels required for a disaster warning satellite without degradation in performance. Results predicted by the queuing model were compared with the results of digital simulation.

Hein, G. F.↗

DOE OSTI2026

Integral Kernel Methods for Nonlinear Parabolic-Elliptic Systems

Nonlinear parabolic-elliptic systems arise in many physical, biological, and chemical phenomena such as chemotaxis, ion transport, self-gravitating particles, and Brownian vortices. Existing methods struggle with the strong coupling and high nonlinearity and nonlocality of some of these systems, especially the ill-conditioned, convection-dominated problems. To overcome numerical difficulties, current approaches rely on initial guesses, preconditioning, or iterative techniques with no convergence guarantees. They might suffer from poor scalability, large memory usage, and difficulty to parallelize. Inspired by the connection of parabolic-elliptic systems to stochastic processes, we introduce a novel meshless, monolithic, and fully explicit method that naturally encapsulates the elliptic and parabolic operators into a single step which updates each node deterministically with global information. By being fully quadrature-based, it avoids solving systems of discretized equations and does not utilize initial guesses or preconditioning, while requiring little memory and being easy to parallelize. We first derive the method in an integral kernel formulation with quadratic complexity in the number of integration nodes and then leverage kernel-independent fast multipole methods (FMM) to present a scalable algorithm with linear complexity. We provide numerical examples for the Poisson-Nernst-Planck equations in one, two, and three dimensions, together with the derivation of the integral kernel for each case. Furthermore, the examples demonstrate the fast convergence and scalability of the FMM-accelerated algorithm, as well as its suitability for convection-dominated problems, making it competitive against traditional PDE solvers.

PDE systems↗

NASA NTRS1991

Charge collected by diffusion from an ion track under mixed boundary conditions

Charge-carrier diffusion from an ion track in a silicon substrate at least a few hundred microns thick is analyzed. The substrate upper surface is treated as reflective except for a small section, intended to represent a reverse-biased junction, which is treated as a sink. Total charge collected by the sink is calculated by assuming transport to be governed by an ambipolar diffusion equation with temporally constant and spatially uniform carrier lifetime and diffusion coefficient. Present results apply to a normally incident track but could easily be generalized to arbitrary track direction. The collected charge is found to depend on track length and on the electrostatic capacitance, rather than the area, of the sink. Theoretical predictions are compared to the results of a numerical simulation called the Poisson and Continuity Equation Solver (PISCES) for three cases and are found to agree within a factor of two in the worst case.

Edmonds, Larry D.↗

NASA NTRS1995

Multiparameter linear least-squares fitting to Poisson data one count at a time

A standard problem in gamma-ray astronomy data analysis is the decomposition of a set of observed counts, described by Poisson statistics, according to a given multicomponent linear model, with underlying physical count rates or fluxes which are to be estimated from the data. Despite its conceptual simplicity, the linear least-squares (LLSQ) method for solving this problem has generally been limited to situations in which the number n(sub i) of counts in each bin i is not too small, conventionally more than 5-30. It seems to be widely believed that the failure of the LLSQ method for small counts is due to the failure of the Poisson distribution to be even approximately normal for small numbers. The cause is more accurately the strong anticorrelation between the data and the wieghts w(sub i) in the weighted LLSQ method when square root of n(sub i) instead of square root of bar-n(sub i) is used to approximate the uncertainties, sigma(sub i), in the data, where bar-n(sub i) = E(n(sub i)), the expected value of N(sub i). We show in an appendix that, avoiding this approximation, the correct equations for the Poisson LLSQ (PLLSQ) problems are actually identical to those for the maximum likelihood estimate using the exact Poisson distribution. We apply the method to solve a problem in high-resolution gamma-ray spectroscopy for the JPL High-Resolution Gamma-Ray Spectrometer flown on HEAO 3. Systematic error in subtracting the strong, highly variable background encountered in the low-energy gamma-ray region can be significantly reduced by closely pairing source and background data in short segments. Significant results can be built up by weighted averaging of the net fluxes obtained from the subtraction of many individual source/background pairs. Extension of the approach to complex situations, with multiple cosmic sources and realistic background parameterizations, requires a means of efficiently fitting to data from single scans in the narrow (approximately = 1.2 keV, HEAO 3) energy channels of a Ge spectrometer, where the expected number of counts obtained per scan may be very low. Such an analysis system is discussed and compared to the method previously used.

Wheaton, Wm. A.↗

DOE OSTI2021

Microstructural changes and mass transport in cement-based materials: A modeling approach

Highlights: • Concept to relate microstructural changes, mass transport, and chemical reactions in cement-based materials is presented. • Multi-physics modeling approach utilizing an extended version of the PNP equation and the geochemical code IPHREEQC. • Case studies indicate significant changes in mass transport properties due to alterations in the microstructure. • Microstructural changes depend on exposure and interactions between chemical reactions, pore structure, and mass transport. A generic modeling framework is presented to relate microstructural changes, i.e. changes in the pore volume distribution and phase assemblage, and mass transport in cement-based materials. The modeling framework accounts for mass transport and chemical equilibrium between ions in the pore solution and solid hydrates by means of an extended version of the Poisson-Nernst-Planck equation, in which the chemical equilibrium is solved by the external geochemical code IPHREEQC. Results of numerical studies concerning carbonation and leaching of a cement-based material indicate significant changes in mass transport properties due to chemical reactions leading to alterations in the microstructure. The resulting highly non-linear (both spatial and temporal) microstructural changes are thereby depending on the exposure and indicate that interactions between chemical reactions, pore structure changes, and mass transport have to be accounted for as a whole.

36 MATERIALS SCIENCE↗

NASA NTRS2005

Full Multigrid Flow Solver

FMG3D (full multigrid 3 dimensions) is a pilot computer program that solves equations of fluid flow using a finite difference representation on a structured grid. Infrastructure exists for three dimensions but the current implementation treats only two dimensions. Written in Fortran 90, FMG3D takes advantage of the recursive subroutine feature, dynamic memory allocation, and structured-programming constructs of that language. FMG3D supports multi-block grids with three types of block-to-block interfaces: periodic, C-zero, and C-infinity. For all three types, grid points must match at interfaces. For periodic and C-infinity types, derivatives of grid metrics must be continuous at interfaces. The available equation sets are as follows: scalar elliptic equations, scalar convection equations, and the pressure-Poisson formulation of the Navier-Stokes equations for an incompressible fluid. All the equation sets are implemented with nonzero forcing functions to enable the use of user-specified solutions to assist in verification and validation. The equations are solved with a full multigrid scheme using a full approximation scheme to converge the solution on each succeeding grid level. Restriction to the next coarser mesh uses direct injection for variables and full weighting for residual quantities; prolongation of the coarse grid correction from the coarse mesh to the fine mesh uses bilinear interpolation; and prolongation of the coarse grid solution uses bicubic interpolation.

Mineck, Raymond E.↗