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 91 records · Page 5

Improved Pyrotechnic Shock Scaling Method Based on Shock Response Spectrum and Statistical Energy Analysis

It has been more than forty years since the pyrotechnic shock scaling method was introduced. The scaling method estimates the attenuation of the Shock Response Spectrum (SRS) based on the distance from the source, structural configurations, types of structural joints and interfaces, and intervening structure [1]. The method has been successfully used in the spacecraft community countless times and still is frequently used to develop pyrotechnic shock requirements at various levels of assembly of a spacecraft. However, since the method was derived empirically from a limited set of shock test data [2], the aerospace community has been looking for an alternative approach. In this paper, a computational scaling method based on Statistical Energy Analysis (SEA) is re-visited. SEA is traditionally used as a method for investigating the diffusion of acoustic and vibratory energies of a system at steady-state [3]. Previously, an approach using the SRS as an acceleration constraint condition to the SEA model was introduced to estimate the attenuation from a shock source [4,5]. In the current investigation, additional examples are provided to further validate the approach. In addition, responses in the time domain are produced by the Local Modal Phase Reconstruction (LMPR) approach which was recently introduced to the community [6].

Fernandez, Juan P.↗

Axial and centrifugal pump meanline performance analysis

A meanline pump flow modeling method has been developed to provide a fast capability for modeling pumps of cryogenic rocket engines. Based on this method, a meanline pump flow code (PUMPA) has been written that can predict the performance of pumps at off-design operating conditions, given the loss of the diffusion system at the design point. The design point rotor efficiency is obtained from empirically derived correlations of loss to rotor specific speed. The rapid input setup and computer run time for the meanline pump flow code makes it an effective analysis and conceptual design tool. The map generation capabilities of the PUMPA code provide the information needed for interfacing with a rocket engine system modeling code.

Veres, Joseph P.↗

Semi-Analytic Reconstruction of Flux in Finite Volume Formulations

Semi-analytic reconstruction uses the analytic solution to a second-order, steady, ordinary differential equation (ODE) to simultaneously evaluate the convective and diffusive flux at all interfaces of a finite volume formulation. The second-order ODE is itself a linearized approximation to the governing first- and second- order partial differential equation conservation laws. Thus, semi-analytic reconstruction defines a family of formulations for finite volume interface fluxes using analytic solutions to approximating equations. Limiters are not applied in a conventional sense; rather, diffusivity is adjusted in the vicinity of changes in sign of eigenvalues in order to achieve a sufficiently small cell Reynolds number in the analytic formulation across critical points. Several approaches for application of semi-analytic reconstruction for the solution of one-dimensional scalar equations are introduced. Results are compared with exact analytic solutions to Burger s Equation as well as a conventional, upwind discretization using Roe s method. One approach, the end-point wave speed (EPWS) approximation, is further developed for more complex applications. One-dimensional vector equations are tested on a quasi one-dimensional nozzle application. The EPWS algorithm has a more compact difference stencil than Roe s algorithm but reconstruction time is approximately a factor of four larger than for Roe. Though both are second-order accurate schemes, Roe s method approaches a grid converged solution with fewer grid points. Reconstruction of flux in the context of multi-dimensional, vector conservation laws including effects of thermochemical nonequilibrium in the Navier-Stokes equations is developed.

Gnoffo, Peter A.↗

Diagnostic study of BSF silicon solar cells

Solar cells equipped with back surface fields (BSF) were fabricated by means of an aluminum alloy, boron diffusion, and boron ion implantation. The importance of initial thickness, resistivity, orientation, and the crystal growth method were examined. Aluminum in paste form was screen printed on the cell, followed by alloy formation at 800 C for a minute. Application of the BSFs resulted in open circuit voltage improvements dependent on the substrate resistivity, ranging from 10 mV at 0.3 ohm-cm to 50 mV at 15 ohm-cm. The gains in performance were attributed to a longer penetration depth, a heightened concentration profile at the interface, and higher concentrations in the p+ layer.

Yoo, H.↗

A Gas-Kinetic Method for Hyperbolic-Elliptic Equations and Its Application in Two-Phase Fluid Flow

A gas-kinetic method for the hyperbolic-elliptic equations is presented in this paper. In the mixed type system, the co-existence and the phase transition between liquid and gas are described by the van der Waals-type equation of state (EOS). Due to the unstable mechanism for a fluid in the elliptic region, interface between the liquid and gas can be kept sharp through the condensation and evaporation process to remove the "averaged" numerical fluid away from the elliptic region, and the interface thickness depends on the numerical diffusion and stiffness of the phase change. A few examples are presented in this paper for both phase transition and multifluid interface problems.

Xu, Kun↗

Diffusion analysis for two-phase metal-matrix composite

Diffusion controlled filament matrix interaction in a metal matrix composite, where the filaments and matrix comprise a two phase binary alloy system, was mathematically modeled. The problem of a diffusion controlled, two phase moving interface by means of a one dimensional, variable grid, finite difference technique was analyzed. Concentration dependent diffusion coefficients and equilibrium solubility limits were used, and the change in filament diameter and compositional changes in the matrix were calculated as a function of exposure time at elevated temperatures. With the tungsten nickel (W-Ni) system as a model composite system, unidirectional composites containing from 0.06 to 0.44 initial filament volume fraction were modeled. Compositional changes in the matrix were calculated by superposition of the contributions from neighboring filaments. Alternate methods for determining compositional changes between first and second nearest neighbor filaments were also considered. The results show the relative importance of filament volume fraction, filament diameter, exposure temperature, and exposure time as they affect the rate and extent of filament matrix interaction.

Tenney, D. R.↗

Mean Line Pump Flow Model in Rocket Engine System Simulation

A mean line pump flow modeling method has been developed to provide a fast capability for modeling turbopumps of rocket engines. Based on this method, a mean line pump flow code PUMPA has been written that can predict the performance of pumps at off-design operating conditions, given the loss of the diffusion system at the design point. The pump code can model axial flow inducers, mixed-flow and centrifugal pumps. The code can model multistage pumps in series. The code features rapid input setup and computer run time, and is an effective analysis and conceptual design tool. The map generation capability of the code provides the map information needed for interfacing with a rocket engine system modeling code. The off-design and multistage modeling capabilities of the code permit parametric design space exploration of candidate pump configurations and provide pump performance data for engine system evaluation. The PUMPA code has been integrated with the Numerical Propulsion System Simulation (NPSS) code and an expander rocket engine system has been simulated. The mean line pump flow code runs as an integral part of the NPSS rocket engine system simulation and provides key pump performance information directly to the system model at all operating conditions.

Veres, Joseph P.↗

Convection and segregation in directional solidification of dilute and non-dilute binary alloys - Effects of ampoule and furnace design

A Petrov-Galerkin/finite-element method is used to analyze the effect of furnace configuration and ampoule design on the temperature field, the convection in the melt, the shape of the melt-solid interface, and the segregation of solute in the crystal, in the directional solidification of several dilute and nondilute binary semiconductor alloys. The vertical Bridgman-Stockbarger system leads to a two-cell flow structure, and the constant gradient furnace has only a single cell near the interface. Radial temperature gradients interact with the solute field to cause the previously predicted sideways diffusive instability, and addition of a solute that is less dense than the bulk melt and that is incorporated upon solidification decreases the intensity of the flow near the interface by increasing the melt density there. The present results have application to understanding of the effects on solute segregation of microgavity solidification and of applied magnetic fields.

Adornato, Peter M.↗

The study of flow pattern and phase-change problem in die casting process

The flow pattern and solidification phenomena in die casting process have been investigated in the first phase study. The flow pattern in filling process is predicted by using a VOF (volume of fluid) method. A good agreement with experimental observation is obtained for filling the water into a die cavity with different gate geometry and with an obstacle in the cavity. An enthalpy method has been applied to solve the solidification problem. By treating the latent heat implicitly into the enthalpy instead of explicitly into the source term, the CPU time can be reduced at least 20 times. The effect of material properties on solidification fronts is tested. It concludes that the dependence of properties on temperature is significant. The influence of the natural convection over the diffusion has also been studied. The result shows that the liquid metal solidification phenomena is diffusion dominant, and the natural convection can affect the shape of the interface. In the second phase study, the filling and solidification processes will be considered simultaneously.

Wang, T. S.↗

The Effects of Ground and Space Processing on the Properties of Organic, Polymeric, and Colloidal Materials

In recent years, a great deal of interest has been directed toward the use of organic materials in the development of high-efficiency optoelectronic and phototonic devices. There is a myriad of possibilities among organic materials which allow flexibility in the design of unique structures with a variety of functional groups. The use of nonlinear optical (NLO) organic materials as thin film wave-guides allows full exploitation of their desirable qualifies by permitting long interaction lengths and large susceptibilities allowing modest power input. There are several methods in use to prepare thin films such as Langmuir-Blodgett (LB) and self-assembly techniques, vapor deposition, growth from sheared solution or melt, and melt growth between glass plates. Organic-based materials have many features that make them desirable for use in optical devices, such as high second-and third-order nonlinearity, flexibility of molecular design, and damage resistance to optical radiation. However, their use in devices has been hindered by processing difficulties for crystals and thin films. We discuss the potential role of microgravity processing of a few organic and polymeric materials. It is of interest to note how materials with second-and third-order NLO behavior may be improved in a diffusion-limited environment and ways in which convection may be detrimental to these materials. We focus our discussion on third-order materials for all-optical switching, and second-order materials for frequency conversion and electrooptics. The goal of minimizing optical loss obviously depends on processing methods. For solution-based processes, such as solution crystal growth and solution photopolymerization, it is well known that thermal and solutal density gradients can initiate buoyancy-driven convection. Resultant fluid flows can affect transport of material to and from growth interfaces and become manifest in the morphology and homogeneity of the growing film or crystal. Likewise, buoyancy-driven convection can hinder production of defect-free, high-quality crystals or films during crystal and film growth by vapor deposition.

Frazier, Donald O.↗

Increasing Durability of Flame-Sprayed Strain Gauges

Thermally sprayed dielectric ceramic coatings are the primary means of attaching strain and temperature gauges to hot-section rotating parts of turbine engines. As hot-section temperatures increase, lifetimes of installed gauges decrease, and seldom exceed one hour above 2,000 F (approx.1,100 C). Advanced engine components are expected to operate at temperatures approaching 2,200 F (approx.1,200 C), and the required high-temperature lifetime is 10 hours minimum. Typically, to enable a ceramic coating to adhere to the smooth surface of an engine component, a thermally sprayed NiCrAlY or NiCoCrAlY bond coat is applied to the smooth surface, thereby providing a textured surface to which the ceramic coat can adhere. The main failure mechanism of this system is decohesion and/or delamination at the interface between the ceramic top coat and the bond coat, caused by oxidation of the bond coat and stresses from the mismatch between the coefficients of thermal expansion of the ceramic top coat and the metallic bond coat. The approach taken to increase the high-temperature lifetime of a gauge attached to an engine component by the method described above involves (1) selective oxidation of the bond coat by means of a heat treatment in reduced oxygen partial pressure followed by (2) the application of a noble-metal diffusion barrier. In experiments to test this approach, heat treatments of NiCoCrAlY bond coats were carried out in a tube furnace in which, in each case, the temperature was alternately (1) increased at a rate of 3 C per minute and (2) held steady for one hour until the desired temperature was reached. The tube furnace was continuously purged with dry nitrogen gas. A final heat-treatment temperature range of 1,600 to 1,800 F (871 to 982 C) proved most beneficial.

Gregory, Otto J.↗

Cooling rates of group IVA iron meteorites

Cooling rates of six group IVA iron meteorites were estimated by a taenite central Ni concentration-taenite half-width method. Calculated cooling rates range from 13 to 25 C/Myr, with an average of 20 C/Myr. No correlation between cooling rate and bulk Ni content is observed, and the data appear to be consistent with a uniform cooling rate as expected from an igneous core origin. This result differs from previous studies reporting a wide range in cooling rates that were strongly correlated with bulk Ni content. The differences result mainly from differences in the phase diagram and the selected diffusion coefficients. Cooling rates inferred from taenite Ni concentrations at the interface with kamacite are consistent with those based on taenite central Ni content.

Willis, J.↗

An Interactive, Design and Educational Tool for Supersonic External-Compression Inlets

A workstation-based interactive design tool called VU-INLET was developed for the inviscid flow in rectangular, supersonic, external-compression inlets. VU-INLET solves for the flow conditions from free stream, through the supersonic compression ramps, across the terminal normal shock region and the subsonic diffuser to the engine face. It calculates the shock locations, the capture streamtube, and the additive drag of the inlet. The inlet geometry can be modified using a graphical user interface and the new flow conditions recalculated interactively. Free stream conditions and engine airflow can also be interactively varied and off-design performance evaluated. Flow results from VU-INLET can be saved to a file for a permanent record, and a series of help screens make the simulator easy to learn and use. This paper will detail the underlying assumptions of the models and the numerical methods used in the simulator.

Benson, Thomas J.↗

Advanced methods for preparation and characterization of infrared detector materials

Crystals were prepared by the Bridgman-Stockbarger method with a wide range of crystal growth rates and temperature gradients adequate to prevent constitutional supercooling under diffusion-limited, steady-state, growth conditions. The longitudinal compositional gradients for different growth conditions and alloy compositions were calculated and compared with experimental data to develop a quantitative model of solute redistribution during the crystal growth of the alloys. Measurements were performed to ascertain the effect of growth conditions on radial compositional gradients. The pseudobinary HgTe-CdTe constitutional phase diagram was determined by precision differential-thermal-analysis measurements and used to calculate the segregation coefficient of Cd as a function of x and interface temperature. Experiments were conducted to determine the ternary phase equilibria in selected regions of the Hg-Cd-Te constitutional phase diagram. Electron and hole mobilities as functions of temperature were analyzed to establish charge-carrier scattering probabilities. Computer algorithms specific to Hg(1-x)CdxTe were developed for calculations of the charge-carrier concentration, charge-carrier mobilities, Hall coefficient, and Dermi Fermi energy as functions of x, temperature, ionized donor and acceptor concentrations, and neutral defect concentrations.

Broerman, J. G.↗

Advanced methods for preparation and characterization of infrared detector materials

Mercury cadmium telluride crystals were prepared by the Bridgman method with a wide range of crystal growth rates and temperature gradients adequate to prevent constitutional supercooling under diffusion-limited, steady state, growth conditions. The longitudinal compositional gradients for different growth conditions and alloy compositions were calculated and compared with experimental data to develop a quantitative model of the crystal growth kinetics for the Hg(i-x)CdxTe alloys, and measurements were performed to ascertain the effect of growth conditions on radial compositional gradients. The pseudobinary HgTe-CdTe constitutional phase diagram was determined by precision differential thermal analysis measurements and used to calculate the segregation coefficient of Cd as a function of x and interface temperature. Computer algorithms specific to Hg(1-x)CdxTe were developed for calculations of the charge carrier concentrations, charge carrier mobilities, Hall coefficient, optical absorptance, and Fermi energy as functions of x, temperature, ionized donor and acceptor concentrations, and neutral defect concentrations.

Lehoczky, S. L.↗

A quantitative study of factors influencing lamellar eutectic morphology during solidification

The factors that influence the shape of the solid-liquid interface of a lamellar binary eutectic alloy are evaluated. Alloys of carbon tetrabromide and hexachloroethane which serve as a transparent analogue of lamellar metallic eutectics are used. The observed interface shapes are analyzed by computer-aided methods. The solid-liquid interfacial free energies of each of the individual phases comprising the eutectic system are measured as a function of composition using a 'grain boundary groove' technique. The solid-liquid interfacial free energy of the two phases are evaluated directly from the eutectic interface. The phase diagram for the system, the heat of fusion as a function of composition, and the density as a function of composition are measured. The shape of the eutectic interface is controlled mainly by the solid-liquid and solid-solid interfacial free energy relationships at the interface and by the temperature gradient present, rather than by interlamellar diffusion in the liquid at the interface, over the range of growth rates studied.

Kaukler, W. F. S.↗

An Investigation of Wave Propagations in Discontinuous Galerkin Method

Analysis of the discontinuous Galerkin method has been carried out for one- and two-dimensional system of hyperbolic equations. Analytical, as well as numerical, properties of wave propagation in a DGM scheme are derived and verified with direct numerical simulations. In addition to a systematic examination of the dissipation and dispersion errors, behaviours of a DG scheme at an interface of two different grid topologies are also studied. Under the same framework, a quantitative discrete analysis of various artificial boundary conditions is also conducted. Progress has been made in numerical boundary condition treatment that is closely related to the application of DGM in aeroacoustics problems. Finally, Fourier analysis of DGM for the Convective diffusion equation has also be studied in connection with the application of DG schemes for the Navier-Stokes equations. This research has resulted in five(5) publications, plus one additional manuscript in preparation, four(4) conference presentations, and three(3) departmental seminars, as summarized in part II. Abstracts of papers are given in part 111 of this report.

Hu, Fang Q.↗

Centrifugal and Axial Pump Design and Off-Design Performance Prediction

A meanline pump-flow modeling method has been developed to provide a fast capability for modeling pumps of cryogenic rocket engines. Based on this method, a meanline pump-flow code PUMPA was written that can predict the performance of pumps at off-design operating conditions, given the loss of the diffusion system at the design point. The design-point rotor efficiency and slip factors are obtained from empirical correlations to rotor-specific speed and geometry. The pump code can model axial, inducer, mixed-flow, and centrifugal pumps and can model multistage pumps in series. The rapid input setup and computer run time for this meanline pump flow code make it an effective analysis and conceptual design tool. The map-generation capabilities of the code provide the information needed for interfacing with a rocket engine system modeling code. The off-design and multistage modeling capabilities of PUMPA permit the user to do parametric design space exploration of candidate pump configurations and to provide head-flow maps for engine system evaluation.

Veres, Joseph P.↗