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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Sound Velocities in Vanadium Reveal Complex Elastic Behavior at High Pressures

Compressional (VP) and shear (VS) wave velocities of polycrystalline vanadium were measured simultaneously up to 11.5 GPa at room temperature using ultrasonic interferometry in a multi-anvil press. Complex softening behavior in VS and resulting shear moduli are discovered, possibly revealing a precursor to the reported phase transition within 30–60 GPa. The current data enables a comprehensive assessment of the elastic and mechanical properties of vanadium at high pressures, including bulk and shear moduli, Young’s modulus, Poisson’s ratio, and Pugh’s ratio. Through fitting to the 3rd-order finite strain equations, the elastic moduli and their pressure derivatives were determined to be K S0 = 151 (2) GPa, G 0 = 46.9 (8) GPa, K$^{’}_{S0}$ = 3.47 (5), and G$^{’}_{0}$ = 0.62 (1). These experimental results allow us to compare with and benchmark the existing Steinberg–Guinan models for extrapolations to extreme pressure and temperature conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Theoretical investigation on improving the ductility of Rh{sub 3}V by ternary addition

The band structure calculations are performed on the intermetallic compounds Rh{sub 3}V and Rh{sub 3}V{sub x}(Al){sub 1-x} by Full Potential Linearized Augmented Plane Wave(FP-LAPW) method. Total energies are calculated as a function of volume and fitted to Birch-Murnaughan equation of state to find the lattice parameter and the other ground state properties. The ductility of these compounds has been analyzed using Cauchy’s Pressure, Pugh Rule and Poisson’s ratio. From this study we observed that Rh{sub 3}V is brittle while its alloy Rh{sub 3}V{sub 0.875}Al{sub 0.125} is ductile. Electron density plots are illustrated and compared.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Crossing the dividing surface of transition state theory. III. Once and only once. Selecting reactive trajectories

The purpose of the present work is to determine initial conditions that generate reacting, recrossing-free trajectories that cross the conventional dividing surface of transition state theory (i.e., the plane in configuration space passing through a saddle point of the potential energy surface and perpendicular to the reaction coordinate) without ever returning to it. Local analytical equations of motion valid in the neighborhood of this planar surface have been derived as an expansion in Poisson brackets. We show that the mere presence of a saddle point implies that reactivity criteria can be quite simply formulated in terms of elements of this series, irrespective of the shape of the potential energy function. Some of these elements are demonstrated to be equal to a sum of squares and thus to be necessarily positive, which has a profound impact on the dynamics. The method is then applied to a three-dimensional model describing an atom-diatom interaction. A particular relation between initial conditions is shown to generate a bundle of reactive trajectories that form reactive cylinders (or conduits) in phase space. This relation considerably reduces the phase space volume of initial conditions that generate recrossing-free trajectories. Loci in phase space of reactive initial conditions are presented. Reactivity is influenced by symmetry, as shown by a comparative study of collinear and bent transition states. Finally, it is argued that the rules that have been derived to generate reactive trajectories in classical mechanics are also useful to build up a reactive wave packet.

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

Helically symmetric extended magnetohydrodynamics: Hamiltonian formulation and equilibrium variational principles

Hamiltonian extended magnetohydrodynamics (XMHD) is restricted to respect helical symmetry by reducing the Poisson bracket for the three-dimensional dynamics to a helically symmetric one, as an extension of the previous study for translationally symmetric XMHD (Kaltsas et al. , Phys. Plasmas , vol. 24, 2017, 092504). Four families of Casimir invariants are obtained directly from the symmetric Poisson bracket and they are used to construct Energy–Casimir variational principles for deriving generalized XMHD equilibrium equations with arbitrary macroscopic flows. The system is then cast into the form of Grad–Shafranov–Bernoulli equilibrium equations. The axisymmetric and the translationally symmetric formulations can be retrieved as geometric reductions of the helically symmetric one. As special cases, the derivation of the corresponding equilibrium equations for incompressible plasmas is discussed and the helically symmetric equilibrium equations for the Hall MHD system are obtained upon neglecting electron inertia. An example of an incompressible double-Beltrami equilibrium is presented in connection with a magnetic configuration having non-planar helical magnetic axis.

Physics↗
NASA NTRS1988

An exactly solvable model for calculating critical misfit and thickness in epitaxial superlattices - Layers of equal elastic constants and thicknesses

A parabolic interaction potential has been used to develop a model for calculating the misfit dislocation (MD) energy in the case of a superlattice of alternating layers of materials with equal elastic constants and thicknesses. The model, which is believed to be a good one for small misfits and to have some merit for covalent bonded materials, is exactly solvable for the critical thickness above which it is energetically favorable to lose coherency by the introduction of MDs into the interfaces. It was found, for a given misfit f, that the critical thickness for epitaxial superlattices free from their substrate is somewhat more than four times that for a single epilayer on a thick substrate. Furthermore, the critical thickness varies almost inversely with misfit to the power 1.22 when Poisson's ratio is 1/3. It was also shown that the critical misfit f(c) obtained by equating maximal misfit strain and MD energies is a significant overestimate of f(c). The results for a superlattice are compared with those of a thin layer on a thick substrate.

Van Der Merwe, Jan H.↗
NASA NTRS2008

Application of an Elongated Kelvin Model to Space Shuttle Foams

Spray-on foam insulation is applied to the exterior of the Space Shuttle s External Tank to limit propellant boil-off and to prevent ice formation. The Space Shuttle foams are rigid closed-cell polyurethane foams. The two foams used most extensively on the Space Shuttle External Tank are BX-265 and NCFI24-124. Since the catastrophic loss of the Space Shuttle Columbia, numerous studies have been conducted to mitigate the likelihood and the severity of foam shedding during the Shuttle s ascent to space. Due to the foaming and rising process, the foam microstructures are elongated in the rise direction. As a result, these two foams exhibit a non-isotropic mechanical behavior. In this paper, a detailed microstructural characterization of the two foams is presented. The key features of the foam cells are summarized and the average cell dimensions in the two foams are compared. Experimental studies to measure the room temperature mechanical response of the two foams in the two principal material directions (parallel to the rise and perpendicular to the rise) are also reported. The measured elastic modulus, proportional limit stress, ultimate tensile stress and the Poisson s ratios for the two foams are compared. The generalized elongated Kelvin foam model previously developed by the authors is reviewed and the equations which result from this model are presented. The resulting equations show that the ratio of the elastic modulus in the rise direction to that in the perpendicular-to-rise direction as well as the ratio of the strengths in the two material directions is only a function of the microstructural dimensions. Using the measured microstructural dimensions and the measured stiffness ratio, the foam tensile strength ratio and Poisson s ratios are predicted for both foams. The predicted tensile strength ratio is in close agreement with the measured strength ratios for both BX-265 and NCFI24-124. The comparison between the predicted Poisson s ratios and the measured values is not as favorable.

Sullivan, Roy M.↗
DOE OSTI2017

Canonical symplectic structure and structure-preserving geometric algorithms for Schrödinger–Maxwell systems

An infinite dimensional canonical symplectic structure and structure-preserving geometric algorithms are developed for the photon–matter interactions described by the Schrödinger–Maxwell equations. The algorithms preserve the symplectic structure of the system and the unitary nature of the wavefunctions, and bound the energy error of the simulation for all time-steps. This new numerical capability enables us to carry out first-principle based simulation study of important photon–matter interactions, such as the high harmonic generation and stabilization of ionization, with long-term accuracy and fidelity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Accurate derivative evaluation for any Grad–Shafranov solver

We present a numerical scheme that can be combined with any fixed boundary finite element based Poisson or Grad–Shafranov solver to compute the first and second partial derivatives of the solution to these equations with the same order of convergence as the solution itself. At the heart of our scheme is an efficient and accurate computation of the Dirichlet to Neumann map through the evaluation of a singular volume integral and the solution to a Fredholm integral equation of the second kind. Our numerical method is particularly useful for magnetic confinement fusion simulations, since it allows the evaluation of quantities such as the magnetic field, the parallel current density and the magnetic curvature with much higher accuracy than has been previously feasible on the affordable coarse grids that are usually implemented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Stochastic simulation of reaction-diffusion systems: A fluctuating-hydrodynamics approach

Here, we develop numerical methods for stochastic reaction-diffusion systems based on approaches used for fluctuating hydrodynamics (FHD). For hydrodynamic systems, the FHD formulation is formally described by stochastic partial differential equations (SPDEs). In the reaction-diffusion systems we consider, our model becomes similar to the reaction-diffusion master equation (RDME) description when our SPDEs are spatially discretized and reactions are modeled as a source term having Poisson fluctuations. However, unlike the RDME, which becomes prohibitively expensive for an increasing number of molecules, our FHD-based description naturally extends from the regime where fluctuations are strong, i.e., each mesoscopic cell has few (reactive) molecules, to regimes with moderate or weak fluctuations, and ultimately to the deterministic limit. By treating diffusion implicitly, we avoid the severe restriction on time step size that limits all methods based on explicit treatments of diffusion and construct numerical methods that are more efficient than RDME methods, without compromising accuracy. Guided by an analysis of the accuracy of the distribution of steady-state fluctuations for the linearized reaction-diffusion model, we construct several two-stage (predictor-corrector) schemes, where diffusion is treated using a stochastic Crank-Nicolson method, and reactions are handled by the stochastic simulation algorithm of Gillespie or a weakly second-order tau leaping method. We find that an implicit midpoint tau leaping scheme attains second-order weak accuracy in the linearized setting and gives an accurate and stable structure factor for a time step size of an order of magnitude larger than the hopping time scale of diffusing molecules. We study the numerical accuracy of our methods for the Schlögl reaction-diffusion model both in and out of thermodynamic equilibrium. We demonstrate and quantify the importance of thermodynamic fluctuations to the formation of a two-dimensional Turing-like pattern and examine the effect of fluctuations on three-dimensional chemical front propagation. Furthermore, by comparing stochastic simulations to deterministic reaction-diffusion simulations, we show that fluctuations accelerate pattern formation in spatially homogeneous systems and lead to a qualitatively different disordered pattern behind a traveling wave.

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

Conservative fourth-order finite-volume Vlasov–Poisson solver for axisymmetric plasmas in cylindrical (r,v r ,v θ ) phase space coordinates

In this paper, a fourth-order finite-volume Vlasov–Poisson algorithm is developed for simulating axisymmetric plasma configurations in ( r , v r , v θ ) phase space coordinates. The Vlasov equation for cylindrical phase space coordinates is cast into conservation-law form and is discretized on a structured grid. The conservative finite-volume discretization is based on fifth-order upwind reconstructions of the distribution function and a fourth-order quadrature rule that accounts for transverse variations of fluxes along control volume surfaces. High-order specular reflection boundary conditions enable high-fidelity treatment of plasma distribution functions at the axis and at wall boundaries. The numerical method is applied to simulate a confined uniform neutral gas to assess convergence properties for an equilibrium system in which effects of finite-temperature and acceleration due to centrifugal and Coriolis forces are present. The discretization is also applied to a Z-pinch configuration to study electrostatic ion confinement and the dynamics of the associated breathing mode. Simulations show that the ion distribution function exhibits non-Maxwellian features and that a temperature anisotropy develops and is sustained. Finally, the finite-volume implementation is demonstrated to converge at fourth-order for both applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1999

A High Order Discontinuous Galerkin Method for 2D Incompressible Flows

In this paper we introduce a high order discontinuous Galerkin method for two dimensional incompressible flow in vorticity streamfunction formulation. The momentum equation is treated explicitly, utilizing the efficiency of the discontinuous Galerkin method The streamfunction is obtained by a standard Poisson solver using continuous finite elements. There is a natural matching between these two finite element spaces, since the normal component of the velocity field is continuous across element boundaries. This allows for a correct upwinding gluing in the discontinuous Galerkin framework, while still maintaining total energy conservation with no numerical dissipation and total enstrophy stability The method is suitable for inviscid or high Reynolds number flows. Optimal error estimates are proven and verified by numerical experiments.

Liu, Jia-Guo↗
NASA NTRS1989

A dual potential formulation of the Navier-Stokes equations

A dual potential formulation for numerically solving the Navier-Stokes equations is developed and presented. The velocity field is decomposed using a scalar and vector potential. Vorticity and dilatation are used as the dependent variables in the momentum equations. Test cases in two dimensions verify the capability to solve flows using approximations from potential flow to full Navier-Stokes simulations. A three-dimensional incompressible flow formulation is also described. An interesting feature of this approach to solving the Navier-Stokes equations is the decomposition of the velocity field into a rotational part (vector potential) and an irrotational part (scalar potential). The Helmholtz decomposition theorem allows this splitting of the velocity field. This approach has had only limited use since it increases the number of dependent variables in the solution. However, it has often been used for incompressible flows where the solution scheme is known to be fast and accurate. This research extends the usage of this method to fully compressible Navier-Stokes simulations by using the dilatation variable along with vorticity. A time-accurate, iterative algorithm is used for the uncoupled solution of the governing equations. Several levels of flow approximation are available within the framework of this method. Potential flow, Euler and full Navier-Stokes solutions are possible using the dual potential formulation. Solution efficiency can be enhanced in a straightforward way. For some flows, the vorticity and/or dilatation may be negligible in certain regions (e.g., far from a viscous boundary in an external flow). It is possible to drop the calculation of these variables then and optimize the solution speed. Also, efficient Poisson solvers are available for the potentials. The relative merits of non-primitive variables versus primitive variables for solution of the Navier-Stokes equations are also discussed.

Gegg, S. G.↗
DOE OSTI2020

Uncertainty-based weighted least squares density integration for background-oriented schlieren

We propose an improved density integration methodology for Background-Oriented Schlieren (BOS) measurements that overcomes the noise sensitivity of the commonly used Poisson solver. Here, the method employs a weighted least-squares (WLS) optimization of the 2D integration of the density gradient field by solving an over-determined system of equations. Weights are assigned to the grid points based on density gradient uncertainties to ensure that a less reliable measurement point has less effect on the integration procedure. Synthetic image analysis with a Gaussian density field shows that WLS constrains the propagation of random error and reduces it by 80% in comparison to Poisson for the highest noise level. Using WLS with experimental BOS measurements of flow induced by a spark plasma discharge shows a 30% reduction in density uncertainty in comparison to Poisson, thereby increasing the overall precision of the BOS density measurements.

42 ENGINEERING↗
DOE OSTI2016

The classical limit of minimal length uncertainty relation: revisit with the Hamilton-Jacobi method

The existence of a minimum measurable length could deform not only the standard quantum mechanics but also classical physics. The effects of the minimal length on classical orbits of particles in a gravitation field have been investigated before, using the deformed Poisson bracket or Schwarzschild metric. In this paper, we first use the Hamilton-Jacobi method to derive the deformed equations of motion in the context of Newtonian mechanics and general relativity. We then employ them to study the precession of planetary orbits, deflection of light, and time delay in radar propagation. We also set limits on the deformation parameter by comparing our results with the observational measurements. Finally, comparison with results from previous papers is given at the end of this paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Cavitation in a metallic liquid: Homogeneous nucleation and growth of nanovoids

Large-scale molecular dynamics (MD) simulations are performed to investigate homogeneous nucleation and growth of nanovoids during cavitation in liquid Cu. We characterize in detail the atomistic cavitation processes by following the temporal evolution of cavities or voids, analyze the nucleation behavior with the mean first-passage time (MFPT) and survival probability (SP) methods, and discuss the results against classical nucleation theory (CNT), the Tolman equation for surface energy, independent calculation of surface tension via integrating the stress profiles, the Johnson-Mehl-Avrami (JMA) growth law, and the power law for nucleus size distributions. Cavitation in this representative metallic liquid is a high energy barrier Poisson processes, and the steady-state nucleation rates obtained from statistical runs with the MFPT and SP methods are in agreement. The MFPT method also yields the critical nucleus size and the Zeldovich factor. Fitting with the Tolman's equation to the MD simulations yields the surface energy of a planar interface (∼0.9 J m{sup −2}) and the Tolman length (0.4–0.5 Å), and those values are in accord with those from integrating the stress profiles of a planar interface. Independent CNT predictions of the nucleation rate (10{sup 33−34} s{sup −1} m{sup −3}) and critical size (3–4 Å in radius) are in agreement with the MFPT and SP results. The JMA law can reasonably describe the nucleation and growth process. The size distribution of subcritical nuclei appears to follow a power law with an exponent decreasing with increasing tension owing to coupled nucleation and growth, and that of the supercritical nuclei becomes flattened during further stress relaxation due to void coalescence.

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

SU(2s+1) symmetry and nonlinear dynamics of high spin magnets

The article is devoted to the description of dynamics of magnets with arbitrary spin on the basis of the Hamiltonian formalism. The relationship of quantum states and magnetic degrees of freedom has been considered. Subalgebras of Poisson bracket of magnetic values for spin s=1/2; 1; 3/2 have been established. We have obtained non-linear dynamic equations for the normal and degenerate non-equilibrium states of high-spin magnets with the SO(3), SU(4), SU(2)×SU(2), SU(3), SO(4), SO(5) symmetries of exchange interaction. The connection between models of magnetic exchange energy and the Casimir invariants has been discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

MechElastic: A Python library for analysis of mechanical and elastic properties of bulk and 2D materials

We report the MechElastic Python package evaluates the mechanical and elastic properties of bulk and 2D materials using the elastic coefficient matrix ( C ij ) obtained from any ab-initio density-functional theory (DFT) code. The current version of this package reads the output of VASP, ABINIT, and Quantum Espresso codes (but it can be easily generalized to any other DFT code) and performs the appropriate post-processing of elastic constants as per the requirement of the user. This program can also detect the input structure's crystal symmetry and test the mechanical stability of all crystal classes using the Born-Huang criteria. Various useful material-specific properties such as elastic moduli, longitudinal and transverse elastic wave velocities, Debye temperature, elastic anisotropy, 2D layer modulus, hardness, Pugh's ratio, Cauchy's pressure, Kleinman's parameter, and Lame's coefficients, can be estimated using this program. Another existing feature of this program is to employ the ELATE package (2016) [29] and plot the spatial variation of several elastic properties such as Poisson's ratio, linear compressibility, shear modulus, and Young's modulus in three dimensions. Further, the MechElastic package can plot the equation of state (EOS) curves for energy and pressure for a variety of EOS models such as Murnaghan, Birch, Birch-Murnaghan, and Vinet, by reading the inputted energy/pressure versus volume data obtained via numerical calculations or experiments. This package is particularly useful for the high-throughput analysis of elastic and mechanical properties of materials.

2D materials↗
DOE OSTI2024

Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds

Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.

97 MATHEMATICS AND COMPUTING↗