Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion equation”

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,477 records · Page 82

Dissolution of multicomponent bubbles

The behavior of an isolated, stationary, multicomponent gas bubble in a glassmelt containing several dissolved gases is considered. The relevant mass-transport equations are formulated and calculations are performed for the case of two diffusing gases using a quasi-stationary model and a numerical solution of the exact mass-transfer equations. The results obtained from these two approaches are compared. The factors which govern the dissolution or growth of a bubble are thermodynamic and kinetic in origin. The tendency of a bubble to grow or shrink at long times is controlled by departure from overall equilibrium, whereas the short-time bubble dynamics may be dominated by kinetic effects. As a result of the existence of these dual influences, maxima and/or minima occur in the functional dependence of the bubble radius on time.

Weinberg, M. C.↗

Multi-Scale Modeling of Magnetospheric Reconnection

One of the major challenges in modeling the magnetospheric magnetic reconnection is to quantify the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions near reconnection sites. There is still considerable debate as to what degree microphysical processes on kinetic scales affect the global evolution and how important it is to substitute numerical dissipation and/or ad hoc anomalous resistivity by a physically motivated model of dissipation. Comparative studies of magnetic reconnection in small scale geometries demonstrated that MHD simulations that included non-ideal processes in terms of a resistive term $\eta J$ did not produce the fast reconnection rates observed in kinetic simulations. For a broad range of physical parameters in collisionless magnetospheric plasma, the primary mechanism controlling the dissipation in the vicinity of the reconnection site is non-gyrotropic effects with spatial scales comparable with the particle Larmor radius. We utilize the global MHD code BATSRUS and incorporate nongyrotropic effects in diffusion regions in terms of corrections to the induction equation. We developed an algorithm to search for magnetotail reconnection sites, specifically where the magnetic field components perpendicular to the local current direction approaches zero and form an X-type configuration. Spatial scales of the diffusion region and magnitude of the reconnection electric field are calculated selfconsistently using MHD plasma and field parameters in the vicinity of the reconnection site. The location of the reconnection sites is updated during the simulations. To clarify the role of nongyrotropic effects in diffusion region on the global magnetospheric dynamic we perform simulations with steady southward IMF driving of the magnetosphere. Ideal MHD simulations with magnetic reconnection supported by numerical resistivity produce steady configuration with almost stationary near-earth neutral line (NENL). Simulations with non-gyrotropic corrections demonstrate dynamic quasi-periodic response to the steady driving condition. The loading/unloading cycle in non-gyrotropic MHD results has a non-stationary reconnection site in the magnetotail, with the retreating during the stretching phase and then a new NENL forming in the resulting thin plasma sheet. We expect that this model will lead to improved representations of space weather event in the magnetosphere.

Kuznetsova, M. M.↗

Spiking Neurons for Analysis of Patterns

Artificial neural networks comprising spiking neurons of a novel type have been conceived as improved pattern-analysis and pattern-recognition computational systems. These neurons are represented by a mathematical model denoted the state-variable model (SVM), which among other things, exploits a computational parallelism inherent in spiking-neuron geometry. Networks of SVM neurons offer advantages of speed and computational efficiency, relative to traditional artificial neural networks. The SVM also overcomes some of the limitations of prior spiking-neuron models. There are numerous potential pattern-recognition, tracking, and data-reduction (data preprocessing) applications for these SVM neural networks on Earth and in exploration of remote planets. Spiking neurons imitate biological neurons more closely than do the neurons of traditional artificial neural networks. A spiking neuron includes a central cell body (soma) surrounded by a tree-like interconnection network (dendrites). Spiking neurons are so named because they generate trains of output pulses (spikes) in response to inputs received from sensors or from other neurons. They gain their speed advantage over traditional neural networks by using the timing of individual spikes for computation, whereas traditional artificial neurons use averages of activity levels over time. Moreover, spiking neurons use the delays inherent in dendritic processing in order to efficiently encode the information content of incoming signals. Because traditional artificial neurons fail to capture this encoding, they have less processing capability, and so it is necessary to use more gates when implementing traditional artificial neurons in electronic circuitry. Such higher-order functions as dynamic tasking are effected by use of pools (collections) of spiking neurons interconnected by spike-transmitting fibers. The SVM includes adaptive thresholds and submodels of transport of ions (in imitation of such transport in biological neurons). These features enable the neurons to adapt their responses to high-rate inputs from sensors, and to adapt their firing thresholds to mitigate noise or effects of potential sensor failure. The mathematical derivation of the SVM starts from a prior model, known in the art as the point soma model, which captures all of the salient properties of neuronal response while keeping the computational cost low. The point-soma latency time is modified to be an exponentially decaying function of the strength of the applied potential. Choosing computational efficiency over biological fidelity, the dendrites surrounding a neuron are represented by simplified compartmental submodels and there are no dendritic spines. Updates to the dendritic potential, calcium-ion concentrations and conductances, and potassium-ion conductances are done by use of equations similar to those of the point soma. Diffusion processes in dendrites are modeled by averaging among nearest-neighbor compartments. Inputs to each of the dendritic compartments come from sensors. Alternatively or in addition, when an affected neuron is part of a pool, inputs can come from other spiking neurons. At present, SVM neural networks are implemented by computational simulation, using algorithms that encode the SVM and its submodels. However, it should be possible to implement these neural networks in hardware: The differential equations for the dendritic and cellular processes in the SVM model of spiking neurons map to equivalent circuits that can be implemented directly in analog very-large-scale integrated (VLSI) circuits.

Huntsberger, Terrance↗

Computation of three-dimensional viscous flows using a space-marching method

A space-marching method, developed to compute three-dimensional flows for internal geometries, has been utilized to predict viscous flows through a curved duct and over a swept wing. The Navier-Stokes equations have been posed as an initial value problem by neglecting the streamwise viscous diffusion terms and by treating the pressure gradient as a known source term. The resulting equations have been solved by a non-iterative (single pass) algorithm at each streamwise step. The results are compared with earlier computations (based on iterative methods) and the experimental data. The agreement between the present predictions, the experimental data, and the earlier predictions is good for the cases computed. The computation time is only a fraction of the iterative methods.

Murthy, K. N. S.↗

TWO-DIMENSIONAL CORE-COLLAPSE SUPERNOVA SIMULATIONS WITH THE ISOTROPIC DIFFUSION SOURCE APPROXIMATION FOR NEUTRINO TRANSPORT

The neutrino mechanism of core-collapse supernova is investigated via non-relativistic, two-dimensional (2D), neutrino radiation–hydrodynamic simulations. For the transport of electron flavor neutrinos, we use the interaction rates defined by Bruenn and the isotropic diffusion source approximation (IDSA) scheme, which decomposes the transported particles into trapped-particle and streaming-particle components. Heavy neutrinos are described by a leakage scheme. Unlike the “ray-by-ray” approach in some other multidimensional supernova models, we use cylindrical coordinates and solve the trapped-particle component in multiple dimensions, improving the proto-neutron star resolution and the neutrino transport in angular and temporal directions. We provide an IDSA verification by performing one-dimensional (1D) and 2D simulations with 15 and 20 M{sub ⊙} progenitors from Woosley et al. and discuss the difference between our IDSA results and those existing in the literature. Additionally, we perform Newtonian 1D and 2D simulations from prebounce core collapse to several hundred milliseconds postbounce with 11, 15, 21, and 27 M{sub ⊙} progenitors from Woosley et al. with the HS(DD2) equation of state. General-relativistic effects are neglected. We obtain robust explosions with diagnostic energies E{sub dia} ≳ 0.1–0.5 B (1 B ≡ 10{sup 51} erg) for all considered 2D models within approximately 100–300 ms after bounce and find that explosions are mostly dominated by the neutrino-driven convection, although standing accretion shock instabilities are observed as well. We also find that the level of electron deleptonization during collapse dramatically affects the postbounce evolution, e.g., the neglect of neutrino–electron scattering during collapse will lead to a stronger explosion.

79 ASTRONOMY AND ASTROPHYSICS↗

Improved Prediction of Cold-Air Pools in the Weather Research and Forecasting Model Using a Truly Horizontal Diffusion Scheme for Potential Temperature

The terrain-following vertical coordinate system used by many atmospheric models, including the Weather Research and Forecasting (WRF) Model, is prone to errors in regions of complex terrain. These errors stem, in part, from the calculation of horizontal gradients within the diffusion term of the momentum or scalar evolution equations. In WRF, such gradients can be calculated along coordinate surfaces, or using metric terms that help account for grid skewness. However, neither of these options ensures a truly horizontal gradient calculation, especially if a grid cell is skewed enough that the heights of the neighboring grid points used in the calculation fall outside the vertical range of the cell. In this work, an improved scheme that uses Taylor series approximations to vertically interpolate variables to the level necessary for a truly horizontal gradient calculation is implemented in WRF for the diffusion of potential temperature. The scheme is validated using an atmosphere-at-rest configuration, in which spurious flows develop only as a result of numerical errors and can thus be used as a proxy for model performance. Following validation, the method is applied to the simulation of cold-air pools (CAPs), which occur in regions of complex terrain and are characterized by strong near-surface temperature gradients. Using the truly horizontal scheme, idealized simulations demonstrate reduced numerical mixing in a quiescent CAP, and a realistic case study in the Columbia River basin shows a reduction in positive wind speed bias by up to roughly 20% compared to observations from the Second Wind Forecast Improvement Project.

54 ENVIRONMENTAL SCIENCES↗

Self-consistent drift-diffusion-reaction model for the electron beam interaction with dielectric samples

The charging of insulating samples degrades the quality and complicates the interpretation of images in scanning electron microscopy and is important in other applications, such as particle detectors. In this paper, we analyze this nontrivial phenomenon on different time scales employing the drift-diffusion-reaction approach augmented with the trapping rate equations and a realistic semi-empirical source function describing the pulsed nature of the electron beam. We consider both the fast processes following the impact of a single primary electron, the slower dynamics resulting from the continuous bombardment of a sample, and the eventual approach to the steady-state regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quasielastic neutron scattering measurements and ab initio MD-simulations on single ion motions in molten NaF

The ionic stochastic motions in the molten alkali halide NaF are investigated by quasielastic neutron scattering and first principles molecular dynamics simulation. Quasielastic neutron scattering was employed to extract the diffusion behavior of the sodium ions in the melt. An extensive first principles based simulation on a box of up to 512 particles has been performed to complement the experimental data. From that large box, a smaller 64-particle box has then been simulated over a runtime of 60 ps. A good agreement between calculated and neutron data on the level of spectral shape has been obtained. The obtained sodium diffusion coefficients agree very well. The simulation predicts a fluorine diffusion coefficient similar to the sodium one. Applying the Nernst-Einstein equation, a remarkable large cross correlation between both ions can be deduced. The velocity cross correlations demonstrate a positive correlation between the ions over a period of 0.1 ps. That strong correlation is evidence that the unlike ions do not move completely statistically independent and have a strong association over a short period of time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An elliptic representation of coupled boundary layers and inviscid core for computation of separated internal flows

This report presents a viscous-inviscid calculation method for the prediction of turbulent incompressible flows in diffusers with small regions of stall. Integral turbulent boundary layer equations (BLE) are used to provide boundary conditions for a finite difference representation of the inviscid core. The BLE are applied at the floating displacement-thickness line, resulting in a simulataneous set of nonlinear block tridiagonal equations which are solved iteratively using a sucessive line-relaxation technique. Results of the method have been applied to calculate the performance of a variety of asymmetric stalled diffusers.

Ghose, S.↗

Simulations of surface winds at the Viking Lander sites using a one-level model

The one-level model developed by Mass and Dempsey for use in predicting surface flows in regions of complex terrain was adapted to simulate surface flows at the Viking lander sites on Mars. In the one-level model, prediction equations for surface winds and temperatures are formulated and solved. Surface temperatures change with time in response to diabatic heating, horizontal advection, adiabatic heating and cooling effects, and horizontal diffusion. Surface winds can change in response to horizontal advection, pressure gradient forces, Coriolis forces, surface drag, and horizontal diffusion. Surface pressures are determined by integration of the hydrostatic equation from the surface to some reference level. The model has successfully simulated surface flows under a variety of conditions in complex-terrain regions on Earth.

Bridger, Alison F. C.↗

DISPERSAL OF G-BAND BRIGHT POINTS AT DIFFERENT LONGITUDINAL MAGNETIC FIELD STRENGTHS

G-band bright points (GBPs) are thought to be the foot-points of magnetic flux tubes. The aim of this paper is to investigate the relation between the diffusion regimes of GBPs and the associated longitudinal magnetic field strengths. Two high resolution observations of different magnetized environments were acquired with the Hinode/Solar Optical Telescope. Each observation was recorded simultaneously with G-band filtergrams and Narrow-band Filter Imager (NFI) Stokes I and V images. GBPs are identified and tracked automatically, and then categorized into several groups by their longitudinal magnetic field strengths, which are extracted from the calibrated NFI magnetograms using a point-by-point method. The Lagrangian approach and the distribution of diffusion indices approach are adopted separately to explore the diffusion regime of GBPs for each group. It is found that the values of diffusion index and diffusion coefficient both decrease exponentially with the increasing longitudinal magnetic field strengths whichever approach is used. The empirical formulas deduced from the fitting equations are proposed to describe these relations. Stronger elements tend to diffuse more slowly than weak elements, independently of the magnetic flux of the surrounding medium. This may be because the magnetic energy of stronger elements is not negligible compared with the kinetic energy of the gas, and therefore the flows cannot perturb them so easily.

79 ASTRONOMY AND ASTROPHYSICS↗

Studies on structural and electrical properties of nanostructured GdMnO{sub 3}

In this communication, we report the results of the studies on structural and electrical properties of cost effective sol-gel grown nanostructured GdMnO{sub 3} manganites. Samples were sintered at different temperatures. Structural phases have been analyzed by performing the X-ray diffraction (XRD) measurements while crystallite size has been calculated using Scherer’s formula. AC conductivity (σ{sub AC}) has been measured for all the samples and fitted theoretically by using power law equation. Variation in σ{sub AC} has been understood on the basis of crystallite boundaries and oxygen ordering in all the samples. Nature of conduction in the samples has been discussed in the light of correlated barrier hopping mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Influence of a Simple Heat Loss Profile on a Pure Diffusion Flame

The presence of soot on the fuel side of a diffusion flame results in significant radiative heat losses. The influence of a fuel side heat loss zone on a pure diffusion flame established between a fuel and an oxidizer wall is investigated by assuming a hypothetical sech(sup 2) heat loss profile. The intensity and width of the loss zone are parametrically varied. The loss zone is placed at different distances from the Burke-Schumann flame location. The migration of the temperature and reactivity peaks are examined for a variety of situations. For certain cases the reaction zone breaks through the loss zone and relocates itself on the fuel side of the loss zone. In all cases the temperature and reactivity peaks move toward the fuel side with increased heat losses. The flame structure reveals that the primary balance for the energy equation is between the reaction term and the diffusion term. Extinction plots are generated for a variety of situations. The heat transfer from the flame to the walls and the radiative fraction is also investigated, and an analytical correlation formula, derived in a previous study, is shown to produce excellent predictions of our numerical results when an O(l) numerical multiplicative constant is employed.

Ray, Anjan↗

The diffusion of stars through phase space

An orbit-averaged Fokker-Planck equation has been derived to study the secular evolution of stellar systems with regular orbits and the heating of stellar disks. It is shown that a population of stars with an initially Maxwellian peculiar-velocity distribution will remain Maxwellian as it diffuses through orbit space only if: (1) a second-order diffusion tensor is proportional to epicycle energy; and (2) the population's velocity dispersion grows as the square root of time. Scattering by ephemeral spiral waves is able to account for the observed kinematics of the solar neighborhood only if the waves have wavelengths in excess of 9 kpc and constantly drifting pattern speeds.

Binney, James↗

Dissipative particle dynamics of diffusion-NMR requires high Schmidt-numbers

We present an efficient mesoscale model to simulate the diffusion measurement with nuclear magnetic resonance (NMR). On the level of mesoscopic thermal motion of fluid particles, we couple the Bloch equations with dissipative particle dynamics (DPD). Thereby we establish a physically consistent scaling relation between the diffusion constant measured for DPD-particles and the diffusion constant of a real fluid. The latter is based on a splitting into a centre-of-mass contribution represented by DPD, and an internal contribution which is not resolved in the DPD-level of description. As a consequence, simulating the centre-of-mass contribution with DPD requires high Schmidt numbers. After a verification for fundamental pulse sequences, we apply the NMR-DPD method to NMR diffusion measurements of anisotropic fluids, and of fluids restricted by walls of microfluidic channels. For the latter, the free diffusion and the localisation regime are considered.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Probability and Cumulative Density Function Methods for the Stochastic Advection-Reaction Equation

We present a cumulative density function (CDF) method for the probabilistic analysis of $d$-dimensional advection-dominated reactive transport in heterogeneous media. We employ a probabilistic approach in which epistemic uncertainty on the spatial heterogeneity of Darcy-scale transport coefficients is modeled in terms of random fields with given correlation structures. Our proposed CDF method employs a modified Large-Eddy-Diffusivity (LED) approach to close and localize the nonlocal equations governing the one-point PDF and CDF of the concentration field, resulting in a $(d + 1)$ dimensional PDE. Compared to the classsical LED localization, the proposed modified LED localization explicitly accounts for the mean-field advective dynamics over the phase space of the PDF and CDF. To illustrate the accuracy of the proposed closure, we apply our CDF method to one-dimensional single-species reactive transport with uncertain, heterogeneous advection velocities and reaction rates modeled as random fields.

advection-reaction equation↗

On the surface photovoltaic effect in a multivalley semiconductor in an external magnetic field

A theory of the photovoltaic effect in a semi-infinite multvalley semiconductor upon the absorption of polarized light at free carriers caused by the specular reflection and diffuse scattering of electrons at a film surface is developed by us. The kinetic Boltzmann equation in the approximation of time relaxation and boundary conditions, which determine the correlation between the distribution function of electrons reflected from a semi-infinite crystal surface and the distribution function of electrons incident onto a surface for the case of both specular reflection and diffuse scattering is used. The aforementioned takes into account the fact that the distribution function of electrons diffusely scattered from a surface depends only on their energy and is determined from the condition of total electron flux vanishing at the surface. Expressions for analyzing the spectral dependence of the current which is linearly dependent on the magnetic-field strength are obtained.

36 MATERIALS SCIENCE↗

Langevin and Fokker-Planck analyses for diffusion-mediated passing of circular and discorectangular species in two-dimensional channels

The propensity for pairs of diffusing species to pass each other within narrow channels or pores is of basic interest as a first-passage-type problem. It is also of relevance for solution-phase transport in nanoporous materials, and in particular for catalytic conversion reactions where high yield requires that product species can efficiently pass reactant species to exit the pores. Here, we analyze a two-dimensional model with nonoverlapping circular and discorectangular species confined to a rectangular channel, and where passing is mediated by Brownian dynamics in an implicit solvent. For narrower channels where passing is still possible, the discorectangle must align with the channel to pass the circular species. Behavior of the passing propensity, 𝑃, can be assessed by strongly damped Langevin simulations, or within an equivalent Fokker-Planck equation (FPE) formalism. The latter corresponds to a diffusion problem in a “higher-dimensional channel” with a constriction. We assess the variation of the passing propensity, 𝑃, for a broad range of channel width including its scaling just above the threshold where passing is sterically blocked. Analysis of 𝑃 versus the rotational diffusion coefficient 𝐷 𝑟 of the discorectangle reveals a significant decrease in 𝑃 for lower 𝐷 𝑟 for moderate channel width. This prompts a direct analysis of the regime where 𝐷 𝑟 → 0 , for which the FPE can be reduced to a three-dimensional diffusion problem, precise analysis of which is facilitated by adaptive-mesh finite element methods. The dependence of 𝑃 on the aspect ratio of the discorectangle is also assessed.

Rahman, Md Khaledur [Ames Laboratory (AMES), Ames,↗