Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Green 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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

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↗
NASA NTRS2001

2D Quantum Transport Modeling in 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, oxide tunneling and phase-breaking scattering are treated on equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. Quantum simulations are focused on MIT 25, 50 and 90 nm "well- tempered" MOSFETs and compared to 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 quantitatively consistent with I D Schroedinger-Poisson calculations. The effect of gate length on gate-oxide leakage and sub-threshold 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↗
NASA NTRS1998

Predictive Models for Semiconductor Device Design and Processing

The device feature size continues to be on a downward trend with a simultaneous upward trend in wafer size to 300 mm. Predictive models are needed more than ever before for this reason. At NASA Ames, a Device and Process Modeling effort has been initiated recently with a view to address these issues. Our activities cover sub-micron device physics, process and equipment modeling, computational chemistry and material science. This talk would outline these efforts and emphasize the interaction among various components. The device physics component is largely based on integrating quantum effects into device simulators. We have two parallel efforts, one based on a quantum mechanics approach and the second, a semiclassical hydrodynamics approach with quantum correction terms. Under the first approach, three different quantum simulators are being developed and compared: a nonequlibrium Green's function (NEGF) approach, Wigner function approach, and a density matrix approach. In this talk, results using various codes will be presented. Our process modeling work focuses primarily on epitaxy and etching using first-principles models coupling reactor level and wafer level features. For the latter, we are using a novel approach based on Level Set theory. Sample results from this effort will also be presented.

Meyyappan, Meyya↗
NASA NTRS2002

Influence of Scattering on Ballistic Nanotransistor Design

Importance of this work: (1) This is the first work to model electron-phonon scattering within a quantum mechanical approach to nanotransistors. The simulations use the non equilibrium Green's function method. (2) A simple equation which captures the importance of scattering as a function of the spatial location from source to drain is presented. This equation helps interpret the numerical simulations. (3) We show that the resistance per unit length in the source side is much larger than in the drain side. Thus making scattering in the source side of the device much more important than scattering in the drain side. Numerical estimates of ballisticity for 10nm channel length devices in the presence of of electron-phonon scattering are given. Based on these calculations, we propose that to achieve a larger on-current in nanotransistors, it is crucial to keep the highly doped source extension region extremely small, even if this is at the cost of making the highly doped drain extension region longer.

Anantram, M. P.↗
NASA NTRS2002

Slat Cove Noise Modeling: A Posteriori Analysis of Unsteady RANS Simulations

A companion paper by Khorrami et al demonstrates the feasibility of simulating the (nominally) self-sustained, large-scale unsteadiness within the leading-edge slat-cove region of multi-element airfoils using unsteady Reynolds-Averaged Navier-Stokes (URANS) equations, provided that the turbulence production term in the underlying two-equation turbulence model is switched off within the cove region. In conjunction with a FfowesWilliams-Hawkings solver, the URANS computations were shown to capture the dominant portion of the acoustic spectrum attributed to slat noise, as well as reproducing the increased intensity of slat cove motions (and, correspondingly, far-field noise as well) at the lower angles of attack. This paper examines that simulation database, augmented by additional simulations, with the objective of transitioning this apparent success to aeroacoustic predictions in an engineering context. As a first step towards this goal, the simulated flow and acoustic fields are compared with experiment and simplified analytical model. Rather intense near-field fluctuations in the simulated flow are found to be associated with unsteady separation along the slat bottom surface, relatively close to the slat cusp. Accuracy of the laminar-cove simulations in this near-wall region is raised to be an open issue. The adjoint Green's function approach is also explored in an attempt to identify the most efficient noise source locations.

Choudhari, Meelan↗
NASA NTRS2002

Where is Scattering Important in Nanotransistors?

This viewgraph presentation provides information on where scattering is most detrimental to on-current in a transistor nanostructure. It also provides information on series resistance, potential, and current. The presentation includes algorithms and models for two dimensional (quasi-one dimensional) simulation of nanostructures, and equations for retarded and less-than Green's functions.

Anantram, M. P.↗
NASA NTRS2003

Where is scattering important in nanotransistors?

This viewgraph presentation provides an overview of scattering in nanotransistors. Topics covered include: scattering mechanisms, devices used and series resistance. The authors conclude: (1) the influence of electron-photon scattering on the on-current of dual gate MOSFET's using the non-equilibrium Green's function method and (2) the influence of the drain extension region can no longer be modeled as a series resistance.

Anantram, M. P.↗
NASA NTRS2003

A New Look into the Effect of Large Drops on Radiative Transfer Process

Recent studies indicate that a cloudy atmosphere absorbs more solar radiation than any current 1D or 3D radiation model can predict. The excess absorption is not large, perhaps 10-15 W/sq m or less, but any such systematic bias is of concern since radiative transfer models are assumed to be sufficiently accurate for remote sensing applications and climate modeling. The most natural explanation would be that models do not capture real 3D cloud structure and, as a consequence, their photon path lengths are too short. However, extensive calculations, using increasingly realistic 3D cloud structures, failed to produce photon paths long enough to explain the excess absorption. Other possible explanations have also been unsuccessful so, at this point, conventional models seem to offer no solution to this puzzle. The weakest link in conventional models is the way a size distribution of cloud particles is mathematically handled. Basically, real particles are replaced with a single average particle. This "ensemble assumption" assumes that all particle sizes are well represented in any given elementary volume. But the concentration of larger particles can be so low that this assumption is significantly violated. We show how a different mathematical route, using the concept of a cumulative distribution, avoids the ensemble assumption. The cumulative distribution has jumps, or steps, corresponding to the rarer sizes. These jumps result in an additional term, a kind of Green's function, in the solution of the radiative transfer equation. Solving the cloud radiative transfer equation with the measured particle distributions, described in a cumulative rather than an ensemble fashion, may lead to increased cloud absorption of the magnitude observed.

Marshak, Alexander↗
NASA NTRS2003

Three-Gorge Reservoir: A 'Controlled Experiment' for Calibration/Validation of Time-Variable Gravity Signals Detected from Space

With the advances of measurements, modern space geodesy has become a new type of remote sensing for the Earth dynamics, especially for mass transports in the geophysical fluids on large spatial scales. A case in point is the space gravity mission GRACE (Gravity Recovery And Climate Experiment) which has been in orbit collecting gravity data since early 2002. The data promise to be able to detect changes of water mass equivalent to sub-cm thickness on spatial scale of several hundred km every month or so. China s Three-Gorge Reservoir has already started the process of water impoundment in phases. By 2009,40 km3 of water will be stored behind one of the world s highest dams and spanning a section of middle Yangtze River about 600 km in length. For the GRACE observations, the Three-Gorge Reservoir would represent a geophysical controlled experiment , one that offers a unique opportunity to do detailed geophysical studies. -- Assuming a complete documentation of the water level and history of the water impoundment process and aided with a continual monitoring of the lithospheric loading response (such as in area gravity and deformation), one has at hand basically a classical forwardinverse modeling problem of surface loading, where the input and certain output are known. The invisible portion of the impounded water, i.e. underground storage, poses either added values as an observable or a complication as an unknown to be modeled. Wang (2000) has studied the possible loading effects on a local scale; we here aim for larger spatial scales upwards from several hundred km, with emphasis on the time-variable gravity signals that can be detected by GRACE and follow-on missions. Results using the Green s function approach on the PREM elastic Earth model indicate the geoid height variations reaching several millimeters on wavelengths of about a thousand kilometers. The corresponding vertical deformations have amplitude of a few centimeters. In terms of long-wavelength spherical harmonics, the induced geoid height variations are very close to the accuracy of GRACE- recoverable gravity field, while the low-degree (2 to 5) harmonics should be detectable. With a large regional time-variable gravity signal, the Three-Gorge experiment can serve as a useful calibration/verification for GRACE (including the elastic loading effects), and future gravity missions (especially for visco-elastic yielding as well as underground water variations).

Chao, Benjamin F.↗
NASA NTRS2003

Spectral Characteristics of Wake Vortex Sound During Roll-Up

This report presents an analysis of the sound spectra generated by a trailing aircraft vortex during its rolling-up process. The study demonstrates that a rolling-up vortex could produce low frequency (less than 100 Hz) sound with very high intensity (60 dB above threshold of human hearing) at a distance of 200 ft from the vortex core. The spectrum then drops o rapidly thereafter. A rigorous analytical approach has been adopted in this report to derive the spectrum of vortex sound. First, the sound pressure was solved from an alternative treatment of the Lighthill s acoustic analogy approach [1]. After the application of Green s function for free space, a tensor analysis was applied to permit the removal of the source term singularity of the wave equation in the far field. Consequently, the sound pressure is expressed in terms of the retarded time that indicates the time history and spacial distribution of the sound source. The Fourier transformation is then applied to the sound pressure to compute its spectrum. As a result, the Fourier transformation greatly simplifies the expression of the vortex sound pressure involving the retarded time, so that the numerical computation is applicable with ease for axisymmetric line vortices during the rolling-up process. The vortex model assumes that the vortex circulation is proportional to the time and the core radius is a constant. In addition, the velocity profile is assumed to be self-similar along the aircraft flight path, so that a benchmark vortex velocity profile can be devised to obtain a closed form solution, which is then used to validate the numerical calculations for other more realistic vortex profiles for which no closed form solutions are available. The study suggests that acoustic sensors operating at low frequency band could be profitably deployed for detecting the vortex sound during the rolling-up process.

Booth, Earl R., Jr.↗
NASA NTRS2004

A-Posteriori Error Estimation for Hyperbolic Conservation Laws with Constraint

This lecture considers a-posteriori error estimates for the numerical solution of conservation laws with time invariant constraints such as those arising in magnetohydrodynamics (MHD) and gravitational physics. Using standard duality arguments, a-posteriori error estimates for the discontinuous Galerkin finite element method are then presented for MHD with solenoidal constraint. From these estimates, a procedure for adaptive discretization is outlined. A taxonomy of Green's functions for the linearized MHD operator is given which characterizes the domain of dependence for pointwise errors. The extension to other constrained systems such as the Einstein equations of gravitational physics are then considered. Finally, future directions and open problems are discussed.

Barth, Timothy↗
NASA NTRS2000

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.↗
NASA NTRS1994

HZE beam transport in multilayered materials

A nonperturbative analytic solution of the high charge and energy (HZE) Green's function is used to implement a computer code for laboratory ion beam transport in multiple-layered materials. The code is established to operate on the Langley nuclear fragmentation model used in space engineering applications. Computational procedures are established to generate linear energy transfer (LET) distributions for a specified ion beam and target for comparison with experimental measurement. Comparison with 56Fe ion with Pb-Al and Pb-(CH2)x targets shows reasonable agreement.

NASA Discipline Radiation Health↗
NASA NTRS2004

Modified Coulomb-Dipole Theory for 2e Photoionization

In the light of recent experiment on 2e photoionization of Li near threshold, we have considered a modification of the Coulomb-dipole theory, retaining the basic assumption that the threshold is dominated by asymmetric events in phase space [implies r(sub 1), k(sub 1)) greater than or equal to 2(r(sub 2), k(sub )]. In this region [in a collinear model, 2/r(sub 12) approached + 2/(r(sub 1)+r(sub 2)] the interaction reduces to V(rIsub 1) is greater than or equal to 2r(sub 2) is identically equal to [(-Z/r(sub 2)-(A-1)/r(sub 1)] + [(-2r(sub 2)/r(sub 1 exp 2)] is identically equal to V(sub c)+[V(sub d)]. For two electron emission Z = 2, thus both electrons see a Coulomb potential (V(sub c)) asymptotically, albeit each seeing a different charge. The residual potential (V(sub d)) is dipole in character. Writing the total psi = psi (sub c) + psi(sub d) = delta psi, and noting that. (T+V(sub c)-E)psy(sub c) = 0 and (T+V(sub c))psi(sub d) = 0 can be solved exactly, we find, substituting psi into the complete Schrod. Eq., that delta psi = -(H-E)(exp -1)(V(sub d) psi(sub 0)+V(sub c psi (sub 1). Using the fact that the absolute value of V(sub c) is much more than the absolute value of V(sub d) in almost all of configuration space, we can replace H by H(sub 0) in 9H-E)(exp -1) to obtain an improved approximation psi (improved) = psi(sub c) + psi(sub d) -(H(sub 0)-E)(exp -1) (V(sub c) psi (sub 0) + V(sub c) psi(sub 1). Here's the Green's function (H(sub 0)-E)(exp -1), can be exhibited explicitly, but the last term in psi (improved) is small, compared to the first two terms. Inserting them into the transition matrix element, which one handles in the usual way, we obtain in the limit E approaches 0, the threshold law: Q(E) alpha E + M(E)E(exp 5/4) + higher order (Eq. 1a). The modulation function, M(E), is a well-defined (but very non-trivial integral, but it is expected to be well approximated by a sinusoidal function containing a dipole phase term (M(E) = c sin[alpha log (E) + micron] (Eq. 1b). Experimental results show definite modulations, and are well fitted by Eqs (1).

Source record↗
NASA NTRS2004

The Role of Instability Waves in Predicting Jet Noise

There has been an ongoing debate about the role of linear instability waves in the prediction of jet noise. Parallel mean flow models, such as the one proposed by Lilley, usually neglect these waves because they cause the solution to become infinite. The resulting solution is then non-causal and can, therefore, be quite different from the true causal solution for the chaotic flows being considered here. The present paper solves the relevant acoustic equations for a non-parallel mean flow by using a vector Green s function approach and assuming the mean flow to be weakly non-parallel, i.e., assuming the spread rate to be small. It demonstrates that linear instability waves must be accounted for in order to construct a proper causal solution to the jet noise problem. . Recent experimental results (e.g., see Tam, Golebiowski, and Seiner,1996) show that the small angle spectra radiated by supersonic jets are quite different from those radiated at larger angles (say, at 90deg) and even exhibit dissimilar frequency scalings (i.e., they scale with Helmholtz number as opposed to Strouhal number). The present solution is (among other things )able to explain this rather puzzling experimental result.

Goldstein, M. E.↗
NASA NTRS1998

A Hybrid SBR/MoM Technique for Analysis of Scattering from Small Protrusions on a Large Conducting Body

For analysis of large-scale electromagnetic scattering problems, high-frequency asymptotic methods are fast but approximate, whereas low-frequency numerical methods are accurate but slow. Neither can produce an efficient and accurate solution to scattering by large bodies containing small structures. A promising approach is to combine the best features of both types of methods to produce a hybrid technique that is sufficiently fast, reasonably accurate, and applicable to a class of unsolvable problems such as the scatterers mentioned above. There are two extremes for this type of hybridization. One is simply to superimpose solutions from asymptotic and numerical methods. While this approach is most widely used in practical applications, it neglects the interactions between the two solutions, which can be significant in many problems. The other extreme is to combine an asymptotic and a numerical method in an exact manner. In this approach, the effect of a large body is included by incorporating its diffraction into the Green's function in the integral equation for the small structures, which accounts for all interactions. While this approach is accurate, it is difficult to be implemented in a general-purpose computer code because of its complex nature. A more practical approach is to develop a technique that can include all significant interactions and neglect all trivial interactions. The resulting hybrid technique can produce sufficient accuracy and can be implemented in a general-purpose computer code. In this paper, we develop a technique that combines the shooting and- bouncing-ray (SBR) method and the method of moments (MoM) to solve for the scattering by large conducting bodies with small structures mounted on their surfaces.

Jin, J. M.↗
NASA NTRS2007

Derivation of Formulations 1 and 1A of Farassat

Formulations 1 and 1A are the solutions of the Ffowcs Williams-Hawkings (FW-H) equation with surface sources only when the surface moves at subsonic speed. Both formulations have been successfully used for helicopter rotor and propeller noise prediction for many years although we now recommend using Formulation 1A for this purpose. Formulation 1 has an observer time derivative that is taken numerically, and thus, increasing execution time on a computer and reducing the accuracy of the results. After some discussion of the Green's function of the wave equation, we derive Formulation 1 which is the basis of deriving Formulation 1A. We will then show how to take this observer time derivative analytically to get Formulation 1A. We give here the most detailed derivation of these formulations. Once you see the whole derivation, you will ask yourself why you did not do it yourself!

Farassat, F.↗