Generalized Landau kinetic equation describing the transverse diffusion of a nongyrotropic plasma column
Plasma kinetic theory application to particle distribution function for diffusion effects
PLASMA PHYSICS↗
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.
A High-Order Nonlinear Diffusion Acceleration for the SN Equations Discretized with the Discontinuous FEM II: Fourier Analysis
In this paper, a Fourier Analysis (FA) of a new high order Nonlinear Diffusion Acceleration Method (NDA) for the $S_N$ equations discretized with the Discontinuous Finite Element Method (DFEM) is presented. In contrast to coarse mesh finite difference methods (CMFD), the associated diffusion system is discretized using DFEM with a function space identical to the $S_N$ equations. The utilized $S_N$ solver is the latest addition to the Rattlesnake\cite{Wang2013} radiation transport code that is based on the MOOSE\cite{Gaston2009} framework.
97 MATHEMATICS AND COMPUTING↗
A Research Program for Fission Product/Dust Transport in HTGR’s
High and Very High Temperatures Gas Reactors (HTGRs/VHTRs) have five barriers to fission product (FP) release: the TRISO fuel coating, the fuel elements, the core graphite, the primary coolant system, and the reactor building. This project focused on measurements and computations of FP diffusion in graphite, FP adsorption on graphite and FP interactions with dust particles of arbitrary shape. Diffusion Coefficients of Cs and Iodine in two nuclear graphite were obtained by the release method and use of Inductively Coupled Plasma-Mass Spectroscopy (ICP-MS) and Instrumented Neutron Activation Analysis (INAA). A new mathematical model for fission gas release from nuclear fuel was also developed. Several techniques were explored to measure adsorption isotherms, notably a Knudsen Effusion Mass Spectrometer (KEMS) and Instrumented Neutron Activation Analysis (INAA). Some of these measurements are still in progress. The results will be reported in a supplemental report later. Studies of FP interactions with dust and shape factors for both chain-like particles and agglomerates over a wide size range were obtained through solutions of the diffusion and transport equations. The Green's Function Method for diffusion and Monte Carlo technique for transport were used, and it was found that the shape factors are sensitive to the particle arrangements, and that diffusion and transport of FPs can be hindered. Several journal articles relating to the above work have been published, and more are in submission and preparation.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Influence of parent and daughter decay on tracer diffusion profiles
Differential equation for parent and daughter decay influence on radioactive tracer diffusion profiles
PARTICLE DIFFUSION↗
Trap-Assisted Dopant Compensation Prevents Shunting in poly-Si Passivating Interdigitated Back Contact Silicon Solar Cells
Interdigitated back contact (IBC) solar cells achieve the highest efficiencies of single-junction architectures, but complicated patterning of the rear fingers and spreading of dopants during processing inhibit their mainstream adoption due to concerns of shunting between the IBC fingers. One method of simplifying patterning at the rear is by using contact masks combined with plasma-enhanced chemical vapor deposition (PECVD) or ion implantation. However, the intrinsic isolation region becomes contaminated during high-temperature annealing by lateral diffusion of dopants and during masked PECVD by spreading of dopant radicals through region between the mask and the substrate. Despite this contamination, we show through scanning spreading resistance microscopy and Kelvin probe force microscopy that a ~20 µm wide compensating region exists with high enough resistivity to prevent shunting. We model this p-i-n poly-Si system using two simulation models: a simple resistor model considering only the capture of charge carriers by trap defects in poly-Si to reduce the conductivity, and a more refined 1-dimensional finite element model using Poisson’s equation, drift-diffusion equations, and recombination of carriers. Using this model, we show that high defect density significantly decreases the current across the region between the p- and n-type fingers, preventing shunting.
atom probe↗
Implementation of the full viscoresistive magnetohydrodynamic equations in a nonlinear finite element code
Numerical simulations form an indispensable tool to understand the behavior of a hot plasma that is created inside a tokamak for providing nuclear fusion energy. Various aspects of tokamak plasmas have been successfully studied through the reduced magnetohydrodynamic (MHD) model. The need for more complete modeling through the full MHD equations is addressed here. Our computational method is presented along with measures against possible problems regarding pollution, stability, and regularity. The problem of ensuring continuity of solutions in the center of a polar grid is addressed in the context of a finite element discretization of the full MHD equations. A rigorous and generally applicable solution is proposed here. Useful analytical test cases are devised to verify the correct implementation of the momentum and induction equation, the hyperdiffusive terms, and the accuracy with which highly anisotropic diffusion can be simulated. A striking observation is that highly anisotropic diffusion can be treated with the same order of accuracy as isotropic diffusion, even on non-aligned grids, as long as these grids are generated with sufficient care. This property is shown to be associated with our use of a magnetic vector potential to describe the magnetic field. Several well-known instabilities are simulated to demonstrate the capabilities of the new method. The linear growth rate of an internal kink mode and a tearing mode are benchmarked against the results of a linear MHD code. The evolution of a tearing mode and the resulting magnetic islands are simulated well into the nonlinear regime. The results are compared with predictions from the reduced MHD model. Finally, a simulation of a ballooning mode illustrates the possibility to use our method as an ideal MHD method without the need to add any physical dissipation.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Optimization of four terminal rear heterojunction GaAs on Si interdigitated back contact tandem solar cells
High-efficiency, four-terminal tandem solar cells composed of thin GaAs films mechanically stacked onto interdigitated back contact silicon solar cells with a glass interlayer are demonstrated. The optimal thickness of the absorber layer of a rear heterojunction GaAs subcell for use in four terminal tandem solar cells was studied. GaAs top cells with absorber layer thicknesses of 1.5, 1.9, 2.3, 2.8, and 3.5 µm were fabricated on glass and mechanically stacked onto interdigitated back-contact Si bottom cells. All tandem cells were found to have efficiencies above 30% under the AM1.5 G spectrum demonstrating a relatively weak sensitivity to thickness in the four-terminal configuration. We found the 2.8 µm absorber layer cell to have the highest top cell and tandem cell efficiency at 26.38% and 32.57%, respectively. Optical modeling with transfer matrix method for the planar top cell and Lambertian light trapping in the textured Si subcell, along with drift-diffusion Hovel equations, were used to show photon recycling enhancement to the effective diffusion length and VOC of the top cell as a result of the low-index glass interlayer.
14 SOLAR ENERGY↗
A plateau in the sensitivity of a compact optically pumped atomic magnetometer
In a compact optically pumped atomic magnetometer (OPAM), there is a plateau in the sensitivity where the dependence of the sensitivity on pumping power is small compared with that predicted by the uniform polarization model. The mechanism that generates this plateau was explained by numerical analysis. The distribution of spin polarization in the alkali metal cell of an OPAM was modeled using the Bloch equation incorporating a diffusion term and an equation for the attenuation of the pump beam. The model was well-fitted to the experimental results for a module with a cubic cell with 20 mm sides and pump and probe beams with 8 mm diameter. On the plateau, strong magnetic response was generated at the regions that were not illuminated directly by the intense pump beam, while at the same time spin polarization as large as 0.5 was maintained due to diffusion of the spin-polarized atoms. Thus, the sensitivity of the magnetometer monitored with a probe beam decreases only slightly with increasing pump beam intensity because the spin polarization under an intense pump beam is saturated. This plateau, which is characteristic of this type of magnetometer using a narrow pump and probe beams, can be used in arrays of magnetometers because it enables stable operation with little sensitivity fluctuation from changes in pump beam power.
36 MATERIALS SCIENCE↗
Momentum-dependent diffusive particle acceleration in modified shock fronts
In the presently derived analytic solutions of the steady transport equation for diffusive particle acceleration in a modified, planar shock front having free escape boundaries, the fluid velocity profile through the shock transition decreases monotonically between the upstream and downstream boundaries. The spatial diffusion coefficient's spatial dependence is linked to that of the fluid velocity profile. Attention is given to the solution corresponding to monoenergetic particle injection at the shock front, with free particle escape at finite distances both upstream and downstream of the shock. The accelerated particle spectrum is dominated at high energies by an exponential cutoff, due to the competition between acceleration by the first-order Fermi mechanism and particle loss through the free escape boundaries.
Webb, G. M.↗
Mathematical Model of Transfer and Deposition of Finely Dispersed Particles in a Turbulent Flow of Emulsions and Suspensions
The problem of modeling turbulent transfer of finely dispersed particles in liquids has been considered. An approach is used where the transport of particles is represented in the form of a variety of the diffusion process with the coefficient of turbulent transfer to the wall. Differential equations of transfer are written for different cases, and a solution of the cell model is obtained for calculating the efficiency of separation in a channel. Based on the theory of turbulent transfer of particles and of the boundary layer model, an expression has been obtained for calculating the rate of turbulent deposition of finely dispersed particles. The application of this expression in determining the efficiency of physical coagulation of emulsions in different channels and on the surface of chaotic packings is shown.
42 ENGINEERING↗
Anomalous Transport of Cosmic Rays in a Nonlinear Diffusion Model
We investigate analytically and numerically the transport of cosmic rays following their escape from a shock or another localized acceleration site. Observed cosmic-ray distributions in the vicinity of heliospheric and astrophysical shocks imply that anomalous, superdiffusive transport plays a role in the evolution of the energetic particles. Several authors have quantitatively described the anomalous diffusion scalings, implied by the data, by solutions of a formal transport equation with fractional derivatives. Yet the physical basis of the fractional diffusion model remains uncertain. We explore an alternative model of the cosmic-ray transport: a nonlinear diffusion equation that follows from a self-consistent treatment of the resonantly interacting cosmic-ray particles and their self-generated turbulence. The nonlinear model naturally leads to superdiffusive scalings. In the presence of convection, the model yields a power-law dependence of the particle density on the distance upstream of the shock. Although the results do not refute the use of a fractional advection–diffusion equation, they indicate a viable alternative to explain the anomalous diffusion scalings of cosmic-ray particles.
79 ASTRONOMY AND ASTROPHYSICS↗
Discovering a reaction–diffusion model for Alzheimer’s disease by combining PINNs with symbolic regression
Misfolded tau proteins play a critical role in the progression and pathology of Alzheimer's disease. Recent studies suggest that the spatio-temporal pattern of misfolded tau follows a reaction-diffusion type equation. However, the precise mathematical model and parameters that characterize the progression of misfolded protein across the brain remain incompletely understood. Here, we use deep learning and artificial intelligence to discover a mathematical model for the progression of Alzheimer's disease using longitudinal tau positron emission tomography from the Alzheimer's Disease Neuroimaging Initiative database. Specifically, we integrate physics informed neural networks (PINNs) and symbolic regression to discover a reaction-diffusion type partial differential equation for tau protein misfolding and spreading. First, we demonstrate the potential of our model and parameter discovery on synthetic data. Then, we apply our method to discover the best model and parameters to explain tau imaging data from 46 individuals who are likely to develop Alzheimer's disease and 30 healthy controls. Our symbolic regression discovers different misfolding models f(c) for two groups, with a faster misfolding for the Alzheimer's group, f(c) = 0.23c 3 – 1.34c 2 + 1.11c, than for the healthy control group, f(c) = –c 3 + 0.62c 2 + 0.39c. Our results suggest that PINNs, supplemented by symbolic regression, can discover a reaction-diffusion type model to explain misfolded tau protein concentrations in Alzheimer's disease. Furthermore, we expect our study to be the starting point for a more holistic analysis to provide image-based technologies for early diagnosis, and ideally early treatment of neurodegeneration in Alzheimer's disease and possibly other misfolding-protein based neurodegenerative disorders.
60 APPLIED LIFE SCIENCES↗
The Poisson Equation and Estimates for Distances Between Stationary Distributions of Diffusions
We estimate distances between stationary solutions to Fokker–Planck–Kolmogorov equations with different diffusion and drift coefficients. To this end we study the Poisson equation on the whole space. We have obtained sufficient conditions for stationary solutions to satisfy the Poincaré and logarithmic Sobolev inequalities.
97 MATHEMATICS AND COMPUTING↗
A computationally-efficient method for flamelet calculations
A new open-source code for the simulation of the diffusion flamelet equations is proposed. Emphasis is placed on using an approximate Jacobian to reduce the computational cost of the matrix operations. Performance of the proposed solvers is tested by performing flamelet calculations with kinetic mechanisms of varying sizes. For the unity Lewis number equations, the present iterative Newton solver using an approximate Jacobian greatly outperforms direct Newton solvers using exact Jacobians. The computation cost scales linearly with the number of species, leading to a reduction in solution times by two orders of magnitude for mechanisms containing thousands of species. The applicability of the Jacobian approximations to the solution of the non-unity Lewis number flamelet equations is assessed. The approximations are generally inadequate to solve the full non-unity Lewis number equations but can be used in some applications depending on the balance of terms in the flamelet equations. As an example, the flamelet solver is applied to the study of sooting tendencies in laminar co-flow diffusion flames where modified non-unity Lewis number flamelet equations, previously shown to accurately reproduce experimentally-measured Yield Sooting Indices (YSI), are solved. Here, the accelerated flamelet solver is well suited for sensitivity analysis and uncertainty quantification with large detailed kinetic mechanisms, tasks for which the computational cost was previously prohibitive.
42 ENGINEERING↗
Reaction rates for mesoscopic reaction-diffusion kinetics
The mesoscopic reaction-diffusion master equation (RDME) is a popular modeling framework frequently applied to stochastic reaction-diffusion kinetics in systems biology. The RDME is derived from assumptions about the underlying physical properties of the system, and it may produce unphysical results for models where those assumptions fail. In that case, other more comprehensive models are better suited, such as hard-sphere Brownian dynamics (BD). Although the RDME is a model in its own right, and not inferred from any specific microscale model, it proves useful to attempt to approximate a microscale model by a specific choice of mesoscopic reaction rates. In this paper we derive mesoscopic scale-dependent reaction rates by matching certain statistics of the RDME solution to statistics of the solution of a widely used microscopic BD model: the Smoluchowski model with a Robin boundary condition at the reaction radius of two molecules. We also establish fundamental limits on the range of mesh resolutions for which this approach yields accurate results and show both theoretically and in numerical examples that as we approach the lower fundamental limit, the mesoscopic dynamics approach the microscopic dynamics. Finally, we show that for mesh sizes below the fundamental lower limit, results are less accurate. Thus, the lower limit determines the mesh size for which we obtain the most accurate results.
74 ATOMIC AND MOLECULAR PHYSICS↗
Neutron diffusion calculation in heterogeneous geometry based on local/global iteration using proper orthogonal decomposition
This study newly proposes a heterogeneous core calculation method based on local/global iteration using proper orthogonal decomposition (POD). By using the singular value decomposition (SVD) and the low-rank approximation, appropriate POD bases for expanding the neutron flux can be obtained from snapshot data of the neutron flux obtained by fine mesh calculations. By projection using the POD bases, the dimension of the target equation (e.g., discretized neutron diffusion equation) can be dramatically reduced. In the proposed method, POD is effectively applied to each single assembly calculation (local calculation). Furthermore, using the local/global iteration, the effective neutron multiplication factor and the neutron flux distribution in the whole core geometry can be obtained by combining the numerical results of the local calculation for each fuel assembly and the global calculation for the whole core. As a feasibility study, the proposed method is applied to a one-dimensional heterogeneous core analysis, and the accuracy is investigated by changing the total number of POD bases. (authors)
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Large-Time Behavior of Solutions to Vlasov-Poisson-Fokker-Planck Equations: From Evanescent Collisions to Diffusive Limit
The present contribution investigates the dynamics generated by the two-dimensional Vlasov-Poisson-Fokker-Planck equation for charged particles in a steady inhomogeneous background of opposite charges. We provide global in time estimates that are uniform with respect to initial data taken in a bounded set of a weighted L{sup 2} space, and where dependencies on the mean-free path τ and the Debye length δ are made explicit. In our analysis the mean free path covers the full range of possible values: from the regime of evanescent collisions τ→∞ to the strongly collisional regime τ→0. As a counterpart, the largeness of the Debye length, that enforces a weakly nonlinear regime, is used to close our nonlinear estimates. Accordingly we pay a special attention to relax as much as possible the τ-dependent constraint on δ ensuring exponential decay with explicit τ-dependent rates towards the stationary solution. In the strongly collisional limit τ→0, we also examine all possible asymptotic regimes selected by a choice of observation time scale. Here also, our emphasis is on strong convergence, uniformity with respect to time and to initial data in bounded sets of a L{sup 2} space. Our proofs rely on a detailed study of the nonlinear elliptic equation defining stationary solutions and a careful tracking and optimization of parameter dependencies of hypocoercive/hypoelliptic estimates.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
PPINN: Parareal physics-informed neural network for time-dependent PDEs
Physics-informed neural networks (PINNs) encode physical conservation laws and prior physical knowledge into the neural networks, ensuring the correct physics is represented accurately while alleviating the need for supervised learning to a great degree. While effective for relatively short-term time integration, when long time integration of the time-dependent PDEs is sought, the time–space domain may become arbitrarily large and hence training of the neural network may become prohibitively expensive. To this end, we develop a parareal physics-informed neural network (PPINN), hence decomposing a long-time problem into many independent short-time problems supervised by an inexpensive/fast coarse-grained (CG) solver. In particular, the serial CG solver is designed to provide approximate predictions of the solution at discrete times, while initiate many fine PINNs simultaneously to correct the solution iteratively. There is a two-fold benefit from training PINNs with small-data sets rather than working on a large-data set directly, i.e., training of individual PINNs with small-data is much faster, while training the fine PINNs can be readily parallelized. Consequently, compared to the original PINN approach, the proposed PPINN approach may achieve a significant speed-up for long-time integration of PDEs, assuming that the CG solver is fast and can provide reasonable predictions of the solution, hence aiding the PPINN solution to converge in just a few iterations. To investigate the PPINN performance on solving time-dependent PDEs, we first apply the PPINN to solve the Burgers equation, and subsequently we apply the PPINN to solve a two-dimensional nonlinear diffusion–reaction equation. Furthermore, our results demonstrate that PPINNs converge in a few iterations with significant speed-ups proportional to the number of time-subdomains employed.
42 ENGINEERING↗