Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffuse interface methods”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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 NTRS2016

Modification of Surface Density of a Porous Medium

A method for increasing density of a region of a porous, phenolic bonded ("PPB") body adjacent to a selected surface to increase failure tensile strength of the adjacent region and/or to decrease surface recession at elevated temperatures. When the surface-densified PPB body is brought together with a substrate, having a higher failure tensile strength, to form a composite body with a PPB body/substrate interface, the location of tensile failure is moved to a location spaced apart from the interface, the failure tensile strength of the PPB body is increased, and surface recession of the material at elevated temperature is reduced. The method deposits and allows diffusion of a phenolic substance on the selected surface. The PPB body and the substrate may be heated and brought together to form the composite body. The phenolic substance is allowed to diffuse into the PPB body, to volatilize and to cure, to provide a processed body with an increased surface density.

Stackpoole, Margaret M.↗
NASA NTRS2023

Solidification Furnace for In Situ Observation of Bulk Transparent Systems and Image Analysis Methods

This paper aims to describe the experimental framework of the Directional Solidification Insert (DSI), installed onboard the International Space Station (ISS), dedicated to the in situ and real-time characterization of the dynamic selection of the solid-liquid interface morphology in bulk samples of transparent materials under diffusive growth conditions. The in situ observation of the solid-liquid interface is an invaluable tool for gaining knowledge on the time-evolution of the interface pattern because the initial morphological instability evolves nonlinearly and undergoes a reorganization process. The result of each experiment, characterized by the sample concentration, a thermal gradient and a pulling rate, is a large number of images. Interpretation of these images necessitates a robust identification of each cell/dendrite position and size during the entire solidification. Several image analysis methods have been developed to achieve this goal reliably despite varying contrast and noise levels, and are described in details. Typical solidification experiments are presented and the dynamics of the pattern formation are analyzed to illustrate the application of the image analysis methods.

Programming languages↗
NASA NTRS2005

X-Ray Fluorescence Solvent Detection at the Substrate-Adhesive Interface

With environmental regulations limiting the use of volatile organic compounds, low-vapor pressure solvents have replaced traditional degreasing solvents for bond substrate preparation. When used to clean and prepare porous bond substrates such as phenolic composites, low vapor pressure solvents can penetrate deep into substrate pore networks and remain there for extended periods. Trapped solvents can interact with applied adhesives either prior to or during cure, potentially compromising bond properties. Currently, methods for characterizing solvent time-depth profiles in bond substrates are limited to bulk gravimetric or sectioning techniques. While sectioning techniques such as microtome allow construction of solvent depth profiles, their depth resolution and reliability are limited by substrate type. Sectioning techniques are particularly limited near the adhesive-substrate interface where depth resolution is further limited by adhesive-substrate hardness and, in the case of a partially cured adhesive, mechanical properties differences. Additionally, sectioning techniques cannot provide information about lateral solvent diffusion. Cross-section component mapping is an alternative method for measuring solvent migration in porous substrates that eliminates the issues associated with sectioning techniques. With cross-section mapping, the solvent-wiped substrate is sectioned perpendicular rather than parallel to the wiped surface, and the sectioned surface is analyzed for the solvent or solvent components of interest using a two-dimensional mapping or imaging technique. Solvent mapping can be performed using either direct or indirect methods. With a direct method, one or more solvent components are mapped using red or Raman spectroscopy together with a moveable sample stage and/or focal plane array detector. With an indirect method, an elemental "tag" not present in the substrate is added to the solvent before the substrate is wiped. Following cross sectioning, the tag element can then be mapped by its characteristic x-ray emission using either x-ray fluorescence, or electron-beam energy-and wavelength-dispersive x-ray spectrometry. The direct mapping techniques avoid issues of different diffusion or migration rates of solvents and elemental tags, while the indirect techniques avoid spectral resolution issues in cases where solvents and substrates have adjacent or overlapping peaks. In this study, cross-section component indirect mapping is being evaluated as a method for measuring migration of d-limonene based solvents in glass-cloth phenolic composite (GCP) prior to and during subsequent bonding and epoxy adhesive cure.

Wurth, Laura↗
NASA NTRS1972

Diffusion bonding and its application to manufacturing

In its simplest form diffusion bonding is accomplished by placing clean metal surfaces together under a sufficient load and heating. The natural interatomic attractive force between atoms transforms the interface into a natural grain boundary. Therefore, in principle, the properties of the bond area are identical to those of the parent metal. Other advantages of diffusion bonding over conventional methods of bonding include freedom from residual stresses, excessive deformation, foreign metals, or changed crystal structures. Stainless steels, nickel-base superalloys, and aluminum alloys have all been successfully joined. Complex hardware, including integrated flueric devices, jet engine servovalves, and porous woven structures have been fabricated. The processing involved is discussed, along with such theoretical considerations as the role of metal surfaces, the formation of metal contact junctions, and the mechanisms of material transport in diffusion bonding.

Spurgeon, W. M.↗
NASA NTRS1993

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. The Roe approximate Riemann solution scheme or the computationally less expensive Advection Upstream Splitting Method (AUSM) flux-splitting scheme are used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passage and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗
NASA NTRS1994

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. Roe's approximate Riemann solution scheme or the computationally less expensive advection upstream splitting method (AUSM) flux-splitting scheme is used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passages and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗
NASA NTRS2024

Offset-Stoichiometric Reflowable Composite Bonding Method with Adhesive for Mitigating Strict Faying Surface Tolerances

Inherent susceptibility of adhesive bonds to miniscule quantities of contamination can cause undetectable weakened bonds. For this reason, the Federal Aviation Administration (FAA) places strict regulations on adhesively bonded joints in primary aircraft structures. To meet certification requirements aircraft manufactures resort to redundant load paths in the form of fasteners which inherently add weight to the structure and increase manufacturing time. In prior work, a secondary bonding technique called AERoBOND was developed, which utilized off-stoichiometric epoxymatrix resins to facilitate reflow and diffusion of the resin within the joint interface during a secondary bonding/cure process, thus achieving a bond similar to a co-cured joint. However, the AERoBOND process required tight spatial tolerances between the two parts being joined. This study examined the utilization of conventional adhesive with the AERoBOND method to act as a filler in the joint line, effectively reducing the need for tight tolerances on the joining parts and serving as a flexible alternative for existing manufacturing processes. Ultrasonic inspection, optical microscopy, and ASTM International standard tests were performed to analyze the joints for defects and to quantify the mode-I and mode-II interlaminar fracture toughness and short beam strength of the proposed methodology with varying manufacturing parameters. The comprehensive results indicate that the AERoBOND+ method with and without surface preparation performs comparably to co-cured and conventional, adhesively bonded joints when secondary cured in an autoclave with 791 kPa of pressure. As an example, the AERoBOND+ panel without surface preparation bonded with 791 kPa of pressure (referred to as AB+3 throughout paper) had a mode-I fracture toughness (G Ic ) of 0.643 kJ/m 2 and a mode-II fracture toughness (G IIc ) of 4.000 kJ/m 2 in the non-precracked condition and 4.218 kJ/m 2 in the precracked condition at the adhesive-to-prepreg interface. These results were 96%, 142%, and 217%, respectively, of a co-cured baseline panel (referred to as C1 throughout paper).

Composites↗
NASA NTRS1989

Similarity solutions for phase-change problems

A modification of Ivantsov's (1947) similarity solutions is proposed which can describe phase-change processes which are limited by diffusion. The method has application to systems that have n-components and possess cross-diffusion and Soret and Dufour effects, along with convection driven by density discontinuities at the two-phase interface. Local thermal equilibrium is assumed at the interface. It is shown that analytic solutions are possible when the material properties are constant.

Canright, D.↗
NASA NTRS2015

Thermal Pyrolytic Graphite Enhanced Components

A thermally conductive composite material, a thermal transfer device made of the material, and a method for making the material are disclosed. Apertures or depressions are formed in aluminum or aluminum alloy. Plugs are formed of thermal pyrolytic graphite. An amount of silicon sufficient for liquid interface diffusion bonding is applied, for example by vapor deposition or use of aluminum silicon alloy foil. The plugs are inserted in the apertures or depressions. Bonding energy is applied, for example by applying pressure and heat using a hot isostatic press. The thermal pyrolytic graphite, aluminum or aluminum alloy and silicon form a eutectic alloy. As a result, the plugs are bonded into the apertures or depressions. The composite material can be machined to produce finished devices such as the thermal transfer device. Thermally conductive planes of the thermal pyrolytic graphite plugs may be aligned in parallel to present a thermal conduction path.

Hardesty, Robert E.↗
NASA NTRS1993

Ultrasonic inspection of a diffusion-bonded platelet rocket chamber liner

An ultrasonic inspection technique, which applies pulse-echo technology for examining the rocket chamber's liner-to-jacket interface was developed, which makes it possible to characterize disbonded layers in intimate contact. The method was used to inspect a simulated flat-geometry rocket chamber liner structure with debonded liner layers and to validate the use of boron nitride powder as a flaw introduction technique.

Matson, Douglas M.↗
NASA NTRS2000

Thermal Edge-Effects Model for Automated Tape Placement of Thermoplastic Composites

Two-dimensional thermal models for automated tape placement (ATP) of thermoplastic composites neglect the diffusive heat transport that occurs between the newly placed tape and the cool substrate beside it. Such lateral transport can cool the tape edges prematurely and weaken the bond. The three-dimensional, steady state, thermal transport equation is solved by the Green's function method for a tape of finite width being placed on an infinitely wide substrate. The isotherm for the glass transition temperature on the weld interface is used to determine the distance inward from the tape edge that is prematurely cooled, called the cooling incursion Delta a. For the Langley ATP robot, Delta a = 0.4 mm for a unidirectional lay-up of PEEK/carbon fiber composite, and Delta a = 1.2 mm for an isotropic lay-up. A formula for Delta a is developed and applied to a wide range of operating conditions. A surprise finding is that Delta a need not decrease as the Peclet number Pe becomes very large, where Pe is the dimensionless ratio of inertial to diffusive heat transport. Conformable rollers that increase the consolidation length would also increase Delta a, unless other changes are made, such as proportionally increasing the material speed. To compensate for premature edge cooling, the thermal input could be extended past the tape edges by the amount Delta a. This method should help achieve uniform weld strength and crystallinity across the width of the tape.

Costen, Robert C.↗
NASA NTRS1988

Domain decomposition algorithms and computation fluid dynamics

In the past several years, domain decomposition was a very popular topic, partly motivated by the potential of parallelization. While a large body of theory and algorithms were developed for model elliptic problems, they are only recently starting to be tested on realistic applications. The application of some of these methods to two model problems in computational fluid dynamics are investigated. Some examples are two dimensional convection-diffusion problems and the incompressible driven cavity flow problem. The construction and analysis of efficient preconditioners for the interface operator to be used in the iterative solution of the interface solution is described. For the convection-diffusion problems, the effect of the convection term and its discretization on the performance of some of the preconditioners is discussed. For the driven cavity problem, the effectiveness of a class of boundary probe preconditioners is discussed.

Chan, Tony F.↗
NASA NTRS1985

Diffusion welding of MA 6000 and a conventional nickel-base superalloy

A feasibility study of diffusion welding the oxide dispersion strengthened (ODS) alloy MA 6000 to itself and to conventional Ni-base superalloy Udimet 700 was conducted. Butt joints between MA 6000 pieces and lap joints between Udimet 700 and the ODS alloy were produced by hot pressing for 1.25 hr at temperatures ranging from 1000 to 1200 C (1832-2192 F) in vacuum. Following pressing, all weldments were heat treated and machined into mechanical property test specimens. While three different combinations of recrystallized and unrecrystallized MA 6000 butt joints were produced, the unrecrystallized to unrecrystallized joint was most successful as determined by mechanical properties and microstructural examination. Failure to weld the recrystallized material probably related to a lack of adequate deformation at the weld interface. While recrystallized MA 6000 could be diffusion welded to Udimet 700 in places, complete welding over the entire lap joint was not achieved, again due to the lack of sufficient deformation at the faying surfaces. Several methods are proposed to promote the intimate contact necessary for diffusion welding MA 6000 to itself and to superalloys.

Moore, T. J.↗
NASA NTRS1992

Convection of tin in a Bridgman system. I - Flow characterization by effective diffusivity measurements

An electrochemical titration method was used to investigate the dynamic states in a cylindrical layer of convecting tin. The liquid tin was contained in a cell, with curved boundaries made of quartz and flat boundaries made of a solid state electrolyte - yttria-stabilized zirconia (YSZ). The electrolyte acted as a window through which a trace amount of oxygen could be pumped in or out by the application of a constant voltage. The concentration at the YSZ interface was monitored by operating the electrochemical cell in the galvanic mode. Experimentally determined effective diffusivities of oxygen were compared with the molecular diffusivity. Dynamic states in the convective flow were thus inferred. Temperature measurements were simultaneously made in order to identify the onset of oscillations from a steady convective regime. The experiments were conducted for two different aspect ratios for various imposed temperature gradients and two different orientations with respect to gravity. Transcritical states were identified and comparison to two-dimensional numerical models were made.

Sears, B.↗
NASA NTRS1994

Process modelling for materials preparation experiments

The main goals of the research under this grant consist of the development of mathematical tools and measurement techniques for transport properties necessary for high fidelity modelling of crystal growth from the melt and solution. Of the tasks described in detail in the original proposal, two remain to be worked on: development of a spectral code for moving boundary problems, and development of an expedient diffusivity measurement technique for concentrated and supersaturated solutions. We have focused on developing a code to solve for interface shape, heat and species transport during directional solidification. The work involved the computation of heat, mass and momentum transfer during Bridgman-Stockbarger solidification of compound semiconductors. Domain decomposition techniques and preconditioning methods were used in conjunction with Chebyshev spectral methods to accelerate convergence while retaining the high-order spectral accuracy. During the report period we have further improved our experimental setup. These improvements include: temperature control of the measurement cell to 0.1 C between 10 and 60 C; enclosure of the optical measurement path outside the ZYGO interferometer in a metal housing that is temperature controlled to the same temperature setting as the measurement cell; simultaneous dispensing and partial removal of the lower concentration (lighter) solution above the higher concentration (heavier) solution through independently motor-driven syringes; three-fold increase in data resolution by orientation of the interferometer with respect to diffusion direction; and increase of the optical path length in the solution cell to 12 mm.

Rosenberger, Franz↗
NASA NTRS1996

MGLab: An Interactive Multigrid Environment

MGLab is a set of Matlab functions that defines an interactive environment for experimenting with multigrid algorithms. The package solves two-dimensional elliptic partial differential equations discretized using either finite differences or finite volumes, depending on the problem. Built-in problems include the Poisson equation, the Helmholtz equation, a convection-diffusion problem, and a discontinuous coefficient problem. A number of parameters controlling the multigrid V-cycle can be set using a point-and-click mechanism. The menu-based user interface also allows a choice of several Krylov subspace methods, including CG, GMRES(k), and Bi-CGSTAB, which can be used either as stand-alone solvers or as multigrid acceleration schemes. The package exploits Matlab's visualization and sparse matrix features and has been structured to be easily extensible.

Bordner, James↗
NASA NTRS1982

On the cause of the flat-spot phenomenon observed in silicon solar cells at low temperatures and low intensities

A model that explains the flat-spot power loss phenomenon is presented. Evidence suggests that the effect is due to localized metallurgical interactions between the silicon substrate and the contact metallization. These reactions are shown to result in localized regions in which the P-N junction is destroyed and replaced with a metal semiconductor-like interface. The effects of thermal treatment, crystallographic orientation, junction depth, and metallization are shown along with a method of preventing the effect through the suppression of vacancy formation at the free surface of the contact metallization. Data indicating the effectiveness of a TiN diffusion barrier in preventing the effect are also given.

Weizer, V. G.↗