Search NASA⌕ Search

SEARCH · Search NASA

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Spherical Diffusion Flames: Structure and Dynamics

The spherical geometry is the most suitable one to study fundamental issues concerning the structure and the dynamics of diffusion flames. From a theoretical point of view, it is the only geometry that permits the existence of a truly one dimensional stationary diffusion flame. A stationary planar diffusion flame with the fuel supplied upstream, at x approaches -infinity say, and the oxidant downstream, at x approaches +infinity is not possible. For a steady diffusion flame to exist, one must have nonzero fluxes of fuel and oxidant towards the flame. However, in the unlimited region behind the planar flame the only bounded solutions to the reaction-free convective-diffusive operator are constants. Hence the oxidant concentration behind the flame remains constant and there is no mechanism to generate the necessary flux towards the flame. A one-dimensional problem can be formulated if the reactants are supplied at finite locations; but the boundary conditions in this case introduce unnecessary complications and do not appropriately model the physical reality. Indeed, a planar diffusion flame can be established in the stagnation-point flow of two opposed jets but the flame in this case is stretched and the flow is essentially two-dimensional. The only stationary one-dimensional diffusion flame in an unlimited environment is therefore the spherical flame.

Matalon, Moshe↗
NASA NTRS1993

Applying and validating the RANS-3D flow-solver for evaluating a subsonic serpentine diffuser geometry

Subsonic inlet ducts for advanced, high-performance aircraft are evolving towards complex three-dimensional shapes for reasons of overall integration and weight. These factors lead to diffuser geometries that may sacrifice inlet performance, unless careful attention to design details and boundary layer management techniques are employed. The ability of viscous computational fluid dynamic (CFD) analysis of such geometries to aid the aircraft configurator in this complex design problem is herein examined. The RANS-3D Reynolds-Averaged Navier-Stokes solver is applied to model the complex flowfield occurring in a representative diffuser geometry and the solutions are compared to experimental results from a static test of the inlet duct. The computational results are shown to compare very favorably with experimental results over a range of mass flow rates, including those involving large amounts of separation in the diffuser. In addition, a novel grid topology is presented, and two turbulence models are evaluated in this study as part of the RANS-3D code.

Fletcher, Michael J.↗
NASA NTRS1973

A review of methods for predicting air pollution dispersion

Air pollution modeling, and problem areas in air pollution dispersion modeling were surveyed. Emission source inventory, meteorological data, and turbulent diffusion are discussed in terms of developing a dispersion model. Existing mathematical models of urban air pollution, and highway and airport models are discussed along with their limitations. Recommendations for improving modeling capabilities are included.

Mathis, J. J., Jr.↗
NASA NTRS1987

An Exponential Finite Difference Technique for Solving Partial Differential Equations

An exponential finite difference algorithm, as first presented by Bhattacharya for one-dimensianal steady-state, heat conduction in Cartesian coordinates, has been extended. The finite difference algorithm developed was used to solve the diffusion equation in one-dimensional cylindrical coordinates and applied to two- and three-dimensional problems in Cartesian coordinates. The method was also used to solve nonlinear partial differential equations in one (Burger's equation) and two (Boundary Layer equations) dimensional Cartesian coordinates. Predicted results were compared to exact solutions where available, or to results obtained by other numerical methods. It was found that the exponential finite difference method produced results that were more accurate than those obtained by other numerical methods, especially during the initial transient portion of the solution. Other applications made using the exponential finite difference technique included unsteady one-dimensional heat transfer with temperature varying thermal conductivity and the development of the temperature field in a laminar Couette flow.

Handschuh, Robert F.↗
NASA NTRS1989

Effects of ion two-stream instability on auroral ion heating

The ion two-stream instability develops between upflowing H(+) and O(+) in the auroral acceleration region, if the two species are accelerated to different velocities in a semistationary electrostatic potential. Parallel propagating modes are unstable only for small relative drifts, while oblique modes remain unstable at larger drifts and, consequently, higher altitudes. The effects of oblique modes on ion heating are investigated, and it is found that O(+) heats primarily by trapping in the direction of wave propagation during the early, coherent phase of wave growth in the temporal evolution problem. After undergoing trapping along the magnetic field direction, H(+) heats primarily by quasi-linear diffusion when the spectrum broadens nonlinearly. The parallel ion distributions evolve by formation of a high-energy (low-energy) tail on oxygen (hydrogen) along with bulk slowing of hydrogen, which results in more energetic oxygen than hydrogen. Total heating is greatest in the direction of wave propagation, and this effect is greater for hydrogen than for oxygen.

Roth, I.↗
NASA NTRS1988

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.↗
NASA NTRS1993

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.↗
NASA NTRS1993

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.↗
NASA NTRS1990

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.↗
NASA NTRS1986

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.↗
NASA NTRS2011

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↗
NASA NTRS1990

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↗
NASA NTRS1993

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.↗
NASA NTRS1997

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↗
NASA NTRS1995

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↗
NASA NTRS1998

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↗
NASA NTRS1984

Vibration characteristics of the HPOTP (High-Pressure Oxygen Turbopump) of the SSME (Space Shuttle Main Engine)

Attention is given to rotor dynamic problems that have been encountered and eliminated in the course of Space Shuttle Main Engine (SSME) development, as well as continuing, subsynchronous problems which are being encountered in the development of a 109-percent power level engine. The basic model for the SSME's High Pressure Oxygen Turbopump (HPOTP) encompasses a structural dynamic model for the rotor and housing, and component models for the liquid and gas seals, turbine clearance excitation forces, and impeller diffuser forces. Linear model results are used to examine the synchronous response and stability characteristics of the HPOTP, with attention to bearing load and stability problems associated with the second critical speed. Differences between linear and nonlinear model results are discussed and explained in terms of simple models. Simulation results indicate that while synchronous bearing loads can be reduced, subsynchronous motion is not eliminated by seal modifications.

Childs, D. W.↗
DOE OSTI2024

A novel conditional generative model for efficient ensemble forecasts of state variables in large-scale geological carbon storage

Integrating monitoring data to efficiently update reservoir pressure and CO 2 plume distribution forecasts presents a significant challenge in geological carbon storage (GCS) applications. Inverse modeling techniques are commonly used to fuse observational data and refine reservoir model parameters, thereby improving state variable forecasts. However, these techniques often rely on linear or Gaussian assumptions, which can limit their effectiveness in accurately predicting state variables. Moreover, simulating large-scale three-dimensional (3D) GCS problems is computationally expensive, making iterative runs in inverse problems prohibitive. To address these challenges, we propose a conditional generative model utilizing the score-based diffusion method for real-time 3D pressure and saturation field distribution predictions. Our approach involves solving the score function with a mini-batch-based Monte Carlo estimator to generate labeled data. This data is subsequently employed to train a fully connected neural network, enabling it to learn the conditional sample generator within a supervised learning framework. This method enables the rapid generation of a large ensemble of predictions, facilitating comprehensive uncertainty quantification of state variables. Here we applied our method to forecast the dynamic 3D distributions of pressure and saturation fields over a 30-year injection period. The statistical assessment with low root mean square error (RMSE) values demonstrates that our method can accurately predict the spatiotemporal distributions of both pressure and saturation fields. Moreover, the developed conditional generative model shows high computational efficiency by generating 100 ensemble forecasts of 3D state variables in less than 10 min. The consistency between ensemble averages and ground truth values further illustrates the model’s capability to capture state variable dynamics during the CO 2 plume injection process. Notably, the ground truth values fall within the ensemble forecasts, indicating that our uncertainty quantification effectively captures variability and potential noise in the observations. Thus, the developed conditional generative model proves to be a more efficient, accurate, and practical tool for GCS applications, facilitating timely risk analysis and informed decision-making.

58 GEOSCIENCES↗