Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Green's function methods”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Unified wavefront singularity characterization of three-dimensional elastodynamic time-domain half-space Green's function under impulsive boundary and internal loads

Founded on a novel analytical formulation that led to a rigorous yet compact path-integral representation of the time-domain elastodynamic half-space Green's function, a unified analysis of the possible occurrence of different singular wavefront behaviour in the response under arbitrary impulsive internal or surface point loads at arbitrary source-receiver locations is presented. With the decomposition of the general solution into distinct initiating and reflected wave group integrals that share a common factored format and simple contour definitions, the mathematical framework is shown to allow a straightforward identification of the specific conditions and the particular wave groups that are responsible for the singular wavefront phenomena without resorting to advanced analytic function theories or asymptotic methods. Analytic characterizations of the nature, strength and direction of all intrinsic singular wavefront behaviours of the three-dimensional Green's function in three canonical cases of source-receiver configurations are given in a dual integral-closed form format to facilitate their theoretical understanding as well as computational applications. Graphical illustrations of their variation with the source-receiver configuration and the medium's Poisson's ratio together with relevant comparison and clarifications of some classical treatments are included.

Green’s functions↗
DOE OSTI2017

Quantum Monte Carlo Calculations in Nuclear Theory

To a remarkable extent, atomic nuclei can be described as collections of point-like particles whose dynamics are dictated by a nonrelativistic Hamiltonian involving two- and three-nucleon potentials.Ab initio approaches are aimed at solving the many-body Schr$ \overline{o} $dinger equation associated with the nuclear Hamiltonian, which is a highly non-trivial problem because of the nonperturbative nature and the strong spin-isospin dependence of the nuclear forces. Within this framework it is possible to disentangle the theoretical uncertainty coming from modeling the nuclear potential and currents from that due to the approximations that are usually inherent in many-body techniques.Our many-body method of choice is quantum Monte Carlo (QMC), in particular Green’s function Monte Carlo (GFMC), which allows solving the nuclear Schrodinger equation with the required 1%accuracy level for both the ground- and the low-lying excited states of A ≤ 12 (A is the number of nucleons) nuclei. Since 2000 we have been achieving excellent agreement with experiment for nuclei with increasing number of nucleons by using the Argonne V18 two-nucleon interaction and the Illinois three-nucleon interactions (AV18+IL7). During this period many other groups have been developing chiral effective field theory (χEFT) potentials; these are generally nonlocal and notsuitable for QMC methods. We and collaborators have recently developed a local χEFT potential including the Δ-isobar degrees of freedom which can be used with QMC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

Non-equilibrium magnetic interactions in strongly correlated systems

We formulate a low-energy theory for the magnetic interactions between electrons in the multi-band Hubbard model under non-equilibrium conditions determined by an external time-dependent electric field which simulates laser-induced spin dynamics. We derive expressions for dynamical exchange parameters in terms of non-equilibrium electronic Green functions and self-energies, which can be computed, e.g., with the methods of time-dependent dynamical mean-field theory. Moreover, we find that a correct description of the system requires, in addition to exchange, a new kind of magnetic interaction, that we name twist exchange, which formally resembles Dzyaloshinskii–Moriya coupling, but is not due to spin–orbit, and is actually due to an effective three-spin interaction. Our theory allows the evaluation of the related time-dependent parameters as well. -- Highlights: •We develop a theory for magnetism of strongly correlated systems out of equilibrium. •Our theory is suitable for laser-induced ultrafast magnetization dynamics. •We write time-dependent exchange parameters in terms of electronic Green functions. •We find a new magnetic interaction, a “twist exchange”. •We give general expressions for magnetic noise in itinerant-electron systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Quantum conductance of zigzag graphene oxide nanoribbons

The electronic properties of zigzag graphene oxide nanoribbons (ZGOR) are presented. The results show interesting behaviors which are considerably different from the properties of the perfect graphene nanoribbons (GNRs). The theoretical methods include a Huckel-tight binding approach, a Green's function methodology, and the Landauer formalism. The presence of oxygen on the edge results in band bending, a noticeable change in density of states and thus the conductance. Consequently, the occupation in the valence bands increase for the next neighboring carbon atom in the unit cell. Conductance drops in both the conduction and valence band regions are due to the reduction of allowed k modes resulting from band bending. The asymmetry of the energy band structure of the ZGOR is due to the energy differences of the atoms. The inclusion of a foreign atom's orbital energies changes the dispersion relation of the eigenvalues in energy space. These novel characteristics are important and valuable in the study of quantum transport of GNRs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

The Hubbard dimer: a density functional case study of a many-body problem

This review explains the relationship between density functional theory and strongly correlated models using the simplest possible example, the two-site Hubbard model. The relationship to traditional quantum chemistry is included. Even in this elementary example, where the exact ground-state energy and site occupations can be found analytically, there is much to be explained in terms of the underlying logic and aims of density functional theory. Although the usual solution is analytic, the density functional is given only implicitly. We overcome this difficulty using the Levy–Lieb construction to create a parametrization of the exact function with negligible errors. The symmetric case is most commonly studied, but we find a rich variation in behavior by including asymmetry, as strong correlation physics vies with charge-transfer effects. We explore the behavior of the gap and the many-body Green's function, demonstrating the 'failure' of the Kohn–Sham (KS) method to reproduce the fundamental gap. We perform benchmark calculations of the occupation and components of the KS potentials, the correlation kinetic energies, and the adiabatic connection. We test several approximate functionals (restricted and unrestricted Hartree–Fock and Bethe ansatz local density approximation) to show their successes and limitations. Here, we also discuss and illustrate the concept of the derivative discontinuity. Useful appendices include analytic expressions for density functional energy components, several limits of the exact functional (weak- and strong-coupling, symmetric and asymmetric), various adiabatic connection results, proofs of exact conditions for this model, and the origin of the Hubbard model from a minimal basis model for stretched H 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

The impact of including spatially longitudinal heterogeneities of vessel oxygen content and vascular fraction in 3D tumor oxygenation models on predicted radiation sensitivity

Purpose: Oxygen distribution models have been used to analyze the influences of oxygen tensions on tissue response after radiotherapy. These distributions are often generated assuming constant oxygen tension in the blood vessels. However, as red blood cells progress through the vessels, oxygen is continuously released into the plasma and the surrounding tissue, resulting in longitudinally varying oxygen levels in the blood vessels. In the present study, the authors investigated whether a tumor oxygenation model that incorporated longitudinally varying oxygen levels would provide different predictions of necrotic fractions and radiosensitivity compared to commonly used models with a constant oxygen pressure. Methods: Our models simulated oxygen diffusion based on a Green's function approach and oxygen consumption according to the Michaelis-Menten equation. The authors constructed tumor models with different vascular fractions (VFs), from which they generated depth oxygenation curves and a look-up table of oxygen pressure gradients. The authors evaluated models of spherical tumors of various sizes, from 1 to 10{sup 4} mg. The authors compared the results from a model with constant vessel oxygen (CVO) pressure to those from models with longitudinal variations in oxygen saturation and either a constant VF (CVF) or variable VF (VVF) within the tumor tissue. The authors monitored the necrotic fractions, defined as tumor regions with an oxygen pressure below 1 mmHg. Tumor radiation sensitivity was expressed as D{sub 99,} the homogeneous radiation dose required for a tumor control probability of 0.99. Results: In the CVO saturation model, no necrosis was observed, and decreasing the VF could only decrease the D{sub 99} by up to 10%. Furthermore, the D{sub 99} vs VF dependence was similar for different tumor masses. Compared to the CVO model, the extended CVF and VVF models provided clearly different results, including pronounced effects of VF and tumor size on the necrotic fraction and D{sub 99}, necrotic fractions ranging from 0% to 97%, and a maximal D{sub 99} increment of 57%. Only minor differences were observed between different vessel architectures, i.e., CVF vs VVF. In the smallest tumor with a low necrotic fraction, the D{sub 99} strictly decreased with increasing blood velocity. Increasing blood velocity also decreased the necrotic fraction in all tumor sizes. VF had the most profound influence on both the necrotic fraction and on D{sub 99}. Conclusions: Our present analysis of necrotic formation and the impact of tumor oxygenation on D{sub 99} demonstrated the importance of including longitudinal variations in vessel oxygen content in tumor models. For small tumors, radiosensitivity was particularly dependent on VF and slightly dependent on the blood velocity and vessel arrangement. These dependences decreased with increasing tumor size, because the necrotic fraction also increased, thereby decreasing the number of viable tumor cells that required sterilization. The authors anticipate that the present model will be useful for estimating tumor oxygenation and radiation response in future detailed studies.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS1975

Turbulent Chemically Reacting Flows According to a Kinetic Theory

A review of various methods of calculating turbulent chemically reacting flow such as the Green Function, Navier-Stokes equation, and others is presented. Nonequilibrium degrees of freedom were employed to study the mixing behavior of a multiscale turbulence field. Classical and modern theories are discussed.

Hong, Z. C.↗
NASA NTRS1995

Potential theoretic methods for far field sound radiation calculations

In the area of computational acoustics, procedures which accurately predict the far-field sound radiation are much sought after. A systematic development of such procedures are found in a sequence of papers by Atassi. The method presented here is an alternate approach to predicting far field sound based on simple layer potential theoretic methods. The main advantages of this method are: it requires only a simple free space Green's function, it can accommodate arbitrary shapes of Kirchoff surfaces, and is readily extendable to three-dimensional problems. Moreover, the procedure presented here, though tested for unsteady lifting airfoil problems, can easily be adapted to other areas of interest, such as jet noise radiation problems. Results are presented for lifting airfoil problems and comparisons are made with the results reported by Atassi. Direct comparisons are also made for the flat plate case.

Hariharan, S. I.↗
DOE OSTI2019

Comparing EGF Methods for Estimating Corner Frequency and Stress Drop from P-wave Spectra

Empirical Green's functions (EGF) are widely applied to correct earthquake spectra for attenuation and other path effects in order to estimate corner frequencies and stress drops, but these source parameter estimates often exhibit poor agreement between different studies. We investigate this issue by analyzing a compact cluster of over 3000 aftershocks of the 1992 Landers earthquake. We apply and compare two different analysis and modeling methods: (1) the spectral decomposition and global EGF fitting approach, and (2) a more traditional EGF method of modeling spectral ratios. We find that spectral decomposition yields event terms that are consistent with stacks of spectral ratios for individual events, but source-parameter estimates nonetheless vary between the methods. The main source of differences comes from the modeling approach used to estimate the EGF. The global EGF-fitting approach suffers from parameter tradeoffs among the absolute stress drop, the stress-drop scaling with moment, and the high-frequency falloff rate, but has the advantage that the relative spectral shapes and stress drops among the different events in the cluster are well-resolved even if their absolute levels are not. The spectral-ratio approach solves for a different EGF for each target event without imposing any constraint on the corner frequency, f c , of the smaller events, and so can produce biased results for target-event f c . Placing constraints on the small-event f c improves the performance of the spectral-ratio method and allows the two methods to yield very similar results.

58 GEOSCIENCES↗
NASA NTRS1973

Geodesic synchrotron radiation

This paper presents the results and methods of computing the high-frequency radiation emitted by freely falling particles moving in circular geodesic orbits in a spherically symmetric gravitational field. The high-frequency radiation, to which the methods of this paper apply, is the principal part of radiated energy only in the case of a particle moving in a highly relativistic, and therefore unstable, circular geodesic. The geodesic synchrotron radiation emitted in this case shows excitation of high-frequency harmonics and a narrow angular distribution. A Green's-function solution of the scalar wave equation is obtained using WKB methods. For application to relativistic circular orbits, a parabolic WKB approximation is required and yields solutions in terms of parabolic cylinder functions.

Breuer, R. A.↗
DOE OSTI2019

Acoustic Source Parameters and Yield Estimation for SPE Phase I and Forensics Surface Events

Acoustic signals recorded for explosive events are sensitive to pressure changes in the source region and can provide valuable information about source characteristics. However, acoustic wave propagation in the atmosphere is significantly affected by weather conditions (e.g., temperature, pressure, and wind), and the atmospheric propagation effects should be considered properly to obtain accurate source information. Acoustic waveform inversion method has been increasing useful to determine quantitative source parameters by incorporating accurate Green’s function into the inversion and by effectively removing propagation path effects. This paper reviews recent progress in the acoustic waveform inversion and its application to explosion yield/emplacement (e.g., height-of-burst or depth-of-burial) estimation. The inversion method is applied to acoustic signals recorded during Source Physics Experiment Phase I and Forensics Surface Event conducted at the Nevada National Security Site (NNSS). Inverted acoustic source strengths are related to the explosion yield and emplacement, and their estimation uncertainties are explored by the comparison with the ground-truth data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS2017

Imaging a Free-Geometry Magma Reservoir Using FEM-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗
NASA NTRS2017

Imaging the Complex Geometry of a Magma Reservoir Using Fem-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluid-filled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗
DOE OSTI2023

Phase-field modeling of dislocation–interstitial interactions

The mechanical behavior of body-centered cubic (BCC) materials can be dramatically affected by the presence of interstitial solute atoms. Here we present a new phase-field dislocation dynamics formulation to include the diffusion of interstitials. Short-range interactions are accounted for by a concentration-dependent lattice energy, and long-range interactions are accounted for by modifications to the elastic energy. The interstitial diffusion law introduces gradients that require methods for minimizing Gibbs oscillations, which is done via a modified Green’s function. The formulation is general to any solute-solvent system and is applied here to Nb-O as a model system, whose interstitial parameters are obtained from ab initio calculations. The effect of O on the core structures of Nb edge and screw dislocations is calculated. The diffusion of O to form interstitial atmospheres around dislocation cores is simulated, as well as the critical stresses required for dislocations to break away or cross slip from these atmospheres. Future applications of the method to simulate complex interstitial embrittlement mechanisms are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Improved quasiparticle self-consistent electronic band structure and excitons in β – LiGaO 2

The band structure of β–LiGaO 2 is calculated using the quasiparticle self-consistent QSGWˆ method where the screened Coulomb interaction Wˆ is evaluated including electron-hole interaction ladder diagrams and G is the one-electron Green's function. Improved convergence compared to previous calculations leads to a significantly larger band gap of about 7.0 eV. However, exciton binding energies are found to be large and lead to an exciton gap of about 6.0 eV if also a zero-point-motion correction of about –0.4 eV is included. Furthermore, these results are in excellent agreement with recent experimental results on the onset of absorption. Besides the excitons observed thus far, the calculations indicate the existence of a Rydberg-like series of exciton excited states, which is however modified from the classical Wannier exciton model by the anisotropies of the material and the more complex mixing of Bloch states in the excitons resulting from the Bethe-Salpeter equation. The exciton fine structure and the exciton wave functions are visualized and analyzed in various ways.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Time-dependent integral transport in one-dimensional infinite media using dimensionless variables and the reduced collision formulation

A new method for calculating time-dependent integral neutron transport single collision fluxes has been developed for one-dimensional infinite media problems using dimensionless variables and the reduced collision formulation. By transforming the integral scalar flux equation into dimensionless variables, one may immediately integrate over the time domain and reduce the problem complexity to integrations strictly over the spatial domain, significantly reducing the computational work required to obtain a solution. The results have been compared to published benchmark solutions for validation. Further, reformulating the scalar flux as a Green’s function has allowed for the determination of an analytical expression for the relative error of any collision level scalar flux to the summed scalar flux. This provides the absolute number of collisions required for convergence at any given mean free time, using previous work in one-dimensional Cartesian transport and derived work herein for the one-dimensional spherical transport.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Optimization performance, fidelity, and cost: SIAM VQE

This dataset contains files storing results from classically-simulated quantum subroutines within a dynamical mean-field theory workflow, and jupyter notebooks processing the data in these files to generate plots. The files store: (1) Results from variational quantum eigensolver (VQE) simulations searching for optimal parameters allowing parametrized quantum circuits to prepare approximations to ground states of different Anderson impurity models (AIMs) (2) Results from simulations of a quantum Lanczos algorithm (QLA) estimating the Lanczos coefficients defining the continued-fraction representation of an (AIM) Green’s function Description: Any file named vqe_gs_results* stores approximations to the ground state and energy of a given AIM estimated using three different methods: (1) Numerical diagonalization (2) Ideal VQE simulation (3) VQE simulation with sampling noise For each VQE simulations metadata about the optimization (optimization results plus number of quantum circuits that would have been executed on real hardware) is also stored. Any file named qla_dos_results* estimations for the Lanczos coefficients defining the Green’s function of an AIM. The stored estimations are achieved using different methods: (1) Numerical Lanczos algorithm from initial states obtained from numerical diagonalization (2) Simulated quantum Lanczos algorithm from initial states prepared from parametrized quantum circuits yielded by corresponding ideal and noisy VQE subroutines. The dataset is used and described in M. Karabin et al., "Quantum solver for single-impurity Anderson models with particle-hole symmetry", Phys. Rev. Research 8, 033066 (2026). DOI: https://doi.org/10.1103/7ys3-tl4l

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Molecular Junctions: Control of the Energy Gap Achieved by a Pinning Effect

Single-molecule junctions are the constitutive components of molecular electronics circuits. For any potential application, the energy gap in the junction, i.e., the accumulated energy difference between the electrode Fermi level and the two frontier energy levels of the molecule, is a key property. In this study, using the nonequilibrium Green’s function coupled to the density functional theory framework (NEGF-DFT) method, we show that the gap of the molecule inserted between electrodes can differ largely from the gap of the same molecule, at the isolated level. It can be widely compressed by tuning the alignment mechanism at each metal/molecule interface. In the context of molecular rectification, we show that this mechanism relates to the pinning effect. We discuss the different parameters affecting the compression of the gap and its efficiency. Interestingly, we find that the structure both of the molecule and of the anchoring group plays an important role. Finally, we investigate the evolution of these features out-of-equilibrium.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗