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 1,009 records · Page 56

Finite element procedures for time-dependent convection-diffusion-reaction systems

New finite element procedures based on the streamline-upwind/Petrov-Galerkin formulations are developed for time-dependent convection-diffusion-reaction equations. These procedures minimize spurious oscillations for convection-dominated and reaction-dominated problems. The results obtained for representative numerical examples are accurate with minimal oscillations. As a special application problem, the single-well chemical tracer test (a procedure for measuring oil remaining in a depleted field) is simulated numerically. The results show the importance of temperature effects on the interpreted value of residual oil saturation from such tests.

Tezduyar, T. E.↗

Recent developments in high order K-exact reconstruction on unstructured meshes

This paper presents recent improvements in high-order K-exact reconstruction on unstructured meshes. The new reconstruction procedures are incorporated into a basic upwind finite-volume scheme suitable for solving scalar advection-diffusion equations as well as the Euler and Navier-Stokes equations. Numerical calculations are performed comparing the present method with lower order accurate reconstruction procedures (piecewise constant and piecewise linear) and various competing technologies such as the fluctuation splitting method of Roe (1987) and Deconinck et al. (1992) and a system-variant of the streamline diffusion Petrov-Galerkin method developed by Hansbo (1991) and Hansbo and Johnson (1991). Five test problems are used in the numerical comparisons: scalar circular advection, transonic and supersonic Euler flow, laminar boundary-layer flow, and general compressible Navier-Stokes flow.

Barth, Timothy J.↗

Morphological instability in epitaxially strained dislocation-free solid films - Linear stability theory

The morphological instability of a growing epitaxially strained dislocation-free solid film is analyzed. An evolution equation for the film surface is derived in the dilute limit of vacancies based on surface diffusion driven by a stress-dependent chemical potential. From the time-dependent linear stability problem the conditions for which a growing film is unstable are determined. It is found that the instability is driven by the lattice mismatch between the film and the substrate; however, low temperatures as well as elastically stiff substrates are stabilizing influences. The results also reveal that the critical film thickness for instability depends on the growth rate of the film itself. Detailed comparison with experimental observations indicates that the instability described exhibits many of the observed features of the onset of the 'island instability'.

Spencer, B. J.↗

Micromechanical thermal analysis of interphase region in a titanium aluminide MMC

The high reactivity between the fiber and matrix in silicon carbide/titanium aluminide MMCs leads to the formation of brittle reaction products at the fiber/matrix (F/M) interface. Also, the high thermal expansion coefficient mismatch between the fiber and matrix leads to high tensile residual stresses at the F/M interface, and this can lead to premature cracking during cooldown. One solution to these problems is the use of a metallic fiber coating like Ta which acts as an F/M diffusion barrier and reacts with the matrix to form a beta stabilized compliant layer. A finite element micromechanics analysis was performed to study the effects of Ta and beta interphase layers on the thermal residual stresses during consolidation. A 5-micron-thick beta layer reduced cool-down stresses by 8 percent compared to a 2 percent reduction computed for a Ta layer of the same thickness. Plastic yielding in the Ta was not effective in reducing cool-down stresses. Compliant alpha-2 particles next to the stiffer gamma particles reduced stresses in the gamma particles by less than 2 percent. A simple closed form analysis was developed to calculate thermal residual stresses in a fiber/interphase/matrix system. A 2-micron-thick Ag interphase layer was found to reduce residual stresses by about 11 percent.

Naik, R. A.↗

Multigrid Algorithms with Projection and Prolongation over Elements of the Phase Space for K-Eigenvalue Transport Problems

This paper describes new multilevel acceleration methods for solving the multigroup neu- tron transport eigenvalue problems. These multilevel algorithms use different projection and prolongation operators in the phase space. The Nonlinear Diffusion Acceleration (NDA) method with multiple grids in energy is formulated with the prolongation oper- ator based on multiplication iterative correction and linear-in-energy mapping. Another multilevel NDA method uses the projection operator with coarsening in energy between the high-order transport and low-order NDA equations. The third algorithm is formu- lated with the partial-current based CMFD low-order equations and applies projection operators in space and energy. The numerical results are presented.

Cornejo, Luke↗

Transient MOC with frequency transform and DSA on unstructured mesh

We present an implementation of the transient method of characteristics (MOC) with isotropic time derivatives, accelerated by diffusion synthetic acceleration (DSA). The fully implicit frequency transform method is used to solve the transient problem with analytic precursor integration. The code works on meshes composed of almost any of the commonly used non-curvilinear finite element types, and can handle the deformation of geometry in time-dependent transport calculations. We present results of a continuous Fourier analysis for the transient multigroup DSA problem, and representative benchmarking results are presented for the C5G7-TD benchmark in 2D showing reasonable performance and agreement compared to other codes. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Improvement of the second- and third-moment modeling of turbulence: A study of Reynolds-stress closure model

Four parts of the Reynolds-stress closure modeling are reported: (1) improvement of the k and epsilon equaitons; (2) development of the third-moment transport equation; (3) formulation of the diffusion coefficient of the momentum equation by using the algebraic-stress model of turbulence; and (4) the application of the Reynolds-stress model to a heat exchanger problem. It was demonstrated that the third-moment transport model improved the prediction of the triple-velocity products in the recirculating and reattaching flow regions in comparison with the existing algebraic models for the triple-velocity products. Optimum values for empirical coefficients are obtained for the prediction of the backward-facing step flows. A functional expression is derived for the coefficient of the momentum diffusion by employing the algebraic-stress model. The second-moment closure is applied to a heat transfer problem. The computations for the flow in a corrugated-wall channel show that the second-moment closure improves the prediction of the heat transfer rates by 30% over the k - epsilon model.

Amano, R. S.↗

The frequency-dependent response of single aerosol particles to vapour phase oscillations and its application in measuring diffusion coefficients

A new method for measuring diffusion in the condensed phase of single aerosol particles is proposed and demonstrated. The technique is based on the frequency-dependent response of a binary particle to oscillations in the vapour phase of one of its chemical components. Here, we discuss how this physical situation allows for what would typically be a non-linear boundary value problem to be approximately reduced to a linear boundary value problem. For the case of aqueous aerosol particles, we investigate the accuracy of the closed-form analytical solution to this linear problem through a comparison with the numerical solution of the full problem. Then, using experimentally measured whispering gallery modes to track the frequency-dependent response of aqueous particles to relative humidity oscillations, we determine diffusion coefficients as a function of water activity. The measured diffusion coefficients are compared to previously reported values found using the two common experiments: (i) the analysis of the sorption/desorption of water from a particle after a step-wise change to the surrounding relative humidity and (ii) the isotopic exchange of water between a particle and the vapour phase. The technique presented here has two main strengths: first, when compared to the sorption/desorption experiment, it does not require the numerical evaluation of a boundary value problem during the fitting process as a closed-form expression is available. Second, when compared to the isotope exchange experiment, it does not require the use of labeled molecules. Therefore, the frequency-dependent experiment retains the advantages of these two commonly used methods but does not suffer from their drawbacks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Numerical Treatment of Degenerate Diffusion Equations via Feller's Boundary Classification, and Applications

A numerical method is devised to solve a class of linear boundary-value problems for one-dimensional parabolic equations degenerate at the boundaries. Feller theory, which classifies the nature of the boundary points, is used to decide whether boundary conditions are needed to ensure uniqueness, and, if so, which ones they are. The algorithm is based on a suitable preconditioned implicit finite-difference scheme, grid, and treatment of the boundary data. Second-order accuracy, unconditional stability, and unconditional convergence of solutions of the finite-difference scheme to a constant as the time-step index tends to infinity are further properties of the method. Several examples, pertaining to financial mathematics, physics, and genetics, are presented for the purpose of illustration.

Cacio, Emanuela↗

A high temperature thermal diffusivity determination procedure for solids and liquids

A new method for measuring the thermal diffusivity of materials at high temperatures is presented. The method is applicable to solids on Earth, and to liquids in the reduced gravity environment of space. It is especially suited to levitated liquid metals at elevated temperatures where thermal diffusivity data is not available. The method is applied in two parts, such that lumped analysis is valid in the first part, and Fourier's law of conduction is valid in the second. In both parts, the spherical specimen is assumed to have been heated to a desired temperature and cooled. An inverse conduction problem is then formulated and solved using Laplace transformation techniques. Using this solution, sample sizes, and experimentally obtained surface temperature history, the thermal diffusivity is determined by minimizing a function that satisifies the heat balance at the surface. Minimization is performed using a modified quasilinearization algorithm. Accuracy is very sensitive to error in the temperature data and increases with better curve fits to the temperature data. An error analysis is also performed, and the effect of errors in the various parameters on the evaluated thermal diffusivity is determined. An experimental study for solids on Earth is suggested, before development for implementation in space.

Bayazitoglu, Yildiz↗

Numerical vorticity capturing for vortex-solid body interaction problems

A new numerical method based on adding a term to the Euler/Navier-Stokes equations has been demonstrated it can effectively treat vortex-dominated flows using low-order numerical schemes and coarse grids. The numerical diffusion introduced by convection schemes can be eliminated by modifying the velocity field as a result of solving the modified flow governing equations. The modification of the velocity field, which conserves the total vorticity, essentially convects the vorticity toward its local extreme to offset the numerical diffusion. This method preserves the vortex structure even when vortices travel on coarse grid region. The method is implemented with two distinctive codes: a Navier-Stokes flow solver based on vorticity-velocity formulation and an Euler code based on velocity-pressure formulation. Problems presented in this paper include airfoil dynamic stall, vortex airfoil interaction, and vortex-fuselage interactions.

Wang, Clin M.↗

Verification of RadCAD: Specular Capabilities

As part of the RadCAD's development process, it is necessary to compare RadCAD's results with other radiation tools and exact solutions when and where possible. Form factor algorithms have been previously verified with exact solutions. This paper will consider RadCAD's specular capabilities. First, radiation exchange factors will be compared against exact solutions and results from TRASYS for various geometries. Critical dimensions and optical properties are changed for each geometry. Second, a specular adjunct plate system will be used to verify absorbed heat fluxes. This particular geometric problem has had some attention in the literature. Previous authors have used this problem to validate software results with exact analytical solution. This paper will compare absorbed heat rates against the exact solution and other published results from other thermal radiation tools. The agreement between RadCAD and the exact solutions is good. The maximum error for both specular and diffuse exchange factors for both geometries and all optical properties was 3%. The absorbed fluxes differed by a maximum of 4% for the adjunct plate problem.

Lucas, Sam↗

Hydrogen Effects on GaAs, Status and Progress

It has been observed that GaAs devices in hermetically sealed packages when exposed to temperatures as low as 125oC can exhibit unacceptable degradation in both RF and DC characteristics. The source of the problem was found to be hydrogen that has been absorbed in the package's metals (Kovar, plating, etc.) and converted into atomic hydrogen within the Pt metallization of the gate structure. Subsequently, atomic hydrogen diffuses into the channel region of the FET structure and neutralizes the Si donors, resulting in a degradation of the device characteristics. Most of the data indicates the onslaught of the problem to occur after 500 hours at 125oC. This onslaught has been observed to be dependent on the thickness and processing conditions of the passivation layer, the sealing environment, and the amount of Pt or Pd in the gate structure of the device.

Hydrogen GaAs↗

Implementation of the full viscoresistive magnetohydrodynamic equations in a nonlinear finite element code

Numerical simulations form an indispensable tool to understand the behavior of a hot plasma that is created inside a tokamak for providing nuclear fusion energy. Various aspects of tokamak plasmas have been successfully studied through the reduced magnetohydrodynamic (MHD) model. The need for more complete modeling through the full MHD equations is addressed here. Our computational method is presented along with measures against possible problems regarding pollution, stability, and regularity. The problem of ensuring continuity of solutions in the center of a polar grid is addressed in the context of a finite element discretization of the full MHD equations. A rigorous and generally applicable solution is proposed here. Useful analytical test cases are devised to verify the correct implementation of the momentum and induction equation, the hyperdiffusive terms, and the accuracy with which highly anisotropic diffusion can be simulated. A striking observation is that highly anisotropic diffusion can be treated with the same order of accuracy as isotropic diffusion, even on non-aligned grids, as long as these grids are generated with sufficient care. This property is shown to be associated with our use of a magnetic vector potential to describe the magnetic field. Several well-known instabilities are simulated to demonstrate the capabilities of the new method. The linear growth rate of an internal kink mode and a tearing mode are benchmarked against the results of a linear MHD code. The evolution of a tearing mode and the resulting magnetic islands are simulated well into the nonlinear regime. The results are compared with predictions from the reduced MHD model. Finally, a simulation of a ballooning mode illustrates the possibility to use our method as an ideal MHD method without the need to add any physical dissipation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Origin and Time Dependence of the Amount and Composition of Non-Constituent Gases Present in Crystal Growth Systems

Presence of different, non-constituent gases may be a critical factor in crystal growth systems. In Physical Vapor Transport processes the cras(es) can be used intentionally (to prevent excessively high, unstable growth conditions), or can evolve unintentionally during the course of the process (which may lead to undesired reduction in the -rowth rate). In melt growth, particularly under low gravity conditions (reduced hydrostatic pressure) the gas present in the system may contribute to formation of voids in the growing crystals and even to a separation of the crystal and the liquid phase [1]. On the other hand, some amount of gas may facilitate 'contactless' crystal growth particularly under reduced gravity conditions [2 - 6]. Different non-constituent gases may be present in growth ampoules, and their amount and composition may change during the crystallization process. Some gases can appear even in empty ampoules sealed originally under high vacuum: they may diffuse in from the outside, and/or desorb from the ampoule walls. Residual gases can also be generated by the source materials: even very high purity commercial elements and compounds may contain trace amounts of impurities, particularly oxides. The oxides may have low volatilities themselves but their reaction with other species, particularly carbon and hydrogen, may produce volatile compounds like water or carbon oxides. The non-constituent gases, either added initially to the system or evolved during the material processing, may diffuse out of the ampoule during the course of the experiment. Gases present outside (e.g. as a protective atmosphere or thermal conductor) may diffuse into the ampoule. In either case the growth conditions and the quality of the crystals may be affected. The problem is of a particular importance in sealed systems where the amount of the gases cannot be directly controlled. Therefore a reasonable knowledge and understanding of the origin, composition, magnitude, and change with time of gases present in sealed ampoules may be important for a meaningful control and interpretation of crystal growth processes. This problem is of a particular importance for processing of electronic materials in space because (i) safety considerations require using sealed systems only, and (ii) high cost of crystal growth experiments in microgravity calls for a throughout, accurate description of the processing conditions necessary for a meaningful, efficient, and conclusive interpretation of the space results. In this paper we present the results of our extensive studies on gases in closed crystal growth systems which include: (a) Degassing properties of fused silica; (b) Generation of inert gases by source materials (CdTe, ZnTe, CdZnTe, ZnSe, PbTe, PbSe, PbSeTe); (c) Diffusive cas losses from silica glass ampoules.

Palosz, Witold↗

COSMIC RAY TRANSPORT WITH MAGNETIC FOCUSING AND THE “TELEGRAPH” MODEL

Cosmic rays (CR), constrained by scattering on magnetic irregularities, are believed to propagate diffusively. However, a well-known defect of diffusive approximation, whereby some of the particles propagate unrealistically fast, has directed interest toward an alternative CR transport model based on the “telegraph” equation. Though, its derivations often lack rigor and transparency leading to inconsistent results. We apply the classic Chapman–Enskog method to the CR transport problem. We show that no “telegraph” (second order time derivative) term emerges in any order of a proper asymptotic expansion with systematically eliminated short timescales. Nevertheless, this term may formally be converted from the fourth order hyper-diffusive term of the expansion. However, both the telegraph and hyperdiffusive terms may only be important for a short relaxation period associated with either strong pitch-angle anisotropy or spatial inhomogeneity of the initial CR distribution. Beyond this period the system evolves diffusively in both cases. The term conversion, that makes the telegraph and Chapman–Enskog approaches reasonably equivalent, is possible only after this relaxation period. During this period, the telegraph solution is argued to be unphysical. Unlike the hyperdiffusion correction, it is not uniformly valid and introduces implausible singular components to the solution. These dominate the solution during the relaxation period. Because they are shown not to be inherent in the underlying scattering problem, we argue that the telegraph term is involuntarily acquired in an asymptotic reduction of the problem.

79 ASTRONOMY AND ASTROPHYSICS↗

Partitioned coupling strategies for multi-physically coupled radiative heat transfer problems

This article aims to propose new aspects concerning a partitioned solution strategy for multi-physically coupled fields including the physics of thermal radiation. Particularly, we focus on the partitioned treatment of electro–thermo-mechanical problems with an additional fourth thermal radiation field. One of the main goals is to take advantage of the flexibility of the partitioned approach to enable combinations of different simulation software and solvers. Within the frame of this article, we limit ourselves to the case of nonlinear thermoelasticity at finite strains, using temperature-dependent material parameters. For the thermal radiation field, diffuse radiating surfaces and gray participating media are assumed. Moreover, we present a robust and fast partitioned coupling strategy for the fourth field problem. Stability and efficiency of the implicit coupling algorithm are improved drawing on several methods to stabilize and to accelerate the convergence. To conclude and to review the effectiveness and the advantages of the additional thermal radiation field several numerical examples are considered to study the proposed algorithm. In particular we focus on an industrial application, namely the electro–thermo-mechanical modeling of the field-assisted sintering technology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Results of Oklo Diffusion Couple Experiments

Through the Gateway for Accelerated Innovation in Nuclear (GAIN) initiative, the U.S. Department of Energy funded a GAIN Nuclear Energy Voucher to perform out-of-pile experiments related to the fuel design for Oklo Inc., headquartered in Santa Clara, CA. GAIN’s goal is to accelerate the innovation and application of advanced nuclear technologies, and voucher recipients do not receive direct financial awards but are provided access to national laboratory capabilities at no cost. This report documents the results of diffusion couple experiments performed as part of the GAIN-funded project for Oklo. This work addresses the specific problem of filling the gap in understanding of the fuel-cladding chemical interaction (FCCI) that result from the past focus on a particular fuel design. These knowledge and mechanistic modeling gaps make it difficult to extrapolate from prior fuel performance data to advance an economically competitive fuel design. Specifically, the current state-of-the-art fuel performance modeling cannot separate the entangled effects of temperature ranges, burnups, materials, and geometries.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗