Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion equations”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Comparative experimental investigation on the actuation mechanisms of ionic polymer–metal composites with different backbones and water contents

Water-based ionic polymer–metal composites (IPMCs) exhibit complex deformation properties, especially when the water content changes. To explore the general actuation mechanisms, both Nafion and Flemion membranes are used as the polymer backbones. IPMC deformation includes three stages: fast anode deformation, relaxation deformation, and slow anode deformation, which is mainly dependent on the water content and the backbone. When the water content decreases from 21 to 14 wt. %, Nafion–IPMC exhibits a large negative relaxation deformation, zero deformation, a positive relaxation deformation, and a positive steady deformation without relaxation in sequence. Despite the slow anode deformation, Flemion–IPMC also shows a slight relaxation deformation, which disappears when the water content is less than 13 wt. %. The different water states are investigated at different water contents using nuclear magnetic resonance spectroscopy. The free water, which decreases rapidly at the beginning through evaporation, is proven to be critical for relaxation deformation. For the backbone, indirect evidence from the steady current response is correlated with the slow anode deformation of Flemion-IPMC. The latter is explained by the secondary dissociation of the weak acid group –COOH. Finally, we thoroughly explain not only the three deformations by swelling but also their evolvement with decreasing water content. A fitting model is also presented based on a multi-diffusion equation to reveal the deformation processes more clearly, the results from which are in good agreement with the experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Identification of the Diffusion Parameter in Nonlocal Steady Diffusion Problems

The problem of identifying the diffusion parameter appearing in a nonlocal steady diffusion equation is considered. The identification problem is formulated as an optimal control problem having a matching functional as the objective of the control and the parameter function as the control variable. The analysis makes use of a nonlocal vector calculus that allows one to define a variational formulation of the nonlocal problem. In a manner analogous to the local partial differential equations counterpart, we demonstrate, for certain kernel functions, the existence of at least one optimal solution in the space of admissible parameters. We introduce a Galerkin finite element discretization of the optimal control problem and derive a priori error estimates for the approximate state and control variables. Using one-dimensional numerical experiments, we illustrate the theoretical results and show that by using nonlocal models it is possible to estimate non-smooth and discontinuous diffusion parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Improved Modelling for the Assembly Discontinuity Factors in the APEC Method

Modern reactor analysis of light water reactor depends largely on a two-step procedure based on the simplified equivalence theory (SET). Consequently, the resulting homogenized two-group cross-section (XS) and associated discontinuity factor (DF) are subject to noticeable, intrinsic uncertainties. In our previous study, we introduced the albedo-corrected parameterized equivalence constants (APEC) method, which is a new leakage correction method applicable to the fuel assembly homogenized two-group XS of the standard two-step procedure. In the APEC method, fuel assembly two-group XS is expressed as a simple polynomial function of assembly-wise current-to-flux ratio (CFR), and it was demonstrated that two-group XSs can be well corrected during a nodal calculation and the nodal accuracy is improved a lot. However, since the conventional assembly discontinuity factor (ADF) is still used in the current APEC method, a small but non-negligible error is unavoidable. Therefore, it is necessary to include improved DF model in the APEC method for a higher accuracy. The functionalized discontinuity factor (FDF) using the definition of CFR was investigated by authors before the development of the APEC method, where it was motivated by the boundary perturbation theory by F. Rahnema. However, in our previous study on FDF, lattice calculations were performed only with one-dimensional two-group diffusion equation and could not be applied to actual core analysis. In this paper, feasibility study of DF functionalization is performed based on a 2-D method of characteristics (MOC) lattice code, DeCART2D.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2014

Comprehensive studies of response characteristics of organic photodetectors based on rubrene and C{sub 60}

We studied the transient response characteristics of organic photodetectors composing of high mobility materials of rubrene and C{sub 60}, respectively, as donor and acceptor. It was found that the response speed was limited by the delay time of both exciton diffusion and transit, and an anomalous phenomenon that the bandwidth decreases as the reverse bias increases was found for the first time. The investigation of frequency dependence at different device structures and light excitations demonstrated that the light absorption of rubrene prevents the photodetector from obtaining a fast response. With the help of magnetic field effect study, it was clearly elucidated that the slow diffusion time of the long lifetime triplet excitons generated from singlet fission in rubrene limited the bandwidth of the device. Moreover, the simulation of the response of photodetector under transient and steady state by exciton transport-diffusion equation showed that the exciton dissociation efficiency in rubrene increases more quickly than that in C{sub 60}, which should account for the negative dependence of bandwidth on the reverse bias in rubrene-based device.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Effects of spin diffusion on electron spin relaxation time measured with a time-resolved microscopic photoluminescence technique

We performed measurements at room temperature for a GaAs/AlGaAs multiple quantum well grown on GaAs(110) using a time-resolved microscopic photoluminescence (micro-PL) technique to find what effects spin diffusion had on the measured electron spin relaxation time, τ{sub s}, and developed a method of estimating the spin diffusion coefficient, D{sub s}, using the measured data and the coupled drift-diffusion equations for spin polarized electrons. The spatial nonuniformities of τ{sub s} and the initial degree of electron spin polarization caused by the pump intensity distribution inside the focal spot were taken into account to explain the dependence of τ{sub s} on the measured spot size, i.e., a longer τ{sub s} for a smaller spot size. We estimated D{sub s} as ∼100 cm{sup 2}/s, which is similar to a value reported in the literature. We also provided a qualitative understanding on how spin diffusion lengthens τ{sub s} in micro-PL measurements.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Charge sharing in pixelated semiconductor sensors

The charge sharing between neighboring pixels in pixelated sensors can be used to measure particle or x-ray coordinates with accuracy better than the pixel pitch. The accurate model of the charge distribution shape is essential to achieve ultimate coordinate accuracy. The charge sharing is caused by charge carriers' diffusion on the path from the generation point to pixels. This paper is focused on the diffusion of the initially compact charge cloud in the field free region. The diffusion equation solutions are obtained using separation of variable and Fourier synthesis method for different initial conditions and resulting charge distributions are integrated over pixel areas. In conclusion, the look up table containing pre-calculated values for pixel charge fractions is proposed to speed up numerical calculations.

47 OTHER INSTRUMENTATION↗
DOE OSTI2017

Growth of Silicon Carbide Nanolayers on Contact of Porous Carbon with Molten Silicon

The mathematical model of the growth of a silicon carbide layer includes two diffusion equations for silicon and carbon atoms and an equation describing the motion of the chemical conversion front. Analytical estimates and results of calculations are presented. Parametric calculations for different values of silicon solubility in silicon carbide are carried out. It is shown that the effective density of carbon exerts its influence on the duration of the process of silicon impregnation and the composite morphology.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2015

Compression induced phase transition of nematic brush: A mean-field theory study

Responsive behavior of polymer brush to the external compression is one of the most important characters for its application. For the flexible polymer brush, in the case of low grafting density, which is widely studied by the Gaussian chain model based theory, the compression leads to a uniform deformation of the chain. However, in the case of high grafting density, the brush becomes anisotropic and the nematic phase will be formed. The normal compression tends to destroy the nematic order, which leads to a complex responsive behaviors. Under weak compression, chains in the nematic brush are buckled, and the bending energy and Onsager interaction give rise to the elasticity. Under deep compression, the responsive behaviors of the nematic polymer brush depend on the chain rigidity. For the compressed rigid polymer brush, the chains incline to re-orientate randomly to maximize the orientational entropy and its nematic order is destroyed. For the compressed flexible polymer brush, the chains incline to fold back to keep the nematic order. A buckling-folding transition takes place during the compressing process. For the compressed semiflexible brush, the chains are collectively tilted to a certain direction, which leads to the breaking of the rotational symmetry in the lateral plane. These responsive behaviors of nematic brush relate to the properties of highly frustrated worm-like chain, which is hard to be studied by the traditional self-consistent field theory due to the difficulty to solve the modified diffusion equation. To overcome this difficulty, a single chain in mean-field theory incorporating Monte Carlo simulation and mean-field theory for the worm-like chain model is developed in present work. This method shows high performance for entire region of chain rigidity in the confined condition.

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

An information theoretic approach to use high-fidelity codes to calibrate low-fidelity codes

For many simulation models, it can be prohibitively expensive or physically infeasible to obtain a complete set of experimental data to calibrate model parameters. In such cases, one can alternatively employ validated higher-fidelity codes to generate simulated data, which can be used to calibrate the lower-fidelity code. In this paper, we employ an information-theoretic framework to determine the reduction in parameter uncertainty that is obtained by evaluating the high-fidelity code at a specific set of design conditions. These conditions are chosen sequentially, based on the amount of information that they contribute to the low-fidelity model parameters. The goal is to employ Bayesian experimental design techniques to minimize the number of high-fidelity code evaluations required to accurately calibrate the low-fidelity model. We illustrate the performance of this framework using heat and diffusion examples, a 1-D kinetic neutron diffusion equation, and a particle transport model, and include initial results from the integration of the high-fidelity thermal-hydraulics code Hydra-TH with a low-fidelity exponential model for the friction correlation factor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

The loss rates of O{sup +} in the inner magnetosphere caused by both magnetic field line curvature scattering and charge exchange reactions

With consideration of magnetic field line curvature (FLC) pitch angle scattering and charge exchange reactions, the O{sup +} (>300 keV) in the inner magnetosphere loss rates are investigated by using an eigenfunction analysis. The FLC scattering provides a mechanism for the ring current O{sup +} to enter the loss cone and influence the loss rates caused by charge exchange reactions. Assuming that the pitch angle change is small for each scattering event, the diffusion equation including a charge exchange term is constructed and solved; the eigenvalues of the equation are identified. The resultant loss rates of O{sup +} are approximately equal to the linear superposition of the loss rate without considering the charge exchange reactions and the loss rate associated with charge exchange reactions alone. The loss time is consistent with the observations from the early recovery phases of magnetic storms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Spin drift in highly doped n-type Si

A quantitative estimation of spin drift velocity in highly doped n-type silicon (Si) at 8 K is presented in this letter. A local two-terminal Hanle measurement enables the detection of a modulation of spin signals from the Si as a function of an external electric field, and this modulation is analyzed by using a spin drift-diffusion equation and an analytical solution of the Hanle-type spin precession. The analyses reveal that the spin drift velocity is linearly proportional to the electric field. The contribution of the spin drift effect to the spin signals is crosschecked by introducing a modified nonlocal four-terminal method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

The small-voxel tracking algorithm for simulating chemical reactions among diffusing molecules

Simulating the evolution of a chemically reacting system using the bimolecular propensity function, as is done by the stochastic simulation algorithm and its reaction-diffusion extension, entails making statistically inspired guesses as to where the reactant molecules are at any given time. Those guesses will be physically justified if the system is dilute and well-mixed in the reactant molecules. Otherwise, an accurate simulation will require the extra effort and expense of keeping track of the positions of the reactant molecules as the system evolves. One molecule-tracking algorithm that pays careful attention to the physics of molecular diffusion is the enhanced Green's function reaction dynamics (eGFRD) of Takahashi, Tănase-Nicola, and ten Wolde [Proc. Natl. Acad. Sci. U.S.A. 107, 2473 (2010)]. We introduce here a molecule-tracking algorithm that has the same theoretical underpinnings and strategic aims as eGFRD, but a different implementation procedure. Called the small-voxel tracking algorithm (SVTA), it combines the well known voxel-hopping method for simulating molecular diffusion with a novel procedure for rectifying the unphysical predictions of the diffusion equation on the small spatiotemporal scale of molecular collisions. Indications are that the SVTA might be more computationally efficient than eGFRD for the problematic class of non-dilute systems. A widely applicable, user-friendly software implementation of the SVTA has yet to be developed, but we exhibit some simple examples which show that the algorithm is computationally feasible and gives plausible results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Comment on “Towards high efficiency thin-film crystalline silicon solar cells: The roles of light trapping and non-radiative recombinations” [J. Appl. Phys. 115, 094501 (2014)]

In the above paper, an analytical approach including a new solution to the differential diffusion equation in illuminated quasi-neutral regions (QNR) is exploited to calculate the short-circuit current density (J{sub sc}), open-circuit voltage (V{sub oc}), fill factor (FF), and efficiency (η) of light-trapping (LT) c-Si solar cells with a given structure. Comparisons with numerical results calculated by the Silvaco ATLAS device simulator in the same LT cells show that the analytical results are systematically overestimated. According to the authors, the inaccuracies in J{sub sc}, V{sub oc}, and η are due to the fact that assuming ideal collection from space-charge region (SCR) and using the superposition approximation introduce systematic errors into analytical models. In this comment, an analytical approach using reported solutions to the transport equations in QNR and SCR, where ideal collection from SCR is assumed and the superposition approximation is used, is shown to agree with both the Silvaco and PC1d numerical approaches in calculating J{sub sc}, V{sub oc}, and η, in the same LT devices as considered in the commented paper. Reasons for the inaccuracies detected in the commented paper are suggested.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

POLCA8 - modelling of cross section variations inside hexagonal assemblies

This paper presents the POLCA8 approach for modelling non-constant cross section distributions inside hexagonal fuel assemblies. The multigroup diffusion equation is modified to account for intranodal cross section variations. The obtained equation is solved in a node-wise manner based on the Fourier expansion method. As a result of varying cross sections, the solution includes a particular part additionally to the homogeneous one. A method for obtaining the particular solution is derived. Numerical tests on a VVER-1000 core are presented showing the impact of cross-section variations to some key parameters for reactor operation. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2014

Analysis of Hanle-effect signals observed in Si-channel spin accumulation devices

We reexamined curve-fitting analysis for spin-accumulation signals observed in Si-channel spin-accumulation devices, employing widely-used Lorentz functions and a new formula developed from the spin diffusion equation. A Si-channel spin-accumulation device with a high quality ferromagnetic spin injector was fabricated, and its observed spin-accumulation signals were verified. Experimentally obtained Hanle-effect signals for spin accumulation were not able to be fitted by a single Lorentz function and were reproduced by the newly developed formula. Our developed formula can represent spin-accumulation signals and thus analyze Hanle-effect signals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Simulation of stochastic diffusion via first exit times

In molecular biology it is of interest to simulate diffusion stochastically. In the mesoscopic model we partition a biological cell into unstructured subvolumes. In each subvolume the number of molecules is recorded at each time step and molecules can jump between neighboring subvolumes to model diffusion. The jump rates can be computed by discretizing the diffusion equation on that unstructured mesh. If the mesh is of poor quality, due to a complicated cell geometry, standard discretization methods can generate negative jump coefficients, which no longer allows the interpretation as the probability to jump between the subvolumes. We propose a method based on the mean first exit time of a molecule from a subvolume, which guarantees positive jump coefficients. Two approaches to exit times, a global and a local one, are presented and tested in simulations on meshes of different quality in two and three dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Migration Rates of Some Cations in Unsaturated Layer around the Near-Surface Disposal Facilities of Radioactive Wastes - 16116

In the porous layer around the near-surface underground, the intermittent recharge process of rainwater forms a layer perfectly not saturated with liquid phase. Here, such a layer is simply called 'unsaturated layer'. In this study, the retardation effects of radionuclide migration were examined under an unsaturated condition by using a column packed with silica sand particles. As tracer ions, sodium (Na), europium (Eu) and cesium (Cs) were selected for estimating the retardation effects under both the saturated and unsaturated condition. Particularly, Cs is one of the key elements included in the contaminated soil due to Fukushima Daiichi Nuclear Power Plant Accident. These tracer ions were injected into column with the condition of constant flow rate or constant pressure flow. To quantify the influence of the unsaturated condition on the migration of radionuclides, retardation coefficient (R{sub d}) and Peclet number (P{sub e}) were estimated by applying advection diffusion equation to the experimental results. As a result, the conflicting effects on the retardation effect in the unsaturated layer were suggested; one is the decrease in the retardation effect by interrupting the sorption of ions with gas phase, the other is the increase in the retardation effect by the clogging with gas phase in flow paths. For the assessment of the radionuclide migration in the unsaturated layer, it is important to clarify such complicated behaviors with considering key parameters related to the migration of ions. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗