Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffractive computing”

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.

340 records · Page 19

A study of fracture mechanisms in ATD roller bearing

The purpose was to investigate how microstructures, especially anisotropy, affects internal stresses and the overall mechanical response of bearings. Samples with the stress axis along the aligned carbide direction possessed high modulus values compared to those with their axis perpendicular to the carbide stringers. The difference in the modulus was found to be more than a factor of two. A series of experiments was conducted on rolled samples to further investigate this effect; the two sets of results were consistent with each other. The degree of anisotropy of the microstructure in terms of the carbide and matrix orientations was determined using x-ray diffraction. The stress state determination was conducted using neutron diffraction. It was found that there was little variation in the distribution of the internal stresses amongst different samples, indicating that small changes in the processing and geometrical conditions did not result in significant variations in the internal stress. A nominal tensile hoop stress of 39 ksi was obtained for the inner raceway. Furthermore, during the course of the investigation, it was apparent that there was a need to determine the stress state induced by the shrink fitting process. Therefore, a series of experiments was conducted using strain gages to identify the stress distribution in a shrink fitting process in three different types of geometries. Correlations were obtained to estimate the highest stress values in the outer and inner groove geometry. A finite element program based on the ANSYS system was developed to compute the stress distribution in the inner raceway geometry. This analysis indicates that the highest tensile stress in the system occurs at the ID of the ring with a stress value of over 5 times that of the applied radial stress. Results from all these facets were correlated with one another. It appears that the material does not fail as a result of any one single factor, but results from a combination of the various factors investigated.

Zee, Ralph H.↗

Capture of Hydrogen Using ZrNi

Water, as ice, is thought to reside in craters at the lunar poles along with CH4 and H2 . A proposed robotic mission for 2012 will utilize metal/metal hydrides for H2 recovery. Specifications are 99% capture of H2 initially at 5 bar and 100C (or greater), and degassing completely at 300C. Of 47-systems examined using the van't Hoff equation, 4 systems, Mg/MgH2, Mg2Ni/Mg2NiH4, ZrNi/ZrNiH2.8, and Pd/PdH0.77, were considered likely candidates for further examination. It is essential, when selecting a system, to also examine questions regarding activation, kinetics, cyclic stability, and gas impurity effects. After considering those issues, ZrN1 was selected as the most promising candidate, as it is easily activated and rapidly forms ZrNiH 2.8 . In addition, it resists oxide poisoning by CO2, and H2O, while some oxidation by O2 is recommended for improved activation . The presence of hydrogen in the as received Zr-Ni alloy from Alfa Aesar posed additional technical problems. X-ray diffraction of the Zr-Ni powder (-325 mesh), with a Zr:Ni wt% ratio of 70:30, was found to consist of ZrH2, ZrNiH2.8, and ZrNi. ZrH2 in the alloy presented the risk that after degassing that both Zr and ZrNi would be present, and thus lead to erroneous results regarding the reactivity of ZrNi with H2 . Fortunately, ZrH2 is a highly stable hydride that does not degas H2 to any significant extent at temperatures below 300C. Based on equilibrium calculations for the decomposition of ZrH2, only 1 millionth of the hydride decomposed at 300C under a N2 atmosphere flowing at 25 ccm for 64 hours, the longest time for pretreatment employed in the investigation. It was possible, from the X-ray results and knowledge of the Zr:Ni ratio, to compute the composition of a pretreated specimen as being 76 wt% ZrNi and the balance ZrH2.

Patton, Lisa↗

Machine Learning Approaches for Rare-Earth Silicate Environmental Barrier Coating Thermochemical and Thermomechanical Property Predictions

Environmental barrier coatings (EBCs) are a necessary enabling technology for the transition from superalloys to silicon carbide (SiC) ceramic matrix composites (CMCs) in gas turbine engines for increased efficiency and decreased fuel costs. SiC-based CMCs are prone to oxidation-based degradation in the engine hot section, and rare-earth (RE) silicates are promising candidates for EBCs due to their close thermal expansion match to the composite substrate and oxidation resistance. However, the design of EBCs is hindered by the large chemical space of candidate materials and the difficulty in obtaining material properties for engineering optimization. This is especially difficult as research continues into mixed-cation or “high-entropy” RE silicates. First-principles computational methods such as density functional theory (DFT) are highly effective at calculating material properties to guide coating design but are limited by their computational cost. Atomistic simulations have the potential to both accelerate property calculations and expand the properties able to be calculated due to their lower computational compared to DFT. However, they require interatomic potentials (IAPs) specific to the material system of interest, and, to our knowledge, there are no suitable IAPs for RE silicates. Machine learning (ML) is a promising technique to accelerate material property predictions indirectly by generating IAPs for atomistic simulations or via direct prediction. In this work, we present two ML approaches to accelerate the calculation of RE silicate properties relevant to EBC design: 1) a ML-derived interatomic potential (IAP) for atomistic simulations of yttrium disilicate (Y2Si2O7) from DFT training data, and 2) a neural network (NN) model to directly predict thermochemical properties of RE silicates and oxides directly from easily obtainable unit cell parameters. Classical MD simulations using the IAP yield lattice properties and bond lengths in good agreement with both DFT and experimental results from x-ray diffraction. Thermodynamic properties calculated using the finite-displacement phonon method and quasi-harmonic approximation were orders of magnitude faster than DFT with good agreement to the DFT results. The IAP was also used to calculate properties such as coefficient of thermal expansion (CTE) that require large simulation supercells and are therefore difficult with DFT. The IAP correctly predicted the anisotropic nature of the CTE in three different phases of Y2Si2O7. The NN model predicts constant pressure heat capacity, Cp, orders of magnitude faster than DFT calculations, which can enable its use as a surrogate model for multiscale simulations. The two methods presented in this work demonstrate the utility of ML for accelerating the prediction of RE silicate properties, which can in turn accelerate EBC design and optimization.

machine learning↗

Blocking Losses With a Photon Counter

It was not known how to assess accurately losses in a communications link due to photodetector blocking, a phenomenon wherein a detector is rendered inactive for a short time after the detection of a photon. When used to detect a communications signal, blocking leads to losses relative to an ideal detector, which may be measured as a reduction in the communications rate for a given received signal power, or an increase in the signal power required to support the same communications rate. This work involved characterizing blocking losses for single detectors and arrays of detectors. Blocking may be mitigated by spreading the signal intensity over an array of detectors, reducing the count rate on any one detector. A simple approximation was made to the blocking loss as a function of the probability that a detector is unblocked at a given time, essentially treating the blocking probability as a scaling of the detection efficiency. An exact statistical characterization was derived for a single detector, and an approximation for multiple detectors. This allowed derivation of several accurate approximations to the loss. Methods were also derived to account for a rise time in recovery, and non-uniform illumination due to diffraction and atmospheric distortion of the phase front. It was assumed that the communications signal is intensity modulated and received by an array of photon-counting photodetectors. For the purpose of this analysis, it was assumed that the detectors are ideal, in that they produce a signal that allows one to reproduce the arrival times of electrons, produced either as photoelectrons or from dark noise, exactly. For single detectors, the performance of the maximum-likelihood (ML) receiver in blocking is illustrated, as well as a maximum-count (MC) receiver, that, when receiving a pulse-position-modulated (PPM) signal, selects the symbol corresponding to the slot with the largest electron count. Whereas the MC receiver saturates at high count rates, the ML receiver may not. The loss in capacity, symbol-error-rate (SER), and count-rate were numerically computed. It was shown that the capacity and symbol-error-rate losses track, whereas the count-rate loss does not generally reflect the SER or capacity loss, as the slot-statistics at the detector output are no longer Poisson. It is also shown that the MC receiver loss may be accurately predicted for dead times on the order of a slot.

Moision, Burce E.↗

Determining Sizes of Particles in a Flow from DPIV Data

A proposed method of measuring the size of particles entrained in a flow of a liquid or gas would involve utilization of data from digital particle-image velocimetry (DPIV) of the flow. That is to say, with proper design and operation of a DPIV system, the DPIV data could be processed according to the proposed method to obtain particle sizes in addition to particle velocities. As an additional benefit, one could then compute the mass flux of the entrained particles from the particle sizes and velocities. As in DPIV as practiced heretofore, a pulsed laser beam would be formed into a thin sheet to illuminate a plane of interest in a flow field and the illuminated plane would be observed by means of a charge-coupled device (CCD) camera aimed along a line perpendicular to the illuminated plane. Unlike in DPIV as practiced heretofore, care would be taken to polarize the laser beam so that its electric field would lie in the illuminated plane, for the reason explained in the next paragraph. The proposed method applies, more specifically, to transparent or semitransparent spherical particles that have an index of refraction different from that of the fluid in which they are entrained. The method is based on the established Mie theory, which describes the scattering of light by diffraction, refraction, and specular reflection of light by such particles. In the case of a particle illuminated by polarized light and observed in the arrangement described in the preceding paragraph, the Mie theory shows that the image of the particle on the focal plane of the CCD camera includes two glare spots: one attributable to light reflected toward the camera and one attributable to light refracted toward the camera. The distance between the glare spots is a known function of the size of the particle, the indices of refraction of the particle material, and design parameters of the camera optics. Hence, the size of a particle can be determined from the distance between the glare spots. The proposed method would be implemented in an algorithm that would automatically identify, and measure the distance between, the glare spots for each particle for which a suitable image has been captured in a DPIV image frame. The algorithm (see figure) would begin with thresholding of data from the entire image frame to reduce noise, thereby facilitating discrimination of particle images from the background and aiding in the separation of overlapping particles. It is important not to pick a threshold level so high that the light intensity between a given pair of glare spots does not fall below the threshold value, leaving the glare spots disconnected. The image would then be scanned in a sequence of rows and columns of pixels to identify groups of adjacent pixels that contain nonzero brightnesses and that are surrounded by pixels of zero brightness. Each such group would be assumed to constitute the image of one particle. Each such group would be further analyzed to determine whether the image was saturated; saturated particle images must be rejected because the locations of glare spots in saturated images cannot accurately be determined. Within each unsaturated particle image, the centroids (deemed to be the locations) of the glare spots would be determined by means of gradients of brightness distributions and three-point horizontal and three-point vertical Gaussian estimates based on the brightness values of the brightest pixels and the pixels adjacent to them. If the brightness of a given particle image contained only one peak, then it would be assumed that a second glare spot did not exist and that image would be rejected.

Wernet, M. P.↗

Interference from the Deep Space Network's 70-m High Power Transmitter in Goldstone, CA to 3G Mobile Users Operating in the Surrounding Area

The International Telecommunications Union (ITU) has allocated 2110-2200 MHz for the third generation (3G) mobile services. Part of the spectrum (2110-2120 MHz) is allocated for space research service and has been used by the DSN for years for sending command uplinks to deep space missions. Due to the extremely high power transmitted, potential interference to 3G users in areas surrounding DSN Goldstone exists. To address this issue, a preliminary analytical study has been performed and computer models have been developed. The goal is to provide theoretical foundation and tools to estimate the strength of interference as a function of distance from the transmitter for various interference mechanisms, (or propagation modes), and then determine the size of the area in which 3G users are susceptible to interference from the 400-kW transmitter in Goldstone. The focus is non-line-of-sight interference, taking into account of terrain shielding, anomalous propagation mechanisms, and technical and operational characteristics of the DSN and the 3G services.

propagation↗

Efficiency measurements of reflection gratings in the 100-300-A band

Zero-order and first-order efficiencies of plane diffraction grazings over the 100-300-A band have been measured using the NRL beamline facility at the National Synchrotron Light Source. Measurements were taken at grating angles of 5, 10, and 15 deg, and the incident radiation was polarized primarily in the transverse magnetic mode. Four (three laminar, one blazed) of the gratings were ion etched in high quality quartz blanks, overcoated with 200 A of gold. A fifth (blazed) grating was ruled in gold. The line densities are 1000, 2000, and 3600 grooves/mm for those with laminar profiles. Both blazed gratings have a blaze angle of 8 deg and line densities of 3600 grooves/mm. Measured first-order efficiencies up to 15 percent were obtained for the laminar gratings and up to 8 percent was obtained for the ion-etched blazed grating. Much lower efficiencies (not greater than 2 percent) were measured for the ruled grating. A computer program was written to calculate grating efficiencies using a full electromagnetic model. Measured efficiencies of the ion-etched gratings agree well with predicted values, but the measured first-order efficiency of the ruled grating is much lower than predicted.

Meekins, John F.↗

Radio and optical interferometric imaging

Since diffraction-limited imaging with a single aperture yields angular resolution approx. lambda/D, the attainment of high angular resolution with single apertures requires the construction of correspondingly large monolithic apertures, the whole surface of which must be figured to much less than a wavelength. At the longer wavelengths, it is impossible to build a sufficiently large single aperture: for example, at lambda 21 cm, arcsec resolution requires an aperture of diameter approx. 50 km. At the shorter wavelengths, the atmosphere imposes a natural limit in resolution of about one arcsec. However, another route is possible; that is, using synthetic apertures to image the sky. Synthetic apertures are now in use in many fields, e.g., radio interferometry, radar imaging, and magnetic-resonance imaging. Radio-interferometric techniques developed in radio astronomy over the past 40 years are now being applied to optical and IR astronomical imaging by a number of groups. Furthermore, the problem of figuring synthetic apertures is considerably simpler, and can be implemented in a computer: new 'self-calibration' techniques allow imaging even in the presence of phase errors due to the atmosphere.

Cornwell, Tim J.↗

Building and Breaking Carbon Composites with REACTER

Carbon composites have become indispensable for aerospace and other high-performance applications, and a detailed picture of their morphology and failure mechanisms remains difficult to obtain through experiment. REACTER is a versatile computational modeling tool for atomistic molecular dynamics designed to model chemical reactions at the speed and length scales of classical force fields.1 In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites, carbon-fiber reinforced polymers (CFRP) and carbon nanotube (CNT) composites. Carbon fiber core morphologies were created by the method of Desai et al.,2 but using the advanced reaction constraints framework of REACTER, their proposed multistep procedure was reduced to a single uninterrupted molecular dynamics simulation. The carbon fiber filler was embedded into a polymer matrix by simulated in situ polymerization of several thermosetting resins, including bismaleimide and polyarylacetylene, to obtain the final CFRP model. To generate the second class of carbon composite, CNT networks were grown dynamically using the new ‘create atoms’ feature of REACTER, and similarly infiltrated with resin to obtain CNT composites. The resulting models were compared directly to experiment using simulated high resolution transmission electron microscopy and x-ray diffraction. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third-order DFT-based tight-binding (DFTB3) simulations, via a reaction constraint on the total potential energy of the involved atoms.

Molecular Dynamics↗

Evaluation of control laws and actuator locations for control systems applicable to deformable astronomical telescope mirrors

Some of the major difficulties associated with large orbiting astronomical telescopes are the cost of manufacturing the primary mirror to precise tolerances and the maintaining of diffraction-limited tolerances while in orbit. One successfully demonstrated approach for minimizing these problem areas is the technique of actively deforming the primary mirror by applying discrete forces to the rear of the mirror. A modal control technique, as applied to active optics, has previously been developed and analyzed. The modal control technique represents the plant to be controlled in terms of its eigenvalues and eigenfunctions which are estimated via numerical approximation techniques. The report includes an extension of previous work using the modal control technique and also describes an optimal feedback controller. The equations for both control laws are developed in state-space differential form and include such considerations as stability, controllability, and observability. These equations are general and allow the incorporation of various mode-analyzer designs; two design approaches are presented. The report also includes a technique for placing actuator and sensor locations at points on the mirror based upon the flexibility matrix of the uncontrolled or unobserved modes of the structure. The locations selected by this technique are used in the computer runs which are described. The results are based upon three different initial error distributions, two mode-analyzer designs, and both the modal and optimal control laws.

Ostroff, A. J.↗

Resonant Transmissive Modulator Construction for Use in Beam Steering Array

An agile electronically configurable optical beam steering device for use in directional free-space transmission of optical signals is described. The device design employs an array of tunable resonant transmissive modulators constructed from multi-layered stacks of electro-optically active dielectric materials. Each modulator may be individually configured to transmit an optical signal with a specific amount of phase and group velocity modulation. The resulting diffractive interactions between modulator outputs provide a method for directional optical signal transmission. Of particular focus within this presentation is the optimized design of the individual modulator. It is widely known that periodic multi-layered stacks of dielectric materials exhibit reflective and transmissive resonances as described by the stack's transmission function. If an electric field is applied normal to the stack's layers, phenomena such as the electro-optic effect and/or the quantum-confined Stark effect may be employed to change the effective refractive index within the stack. This refractive index change alters the stack's transmission function and modulates the phase and group velocity of the transmitted optical signal. By varying the parameters that describe the stack's specific construction and operation, such as number of layers within the stack, refractive indexes of stack materials, layer thicknesses, and combinations of periodic verses non-periodic layer repetitions, custom transmittance functions may be generated. The transmitted optical signal carrier frequency and bandwidth must be carefully selected to maximize transmission and to minimize absorption. A computational optimization of the variables describing the stack's construction strives to maximize the amount of optical signal modulation verses applied voltage. Trade-offs between methods of increasing device performance verses the limitations of fabrication technologies are considered.

Keys, Andrew S.↗

Delta L: An Apparatus for Measuring Macromolecule Crystal Growth Rates in Microgravity

Strongly diffracting high quality macromolecule crystals of suitable volume are keenly sought for X-ray diffraction analysis so that high-resolution molecular structure data can be obtained. Such data is of tremendous value to medical research, agriculture and commercial biotechnology. In previous studies by many investigators microgravity has been reported in some instances to improve biological macromolecule X-ray crystal quality while little or no improvement was observed in other cases. A better understanding of processes effecting crystal quality improvement in microgravity will therefore be of great benefit in optimizing crystallization success in microgravity. In ground based research with the protein lysozyme we have previously shown that a population of crystals grown under the same solution conditions, exhibit a variation in X-ray diffraction properties (Judge et al., 1999). We have also observed that under the same solution conditions, individual crystals will grow at slightly different growth rates. This phenomenon is called growth rate dispersion. For small molecule materials growth rate dispersion has been directly related to crystal quality (Cunningham et al., 1991; Ristic et al., 1991). We therefore postulate that microgravity may act to improve crystal quality by reducing growth rate dispersion. If this is the case then as different, Materials exhibit different degrees of growth rate dispersion on the ground then growth rate dispersion could be used to screen which materials may benefit the most from microgravity crystallization. In order to assess this theory the Delta L hardware is being developed so that macromolecule crystal growth rates can be measured in microgravity. Crystal growth rate is defined as the change or delta in crystal size (defined as a characteristic length, L) over time; hence the name of the hardware. Delta L will consist of an optics, a fluids, and a data acquisition sub-assemblies. The optics assembly will consist of a video microscope camera mounted on three axis computer controlled translation stages. The fluids assembly consists of macromolecule and precipitant reservoirs, a temperature controlled growth cell and waste container, The data acquisition is achieved by using a frame-gabber, with images being stored on a hard drive. In operation, macromolecule and precipitant solution will be injected into the temperature controlled growth cell. As macromolecule crystals grow, the video microscope camera controlled by the translation stages, will be used to locate and record images of individual crystals, returning to the same crystals at specific time intervals. The images will be stored on the hard drive and used to calculate the crystal growth rate. To prevent vibrations interfering in the crystal growth rate measurements (Snell et al., 1997) Delta L will be used in connection with the Glovebox Integrated Microgravity Isolation Technology (g-LIMIT) inside the Microgravity Science Glovebox (MSG), onboard the International Space Station (ISS).

Judge, Russell A.↗

Materials Science

The Materials Science Program is structured so that NASA s headquarters is responsible for the program content and selection, through the Enterprise Scientist, and MSFC provides for implementation of ground and flight programs with a Discipline Scientist and Discipline Manager. The Discipline Working Group of eminent scientists from outside of NASA acts in an advisory capacity and writes the Discipline Document from which the NRA content is derived. The program is reviewed approximately every three years by groups such as the Committee on Microgravity Research, the National Materials Advisory Board, and the OBPR Maximization and Prioritization (ReMaP) Task Force. The flight program has had as many as twenty-six principal investigators (PIs) in flight or flight definition stage, with the numbers of PIs in the future dependent on the results of the ReMaP Task Force and internal reviews. Each project has a NASA-appointed Project Scientist, considered a half-time job, who assists the PI in understanding and preparing for internal reviews such as the Science Concept Review and Requirements Definition Review. The Project Scientist also insures that the PI gets the maximum science support from MSFC, represents the PI to the MSFC community, and collaborates with the Project Manager to insure the project is well-supported and remains vital. Currently available flight equipment includes the Materials Science Research Rack (MSRR-1) and Microgravity Science Glovebox. Ground based projects fall into one or more of several categories. Intellectual Underpinning of Flight Program projects include theoretical studies backed by modeling and computer simulations; bring to maturity new research, often by young researchers, and may include preliminary short duration low gravity experiments in the KC-135 aircraft or drop tube; enable characterization of data sets from previous flights; and provide thermophysical property determinations to aid PIs. Radiation Shielding and preliminary In Situ Resource Utilization (ISRU) studies work towards future long duration missions. Biomaterials support materials issues affecting crew health. Nanostructured Materials are currently considered to be maturing new research, and Advanced Materials for Space Transportation has as yet no PIs. PIs are assigned a NASA Technical Monitor to maintain contact, a position considered to be a 5 percent per PI effort. Currently 33 PIs are supported on the 1996 NRA, which is about to expire, and 59 on the 1998 NRA. Two new NRAs, one for Radiation Shielding and one for Materials Science for Advanced Space Propulsion are due to be announced by the 2003 fiscal year. MSFC has a number of facilities supporting materials science. These include the Microgravity Development Laboratory/SD43; Electrostatic Levitator Facility; SCN Purification Facility; Electron Microscope/Microprobe Facility; Static and Rotating Magnetic Field Facility; X-Ray Diffraction Facility; and the Furnace Development Laboratory.

Source record↗

Fabricating Radial Groove Gratings Using Projection Photolithography

Projection photolithography has been used as a fabrication method for radial grove gratings. Use of photolithographic method for diffraction grating fabrication represents the most significant breakthrough in grating technology in the last 60 years, since the introduction of holographic written gratings. Unlike traditional methods utilized for grating fabrication, this method has the advantage of producing complex diffractive groove contours that can be designed at pixel-by-pixel level, with pixel size currently at the level of 45 45 nm. Typical placement accuracy of the grating pixels is 10 nm over 30 nm. It is far superior to holographic, mechanically ruled or direct e-beam written gratings and results in high spatial coherence and low spectral cross-talk. Due to the smooth surface produced by reactive ion etch, such gratings have a low level of randomly scattered light. Also, due to high fidelity and good surface roughness, this method is ideally suited for fabrication of radial groove gratings. The projection mask is created using a laser writer. A single crystal silicon wafer is coated with photoresist, and then the projection mask, with its layer of photoresist, is exposed for patterning in a stepper or scanner. To develop the photoresist, the fabricator either removes the exposed areas (positive resist) of the unexposed areas (negative resist). Next, the patterned and developed photoresist silicon substrate is subjected to reactive ion etching. After this step, the substrate is cleaned. The projection mask is fabricated according to electronic design files that may be generated in GDS file format using any suitable CAD (computer-aided design) or other software program. Radial groove gratings in off-axis grazing angle of incidence mount are of special interest for x-ray spectroscopy, as they allow achieving higher spectral resolution for the same grating area and have lower alignment tolerances than traditional in-plane grating scheme. This is especially critical for NASA Constellation- X project that will utilize hundreds of gratings all of which need to be precisely aligned for x-ray observation of space.

Iazikov, Dmitri↗

Phase holograms in PMMA with proximity effect correction

Complex computer generated phase holograms (CGPH's) have been fabricated in PMMA by partial e-beam exposure and subsequent partial development. The CGPH was encoded as a sequence of phase delay pixels and written by the JEOL JBX-5D2 E-beam lithography system, a different dose being assigned to each value of phase delay. Following carefully controlled partial development, the pattern appeared rendered in relief in the PMMA, which then acts as the phase-delay medium. The exposure dose was in the range 20-200 micro-C/sq cm, and very aggressive development in pure acetone led to low contrast. This enabled etch depth control to better than plus or minus lambda(sub vis)/60. That result was obtained by exposing isolated 50 micron square patches and measuring resist removal over the central area where the proximity effect dose was uniform and related only to the local exposure. For complex CGPH's with pixel size of the order of the e-beam proximity effect radius, the patterns must be corrected for the extra exposure caused by electrons scattered back up out of the substrate. This has been accomplished by deconvolving the two-dimensional dose deposition function with the desired dose pattern. The deposition function, which plays much the same role as an instrument response function, was carefully measured under the exact conditions used to expose the samples. The devices fabricated were designed with 16 equal phase steps per retardation cycle, were up to 1 cm square, and consisted of up to 100 million 0.3-2.0 micron square pixels. Data files were up to 500 MB long and exposure times ranged to tens of hours. A Fresnel phase lens was fabricated that had diffraction limited optical performance with better than 85 percent efficiency.

Maker, Paul D.↗

A model predicting the evolution of ice particle size spectra and radiative properties of cirrus clouds. Part 2: Dependence of absorption and extinction on ice crystal morphology

This study builds upon the microphysical modeling described in Part 1 by deriving formulations for the extinction and absorption coefficients in terms of the size distribution parameters predicted from the micro-physical model. The optical depth and single scatter albedo of a cirrus cloud can then be determined, which, along with the asymmetry parameter, are the input parameters needed by cloud radiation models. Through the use of anomalous diffraction theory, analytical expressions were developed describing the absorption and extinction coefficients and the single scatter albedo as functions of size distribution parameters, ice crystal shapes (or habits), wavelength, and refractive index. The extinction coefficient was formulated in terms of the projected area of the size distribution, while the absorption coefficient was formulated in terms of both the projected area and mass of the size distribution. These properties were formulated as explicit functions of ice crystal geometry and were not based on an 'effective radius.' Based on simulations of the second cirrus case study described in Part 1, absorption coefficients predicted in the near infrared for hexagonal columns and rosettes were up to 47% and 71% lower, respectively, than absorption coefficients predicted by using equivalent area spheres. This resulted in single scatter albedos in the near-infrared that were considerably greater than those predicted by the equivalent area sphere method. Reflectances in this region should therefore be underestimated using the equivalent area sphere approach. Cloud optical depth was found to depend on ice crystal habit. When the simulated cirrus cloud contained only bullet rosettes, the optical depth was 142% greater than when the cloud contained only hexagonal columns. This increase produced a doubling in cloud albedo. In the near-infrared (IR), the single scatter albedo also exhibited a significant dependence on ice crystal habit. More research is needed on the geometrical properties of ice crystals before the influence of ice crystal shape on cirrus radiative properties can be adequately understood. This study provides a way of coupling the radiative properties of absorption, extinction, and single scatter albedo to the microphysical properties of cirrus clouds. The dependence of extinction and absorption on ice crystal shape was not just due to geometrical differences between crystal types, but was also due to the effect these differences had on the evolution of ice particle size spectra. The ice particle growth model in Part 1 and the radiative properties treated here are based on analytical formulations, and thus represent a computationally efficient means of modeling the microphysical and radiative properties of cirrus clouds.

Mitchell, David L.↗