Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CaI”

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 415 records · Page 23

Correction of thin cirrus effects in AVIRIS images using the sensitive 1.375-micron cirrus detecting channel

Using spectral imaging data acquired with the Airborne Visible Infrared Imaging Spectrometer (AVIRIS) from an ER-2 aircraft at 20 km altitude during various field programs, it was found that narrow channels near the center of the strong 1.38-micrometer water vapor band are very effective in detecting think cirrus clouds. Based on this observation from AVIRIS data, Gao and Kaufman proposed to put a channel centered at 1.375 micrometers with a width of 30 nm on the Moderate Resolution Imaging Spectrometer (MODIS) for remote sensing of cirrus clouds from space. The sensitivity of the 1.375-micrometer MODIS channel to detect thin cirrus clouds during the day time is expected to be one to two orders of magnitude better than the current infrared emission techniques. As a result, much larger fraction of the satellite data is expected to be identified as being covered by cirrus clouds, some of them so thin that their obscuration of the surface is very small. In order to make better studies of surface reflectance properties, thin cirrus effects must be removed from satellite images. Therefore, there is a need to study radiative properties of thin cirrus clouds, so that a strategy for correction or removal of the thin cirrus effects, similar to the correction of atmospheric aerosol effect, can be formed. In this extended abstract, we describe an empirical approach for removing/correcting thin cirrus effects in AVIRIS images using channels near 1.375 microns - one step beyond the detection of cirrus clouds using these channels.

Gao, Bo-Cai↗

Remote sensing of smoke, clouds, and radiation using AVIRIS during SCAR experiments

During the past two years, researchers from several institutes joined together to take part in two SCAR experiments. The SCAR-A (Sulfates, Clouds And Radiation - Atlantic) took place in the mid-Atlantic region of the United States in July, 1993. remote sensing data were acquired with the Airborne Visible Infrared Imaging Spectrometer (AVIRIS), the MODIS Airborne Simulator (MAS), and a RC-10 mapping camera from an ER-2 aircraft at 20 km. In situ measurements of aerosol and cloud microphysical properties were made with a variety of instruments equipped on the University of Washington's C-131A research aircraft. Ground based measurements of aerosol optical depths and particle size distributions were made using a network of sunphotometers. The main purpose of SCAR-A experiment was to study the optical, physical and chemical properties of sulfate aerosols and their interaction with clouds and radiation. Sulfate particles are believed to affect the energy balance of the earth by directly reflecting solar radiation back to space and by increasing the cloud albedo. The SCAR-C (Smoke, Clouds And Radiation - California) took place on the west coast areas during September - October of 1994. Sets of aircraft and ground-based instruments, similar to those used during SCAR-A, were used during SCAR-C. Remote sensing of fires and smoke from AVIRIS and MAS imagers on the ER-2 aircraft was combined with a complete in situ characterization of the aerosol and trace gases from the C-131A aircraft of the University of Washington and the Cesna aircraft from the U.S. Forest Service. The comprehensive data base acquired during SCAR-A and SCAR-C will contribute to a better understanding of the role of clouds and aerosols in global change studies. The data will also be used to develop satellite remote sensing algorithms from MODIS on the Earth Observing System.

Gao, Bo-Cai↗

Surface-induced brightness temperature variations and their effects on detecting thin cirrus clouds using IR emission channels in the 8-12 microns region

A method for detecting cirrus clouds in terms of brightness temperature differences between narrowbands at 8, 11, and 12 microns has been proposed by Ackerman et al. In this method, the variation of emissivity with wavelength for different surface targets was not taken into consideration. Based on state-of-the-art laboratory measurements of reflectance spectra of terrestrial materials by Salisbury and D'Aria, it is found that the brightness temperature differences between the 8- and 11-microns bands for soils, rocks, and minerals, and dry vegetation can vary between approximately -8 and +8 K due solely to surface emissivity variations. The large brightness temperature differences are sufficient to cause false detection of cirrus clouds from remote sensing data acquired over certain surface targets using the 8-11-12-microns method directly. It is suggested that the 8-11-12-microns method should be improved to include the surface emissivity effects. In addition, it is recommended that in the future the variation of surface emissivity with wavelength should be taken into account in algorithms for retrieving surface temperatures and low-level atmospheric temperature and water vapor profiles.

Gao, Bo-Cai↗

Momentum transport near a magnetic X line in collisionless reconnection

Plasma dynamics and momentum transport near an X line during time-dependent magnetic reconnection in a collisionless plasma are investigated based on two-dimensional particle simulations. We find that a weakly skewed velocity distribution is formed near the magnetic X line, leading to the presence of off-diagonal elements of the plasma pressure tensor. Let the reconnection electric field be in the y direction. The gradients of the off-diagonal elements of the pressure tensor can provide a transport of the y momentum. During the normal magnetic reconnection, the momentum transport associated with the off-diagonal terms of the pressure tensor mediates a transfer of the y momentum from the region near the X line to regions outside the X line. A period of 'reverse magnetic reconnection,' during which the plasma kinetic energy is converted into magnetic energy, is also observed in the simulation. When reverse reconnection occurs, the gradients of the off-diagonal pressure tensor elements can mediate a transfer of y momentum into the X line. It is found that the inertial term also plays a significant role in the force balance near the magnetic X line. An explanation for the origin of the off-diagonal pressure terms is also given in this paper.

Cai, H. J.↗

Formation of a very thin current sheet in the near-earth magnetotail and the explosive growth phase of substorms

A magnetofricional method is used to construct two-dimensional MHD equilibria of the Earth's magnetosphere for a given distribution of entropy functions(S = pV(exp gamma), where p is the plasma pressure and V is the tube volume per unit magnetic flux. It is found that a very thin current sheet with B (sub zeta) is less than 0.5 nu T and thickness less than 1000 km can be formed in the near-earth magnetotail (x is approximately -8 to -20R(sub e) during the growth phase of substorm. The tail current sheets are found to become thinner as the entropy or the entropy gradient increases. It is suggested that the new entropy anti-diffusion instability associated with plasma transport across field lines leads to magnetic field dipolarization and accelerates the formation of thin current sheet, which may explain the observed explosive growth phase of substorms.

Lee, L. C.↗

Exploring the use of I/O nodes for computation in a MIMD multiprocessor

As parallel systems move into the production scientific-computing world, the emphasis will be on cost-effective solutions that provide high throughput for a mix of applications. Cost effective solutions demand that a system make effective use of all of its resources. Many MIMD multiprocessors today, however, distinguish between 'compute' and 'I/O' nodes, the latter having attached disks and being dedicated to running the file-system server. This static division of responsibilities simplifies system management but does not necessarily lead to the best performance in workloads that need a different balance of computation and I/O. Of course, computational processes sharing a node with a file-system service may receive less CPU time, network bandwidth, and memory bandwidth than they would on a computation-only node. In this paper we begin to examine this issue experimentally. We found that high performance I/O does not necessarily require substantial CPU time, leaving plenty of time for application computation. There were some complex file-system requests, however, which left little CPU time available to the application. (The impact on network and memory bandwidth still needs to be determined.) For applications (or users) that cannot tolerate an occasional interruption, we recommend that they continue to use only compute nodes. For tolerant applications needing more cycles than those provided by the compute nodes, we recommend that they take full advantage of both compute and I/O nodes for computation, and that operating systems should make this possible.

Kotz, David↗

Newton-Krylov-Schwarz: An implicit solver for CFD

Newton-Krylov methods and Krylov-Schwarz (domain decomposition) methods have begun to become established in computational fluid dynamics (CFD) over the past decade. The former employ a Krylov method inside of Newton's method in a Jacobian-free manner, through directional differencing. The latter employ an overlapping Schwarz domain decomposition to derive a preconditioner for the Krylov accelerator that relies primarily on local information, for data-parallel concurrency. They may be composed as Newton-Krylov-Schwarz (NKS) methods, which seem particularly well suited for solving nonlinear elliptic systems in high-latency, distributed-memory environments. We give a brief description of this family of algorithms, with an emphasis on domain decomposition iterative aspects. We then describe numerical simulations with Newton-Krylov-Schwarz methods on aerodynamics applications emphasizing comparisons with a standard defect-correction approach, subdomain preconditioner consistency, subdomain preconditioner quality, and the effect of a coarse grid.

Cai, Xiao-Chuan↗

Local multiplicative Schwarz algorithms for convection-diffusion equations

We develop a new class of overlapping Schwarz type algorithms for solving scalar convection-diffusion equations discretized by finite element or finite difference methods. The preconditioners consist of two components, namely, the usual two-level additive Schwarz preconditioner and the sum of some quadratic terms constructed by using products of ordered neighboring subdomain preconditioners. The ordering of the subdomain preconditioners is determined by considering the direction of the flow. We prove that the algorithms are optimal in the sense that the convergence rates are independent of the mesh size, as well as the number of subdomains. We show by numerical examples that the new algorithms are less sensitive to the direction of the flow than either the classical multiplicative Schwarz algorithms, and converge faster than the additive Schwarz algorithms. Thus, the new algorithms are more suitable for fluid flow applications than the classical additive or multiplicative Schwarz algorithms.

Cai, Xiao-Chuan↗

Parallel Newton-Krylov-Schwarz algorithms for the transonic full potential equation

We study parallel two-level overlapping Schwarz algorithms for solving nonlinear finite element problems, in particular, for the full potential equation of aerodynamics discretized in two dimensions with bilinear elements. The overall algorithm, Newton-Krylov-Schwarz (NKS), employs an inexact finite-difference Newton method and a Krylov space iterative method, with a two-level overlapping Schwarz method as a preconditioner. We demonstrate that NKS, combined with a density upwinding continuation strategy for problems with weak shocks, is robust and, economical for this class of mixed elliptic-hyperbolic nonlinear partial differential equations, with proper specification of several parameters. We study upwinding parameters, inner convergence tolerance, coarse grid density, subdomain overlap, and the level of fill-in in the incomplete factorization, and report their effect on numerical convergence rate, overall execution time, and parallel efficiency on a distributed-memory parallel computer.

Cai, Xiao-Chuan↗

First-Order System Least-Squares for the Navier-Stokes Equations

This paper develops a least-squares approach to the solution of the incompressible Navier-Stokes equations in primitive variables. As with our earlier work on Stokes equations, we recast the Navier-Stokes equations as a first-order system by introducing a velocity flux variable and associated curl and trace equations. We show that the resulting system is well-posed, and that an associated least-squares principle yields optimal discretization error estimates in the H(sup 1) norm in each variable (including the velocity flux) and optimal multigrid convergence estimates for the resulting algebraic system.

Bochev, P.↗

First-Order System Least Squares for Velocity-Vorticity-Pressure Form of the Stokes Equations, with Application to Linear Elasticity

In this paper, we study the least-squares method for the generalized Stokes equations (including linear elasticity) based on the velocity-vorticity-pressure formulation in d = 2 or 3 dimensions. The least squares functional is defined in terms of the sum of the L(exp 2)- and H(exp -1)-norms of the residual equations, which is weighted appropriately by by the Reynolds number. Our approach for establishing ellipticity of the functional does not use ADN theory, but is founded more on basic principles. We also analyze the case where the H(exp -1)-norm in the functional is replaced by a discrete functional to make the computation feasible. We show that the resulting algebraic equations can be uniformly preconditioned by well-known techniques.

Cai, Zhiqiang↗

First-Order System Least Squares for the Stokes Equations, with Application to Linear Elasticity

Following our earlier work on general second-order scalar equations, here we develop a least-squares functional for the two- and three-dimensional Stokes equations, generalized slightly by allowing a pressure term in the continuity equation. By introducing a velocity flux variable and associated curl and trace equations, we are able to establish ellipticity in an H(exp 1) product norm appropriately weighted by the Reynolds number. This immediately yields optimal discretization error estimates for finite element spaces in this norm and optimal algebraic convergence estimates for multiplicative and additive multigrid methods applied to the resulting discrete systems. Both estimates are uniform in the Reynolds number. Moreover, our pressure-perturbed form of the generalized Stokes equations allows us to develop an analogous result for the Dirichlet problem for linear elasticity with estimates that are uniform in the Lame constants.

Cai, Z.↗

Removal of Thin Cirrus Path Radiances in the 0.4-1.0 micron Spectral Region Using the 1.375-micron Strong Water Vapor Absorption Channel

Through analysis of spectral imaging data acquired with the Airborne Visible Infrared Imaging Spectrometer (AVIRIS) from an ER-2 aircraft at 20 km altitude during several field programs, it was found that narrow channels near the center of the strong 1.38-micron water vapor band are very sensitive in detecting thin cirrus clouds. Based on this observation from AVIRIS data, a channel centered at 1.375 microns with a width of 30 nm was selected for the Moderate Resolution Imaging Spectrometer (MODIS) for remote sensing of cirrus clouds from space. The sensitivity of the 1.375-micron MODIS channel to detect thin cirrus clouds during the day time is expected to be one to two orders of magnitude better than the current infrared emission techniques. As a result, a larger fraction of the satellite data will likely be identified as containing cirrus clouds. In order to make better studies of surface reflectance properties, thin cirrus effects must be removed from satellite images. We have developed an empirical approach for removing/correcting thin cirrus effects in the 0.4 - 1.0 micron region using channels near 1.375 microns. This algorithm will be incorporated into the present MODIS atmospheric correction algorithms for ocean color and land applications and will yield improved MODIS atmospheric aerosol, land surface, and ocean color products.

Gao, Bo-Cai↗

A New and Fast Method for Smoothing Spectral Imaging Data

The Airborne Visible Infrared Imaging Spectrometer (AVIRIS) acquires spectral imaging data covering the 0.4 - 2.5 micron wavelength range in 224 10-nm-wide channels from a NASA ER-2 aircraft at 20 km. More than half of the spectral region is affected by atmospheric gaseous absorption. Over the past decade, several techniques have been used to remove atmospheric effects from AVIRIS data for the derivation of surface reflectance spectra. An operational atmosphere removal algorithm (ATREM), which is based on theoretical modeling of atmospheric absorption and scattering effects, has been developed and updated for deriving surface reflectance spectra from AVIRIS data. Due to small errors in assumed wavelengths and errors in line parameters compiled on the HITRAN database, small spikes (particularly near the centers of the 0.94- and 1.14-micron water vapor bands) are present in this spectrum. Similar small spikes are systematically present in entire ATREM output cubes. These spikes have distracted geologists who are interested in studying surface mineral features. A method based on the "global" fitting of spectra with low order polynomials or other functions for removing these weak spikes has recently been developed by Boardman (this volume). In this paper, we describe another technique, which fits spectra "locally" based on cubic spline smoothing, for quick post processing of ATREM apparent reflectance spectra derived from AVIRIS data. Results from our analysis of AVIRIS data acquired over Cuprite mining district in Nevada in June of 1995 are given. Comparisons between our smoothed spectra and those derived with the empirical line method are presented.

Gao, Bo-Cai↗

Asymptotic Solutions for Optical Properties of Large Particles with Strong Absorption

For scattering calculations involving nonspherical particles such as ice crystals, we show that the transverse wave condition is not applicable to the refracted electromagnetic wave in the context of geometric optics when absorption is involved. Either the TM wave condition (i.e., where the magnetic field of the refracted wave is transverse with respect to the wave direction) or the TE wave condition (i.e., where the electric field is transverse with respect to the propagating direction of the wave) may be assumed for the refracted wave in an absorbing medium to locally satisfy the electromagnetic boundary condition in the ray tracing calculation. The wave mode assumed for the refracted wave affects both the reflection and refraction coefficients. As a result, a nonunique solution for these coefficients is derived from the electromagnetic boundary condition. In this study we have identified the appropriate solution for the Fresnel reflection/refraction coefficients in light scattering calculation based on the ray tracing technique. We present the 3 x 2 refraction or transmission matrix that completely accounts for the inhomogeneity of the refracted wave in an absorbing medium. Using the Fresnel coefficients for an absorbing medium, we derive an asymptotic solution in an analytical format for the scattering properties of a general polyhedral particle. Numerical results are presented for hexagonal plates and columns with both preferred and random orientations. The asymptotic theory can produce reasonable accuracy in the phase function calculations in the infrared window region (wavelengths near 10 micron) if the particle size (in diameter) is on the order of 40 micron or larger. However, since strong absorption is assumed in the computation of the single-scattering albedo in the asymptotic theory, the single scattering albedo does not change with variation of the particle size. As a result, the asymptotic theory can lead to substantial errors in the computation of single-scattering albedo for small and moderate particle sizes. However, from comparison of the asymptotic results with the FDTD solution, it is expected that a convergence between the FDTD results and the asymptotic theory results can be reached when the particle size approaches 200 micron. We show that the phase function at side-scattering and backscattering angles is insensitive to particle shape if the random orientation condition is assumed. However, if preferred orientations are assumed for particles, the phase function has a strong dependence on scattering azimuthal angle. The single-scattering albedo also shows very strong dependence on the inclination angle of incident radiation with respect to the rotating axis for the preferred particle orientations.

Yang, Ping↗

Radiative Properties of Cirrus Clouds in the Infrared (8-13 microns) Spectral Region

Atmospheric radiation in the infrared (IR) 8-13 microns spectral region contains a wealth of information that is very useful for the retrieval of ice cloud properties from aircraft or space-borne measurements. To provide the scattering and absorption properties of nonspherical ice crystals that are fundamental to the IR retrieval implementation, we use the finite-difference time domain (FDTD) method to solve for the extinction efficiency, single-scattering albedo, and the asymmetry parameter of the phase function for ice crystals smaller than 40 microns. For particles larger than this size, the improved geometric optics method (IGOM) can be employed to calculate the asymmetry parameter with an acceptable accuracy, provided that we properly account for the inhomogeneity of the refracted wave due to strong absorption inside the ice particle. A combination of the results computed from the two methods provides the asymmetry parameter for the entire practical range of particle sizes between 1 micron and 10000 microns over wavelengths ranging from 8 microns to 13 microns. For the extinction and absorption efficiency calculations, several methods including the IGOM, Mie solution for equivalent spheres (MSFES), and the anomalous diffraction theory (ADT) can lead to a substantial discontinuity in comparison with the FDTD solutions for particle sizes on the order of 40 microns. To overcome this difficulty, we have developed a novel approach called the stretched scattering potential method (SSPM). For the IR 8-13 microns spectral region, we show that SSPM is a more accurate approximation than ADT, MSFES, and IGOM. The SSPM solution can be further refined numerically. Through a combination of the FDTD and SSPM, we have computed the extinction and absorption efficiency for hexagonal ice crystals with sizes ranging from 1 to 10000 microns at 12 wavelengths between 8 and 13 microns Calculations of the cirrus bulk scattering and absorption properties are performed for 30 size distributions obtained from various field campaigns for midlatitude and tropical cirrus cloud systems. Parameterization of these bulk scattering properties is carried out by using second-order polynomial functions for the extinction efficiency and the single-scattering albedo and the power law expression for the asymmetry parameter. We note that the volume-normalized extinction coefficient can be separated into two parts: one is inversely proportional to effective size and is independent of wavelength, and the other is the wavelength-dependent effective extinction efficiency. Unlike conventional parameterization efforts, the present parameterization scheme is more accurate because only the latter part of the volume-normalized extinction coefficient is approximated in terms of an analytical expression. After averaging over size distribution, the single-scattering albedo is shown to decrease with an increase in effective size for wavelengths shorter than 10.0 microns whereas the opposite behavior is observed for longer wavelengths. The variation of the asymmetry parameter as a function of effective size is substantial when the effective size is smaller than 50 microns. For effective sizes larger than 100 microns, the asymmetry parameter approaches its asymptotic value. The results derived in this study can be useful to remote sensing applications involving IR window bands under cirrus cloud conditions.

Yang, Ping↗

The Effect of the Wall Contact and Post-Growth, Cool-Down on Defects in CdTe Crystals Grown By 'Contactless' PVT

To take a maximum advantage of materials processing in microgravity for understanding the effects of gravity, gravity-independent effects should be minimized. In crystal growth, the quality of the grown crystals may depend, among other factors, on their interaction with the walls of the processing container during and after growth, and on the rate of the crystal cool-down at the end of the process. To investigate the above phenomena, a series of CdTe crystal growth processes was carried out. The crystals were grown by physical vapor transport without contact with the side walls of the silica glass ampoules. To eliminate the effect of the seed quality, and to reduce the number of nuclei and related crystal grains, the Low Supersaturation Nucleation technique was applied. The source temperature was 930 C, the undercooling was a few degrees. The crystals, having the diameter of 25 mm, grew at the rate of a few mm per day. The post-growth cool-down to the room temperature was conducted at different rates, and lasted from a few minutes to four days. The crystals were characterized using chemical etching, low temperature luminescence, and Synchrotron White Beam X-ray Topography techniques. The dislocation (etch pit) density was measured and its distribution was analyzed by comparison with Poisson curves and with the Normalized Radial Distribution Correlation Function. In the regions where the crystal is in contact with silica, the materials show a considerable strain field which extends for a few mm or more from the silica-crystal interface. In the reference crystal grown with contact with the ampoule walls, and when the crystals are cooled at the highest rates, the etch pit/dislocation density is in the high 10(exp 5) per square centimeter region. Typical EPD values for lower cool-down rates are in the lower 10(exp 4) per square centimeter region. In some areas the actual dislocation density was about 10(exp 3) per square centimeter or even less. No apparent effect of the cool-down rate on polygonization was observed. Low temperature PL spectra show, that the dominant peak is (D(sup 0), h) and (A(sup 0), e) for samples with low and high dislocation densities, respectively. For low EPD crystals a peak at 1.45 eV with 21 meV phonon replicas was observed and attributed to donor-acceptor pair to neutral copper-acceptor transition. In high EPD crystals this PL structure was not observed, apparently due to the masking effect of the strong contribution from the dislocation band.

Palosz, W.↗

The Effect of the Wall Contact and Post-Growth C001-Down on Defects in CdTe Crystals Grown by Contactless PVT

In crystal growth, the quality of the final material may depend, among other factors, on its interaction with the walls of the ampoule during and after the growth, and on the rate of the crystal cool-down at the end of ate the process. To investigate the above phenomena, a series of CdTe crystal growth processes was carried out, The crystals were grown by physical vapor transport without contact with the side walls of the silica glass ampoules, applying the Low Supersaturation Nucleation technique. The source temperature was 930 C, the undercooling was a few degrees. The crystals, having the diameter of 25 mm, grew at the rate of a few mm per day. The post-growth cool-down to the room temperature was conducted at different rates, and lasted from a few minutes to four days. The crystals were characterized using chemical etching low temperature luminescence, and Synchrotron White Beam X-ray Topography techniques. The dislocation (etch pit) density was measured and its distribution was analyzed by comparison with Poisson curves and with the Normalized Radial Distribution Correlation Function. It was found that the contact of the crystal with silica leads to a strain field and high (in the 105 sq cm range) dislocation (etch pit) density. Similar defect concentrations were found in crystals subjected to fast post-growth cool-down. Typical EPD values for lower cool-down rates and in regions not affected by wall interactions are in the lower 10(exp 4) sq cm range. In some areas the actual dislocation density was about 10(exp 3) sq cm or even less. No apparent effect of the cool-down rate on polygonization was observed. A fine structure could be discerned in low-temperature PL spectra of crystals with low dislocation density.

Palosz, W.↗