Search NASASearch

SEARCH · Search NASA

Results for “Mass diffusivity”

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 73 records · Page 4

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.

Strain-Based Assessment of Shale Caprock during Cyclic Underground Hydrogen Storage

Successful large-scale underground hydrogen storage (UHS) in depleted gas reservoirs depends on the integrity of the overlying caprock to prevent hydrogen loss during cyclic injection and depletion. Prior studies on crushed Marcellus shale, a potential caprock, indicate that cyclic hydrogen injection and depletion induces microstructural changes, increasing porosity and permeability. However, the extent of these changes in intact shale remains unclear. This study presents a strain-based experimental approach to quantify volumetric strain evolution in intact Marcellus shale matrix under unconstrained stress conditions. A quadrant-shaped shale sample without visible fractures underwent eight hydrostatic pore-pressure cycles (injection to 1500 psi and depletion to 500 psi in 250 psi steps). Linear strain gauges measured strain in three orthogonal directions. Results indicate progressive plastic strain accumulation, leading to an ∼12% increase in matrix porosity after eight cycles, with an estimated 19% increase after 30 cycles. This porosity increase follows a logarithmic trend, suggesting a diminishing effect in later cycles. Additionally, permeability and diffusive mass flux are projected to rise by ∼70% over 30 cycles, enhancing hydrogen migration risk. The shale matrix also exhibited mechanical stiffening over successive cycles, limiting large-scale deformation but not preventing porosity enhancement. A new parameter, α, was introduced to characterize shale sensitivity to cyclic loading, aiding UHS caprock assessments. These findings underscore the necessity of incorporating cyclic loading effects in UHS site selection and operational strategies to ensure long-term storage integrity.

08 HYDROGEN

Comparative Evaluation of the Ability of the MYNN‐EDMF PBL Scheme in WRF Model to Reproduce Near Surface Wind Speed Over Different Topographical Types

Abstract This study systematically evaluates the performance of the Mellor‐Yamada‐Nakanishi‐Niino‐Eddy‐Diffusion‐Mass‐Flux planetary boundary layer (PBL) scheme within the Weather Research and Forecasting (WRF) model in simulating near‐surface wind speeds across various topographies in New York State (NYS). Simulated wind speeds are compared with in‐situ measurements from 22 surface sites, grouped into six topographic categories: continental plain (CT), lakeside (LS), river valley (RV), Long Island (LI), Block Island (BI), and offshore ocean (OO). A quantitative evaluation based on Relative Euclidean Distance shows that wind speeds at the OO site are the most accurately reproduced, followed by those at LI sites, while the model performs less accurately for the remaining topographic groups. Wind speeds over CT sites tend to be overestimated by approximately 1 m/s, although their diurnal variability (DV) is well captured. In contrast, the model underestimates wind DV at LS, RV, LI, and BI sites, with the largest biases occurring at LI and BI, resulting in underestimated daytime wind speed and/or overestimated nighttime wind speed. The OO winds exhibit minimal diurnal variation, accurately captured by our WRF model. The surface wind diurnal variation is closely linked to PBL development. Among the indicators of PBL development, surface potential temperature biases most strongly correlate with wind speed biases. Our WRF model faces challenges in capturing the distinctions between winds influenced by local circulations and those over continental plains, and the significantly stronger winds at OO compared to BI. Potential causes for these biases are discussed, offering pathways for improving surface wind simulations in future.

54 ENVIRONMENTAL SCIENCES

Venus - Vertical transport rates in the visible atmosphere

The cloud particles on Venus are sufficiently small for their vertical distribution to be strongly affected by atmospheric turbulence. Reasonably firm estimates of the vertical distribution of cloud particles can be made from current interpretations of refraction, polarization, and absorption band data. These enable very firm upper limits to be placed on average vertical mass diffusion coefficients in the visible atmosphere. These upper limits are independent of the cloud particle composition or mode of formation, and of changes in the particle size distribution with altitude.

Prinn, R. G.

Magnetic storm dynamics of the thermosphere

A theoretical study of the Dst component of magnetic storms is presented. The dynamic characteristics are found significantly different for Joule dissipation and electron precipitation, leading to the conclusion that the former is probably the predominant heat source for the upper thermosphere. Composition measurements on OGO-6, which reveal markedly different characteristics in N2, O and He, can be explained on the basis of energy advection and diffusive mass transport by thermospheric winds. Essential features in the F2-region response are explicable in terms of these dynamic processes. Electric field induced motions are estimated and it is concluded that resultant adiabatic heating could be significant.

Mayr, H. G.

Some properties of upper atmosphere dynamics

A unifying description is provided of some important dynamic properties of the upper atmosphere in composition and temperature characteristic of a variety of phenomena, including diurnal and seasonal tides, magnetic storms, and momentum coupling with the magnetosphere. A theoretical multiconstituent model is used which can link the large-scale variations of composition and temperature to the dynamics and energetics of the thermosphere. Global mean properties of the thermosphere are reviewed, and an attempt is made to convey some understanding of the dynamic properties of energy and diffusive mass transport in the thermosphere. Attention is given to sources of energy for the thermosphere, the transport processes involved in the solar diurnal tide of the thermosphere, energy and particle sources for the annual tide, feedback from composition changes to wind-field and temperature variations, energy deposition in the thermosphere during magnetic storms and substorms, and momentum source signatures in the thermosphere.

Mayr, H. G.

Theory of the low-latitude boundary layer

A one-dimensional steady state fluid mechanical model is developed of the low-latitude plasma boundary layer inside the dawn and dusk magnetopause. Momentum transfer in the layer is produced by viscosity and/or mass diffusion. Coupling to the ionosphere is achieved via field-aligned currents, the magnitude of which is limited by parallel potential drops. These currents flow into and out of the ionosphere in the manner described by Iijima and Potemra. The higher-latitude (region 1) currents are associated with the boundary layer proper, while the lower-latitude (region 2) ones are associated with a region of sunward return flow adjacent to the boundary layer. The parallel potential drops have a magnitude of typically 2-3 kV and a north-south extent of 100-200 km. The calculated potential profile corresponds reasonably well to observed inverted V precipitation events.

Sonnerup, B. U. OE.

Modeling of Fluctuating Mass Flux in Variable Density Flows

The approach solves for both Reynolds and Favre averaged quantities and calculates the scalar pdf. Turbulent models used to close the governing equations are formulated to account for complex mixing and variable density effects. In addition, turbulent mass diffusivities are not assumed to be in constant proportion to turbulent momentum diffusivities. The governing equations are solved by a combination of finite-difference technique and Monte-Carlo simulation. Some preliminary results on simple variable density shear flows are presented. The differences between these results and those obtained using conventional models are discussed.

So, R. M. C.

Analytical and Numerical Solution for a Solidifying Liquid Alloy Slab

Numerical and analytical solutions are presented for the temperature and concentration distributions during the solidification of a binary liquid alloy slab. The slab is taken to be of a finite depth but infinite in the horizontal direction. The solidification process is started by withdrawing a fixed amount of heat from the lower surface of the slab. The upper surface of the slab is subjected to both radiation and convective conditions. The solution gives the concentration and temperature profiles and the interface position as a function of time. Due to the smallness of the mass diffusion coefficient in the solid, the numerical solution method breaks down whenever the ratio of the diffusivities in the solid and the liquid falls below a certain value. An analytical method is developed which gives accurate solution for any value of the diffusivity ratio.

Antar, B. N.

Dynamic behaviour of the K-corona above a type I radio source

Coronagraph images can be greatly enhanced by subtracting from the brightness of each picture point, the running median brightness over a square surrounding each picture point. The application of this technique to coronal analysis appears to be new. After such enhancement, images recorded by the SMM Coronagraph/Polarimeter reveal a coronal transient with previously undescribed characteristics. The transient began with bright ejecta moving up pre-existing rays and weaker, more diffuse, mass injection in the area between rays. At the sudden onset of this mass ejection the rays began to move apart from one another: during the ensuing 22 hr the separation of one pair of rays increased by 37 deg of latitude. Towards the end of the event these moving rays looked much like the legs of previously described loop transients. A type I radio source lay close to the center of symmetry of the transient.

Duncan, R. A.

The mechanisms of flame holding in the wake of a bluff body

The flame holding mechanism for lean methane- and lean propane-air flames is examined under conditions where the recirculation zone is absent. The main objective of this work is to study the holding process in detail in an attempt to determine the mechanism of flame holding and also the conditions where this mechanism is viable and when it fails and blow-off occurs. Inverted flames held in the wake of a flat strip were studied. Experiments with different sizes of flame holders were performed. The velocity flow field was determined using a laser Doppler velocimetry technique. Equation of continuity was used to calculate the flame temperature from the change in area of flow streamlines before and after the flame. Observations of the inverted flame itself were obtained using schlieren and direct photography. Results show that there are different mechanisms operative at the time of blow-off for lean propane and methane flames. Blow-off or extinction occurs for lean propane-air flame in spite of the reaction going to completion and the disparity between the heat loss and the gain in mass diffusion in the reaction zone i.e., Le 1.0 causes the flame to blow-off. For methane-air flame the controlling factor or blow-off is incomplete reaction due to higher blowing rate leading to reduced residence time in the reaction zone.

Strehlow, R. A.

Solution Growth of Crystals in Zero-Gravity

In a low-g environment, buoyancy driven convection effects in solution crystal growth are greatly reduced and, thus, one can study diffusion mass transport which in 1-g is masked by convective phenomena. Crystals of triglycine sulfate (TGS) will be grown aboard the Spacelab 3 mission, using a specially developed Fluids Experiments System (FES). The objectives of the experiment are: (1) to develop a technique for solution crystal growth in a low-g environment, (2) to characterize the growth environment provided by an orbiting spacecraft and to determine the influence of the environment on the growth behavior, and (3) to determine how gravity in a low-gravity environment influences the properties of a resulting TGS crystal. Single crystals of TGS have been grown using conventional low-temperature solution crystal growth method and the growth process has been extensively characterized. Various physical properties of TGS solution have been measured. Also, a unique technique of growing solution growth crystals by extracting heat at a programmed rate from the crystal through a semi-insulating sting has been developed and tested in 1-g environment.

Lal, R. B.

Results and further experiments using Spacelab holography

Spacelab 3 provided the platform for the first extensive application of holography in space. Holograms of excellent quality were recorded during the flight, providing information on the solution growth of triglycine sulfate crystals under microgravity conditions. The suppression of convection and the development of a symmetric concentration-depletion region around each crystal were clearly shown. Crystal growth was governed by diffusion mass transport across the depletion region. The successful performance of the holographic unit has prompted NASA to schedule a slightly modified version for reflight on IML-1; opportunities now exist for other investigators to use this facility.

Owen, R. B.

Numerical Modeling Of Two-Phase, Reactive Flows

Capabilities of several computational schemes combined and enhanced. Advanced Rocket Injector-Combustor Code is computer program for numerical simulation of two-phase flows of chemically reacting or nonreacting components. Applied to sprays of fuel in injectors and combustors, flows in heat exchangers, evaporators, and mass-diffusion systems; and two-dimensional or axisymmetric flows containing liquid or solid suspensions. Used in research for development and verification of various physical models.

Liang, Pak-Yan

Catalytic ignition model in a monolithic reactor with in-depth reaction

Two transient models have been developed to study the catalytic ignition in a monolithic catalytic reactor. The special feature in these models is the inclusion of thermal and species structures in the porous catalytic layer. There are many time scales involved in the catalytic ignition problem, and these two models are developed with different time scales. In the full transient model, the equations are non-dimensionalized by the shortest time scale (mass diffusion across the catalytic layer). It is therefore accurate but is computationally costly. In the energy-integral model, only the slowest process (solid heat-up) is taken as nonsteady. It is approximate but computationally efficient. In the computations performed, the catalyst is platinum and the reactants are rich mixtures of hydrogen and oxygen. One-step global chemical reaction rates are used for both gas-phase homogeneous reaction and catalytic heterogeneous reaction. The computed results reveal the transient ignition processes in detail, including the structure variation with time in the reactive catalytic layer. An ignition map using reactor length and catalyst loading is constructed. The comparison of computed results between the two transient models verifies the applicability of the energy-integral model when the time is greater than the second largest time scale of the system. It also suggests that a proper combined use of the two models can catch all the transient phenomena while minimizing the computational cost.

Tien, Ta-Ching