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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

OpenSn: A massively parallel, open-source simulation environment for discrete ordinates radiation transport

OpenSn is an open-source, massively parallel deterministic radiation transport code for solving the discrete-ordinates ( S N ) form of the Boltzmann transport equation on unstructured, arbitrary polyhedral meshes. It supports high-fidelity simulations involving steady-state, eigenvalue, and adjoint problems for neutral particles (e.g., neutrons, photons, multi-particles), using the multigroup approximation in energy. OpenSn combines angular discretization via discrete ordinates with a discontinuous Galerkin finite element method (DGFEM) in space, enabling accurate resolution of transport physics on arbitrary polyhedral cells, included locally refined spatial grids. It includes multiple angular quadrature types, including locally refined angular quadratures. Written in modern C++ with a Python API, OpenSn runs efficiently on platforms ranging from laptops to supercomputers. The transport sweep algorithm is implemented using a task-based, directed-acyclic-graph (DAG) approach for each angle and supports asynchronous parallelism across thousands of MPI ranks. Group-set aggregation improves compute intensity, and synthetic acceleration techniques (e.g., diffusion synthetic acceleration, second-moment method) enhance solver convergence. OpenSn has been verified on reactor physics problems and demonstrated excellent weak and strong scaling performance on more than 32,768 processes, making it a versatile and robust platform for large-scale transport simulations in complex geometries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Synthesizing High-Capacity Oxyfluoride Conversion Anodes by Direct Fluorination of Molybdenum Dioxide (MoO 2 )

High-capacity metal oxide conversion anodes for lithium-ion batteries (LIBs) are primarily limited by their poor reversibility and cycling stability. In this study, a promising approach has been developed to improve the electrochemical performance of a MoO 2 anode by direct fluorination of the prelithiated MoO 2 . The fluorinated anode contains a mixture of crystalline MoO 2 and amorphous molybdenum oxyfluoride phases, as determined from a suite of characterization methods including X-ray diffraction, Raman spectroscopy, and X-ray photoelectron spectroscopy, and scanning transmission electron microscopy. Electrochemical measurements indicate that fluorination facilitates the conversion reaction kinetics, which leads to increased capacity, higher coulombic efficiency, and better cycling stability as compared to the nonfluorinated samples. Overall, these results suggest that fluorination after prelithiation not only favors formation of the oxyfluoride phase but also improves the lithium-ion diffusivity and reversibility of the conversion reaction, making it an attractive approach to address the problems of conversion electrodes. These findings provide a new route to design high-capacity negative electrodes for LIBs.

25 ENERGY STORAGE↗
DOE OSTI2022

Verification of Bison fission product species conservation under TRISO reactor conditions

When assessing the reliability and predictive capabilities of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. During this process, numerical results in a discrete solution are compared to the analytical solution of the mathematical model. Here, in this paper, the code verification process is applied to one-dimensional spatiotemporal problems that exercise partial differential equation governing the conservation of fission product species (or mass diffusion). Numerical experiments were performed in the Bison fuel performance code to evaluate its predictive capability under various TRISO reactor conditions such as base irradiation and safety heating test conditions for either short- or long-lived fission product species, as well as a case concerning evaporation from the outer surface of a particle. The code predictions were compared with the expected exact results obtained from the analytical expressions, and the fact that they demonstrate the correct analytical behavior provides strong evidence of proper numerical algorithm implementation.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2019

Observing crystal nucleation in four dimensions using atomic electron tomography

Nucleation plays a critical role in many physical and biological phenomena that range from crystallization, melting and evaporation to the formation of clouds and the initiation of neurodegenerative diseases. However, nucleation is a challenging process to study in experiments, especially in the early stage when several atoms or molecules start to form a new phase from a parent phase. A number of experimental and computational methods have been used to investigate nucleation processes, but experimental determination of the three-dimensional atomic structure and the dynamics of early stage nuclei has been unachievable. Here we use atomic electron tomography to study early stage nucleation in four dimensions (4D: that is, including time) at atomic resolution. Using FePt nanoparticles as a model system, we find that early stage nuclei are irregularly shaped, each has a core of one to a few atoms with the maximum order parameter, and the order parameter gradient points from the core to the boundary of the nucleus. We capture the structure and dynamics of the same nuclei undergoing growth, fluctuation, dissolution, merging and/or division, which are regulated by the order parameter distribution and its gradient. These experimental observations are corroborated by molecular dynamics simulations of heterogeneous and homogeneous nucleation in liquid–solid phase transitions of Pt. Our experimental and molecular dynamics results differ from classical nucleation theory, indicating that a theory beyond this is needed to describe early stage nucleation at the atomic scale. Looking forward, we anticipate that the reported approach will open the door to the study of many fundamental problems in materials science, nanoscience, condensed matter physics and chemistry, such as phase transition, atomic diffusion, grain boundary dynamics, interface motion, defect dynamics and surface reconstruction with 4D atomic resolution.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2015

Quantum Monte Carlo calculation of the binding energy of the beryllium dimer

The accurate calculation of the binding energy of the beryllium dimer is a challenging theoretical problem. In this study, the binding energy of Be{sub 2} is calculated using the diffusion Monte Carlo (DMC) method, using single Slater determinant and multiconfigurational trial functions. DMC calculations using single-determinant trial wave functions of orbitals obtained from density functional theory calculations overestimate the binding energy, while DMC calculations using Hartree-Fock or CAS(4,8), complete active space trial functions significantly underestimate the binding energy. In order to obtain an accurate value of the binding energy of Be{sub 2} from DMC calculations, it is necessary to employ trial functions that include excitations outside the valence space. Our best estimate DMC result for the binding energy of Be{sub 2}, obtained by using configuration interaction trial functions and extrapolating in the threshold for the configurations retained in the trial function, is 908 cm{sup −1}, only slightly below the 935 cm{sup −1} value derived from experiment.

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

Verification of MOOSE/Bison's Heat Conduction Solver Using Combined Spatiotemporal Convergence Analysis

Bison is a computational physics code that uses the finite element method to model the thermo-mechanical response of nuclear fuel. Since Bison is used to inform high-consequence decisions, it is important that its computational results are reliable and predictive. One important step in assessing the reliability and predictive capabilities of a simulation tool is the verification process, which quantifies numerical errors in a discrete solution relative to the exact solution of the mathematical model. One step in the verification process—called code verification—ensures that the implemented numerical algorithm is a faithful representation of the underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. In this paper, the code verification process is applied to spatiotemporal heat conduction problems in Bison. Simultaneous refinement of the discretization in space and time is employed to reveal any potential mistakes in the numerical algorithms for the interactions between the spatial and temporal components of the solution. For each verification problem, the correct spatial and temporal order of accuracy is demonstrated for both first- and second-order accurate finite elements and a variety of time-integration schemes. Furthermore, these results provide strong evidence that the Bison numerical algorithm for solving spatiotemporal problems reliably represents the underlying mathematical model in MOOSE. The selected test problems can also be used in other simulation tools that numerically solve for conduction or diffusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Rapid Simulation of Unconventional Reservoirs by Multidomain Multiresolution Modeling Based on the Diffusive Time of Flight

Summary The fast marching method (FMM)-based rapid flow simulation has been shown to accelerate simulation efficiency by orders of magnitude by transforming 3D simulation to equivalent 1D simulation using the concept of the “diffusive time of flight” (DTOF). However, the 1D transformation does not directly apply to multiwell problems. In this paper, we propose a novel DTOF-based multidomain multiresolution discretization scheme to accelerate multiwell simulation of unconventional reservoirs. Our method formulates multiwell simulation problems based on the DTOF which displays the pressure front propagation in unconventional reservoirs. The DTOF contours are used to partition the reservoir into local and shared domains. A local domain is where the flow is dominated by a single well, and the shared domain is where the fluid flow is influenced by multiple wells. The DTOF contours expand independently in local domains and interfere in the shared domain. After the partitioning, each domain is discretized using a multiresolution scheme whereby the original 3D fine mesh is preserved near the wells to account for detailed physics including gravity, and the rest of the domain is discretized into 1D mesh based on the DTOF contours to alleviate the simulation workload. The power and efficacy of our approach are demonstrated using synthetic and field-scale simulation models with different degrees of geologic and well-completion complexity. The simulation results, number of active cells, and computation time for the proposed discretization scheme are compared with the original high-fidelity 3D model for each case. The results show that the proposed method is suitable for multiwell simulation problems in unconventional reservoirs and can accelerate flow simulations by orders of magnitude with minimal loss of accuracy. The novelty of this work is the creation of DTOF-derived multiresolution discretization with local and shared domains to simplify and accelerate the calculation of subsurface flow problems, especially in unconventional reservoirs. Our workflow can be easily interfaced with commercial simulators, making it suitable for large-scale field applications.

Engineering↗
NASA NTRS1974

Technology application at Rockwell International

Technology diffusion at Rockwell International has progressed from the application of special skills to solve specific problems to the transfer of the required people and skills to allow the commercial divisions to solve their own problems and develop or improve their own products. Our prime effort is concentrated on commercial industrial applications. Rockwell's major emphasis on advanced technology utilization is directed through three operational modes, namely, transfer from technologically developed to underdeveloped organizations within the company; transfer between technologically developed organizations to form new operations within the company; and formation of high technology spin-off organizations beyond the corporate entity.

Meechan, C. J.↗
NASA NTRS1982

Silicon solar cell process development, fabrication and analysis

For UCP Si, randomly selected wafers and wafers cut from two specific ingots were studied. For the randomly selected wafers, a moderate gettering diffusion had little effect. Moreover, an efficiency up to 14% AMI was achieved with advanced processes. For the two specific UCP ingots, ingot #5848-13C displayed severe impurity effects as shown by lower 3sc in the middle of the ingot and low CFF in the top of the ingot. Also the middle portions of this ingot responded to a series of progressively more severe gettering diffusion. Unexplained was the fact that severely gettered samples of this ingot displayed a negative light biased effect on the minority carrier diffusion length while the nongettered or moderately gettered ones had the more conventional positive light biased effect on diffusion length. On the other hand, ingot C-4-21A did not have the problem of ingot 5848-13C and behaved like to the randomly selected wafers. The top half of the ingot was shown to be slightly superior to the bottom half, but moderate gettering helped to narrow the gap.

Iles, P. A.↗
NASA NTRS1985

Multigrid techniques for the solution of the passive scalar advection-diffusion equation

The solution of elliptic passive scalar advection-diffusion equations is required in the analysis of many turbulent flow and convective heat transfer problems. The accuracy of the solution may be affected by the presence of regions containing large gradients of the dependent variables. The multigrid concept of local grid refinement is a method for improving the accuracy of the calculations in these problems. In combination with the multilevel acceleration techniques, an accurate and efficient computational procedure is developed. In addition, a robust implementation of the QUICK finite-difference scheme is described. Calculations of a test problem are presented to quantitatively demonstrate the advantages of the multilevel-multigrid method.

Phillips, R. E.↗
NASA NTRS1985

Excimer laser annealing for fabrication of low-cost solar cells

Pulsed excimer laser annealing was successfully performed using a 50 w laser. Both polished and texturized cells were tried, however, there are serious problems with nonuniformity on texturized cells. A number of cells were produced and compared to diffusion furnace annealed cells. There was no clear economic advantage in using an excimer laser and there was a small penalty on average efficiency. The conclusion was that the excimer laser anneal process must be able to produce superior cells to be considered as a viable process option.

Spitzer, M. B.↗
NASA NTRS1987

Optimal cure cycle design of a resin-fiber composite laminate

A unified computed aided design method was studied for the cure cycle design that incorporates an optimal design technique with the analytical model of a composite cure process. The preliminary results of using this proposed method for optimal cure cycle design are reported and discussed. The cure process of interest is the compression molding of a polyester which is described by a diffusion reaction system. The finite element method is employed to convert the initial boundary value problem into a set of first order differential equations which are solved simultaneously by the DE program. The equations for thermal design sensitivities are derived by using the direct differentiation method and are solved by the DE program. A recursive quadratic programming algorithm with an active set strategy called a linearization method is used to optimally design the cure cycle, subjected to the given design performance requirements. The difficulty of casting the cure cycle design process into a proper mathematical form is recognized. Various optimal design problems are formulated to address theses aspects. The optimal solutions of these formulations are compared and discussed.

Hou, Jean W.↗
NASA NTRS1989

Numerical simulations in combustion

This paper reviews numerical simulations in reacting flows in general and combustion phenomena in particular. It is shown that use of implicit schemes and/or adaptive mesh strategies can improve convergence, stability, and accuracy of the solution. Difficulties increase as turbulence and multidimensions are considered, particularly when finite-rate chemistry governs the given combustion problem. Particular attention is given to the areas of solid-propellant combustion dynamics, turbulent diffusion flames, and spray droplet vaporization.

Chung, T. J.↗
NASA NTRS1990

Is there a classical Hyades lithium problem?

The distribution of surface lithium abundances in the Hyades cluster, as a function of effective temperature, has resisted a satisfactory theoretical explanation for more than two decades. Trends in the historically recomputed values of stellar opacity have stimulated a systematic investigation of the effects of larger opacities on lithium depletion during the premain-sequence evolution of low-mass stars. It is shown that the calculated depletion depends quite strongly on the opacities employed, and that both surface and interior opacities play an important role. By performing standard evolutionary calculations with parameterized increments in the opacities, the observed Hyades Li and T(eff) relation for G and K dwarfs have been determined in a manner consistent with the other observational constraints. The solution presented herein requires an increase in interior opacities at about 4 million K over older values by approximately less than 37 percent. Such a change is not incompatible with the historical trends in opacity computations and suggests that the resolution of this problem lies solely in the opacities, without recourse to other physical mechanisms such as turbulent diffusion or 'extra mixing'.

Swenson, Fritz J.↗
NASA NTRS1992

Process modelling for materials preparation experiments

The development is examined of mathematical tools and measurement of transport properties necessary for high fidelity modeling of crystal growth from the melt and solution, in particular for the Bridgman-Stockbarger growth of mercury cadmium telluride (MCT) and the solution growth of triglycine sulphate (TGS). The tasks include development of a spectral code for moving boundary problems, kinematic viscosity measurements on liquid MCT at temperatures close to the melting point, and diffusivity measurements on concentrated and supersaturated TGS solutions. A detailed description is given of the work performed for these tasks, together with a summary of the resulting publications and presentations.

Rosenberger, Franz↗
NASA NTRS2011

Analysis of the Radiative Transfer Equation with Highly Asymmetric Phase Function

This paper considers a scalar radiative transfer problem with high scattering anisotropy, Two computational methods are presented based on decomposition of the diffuse light field into a regular and anisotropic part. The first algorithm (DOMAS) singles out the anisotropic radiance in the forward scattering peak using the Small-Angle Modification of RTE. The second algorithm (DOM2+) separates the single scattering radiance as an anisotropic part, which largely defines the fine detail of the total radiance in the backscattering directions. In both cases, the anisotropic part is represented analytically. With anisotropy subtraction, the regular part of the signal. which requires a numerical solution, is essentially smoothed as a function of angles. Further, the transport equation is obtained for the regular part that contains an additional source function from the anisotropic part of the signal. This equation is solved with the discrete ordinates method. A conducted numerical analysis of this work showed that algorithm DOMAS has a strong advantage as compared to the standard discrete ordinates method for simulation of the radiance transmission, and DOM2 + is the best of the three for the reflection computations. Both algorithms offer at least a factor of three acceleration of convergence of the azimuthal series for highly anisotropic phase functions.

Korkin, Sergey V.↗
NASA NTRS1987

An unconditionally-stable central differencing scheme for high Reynolds number flows

The central difference scheme (CDS) is a second order accurate scheme which is free of numerical diffusion (in the second order sense) and is simple to implement: However, for grid Peclet numbers larger than 2, the CDS leads to over- and undershoots and is unstable. The present paper describes a method, called CONDIF, which retains the essential nature of the CDS but eliminates the over- and under-shoots. It leads to unconditionally positive coefficients and, in the limit, approaches the CDS for all values of grid Peclet numbers. The CONDIF modifies the CDS by introducing a controlled amount of numerical diffusion based on the local gradients. In the worst case the scheme yields results similar to those of the hybrid scheme. This paper reports the results obtained from CONDIF for a number of test problems which have been widely used for comparative study of numerical schemes in the published literature. For most of these problems, the CONDIF results are significantly more accurate than the hybrid scheme at high Peclet numbers. In particular, the CONDIF scheme depicts much lower level of numerical diffusion than the hybrid scheme even when the Peclet number is very high and the flow is at large angles to the grid.

Runchal, Akshai K.↗
DOE OSTI2014

The exit-time problem for a Markov jump process

The purpose of this paper is to consider the exit-time problem for a finite-range Markov jump process, i.e, the distance the particle can jump is bounded independent of its location. Additionally, such jump diffusions are expedient models for anomalous transport exhibiting super-diffusion or nonstandard normal diffusion. We refer to the associated deterministic equation as a volume-constrained nonlocal diffusion equation. The volume constraint is the nonlocal analogue of a boundary condition necessary to demonstrate that the nonlocal diffusion equation is well-posed and is consistent with the jump process. A critical aspect of the analysis is a variational formulation and a recently developed nonlocal vector calculus. Finally, this calculus allows us to pose nonlocal backward and forward Kolmogorov equations, the former equation granting the various moments of the exit-time distribution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗