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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

First results of a general circulation model applied to the SST-NOx problem

Results of model runs, one using two-dimensional distribution of NO2 in the unperturbed stratosphere and another including an additional localized source of NO2 to approximate the effect of SSTs, are reported. The general circulation model and chemical reactions are described, and corrections in the previous latitudinal and seasonal gradients of total columnar ozone and the diffusion coefficient in the neighborhood of the tropopause are noted. Excellent agreement with previous observations was obtained. Global distributions of ozone and NO2 are described and represented in graphs. Results indicate ozone reduction of approximately 16% in the Northern Hemisphere and approximately 8% in the Southern Hemisphere, with the mid-latitude source of NO2 apparently having a blocking effect on the horizontal transport of ozone, resulting in larger reductions of ozone at high latitudes than at low ones.

Cunnold, D. M.↗
NASA NTRS1980

Design studies of Laminar Flow Control (LFC) wing concepts using superplastics forming and diffusion bonding (SPF/DB)

Alternate concepts and design approaches were developed for suction panels and techniques were defined for integrating these panel designs into a complete LFC 200R wing. The design concepts and approaches were analyzed to assure that they would meet the strength, stability, and internal volume requirements. Cost and weight comparisions of the concepts were also made. Problems of integrating the concepts into a complete aircraft system were addressed. Methods for making splices both chordwise and spanwise, fuel light joints, and internal duct installations were developed. Manufacturing problems such as slot aligment, tapered slot spacing, production methods, and repair techniques were addressed. An assessment of the program was used to developed recommendations for additional research in the development of SPF/DB for LFC structure.

Wilson, V. E.↗
NASA NTRS1983

Doubly diffusive magnetic buoyancy instability in the solar interior

An investigation of the buoyancy of diffuse magnetic fields has shown that in the presence of rotation, static equilibrium configurations of the toroidal magnetic field and ambient plasma can exist. In that case, the escape of toroidal magnetic flux from the solar interior may be determined by the growth of instabilities which the equilibrium configuration may be subject to. In connection with the present investigation, it is assumed that in the region of toroidal magnetic flux amplification, the magnetic field has not as yet filamented into flux ropes, and is therefore 'diffuse'. A study is conducted of the MHD stability of an electrically conducting and differentially rotating gas in the presence of a toroidal magnetic field, an external constant gravitational field, and radiance pressure. The full dispersion relation for the magnetic buoyancy problem is developed, and the solutions of the dispersion relation are discussed.

Schmitt, J. H. M. M.↗
NASA NTRS1984

On the acceleration of energetic ions in Jupiter's magnetosphere

Several aspects of the problem of high-energy ions in the Jovian magnetosphere are addressed. Voyager observations pertaining to the problem of high-energy ions in the magnetosphere are summarized, and the charge exchange emission of fast neutral sulfur and oxygen atoms and their subsequent recapture by electron impact, charge exchange, and photoionization is considered. Solutions are given to the diffusion equation assuming a source of ions injected with a gyroenergy corresponding to pickup in the middle and outer magnetosphere. It is concluded that no reasonable model parameters exist to produce the required steep spectra of the particle observations with only pickup and adiabatic radial diffusion included. A local acceleration mechanism based on nonadiabatic wave-particle interactions is needed. The assumptions and model predictions of stochastic acceleration by MHD turbulence for the Jovian magnetosphere are described. The model makes a specific correspondence between MHD wave spectrum properties and particle spectrum properties at energies above the Alfven energy.

Barbosa, D. D.↗
NASA NTRS1985

Calculation of three-dimensional, viscous flow through turbomachinery blade passages by parabolic marching

The three-dimensional compressible Navier-Stokes equations are formulated in a rotating coordinate system, so as to include centrifugal and Coriolis forces. The equations are parabolized by using a previously calculated inviscid static pressure field. The thin layer Navier-Stokes approximation, which neglects streamwise diffusion, is used. A body-fitted coordinate system is used. The streamwise momentum equation is uncoupled from the cross-stream momentum equation by using contravariant momentum components, and then using the contravariant velocity components as primary unknowns. To reduce problems with small separated regions, the Reyhner and Flugge-Lotz approximation is used. The energy equation is included to allow for calculation of heat transfer. The flow may be laminar, or a simple eddy-viscosity turbulence may be used. A number of curved ducts and an axial stator were analyzed, including cases for which experimental data are available.

Katsanis, T.↗
NASA NTRS1985

Calculation of three-dimensional, viscous flow through turbomachinery blade passages by parabolic marching

The three-dimensional compressible Navier-Stokes equations are formulated in a rotating coordinate system, so as to include centrifugal and Coriolis forces. The equations are parabolized by using a previously calculated inviscid static pressure field. The thin layer Navier-Stokes approximation, which neglects streamwise diffusion, is used. A body-fitted coordinate system is used. The streamwise momentum equation is uncoupled from the cross-stream momentum equation by using contravariant momentum components, and then using the contravariant velocity components as primary unknowns. To reduce problems with small separating regions, the Reyhner and Flugge-Lotz approximation is used. The energy equation is included to allow for calculation of heat transfer. The flow may be laminar, or a simple eddy-viscosity turbulence may be used. A number of curved ducts and an axial stator were analyzed, including cases for which experimental data are available.

Katsanis, T.↗
NASA NTRS1987

Thermochemical erosion

A quasi-steady multidimensional Stefan problem is treated by means of similarity, asymptotic methods, and numerical analysis. The continuous wearing away of a solid 'erodee' is effected by gradual consumption of a solid 'eroder'. In the process under study, the erodee melts and the (liquid) melt diffuses to the eroder, where an exothermic, indefinitely rapid, diffusion-controlled surface reaction occurs. The products of reaction convey some of the chemically-generated heat to the surface of the erodee, so that further melting may occur, and the process may proceed. The analysis seeks to estimate the quasi-steady rates of consumption of the erodee and eroder, for the (convenient) self-similar case of confocal parabolas defining the boundaries of the melt layer. The analysis is simplified by observing that the flow is of lubrication-theory type: creeping, primarily unidirectional motion owing to the thinness of the melt layer between the erodee and eroder, with diffusion of heat and mass predominantly perpendicular to the principal motion. The analysis is completed by consideration of the force balance on the eroder.

Butler, Gerald↗
NASA NTRS2017

Decoupled Method for Reconstruction of Surface Conditions From Internal Temperatures On Ablative Materials With Uncertain Recession Model

Obtaining measurements of flight environments on ablative heat shields is both critical for spacecraft development and extremely challenging due to the harsh heating environment and surface recession. Thermocouples installed several millimeters below the surface are commonly used to measure the heat shield temperature response, but an ill-posed inverse heat conduction problem must be solved to reconstruct the surface heating environment from these measurements. Ablation can contribute substantially to the measurement response making solutions to the inverse problem strongly dependent on the recession model, which is often poorly characterized. To enable efficient surface reconstruction for recession model sensitivity analysis, a method for decoupling the surface recession evaluation from the inverse heat conduction problem is presented. The decoupled method is shown to provide reconstructions of equivalent accuracy to the traditional coupled method but with substantially reduced computational effort. These methods are applied to reconstruct the environments on the Mars Science Laboratory heat shield using diffusion limit and kinetically limited recession models.

Oliver, A. Brandon↗
NASA NTRS1984

Computation of viscous flow in planar and axisymmetric ducts by an implicit marching procedure

A streamwise marching procedure, approximately 200 times faster than a full Navier-Stokes procedure with comparable accuracy, is presented for solving problems of compressible viscous subsonic flow. Results are presented and compared with experimental data for the cases of developing turbulent flow in a circular pipe; turbulent flow in a two-dimensional S-duct; and turbulent flow in a typical subsonic diffuser. Prior to each main marching step, a preliminary marching step is taken in which the integral continuity equation and an uncoupled form of the streamwise momentum equation are solved simultaneously to obtain the viscous pressure correction. During the main marching step the equations for continuity, streamwise momentum, cross-flow momentum, and energy are solved simultaneously as a coupled system using an implicit finite-difference method, with the viscous pressure correction treated as a source term. The analysis may be used for flows with both favorable and adverse pressure gradients and to predict the location of flow separation.

Towne, C. E.↗
NASA NTRS1986

Magnetized thermal conduction fronts

The evolution of planar thermal conduction fronts in the presence of a dynamically weak, but otherwise self-consistent, magnetic field is considered. The field is assumed to be connected and untangled. In the diffusion limit for the thermal conductivity, these fronts exhibit self-similar behavior, even in the presence of a field. The role of the field is restricted to channeling the heat flux along its lines of force, and it enters into the problem as a dimensionless angle variable. 'Combing' (or opening) of insulating field lines by the evaporative flow is explicitly demonstrated. Unless the field is nearly perpendicular to the front normal in the hot gas, insulating effects are not profound. Self-similarity breaks down if the front becomes saturated, and under certain conditions magnetized saturated conduction fronts cannot propagate: the solution characteristics of the wave equation form caustics. The physical resolution is the advent of two-fluid (nonlocal) heating. Such Coulomb-heated fronts are expected to be relatively rare in typical astrophysical systems. The large-scale effects of a magnetic field on cloud evaporation in the interstellar medium are briefly discussed, and it is suggested that these fields preclude the presence of time-independent evaporative solutions. Thermal interfaces may then continue to evolve until radiative cooling halts their development; large tracts of warm 10,000 K gas may result.

Balbus, S. A.↗
NASA NTRS1993

An adaptive grid with directional control

An adaptive grid generator for adaptive node movement is here derived by combining a variational formulation of Winslow's (1981) variable-diffusion method with a directional control functional. By applying harmonic-function theory, it becomes possible to define conditions under which there exist unique solutions of the resulting elliptic equations. The results obtained for the grid generator's application to the complex problem posed by the fluid instability-driven magnetic field reconnection demonstrate one-tenth the computational cost of either a Eulerian grid or an adaptive grid without directional control.

Brackbill, J. U.↗
NASA NTRS1995

The Method of Space-time Conservation Element and Solution Element: Development of a New Implicit Solver

The method of space-time conservation element and solution element is a nontraditional numerical method designed from a physicist's perspective, i.e., its development is based more on physics than numerics. It uses only the simplest approximation techniques and yet is capable of generating nearly perfect solutions for a 2-D shock reflection problem used by Helen Yee and others. In addition to providing an overall view of the new method, we introduce a new concept in the design of implicit schemes, and use it to construct a highly accurate solver for a convection-diffusion equation. It is shown that, in the inviscid case, this new scheme becomes explicit and its amplification factors are identical to those of the Leapfrog scheme. On the other hand, in the pure diffusion case, its principal amplification factor becomes the amplification factor of the Crank-Nicolson scheme.

Chang, S. C.↗
NASA NTRS1995

An Extended Lagrangian Method

A unique formulation of describing fluid motion is presented. The method, referred to as 'extended Lagrangian method,' is interesting from both theoretical and numerical points of view. The formulation offers accuracy in numerical solution by avoiding numerical diffusion resulting from mixing of fluxes in the Eulerian description. The present method and the Arbitrary Lagrangian-Eulerian (ALE) method have a similarity in spirit-eliminating the cross-streamline numerical diffusion. For this purpose, we suggest a simple grid constraint condition and utilize an accurate discretization procedure. This grid constraint is only applied to the transverse cell face parallel to the local stream velocity, and hence our method for the steady state problems naturally reduces to the streamline-curvature method, without explicitly solving the steady stream-coordinate equations formulated a priori. Unlike the Lagrangian method proposed by Loh and Hui which is valid only for steady supersonic flows, the present method is general and capable of treating subsonic flows and supersonic flows as well as unsteady flows, simply by invoking in the same code an appropriate grid constraint suggested in this paper. The approach is found to be robust and stable. It automatically adapts to flow features without resorting to clustering, thereby maintaining rather uniform grid spacing throughout and large time step. Moreover, the method is shown to resolve multi-dimensional discontinuities with a high level of accuracy, similar to that found in one-dimensional problems.

Liou, Meng-Sing↗
NASA NTRS2002

Method for improving the performance of oxidizable ceramic materials in oxidizing environments

Improved adhesion of thermal barrier coatings to nonmetallic substrates using a dense layer of ceramic on an underlying nonmetallic substrate that includes at least one oxidizable component. The improved adhesion occurs because the application of the dense ceramic layer forms a diffusion barrier for oxygen. This diffusion barrier prevents the oxidizable component of the substrate from decomposing. The present invention applies ceramic by a process that deposits a relatively thick and dense ceramic layer on the underlying substrate. The formation of the dense layer of ceramic avoids the problem of void formation associated with ceramic formation by most prior art thermal decomposition processes. The formation of voids has been associated with premature spalling of thermal barrier layers and other protective layers applied to substrates.

Nagaraj, Bangalore A.↗
NASA NTRS1989

Stochastic histories of dust grains in the interstellar medium

The purpose is to study an evolving system of refractory dust grains within the Interstellar Medium (ISM). This is done via a combination of Monte Carlo processes and a system of partial differential equations, where refractory dust grains formed within supernova remnants and ejecta from high mass loss stars are subjected to the processes of sputtering and collisional fragmentation in the diffuse media and accretion within the cold molecular clouds. In order to record chemical detail, the authors take each new particle to consist of a superrefractory core plus a more massive refractory mantle. The particles are allowed to transfer to and fro between the different phases of the interstellar medium (ISM) - on a time scale of 10(exp 8) years - until either the particles are destroyed or the program finishes at a Galaxy time of 6x10(exp 9) years. The resulting chemical and size spectrum(s) are then applied to various astrophysical problems with the following results. For an ISM which has no collisional fragmentation of the dust grains, roughly 10 percent by mass of the most refractory material survives the rigors of the ISM intact, which leaves open the possibility that fossilized isotopically anomalous material may have been present within the primordial solar nebula. Stuctured or layered refractory dust grains within the model cannot explain the observed interstellar depletions of refractory material. Fragmentation due to grain-grain collisions in the diffuse phase plus the accretion of material in the molecular cloud phase can under certain circumstances cause a bimodal distribution in grain size.

Liffman, Kurt↗
NASA NTRS2010

Fabrication of Lanthanum Telluride 14-1-11 Zintl High-Temperature Thermoelectric Couple

The development of more efficient thermoelectric couple technology capable of operating with high-grade heat sources up to 1,275 K is key to improving the performance of radioisotope thermoelectric generators. Lanthanum telluride La3-xTe4 and 14-1-11 Zintls (Yb14MnSb11) have been identified as very promising materials. The fabrication of advanced high-temperature thermoelectric couples requires the joining of several dissimilar materials, typically including a number of diffusion bonding and brazing steps, to achieve a device capable of operating at elevated temperatures across a large temperature differential (up to 900 K). A thermoelectric couple typically comprises a heat collector/ exchanger, metallic interconnects on both hot and cold sides, n-type and ptype conductivity thermoelectric elements, and cold-side hardware to connect to the cold-side heat rejection and provide electrical connections. Differences in the physical, mechanical, and chemical properties of the materials that make up the thermoelectric couple, especially differences in the coefficients of thermal expansion (CTE), result in undesirable interfacial stresses that can lead to mechanical failure of the device. The problem is further complicated by the fact that the thermoelectric materials under consideration have large CTE values, are brittle, and cracks can propagate through them with minimal resistance. The inherent challenge of bonding brittle, high-thermal-expansion thermoelectric materials to a hot shoe material that is thick enough to carry the requisite electrical current was overcome. A critical advantage over prior art is that this device was constructed using all diffusion bonds and a minimum number of assembly steps. The fabrication process and the materials used are described in the following steps: (1) Applying a thin refractory metal foil to both sides of lanthanum telluride. To fabricate the n-type leg of the advanced thermoelectric couple, the pre-synthesized lanthanum telluride coupon was diffusion bonded to the metal foil using a thin adhesion layer. (2) Repeating a similar process for the 14-1-11 Zintl p-type leg of the advanced thermoelectric couple. (3) Bonding thick CTE-matched metal plates on the metallized lanthanum telluride and Yb14MnSb11 to form the hot and cold sides of the thermoelectric couple. The calculated conversion efficiency of such an advanced couple would be about 10.5 percent, about 35 percent better than heritage radioisotope thermoelectric technology that relies on Si-Ge alloys. In addition, unlike Si-Ge alloys, these materials can be combined with many other thermoelectric materials optimized for operation at lower temperatures to achieve conversion efficiency in excess of 15 percent (a factor of 2 increase over heritage technology).

Ravi, Vilupanur A.↗
NASA NTRS1994

Major unresolved issues in the formation of chondrules and CAIs

So much has been learned and published about chondrites and so many models have been floated (but neither proven nor refuted) that the subject has become extremely diffuse and amorphous. It is impossible to remain simultaneously aware of all the many constraints, from many different subdisciplines, that bear on nebular processes. This paper attempts to impose some structure on the situation and loosen the interpretational logjam by breaking the chondrite problem into a series of well-defined questions. Most are open questions; those that some might consider closed by now may deserve formal reexamination.

Wood, J. A.↗
NASA NTRS1999

On Multi-Dimensional Unstructured Mesh Adaption

Anisotropic unstructured mesh adaption is developed for a truly multi-dimensional upwind fluctuation splitting scheme, as applied to scalar advection-diffusion. The adaption is performed locally using edge swapping, point insertion/deletion, and nodal displacements. Comparisons are made versus the current state of the art for aggressive anisotropic unstructured adaption, which is based on a posteriori error estimates. Demonstration of both schemes to model problems, with features representative of compressible gas dynamics, show the present method to be superior to the a posteriori adaption for linear advection. The performance of the two methods is more similar when applied to nonlinear advection, with a difference in the treatment of shocks. The a posteriori adaption can excessively cluster points to a shock, while the present multi-dimensional scheme tends to merely align with a shock, using fewer nodes. As a consequence of this alignment tendency, an implementation of eigenvalue limiting for the suppression of expansion shocks is developed for the multi-dimensional distribution scheme. The differences in the treatment of shocks by the adaption schemes, along with the inherently low levels of artificial dissipation in the fluctuation splitting solver, suggest the present method is a strong candidate for applications to compressible gas dynamics.

Wood, William A.↗