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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Prediction of unsteady transonic flow around missile configurations

This paper describes the preliminary development of a method for predicting the unsteady transonic flow around missiles at transonic and supersonic speeds, with the final goal of developing a computer code for use in aeroelastic calculations or during maneuvers. The basic equations derived for this method are an extension of those derived by Klopfer and Nixon (1989) for steady flow and are a subset of the Euler equations. In this approach, the five Euler equations are reduced to an equation similar to the three-dimensional unsteady potential equation, and a two-dimensional Poisson equation. In addition, one of the equations in this method is almost identical to the potential equation for which there are well tested computer codes, allowing the development of a prediction method based in part on proved technology.

Nixon, D.↗
NASA NTRS1994

On the prediction of velocity fields from redshift space galaxy samples

We present a new method for recovering the underlying velocity field from an observed distribution of galaxies in redshift space. The method is based on a kinematic Zel'dovich relation between the velocity and density fields in redshift space. This relation is expressed in a differential equation slightly modified from the usual Poisson equation and which depends nontrivially on Beta identically equal to Omega(exp 0.6)/b. The linear equation can be readily solved by standard techniques of separation of variables by means of spherical harmonics. One can also include a term describing the 'rocket effect' discussed by Kaiser (1987). From this redshift space information alone, one can generate a prediction of the peculiar velocity field for each harmonic (l, m) as a function of distance. We note that for the quadrupole and higher order moments, the equation is a boundary value problem with solutions dependent on both the interior and exterior mass distribution. However, for a shell at distance r, the dipole, as well as the monopole, of the velocity field in the Local Group frame is fully determined by the interior mass distribution. This implies that the shear of the measured velocity field, when fitted to a dipole distortion, should be aligned and consistent with the gravity field inferred from the well determined local galaxy distribution. As a preliminary application we compute the velocity dipole of distant shells as predicted from the 1.2 Jy IRAS survey compared to the measured velocity dipole on shells, as inferred from a recent POTENT analysis. The coherence between the two fields is good, yielding a best estimate of Beta = 0.6 +/- 0.2.

Nusser, Adi↗
NASA NTRS1995

Fractional Step and Pseudocompressibility Methods for the Solution of Incompressible Navier-Stokes Equations

The fractional step and the pseudocompressibility methods for the solution of the incompressible Navier-Stokes equations are outlined. The fractional step method is based on finite-volume formulation and uses the pressure and the volume fluxes across the faces of each cell as dependent variables. The momentum equations are solved implicitly and the Poisson equation for the pressure is solved by using the multigrid method. The pseudocompressibility approach uses an implicit-higher-order-upwind differencing scheme for the convective terms together with the Gauss-Seidel line relaxation method. The dependent variables in the pseudocompressibility approach are the pressure and the cartesian velocity components in unstaggered mesh orientation. The 90-degree square duct flow, the wing-tip vortex wake flow and unsteady turbulent flows over an oscillating NACA 0015 airfoil are computed using both the fractional step and the pseudocompressibility methods. The results obtained from two different schemes are compared against experimental measurements.

Kiris, Cetin↗
DOE OSTI2017

Screening in a multicomponent plasma by the example of a wet air plasma

A general theory of screening of a dust-particle charge in a multicomponent plasma is developed based on the asymptotic screening theory. In the asymptotic theory, the balance equations for the number of charged plasma particles in the diffusion-drift approximation supplemented with point sinks of plasma particles on a dust particle and the Poisson equation for the self-consistent electric field potential are solved by the linearization method using the three-dimensional Fourier transform. The coefficients of a secular or a characteristic polynomial, whose zeroes determine screening constants in a multicomponent plasma, were calculated by the Leverrier–Faddeev method, allowing one to find these coefficients for any number of initial balance equations for the number of charged plasma particles. By the example of a wet air plasma produced by an external gas-ionization source, the procedure is considered for constructing the truncated system of balance equations for main plasma ions and screening constants are determined. The influence of the number of main ions on screening constants is considered and the physical meaning of all screening constants is found.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Computational modeling and analysis of thermoelectric properties of nanoporous silicon

In this paper, thermoelectric properties of nanoporous silicon are modeled and studied by using a computational approach. The computational approach combines a quantum non-equilibrium Green's function (NEGF) coupled with the Poisson equation for electrical transport analysis, a phonon Boltzmann transport equation (BTE) for phonon thermal transport analysis and the Wiedemann-Franz law for calculating the electronic thermal conductivity. By solving the NEGF/Poisson equations self-consistently using a finite difference method, the electrical conductivity σ and Seebeck coefficient S of the material are numerically computed. The BTE is solved by using a finite volume method to obtain the phonon thermal conductivity k{sub p} and the Wiedemann-Franz law is used to obtain the electronic thermal conductivity k{sub e}. The figure of merit of nanoporous silicon is calculated by ZT=S{sup 2}σT/(k{sub p}+k{sub e}). The effects of doping density, porosity, temperature, and nanopore size on thermoelectric properties of nanoporous silicon are investigated. It is confirmed that nanoporous silicon has significantly higher thermoelectric energy conversion efficiency than its nonporous counterpart. Specifically, this study shows that, with a n-type doping density of 10{sup 20} cm{sup –3}, a porosity of 36% and nanopore size of 3 nm × 3 nm, the figure of merit ZT can reach 0.32 at 600 K. The results also show that the degradation of electrical conductivity of nanoporous Si due to the inclusion of nanopores is compensated by the large reduction in the phonon thermal conductivity and increase of absolute value of the Seebeck coefficient, resulting in a significantly improved ZT.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Dynamic mesoscale model of dipolar fluids via fluctuating hydrodynamics

Fluctuating hydrodynamics (FHD) is a general framework of mesoscopic modeling and simulation based on conservational laws and constitutive equations of linear and nonlinear responses. However, explicit representation of electrical forces in FHD has yet to appear. In this work, we devised an Ansatz for the dynamics of dipole moment densities that is linked with the Poisson equation of the electrical potential ϕ in coupling to the other equations of FHD. The resulting ϕ-FHD equations then serve as a platform for integrating the essential forces, including electrostatics in addition to hydrodynamics, pressure-volume equation of state, surface tension, and solvent-particle interactions that govern the emergent behaviors of molecular systems at an intermediate scale. This unique merit of ϕ-FHD is illustrated by showing that the water dielectric function and ion hydration free energies in homogeneous and heterogenous systems can be captured accurately via the mesoscopic simulation. Furthermore, we show that the field variables of ϕ-FHD can be mapped from the trajectory of an all-atom molecular dynamics simulation such that model development and parametrization can be based on the information obtained at a finer-grained scale. With the aforementioned multiscale capabilities and a spatial resolution as high as 5 Å, the ϕ-FHD equations represent a useful semi-explicit solvent model for the modeling and simulation of complex systems, such as biomolecular machines and nanofluidics.

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

Effect of electron–electron interactions on optical properties of GaN/AlN quantum wells: a nonlinear Schrödinger equation approach

Researchers usually neglect the electron–electron interaction effect when they study the optical properties of semiconducting nanostructures through the compact density matrix approach. In the existing papers, this work has also been done through self-consistent solution of the Schrödinger and Poisson equations. For the first time, we have investigated the electronic and optical properties of wurtzite GaN/AlN quantum wells in the presence of electron–electron interactions assuming weak electron–electron scattering through a numerical solution of a nonlinear Schrödinger equation. It is shown that the electron–electron interaction effect leads to the removal of the N-fold degeneracy in the energy levels of N-well multiple quantum well systems. Effects of different parameters such as electric field, the strength of the electron–electron interaction, a total length of the system and also two different numbers of wells (1 and 2) on the mentioned optical properties have been investigated. Also, we observed redshifts in the optical rectification coefficient peak positions when the nonlinear parameter increased. Finite difference self-consistent method was employed to solve the nonlinear Schrödinger equation and we demonstrated the method convergence. Finally, we illustrated the accuracy of the method by comparing our energy values with those of the Sinc method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Extending cosmological tests of general relativity with the square kilometre array

Tests of general relativity (GR) are still in their infancy on cosmological scales, but forthcoming experiments promise to greatly improve their precision over a wide range of distance scales and redshifts. One such experiment, the Square Kilometre Array (SKA), will carry out several wide and deep surveys of resolved and unresolved neutral hydrogen (H i) 21 cm line-emitting galaxies, mapping a significant fraction of the sky from 0⩽z≲6. I present forecasts for the ability of a suite of possible SKA H i surveys to detect deviations from GR by reconstructing the cosmic expansion and growth history. SKA Phase 1 intensity mapping surveys can achieve sub-1% measurements of fσ{sub 8} out to z≈1, with an SKA1-MID Band 2 survey out to z ≲ 0.6 able to surpass contemporary spectroscopic galaxy surveys such as DESI and Euclid in terms of constraints on modified gravity parameters if challenges such as foreground contamination can be tackled effectively. A more futuristic Phase 2 H i survey of ∼10{sup 9} spectroscopic galaxy redshifts would be capable of detecting a ∼2% modification of the Poisson equation out to z ≈ 2.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Study on longitudinal dispersion relation in one-dimensional relativistic plasma: Linear theory and Vlasov simulation

The dispersion relation of one-dimensional longitudinal plasma waves in relativistic homogeneous plasmas is investigated with both linear theory and Vlasov simulation in this paper. From the Vlasov-Poisson equations, the linear dispersion relation is derived for the proper one-dimensional Jüttner distribution. Numerically obtained linear dispersion relation as well as an approximate formula for plasma wave frequency in the long wavelength limit is given. The dispersion of longitudinal wave is also simulated with a relativistic Vlasov code. The real and imaginary parts of dispersion relation are well studied by varying wave number and plasma temperature. Simulation results are in agreement with established linear theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Operation of a semiconductor microcavity under electric excitation

We present a microscopic theory for the description of the bias-controlled operation of an exciton-polariton-based heterostructure, in particular, the polariton laser. Combining together the Poisson equations for the scalar electric potential and Fermi quasi-energies of electrons and holes in a semiconductor heterostructure, the Boltzmann equation for the incoherent excitonic reservoir and the Gross-Pitaevskii equation for the exciton-polariton mean field, we simulate the dynamics of the system minimising the number of free parameters and build a theoretical threshold characteristic: number of particles vs applied bias. This approach, which also accounts for the nonlinear (exciton-exciton) interaction, particle lifetime, and which can, in principle, account for any relaxation mechanisms for the carriers of charge inside the heterostructure or polariton loss, allows to completely describe modern experiments on polariton transport and model devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Cosmological constraints on Lorentz violating dark energy

The role of Lorentz invariance as a fundamental symmetry of nature has been lately reconsidered in different approaches to quantum gravity. It is thus natural to study whether other puzzles of physics may be solved within these proposals. This may be the case for the cosmological constant problem. Indeed, it has been shown that breaking Lorentz invariance provides Lagrangians that can drive the current acceleration of the universe without experiencing large corrections from ultraviolet physics. In this work, we focus on the simplest model of this type, called ΘCDM, and study its cosmological implications in detail. At the background level, this model cannot be distinguished from ΛCDM. The differences appear at the level of perturbations. We show that in ΘCDM, the spectrum of CMB anisotropies and matter fluctuations may be affected by a rescaling of the gravitational constant in the Poisson equation, by the presence of extra contributions to the anisotropic stress, and finally by the existence of extra clustering degrees of freedom. To explore these modifications accurately, we modify the Boltzmann code class. We then use the parameter inference code Monte Python to confront ΘCDM with data from WMAP-7, SPT and WiggleZ. We obtain strong bounds on the parameters accounting for deviations from ΛCDM. In particular, we find that the discrepancy between the gravitational constants appearing in the Poisson and Friedmann equations is constrained at the level of 1.8%.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Shielding effect and wakefield pattern of a moving test charge in a non-Maxwellian dusty plasma

By using the Vlasov-Poisson equations, we calculate an expression for the electrostatic potential caused by a test charge in an unmagnetized non-Maxwellian dusty plasma, whose constituents are the superthermal hot-electrons, the mobile cold-electrons with a neutralizing background of cold ions, and charge fluctuating isolated dust grains. The superthermality effects due to hot electrons not only modify the dielectric constant of the electron-acoustic waves but also significantly affect the electrostatic potential. The latter can be decomposed into the Debye-Hückel and oscillatory wake potentials. Analytical and numerical results reveal that the Debye-Hückel and wakefield potentials converge to the Maxwellian case for large values of superthermality parameter. Furthermore, the plasma parameters play a vital role in the formation of shielding and wakefield pattern in a two-electron temperature plasma. The present results should be important for laboratory and space dusty plasmas, where hot-electrons can be assumed to follow the non-Maxwellian distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Kinetic study of ion acoustic twisted waves with kappa distributed electrons

The kinetic theory of Landau damping of ion acoustic twisted modes is developed in the presence of orbital angular momentum of the helical (twisted) electric field in plasmas with kappa distributed electrons and Maxwellian ions. The perturbed distribution function and helical electric field are considered to be decomposed by Laguerre-Gaussian mode function defined in cylindrical geometry. The Vlasov-Poisson equation is obtained and solved analytically to obtain the weak damping rates of the ion acoustic twisted waves in a non-thermal plasma. The strong damping effects of ion acoustic twisted waves at low values of temperature ratio of electrons and ions are also obtained by using exact numerical method and illustrated graphically, where the weak damping wave theory fails to explain the phenomenon properly. The obtained results of Landau damping rates of the twisted ion acoustic wave are discussed at different values of azimuthal wave number and non-thermal parameter kappa for electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Twisted Landau damping rates in multi-component dusty plasmas

Keeping in view the kinetic treatment for plasma particles, the electrostatic twisted dust-acoustic (DA) and dust-ion-acoustic (DIA) waves are investigated in a collisionless unmagnetized multi-component dusty plasma, whose constituents are the electrons, singly ionized positive ions, and negatively charged massive dust particulates. With this background, the Vlasov–Poisson equations are coupled together to derive a generalized dielectric constant by utilizing the Laguerre-Gaussian perturbed distribution function and electrostatic potential in the paraxial limit. The dispersion and damping rates of twisted DA and DIA waves are analyzed with finite orbital angular momentum states in a multi-component dusty plasma. Significant modifications concerning the real wave frequencies and damping rates appeared with varying twisted dimensionless parameter and dust concentration. In particular, it is shown that dust concentration enhances the phase speed of the DIA waves in contrary to DA waves, whereas the impact of twisted parameter reduces the frequencies of both DA and DIA waves. The results should be useful for the understanding of particle transport and trapping phenomena caused by wave excitation in laboratory dusty plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Maxwell-Garnett effective medium theory: Quantum nonlocal effects

We develop the Maxwell-Garnett theory for the effective medium approximation of composite materials with metallic nanoparticles by taking into account the quantum spatial dispersion effects in dielectric response of nanoparticles. We derive a quantum nonlocal generalization of the standard Maxwell-Garnett formula, by means the linearized quantum hydrodynamic theory in conjunction with the Poisson equation as well as the appropriate additional quantum boundary conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1988

Three-dimensional elliptic grid generation for an F-16

A case history depicting the effort to generate a computational grid for the simulation of transonic flow about an F-16 aircraft at realistic flight conditions is presented. The flow solver for which this grid is designed is a zonal one, using the Reynolds averaged Navier-Stokes equations near the surface of the aircraft, and the Euler equations in regions removed from the aircraft. A body conforming global grid, suitable for the Euler equation, is first generated using 3-D Poisson equations having inhomogeneous terms modeled after the 2-D GRAPE code. Regions of the global grid are then designated for zonal refinement as appropriate to accurately model the flow physics. Grid spacing suitable for solution of the Navier-Stokes equations is generated in the refinement zones by simple subdivision of the given coarse grid intervals. That grid generation project is described, with particular emphasis on the global coarse grid.

Reese L. Sorenson↗
NASA NTRS1991

Convection in vertical Bridgman configurations

We are interested in tracing the convective profiles of vertical Bridgman growth in two configurations, the pure Rayleigh convective mode and the combined Rayleigh-Marangoni mode. In order to do so, we conducted a numerical investigation that involved a finite volume calculation. The governing equations were integrated about a cell volume, using the Gauss Theorem and the volume variables like temperature and velocity were related to the surface variables. In order to solve for the pressure field, we employed the continuity equation and the residuals resulted in a Poisson equation. Results and comments for the Rayleigh and Marangoni problems in a vertical cylinder or Bridgman configuration are given.

Narayanan, Ranga↗