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 469 records · Page 26

Real-time equation-of-motion CC cumulant and CC Green’s function simulations of photoemission spectra of water and water dimer

Newly developed coupled-cluster (CC) methods enable simulations of ionization potentials and spectral functions of molecular systems in a wide range of energy scales ranging from core-binding to valence. This paper discusses the results obtained with the real-time equation-of-motion CC cumulant (RT-EOM-CC) approach and CC Green’s function (CCGF) approaches in applications to the water and water dimer molecules. Here, we compare the ionization potentials obtained with these methods for the valence region with the results obtained with the coupled-cluster with singles, doubles, and perturbative triples formulation as a difference of energies for N and N – 1 electron systems. All methods show good agreement with each other. They also agree well with the experiment with errors usually below 0.1 eV for the ionization potentials. We also analyze unique features of the spectral functions, associated with the position of satellite peaks, obtained with the RT-EOM-CC and CCGF methods employing single and double excitations, as a function of the monomer OH bond length and the proton transfer coordinate in the dimer. Finally, we analyze the impact of the basis set effects on the quality of calculated ionization potentials and find that the basis set effects are less pronounced for the augmented-type sets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Josephson current in Fe-based superconducting junctions: Theory and experiment

Here we present a theory of the dc Josephson effect in contacts between Fe-based and spin-singlet s -wave superconductors. The method is based on the calculation of temperature Green's function in the junction within the tight-binding model. We calculate the phase dependencies of the Josephson current for different orientations of the junction relative to the crystallographic axes of Fe-based superconductor. Further, we consider the dependence of the Josephson current on the thickness of an insulating layer and on temperature. Experimental data for PbIn / Ba 1 – x K x ( FeAs ) 2 point-contact Josephson junctions are consistent with theoretical predictions for s ± symmetry of an order parameter in this material. The proposed method can be further applied to calculations of the dc Josephson current in contacts with other new unconventional multiorbital superconductors, such as Sr 2 RuO 4 and the superconducting topological insulator Cu x Bi 2 Se 3 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Discrete Green’s functions and spectral graph theory for computationally efficient thermal modeling

Here, this work concerns solutions of the heat equation with the spectral graph method, for which the temperature is defined at discrete points in the domain and the spatial relationship among the points is described by a graph. The heat equation on the graph is solved using matrix techniques involving the eigenvectors and eigenvalues of the Laplacian matrix. The spectral graph approach precludes the computationally intensive meshing and numerous time-integration steps of the finite element method. In the present work, the spectral graph method is extended to include heat loss at the boundaries with a generalized boundary condition, and physics-based edge weights are introduced which simplify the calibration process. From this approach a discrete Green’s function is defined which allows for solutions under a variety of heating conditions including: space-varying initial conditions; time-and-space varying internal heating; and, time-and-space-varying heating at boundaries of type 1 (Dirichlet), type 2 (Neumann) and type 3 (Robin). Results are provided for benchmark heat transfer problems in one spatial dimension and in three spatial dimensions, and verification is provided by comparison with exact analytical solutions and finite difference solutions. The spectral graph method converges within 0.4% error of the analytical solution. The practical utility of the approach is demonstrated by thermal simulation of a multilayer additive manufacturing process. The spectral graph results are compared to experimentally-obtained temperature data for two metal parts, with error less than 5% of the experimental measurements, with computation time less than one minute on a desktop computer.

36 MATERIALS SCIENCE↗

Multiscale Neural Networks for Approximating Green’s Functions

Neural networks (NNs) have been widely used to solve partial differential equations (PDEs) in the applications of physics, biology, and engineering. One effective approach for solving PDEs with a fixed differential operator is learning Green’s functions. However, Green’s functions are notoriously difficult to learn due to their poor regularity, which typically requires larger NNs and longer training times. In this work, we address these challenges by leveraging multiscale NNs to learn Green’s functions. Through theoretical analysis using multiscale Barron space methods and experimental validation, we show that the multiscale approach significantly reduces the necessary NN size and accelerates training.

97 MATHEMATICS AND COMPUTING↗

Application of boundary integral method to elastic analysis of V-notched beams

A semidirect boundary integral method, using Airy's stress function and its derivatives in Green's boundary integral formula, is used to obtain an accurate numerical solution for elastic stress and strain fields in V-notched beams in pure bending. The proper choice of nodal spacing on the boundary is shown to be necessary to achieve an accurate stress field in the vicinity of the tip of the notch. Excellent agreement is obtained with the results of the collocation method of solution.

Rzasnicki, W.↗

Efficient three-dimensional Poisson solvers in open rectangular conducting pipe

Three-dimensional (3D) Poisson solver plays an important role in the study of space-charge effects on charged particle beam dynamics in particle accelerators. In this paper, we propose three new 3D Poisson solvers for a charged particle beam in an open rectangular conducting pipe. These three solvers include a spectral integrated Green function (IGF) solver, a 3D spectral solver, and a 3D integrated Green function solver. These solvers effectively handle the longitudinal open boundary condition using a finite computational domain that contains the beam itself. This saves the computational cost of using an extra larger longitudinal domain in order to set up an appropriate finite boundary condition. Using an integrated Green function also avoids the need to resolve rapid variation of the Green function inside the beam. The numerical operational cost of the spectral IGF solver and the 3D IGF solver scales as , where is the number of grid points. The cost of the 3D spectral solver scales as , where is the maximum longitudinal mode number. We compare these three solvers using several numerical examples and discuss the advantageous regime of each solver in the physical application.

Computer Science↗

Large Scale GW Calculations

We present GW calculations of molecules, ordered and disordered solids and interfaces, which employ an efficient contour deformation technique for frequency integration and do not require the explicit evaluation of virtual electronic states nor the inversion of dielectric matrices. We also present a parallel implementation of the algorithm which takes advantage of separable expressions of both the single particle Green's function and the screened Coulomb interaction. The method can be used starting from density functional theory calculations performed with semilocal or hybrid functionals. We applied the newly developed technique to GW calculations of systems of unprecedented size, including water/semiconductor interfaces with thousands of electrons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The effect of Coulomb interactions on thermoelectric properties of quantum dots

Thermoelectric effects in a quantum dot coupled to the source and drain charge reservoirs are explored using a nonequilibrium Green's functions formalism beyond the Hartree-Fock approximation. Thermal transport is analyzed within a linear response regime. A transition from Coulomb blockade regime to Kondo regime in thermoelectric transport through a single-level quantum dot is traced using unified approximations for the relevant Green's functions.

36 MATERIALS SCIENCE↗

The electronic heat capacity of YBa{sub 2}Cu{sub 3}O{sub 7−δ} superconductor

The contributions due to the point defects or disorder and anharmonicities which play deterministic role in the understanding of electronic heat capacity (EHC) of high temperature superconductors (HTS) have been investigated via electron density of states (EDOS) approach on the basis of quantum dynamical many body theory. The evaluation of EDOS has been carried out with the help of most versatile method of double time temperature dependent electron Green’s functions (GF) via a Hamiltonian (non BCS type) which includes the effects of electrons, phonons, defects, anharmonicity, and electron-phonon interactions which enables to account the effects of cubic anharmonicity besides with both the force constant changes and mass difference caused by the impurities in developing the results for EDOS and EHC. The new results reveal some striking features of EHC of HTS.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Global centroid moment tensor solutions in a heterogeneous earth: the CMT3D catalogue

SUMMARY For over 40 yr, the global centroid-moment tensor (GCMT) project has determined location and source parameters for globally recorded earthquakes larger than magnitude 5.0. The GCMT database remains a trusted staple for the geophysical community. Its point-source moment-tensor solutions are the result of inversions that model long-period observed seismic waveforms via normal-mode summation for a 1-D reference earth model, augmented by path corrections to capture 3-D variations in surface wave phase speeds, and to account for crustal structure. While this methodology remains essentially unchanged for the ongoing GCMT catalogue, source inversions based on waveform modelling in low-resolution 3-D earth models have revealed small but persistent biases in the standard modelling approach. Keeping pace with the increased capacity and demands of global tomography requires a revised catalogue of centroid-moment tensors (CMT), automatically and reproducibly computed using Green's functions from a state-of-the-art 3-D earth model. In this paper, we modify the current procedure for the full-waveform inversion of seismic traces for the six moment-tensor parameters, centroid latitude, longitude, depth and centroid time of global earthquakes. We take the GCMT solutions as a point of departure but update them to account for the effects of a heterogeneous earth, using the global 3-D wave speed model GLAD-M25. We generate synthetic seismograms from Green's functions computed by the spectral-element method in the 3-D model, select observed seismic data and remove their instrument response, process synthetic and observed data, select segments of observed and synthetic data based on similarity, and invert for new model parameters of the earthquake’s centroid location, time and moment tensor. The events in our new, preliminary database containing 9382 global event solutions, called CMT3D for ‘3-D centroid-moment tensors’, are on average 4 km shallower, about 1 s earlier, about 5 per cent larger in scalar moment, and more double-couple in nature than in the GCMT catalogue. We discuss in detail the geographical and statistical distributions of the updated solutions, and place them in the context of earlier work. We plan to disseminate our CMT3D solutions via the online ShakeMovie platform.

58 GEOSCIENCES↗

Lattice vibrations in the Frenkel-Kontorova model. I. Phonon dispersion, number density, and energy

We studied the lattice vibrations of two inter-penetrating atomic sublattices via the Frenkel-Kontorova (FK) model of a linear chain of harmonically interacting atoms subjected to an on-site potential, using the technique of thermodynamic Green's functions based on quantum field-theoretical methods. General expressions were deduced for the phonon frequency-wave-vector dispersion relations, number density, and energy of the FK model system. In addition, as the application of the theory, we investigated in detail cases of linear chains with various periods of the on-site potential of the FK model. Some unusual but interesting features for different amplitudes of the on-site potential of the FK model are discussed. In the commensurate structure, the phonon spectrum always starts at a finite frequency, and the gaps of the spectrum are true ones with a zero density of modes. In the incommensurate structure, the phonon spectrum starts from zero frequency, but at a non-zero wave vector; there are some modes inside these gap regions, but their density is very low. In our approximation, the energy of a higher-order commensurate state of the one-dimensional system at a finite temperature may become indefinitely close to the energy of an incommensurate state. This finding implies that the higher-order incommensurate-commensurate transitions are continuous ones and that the phase transition may exhibit a “devil's staircase” behavior at a finite temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Radiation from highly relativistic geodesics

A number of recent works are reviewed concerning the generation and emission of gravitational waves. It is shown that at high frequencies, the generation of gravitational radiation is a local phenomenon. Two examples are described illustrating this generation when a high-energy particle collides against the space-time curvature. One, after Matzner and Nutku, uses a method of virtual photons; the other, after Chrzanowski and Misner, is based on the W.K.B. approximation, corresponding to geometric optics, for the inhomogeneous wave equation. This method uses a factorized integral representation of the Green function which is valid asymptotically to infinity in space.

Misner, C. W.↗

Structural sensitivity analysis: Methods, applications and needs

Innovative techniques applicable to sensitivity analysis of discretized structural systems are reviewed. The techniques include a finite difference step size selection algorithm, a method for derivatives of iterative solutions, a Green's function technique for derivatives of transient response, simultaneous calculation of temperatures and their derivatives, derivatives with respect to shape, and derivatives of optimum designs with respect to problem parameters. Computerized implementations of sensitivity analysis and applications of sensitivity derivatives are also discussed. Some of the critical needs in the structural sensitivity area are indicated along with plans for dealing with some of those needs.

Adelman, H. M.↗

Structural sensitivity analysis: Methods, applications, and needs

Some innovative techniques applicable to sensitivity analysis of discretized structural systems are reviewed. These techniques include a finite-difference step-size selection algorithm, a method for derivatives of iterative solutions, a Green's function technique for derivatives of transient response, a simultaneous calculation of temperatures and their derivatives, derivatives with respect to shape, and derivatives of optimum designs with respect to problem parameters. Computerized implementations of sensitivity analysis and applications of sensitivity derivatives are also discussed. Finally, some of the critical needs in the structural sensitivity area are indicated along with Langley plans for dealing with some of these needs.

Adelman, H. M.↗

Electronic characterization of defects in narrow gap semiconductors

The study of point defects in semiconductors has a long and honorable history. In particular, the detailed understanding of shallow defects in common semiconductors traces back to the classic work of Kohn and Luttinger. However, the study of defects in narrow gap semiconductors represents a much less clear story. Here, both shallow defects (caused by long range potentials) and deep defects (from short range potentials) are far from being completely understood. In this study, all results are calculational and our focus is on the chemical trend of deep levels in narrow gap semiconductors. We study substitutional (including antisite), interstitial and ideal vacancy defects. For substitutional and interstitial impurities, the efects of relaxation are included. For materials like Hg(1-x)Cd(x)Te, we study how the deep levels vary with x, of particular interest is what substitutional and interstitial atoms yield energy levels in the gap i.e. actually produce deep ionized levels. Also, since the main technique utilized is Green's functions, we include some summary of that method.

Patterson, James D.↗

Unsteady three-dimensional thermal field prediction in turbine blades using nonlinear BEM

A time-and-space accurate and computationally efficient fully three dimensional unsteady temperature field analysis computer code has been developed for truly arbitrary configurations. It uses boundary element method (BEM) formulation based on an unsteady Green's function approach, multi-point Gaussian quadrature spatial integration on each panel, and a highly clustered time-step integration. The code accepts either temperatures or heat fluxes as boundary conditions that can vary in time on a point-by-point basis. Comparisons of the BEM numerical results and known analytical unsteady results for simple shapes demonstrate very high accuracy and reliability of the algorithm. An example of computed three dimensional temperature and heat flux fields in a realistically shaped internally cooled turbine blade is also discussed.

Martin, Thomas J.↗

Optical Control of Intersubband Absorption in a Multiple Quantum Well-Embedded Semiconductor Microcravity

Optical intersubband response of a multiple quantum well (MQW)-embedded microcavity driven by a coherent pump field is studied theoretically. The n-type doped MQW structure with three subbands in the conduction band is sandwiched between a semi-infinite medium and a distributed Bragg reflector (DBR). A strong pump field couples the two upper subbands and a weak field probes the two lower subbands. To describe the optical response of the MQW-embedded microcavity, we adopt a semi-classical nonlocal response theory. Taking into account the pump-probe interaction, we derive the probe-induced current density associated with intersubband transitions from the single-particle density-matrix formalism. By incorporating the current density into the Maxwell equation, we solve the probe local field exactly by means of Green's function technique and the transfer-matrix method. We obtain an exact expression for the probe absorption coefficient of the microcavity. For a GaAs/Al(sub x)Ga(sub 1-x)As MQW structure sandwiched between a GaAs/AlAs DBR and vacuum, we performed numerical calculations of the probe absorption spectra for different parameters such as pump intensity, pump detuning, and cavity length. We find that the probe spectrum is strongly dependent on these parameters. In particular, we find that the combination of the cavity effect and the Autler-Townes effect results in a triplet in the optical spectrum of the MQW system. The optical absorption peak value and its location can be feasibly controlled by varying the pump intensity and detuning.

Liu, Ansheng↗

Quadratic time dependent Hamiltonians and separation of variables

Time dependent quantum problems defined by quadratic Hamiltonians are solved using canonical transformations. The Green’s function is obtained and a comparison with the classical Hamilton–Jacobi method leads to important geometrical insights like exterior differential systems, Monge cones and time dependent Gaussian metrics. The Wei–Norman approach is applied using unitary transformations defined in terms of generators of the associated Lie groups, here the semi-direct product of the Heisenberg group and the symplectic group. A new explicit relation for the unitary transformations is given in terms of a finite product of elementary transformations. The sequential application of adequate sets of unitary transformations leads naturally to a new separation of variables method for time dependent Hamiltonians, which is shown to be related to the Inönü–Wigner contraction of Lie groups. The new method allows also a better understanding of interacting particles or coupled modes and opens an alternative way to analyze topological phases in driven systems. - Highlights: • Exact unitary transformation reducing time dependent quadratic quantum Hamiltonian to zero. • New separation of variables method and simultaneous uncoupling of modes. • Explicit examples of transformations for one to four dimensional problems. • New general evolution equation for quadratic form in the action, respectively Green’s function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗