Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffused interface model”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Variable Property Flow in Rectangular Ducts with Repeated Rectangular Rib Roughness

PHOENICS was used to simulate turbulent flow of supercritical hydrogen through asymmetrically heated rectangular ducts with one side rib-roughened, in order to provide design data for active cooling systems of hypersonic aircraft where the hydrogen fuel is used as coolant. Variable property and conjugate heat transfer effects were considered. The computations using PHOENICS required provision of property sub-routines for supercritical hydrogen, a low-Reynolds number form of the k-epsilon turbulence model, a rough-wall wall function and a procedure for calculating proper diffusion coefficients at the fluid-solid interface. The effect of each design variable such as width and height of duct and wall thickness, on pressure drop and maximum wall temperature was investigated. PHOENICS version 1.4 was used on an IBM 3090 computer.

Youn, B.↗
NASA NTRS1993

A numerical treatment of geodynamic viscous flow problems involving the advection of material interfaces

Effective numerical treatment of multicomponent viscous flow problems involving the advection of sharp interfaces between materials of differing physical properties requires correction techniques to prevent spurious diffusion and dispersion. We develop a particular algorithm, based on modern shock-capture techniques, employing a two-step nonlinear method. The first step involves the global application of a high-order upwind scheme to a hyperbolic advection equation used to model the distribution of distinct material components in a flow field. The second step is corrective and involves the application of a global filter designed to remove dispersion errors that result from the advection of discontinuities (e.g., material interfaces) by high-order, minimally dissipative schemes. The filter introduces no additional diffusion error. Nonuniform viscosity across a material interface is allowed for by the implementation of a compositionally weighted-inverse interface viscosity scheme. The combined method approaches the optimal accuracy of modern shock-capture techniques with a minimal increase in computational time and memory. A key advantage of this method is its simplicity to incorporate into preexisting codes be they finite difference, element, or volume of two or three dimensions.

Lenardic, A.↗
NASA NTRS2000

The Turbulent/Non-Turbulent Interface Bounding a Far-Wake

The velocity fields of a turbulent wake behind a flat plate obtained from the direct numerical simulations of Moser et al. are used to study the structure of the flow in the intermittent zone where there are, alternately, regions of fully turbulent flow and non-turbulent velocity fluctuations either side of a thin randomly moving interface. Comparisons are made with a wake that is 'forced' by amplifying initial velocity fluctuations. There is also a random temperature field T in the flow; T varies between constant values of 0.0 and 1.0 on the sides of the wake. The value of the Reynolds number based on the centreplane mean velocity defect and halfwidth b of the wake is Re approx. = 2000. It is found that the thickness of the continuous interface is about equal to 0.07b, whereas the amplitude of fluctuations of the instantaneous interface displacement y(sub I)(t) is an order of magnitude larger, being about 0.5b. This explains why the mean statistics of vorticity in the intermittent zone can be calculated in terms of the probability distribution of y(sub I) and the instantaneous discontinuity in vorticity across the interface. When plotted as functions of y - y(sub I), the conditional mean velocity (U) and temperature (T) profiles show sharp jumps Delta(U) and Delta(T) at the interface adjacent to a thick zone where (U) and (T) vary much more slowly. Statistics for the vorticity and velocity variances, available in such detail only from DNS data, show how streamwise and spanwise components of vorticity are generated by vortex stretching in the bulges of the interface. Flow fields around the interface, analyzed in terms of the local streamline pattern, confirm previous results that the advancement of the vortical interface into the irrotational flow is driven by large-scale eddy motion. It is argued that because this is an inviscid mechanism the entrainment process is not sensitive to the value of Re, and that small-scale nibbling only plays a subsidiary role. While mean Reynolds stresses decrease gradually in the intermittent zone, conditional stresses are found to decrease sharply towards zero at the interface. Using one-point turbulence models applied to either unconditional or conditional statistics for the turbulent region and then averaged, the entrainment rate E(sub b) would, if calculated exactly, be zero. But if computed with standard computational methods, E(sub b) would be non-zero because of numerical diffusion. It is concluded that the current practice in statistical models of approximating entrainment by a diffusion process is computationally arbitrary and physically incorrect. An analysis shows how E(sub b) is related to Delta(U) and the jump in shear stress at the interface, and correspondingly to Delta(T) and the heat flux.

Bisset, David K.↗
NASA NTRS1985

Study of the kinetics and mechanism of the thermal nitridation of SiO2

X-ray photoelectron spectroscopy (XPS) has been used to study the nitridation time and temperature dependence of the nitrogen distribution in thermally nitrided SiO2 films. The XPS data show that the maximum nitrogen concentration near the (SiO(x)N(y)/Si interface is initially at the interface, but moves 20-25 A away from the interface with increasing nitridation time. Computer modeling of the kinetic processes involved is carried out and reveals a mechanism in which diffusing species, initially consisting primarily of nitrogen, react with the substrate, followed by formation of the oxygen-rich oxynitride due to reaction of the diffusing oxygen displaced by the slower nitridation of the SiO2. The data are consistent with this mechanism provided the influence of the interfacial strain on the nitridation and oxidation kinetics is explicitly accounted for.

Vasquez, R. P.↗
NASA NTRS1990

Growth morphology with anisotropic surface kinetics

The morphological evolution of crystals growing from an incongruent vapor phase is studied using a Monte Carlo model, and the full range of growth morphologies is recovered. The diffusion in the bulk nutrient and the anisotropy in the interface kinetics are morphologically destabilizing and stabilizing, respectively. For a given set of simulation parameters and lattice symmetries there is a critical size, which scales linearly with the mean free path in the vapor, beyond which a crystal cannot retain its stable, macroscopically faceted growth shape. Surface diffusion stabilizes faceted growth on the shorter scale of the mean surface diffusion length. In simulations with a uniform drift superimposed on the random walk nutrient transport, crystal faces oriented toward the drift show enhanced morphological stability compared to the purely diffusive situation. Rotational drifts with periodic reversal of direction are morphologically stabilizing for all crystal facets.

Xiao, Rong-Fu↗
NASA NTRS1999

The Dissolution of an Interfween Miscible Liquids

The disappearance of the surface tension of the interface of a binary mixture, measured using the dynamic surface light scattering technique, is slower for a binary mixture of higher density contrast. A comparison with a naive diffusion model, expected to provide a lower limit for the speed of dissolution in the absence of gravity shows that the interfacial surface tension disappears much slower than even by diffusion with the effect becoming much more pronounced when density contrast between the liquid phases is increased. Thus, the factor most likely to be responsible for this anomalously slow dissolution is gravity. A mechanism could be based on the competition between diffusive relaxation and sedimentation at the dissolving interface.

Vlad, D.H.↗
NASA NTRS1998

3-D Modeling of Double-Diffusive Convection During Directional Solidification of a Non-Dilute Alloy with Application to the HgCdTe Growth Under Microgravity Conditions

A numerical calculation for a non-dilute alloy solidification was performed using the FIDAP finite element code. For low growth velocities plane front solidification occurs. The location and the shape of the interface was determined using melting temperatures from the HgCdTe liquidus curve. The low thermal conductivity of the solid HgCdTe causes thermal short circuit through the ampoule walls, resulting in curved isotherms in the vicinity of the interface. Double-diffusive convection in the melt is caused by radial temperature gradients and by material density inversion with temperature. Cooling from below and the rejection at the solid-melt interface of the heavier HgTe-rich solute each tend to reduce convection. Because of these complicating factors dimensional rather then non-dimensional modeling was performed. Estimates of convection contributions for various gravity conditions was performed parametrically. For gravity levels higher then 1 0 -7 of earth's gravity it was found that the maximum convection velocity is extremely sensitive to gravity vector orientation and can be reduced at least by factor of 50% for precise orientation of the ampoule in the microgravity environment. The predicted interface shape is in agreement with one obtained experimentally by quenching. The results of 3-D modeling are compared with previous 2-D finding. A video film featuring melt convection will be presented.

Bune, Andris V.↗
NASA NTRS1999

Effects of Gravity on the Double-Diffusive Convection During Directional Solidification of a Non-Dilute Alloy with Application to HgCdTe

A general 2-D and 3-D finite element model of non-dilute alloy solidification was used to simulate growth of HgCdTe in terrestrial and microgravity conditions. Verification of the 3-D model was undertaken by comparison with previously published results on convection in an inclined cylinder. For low growth velocities, plane front solidification occurs. The location and the shape of the interface were determined using melting temperatures obtained from the HgCdTe liquidus curve. The low thermal conductivity of the solid HgCdTe causes a thermal short circuit through the ampoule walls, resulting in curved isotherms in the vicinity of the interface. Double-diffusive convection in the melt is caused by radial temperature gradients and by material density inversion due to the combined effects of composition and temperature. Cooling from below and the rejection at the solid-melt interface of the heavier HgTe-rich solute each tend to reduce convection. Because of these complicating factors, dimensional rather than non-dimensional modeling was performed. the predicted interface shape is in agreement with one obtained experimentally by quenching.

Bune, Andris V.↗
NASA NTRS2006

A Theoretical Solid Oxide Fuel Cell Model for System Controls and Stability Design

As the aviation industry moves towards higher efficiency electrical power generation, all electric aircraft, or zero emissions and more quiet aircraft, fuel cells are sought as the technology that can deliver on these high expectations. The Hybrid Solid Oxide Fuel Cell system combines the fuel cell with a microturbine to obtain up to 70 percent cycle efficiency, and then distributes the electrical power to the loads via a power distribution system. The challenge is to understand the dynamics of this complex multi-discipline system, and design distributed controls that take the system through its operating conditions in a stable and safe manner while maintaining the system performance. This particular system is a power generation and distribution system and the fuel cell and microturbine model fidelity should be compatible with the dynamics of the power distribution system in order to allow proper stability and distributed controls design. A novel modeling approach is proposed for the fuel cell that will allow the fuel cell and the power system to be integrated and designed for stability, distributed controls, and other interface specifications. This investigation shows that for the fuel cell, the voltage characteristic should be modeled, but in addition, conservation equation dynamics, ion diffusion, charge transfer kinetics, and the electron flow inherent impedance should also be included.

Kopasakis, George↗
NASA NTRS2008

A Theoretical Solid Oxide Fuel Cell Model for Systems Controls and Stability Design

As the aviation industry moves toward higher efficiency electrical power generation, all electric aircraft, or zero emissions and more quiet aircraft, fuel cells are sought as the technology that can deliver on these high expectations. The hybrid solid oxide fuel cell system combines the fuel cell with a micro-turbine to obtain up to 70% cycle efficiency, and then distributes the electrical power to the loads via a power distribution system. The challenge is to understand the dynamics of this complex multidiscipline system and the design distributed controls that take the system through its operating conditions in a stable and safe manner while maintaining the system performance. This particular system is a power generation and a distribution system, and the fuel cell and micro-turbine model fidelity should be compatible with the dynamics of the power distribution system in order to allow proper stability and distributed controls design. The novelty in this paper is that, first, the case is made why a high fidelity fuel cell mode is needed for systems control and stability designs. Second, a novel modeling approach is proposed for the fuel cell that will allow the fuel cell and the power system to be integrated and designed for stability, distributed controls, and other interface specifications. This investigation shows that for the fuel cell, the voltage characteristic should be modeled but in addition, conservation equation dynamics, ion diffusion, charge transfer kinetics, and the electron flow inherent impedance should also be included.

Kopasakis, George↗
NASA NTRS1980

Liquid-phase electroepitaxy - Dopant segregation

A theoretical model is presented which accounts for the dopant segregation in liquid-phase electroepitaxy in terms of dopant transport in the liquid phase (by electromigration and diffusion), the growth velocity, and the Peltier effect at the substrate-solution interface. The contribution of dopant electromigration to the magnitude of the effective segregation coefficient is dominant in the absence of convection; the contribution of the Peltier effect becomes significant only in the presence of pronounced convection. Quantitative expressions which relate the segregation coefficient to the growth parameters also permit the determination of the diffusion constant and electromigration mobility of the dopant in the liquid phase. The model was found to be in good agreement with the measured segregation characteristics of Sn in the electroepitaxial growth of GaAs from Ga-As solutions. For Sn in Ga-As solution at 900 C the diffusion constant was found to be 4 x 10 to the -5 sq cm/s and the electromigration velocity (toward the substrate with a positive polarity 2 x 10 to the -5 cm/s current density of 10 A/sq cm.

Lagowski, J.↗
NASA NTRS1980

Electroepitaxy of multicomponent systems - Ternary and quarternary compounds

A theoretical model is presented which accounts for the electroepitaxial growth kinetics and composition of multicomponent compounds in terms of mass transport in the liquid and phase diagram relationships. The mass transport in the interface is dominated by electromigration in the absence of convection and by diffusion in the presence of convection. The composition of the solid is controlled by the Peltier effect at the growth interface and by the diffusion and mobility constants of the solute components and the growth velocity (current density). Thus, for a given solution composition, the composition of the solid can be varied by varying the current density. For a given current density the composition remains constant even in the case of relatively thick epitaxial layers. All aspects of the model were found to be in good agreement with the growth and composition characteristics of Ga/x-1/Al/x/As layers.

Bryskiewicz, T.↗
NASA NTRS1989

Finite element modeling of frictionally restrained composite interfaces

The use of special interface finite elements to model frictional restraint in composite interfaces is described. These elements simulate Coulomb friction at the interface, and are incorporated into a standard finite element analysis of a two-dimensional isolated fiber pullout test. Various interfacial characteristics, such as the distribution of stresses at the interface, the extent of slip and delamination, load diffusion from fiber to matrix, and the amount of fiber extraction or depression are studied for different friction coefficients. The results are compared to those obtained analytically using a singular integral equation approach, and those obtained by assuming a constant interface shear strength. The usefulness of these elements in micromechanical modeling of fiber-reinforced composite materials is highlighted.

Ballarini, Roberto↗
NASA NTRS2001

The Isothermal Dendritic Growth Experiment (IDGE)

Dendritic solidification is one of the simplest examples of pattern formation where a structureless melt evolves into a ramified crystalline microstructure; it is a common mode of solidification in many materials, but especially so in metals and alloys. There is considerable engineering interest in dendrites because of the role dendrites play in the determination of microstructure, and thereby in influencing the physical properties of cast metals and alloys. Dendritic solidification provides important examples of non-equilibrium physics, pattern formation dynamics, and models for computational condensed matter and material physics. Current theories of dendritic growth generally couple diffusion effects in the melt with the physics introduced by the interface. Unfortunately, in terrestrial based experiments, convective effects in the melt alter the growth process in such a manner as to prevent definitive analysis of convective, diffusive or interfacial effects. Thus, the effective elimination of convection in the melt by operating experiments on orbit were required to produce high-fidelity data needed for achieving further progress. This simple fact comprised the scientific justification for the IDGE.

Glicksman, M. E.↗
NASA NTRS1985

Thermoplastic matrix composite processing model

The effects the processing parameters pressure, temperature, and time have on the quality of continuous graphite fiber reinforced thermoplastic matrix composites were quantitatively accessed by defining the extent to which intimate contact and bond formation has occurred at successive ply interfaces. Two models are presented predicting the extents to which the ply interfaces have achieved intimate contact and cohesive strength. The models are based on experimental observation of compression molded laminates and neat resin conditions, respectively. Identified as the mechanism explaining the phenomenon by which the plies bond to themselves is the theory of autohesion (or self diffusion). Theoretical predictions from the Reptation Theory between autohesive strength and contact time are used to explain the effects of the processing parameters on the observed experimental strengths. The application of a time-temperature relationship for autohesive strength predictions is evaluated. A viscoelastic compression molding model of a tow was developed to explain the phenomenon by which the prepreg ply interfaces develop intimate contact.

Dara, P. H.↗
NASA NTRS1990

Creep and creep rupture of strongly reinforced metallic composites

A creep and creep damage theory is presented for metallic composites with strong fibers. Application is to reinforced structures in which the fiber orientation may vary throughout but a distinct fiber direction can be identified locally (local transverse isotropy). The creep deformation model follows earlier work and is based on a flow potential function that depends on invariants reflecting stress and the material symmetry. As the focus is on the interaction of creep and damage, primary creep is ignored. The creep rupture model is an extension of continuum damage mechanics and includes an isochronous damage function that depends on invariants specifying the local maximum transverse tension and the maximum longitudinal shear stress. It is posited that at high temperature and low stress, appropriate to engineering practice, these stress components damage the fiber/matrix interface through diffusion controlled void growth, eventually causing creep rupture. Experiments are outlined for characterizing a composite through creep rupture tests under transverse tension and longitudinal shear. Application is made to a thin-walled pressure vessel with reinforcing fibers at an arbitrary helical angle. The results illustrate the usefulness of the model as a means of achieving optimal designs of composite structures where creep and creep rupture are life limiting.

Robinson, D. N.↗
NASA NTRS1980

Frequency response of charge transfer in MOS inversion layers

The dynamics of charge transfer from a reservoir into an MOS inversion layer, which limits the frequency response of an MOS transistor or a charge-coupled device, is investigated. Using Berman and Kerr's model of space-charge capacitance in the semiconductor, a small-signal distributed model is developed for an MOS structure which transfers charge in an inversion channel due to a variation in the gate voltage. The dynamics of the charge transfer is characterized by a time constant which is determined by the length of the inversion channel and its mobility. Experimental data of gate capacitance vs frequency, taken from a test structure with a diffused source/drain well, are satisfactorily fitted by theoretical curves derived from the model. The channel mobility is precisely determined from the adjusted time constant. The influence of interface states on the capacitance-frequency relationship is also briefly discussed.

Lieneweg, U.↗
NASA NTRS1986

Oxidation and embrittlement of Ti-6Al-2Sn-4Zr-2Mo alloy

A thermogravimetric analysis monitoring study has been conducted for Ti-6Al-2Sn-4Zr-2Mo in the 593-760 C range, with a view to its total oxidation kinetics in laboratory air. In addition, tensile elongation was obtained for selected foil specimens after their exposure to simulated Space Shuttle reentry conditions. The total weight gain with time exhibited two parabolic stages separated by a transient region which is the result of a simultaneous increase in the parabolic growth rate constant for the oxide and in oxygen solubility at the oxide-metal interface. A model for the total oxidation kinetics accounting for the two individual components (oxide growth and solid solution formation) is proposed; the activation energy for oxygen diffusion in the alpha-solid solution is shown to approximate the activation energy for the degradation of tensile elongation of the alloy in its foil gage condition.

Shenoy, R. N.↗