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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Quantum Transmission Conditions for Diffusive Transport in Graphene with Steep Potentials

We present a formal derivation of a drift-diffusion model for stationary electron transport in graphene, in presence of sharp potential profiles, such as barriers and steps. Assuming the electric potential to have steep variations within a strip of vanishing width on a macroscopic scale, such strip is viewed as a quantum interface that couples the classical regions at its left and right sides. In the two classical regions, where the potential is assumed to be smooth, electron and hole transport is described in terms of semiclassical kinetic equations. The diffusive limit of the kinetic model is derived by means of a Hilbert expansion and a boundary layer analysis, and consists of drift-diffusion equations in the classical regions, coupled by quantum diffusive transmission conditions through the interface. The boundary layer analysis leads to the discussion of a four-fold Milne (half-space, half-range) transport problem.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Fully quantum algorithm for mesoscale fluid simulations with application to partial differential equations

Fluid flow simulations marshal our most powerful computational resources. In many cases, even this is not enough. Quantum computers provide an opportunity to speed up traditional algorithms for flow simulations. We show that lattice-based mesoscale numerical methods can be executed as efficient quantum algorithms due to their statistical features. This approach revises a quantum algorithm for lattice gas automata to reduce classical computations and state preparation at every time step. For this, the algorithm approximates the qubit relative phases and subtracts them at the end of each time step. Phases are evaluated using the iterative phase estimation algorithm and subtracted using single-qubit rotation phase gates. Further, this method optimizes the quantum resource required and makes it more appropriate for near-term quantum hardware. We also demonstrate how the checkerboard deficiency that the D1Q2 scheme presents can be resolved using the D1Q3 scheme. The algorithm is validated by simulating two canonical partial differential equations: the diffusion and Burgers' equations on different quantum simulators. We find good agreement between quantum simulations and classical solutions for the presented algorithm.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1973

Matrix equations of radiation diffusion, 4

The matrices A sup (-1) and T (tau sub O = 0.1 - 1.0) are tabulated, permitting a source function and emission line profile to be easily calculated. It is assumed that the radiation diffuses with complete redistribution in frequency.

Kostyk, R. I.↗
NASA NTRS1979

Numerical methods of solving a system of multi-dimensional nonlinear equations of the diffusion type

The principles of conservation and stability of difference schemes achieved using the iteration control method were examined. For the schemes obtained of the predictor-corrector type, the conversion was proved for the control sequences of approximate solutions to the precise solutions in the Sobolev metrics. Algorithms were developed for reducing the differential problem to integral relationships, whose solution methods are known, were designed. The algorithms for the problem solution are classified depending on the non-linearity of the diffusion coefficients, and practical recommendations for their effective use are given.

Agapov, A. V.↗
NASA NTRS2000

Two-Dimensional Array Beam Scanning Via Externally and Mutually Injection Locked Coupled Oscillators

Some years ago, Stephan proposed an approach to one dimensional (linear) phased array beam steering which requires only a single phase shifter. This involves the use of a linear array of voltage-controlled electronic oscillators coupled to nearest neighbors. The oscillators are mutually injection locked by controlling their coupling and tuning appropriately. Stephan's approach consists of deriving two signals from a master oscillator, one signal phase shifted with respect to the other by means of a single phase shifter. These two signals are injected into the end oscillators of the array. The result is a linear phase progression across the oscillator array. Thus, if radiating elements are connected to each oscillator and spaced uniformly along a line, they will radiate a beam at an angle to that line determined by the phase gradient which is, in turn, determined by the phase difference between the injection signals.The beam direction is therefore controlled by adjusting this phase difference. Recently, Pogorzelski and York presented a formulation which facilitates theoretical analysis of the above beam steering technique. This was subsequently applied by Pogorzelski in analysis of two dimensional beam steering using perimeter detuning of a coupled oscillator array. The formulation is based on a continuum model in which the oscillator phases are represented by a continuous function satisfying a partial differential equation of diffusion type. This equation can be solved via the Laplace transform and the resulting solution exhibits the dynamic behavior of the array as the beam is steered. Stephan's beam steering technique can be similarly generalized to two-dimensional arrays in which the beam control signals are applied to the oscillators on the perimeter of the array. In this paper the continuum model for this two-dimensional case is developed and the dynamic solution for the corresponding aperture phase function is obtained. The corresponding behavior of the resulting far-zone radiation pattern is displayed as well.

Pogorzelski, Ronald J.↗
DOE OSTI2017

Reaction rates for reaction-diffusion kinetics on unstructured meshes

The reaction-diffusion master equation is a stochastic model often utilized in the study of biochemical reaction networks in living cells. It is applied when the spatial distribution of molecules is important to the dynamics of the system. A viable approach to resolve the complex geometry of cells accurately is to discretize space with an unstructured mesh. Diffusion is modeled as discrete jumps between nodes on the mesh, and the diffusion jump rates can be obtained through a discretization of the diffusion equation on the mesh. Reactions can occur when molecules occupy the same voxel. Here in this paper, we develop a method for computing accurate reaction rates between molecules occupying the same voxel in an unstructured mesh. For large voxels, these rates are known to be well approximated by the reaction rates derived by Collins and Kimball, but as the mesh is refined, no analytical expression for the rates exists. We reduce the problem of computing accurate reaction rates to a pure preprocessing step, depending only on the mesh and not on the model parameters, and we devise an efficient numerical scheme to estimate them to high accuracy. We show in several numerical examples that as we refine the mesh, the results obtained with the reaction-diffusion master equation approach those of a more fine-grained Smoluchowski particle-tracking model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Discrete generative diffusion models without stochastic differential equations: A tensor network approach

Diffusion models (DMs) are a class of generative machine learning methods that sample a target distribution by transforming samples of a trivial (often Gaussian) distribution using a learned stochastic differential equation. In standard DMs, this is done by learning a “score function” that reverses the effect of adding diffusive noise to the distribution of interest. Here we consider the generalisation of DMs to lattice systems with discrete degrees of freedom, and where noise is added via Markov chain jump dynamics. We show how to use tensor networks (TNs) to efficiently define and sample such “discrete diffusion models” (DDMs) without explicitly having to solve a stochastic differential equation. We show the following: (i) by parametrising the data and evolution operators as TNs, the denoising dynamics can be represented exactly; (ii) the auto-regressive nature of TNs allows to generate samples efficiently and without bias; (iii) for sampling Boltzmann-like distributions, TNs allow to construct an efficient learning scheme that integrates well with Monte Carlo. We illustrate this approach to study the equilibrium of two models with non-trivial thermodynamics, the d = 1 constrained Fredkin chain and the d = 2 Ising model. Published by the American Physical Society 2025

Causer, Luke (ORCID:0000000194243473)↗
DOE OSTI2024

A Second Moment Method for k -Eigenvalue Acceleration with Continuous Diffusion and Discontinuous Transport Discretizations

The second moment method is a linear acceleration technique that couples the transport equation to a diffusion equation with transport-dependent additive closures. The resulting low-order diffusion equation can be discretized independent of the transport discretization, unlike diffusion synthetic acceleration, and is symmetric positive definite, unlike quasidiffusion. While this method has been shown to be comparable to quasidiffusion in iterative performance for fixed source and time-dependent problems, it is largely unexplored as an eigenvalue problem acceleration scheme due to the belief that the resulting inhomogeneous source makes the problem ill posed. Recently, a preliminary feasibility study was performed on the second moment method for eigenvalue problems. The results suggested comparable performance to quasidiffusion and more robust performance than diffusion synthetic acceleration. This work extends the initial study to more realistic reactor problems using state-of-the-art discretization techniques. Finally, the results in this paper show that the second moment method is more computationally efficient than its alternatives on complex reactor problems with unstructured meshes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

A flexible nonlinear diffusion acceleration method for the S N transport equations discretized with discontinuous finite elements

This paper presents a flexible nonlinear diffusion acceleration (NDA) method that discretizes both the S N transport equation and the diffusion equation using the discontinuous finite element method (DFEM). The method is flexible in that the diffusion equation can be discretized on a coarser mesh with the only restriction that it is nested within the transport mesh and the FEM shape function orders of the two equations can be different. The consistency of the transport and diffusion solutions at convergence is defined by using a projection operator mapping the transport into the diffusion FEM space. The diffusion weak form is based on the modified incomplete interior penalty (MIP) diffusion DFEM discretization that is extended by volumetric drift, interior face, and boundary closure terms. In contrast to commonly used coarse mesh finite difference (CMFD) methods, the presented NDA method uses a full FEM discretized diffusion equation for acceleration. Suitable projection and prolongation operators arise naturally from the FEM framework. Via Fourier analysis and numerical experiments for a one-group, fixed source problem the following properties of the NDA method are established for structured quadrilateral meshes: (1) the presented method is unconditionally stable and effective in the presence of mild material heterogeneities if the same mesh and identical shape functions either of the bilinear or biquadratic type are used, (2) the NDA method remains unconditionally stable in the presence of strong heterogeneities, (3) the NDA method with bilinear elements extends the range of effectiveness and stability by a factor of two when compared to CMFD if a coarser diffusion mesh is selected. In addition, the method is tested for solving the C5G7 multigroup, eigenvalue problem using coarse and fine mesh acceleration. Finally, while NDA does not offer an advantage over CMFD for fine mesh acceleration, it reduces the iteration count required for convergence by almost a factor of two in the case of coarse mesh acceleration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

DG-IMEX Stochastic Galerkin Schemes for Linear Transport Equation with Random Inputs and Diffusive Scalings

In this paper, we consider the linear transport equation under diffusive scaling and with random inputs. The method is based on the generalized polynomial chaos approach in the stochastic Galerkin framework. Several theoretical aspects will be addressed. Additionally, a uniform numerical stability with respect to the Knudsen number ϵ, and a uniform in ϵ error estimate is given. For temporal and spatial discretizations, we apply the implicit–explicit scheme under the micro–macro decomposition framework and the discontinuous Galerkin method, as proposed in Jang et al. (SIAM J Numer Anal 52:2048–2072, 2014) for deterministic problem. Lastly, we provide a rigorous proof of the stochastic asymptotic-preserving (sAP) property. Extensive numerical experiments that validate the accuracy and sAP of the method are conducted.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

A two-dimensional kinematic dynamo model of the ionospheric magnetic field at Venus

The results of a high-resolution, two-dimensional, time dependent, kinematic dynamo model of the ionospheric magnetic field of Venus are presented. Various one-dimensional models are considered and the two-dimensional model is then detailed. In this model, the two-dimensional magnetic induction equation, the magnetic diffusion-convection equation, is numerically solved using specified plasma velocities. Origins of the vertical velocity profile and of the horizontal velocities are discussed. It is argued that the basic features of the vertical magnetic field profile remain unaltered by horizontal flow effects and also that horizontal plasma flow can strongly affect the magnetic field for altitudes above 300 km.

Cravens, T. E.↗
DOE OSTI2018

Equilibrium, Kinetic and Thermodynamic Studies Of Lead (II) Sorption on Hydrolyzed Starch Graft Copolymers

Hydrolyzed starch graft poly (acrylic acid) and starch graft poly(acrylonitrile) copolymers were used as a sorbents for the removal of Pb(II) from aqueous solution in batch process. The sorbents were characterized by scanning electron microscopy, Fourier transform infrared spectroscopy and thermogravimetric analysis. The effect of contact time, initial concentration and temperature were assessed to optimize the conditions for maximum sorption of the metal ions. The experimental data were analyzed by Langmuir, Freundlich, Temkin, Dubinin Radushkevich and Sips isotherm models. The Freundlich and the Sips isotherms confirmed the applicability of the models. The sorption capacity of the sorbent was found to be 118.61 mg/g for starch graft poly (acrylic) acid and 115.83 mg/g for starch graft poly (acrylonitrile). The Sorption kinetics was assessed by Lagergren pseudo first order, pseudo second order, Elovich equation, Intraparticle diffusion and the liquid film diffusion models. The experimental data fitted very well with the Elovich equation and the Intraparticle kinetic models. The thermodynamic analysis showed that the sorption was a spontaneous and endothermic process. The results indicated that starch graft copolymers can be used as an effective sorbents for Pb (II) removal from aqueous media.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Existence in the Sense of Sequences of Stationary Solutions for Some Non-Fredholm Integro-Differential Equations

We establish the existence in the sense of sequences of stationary solutions for some reaction-diffusion type equations in appropriate H{sup 2} spaces. It is shown that, under reasonable technical conditions, the convergence in L{sup 1} of the integral kernels implies the existence and convergence in H{sup 2} of solutions. The nonlocal elliptic equations involve second order differential operators with and without the Fredholm property. Bibliography: 21 titles.

97 MATHEMATICS AND COMPUTING↗