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 163 records · Page 9

Biharmonic spline interpolation of GEOS-3 and Seasat altimeter data

An algorithm is presented for determining the minimum curvature surface passing through a set of nonuniformly spaced data points. The curve is generated as a linear combination of Green functions for the biharmonic operator at each data point, with the amplitudes of the functions adjusted so that the interpolating surfaces passes through each point. The function passing through the points can be regarded as a spline to which point forces are applied, defining the minimum curvature between the points. The technique was used to combine the along track slopes of the GEOS-3 and Seasat altimeter data into a consistent geoid height map of the Caribbean area, covering 0.5 million data points in the process. Sample images are provided and new topographic features that are revealed are discussed.

Sandwell, David T.↗

Design of magnets inside cylindrical superconducting shields

The design of magnets inside closed, cylindrical, superconducting shields is discussed. The Green function is given for the magnetic vector potential for cylindrically symmetric currents inside such a shield. The magnetic field everywhere inside the shield can be obtained from this function, which includes the effects of the induced shield currents exactly. The field is given for a thin solenoid as an example and the convergence of the series solution for this case is discussed. The shield can significantly reduce the strength and improve the homogeneity of a magnet. The improvement in homogeneity is of particular importance in the design of correction coils. These effects, and the maximum field on the shield, are examined for a typical solenoid. The results given are also useful, although not exact, for long shields with one or two open ends.

Rigby, K. W.↗

On the integration of image sources in exact image method of field analysis

Convergence conditions of image integration in the exact image method of field calculation have been investigated, and the method is extended to include more general media than previously considered. It is demonstrated that the integral is well behaved and the method works best when the field is calculated in a medium with less loss. If the medium has more loss, it is shown that the image line may enter the improper half-space and still produce valid results. Means of correctly selecting the integration path branch for the case of crossing branch cuts of the Green function in complex integration planes are proposed.

Lindell, I. V.↗

The use of the FFT for the efficient solution of the problem of electromagnetic scattering by a body of revolution

The enhancement of the computational efficiency of the body of revolution (BOR) scattering problem is discused with a view to making it practical for solving large-body problems. The problem of EM scattering by a perfectly conducting BOR is considered, although the methods can be extended to multilayered dielectric bodies as well. Typically, the generation of the elements of the moment method matrix consumes a major portion of the computational time. It is shown how this time can be significantly reduced by manipulating the expression for the matrix elements to permit efficient FFT computation. A technique for extracting the singularity of the Green function that appears within the integrands of the matrix diagonal is also presented, further enhancing the usefulness of the FFT. The computation time can thus be improved by at least an order of magnitude for large bodies in comparison to that for previous algorithms.

Gedney, Stephen D.↗

A boundary integral approach to the scattering of nonplanar acoustic waves by rigid bodies

The acoustic scattering of an incident wave by a rigid body can be described by a singular Fredholm integral equation of the second kind. This equation is derived by solving the wave equation using generalized function theory, Green's function for the wave equation in unbounded space, and the acoustic boundary condition for a perfectly rigid body. This paper will discuss the derivation of the wave equation, its reformulation as a boundary integral equation, and the solution of the integral equation by the Galerkin method. The accuracy of the Galerkin method can be assessed by applying the technique outlined in the paper to reproduce the known pressure fields that are due to various point sources. From the analysis of these simpler cases, the accuracy of the Galerkin solution can be inferred for the scattered pressure field caused by the incidence of a dipole field on a rigid sphere. The solution by the Galerkin technique can then be applied to such problems as a dipole model of a propeller whose pressure field is incident on a rigid cylinder. This is the groundwork for modeling the scattering of rotating blade noise by airplane fuselages.

Gallman, Judith M.↗

Energetics of molecular-beam epitaxy models

The removal energies of constituent atoms from various unreconstructed semiconductor surfaces are calculated using a Green function method. An efficient difference-equation approach within the second-neighbor tight-binding model is employed. For a compound AB, binding energies for the A and B atoms on the (111), (-1 -1 -1), (100), and (110) surfaces are calculated. Analyses are made of the energy to remove an atom from the nearly full surface and from the nearly empty surface. Results are presented for Si, GaAs, CdTe, and HgTe; and the surface sublimation energies are found to depend on surface coverage and do not display a simple linear relationship to the number of bonds broken, as is often assumed in modeling growth by MBE.

Krishnamurthy, Srinivasan↗

Solutions of the heat conduction equation in multilayers for photothermal deflection experiments

Analytical expressions for temperature and laser beam deflection in multilayer medium is derived using Green function techniques. The approach is based on calculation of the normal component of heat fluxes across the boundaries, from which either the beam deflections or the temperature anywhere in space can be found. A general expression for the measured signals for the case of four-quadrant detection is also presented and compared with previous calculations of detector response for finite probe beams.

Mcgahan, William A.↗

Nonlinear mechanics of composite materials with periodic microstructure

This report summarizes the result of research done under NASA NAG3-882 Nonlinear Mechanics of Composites with Periodic Microstructure. The effort involved the development of non-finite element methods to calculate local stresses around fibers in composite materials. The theory was developed and some promising numerical results were obtained. It is expected that when this approach is fully developed, it will provide an important tool for calculating local stresses and averaged constitutive behavior in composites. NASA currently has a major contractual effort (NAS3-24691) to bring the approach developed under this grant to application readiness. The report has three sections. One, the general theory that appeared as a NASA TM, a second section that gives greater details about the theory connecting Greens functions and Fourier series approaches, and a final section shows numerical results.

Jordan, E. H.↗

Radiation and scattering by cavity-backed antennas on a circular cylinder

Conformal arrays are popular antennas for aircraft and missile platforms due to their inherent low weight and drag properties. However, to date there has been a dearth of rigorous analytical or numerical solutions to aid the designer. In fact, it has been common practice to use limited measurements and planar approximations in designing such non-planar antennas. The finite element-boundary integral method is extended to scattering and radiation by cavity-backed structures in an infinite, metallic cylinder. In particular, the formulation specifics such as weight functions, dyadic Green's function, implementation details, and particular difficulties inherent to cylindrical structures are discussed. Special care is taken to ensure that the resulting computer program has low memory demand and minimal computational requirements. Both scattering and radiation parameters are computed and validated as much as possible.

Kempel, Leo C.↗

Scattering by cavity-backed antennas on a circular cylinder

Conformal arrays are popular antennas for aircraft, spacecraft, and land vehicle platforms due to their inherent low weight and drag properties. However, to date there has been a dearth of rigorous analytical or numerical solutions to aid the designer. In fact, it has been common practice to use limited measurements and planar approximations in designing such non-planar antennas. The finite element-boundary integral method is extended to scattering by cavity-backed structures in an infinite, metallic cylinder. In particular, the formulation specifics such as weight functions, dyadic Green's function, implementation details and particular difficulties inherent to cylindrical structures are discussed. Special care is taken to ensure that the resulting computer program has low memory demand and minimal computational requirements. Scattering results are presented and validated as much as possible.

Kempel, Leo C.↗

Computational modeling of properties

A simple model was developed to calculate the electronic transport parameters in disordered semiconductors in strong scattered regime. The calculation is based on a Green function solution to Kubo equation for the energy-dependent conductivity. This solution together with a rigorous calculation of the temperature-dependent chemical potential allows the determination of the dc conductivity and the thermopower. For wide-gap semiconductors with single defect bands, these transport properties are investigated as a function of defect concentration, defect energy, Fermi level, and temperature. Under certain conditions the calculated conductivity is quite similar to the measured conductivity in liquid 2-6 semiconductors in that two distinct temperature regimes are found. Under different conditions the conductivity is found to decrease with temperature; this result agrees with measurements in amorphous Si. Finally the calculated thermopower can be positive or negative and may change sign with temperature or defect concentration.

Franz, Judy R.↗

Computational modeling of properties

A simple model was developed to calculate the electronic transport parameters in disordered semiconductors in strong scattered regime. The calculation is based on a Green function solution to Kubo equation for the energy-dependent conductivity. This solution together with a rigorous calculation of the temperature-dependent chemical potential allows the determination of the dc conductivity and the thermopower. For wise-gap semiconductors with single defect bands, these transport properties are investigated as a function of defect concentration, defect energy, Fermi level, and temperature. Under certain conditions the calculated conductivity is quite similar to the measured conductivity in liquid II-VI semiconductors in that two distinct temperature regimes are found. Under different conditions the conductivity is found to decrease with temperature; this result agrees with measurements in amorphous Si. Finally the calculated thermopower can be positive or negative and may change sign with temperature or defect concentration.

Franz, Judy R.↗

Cloud Structure and Cloud-Radiation Interaction: From Smallest Observable Scales (Less than 10 m) to GCM Grid Scales (More than 100 km)

In this presentation we review the fractal nature of internal cloud structure from cm- to km-scales as captured by in-situ probes during long horizontal penetrations by aircraft. We uncover the non-Poissonian spatial distribution of cloud droplets at submeter scales and confirm scale-invariant behavior for large scales. Based on these structural characteristics, we generate simple fractal cloud models that reproduce statistical scaling properties of real clouds. These stochastic models represent a link between nonlinear science, in general, and cloud-radiation interaction, in particular. Next we run three-dimensional radiative transfer computations on these synthetic fractal clouds and compare the structure of the resulting radiation fields with the known structure of the cloud model and with satellite images of real clouds. The different behaviors observed for small and large-scale variabilities will be discussed in detail. We find that while the large-scale fluctuations of the resulting radiation fields resemble those in the original scale-invariant cloud structure, the radiation at small scales is much smoother than its cloud liquid water counterpart. This violates scale-invariance and produces a scale-break at 0.2-0.5 km that is clearly observed in high-resolution satellite data such as from Landsat. Finally, we show how radiative transfer Green function theory in the photon diffusion limit explains (and predicts) the above phenomena of "radiative smoothing."

Marshak, Alexander↗

Dynamics of Quantum Adiabatic Evolution Algorithm for Number Partitioning

We have developed a general technique to study the dynamics of the quantum adiabatic evolution algorithm applied to random combinatorial optimization problems in the asymptotic limit of large problem size n. We use as an example the NP-complete Number Partitioning problem and map the algorithm dynamics to that of an auxiliary quantum spin glass system with the slowly varying Hamiltonian. We use a Green function method to obtain the adiabatic eigenstates and the minimum exitation gap, gmin = O(n2(sup -n/2)), corresponding to the exponential complexity of the algorithm for Number Partitioning. The key element of the analysis is the conditional energy distribution computed for the set of all spin configurations generated from a given (ancestor) configuration by simultaneous flipping of a fixed number of spins. For the problem in question this distribution is shown to depend on the ancestor spin configuration only via a certain parameter related to the energy of the configuration. As the result, the algorithm dynamics can be described in terms of one-dimensional quantum diffusion in the energy space. This effect provides a general limitation of a quantum adiabatic computation in random optimization problems. Analytical results are in agreement with the numerical simulation of the algorithm.

Smelyanskiy, Vadius↗

Dynamics of Quantum Adiabatic Evolution Algorithm for Number Partitioning

We have developed a general technique to study the dynamics of the quantum adiabatic evolution algorithm applied to random combinatorial optimization problems in the asymptotic limit of large problem size n. We use as an example the NP-complete Number Partitioning problem and map the algorithm dynamics to that of an auxiliary quantum spin glass system with the slowly varying Hamiltonian. We use a Green function method to obtain the adiabatic eigenstates and the minimum excitation gap. g min, = O(n 2(exp -n/2), corresponding to the exponential complexity of the algorithm for Number Partitioning. The key element of the analysis is the conditional energy distribution computed for the set of all spin configurations generated from a given (ancestor) configuration by simultaneous flipping of a fixed number of spins. For the problem in question this distribution is shown to depend on the ancestor spin configuration only via a certain parameter related to 'the energy of the configuration. As the result, the algorithm dynamics can be described in terms of one-dimensional quantum diffusion in the energy space. This effect provides a general limitation of a quantum adiabatic computation in random optimization problems. Analytical results are in agreement with the numerical simulation of the algorithm.

Smelyanskiy, V. N.↗

The Role of Instability Waves in Predicting Jet Noise

Debate over whether linear instability waves play a role in the prediction of jet noise has been going on for many years. Parallel mean flow models, such as the one proposed by Lilley, usually neglect these waves because they cause the solution to become infinite. The present paper solves the true non-parallel acoustic equations for a two-dimensional shear layer by using a vector Greens function and assuming small mean flow spread rate. The results show that linear instability waves must be accounted for in order to construct a proper causal solution to the problem.

Goldstein, Marvin E.↗

Constraining Oceanic Sources and Sinks of CO2 by Inverse Modeling: A First Step Toward Assimilation of In Situ and Remote Data

The primary goal of this collaborative project, lead by Nicholas Gruber at UCLA, is to better quantify the air-sea fluxes of pre-industrial and anthropogenic carbon dioxide CO2 for the global ocean with a newly developed inverse method. The method involves using ocean general circulation models (OGCMs) to characterize the subsurface distribution of tracers injected for different regions of the ocean surface. By linearly combining a series of these basis (or Greens) functions, one should be able to recreate any arbitrary ocean tracer field. In our case, we are using the global database of ocean dissolved inorganic carbon (DIC) from the WOCE/JGOFS global CO2 survey from the 1990s. The observed DIC field has been partitioned into pre-industrial and anthropogenic components using the empirical C* technique. We then solve for an optimal combination of the regional basis functions that best match the two sets observations in a least-squares sense with singular value decomposition. This in turn provides solutions for the regional air-sea CO2 fluxes (pre-industrial and anthropogenic) that are independent of more traditional estimates derived from surface water disequilibria in the partial pressure of CO2 and wind speed gas exchange parameterizations. We are apply the method to a suite of nine different OGCMs to test the robustness of the derived surface fluxes fields to differences in ocean circulation. Together with a collegue Keith Lindsay at NCAR, S. Doney's role in this project was to create the appropriate steady-state and temporally evolving basis functions from the NCAR CSM-1 ocean circulation model. Because of a slow start across the project in defining a common set of experimental protocols, we requested and were given a one year, no cost extension. Both the steady-state and temporally evolving basis functions were delivered to the UCLA group for further analysis. Additional information is included in the original extended abstract.

Doney, Scott↗

Global Analysis of Empirical Relationships Between Annual Climate and Seasonality of NDVI

This study describes the use of satellite data to calibrate a new climate-vegetation greenness function for global change studies. We examined statistical relationships between annual climate indexes (temperature, precipitation, and surface radiation) and seasonal attributes of the AVHRR Normalized Difference Vegetation Index (NDVI) time series for the mid-1980s in order to refine our empirical understanding of intraannual patterns and global abiotic controls on natural vegetation dynamics. Multiple linear regression results using global l(sup o) gridded data sets suggest that three climate indexes: growing degree days, annual precipitation total, and an annual moisture index together can account to 70-80 percent of the variation in the NDVI seasonal extremes (maximum and minimum values) for the calibration year 1984. Inclusion of the same climate index values from the previous year explained no significant additional portion of the global scale variation in NDVI seasonal extremes. The monthly timing of NDVI extremes was closely associated with seasonal patterns in maximum and minimum temperature and rainfall, with lag times of 1 to 2 months. We separated well-drained areas from l(sup o) grid cells mapped as greater than 25 percent inundated coverage for estimation of both the magnitude and timing of seasonal NDVI maximum values. Predicted monthly NDVI, derived from our climate-based regression equations and Fourier smoothing algorithms, shows good agreement with observed NDVI at a series of ecosystem test locations from around the globe. Regions in which NDVI seasonal extremes were not accurately predicted are mainly high latitude ecosystems and other remote locations where climate station data are sparse.

Potter, C. S.↗