Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical diffusion”

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,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

MIGRATION OF SMALL MOONS IN SATURN's RINGS

The motions of small moons through Saturn's rings provide excellent tests of radial migration models. In theory, torque exchange between these moons and ring particles leads to radial drift. We predict that moons with Hill radii r {sub H} {approx} 2-24 km should migrate through the A ring in 1000 yr. In this size range, moons orbiting in an empty gap or in a full ring eventually migrate at the same rate. Smaller moons or moonlets-such as the propellers-are trapped by diffusion of disk material into corotating orbits, creating inertial drag. Larger moons-such as Pan or Atlas-do not migrate because of their own inertia. Fast migration of 2-24 km moons should eliminate intermediate-size bodies from the A ring and may be responsible for the observed large-radius cutoff of r {sub H} {approx} 1-2 km in the size distribution of the A ring's propeller moonlets. Although the presence of Daphnis (r {sub H} Almost-Equal-To 5 km) inside the Keeler gap challenges this scenario, numerical simulations demonstrate that orbital resonances and stirring by distant, larger moons (e.g., Mimas) may be important factors. For Daphnis, stirring by distant moons seems the most promising mechanism to halt fast migration. Alternatively, Daphnis may be a recent addition to the ring that is settling into a low inclination orbit in {approx}10{sup 3} yr prior to a phase of rapid migration. We provide predictions of observational constraints required to discriminate among possible scenarios for Daphnis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Dimerization Induced Deprotonation of Water on RuO2(110)

RuO2 has proven to be indispensable as a co-catalyst in numerous systems designed for photocatalytic water splitting. In this study we have carried out a detailed mechanistic study of water behavior on the most stable RuO2 face, RuO2(110), by employing variable temperature scanning tunneling microscopy and density functional theory calculations. We show that water monomers adsorb molecularly on Ru sites, become mobile above 238 K, diffuse along the Ru rows and form water dimers. The onset for dimer diffusion is observed at ~277 K indicating significantly higher diffusion barrier than that for monomers. More importantly, we find that water dimers deprotonate readily to form Ru-bound H3O2 and bridging OH species. The observed behavior is compared and contrasted with that observed for water on isostructural rutile TiO2(110).

Environmental Molecular Sciences Laboratory↗
DOE OSTI2015

On the pressure effect in energetic deposition of Cu thin films by modulated pulsed power magnetron sputtering: A global plasma model and experiments

The modulated pulsed power magnetron sputtering (MPPMS) discharge processes are numerically modeled and experimentally investigated, in order to explore the effect of the pressure on MPPMS discharges as well as on the microstructure of the deposited thin films. A global plasma model has been developed based on a volume-averaged global description of the ionization region, considering the loss of electrons by cross-B diffusion. The temporal variations of internal plasma parameters at different pressures from 0.1 to 0.7 Pa are obtained by fitting the model to duplicate the experimental discharge data, and Cu thin films are deposited by MPPMS at the corresponding pressures. The surface morphology, grain size and orientation, and microstructure of the deposited thin films are investigated by scanning electron microscopy, transmission electron microscopy, and x-ray diffraction. By increasing the pressure from 0.1 to 0.7 Pa, both the ion bombardment energy and substrate temperature which are estimated by the modeled plasma parameters decrease, corresponding to the observed transition of the deposited thin films from a void free structure with a wide distribution of grain size (zone T) into an underdense structure with a fine fiber texture (zone 1) in the extended structure zone diagram (SZD). The microstructure and texture transition of Cu thin films are well-explained by the extended SZD, suggesting that the primary plasma processes are properly incorporated in the model. The results contribute to the understanding of the characteristics of MPPMS discharges, as well as its correlation with the microstructure and texture of deposited Cu thin films.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Acceleration of relativistic electrons by magnetohydrodynamic turbulence: Implications for non-thermal emission from black hole accretion disks

We use analytic estimates and numerical simulations of test particles interacting with magnetohydrodynamic (MHD) turbulence to show that subsonic MHD turbulence produces efficient second-order Fermi acceleration of relativistic particles. This acceleration is not well described by standard quasi-linear theory but is a consequence of resonance broadening of wave-particle interactions in MHD turbulence. We provide momentum diffusion coefficients that can be used for astrophysical and heliospheric applications and discuss the implications of our results for accretion flows onto black holes. In particular, we show that particle acceleration by subsonic turbulence in radiatively inefficient accretion flows can produce a non-thermal tail in the electron distribution function that is likely important for modeling and interpreting the emission from low-luminosity systems such as Sgr A* and M87.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Forecasting fluid-injection induced seismicity to choose the best injection strategy for safety and efficiency

Induced seismicity poses a challenge to the development of Enhanced Geothermal Systems (EGS). Improving monitoring and forecasting techniques is essential to mitigate induced seismicity and thereby fostering a positive perception of EGS projects among local authorities and population. Induced seismicity is the result of complex and coupled thermo-hydro-mechanical-chemical mechanisms. Injection flux and pressure are crucial controlling parameters for both hydraulic stimulation and circulation protocols. We develop a methodology combining a hydro-mechanical model with a seismicity rate model to estimate the magnitude and frequency of mainshocks and aftershocks induced by fluid injection. We apply the methodology to the case of the Basel EGS (2006, Switzerland) to compare the effects of progressive, cyclic and constant injections on the mechanical response of discrete faults. Results from the coupled hydro-mechanical models show that the pore pressure diffusion and consequent enhancement of fault permeability are limited to the vicinity of the injection well during cyclic injection. Additionally, constant injection induces seismicity from the start of the injection but enhances the permeability of most of the faults within a shorter duration, inducing less post-injection seismicity. The methodology can be adapted to any numerical model and allows new projects to be developed by anticipating the safest injection protocol. This article is part of the theme issue ‘Induced seismicity in coupled subsurface systems’.

58 GEOSCIENCES↗
DOE OSTI2014

Modeling Damage Complexity-Dependent Non-Homologous End-Joining Repair Pathway

Non-homologous end joining (NHEJ) is the dominant DNA double strand break (DSB) repair pathway and involves several repair proteins such as Ku, DNA-PKcs, and XRCC4. It has been experimentally shown that the choice of NHEJ proteins is determined by the complexity of DSB. In this paper, we built a mathematical model, based on published data, to study how NHEJ depends on the damage complexity. Under an appropriate set of parameters obtained by minimization technique, we can simulate the kinetics of foci track formation in fluorescently tagged mammalian cells, Ku80-EGFP and DNA-PKcs-YFP for simple and complex DSB repair, respectively, in good agreement with the published experimental data, supporting the notion that simple DSB undergo fast repair in a Ku-dependent, DNA-PKcs-independent manner, while complex DSB repair requires additional DNA-PKcs for end processing, resulting in its slow repair, additionally resulting in slower release rate of Ku and the joining rate of complex DNA ends. Based on the numerous experimental descriptions, we investigated several models to describe the kinetics for complex DSB repair. An important prediction of our model is that the rejoining of complex DSBs is through a process of synapsis formation, similar to a second order reaction between ends, rather than first order break filling/joining. The synapsis formation (SF) model allows for diffusion of ends before the synapsis formation, which is precluded in the first order model by the rapid coupling of ends. Therefore, the SF model also predicts the higher number of chromosomal aberrations observed with high linear energy transfer (LET) radiation due to the higher proportion of complex DSBs compared to low LET radiation, and an increased probability of misrejoin following diffusion before the synapsis is formed, while the first order model does not provide a mechanism for the increased effectiveness in chromosomal aberrations observed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

MIDA: A Multimodal Imaging-Based Detailed Anatomical Model of the Human Head and Neck

Computational modeling and simulations are increasingly being used to complement experimental testing for analysis of safety and efficacy of medical devices. Multiple voxel- and surface-based whole- and partial-body models have been proposed in the literature, typically with spatial resolution in the range of 1–2 mm and with 10–50 different tissue types resolved. We have developed a multimodal imaging-based detailed anatomical model of the human head and neck, named “MIDA”. The model was obtained by integrating three different magnetic resonance imaging (MRI) modalities, the parameters of which were tailored to enhance the signals of specific tissues: i) structural T1- and T2-weighted MRIs; a specific heavily T2-weighted MRI slab with high nerve contrast optimized to enhance the structures of the ear and eye; ii) magnetic resonance angiography (MRA) data to image the vasculature, and iii) diffusion tensor imaging (DTI) to obtain information on anisotropy and fiber orientation. The unique multimodal high-resolution approach allowed resolving 153 structures, including several distinct muscles, bones and skull layers, arteries and veins, nerves, as well as salivary glands. The model offers also a detailed characterization of eyes, ears, and deep brain structures. A special automatic atlas-based segmentation procedure was adopted to include a detailed map of the nuclei of the thalamus and midbrain into the head model. The suitability of the model to simulations involving different numerical methods, discretization approaches, as well as DTI-based tensorial electrical conductivity, was examined in a case-study, in which the electric field was generated by transcranial alternating current stimulation. The voxel- and the surface-based versions of the models are freely available to the scientific community.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

Effects of anisotropies in turbulent magnetic diffusion in mean-field solar dynamo models

We study how anisotropies of turbulent diffusion affect the evolution of large-scale magnetic fields and the dynamo process on the Sun. The effect of anisotropy is calculated in a mean-field magnetohydrodynamics framework assuming that triple correlations provide relaxation to the turbulent electromotive force (so-called the 'minimal τ-approximation'). We examine two types of mean-field dynamo models: the well-known benchmark flux-transport model and a distributed-dynamo model with a subsurface rotational shear layer. For both models, we investigate effects of the double- and triple-cell meridional circulation, recently suggested by helioseismology and numerical simulations. To characterize the anisotropy effects, we introduce a parameter of anisotropy as a ratio of the radial and horizontal intensities of turbulent mixing. It is found that the anisotropy affects the distribution of magnetic fields inside the convection zone. The concentration of the magnetic flux near the bottom and top boundaries of the convection zone is greater when the anisotropy is stronger. It is shown that the critical dynamo number and the dynamo period approach to constant values for large values of the anisotropy parameter. The anisotropy reduces the overlap of toroidal magnetic fields generated in subsequent dynamo cycles, in the time-latitude 'butterfly' diagram. If we assume that sunspots are formed in the vicinity of the subsurface shear layer, then the distributed dynamo model with the anisotropic diffusivity satisfies the observational constraints from helioseismology and is consistent with the value of effective turbulent diffusion estimated from the dynamics of surface magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Lithiation Mechanism of Tunnel-Structured MnO 2 Electrode Investigated by In Situ Transmission Electron Microscopy

Manganese oxide (α-MnO 2 ) has been considered as a promising energy material, including as a lithium-based battery electrode candidate, due to its environmental friendliness. Thanks to its unique 1D [2 × 2] tunnel structure, α-MnO 2 can be applied to a cathode by insertion reaction and to an anode by conversion reaction in corresponding voltage ranges, in a lithium-based battery. Numerous reports have attributed its remarkable performance to its unique tunnel structure; however, the precise electrochemical reaction mechanism remains unknown. In this study, finding of the lithiation mechanism of α-MnO 2 nanowire by in situ transmission electron microscopy (TEM) is reported. By elaborately modifying the existing in situ TEM experimental technique, rapid lithium-ion diffusion through the tunnels is verified. Furthermore, by tracing the full lithiation procedure, the evolution of the MnO intermediate phase and the development of the MnO and Li 2 O phases with preferred orientations is demonstrated, which explains how the conversion reaction occurs in α-MnO 2 material. This study provides a comprehensive understanding of the electrochemical lithiation process and mechanism of α-MnO 2 material, in addition to the introduction of an improved in situ TEM biasing technique.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Effects of unburnt reaction progress on stretch flame dynamics under elevated temperatures

Flames in practical combustors and engines are inevitably subject to the effects of upstream chemical reaction progress and stretch, induced by elevated thermodynamic conditions and flow non-uniformities, respectively. Recent shock tube experiments and simulation studies on flame propagation have shown that flame speed under engine-relevant conditions can be enhanced with non-negligible upstream chemical reaction progress, especially when low-temperature heat release is involved in the unburnt mixture. On the other hand, depending on the mixture equivalence ratio and diluents, nonequidiffusion (including the non-unity Lewis number effect and the preferential diffusion effect) can couple with the flame stretch to fundamentally affect the flame propagation, which can manifest as either facilitation or suppression. Depending on the transport property of a reacting mixture, there hence can be either inhibition or promotion from unburnt reaction progress and stretch effects on flame propagation. In the current work, through one-dimensional numerical simulations of transient planar and spherical flames of n-heptane/air under elevated thermodynamic conditions, the combined effects of upstream chemical reaction progress and stretch on flame propagation are investigated. Results show that for both lean and rich n-heptane/air mixtures, flame speed can be substantially promoted with reaction progress, while the rich mixtures can exhibit opposite stretch dependence after first-stage ignition. Different definitions of effective Lewis number are adopted to explain the change in Markstein length for spherical flame in lean and rich mixtures reformed by the low temperature chemistry. Further, this work fills an important gap in laminar premixed flame research relevant to practical combustion systems and can provide useful insight into local turbulent flame behaviors and phi-sensitivity in engine combustion.

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

Modeling transient edge plasma transport with dynamic recycling

The work presents numerical simulation studies of the role that dynamic plasma recycling on the main wall and divertor target surfaces plays in transient edge plasma transport phenomena, such as edge localized modes (ELMs). The studies are performed by coupling the edge plasma transport code UEDGE [Rognlien et al., J. Nucl. Mater. 196–198, 347 (1992)] and the wall reaction–diffusion transport code FACE [Smirnov et al., Fusion Sci. Technol. 71, 75 (2017)]. The two-dimensional, time-dependent, two-way coupling of the codes, in a realistic tokamak geometry, is accomplished using the Integrated Plasma Simulator framework [Elwasif et al., in 18th Euromicro Conference on Parallel, Distributed and Network-Based Processing (PDP 2010), Pisa, Italy (IEEE, 2010), pp. 419–427] for all modeled material plasma boundaries. The simulations show that dynamic plasma recycling has substantially different characteristics on the main wall and on the divertor plates. It is demonstrated that during an ELM cycle the outer wall can dynamically absorb and release a number of particles comparable to that expelled by the ELM from the core plasma, by far exceeding the dynamic retention capacity of the divertor surfaces. The resulting evolution of the edge and divertor plasma conditions during an ELM cycle is analyzed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Designing a validation experiment for radio frequency condensation

Abstract Theoretical studies have suggested that nonlinear effects can lead to ‘radio frequency (RF) condensation’, where an initially broad current profile can coalesce in islands when they reach sufficient width. In suitable conditions, RF condensation can ‘self-focus’ the driven current to the center of an island, improving stabilization efficiency and reducing control complexity. In unsuitable conditions, the effect can prematurely deplete the RF energy before it reaches the island center, impairing stabilization. It is predicted that the RF condensation effect can significantly impact reactor-scale tokamaks. This paper presents a set of simulations investigating the conditions under which RF condensation might be encountered in present-day tokamaks. For concreteness, the calculations use equilibrium reconstructions for two shots from DIII-D and AUG. The Current Condensation Amid Magnetic Islands (OCCAMI) simulation code has been used for this investigation. The code takes as its input a numerically specified axisymmetric EFIT equilibrium solution, and it perturbatively constructs a 3D field with an island embedded at the appropriate rational surface. In the OCCAMI code, the GENRAY code is used for ray tracing and for calculating the power deposition along a ray trajectory, and GENRAY is coupled self-consistently to a solution of the thermal diffusion equation in the island. The simulation results described in the paper illuminate the conditions required for experimental validation of the theory of RF condensation. The simulations also provide an explanation of why the effect was not noticed in experiments prior to the publication of theoretical papers on the subject.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Application of Scalar Filtered Density Function to Turbulent Flows Under Supercritical Condition

Abstract The scalar filtered density function (FDF) methodology is extended and employed for large eddy simulation (LES) of turbulent flows under supercritical condition. To describe real fluid behavior, the extended methodology incorporates the generalized heat and mass diffusion models along with real fluid thermodynamic relations which are derived using the cubic Peng–Robinson equation of state. These models are implemented within the stochastic differential equations comprising the scalar FDF transport. Simulations are conducted of a temporally developing mixing layer under supercritical condition and the results are assessed by comparing with data generated by direct numerical simulation (DNS) of the same layer. The consistency of the proposed FDF methodology is assessed. The LES-FDF predictions are shown to agree favorably with the DNS data and exhibit several key features pertaining to supercritical turbulent flows.

Energy & Fuels↗
DOE OSTI2018

A Low-Rank Multigrid Method for the Stochastic Steady-State Diffusion Problem

We study a multigrid method for solving large linear systems of equations with tensor product structure here. Such systems are obtained from stochastic finite element discretization of stochastic partial differential equations such as the steady-state diffusion problem with random coefficients. When the variance in the problem is not too large, the solution can be well approximated by a low-rank object. In the proposed multigrid algorithm, the matrix iterates are truncated to low rank to reduce memory requirements and computational effort. The technique is proved convergent with an analytic error bound. Numerical experiments show its effectiveness in solving the Galerkin systems compared to the original multigrid solver, especially when the number of degrees of freedom associated with the spatial discretization is large.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Updated Phase Calculation Method for Metallic Nuclear Fuel in BISON

Zirconium based metallic nuclear fuel undergoes significant redistribution of its constituents, resulting in several different phases present within a single fuel pin. Primarily dominated by the movement of zirconium in U-Pu-Zr and UZr fuels, phase changes result in a three (or more) zoned fuel structure clearly visible during post-irradiation examination (PIE). Each of these regions may posses differing material properties, requiring knowledge of the relative concentration of each phase in order for fuel performance codes to capture the relevant fuel behavior. The model described here was implemented into the fuel performance code BISON, similar to previous work for modeling zirconium redistribution in U-Pu-Zr fuel. However, the modules have since been updated with several key assumptions for easier use in future zirconium redistribution work and other phase-dependent material properties implementation, as discussed in the following sections. A new material model has been implemented in BISON to calculate the fractional phase concentrations as a function of temperature and zirconium concentration. Instead of smoothing the material properties to ensure numerical convergence, the new model smoothes the phase contributions directly, allowing easier implementation of material properties. By comparing the zirconium distribution for different smoothing parameters, it was shown that the smoothing parameters utilized here was adequate for the parameters of interest, while decreasing computational expense. The calculated phase information has been used to determine phase-dependent diffusivities, and is available for other phase-dependent properties. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2013

STELLAR DYNAMOS AND CYCLES FROM NUMERICAL SIMULATIONS OF CONVECTION

We present a series of kinematic axisymmetric mean-field αΩ dynamo models applicable to solar-type stars, for 20 distinct combinations of rotation rates and luminosities. The internal differential rotation and kinetic helicity profiles required to calculate source terms in these dynamo models are extracted from a corresponding series of global three-dimensional hydrodynamical simulations of solar/stellar convection, so that the resulting dynamo models end up involving only one free parameter, namely, the turbulent magnetic diffusivity in the convecting layers. Even though the αΩ dynamo solutions exhibit a broad range of morphologies, and sometimes even double cycles, these models manage to reproduce relatively well the observationally inferred relationship between cycle period and rotation rate. On the other hand, they fail in capturing the observed increase of magnetic activity levels with rotation rate. This failure is due to our use of a simple algebraic α-quenching formula as the sole amplitude-limiting nonlinearity. This suggests that α-quenching is not the primary mechanism setting the amplitude of stellar magnetic cycles, with magnetic reaction on large-scale flows emerging as the more likely candidate. This inference is coherent with analyses of various recent global magnetohydrodynamical simulations of solar/stellar convection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

A mechanistic modeling framework for gas-phase adsorption kinetics and fixed-bed transport

Adsorption is a complex physicochemical process involving interparticle transport, interphase mass-transfer, intraparticle diffusion, and surface reactions. Although the exact description of the adsorption process will inevitably vary from system to system, it will always be governed by those primary mechanisms. Thus, by devising a model framework that can inherently include those mechanisms, it would be possible to create a modeling platform on which many different adsorption problems could be solved numerically. In order to accomplish this task, a generalized 1-D conservation law model was created to include the necessary mechanisms of adsorption on several different geometrical domains. Specific model applications for adsorption were developed under that framework and validated using experimental data available in literature or obtained in this work. This modeling platform makes it easier to model various adsorption problems and develop new adsorption models because of the common treatment of the mathematics governing the physical processes.

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

A splice method for local-to–nonlocal coupling of weak forms

Here, we propose a method to couple local and nonlocal diffusion models. By inheriting desirable properties such as patch tests, asymptotic compatibility and unintrusiveness from related splice and optimization-based coupling schemes, it enables the use of weak (or variational) formulations, is computationally efficient and straightforward to implement. We prove well-posedness of the coupling scheme and demonstrate its properties and effectiveness in a variety of numerical examples.

97 MATHEMATICS AND COMPUTING↗