Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Greens function”

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 595 records · Page 33

Propagation of transients in a random medium

The propagation of transient scalar waves in a three-dimensional random medium is considered. The analysis is based on the smoothing method. An integro-differential equation for the coherent (or average) wave is derived and solved for the case of a statistically homogeneous and isotropic medium and a delta-function source. This yields the coherent Green's function of the medium. It is found that the waveform of the coherent wave depends generally on the distance from the source measured in terms of a certain dimensionless parameter. Based on the magnitude of this parameter, three propagation zones, called the near zone, the far zone, and the intermediate zone, are defined.

Wenzel, A. R.↗

Boundary-Layer Receptivity and Integrated Transition Prediction

The adjoint parabold stability equations (PSE) formulation is used to calculate the boundary layer receptivity to localized surface roughness and suction for compressible boundary layers. Receptivity efficiency functions predicted by the adjoint PSE approach agree well with results based on other nonparallel methods including linearized Navier-Stokes equations for both Tollmien-Schlichting waves and crossflow instability in swept wing boundary layers. The receptivity efficiency function can be regarded as the Green's function to the disturbance amplitude evolution in a nonparallel (growing) boundary layer. Given the Fourier transformed geometry factor distribution along the chordwise direction, the linear disturbance amplitude evolution for a finite size, distributed nonuniformity can be computed by evaluating the integral effects of both disturbance generation and linear amplification. The synergistic approach via the linear adjoint PSE for receptivity and nonlinear PSE for disturbance evolution downstream of the leading edge forms the basis for an integrated transition prediction tool. Eventually, such physics-based, high fidelity prediction methods could simulate the transition process from the disturbance generation through the nonlinear breakdown in a holistic manner.

Chang, Chau-Lyan↗

Analytical continuation of matrix-valued functions: Carathéodory formalism

Finite-temperature quantum field theories are formulated in terms of Green's functions and self-energies on the Matsubara axis. In multiorbital systems, these quantities are related to positive semidefinite matrix-valued functions of the Carathéodory and Schur class. Analysis, interpretation, and evaluation of derived quantities such as real-frequency response functions requires analytic continuation of the off-diagonal elements to the real axis. Here we derive the criteria under which such functions exist for given Matsubara data and present an interpolation algorithm that intrinsically respects their mathematical properties. For small systems with precise Matsubara data, we find that the continuation exactly recovers all off-diagonal and diagonal elements. In real-materials systems, we show that the precision of the continuation is sufficient for the analytic continuation to commute with the Dyson equation, and we show that the truncation of the off-diagonal self-energy elements leads to considerable approximation artifacts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electric response of a metal-molecule-metal junction to laser pulse by solving hierarchical equations of motion

We have combined the quantum dissipative theory and the time dependent density functional theory to perform the first principle calculation of laser induced quantum dynamical electron transport through a molecule weak bridged to two electrodes. The formalism of hierarchical equations of motion based on non-equilibrium Green’s function theory has been taken in this work. Numerical simulations of optical absorption spectra of benzene, laser induced transient current without and with bias, charge pumping effect, as well as the spectrum analysis from the current in Au-benzene-Au molecular junction are presented and discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Algebraic Structures in the Coupling of Gravity to Gauge Theories

Highlights: • Geometric foundation for Quantum General Relativity coupled to Quantum Electrodynamics. • Rigorous introduction to Hopf algebraic renormalization. • Appropriate modifications in case renormalization Hopf algebra is ill-defined. • Generalization of Furry’s Theorem to include gravitons and ghosts. • Discussion of obstructions to multiplicative renormalization. This article is an extension of the author’s second master thesis, Prinz (2017). It aims to introduce to the theory of perturbatively quantized General Relativity coupled to Spinor Electrodynamics, provide the results thereof and set the notation to serve as a starting point for further research in this direction. It includes the differential geometric and Hopf algebraic background, as well as the corresponding Lagrange density and some renormalization theory. Then, a particular problem in the renormalization of Quantum General Relativity coupled to Quantum Electrodynamics is addressed and solved by a generalization of Furry’s Theorem. Next, the restricted combinatorial Green’s functions for all two-loop propagators and all one-loop divergent subgraphs thereof are presented. Finally, relations between these one-loop restricted combinatorial Green’s functions necessary for multiplicative renormalization are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nonequilibrium electron dynamics in pump-probe spectroscopy: Role of excited phonon populations

We study the role of excited phonon populations in the relaxation rates of nonequilibrium electrons using a nonequilibrium Green's function formalism. The transient modi cations in the phononic properties are accounted for by self-consistently solving the Dyson equation for the electron and phonon Green's functions. The pump induced changes manifest in both the electronic and phononic spectral functions. We nd that the excited phonon populations suppress the decay rates of nonequilibrium electrons due to enhanced phonon absorption. The increased phonon occupation also sets the nonequilibrium decay rates and the equilibrium scattering rates apart. The decay rates are found to be time-dependent, and this is illustrated in the experimentally observed population decay of photoexcited Bi 1.5 Sb 0.5 Te 1.7 Se 1.3 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effect of Free Jet on Refraction and Noise

This article investigates the role of a free jet on the sound radiated from a jet. In particular, the role of an infinite wind tunnel, which simulates the forward flight condition, is compared to that of a finite wind tunnel. The second configuration is usually used in experiments, where the microphones are located in a static ambient medium far outside the free jet. To study the effect of the free jet on noise, both propagation and source strength need to be addressed. In this work, the exact Green's function in a locally parallel flow is derived for a simulated flight case. Numerical examples are presented that show a reduction in the magnitude of the Green's function in the aft arc and an increase in the forward arc for the simulated flight condition. The effect of finite wind tunnel on refraction is sensitive to the source location and is most pronounced in the aft arc. A Reynolds-averaged Navier-Stokes solution (RANS) yields the required mean flow and turbulence scales that are used in the jet mixing noise spectrum calculations. In addition to the sound/flow interaction, the separate effect of source strength and elongation of the noise-generating region of the jet in a forward flight is studied. Comparisons are made with experiments for the static and finite tunnel cases. Finally, the standard free-jet shear corrections that convert the finite wind tunnel measurements to an ideal wind tunnel arrangement are evaluated.

Khavaran, Abbas↗

Nonlocal, Pattern-Aware Response and Feedback Framework for Regional Climate Response

We devise a pattern-aware feedback framework for representing the forced climate response using a suite of Green’s function experiments with solar radiation perturbations. By considering the column energy balance, a comprehensive linear response function (CLRF) for important climate variables and feedback quantities such as moist static energy, sea surface temperature, albedo, cloud optical depth, and lapse rate is learned from Green’s function data. The learned CLRF delineates the effects of the energy diffusion in both the ocean and atmosphere and the pattern-aware feedbacks from the aforementioned radiatively active processes. The CLRF can then be decomposed into forcing–response mode pairs, which are in turn used to construct a reduced-order model describing the dominant dynamics of climate responses. These mode pairs capture nonlocal effects and teleconnections in the climate and thus make the reduced-order model apt for capturing regional features of climate response. A key observation is that the CLRF captures the polar-amplified response as the most excitable mode of the climate system, and this mode is explainable in the data-learned pattern-aware feedback framework. The reduced-order model can be used for predicting the response for a given forcing and for reconstructing the forcing from a given response; we demonstrate these capabilities for multiple independent forcing scenarios.

Feedback↗

Smooth indentation of orthotropic beams

An approximate Green's function for surface displacements in an orthotropic beam is derived as the superposition of the half-plane solution for displacements and the beam theory deflections. The Green's function is used to formulate the integral equation for the problem of smooth contact between a rigid-cylinder and a simply supported orthotropic beam. The integral equation is solved using a least-squares approximation procedure. The contact stress distribution is presented for three materials with different degrees of orthotropy. Numerical results are given for contact force-contact length and contact force-indentation relations. The effects of curvature of the deflected beam, shear deformation and length-to-thickness ratio of the beam on the contact behavior are discussed.

Sankar, B. V.↗

Density functional theory based simulations of silicon nanowire field effect transistors

First-principles density functional theory (DFT) based, atomistic, self-consistent device simulations are performed for realistically sized Si nanowire field effect transistors (NW FETs) having tens of thousands of atoms. Through mode space transformation, DFT Hamiltonian and overlap matrices are reduced in size from a few thousands to around one hundred. Ultra-efficient quantum-mechanical transport calculations in the non-equilibrium Green's function formalism in a non-orthogonal basis are therefore made possible. The n-type and p-type Si NW FETs are simulated and found to exhibit similar device performance in the nanoscale regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bayesian Inference for the Seismic Moment Tensor Using Regional Waveforms and Teleseismic- P Polarities with a Data-Derived Distribution of Velocity Models and Source Locations

The largest source of uncertainty in any source inversion is the velocity model used in the transfer function that relates observed ground motion to the seismic moment tensor. However, standard inverse procedure often does not quantify uncertainty in the seismic moment tensor due to error in the Green’s functions from uncertain event location and Earth structure. Here, we incorporate this uncertainty into an estimation of the seismic moment tensor using a data-derived distribution of velocity models based on complementary geophysical data sets, including thickness constraints, velocity profiles, gravity data, surface-wave group velocities, and regional body-wave travel times. The data-derived distribution of velocity models is then used as a prior distribution of Green’s functions for use in Bayesian inference of an unknown seismic moment tensor using regional and teleseismic-P waveforms. The use of multiple data sets is important for gaining resolution to different components of the moment tensor. The combined likelihood is estimated using data-specific error models and the posterior of the seismic moment tensor is estimated and interpreted in terms of the most probable source type.

58 GEOSCIENCES↗

Fully self-consistent solution of the Dyson equation using a plane-wave basis set

Self-consistent solutions of the Dyson equation are obtained using a plane-wave basis set for seven small molecules. Such self-consistent solutions can help to unify the different GW self-consistent schemes, reduce the scatter of results in current GW calculations, and shed light on the true effects of GW self-consistency. Unlike other works of self-consistent GW calculations, in the present work the Green's function is expressed as a matrix under the plane-wave basis set. The algorithmic details which enable such calculations are presented. Finally, the ability to solve the full Green's function using a plane-wave basis set may open the door for future beyond- GW many-body perturbation theory calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A simulation technique for 3D MR-guided acoustic radiation force imaging

Purpose: In magnetic resonance-guided focused ultrasound (MRgFUS) therapies, the in situ characterization of the focal spot location and quality is critical. MR acoustic radiation force imaging (MR-ARFI) is a technique that measures the tissue displacement caused by the radiation force exerted by the ultrasound beam. This work presents a new technique to model the displacements caused by the radiation force of an ultrasound beam in a homogeneous tissue model. Methods: When a steady-state point-source force acts internally in an infinite homogeneous medium, the displacement of the material in all directions is given by the Somigliana elastostatic tensor. The radiation force field, which is caused by absorption and reflection of the incident ultrasound intensity pattern, will be spatially distributed, and the tensor formulation takes the form of a convolution of a 3D Green’s function with the force field. The dynamic accumulation of MR phase during the ultrasound pulse can be theoretically accounted for through a time-of-arrival weighting of the Green’s function. This theoretical model was evaluated experimentally in gelatin phantoms of varied stiffness (125-, 175-, and 250-bloom). The acoustic and mechanical properties of the phantoms used as parameters of the model were measured using independent techniques. Displacements at focal depths of 30- and 45-mm in the phantoms were measured by a 3D spin echo MR-ARFI segmented-EPI sequence. Results: The simulated displacements agreed with the MR-ARFI measured displacements for all bloom values and focal depths with a normalized RMS difference of 0.055 (range 0.028–0.12). The displacement magnitude decreased and the displacement pattern broadened with increased bloom value for both focal depths, as predicted by the theory. Conclusions: A new technique that models the displacements caused by the radiation force of an ultrasound beam in a homogeneous tissue model theory has been rigorously validated through comparison with experimentally obtained 3D displacement data in homogeneous gelatin phantoms using a 3D MR-ARFI sequence. The agreement of the experimentally measured and simulated results demonstrates the potential to use MR-ARFI displacement data in MRgFUS therapies.

60 APPLIED LIFE SCIENCES↗

Beyond the Rigid Lid: Baroclinic Modes in a Structured Atmosphere

Abstract The baroclinic-mode decomposition is a fixture of the tropical-dynamics literature because of its simplicity and apparent usefulness in understanding a wide range of atmospheric phenomena. However, its derivation relies on the assumption that the tropopause is a rigid lid that artificially restricts the vertical propagation of wave energy. This causes tropospheric buoyancy anomalies of a single vertical mode to remain coherent for all time in the absence of dissipation. Here, the authors derive the Green’s functions for these baroclinic modes in a two-dimensional troposphere (or, equivalently, a three-dimensional troposphere with one translational symmetry) that is overlain by a stratosphere. These Green’s functions quantify the propagation and spreading of gravity waves generated by a horizontally localized heating, and they can be used to reconstruct the evolution of any tropospheric heating. For a first-baroclinic two-dimensional right-moving or left-moving gravity wave with a characteristic width of 100 km, its initial horizontal shape becomes unrecognizable after 4 h, at which point its initial amplitude has also been reduced by a factor of 1/π. After this time, the gravity wave assumes a universal shape that widens linearly in time. For gravity waves on a periodic domain the length of Earth’s circumference, it takes only 10 days for the gravity waves to spread their buoyancy throughout the entire domain.

54 ENVIRONMENTAL SCIENCES↗

Shielding from space radiations

This Progress Report covering the period of 1 June 1993 to 1 Dec. 1993 presents the development of an analytical solution to the heavy ion transport equation in terms of a one-layer Green's function formalism. The mathematical developments are recasted into an efficient computer code for space applications. The efficiency of this algorithm is accomplished by a nonperturbative technique of extending the Green's function over the solution domain. The code may also be applied to accelerator boundary conditions to allow code validation in laboratory experiments. Results from the isotopic version of the code with 80 isotopes present for a single layer target material, for the case of an Iron beam projectile at 600 MeV/nucleon in water is presented.

Chang, C. Ken↗

Shielding from space radiations

This Progress Report covering the period of 1 June 1993 to 1 Dec. 1993 presents the development of an analytical solution to the heavy ion transport equation in terms of a one-layer Green's function formalism. The mathematical developments are recasted into an efficient computer code for space applications. The efficiency of this algorithm is accomplished by a nonperturbative technique of extending the Green's function over the solution domain. The code may also be applied to accelerator boundary conditions to allow code validation in laboratory experiments. Results from the isotopic version of the code with 80 isotopes present for a single layer target material, for the case of an iron beam projectile at 600 MeV/nucleon in water is presented.

Chang, C. Ken↗

The Use of Vegetation For Estimating Broken-Cloud Optical Properties From Surface Measurements

It is known that because of complex three-dimensional (3D) radiative effects of broken clouds, the retrieval of cloud optical properties from upward measurements based on a one-dimensional (1D) inversion technique almost surely fails. To remove radiative effects of 3D cloud structure, we have developed a new technique that retrieves cloud optical thickness for broken clouds above green vegetation from simultaneous surface measurements in the VIS and Near Infrared (NIR) spectral regions. The theoretical basis of the method is the very different spectral behavior of cloud liquid water drops and green vegetation. For example, cloud optical properties, and hence cloud reflectivities, change little between 650 and 860 nm, while the vegetated surface albedo changes from 0.05 to 0.5 between the same two wavelengths. This spectral contrast in surface albedo suggests using ground measurements at both wavelengths not independently, but as an algebraic combination (a spectral index). For a spectral band in the NIR region, the green vegetation acts as a powerful reflector that "illuminates" horizontally inhomogeneous clouds from below. This provides the extra information needed to largely remove the 3D radiative effects, especially in the case of broken clouds; this in turn allows the retrieval of cloud optical depth using traditional 1D radiative transfer theory. This approach is similar to the so-called Green's function problem for radiative transfer where a laser beam illuminates clouds and the resulting "spot-size" of the reflected light around the beam characterizes cloud properties. We generalize Green's function theory to surf ace-cloud interaction and develop new spectral indices from which broken-cloud optical depth can be retrieved.

Marshak, A.↗

Bayesian inference for the seismic moment tensor using regional waveforms and a data-derived distribution of velocity models

The largest source of uncertainty in any source inversion is the velocity model used to construct the transfer function employed in the forward model that relates observed ground motion to the seismic moment tensor. However, standard inverse procedures often does not quantify uncertainty in the seismic moment tensor due to error in the Green’s functions from uncertain event location and Earth structure. We attempt to incorporate this uncertainty into an estimation of the seismic moment tensor using a distribution of velocity models calculated in a prior effort based on different and complementary data sets. The posterior distribution of velocity models is then used to construct Green’s functions for use in Bayesian inference of an unknown seismic moment tensor using regional waveform data. The combined likelihood is estimated using data-specific error models and the posterior of the seismic moment tensor is estimated and can be interpreted in terms of most-probable source-type.

58 GEOSCIENCES↗