Search NASA⌕ Search

SEARCH · Search NASA

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 73 records · Page 4

Thermal Edge-Effects Model for Automated Tape Placement of Thermoplastic Composites

Two-dimensional thermal models for automated tape placement (ATP) of thermoplastic composites neglect the diffusive heat transport that occurs between the newly placed tape and the cool substrate beside it. Such lateral transport can cool the tape edges prematurely and weaken the bond. The three-dimensional, steady state, thermal transport equation is solved by the Green's function method for a tape of finite width being placed on an infinitely wide substrate. The isotherm for the glass transition temperature on the weld interface is used to determine the distance inward from the tape edge that is prematurely cooled, called the cooling incursion Delta a. For the Langley ATP robot, Delta a = 0.4 mm for a unidirectional lay-up of PEEK/carbon fiber composite, and Delta a = 1.2 mm for an isotropic lay-up. A formula for Delta a is developed and applied to a wide range of operating conditions. A surprise finding is that Delta a need not decrease as the Peclet number Pe becomes very large, where Pe is the dimensionless ratio of inertial to diffusive heat transport. Conformable rollers that increase the consolidation length would also increase Delta a, unless other changes are made, such as proportionally increasing the material speed. To compensate for premature edge cooling, the thermal input could be extended past the tape edges by the amount Delta a. This method should help achieve uniform weld strength and crystallinity across the width of the tape.

Costen, Robert C.↗

Coupling of Multiple Coulomb Scattering with Energy Loss and Straggling in HZETRN

The new version of the HZETRN deterministic transport code based on Green's function methods, and the incorporation of ground-based laboratory boundary conditions, has lead to the development of analytical and numerical procedures to include off-axis dispersion of primary ion beams due to small-angle multiple Coulomb scattering. In this paper we present the theoretical formulation and computational procedures to compute ion beam broadening and a methodology towards achieving a self-consistent approach to coupling multiple scattering interactions with ionization energy loss and straggling. Our initial benchmark case is a 60 MeV proton beam on muscle tissue, for which we can compare various attributes of beam broadening with Monte Carlo simulations reported in the open literature.

Mertens, Christopher J.↗

An Empirical Temperature Variance Source Model in Heated Jets

An acoustic analogy approach is implemented that models the sources of jet noise in heated jets. The equivalent sources of turbulent mixing noise are recognized as the differences between the fluctuating and Favre-averaged Reynolds stresses and enthalpy fluxes. While in a conventional acoustic analogy only Reynolds stress components are scrutinized for their noise generation properties, it is now accepted that a comprehensive source model should include the additional entropy source term. Following Goldstein s generalized acoustic analogy, the set of Euler equations are divided into two sets of equations that govern a non-radiating base flow plus its residual components. When the base flow is considered as a locally parallel mean flow, the residual equations may be rearranged to form an inhomogeneous third-order wave equation. A general solution is written subsequently using a Green s function method while all non-linear terms are treated as the equivalent sources of aerodynamic sound and are modeled accordingly. In a previous study, a specialized Reynolds-averaged Navier-Stokes (RANS) solver was implemented to compute the variance of thermal fluctuations that determine the enthalpy flux source strength. The main objective here is to present an empirical model capable of providing a reasonable estimate of the stagnation temperature variance in a jet. Such a model is parameterized as a function of the mean stagnation temperature gradient in the jet, and is evaluated using commonly available RANS solvers. The ensuing thermal source distribution is compared with measurements as well as computational result from a dedicated RANS solver that employs an enthalpy variance and dissipation rate model. Turbulent mixing noise predictions are presented for a wide range of jet temperature ratios from 1.0 to 3.20.

Khavaran, Abbas↗

Angle-dependent Green's functions for relativistic Compton reflection

Using analytic methods, we derived the first-order Green's functions, in terms of the observing angle, for the reflected spectra after a single Compton scattering of relativistic photons by a plane-parallel semiinfinite medium of cold electrons. We also calculated the total Green's functions for incident photons of 511 keV, using Monte Carlo methods. We show that both the first-order and total Green's functions have strong angular dependences and that the first-order Green's functions play a dominant role in determining the properties of the total reflected spectra.

Hua, Xin-Min↗

Elementary introduction to the Green's function

A technique, using the method of variation of parameters for solving differential equations, is developed for introducing Green's functions early in an undergraduate curriculum. Various examples are presented.

Whitten, R. C.↗

Full CI benchmark calculations for several states of the same symmetry

Full CI (FCI) wave functions are used to compute energies for several electronic states of the same symmetry for SiH2, CH2, and CH2(+). It is found that CASSCF/multireference CI wave functions yield results very similar to FCI, irrespective of whether the CASSCF MOs are optimized independently for each state or using an average of the CASSCF energies for all desired states. The ionization potentials and excitation energies obtained from the FCI calculations should help calibrate methods (such as Green's function approaches, equations of motion and propagator methods, and cluster expansions) in which energy differences are computed directly.

Bauschlicher, Charles W., Jr.↗

Note on the eigensolution of a homogeneous equation with semi-infinite domain

The 'variation-iteration' method using Green's functions to find the eigenvalues and the corresponding eigenfunctions of a homogeneous Fredholm integral equation is employed for the stability analysis of fluid hydromechanics problems with a semiinfinite (infinite) domain of application. The objective of the study is to develop a suitable numerical approach to the solution of such equations in order to better understand the full set of equations for 'real-world' flow models. The study involves a search for a suitable value of the length of the domain which is a fair finite approximation to infinity, which makes the eigensolution an approximation dependent on the length of the interval chosen. In the examples investigated y = 1 = a seems to be the best approximation of infinity; for y greater than unity this method fails due to the polynomial nature of Green's functions.

Wadia, A. R.↗

Scattering from perfectly conducting and resistive strips on a grounded dielectric slab

The scattering properties of perfectly conducting and resistive strips are predicted for strips which are located on a dielectric slab backed by a perfectly conducting ground plane. The spectral domain Green's function is used to relate the currents and fields on the strip, and the resulting integral equation is solved using the method of moments. Both TE and TM strips are examined using piecewise linear and pulse subdomain basis functions, respectively, to model the current on the strip. Calculated results are compared with results measured at the NASA Langley Research Center.

PULSE SUBDOMAIN BASIS FUNCTION↗

Acoustic loading in planar networks

The acoustic loading in a complex planar network of ducts is determined by a method in which Green function surface elements are used. The network consists of straight ducts, elbows and branched ducts. A transfer matrix technique is developed in which each duct is treated separately and the matrix of the influence coefficients is transformed to tri-diagonal form allowing efficient inversion.

El-Raheb, M.↗

Sensitivity analysis of discrete structural systems: A survey

Methods for calculating sensitivity derivatives for discrete structural systems are surveyed, primarily covering literature published during the past two decades. Methods are described for calculating derivatives of static displacements and stresses, eigenvalues and eigenvectors, transient structural response, and derivatives of optimum structural designs with respect to problem parameters. The survey is focused on publications addressed to structural analysis, but also includes a number of methods developed in nonstructural fields such as electronics, controls, and physical chemistry which are directly applicable to structural problems. Most notable among the nonstructural-based methods are the adjoint variable technique from control theory, and the Green's function and FAST methods from physical chemistry.

Adelman, H. M.↗

Measurements of large scale disturbances in rectangular wall jets

Detailed measurements of the velocity field and wall surface pressure distribution of a tone-generating wall jet were obtained in the Ames 7- by 10-foot Wind Tunnel and anechoeic chamber. The jet issued from a 0.508 cm by 10.16 cm nozzle adjacent to a 3.81 cm wall, with an initially parabolic velocity profile at the exit. The jet was operated at Reynolds numbers, R, of 4910, 7740, and 9410. At R = 4910 and 7740, the flow appeared as a double-rowed vortex array which formed upstream of the edge and convected downstream without change in structure. At R = 9410, apparent vortex pairing occurred at the edge location. These flow characteristics were observed in both velocity field measurements and flow visualization records. The unsteady vorticity field of the jet was computed from phase-averaged hot-wire anemometer measurements of the velocity field. The associated acoustic field was estimated by the method of Green's functions, and was found to be similar to measurements of the sound field. These observations tend to support theoretical predictions that vortex shedding at the edge tends to cancel the effect of acoustic sources in other parts of the flow field.

Horne, W. C.↗

Sensitivity analysis of discrete structural systems

Methods for calculating sensitivity derivatives for discrete structural systems are surveyed, primarily covering literature published during the past two decades. Methods are described for calculating derivatives of static displacements and stresses, eigenvalues and eigenvectors, transient structural response, and derivatives of optimum structural designs with respect to problem parameters. The survey is focused on publications developed in nonstructural fields such as electronics, controls, and physical chemistry which are directly applicable to structural problems. Most notable among the nonstructural-based methods are the adjoint variable technique from control theory, and the Green's function and FAST methods from physical chemistry.

Adelman, H. M.↗

Vapor concentration measurement with photothermal deflectometry

Theoretical and experimental results for using the photothermal deflection technique to measure vapor species concentration, while minimizing the disturbance of the transport (material) parameters due to vapor heating, are developed and described. In contrast to common practice, the above constraints require using a pump-beam duty cycle of less than 50 percent. The theoretical description of the shortened heating time is based on a step-function formulation of the pumping cycle. The results are obtained as closed-form solutions of the energy equation for many chopping cycles until steady state is reached, by use of a Green's-function method. The Euler formulation of the Fermat principle is used to calculate the deflection angle. The equations are expanded to include the effects of vapor velocity on both the temperature and temperature gradient profiles. The effects of finite (unfocused) pump and probe beams and thermal (Soret) diffusion are also accounted for. Excellent agreement between theory and experiment is obtained.

Banish, R. Michael↗

Electronic structures of HgTe and CdTe surfaces and HgTe/CdTe interfaces

A Green's-function method has been used to study the surface and interface electronic structures of the II-VI compounds HgTe and CdTe. Localized surface and resonance states near the cation-terminated (100) surface of CdTe and the anion-terminated surface of HgTe have been found for the ideal surfaces. The energies and strengths of these surface states are altered by surface perturbations. The bulk states near the surface are drastically modified by the creation of the surface, but the band gaps remain unchanged. Numerical evaluation of the local densities of states at the Gamma and J points shows that, at the (100) interface of HgTe/CdTe, the previously observed surface states are no longer present. However, in the interface region, bulk states of one material penetrate some distance into the other material.

Schick, J. T.↗

Crustal deformation along the San Andreas, California

The goal is to achieve a better understanding of the regional and local deformation and crustal straining processes in western North America, particularly the effects of the San Andreas and nearby faults on the spatial and temporal crustal deformation behavior. Construction of theoretical models based on the mechanics of coupled elastic plate, viscoelastic foundation and large scale crack mechanics provide a rational basis for the interpretation of seismic and aseismic anomalies and expedite efforts in forecasting the stability of plate boundary deformation. Special focus is placed on the three dimensional time dependent surface deformation due to localized slippage in a elastic layer coupled to a visco-elastic substrate. The numerical analysis is based on a 3-D boundary element technique. Extension to visco-elastic coupling demands the derivation of 3-D time dependent Green's function. This method was applied to analyze the viscoelastic surface displacements due to a dislocated embedded patch. Surface uplift as a function of time and position are obtained. Comparisons between surface uplift for long and short dislocated patches are made.

Li, Victor C.↗

A numerical solution method for acoustic radiation from axisymmetric bodies

A new and very efficient numerical method for solving equations of the Helmholtz type is specialized for problems having axisymmetric geometry. It is then demonstrated by application to the classical problem of acoustic radiation from a vibrating piston set in a stationary infinite plane. The method utilizes 'Green's Function Discretization', to obtain an accurate resolution of the waves using only 2-3 points per wave. Locally valid free space Green's functions, used in the discretization step, are obtained by quadrature. Results are computed for a range of grid spacing/piston radius ratios at a frequency parameter, omega R/c(sub 0), of 2 pi. In this case, the minimum required grid resolution appears to be fixed by the need to resolve a step boundary condition at the piston edge rather than by the length scale imposed by the wave length of the acoustic radiation. It is also demonstrated that a local near-field radiation boundary procedure allows the domain to be truncated very near the radiating source with little effect on the solution.

Caruthers, John E.↗

2D Quantum Transport Modeling in Nanoscale MOSFETs

We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions, oxide tunneling and phase-breaking scattering are treated on an equal footing. Electron bandstructure is treated within the anisotropic effective mass approximation. We present the results of our simulations of MIT 25 and 90 nm "well-tempered" MOSFETs and compare them to those of classical and quantum corrected models. The important feature of quantum model is smaller slope of Id-Vg curve and consequently higher threshold voltage. These results are consistent with 1D Schroedinger-Poisson calculations. The effect of gate length on gate-oxide leakage and subthreshold current has been studied. The shorter gate length device has an order of magnitude smaller leakage current than the longer gate length device without a significant trade-off in on-current.

Svizhenko, Alexei↗

2D Quantum Mechanical Study of Nanoscale MOSFETs

With the onset of quantum confinement in the inversion layer in nanoscale MOSFETs, behavior of the resonant level inevitably determines all device characteristics. While most classical device simulators take quantization into account in some simplified manner, the important details of electrostatics are missing. Our work addresses this shortcoming and provides: (a) a framework to quantitatively explore device physics issues such as the source-drain and gate leakage currents, DIBL, and threshold voltage shift due to quantization, and b) a means of benchmarking quantum corrections to semiclassical models (such as density-gradient and quantum-corrected MEDICI). We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions and oxide tunneling are treated on an equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. We present the results of our simulations of MIT 25, 50 and 90 nm "well-tempered" MOSFETs and compare them to those of classical and quantum corrected models. The important feature of quantum model is smaller slope of Id-Vg curve and consequently higher threshold voltage. Surprisingly, the self-consistent potential profile shows lower injection barrier in the channel in quantum case. These results are qualitatively consistent with ID Schroedinger-Poisson calculations. The effect of gate length on gate-oxide leakage and subthreshold current has been studied. The shorter gate length device has an order of magnitude smaller current at zero gate bias than the longer gate length device without a significant trade-off in on-current. This should be a device design consideration.

Svizhenko, Alexei↗