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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Thermal Transport Model for Heat Sink Design

A document discusses the development of a finite element model for describing thermal transport through microcalorimeter arrays in order to assist in heat-sinking design. A fabricated multi-absorber transition edge sensor (PoST) was designed in order to reduce device wiring density by a factor of four. The finite element model consists of breaking the microcalorimeter array into separate elements, including the transition edge sensor (TES) and the silicon substrate on which the sensor is deposited. Each element is then broken up into subelements, whose surface area subtends 10 10 microns. The heat capacity per unit temperature, thermal conductance, and thermal diffusivity of each subelement are the model inputs, as are the temperatures of each subelement. Numerical integration using the Finite in Time Centered in Space algorithm of the thermal diffusion equation is then performed in order to obtain a temporal evolution of the subelement temperature. Thermal transport across interfaces is modeled using a thermal boundary resistance obtained using the acoustic mismatch model. The document concludes with a discussion of the PoST fabrication. PoSTs are novel because they enable incident x-ray position sensitivity with good energy resolution and low wiring density.

Chervenak, James A.↗
NASA NTRS2009

Integration of the Radiation Belt Environment Model Into the Space Weather Modeling Framework

We have integrated the Fok radiation belt environment (RBE) model into the space weather modeling framework (SWMF). RBE is coupled to the global magnetohydrodynamics component (represented by the Block-Adaptive-Tree Solar-wind Roe-type Upwind Scheme, BATS-R-US, code) and the Ionosphere Electrodynamics component of the SWMF, following initial results using the Weimer empirical model for the ionospheric potential. The radiation belt (RB) model solves the convection-diffusion equation of the plasma in the energy range of 10 keV to a few MeV. In stand-alone mode RBE uses Tsyganenko's empirical models for the magnetic field, and Weimer's empirical model for the ionospheric potential. In the SWMF the BATS-R-US model provides the time dependent magnetic field by efficiently tracing the closed magnetic field-lines and passing the geometrical and field strength information to RBE at a regular cadence. The ionosphere electrodynamics component uses a two-dimensional vertical potential solver to provide new potential maps to the RBE model at regular intervals. We discuss the coupling algorithm and show some preliminary results with the coupled code. We run our newly coupled model for periods of steady solar wind conditions and compare our results to the RB model using an empirical magnetic field and potential model. We also simulate the RB for an active time period and find that there are substantial differences in the RB model results when changing either the magnetic field or the electric field, including the creation of an outer belt enhancement via rapid inward transport on the time scale of tens of minutes.

Glocer, A.↗
NASA NTRS2010

The Evolution of Remnant Ice at the Lunar South Pole from Diviner Surface Temperature Results

The Diviner lunar radiometer instrument aboard the Lunar Reconnaissance Orbiter mission has revealed large areas of lunar polar terrain with surface temperatures well below 100K. At these temperatures, the sublimation rate of water ice is well below 1 mm per billion years. In contrast, the loss rate at 120K is more than 1 meter of ice in that time consequently volatiles delivered to the coldest locations can be trapped for over 1 Ga, but will be quickly lost from warmer locales. Here we investigate the loss or retention of a layer of ice-bearing regolith at the lunar south poe, assuming contemporary surface temperature conditions and no other loss processes. We use an analytic solution for the one-dimensional diffusion equation of water ice, assuming an isothermal regolith with pore space comparable to mean grain size, 75 micrometers. Only the top meter of soil is assumed to be ice-bearing. We can then calculate the history of ice content with time based on local temperature, and predict what the epithermal neutron output would be in the presence of such a concentration of hydrogen. We compare the present, observed distribution of hydrogen with what one would expect from the temperature-dependent loss or retention of ice for various times since emplacement.

Elphic, Richard C.↗
NASA NTRS2012

A Review of the Ginzburg-Syrovatskii's Galactic Cosmic-Ray Propagation Model and its Leaky-Box Limit

Phenomenological models of galactic cosmic-ray propagation are based on a diffusion equation known as the Ginzburg-Syrovatskii s equation, or variants (or limits) of this equation. Its one-dimensional limit in a homogeneous volume, known as the leaky-box limit or model, is sketched here. The justification, utility, limitations, and a typical numerical implementation of the leaky-box model are examined in some detail.

Barghouty, A. F.↗
NASA NTRS2014

RITRACKS: A Software for Simulation of Stochastic Radiation Track Structure, Micro and Nanodosimetry, Radiation Chemistry and DNA Damage for Heavy Ions

The code RITRACKS (Relativistic Ion Tracks) has been developed over the last few years at the NASA Johnson Space Center to simulate the effects of ionizing radiations at the microscopic scale, to understand the effects of space radiation at the biological level. The fundamental part of this code is the stochastic simulation of radiation track structure of heavy ions, an important component of space radiations. The code can calculate many relevant quantities such as the radial dose, voxel dose, and may also be used to calculate the dose in spherical and cylindrical targets of various sizes. Recently, we have incorporated DNA structure and damage simulations at the molecular scale in RITRACKS. The direct effect of radiations is simulated by introducing a slight modification of the existing particle transport algorithms, using the Binary‐Encounter‐Bethe model of ionization cross sections for each molecular orbitals of DNA. The simulation of radiation chemistry is done by a step‐by‐step diffusion‐reaction program based on the Green's functions of the diffusion equation]. This approach is also used to simulate the indirect effect of ionizing radiation on DNA. The software can be installed independently on PC and tablets using the Windows operating system and does not require any coding from the user. It includes a Graphic User Interface (GUI) and a 3D OpenGL visualization interface. The calculations are executed simultaneously (in parallel) on multiple CPUs. The main features of the software will be presented.

Plante, I↗
NASA NTRS2015

Free-Energy Calculations. A Mathematical Perspective

Ion channels are pore-forming assemblies of transmembrane proteins that mediate and regulate ion transport through cell walls. They are ubiquitous to all life forms. In humans and other higher organisms they play the central role in conducting nerve impulses. They are also essential to cardiac processes, muscle contraction and epithelial transport. Ion channels from lower organisms can act as toxins or antimicrobial agents, and in a number of cases are involved in infectious diseases. Because of their important and diverse biological functions they are frequent targets of drug action. Also, simple natural or synthetic channels find numerous applications in biotechnology. For these reasons, studies of ion channels are at the forefront of biophysics, structural biology and cellular biology. In the last decade, the increased availability of X-ray structures has greatly advanced our understanding of ion channels. However, their mechanism of action remains elusive. This is because, in order to assist controlled ion transport, ion channels are dynamic by nature, but X-ray crystallography captures the channel in a single, sometimes non-native state. To explain how ion channels work, X-ray structures have to be supplemented with dynamic information. In principle, molecular dynamics (MD) simulations can aid in providing this information, as this is precisely what MD has been designed to do. However, MD simulations suffer from their own problems, such as inability to access sufficiently long time scales or limited accuracy of force fields. To assess the reliability of MD simulations it is only natural to turn to the main function of channels - conducting ions - and compare calculated ionic conductance with electrophysiological data, mainly single channel recordings, obtained under similar conditions. If this comparison is satisfactory it would greatly increase our confidence that both the structures and our computational methodologies are sufficiently accurate. Channel conductance, defined as the ratio of ionic current through the channel to applied voltage, can be calculated in MD simulations by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. If the current is small, a voltage significantly higher than the experimental one needs to be applied to collect sufficient statistics of ion crossing events. Then, the calculated conductance has to be extrapolated to the experimental voltage using procedures of unknown accuracy. Instead, we propose an alternative approach that applies if ion transport through channels can be described with sufficient accuracy by the one-dimensional diffusion equation in the potential given by the free energy profile and applied voltage. Then, it is possible to test the assumptions of the equation, recover the full voltage/current dependence, determine the reliability of the calculated conductance and reconstruct the underlying (equilibrium) free energy profile, all from MD simulations at a single voltage. We will present the underlying theory, model calculations that test this theory and simulations on ion conductance through a channel that has been extensively studied experimentally. To our knowledge this is the first case in which the complete, experimentally measured dependence of the current on applied voltage has been reconstructed from MD simulations.

free energy↗
DOE OSTI2021

Katana

Katana is built on the pre-existing MPACT and Futility libraries. Combining those solves with several new ones that are part of Katana, nodal diffusion equations are solved in 3D using pre-tabulated nodal cross sections. 3D coarse mesh finite difference acceleration is applied to accelerate convergence. Feedback effects are accounted for during the course of the solution.

Graham, Aaron M↗
DOE OSTI2025

Dark Matter-Electron Interactions Alter the Luminosity and Spectral Index of M87

We investigate the possibility that cosmic-ray electron cooling through dark matter-electron scatterings contributes to the low radiative efficiency observed in radio-loud galaxies such as M87. Light dark matter can scatter efficiently off electrons in M87, lowering the observed bolometric luminosity compared to astrophysical expectations. This consideration allows us to probe previously unexplored regions of the parameter space of dark matter-electron interactions. We further model the cosmic-ray electron distribution by numerically solving a diffusion equation along the jet and find that efficient dark matter-electron interactions can induce a flattening of the spectral index at different distances from the central supermassive black hole, in better alignment with radio observations from M87.

Hussein, Abdelaziz [MIT] (ORCID:0000000255608668)↗
NASA NTRS1992

Series integration of the diaphragm cell transport equation when the diffusion coefficient is a function of concentration

An integral of the form, t = B0 + BL ln(Delta-c) + B1(Delta-c) + B2(Delta-c)-squared + ..., where t is the time and Delta-c is the concentration difference across the frit, is derived in the case of the diaphragm cell transport equation where the interdiffusion coefficient is a function of concentration. The coefficient, B0, is a constant of the integration, while the coefficients, BL, B1, B2,..., depend in general upon the constant, the compartment volumes, and the interdiffusion coefficient and various of its concentration derivatives evaluated at the mean concentration for the cell. Explicit formulas for BL, B1, B2,... are given.

Cain, Judith B.↗
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 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↗