Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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 1,027 records · Page 57

Diffusion Of Mass In Evaporating Multicomponent Drops

Report summarizes study of diffusion of mass and related phenomena occurring in evaporation of dense and dilute clusters of drops of multicomponent liquids intended to represent fuels as oil, kerosene, and gasoline. Cluster represented by simplified mathematical model, including global conservation equations for entire cluster and conditions on boundary between cluster and ambient gas. Differential equations of model integrated numerically. One of series of reports by same authors discussing evaporation and combustion of sprayed liquid fuels.

Bellan, Josette↗

Modeling of turbulent chemical reaction

Viewgraphs are presented on modeling turbulent reacting flows, regimes of turbulent combustion, regimes of premixed and regimes of non-premixed turbulent combustion, chemical closure models, flamelet model, conditional moment closure (CMC), NO(x) emissions from turbulent H2 jet flames, probability density function (PDF), departures from chemical equilibrium, mixing models for PDF methods, comparison of predicted and measured H2O mass fractions in turbulent nonpremixed jet flames, experimental evidence of preferential diffusion in turbulent jet flames, and computation of turbulent reacting flows.

Chen, J.-Y.↗

The Solar-Sail Launched Interstellar Probe: Pre-Perihelion Trajectories and Application of Holography

Design of missions beyond our solar system presents many challenges. Here, we consider certain aspects of the solar-sail launched interstellar probe (ISP), a spacecraft slated for launch in the 2010 time period that is planned to reach the heliopause, at 200 Astronomical Units (AU) from the Sun after a flight of about 20-years duration. The baseline mission under consideration by NASA / JPL has a sail radius of 200 m, a science payload of 25 kg, a spacecraft areal mass thickness of about two grams per square meter and is accelerated out of the solar system at about 14 AU per year after performing a perihelion pass of about 0.25 AU. In current plans, the sail is to be dropped near Jupiter's orbit (5.2 AU from the Sun) on the outbound trajectory leg. One aspect of this study is application of a realistic model of sail thermo-optics to sail kinematics that includes diffuse / specular reflectance and sail roughness. The effects of solar-wind degradation of sail material, based on recent measurements at the NASA MSFC (Marshall Space Flight Center) Space Environment Facility were incorporated in the kinematical model. After setting initial and final conditions for the spacecraft, trajectory was optimized using the provision of variable sail aspect angle. The second phase of the study included consideration of rainbow holography as a medium for a message plaque that would be carried aboard the ISP in the spirit of the message plaques aboard Pioneer 10 /11 and Voyager 1 /2. A prototype holographic message plaque was designed and created by artist C. Bangs with the assistance of Ana Maria Nicholson and Dan Schweitzer of the Center for Holographic Arts in Long Island City, NY. The piece was framed by Simon Liu Inc. of Brooklyn, NY. Concurrent to the creation of the prototype message plaque, we explored the potential of this medium to transmit large amounts of visual information to any extraterrestrial civilization that might detect and intercept ISP. It was also necessary to investigate possible degradation of holograms by the space environment. We developed a new way of characterizing the optical quality of holograms.

Gregory L. Matloff↗

Characterization and Computational Modeling of Minor Phases in Alloy LSHR

The minor phases of powder metallurgy disk superalloy LSHR were studied. Samples were consistently heat treated at three different temperatures for long times to approach equilibrium. Additional heat treatments were also performed for shorter times, to assess minor phase kinetics in non-equilibrium conditions. Minor phases including MC carbides, M23C6 carbides, M3B2 borides, and sigma were identified. Their average sizes and total area fractions were determined. CALPHAD thermodynamics databases and PrecipiCalc(TradeMark), a computational precipitation modeling tool, were employed with Ni-base thermodynamics and diffusion databases to model and simulate the phase microstructural evolution observed in the experiments with an objective to identify the model limitations and the directions of model enhancement.

Jou, Herng-Jeng↗

Sodium Pumping via Condensation within a Non-Wetting Porous Structure

The transfer of liquid sodium from a low-pressure condenser to a high-pressure evaporator is required to complete the thermodynamic cycle of a sodium thermal electrochemical converter (Na-TEC). A unique capillary pump is explored in this work, whereby low-pressure sodium vapor is condensed within a non-wetting (i.e. contact angle > pi/2) stainless steel porous structure. Due to the curvature of the non-wetting liquid-vapor interface, the liquid adjacent to the interface is at a higher pressure than the vapor. This is in contrast to traditional wicks, where the liquid adjacent to the interface has a lower pressure. Conjugate heat transfer techniques are used to model the coupled momentum and thermal transport processes within this non-wetting porous structure. Due to the low mass flowrates in the Na-TEC (< 0.1 mg/s), higher order friction terms are neglected and the classical Darcy’s law is used for the non-isothermal, compressible vapor flow. At lower temperatures, the sodium vapor density becomes sufficiently small and there is a transition from classic Poiseuille flow (Kn < 0.1) to Knudsen flow (Kn > 10). Furthermore, transport due to binary diffusion is also considered for situations where there is a substantial presence of non-condensable argon in the sodium vapor plenum. Initial modeling results showing expected mass flowrates, temperature profiles, and liquid-vapor interface locations are presented for various physical conditions. The main design features of an experimental set-up used to demonstrate this capillary pumping is also described in detail, including the method used to measure the low sodium mass flowrates. Initial experimental results demonstrating sodium pumping for a range of temperatures (800K – 1100K) will be presented and compared against the model predictions. Finally, several of the challenges with this experiment will be discussed, and specific strategies used to address these are presented.

14 SOLAR ENERGY↗

Numerical analysis of flows of rarefied gases in long channels with octagonal cross section shapes

Isothermal, pressure driven rarefied gas flows through long channels with octagonal cross section shapes are analyzed computationally. The capillary is between inlet and outlet reservoirs. The cross section is constant along the axial direction. The boundary condition at the solid-gas interface is assumed to be diffuse reflection. Since the channel is long, the gaseous velocity is small compared to the average molecular speed. Consequently, a linearized description can be used. The flow is described by the linearized Bhatnagar-Gross-Krook kinetic model. The solution of the problem is divided into two stages. First, the local flow field is determined by assuming the local pressure gradient. Secondly, the global flow behavior is deduced by the consideration of the conservation of the mass along the axis of the capillary. The kinetic equation is solved by the discrete velocity method on the cross section. Both spatial and velocity spaces are discretized. A body fitted rectangular grid is used for the spatial space. Near the boundary, first-order, while in the interior part of the flow domain, second-order finite-differences are applied to approximate the spatial derivatives. This combination results into an efficient and straightforward numerical treatment. The velocity space is represented by a Gauss-Legendre quadrature. The kinetic equation is solved in an iterative manner. The local dimensionless flow rate is calculated and tabulated for a wide range of the gaseous rarefaction for octagonal cross sections with various geometrical parameters. It exhibits the Knudsen minimum phenomenon. The flow rates in the octagonal channel are compared to those through capillaries with circular and square cross sections. Typical velocity profiles are also shown. The mass flow rate and the distribution of the pressure are determined and presented for global pressure driven flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Equilibration in fermionic systems

Highlights: •Local equilibration in finite fermion systems is based on a nonlinear diffusion equation. •The fermionic diffusion equation is solved analytically in closed form. •The exact solutions agree with the numerical results. •Particle-number and energy conservation are shown to be fulfilled. •The model is applied to hadrons at low energies, and relativistic quarks. -- Abstract: The time evolution of a finite fermion system towards local statistical equilibrium is investigated using analytical solutions of a nonlinear partial differential equation that had been derived earlier from the Boltzmann collision term. The solutions of this fermionic diffusion equation are rederived in closed form, evaluated exactly for simplified initial conditions, and applied to hadron systems at low energies in the MeV-range, as well as to quark systems at relativistic energies in the TeV-range where antiparticle production is abundant. Conservation laws for particle number including created antiparticles, and for the energy are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A theoretical analysis of steady-state charge collection in simple diodes under high-injection conditions

A previous rigorous mathematical analysis of drift-diffusion equations was used to investigate collected charge in a simple reverse-biased p-n junction diode exposed to an ionization source that liberates carriers (electron-hole pairs) in a quasi-neutral region within the diode. Each of two simple models was found to agree with the more rigorous analysis when carrier liberation is sufficiently intense. One is the sensitive volume (SV) model, and the other was called “ambipolar diffusion with a cutoff” (ADC). The earlier rigorous analysis was worked out in detail only for a localized source, i.e., a point source of carrier liberation, so it was able to validate the applicability of each simple model only for that case. The present paper treats an arbitrary spatial distribution of carrier generation and concludes that the ADC model remains valid for this more general case, but the SV model does

Edmonds, Larry D.↗

Ni coarsening under humid atmosphere in the electrode of solid oxide cells: A combined study of density-functional theory and phase-field modeling

Ni coarsening in solid oxide cell (SOC) electrodes is known to be significantly faster under a humid atmosphere. The underlying mechanisms, though, have not been fully understood. In this work, we examine the surface diffusion of Ni(OH) x by a combination of density-functional theory and phase-field modeling. Here, density-functional theory is used to evaluate the adsorption and diffusion of the Ni(OH) x species on Ni (111) surface, and the results are used to obtain the effective surface diffusivity of Ni versus steam and hydrogen gas partial pressures. This diffusivity is used as an input for a phase-field model to investigate the Ni coarsening in SOC electrodes. It is found that Ni(OH) x formed through chemical reaction cannot accelerate coarsening unless an extremely high steam to hydrogen ratio is reached. However, if Ni(OH) x is formed through electrochemical reactions near triple phase boundaries (TPBs), surface diffusion of Ni(OH) x may cause faster Ni coarsening in fuel cell mode under a large overpotential. Specifically, surface diffusion of Ni(OH) may be comparable to or faster than that of Ni under an overpotential that is large but still possible under fuel cell operating conditions.

25 ENERGY STORAGE↗

An isogeometric variational multiscale method for large-eddy simulation of coupled multi-ion transport in turbulent flow

Highlights: •We present a computational method for coupled multi-ion transport in turbulent flow. •The underlying formulation is a variational multiscale finite element method. •It is combined with the isogeometric concept for electrochemical systems. •Coupled multi-ion transport in fully turbulent Taylor–Couette flow is simulated. •This example is an important model problem for rotating cylinder electrodes. -- Abstract: Electrochemical processes, such as electroplating of large items in galvanic baths, are often coupled to turbulent flow. In this study, we propose an isogeometric residual-based variational multiscale finite element method for multi-ion transport in dilute electrolyte solutions under turbulent flow conditions. In other words, this means that the concepts of isogeometric discretization and variational multiscale methods are successfully combined for developing a method capable of simulating the challenging problem of coupled multi-ion transport in turbulent flow. We present a comprehensive three-dimensional computational method taking into account, among others, coupled convection–diffusion-migration equations subject to an electroneutrality constraint in combination with phenomenological electrode-kinetics modeling. The electrochemical subproblem is one-way coupled to turbulent incompressible flow via convection. Ionic mass transfer in turbulent Taylor–Couette flow is investigated, representing an important model problem for rotating-cylinder-electrode configurations. Multi-ion transport as considered here is an example for mass transport at high Schmidt number (Sc=1389). An isogeometric discretization is especially advantageous for the present problem, since (i) curved boundaries can be represented exactly, and (ii) it has been proven to provide very accurate solutions for flow quantities when being applied in combination with residual-based variational multiscale modeling. We demonstrate that the method is robust and provides results which are in good agreement with direct numerical simulation results as well as empirical mass-transfer correlations reported in literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Adaptive Sampling of Large Deviations

We introduce and test an algorithm that adaptively estimates large deviation functions characterizing the fluctuations of additive functionals of Markov processes in the long-time limit. These functions play an important role for predicting the probability and pathways of rare events in stochastic processes, as well as for understanding the physics of nonequilibrium systems driven in steady states by external forces and reservoirs. The algorithm uses methods from risk-sensitive and feedback control to estimate from a single trajectory a new process, called the driven process, known to be efficient for importance sampling. Its advantages compared to other simulation techniques, such as splitting or cloning, are discussed and illustrated with simple equilibrium and nonequilibrium diffusion models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Wall loss of atomic nitrogen determined by ionization threshold mass spectrometry

In the afterglow of an inductively coupled N{sub 2} plasma, relative N atom densities are measured by ionization threshold mass spectrometry as a function of time in order to determine the wall loss time t{sub wN} from the exponential decay curves. The procedure is performed with two mass spectrometers on different positions in the plasma chamber. t{sub wN} is determined for various pressures, i.e., for 3.0, 5.0, 7.5, and 10 Pa. For this conditions also the internal plasma parameters electron density n{sub e} and electron temperature T{sub e} are determined with the Langmuir probe and the rotational temperature T{sub rot}{sup N{sub 2}} of N{sub 2} is determined with the optical emission spectroscopy. For T{sub rot}{sup N{sub 2}}, a procedure is presented to evaluate the spectrum of the transition υ{sup ′}=0→υ{sup ″}=2 of the second positive system (C{sup 3}Π{sub u}→B{sup 3}Π{sub g}) of N{sub 2}. With this method, a gas temperature of 610 K is determined. For both mass spectrometers, an increase of the wall loss times of atomic nitrogen with increasing pressure is observed. The wall loss time measured with the first mass spectrometer in the radial center of the cylindrical plasma vessel increases linearly from 0.31 ms for 3 Pa to 0.82 ms for 10 Pa. The wall loss time measured with the second mass spectrometer (further away from the discharge) is about 4 times higher. A model is applied to describe the measured t{sub wN.} The main loss mechanism of atomic nitrogen for the considered pressure is diffusion to the wall. The surface loss probability β{sub N} of atomic nitrogen on stainless steel was derived from t{sub wN} and is found to be 1 for the present conditions. The difference in wall loss times measured with the mass spectrometers on different positions in the plasma chamber is attributed to the different diffusion lengths.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The SMART Theory and Modeling Team: An Integrated Element of Mission Development and Science Analysis

When targeting physical understanding of space plasmas, our focus is gradually shifting away from discovery-type investigations to missions and studies that address our basic understanding of processes we know to be important. For these studies, theory and models provide physical predictions that need to be verified or falsified by empirical evidence. Within this paradigm, a tight integration between theory, modeling, and space flight mission design and execution is essential. NASA's Magnetospheric MultiScale (MMS) mission is a pathfinder in this new era of space research. The prime objective of MMS is to understand magnetic reconnection, arguably the most fundamental of plasma processes. In particular, MMS targets the microphysical processes, which permit magnetic reconnection to operate in the collisionless plasmas that permeate space and astrophysical systems. More specifically, MMS will provide closure to such elemental questions as how particles become demagnetized in the reconnection diffusion region, which effects determine the reconnection rate, and how reconnection is coupled to environmental conditions such as magnetic shear angles. Solutions to these problems have remained elusive in past and present spacecraft missions primarily due to instrumental limitations - yet they are fundamental to the large-scale dynamics of collisionless plasmas. Owing to the lack of measurements, most of our present knowledge of these processes is based on results from modern theory and modeling studies of the reconnection process. Proper design and execution of a mission targeting magnetic reconnection should include this knowledge and have to ensure that all relevant scales and effects can be resolved by mission measurements. The SMART mission has responded to this need through a tight integration between instrument and theory and modeling teams. Input from theory and modeling is fed into all aspects of science mission design, and theory and modeling activities are tailored to SMART needs during mission development and science analysis. In this presentation, we will present an overview of SMART theory and modeling team activities. In particular, we will provide examples of science objectives derived from state-of-the art models, and of recent research results that continue to be utilized in SMART mission development.

Hesse, Michael↗

Leaching behavior of copper from waste printed circuit boards with Brønsted acidic ionic liquid

Highlights: • A Brønsted acidic ILs was used to leach Cu from WPCBs for the first time. • The particle size of WPCBs has significant influence on Cu leaching rate. • Cu leaching rate was higher than 99% under the optimum leaching conditions. • The leaching process can be modeled with shrinking core model, and the E{sub a} was 25.36 kJ/mol. - Abstract: In this work, a Brønsted acidic ionic liquid, 1-butyl-3-methyl-imidazolium hydrogen sulfate ([bmim]HSO{sub 4}), was used to leach copper from waste printed circuit boards (WPCBs, mounted with electronic components) for the first time, and the leaching behavior of copper was discussed in detail. The results showed that after the pre-treatment, the metal distributions were different with the particle size: Cu, Zn and Al increased with the increasing particle size; while Ni, Sn and Pb were in the contrary. And the particle size has significant influence on copper leaching rate. Copper leaching rate was higher than 99%, almost 100%, when 1 g WPCBs powder was leached under the optimum conditions: particle size of 0.1–0.25 mm, 25 mL 80% (v/v) ionic liquid, 10 mL 30% hydrogen peroxide, solid/liquid ratio of 1/25, 70 °C and 2 h. Copper leaching by [bmim]HSO{sub 4} can be modeled with the shrinking core model, controlled by diffusion through a solid product layer, and the kinetic apparent activation energy has been calculated to be 25.36 kJ/mol.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Factors controlling surface oxygen exchange in oxides

Decreasing the working temperature of solid oxide fuel cells is critical to their increased commercialization but is inhibited by the slow oxygen exchange kinetics at the cathode, which limits the overall rate of the oxygen reduction reaction. We use ab initio methods to develop a quantitative elementary reaction model of oxygen exchange in a representative cathode material, La 0.5 Sr 0.5 CoO 3-δ , and predict that under operating conditions the rate-limiting step for oxygen incorporation from O 2 gas on the stable, (001)-SrO surface is lateral (surface) diffusion of O-adatoms and oxygen surface vacancies. We predict that a high vacancy concentration on the metastable CoO 2 termination enables a vacancy-assisted O-2 dissociation that is 10 2 -10 3 times faster than the rate limiting step on the Sr-rich (La,Sr)O termination. This result implies that dramatically enhanced oxygen exchange performance could potentially be obtained by suppressing the (La,Sr)O termination and stabilizing highly active CoO 2 termination.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stratification of drying particle suspensions: Comparison of implicit and explicit solvent simulations

We report that large scale molecular dynamics simulations are used to study drying suspensions of a binary mixture of large and small particles in explicit and implicit solvents. The solvent is first modeled explicitly and then mapped to a uniform viscous medium by matching the diffusion coefficients and the pair correlation functions of the particles. “Small-on-top” stratification of the particles, with an enrichment of the smaller ones at the receding liquid-vapor interface during drying, is observed in both models under the same drying conditions. With the implicit solvent model, we are able to model much thicker films and study the effect of the initial film thickness on the final distribution of particles in the dry film. Laslty, our results show that the degree of stratification is controlled by the Péclet number defined using the initial film thickness as the characteristic length scale. When the Péclet numbers of large and small particles are much larger than 1, the degree of “small-on-top” stratification is first enhanced and then weakens as the Péclet numbers are increased.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The origin of the Stokes–Einstein relation in simple dense liquids

Here, we investigate the origin of the universal relation between structural relaxation and diffusion in simple dense liquids, known as the Stokes–Einstein (SE) relation. The fact that this relation, originally derived from a hydrodynamic model of a macroscopic particle in a viscous medium, can describe the microscopic-scale liquid dynamics still eludes understanding. We introduce a new universal measure of structural relaxation in a system of N identical particles based on an explicit decomposition of the configuration space into N! congruent convex polyhedra. This measure makes it possible to quantify the correlation between two distinct particle configurations in terms of their minimal Euclidean distance, optimized with respect to particle permutations. Using this measure alongside a model of independent random walkers under the single-occupancy constraint, we derive a master equation that quantifies the SE relation. It allows us to demonstrate that the universal relation between structural relaxation and diffusion in simple dense liquids is caused by two conditions: (a) the confinement of the dominant density fluctuations to the first coordination shell, manifested by de Gennes narrowing, and (b) Gaussianity of the diffusion process; the former is shown to be violated in low-density fluids, and the latter is known to be violated in supercooled liquids.

Physics - Condensed matter physics↗

Ultrasonic determination of thermodynamic threshold parameters for irreversible cutaneous burns

In vivo ultrasonic measurements of the depth of conductive cutaneous burns experimentally induced in anesthetized Yorkshire pigs are reported as a function of burn time for the case in which the skin surface temperature is maintained at 100 C. The data are used in the solution of the one-dimensional heat diffusion equation with time-dependent boundary conditions to obtain the threshold temperature and the energy of transformation per unit mass associated with the transition of the tissue from the state of viability to the state of necrosis. The simplicity of the mathematical model and the expediency of the ultrasonic measurements in studies of thermal injury are emphasized.

Cantrell, J. H., Jr.↗