Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave stretching”

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 91 records · Page 5

Transport out of the lower stratospheric Arctic vortex by Rossby wave breaking

The fine-scale structure in lower stratospheric tracer transport during the period of the two Arctic Airborne Stratospheric Expeditions (January and February 1989; December 1991 to March 1992) is investigated using contour advection with surgery calculations. These calculations show that Rossby wave breaking is an ongoing occurrence during these periods and that air is ejected from the polar vortex in the form of long filamentary structures. There is good qualitative agreement between these filaments and measurements of chemical tracers taken aboard the NASA ER-2 aircraft. The ejected air generally remains filamentary and is stretched and mixed with midlatitude air as it is wrapped around the vortex. This process transfers vortex air into midlatitudes and also produces a narrow region of fine-scale filaments surrounding the polar vortex. Among other things, this makes it difficult to define a vortex edge. The calculations also show that strong stirring can occur inside as well as outside the vortex.

Waugh, D. W.↗

Solution of Third Computational Workshop Internal Propagation Problems using Low Order Schemes

The problems solved are propagation of sound waves through a transonic nozzle and shock-sound interaction. The first problem is solved using a MacCormack scheme. The second problem is solved by perturbation of a conservative Euler equation solution to the steady state. The problems are solved on evenly spaced grids. While solutions were found, the methods selected and not using stretched grids lead to solutions which do not compare well with those found using more accurate schemes. However, the perturbation of the mean flow scheme used to solve the second problem shows promise.

Miles, Jeffrey Hilton↗

Soleus H-reflex gain in humans walking and running under simulated reduced gravity

The Hoffmann (H-) reflex is an electrical analogue of the monosynaptic stretch reflex, elicited by bypassing the muscle spindle and directly stimulating the afferent nerve. Studying H-reflex modulation provides insight into how the nervous system centrally modulates stretch reflex responses.A common measure of H-reflex gain is the slope of the relationship between H-reflex amplitude and EMG amplitude. To examine soleus H-reflex gain across a range of EMG levels during human locomotion, we used simulated reduced gravity to reduce muscle activity. We hypothesised that H-reflex gain would be independent of gravity level.We recorded EMG from eight subjects walking (1.25 m s-1) and running (3.0 m s-1) at four gravity levels (1.0, 0.75, 0.5 and 0.25 G (Earth gravity)). We normalised the stimulus M-wave and resulting H-reflex to the maximal M-wave amplitude (Mmax) elicited throughout the stride to correct for movement of stimulus and recording electrodes relative to nerve and muscle fibres. Peak soleus EMG amplitude decreased by ~30% for walking and for running over the fourfold change in gravity. As hypothesised, slopes of linear regressions fitted to H-reflex versus EMG data were independent of gravity for walking and running (ANOVA, P > 0.8). The slopes were also independent of gait (P > 0.6), contrary to previous studies. Walking had a greater y-intercept (19.9% Mmax) than running (-2.5% Mmax; P < 0.001). At all levels of EMG, walking H-reflex amplitudes were higher than running H-reflex amplitudes by a constant amount. We conclude that the nervous system adjusts H-reflex threshold but not H-reflex gain between walking and running. These findings provide insight into potential neural mechanisms responsible for spinal modulation of the stretch reflex during human locomotion.

Non-NASA Center↗

A small-scale turbulence model

A model for the small-scale structure of turbulence is reformulated in such a way that it may be conveniently computed. The model is an ensemble of randomly oriented structured two dimensional vortices stretched by an axially symmetric strain flow. The energy spectrum of the resulting flow may be expressed as a time integral involving only the enstrophy spectrum of the time evolving two-dimensional cross section flow, which may be obtained numerically. Examples are given in which a k(exp -5/3) spectrum is obtained by this method without using large wave number asymptotic analysis. The k(exp -5/3) inertial range spectrum is shown to be related to the existence of a self-similar enstrophy preserving range in the two-dimensional enstrophy spectrum. The results are insensitive to time dependence of the strain-rate, including even intermittent on-or-off strains.

Lundgren, T. S.↗

Current and Future Plans of the NASA Data Assimilation Office (DAO)

The mission of the Data Assimilation Office (DAO) is to advance the state of the art of data assimilation and produce research-quality assimilated data sets which make optimal use of space-based observations. Development efforts over the last few years have focused on delivering a production data assimilation system in support of NASA's Terra platform. That system, called the Goddard Earth Observing System - version 2 or GEOS-2, represents a major upgrade to the baseline GEOS-1 system employed in NASA's first reanalysis effort. GEOS-2 includes a physical-space three dimensional variational analysis algorithm (the Physical-space Statistical Analysis System or PSAS) and numerous improvements to the general circulation model. The latter include a Soil-Vegetation-Atmosphere Transfer (SVAT) land surface scheme, a level 2.5 moist turbulence scheme and new Short Wave (SW) and Long Wave (LW) radiation code. The system also includes an off-line ozone assimilation system, and the capability to assimilate scatterometer surface winds, and TIROS Operational Vertical Sounder (TOVS) and Special Sensor Microwave Imager (SSM/I) moisture data. GEOS-2 is currently run at 1 degree horizontal resolution and 48 levels extending to O.Olmb. Experimental versions of GEOS-2 are run with a global stretched grid allowing enhanced (e.g. 1/4 deg) regional resolution. Other capabilities being developed include, the assimilation of Tropical Rainfall Measuring Mission (TRMM) precipitation and Global Positioning System (GPS) data, an off-line land surface assimilation system, and a retrospective analysis scheme. The DAO is also engaged in a number of collaborative efforts to help accelerate the development of the next generation data assimilation system. These include, a joint modeling effort between the DAO and NCAR/CGDD to develop a new Global Circulation Model (GCM), and a Department of Energy Lawrence Livermore National Laboratory (DOE/LLNL) collaboration on model parallelization. Plans for the next reanalysis will be discussed in the context of current and near term system quality and computing capabilities, and the need for multiple reanalysis products.

Atlas, Robert↗

Stereo Imaging Velocimetry of Mixing Driven by Buoyancy Induced Flow Fields

Mixing of two fluids generated by steady and particularly g-jitter acceleration is fundamental towards the understanding of transport phenomena in a microgravity environment. We propose to carry out flight and ground-based experiments to quantify flow fields due to g-jitter type of accelerations using Stereo Imaging Velocimetry (SIV), and measure the concentration field using laser fluorescence. The understanding of the effects of g-jitter on transport phenomena is of great practical interest to the microgravity community and impacts the design of experiments for the Space Shuttle as well as the International Space Station. The aim of our proposed research is to provide quantitative data to the community on the effects of g-jitter on flow fields due to mixing induced by buoyancy forces. The fundamental phenomenon of mixing occurs in a broad range of materials processing encompassing the growth of opto-electronic materials and semiconductors, (by directional freezing and physical vapor transport), to solution and protein crystal growth. In materials processing of these systems, crystal homogeneity, which is affected by the solutal field distribution, is one of the major issues. The understanding of fluid mixing driven by buoyancy forces, besides its importance as a topic in fundamental science, can contribute towards the understanding of how solutal fields behave under various body forces. The body forces of interest are steady acceleration and g-jitter acceleration as in a Space Shuttle environment or the International Space Station. Since control of the body force is important, the flight experiment will be carried out on a tunable microgravity vibration isolation mount, which will permit us to precisely input the desired forcing function to simulate a range of body forces. To that end, we propose to design a flight experiment that can only be carried out under microgravity conditions to fully exploit the effects of various body forces on fluid mixing. Recent flight experiments, by the P.I. through collaboration with the Canadian Space Agency (STS-85, August 1997), aimed at determining the stability of the interface between two miscible liquids inside an enclosure show that a long liquid column (5 cm) under microgravity isolation conditions can be stable, i.e. the interface remains sharp and vertical over a short time scale; thus transport occurs by molecular mass diffusion. On the other hand, when the two liquids were excited from a controlled vibration source (Microgravity Vibration Isolation Mount) two to four mode large amplitude quasi-stationary waves were observed. The data was limited to CCD recording of the dynamics of the interface between the two fluids. We propose to carry out flight experiments to quantify the dynamics of the flow field using Stereo Imaging Velocimetry and measure the concentration field using laser fluorescence. The results will serve as a basis to understand effects of g-jitter on transport phenomena, in this case mass diffusion. As the measurement of the kinematics of the flow field will shed light on the instability mechanism. The research will allow measurement of the flow field in microgravity environment to prove two hypotheses: (1) Maxwell's hypothesis: finite convection always exists in diffusing systems, and (2) Quasi-stationary waves inside a bounded enclosure in a microgravity environment is generated by Kelvin-Helmholtz instability; resonance of the interface which produces incipient mixing is due to Rayleigh-Taylor instability. The first hypothesis can be used as a benchmark experiment to illustrate diffusive mixing. The second hypothesis will lead to the understanding of g-jitter effects on buoyancy driven flow fields which occur in many situations involving materials processing, and other basic fluid physics phenomena. In addition, the second hypothesis will also provide insight in how Rayleigh-Taylor and Kelvin-Helmholtz instabilities propagate concentration fronts during mixing. Measurement of the flow field using SIV is important because it is the flow field which causes instability at the interface between the two fluids. Mixing driven by buoyancy induced flow fields will be addressed both experimentally and computationally. The experimental effort will address the kinematics of mixing: stretching, transport and chaos. Quantification of the mechanisms of mixing will consists of measuring the flow field using the SIV system at Glenn and capturing the dynamics of the interface, to measure mass transport, using a CCD camera. These experiments will be carried out within the framework of Earth's gravity and g-jitter microgravity acceleration as in a Space Shuttle environment or the International Space Station. The g-jitter will be induced and controlled using a tunable vibration isolation platform to isolate against vibration as well as input periodic and random vibration to the system. The parametric range of the microgravity experiment will be extended from the experiments on STS-85 to investigate higher mode quasi-stationary waves (8 to 12), as well as resonance regions which leads to chaos and turbulence. Ground-based experiments will focus on effects of vibration on stably stratified fluid layers in order to scale for possible scenarios in a microgravity environment. These vibrations will be subjected perpendicular to the concentration field on the ground since the parallel case can only be carried out in a microgravity environment. The concept of dynamical similarity will be applied to tune the experiments as closely as possible to a Space Shuttle environment or the International Space Station. The computational effort will take advantage of the Computational Laboratory at Glenn to corroborate the experimental findings with predictions of the dynamics of the flow field using the codes FLUENT (finite difference based) and FIDAP (finite element based). We will investigate two important cases, single-fluid model to address dilute systems with negligible jump in viscosity and the more general two-fluid model which accounts for finite jump in viscosity. Apart from its microgravity relevance, this experiment is well suited to study dynamics in nonlinear systems.

Duval, W. M. B.↗

Enstrophy budget in decaying compressible turbulence

The interaction of vorticity and shock waves in decaying compressible turbulence is investigated by direct numerical simulation of the Navier-Stokes equation using a spectral method with 64-cubed collocation points. The Taylor microscale Reynolds number decreases from 36 to 17 and the rms Mach number from 0.8 to 0.4 during the decay. It is observed that vorticity is created inside curved shocks through baroclinic interaction and intensified by the compression of fluid elements. Vortex stretching is, on the other hand, dominant outside shocks and is responsible for enhancing the strength of vorticity. A dipole of vorticity is typically created at the vertex of curved shocks, the structure of which is explained in terms of jump relations of velocity derivatives across a curved shock.

Kida, Shigeo↗

Dynamics and control of coherent structures in the turbulent wall layer: An overview

We expand the velocity field in the vicinity of the wall in empirical eigenfunctions obtained from experiment. Truncating our system and using Galerkin projection, we obtain a closed set of non-linear ordinary differential equations with ten degrees of freedom. We find a rich dynamical behavior, including in particular a heteroclinic attracting orbit giving rise to intermittency. The intermittent jump from one attracting point to the other resembles in many respects the bursts observed in experiments. Specifically, the time between jumps and the duration of the jumps, is approximately that observed in a burst; the jump begins with the formation of a narrowed and intensified updraft, like the ejection phase of a burst, and is followed by a gentle, diffuse downdraft like the sweep phase of a burst. The magnitude of the Reynolds stress spike produced during a burst is limited by our truncation. The behavior is quite robust, much of it being due to the symmetries present (Aubry's group has examined dimensions up to 128 with persistence of the global behavior). We have examined eigenvalues and coefficients obtained from experiment, and from exact simulation, which differ in magnitude. Similar behavior is obtained in both cases; in the latter case, the heteroclinic orbits connect limit cycles instead of fixed points, corresponding to cross-stream waving of the streamwise rolls. The bifurcation diagram remains structurally similar, but somewhat distorted. The role of the pressure term is made clear - it triggers the intermittent jumps, which otherwise would occur at longer and longer intervals, as the system trajectory is attracted closer and closer to the heteroclinic cycle. The pressure term results in the jumps occurring at essentially random times, and the magnitude of the signal determines the average timing. Stretching of the wall region shows that the model is consistent with observations of polymer drag reduction. Change of the third order coefficients, corresponding to acceleration or deceleration of the mean flow, changes the heteroclinic cycles from attracting to repelling, increasing or decreasing the stability, in agreement with observations. The existence of fixed points is an artifact introduced by the projection; however, a decoupled model still displays the rich dynamics. Numerous assumptions made in Aubry et al. (1988) can now be proved exactly. Feeding back eigenfuncitons with the proper phase can delay the bursting, (the heteroclinic jump to the other fixed point), decreasing the drag. It is also possible to speed up the bursting, increasing mixing to control separation. Our approach is optimal for short time tracking in control.

Berkooz, Gal↗

Lightning induced currents in aircraft wiring using low level injection techniques

Various techniques were studied to predict the transient current induced into aircraft wiring bundles as a result of an aircraft lightning strike. A series of aircraft measurements were carried out together with a theoretical analysis using computer modeling. These tests were applied to various aircraft and also to specially constructed cylinders installed within coaxial return conductor systems. Low level swept frequency CW (carrier waves), low level transient and high level transient injection tests were applied to the aircraft and cylinders. Measurements were made to determine the transfer function between the aircraft drive current and the resulting skin currents and currents induced on the internal wiring. The full threat lightning induced transient currents were extrapolated from the low level data using Fourier transform techniques. The aircraft and cylinders used were constructed from both metallic and CFC (carbon fiber composite) materials. The results show the pulse stretching phenomenon which occurs for CFC materials due to the diffusion of the lightning current through carbon fiber materials. Transmission Line Matrix modeling techniques were used to compare theoretical and measured currents.

Stevens, E. G.↗

An evaluation of short-wave-infrared (SWIR) data from the AVIRIS and GEOSCAN instruments for mineralogical mapping at Cuprite, Nevada

An evaluation was performed on SWIR (2000-2400 nm) data from two airborne remote sensing systems for discriminating and identifying alteration minerals at Cuprite, Nevada. The data were acquired by the NASA Airborne Visible/InfraRed Imaging Spectrometer (AVIRIS) and the GEOSCAN Mk II multispectral scanner. The evaluation involved comparison of processed imagery and image-derived spectra with existing alteration maps and laboratory spectra of rock samples from Cuprite. Results indicate that both the AVIRIS and GEOSCAN data permit the discrimination of areas of alunite, buddingtonite, kaolinite, and silicification using color composite images formed from three SWIR bands processed with either the decorrelation stretch or a log residual algorithm. The laboratory spectral features alunite, kaolinite and buddingtonite could be seen clearly only in the log residual processed AVIRIS data. However, this does not preclude their identification with the GEOSCAN data.

Hook, Simon J.↗

Premixed Gas Combustion: An Excitable System

Rotating spiral and target patterns have been observed experimentally on freely-propagating premixed gas flames in large diameter tubes at normal gravity (1-g). These modes of propagation occur in near-limit mixtures which have a Lewis number (Le, defined as the ratio of the thermal diffusivity of the cold mixture to the mass diffusivity of the scarce component into the mixture) sufficiently greater than one. However, at 1-g, buoyant flows strongly distort the flame curvature, hydrodynamics (thus stretch) and convective transport of species and heat. In turn, these alter the critical Le required for onset of instability. To isolate and better understand the mechanisms which drive the observed patterns and their dynamics, 1-g and microgravity (micro-g) experiments are being conducted to determine: (1) the structure and dynamics of the patterns, (2) a map of the critical Le and heat loss for their occurrence, (3) the relative significance of the chemical kinetics, and (4) the effect of curvature (local wave and global flame front) on wave propagation. With this in hand, we will be better prepared to discuss an additional mode, a state of 'chemical turbulence,' which seems to be the ultimate fate of many of these near-limit flames prior to extinction.

Pearlman, Howard↗

A Model for Flux Rope Formation and Disconnection in Pseudostreamer Coronal Mass Ejections

Coronal mass ejections (CMEs) from pseudostreamers represent a significant fraction of large-scale eruptions from the Sun. In some cases, these CMEs take a narrow jet-like form reminiscent of coronal jets; in others, they have a much broader fan-shaped morphology like CMEs from helmet streamers. We present results from a magnetohydrodynamic simulation of a broad pseudostreamer CME. The early evolution of the eruption is initiated through a combination of breakout interchange reconnection at the overlying null point and ideal instability of the flux rope that forms within the pseudostreamer. This stage is characterised by a rolling motion and deflection of the flux rope toward the breakout current layer. The stretching out of the strapping field forms a flare current sheet below the flux rope; reconnection onset there forms low-lying flare arcade loops and the two-ribbon flare footprint. Once the CME flux rope breaches the rising breakout current layer, interchange reconnection with the external open field disconnects one leg from the Sun. This induces a whip-like rotation of the flux rope, generating the unstructured fan shape characteristic of pseudostreamer CMEs. Interchange reconnection behind the CME releases torsional Alfv´en waves and bursty dense outflows into the solar wind. Our results demonstrate that pseudostreamer CMEs follow the same overall magnetic evolution as coronal jets, although they present different morphologies of their ejecta. We conclude that pseudostreamer CMEs should be considered a class of eruptions that are distinct from helmet streamer CMEs, in agreement with previous observational studies.

Sun: corona↗

Observations of field-aligned density fluctuations in the inner solar wind

This paper reports the results of radio wave scattering observations which show highly anisotropic density microstructure in the inner solar wind. These observations were made in October 1983 and 1985, using the National Radio Astronomy Observatory's Very Large Array. Heliocentric distances of the observations ranged from 2.2. to about 13 solar radii. The axial ratio of the density fluctuations, projected onto the plane of the sky, increased from about 4 at 10 solar radii to about 14 at 2.2 solar radii. The major axis appeared to be field aligned; that is, the irregularities were stretched out approximately along the radial. These results are, in general, consistent with the results of earlier observations. The present observations differ in that they were taken closer to the sun, they were sensitive to larger scale structures (up to about 35 km), and showed much higher anisotropy. Combining these data with previously taken data strongly indicates that the anisotropy is scale dependent; scales greater than 10 km appear to be more anisotropic than those less than 2 km.

Armstrong, John W.↗

The vibrational spectrum of H2O3: An ab initio investigation

Theoretically determined frequencies and absorption intensities are reported for the vibrational spectrum of the covalent HOOOH and hydrogen bonded HO---HOO intermediates that may form in the reaction of the hydroxyl and hydroperoxyl radicals. Basis sets of DZP quality, augmented by diffuse and second sets of polarization functions have been used with CASSCF wave functions. The calculated harmonic vibrational frequencies of HOOOH have been corrected with empirical factors and presented in the form of a 'stick' spectrum. The oxygen backbone vibrations, predicted to occur at 519, 760, and 870 cm(exp -1), are well separated from most interferences, and may be the most useful for the species' identification. In the case of the hydrogen bonded isomer, emphasis has been placed upon prediction of the shifts in the intramolecular vibrational frequencies that take place upon formation of the complex. In particular, the HO stretch and HOO bend of HO2 are predicted to have shifts of -59 and 53 cm(exp -1), respectively, which should facilitate their identification. It is also noted that the antisymmetric stretching frequency of the oxygen backbone in HOOOH exhibits a strong sensitivity to the degree of electron correlation, such as has been previously observed for the same mode in ozone.

Jackels, Charles F.↗

Magnetic Fluctuations in the Jovian Magnetosphere

The engine that drives the jovian magnetosphere is the mass added to the Io ion torus, accelerated to corotational velocities by field-aligned currents that couple the Io torus to the jovian ionosphere. The mass of the torus builds up to an amount that the magnetic forces cannot contain and the plasma, first slowly and then more rapidly, drifts outward. Numerous authors have treated this problem based first on the observations of the Pioneer 10 and 11 flybys; then on Voyager 1 and 2, and Ulysses; and finally most recently the Galileo orbiter. The initial observations revealed the now familiar magnetodisk, in which the field above and below the magnetic equator became quite radial in orientation and much less dipolar. The Galileo observations show this transformation to occur on average at 24 R(sub J) and to often be quite abrupt. These observations are consistent with outward transport of magnetized plasma that moves ever faster radially until about 50 R(sub J) on the nightside where the field lines stretch to the breaking point, reconnection occurs, and plasma and field islands are transported down the tail ultimately removing the mass from the magnetosphere that Io had deposited deep in the inner torus. The reconnection process creates empty flux tubes connected to Jupiter that are buoyant and thought to float inward and replace the flux carried out with the torus plasma. As described above, the jovian magnetosphere could very well be in a state of steady laminar circulation, but indeed it is not. The process is very unsteady and the wave levels can be very intense. The existence of these waves in turn can lead to processes that compete with the radial circulation pattern in removing plasma from the system. These waves can scatter particles so that they precipitate into the ionosphere. This process should be important in the Io torus where the atmospheric loss cone is relatively large and becomes less important as the loss cone decreases in size with radial distance. However, the Io torus is relatively quiet compared to the region outside the torus and it is not obvious without studying this scattering carefully whether the loss in the torus or out of the torus is greater and whether it can act rapidly enough to compete with the radial transport of ions to the tail in the life cycle of the mass added at Io. Closer to Io the ion cyclotron waves are most intense and possibly are associated with the losses in the Io flux tube. The waves are also diagnostic of both the Io atmospheric composition and the size and strength of the massloading process.

Russell, Christopher T.↗

A Meteorological Overview of the TC4 Mission

The TC4 mission in Central America during summer 2007 examined convective transport into the tropical Upper Troposphere/Lower Stratosphere (UTLS) and the evolution of cirrus clouds. The tropical tropopause layer (TTL) circulation is dominated by the Asian monsoon anticyclone and westward winds that stretch from the western Pacific into the Atlantic. During TC4, TTL westward flow over Central America was stronger than normal. Incidence of cold clouds over the Central American region was the third lowest out of 34 years sampled. The major factor was an incipient La Nina, specifically anomalously cold temperatures off the Pacific Coast of South America. Weakness in the low level Caribbean jet caused a shift in the coldest clouds from the Caribbean to the Pacific side of Central America. The character of tropopause temperature variability was that of upward propagating waves generated by local and nonlocal convection. These waves produced tropopause temperature variations of 3 K, with peak-to-peak variations of 8 K. At low levels in Central America, flow from the Sahara desert predominated; further south, the air came from the Amazon region. Convectively influenced air in the upper troposphere came from Central America, the northern Amazon region, the Atlantic ITCZ, and the North American monsoon. In the TTL, Asian and African convection affected the observed air masses. North of 10N in the Central American TTL, African and Asian convection may have contributed as much to the air masses as Central and South American convection. South of 8N, Asian and African convection had far less impact.

Pfister, L.↗

The terminal area simulation system. Volume 1: Theoretical formulation

A three-dimensional numerical cloud model was developed for the general purpose of studying convective phenomena. The model utilizes a time splitting integration procedure in the numerical solution of the compressible nonhydrostatic primitive equations. Turbulence closure is achieved by a conventional first-order diagnostic approximation. Open lateral boundaries are incorporated which minimize wave reflection and which do not induce domain-wide mass trends. Microphysical processes are governed by prognostic equations for potential temperature water vapor, cloud droplets, ice crystals, rain, snow, and hail. Microphysical interactions are computed by numerous Orville-type parameterizations. A diagnostic surface boundary layer is parameterized assuming Monin-Obukhov similarity theory. The governing equation set is approximated on a staggered three-dimensional grid with quadratic-conservative central space differencing. Time differencing is approximated by the second-order Adams-Bashforth method. The vertical grid spacing may be either linear or stretched. The model domain may translate along with a convective cell, even at variable speeds.

Proctor, F. H.↗

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↗