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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Kinetic modelling for zinc (II) ions biosorption onto Luffa cylindrica

The biosorption of Zinc (II) ions onto a biomaterial - Luffa cylindrica has been studied. This biomaterial was characterized by elemental analysis, surface area, pore size distribution, scanning electron microscopy, and the biomaterial before and after sorption, was characterized by Fourier Transform Infra Red (FTIR) spectrometer. The kinetic nonlinear models fitted were Pseudo-first order, Pseudo-second order and Intra-particle diffusion. A comparison of non-linear regression method in selecting the kinetic model was made. Four error functions, namely coefficient of determination (R{sup 2}), hybrid fractional error function (HYBRID), average relative error (ARE), and sum of the errors squared (ERRSQ), were used to predict the parameters of the kinetic models. The strength of this study is that a biomaterial with wide distribution particularly in the tropical world and which occurs as waste material could be put into effective utilization as a biosorbent to address a crucial environmental problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Competitive autocatalytic reactions in chaotic flows with diffusion: Prediction using finite-time Lyapunov exponents

We investigate chaotic advection and diffusion in autocatalytic reactions for time-periodic sine flow computationally using a mapping method with operator splitting. We specifically consider three different autocatalytic reaction schemes: a single autocatalytic reaction, competitive autocatalytic reactions, which can provide insight into problems of chiral symmetry breaking and homochirality, and competitive autocatalytic reactions with recycling. In competitive autocatalytic reactions, species B and C both undergo an autocatalytic reaction with species A such that A+B→2B and A+C→2C. Small amounts of initially spatially localized B and C and a large amount of spatially homogeneous A are advected by the velocity field, diffuse, and react until A is completely consumed and only B and C remain. We find that local finite-time Lyapunov exponents (FTLEs) can accurately predict the final average concentrations of B and C after the reaction completes. The species that starts in the region with the larger FTLE has, with high probability, the larger average concentration at the end of the reaction. If B and C start in regions with similar FTLEs, their average concentrations at the end of the reaction will also be similar. When a recycling reaction is added, the system evolves towards a single species state, with the FTLE often being useful in predicting which species fills the entire domain and which is depleted. The FTLE approach is also demonstrated for competitive autocatalytic reactions in journal bearing flow, an experimentally realizable flow that generates chaotic dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Convergence Study of LR-NDA Using Fourier Analysis

Numerical solution of neutron transport equations is challenging. It often requires acceleration techniques to improve computational efficiency. In the last two decades, various acceleration methods have been developed. Among those methods, the coarse-mesh nonlinear diffusion acceleration (CM-NDA) technique, or commonly known as coarse-mesh finite difference (CMFD), becomes one of the most popular methods due to its effectiveness and simplicity in recent years. However, CM-NDA will lose its effectiveness and could become unstable when the mesh optical thickness becomes large. There exists a number of stabilization techniques to improve the stability of CM-NDA, such as pCMFD and most recent odCMFD. Wang et al., proposed a new CM-NDA scheme, LR-NDA, which employs local refinement on the coarse mesh. The new scheme can significantly improve the effectiveness and stability of CM-NDA. In this paper, we perform a convergence study of LR-NDA using Fourier analysis. The spectral radius calculated from Fourier analysis is consistent with the numerical results for a 1-D fixed source model problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Excitation Energies of Localized Correlated Defects via Quantum Monte Carlo: A Case Study of Mn 4+ -Doped Phosphors

Accurate excitation energies of localized defects have been a long-standing problem for electronic structure calculation methods. Using Mn 4+ -doped solids as our proof of principle, we show that diffusion quantum Monte Carlo (DMC) is able to predict phosphorescence emission energies within statistical error. To demonstrate the generality of our DMC approach for other possible localized defects, we conduct charge density analyses using DMC and density functional theory (DFT). We also identify a new material with an emission energy of 1.97(8) eV, which is close to the optimum of 2.03 eV for a red-emitting phosphor. To our knowledge, our work is the first report on studying excitation energies of a transition metal impurity using an ab initio many-body electronic structure method. In contrast, semilocal and hybrid-DFT largely underestimates, and fails to reproduce, some of the trends in the emission energies. In conclusion, our work underscores the importance of an accurate account of exchange, correlation, and excitonic effects for localized excitations in defective solids.

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

ANALYTICAL MODELS OF EXOPLANETARY ATMOSPHERES. II. RADIATIVE TRANSFER VIA THE TWO-STREAM APPROXIMATION

We present a comprehensive analytical study of radiative transfer using the method of moments and include the effects of non-isotropic scattering in the coherent limit. Within this unified formalism, we derive the governing equations and solutions describing two-stream radiative transfer (which approximates the passage of radiation as a pair of outgoing and incoming fluxes), flux-limited diffusion (which describes radiative transfer in the deep interior), and solutions for the temperature-pressure profiles. Generally, the problem is mathematically underdetermined unless a set of closures (Eddington coefficients) is specified. We demonstrate that the hemispheric (or hemi-isotropic) closure naturally derives from the radiative transfer equation if energy conservation is obeyed, while the Eddington closure produces spurious enhancements of both reflected light and thermal emission. We concoct recipes for implementing two-stream radiative transfer in stand-alone numerical calculations and general circulation models. We use our two-stream solutions to construct toy models of the runaway greenhouse effect. We present a new solution for temperature-pressure profiles with a non-constant optical opacity and elucidate the effects of non-isotropic scattering in the optical and infrared. We derive generalized expressions for the spherical and Bond albedos and the photon deposition depth. We demonstrate that the value of the optical depth corresponding to the photosphere is not always 2/3 (Milne's solution) and depends on a combination of stellar irradiation, internal heat, and the properties of scattering in both the optical and infrared. Finally, we derive generalized expressions for the total, net, outgoing, and incoming fluxes in the convective regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Fluctuating hydrodynamics of multispecies nonreactive mixtures

In this study we discuss the formulation of the fluctuating Navier-Stokes equations for multispecies, nonreactive fluids. In particular, we establish a form suitable for numerical solution of the resulting stochastic partial differential equations. An accurate and efficient numerical scheme, based on our previous methods for single species and binary mixtures, is presented and tested at equilibrium as well as for a variety of nonequilibrium problems. These include the study of giant nonequilibrium concentration fluctuations in a ternary mixture in the presence of a diffusion barrier, the triggering of a Rayleigh-Taylor instability by diffusion in a four-species mixture, as well as reverse diffusion in a ternary mixture. Finally, good agreement with theory and experiment demonstrates that the formulation is robust and can serve as a useful tool in the study of thermal fluctuations for multispecies fluids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Analytical Sensitivity Analysis of a Spent Nuclear Fuel Cask

Here, we report nuclear science and engineering is a field increasingly dominated by computational studies resulting from increasingly powerful computational tools. As a result, analytical studies, which previously pioneered nuclear engineering, are increasingly viewed as secondary or unnecessary. However, analytical solutions to reduced-fidelity models can provide important information concerning the underlying physics of a problem and aid in guiding computational studies. Similarly, there is increased interest in sensitivity analysis studies. These studies commonly use computational tools. However, providing a complementary sensitivity study of relevant analytical models can lead to a deeper analysis of a problem. This work provides the analytical sensitivity analysis of the one-dimensional (1D) cylindrical mono-energetic neutron diffusion equation using the forward sensitivity analysis procedure (FSAP) developed by Cacuci. Further, these results are applied to a reduced-fidelity model of a spent nuclear fuel cask, demonstrating how computational analysis might be improved with a complementary analytic sensitivity analysis.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2019

Theoretical Analysis of the Entrainment–Mixing Process at Cloud Boundaries. Part II: Motion of Cloud Interface

The entrainment–mixing process leads to formation of a cloud–environment interface zone, consisting of the cloud dilution zone, where the liquid water content decreases as the distance from the cloud interior increases, and a humid shell with higher humidity in comparison with the dry environment. The time evolution of the location and width of these zones is analyzed by solving the diffusion–evaporation equation for the open region in the vicinity of the cloud–dry air interface. Upon normalization, the problem is reduced to a one-parametric one, the governing parameter being the potential evaporation parameter R < 0, which is proportional to the ratio of saturation deficit in the dry air to the available liquid water content in the cloud air. It is shown that the widths of the dilution zone and the humid shell increase with time. At R < −1, the interface between the dilution zone and the humid shell (i.e., the cloud edge) moves toward the cloud interior (i.e., the cloud dissipates). If R > −1, the cloud edge moves outward; that is, the mixing leads to an increase in the cloud volume. It is shown that the relative humidity remains high within most of the cloud dilution zone, and the effective radius changes only slightly, especially in growing clouds, where the effective radius changes by less than 10% within more than 80% of the cloud dilution zone. The strong decrease in the LWC is accompanied by a corresponding decrease in the droplet concentration. The study is illustrated by means of two examples of clouds: a small cumulus (Cu) typical of the trade wind zone and a deep convective cloud typical of the monsoon period over India. It is shown that while the small Cu tends to dissipate, the deep convective cloud expands. The characteristic sizes of the dilution zone and the humid shell are several tens of meters for the small Cu and several hundred meters for the deep convective cloud.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Polarization Tomography of Residual Stresses in Cylindrical Gradient-Index Lenses

In most cases, cylindrical structures with a gradient refractive index (GRIN) are formed on the basis of diffusion technology. The refractive index and the residual stresses in them are axially symmetrical. The distortion determined by diffusion and the thermal internal stresses in them are simulated by the fictive temperature of residual stresses. The use of the integral photoelasticity for the reconstruction of residual stresses in GRIN lenses is usually limited to the case of a planar deformed state. In this study, this algorithm is extended to the case of the axial change of residual stresses. It is assumed that the illumination procedure is carried out in the sample cross section. The optical problem is solved within the quasi-isotropic approximation. The reconstruction is based on measuring the characteristic parameters of polarized light passed through the sample and on the combined solution of the thermoelasticity problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Mesoscopic-microscopic spatial stochastic simulation with automatic system partitioning

The reaction-diffusion master equation (RDME) is a model that allows for efficient on-lattice simulation of spatially resolved stochastic chemical kinetics. Compared to off-lattice hard-sphere simulations with Brownian dynamics or Green’s function reaction dynamics, the RDME can be orders of magnitude faster if the lattice spacing can be chosen coarse enough. However, strongly diffusion-controlled reactions mandate a very fine mesh resolution for acceptable accuracy. It is common that reactions in the same model differ in their degree of diffusion control and therefore require different degrees of mesh resolution. This renders mesoscopic simulation inefficient for systems with multiscale properties. Mesoscopic-microscopic hybrid methods address this problem by resolving the most challenging reactions with a microscale, off-lattice simulation. However, all methods to date require manual partitioning of a system, effectively limiting their usefulness as “black-box” simulation codes. In this paper, we propose a hybrid simulation algorithm with automatic system partitioning based on indirect a priori error estimates. We demonstrate the accuracy and efficiency of the method on models of diffusion-controlled networks in 3D.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Physical and mechanical metallurgy of zirconium alloys for nuclear applications: a multi-scale computational study

In the post-Fukushima world, the stability of materials under extreme conditions is an important issue for the safety of nuclear reactors. Because the nuclear industry is going to continue using advanced zirconium cladding materials in the foreseeable future, it become critical to gain fundamental understanding of the several interconnected problems. First, what are the thermodynamic and kinetic factors affecting the oxidation and hydrogen pick-up by these materials at normal, off-normal conditions, and in long-term storage? Secondly, what protective coatings (if any) could be used in order to gain extremely valuable time at off-normal conditions, e.g., when temperature exceeds the critical value of 2200°F? Thirdly, the kinetics of oxidation of such protective coating or braiding needs to be quantified. Lastly, even if some degree of success is achieved along this path, it is absolutely critical to have automated inspection algorithms allowing identifying defects of cladding as soon as possible. This work strives to explore these interconnected factors from the most advanced computational perspective, utilizing such modern techniques as first-principles atomistic simulations, computational thermodynamics of materials, diffusion modeling, and the morphological algorithms of image processing for defect identification. Consequently, it consists of the four parts dealing with these four problem areas preceded by the introduction and formulation of the studied problems. In the 1st part an effort was made to employ computational thermodynamics and ab initio calculations to shed light upon the different stages of oxidation of ziraloys (2 and 4), the role of microstructure optimization in increasing their thermal stability, and the process of hydrogen pick-up, both in normal working conditions and in long-term storage. The 2nd part deals with the need to understand the influence and respective roles of the two different plasticity mechanisms in Zr nuclear alloys: twinning (at low T) and crystallographic slip (higher T’s). For that goal, a description of the advanced plasticity model is outlined featuring the non-associated flow rule in hcp materials including Zr. The 3rd part describes the kinetic theory of oxidation of the several materials considered to be perspective coating materials for Zr alloys: SiC and ZrSiO 4 . In the 4th part novel and advanced projectional algorithms for defect identification in zircaloy coatings are described. In so doing, the author capitalized on some 12 years of his applied industrial research in this area. Our conclusions and recommendations are presented in the 5th part of this work, along with the list of used literature and the scripts for atomistic, thermodynamic, kinetic, and morphological computations.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Acceleration of diffusive molecular dynamics simulations through mean field approximation and subcycling time integration

Highlights: • Two numerical methods are presented to accelerate DMD simulations. • The methods greatly reduce the computation time without much loss of accuracy. • The accelerated DMD model is capable of simulating long-term mass transport at atomistic length scales. • Demonstrated for the propagation of an atomically sharp phase boundary over 30 seconds. - Abstract: Diffusive Molecular Dynamics (DMD) is a class of recently developed computational models for the simulation of long-term diffusive mass transport at atomistic length scales. Compared to previous atomistic models, e.g., transition state theory based accelerated molecular dynamics, DMD allows the use of larger time-step sizes, but has a higher computational complexity at each time-step due to the need to solve a nonlinear optimization problem at every time-step. This paper presents two numerical methods to accelerate DMD based simulations. First, we show that when a many-body potential function, e.g., embedded atom method (EAM), is employed, the cost of DMD is dominated by the computation of the mean of the potential function and its derivatives, which are high-dimensional random variables. To reduce the cost, we explore both first- and second-order mean field approximations. Specifically, we show that the first-order approximation, which uses a point estimate to calculate the mean, can reduce the cost by two to three orders of magnitude, but may introduce relatively large error in the solution. We show that adding an approximate second-order correction term can significantly reduce the error without much increase in computational cost. Second, we show that DMD can be significantly accelerated through subcycling time integration, as the cost of integrating the empirical diffusion equation is much lower than that of the optimization solver. To assess the DMD model and the numerical approximation methods, we present two groups of numerical experiments that simulate the diffusion of hydrogen in palladium nanoparticles. In particular, we show that the computational framework is capable of capturing the propagation of an atomically sharp phase boundary over a time window of more than 30 seconds. The effects of the proposed numerical methods on solution accuracy and computation time are also assessed quantitatively.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Turbulence-induced relative velocity of dust particles. III. The probability distribution

Motivated by its important role in the collisional growth of dust particles in protoplanetary disks, we investigate the probability distribution function (PDF) of the relative velocity of inertial particles suspended in turbulent flows. Using the simulation from our previous work, we compute the relative velocity PDF as a function of the friction timescales, τ{sub p1} and τ{sub p2}, of two particles of arbitrary sizes. The friction time of the particles included in the simulation ranges from 0.1τ{sub η} to 54T {sub L}, where τ{sub η} and T {sub L} are the Kolmogorov time and the Lagrangian correlation time of the flow, respectively. The relative velocity PDF is generically non-Gaussian, exhibiting fat tails. For a fixed value of τ{sub p1}, the PDF shape is the fattest for equal-size particles (τ{sub p2} = τ{sub p1}), and becomes thinner at both τ{sub p2} < τ{sub p1} and τ{sub p2} > τ{sub p1}. Defining f as the friction time ratio of the smaller particle to the larger one, we find that, at a given f in (1/2) ≲ f ≲ 1, the PDF fatness first increases with the friction time τ{sub p,h} of the larger particle, peaks at τ{sub p,h} ≅ τ{sub η}, and then decreases as τ{sub p,h} increases further. For 0 ≤ f ≲ (1/4), the PDF becomes continuously thinner with increasing τ{sub p,h}. The PDF is nearly Gaussian only if τ{sub p,h} is sufficiently large (>>T {sub L}). These features are successfully explained by the Pan and Padoan model. Using our simulation data and some simplifying assumptions, we estimated the fractions of collisions resulting in sticking, bouncing, and fragmentation as a function of the dust size in protoplanetary disks, and argued that accounting for non-Gaussianity of the collision velocity may help further alleviate the bouncing barrier problem.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

An effective numerical algorithm for intra-granular fission gas release during non-equilibrium trapping and resolution

Both fission gas release and gaseous swelling in nuclear fuel are driven by the transport of fission gas from within the fuel grains to grain boundaries (intra-granular fission gas release). The process involves gas atom diffusion in conjunction with trapping in and resolution from intra-granular bubbles, and is described mathematically by a system of two partial differential equations (PDE). Following Speight, under the assumption of equilibrium between trapping and resolution (quasi-stationary approximation) the system can be reduced to a single diffusion equation with an effective diffusion coefficient. Numerical solutions used in engineering fuel performance calculations invariably rely on this simplification. First, we investigate the validity of the quasi-stationary approximation compared to the solution of the general system of PDEs. Results demonstrate that the approximation is valid under most conditions of practical interest, but is inadequate to model intra-granular fission gas release during rapid transients to relatively high temperatures such as postulated reactivity initiated accidents (RIA). Then, we develop a novel numerical algorithm for the solution of the general problem in time-varying conditions. We verify the PolyPole-2 algorithm against a reference finite difference solution for a large number of randomly generated operation histories including prototypical RIA transients. Results demonstrate that PolyPole-2 captures the solution of the general PDE system with a high accuracy and a low computational cost. As a result, the PolyPole-2 algorithm overcomes the quasi-stationary approximation and the concept of an effective diffusion coefficient for the solution of the intra-granular fission gas release problem in nuclear fuel analysis.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Automatically adaptive stabilized finite elements and continuation analysis for compaction banding in geomaterials

Under compressive creep, viscoplastic solids experiencing internal mass transfer processes can accommodate singular cnoidal wave solutions as material instabilities at the stationary wave limit. These instabilities appear when the loading rate is significantly faster than the material's capacity to diffusive internal perturbations, leading to localized failure features (e.g., cracks and compaction bands). These cnoidal waves, generally found in fluids, have strong nonlinearities that produce periodic patterns. Due to the singular nature of the solutions, the applicability of the theory is currently limited. Additionally, practical simulation tools require proper regularization to overcome the challenges that singularity induces. We focus on the numerical treatment of the governing equation using a nonlinear approach building on a recent adaptive stabilized finite element method. This automatic refinement method provides an error estimate that drives mesh adaptivity, a crucial feature for the problem at hand. We compare the performance of this adaptive strategy against analytical and standard finite element solutions. We then investigate the sensitivity of the diffusivity ratio, the parameter controlling the process, and identify multiple possible solutions with several stress peaks. We also show the evolution of the spacing between peaks for all solutions as a function of that parameter.

42 ENGINEERING↗
DOE OSTI2022

Verification of the AGREE and Serpent for the steady state HTR-10 benchmark problems

This paper provides results from steady state HTR-10 benchmark calculations performed using AGREE and Serpent that are compared to experimental results as well as calculations performed by INET. The purpose of completing this benchmark is to validate AGREE and Serpent for the prediction of HTGR operation so that they may ultimately be used to support licensing and deployment efforts of advanced reactors. The benchmark consists of several problems ranging from control rod worth calculations to k{sub eff} calculations at various temperatures. Overall, AGREE and Serpent show good agreement with the reference solutions and are effectively able to predict HTGR operation for a variety of steady state cases. The largest difference from the reference was for the initial core single rod worth, which is possibly due to the larger core helium cavity causing inaccuracies in the diffusion calculation whereas the control rod worth determined by Serpent and AGREE is much closer to the reference result for the full core loading. Moreover, in every benchmark problem, increasing the number of energy groups in the cross sections results in improved agreement of the AGREE result with the reference, although it also causes an increase in computation time. For problems with experimental data available, the accuracy of results generated by Serpent and AGREE is comparable to results obtained by other benchmark participants. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Connectivity‐based parcellation of normal and anatomically distorted human cerebral cortex

Abstract For over a century, neuroscientists have been working toward parcellating the human cortex into distinct neurobiological regions. Modern technologies offer many parcellation methods for healthy cortices acquired through magnetic resonance imaging. However, these methods are suboptimal for personalized neurosurgical application given that pathology and resection distort the cerebrum. We sought to overcome this problem by developing a novel connectivity‐based parcellation approach that can be applied at the single‐subject level. Utilizing normative diffusion data, we first developed a machine‐learning (ML) classifier to learn the typical structural connectivity patterns of healthy subjects. Specifically, the Glasser HCP atlas was utilized as a prior to calculate the streamline connectivity between each voxel and each parcel of the atlas. Using the resultant feature vector, we determined the parcel identity of each voxel in neurosurgical patients ( n = 40) and thereby iteratively adjusted the prior. This approach enabled us to create patient‐specific maps independent of brain shape and pathological distortion. The supervised ML classifier re‐parcellated an average of 2.65% of cortical voxels across a healthy dataset ( n = 178) and an average of 5.5% in neurosurgical patients. Our patient dataset consisted of subjects with supratentorial infiltrating gliomas operated on by the senior author who then assessed the validity and practical utility of the re‐parcellated diffusion data. We demonstrate a rapid and effective ML parcellation approach to parcellation of the human cortex during anatomical distortion. Our approach overcomes limitations of indiscriminately applying atlas‐based registration from healthy subjects by employing a voxel‐wise connectivity approach based on individual data.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

A new continuum approach for nonlinear kinetic simulation and transport analysis

A numerical code based on the I-transform approach is developed to solve the nonlinear Vlasov equation and carry out the transport analysis. The numerical results given by the I-transform approach agree with the conservative semi-Lagrangian approach in the Landau damping case and the bump-on-tail instability case. The diffusivities induced by the random fields and the quasilinear transport are also successfully demonstrated by using the new approach. It is found that the nonlinear transport in the one-dimensional Langmuir turbulence cannot be well-described by a simple diffusion model, due to the strong particle trapping at the nonlinear stage.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗