Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion problems”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

A low-communication, parallel algorithm for solving PDEs based on range decomposition

This work proposes a new, low-communication algorithm for solving PDEs on massively parallel computers. The range decomposition (RD) algorithm exposes coarse-grain parallelism by applying nested iteration and adaptive mesh refinement locally before performing a global communication step. Just a few such steps are observed to be sufficient to obtain accuracy within a small multiple of discretization error. The target applications are petascale and exascale machines, where hierarchical parallelism is required and traditional parallel numerical PDE communication patterns are costly because of message latency. The RD algorithm uses a partition of unity to equally distribute the error, and thus, the work. The computational advantages of this approach are that the decomposed problems can be solved in parallel without any communication until the partitioned solutions are summed. This offers potential advantages in the paradigm of expensive communication but very cheap computation. This paper introduces the method and explains the details of the communication step. Two performance models are developed, showing that the latency cost associated with a traditional parallel implementation of nested iteration is proportional to log(P) 2 , whereas the RD method reduces the communication latency to log(P), while maintaining similar bandwidth costs. Numerical results for two problems, Laplace and advection diffusion, demonstrate the enhanced performance, and a heuristic argument explains why the method converges quickly.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

On the fractional Harry Dym equation

We present and discuss the general case of a fractional nonlinear partial differential equation using similarity reductions and recover results associated with Harry Dym-type equations. Particular cases are discussed and some plots illustrate the results. The fractional derivative is considered in the Caputo sense.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Fast multiscale contrast independent preconditioners for linear elastic topology optimization problems

The goal of this work is to present a fast and viable approach for the numerical solution of the high-contrast state problems arising in topology optimization. The optimization process is iterative, and the gradients are obtained by an adjoint analysis, which requires the numerical solution of large high-contrast linear elastic problems with features spanning several length scales. The size of the discretized problems forces the utilization of iterative linear solvers with solution time dependent on the quality of the preconditioner. The lack of clear separation between the scales, as well as the high-contrast, imposes severe challenges on the standard preconditioning techniques. Thus, here we propose new methods for the high-contrast elasticity equation with performance independent of the high-contrast and the multi-scale structure of the elasticity problem. The solvers are based on two-levels domain decomposition techniques with a carefully constructed coarse level to deal with the high-contrast and multi-scale nature of the problem. The construction utilizes spectral equivalence between scalar diffusion and each displacement block of the elasticity problems and, in contrast to previous solutions proposed in the literature, is able to select the appropriate dimension of the coarse space automatically. The new methods inherit the advantages of domain decomposition techniques, such as easy parallelization and scalability. Finally, the presented numerical experiments demonstrate the excellent performance of the proposed methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A continuum dislocation-based model of wedge microindentation of single crystals

Recent Electron Backscatter Diffraction (EBSD) experiments have revealed the emergence of heterogeneous dislocation microstructures forming under a wedge indenter in fcc crystals, where micro-meter dislocation patterns challenge the predictions of traditional models of plasticity. In order to explain the formation of these features and develop a relationship between the force-displacement curve and the dislocation substructure, here we present a model of wedge indentation based on the continuum theory of dislocations. The model accounts for large deformation kinematics through the multiplicative split of the deformation gradient tensor, where the incompatible plastic component of deformation results from the flux of dislocations on different and interacting slips systems. Constitutive equations for dislocation fluxes are determined from a dissipative variational principle. As a result, each dislocation density satisfies an initial-boundary value problem with convective-diffusive character, which is coupled to the macroscopic stress and displacement fields governing the deformation process. Solution to the self-consistent continuum formulation is found using the finite element method. Computer simulations mimic the experimental conditions of wedge micro-indentation experiments into Ni single-crystals used by Kysar et al. (2010a). A comparison of overall dislocation density distribution and macroscopic mechanical response shows good overall agreement with the experimental results in terms of the detailed features of dislocation patterns and lattice rotations as well as the macroscopic force-displacement response.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Energy and momentum preserving Coulomb collision model for kinetic Monte Carlo simulations of plasma steady states in toroidal fusion devices

A kinetic Monte Carlo model suited for self-consistent transport studies is proposed and tested. The Monte Carlo collision operator is based on a widely used model of Coulomb scattering by a drifting Maxwellian and a new algorithm enforcing the momentum and energy conservation laws. The difference to other approaches consists in a specific procedure of calculating the background Maxwellian parameters, which does not require ensemble averaging and, therefore, allows for the use of single-particle algorithms. This possibility is useful in transport balance (steady state) problems with a phenomenological diffusive ansatz for the turbulent transport, because it allows a direct use of variance reduction methods well suited for single particle algorithms. In addition, a method for the self-consistent calculation of the electric field is discussed. Results of testing of the new collision operator using a set of 1D examples, and preliminary results of 2D modelling in realistic tokamak geometry, are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

On the use of reverse Brownian motion to accelerate hybrid simulations

Multiscale and multiphysics simulations are two rapidly developing fields of scientific computing. Efficient coupling of continuum (deterministic or stochastic) constitutive solvers with their discrete (stochastic, particle-based) counterparts is a common challenge in both kinds of simulations. We focus on interfacial, tightly coupled simulations of diffusion that combine continuum and particle-based solvers. The latter employs the reverse Brownian motion (rBm), a Monte Carlo approach that allows one to enforce inhomogeneous Dirichlet, Neumann, or Robin boundary conditions and is trivially parallelizable. We discuss numerical approaches for improving the accuracy of rBm in the presence of inhomogeneous Neumann boundary conditions and alternative strategies for coupling the rBm solver with its continuum counterpart. Numerical experiments are used to investigate the convergence, stability, and computational efficiency of the proposed hybrid algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Parametric dynamic mode decomposition for reduced order modeling

Dynamic Mode Decomposition (DMD) is a model-order reduction approach, whereby spatial modes of fixed temporal frequencies are extracted from numerical or experimental data sets. The DMD low-rank or reduced operator is typically obtained by singular value decomposition of the temporal data sets. For parameter-dependent models, as found in many multi-query applications such as uncertainty quantification or design optimization, the only parametric DMD technique developed was a stacked approach, with data sets at multiple parameter values were aggregated together, increasing the computational work needed to devise low-rank dynamical reduced-order models. Here in this paper, we present two novel approach to carry out parametric DMD: one based on the interpolation of the reduced-order DMD eigen-pair and the other based on the interpolation of the reduced DMD (Koopman) operator. Numerical results are presented for diffusion-dominated nonlinear dynamical problems, including a multiphysics radiative transfer example. All three parametric DMD approaches are compared.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Regularizing INR with Diffusion Prior for Self-Supervised 3D Reconstruction OF Neutron Computed Tomography Data

Recently, generative diffusion priors have made huge strides as inverse problem solvers, including the ability to be adapted for inference on out-of-distribution data. Concurrently, implicit neural representations (INRs) have emerged as fast and lightweight inverse imaging solvers that are amenable to hybrid approaches that combine learned priors with traditional inverse problem formulations. In this paper, we present a diffusive computed tomography (CT) inversion framework for regularizing INRs called Diffusive INR (DINR), designed to enable high-quality reconstruction from sparse-view neutron CT. Pretrained purely on synthetic data, DINR is evaluated on simulated and experimentally obtained observations of concrete microstructures, where traditional reconstruction methods suffer substantial degradation when the number of views is reduced. Our approach delivers superior performance, reduces reconstruction artifacts, and achieves gains in PSNR and SSIM, enabling accurate micro-structural characterization even under extreme data limitations compared to state-of-the-art sparse-view reconstruction techniques.

Hossain, Maliha [ORNL]↗
DOE OSTI2018

An Adaptive Reduced Basis Collocation Method Based on PCM ANOVA Decomposition for Anisotropic Stochastic PDEs

The combination of reduced basis and collocation methods enables efficient and accurate evaluation of the solutions to parameterized PDEs. Here, we study the stochastic collocation methods that can be combined with reduced basis methods to solve high-dimensional parameterized stochastic PDEs. We also propose an adaptive algorithm using a probabilistic collocation method (PCM) and ANOVA decomposition. This procedure involves two stages. First, the method employs an ANOVA decomposition to identify the effective dimensions, i.e., subspaces of the parameter space in which the contributions to the solution are larger, and sort the reduced basis solution in a descending order of error.Then, the adaptive search refines the parametric space by increasing the order of polynomials until the algorithm is terminated by a saturation constraint. We demonstrate the effectiveness of the proposed algorithm for solving a stationary stochastic convection-diffusion equation, a benchmark problem chosen because solutions contain steep boundary layers and anisotropic features. Lastly, we show that two stages of adaptivity are critical in a benchmark problem with anisotropic stochasticity.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1977

Systems improved numerical differencing analyzer

Analyzer is software system that possesses capabilities which make it well suited for solving lumped-parameter representations of physical problems governed by diffusion-type equations, such as Fourier, Poisson, or La Place.

Fink, L. C.↗
NASA NTRS1980

Systems improved numerical differencing analyzer

Program solves physical problems governed by diffusion-type equations, provided that equations can be modeled by lumped-parameter representation. Program is used for thermal analysis, and could be adapted to solve Fourier, Poisson, and Laplace differential equations. Program is in FORTRAN IV and Assembler for execution on UNIVAC 1100-series or CYBER 175.

Source record↗
NASA NTRS1985

Spatial distribution of cosmic rays in magnetic cycles of the Sun

The space-time distribution of cosmic ray (CR) density is determined from the solution of the anisotropic diffusion equation. The boundary problem is solved in a quasistationary axisymmetric case, taking into account the delay of cosmic rays relative to the cyclic variations of solar activity. The specificty of the modulation model is the use of the empirical dependence of the free path of particles for scattering on the observed latitudinal distribution of solar activity. Another specific feature of the model is an account of an additional effect due to the drift of charged particles towards the helio-equator plane or backward, depending on the direction of the total magnetic field of the Sun.

Alania, M. V.↗
NASA NTRS1987

Excimer laser annealing: A gold process for CZ silicon junction formation

A cold process using an excimer laser for junction formation in silicon has been evaluated as a way to avoid problems associated with thermal diffusion. Conventional thermal diffusion can cause bulk precipitation of SiOx and SiC or fail to completely activate the dopant, leaving a degenerate layer at the surface. Experiments were conducted to determine the feasibility of fabricating high quality p-n junctions using a pulsed excimer laser for junction formation at remelt temperature with ion-implanted surfaces. Solar-cell efficiency exceeding 16 percent was obtained using Czochralski single-crystal silicon without benefit of back surface field or surface passivation. Characterization shows that the formation of uniform, shallow junctions (approximately 0.25 micron) by excimer laser scanning preserves the minority carrier lifetime that leads to high current collection. However, the process is sensitive to initial surface conditions and handling parameters that drive the cost up.

Wong, David C.↗
NASA NTRS1988

Wave-particle interactions in the magnetosphere of Uranus

The Voyager 2 encounter of Uranus has provided observations of plasma waves in and near the magnetosphere. These data, while the first from Uranus, will also be the only direct information on wave-particle interactions at this planet for many years to come. The observations include electrostatic waves upstream of the bow shock, turbulence in the shock, Bernstein emissions and whistler mode waves in the magnetosphere, broadband electrostatic noise in the magnetotail, and a number of the other types of plasma waves which have yet to be clearly identified. Each of these types of waves exist in a plasma environment which both supports the growth of the waves and is modified by interactions with the waves. Wave-particle interactions provide the channels through which the waves can accelerate, scatter, or thermalize the plasmas. The most spectacular example in the case of Uranus is the extremely intense whistler mode activity in the inner magnetosphere which is the source of strong pitch angle diffusion. The resulting electron precipitation is sufficient to produce the auroral emissions observed by Voyager. The strong diffusion, however, presents the problem of supplying electrons in the range of 5 to 40 keV in order to support the losses to the atmosphere.

Kurth, W. S.↗
NASA NTRS1991

Wave-particle interactions in the magnetosphere of Uranus

The Voyager 2 encounter of Uranus has provided observations of plasma waves in and near the magnetosphere. These data, while the first from Uranus, will also be the only direct information on wave-particle interactions at this planet for many years to come. The observations include electrostatic waves upstream of the bow shock, turbulence in the shock Bernstein emissions and whistler mode waves in the magnetosphere, broadband electrostatic noise in the magnetotail, and a number of the other types of plasma waves which have yet to be clearly identified. Each of these types of waves exist in a plasma environment which both supports the growth of the waves and is modified by interactions with the waves. Wave-particle interactions provide the channels through which the waves can accelerate, scatter, or thermalize the plasmas. The most spectacular example in the case of Uranus is the extremely intense whistler mode activity in the inner magnetosphere which is the source of strong pitch angle diffusion. The resulting electron precipitation is sufficient to produce the auroral emissions observed by Voyager. The strong diffusion, however, presents the problem of supplying electrons in the range of 5 to 40 keV in order to support the losses to the atmosphere.

Kurth, W. S.↗
NASA NTRS1994

Counting particles emitted by stratospheric aircraft and measuring size of particles emitted by stratospheric aircraft

The ER-2 condensation nuclei counter (CNC) has been modified to reduce the diffusive losses of particles within the instrument. These changes have been successful in improving the counting efficiency of small particles at low pressures. Two techniques for measuring the size distributions of particles with diameters less than 0.17 micrometers have been evaluated. Both of these methods, the differential mobility analyzer (DMA) and the diffusion battery, have fundamental problems that limit their usefulness for stratospheric applications. We cannot recommend either for this application. Newly developed, alternative methods for measuring small particles include inertial separation with a low-loss critical orifice and thin-plate impactor device. This technique is now used to collect particles in the multisample aerosol collector housed in the ER-2 CNC-2, and shows some promise for particle size measurements when coupled with a CNC as a counting device. The modified focused-cavity aerosol spectrometer (FCAS) can determine the size distribution of particles with ambient diameters as small as about 0.07 micrometers. Data from this instrument indicates the presence of a nuclei mode when CNC-2 indicates high concentrations of particles, but cannot resolve important parameters of the distribution.

Wilson, James Charles↗
NASA NTRS1995

Energization in regions of CIRs unconnected to shocks are probably not the result of cross-field transport

Corotating energetic ion populations (CEIPs) associated with the forward and reverse shocks of corotating interaction regions (CIRs) are observed in CIRs at places where models say are magnetically unconnected to either shock. Such disconnections between CEIPs and shocks are common and have been documented with data from Pioneers 10 and 11 and confirmed with data from Ulysses. They pose a problem for models that account for these CEIPs in terms of ion energization at the shocks followed by ion propagation along field lines. Two possible resolutions to this problem have been suggested: diffusion of the ions across field lines and extension of the ion energization process to regions beyond the shock waves. Here we quantitatively examine the first of these possibilities. We give the Green's function solution to the convection-diffusion equation applied to idealized CIR geometry, with a source at the reverse shock -- the main producer of CEIPs. Two kinds of diffusion are considered: resonant diffusion and stochastic field line diffusion. We find that for resonant diffusion the computed ratio is many orders of magnitudes below the observed ratio. For stochastic field line diffusion, the computed ratio approximately equals the observed ratio if a diffusion coefficient appropriate to the free solar wind is used. It is several orders of magnitude below the observed ratio, however, if a diffusion coefficient appropriate to CIRs is used. We conclude that cross-field diffusion probably does not account for the presence of energetic ions in regions of CIRs that are magnetically unconnected to its shock waves. We suggest that the alternative possibility -- the energetic ions in regions magnetically unconnected to shocks result from an acceleration process that is independent of shocks -- be pursued to the point where quantitative tests can be performed.

Intriligator, Devrie S.↗