Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion couples”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

A Meteoroid Handbook for Aerospace Engineers and Managers

At the beginning of the Space Age, spacecraft designers and mission planners were very concerned about meteoroids. They envisioned vehicles being ripped to pieces by streams of fast-moving space rocks, a notion promoted by the science fiction novels and movies of the time. The reality is, of course, different—the meteoroid streams that produce meteor showers are not dense by laypeople’s standards, having spatial densities of just a handful of particles per cubic kilometer, even during meteor outbursts. The ever-present, diffuse, sporadic background, which produces observed meteor rates of only 5 to 8 meteors per hour, makes up 90% of the meteoroid risk to spacecraft that spend at least a year in low Earth orbit (LEO), whereas the visually spectacular but short-lived meteor showers make up the other 10%. Still, meteoroids do pose a significant risk to spacecraft. At Earth, they can travel 12 to 72 km/s. These high speeds cause even small meteoroids to carry enormous kinetic energy, making them capable of doing serious damage to spacecraft. For example, a 1-mm-diameter meteoroid moving at 25 km/s can inflict the same damage as a bullet fired from a 0.357 Magnum pistol. An exterior wire can be severed by a 0.1-mm (100 mm) particle, a spacesuit can be penetrated by a 0.5-mm meteoroid, and an unshielded pressure wall (like the cabin of the Space Shuttle) can be perforated by centimeter-sized particles. Along with mechanical damage, meteoroids can also cause other types of spacecraft anomalies. Meteoroids can transfer their momentum to the spacecraft, which can destroy or damage equipment such as shunt resistors and charge-coupled device (CCD) detectors with a clear view of space. Meteoroid impacts can also generate plasma. The impact vaporizes material, producing a crater and an expanding plasma, which can in turn provide a conductive path for any charge accumulated on the spacecraft. This effect is thought to be responsible for the demise of a satellite in one case: the OLYMPUS communications satellite was sent tumbling out of control during the 1993 Perseid outburst, and a Perseid meteoroid strike has been posited as a possible cause (McDonnell et al. 1993; Caswell et al. 1995). Other researchers have suggested that very fast meteoroids could produce a small electromagnetic pulse capable of disrupting spacecraft function (Close et al. 2010).

Moorhead, A.↗
NASA NTRS2011

Ion-Neutral Coupling in Solar Prominences

Interactions between ions and neutrals in a partially ionized plasma are important throughout heliophysics, including near the solar surface in prominences. Understanding how ion-neutral coupling affects formation, support, structure, and dynamics of prominences will advance our physical understanding of magnetized systems involving a transition from a weakly ionized dense gas to a fully ionized tenuous plasma. We address the fundamental physics of prominence support, which is normally described in terms of a magnetic force on the prominence plasma that balances the solar gravitational force, and the implications for observations. Because the prominence plasma is only partially ionized, it is necessary to consider the support of the both the ionized and neutral components. Support of the neutrals is accomplished through a frictional interaction between the neutral and ionized components of the plasma, and its efficacy depends strongly on the degree of ionization of the plasma. More specifically, the frictional force is proportional to the relative flow of neutral and ion species, and for a sufficiently weakly ionized plasma, this flow must be relatively large to produce a frictional force that balances gravity. A large relative flow, of course, implies significant draining of neutral particles from the prominence. We evaluate the importance of this draining effect for a hydrogen-helium plasma, and consider the observational evidence for cross-field diffusion of neutral prominence material.

Gilbert, Holly↗
NASA NTRS2004

Coastal Research Imaging Spectrometer

The Coastal Research Imaging Spectrometer (CRIS) is an airborne remote sensing system designed specifically for research on the physical, chemical, and biological characteristics of coastal waters. The CRIS includes a visible-light hyperspectral imaging subsystem for measuring the color of water, which contains information on the biota, sediment, and nutrient contents of the water. The CRIS also includes an infrared imaging subsystem, which provides information on the temperature of the water. The combination of measurements enables investigation of biological effects of both natural and artificial flows of water from land into the ocean, including diffuse and point-source flows that may contain biological and/or chemical pollutants. Temperature is an important element of such measurements because temperature contrasts can often be used to distinguish among flows from different sources: for example, a sewage outflow could manifest itself in spectral images as a local high-temperature anomaly. Both the visible and infrared subsystems scan in pushbroom mode: that is, an aircraft carrying the system moves along a ground track, the system is aimed downward, and image data are acquired in across-track linear arrays of pixels. Both subsystems operate at a frame rate of 30 Hz. The infrared and visible-light optics are adjusted so that both subsystems are aimed at the same moving swath, which has across-track angular width of 15 . Data from the infrared and visible imaging subsystems are stored in the same file along with aircraft- position data acquired by a Global Positioning System receiver. The combination of the three sets of data is used to construct infrared and hyperspectral maps of scanned areas (see figure). The visible subsystem is based on a grating spectrograph and a rapid-readout charge-coupled-device camera. Images of the swatch are acquired in 256 spectral bands at wavelengths from 400 to 800 nm. The infrared subsystem, which is sensitive in a single wavelength band of 8 to 10 m, is based on a focal-plane array of HgCdTe photodetectors that are cooled to an operating temperature of 77 K by use of a closed-Stirling-cycle mechanical cooler. The nonuniformities of the HgCdTe photodetector array are small enough that the raw pixel data from the infrared subsystem can be used to recognize temperature differences on the order of 1 C. By use of a built-in blackbody calibration source that can be switched into the field of view, one can obtain bias and gain offset terms for individual pixels, making it possible to offset the effects of nonuniformities sufficiently to enable the measurement of temperature differences as small as 0.1 C.

Lucey, Paul G.↗
NASA NTRS1996

Self-Consistent Simulation of the Brownian Stage of Dust Growth

It is a widely accepted view that in proto-planetary accretion disks the collision and following sticking of dust particles embedded in the gas eventually leads to the formation of planetesimals (coagulation). For the smallest dust grains, Brownian motion is assumed to be the dominant source of their relative velocities leading to collisions between these dust grains. As the dust grains grow they eventually couple to the turbulent motion of the gas which then drives the coagulation much more efficiently. Many numerical coagulation simulations have been carried out to calculate the fractal dimension of the aggregates, which determines the duration of the ineffective Brownian stage of growth. Predominantly on-lattice and off-lattice methods were used. However, both methods require simplification of the astrophysical conditions. The aggregates found by those methods had a fractal dimension of approximately 2 which is equivalent to a constant, mass-independent friction time. If this value were valid for the conditions in an accretion disk, this would mean that the coagulation process would finally 'freeze out' and the growth of a planetesimal would be impossible within the lifetime of an accretion disk. In order to investigate whether this fractal dimension is model independent, we simulate self-consistently the Brownian stage of the coagulation by an N-particle code. This method has the advantage that no further assumptions about homogeneity of the dust have to be made. In our model, the dust grains are considered as aggregates built up of spheres. The equation of motion of the dust grains is based on the probability density for the diffusive transport within the gas atmosphere. Because of the very low number density of the dust grains, only 2-body-collisions have to be considered. As the Brownian stage of growth is very inefficient, the system is to be simulated over long periods of time. In order to find close particle pairs of the system which are most likely to undergo a collision, we use a particle-in-cell (PIC) method for the early stages of the simulation where the system is still very homogeneous and a tree method later when the particles are more clustered.

Kempf, S.↗
NASA NTRS2018

Quantifying the Effect of Mixing on the Mean Age of Air in CCMVal-2 and CCMI-1 Models

The stratospheric age of air (AoA) is a useful measure of the overall capabilities of a general circulation model (GCM) to simulate stratospheric transport. Previous studies have reported a large spread in the simulation of AoA by GCMs and coupled chemistry-climate models (CCMs). Compared to observational estimates, simulated AoA is mostly too low. Here we attempt to untangle the processes that lead to the AoA differences between the models and between models and observations. AoA is influenced by both mean transport by the residual circulation and two-way mixing; we quantify the effects of these processes using data from the CCM inter-comparison projects CCMVal-2 (Chemistry-Climate Model Validation Activity 2) and CCMI-1 (Chemistry-Climate Model Initiative, phase 1). Transport along the residual circulation is measured by the residual circulation transit time (RCTT). We interpret the difference between AoA and RCTT as additional aging by mixing. Aging by mixing thus includes mixing on both the resolved and subgrid scale. We find that the spread in AoA between the models is primarily caused by differences in the effects of mixing and only to some extent by differences in residual circulation strength. These effects are quantified by the mixing efficiency, a measure of the relative increase in AoA by mixing. The mixing efficiency varies strongly between the models from 0.24 to 1.02. We show that the mixing efficiency is not only controlled by horizontal mixing, but by vertical mixing and vertical diffusion as well. Possible causes for the differences in the models' mixing efficiencies are discussed. Differences in subgrid-scale mixing (including differences in advection schemes and model resolutions) likely contribute to the differences in mixing efficiency. However, differences in the relative contribution of resolved versus parameterized wave forcing do not appear to be related to differences in mixing efficiency or AoA.

CCMI↗
NASA NTRSDate not supplied

Implementation of Charged Particle Behavior in Discrete Element Method (DEM) Simulations

Lunar landers will agitate the surface of the Moon with an exhaust plume during descent which will, due to the particulate nature of the lunar regolith, loosen and eject grains from the surface. This ejection is not only coupled with the charged plume gas, but also results in significant particle-particle interactions. Settling of these grains after plume effects have subsided takes much longer than expected in a ballistic trajectory. The prevailing hypothesis attributes this behavior to the accumulated charge on the particles. We are thus developing a discrete element method (DEM) approach to explore these charged particle interactions on the lunar surface. The Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) Improved for General Granular and Granular Heat Transfer Simulations (LIGGGHTS) software package provides a DEM modeling framework for granular interactions. It includes many complexities such as non-spherical particle shapes, cohesion and frictional forces, and heat transfer, but has no provision for inter-particle electrostatic forces and charge transfer that are important to examine in the lunar environment. In this work, a standard Coulomb potential and a Yukawa potential are integrated into the LIGGGHTS framework to provide a basis for particle-particle electrostatic interactions, as well as a gravitational potential to enable inter-grain gravitational interactions. A preliminary approach to charge transfer between grains incorporating properties such as work function and electrical conductivity to the library of available material characteristics will be presented. Several scenarios have been simulated that include charged particle interactions within a diffuse granular gas, settling of charged grains into a regolith bed, sliding of granular material along an incline, and vibration of settled grains to produce a behavior similar to fluidization. There are numerous challenges to incorporate realistic interactions between complex lunar particles. Currently, grains are modeled to behave as if the entirety of the charge acts at the center of mass, such as conductors with spherical symmetry and insulators with homogeneously distributed charge. We are developing improvements that will include the use of non-spherical particle geometries, as well as reasonable approximations of insulating/dielectric materials that have non-uniform charge distributions. The cases simulated thus far will be examined in a relevant environment within a vacuum chamber to validate the simulations. These simulations will be bounded by experiments utilizing high-speed camera observations of the motion for validation. The grains in the experiment will exchange charge during their motion and this can be quantified by collection within a charge measurement device such as a Faraday cup. Such a device may be modeled within the software by defining an integration region and computing the contained charge as a function of simulation time, allowing for side-by-side comparison of simulated and measured bulk charging results. Any differences will be reconciled by updating the mathematical mechanisms described within the simulation suite. Successfully combining results from experiments within a relevant environment into the LIGGGHTS framework will improve modeling of the charged grain dynamics experienced on the Moon to provide insights into dust behavior for future lunar exploration missions.

Electrostatics↗
NASA NTRS2004

Synthesis and Characterization of Ru(II) Tris(1,1O-phenanthroline)-Electron Acceptor Dyads Incorporating the 4-benzoyl-N-methylpyridinium Cation or N-Benzyl-N'-methyl-viologen. Improving the Dynamic Range, Sensitivity and Response Time of Sol-Gel Based Optical Oxygen Sensors

The title compounds (1 and 2, above) were synthesized by Sonogashira coupling reactions of appropriate Ru(1I) complexes with the electron a cceptors. Characterization was conducted in solution and in frozen ma trices. Finally, the title compounds were evaluated as dopants of sol-gel materials. It was found that the intramolecular quenching efficie ncy of 4-benzoyl-Nmethylpyridinium cation in solution depends on the solvent: photoluminescence is quenched completely in CH,CN, but not i n methanol or ethanol. On the other hand, intramolecular emission que nching by 4-benzyl-N-methyl viologen is complete in all solvents. The difference between the two quenchers is traced electrochemically to t he solvation of the 4-benzoyl-Nmethylpyridiniums by alcohol. In froze n matrices or adsorbed on the surfaces of silica aerogel, both Ru(I1) complex/electron acceptor dyads of this study are photoluminescent, and the absence of quenching has been traced to the environmental rigi dity. When doped aerogels are cooled at 77 K, the emission intensity increases by approximately 4x, and the spectra shift to the blue, analogous to what is observed with Ru(I1) complexes in solutions undergoi ng fluid-to-rigid transition. However, in contrast to frozen solution s, the luminescent moieties in the bulk of aerogels kept at low tempe ratures are still accessible to gas-phase quenchers diffusing through the mesopores, leading to more sensitive platforms for sensors than o ther room-temperature configurations. Thus the photoluminescence of o ur Ru(I1) complex dyads adsorbed on aerogel is quenchable by O2 both at room temperature and at 77 K. Furthermore, it was also found that O 2 modulates the photoluminescence of aerogels doped with 4-benzoyl -N -methylpyridinium-based dyads over a wider dynamic range compared wi th aerogels doped with either our vislogen-based dyads or with Ru(I1) tris(1,lO-phenanthroline) itself.

Leventis, Nicholas↗
NASA NTRS1996

Soot and Radiation Measurements in Microgravity Jet Diffusion Flames

The subject of soot formation and radiation heat transfer in microgravity jet diffusion flames is important not only for the understanding of fundamental transport processes involved but also for providing findings relevant to spacecraft fire safety and soot emissions and radiant heat loads of combustors used in air-breathing propulsion systems. Our objectives are to measure and model soot volume fraction, temperature, and radiative heat fluxes in microgravity jet diffusion flames. For this four-year project, we have successfully completed three tasks, which have resulted in new research methodologies and original results. First is the implementation of a thermophoretic soot sampling technique for measuring particle size and aggregate morphology in drop-tower and other reduced gravity experiments. In those laminar flames studied, we found that microgravity soot aggregates typically consist of more primary particles and primary particles are larger in size than those under normal gravity. Comparisons based on data obtained from limited samples show that the soot aggregate's fractal dimension varies within +/- 20% of its typical value of 1.75, with no clear trends between normal and reduced gravity conditions. Second is the development and implementation of a new imaging absorption technique. By properly expanding and spatially-filtering the laser beam to image the flame absorption on a CCD camera and applying numerical smoothing procedures, this technique is capable of measuring instantaneous full-field soot volume fractions. Results from this technique have shown the significant differences in local soot volume fraction, smoking point, and flame shape between normal and reduced gravity flames. We observed that some laminar flames become open-tipped and smoking under microgravity. The third task we completed is the development of a computer program which integrates and couples flame structure, soot formation, and flame radiation analyses together. We found good agreements between model predictions and experimental data for laminar and turbulent flames under both normal and reduced gravity. We have also tested in the laboratory the techniques of rapid-insertion fine-wire thermocouples and emission pyrometry for temperature measurements. These techniques as well as laser Doppler velocimetry and spectral radiative intensity measurement have been proposed to provide valuable data and improve the modeling analyses.

Ku, Jerry C.↗
NASA NTRS2014

An Experimental and Computational Study on Soot Formation in a Coflow Jet Flame Under Microgravity and Normal Gravity

Upon the completion of the Structure and Liftoff in Combustion Experiment (SLICE) in March 2012, a comprehensive and unique set of microgravity coflow diffusion flame data was obtained. This data covers a range of conditions from weak flames near extinction to strong, highly sooting flames, and enabled the study of gravitational effects on phenomena such as liftoff, blowout and soot formation. The microgravity experiment was carried out in the Microgravity Science Glovebox (MSG) on board the International Space Station (ISS), while the normal gravity experiment was performed at Yale utilizing a copy of the flight hardware. Computational simulations of microgravity and normal gravity flames were also carried out to facilitate understanding of the experimental observations. This paper focuses on the different sooting behaviors of CH4 coflow jet flames in microgravity and normal gravity. The unique set of data serves as an excellent test case for developing more accurate computational models.Experimentally, the flame shape and size, lift-off height, and soot temperature were determined from line-of-sight flame emission images taken with a color digital camera. Soot volume fraction was determined by performing an absolute light calibration using the incandescence from a flame-heated thermocouple. Computationally, the MC-Smooth vorticity-velocity formulation was employed to describe the chemically reacting flow, and the soot evolution was modeled by the sectional aerosol equations. The governing equations and boundary conditions were discretized on an axisymmetric computational domain by finite differences, and the resulting system of fully coupled, highly nonlinear equations was solved by a damped, modified Newtons method. The microgravity sooting flames were found to have lower soot temperatures and higher volume fraction than their normal gravity counterparts. The soot distribution tends to shift from the centerline of the flame to the wings from normal gravity to microgravity.

microgravity↗
NASA NTRS2000

The Sensitivity of Model Ozone to Advective and Photochemical Processes in the High Latitude Winter Lower Stratosphere

Three dimensional chemistry and transport models (CTMs) contain a set of coupled continuity equations which describe the evolution of constituents such as ozone and other minor species which affect ozone. Both advection and photochemical processes contribute to constituent evolution, and a CTM provides a means to evaluate these contributions separately. Such evaluation is particularly useful when both terms are important to the modeled tendency. An example is the ozone tendency in the high latitude winter lower stratosphere, where advection tends to increase ozone, and catalytic processes involving chlorine radicals tend to decrease ozone. The Goddard three dimensional chemistry and transport model uses meteorological fields from the Goddard Earth Observing System Data Assimilation System, thus the modeled ozone evolution may reproduce the observed evolution and provide a test of the model representation of photochemical processes if the transport is shown to be modeled appropriately. We have investigated the model advection further using diabatic trajectory calculations. For long lived constituents such as N2O, the model field for a particular time on a potential temperature surface is compared with a field produced by calculating 15 day back trajectories for a fixed latitude longitude grid, and mapping model N2O at the terminus of the back trajectories onto the initial grid. This provides a quantitative means to evaluate two aspects of the CTM transport: one, the model horizontal gradient between middle latitudes and the polar vortex is compared with the gradient produced using the non-diffusive trajectory calculation; two, the model vertical advection, which is produced by the divergence of the horizontal winds, is compared with the vertical transport expected from diabatic cooling.

Douglass, A.↗
NASA NTRS2022

A Coupled Evaluation of Operational MODIS and Model Aerosol Products for Maritime Environments Using Sun Photometry: Evaluation of the Fine and Coarse Mode

Although satellite retrievals and data assimilation have progressed to where there is a good skill for monitoring maritime Aerosol Optical Depth (AOD), there remains uncertainty in achieving further degrees of freedom, such as distinguishing fine and coarse mode dominated species in maritime environments (e.g., coarse mode sea salt and dust versus fine mode terrestrial anthropogenic emissions, biomass burning, and maritime secondary production). For the years 2016 through 2019, we performed an analysis of 550 nm total AOD550, fine mode AOD (FAOD550; also known as FM AOD in the literature), coarse mode AOD (CAOD550), and fine mode fraction (η550) between Moderate Resolution Spectral Imaging Radiometer (MODIS) V6.1 MOD/MYD04 dark target aerosol retrievals and the International Cooperative for Aerosol Prediction (ICAP) core four multi-model consensus (C4C) of analyses/short term forecasts that assimilate total MODIS AOD550. Differences were adjudicated by the global shipboard Maritime Aerosol Network (MAN) and selected island AERONET sun photometer observations with the application of the spectral deconvolution algorithm (SDA). Through a series of conditional and regional analyses, we found divergence included regions of terrestrial influence and latitudinal dependencies in the remote oceans. Notably, MODIS and the C4C and its members, while having good correlations overall, have a persistent +0.04 to +0.02 biases relative to MAN and AERONET for typical AOD550 values (84th% < 0.28), with the C4C underestimating significant events thereafter. Second, high biases in AOD550 are largely associated with the attribution of the fine mode in satellites and models alike. Thus, both MODIS and C4C members are systematically overestimating AOD550 and FAOD550 but perform better in characterizing the CAOD550. Third, for MODIS, findings are consistent with previous reports of a high bias in the retrieved Ångström Exponent, and we diagnosed both the optical model and cloud masking as likely causal factors for the AOD550 and FAOD550 high bias, whereas for the C4C, it is likely from secondary overproduction and perhaps numerical diffusion. Fourth, while there is no wind-speed-dependent bias for surface winds <12 m s−1, the C4C and MODIS AOD550s also overestimate CAOD550 and FAOD550, respectively, for wind speeds above 12 m/s. Finally, sampling bias inherent in MAN, as well as other circumstantial evidence, suggests biases in MODIS are likely MODIS and the C4C products have their own strengths and challenges for a given climate application and discuss needed research. even larger than what was diagnosed here. We conclude with a discussion on how

Aerosols↗
NASA NTRS2003

The Chemical Composition and Gas-to-Dust Mass Ratio of Nearby Interstellar Matter

We use recent results on interstellar gas toward nearby stars and interstellar by-products within the solar system to select among the equilibrium radiative transfer models of the nearest interstellar material presented in Slavin & Frisch. For the assumption that O/H - 400 parts per million, models 2 and 8 are found to yield good fits to available data on interstellar material inside and outside of the heliosphere, with the exception of the Ne abundance in the pickup ion and anomalous cosmic-ray populations. For these models, the interstellar medium (ISM) at the entry point to the heliosphere has n(H(sup 0)) = 0.202-0.208/cu cm, n(He(sup 0) = 0.0137-0.0152/cu cm, and ionizations X(H) = 0.29-0.30, X(He) = 0.47-0.51. These best models suggest that the chemical composition of the nearby ISM is approx.60%-70% subsolar if S is undepleted. Both H(0) and H(+) need to be included when evaluating abundances of ions found in warm diffuse clouds. Models 2 and 8 yield an H filtration factor of approx.0.46. Gas-to-dust mass ratios for the ISM toward epsilon CMa are R(sub gd) = 178-183 for solar abundances of Holweger or R(sub gd) = 611-657 for an interstellar abundance standard 70% solar. Direct observations of dust grains in the solar system by Ulysses and Galileo yield R(sub gd) appr0x. 115 for models 2 and 8, supporting earlier results (Frisch and coworkers). If the local ISM abundances are subsolar, then gas and dust are decoupled over small spatial scales. The inferred variation in R(sub gd) over parsec length scales is consistent with the fact that the ISM near the Sun is part of a dynamically active cluster of cloudlets flowing away from the Sco-Cen association. Observations toward stars within approx.500 pc show that R(sub gd) correlates with the percentage of the dust mass that is carried by iron, suggesting that an Fe-rich grain core (by mass) remains after grain destruction. Evidently large dust grains (>10(exp -13) g) and small dust grains (<10(exp -13) g) are not well mixed over parsec length spatial scales in the ISM. It also appears that very small C-dominated dust grains have been destroyed in the ISM within several parsecs of the Sun, since C appears to be essentially undepleted. However, if gas-dust coupling breaks down over the cloud lifetime, the missing mass arguments applied here to determine R(sub gd) and dust grain mineralogy are not appropriate.

Frisch, Priscilla C.↗
NASA NTRS1993

Study and modeling of finite rate chemistry effects in turbulent non-premixed flames

The development of numerical models that reflect some of the most important features of turbulent reacting flows requires information about the behavior of key quantities in well defined combustion regimes. In turbulent flames, the coupling between turbulent and chemical processes is so strong that it is extremely difficult to isolate the role played by one individual physical phenomenon. Direct numerical simulation (hereafter DNS) allows us to study in detail the turbulence-chemistry interaction in some restricted but completely defined situations. Globally, non-premixed flames are controlled by two limiting regimes: the fast chemistry case, where the turbulent flame can be pictured as a random distribution of local chemical equilibrium problems; and the slow chemistry case, where the chemistry integrates in time the turbulent fluctuations. The Damkoehler number, ratio of a mechanical time scale to chemical time scale, is used to distinguish between these regimes. Today most of the industrial computer codes are able to perform predictions in the hypothesis of local equilibrium chemistry using a presumed shape for the probability density function (pdt) of the conserved scalar. However, the finite rate chemistry situation is of great interest because industrial burners usually generate regimes in which, at some points, the flame is undergoing local extinction or at least non-equilibrium situations. Moreover, this variety of situations strongly influences the production of pollutants. To quantify finite rate chemistry effect, the interaction between a non-premixed flame and a free decaying turbulence is studied using DNS. The attention is focused on the dynamic of extinction, and an attempt is made to quantify the effect of the reaction on the small scale mixing process. The unequal diffusivity effect is also addressed. Finally, a simple turbulent combustion model based on the DNS observations and tractable in real flow configurations is proposed.

Vervisch, Luc↗
NASA NTRS2002

Shadowgraph Study of Gradient Driven Fluctuations

A fluid or fluid mixture, subjected to a vertical temperature and/or concentration gradient in a gravitational field, exhibits greatly enhanced light scattering at small angles. This effect is caused by coupling between the vertical velocity fluctuations due to thermal energy and the vertically varying refractive index. Physically, small upward or downward moving regions will be displaced into fluid having a refractive index different from that of the moving region, thus giving rise to the enhanced scattering. The scattered intensity is predicted to vary with scattering wave vector q, as q(sup -4), for sufficiently large q, but the divergence is quenched by gravity at small q. In the absence of gravity, the long wavelength fluctuations responsible for the enhanced scattering are predicted to grow until limited by the sample dimensions. It is thus of interest to measure the mean-squared amplitude of such fluctuations in the microgravity environment for comparison with existing theory and ground based measurements. The relevant wave vectors are extremely small, making traditional low-angle light scattering difficult or impossible because of stray elastically scattered light generated by optical surfaces. An alternative technique is offered by the shadowgraph method, which is normally used to visualize fluid flows, but which can also serve as a quantitative tool to measure fluctuations. A somewhat novel shadowgraph apparatus and the necessary data analysis methods will be described. The apparatus uses a spatially coherent, but temporally incoherent, light source consisting of a super-luminescent diode coupled to a single-mode optical fiber in order to achieve extremely high spatial resolution, while avoiding effects caused by interference of light reflected from the various optical surfaces that are present when using laser sources. Results obtained for a critical mixture of aniline and cyclohexane subjected to a vertical temperature gradient will be presented. The sample was confined between two horizontal parallel sapphire plates with a vertical spacing of 1 mm. The temperatures of the sapphire plates were controlled by independent circulating water loops that used Peltier devices to add or remove heat from the room air as required. For a mixture with a temperature gradient, two effects are involved in generating the vertical refractive index gradient, namely thermal expansion and the Soret effect, which generates a concentration gradient in response to the applied temperature gradient. For the aniline/cyclohexane system, the denser component (aniline) migrates toward the colder surface. Consequently, when heating from above, both effects result in the sample density decreasing with altitude and are stabilizing in the sense that no convective motion occurs regardless of the magnitude of the applied temperature gradient. The Soret effect is strong near a binary liquid critical point, and thus the dominant effect is due to the induced concentration gradient. The results clearly show the divergence at low q and the predicted gravitational quenching. Results obtained for different applied temperature gradients at varying temperature differences from the critical temperature, clearly demonstrate the predicted divergence of the thermal diffusion ratio. Thus, the more closely the critical point is approached, the smaller becomes the temperature gradient required to generate the same signal. Two different methods have been used to generate pure concentration gradients. In the first, a sample cell was filled with a single fluid, ethylene glycol, and a denser miscible fluid, water, was added from below thus establishing a sharp interface to begin the experiment. As time went on the two fluids diffused into each other, and large amplitude fluctuations were clearly observed at low q. The effects of gravitational quenching were also evident. In the second method, the aniline/cyclohexane sample was used, and after applying a vertical temperature gradient for several hours, the top and bottom temperatures were set equal and the thermal gradient died on a time scale of seconds, leaving the Soret induced concentration gradient in place. Again, large-scale fluctuations were observed and died away slowly in amplitude as diffusion destroyed the initial concentration gradient.

Cannell, David↗
NASA NTRS1999

Dynamics of Single Chains of Suspended Ferrofluid Particles

We present an experimental study of the dynamics of isolated chains made of super-paramagnetic particles under the influence of a magnetic field. The motivation of this work is to understand if the chain fluctuations exist and, if it does, how does the fluctuation affect chain aggregation. We find that single chains strongly fluctuate and that the characteristic frequency of their fluctuations is inversely proportional to the magnetic field strength. The higher the field the lower the characteristic frequency of the chain fluctuations. In the high magnetic field limit, chains behave like rigid rods without any internal motions. In this work, we used ferrofluid particles suspended in water. These particles do not have any intrinsic magnetization. Once a magnetic field is applied, a dipole moment is induced in each particle, proportional to the magnetic field. A dipolar magnetic interaction then occurs between particles. If dipole-dipole magnetic energy is higher than the thermal energy, the result is a structure change inside the dipolar fluid. The ratio of these two energies is expressed by a coupling constant lambda as: lambda = (pi(a(exp 3))(chi(exp 2))(mu(sub 0))(H(sub 0))(exp 2))/18kT Where a is the particle radius, mu(sub 0) is the vacuum magnetic permeability, H(sub 0) the applied magnetic field, k the Boltzmann constant and T the absolute temperature. If lambda > 1, magnetic particles form chains along the field direction. The lateral coalescence of several chains may form bigger aggregates especially if the particle volume fraction is high. While many studies and applications deal with the rheological properties and the structural changes of these dipolar fluids, this work focuses on the understanding of the chain dynamics. In order to probe the chain dynamics, we used dynamic light scattering (DLS) in self-beating mode as our experimental technique. The experimental geometry is such that the scattering plane is perpendicular to the magnetic field. Therefore, only motions in this plane are probed. A very dilute sample of a ferrofluid emulsion with a particle volume fraction of 10(exp -5) is used in this experiment. We chose such a low volume fraction to avoid multiple light scattering as well as lateral chain-chain aggregation. DLS measures the dynamic structure factor S(q,t) of the sample (q is the scattering wave vector, t is the time). In the absence of the magnetic field, identical particles of ferrofluid droplets are randomly distributed and S(q,t) reduces to exp(-q(exp 2)2D(sub 0)t). D(sub 0)=(kT/(6(pi)(eta)(a)) is the diffusion coefficient of Brownian particles (where Xi = (6(pi)(eta)(a)) is the Stokes frictional coefficient of a spherical particle in a fluid of viscosity eta). If interactions or polydispersity can not be ignored, an effective diffusion coefficient is introduced. Formally, D(sub eff) is defined as: D(sub eff) = - q(exp -2) partial derivative of (ln(S(q,t)) with respect to time, as t goes to 0. D(sub eff) reduces to D(sub 0) if no interactions and only a few particles size are present. Therefore, we can use DLS to measure particle size. The particle radius was found to be a=0.23 mu m with 7% of polydispersity. In this case, if we vary the scattering angle theta (and so q) we do not have any change in the measured diffusion coefficient: it is q-independent. When a magnetic field is applied, particles aggregate into chains if lambda > 1. We first studied the kinetics of the chain formation when lambda = 406. At a fixed scattering angle, we measured diffusion coefficient D(sub eff) as a function of time. Experimentally, we find that D(sub eff) decreases monotonously with time. Physically, this means that chains are becoming longer and longer. Since we are only sensitive to motions in the scattering plane and since chains have their main axis perpendicular to this plane, the measured diffusion coefficient is the trans-verse diffusion coefficient. We can relate D(sub eff) to the mean number of particles per chain N(t) at a given time and to the diffusion coefficient of an isolated particle D(sub 0) as D(sub eff)=f(N(t))D(sub 0). Since f(N) is known from other recent work, N can be expressed as a function of the time. We found a square root dependency: N(t) proportional to the square root of t. As expected for very low volume fraction, this behavior is characteristic of a diffusion-limited aggregation as suggested by several authors and by our previous work. In this study, we focus on the dependence of the effective diffusion coefficient on the scattering angle and the magnetic field strength. After the magnetic field is applied (lambda = 406) for a long time, typically 6 hours, kinetics of chain formation becomes very slow. Chain size does not vary much over the next hour period. Thus, we can perform different interesting experiments. First, at a fixed magnetic field, we measure the effective diffusion coefficient as a function of the scattering angle (from 5 to 130 deg). Our results show that the measured diffusion coefficient increases linearly with the scattering angle: D(sub eff) proportional to q. If we do the same experiment for different lambda values, D(sub eff) depends on lambda as D(sub eff) proportional to lambda(exp -1/2). We also find for different lambda values that the same asymptotic D(sub eff) value is obtained when q approaches zero. The angle dependency of D(sub eff) suggests that an additional motion exists besides chain drifting. Chain size is constant during experiment, which was verified by measuring the same diffusion coefficient at the beginning and at the end of the angle switching. If chains are rigid, D(sub eff) is independent of q. Therefore, we found that D(sub eff) not only measures the motion of the entire chain but also its internal fluctuations. These internal motions are the fluctuations of the particles in the chain. To understand the q dependency of D(sub eff), let us look at the probing length used. In our study, the characteristic length scale probed is l=2pi/q which is in the range of 0.9<l/a<20. When l is much larger than the particle radius (l/a going to infinity) we are mainly sensitive to the center of mass diffusion of the chain. D(sub eff) is then the diffusion of the entire chain and depends only on N but not on q and lambda. The value D(sub eff) allows us to obtain the number of particles per chain N. In the opposite limit (i.e. l/a < 1), we are sensitive to motions on the size of individual particles. The main contributions to the measured diffusion coefficient come from internal motions of the chains (i.e. particles' fluctuations inside the chain). We investigated also the effect of the chain size on D(sub eff). For the same value of the magnetic interaction lambda and different chain sizes, we found that the slope of D(sub eff) versus q decreases when N increases as a power law D(sub eff) proportional to N(exp -0.7). This is an important behavior because this means that the longer the chain the slower its internal fluctuations besides the whole chain motion.

Cutillas, S.↗
NASA NTRS2014

The Comprehensive Inner Magnetosphere-Ionosphere Model

Simulation studies of the Earth's radiation belts and ring current are very useful in understanding the acceleration, transport, and loss of energetic particles. Recently, the Comprehensive Ring Current Model (CRCM) and the Radiation Belt Environment (RBE) model were merged to form a Comprehensive Inner Magnetosphere-Ionosphere (CIMI) model. CIMI solves for many essential quantities in the inner magnetosphere, including ion and electron distributions in the ring current and radiation belts, plasmaspheric density, Region 2 currents, convection potential, and precipitation in the ionosphere. It incorporates whistler mode chorus and hiss wave diffusion of energetic electrons in energy, pitch angle, and cross terms. CIMI thus represents a comprehensive model that considers the effects of the ring current and plasmasphere on the radiation belts. We have performed a CIMI simulation for the storm on 5-9 April 2010 and then compared our results with data from the Two Wide-angle Imaging Neutral-atom Spectrometers and Akebono satellites. We identify the dominant energization and loss processes for the ring current and radiation belts. We find that the interactions with the whistler mode chorus waves are the main cause of the flux increase of MeV electrons during the recovery phase of this particular storm. When a self-consistent electric field from the CRCM is used, the enhancement of MeV electrons is higher than when an empirical convection model is applied. We also demonstrate how CIMI can be a powerful tool for analyzing and interpreting data from the new Van Allen Probes mission.

Inner Magnetosphere↗
NASA NTRS2000

Convective Melting of Particles in Flow Under Microgravity Conditions

Study of melting of dispersed or packed solid particles in a fluid under gravity and microgravity conditions provides benchmark information for many engineering applications such as material processing, environmental assessment and protection and space fire protection. During such processes, packed or dispersed solid particles are interacting with fluid flow at above-melting temperatures. By unmasking the buoyancy effects in coupled flow, phase change and heat transport phenomena, a better understanding of melting rate in non-thermal equilibrium, convective conditions can be studied. A series of flight experiments were conducted onboard the NASA KC-135 Microgravity Research Airplane. The Particle Melting in Plow (PMF) module was designed to allow flow through the initially packed ice particles at controlled temperature and velocity. To achieve this, a close-loop flow system was designed. Video images were taken to record the visualization of the melting process, from which a time variation of packed particle thickness distribution at different times can be obtained by the image analysis method. The fluid temperature distribution within the melting zone is measured by thermocouples. An infrared camera was mounted from the top of the test section to record the ice-water thermal images at a given location. The results from thermal images yield local temperature variation between melting solid and liquid and local Stephan number. Typical results for a number of cases are presented. The mathematical model, describing mass, energy, and momentum balance equations for the liquid and solid phases, is presented. It is found that melting rate is influenced mainly by the ratio of Reynolds number (based on the initial particle diameter) to the Froud number, and Stephan number. At the absence of gravity, Froud number approaches zero, Reynolds number and Stephan number become dominant factors governing the melting rate. The numerically determined results are compared with the experimental ones. It is found that the discrepancy between the predicted and measured melting rate is largely due to the inaccuracy in the constitutive equations for effective thermal physical properties, such as effective thermal conductivity and diffusivity, and transport properties, such as particle interaction coefficient, and local heat transfer coefficient of particles.

Stewart, C.↗
NASA NTRS2011

Acoustic Testing of Flight Hardware Using Loudspeakers: How Much do We Know About This Method of Testing?

Loudspeakers have been used for acoustic qualification of spacecrafts, reflectors, solar panels, and other acoustically responsive structures for more than a decade. Even though a lot of hardware has been acoustic tested using this method, the nature of the acoustic field generated by controlling an ensemble of speakers with and without the hardware in the test volume has not been thoroughly investigated. Limited measurements from some of the recent speaker tests used to qualify flight hardware have indicated significant spatial variation of the acoustic field within the test volume. Also structural responses have been reported to differ when similar tests were performed using reverberant chambers. Unlike the reverberant chamber acoustic test, for which the acoustic field in most chambers is known to be diffuse except below several tens of Hz where acoustic standing waves and large spatial variations exist, the characteristics of the acoustic field within the speaker test volume has not been quantified. It has only been recently that a detailed acoustic field characterization of speaker testing has been made at Jet Propulsion Laboratory (JPL) with involvement of various organizations. To address the impact of non-uniform acoustic field on structures, a series of acoustic tests were performed using a flat panel and a 3-ft cylinder exposed to the field controlled by speakers and repeated in a reverberant chamber. The analysis of the data from this exercise reveals that there are significant differences both in the acoustic field and in the structural responses. In this paper the differences between the two methods are reviewed in some detail and the over- or under-testing of articles that could pose un-anticipated structural and flight qualification issues are discussed. A framework for discussing the validity of the speaker acoustic testing method with the current control system and a path forward for improving it will be provided.

Acoustics↗