Search NASA⌕ Search

SEARCH · Search NASA

Results for “height 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 955 records · Page 53

Interaction with the lower ionosphere of electromagnetic pulses from lightning - Heating, attachment, and ionization

A Boltzmann formulation of the electron distribution function and Maxwell's equations for the EM fields are used to simulate the interaction of lightning radiated EM pulses with the lower ionosphere. Ionization and dissociative attachment induced by the heated electrons cause significant changes in the local electron density, N(e). Due to 'slow' field changes of typical lightning EM pulses over time scales of tens of microsec, the distribution function follows the quasi-equilibrium solution of the Boltzmann equation in the altitude range of interest (70 to 100 km). The EM pulse is simulated as a planar 100 microsec long single period oscillation of a 10 kHz wave injected at 70 km. Under nighttime conditions, individual pulses of intensity 10-20 V/m (normalized to 100 km horizontal distance) produce changes in N(e) of 1-30 percent while a sequence of pulses leads to strong modification of N(e) at altitudes less than 95 km. The N(e) changes produce a 'sharpening' of the lower ionospheric boundary by causing a reduction in electron density at 75-85 km (due to attachment) and a substantial increase at 85-95 km (due to ionization) (e.g., the scale height decreases by a factor of about 2 at about 85 km for a single 20 V/m EM pulse). No substantial N(e) changes occur during daytime.

Taranenko, Y. N.↗

Linear operating region in the ozone dial photon counting system

Ozone is a relatively unstable molecule found in Earth's atmosphere. An ozone molecule is made up of three atoms of oxygen. Depending on where ozone resides, it can protect or harm life on Earth. High in the atmosphere, about 15 miles up, ozone acts as a shield to protect Earth's surface from the sun's harmful ultraviolet radiation. Without this shield, we would be more susceptible to skin cancer, cataracts, and impaired immune systems. Closer to Earth, in the air we breathe, ozone is a harmful pollutant that causes damage to lung tissue and plants. Since the early 1980's, airborne lidar systems have been used for making measurements of ozone. The differential absorption lidar (DIAL) technique is used in the remote measurement of O3. This system allows the O3 to be measured as function of the range in the atmosphere. Two frequency-doubled Nd:YAG lasers are used to pump tunable dye lasers. The lasers are operating at 289 nm for the DIAL on-line wavelength of O3, and the other one is operated at 300 nm for the off-line wavelength. The DIAL wavelengths are produced in sequential laser pulses with a time separation of 300 micro s. The backscattered laser energy is collected by telescopes and measured using photon counting systems. The photon counting system measures the light signal by making use of the photon nature of light. The output pulse from the Photo-Multiplier Tube (PE), caused by a photon striking the PMT photo-cathode, is amplified and passed to a pulse height discriminator. The peak value of the pulse is compared to a reference voltage (discrimination level). If the pulse amplitude exceeds the discrimination level, the discriminator generates a standard pulse which is counted by the digital counter. Non-linearity in the system is caused by the overlapping of pulses and the finite response time of the electronics. At low count rates one expects the system to register one event for each output pulse from the PMT corresponding to a photon incident upon the photocathode, however, at higher rates the limitations of the discrimination/counting system will cause the observed count rate to be non-linear with respect to the true count rate. Depending on the pulse height distribution and the discriminator level, the overlapping of pulses (pulse pile-up) can cause count loss or even an additional apparent count gain as the signal levels increase. Characterization of the system, including the pulse height distribution, the signal to noise ratio, and the effect of the discriminator threshold level, is critical in maximizing the linear operating region of the system, thus greatly increasing the useful dynamic range of the system.

Andrawis, Madeleine↗

Morphology and Density Structure of Post-CME Current Sheets

Eruption of a coronal mass ejection (CME) is believed to drag and open the coronal magnetic field, presumably leading to the formation of a large-scale current sheet and field relaxation by magnetic reconnection. This paper analyzes the physical characteristics of ray-like coronal features formed in the aftermath of CMEs, to confirm whether interpreting such phenomena in terms of a reconnecting current sheet is consistent with observations. Methods: The study focuses on UVCS/SOHO and LASCO/SOHO measurements of the ray width, density excess, and coronal velocity field as a function of the radial distance. The morphology of the rays implies that they are produced by Petschek-like reconnection in the large-scale current sheet formed in the wake of CME. The hypothesis is supported by the flow pattern, often showing outflows along the ray, and sometimes also inflows into the ray. The inferred inflow velocities range from 3 to 30 km/s, and are consistent with the narrow opening-angle of rays, which add up to a few degrees. The density of rays is an order of magnitude higher than in the ambient corona. The model results are consistent with the observations, revealing that the main cause of the density excess in rays is a transport of the dense plasma from lower to higher heights by the reconnection outflow.

Vrsnak, B.↗

Influence of Satellite-Based Heterogeneous Vegetation Momentum Roughness on Mesoscale Model Dynamics During IHOP 2002

The sensitivity of mesoscale weather prediction model to a vegetation roughness initialization is investigated for the south central United States. Three different roughness databases are employed: i) a control or standard lookup table roughness that is a function only of land cover type, ii) a spatially heterogeneous roughness database previously derived using a physically based procedure and MODIS imagery, and iii) a MODIS climatologic roughness database that possesses the same spatial heterogeneity as (i) but with mean land class values from (ii). The model used is the Weather Research and Forecast Model (WRF) coupled to the Community Land Model within the Land Information System (LIS). For each simulation, a statistical comparison is made between modeled results and ground observations from meteorological stations within the Oklahoma mesonet and surrounding region during IHOP20O2. A sensitivity analysis on the impact the MODIS-based roughness fields is also made through a time-series intercomparison of temperature bias, probability of detection (POD), average wind speed, boundary layer height, and turbulent kinetic energy (TKE) the results that, for the current replacement of the standard land-cover type based roughness values with the satellite-derived fields statistically improves model performance for most of the observed variables. Further, the satellite-based roughness enhances the surface wind speed, PBL height and TKE production on the order of 3 to l0 percent, with a lesser effect over grassland and cropland domains, and the greater effect over mixed land cover domains

Jasinski, Michael↗

Numerical methods for reduction of topside ionograms

Several alternative methods for solving the group height equation are presented. Three of these are now in operation at Ames Research Center and use data contained in a single ionogram trace. From the data an electron density profile N(h) is computed. If the ionogram also exhibits other traces, reverse ionogram traces are computed, using the N(h) profile, for comparison with the redundant data. When agreement is poor, the initial data trace is reinterpreted, another N(h) profile computed, and the reverse traces generated once again. This process is repeated until a desired degree of consistency is achieved. To reduce the necessity for human intervention and eliminate decision making required in conjunction with the preceding methods, a method is proposed that accepts as input, all data from a single ionogram. In general, no electron density function will satisfy these data exactly, but a best N(h) profile can be computed. Finally, a method is described that eliminates the need to assume that the ionosphere is spherically stratified. Horizontal gradients in electron density are detected and accounted for by processing several ionograms from the same satellite pass simultaneously. This idea is derived as an extension of one of the basic methods.

Mcculley, L.↗

Enhanced Mixing in a Rectangular Duct

An experimental investigation of the mixing of non-reacting opposed rows of jets injected normal to a confined rectangular crossflow has been conducted. Planar Mie-scattering was used to measure the time-average concentration distribution of the jet fluid in planes perpendicular to the duct axis. Particular emphasis was placed on the study of closely spaced orifice configurations applicable to the mixing zone of an RQL combustor. Baseline studies were performed of mixing under "ideal" conditions, i.e., plenum fed jets injecting into a crossflow uniform in velocity and turbulence intensity. In addition, more practical ("non-ideal") issues encountered during hardware design were also studied. As in other studies, mixing effectiveness, determined using a spatial unmixedness parameter based on the variance of mean jet concentration distributions, was found to be optimum when the spacing-to-duct-height ratio was inversely proportional to the square root of the jet-to-mainstream momentum-flux ratio. This relationship is suitable for design under ideal flow conditions. Inlet flow boundary conditions of the jet and approach flow (mainstream) were found to strongly influence mixing performance, but no attempt was made to determine optimum performance under non-ideal conditions. The tests performed do offer some guidance as to expected mixing behavior for several common variables likely to be imposed by hardware constraints. Additionally, in this study it was found that for rows of orifices with opposite centerlines inline, mixing was similar for blockages up to 89 percent (previous crossflow mixing studies concerned with dilution zone configurations, blockages were typically less than 50 percent). Lower levels of unmixedness were obtained as a function of downstream location when axial injection length was minimized. Mixing may be enhanced if orifice centerlines of opposed rows are staggered, but blockage must be =50 percent in this configuration. Round hole and "square" orifice shapes had similar performance. Other variations in orifice shape did not substantially augment overall mixing performance. Furthermore an isothermal mixing data set was generated and used by CFDRC as input to a NO(x) inference code.

Liscinsky, D. S.↗

SOHO/CDS Measurements of Coronal EUV Polarization above the Limb

Attempts to measure polarization in coronal EUV emission above the solar limb have been made using the SOHO/CDS normal-incidence spectrometer which has a polarization sensitivity of about 50%, a property that causes variations in intensity response as a function of the spacecraft's roll angle for polarized light. Such observations were made on the disk and up to 0.22 solar radii above the solar limb in a number of EUV lines during two special roll-maneuvers of the SOHO spacecraft. Measurements of intensity gradients were made above a modestly active equatorial region in 1997 and above a relatively cool polar region in 2001. Observed emission lines include He I 584A, He II 304A, 0 IV 555+610A, 0 V 630A, Mg IX 368A, Mg X 610+625A, and Si XI 303A, formed at temperatures that evenly cover the range in logT from 4.1 to 6.2. Near the disk, measured intensities of all lines fall off exponentially at different rates that can be used to determine the density scale-heights of the emitting plasma, since this emission is dominated by collisional excitation with an Ne-squared dependence. Assuming hydrostatic equilibrium, the intensity gradient for each line can then be converted into a 'scale-height temperature', which is found to be closely related to the ionization temperature of each line over the wide range of lines and solar conditions observed. Beyond a certain distance, intensity gradients of the cooler lines switch over to a flatter exponential slope, suggesting that this radiation is dominated by resonance scattering which varies as Ne to the first power. Such radiation should also be linearly polarized in the plane containing the line-of-sight and the solar center, a signature that would strongly confirm this interpretation.

Thomas, Roger J.↗

Bayesian Estimation of Precipitation from Satellite Passive Microwave Observations Using Combined Radar-Radiometer Retrievals

Precipitation estimation from satellite passive microwave radiometer observations is a problem that does not have a unique solution that is insensitive to errors in the input data. Traditionally, to make this problem well posed, a priori information derived from physical models or independent, high-quality observations is incorporated into the solution. In the present study, a database of precipitation profiles and associated brightness temperatures is constructed to serve as a priori information in a passive microwave radiometer algorithm. The precipitation profiles are derived from a Tropical Rainfall Measuring Mission (TRMM) combined radar radiometer algorithm, and the brightness temperatures are TRMM Microwave Imager (TMI) observed. Because the observed brightness temperatures are consistent with those derived from a radiative transfer model embedded in the combined algorithm, the precipitation brightness temperature database is considered to be physically consistent. The database examined here is derived from the analysis of a month-long record of TRMM data that yields more than a million profiles of precipitation and associated brightness temperatures. These profiles are clustered into a tractable number of classes based on the local sea surface temperature, a radiometer-based estimate of the echo-top height (the height beyond which the reflectivity drops below 17 dBZ), and brightness temperature principal components. For each class, the mean precipitation profile, brightness temperature principal components, and probability of occurrence are determined. The precipitation brightness temperature database supports a radiometer-only algorithm that incorporates a Bayesian estimation methodology. In the Bayesian framework, precipitation estimates are weighted averages of the mean precipitation values corresponding to the classes in the database, with the weights being determined according to the similarity between the observed brightness temperature principal components and the brightness temperature principal components of the classes. Because the classes are stratified by the sea surface temperature and the echo-top-height estimator, the number of classes that are considered for retrieval is significantly smaller than the total number of classes, making the algorithm computationally efficient. The radiometer-only algorithm is applied to TMI observations, and precipitation estimates are compared with combined TRMM precipitation radar (PR) TMI reference estimates. The TMI-only algorithm, supported by the empirically derived database, produces estimates that are more consistent with the reference values than the precipitation estimates from the version-6 TRMM facility TMI algorithm. Cloud-resolving model simulations are used to assign a latent heating profile to each precipitation profile in the empirically derived database, making it possible to estimate latent heating using the radiometer-only algorithm. Although the evaluation of latent heating estimates in this study is preliminary, because realistic conditional probability distribution functions are attached to latent heating structures in the algorithm s database, a generally positive impact on latent heating estimation from passive microwave observations is expected.

Grecu, Mircea↗

Dependence of Lunar Surface Charging on Solar Wind Plasma Conditions and Solar Irradiation

The surface of the Moon is electrically charged by exposure to solar radiation on its dayside, as well as by the continuous flux of charged particles from the various plasma environments that surround it. An electric potential develops between the lunar surface and ambient plasma, which manifests itself in a near-surface plasma sheath with a scale height of order the Debye length. This study investigates surface charging on the lunar dayside and near-terminator regions in the solar wind, for which the dominant current sources are usually from the pohotoemission of electrons, J(sub p), and the collection of plasma electrons J(sub e) and ions J(sub i). These currents are dependent on the following six parameters: plasma concentration n(sub 0), electron temperature T(sub e), ion temperature T(sub i), bulk flow velocity V, photoemission current at normal incidence J(sub P0), and photo electron temperature T(sub p). Using a numerical model, derived from a set of eleven basic assumptions, the influence of these six parameters on surface charging - characterized by the equilibrium surface potential, Debye length, and surface electric field - is investigated as a function of solar zenith angle. Overall, T(sub e) is the most important parameter, especially near the terminator, while J(sub P0) and T(sub p) dominate over most of the dayside.

Photoemission↗

Causes of Long-Term Drought in the United States Great Plains

The United States Great Plains (USGP) experienced a number of multi-year droughts during the last century, most notably the droughts of the 1930s and 1950s. This study examines the causes of such droughts using ensembles of long term (1930-1999) simulations carried out with the NASA Seasonal-to-Interannual Prediction Project (NSIPP-1) atmospheric general circulation model (AGCM) forced with observed sea surface temperatures (SSTs). The results show that the model produces long-term (multi-year) variations in the USGP precipitation that are similar to those observed. A correlative analysis suggests that the ensemble mean low frequency (time scales longer than about 6 years) rainfall variations in the USGP are linked to a pan-Pacific pattern of SST variability that is the leading empirical orthogonal function (EOF) in the low frequency SST data. The link between the SST and the Great Plains precipitation is confirmed in idealized AGCM simulations, in which the model is forced by the 2 polarities of the pan-Pacific SST pattern. The idealized simulations further show that it is primarily the tropical part of the SST anomalies that influence the USGP. As such, the USGP tend to have above normal precipitation when the tropical Pacific SSTs are above normal, while there is a tendency for drought when the tropical SSTs are cold. The upper tropospheric response to the pan-Pacific SST EOF shows a global-scale pattern with a strong wave response in the Pacific and a substantial zonally-symmetric component in which USGP pluvial (drought) conditions are associated with reduced (enhanced) heights throughout the extra-tropics. The potential predictability of rainfall in the USGP associated with SSTs is rather modest, with on average about 1/3 of the total low frequency rainfall variance forced by SST anomalies. Further idealized experiments with climatological SST, suggest that the remaining low frequency variance in the USGP precipitation is the result of interactions with soil moisture. In particular, simulations with soil moisture feedback show a six-fold increase in the variance in annual USGP precipitation compared with simulations in which the soil feedback is excluded. In addition to increasing variance, the interactions with the soil introduce year-to-year memory in the hydrological cycle that is consistent with a red noise process, in which the low frequencies in the deep soil are the result of integrating a net forcing (precipitation-evaporation-runoff) that is white noise on interannual time scales. As such, the role of low frequency SST variability is to introduce a bias to the net forcing on the soil moisture that drives the random process preferentially to either wet or dry conditions.

Schubert, Siegfried D.↗

Angular clustering properties of gamma-ray bursts and quantitative constraints on their distances

Upper limits on the angular correlation function, w(theta), of gamma-ray bursts are obtained which constrain the sources to be very nearby or very distant. The angular clustering properties of the sources in two recent gamma-ray burst catalogs (Golenetskii et al., 1986 and Atteia et al., 1987) have been derived by determining their two-point w(theta) values, and the correlations of both data sets are consistent with w(theta) = 0. It is shown that if the spatial correlation of gamma-ray bursts resembles that of galaxies or galactic clusters, present gamma-ray burst catalogs must be complete to distances in excess of 100 Mpc, while if gamma-ray bursts are associated with a Galactic disk population, then the upper limits on w(theta) constrain the sampling depth to less than two disk scale heights.

Hartmann, Dieter↗

Wake Vortex Transport and Decay in Ground Effect: Vortex Linking with the Ground

Numerical simulations are carried out with a three-dimensional Large-Eddy Simulation (LES) model to explore the sensitivity of vortex decay and transport in ground effect (IGE). The vortex decay rates are found to be strongly enhanced following maximum descent into ground effect. The nondimensional decay rate is found to be insensitive to the initial values of circulation, height, and vortex separation. The information gained from these simulations is used to construct a simple decay relationship. This relationship compares well with observed data from an IGE case study. Similarly, a relationship for lateral drift due to ground effect is constructed from the LES data. In the second part of this paper, vortex linking with the ground is investigated. Our numerical simulations of wake vortices for IGE show that a vortex may link with its image beneath the ground, if the intensity of the ambient turbulence is moderate to high. This linking with the ground (which is observed in real cases)gives the appearance of a vortex tube that bends to become vertically oriented and which terminates at the ground. From the simulations conducted, the linking time for vortices in the free atmosphere; i.e., a function of ambient turbulence intensity.

Proctor, Fred H.↗

Sensitivity of Cirrus Bidirectional Reflectance at MODIS Bands to Vertical Inhomogeneity of Ice Crystal Habits and Size Distribution

A common assumption in satellite imager-based cirrus retrieval algorithms is that the radiative properties of a cirrus cloud may be represented by those associated with a specific ice crystal shape (or habit) and a single particle size distribution. However, observations of cirrus clouds have shown that the shapes and sizes of ice crystals may vary substantially with height within the clouds. In this study we investigate the sensitivity of the top-of-atmosphere bidirectional reflectances at two MODIS bands centered at 0.65 micron and 2.11 micron to the cirrus models assumed to be either a single homogeneous layer or three distinct but contiguous, layers. First, we define the single- and three-layer cirrus cloud models with respect to ice crystal habit and size distribution on the basis of in situ replicator data acquired during the First ISCCP Regional Experiment (FIRE-II), held in Kansas during the fall of 1991. Subsequently, fundamental light scattering and radiative transfer theory is employed to determine the single scattering and the bulk radiative properties of the cirrus cloud. Regarding the radiative transfer computations, we present a discrete form of the adding/doubling principle by introducing a direct transmission function, which is computationally straightforward and efficient an improvement over previous methods. For the 0.65 micron band, at which absorption by ice is negligible, there is little difference between the bidirectional reflectances calculated for the one- and three-layer cirrus models, suggesting that the vertical inhomogeneity effect is relatively unimportant. At the 2.11 micron band, the bidirectional reflectances computed for both optically thin (tau = 1) and thick (tau = 10) cirrus clouds show significant differences between the results for the one- and three-layer models. The reflectances computed for the three-layer cirrus model are substantially larger than those computed for the single-layer cirrus. Finally, we find that cloud reflectance is very sensitive to the optical properties of the small crystals that predominate in the top layer of the three-layer cirrus model. It is critical to define the most realistic geometric shape for the small "quasi-spherical" ice crystals in the top layer for obtaining reliable single-scattering parameters and bulk radiative properties of cirrus.

Yang, P.↗

Potential Scientific Output from the Visible Channels of EPIC Spectroradiometer as Part of the DSCOVR L1 Mission

In addition to 4 UV channels, the EPIC spectroradiometer will provide measurements at 6 visible (and near IR) channels (443,551,680,687.75, 764 and 779.9 nm) at roughly 10 km spatial resolution. The scattering angles are near backscattering and vary between 165 and 176 degree. Two pairs {680 and 687.75 nm} and {764 and 779.9 nm} represent 02 B- and A-bands (and their references) channels; they will be used for cloud height measurements over land and ocean. The B-band channel will contribute to the more absorbing A-band measurements over bright vegetation. A pair {680 and 779.9 nm} will be used for retrieving vegetation properties. Due to the special EPIC geometry, the illuminated part of the leaves will be always observed. As a result, in additional to the traditional Leaf Area Index (LAI), these observations for the first time will provide the sunlit fraction of LAI. Since the sunlit and shaded leaves exhibit different photosynthetic response to incident radiation, these measurements will help to improve global ecological and biogeochemistry models. Finally, a pair {443 and 551 nm} will be used for atmospheric correction. As a by-product of the atmospheric correction algorithm, we also expect to get aerosol optical thickness (AOT) and surface bidirectional reflection function (BRF). The presentation will briefly overview the proposed science algorithms.

Marshak, Alexander↗

A study of marine stratocumulus using lidar and other FIRE aircraft observations

The National Center for Atmospheric Research (NCAR) airborne infrared lidar system (NAILS) used in the 1987 First ISCCP Regional Experiment (FIRE) off the coast of California is a 10.6 microns wavelength carbon dioxide lidar system constructed by Ron Schwiesow and co-workers at NCAR. The lidar is particularly well suited for detailed observations of cloud shapes; i.e., height of cloud top (when flying above cloud and looking down) and cloud base (when flying below cloud and looking up) along the flight path. A brief summary of the lidar design characteristics is given. The lidar height resolution of plus or minus 3 m allows for the distance between the aircraft and cloud edge to be determined with this accuracy; however, the duration of the emitted pulse is approximately 3 microseconds, which corresponds to a 500 m pulse length. Therefore, variations in backscatter intensities within the clouds can normally not be resolved. Hence the main parameter obtainable from the lidar is distance to cloud; in some cases the cloud depth can also be determined. During FIRE the lidar was operational on 7 of the 10 Electra flights, and data were taken when the distance between cloud and aircraft (minimum range) was at least 500 m. The lidar was usually operated at 8 Hz, which at a flight speed of 100 m s(-1) translates into a horizontal resolution of about 12 m. The backscatter as function of time (equivalent to distance) for each laser pulse is stored in digital form on magnetic tape. Currently, three independent variables are available to the investigators on the FIRE Electra data tapes: lidar range to cloud, strength of return (relative power), and pulse width of return, which is related to penetration depth.

Jensen, Jorgen B.↗

Moment method analysis of linearly tapered slot antennas

A method of moments (MOM) model for the analysis of the Linearly Tapered Slot Antenna (LTSA) is developed and implemented. The model employs an unequal size rectangular sectioning for conducting parts of the antenna. Piecewise sinusoidal basis functions are used for the expansion of conductor current. The effect of the dielectric is incorporated in the model by using equivalent volume polarization current density and solving the equivalent problem in free-space. The feed section of the antenna including the microstripline is handled rigorously in the MOM model by including slotline short circuit and microstripline currents among the unknowns. Comparison with measurements is made to demonstrate the validity of the model for both the air case and the dielectric case. Validity of the model is also verified by extending the model to handle the analysis of the skew-plate antenna and comparing the results to those of a skew-segmentation modeling results of the same structure and to available data in the literature. Variation of the radiation pattern for the air LTSA with length, height, and taper angle is investigated, and the results are tabulated. Numerical results for the effect of the dielectric thickness and permittivity are presented.

Koeksal, Adnan↗

An Einstein survey of the 1 keV soft X-ray background in the Galactic plane

We have analyzed 56 Einstein Observatory Imaging Proportional Counter (IPC) observations within +/- 3 deg of the Galactic plane in order to determine the low-latitude soft X-ray background flux in the 0.56-1.73 keV band. Any detected X-ray point source which fell within our regions of study was removed from the image, enabling us to present maps of the background flux as a function of Galactic latitude along 18 meridians. These maps reveal considerable structure to the background in the Galactic plane on an angular scale of approximately 1 deg. Our results are compared with those of an earlier study of the 1 keV X-ray background along l = 25 deg by Kahn & Caillault. The double-peaked structure they found is not discernible in our results, possibly because of the presence of solar backscattered flux in their data. A model which takes into account contributions to the background by extragalactic and stellar sources, the distribution of both atomic and molecular absorbing material with the Galaxy, the energy dependence of the cross section for absorption of X-rays, and the energy dependence of the detector has been constructed and fitted to these new data to derive constraints on the scale height, temperature, and volume emissivity of the unaccounted-for X-ray-emitting material. The results of this model along l = 25 deg are roughly similar to those of the model of Kahn & Caillault along the same meridian.

Stanford, John M.↗

Shaping Diffraction-Grating Grooves to Optimize Efficiency

A method of shaping diffraction-grating grooves to optimize the spectral efficiency, spectral range, and image quality of a spectral imaging instrument is under development. The method is based on the use of an advanced design algorithm to determine the possibly complex shape of grooves needed to obtain a desired efficiency-versus-wavelength response (see figure). Then electron- beam fabrication techniques are used to realize the required groove shape. The method could be used, for example, to make the spectral efficiency of the grating in a given wavelength range proportional to the inverse of the spectral efficiency of a photodetector array so that the overall spectral efficiency of the combination of the grating and the photodetector array would be flat. The method has thus far been applied to one-dimensional gratings only, but in principle, it is also applicable to two-dimensional gratings. The algorithm involves calculations in the spatial-frequency domain. The spatial-frequency spectrum of a grating is represented as a diffraction-order spectral-peak-width function multiplied by an efficiency function for a single grating groove. This representation affords computational efficiency and accuracy by making it possible to consider only the response from one grating groove (one period of the grating), instead of from the whole grating area, in determining the response from the entire grating. This combination of efficiency and accuracy is crucial for future extensions of the algorithm to two-dimensional designs and to designs in which polarization must also be taken into account. The algorithm begins with the definition of target values of relative efficiency that represent the desired spectral response of the grating in certain spectral frequencies calculated from the diffraction order and wavelength. The grating period is divided into a number of cells - typically, 100. The phase contribution from each cell is determined from the phase of the incident electromagnetic wave and the height of the grating surface in the cell. The total contribution from all cells to each target value is then calculated. Then a method known to specialists as the optimum-rotation-angle method is used to adjust the height of each cell so that the total response from all cells is optimized. The computation is iterative and continues until the desired response is obtained. In the event that the desired response is unphysical, the algorithm nevertheless strives to generate a grating-grove profile for which the response approximates the desired one as closely as possible.

Backlund, John↗