Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “strong coupling”

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 199 records · Page 11

The deformation of flux tubes in the solar wind with applications to the structure of magnetic clouds and CMEs

Two dimensional magnetohydrodynamic simulations of the distortion of a magnetic flux tube, accelerated through ambient solar wind plasma, are presented. Vortices form on the trailing edge of the flux tube, and couple strongly to its interior. If the flux tube azimuthal field is weak, it deforms into an elongated banana-like shape after a few Alfven transit times. A significant azimuthal field component inhibits this distortion. In the case of magnetic clouds in the solar wind, it is suggested that the shape observed at 1 AU was determined by distortion of the cloud in the inner heliosphere. Distortion of the cloud beyond 1 AU takes many days. It is estimated that effective drag coefficients slightly greater than unity are appropriate for modeling flux tube propagation. Synthetic magnetic field profiles as would be seen by a spacecraft traversing the cloud are presented.

Cargill, Peter J.↗

Air pollutant transport in a coastal environment. Part 1: Two-dimensional simulations of sea-breeze and mountain effects

Over the southern California coastal region, observations of the vertical distributions of pollutants show that maximum concentrations can occur within temperature inversion layers well above the surface. A mesoscale model is used to study the dynamical phenomena that cause such layers, including sea breezes and mountain flows, and to study the characteristics of air pollutant transport in a coastal environment capped by a temperature inversion. The mathematical and physical structure of the model is described. Two-dimensional simulations corresponding to four configurations of coastal plains and mountains are discussed. The simulations reveal that pollutant transport over a coastal plain is strongly influenced by the topographic configuration, including the height of coastal mountains and their distance from the coastline. Sea breezes induced by land-sea thermal contrasts, as well as upslope winds induced along mountain flanks, both create vertical transport that can lead to the formation of elevated pollution layers. The sea-breeze circulation generates pollution layers by undercutting the mixed layer and lofting pollutants into the stable layer. Heating of mountain slopes acts to vent pollutants above the mountain ridge during the day; during the evening, pollutants can be injected directly into the inversion layer from the decaying upslope flows. In a land-sea configuration with mountains close to the coastline, the sea breeze and heated-mountain flow are strongly coupled. In the afternoon, this interaction can produce upslope flow from which polluted air is detrained into the inversion layer as a return circulation. When the mountains lie farther inland, however, pollutants may be trapped aloft when the mixed layer stabilizes in the late afternoon. As the nocturnal boundary layer forms over the coast in the evening, polluted mixed-layer air is effectively left behind in the inversion layer. In the Los Angeles Basin, the formation mechanism for elevated polluted layers is most similar to our cases with inland mountains.

Lu, Rong↗

Observations of cloud liquid water path over oceans: Optical and microwave remote sensing methods

Published estimates of cloud liquid water path (LWP) from satellite-measured microwave radiation show little agreement, even about the relative magnitudes of LWP in the tropics and midlatitudes. To understand these differences and to obtain more reliable estimate, optical and microwave LWP retrieval methods are compared using the International Satellite Cloud Climatology Project (ISCCP) and special sensor microwave/imager (SSM/I) data. Errors in microwave LWP retrieval associated with uncertainties in surface, atmosphere, and cloud properties are assessed. Sea surface temperature may not produce great LWP errors, if accurate contemporaneous measurements are used in the retrieval. An uncertainty of estimated near-surface wind speed as high as 2 m/s produces uncertainty in LWP of about 5 mg/sq cm. Cloud liquid water temperature has only a small effect on LWP retrievals (rms errors less than 2 mg/sq cm), if errors in the temperature are less than 5 C; however, such errors can produce spurious variations of LWP with latitude and season. Errors in atmospheric column water vapor (CWV) are strongly coupled with errors in LWP (for some retrieval methods) causing errors as large as 30 mg/sq cm. Because microwave radiation is much less sensitive to clouds with small LWP (less than 7 mg/sq cm) than visible wavelength radiation, the microwave results are very sensitive to the process used to separate clear and cloudy conditions. Different cloud detection sensitivities in different microwave retrieval methods bias estimated LWP values. Comparing ISCCP and SSM/I LWPs, we find that the two estimated values are consistent in global, zonal, and regional means for warm, nonprecipitating clouds, which have average LWP values of about 5 mg/sq cm and occur much more frequently than precipitating clouds. Ice water path (IWP) can be roughly estimated from the differences between ISCCP total water path and SSM/I LWP for cold, nonprecipitating clouds. IWP in the winter hemisphere is about 3 times the LWP but only half the LWP in the summer hemisphere. Precipitating clouds contribute significantly to monthly, zonal mean LWP values determined from microwave, especially in the intertropical convergence zone (ITCZ), because they have almost 10 times the liquid water (cloud plus precipitation) of nonprecipitating clouds on average. There are significant differences among microwave LWP estimates associated with the treatment of precipitating clouds.

Lin, Bing↗

Relationships between evaprorative fraction and remotely sensed vegetation index and microwave brightness temperature for semiarid rangelands

Measurements of the microwave brightness temperature (TB) with the Pushbroom Microwave Radiometer (PBMR) over the Walnut Gulch Experiment Watershed were made on selected days during the MONSOON 90 field campaign. The PBMR is an L-band instrument (21-cm wavelength) that can provide estimates of near-surface soil moisture over a variety of surfaces. Aircraft observations in the visible and near-infrared wavelengths collected on selected days also were used to compute a vegetation index. Continuous micrometeorological measurements and daily soil moisture samples were obtained at eight locations during experimental period. Two sites were instrumented with time domain reflectometry probes to monitor the soil moisture profile. The fraction of available energy used for evapotranspiration was computed by taking the ratio of latent heat flux (LE) to the sum of net radiation (Rn) and soil heat flux (G). This ratio is commonly called the evaporative fraction (EF) and normally varies between 0 and 1 under daytime convective conditions with minimal advection. A wide range of environmental conditions existed during the field campaign, resulting in average EF values for the study area varying from 0.4 to 0.8 and values of TB ranging from 220 to 280 K. Comparison between measured TB and EF for the eight locations showed an inverse relationship. Other days were included in the analysis by estimating TB with the soil moisture data. Because transpiration from the vegetation is more strongly coupled to root zone soil moisture, significant scatter in this relationship existed at high values of TB or dry near-surface soil moisture conditions. The variation in EF under dry near-surface soil moisture conditions was correlated to the amount of vegetation cover estimated with a remotely sensed vegetation index. These findings indicate that information obtained from optical and microwave data can be used for quantifying the energy balance of semiarid areas. The microwave data can indicate when soil evaporation is significantly contributing to EF, while the optical data is helpful for quantifying the spatial variation in EF due to the distribution of vegetation cover.

Kustas, W. P.↗

Critical review of electrical conductivity measurements and charge distribution analysis of magnesium oxide

The electrical conductivity sigma of MgO single crystals shows a sharp increase at 500-800 C, in particular of sigma surface, generally attributed to surface contamination. Charge Distribution Analysis (CDA), a new technique providing information on fundamental properties that was previously unavailable, allows for the determination of surface charges, their sign and associated internal electric field. Data on 99.99% purity, arc-fusion grown MgO crystals show that mobile charge carriers start to appear in the bulk of the MgO crystals between 200 and 400 C when sigma (measured by conventional techniques) is in t he 10(exp -14) to 10(exp -16) /omega/cm range. Above 500 C, as sigma increases to 10(exp -6) to 10(exp -7)/omega/cm, more charges appear giving rise to a strong positive surface charge supported by a strong internal field. This indicates that charges are generated in the bulk and diffuse to the surface by an internally controlled process. On the basis of their positive sign they are identified as holes (defect electrons). Because of the low cation content of these very pure MgO crystals, theses holes cannnot be associated with transition metal impurties. Instead, they are associated with the O(2-) sublattice, e.g. consist of O(-) states or positive holes. This conclusion is supported by magnetic susceptibility data showing the appearance of 1000 +/- 500 ppm paramagnetic species between 200-500 C. The magnetic data are consistent with strongly coupled, diamagnetic O(-) pairs below 200-500 C, chemically equivalent to peroxy anions, O2(2-), and probably associated with cation vacancies in the MgO matrix. The formation of O2(2-) in arc-fusion grown MgO crystals is very unexpected because of the highly reducing growth conditions. Their presence implies an internal redox reaction involving dissolved 'water' by which OH(-) pairs convert to O2(2-) plus H2 molecules. This redox conversion is supported by mass spectroscopic measurements of the H2 release from highly OH(-)-doped, finely divided MgO and by wet-chemical analysis of its oxidant concentration.

Freund, Friedemann↗

Geometry of interplanetary magnetic clouds

Two dimensional magnetohydrodynamic simulations are presented of the distortion of a magnetic flux rope that is being accelerated through ambient solar wind plasma. The flux rope magnetic field has an axial component parallel to the solar wind field and an azimuthal component, which lies in the simulation plane. As the flux rope moves through the solar wind plasma, vortices form on its trailing edge and couple strongly to its interior. If the flux rope azimuthal field is weak, it deforms into an elongated banana-like shape a few Alfven transit times. A strong azimuthal field component tends to inhibit this distortion. If the flux rope is taken to model a magnetic cloud, it is suggested that the shape of the cloud at 1 AU is determined by its distortion in the inner solar wind. Distortion timescales beyond 1 AU are estimated as many days. It is estimated that effective drag coefficients somewhat greater than unity are appropriate for modelling flux rope propagation.

Cargill, P. J.↗

Seasonal variations in the atmospheric distribution of a reactive chlorine compound, tetrachloroethene (CCl2 = CCl2)

Tropospheric mixing ratios of CCl2 = CCl2 were measured at remote surface locations in the Pacific between 71 deg N and 47 deg S during September and December of 1989, and March and June of 1990. The observed gradient of decreasing concentrations from the northern to the southern hemisphere, and very low concentrations in the southern hemisphere throughout the year, indicates a predominant input from the northern hemisphere. Seasonal measurements in the northern hemisphere showed maximum CCl2 = CCl2 concentrations occurring in the late summer. This distinct seasonal variation is strongly coupled to the atmompheric abundance of hydroxyl radical, the only important species responsible for CCl2 = CCl2 removal. With the estimated global CCl2 =CCl2 emissions, the lifetime is calculated to be about 5.4 months which is in good agreement with the 4.0 month estimate obtained from the inverse ratio of its measured hydroxyl reaction rate constant compared with that of methyl-chloroform CH3CCl3.

Wang, Charles J.-L.↗

Electronic properties of Y-Ba-Cu-O superconductors as seen by Cu and O NMR/NQR

Nuclear Magnetic Resonance (NMR) and Nuclear Quadrupole Resonance (NQR) allow the investigation of electronic properties at the atomic level. We will report on such studies of typical members of the the Y-Ba-Cu-O family such as YBa2Cu30(6 + x) (1-2-3-(6 + x)), YBa2Cu4O8 (1-2-4) and Y2Ba4Cu7015 (2-4-7) with many examples of recent work performed in our laboratory. In particular, we will deal with Knight shift and relaxation studies of copper and oxygen. We will discuss important issues of current studies such as: (1) Existence of a common electronic spin-susceptibility in the planes (and perhaps in the chains) of 1-2-4; (2) Strong evidence for the existence of a pseudo spin-gap of the antiferromagnetic fluctuations in 1-2-4 and 2-4-7; (3) Evidence for d-wave pairing in 1-2-4; (4) Strong coupling of inequivalent Cu-O planes in 2-4-7 and possible origin for the high Tc value of this compound; and (5) The possibility to describe NMR data in the framework of a charge-excitation picture.

Brinkmann, D.↗

Quaternary borocarbides: New class of intermetallic superconductors

Our recent discovery of superconductivity (SC) in the four-element multiphase Y-Ni-B-C system at an elevated temperature (TC approximately 12 K) has opened up great possibilities of identifying new superconducting materials and generating new physics. Superconductivity with Tc (greater than 20 K) higher than that known so far in bulk intermetallics has been observed in multiphase Y-Pd-B-C and Th-Pd-B-C systems and a family of single phase materials RENi2B2C (RE= Y, rare earth) have been found. Our investigations show YNi2B2C to be a strong coupling hard type-II SC. HC2(T) exhibits an unconventional temperature dependence. Specific heat and magnetization studies reveal coexistence of SC and magnetism in RNi2B2C (R = Ho, Er, Tm) with magnetic ordering temperatures (Tc approximately 8 K, 10.5 K, 11 K and Tm approximately 5 K, approximately 7K, approximately 4 K respectively) that are remarkably higher than those in known magnetic superconductors . Mu-SR studies suggest the possibility of Ni atoms carrying a moment in TmNi2B2C. Resistivity results suggests a double re-entrant transition (SC-normal-SC) in HoNi2B2C. RENi2B2C (RE = Ce, Nd, Gd) do not show SC down to 4.2 K. The Nd- and Gd-compounds order magnetically at approximately 4.5 K and approximately 19.5 K, respectively. Two SC transitions are observed in Y-Pd-B-C (Tc approximately 22 K, approximately 10 K) and in Th-Pd-B-C (Tc approximately 20 K, approximately 14 K) systems, which indicate that there are at least two structures which support SC in these borocarbides. In our multiphase ThNi2B2C we observe SC at approximately 6 K. No SC was seen in multiphase UNi2B2C, UPd2B2C, UOs2Ge2C and UPd5B3C(0.35) down to 4.2 K. Tc in YNi2B2C is depressed by substitutions (Gd, Th and U at Y-sites and Fe, Co at Ni-sites).

Nagarajan, R.↗

Structural, dynamical and electronic properties of CaCuO2

The scalar relativistic version of an accurate first principles full potential self-consistent linearized muffin tin orbital (LMTO) method has been employed for describing the physical properties of the parent system of the high-T(sub c) oxide superconductors, i.e., CaCuO2. The presently employed modified version of the LMTO method is quite fast and goes beyond the usual LMTO ASA method in the sense that it permits a completely general shape of the potential and the charge density. Also, in contrast to LMTO ASA, the present method is also capable of treating distorted lattice structures accurately. The calculated values of the lattice parameters of pure CaCuO2 lie within 3% of the experimentally measured values for the Sr-doped system Ca(0.86)Sr(0.14)CuO(2). The computed electronic structures and the density of states is quite similar to those of the other oxide superconductors, except of their three- dimensional character because of the presence of strong coupling between the closely spaced CuO2 layers. The van Hove singularity peak appears slightly below the Fermi level and a small concentration of oxygenation /or/ substitutional doping may pin it at the Fermi level. The calculated frequencies for some symmetric frozen phonons for undoped CaCuO2 are quite near to the measured data for the Sr-doped CaCuO2.

Agrawal, Bal K.↗

On the Nonlinear Evolution of a Stationary Cross-Flow Vortex in a Fully Three-Dimensional Boundary Layer Flow

We consider the nonlinear stability of a fully three-dimensional boundary layer flow in an incompressible fluid and derive an equation governing the nonlinear development of a stationary cross-flow vortex. The amplitude equation is a novel integro-differential equation which has spatial derivatives of the amplitude occurring in the kernal function. It is shown that the evolution of the cross-flow vortex is strongly coupled to the properties of an unsteady wall layer which is in fact driven by an unknown slip velocity, proportional to the amplitude of the cross-flow vortex. The work is extended to obtain the corresponding equation for rotating disk flow. A number of special cases are examined and the numerical solution for one of cases, and further analysis, demonstrates the existence of finite-distance as well as focussing type singularities. The numerical solutions also indicate the presence of a new type of nonlinear wave solution for a certain set of parameter values.

Gajjar, J. S. B.↗

Two Carrier Analysis of Persistent Photoconductivity in Modulation-Doped Structures

A simultaneous fit of Hall and conductivity data gives quantitative results on the carrier concentration and mobility in both the quantum well and the parallel conduction channel. In this study this method was applied to reveal several new findings on the effect of persistent photoconductivity (PPC) on free-carrier concentrations and mobilities. The increase in the two-dimensional electron-gas (2DEG) concentration is significantly smaller than the apparent one derived from single carrier analysis of the Hall coefficient. In the two types of structures investigated, delta doped and continuously doped barrier, the apparent concentration almost doubles following illumination, while analysis reveals an increase of about 20% in the 2DEG. The effect of PPC on mobility depends on the structure. For the sample with a continuously doped barrier the mobility in the quantum well more than doubles. This increase is attributed to the effective screening of the ionized donors by the large electron concentration in the barrier. In the delta doped barrier sample the mobility is reduced by almost a factor of 2. This decrease is probably caused by strong coupling between the two wells, as is demonstrated by self-consistent analysis.

Schacham, S. E.↗

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗

Computation of Reacting Flows in Combustion Processes

The main objective of this research was to develop an efficient three-dimensional computer code for chemically reacting flows. The main computer code developed is ALLSPD-3D. The ALLSPD-3D computer program is developed for the calculation of three-dimensional, chemically reacting flows with sprays. The ALL-SPD code employs a coupled, strongly implicit solution procedure for turbulent spray combustion flows. A stochastic droplet model and an efficient method for treatment of the spray source terms in the gas-phase equations are used to calculate the evaporating liquid sprays. The chemistry treatment in the code is general enough that an arbitrary number of reaction and species can be defined by the users. Also, it is written in generalized curvilinear coordinates with both multi-block and flexible internal blockage capabilities to handle complex geometries. In addition, for general industrial combustion applications, the code provides both dilution and transpiration cooling capabilities. The ALLSPD algorithm, which employs the preconditioning and eigenvalue rescaling techniques, is capable of providing efficient solution for flows with a wide range of Mach numbers. Although written for three-dimensional flows in general, the code can be used for two-dimensional and axisymmetric flow computations as well. The code is written in such a way that it can be run in various computer platforms (supercomputers, workstations and parallel processors) and the GUI (Graphical User Interface) should provide a user-friendly tool in setting up and running the code.

Keith, Theo G., Jr.↗

Prediction of Transitional Flows in the Low Pressure Turbine

Current turbulence models tend to give too early and too short a length of flow transition to turbulence, and hence fail to predict flow separation induced by the adverse pressure gradients and streamline flow curvatures. Our discussion will focus on the development and validation of transition models. The baseline data for model comparisons are the T3 series, which include a range of free-stream turbulence intensity and cover zero-pressure gradient to aft-loaded turbine pressure gradient flows. The method will be based on the conditioned N-S equations and a transport equation for the intermittency factor. First, several of the most popular 2-equation models in predicting flow transition are examined: k-e [Launder-Sharina], k-w [Wilcox], Lien-Leschiziner and SST [Menter] models. All models fail to predict the onset and the length of transition, even for the simplest flat plate with zero-pressure gradient(T3A). Although the predicted onset position of transition can be varied by providing different inlet turbulent energy dissipation rates, the appropriate inlet conditions for turbulence quantities should be adjusted to match the decay of the free-stream turbulence. Arguably, one may adjust the low-Reynolds-number part of the model to predict transition. This approach has so far not been very successful. However, we have found that the low-Reynolds-number model of Launder and Sharma [1974], which is an improved version of Jones and Launder [1972] gave the best overall performance. The Launder and Sharma model was designed to capture flow re-laminarization (a reverse of flow transition), but tends to give rise to a too early and too fast transition in comparison with the physical transition. The three test cases were for flows with zero pressure gradient but with different free-stream turbulent intensities. The same can be said about the model when considering flows subject to pressure gradient(T3C1). To capture the effects of transition using existing turbulence models, one approach is to make use of the concept of the intermittency to predict the flow transition. It was originally based on the intermittency distribution of Narasimha [1957], and then gradually evolved into a transport equation for the intermittency factor. Gostelow and associates [1994,1995] have made some improvements to Narasimha's method in an attempt to account for both favorable and adverse pressure gradients. Their approach is based on a linear, explicit combination of laminar and turbulent solutions. This approach fails to predict the overshoot of the skin friction on a flat plate near the end of transition zone, even though the length of transition is well predicted. The major flaw of Gostelow's approach is that it assumes the non-turbulent part being the laminar solution and the turbulent part being the turbulent solution and they do not interact across the transitional region. The technique in condition averaging the flow equations in intermittent flows was first introduced by Libby [1975] and Dopazo [1977] and further refined by Dick and associates [1988, 1996]. This approach employs two sets of transport equations for the non-turbulent part and the other for the turbulent part. The advantage of this approach is that it allows the interaction of non-turbulent and turbulent velocities through the introduction of additional source terms in the continuity and momentum equations for the non-turbulent and turbulent velocities. However, the strong coupling of the two sets of equations has caused some numerical difficulties, which requires special attention. The prediction of the skin friction can be improved by this approach via the implicit coupling of non-turbulent and turbulent velocity flelds. Another improvement of the interrmittency model can be further made by allowing the intermittency to vary in the cross-stream direction. This is one step prior to testing any proposal for the transport equation for the intermittency factor. Instead of solving the transport equation for the intermittency factor, the distribution for the intermittency factor is prescribed by Klebanoff's empirical formula [1955]. The skin friction is very well predicted by this new modification, including the overshoot of the profile near the end of the transition zone. The outcome of this study is very encouraging since it indicates that the proper description of the intermittency distribution is the key to the success of the model prediction. This study will be used to guide us on the modelling of the intermittency transport equation.

Huang, George↗

Magnetic Roots and the Driving of Extended Coronal Heating

We report results from a continuation of a previous study, in which we found large bright coronal loops within active regions and extending from active regions that have one end rooted near an island of included magnetic polarity that is a site of enhanced coronal heating and microflares. This suggested that magnetic activity such as microflaring results in enhanced heating in both the compact core field around the island and in the large loops extending from it. We might expect that the intensity variations due to enhanced heating in the compact and extended structures would be correlated. However, although some ex- tended loops do respond to the largest events taking place in the core fields near their feet, they do not show a clear response to most smaller individual events nor to the overall envelope of coronal heating activity in the core fields at their feet as determined from longer-term observations. Thus, while it is clear that the extended loops' heating is being driven from their ends at the magnetic islands, much of this heating is apparently by some form of footpoint activity that is not strongly coupled to the heating in the footpoint core fields. One possibility is that the remote heating in the extended loops is driven by reconnection at the magnetic null over the island, and that this reconnection is driven mainly by core-field activity that produces little coronal heating within the core field itself, perhaps in the manner of the numerical simulations by Karpen, Antiochos, and DeVore.

Porter, Jason G.↗

Detailed Studies on the Structure and Dynamics of Reacting Dusty Flows at Normal and Microgravity

Dusty reacting flows are of particular interest for a wide range of applications. Inert particles can alter the flammability and extinction limits of a combustible mixture. Reacting particles can release substantial amount of heat and can be used either for power generation or propulsion. Accumulation of combustible particles in air can result in explosions which, for example, can occur in grain elevators, during lumber milling and in mine galleries. Furthermore, inert particles are used as flow velocity markers in reacting flows, and their velocity is measured by non-intrusive laser diagnostic techniques. Despite their importance, dusty reacting flows have been less studied and understood compared to gas phase as well as sprays. The addition of solid particles in a flowing gas stream can lead to strong couplings between the two phases, which can be of dynamic, thermal, and chemical nature. The dynamic coupling between the two phases is caused by the inertia that causes the phases to move with different velocities. Furthermore, gravitational, thermophoretic, photophoretic, electrophoretic, diffusiophoretic, centrifugal, and magnetic forces can be exerted on the particles. In general, magnetic, electrophoretic, centrifugal, photophoretic, and diffusiophoretic can be neglected. On the other hand, thermophoretic forces, caused by steep temperature gradients, can be important. The gravitational forces are almost always present and can affect the dynamic response of large particles. Understanding and quantifying the chemical coupling between two phases is a challenging task. However, all reacting particles begin this process as inert particles, and they must be heated before they participate in the combustion process. Thus, one must first understand the interactions of inert particles in a combustion environment. The in-detail understanding of the dynamics and structure of dusty flows can be only advanced by considering simple flow geometries such as the opposed-jet, stagnation-type. In such configurations the imposed strain rate is well characterized, and the in-depth understanding of the details of the physico-chemical processes can be systematically obtained. A number of computational and experimental studies on spray and particle flows have been conducted in stagnation-type configurations. Numerically, the need for a hybrid Eulerian-Lagrangian approach has been identified by Continillo and Sirignano, and the use of such approach has allowed for the prediction of the phenomenon of droplet flow reversal. Gomez and Rosner have conducted a detailed study on the particle response in the opposed-jet configuration, and the particle thermophoretic diffusivities were determined experimentally. Sung, Law and co-workers have conducted numerical studies on the effect of strain rate and temperature gradients on the dynamics of inert particles, as a way of understanding potential errors in experimental LDV data that may arise from thermophoretic forces. This investigation is a combined experimental and numerical study on the details of reacting dusty flows. The specific tasks are: (1) Experimental determination of laminar flame speeds, and extinction strain rates of dusty flows at normal- and micro-gravity as functions of the particle type, particle initial diameter, particle initial number density, and gas phase chemical composition; (2) Detailed numerical simulation of the experiments. Results are compared with experiments and the adequacy of theoretical models is assessed; and (3) Provision of enhanced insight into the thermo-chemical coupling between the two phases.

Andac, M. Gurhan↗

Performance of the GEOS-3/Terra Data Assimilation System in the Northern Stratospheric Winter 1999/2000

As part of NASA's support for the Terra satellite, which became operational in January 2000, the Data Assimilation Office introduced a new version of the GEOS data assimilation system (DAS) in November 1999. This system, GEOS-3/Terra, differs from its predecessor in several ways, notably through an increase in horizontal resolution (from 2-by-2.5 degrees to 1-by-1 degree), a slightly lower upper boundary (0.1 instead of 0.01hPa) with fewer levels (48 as opposed to 70), and substantial changes to the tropospheric physics package. This paper will address the performance of the GEOS-3/Terra DAS in the stratosphere. it focusses on the analyses (produced four times daily) and the five-day forecasts (produced twice daily). These were important for the meteorological support of the SAGE-3 Ozone Loss and Validation Experiment, based in Kiruna, Northern Sweden, in the winter of 1999/2000. It is shown that the analyses of basic meteorological fields (temperature, geopotential height, and horizontal wind) are in good agreement with those from other centers. The analyses captured the cold polar vortex which persisted through most of the winter. It is shown that forecasts (up to five days) tend to have a warm bias, which is important for the prediction of polar stratospheric clouds, which are triggered by temperatures of 195K (or lower). The importance of accurate upper tropospheric forecasts in predicting the stratospheric flow is highlighted in the context of the evolution of the shape of the stratospheric polar vortex. A prominent blocking high in the Atlantic region in January was an important factor determining the shape of the distorted lower stratospheric vortex; the predictive skill of these features was strongly coupled in the GEOS-3/Terra system.

Pawson, S.↗