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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Development of deterministic transport methods for low energy neutrons for shielding in space

Transport of low energy neutrons associated with the galactic cosmic ray cascade is analyzed in this dissertation. A benchmark quality analytical algorithm is demonstrated for use with BRYNTRN, a computer program written by the High Energy Physics Division of NASA Langley Research Center, which is used to design and analyze shielding against the radiation created by the cascade. BRYNTRN uses numerical methods to solve the integral transport equations for baryons with the straight-ahead approximation, and numerical and empirical methods to generate the interaction probabilities. The straight-ahead approximation is adequate for charged particles, but not for neutrons. As NASA Langley improves BRYNTRN to include low energy neutrons, a benchmark quality solution is needed for comparison. The neutron transport algorithm demonstrated in this dissertation uses the closed-form Green's function solution to the galactic cosmic ray cascade transport equations to generate a source of neutrons. A basis function expansion for finite heterogeneous and semi-infinite homogeneous slabs with multiple energy groups and isotropic scattering is used to generate neutron fluxes resulting from the cascade. This method, called the FN method, is used to solve the neutral particle linear Boltzmann transport equation. As a demonstration of the algorithm coded in the programs MGSLAB and MGSEMI, neutron and ion fluxes are shown for a beam of fluorine ions at 1000 MeV per nucleon incident on semi-infinite and finite aluminum slabs. Also, to demonstrate that the shielding effectiveness against the radiation from the galactic cosmic ray cascade is not directly proportional to shield thickness, a graph of transmitted total neutron scalar flux versus slab thickness is shown. A simple model based on the nuclear liquid drop assumption is used to generate cross sections for the galactic cosmic ray cascade. The ENDF/B V database is used to generate the total and scattering cross sections for neutrons in aluminum. As an external verification, the results from MGSLAB and MGSEMI were compared to ANISN/PC, a routinely used neutron transport code, showing excellent agreement. In an application to an aluminum shield, the FN method seems to generate reasonable results.

Ganapol, Barry↗
NASA NTRS1993

Scattering by a groove in an impedance plane

An analysis of two-dimensional scattering from a narrow groove in an impedance plane is presented. The groove is represented by a impedance surface and the problem reduces to that of scattering from an impedance strip in an otherwise uniform impedance plane. On the basis of this model, appropriate integral equations are constructed using a form of the impedance plane Green's functions involving rapidly convergent integrals. The integral equations are solved by introducing a single basis representation of the equivalent current on the narrow impedance insert. Both transverse electric (TE) and transverse magnetic (TM) polarizations are treated. The resulting solution is validated by comparison with results from the standard boundary integral method (BIM) and a high frequency solution. It is found that the presented solution for narrow impedance inserts can be used in conjunction with the high frequency solution for the characterization of impedance inserts of any given width.

Bindiganavale, Sunil↗
NASA NTRS1993

Research on Earth's rotation and the effect of atmospheric pressure on vertical deformation and sea level variability

The work done under NASA grant NAG5-485 included modelling the deformation of the earth caused by variations in atmospheric pressure. The amount of deformation near coasts is sensitive to the nature of the oceanic response to the pressure. The PSMSL (Permanent Service for Mean Sea Level) data suggest the response is inverted barometer at periods greater than a couple months. Green's functions were constructed to describe the perturbation of the geoid caused by atmospheric and oceanic loading and by the accompanying load-induced deformation. It was found that perturbation of up to 2 cm are possible. Ice mass balance data was used for continental glaciers to look at the glacial contributions to time-dependent changes in polar motion, the lod, the earth's gravitational field, the position of the earth's center-of-mass, and global sea level. It was found that there can be lateral, non-hydrostatic structure inside the fluid core caused by gravitational forcing from the mantle, from the inner core, or from topography at the core/mantle or inner core/outer core boundaries. The nutational and tidal response of a non-hydrostatic earth with a solid inner core was modeled. Monthly, global tide gauge data from PSMSL was used to look at the 18.6-year ocean tide, the 14-month pole tide, the oceanic response to pressure, the linear trend and inter-annual variability in the earth's gravity field, the global sea level rise, and the effects of post glacial rebound. The effects of mantle anelasticity on nutations, earth tides, and tidal variation in the lod was modeled. Results of this model can be used with Crustal Dynamics observations to look at the anelastic dissipation and dispersion at tidal periods. The effects of surface topography on various components of crustal deformation was also modeled, and numerical models were developed of post glacial rebound.

Wahr, John↗
NASA NTRS1993

A theoretical and numerical study of polarimetric scattering and emission from random rough surfaces with anisotropic directional spectrum

In this paper, theoretical and numerical results of the polarimetric scattering and emission from random rough surfaces with anisotropic directional spectrum are presented for the remote sensing of ocean and soil surfaces. The polarimetric scattered field for rough dielectric surfaces is derived to the second order by the small perturbation method (SPM). It is found that the second-order scattered field is coherent in nature, and its coefficients for different polarizations present the lowest-order corrections to the Fresnel reflection coefficients of the surfaces. In addition, the cross-polarized (HV and VH) components of the coherent fields are reciprocal and not zero for surfaces with anisotropic directional spectrum when the azimuth angle of the incident direction is not aligned with the symmetry directions of surfaces. In order to verify the energy conservation condition of the theoretical results, which is important if the theory is to be applied to the passive polarimetry of rough surfaces, a Monte Carlo simulation is performed to numerically calculate the polarimetric reflectivities of one-dimensional random rough surfaces which are generated with a prescribed power-law spectrum in the spectral domain and transformed to the spatial domain by the FFT. The surfaces simulated by this approach are periodic with the period corresponding to the low-wavenumber cutoff. To calculate the scattering from periodic dielectric surfaces, the authors present a new numerical technique which applies the Floquet theorem to reduce the problem to one period and does not require the evaluation of one-dimensional periodic Green's function used in the conventional method of moment formulation. Once the scattering coefficients are obtained, the polarimetric Stokes vectors for the emission from the random surfaces are then calculated according to the Kirchhoff's law and are illustrated as functions of relative azimuth observation and row directions. The second-order SPM is also reduced to this one-dimensional case, and theoretical results for all four Stokes parameters are in excellent agreement with the numerical results obtained from the Monte Carlo simulation. Moreover, the second order coherent fields are found to be necessary for the theoretical evaluation of the third and fourth Stokes parameters. Without that, the reflectivities of random rough surfaces are significantly over-estimated, and the signs of the third and the fourth Stokes parameters are found to be incorrect, if calculated by using the first-order incoherent fields alone.

Yueh, S. H.↗
NASA NTRS1994

Radiation and scattering from printed antennas on cylindrically conformal platforms

The goal was to develop suitable methods and software for the analysis of antennas on cylindrical coated and uncoated platforms. Specifically, the finite element boundary integral and finite element ABC methods were employed successfully and associated software were developed for the analysis and design of wraparound and discrete cavity-backed arrays situated on cylindrical platforms. This work led to the successful implementation of analysis software for such antennas. Developments which played a role in this respect are the efficient implementation of the 3D Green's function for a metallic cylinder, the incorporation of the fast Fourier transform in computing the matrix-vector products executed in the solver of the finite element-boundary integral system, and the development of a new absorbing boundary condition for terminating the finite element mesh on cylindrical surfaces.

Kempel, Leo C.↗
NASA NTRS1993

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

LET spectra of fragmented 600-MeV/A iron beam on aluminum and polyethylene: A preliminary analysis

An iron beam experiment recently conducted at the Lawrence Berkeley Laboratory Bevalac by Benton et al. provides an opportunity for verifying the new Green's function computer code (GRNTRN) and assessing the related nuclear database. The iron beam with 600 MeV/A at extraction traversed a series of beam transport elements, lead foil, and several triggering devices before impacting the target. Of these, only the 2.24 g/sq cm lead foil and target are considered in the transport analysis with an assumed (inferred from calibration) beam energy of 557 MeV/A. A thin layer of CR-39 plastic nuclear track detector (PNTD) was placed in front of the target to monitor the incident flux and a stack of four PNTD's placed behind to measure the linear energy transfer (LET) distribution of the transported beam. Test data are analyzed for three separate targets: 2 g/sq cm aluminum, 5-cm polyethylene, and 8-cm polyethylene. The two-layer GRNTRN results were mapped into the detector response function for comparison with the measured LET spectra. Reasonable agreement is obtained. Future research and analysis can be improved by using a more accurate isotope set or including other important media which significantly alter the beam. The assumption that the fragmentation cross sections are too small for aluminum is inconclusive because the production of fragments by 2 g/sq cm aluminum can be more significantly affected by the fragment contribution originating from the media in front of the target than that currently considered in the analysis.

Shinn, Judy L.↗
NASA NTRS1994

Strength characteristics and crack growth behavior of a composite with well aligned fibers

Continuous fiber composites have shown tremendous promise in industrial applications. Their microstructures, however, are very complex and in many instances difficult to characterize. In this project, the fracture characteristics of a specially made fiber reinforced composite with different fiber spacing are investigated. The experimental results so far have shown that after an initial transient phase the crack speed reaches a steady phase, i.e., independent of the crack length. Within the steady crack growth phase debonding along the fibers in the bridging zone grows in a self-similar manner. During the steady phase the energy dissipation per cycle is constant. Afterwards, an increase of the energy dissipation is observed that is accompanied by a decrease in crack speed. This latter trend is presumed to be the result of relatively large amounts of energy dissipated in the bulk of the specimen. Using appropriate Green's function and computer simulations, the stress intensity factor at the crack tip is evaluated for various cases of bridging stresses. In this way the effects of specimen size and fiber spacing on the overall fracture behavior of the composite system are analyzed. The steady crack speed and the steady rate of debonding have a similar power dependence on stress level. Dimensional analysis demonstrates that the particular fracture process is not governed by dimensional invariance but on the detailed micromechanisms in the bridging zone.

Botsis, John↗
NASA NTRS1994

Analysis of cylindrical wrap-around and doubly conformal patch antennas by way of the finite element-artificial absorber method

The goal of this project was to develop analysis codes for computing the scattering and radiation of antennas on cylindrically and doubly conformal platforms. The finite element-boundary integral (FE-BI) method has been shown to accurately model the scattering and radiation of cavity-backed patch antennas. Unfortunately extension of this rigorous technique to coated or doubly curved platforms is cumbersome and inefficient. An alternative approximate approach is to employ an absorbing boundary condition (ABC) for terminating the finite element mesh thus avoiding use of a Green's function. A FE-ABC method is used to calculate the radar cross section (RCS) and radiation pattern of a cavity-backed patch antenna which is recessed within a metallic surface. It is shown that this approach is accurate for RCS and antenna pattern calculations with an ABC surface displaced as little as 0.3 lambda from the cavity aperture. These patch antennas may have a dielectric overlay which may also be modeled with this technique.

Volakis, J. L.↗
NASA NTRS1994

Electronic characterization of defects in narrow gap semiconductors

We use a Green's function technique to calculate the position of deep defects in narrow gap semiconductors. We consider substitutional (including antisite), vacancy, and interstitial (self and foreign) deep defects. We also use perturbation theory to look at the effect of nonparabolic bands on shallow defect energies and find nonparabolicity can increase the binding by 10 percent or so. We consider mercury cadmium telluride (MCT), mercury zinc telluride (MZT), and mercury zinc selenide (MZS). For substitutional and interstitial defects we look at the situation with and without relaxation. For substitutional impurities in MCT, MZT, and MZS, we consider x (the concentration of Cd or Zn) in the range 0.1 less than x less than 0.3 and also consider appropriate x so E(sub g) = 0.1 eV for each of the three compounds. We consider several cation site s-like deep levels and anion site p-like levels. For E(sub g) = 0.1 eV, we also consider the effects of relaxation. Similar comments apply to the interstitial deep levels whereas no relaxation is considered for the ideal vacancy model. Relaxation effects can be greater for the interstitial than the substitutional cases. Specific results are given in figures and tables and comparison to experiment is made in a limited number of cases. We find, for example, that I, Se, S, Rn, and N are possible cation site, s-like deep levels in MCT and Zn and Mg are for anion site, p-like levels (both levels for substitutional cases). The corresponding cation and anion site levels for interstitial deep defects are (Au, Ag, Hg, Cd, Cu, Zn) and (N, Ar, O, F). For the substitutional cases we have some examples of relaxation moving the levels into the band gap, whereas for the interstitial case we have examples where relaxation moves it out of the band gap. Future work involves calculating the effects of charge state interaction and seeing the effect of relaxation on vacancy levels.

Patterson, James D.↗
NASA NTRS1995

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 Fe-56 ion with Pb-Al and Pb-(CH2)(x) targets shows reasonable agreement.

Shinn, J. L.↗
NASA NTRS1995

Electronic characterization of defects in narrow gap semiconductors: Comparison of electronic energy levels and formation energies in mercury cadmium telluride, mercury zinc telluride, and mercury zinc selenide

The project has evolved to that of using Green's functions to predict properties of deep defects in narrow gap materials. Deep defects are now defined as originating from short range potentials and are often located near the middle of the energy gap. They are important because they affect the lifetime of charge carriers and hence the switching time of transistors. We are now moving into the arena of predicting formation energies of deep defects. This will also allow us to make predictions about the relative concentrations of the defects that could be expected at a given temperature. The narrow gap materials mercury cadmium telluride (MCT), mercury zinc telluride (MZT), and mercury zinc selenide (MZS) are of interest to NASA because they have commercial value for infrared detecting materials, and because there is a good possibility that they can be grown better in a microgravity environment. The uniform growth of these crystals on earth is difficult because of convection (caused by solute depletion just ahead of the growing interface, and also due to thermal gradients). In general it is very difficult to grow crystals with both radial and axial homogeneity.

Patterson, James D.↗
NASA NTRS1994

Detection of atmospheric pressure loading using very long baseline interferometry measurements

Loading of the Earth by the temporal redistribution of global atmospheric mass is likely to displace the positions of geodetic monuments by tens of millimeters both vertically and horizontally. Estimates of these displacements are determined by convolving National Meteorological Center (NMC) global values of atmospheric surface pressure with Farrell's elastic Green's functions. An analysis of the distances between radio telescopes determined by very long baseline interferometry (VLBI) between 1984 and 1992 reveals that in many of the cases studied there is a significant contribution to baseline length change due to atmospheric pressure loading. Our analysis covers intersite distances of between 1000 and 10,000 km and is restricted to those baselines measured more than 100 times. Accounting for the load effects (after first removing a best fit slope) reduces the weighted root-mean-square (WRMS) scatter of the baseline length residuals on 11 of the 22 baselines investigated. The slight degradation observed in the WRMS scatter on the remaining baselines is largely consistent with the expected statistical fluctuations when a small correction is applied to a data set having a much larger random noise. The results from all baselines are consistent with approximately 60% of the computed pressure contribution being present in the VLBI length determinations. Site dependent coefficients determined by fitting local pressure to the theoretical radial displacement are found to reproduce the deformation caused by the regional pressure to within 25% for most inland sites. The coefficients are less reliable at near coastal and island stations.

Vandam, T. M.↗
NASA NTRS1994

Transverse particle acceleration and diffusion in a planetary magnetic field

A general model of particle acceleration by plasma waves coupled with adiabatic radial diffusion in a planetary magnetic field is developed. The model assumes that a spectrum of lower hybird waves is present to resonantly accelerate ions transverse to the magnetic field. The steady state Green's function for the combined radial diffusion and wave acceleration equation is found in terms of a series expansion. The results provide a rigorous demonstration of how a quasi-Maxwellian distribution function is formed in the absence of particle collisons and elucidate the nature of turbulent heating of magnetospheric plasmas. The solution is applied to the magnetosphere of Neptune for which a number of examples are given illustrating how the spectrum of pickup N(+) ions from Triton evolves.

Barbosa, D. D.↗
NASA NTRS1995

Theoretical overview and modeling of the sodium and potassium atmospheres of mercury

A general theoretical overview for the sources, sinks, gas-surface interactions, and transport dynamics of sodium and potassium in the exospheric atmsophere of Mercury is given. Information for these four factors, which control the spatial distribution of these two alkali-group gases about the planet, is incorporated in numerical models. The spatial nature and relative importance of the initial source atom atmosphere and the ambient (ballistic hopping) atom atmosphere are then examined and are shown to be controlled and coupled to a great extent by the extremely large and variable solar radiation acceleration experienced by sodium and potassium as they resonantly scatter solar photons. The lateral (antisunward) transport rate of thermally accommodated sodium and potassium ambient atoms is shown to be driven by the solar radiation acceleration and, over a significant portion of Mercury's orbit about the Sun, is sufficiently rapid to be competitive with the short photoionization lifetimes for these atoms when they are located on the summit surface near or within about 30 deg of the terminator. The lateral transport rate is characterized by a migration time determined by model calculations for an ensemble of atoms initially starting at a point source on the surface (i.e., a numerical spacetime dependent Green's function). Four animations for the spacetime evolution of the sodium (or potassium) atmosphere produced by a point source on the surface are presented on a videotape format. For extended surface sources for sodium and potassium, the local column density is determined by competition between the photoionization lifetimes and the lateral transport times of atoms originating from different surface source locations. Sodium surface source fluxes (referenced to Mercury at perihelion) that are required on the sunlit hemisphere to reproduce the typically observed several megarayleighs of D2 emission-line brightness and the inferred column densities of 1-2 x 10(exp 11) atoms per sq cm range from approximately 2-5 x 10(exp 7) atoms/sq cm/sec. The sodium model is applied to study observational data that document an anticorrelation in the average sodium column density and solar radiation acceleration. Lateral transport driven by the solar radiation acceleration is shown to produce this behavior for combinations of different sources and surface accomodation coefficients. The best fit model fits to the observational data require a significant degree of thermal accommodation of the ambient sodium atoms to the surface and a source rate that decreases as an inverse power of 1.5 to 2 in heliocentric distance.

Smyth, William H.↗
NASA NTRS1994

Atmospheric pressure loading parameters from very long baseline interferometry observations

Atmospheric mass loading produces a primarily vertical displacement of the Earth's crust. This displacement is correlated with surface pressure and is large enough to be detected by very long baseline interferometry (VLBI) measurements. Using the measured surface pressure at VLBI stations, we have estimated the atmospheric loading term for each station location directly from VLBI data acquired from 1979 to 1992. Our estimates of the vertical sensitivity to change in pressure range from 0 to -0.6 mm/mbar depending on the station. These estimates agree with inverted barometer model calculations (Manabe et al., 1991; vanDam and Herring, 1994) of the vertical displacement sensitivity computed by convolving actual pressure distributions with loading Green's functions. The pressure sensitivity tends to be smaller for stations near the coast, which is consistent with the inverted barometer hypothesis. Applying this estimated pressure loading correction in standard VLBI geodetic analysis improves the repeatability of estimated lengths of 25 out of 37 baselines that were measured at least 50 times. In a root-sum-square (rss) sense, the improvement generally increases with baseline length at a rate of about 0.3 to 0.6 ppb depending on whether the baseline stations are close to the coast. For the 5998-km baseline from Westford, Massachusetts, to Wettzell, Germany, the rss improvement is about 3.6 mm out of 11.0 mm. The average rss reduction of the vertical scatter for inland stations ranges from 2.7 to 5.4 mm.

Macmillan, D. S.↗
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.↗
NASA NTRS1995

Electromagnetic Scattering From a Rectangular Cavity Recessed in a 3-D Conducting Surface

The problem of electromagnetic (EM) scattering from an aperture backed by a rectangular cavity recessed in a three-dimensional conducting body is analyzed using the coupled field integral equation approach. Using the free-space Green's function, EM fields scattered outside the cavity are determined in terms of (1) an equivalent electric surface current density flowing on the three-dimensional conducting surface of the object including the cavity aperture and (2) an equivalent magnetic surface current density flowing over the aperture only. The EM fields inside the cavity are determined using the waveguide modal expansion functions. Making the total tangential electric and magnetic fields across the aperture continuous and subjecting the total tangential electric field on the outer conducting three-dimensional surface of the object to zero, a set of coupled integral equations is obtained. The equivalent electric and magnetic surface currents are then obtained by solving the coupled integral equation using the Method of Moments. The numerical results on scattering from rectangular cavities embedded in various three-dimensional objects are compared with the results obtained by other numerical techniques.

Deshpande, M. D.↗