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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Finite-Difference Simulation of Eddy Currents in Nonmagnetic Sheets via Electric Vector Potential

The magnetic excitation of eddy currents in thin metal sheets is a mathematically challenging problem with many practical applications. Since the eddy current density is a solenoidal field on a 2-D plane, we may express it as the curl of a separate vector field with only one scalar component. This new field is governed by a unique form of the vector-Poisson equation that is readily solved numerically through standard finite-difference techniques. For the special case of a weakly induced eddy current, the self-inductance terms can be neglected to produce a sparse system matrix that is easily inverted. For a strongly induced eddy current, self-inductance must be included at the cost of a more complex and dense system matrix. The method is validated against the CST EM Studio software suite by producing nearly identical results in only a tiny fraction of the computational time.

42 ENGINEERING↗
DOE OSTI2021

Sparse Approximate Multifrontal Factorization with Butterfly Compression for High-Frequency Wave Equations

In this work, we present a fast and approximate multifrontal solver for large-scale sparse linear systems arising from finite-difference, finite-volume or finite-element discretization of high-frequency wave equations. The proposed solver leverages the butterfly algorithm and its hierarchical matrix extension for compressing and factorizing large frontal matrices via graph-distance guided entry evaluation or randomized matrix-vector multiplication-based schemes. Complexity analysis and numerical experiments demonstrate $\mathcal{O}(N\log^2 N)$ computation and $\mathcal{O}(N)$ memory complexity when applied to an $N\times N$ sparse system arising from 3D high-frequency Helmholtz and Maxwell problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Learning Robust Marking Policies for Adaptive Mesh Refinement

Here in this work, we revisit the marking decisions made in the standard adaptive finite element method (AFEM). Experience shows that a naïve marking policy leads to inefficient use of computational resources for adaptive mesh refinement (AMR). Consequently, using AMR in practice often involves ad-hoc or time-consuming offline parameter tuning to set appropriate parameters for the marking subroutine. To address these practical concerns, we recast AMR as a Markov decision process in which refinement parameters can be selected on-the-fly at run time, without the need for pre-tuning by expert users. In this new paradigm, the refinement parameters are also chosen adaptively via a marking policy that can be optimized using methods from reinforcement learning. We use the Poisson equation to demonstrate our techniques on h- and hp-refinement benchmark problems, and our experiments suggest that superior marking policies remain undiscovered for many classical AFEM applications. Furthermore, an unexpected observation from this work is that marking policies trained on one family of PDEs are sometimes robust enough to perform well on problems far outside the training family. For illustration, we show that a simple hp-refinement policy trained on 2D domains with only a single re-entrant corner can be deployed on far more complicated 2D domains, and even 3D domains, without significant performance loss. For reproduction and broader adoption, we accompany this work with an open-source implementation of our methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

FLUXNET v. 1.0

Solves the Poisson equation on three-dimensional network of line segments, each of uniform material properties, subject to boundary condition forcing.

Weiss, Chester↗
DOE OSTI2019

Adaptive Sparse Grid Discretization (ASGarD)

This code solves the Vlasov-Poisson equations in 6-dimensions without approximation. It uses sparse grid discretization to solve the high-dimensional PDEs. The code functions as a discontinuous Galerkin finite element solver. Advantages/benefits include no approximations, supports high-dimensionality, leadership-scale capability, and implicit time advance.

Hahn, Steven [Oak Ridge National Laboratory (ORNL)↗
DOE OSTI2023

Scalable Reduced Order Model with Discontinuous Galerkin Domain Decomposition

scaleupROM is a scalable, physics-constrained reduced order model (ROM). It aims to provide robust, accelerated physics predictions at extrapolated scales, based on the small, component-level data. This is implemented by combining projection-based ROM with discontinuous Galerkin domain decomposition, in the framework of MFEM and libROM. It currently supports the Poisson equation and Stokes flow equation, and more work is in progress toward general, nonlinear physics systems.

Chung, Seung Whan↗
DOE OSTI2018

Band Bending/Contact Resistivity Codes v. 1

Code #1 solves the one-dimensional Poisson equation and outputs the band bending poteatial near a metal-semiconductor interface. Code #2 uses the output of code #1 and calculates the contact resistivity of the interface.

Spataru, Dan↗
DOE OSTI2024

NeuroFEM

SAND2025-00525O NeuroFEM is a software tool that demonstrates a neuromorphic algorithm for solving finite element problems. It sets up a 2D finite element problem for the Poisson equation on a disk, constructs synaptic matrices, and simulates neural dynamics to solve the resulting sparse linear system. The software illustrates how the algorithm converges to the solution and plots the results, showcasing a neuromorphic counterpart to traditional methods like Conjugate Gradient or GMRES. This tool is designed to highlight the potential of neuromorphic algorithms for solving sparse linear systems, which are prevalent in various computational applications. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗
DOE OSTI2023

Sociodemographic effects on pandemic fatigue are multifaceted and context-specific: A longitudinal analysis of physical distancing adherence

Pandemic fatigue emerged early during the COVID-19 pandemic and remains a concern as new variants emerge and ongoing public health measures are needed to control them. A wide range of factors can affect pandemic fatigue, but empiric research indicating which may be most important to adherence in specific populations is lacking. We conducted a longitudinal study of changes in physical distancing in two cohorts: adults living with children <18 years and adults ≥50 years old. Six types of non-work, non-household contacts were ascertained at six times from April to October 2020. We used generalized estimating equations Poisson regression to estimate the one-week change in contact rate and how this differed based on sociodemographic characteristics. The rate of all contact types increased during the middle of the study period and decreased toward the end. Changes in contact rates over time differed according to several sociodemographic characteristics, including age, gender, race/ethnicity, education, household composition, and access to transportation. Furthermore, the factors influencing the rate of change in contact rates differed by the type or setting of the contact, for example contacts as a result of visiting another person’s home versus during a retail outing. These results provide evidence for potential mechanisms by which pandemic fatigue has resulted in lower physical distancing adherence.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Thermodynamic Bounds on Nonlinear Electrostatic Perturbations in Intense Charged Particle Beams

This paper places a lowest upper bound on the field energy in electrostatic perturbations in single-species charged particle beams with initial temperature anisotropy ( Tll T ⊥ < 1). The result applies to all electrostatic perturbations driven by the natural anisotropies that develop in accelerated particle beams, including Harris-type electrostatic instabilities, known to limit the luminosity and minimum spot size attainable in experiments. The thermodynamic bound on the field perturbation energy of the instabilities is obtained from the nonlinear Vlasov-Poisson equations for an arbitrary initial distribution function, including the effects of intense self-fields, finite geometry and nonlinear processes. This paper also includes analytical estimates of the nonlinear bounds for space-charge-dominated and emittance-dominated anisotropic bi-Maxwellian distributions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Global Field Aligned Mesh and Gyrokinetic Field Solver in a Tokamak Edge Geometry (Final Technical Report)

This project was a collaboration between researchers at the California Institute of Technology and the University of California, Irvine to investigate the utility of a global field-aligned mesh and gyrokinetic field solver for simulations of the tokamak plasma edge region. Mesh generation software from UC Irvine was tested with specific tokamak edge magnetic geometry scenarios and the quality of the meshes and the solutions to the gyrokinetic Poisson equation were evaluated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Boltzmann-Electron Model in Aleph.

We apply the Boltzmann-electron model in the electrostatic, particle-in-cell, finite- element code Aleph to a plasma sheath. By assuming a Boltzmann energy distribution for the electrons, the model eliminates the need to resolve the electron plasma fre- quency, and avoids the numerical "grid instability" that can cause unphysical heating of electrons. This allows much larger timesteps to be used than with kinetic electrons. Ions are treated with the standard PIC algorithm. The Boltzmann-electron model re- quires solution of a nonlinear Poisson equation, for which we use an iterative Newton solver (NOX) from the Trilinos Project. Results for the spatial variation of density and voltage in the plasma sheath agree well with an analytic model

Hughes, Thomas Patrick↗
DOE OSTI2017

Final Technical Report on “Center for Nonlinear Simulation of Energetic Particles in Burning Plasmas (CSEP)”

The gyrokinetic turbulence code GEM is used to study reverse shear Alfven eigenmodes in the DIII-D tokamak, cylindrical tearing mode, the nonlinear saturation mechanism of Alfven eigenmodes, and verify other gyrokinetic codes. It is used to study the stability of Alfven waves in the presence of fusion product alpha particles in ITER. New algorithm features and physics capabilities are added to GEM, including a novel discretization method for the Poisson equation, a vorticity approach to the fluid electron model and a low-n solver in cylindrical geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Crosscutting Recombination Metrology for Expediting VOC Engineering

Work carried out under this grant addressed the need for improved materials and device characterization of solar cell designs. This experimental program aimed to develop joint electrical and optical studies of solar cells based on the group-II-VI semiconductor CdTe and the chalcogenide semiconductor Cu(InGa)Se 2 . Electrical measurements focused on temperature (T) and illumination intensity (I) dependence of the open-circuit voltage (VOC), a quantity that plays an important role in conversion efficiency. The VOC depends on temperature and illumination intensity through the energy gap and positions of the quasiFermi levels (QFLs). Optical measurements exploited the illumination and temperature dependence of absolute photoluminescence (PL) intensity (PL-I), a quantity that also depends on the QFL positions. The QFL splitting is a quantity that also depends on recombination rates within the solar cell, including radiative and non-radiative processes. Other factors affect the electrical and optical properties. These include band offsets at both the front and back interfaces of the device active region where electrical contacts are formed. Understanding, locating, and quantifying the recombination properties of the solar cells was the principal goal of this research. Custom simulations were developed to solve continuity and Poisson equations selfconsistently to obtain free-electron and free-hole concentrations. The measurements of temperature and intensity dependent VOC and PL-I allowed the extraction of a number of physical parameters including the conduction band offset, the shunt constant, and the effective capture coefficient for recombination centers, both within the absorber and at the absorber-buffer interface. In addition, the back-contact barrier height is determined from electrical properties. Each of these is discussed in the report. The illumination intensity dependence of VOC and PL-I do not overlap in experiments primarily because PL is not observed in the low-intensity regime where VOC is studied. Absolute PL-I measurements enabled development of an approach directly comparing it with VOC. This in turn allowed development of a simulation capability that would consistently use both measurements, as functions of T and I, for producing a minimization scheme using a differential evolution algorithm. The approach was applied to both CdTe and Cu(InGa)Se 2 solar cells each prepared under various conditions. Shunt conductivity, QFL splitting, recombination were measured. Also extensively studied are current-voltage and frequency-dependent capacitance voltage. Work completed under this grant included studies of magnesium zinc oxide (MZO) buffered CdTe solar cells including the effects of hydroiodic acid (HI) treatment. From these studies the recombination rates were quantified at the buffer/absorber (front) interface. The effect of back-contact barrier height was investigated. CdTe solar cells were prepared with CdSe and CdS buffer layers. Lower grain-boundary potential fluctuation, higher back-contact barrier height, and greater bulk recombination were observed in early devices with CdSe buffer layer. Investigations were also carried out into CdS/Cu(InGa)Se 2 solar cells. A metal-insulatorsemiconductor capacitance-voltage technique was developed to extract recombinationrelevant properties of the absorber interface. From this work the density of interface states D it , their capture cross-section σ, fixed charge at the interface Q f , and surface recombination velocity were all extracted.

14 SOLAR ENERGY↗
DOE OSTI2023

An Infinite Domain 3D Poisson Solver Based on the Barnes-Hut Algorithm

We present a domain decomposition method for the solution of the 3D Poisson equation with infinite domain boundary conditions. The method is based on an application of the Barnes-Hut tree particle scheme adapted to gridded data. Long range interactions are computed using the first two terms in the Cartesian multipole expansion of Green’s function convoluted with the charge while short range computations are performed using Hockney’s domain doubling algorithm. A standard domain decomposition strategy requires O(N 2 ) applications of Hockney’s algorithm, where N is the number of subdomains that intersect that charge support, while in the present approach only O(Nlog 2 N) such computations suffice. The discretization scheme employed is a sixth order Mehrstellen approximation of the 3D Laplace opera tor. The method exhibits satisfactory accuracy at a substantially reduced computational cost compared to the full domain decomposition Hockney’s algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Parallel Simulation of Beam Dynamics in Particle Accelerators [Slides]

Particle accelerators are among the most versatile and important tools of scientific discovery. The Nation's accelerators are responsible for a wealth of advances in materials science, chemistry, the biosciences, particle physics, and nuclear physics. They also have important applications to national security, the environment, energy, medicine, and on the quality of people's lives. LANL has a long history of making pioneering contributions to Accelerator Science including key contributions to the field of Computational Accelerator Physics. These include the development of early beam dynamics codes with space charge (such as PARMILA and PARMELA), the development of rf cavity codes and magnet codes (including Poisson and Superfish), and the development and distribution of codes to the accelerator community through the Los Alamos Accelerator Code Group. LANL researchers also helped pioneer the development of massively parallel space-charge codes. In project t22_accelsim we have moved beyond electrostatic models of collective effects (i.e., solving the Poisson equation in the bunch frame) to fully electromagnetic models based on the Lienard-Wiechert formalism. This approach enables the large-scale simulation of radiation production and collective effects in high brightness electron beams. This is highly relevant to LANL given its future goal of developing an X-ray Free Electron Laser (XFEL). It also directly impacts a LANL LDRD project to develop an undulator-based non-invasive beam profile monitor for beams created in laser-plasma accelerator systems.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Achieving Higher Order Accuracy in Space in Hydrodynamic Simulations of Self-Gravitating Gas

Modern astrophysical simulation codes employ a variety of numerical algorithms capable of achieving higher-order accuracy in both space and time. Albeit they succeed in achieving an effective higher spatial resolution and in suppressing the numerical damping of waves, to our knowledge, all current astrophysical simulations invoking self-gravity are limited to second-order accuracy in space. If we can devise an algorithm to evaluate self-gravity with a higher-order spatial accuracy, we can better the evaluation of the gravitational acceleration and gravitational energy release which dictate the evolution of many astrophysical systems. Herein, we present a numerical algorithm for self-gravitating hydrodynamics capable of achieving fourth order accuracy for a given density distribution on a Cartesian uniform grid. First, we derive the cell-averaged gravitational potential at fourth-order accuracy from the cell-averaged density by solving the Poisson equation. Next, we obtain the cell average of the product of the density and gravitational acceleration, which differs from the cell-averaged density multiplied by the cell-averaged gravitational acceleration. We then show the verification of the algorithm by applying it to critical test problems: (1) maintaining equilibria of self-gravitating slabs, even upon advection, (2) evolving a polytropic sphere with a massive power-law envelope, and (3) conservation of specific entropy during the propagation of a sound wave.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ~1/$\sqrt{N_p}$ with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly nonequilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗