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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

An improved plate theory of order (1,2) for thick composite laminates

A new (1,2)-order theory is proposed for the linear elasto-static analysis of laminated composite plates. The basic assumptions are those concerning the distribution through the laminate thickness of the displacements, transverse shear strains and the transverse normal stress, with these quantities regarded as some weighted averages of their exact elasticity theory representations. The displacement expansions are linear for the inplane components and quadratic for the transverse component, whereas the transverse shear strains and transverse normal stress are respectively quadratic and cubic through the thickness. The main distinguishing feature of the theory is that all strain and stress components are expressed in terms of the assumed displacements prior to the application of a variational principle. This is accomplished by an a priori least-square compatibility requirement for the transverse strains and by requiring exact stress boundary conditions at the top and bottom plate surfaces. Equations of equilibrium and associated Poisson boundary conditions are derived from the virtual work principle. It is shown that the theory is particularly suited for finite element discretization as it requires simple C(sup 0)- and C(sup -1)-continuous displacement interpolation fields. Analytic solutions for the problem of cylindrical bending are derived and compared with the exact elasticity solutions and those of our earlier (1,2)-order theory based on the assumed displacements and transverse strains.

Tessler, A.↗
NASA NTRS2000

Parallel Unsteady Turbopump Flow Simulations for Reusable Launch Vehicles

An efficient solution procedure for time-accurate solutions of Incompressible Navier-Stokes equation is obtained. Artificial compressibility method requires a fast convergence scheme. Pressure projection method is efficient when small time-step is required. The number of sub-iteration is reduced significantly when Poisson solver employed with the continuity equation. Both computing time and memory usage are reduced (at least 3 times). Other work includes Multi Level Parallelism (MLP) of INS3D, overset connectivity for the validation case, experimental measurements, and computational model for boost pump.

Kiris, Cetin↗
NASA NTRS2004

Computations of Wall Distances Based on Differential Equations

The use of differential equations such as Eikonal, Hamilton-Jacobi and Poisson for the economical calculation of the nearest wall distance d, which is needed by some turbulence models, is explored. Modifications that could palliate some turbulence-modeling anomalies are also discussed. Economy is of especial value for deforming/adaptive grid problems. For these, ideally, d is repeatedly computed. It is shown that the Eikonal and Hamilton-Jacobi equations can be easy to implement when written in implicit (or iterated) advection and advection-diffusion equation analogous forms, respectively. These, like the Poisson Laplacian term, are commonly occurring in CFD solvers, allowing the re-use of efficient algorithms and code components. The use of the NASA CFL3D CFD program to solve the implicit Eikonal and Hamilton-Jacobi equations is explored. The re-formulated d equations are easy to implement, and are found to have robust convergence. For accurate Eikonal solutions, upwind metric differences are required. The Poisson approach is also found effective, and easiest to implement. Modified distances are not found to affect global outputs such as lift and drag significantly, at least in common situations such as airfoil flows.

Tucker, Paul G.↗
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 OSTI2013

Lifting particle coordinate changes of magnetic moment type to Vlasov-Maxwell Hamiltonian dynamics

Techniques for coordinate changes that depend on both dependent and independent variables are developed and applied to the Maxwell-Vlasov Hamiltonian theory. Particle coordinate changes with a new velocity variable dependent on the magnetic field, with spatial coordinates unchanged, are lifted to the field theoretic level, by transforming the noncanonical Poisson bracket and Hamiltonian structure of the Vlasov-Maxwell dynamics. Several examples are given including magnetic coordinates, where the velocity is decomposed into components parallel and perpendicular to the local magnetic field, and the case of spherical velocity coordinates. An example of the lifting procedure is performed to obtain a simplified version of gyrokinetics, where the magnetic moment is used as a coordinate and the dynamics is reduced by elimination of the electric field energy in the Hamiltonian.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2018

Acoustic properties and density of polyurea at pressure up to 13.5 GPa through Brillouin scattering spectroscopy

Brillouin scattering was performed on an elastomeric polyurea, a significant technological polymer. Being widely used for impact modification, of particular interest is its response to extreme pressure conditions. We applied pressures up to 13.5 GPa utilizing a diamond anvil cell and measured the longitudinal and transverse sound velocities via Brillouin light scattering. From these data, the equation of state, the elastic moduli, and Poisson's ratio were obtained. By comparison with previous dilatometry measurements up to 1 GPa, we show how viscoelastic effects can be accounted for in order to obtain an accurate equation of state. Due to the extreme strain-rate hardening of vitrifying polyurea, the property changes associated with its solidification are more subtle in the high frequency Brillouin data than observed in conventional mechanical testing and dilatometry.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Existence, uniqueness, and stability of stochastic neutral functional differential equations of Sobolev-type

In this paper, we are mainly concerned with a class of stochastic neutral functional differential equations of Sobolev-type with Poisson jumps. Under two different sets of conditions, we establish the existence of the mild solution by applying the Leray-Schauder alternative theory and the Sadakovskii’s fixed point theorem, respectively. Furthermore, we use the Bihari’s inequality to prove the Osgood type uniqueness. Also, the mean square exponential stability is investigated by applying the Gronwall inequality. Finally, two examples are given to illustrate the theory results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1984

Generic Karman-Rostovstev plate equations in an affine space

The Karman-Rostovstev plate equations are affinely transformed into equations that depend on the generic plate constants generalized rigidity, generalized extensional rigidity, and generalized Poisson ratio associated with the well-known linearized orthotropic plate/slab equations. Furthermore, since the generalized extensional rigidity is a function of generalized rigidity and generalized Poisson ratio only, with generalized Poisson ratio being a weak parameter, the only fundamental constant is generalized rigidity. This means that regardless of what analysis technique is used, the results can be cast in terms of only one strong material parameter (whose range is from zero to one). The affine stretching is performed on the dependent variables as well as the independent variables, so that the dependence of all quantities of interest on generalized rigidity and generalized Poisson ratio are easily calculated (e.g., even auxiliary quantities such as stress, strain, all weighted stress resultants, etc.). Some typical results (static side load response, collapse load limited by rib buckling) are presented by modifying an approximate technique due to Donnell (1976).

Brunelle, E. J.↗
DOE OSTI2012

Relativistic current collection by a cylindrical Langmuir probe

The current I to a cylindrical Langmuir probe with a bias {Phi}{sub p} satisfying {beta}{identical_to}e{Phi}{sub p}/m{sub e}c{sup 2}{approx}O(1) is discussed. The probe is considered at rest in an unmagnetized plasma composed of electrons and ions with temperatures kT{sub e}{approx}kT{sub i} Much-Less-Than m{sub e}c{sup 2}. For small enough radius, the probe collects the relativistic orbital-motion-limited (OML) current I{sub OML}, which is shown to be larger than the non-relativistic result; the OML current is proportional to {beta}{sup 1/2} and {beta}{sup 3/2} in the limits {beta} Much-Less-Than 1 and {beta} Much-Greater-Than 1, respectively. Unlike the non-relativistic case, the electron density can exceed the unperturbed density value. An asymptotic theory allowed to compute the maximum radius of the probe to collect OML current, the sheath radius for probe radius well below maximum and how the ratio I/I{sub OML} drops below unity when the maximum radius is exceeded. A numerical algorithm that solves the Vlasov-Poisson system was implemented and density and potential profiles presented. The results and their implications in a possible mission to Jupiter with electrodynamic bare tethers are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

The expansion of a collisionless plasma into a plasma of lower density

This paper considers the asymptotic and numerical solution of a simple model for the expansion of a collisionless plasma into a plasma of lower density. The dependence on the density ratio of qualitative and quantitative features of solutions of the well-known cold-ion model is explored. In the cold-ion limit, we find that a singularity develops in the ion density in finite time unless the density ratio is zero or close to unity. The classical cold-ion model may cease to be valid when such a singularity occurs and we then regularize the model by the finite ion-temperature Vlasov-Poisson system. Numerical evidence suggests the emergence of a multi-modal velocity distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1974

Special perturbations employing osculating reference states

The concept of employing osculating reference position and velocity vectors in the numerical integration of the equations of motion of a satellite is examined. The choice of the reference point is shown to have a significant effect upon numerical efficiency and the class of trajectories described by the differential equations of motion. For example, when the position and velocity vectors on the osculating orbit at a fixed reference time are chosen, a universal formulation is yielded. For elliptical orbits, however, this formulation is unattractive for numerical integration purposes due to Poisson terms (mixed secular) appearing in the equations of motion. Other choices for the reference point eliminate this problem but usually at the expense of universality. A number of these formulations, including a universal one, are considered here. Comparisons of the numerical characteristics of these techniques with those of the Encke method are presented.

Born, G. H.↗
DOE OSTI2012

Influence of a strong laser field on Coulomb explosion and stopping power of energetic H{sub 3}{sup +} clusters in plasmas

The influence of a high-intensity laser field on the Coulomb explosion and stopping power for a swift H{sub 3}{sup +} cluster ion in a plasma target is studied by means of the molecular dynamic (MD) method based on the linearized Vlasov-Poisson theory. Excitations of the plasma are described by the classical plasma dielectric function. In the presence of the laser field, the general expressions for the induced potential in the target and the interaction force among the ions within the cluster are derived. Based on the numerical solution of the equations of motion for the constituent ions, the Coulomb explosion patterns and the cluster's stopping power are discussed for a range of laser parameters. Numerical results show that the laser field affects the correlation between the ions and contributes to weaken the wake effect and the stopping power as compared to the laser-free case. On the other hand, the stopping power ratio of H{sub 3}{sup +} cluster is higher than the situation of dicluster of H{sub 2}{sup +} due to the vicinage effect in the cluster.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Observation of threshold energy and hysteresis in high current ion beam guiding and transmission through a micro-glass-capillary

The energy dependent guiding of high current density (0.1-3 A/m{sup 2}) argon ion beams through a micro-glass-capillary is studied. It is observed that ion transmission through the capillary takes place only if its energy is greater than the retarding potential barrier, which depends upon the amount of charge deposited on the capillary inner wall. Foremost evidence of the observation that the transmitted current exhibits hysteresis with ion energy is presented. Particle in cell simulations carried out by solving Poisson's and Newton's force equation self-consistently, agree reasonably well with the experimental results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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)↗
NASA NTRS2016

Viscoplastic Characterization of Ti-6-4: Experiments

As part of a continued effort to improve the understanding of material time-dependent response, a series of mechanical tests have been conducted on the titanium alloy, Ti-6Al-4V. Tensile, creep, and stress relaxation tests were performed over a wide range of temperatures and strain rates to engage various amounts of time-dependent behavior. Additional tests were conducted that involved loading steps, overloads, dwell periods, and block loading segments to characterize the interaction between plasticity and time-dependent behavior. These data will be used to characterize a recently developed, viscoelastoplastic constitutive model with a goal toward better estimates of aerospace component behavior, resulting in improved safety.

mechanical properties↗