Search NASA⌕ Search

SEARCH · Search NASA

Results for “Solute diffusion”

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 73 records · Page 4

A Systematic Solution Approach for Neutron Transport Problems in Diffuse Regimes

A systematic solution approach for the neutron transport equation, based on a least-squares finite-element discretization, is presented. This approach includes the theory for the existence and uniqueness of the analytical as well as of the discrete solution, bounds for the discretization error, and guidance for the development of an efficient multigrid solver for the resulting discrete problem. To guarantee the accuracy of the discrete solution for diffusive regimes, a scaling transformation is applied to the transport operator prior to the discretization. The key result is the proof of the V-ellipticity and continuity of the scaled least-squares bilinear form with constants that are independent of the total cross section and the absorption cross section. For a variety of least-squares finite-element discretizations this leads to error bounds that remain valid in diffusive regimes. Moreover, for problems in slab geometry a full multigrid solver is presented with V(1, 1)-cycle convergence rates approximately equal to 0.1, independent of the size of the total cross section and the absorption cross section.

Manteuffel, T. A.↗

Wind-US Unstructured Flow Solutions for a Transonic Diffuser

The Wind-US Computational Fluid Dynamics flow solver computed flow solutions for a transonic diffusing duct. The calculations used an unstructured (hexahedral) grid. The Spalart-Allmaras turbulence model was used. Static pressures along the upper and lower wall agreed well with experiment, as did velocity profiles. The effect of the smoothing input parameters on convergence and solution accuracy was investigated. The meaning and proper use of these parameters are discussed for the benefit of Wind-US users. Finally, the unstructured solver is compared to the structured solver in terms of run times and solution accuracy.

Mohler, Stanley R., Jr.↗

Velocity and concentration distribution in a Stefan diffusion tube

Numerical solutions for the diffusion- and gravity-driven flow in a cylindrical Stefan tube were obtained from the coupled diffusion and Navier-Stokes equations of Peclet numbers 0.3, 1 and 5. Distributions of binary component concentrations and velocities were calculated. The mass average velocity is parabolic in nature, except at high Peclet numbers. The solvent is not stagnant but recirculates, even in the absence of gravity. Radial concentration gradients develop which act convectively destabilizing. Consequences for the deduction of diffusion coefficients from Stefan tube experiments are discussed.

Markham, B. L.↗

Analytic solution of the box diffusion model for a global ocean

The paper presents an analytic solution for the temperature of the mixed layer for a box diffusion model of a global ocean. It is shown that the exact solution approaches the appropriate limits of a mixed layer ocean model and of a diffusive ocean without a mixed layer on top. It is also shown that, over a large range of parameters of the box diffusion model, the difference between the temperature calculated using a box diffusion model and the temperature calculated using a purely diffusive ocean is buried in the noise.

Lebedeff, Sergej A.↗

Growth and characterization of Hg(1-x)Cd(x)Se alloys

A detailed evaluation of the influence of growth conditions on the radial and axial compositional variations in directionally solidified Hg(1-x)Cd(x)Se alloys was performed. The measured axial compositional profiles were fitted to theoretical profiles to determine the effective solute (CdSe) diffusion coefficient (D) for the Hg(1-x)Cd(x)Se system. The value for D was estimated to be 3.0 + or - 0.5 x 10 to the -5th sq cm/s and did not appear to be significantly affected by the rate of crystal growth. The axial crystal uniformity was shown to be highly growth-rate dependent, with the faster growth rates producing crystals of more uniform composition in the axial direction. The magnitude of the radial variations was also shown to be highly growth-rate dependent, with the slower growth rates producing crystals of greater radial uniformity. This translation rate dependence of the radial uniformity is discussed in terms of lateral solute diffusion and convective interfacial fluid flows.

Andrews, R. N.↗

Comparison of the Radiative Two-Flux and Diffusion Approximations

Approximate solutions are sometimes used to determine the heat transfer and temperatures in a semitransparent material in which conduction and thermal radiation are acting. A comparison of the Milne-Eddington two-flux approximation and the diffusion approximation for combined conduction and radiation heat transfer in a ceramic material was preformed to determine the accuracy of the diffusion solution. A plane gray semitransparent layer without a substrate and a non-gray semitransparent plane layer on an opaque substrate were considered. For the plane gray layer the material is semitransparent for all wavelengths and the scattering and absorption coefficients do not vary with wavelength. For the non-gray plane layer the material is semitransparent with constant absorption and scattering coefficients up to a specified wavelength. At higher wavelengths the non-gray plane layer is assumed to be opaque. The layers are heated on one side and cooled on the other by diffuse radiation and convection. The scattering and absorption coefficients were varied. The error in the diffusion approximation compared to the Milne-Eddington two flux approximation was obtained as a function of scattering coefficient and absorption coefficient. The percent difference in interface temperatures and heat flux through the layer obtained using the Milne-Eddington two-flux and diffusion approximations are presented as a function of scattering coefficient and absorption coefficient. The largest errors occur for high scattering and low absorption except for the back surface temperature of the plane gray layer where the error is also larger at low scattering and low absorption. It is shown that the accuracy of the diffusion approximation can be improved for some scattering and absorption conditions if a reflectance obtained from a Kubelka-Munk type two flux theory is used instead of a reflection obtained from the Fresnel equation. The Kubelka-Munk reflectance accounts for surface reflection and radiation scattered back by internal scattering sites while the Fresnel reflection only accounts for surface reflections.

Spuckler, Charles M.↗

Approximate Solutions Of Equations Of Steady Diffusion

Rigorous analysis yields reliable criteria for "best-fit" functions. Improved "curve-fitting" method yields approximate solutions to differential equations of steady-state diffusion. Method applies to problems in which rates of diffusion depend linearly or nonlinearly on concentrations of diffusants, approximate solutions analytic or numerical, and boundary conditions of Dirichlet type, of Neumann type, or mixture of both types. Applied to equations for diffusion of charge carriers in semiconductors in which mobilities and lifetimes of charge carriers depend on concentrations.

Edmonds, Larry D.↗

Normal evaporation of binary alloys

In the study of normal evaporation, it is assumed that the evaporating alloy is homogeneous, that the vapor is instantly removed, and that the alloy follows Raoult's law. The differential equation of normal evaporation relating the evaporating time to the final solute concentration is given and solved for several important special cases. Uses of the derived equations are exemplified with a Ni-Al alloy and some binary iron alloys. The accuracy of the predicted results are checked by analyses of actual experimental data on Fe-Ni and Ni-Cr alloys evaporated at 1600 C, and also on the vacuum purification of beryllium. These analyses suggest that the normal evaporation equations presented here give satisfactory results that are accurate to within an order of magnitude of the correct values, even for some highly concentrated solutions. Limited diffusion and the resultant surface solute depletion or enrichment appear important in the extension of this normal evaporation approach.

Li, C. H.↗

Effect of Water on the Thermo-Mechanical Behavior of Carbon Cloth Phenolic

The results of thermo-mechanical experiments, which were conducted previously by one of the authors, are reviewed. The strain in the direction normal to the fabric plane was measured as a function of temperature for a variety of initial moisture contents and heating rates. In this paper, the general features of the thermo-mechanical response are discussed and the effect of heating rate and initial moisture content are highlighted. The mechanical interaction between the phenolic polymer and water trapped within its free volumes as the polymer is heated to high temperatures is discussed. An equation for the internal stresses which are generated within the polymer due to trapped water is obtained from the total stress expression for a binary mixture of polymer and water. Numerical solutions for moisture diffusion in the thermo-mechanical experiments were performed and the results of these solutions are presented. The results of the moisture diffusion solutions help to explain the effects of heating rate and moisture content on the strain behavior normal to the fabric plane.

Sullivan, Roy M.↗

Unsteady Crystal Growth Due to Step-Bunch Cascading

Based on our experimental findings of growth rate fluctuations during the crystallization of the protein lysozym, we have developed a numerical model that combines diffusion in the bulk of a solution with diffusive transport to microscopic growth steps that propagate on a finite crystal facet. Nonlinearities in layer growth kinetics arising from step interaction by bulk and surface diffusion, and from step generation by surface nucleation, are taken into account. On evaluation of the model with properties characteristic for the solute transport, and the generation and propagation of steps in the lysozyme system, growth rate fluctuations of the same magnitude and characteristic time, as in the experiments, are obtained. The fluctuation time scale is large compared to that of step generation. Variations of the governing parameters of the model reveal that both the nonlinearity in step kinetics and mixed transport-kinetics control of the crystallization process are necessary conditions for the fluctuations. On a microscopic scale, the fluctuations are associated with a morphological instability of the vicinal face, in which a step bunch triggers a cascade of new step bunches through the microscopic interfacial supersaturation distribution.

Vekilov, Peter G.↗

Solutions for Reacting and Nonreacting Viscous Shock Layers with Multicomponent Diffusion and Mass Injection

Numerical solutions are presented for the viscous shocklayer equations where the chemistry is treated as being either frozen, equilibrium, or nonequilibrium. Also the effects of the diffusion model, surface catalyticity, and mass injection on surface transport and flow parameters are considered. The equilibrium calculations for air species using multicomponent: diffusion provide solutions previously unavailable. The viscous shock-layer equations are solved by using an implicit finite-difference scheme. The flow is treated as a mixture of inert and thermally perfect species. Also the flow is assumed to be in vibrational equilibrium. All calculations are for a 45 deg hyperboloid. The flight conditions are those for various altitudes and velocities in the earth's atmosphere. Data are presented showing the effects of the chemical models; diffusion models; surface catalyticity; and mass injection of air, water, and ablation products on heat transfer; skin friction; shock stand-off distance; wall pressure distribution; and tangential velocity, temperature, and species profiles.

Moss, J. N.↗

An integral solution to a nonlinear diffusion problem

The generalized Galerkin method (or the method of integral relations) is applied to the type of problem described by quasilinear parabolic equations. As an example the problem of nonlinear transient slab diffusion with a general reservoir boundary condition is worked out. The integral relations are given for an arbitrary number of strips, and solutions using up to seven strips have been obtained in order to investigate the convergence of the method.

Hussaini, M. Y.↗

Solute redistribution in dendritic solidification with diffusion in the solid

An investigation of solute redistribution during dendritic solidification with diffusion in the solid has been performed using numerical techniques. The extent of diffusion is characterized by the instantaneous and average diffusion parameters. These parameters are functions of the diffusion Fourier number, the partition ratio and the fraction solid. Numerical results are presented as an approximate model, which is used to predict the average diffusion parameter and calculate the composition of the interdendritic liquid during solidification.

Ganesan, S.↗

A Reconstruction Approach to High-Order Schemes Including Discontinuous Galerkin for Diffusion

We introduce a new approach to high-order accuracy for the numerical solution of diffusion problems by solving the equations in differential form using a reconstruction technique. The approach has the advantages of simplicity and economy. It results in several new high-order methods including a simplified version of discontinuous Galerkin (DG). It also leads to new definitions of common value and common gradient quantities at each interface shared by the two adjacent cells. In addition, the new approach clarifies the relations among the various choices of new and existing common quantities. Fourier stability and accuracy analyses are carried out for the resulting schemes. Extensions to the case of quadrilateral meshes are obtained via tensor products. For the two-point boundary value problem (steady state), it is shown that these schemes, which include most popular DG methods, yield exact common interface quantities as well as exact cell average solutions for nearly all cases.

Huynh, H. T.↗