Search NASA⌕ Search

SEARCH · Search NASA

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Nonlinear hydrodynamics of cosmological sheets. 1: Numerical techniques and tests

We present the numerical techniques and tests used to construct and validate a computer code designed to study the multidimensional nonlinear hydrodynamics of large-scale sheet structures in the universe, especially the fragmentation of such structures under various instabilities. This code is composed of two codes, the hydrodynamical code ZEUS-2D and a particle-mesh code. The ZEUS-2D code solves the hydrodynamical equations in two dimensions using explicit Eulerian finite-difference techniques, with modifications made to incorporate the expansion of the universe and the gas cooling due to Compton scattering, bremsstrahlung, and hydrogen and helium cooling. The particle-mesh code solves the equation of motion for the collisionless dark matter. The code uses two-dimensional Cartesian coordinates with a nonuniform grid in one direction to provide high resolution for the sheet structures. A series of one-dimensional and two-dimensional linear perturbation tests are presented which are designed to test the hydro solver and the Poisson solver with and without the expansion of the universe. We also present a radiative shock wave test which is designed to ensure the code's capability to handle radiative cooling properly. And finally a series of one-dimensional Zel'dovich pancake tests used to test the dark matter code and the hydro solver in the nonlinear regime are discussed and compared with the results of Bond et al. (1984) and Shapiro & Struck-Marcell (1985). Overall, the code is shown to produce accurate and stable results, which provide us a powerful tool to further our studies.

Anninos, Wenbo Y.↗
NASA NTRS2004

Implementation of Improved Transverse Shear Calculations and Higher Order Laminate Theory Into Strain Rate Dependent Analyses of Polymer Matrix Composites

A numerical procedure has been developed to investigate the nonlinear and strain rate dependent deformation response of polymer matrix composite laminated plates under high strain rate impact loadings. A recently developed strength of materials based micromechanics model, incorporating a set of nonlinear, strain rate dependent constitutive equations for the polymer matrix, is extended to account for the transverse shear effects during impact. Four different assumptions of transverse shear deformation are investigated in order to improve the developed strain rate dependent micromechanics model. The validities of these assumptions are investigated using numerical and theoretical approaches. A method to determine through the thickness strain and transverse Poisson's ratio of the composite is developed. The revised micromechanics model is then implemented into a higher order laminated plate theory which is modified to include the effects of inelastic strains. Parametric studies are conducted to investigate the mechanical response of composite plates under high strain rate loadings. Results show the transverse shear stresses cannot be neglected in the impact problem. A significant level of strain rate dependency and material nonlinearity is found in the deformation response of representative composite specimens.

Zhu, Lin-Fa↗
DOE OSTI2025

Algebra of invariants for the Vlasov–Maxwell system

The algebra of invariants for both the relativistic and nonrelativistic multispecies Vlasov–Maxwell system is examined, including the case with a fixed ion background. Invariants and their associated fluxes are obtained directly from the Vlasov–Maxwell system. The invariants are shown to Poisson commute with the Hamiltonian and the rest of the Poisson bracket algebra of invariants is identified. Special attention is given to the role played by the monopole condition, ∇ · B.

Fundamental invariants↗
NASA NTRS1994

Quantum statistics of Raman scattering model with Stokes mode generation

The model describing three coupled quantum oscillators with decay of Rayleigh mode into the Stokes and vibration (phonon) modes is examined. Due to the Manley-Rowe relations the problem of exact eigenvalues and eigenstates is reduced to the calculation of new orthogonal polynomials defined both by the difference and differential equations. The quantum statistical properties are examined in the case when initially: the Stokes mode is in the vacuum state; the Rayleigh mode is in the number state; and the vibration mode is in the number of or squeezed states. The collapses and revivals are obtained for different initial conditions as well as the change in time the sub-Poisson distribution by the super-Poisson distribution and vice versa.

Tanatar, Bilal↗
DOE OSTI2025

Numerical simulations of three-dimensional ion crystal dynamics in a Penning trap using the fast multipole method

We simulate the dynamics, including laser cooling, of three-dimensional (3-D) ion crystals confined in a Penning trap using a newly developed molecular dynamics-like code. The numerical integration of the ions’ equations of motion is accelerated using the fast multipole method to calculate the Coulomb interaction between ions, which allows us to efficiently study large ion crystals with thousands of ions. In particular, we show that the simulation time scales linearly with ion number, rather than with the square of the ion number. By treating the ions’ absorption of photons as a Poisson process, we simulate individual photon scattering events to study laser cooling of 3-D ellipsoidal ion crystals. Initial simulations suggest that these crystals can be efficiently cooled to ultracold temperatures, aided by the mixing of the easily cooled axial motional modes with the low frequency planar modes. In our simulations of a spherical crystal of 1000 ions, the planar kinetic energy is cooled to several millikelvin in a few milliseconds while the axial kinetic energy and total potential energy are cooled even further. This suggests that 3-D ion crystals could be well suited as platforms for future quantum science experiments.

Zaris, John (ORCID:0009000196476323)↗
NASA NTRS2001

Rheology of Foam Near the Order-Disorder Transition

The first part of our research results are summarized in the recent journal publication: J. Gregory McDaniel and R. Glynn Holt, 'Measurement of aqueous foam rheology by acoustic levitation', Phys. Rev. E 61, 2204 (2000). This aspect of the work was a combination of experiment and analysis. We built a levitation system capable of acoustically levitating small samples of aqueous foam of arbitrary gas and liquid volume fractions. We then modulated the acoustic field to induce normal mode oscillations of the foam samples. The observables from the experiment were frequency and mode number. For dry (roughly > 70% gas by volume) foams and small deformations, we developed an effective medium, normal-modes analysis which took the frequency and mode number from experiment, and gave us the shear elastic modulus of the foam as a function of Poisson's ratio. The second part of our results may be found in a soon-to-be submitted manuscript 'Dynamics of aqueous foam drops', I.Sh. Akhatov, J.G. McDaniel and R.G. Holt, describing our modeling in the wet foam limit by considering the acoustic problem. This aspect of the research is purely theoretical. Beginning from a mass-conserving mixture law, the fully nonlinear equations of motion for a wet (roughly < 10% gas by volume) foam drop of initially spherical shape were derived. The frequencies for normal mode oscillations were derived in the linear inviscid limit. The nonlinear equations were numerically solved to elicit the motion of a foam drop under acoustic excitation. The role of the time-varying void fraction in breathing-mode oscillations is of particular interest. As of the end of the current (NAG#3-2121) grant, this work was not yet concluded. We continue to work on this aspect in order to extend the analysis to cover the transition regime of gas volume fractions, as well as to compare to experiments in the wet regime.

Holt, R. Glynn↗
NASA NTRS1986

A fourth order accurate finite difference scheme for the computation of elastic waves

A finite difference for elastic waves is introduced. The model is based on the first order system of equations for the velocities and stresses. The differencing is fourth order accurate on the spatial derivatives and second order accurate in time. The model is tested on a series of examples including the Lamb problem, scattering from plane interf aces and scattering from a fluid-elastic interface. The scheme is shown to be effective for these problems. The accuracy and stability is insensitive to the Poisson ratio. For the class of problems considered here it is found that the fourth order scheme requires for two-thirds to one-half the resolution of a typical second order scheme to give comparable accuracy.

Bayliss, A.↗
NASA NTRS2009

Elongated Tetrakaidecahedron Micromechanics Model for Space Shuttle External Tank Foams

The results of microstructural characterization studies and physical and mechanical testing of BX-265 and NCFI24-124 foams are reported. A micromechanics model developed previously by the authors is reviewed, and the resulting equations for the elastic constants, the relative density, and the strength of the foam in the principal material directions are presented. The micromechanics model is also used to derive equations to predict the effect of vacuum on the tensile strength and the strains induced by exposure to vacuum. Using a combination of microstructural dimensions and physical and mechanical measurements as input, the equations for the elastic constants and the relative density are applied and the remaining microstructural dimensions are predicted. The predicted microstructural dimensions are in close agreement with the average measured values for both BX-265 and NCFI24-124. With the microstructural dimensions, the model predicts the ratio of the strengths in the principal material directions for both foams. The model is also used to predict the Poisson s ratios, the vacuum-induced strains, and the effect of vacuum on the tensile strengths. However, the comparison of these predicted values with the measured values is not as favorable.

Sullivan, Roy M.↗
NASA NTRS2020

An Improved Plastically Dilatant Unified Viscoplastic Constitutive Formulation for Multiscale Analysis of Polymer Matrix Composites Under High Strain Rate Loading

Polymer matrix composites are commonly used to fabricate energy-absorbing structures expected to experience impact loading. As such, a detailed understanding of the dynamic response of the constituent materials is necessary. Since the rate, temperature, and pressure dependence of carbon fiber reinforced polymer matrix composites are primarily manifestations of the rate, temperature, and pressure dependence of the polymer matrix, it is crucial that the constitutive behavior of the matrix be accurately characterized. In this work, an existing unified viscoplastic constitutive formulation is extended to ensure thermodynamic consistency and to more accurately account for the tension-compression asymmetry observed in the response of polymeric materials. A new plastic potential function is proposed, and elementary loading conditions are utilized to determine relations between model constants to ensure nonnegative plastic dissipation, a necessary thermodynamic requirement. Expressions for plastic Poisson’s ratios are derived and are bounded by enforcing nonnegative plastic dissipation. The model is calibrated against available experimental data from tests conducted over a range of strain rates, temperatures, and loading cases on a representative thermoset epoxy; good correlation between simulations and experimental data is obtained. Temperature rises due to the conversion of plastic work to heat are computed via the adiabatic heat energy equation. The viscoplastic polymer model is then used as a constitutive model in the generalized method of cells micromechanics theory to investigate the effects of matrix adiabatic heating on the high strain rate response of a unidirectional composite. The thermodynamic consistency of the model ensures plastic dissipation can only cause an increase in temperature. Simulation results indicate that significant thermal softening due to the conversion of plastic work to heat is observed in the composite for matrix dominated deformation modes.

Impact loading↗
NASA NTRS1982

A cylindrical shell with an arbitrarily oriented crack

The general problem of a shallow shell with constant curvatures is considered. It is assumed that the shell contains an arbitrarily oriented through crack and the material is specially orthotropic. The nonsymmetric problem is solved for arbitrary self equilibrating crack surface tractions, which, added to an appropriate solution for an uncracked shell, would give the result for a cracked shell under most general loading conditions. The problem is reduced to a system of five singular integral equations in a set of unknown functions representing relative displacements and rotations on the crack surfaces. The stress state around the crack tip is asymptotically analyzed and it is shown that the results are identical to those obtained from the two dimensional in plane and antiplane elasticity solutions. The numerical results are given for a cylindrical shell containing an arbitrarily oriented through crack. Some sample results showing the effect of the Poisson's ratio and the material orthotropy are also presented.

Yahsi, O. S.↗
NASA NTRS1983

A cylindrical shell with an arbitrarily oriented crack

The general problem of a shallow shell with constant curvatures is considered. It is assumed that the shell contains an arbitrarily oriented through crack and the material is specially orthotropic. The nonsymmetric problem is solved for arbitrary self equilibrating crack surface tractions, which, added to an appropriate solution for an uncracked shell, would give the result for a cracked shell under most general loading conditions. The problem is reduced to a system to five singular integral equations in a set of unknown functions representing relative displacements and rotations on the crack surfaces. The stress state around the crack tip is asymptotically analyzed and it is shown that the results are identical to those obtained from the two dimensional in plane and antiplane elasticity solutions. The numerical results are given for a cylindrical shell containing an arbitrarily oriented through crack. Some sample results showing the effect of the Poisson's ratio and the material orthotropy are also presented. Previously annunced in STAR as N83-16783

Yahsi, O. S.↗
DOE OSTI2025

Randomized Adiabatic Quantum Linear Solver Algorithm with Optimal Complexity Scaling and Detailed Running Costs

Solving linear systems of equations is a fundamental problem with a wide variety of applications across many fields of science, and there is increasing effort to develop quantum linear solver algorithms. Subaşı et al. [Phys. Rev. Lett. 122, 060504 (2019)] proposed a randomized algorithm inspired by adiabatic quantum computing, based on a sequence of random Hamiltonian simulation steps, with suboptimal scaling in the condition number 𝜅 of the linear system and the target error 𝜖. Here we go beyond these results in several ways. Firstly, using filtering [Lin and Tong, Quantum 4, 361 (2020)] and Poissonization techniques [Cunningham and Roland, ArXiv:2406.03972 (2024)], the algorithm complexity is improved to the optimal scaling 𝑂⁡(𝜅⁢log (1/𝜖))—an exponential improvement in 𝜖, and a shaving of a log 𝜅 scaling factor in 𝜅. Secondly, the algorithm is further modified to achieve constant factor improvements, which are vital as we progress towards hardware implementations on fault-tolerant devices. We introduce a cheaper randomized walk operator method replacing Hamiltonian simulation—which also removes the need for potentially challenging classical precomputations; randomized routines are sampled over optimized random variables; circuit constructions are improved. We obtain a closed formula rigorously upper bounding the expected number of times one needs to apply a block-encoding of the linear system matrix to output a quantum state encoding the solution to the linear system. The upper bound is 837⁢𝜅 at 𝜖 = 10 −10 for Hermitian matrices.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Implementing abstract multigrid or multilevel methods

Multigrid methods can be formulated as an algorithm for an abstract problem that is independent of the partial differential equation, domain, and discretization method. In such an abstract setting, problems not arising from partial differential equations can be treated. A general theory exists for linear problems. The general theory was motivated by a series of abstract solvers (Madpack). The latest version was motivated by the theory. Madpack now allows for a wide variety of iterative and direct solvers, preconditioners, and interpolation and projection schemes, including user callback ones. It allows for sparse, dense, and stencil matrices. Mildly nonlinear problems can be handled. Also, there is a fast, multigrid Poisson solver (two and three dimensions). The type of solvers and design decisions (including language, data structures, external library support, and callbacks) are discussed. Based on the author's experiences with two versions of Madpack, a better approach is proposed. This is based on a mixed language formulation (C and FORTRAN + preprocessor). Reasons for not using FORTRAN, C, or C++ (individually) are given. Implementing the proposed strategy is not difficult.

Douglas, Craig C.↗
NASA NTRS1994

Removing cosmic-ray hits from multiorbit HST Wide Field Camera images

We present an optimized algorithm that removes cosmic rays ('CRs') from multiorbit Hubble Space Telescope (HST) Wide Field/Planetary Camera ('WF/PC') images. It computes the image noise in every iteration from the WF/PC CCD equation. This includes all known sources of random and systematic calibration errors. We test this algorithm on WF/PC stacks of 2-12 orbits as a function of the number of available orbits and the formal Poissonian sigma-clipping level. We find that the algorithm needs greater than or equal 4 WF/PC exposures to locate the minimal sky signal (which is noticeably affected by CRs), with an optimal clipping level at 2-2.5 x sigma(sub Poisson). We analyze the CR flux detected on multiorbit 'CR stacks,' which are constructed by subtracting the best CR filtered images from the unfiltered 8-12 orbit average. We use an automated object finder to determine the surface density of CRS as a function of the apparent magnitude (or ADU flux) they would have generated in the images had they not been removed. The power law slope of the CR 'counts' (gamma approximately = 0.6 for N(m) m(exp gamma)) is steeper than that of the faint galaxy counts down to V approximately = 28 mag. The CR counts show a drop off between 28 less than or approximately V less than or approximately 30 mag (the latter is our formal 2 sigma point source sensitivity without spherical aberration). This prevents the CR sky integral from diverging, and is likely due to a real cutoff in the CR energy distribution below approximately 11 ADU per orbit. The integral CR surface density is less than or approximately 10(exp 8)/sq. deg, and their sky signal is V approximately = 25.5-27.0 mag/sq. arcsec, or 3%-13% of our NEP sky background (V = 23.3 mag/sq. arcsec), and well above the EBL integral of the deepest galaxy counts (B(sub J) approximately = 28.0 mag/sq. arcsec). We conclude that faint CRs will always contribute to the sky signal in the deepest WF/PC images. Since WFPC2 has approximately 2.7x lower read noise and a thicker CCD, this will result in more CR detections than in WF/PC, potentially affecting approximately 10%-20% of the pixels in multiorbit WFPC2 data cubes.

Windhorst, Rogier A.↗
NASA NTRS2014

Properties of Desert Sand and CMAS Glass

As-received desert sand from a Middle East country has been characterized for its phase composition and thermal stability. X-ray diffraction analysis showed the presence of quartz (SiO2), calcite (CaCO3), gypsum (CaSO4.2H2O), and NaAlSi3O8 phases in as-received desert sand and showed weight loss of approx. 35 percent due to decomposition of CaCO3 and CaSO4.2H2O when heated to 1400 C. A batch of as-received desert sand was melted into calcium magnesium aluminosilicate (CMAS) glass at approx. 1500 C. From inductively coupled plasma-atomic emission spectrometry, chemical composition of the CMAS glass was analyzed to be 27.8CaO-4MgO-5Al2O3-61.6SiO2-0.6Fe2O3-1K2O (mole percent). Various physical, thermal and mechanical properties of the glass have been evaluated. Bulk density of CMAS glass was 2.69 g/cc, Young's modulus 92 GPa, Shear modulus 36 GPa, Poisson's ratio 0.28, dilatometric glass transition temperature (T (sub g)) 706 C, softening point (T (sub d)) 764 C, Vickers microhardness 6.3 +/- 0.4 GPa, indentation fracture toughness 0.75 +/- 0.15 MPa.m (sup 1/2), and coefficient of thermal expansion (CTE) 9.8 x 10 (exp -6)/degC in the temperature range 25 to 700 C. Temperature dependence of viscosity has also been estimated from various reference points of the CMAS glass using the Vogel-Fulcher-Tamman (VFT) equation. The glass remained amorphous after heat treating at 850 C for 10 hr but crystallized into CaSiO3 and Ca-Mg-Al silicate phases at 900 C or higher temperatures. Crystallization kinetics of the CMAS glass has also been investigated by differential thermal analysis (DTA). Activation energies for the crystallization of two different phases in the glass were calculated to be 403 and 483 kJ/mol, respectively.

Desert sand↗
NASA NTRS1983

Generic buckling curves for specially orthotropic rectangular plates

Using a double affine transformation, the classical buckling equation for specially orthotropic plates and the corresponding virtual work theorem are presented in a particularly simple fashion. These dual representations are characterized by a single material constant, called the generalized rigidity ratio, whose range is predicted to be the closed interval from 0 to 1 (if this prediction is correct then the numerical results using a ratio greater than 1 in the specially orthotropic plate literature are incorrect); when natural boundary conditions are considered a generalized Poisson's ratio is introduced. Thus the buckling results are valid for any specially orthotropic material; hence the curves presented in the text are generic rather than specific. The solution trends are twofold; the buckling coefficients decrease with decreasing generalized rigidity ratio and, when applicable, they decrease with increasing generalized Poisson's ratio. Since the isotropic plate is one limiting case of the above analysis, it is also true that isotropic buckling coefficients decrease with increasing Poission's ratio.

Brunnelle, E. J.↗
NASA NTRS2000

Possible Statistics of Two Coupled Random Fields: Application to Passive Scalar

We use the relativity postulate of scale invariance to derive the similarity transformations between two coupled scale-invariant random elds at different scales. We nd the equations leading to the scaling exponents. This formulation is applied to the case of passive scalars advected i) by a random Gaussian velocity field; and ii) by a turbulent velocity field. In the Gaussian case, we show that the passive scalar increments follow a log-Levy distribution generalizing Kraichnan's solution and, in an appropriate limit, a log-normal distribution. In the turbulent case, we show that when the velocity increments follow a log-Poisson statistics, the passive scalar increments follow a statistics close to log-Poisson. This result explains the experimental observations of Ruiz et al. about the temperature increments.

Dubrulle, B.↗
NASA NTRS2001

The Constitutive Modeling of Thin Films with Randon Material Wrinkles

Material wrinkles drastically alter the structural constitutive properties of thin films. Normally linear elastic materials, when wrinkled, become highly nonlinear and initially inelastic. Stiffness' reduced by 99% and negative Poisson's ratios are typically observed. This paper presents an effective continuum constitutive model for the elastic effects of material wrinkles in thin films. The model considers general two-dimensional stress and strain states (simultaneous bi-axial and shear stress/strain) and neglects out of plane bending. The constitutive model is derived from a traditional mechanics analysis of an idealized physical model of random material wrinkles. Model parameters are the directly measurable wrinkle characteristics of amplitude and wavelength. For these reasons, the equations are mechanistic and deterministic. The model is compared with bi-axial tensile test data for wrinkled Kaptong(Registered Trademark) HN and is shown to deterministically predict strain as a function of stress with an average RMS error of 22%. On average, fitting the model to test data yields an RMS error of 1.2%

Murphey, Thomas W.↗