SEARCH · Search NASA

Results for “Approximated interface conditions”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE 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 NTRS1989

Effects of aging condition on the fracture toughness of 2XXX and 7XXX series aluminum alloy composites

Results are presented on the effects of matrix aging condition (i.e., matrix temper) on the fracture toughness of 2XXX and 7XXX Al matrix alloys reinforced with SiC particulates, and the results are compared with the mechanical behavior. Fracture toughness testing was conducted on fatigue precracked bend specimens, and fracture surfaces were examined using SEM. Results revealed dramatic differences in the effect of matrix microstructure on the fracture properties of the two composite series. In the 7XXX material, the toughness values decreased from the underaged (UA) condition to the overaged (OA) condition by approximately 40 percent, while in the 2XXX series composite, the effect of matrix microstructure was marginal. In the 7XXX series composites, a transition in fracture mode from particle cracking (in UA) to matrix and linear-interface failure (in OA) was observed, while the 2XXX series composite failed predominantly by particle cracking.

Manoharan, M.↗

NASA NTRS1999

Numerical Analysis of Temperature Gradients and Interface Shape During Directional Solidification of Al and Al-Cu Alloy Under Microgravity Conditions

Numerical modeling was undertaken to analyze the influence of radial thermal gradient on solid/liquid interface shape and convection patterns during solidification of pure Al and Al-4 wt. % Cu alloy. Steady state calculations were performed for different gravity levels and orientations. Furthermore, transient modeling was undertaken to investigate effect of the solidification velocity. The furnace configuration used in this analysis is the proposed International Space Station Furnace. Results from a thermal model of the furnace core were used as initial boundary conditions for solidification modeling. The Solidification model was adopted from previous work and was based on the finite element code FIDAP. Thermocouple data and quenched interface shape from a pure Al sample flown on the Life Sciences and Microgravity Spacelab (LMS) mission, July 1996 was used for model validation. Good agreement was obtained between the predicted interface shape and that measured from the quenched LMS sample. It was found that the imposed temperature boundary condition must be sufficiently smooth for the heat fluxes in the model to be self-consistent. The model predicted that alloy sample Al-Cu is more sensitive to variations in the gravity level compared to pure Al. For Al-Cu alloy, solute diffusivity is approximately 10(exp 4) times smaller than thermal diffusivity of pure Al. Hence, in a microgravity environment the weak convection has no measurable effect on the heat fluxes, but is still strong enough to affect concentration distribution. Since the alloy melting temperature is determined by concentration the interface shape depends on the level of convection. This was found not to be the case for pure metal.

Bune, Andris V.↗

NASA NTRS1973

Electrical properties of dry rocks

The mechanism by which atmospheric moisture affects the conductivity and dielectric constant of rock specimens was studied in time and frequency domains. It is suggested that adsorbed water molecules alter the surface conductivity in a manner similar to that observed in semiconductors and insulators. Powdered basalts show a low-frequency dispersion produced by the atmospheric moisture remaining in the pore system of the sample in a high vacuum; this effect is attributed to isolated adsorption centers. Simulated lunar permafrost at 100 K and a vacuum of 10 to the -8th power torr together with data on lunar samples contaminated with atmospheric moisture and the dielectric properties of ice at various temperatures indicate that, if permafrost exists in the moon it should present a relaxation peak at approximately 300 Hz; for temperatures up to 263 K it may go up to 20 KHz. It is concluded that in order to have electrical steady state conditions in rock samples it is necessary to have volume charge accumulations at interfaces within the sample and at the electrode sample interface. A method for measuring heterogeneous dielectrics with non-negligible ohmic and dielectric conductivities is proposed and experimentally verified.

Morrison, H.↗

NASA NTRS1985

Ellipsometric and optical study of some uncommon insulator films on 3-5 semiconductors

Optical properties of three types of insulating films that show promise in potential applications in the 3-4 semiconductor technology were evaluated, namely a-C:H, BN and CaF2. The plasma deposited a-C:H shows an amorphous behavior with optical energy gaps of approximately 2 to 2.4 eV. These a-C:H films have higher density and/or hardness, higher refractive index and lower optical energy gaps with increasing energy of the particles in the plasma, while the density of states remains unchanged. These results are in agreement, and give a fine-tuned positive confirmation to an existing conjecture on the nature of a-C:H films (1). Ion beam deposited BN films show amorphous behavior with energy gap of 5 eV. These films are nonstoichiometric (B/N approximately 2) and have refractive index, density and/or hardness which are dependent on the deposition conditions. The epitaxially grown CaF2 on GaAs films have optical parameters equal to bulk, but evidence of damage was found in the GaAs at the interface.

Alterovitz, S. A.↗

NASA NTRS1992

Segmented multigrid domain decomposition procedure for incompressible viscous flows

In this paper, the Navier-Stokes (NS) equations are approximated with a reduced NS system, that represents the lowest-order terms in an asymptotic Re expansion. This system allows for simplified boundary conditions, more generality in the location of the outflow boundary, and insures mass conservation in all subdomain grid interfaces, and also at the outflow boundary. The higher-order NS diffusion terms are included through a deferred corrector, in chosen subdomains, when required.

Srinivasan, Kumar↗

NASA NTRS1993

Three-dimensional analysis of anisotropic spatially reinforced structures

The material-adaptive three-dimensional analysis of inhomogeneous structures based on the meso-volume concept and application of deficient spline functions for displacement approximations is proposed. The general methodology is demonstrated on the example of a brick-type mosaic parallelepiped arbitrarily composed of anisotropic meso-volumes. A partition of each meso-volume into sub-elements, application of deficient spline functions for a local approximation of displacements and, finally, the use of the variational principle allows one to obtain displacements, strains, and stresses at anypoint within the structural part. All of the necessary external and internal boundary conditions (including the conditions of continuity of transverse stresses at interfaces between adjacent meso-volumes) can be satisfied with requisite accuracy by increasing the density of the sub-element mesh. The application of the methodology to textile composite materials is described. Several numerical examples for woven and braided rectangular composite plates and stiffened panels under transverse bending are considered. Some typical effects of stress concentrations due to the material inhomogeneities are demonstrated.

Bogdanovich, Alexander E.↗

NASA NTRS1993

Accurate upwind methods for the Euler equations

A new class of piecewise linear methods for the numerical solution of the one-dimensional Euler equations of gas dynamics is presented. These methods are uniformly second-order accurate, and can be considered as extensions of Godunov's scheme. With an appropriate definition of monotonicity preservation for the case of linear convection, it can be shown that they preserve monotonicity. Similar to Van Leer's MUSCL scheme, they consist of two key steps: a reconstruction step followed by an upwind step. For the reconstruction step, a monotonicity constraint that preserves uniform second-order accuracy is introduced. Computational efficiency is enhanced by devising a criterion that detects the 'smooth' part of the data where the constraint is redundant. The concept and coding of the constraint are simplified by the use of the median function. A slope steepening technique, which has no effect at smooth regions and can resolve a contact discontinuity in four cells, is described. As for the upwind step, existing and new methods are applied in a manner slightly different from those in the literature. These methods are derived by approximating the Euler equations via linearization and diagonalization. At a 'smooth' interface, Harten, Lax, and Van Leer's one intermediate state model is employed. A modification for this model that can resolve contact discontinuities is presented. Near a discontinuity, either this modified model or a more accurate one, namely, Roe's flux-difference splitting. is used. The current presentation of Roe's method, via the conceptually simple flux-vector splitting, not only establishes a connection between the two splittings, but also leads to an admissibility correction with no conditional statement, and an efficient approximation to Osher's approximate Riemann solver. These reconstruction and upwind steps result in schemes that are uniformly second-order accurate and economical at smooth regions, and yield high resolution at discontinuities.

Huynh, Hung T.↗

NASA NTRS1998

Performance of the Spacecraft Propulsion Research Facility During Altitude Firing Tests of the Delta 3 Upper Stage

The Spacecraft Propulsion Research Facility at the NASA Lewis Research Center's Plum Brook Station was reactivated in order to conduct flight simulation ground tests of the Delta 3 cryogenic upper stage. The tests were a cooperative effort between The Boeing Company, Pratt and Whitney, and NASA. They included demonstration of tanking and detanking of liquid hydrogen, liquid oxygen and helium pressurant gas as well as 12 engine firings simulating first, second, and third burns at altitude conditions. A key to the success of these tests was the performance of the primary facility systems and their interfaces with the vehicle. These systems included the structural support of the vehicle, propellant supplies, data acquisition, facility control systems, and the altitude exhaust system. While the facility connections to the vehicle umbilical panel simulated the performance of the launch pad systems, additional purge and electrical connections were also required which were unique to ground testing of the vehicle. The altitude exhaust system permitted an approximate simulation of the boost-phase pressure profile by rapidly pumping the test chamber from 13 psia to 0.5 psia as well as maintaining altitude conditions during extended steady-state firings. The performance of the steam driven ejector exhaust system has been correlated with variations in cooling water temperature during these tests. This correlation and comparisons to limited data available from Centaur tests conducted in the facility from 1969-1971 provided insight into optimizing the operation of the exhaust system for future tests. Overall, the facility proved to be robust and flexible for vehicle space simulation engine firings and enabled all test objectives to be successfully completed within the planned schedule.

Meyer, Michael L.↗

NASA NTRS1977

Velocity and temperature distributions of coal-slag layers on magnetohydrodynamic generators walls

Approximate analytical expressions are derived for the velocity and temperature distributions in steady state coal slag deposits flowing over MHD generator walls. Effects of slag condensation and Joule heating are included in the analysis. The transport conditions and the slag temperature at the slag-gas interface are taken to be known parameters in the formulation. They are assumed to have been predetermined either experimentally or from the slag properties and the gas dynamic calculations of the free stream flow. The analysis assumes a power law velocity profile for the slag and accounts for the coupling between the energy and momentum conservation equations. Comparisons are made with the more exact numerical solutions to verify the accuracy of the results.

Pian, C. C. P.↗

NASA NTRS2000

Modeling of a Non-Dilute Alloy Solidification Under Terrestrial and Microgravity Conditions

A general finite element model of non-dilute alloy solidification was used to simulate growth of HgCdTe from the melt. 2-D and 3-D models with and without Oberbeck-Boussinesq (O-B) approximation were used. Verification of the 3-D model was based on experimental data. The location and the shape of the interface were determined using melting temperatures obtained from the HgCdTe liquidus curve. Double-diffusive convection in the melt is caused by radial temperature gradients and by rejection at the solid-melt interface of the heavier HgTe-rich solute. Because of these complicating factors, dimensional rather than non-dimensional modeling was performed. Models based on complete set of Navier-Stokes equations for compressible media are superior to the O-B approximation in microgravity, because they provide more accurate estimation of convection.

Bune, Andris V.↗

NASA NTRS2023

Assessment of the Fluid Dynamics Boundary Condition in Ablating or Blowing Flows

Improved models of ablative thermal protection systems have enabled the treatment of materials and fluid behavior in a coupled manner. This paper reports a new approach to modeling the interface between fluid and material, with attention to the conservation of species mass flux and energy on the fluid side of the interface. The general equation is presented and is shown to recover the traditional uncoupled fluid/materials response interface. Including the chemical reaction terms on the CFD side of the interface makes the heat flux exchange independent of the thermodynamic reference state and, therefore, a measurable quantity. Doing so allows the material response solver to take as input the surface heat flux rather than a film coefficient. Removing the film coefficient approximation enables more direct solution of vehicle thermal response but requires consistency in the wall state. The mixing of the shock layer and pyrolysis gas is then computed with finite rate chemistry within the fluid solver. The boundary conditions described have been implemented in the DPLR v4.05.1 code. Char removal is captured using finite rate chemistry in DPLR’s gas surface interaction module. Aspects of coupling these solutions to material response are discussed.

Ablation↗

NASA NTRS2023

Assessment of the Fluid Dynamics Boundary Condition in Ablating or Blowing Flows

Ablation↗

NASA NTRS2022

A Novel High-Performance Mission-Enabling Multi-Purpose Radioisotope Heat Source

Recent studies indicate science mission concepts targeting access to the sub-surface oceans of icy moons require ice-penetrating cryobots powered by advanced Radioisotope Power Systems (RPS). These systems would deliver waste heat for ice-melting in the range of 10 kW. Minimizing the transit time through the kilometers-thick ice shells to just a few years requires these RPS to utilize heat sources having a higher thermal energy volumetric density than the existing flight-qualified General-Purpose Heat Source (GPHS). A Compact Heat Source (CPHS)has been conceptualized in which the graphite impact shells(GIS) of the existing GPHS are rearranged in a hexagonal aeroshell containing seven GIS per module, as opposed to the standard two per module; offering a thermal energy density of 0.57 W/cm3versus 0.29 W/cm3 offered by the GPHS simply from the repackaging of Technology Readiness Level (TRL)9 subassemblies. Preliminary thermal modeling of the CPHS integrated into a notional radioisotope thermoelectric generator(RTG) structure further suggests that centerline temperatures are well within allowable limits during nominal operation. Given the need for the CPHS for a subset of missions, it is worth exploring the applicability of the CPHS for more general RTG purposes. We discuss herein how the CPHS may be implemented with either heritage or in-development thermoelectric converter technologies into a Next-Generation RTG concept. Due to a higher energy density, the legacy heat rejection fin arrangement must be modified to permit a sufficiently low cold-side temperature. Preliminary finite element analysis suggests fin-root temperatures can be kept as low as 520 K while allowing the generator to fit within the usable dimensions of currently available United States Department of Energy shipping containers. Such temperatures would certainly be compatible with the use of high temperature thermoelectric converter technologies.. A prime candidate is the heritage silicon-germanium (SiGe) unicouple, whose design could be adapted by approximately halving the leg-length, but without changes in hot and cold junction interfaces, which are features critical to the proven performance and reliability of these devices. The estimated Beginning of Life power for a SiGe-based CPHS-RTG using 12 CPHS for a thermal inventory of 10.5 kW is greater than 600 W under deep space operating conditions. Using higher performance segmented couples currently in development that are based on skutterudite,La3−xTe4and 14-1-11 Zintl thermoelectric materials in lieu of the SiGe unicouples would increase the power level to more than1 kW. The high specific power (We/kg) attribute of CPHS-RTGs found in this study could potentially enable Radioisotope Electric Propulsion (REP) mission concepts. Past NASA REP mission concept studies identified specific power needs in excess of 6to 8 We/kg. Based on a GPHS-RTG-like system configuration, we show that at fin root temperatures between 530 K and 570K (deep space environment), specific powers exceeding 10 We/kg are achievable using high performance segmented thermoelectric converters. The compact sizing and power density of the CPHS-RTG would constitute a significant step upgrade in specific power when compared to heritage GPHS-RTG (approximately5.1 We/kg) and off-the-shelf Multi-Mission RTG (approximately 2.6 We/kg).

Nesmith, Bill J.↗

NASA NTRS1993

Thermomechanical and isothermal fatigue behavior of a (90)sub 8 titanium matrix composite

An experimental investigation was conducted to analyze the fatigue damage mechanisms operative in a 35 fiber volume percent (90 deg) titanium matrix composite (TMC) under 427 C isothermal and thermomechanical loading conditions. The thermomechanical fatigue (TMF) tests were performed with a temperature cycle from 200 to 427 C under closely controlled conditions to define both the deformation and fatigue life behavior. Degradation of the TMC's isothermal elastic moduli and coefficient of thermal expansion were monitored throughout the TMF tests. The results indicated TMF life trends which contrasted those established in a comparable (0 deg) system, as TMF loading of the (90 deg) TMC was not found to be 'life-limiting' in comparison to maximum temperature isothermal conditions. In-phase lives were very similar to those established under 427 C isothermal conditions. High stress isothermal fatigue at 427 C produced increased strain ratchetting and stiffness degradation relative to TMF conditions. Out-of-phase loadings produced TMF lives approximately an order of magnitude greater than the lives determined under isothermal and in-phase conditions. Extensive fractography and metallography were also performed. Two key issues were identified and appeared to dominate the fatigue damage and life of the (90 deg) TMC, namely, the weak fiber/matrix interface and environmental attack of the fiber/matrix interface via the exposed (90 deg) fiber ends.

Castelli, Michael G.↗

NASA NTRS1983

An inviscid three-dimensional analysis of the Space Shuttle main engine hot-gas manifold

A numerical study using an inviscid three-dimensional Lagrangian fluid dynamics code has been conducted as a part of an overall effort to understand the flow behavior in the SSME fuel side hot-gas manifold. The model simulates flow from the high-pressure fuel turbine exit through the transfer ducts, including the effects of swirl, inlet flow symmetry, and presence of straightening vanes and struts; a separate, more-detailed effort is in progress that includes viscosity and turbulence effects. The simplified model presented is divided into two parts, the first includes the 180-degree turnaround duct downstream of the turbine exit and the spherical fuel bowl section, while the second models the three transfer ducts. The two parts of the model are coupled together with the interface conditions being updated through iteration. Results indicate that a transverse pressure differential of 165 psi would be imposed on the turbine exit and that unstable flow separation occurs around the vanes, struts, and within the transfer ducts. The three transfer ducts show a mass flux split of approximately 41, 21, and 38 percent. Results to date are encouraging that certain flow characteristics can be usefuly represented using a relatively coarse grid inviscid code.

Liang, P. Y.↗

NASA NTRS1990

A transient response analysis of the space shuttle vehicle during liftoff

A proposed transient response method is formulated for the liftoff analysis of the space shuttle vehicles. It uses a power series approximation with unknown coefficients for the interface forces between the space shuttle and mobile launch platform. This allows the equation of motion of the two structures to be solved separately with the unknown coefficients at the end of each step. These coefficients are obtained by enforcing the interface compatibility conditions between the two structures. Once the unknown coefficients are determined, the total response is computed for that time step. The method is validated by a numerical example of a cantilevered beam and by the liftoff analysis of the space shuttle vehicles. The proposed method is compared to an iterative transient response analysis method used by Martin Marietta for their space shuttle liftoff analysis. It is shown that the proposed method uses less computer time than the iterative method and does not require as small a time step for integration. The space shuttle vehicle model is reduced using two different types of component mode synthesis (CMS) methods, the Lanczos method and the Craig and Bampton CMS method. By varying the cutoff frequency in the Craig and Bampton method it was shown that the space shuttle interface loads can be computed with reasonable accuracy. Both the Lanczos CMS method and Craig and Bampton CMS method give similar results. A substantial amount of computer time is saved using the Lanczos CMS method over that of the Craig and Bampton method. However, when trying to compute a large number of Lanczos vectors, input/output computer time increased and increased the overall computer time. The application of several liftoff release mechanisms that can be adapted to the proposed method are discussed.

Brunty, J. A.↗

NASA NTRS2019

Electromagnetic Scattering by Discrete Random Media. II: The Coherent Field

The computation of the coherent field in the case of a plane electromagnetic wave obliquely incident on a discrete random layer with non-scattering boundaries is addressed. For dense media, the analysis is based on a special-form solution for the conditional configuration-averaged exciting field coefficients, and is restricted to the computation of the so-called zeroth-order fields without a special treatment of the boundary regions. In this setting, we calculate the coherent fields reflected and transmitted by the layer, and the coherent field inside the layer. We found that these fields are analytically equivalent to plane electromagnetic waves, and investigated the fulfillment of the boundary conditions for the electric fields at the layer interfaces. The results are then particularized to the cases of normal incidence and a semi-infinite discrete random medium. For sparsely distributed particles, we present a self-consistent derivation of the coherent field and discuss the Twersky and Foldy approximations.

Adrian Doicu↗