BHFEM3D v. 1.0
Finite element solution to the 3D Poisson equation for computation of DC electric potentials in electromagnetically conducting media.
Weiss, Chester J.↗
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.
Isotropic Standard Poisson in 3D v. 1.0
Finite element solution to the Poisson Equation for electrostatic potntials in an isotropic electrically conducting medium.
Weiss, Chester J.↗
Higher-order space-charge stability in anisotropic beams: Vlasov-Poisson derivation, refined dispersion relations, and stability charts
The Hofmann stability chart is used to screen working points in space-charge-dominated linacs. We identify two errors in its published higher-order dispersion relations: missing $(1\mp2\hatη^2/α)$ factors in the third-order $S^4$ coupling residues, and a sign error in the stated isotropic reduction of the fourth-order relation. Both corrections follow from Hofmann's Vlasov-Poisson equations without fitted parameters. They reproduce coherent tune-shift coefficients in the author's later monograph that the printed forms miss by 24% and 127%. Mode-resolved figures from a published application agree with the corrected relations and reject the printed forms, indicating an inconsistency between the 1998 equations and the calculations underlying those tested figures. We quantify the effect on the non-oscillatory stability chart. Inside the adopted $S^2\le10$ comparison domain, printed and corrected forms disagree on 0.73-2.11% of cells, with no preferred direction. Among excluded cells, disagreement reaches 22%, and the printed relation over-predicts instability at every sampled anisotropy. This concentration may help explain why the errors persisted, although it does not establish their historical cause. For PIP-II, the corrected chart flags four of thirty-two evaluable periods, including one on a third-order odd branch missed by a second-order screen. This count covers non-oscillatory modes only and remains conditional on an unresolved factor-five disagreement between two codes on transverse emittance growth.
Pathak, Abhishek [Fermilab] (ORCID:000000021704208↗
Large-Time Behavior of Solutions to Vlasov-Poisson-Fokker-Planck Equations: From Evanescent Collisions to Diffusive Limit
The present contribution investigates the dynamics generated by the two-dimensional Vlasov-Poisson-Fokker-Planck equation for charged particles in a steady inhomogeneous background of opposite charges. We provide global in time estimates that are uniform with respect to initial data taken in a bounded set of a weighted L{sup 2} space, and where dependencies on the mean-free path τ and the Debye length δ are made explicit. In our analysis the mean free path covers the full range of possible values: from the regime of evanescent collisions τ→∞ to the strongly collisional regime τ→0. As a counterpart, the largeness of the Debye length, that enforces a weakly nonlinear regime, is used to close our nonlinear estimates. Accordingly we pay a special attention to relax as much as possible the τ-dependent constraint on δ ensuring exponential decay with explicit τ-dependent rates towards the stationary solution. In the strongly collisional limit τ→0, we also examine all possible asymptotic regimes selected by a choice of observation time scale. Here also, our emphasis is on strong convergence, uniformity with respect to time and to initial data in bounded sets of a L{sup 2} space. Our proofs rely on a detailed study of the nonlinear elliptic equation defining stationary solutions and a careful tracking and optimization of parameter dependencies of hypocoercive/hypoelliptic estimates.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Stochastic nonlinear wave equation with memory driven by compensated Poisson random measures
In this paper, we study a class of stochastic nonlinear wave equation with memory driven by Lévy noise. We first show the existence and uniqueness of global mild solutions using a suitable energy function. Second, under some additional assumptions we prove the exponential stability of the solutions.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
SPECIAL RELATIVISTIC HYDRODYNAMICS WITH GRAVITATION
Special relativistic hydrodynamics with weak gravity has hitherto been unknown in the literature. Whether such an asymmetric combination is possible has been unclear. Here, the hydrodynamic equations with Poisson-type gravity, considering fully relativistic velocity and pressure under the weak gravity and the action-at-a-distance limit, are consistently derived from Einstein’s theory of general relativity. An analysis is made in the maximal slicing, where the Poisson’s equation becomes much simpler than our previous study in the zero-shear gauge. Also presented is the hydrodynamic equations in the first post-Newtonian approximation, now under the general hypersurface condition. Our formulation includes the anisotropic stress.
79 ASTRONOMY AND ASTROPHYSICS↗
Space charge emission in cylindrical diode
In this paper, a mathematical model to describe cylindrical electron current emissions through a physics approximation method is presented. The proposed mathematical approximation consists of analyzing and solving the nonlinear Poisson's equation, with some determined mathematical restrictions. Our findings tackle the problem when charge-space creates potential barrier that disable the steady-state of the beam propagation. In this problem, the potential barrier effects of electron's speed with zero velocity emitted through the virtual cathode happens. The interaction between particles and the virtual cathode have been to find the inter-atomic potentials as boundary conditions from a quantum mechanics perspective. Furthermore, a non-stationary spatial solution of the electrical potential between anode and cathode is presented. The proposed solution is a 2D differential equation that was linearized from the generalized Poisson equation. A single condition was used solely, throughout the radial boundary conditions of the current density formation.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A generalized Poisson and Poisson-Boltzmann solver for electrostatic environments
The computational study of chemical reactions in complex, wet environments is critical for applications in many fields. It is often essential to study chemical reactions in the presence of applied electrochemical potentials, taking into account the non-trivial electrostatic screening coming from the solvent and the electrolytes. As a consequence, the electrostatic potential has to be found by solving the generalized Poisson and the Poisson-Boltzmann equations for neutral and ionic solutions, respectively. In the present work, solvers for both problems have been developed. A preconditioned conjugate gradient method has been implemented for the solution of the generalized Poisson equation and the linear regime of the Poisson-Boltzmann, allowing to solve iteratively the minimization problem with some ten iterations of the ordinary Poisson equation solver. In addition, a self-consistent procedure enables us to solve the non-linear Poisson-Boltzmann problem. Both solvers exhibit very high accuracy and parallel efficiency and allow for the treatment of periodic, free, and slab boundary conditions. The solver has been integrated into the BigDFT and Quantum-ESPRESSO electronic-structure packages and will be released as an independent program, suitable for integration in other codes.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A meshless stochastic method for Poisson–Nernst–Planck equations
A plethora of biological, physical, and chemical phenomena involve transport of charged particles (ions). Its continuum-scale description relies on the Poisson–Nernst–Planck (PNP) system, which encapsulates the conservation of mass and charge. The numerical solution of these coupled partial differential equations is challenging and suffers from both the curse of dimensionality and difficulty in efficiently parallelizing. We present a novel particle-based framework to solve the full PNP system by simulating a drift–diffusion process with time- and space-varying drift. We leverage Green’s functions, kernel-independent fast multipole methods, and kernel density estimation to solve the PNP system in a meshless manner, capable of handling discontinuous initial states. The method is embarrassingly parallel, and the computational cost scales linearly with the number of particles and dimension. We use a series of numerical experiments to demonstrate both the method’s convergence with respect to the number of particles and computational cost vis-à-vis a traditional partial differential equation solver.
Chemistry↗
Multiphysics simulation of corona discharge induced ionic wind
Ionic wind devices or electrostatic fluid accelerators are becoming of increasing interest as tools for thermal management, in particular for semiconductor devices. In this work, we present a numerical model for predicting the performance of such devices; its main benefit is the ability to accurately predict the amount of charge injected from the corona electrode. Our multiphysics numerical model consists of a highly nonlinear, strongly coupled set of partial differential equations including the Navier-Stokes equations for fluid flow, Poisson's equation for electrostatic potential, charge continuity, and heat transfer equations. To solve this system we employ a staggered solution algorithm that generalizes Gummel's algorithm for charge transport in semiconductors. Predictions of our simulations are verified and validated by comparison with experimental measurements of integral physical quantities, which are shown to closely match.
36 MATERIALS SCIENCE↗
Solution of elliptic partial differential equations by fast Poisson solvers using a local relaxation factor. 1: One-step method
An algorithm for solving a large class of two- and three-dimensional nonseparable elliptic partial differential equations (PDE's) is developed and tested. It uses a modified D'Yakanov-Gunn iterative procedure in which the relaxation factor is grid-point dependent. It is easy to implement and applicable to a variety of boundary conditions. It is also computationally efficient, as indicated by the results of numerical comparisons with other established methods. Furthermore, the current algorithm has the advantage of possessing two important properties which the traditional iterative methods lack; that is: (1) the convergence rate is relatively insensitive to grid-cell size and aspect ratio, and (2) the convergence rate can be easily estimated by using the coefficient of the PDE being solved.
Chang, S. C.↗
The scaling of oblique plasma double layers
Strong oblique plasma double layers are investigated using three methods, i.e., electrostatic particle-in-cell simulations, numerical solutions to the Poisson-Vlasov equations, and analytical approximations to the Poisson-Vlasov equations. The solutions to the Poisson-Vlasov equations and numerical simulations show that strong oblique double layers scale in terms of Debye lengths. For very large potential jumps, theory and numerical solutions indicate that all effects of the magnetic field vanish and the oblique double layers follow the same scaling relation as the field-aligned double layers.
Borovsky, J. E.↗
Instability and dynamics of two nonlinearly coupled intense laser beams in a quantum plasma
We consider nonlinear interactions between two relativistically strong laser beams and a quantum plasma composed of degenerate electron fluids and immobile ions. The collective behavior of degenerate electrons is modeled by quantum hydrodynamic equations composed of the electron continuity, quantum electron momentum (QEM) equation, as well as the Poisson and Maxwell equations. The QEM equation accounts the quantum statistical electron pressure, the quantum electron recoil due to electron tunneling through the quantum Bohm potential, electron-exchange, and electron-correlation effects caused by electron spin, and relativistic ponderomotive forces (RPFs) of two circularly polarized electromagnetic (CPEM) beams. The dynamics of the latter are governed by nonlinear wave equations that include nonlinear currents arising from the relativistic electron mass increase in the CPEM wave fields, as well as from the beating of the electron quiver velocity and electron density variations reinforced by the RPFs of the two CPEM waves. Furthermore, nonlinear electron density variations associated with the driven (by the RPFs) quantum electron plasma oscillations obey a coupled nonlinear Schroedinger and Poisson equations. The nonlinearly coupled equations for our purposes are then used to obtain a general dispersion relation (GDR) for studying the parametric instabilities and the localization of CPEM wave packets in a quantum plasma. Numerical analyses of the GDR reveal that the growth rate of a fastest growing parametrically unstable mode is in agreement with the result that has been deduced from numerical simulations of the governing nonlinear equations. Explicit numerical results for two-dimensional (2D) localized CPEM wave packets at nanoscales are also presented. Possible applications of our investigation to intense laser-solid density compressed plasma experiments are highlighted.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Solution of elliptic partial differential equations by fast Poisson solvers using a local relaxation factor. 2: Two-step method
A two-step semidirect procedure is developed to accelerate the one-step procedure described in NASA TP-2529. For a set of constant coefficient model problems, the acceleration factor increases from 1 to 2 as the one-step procedure convergence rate decreases from + infinity to 0. It is also shown numerically that the two-step procedure can substantially accelerate the convergence of the numerical solution of many partial differential equations (PDE's) with variable coefficients.
Chang, S. C.↗
Hierarchical Statistical 3D ' Atomistic' Simulation of Decanano MOSFETs: Drift-Diffusion, Hydrodynamic and Quantum Mechanical Approaches
When MOSFETs are scaled to deep submicron dimensions the discreteness and randomness of the dopant charges in the channel region introduces significant fluctuations in the device characteristics. This effect, predicted 20 year ago, has been confirmed experimentally and in simulation studies. The impact of the fluctuations on the functionality, yield, and reliability of the corresponding systems shifts the paradigm of the numerical device simulation. It becomes insufficient to simulate only one device representing one macroscopical design in a continuous charge approximation. An ensemble of macroscopically identical but microscopically different devices has to be characterized by simulation of statistically significant samples. The aims of the numerical simulations shift from predicting the characteristics of a single device with continuous doping towards estimating the mean values and the standard deviations of basic design parameters such as threshold voltage, subthreshold slope, transconductance, drive current, etc. for the whole ensemble of 'atomistically' different devices in the system. It has to be pointed out that even the mean values obtained from 'atomistic' simulations are not identical to the values obtained from continuous doping simulations. In this paper we present a hierarchical approach to the 'atomistic' simulation of aggressively scaled decanano MOSFETs. A full scale 3D drift-diffusion'atomostic' simulation approach is first described and used for verification of the more economical, but also more restricted, options. To reduce the processor time and memory requirements at high drain voltage we have developed a self-consistent option based on a thin slab solution of the current continuity equation only in the channel region. This is coupled to the Poisson's equation solution in the whole simulation domain in the Gummel iteration cycles. The accuracy of this approach is investigated in comparison with the full self-consistent solution. At low drain voltage only single solution of the nonlinear Poisson equation is sufficient to extract the current with satisfactory accuracy. A pilot version of a hydrodynamic 'atomistic' simulator has been developed in order to study the effect of the nonequilibrium, non local transport in decanano MOSFETs on the random dopant induced current fluctuations. For the first time we have also applied the density gradient approach in 3D to investigate the effect of the quantum confinement on the threshold voltage fluctuations. The developed 'atomistic' simulation techniques have been applied to study various fluctuation resistant MOSFET architectures including epitaxial and delta doped devices.
Asenov, Asen↗
Optical neural engine for solving scientific partial differential equations
Abstract Solving partial differential equations (PDEs) is the cornerstone of scientific research and development. Data-driven machine learning (ML) approaches are emerging to accelerate time-consuming and computation-intensive numerical simulations of PDEs. Although optical systems offer high-throughput and energy-efficient ML hardware, their demonstration for solving PDEs is limited. Here, we present an optical neural engine (ONE) architecture combining diffractive optical neural networks for Fourier space processing and optical crossbar structures for real space processing to solve time-dependent and time-independent PDEs in diverse disciplines, including Darcy flow equation, the magnetostatic Poisson’s equation in demagnetization, the Navier-Stokes equation in incompressible fluid, Maxwell’s equations in nanophotonic metasurfaces, and coupled PDEs in a multiphysics system. We numerically and experimentally demonstrate the capability of the ONE architecture, which not only leverages the advantages of high-performance dual-space processing for outperforming traditional PDE solvers and being comparable with state-of-the-art ML models but also can be implemented using optical computing hardware with unique features of low-energy and highly parallel constant-time processing irrespective of model scales and real-time reconfigurability for tackling multiple tasks with the same architecture. The demonstrated architecture offers a versatile and powerful platform for large-scale scientific and engineering computations.
Tang, Yingheng (ORCID:0009000153622546)↗
Baryon-Interacting Dark Matter: heating dark matter and the emergence of galaxy scaling relations
The empirical scaling relations observed in disk galaxies remain challenging for models of galaxy formation. The most striking among these is the Mass Discrepancy-Acceleration Relation (MDAR), which encodes both a tight baryonic Tully-Fisher relation (BTFR) and the observed diversity of galaxy rotation curves through the central surface density relation (CSDR). Building on our earlier work [1], we propose here that the MDAR is the result of interactions between baryons and ‘Baryon-Interacting Dark Matter’ (BIDM), which heat up the dark matter. Following a bottom-up, hydrodynamical approach, we find that the MDAR follows if: i) the BIDM equation of state approximates that of an ideal gas; ii) the BIDM relaxation time is order the Jeans time; iii) the heating rate is inversely proportional to the BIDM density. Remarkably, under these assumptions the set of hydrodynamical equations together with Poisson’s equation enjoy an anisotropic scaling symmetry. In the BIDM-dominated regime, this gives rise to an enhanced symmetry which fully captures the low-acceleration limit of the MDAR. We then show that, assuming a cored pseudo-isothermal profile at equilibrium, this set of equations gives rise to parameters reproducing the MDAR. Specifically, in the flat part of the rotation curve the asymptotic rotational velocity matches the parametric dependence of the BTFR. Moreover, in the central region of high-surface brightness galaxies, the profile reproduces the CSDR. Finally, by studying the time-dependent approach to equilibrium, we derive a global combination of the BTFR and CSDR, which matches the expectations in low surface-brightness galaxies. The form of the heating rate also makes model-independent predictions for various cosmological observables. We argue that our scenario satisfies existing observational constraints, and, intriguingly, offers a possible explanation to the EDGES anomaly.
79 ASTRONOMY AND ASTROPHYSICS↗
Charge Transport and Space-Charge Formation in Cd 1 - x Zn x Te 1 - y Se y Radiation Detectors
The electron- and hole-transport properties in cadmium zinc telluride selenide (CZTS) crystals are studied using a laser-induced transient-current technique with pulsed and dc bias. The internal electric field profile and velocity of surface recombination are determined by Monte Carlo simulations of electron and hole transient currents combined with a numerical solution of the drift-diffusion equation coupled with Poisson’s equation. Electron and hole drift mobilities of μe = 830 cm 2 /Vs and μh = 40 cm 2 /Vs, respectively, are determined. We also develop a simple technique for evaluating surface recombination directly from measured current waveforms without the need for numerical simulation. The good quality of the prepared detector at pulsed bias, with electron- and hole-mobility-lifetime products of (μτ)e = 1.9 × 10 -3 cm 2 /V and (μτ)h = 1.4 × 10 -4 cm 2 /V, respectively, are observed. The formation of a positive space charge, originating from hole injection combined with a recombination level, is found. We observe a significant position dependence of the lifetime of electrons and holes in dc bias due to hole injection. The experiment is successfully fitted by a simple model dominated by a single deep recombination level with an energy of E t =E C -0.73eV; concentration of 7.3 × 10 11 cm -3 ; and electron- and hole-capture cross sections of 3.5 × 10 -14 cm 2 and 6.5 × 10 -14 cm 2 , respectively.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗