Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Convection lines”

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 883 records · Page 49

On the Ion Precipitation due to Field Line Curvature (FLC) and EMIC Wave Scattering and Their Subsequent Impact on Ionospheric Electrodynamics

Both field line curvature (FLC) and electromagnetic ion cyclotron (EMIC) wave scattering are believed to be associated with energetic ion precipitation, but their relative contributions to ionospheric ion precipitation and subsequent effects are still unclear. Here, in this study, by using a kinetic ring current model, we investigate their impact on the ionosphere from two aspects: the global distribution of ion precipitation and resulting ionospheric conductance. Our results show that the intensity and coverage of ion precipitation due to EMIC waves are larger than that due to FLC scattering, while the latter mostly contributes to ion precipitation in outer regions (L > 4–5). We then estimate the conductance with empirical models using the simulated ion precipitation energy flux. When the EMIC wave associated proton precipitation is included, the conductance is significantly enhanced in the dusk-to-midnight sector and has a wide magnetic latitude (MLAT) range from around 52° to 62°, considerably altering the electric potential in the dusk sector, which further influences particle dynamics in the magnetosphere. On the contrary, the proton precipitation caused by FLC scattering only occurs at higher latitudes above MLAT = 60° and the corresponding conductance is slightly enhanced at midnight, with negligible influence on the convective electric potential. Although electron precipitation-associated conductance is predominant globally, the results show that proton precipitation can also play an important role in ionospheric electrodynamics, especially when EMIC wave-associated precipitating proton energy flux exceeds that of electrons in the dusk sector, a region where many subauroral coupling processes occur.

79 ASTRONOMY AND ASTROPHYSICS↗

An Unusual Transmission Spectrum for the Sub-Saturn KELT-11b Suggestive of a Subsolar Water Abundance

We present an optical-to-infrared transmission spectrum of the inflated sub-Saturn KELT-11b measured with the Transiting Exoplanet Survey Satellite (TESS), the Hubble Space Telescope (HST) Wide Field Camera 3 G141 spectroscopic grism, and the Spitzer Space Telescope (Spitzer) at 3.6 μm, in addition to a Spitzer 4.5 μm secondary eclipse. The precise HST transmission spectrum notably reveals a low-amplitude water feature with an unusual shape. Based on free-retrieval analyses with varying molecular abundances, we find strong evidence for water absorption. Depending on model assumptions, we also find tentative evidence for other absorbers (HCN, TiO, and AlO). The retrieved water abundance is generally ≲0.1× solar (0.001–0.7× solar over a range of model assumptions), several orders of magnitude lower than expected from planet formation models based on the solar system metallicity trend. We also consider chemical-equilibrium and self-consistent 1D radiative-convective equilibrium model fits and find that they, too, prefer low metallicities ([M/H] ≲ -2, consistent with the free-retrieval results). However, all of the retrievals should be interpreted with some caution because they either require additional absorbers that are far out of chemical equilibrium to explain the shape of the spectrum or are simply poor fits to the data. Finally, we find that the Spitzer secondary eclipse is indicative of full heat redistribution from KELT-11b’s dayside to nightside, assuming a clear dayside. These potentially unusual results for KELT-11b’s composition are suggestive of new challenges on the horizon for atmosphere and formation models in the face of increasingly precise measurements of exoplanet spectra.

79 ASTRONOMY AND ASTROPHYSICS↗

Research activities in nuclear astrophysics and related areas

NASA/GRO grant NAG 5-2081, at the University of Chicago, has provided support for a broad program of theoretical research in nuclear astrophysics and related areas, with regard to gamma-ray and hard X-ray emission from classical nova explosions. This research emphasized the possible detection of 22Na gamma-ray line emission from nearby novae involving ONeMg white dwarfs, the detailed examination of 26Al production in novae, and the possible detection of the predicted early gamma ray emission from novae that arises from the decay of the short lived, positron emitting isotopes of CNO elements. Studies of nova related problems have consumed an increasing fraction of the Principal Investigator's research efforts over the past decade. Current research addresses problems associated with the standard model for the outbursts of the classical novae: the occurrence of thermonuclear runaways (TNR) in the accreted hydrogen rich envelopes on white dwarfs in close binary systems (see, e.g., the reviews by Truran 1982; and Shara 1989). Research in progress and planned for the next three years has three main objectives: (1) to gain an improved understanding of the early evolution of the light curves of, particularly, the fastest novae; (2) to gain an improved understanding of the relative importance of the various possible mechanisms of envelope hydrogen depletion (e.g. winds, common envelope driven mass loss, and nuclear burning) to the long term evolution of novae in outburst; and (3) to seek to provide a somewhat more definitive statement of the role of classical novae in nucleosynthesis. Our proposed 2-D studies of convection during the early phases of the TNR and our systematic attempt to incorporate an improved treatment of radiation hydrodynamics into the hydrodynamic code utilized in our calculations, are particularly relevant to the first of these objectives. Further 2-D studies of the effects of common envelope evolution are intended to provide more realistic constraints on the mass depletion mechanisms. Finally, detailed calculations of the thermonuclear history of the matter ejected in novae will be carried out for representative nova configurations involving both carbon-oxygen (CO) and oxygen-neon-magnesium (ONeMg) white dwarfs.

Source record↗

Uncertainty Determination for Aeroheating in Uranus and Saturn Probe Entries by the Monte Carlo Method

The 2013-2022 Decaedal survey for planetary exploration has identified probe missions to Uranus and Saturn as high priorities. This work endeavors to examine the uncertainty for determining aeroheating in such entry environments. Representative entry trajectories are constructed using the TRAJ software. Flowfields at selected points on the trajectories are then computed using the Data Parallel Line Relaxation (DPLR) Computational Fluid Dynamics Code. A Monte Carlo study is performed on the DPLR input parameters to determine the uncertainty in the predicted aeroheating, and correlation coefficients are examined to identify which input parameters show the most influence on the uncertainty. A review of the present best practices for input parameters (e.g. transport coefficient and vibrational relaxation time) is also conducted. It is found that the 2(sigma) - uncertainty for heating on Uranus entry is no more than 2.1%, assuming an equilibrium catalytic wall, with the uncertainty being determined primarily by diffusion and H(sub 2) recombination rate within the boundary layer. However, if the wall is assumed to be partially or non-catalytic, this uncertainty may increase to as large as 18%. The catalytic wall model can contribute over 3x change in heat flux and a 20% variation in film coefficient. Therefore, coupled material response/fluid dynamic models are recommended for this problem. It was also found that much of this variability is artificially suppressed when a constant Schmidt number approach is implemented. Because the boundary layer is reacting, it is necessary to employ self-consistent effective binary diffusion to obtain a correct thermal transport solution. For Saturn entries, the 2(sigma) - uncertainty for convective heating was less than 3.7%. The major uncertainty driver was dependent on shock temperature/velocity, changing from boundary layer thermal conductivity to diffusivity and then to shock layer ionization rate as velocity increases. While radiative heating for Uranus entry was negligible, the nominal solution for Saturn computed up to 20% radiative heating at the highest velocity examined. The radiative heating followed a non-normal distribution, with up to a 3x variation in magnitude. This uncertainty is driven by the H(sub 2) dissociation rate, as H(sub 2) that persists in the hot non-equilibrium zone contributes significantly to radiation.

Palmer, Grant↗

Quiet-time intensifications along the poleward auroral boundary near midnight

Radar and optical measurements from Sondrestrom are combined with satellite and Goose Bay data in a study of the poleward edge of the nightside auroral oval during a quiet period. The B(sub y) and B(sub z) components of the interplanetary magnetic field were close to zero, and the B(sub x) component was approximately 8 nT for more than 24 hours. On a large scale, the convection and precipitation patterns remained almost constant during this period; on a small scale, however, the conditions were quite dynamic. At 10- to 20-min intervals the arc that marked the poleward auroral boundary intensified, and a new arc appeared poleward of it. About once per hour, stronger intensifications were observed. One such event is examined in detail. The auroral arcs first appeared to dim, and then they brightened, with a factor of 10 increase in E region electron density. At the time of the brightening a new arc formed poleward of all the arcs. The arcs then drifted southward at velocities of approximately 270 m/s. A plasma drift disturbance, characterized by a doubling of the southward velocity and a reversal in the east-west component, propagated westward at 900 m/s through the fields of view of the Sondrestrom and Goose Bay radars. A simultaneous satellite overpass close to the radars revealed the presence of an energetic ion event similar to the 'velocity dispersed ion structures' observed on the Aureol satellite and presumed to be the signature of fast ion beams within the plasma sheet boundary layer. The stronger arc intensification events observed by the Sondrestrom radar are associated with an increase in plasma flow across the boundary between open and closed magnetic field lines. We interpret this increased flow as the ionospheric signature of abrupt, localized increases in the reconnection rate in the midnight sector.

De La Beaujardiere, O.↗

Modeling of Inner Magnetosphere Coupling Processes

The Ring Current (RC) is the biggest energy player in the inner magnetosphere. It is the source of free energy for Electromagnetic Ion Cyclotron (EMIC) wave excitation provided by a temperature anisotropy of RC ions, which develops naturally during inward E B convection from the plasmasheet. The cold plasmasphere, which is under the strong influence of the magnetospheric electric field, strongly mediates the RC-EMIC wave-particle-coupling process and ultimately becomes part of the particle and energy interplay. On the other hand, there is a strong influence of the RC on the inner magnetospheric electric and magnetic field configurations and these configurations, in turn, are important to RC dynamics. Therefore, one of the biggest needs for inner magnetospheric research is the continued progression toward a coupled, interconnected system with the inclusion of nonlinear feedback mechanisms between the plasma populations, the electric and magnetic fields, and plasma waves. As we clearly demonstrated in our studies, EMIC waves strongly interact with electrons and ions of energies ranging from approx.1 eV to approx.10 MeV, and that these waves strongly affect the dynamics of resonant RC ions, thermal electrons and ions, and the outer RB relativistic electrons. As we found, the rate of ion and electron scattering/heating in the Earth's magnetosphere is not only controlled by the wave intensity-spatial-temporal distribution but also strongly depends on the spectral distribution of the wave power. The latter is also a function of the plasmaspheric heavy ion content, and the plasma density and temperature distributions along the magnetic field lines. The above discussion places RC-EMIC wave coupling dynamics in context with inner magnetospheric coupling processes and, ultimately, relates RC studies with plasmaspheric and Superthermal Electrons formation processes as well as with outer RB physics.

Khazanov, George V.↗

Optimal gas cooling with minimal aero-opticaldistortion for next generation high-average-power lasers (Progress report, Year 2)

Next-generation lasers that combine high-peak-power ultrashort pulses with high repetition rates have the potential to efficiently deliver the high-average-power outputs required for inertial fusion energy generation, laser-based particle acceleration, and to generate compact secondary x- and 𝛾-ray radiation sources, among others. The combination of long gain length with high surface-to-volume ratio to remove residual heat makes gas-cooled multi-slab amplifier heads an inherently aperture- and average power scalable technology that enables this new class of solid-state lasers, and their thermal management becomes the most pressing technical challenge. As the index-of refraction depends on the gas density, heat removal from the slab surfaces unavoidably entails scattering of propagated light from regions of high temperature fluctuations. The requirement to maximize convective heat extraction from the gain medium under the constraint of minimum optical aberration of the beam renders amplifier slab cooling a technological challenge at the intersection of aero-optics, fluid-mechanics and heat transfer. The goal of this research is to leverage computational fluid dynamics and aero optics simulations of different levels of fidelity to devise a numerical modeling and optimization framework for gas-cooled multi-slab amplifiers. The project will be conducted in close collaboration with the Advanced Photon Technologies Group at Lawrence Livermore National Laboratory, where a gas-cooled test stand for laser materials and thermal management structures will be developed and installed. The final goal of the collaborative effort is to manufacture and benchmark a physical prototype of an numerically optimized amplifier head design. Measurement data on solid-state and thermo-optical distortions will be used for the validation of the numerical models. The key technical challenge in this multi-objective design problem is that of balancing optimal heat removal with minimal optical distortion. To facilitate computationally efficient design optimization and, at the same time, deepen our basic understanding of the interaction of light with turbulence, a complementary study that combines high and low-fidelity flow simulations and aero-optical models is proposed. Along the lines of these two objectives and solution strategies, the project is organized into two thrusts.

42 ENGINEERING↗

Interplanetary magnetic field control of mantle precipitation and associated field-aligned currents

Dayside reconnection, which is particularly effective for a southward interplanetary magnetic field (IMF), allows magnetosheath particles to enter the magnetosphere where they form the plasma mantle. The motions of the reconnected flux tube produce convective flows in the ionosphere. It is known that the convection patterns in the polar cap are skewed to the dawnside for a positive IMF B(sub y) (or duskside for a negative IMF B(sub y)) in the northern polar cap. Correspondingly, one would expect to find asymmetric distributions of mantle particle precipitation, but previous results have been unclear. In this paper the correlation between B(sub y) and the distribution of mantle particle precipitation is studied for steady IMF conditions with southward IMF. Ion and electron data from the Defense Meteorological Satellite Program (DMSP) F6 and F7 satellites are used to identify the mantle region and IMP 8 is used as a solar wind monitor to characterize the IMF. We study the local time extension of mantle precipitation in the prenoon and postnoon regions. We find that, in accordance with theoretical expectations for a positive (negative) IMF B(sub y), mantle particle precipitation mainly appears in the prenoon region of the northern (southern) hemisphere. The mantle particle precipitation can extend to as early as 0600 magnetic local time (MLT) in the prenoon region but extends over a smaller local time region in the postnoon sector (we did not find mantle plasma beyond 1600 MLT in our data set although coverage is scant in this area). Magnetometer data from F7 are used to determine whether part of the region 1 current flows on open field lines. We find that at times part of the region 1 sense current extends into the region of mantle particle precipitation, and is therefore on open field lines. In other cases, region 1 currents are absent on open field lines. Most of the observed features can be readily interpreted in terms of the open magnetosphere model.

Xu, Dingan↗

POLAR-UVI and other Coordinated Observations of a Traveling Convection Vortex Event Observed on 24 July 1996

Coordinated analysis of data from the POLAR UVI instrument, ground magnetometers, incoherent scatter radar, solar wind monitors IMP-8 and WIND, and DMSP satellite is focused on a traveling convection vortex (TCV) event on 24 July 1966. Starting at 10:48 UT, ground magnetometers in Greenland and eastern Canada measure pulsations consistent with the passing overhead of a series of alternating TCV field-aligned current pairs. Sondrestrom incoherent scatter radar measures strong modulation of the strength and direction of ionospheric plasma flow, The magnetometer pulsations grow in magnitude over the next hour, peaking in intensity at 11:39 UT, at which time the UVI instrument measures a localized intensification of auroral emissions over central and western Greenland. Subsequent images show the intensification grow in strength and propagate westward (tailward) until approximately 11:58 UT at which time the emissions fade. These observations are consistent with the westward passage of two pairs of moderately intense TCVs over central Greenland followed by a third very intense TCV pair. The intensification of auroral emissions at 11:39 UT is associated with the trailing vortex of the third TCV pair, thought to be the result of an upward field-aligned current. Measurements of the solar wind suggest that a pressure change may be responsible for triggering the first two pairs of TCVS, and that a subsequent sudden change in orientation of the IMF may have produced the intensification of the third TCV pair and the associated aurora] brightening. DMSP particle data indicate that the TCVs occur on field lines which map to the boundary plasma sheet or outer edge of the low latitude boundary layer.

Clauer, C. R.↗

Photospheric Current Spikes and Their Possible Association with Flares - Results from an HMI Data Driven Model

A data driven, near photospheric magnetohydrodynamic model predicts spikes in the horizontal current density, and associated resistive heating rate per unit volume Q. The spikes appear as increases by orders of magnitude above background values in neutral line regions (NLRs) of active regions (ARs). The largest spikes typically occur a few hours to a few days prior to M or X flares. The spikes correspond to large vertical derivatives of the horizontal magnetic field. The model takes as input the photospheric magnetic field observed by the Helioseismic & Magnetic Imager (HMI) on the Solar Dynamics Observatory (SDO) satellite. This 2.5 D field is used to determine an analytic expression for a 3 D magnetic field, from which the current density, vector potential, and electric field are computed in every AR pixel for 14 ARs. The field is not assumed to be force-free. The spurious 6, 12, and 24 hour Doppler periods due to SDO orbital motion are filtered out of the time series of the HMI magnetic field for each pixel using a band pass filter. The subset of spikes analyzed at the pixel level are found to occur on HMI and granulation scales of 1 arcsec and 12 minutes. Spikes are found in ARs with and without M or X flares, and outside as well as inside NLRs, but the largest spikes are localized in the NLRs of ARs with M or X flares. The energy to drive the heating associated with the largest current spikes comes from bulk flow kinetic energy, not the electromagnetic field, and the current density is highly non-force free. The results suggest that, in combination with the model, HMI is revealing strong, convection driven, non-force free heating events on granulation scales, and that it is plausible these events are correlated with subsequent M or X flares. More and longer time series need to be analyzed to determine if such a correlation exists. Above an AR dependent threshold value of Q, the number of events N(Q) with heating rates greater than or equal to Q obeys a scale invariant power law distribution for each AR given by N(Q) varies Q(sup -s), where 0.40 less than or equal to S less than or equal to 0.53, with a mean and standard deviation across the 14 ARs of 0.47 and 0.045, showing there is little variation of S from one AR to another. These properties of N(Q) are in close agreement with those of the distribution N(E) for the total energy E of solar flares, determined from observations to be N(E) = constant x E(sup -alpha). From observations of nanoflares in the 0.7 to 4 MK range, and from observations of flares in hard X-rays, it is found that 0.51 less than or equal to alpha less than or equal to 0.57, and 0.4 less than or equal to alpha less than or equal to 0.6, respectively (Crosby et al. 1993, Sol. Phys., 143, 275; Aschwanden & Parnell 2002, ApJ, 572, 1048). Observations also show that, as is found here for the exponent S, there is little variation of alpha with AR (Wheatland 2000, ApJ, 532, 1209), indicating N(E) and N(Q) are largely independent of individual properties of ARs such as area, total magnetic flux, and distribution of current density (i.e. non-potentiality). Therefore the power law scaling of the photospheric heating rate Q computed here on granulation scales is essentially identical to that found for coronal observations of flare energies on scales 1-2 orders of magnitude larger. This suggests the physical mechanisms that cause Q and coronal flares are closely related. It seems likely that Q is the signature of a magnetic reconnection process in an energy range and volume orders of magnitude smaller than those of flares. In this context, at least the larger spikes in Q might be signatures of UV photospheric or lower chromospheric bombs in which plasma is heated to temperatures approximately 10(exp -5) K (Peter et al. 2014, Science 346, 1255726; Judge 2015, ApJ, 808, 116). In addition, lattice based avalanche simulations of flare energy release predict 0.4 less than or equal to alpha less than or equal to 0.5, while analytic, fractal-diffusive self-organized criticality models predict 0.4 less than or equal to alpha less than or equal to 0.67, in excellent agreement with observations, and the results presented here (Aschwanden & Parnell 2002, ApJ, 572, 1048; Aschwanden 2012, A&A, 539, A2; Aschwanden 2013, in "Self Organized Criticality Systems"; Aschwanden et al. 2016, SSR, 198, 47).

Goodman, Michael↗

Experimental Investigation of Average Heat-Transfer and Friction Coefficients for Air Flowing in Circular Tubes Having Square-Thread-Type Roughness

An investigation of forced-convection heat transfer and associated pressure drops was conducted with air flowing through electrically heated Inconel tubes having various degrees of square-thread-type roughness, an inside diameter of 1/2 inch, and a length of 24 inches. were obtained for tubes having conventional roughness ratios (height of thread/radius of tube) of 0 (smooth tube), 0.016, 0.025, and 0.037 over ranges of bulk Reynolds numbers up to 350,000, average inside-tube-wall temperatures up to 1950deg R, and heat-flux densities up to 115,000 Btu per hour per square foot. Data The experimental data showed that both heat transfer and friction increased with increase in surface roughness, becoming more pronounced with increase in Reynolds number; for a given roughness, both heat transfer and friction were also influenced by the tube wall-to-bulk temperature ratio. Good correlation of the heat-transfer data for all the tubes investigated was obtained by use of a modification of the conventional Nusselt correlation parameters wherein the mass velocity in the Reynolds number was replaced by the product of air density evaluated at the average film temperature and the so-called friction velocity; in addition, the physical properties of air were evaluated at the average film temperature. The isothermal friction data for the rough tubes, when plotted in the conventional manner, resulted in curves similar to those obtained by other investigators; that is, the curve for a given roughness breaks away from the Blasius line (representing turbulent flow in smooth tubes) at some value of Reynolds number, which decreases with increase in surface roughness, and then becomes a horizontal line (friction coefficient independent of Reynolds number). A comparison of the friction data for the rough tubes used herein indicated that the conventional roughness ratio is not an adequate measure of relative roughness for tubes having a square-thread-type element. The present data, as well as those of other investigators, were used to isolate the influence of ratios of thread height to width, thread spacing to width, and the conventional roughness ratio on the friction coefficient. A fair correlation of the friction data was obtained for each tube with heat addition when the friction coefficient and Reynolds number were defined on the basis of film properties; however, the data for each tube retained the curve characteristic of that particular roughness. The friction data for all the rough tubes could be represented by a single line for the complete turbulence region by incorporating a roughness parameter in the film correlation. No correlation was obtained for the region of incomplete turbulence.

Sams, E. W.↗

Polar Spacecraft Based Comparisons of Intense Electric Fields and Poynting Flux Near and Within the Plasma Sheet-Tail Lobe Boundary to UVI Images: An Energy Source for the Aurora

In this paper, we present measurements from two passes of the Polar spacecraft of intense electric and magnetic field structures associated with Alfven waves at and within the outer boundary of the plasma sheet at geocentric distances of 4-6 R(sub E), near local midnight. The electric field variations have maximum values exceeding 100 mV/m and are typically polarized approximately normal to the plasma sheet boundary. The electric field structures investigated vary over timescales (in the spacecraft frame.) ranging front 1 to 30 s. They are associated with strong magnetic field fluctuations with amplitudes of 10-40 nT which lie predominantly ill the plane of the plasma sheet and are perpendicular to the local magnetic field. The Poynting flux associated with the perturbation fields measured at these altitudes is about 1-2 ergs per square centimeters per second and is directed along the average magnetic field direction toward the ionosphere. If the measured Poynting flux is mapped to ionospheric altitudes along converging magnetic field lines. the resulting energy flux ranges up to 100 ergs per centimeter squared per second. These strongly enhanced Poynting fluxes appear to occur in layers which are observed when the spacecraft is magnetically conjugate (to within a 1 degree mapping accuracy) to intense auroral structures as detected by the Polar UV Imager (UVI). The electron energy flux (averaged over a spatial resolution of 0.5 degrees) deposited in the ionosphere due to auroral electron beams as estimated from the intensity in the UVI Lyman-Birge-Hopfield-long filters is 15-30 ergs per centimeter squared per second. Thus there is evidence that these electric field structures provide sufficient Poynting flux to power the acceleration of auroral electrons (as well as the energization of upflowing ions and Joule heating of the ionosphere). During some events the phasing and ratio of the transverse electric and magnetic field variations are consistent with earthward propagation of Alfven surface waves with phase velocities of 4000-10000 kilometers per second. During other events the phase shifts between electric and magnetic fields suggest interference between upward and downward propagating Alfven waves. The E/B ratios are about an order of magnitude larger than typical values of C/SIGMA(sub p), where SIGMA(sub p), is the height integrated Pedersen conductivity. The contribution to the total energy flux at these altitudes from Poynting flux associated with Alfven waves is comparable to or larger than the contribution from the particle energy flux and 1-2 orders of magnitude larger than that estimated from the large-scale steady state convection electric field and field-aligned current system.

Wygant, J. R.↗

A double-moment multiple-phase four-class bulk ice scheme, Part 1: Description

A detailed ice-phase bulk microphysical scheme has been developed for simulating the hydrometeor distributions of convective and stratiform precipitation in different large-scale environmental conditions. The proposed scheme involves 90 distinct microphysical processes, which predict the mixing ratios and the number concentrations of small ice crystals, snow, graupel, and frozen drops/hail, as well as the mixing ratios of liquid water on wet precipitation ice (snow, graupel, frozen drops). The number of adjustable coefficients has been significantly reduced in comparison with other bulk schemes. Additional improvements have been made to the parameterization in the following areas: (1) representing small ice crystals with nonzero terminal fall velocities and dispersive size distributions, (2) accurate and computationally efficient calculations of precipitation collection processes, (3) reformulating the collection equation to prevent unrealistically large accretion rates, (4) more realistic conversion by riming between different classes of precipitation ice, (5) preventing unrealistically large rates of raindrop freezing and freezing of liquid water on ice, (6) detailed treatment of various rime-splintering ice multiplication mechanisms, (7) a simple representation of the Hobbs-Rangno ice enhancement process, (8) aggregation of small ice crystals and snow, and (9) allowing explicit competition between cloud water condensation and ice deposition rates rather than using saturation adjustment techniques. For the purposes of conserving the higher moments of the particle distributions, preserving the spectral widths (or slopes) of the particle spectra is shown to be more important than strict conservation of particle number concentration when parameterizing changes in ice-particle number concentrations due to melting, vapor transfer processes (sublimation of dry ice, evaporation from wet ice), and conversion between different hydrometeor species. The microphysical scheme is incorporated into a nonhydrostatic cloud model in Part 2 of this study. The model performed well in simulating the radar and microphysical structures of a midlatitude-continental squall lines and a tropical-maritime squall system with minimal tuning of the parameterization, even though the vertical profiles of radar reflectivity differed substantially between these storms.

Ferrier, Brad Shoenberg↗

Large Eddy Simulation of Cirrus Clouds

The Regional Atmospheric Modeling System (RAMS) with mesoscale interactive nested-grids and a Large-Eddy Simulation (LES) version of RAMS, coupled to two-moment microphysics and a new two-stream radiative code were used to investigate the dynamic, microphysical, and radiative aspects of the November 26, 1991 cirrus event. Wu (1998) describes the results of that research in full detail and is enclosed as Appendix 1. The mesoscale nested grid simulation successfully reproduced the large scale circulation as compared to the Mesoscale Analysis and Prediction System's (MAPS) analyses and other observations. Three cloud bands which match nicely to the three cloud lines identified in an observational study (Mace et al., 1995) are predicted on Grid #2 of the nested grids, even though the mesoscale simulation predicts a larger west-east cloud width than what was observed. Large-eddy simulations (LES) were performed to study the dynamical, microphysical, and radiative processes in the 26 November 1991 FIRE 11 cirrus event. The LES model is based on the RAMS version 3b developed at Colorado State University. It includes a new radiation scheme developed by Harrington (1997) and a new subgrid scale model developed by Kosovic (1996). The LES model simulated a single cloud layer for Case 1 and a two-layer cloud structure for Case 2. The simulations demonstrated that latent heat release can play a significant role in the formation and development of cirrus clouds. For the thin cirrus in Case 1, the latent heat release was insufficient for the cirrus clouds to become positively buoyant. However, in some special cases such as Case 2, positively buoyant cells can be embedded within the cirrus layers. These cells were so active that the rising updraft induced its own pressure perturbations that affected the cloud evolution. Vertical profiles of the total radiative and latent heating rates indicated that for well developed, deep, and active cirrus clouds, radiative cooling and latent heating could be comparable in magnitude in the cloudy layer. This implies that latent heating cannot be neglected in the construction of a cirrus cloud model. The probability density function (PDF) of w was analyzed to assist in the parameterization of cloud-scale velocities in large-scale models. For the more radiatively-driven, thin cirrus case, the PDFs are approximately Gaussian. However, in the interior of the deep, convectively unstable case, the PDFs of w are multi-modal and very broad, indicating that parameterizing cloud-scale motions for such clouds can be very challenging. The results of this research are described in detail in a paper submitted to the Journal of Atmospheric Science (Wu and Cotton, 1999), which is enclosed as Appendix 2. Using soundings extracted from a mesoscale simulation of the November 26, 1991 cirrus event, the radiative effects on vapor deposition/sublimation of ice crystals was studied using a two-dimensional cloud-resolving model (CRM) version of RAMS, coupled to an explicit bin-resolving microphysics. The CRM simulations of the November 26, 1991 cirrus event demonstrate that the radiative impact on the diffusional growth (or sublimation) of ice crystals is significant. In this case, the ice particles experienced radiative warming. Model results show that radiative feedbacks in the diffusional growth of ice particles can be very complex. Radiative warming of an ice particle will restrict the particle's diffusional growth. In the case of radiative warming, ice particles larger than a certain size will experience so much radiative warming that surface ice saturation vapor pressures become large enough to cause sublimation of the larger crystals, while smaller crystals are growing by vapor deposition. However, ice mass production can be enhanced in the case of radiative cooling of an ice particle. For the November 26, 1991 cirrus event, radiative feedback results in significant reduction in the total ice mass, especially in the production of large ice crystals, and consequently, both radiative and dynamic properties of the cirrus cloud are significantly affected. A complete description of this research has been submitted as a paper to the Journal of Atmospheric Science (Wu et al., 1999), and included as Appendix 3.

Wu, Ting↗

Multi-Instrument Analysis of a Traveling Convection Vortex Event on July 24, 1996 Coordinated with the Polar UVI

We present the analysis of a coordinated set of observations from the POLAR Ultraviolet Imager (UVI), ground magnetometers, incoherent scatter radar, solar wind monitors, DMSP and GOES satellites, focused on a traveling convection vortex (TCV) event on 24th July 1996. Starting at approximately 10:48 UT, around magnetometers in Greenland and northern Canada observe pulsations consistent with the passing overhead of a series of alternating TCV filed-aligned current pairs. Azimuthal scans by the Sondrestrom incoherent scatter radar located near Kangerlussuaq (formerly Sondrestrom), Greenland, at this time show strong modulation in the strength and direction of ionospheric plasma flow. The magnetometer pulsations grow in magnitude over the next hour, peaking in intensity at 11:39 UT, at which time images form the UVI instrument show a localized intensification of auroral emissions over central and western Greenland. Subsequent images show the intensification grow in strength and propagate westward (tailward) until approximately 11:58 UT at which time the intensification fades. These observations are consistent with the westward passage of two pairs of moderately intense TCVs over central Greenland followed by a third very intense TCV pair. The intensification of auroral emissions at 11:39 UT is associated with the trailing vortex of the third TCV pair, thought to be the result of an upward field aligned current. The modulated flow observed by the radar is the result of the strong electric fields associated with the impulsive TCV related field aligned current systems as they pass through the field of view of the radar. Measurements of the solar wind from the V;IND and IMP-8 spacecraft suggest that a pressure change may be responsible for triggering the first two pairs of TCVS, and that a subsequent sudden change in the orientation of the interplanetary magnetic field may have produced the intensification of the third TCV pair and the associated auroral brightening. Magnetometer data from the GOES satellite located over the eastern United States at geostationary orbit is consistent with a series of field-aligned moving tailward past the satellite. DMSP particle data indicated that the TCVs occur on field lines which map to the boundary plasma sheet (BPS).

Sitar, R. J.↗

Coherence Analysis of the Space Launch System using Unsteady Pressure Sensitive Paint

Transonic buffet during atmospheric ascent is a major source of unsteady loading on launch vehicles that, in the past, has led to structural failures. Thus, determining accurate buffet forcing functions (BFFs) to properly predict the vehicle response to buffet is of vital importance in launch vehicle design. The state of the art for obtaining the BFFs relies on wind-tunnel tests where the fluctuating pressures are measured by pressure transducers (PTs) at sparse locations on a rigid, geometrically-scaled buffet model. To compute the BFFs, the outer mold line (OML) of the vehicle is mapped onto contiguous panels centered at the location of the PTs and the measured fluctuating pressures are integrated over the panels’ areas. To mitigate conservatism due to the assumption that the measured buffet pressures act in phase across each panel, coherence factors are applied, effectively reducing the integration areas and, therefore, the buffet forces. For coherence factors to provide the proper level of BFF attenuation, accurate knowledge of the spatial and temporal correlation of the buffet pressures is paramount. Unfortunately, even with hundreds of unsteady pressure transducers instrumenting the models, compromises must be made on the spatial resolution of the PTs. Typically, the PT layout aims at resolving the pressure correlation along the longitudinal axis of the vehicle, assuming that coherent structures propagate mostly in the longitudinal direction. As a result, the distribution of azimuthal/cross-stream correlation is not well known and its impact on the estimated BFFs is often neglected. To fill this and other knowledge gaps, extremely high spatial-resolution uPSP data were collected for three different production-design configurations of the Space Launch System in the NASA Ames Research Center Unitary Plan Wind Tunnel 11-Foot Transonic Wind Tunnel. Specifically, two cargo configurations, the Block 1 and Block 1B, and one crew configuration, the Block 1B crew, were surveyed at resolutions ranging from 600,000 to over 1 million uPSP measurement locations. To shed light on the temporal and spectral behavior of are presented. Several OML regions and flow features of interest are investigated, from theexpansion/shock on the Orion Multi-Purpose Crew Vehicle, to the terminal shock environment on the core stage, and the Strouhal shedding behind the boosters forward attach. The sensitivity of these environments to the vehicle attitude is examined. Furthermore, for selected panels,coherence factors that accurately capture the azimuthal coherence distribution of the buffetpressures are derived and their impact on the estimated BFFs is discussed. Finally, distributions of the local convection velocity and cross spectrum phase are presented.

buffet↗

Origin, Bulk Chemical Composition and Physical Structure of the Galilean Satellites of Jupiter: A Post-Galileo Analysis

The origin of Jupiter and the Galilean satellite system is examined in the light of the new data that has been obtained by the NASA Galileo Project. In particular, special attention is given to a theory of satellite origin which was put forward at the start of the Galileo Mission and on the basis of which several predictions have now been proven successful. These predictions concern the chemical composition of Jupiter's atmosphere and the physical structure of the satellites. According to the proposed theory of satellite origin, each of the Galilean satellites formed by chemical condensation and gravitational accumulation of solid grains within a concentric family of orbiting gas rings. These rings were cast off equatorially by the rotating proto-Jovian cloud (PJC which contracted gravitationally to form Jupiter some 4 1/2 billion years ago. The PJC formed from the gas and grains left over from the gas ring that had been shed at Jupiter's orbit by the contracting proto-solar cloud (PSC Supersonic turbulent convection provides the means for shedding discrete gas rings. The temperatures T (sub n) of the system of gas rings shed by the PSC and PJC vary with their respective mean orbital radii R (sub n) (n = 0, 1,2,...) according as T (sub n) proportional to R (sub n) (exp -0.9). If the planet Mercury condenses at 1640 K, so accounting for the high density of that planet via a process of chemical fractionation between iron and silicates, then T (sub n) at Jupiter's orbit is 158 K. Only 35% of the water vapour condenses out. Thus fractionation between rock and ice, together with an enhancement in the abundance of solids relative to gas which takes place through gravitational sedimentation of solids onto the mean orbit of the gas ring, ensures nearly equal proportions of rock and ice in each of Ganymede and Callisto. Io and Europa condense above the H20 ice point and consist solely of hydrated rock (h-rock). The Ganymedan condensate consists of h-rock and H20 ice. For Callisto, NH3 ice makes up -5% of the condensate mass next to h-rock (approximately 50%) and H2O ice (approximately 45%). Detailed thermal and structural models for each of Europa, Ganymede and Callisto are constructed on the basis of the above initial bulk chemical compositions. For Europa (E), a predicted 2-zone model consisting of a dehydrated rock core of mass 0.912 M (sub E) and a 150 km thick frozen mantle of salty H2O yields a moment-of-inertia coefficient which matches the Galileo Orbiter gravity measurement. For Ganymede (G), a 3-zone model possessing an inner core of solid FeS and mass approximately 0.116 M (sub G), and an outer H2O ice mantle of mass approximately 0.502 M (sub G) is needed to explain the gravity data. Ganymede's native magnetic field was formed by thermoremanent magnetization of Fe3O4. A new Callisto (C) model is proposed consisting of a core of mass 0.826 M (sub C) containing a uniform mixture of h-rock (60% by mass) and H2O and NH3 ices, and capped by a mantle of pure ice. This model may have the capacity to yield a thin layer of liquid NH3 (raised dot) 2H2O at the core boundary, in line with Galileo's discovery of an induced magnetic field.

Prentice, A. J. R.↗

Rewet Temperature Correlations for Liquid-Nitrogen Boiling Pipe Flows Across Varying Flow Conditions and Orientations

In many convective liquid-vapor phase change heat transfer engineering applications, cryogenic fluids are widely used in industrial processes, spacecraft and cryosurgery systems, and so on. For example, cryogens are usually used as liquid fuels such as liquid hydrogen and oxygen in the rocket industry, liquid nitrogen (LN2) and helium are frequently used to cool superconducting magnetic device for medical applications. In these systems, proper transport, handling, and storage of cryogenic fluids are of extreme importance. Among all the cryogenic transport processes performed in room temperatures, quenching, also termed chilldown, is a unavoidable initial, transient phase-change heat transfer process that brings the system down to the cryogenic condition. The Leidenfrost temperature or rewet temperature that signals the end of film boiling is practically considered the completion point of a quenching process. Therefore, rewet temperature has been considered the most important parameter for the engineering design of cryogenic thermal management systems. As most of the previous correlations for predicting the Leidenfrost temperature and the rewet temperature have been basically developed for water, they are shown to disagree with recent liquid nitrogen pipe chilldown experiments in upward and downward flow directions over a wide range of flow rates, pressures, and degrees of inlet subcooling. In addition to a complete review of the literature, two new correlations are presented in this work, one based on bubble growth and another based on the theoretical maximum limit of superheat. Each correlation performs well over the entire data set.

liquid nitrogen↗