Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion equations”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Derivation of a continuum model and the energy law for moving contact lines with insoluble surfactants

A continuous model is derived for the dynamics of two immiscible fluids with moving contact lines and insoluble surfactants based on thermodynamic principles. The continuum model consists of the Navier-Stokes equations for the dynamics of the two fluids and a convection-diffusion equation for the evolution of the surfactant on the fluid interface. The interface condition, the boundary condition for the slip velocity, and the condition for the dynamic contact angle are derived from the consideration of energy dissipations. Different types of energy dissipations, including the viscous dissipation, the dissipations on the solid wall and at the contact line, as well as the dissipation due to the diffusion of surfactant, are identified from the analysis. A finite element method is developed for the continuum model. Numerical experiments are performed to demonstrate the influence of surfactant on the contact line dynamics. The different types of energy dissipations are compared numerically.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

A Nth-order linear algorithm for extracting diffuse correlation spectroscopy blood flow indices in heterogeneous tissues

Conventional semi-infinite analytical solutions of correlation diffusion equation may lead to errors when calculating blood flow index (BFI) from diffuse correlation spectroscopy (DCS) measurements in tissues with irregular geometries. Very recently, we created an algorithm integrating a Nth-order linear model of autocorrelation function with the Monte Carlo simulation of photon migrations in homogenous tissues with arbitrary geometries for extraction of BFI (i.e., αD{sub B}). The purpose of this study is to extend the capability of the Nth-order linear algorithm for extracting BFI in heterogeneous tissues with arbitrary geometries. The previous linear algorithm was modified to extract BFIs in different types of tissues simultaneously through utilizing DCS data at multiple source-detector separations. We compared the proposed linear algorithm with the semi-infinite homogenous solution in a computer model of adult head with heterogeneous tissue layers of scalp, skull, cerebrospinal fluid, and brain. To test the capability of the linear algorithm for extracting relative changes of cerebral blood flow (rCBF) in deep brain, we assigned ten levels of αD{sub B} in the brain layer with a step decrement of 10% while maintaining αD{sub B} values constant in other layers. Simulation results demonstrate the accuracy (errors < 3%) of high-order (N ≥ 5) linear algorithm in extracting BFIs in different tissue layers and rCBF in deep brain. By contrast, the semi-infinite homogenous solution resulted in substantial errors in rCBF (34.5% ≤ errors ≤ 60.2%) and BFIs in different layers. The Nth-order linear model simplifies data analysis, thus allowing for online data processing and displaying. Future study will test this linear algorithm in heterogeneous tissues with different levels of blood flow variations and noises.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Modeling and sensitivity analysis of mass transfer in active multilayer polymeric film for food applications

The barrier performance of multilayer polymeric films for food applications has been significantly improved by incorporating oxygen scavenging materials. The scavenging activity depends on parameters such as diffusion coefficient, solubility, concentration of scavenger loaded and the number of available reactive sites. These parameters influence the barrier performance of the film in different ways. Virtualization of the process is useful to characterize, design and optimize the barrier performance based on physical configuration of the films. Also, the knowledge of values of parameters is important to predict the performances. Inverse modeling and sensitivity analysis are sole way to find reasonable values of poorly defined, unmeasured parameters and to analyze the most influencing parameters. Thus, the objective of this work was to develop a model to predict barrier properties of multilayer film incorporated with reactive layers and to analyze and characterize their performances. Polymeric film based on three layers of Polyethylene terephthalate (PET), with a core reactive layer, at different thickness configurations was considered in the model. A one dimensional diffusion equation with reaction was solved numerically to predict the concentration of oxygen diffused into the polymer taking into account the reactive ability of the core layer. The model was solved using commercial software for different film layer configurations and sensitivity analysis based on inverse modeling was carried out to understand the effect of physical parameters. The results have shown that the use of sensitivity analysis can provide physical understanding of the parameters which highly affect the gas permeation into the film. Solubility and the number of available reactive sites were the factors mainly influencing the barrier performance of three layered polymeric film. Multilayer films slightly modified the steady transport properties in comparison to net PET, giving a small reduction in the permeability and oxygen transfer rate values. Scavenging capacity of the multilayer film increased linearly with the increase of the reactive layer thickness and the oxygen absorption reaction at short times decreased proportionally with the thickness of the external PET layer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Averaging Principle for the Higher Order Nonlinear Schrödinger Equation with a Random Fast Oscillation

This work concerns the problem associated with averaging principle for a higher order nonlinear Schrödinger equation perturbed by a oscillating term arising as the solution of a stochastic reaction–diffusion equation evolving with respect to the fast time. This model can be translated into a multiscale stochastic partial differential equations. Stochastic averaging principle is a powerful tool for studying qualitative analysis of stochastic dynamical systems with different time-scales. To be more precise, under suitable conditions, we prove that there is a limit process in which the fast varying process is averaged out and the limit process which takes the form of the higher order nonlinear Schrödinger equation is an average with respect to the stationary measure of the fast varying process. Finally, by using the Khasminskii technique we can obtain the rate of strong convergence for the slow component towards the solution of the averaged equation, and as a consequence, the system can be reduced to a single higher order nonlinear Schrödinger equation with a modified coefficient.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

A space-time smooth artificial viscosity method with wavelet noise indicator and shock collision scheme, Part 2: The 2-D case

This is the second part to our companion paper. Herein, we generalize to two space dimensions the C-method developed in for adding localized, space-time smooth artificial viscosity to nonlinear systems of conservation laws that propagate shock waves, rarefaction waves, and contact discontinuities. For gas dynamics, the C-method couples the Euler equations to scalar reaction-diffusion equations, which we call C-equations, whose solutions serve as space-time smooth artificial viscosity indicators for shocks and contacts. We develop a high-order numerical algorithm for gas dynamics in 2-D which can accurately simulate the Rayleigh-Taylor (RT) instability with Kelvin-Helmholtz (KH) roll-up of the contact discontinuity, as well as shock collision and bounce-back. Solutions to our C-equations not only indicate the location of the shocks and contacts, but also track the geometry of the evolving fronts. This allows us to implement both directionally isotropic and anisotropic artificial viscosity schemes, the latter adding diffusion only in directions tangential to the evolving front. We additionally produce a novel shock collision indicator function, which naturally activates during shock collision, and then smoothly deactivates. Moreover, we implement a high-frequency 2-D wavelet-based noise detector together with an efficient and localized noise removal algorithm. To test the methodology, we use a highly simplified WENO-based discretization scheme. We provide numerical results for some classical 2-D test problems, including the RT problem, the Noh problem, a circular explosion problem from the Liska & Wendroff [13] review paper, the Sedov blast wave problem, the double Mach 10 reflection test, and a shock-wall collision problem. In particular, we show that our artificial viscosity method can eliminate the wall-heating phenomenon for the Noh problem, and thereby produce an accurate, non-oscillatory solution, even though our simplified WENO-type scheme fails to run for this problem.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Stochastic acceleration of energetic ions in Jupiter's magnetosphere

An equation governing the combined radial diffusion and stochastic acceleration of super-Alfvenic ions by magnetohydrodynamic waves in Jupiter's outer magnetosphere is derived. The formulation is based upon a total energy invariant of the adiabatic transport which applies to an isotropic distribution undergoing rapid pitch angle scattering by waves. An analytic solution to the double diffusion equation is obtained and numerical results are presented for two models of ion injection. The first model assumes S(+) and O(+) are injected throughout a broad region of space through photoionization of Jupiter's magnetospheric neutral wind and obtain an initial energy corresponding to the local corotation energy after pickup by the planetary magnetic field. The second model assumes a monoenergetic distribution of energetic protons is implanted in the middle magnetosphere by the action of field-aligned potential drops in Jupiter's auroral ionosphere. For both light and heavy ions the injection process creates a seed population of particles which are further accelerated nonadiabatically by the MHD waves and adiabatically through radial diffusion. A comparison of the theoretical results with a recent data analysis of Voyager low-energy charged particle measurements is made with very good agreement, thus providing a rigorous quantitative account of and definitive explanation for the high-energy ion component of Jupiter's magnetosphere.

Barbosa, D. D.↗
DOE OSTI2019

Determination of the Thicknesses and Visualization of Ion-Exchange Waveguides in Glasses by Scanning Electron Microscopy

The formation of a surface layer with silver and sodium concentration gradients in K8-grade glasses as a result of ion exchange is demonstrated by scanning electron microscopy (SEM) and electron probe X-ray microanalysis. The silver ions enriched layers with different thicknesses, depending on the duration of the ion exchange from the silver melt, are visualized in the secondary electrons mode. SEM data on the thickness of the silver enriched layer are in good agreement with the results of calculation from the diffusion equation for silver in glass. The dependence of the silver concentration on the diffusion depth is obtained using elemental composition mapping over the glass cross section. The presence of a layer with a silver concentration gradient leads to formation of a gradient waveguide. Discrete peaks in silver concentration profiles, caused by specific features of mutual diffusion of sodium and silver ions in glasses, are found.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Resistive drift wave turbulence and anomalous transport of multi-species plasma

Anomalous transport of multi-species plasma is considered with the generalized Hasegawa–Wakatani model. It is shown that the transport of all plasma species is described by fractional diffusion equations with the same effective diffusion coefficient. Strongly enhanced perturbations of heavy impurity density are found in long-living plasma flow vortices.

Physics↗
DOE OSTI2018

Lessons from HAWC pulsar wind nebulae observations: The diffusion constant is not a constant; pulsars remain the likeliest sources of the anomalous positron fraction; cosmic rays are trapped for long periods of time in pockets of inefficient diffusion

Recent TeV observations of nearby pulsars with the HAWC telescope have been interpreted as evidence that the diffusion of high-energy electrons and positrons within pulsar wind nebulae is highly inefficient compared to the rest of the interstellar medium. If the diffusion coefficient well outside the nebula is close to the value inferred for the region inside the nebula, high-energy electrons and positrons produced by the two observed pulsars could not contribute significantly to the local measured cosmic-ray flux. The HAWC Collaboration thus concluded that, under the assumption of isotropic and homogeneous diffusion, the two pulsars are ruled out as sources of the anomalous high-energy positron flux. Here, we argue that since the diffusion coefficient is likely not spatially homogeneous, the assumption leading to this conclusion is flawed. We solve the diffusion equation with a radially dependent diffusion coefficient, and show that the pulsars observed by HAWC produce potentially perfect matches to the observed high-energy positron fluxes. We also study the implications of inefficient diffusion within pulsar wind nebulae on Galactic scales, and show that cosmic rays are likely to have very long residence times in regions of inefficient diffusion. We describe how this prediction can be tested with studies of the diffuse Galactic emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Final Report - Subcontract B6183769

During my visit to LLNL during July 5{15, 2016, we worked on linear system solvers. The two level hierarchical solver that initiated our study was developed to solve linear system arising from hp adaptive finite element calculations, and is implemented in the PLTMG software package, version 12 [1]. This preconditioner typically requires 3-20% of the space used by the stiffness matrix for higher order elements. It has multigrid like convergence rates for a wide variety of PDEs (self-adjoint positive definite elliptic equations, convection dominated convection-diffusion equations, and highly inde nite Helmholtz equations, among others). The convergence rate is not independent of the polynomial degree p as p ! 1, but but remains strong for p 9, which is the highest polynomial degree allowed in PLTMG, due to limitations of numerical quadrature formulae. A more complete description of the method and some numerical experiments illustrating its effectiveness appear in [2]. Like traditional geometric multilevel methods, this scheme relies on extensive knowledge of the underlying finite element space in order to construct both the smoother and the coarse grid correction components.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Final Report - Subcontract B623760

During my visit to LLNL during July 17{27, 2017, I worked on linear system solvers. The two level hierarchical solver that initiated our study was developed to solve linear systems arising from hp adaptive finite element calculations, and is implemented in the PLTMG software package, version 12. This preconditioner typically requires 3-20% of the space used by the stiffness matrix for higher order elements. It has multigrid like convergence rates for a wide variety of PDEs (self-adjoint positive de nite elliptic equations, convection dominated convection-diffusion equations, and highly indefinite Helmholtz equations, among others). The convergence rate is not independent of the polynomial degree p as p ! 1, but but remains strong for p 9, which is the highest polynomial degree allowed in PLTMG, due to limitations of the numerical quadrature rules implemented in the software package. A more complete description of the method and some numerical experiments illustrating its effectiveness appear in. Like traditional geometric multilevel methods, this scheme relies on knowledge of the underlying finite element space in order to construct the smoother and the coarse grid correction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Final Report - Subcontract B629114

During my visits to LLNL August 20-24, 2018, and September 10-14, 2018, I worked on linear system solvers. The two level hierarchical solver that initiated our study was developed to solve linear systems arising from $hp$ adaptive finite element calculations, and is implemented in the PLTMG software package, version 13. This preconditioner typically requires 3 – 20% of the space used by the stiffness matrix for higher order elements. It has multigrid like convergence rates for a wide variety of PDEs (self-adjoint positive definite elliptic equations, convection dominated convection-diffusion equations, and highly indefinite Helmholtz equations, among others). The convergence rate is not independent of the polynomial degree $p$ as $p → ∞$, but but remains strong for $p ≤$ 9, which is the highest polynomial degree allowed in PLTMG , due to limitations of the numerical quadrature rules implemented in the software package. A more complete description of the method and some numerical experiments illustrating its effectiveness appear. Like traditional geometric multilevel methods, this scheme relies on knowledge of the underlying finite element space in order to construct the smoother and the coarse grid correction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Subcontract B635525 (Final Report)

During my visits to LLNL July 8–12 and July 22–26, 2019, I worked on linear system solvers. The two level hierarchical solver that initiated our study was developed to solve linear systems arising from $\textit{hp}$ adaptive finite element calculations, and is implemented in the $\textit{PLTMG}$ software package, version 13. This preconditioner typically requires 3–20% of the space used by the stiffness matrix for higher order elements. It has multigrid like convergence rates for a wide variety of PDEs (self-adjoint positive definite elliptic equations, convection dominated convection-diffusion equations, and highly indefinite Helmholtz equations, among others). The convergence rate is not independent of the polynomial degree $\textit{p}$ as $\textit{p} → ∞$, but remains strong for $\textit{p}$ ≤ 9, which is the highest polynomial degree allowed in $\textit{PLTMG}$, due to limitations of the numerical quadrature rules implemented in the soft- ware package. A more complete description of the method and some numerical experiments illustrating its effectiveness appear in. Like traditional geometric multilevel methods, this scheme relies on knowledge of the underlying finite element space in order to construct the smoother and the coarse grid correction.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1981

Mixing of a passive scalar in isotropic and sheared homogeneous turbulence

In order to calculate the velocity and scalar fields, the three dimensional, time-dependent equations of motion and the diffusion equation were solved numerically. The following cases were treated: isotropic, homogeneous turbulence with decay of a passive scalar; and homogeneous turbulent shear flow with a passive scalar whose mean varies linearly in the spanwise direction. The solutions were obtained at relatively low Reynolds numbers so that all of the turbulent scales could be resolved without modeling. Turbulent statistics such as integral length scales, Taylor microscales, Kolmogorov length scale, one- and two-point correlations of velocity-velocity and velocity-scalar, turbulent Prandtl/Schmidt number, r.m.s. values of velocities, the scalar quantity and pressure, skewness, decay rates, and decay exponents were calculated. The results are compared with the available expermental results, and good agreement is obtained.

Shirani, E.↗
NASA NTRS1981

Centrifugally driven diffusion of Iogenic plasma

The plasma distribution around Io as measured by Voyager 1 displays an asymmetric discontinuity at Io's orbit that has been suggested to be the signature of centrifugally driven interchange diffusion fed by plasma derived from Io. This hypothesis is explored further and found to be valid. The particular form for the diffusion coefficient appropriate to centrifugally driven turbulence is derived. The nonlinear character of this kind of diffusion is thereby made explicit. Solutions to the nonlinear, time-independent and linearized, time-dependent diffusion equations are given. These display a markedly conservative behavior. The nonlinear, steady state solutions are identical in form to the solutions of the previously studied equation of linear, atmospherically driven diffusion. The linearized, time-dependent solutions exhibit a negative feed-back quality that buffers the response of the density to changes in the source strength. Estimates of the source strength, the diffusion coefficient, and the signal propagation speed are also given.

Siscoe, G. L.↗
NASA NTRS1986

The propagation of solar flare particles in a coronal loop

A time-dependent diffusion equation with velocity-dependent diffusion and energy-loss coefficients was solved for the case where energetic solar particles are injected into a coronal loop and then diffuse out the ends of the loop into the lower corona/chromosphere. The solution yields for the case of relativistic electrons, precipitation rates and populations which are necessary for calculating thick and thin target X-ray emission. It follows that the thick target emission is necessarily delayed with respect to the particle acceleration on injection by more than the mere travel time of the particle over the loop length. In addition the time-dependent electron population at the top of the loop is calculated. This is useful in estimating the resulting micron-wave emission. The results show relative timing differences in the different emission processes which are functions of particle species, energy and the point of injection of the particles into the loop. Equivalent quantities are calculated for non-relativistic protons.

Ryan, J. M.↗