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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2016

NEUTRON STAR MASS–RADIUS CONSTRAINTS USING EVOLUTIONARY OPTIMIZATION

The equation of state of cold supra-nuclear-density matter, such as in neutron stars, is an open question in astrophysics. A promising method for constraining the neutron star equation of state is modeling pulse profiles of thermonuclear X-ray burst oscillations from hot spots on accreting neutron stars. The pulse profiles, constructed using spherical and oblate neutron star models, are comparable to what would be observed by a next-generation X-ray timing instrument like ASTROSAT , NICER , or a mission similar to LOFT . In this paper, we showcase the use of an evolutionary optimization algorithm to fit pulse profiles to determine the best-fit masses and radii. By fitting synthetic data, we assess how well the optimization algorithm can recover the input parameters. Multiple Poisson realizations of the synthetic pulse profiles, constructed with 1.6 million counts and no background, were fitted with the Ferret algorithm to analyze both statistical and degeneracy-related uncertainty and to explore how the goodness of fit depends on the input parameters. For the regions of parameter space sampled by our tests, the best-determined parameter is the projected velocity of the spot along the observer’s line of sight, with an accuracy of ≤3% compared to the true value and with ≤5% statistical uncertainty. The next best determined are the mass and radius; for a neutron star with a spin frequency of 600 Hz, the best-fit mass and radius are accurate to ≤5%, with respective uncertainties of ≤7% and ≤10%. The accuracy and precision depend on the observer inclination and spot colatitude, with values of ∼1% achievable in mass and radius if both the inclination and colatitude are ≳60°.

79 ASTRONOMY AND ASTROPHYSICS↗
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.↗
DOE OSTI2019

The Shifted Interface Method: A flexible approach to embedded interface computations

We propose a new embedded finite element method to simulate partial differential equations over domains with internal interfaces. Our approach belongs to the family of surrogate/approximate interface methods and relies on the idea of shifting the location and value of jump interface conditions. This choice has the goal of preserving optimal convergence rates while avoiding small cut cells and related problematic issues, typical of traditional embedded methods. The proposed approach is accurate, robust, efficient, and simple to implement. We apply this concept to internal interface computations in the context of the mixed Poisson problem, also known as the Darcy flow problem. An extensive set of numerical tests is provided to demonstrate the performance of the proposed approach.

42 ENGINEERING↗
DOE OSTI2014

High pressure elastic properties of minerals from ab initio simulations: The case of pyrope, grossular and andradite silicate garnets

A computational strategy is devised for the accurate ab initio simulation of elastic properties of crystalline materials under pressure. The proposed scheme, based on the evaluation of the analytical stress tensor and on the automated computation of pressure-dependent elastic stiffness constants, is implemented in the CRYSTAL solid state quantum-chemical program. Elastic constants and related properties (bulk, shear and Young moduli, directional seismic wave velocities, elastic anisotropy index, Poisson's ratio, etc.) can be computed for crystals of any space group of symmetry. We apply such a technique to the study of high-pressure elastic properties of three silicate garnet end-members (namely, pyrope, grossular, and andradite) which are of great geophysical interest, being among the most important rock-forming minerals. The reliability of this theoretical approach is proved by comparing with available experimental measurements. The description of high-pressure properties provided by several equations of state is also critically discussed.

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

A high-order integral solver for scalar problems of diffraction by screens and apertures in three-dimensional space

We present a novel methodology for the numerical solution of problems of diffraction by infinitely thin screens in three-dimensional space. Our approach relies on new integral formulations as well as associated high-order quadrature rules. The new integral formulations involve weighted versions of the classical integral operators related to the thin-screen Dirichlet and Neumann problems as well as a generalization to the open-surface problem of the classical Calderón formulae. The high-order quadrature rules we introduce for these operators, in turn, resolve the multiple Green function and edge singularities (which occur at arbitrarily close distances from each other, and which include weakly singular as well as hypersingular kernels) and thus give rise to super-algebraically fast convergence as the discretization sizes are increased. When used in conjunction with Krylov-subspace linear algebra solvers such as GMRES, the resulting solvers produce results of high accuracy in small numbers of iterations for low and high frequencies alike. We demonstrate our methodology with a variety of numerical results for screen and aperture problems at high frequencies—including simulation of classical experiments such as the diffraction by a circular disc (featuring in particular the famous Poisson spot), evaluation of interference fringes resulting from diffraction across two nearby circular apertures, as well as solution of problems of scattering by more complex geometries consisting of multiple scatterers and cavities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Structure and osmotic pressure of ionic microgel dispersions

We investigate structural and thermodynamic properties of aqueous dispersions of ionic microgels—soft colloidal gel particles that exhibit unusual phase behavior. Starting from a coarse-grained model of microgel macroions as charged spheres that are permeable to microions, we perform simulations and theoretical calculations using two complementary implementations of Poisson-Boltzmann (PB) theory. Within a one-component model, based on a linear-screening approximation for effective electrostatic pair interactions, we perform molecular dynamics simulations to compute macroion-macroion radial distribution functions, static structure factors, and macroion contributions to the osmotic pressure. For the same model, using a variational approximation for the free energy, we compute both macroion and microion contributions to the osmotic pressure. Within a spherical cell model, which neglects macroion correlations, we solve the nonlinear PB equation to compute microion distributions and osmotic pressures. By comparing the one-component and cell model implementations of PB theory, we demonstrate that the linear-screening approximation is valid for moderately charged microgels. By further comparing cell model predictions with simulation data for osmotic pressure, we chart the cell model’s limits in predicting osmotic pressures of salty dispersions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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.↗
DOE OSTI2015

Boundary scattering of phonons: Specularity of a randomly rough surface in the small-perturbation limit

Scattering of normally incident longitudinal and transverse acoustic waves by a randomly rough surface of an elastically isotropic solid is analyzed within the small-perturbation approach. In the limiting case of a large correlation length L compared with the acoustic wavelength, the specularity reduction is given by 4η 2 k 2 ,where η is the rms roughness and k is the acoustic wave vector, which is in agreement with the well-known Kirchhoff approximation result often referred to as Ziman's equation [J. M. Ziman, Electrons and Phonons (Clarendon Press, Oxford, 1960)]. In the opposite limiting case of a small correlation length, the specularity reduction is discovered to be proportional to η 2 k 4 L 2 , with the fourth power dependence on frequency as in Rayleigh scattering. Numerical calculations for a Gaussian autocorrelation function of surface roughness connect these limiting cases and reveal a maximum of diffuse scattering at an intermediate value of L . This maximum becomes increasingly pronounced for the incident longitudinal wave as the Poisson's ratio of the medium approaches 1/2 as a result of increased scattering into transverse and Rayleigh surface waves. The results imply that thermal transport models using Ziman's formula are likely to overestimate the heat flux dissipation due to boundary scattering, whereas modeling interface roughness as atomic disorder is likely to underestimate scattering.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Thermodynamically consistent Cahn–Hilliard–Navier–Stokes equations using the metriplectic dynamics formalism

Cahn–Hilliard–Navier–Stokes (CHNS) systems describe flows with two-phases, e.g., a liquid with bubbles. Obtaining constitutive relations for general dissipative processes for such systems, which are thermodynamically consistent, can be a challenge. We show how the metriplectic 4-bracket formalism (Morrison and Updike, 2024) achieves this in a straightforward, in fact algorithmic, manner. First, from the noncanonical Hamiltonian formulation for the ideal part of a CHNS system we obtain an appropriate Casimir to serve as the entropy in the metriplectic formalism that describes the dissipation (e.g. viscosity, heat conductivity and diffusion effects). General thermodynamics with the concentration variable and its thermodynamics conjugate, the chemical potential, are included. Having expressions for the Hamiltonian (energy), entropy, and Poisson bracket, we describe a procedure for obtaining a metriplectic 4-bracket that describes thermodynamically consistent dissipative effects. The 4-bracket formalism leads naturally to a general CHNS system that allows for anisotropic surface energy effects. Furthermore, this general CHNS system reduces to cases in the literature, to which we can compare.

Cahn–Hilliard↗
DOE OSTI2024

Inferring Stochastic Rates from Heterogeneous Snapshots of Particle Positions

Many imaging techniques for biological systems—like fixation of cells coupled with fluorescence microscopy—provide sharp spatial resolution in reporting locations of individuals at a single moment in time but also destroy the dynamics they intend to capture. In this study, these snapshot observations contain no information about individual trajectories, but still encode information about movement and demographic dynamics, especially when combined with a well-motivated biophysical model. The relationship between spatially evolving populations and single-moment representations of their collective locations is well-established with partial differential equations (PDEs) and their inverse problems. However, experimental data is commonly a set of locations whose number is insufficient to approximate a continuous-in-space PDE solution. Here, motivated by popular subcellular imaging data of gene expression, we embrace the stochastic nature of the data and investigate the mathematical foundations of parametrically inferring demographic rates from snapshots of particles undergoing birth, diffusion, and death in a nuclear or cellular domain. Toward inference, we rigorously derive a connection between individual particle paths and their presentation as a Poisson spatial process. Using this framework, we investigate the properties of the resulting inverse problem and study factors that affect quality of inference. One pervasive feature of this experimental regime is the presence of cell-to-cell heterogeneity. Rather than being a hindrance, we show that cell-to-cell geometric heterogeneity can increase the quality of inference on dynamics for certain parameter regimes. Altogether, the results serve as a basis for more detailed investigations of subcellular spatial patterns of RNA molecules and other stochastically evolving populations that can only be observed for single instants in their time evolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Elastic and magnetic properties of Fe 3 P up to core pressures: Phosphorus in the Earth's core

Phosphorus (P) is considered a possible light element alloying with iron (Fe) in the Earth's core due to its siderophile nature and the ubiquity of P-bearing iron alloys in iron meteorites. The sequestration of P by liquid metals during the core formation possibly results in the relatively low concentration of P in the bulk silicate Earth. Here, we performed single-crystal and powder X-ray diffraction, synchrotron Mossbauer spectroscopy and nuclear resonant inelastic X-ray scattering measurements in diamond anvil cells to investigate the elastic and magnetic properties of Fe 3 P under high pressures. Our X-ray diffraction results suggest that there is no structural phase transition up to 111 GPa. However, a volume collapse was observed at 21.5 GPa in Fe 3 P, ascribed to a magnetic transition as evidenced by synchrotron Mhssbauer spectroscopy results. Fitting the volume-pressure data by the Birch-Murnaghan equation of state gives bulk modulus $K_{T0}$ = 162.4(7) GPa, its first pressure derivative $K'_{T0}$ = 4.0 (fixed) and zeropressure volume $V_0$ = 370.38(6) Å 3 for the magnetic phase and $K_{T0}$ = 220(7) GPa, $K'_{T0}$ = 4.0 (fixed) and $V_0$ = 357(1) Å 3 for the non-magnetic phase. Sound velocities of Fe 3 P were determined up to 152 GPa by nuclear resonant inelastic X-ray scattering, demonstrating that Fe 3 P bears a low shear velocity and high Poisson's ratio at core pressures compared to Fe and Fe 3 S. When forming a solid solution Fe 3 (S,P) with Fe 3 S at core pressures, Fe 3 P may favorably influence the elastic properties of Fe 3 (S,P) to match the seismic observations of the inner core.

58 GEOSCIENCES↗