Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

A Unified SHOC-MF Boundary Layer Parameterization in the NASA Goddard Earth Observing System Model

We present results from a novel extension of the boundary layer parameterization in the NASA Goddard Earth Observing System (GEOS) model. An Eddy Diffusivity-Mass Flux (EDMF) multi-plume updraft scheme was previously implemented in GEOS. Here we replace the K-profile eddy diffusivity scheme of Lock et. al. (2000) with the Simplified Higher Order Closure (SHOC; Bogenschutz and Krueger, 2013). The scheme prognoses turbulent kinetic energy, and derives buoyancy flux, cloud fraction and condensate from a trivariate double-gaussian joint PDF. The original SHOC scheme is modified in two significant respects: (i) contributions from the EDMF updrafts are included in estimates of all higher order moments used to determine the PDF; and (ii) the second and third moments of total water are made prognostic variables, with source terms associated with EDMF updrafts. The predicted higher order moments as well as updraft and environmental properties are compared with large eddy simulations in several case studies. Results indicate that many higher order moments are underestimated when the mass flux contribution is neglected.

Nathan Arnold↗

Numerical evaluation of AGR-2 fission product release

The AGR-2 experiment produced normal operation and accident condition fuel performance data for tri-structural isotropic (TRISO) particles with UCO and UO 2 kernels. Data for compacts with no failed particles and for compacts with one or more failed particles are available. To model the fission product diffusion and release from these compacts, it is important to consider the computational mesh as well as the heat conduction and mass diffusion properties of the materials in a TRISO particle. Code and solution verification studies with the Bison fuel performance code were performed to show that the code is computing correct solutions. Bison comparisons to post-irradiation examination data and PARFUME predictions were made for several sets of compacts. Furthermore, these comparisons considered silver, cesium, strontium, and krypton with and without failed particles. Bison predictions closely match those of PARFUME, while both codes generally overpredict post-irradiation examination data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of a Unified SHOC-EDMF PBL Parameterization in the NASA GEOS Model

We present results from a novel extension of the boundary layer parameterization in the NASA Goddard Earth Observing System (GEOS) model. An Eddy Diffusivity-Mass Flux (EDMF) multi-plume updraft scheme was previously implemented in GEOS. Here we replace the K-profile eddy diffusivity scheme of Lock et. al. (2000) with the Simplified Higher Order Closure (SHOC). The scheme prognoses turbulent kinetic energy, and derives buoyancy flux, cloud fraction and condensate from a trivariate double-gaussian joint PDF. The original SHOC scheme is modified in two significant respects: (i) contributions from the EDMF updrafts are included in estimates of all higher order moments used to determine the PDF; and (ii) the second and third moments of total water are made prognostic variables, with source terms associated with EDMF updrafts. The predicted higher order moments as well as updraft and environmental properties are compared with large eddy simulations in several case studies.

SHOC-EDMF↗

The effect of saturated thermal conduction on clouds in a hot plasma

We numerically investigate the internal evolution of multiphase clouds, which are at rest with respect to an ambient, highly ionized medium (HIM) representing the hot component of the circumgalactic medium. Time-dependent saturated thermal conduction and its implications like condensation rates and mixing efficiency are assessed in multiphase clouds. Our simulations are carried out by using the adaptive mesh refinement code FLASH. The model clouds are initially in both hydrostatic and thermal equilibrium and are in pressure balance with the HIM. Thus, they have steep gradients in both temperature and density at the interface to HIM leading to non-negligible thermal conduction. Several physical processes are considered numerically or semi-analytically: thermal conduction, radiative cooling and external heating of gas, self-gravity, mass diffusion, and dissociation of molecules and ionization of atoms. It turns out that saturated thermal conduction triggers a continuous condensation irrespective of cloud mass. Dynamical interactions with ambient HIM all relate to the radial density gradient in the clouds: (1) mass flux due to condensation is the higher the more homogeneous the clouds are; (2) mixing of condensed gas with cloud gas is easier in low-mass clouds, because of their shallower radial density gradient; and thus (3) accreted gas is distributed more efficiently. A distinct and sub-structured transition zone forms at the interface between cloud and HIM, which starts at smaller radii and is much narrower as deduced from analytical theory.

79 ASTRONOMY AND ASTROPHYSICS↗

Mass Changes the Diffusion Coefficient of Particles with Ligand-Receptor Contacts in the Overdamped Limit

Inertia does not generally affect the long-time diffusion of passive overdamped particles in fluids. Yet a model starting from the Langevin equation predicts a surprising property of particles coated with ligands that bind reversibly to surface receptors: heavy particles diffuse more slowly than light ones of the same size. We show this by simulation and by deriving an analytic formula for the mass-dependent diffusion coefficient in the overdamped limit. We estimate the magnitude of this effect for a range of biophysical ligand-receptor systems, and find it is potentially observable for tailored micronscale DNA-coated colloids.

97 MATHEMATICS AND COMPUTING↗

The crucial role of cooling in the making of molecular clouds and stars

The role that velocity or pressure fluctuations in the H I clouds medium can play in initiating compression of sub-Jeans mass diffuse clouds is outlined. The frequently discussed idea is reviewed that substantial overpressures arising in the warm, medium and/or highly supersonic compressions of H I clouds can push sub-Jeans mass clumps to densities where gravity can take over and complete the star formation process. A nonequilibrium model is used to analyze the impact of external disturbances on gas clouds, and estimates of the required disturbance amplitudes are substantially reduced from previous estimates based on equilibrium arguments. The results reveal that a cloud which cools under compression is particularly sensitive to mild disturbances from its environment. Furthermore, the specific energy required to trigger effective compressions in a cooling medium is nearly independent of the cloud's mass. It is proposed that mildly nonlinear disturbances play a primary role in the formation of molecular cloud and stars.

Tohline, Joel E.↗

Temperature and salt controlled tuning of protein clusters

The formation of molecular assemblies in protein solutions is of strong interest both from a fundamental viewpoint and for biomedical applications. While ordered and desired protein assemblies are indispensable for some biological functions, undesired protein condensation can induce serious diseases. As a common cofactor, the presence of salt ions is essential for some biological processes involving proteins, and in aqueous suspensions of proteins can also give rise to complex phase diagrams including homogeneous solutions, large aggregates, and dissolution regimes. Here, we systematically study the cluster formation approaching the phase separation in aqueous solutions of the globular protein BSA as a function of temperature (T), the protein concentration (c p ) and the concentrations of the trivalent salts YCl 3 and LaCl 3 (c s ). As an important complement to structural, i.e. time-averaged, techniques we employ a dynamical technique that can detect clusters even when they are transient on the order of a few nanoseconds. By employing incoherent neutron spectroscopy, we unambiguously determine the short-time self-diffusion of the protein clusters depending on c p , c s and T. We determine the cluster size in terms of effective hydrodynamic radii as manifested by the cluster center-of-mass diffusion coefficients D. For both salts, we find a simple functional form D(c p , c s , T) in the parameter range explored. The calculated inter-particle attraction strength, determined from the microscopic and short-time diffusive properties of the samples, increases with salt concentration and temperature in the regime investigated and can be linked to the macroscopic behavior of the samples.

59 BASIC BIOLOGICAL SCIENCES↗

Effect of growth parameters on compositional variations in directionally solidified HgCdTe alloys

A series of Hg(1-x)Cd(x)Te alloy crystals was grown by directional solidification with compositions ranging from x = 0.2 to x = 0.4. The measured axial compositional profiles were interpreted in terms of an exact numerical solution of the appropriate diffusion equation that takes into account both the variations of the segregation coefficient and solidification rate with composition. The solutions for all growth rates agree generally well with the experimental data, and confirm earlier observations that employed approximate analytical solutions. The effective mass diffusion coefficients showed no significant composition or melt temperature dependence. Slightly higher diffusion coefficients were obtained, however, for the highest growth rates.

Szofran, F. R.↗

Irradiation-enhanced Interactions at UO 2 /Al 2 O 3 /Al Interfaces

Combustion synthesis is used to prepare thin UO 2 films on aluminum alloy substrates. This simple preparation method involves electrospraying uranyl nitrate + acetylacetone + 2-methoxyethanol solution on the substrate, followed by a short annealing at 350 or 550 °C. The irradiation of films with a 40 Ar 2+ ion beam (energy of 1.7 MeV and fluences of 7.7 × 10 16 and 1.3 × 10 17 ions/cm 2 ) is conducted to investigate irradiation-induced restructuring processes. High-resolution transmission electron microscopy (TEM) and X-ray photoelectron spectroscopy (XPS) investigations show that the annealing temperature significantly influences the crystallinity and stability of materials during irradiation. A small amount of Mg in the alloy substrate diffuses into the amorphous Al 2 O 3 interfacial layer between the film and the substrate. Local thermal spikes from the incoming ions facilitate the irradiation-induced mixing of immiscible Al 2 O 3 and UO 2 for the materials prepared at 350 °C. This mass diffusion produces relatively large cavities at the interface. Selective diffusion of a more significant amount of Mg for the materials prepared at 550 °C suppresses the mixing of the Al 2 O 3 interlayer with the film but forms Mg y U 1–y O 2±x solid solutions during irradiation. Local thermal heating triggers the precipitation of a discontinuous crystalline MgO layer close to the film surface. As a result, the enhanced and selective diffusion of Mg into the film makes the materials prepared at 550 °C more robust and mechanically stable than those prepared at 350 °C.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Time-Varying Model for Predicting Formaldehyde Emission Rates in Homes

Recent studies have succeeded in relating emissions of various volatile organic compounds to material mass diffusion transfer using detailed empirical characteristics of each of the individual emitting materials. While significant, the resulting models are often scenario specific and/or require a host of individual component parameters to estimate emission rates. This study developed an approach to estimate aggregated emissions rates based on a wide number of field measurements. We used a multi-parameter regression model based on previous mass transfer models to predict formaldehyde emission rate for a whole dwelling using field-measured, time-resolved formaldehyde concentrations, air exchange rates, and indoor environmental parameters in 63 California single-family houses built between 2011 and 2017. The resulting model provides time-varying formaldehyde emission rates, normalized by floor area, for each study home, assuming a well-mixed mass balance transport model of the home, and a well-mixed layer transport model of indoor surfaces. The surface layer model asserts an equilibrium concentration within the surface layer of the emitted materials that is a function of temperature and RH; the dwelling ventilation rate serves as a surrogate for indoor concentration. We also developed a more generic emission model that is suitable for broad prediction of emission for a population of buildings. This model is also based on measurements aggregated from 27 homes from the same study. We showed that errors in predicting household formaldehyde concentrations using this approach were substantially less than those using a traditional constant emission rate model, despite requiring less unique building information.

54 ENVIRONMENTAL SCIENCES↗

Control of Meridional Flow by a Non-Uniform Rotational Magnetic Field

The diffusive mass transfer of species during crystal growth in vertical ampoules is significantly affected by fluid flow in the liquid mother phase (melt). For electrically conductive melts, an elegant way of remotely inducing and controlling this flow is by utilizing a uniform rotational magnetic field (RMF) in the transverse direction. It induces an azimuthal flow which tends to homogenize the thermal and solutal fields. The rotating field also reduces the diffusion boundary layer, stabilizes temperature fluctuations, and promotes better overall crystal growth. For moderate strengths of the applied magnetic field (2-20 m Tesla) with frequencies of up to 400 Hz, the induced secondary meridional flow becomes significant. It typically consists of one roll at the bottom of the liquid column and a second roll (vortex) at the top. The flow along the centerline (ampoule axis) is directed from the growing solid (interface) towards the liquid (melt). In case of convex interfaces (e.g. in floating zone crystal growth) such flow behavior is beneficial since it suppresses diffusion at the center. However, for concave interfaces (e.g. vertical Bridgman crystal growth) such a flow tends to exacerbate the situation in making the interface shape more concave. It would be beneficial to have some control of this meridional flow- for example, a single recirculating cell with controllable direction and flow magnitude will make this technique even more attractive for crystal growth. Such flow control is a possibility if a non-uniform PNE field is utilized for this purpose. Although this idea has been proposed earlier, it has not been conclusively demonstrated so far. In this work, we derive the governing equations for the fluid dynamics for such a system and obtain solutions for a few important cases. Results from parallel experimental measurements of fluid flow in a mercury column subjected to non-uniform RMF will also be presented.

Mazuruk, Konstantin↗

Direct numerical simulation of a turbulent methane/air flame impinged by a sub-breakdown electric field

Here, the effects of an external electric field on a turbulent methane/air diffusion flame are analyzed in this work using direct numerical simulations. The analyzed configuration consists of a temporally evolving mixing layer of air and a mixture of methane and nitrogen at 1 atm that is impinged by an electric field in the direction normal to the mean mixing plane. The combustion and chemi-ionization reactions involved in the flow are modeled using finite-rate chemistry and a reduced reaction scheme consisting of 26 species and 134 reactions. The mass diffusion and ion-wind effects are modeled using a detailed description of the diffusion coefficients and electric mobilities based on kinetic theory. The presented calculations show that the turbulence generated within the mixing layer is mostly unaffected by the applied electric field for the configuration under exam. In fact, the electric body force is developed away from the mixing region where the flow is uniform. Conversely, the turbulence is able to introduce very high intermittency in the electrically charged species concentration and, consequently, in the electric body force. Such intermittency will constitute a challenge in the future formulation of combustion models that take into account ion-wind effects.

09 BIOMASS FUELS↗

A model of the evaporation of binary-fuel clusters of drops

A formulation has been developed to describe the evaporation of dense or dilute clusters of binary-fuel drops. The binary fuel is assumed to be made of a solute and a solvent whose volatility is much lower than that of the solute. Convective flow effects, inducing a circulatory motion inside the drops, are taken into account, as well as turbulence external to the cluster volume. Results obtained with this model show that, similar to the conclusions for single isolated drops, the evaporation of the volatile is controlled by liquid mass diffusion when the cluster is dilute. In contrast, when the cluster is dense, the evaporation of the volatile is controlled by surface layer stripping, that is, by the regression rate of the drop, which is in fact controlled by the evaporation rate of the solvent. These conclusions are in agreement with existing experimental observations. Parametric studies show that these conclusions remain valid with changes in ambient temperature, initial slip velocity between drops and gas, initial drop size, initial cluster size, initial liquid mass fraction of the solute, and various combinations of solvent and solute. The implications of these results for computationally intensive combustor calculations are discussed.

Harstad, K.↗

Effects of G-Jitter on Interfacial Dynamics of Two Miscible Liquids: Application of MIM

We designed an experiment to examine the effects of g-jitter on mixing of two miscible liquids using the Microgravity Vibration Isolation Mount (MIM). The global bifurcation of the interface was observed with the MIM operating alternatively to either transmit the g-jitter, isolate from the g-jitter or to provide controlled vibration levels with well defined amplitude and frequency content. With the MIM in isolation mode, the interface remains stationary indicating buoyancy induced convection is negligibly small such that mixing occurs via intrinsic mass diffusion without the masking effect of vibration driven convection. Analytical and computational results are in agreement with the experimental findings. Operation of the MIM in forced mode with conditions typical of g-jitter shows that vibration induced convective flows can excite instability mechanisms such as Kelvin-Helmholtz to generate large amplitude quasi-stationary waves oriented vertically for various cases with Stokes-Reynolds number in the range of 0.003 to 0.5. The two and four mode quasi-stationary waves are also predicted with a mathematical model. Though unplanned, the effect of a primary thruster filing was captured and shown to cause a catastrophic bifurcation, enhancing local mass transport. In light of the findings, experiments planned for the International Space Station should consider the potential effects of g-jitter.

Duval, Walter M. B.↗

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.↗

Matrimid substrates with bicontinuous surface and macrovoids in the bulk: A nearly ideal substrate for composite membranes in CO 2 capture

Scalable fabrication of thin-film composite (TFC) membrane for post-combustion carbon capture is often limited by the availability of high-performance nanoporous substrate. Ideally, the substrate should allow for fast gas transport at the selective layer/substrate interface as well as in the bulk of the substrate. In this study, highly permeable Matrimid substrates were prepared via vapor- and nonsolvent-induced phase separations. The phase separation mechanism of the Matrimid®/N-methyl-2-pyrrolidone casting solution was modulated by the addition of LiCl to control the solution thermodynamic stability, and nanoporous Matrimid substrates were formed with a bicontinuous surface and macrovoids in the bulk. This hierarchically optimized structure resulted in a CO 2 permeance of 2.60 × 10 5 GPU (1 GPU = 1 × 10 –6 cm 3 (STP) cm –2 s –1 cmHg –1 = 3.349 × 10 –10 mol m –2 s –1 Pa –1 ), which was ca. 11 times more permeable than a benchmark polyethersulfone substrate with cellular pores. The benefit of using this new substrate was demonstrated by coating a 170-nm amine-containing polymer to form a TFC facilitated transport membrane. The membrane exhibited a CO 2 permeance of 932 GPU at 57°C, which was 72 GPU higher than the counterpart coated on the benchmark substrate. Meanwhile, the CO 2 /N 2 selectivity was remained at 158. This permeance improvement could be well explained by the resistance-in-series model, where the improved permeance was attributed to the reductions in substrate and lateral diffusion mass transfer resistances. The upper bound (UB) analysis indicates that the Matrimid substrate is nearly ideal for the state-of-the-art polymers for CO 2 /N 2 separation. Here, the substrate improvement could also reduce the parasitic energy associated with the membrane process due to the reduced membrane cost and flue gas compression requirement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗