Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical performance”

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 487 records · Page 27

Measurement of Temperature Fluctuations and Microscopic Growth Rates in a Silicon Floating Zone on TEXUS36

Several earlier (micro)g experiments have shown that time-dependent thermocapillary (Marangoni) convection is the major cause for the formation of dopant striations in floating-zone grown semiconductor crystals, at least in small-scale systems not employing RF heating. To quantify this correlation, a silicon floating-zone experiment was performed on the TEXUS36 flight (February 7, 1 998) in the monoellipsoid mirror furnace TEM02-ELLI. During the experiment, temperature fluctuations in the silicon melt zone and the microscopic growth rate were simultaneously measured. Temperature fluctuations of 0.5 C - 0.7 C with main frequencies between 0.1 Hz and 0.3 Hz were detectable. The microscopic growth rate fluctuated considerably around the average growth rate of 1 mm/min: rates from 4mm/min to negative values (backmelting) were observed. Dopant striations are clearly visible in the Sb-doped crystal. The frequencies associated with the dopant inhomogeneities correspond quite well with those of the temperature fluctuations and microscopic growth rates. 3D numerical simulations were performed to predict the optimum position of the temperature sensor, to evaluate characteristic temperature amplitudes and frequencies, and to give insight into the instability mechanisms of Maran-goni convection in this configuration. The simulations were in good agreement with the experimental values, showing temperature fluctuations with frequencies f? 0.25 Hz and amplitudes up to 1.8 C at a position equivalent to that of the sensor tip in the experiment.

Croell, Arne↗

Numerical analysis of natural convection in liquid droplets by phase change

A numerical analysis is performed on thermocapillary buoyancy convection induced by phase change in a liquid droplet. A finite-difference code is developed using an alternating-direction implicit (ADI) scheme. The intercoupling relation between thermocapillary force, buoyancy force, fluid property, heat transfer, and phase change, along with their effects on the induced flow patterns, are disclosed. The flow is classified into three types: thermocapillary, buoyancy, and combined convection. Among the three mechanisms, the combined convection simulates the experimental observations quite well, and the basic mechanism of the observed convection inside evaporating sessile drops is thus identified. It is disclosed that evaporation initiates unstable convection, while condensation always brings about a stable density distribution which eventually damps out all fluid disturbances. Another numerical model is presented to study the effect of boundary recession due to evaporation, and the 'peeling-off' effect (the removal of the surface layer of fluid by evaporation) is shown to be relevant.

Duh, J. C.↗

Runge-Kutta methods combined with compact difference schemes for the unsteady Euler equations

Recent development using compact difference schemes to solve the Navier-Stokes equations show spectral-like accuracy. A study was made of the numerical characteristics of various combinations of the Runge-Kutta (RK) methods and compact difference schemes to calculate the unsteady Euler equations. The accuracy of finite difference schemes is assessed based on the evaluations of dissipative error. The objectives are reducing the numerical damping and, at the same time, preserving numerical stability. While this approach has tremendous success solving steady flows, numerical characteristics of unsteady calculations remain largely unclear. For unsteady flows, in addition to the dissipative errors, phase velocity and harmonic content of the numerical results are of concern. As a result of the discretization procedure, the simulated unsteady flow motions actually propagate in a dispersive numerical medium. Consequently, the dispersion characteristics of the numerical schemes which relate the phase velocity and wave number may greatly impact the numerical accuracy. The aim is to assess the numerical accuracy of the simulated results. To this end, the Fourier analysis is to provide the dispersive correlations of various numerical schemes. First, a detailed investigation of the existing RK methods is carried out. A generalized form of an N-step RK method is derived. With this generalized form, the criteria are derived for the three and four-step RK methods to be third and fourth-order time accurate for the non-linear equations, e.g., flow equations. These criteria are then applied to commonly used RK methods such as Jameson's 3-step and 4-step schemes and Wray's algorithm to identify the accuracy of the methods. For the spatial discretization, compact difference schemes are presented. The schemes are formulated in the operator-type to render themselves suitable for the Fourier analyses. The performance of the numerical methods is shown by numerical examples. These examples are detailed. described. The third case is a two-dimensional simulation of a Lamb vortex in an uniform flow. This calculation provides a realistic assessment of various finite difference schemes in terms of the conservation of the vortex strength and the harmonic content after travelling a substantial distance. The numerical implementation of Giles' non-refelctive equations coupled with the characteristic equations as the boundary condition is discussed in detail. Finally, the single vortex calculation is extended to simulate vortex pairing. For the distance between two vortices less than a threshold value, numerical results show crisp resolution of the vortex merging.

Yu, Sheng-Tao↗

Measurement of Temperature Fluctuations and Microscopic Growth Rates in a Silicon Floating Zone Under Microgravity

USA Several microgravity experiments on sounding rockets and the Space Shuttle have shown that time-dependent thermocapillary (Marangoni) convection is the major cause for the formation of dopant striations in floating-zone grown semiconductor crystals, at least in small-scale systems not employing RF heating. To quantify this correlation, a silicon floating-zone experiment was performed during the TEXUS36 flight (February 7, 1998) in the monoellipsoid mirror furnace TEM02-ELLI. During the experiment, temperature fluctuations in the silicon melt zone and the microscopic growth rate were simultaneously measured. Temperature fluctuations of 0.5 C - 0.7 C with main frequencies between 0.1 Hz and 0.3 Hz were detectable. The microscopic growth rate fluctuated considerably around the average growth rate of 1 mm/min: rates from 4 mm/min to negative values (backmelting) were observed. Dopant striations are clearly visible in the Sb-doped crystal. They were characterized by Spreading Resistance measurements and Differential Interference Contrast microscopy. The frequencies associated with the dopant inhomogeneities correspond quite well with those of the temperature fluctuations and microscopic growth rates. 3D numerical simulations were performed to predict the optimum position of the temperature sensor, to evaluate characteristic temperature amplitudes and frequencies, and to give insight into the instability mechanisms of Marangoni convection in this configuration. The simulations were in good agreement with the experimental values, showing temperature fluctuations with frequencies ? 0.25 Hz and amplitudes up to 1.8 C at a position equivalent to that of the sensor tip in the experiment.

Croell, Arne↗

Measurement of Temperature Fluctuations and Microscopic Growth Rates in a Silicon Floating Zone on TEXUS36

Several pg experiments on sounding rockets and the Space Shuttle have shown that time-dependent thermocapillary (Marangoni) convection is the major cause for the formation of dopant striations in floating-zone grown semiconductor crystals, at least in small-scale systems not employing RF heating. To quantify this correlation, a silicon floating-zone experiment was performed during the TEXUS36 flight (February 7, 1 998) in the monoellipsoid mirror furnace TEM02-ELLI. During the experiment, temperature fluctuations in the silicon melt zone and the microscopic growth rate were simultaneously measured. Temperature fluctuations of 0.5 C - 0.7 C with main frequencies between 0.1 Hz and 0.3Hz were detectable. The microscopic growth rate fluctuated considerably around the average growth rate of 1 mm/min: rates from 4mm/min to negative values (backmelting) were observed. Dopant striations are clearly visible in the Sb-doped crystal. They were characterized by Spreading Resistance measurements and Differential Interference Contrast microscopy. The frequencies associated with the dopant inhomogeneities correspond quite well with those of the temperature fluctuations and microscopic growth rates. 3D numerical simulations were performed to predict the optimum position of the temperature sensor, to evaluate characteristic temperature amplitudes and frequencies, and to give insight into the instability mechanisms of Marangoni convection in this configuration. The simulations were in good agreement with the experimental values, showing temperature fluctuations with frequencies? 0.25 Hz and amplitudes up to 1.8 C at a position equivalent to that of the sensor tip in the experiment.

Croell, Arne↗

Rocket-in-a-Duct Performance Analysis

An axisymmetric, 110 N class, rocket configured with a free expansion between the rocket nozzle and a surrounding duct was tested in an altitude simulation facility. The propellants were gaseous hydrogen and gaseous oxygen and the hardware consisted of a heat sink type copper rocket firing through copper ducts of various diameters and lengths. A secondary flow of nitrogen was introduced at the blind end of the duct to mix with the primary rocket mass flow in the duct. This flow was in the range of 0 to 10% of the primary massflow and its effect on nozzle performance was measured. The random measurement errors on thrust and massflow were within +/-1%. One dimensional equilibrium calculations were used to establish the possible theoretical performance of these rocket-in-a-duct nozzles. Although the scale of these tests was small, they simulated the relevant flow expansion physics at a modest experimental cost. Test results indicated that lower performance was obtained at higher free expansion area ratios and longer ducts, while, higher performance was obtained with the addition of secondary flow. There was a discernable peak in specific impulse efficiency at 4% secondary flow. The small scale of these tests resulted in low performance efficiencies, but prior numerical modeling of larger rocket-in-a-duct engines predicted performance that was comparable to that of optimized rocket nozzles. This remains to be proven in large-scale, rocket-in-a-duct tests.

Schneider, Steven J.↗

Three-Dimensional Effects in Modeling of Dual-Mode Scramjets

A numerical investigation of an experimental dual-mode scramjet configuration is performed. Both experimental and numerical results indicate significant upstream interaction for this case. Several computational cases are examined: these include the use of jet-to-jet symmetry and entire half-duct modeling. Grid convergence, turbulence modeling, and wall temperature effects are studied in terms of wall pressure predictions and flow-field characteristics. Wall pressure comparisons between CFD and experiment show fair agreement for the jet-to-jet case. However, further computations of the entire half-duct show the development of a large sidewall separation zone extending much further upstream than the separation zone at the duct centerline. This sidewall separation is the dominant feature in the CFD-generated flowfield but is not evident in the experimental data, resulting in a unfavorable comparison between CFD and experimental data. Current work aimed at resolving this issue and at further understanding asymmetric flow-structures in dual-mode flow-fields is discussed.

Rodriguez, C. G.↗

Three-Dimensional Planetary Surface Tracking Based on a Simple Ultra-Wideband Impulse-Radio Infrastructure

Several prototype ultra-wideband (UWB) impulse-radio (IR) tracking systems are currently under development at NASA Johnson Space Center (JSC). These systems are being studied for use in tracking of Lunar/Mars rovers and astronauts during early exploration missions when satellite navigation systems (such as GPS) are not available. To date, the systems that have been designed and tested are intended only for two-dimensional location and tracking, but these designs can all be extended to three-dimensional tracking with only minor modifications and increases in complexity. In this presentation, we will briefly review the design and performance of two of the current 2-D systems: one designed specifically for short-range, extremely high-precision tracking (approximately 1-2 cm resolution) and the other designed specifically for much longer range tracking with less stringent precision requirements (1-2 m resolution). We will then discuss a new multi-purpose system design based on a simple UWB-IR architecture that can be deployed easily on a planetary surface to support arbitrary three-dimensional localization and tracking applications. We will discuss utilization of this system as an infrastructure to provide both short-range and long-range tracking and analyze the localization performance of the system in several different configurations. We will give theoretical performance bounds for some canonical system configurations and compare these performance bounds with both numerical simulations of the system as well as actual experimental system performance evaluations.

Barton, Richard J.↗

Entropy Stable Method for the Euler Equations Revisited: Central Differencing via Entropy Splitting and SBP

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Ols- son & Oliger (1994), Gerritsen & Olsson (1996), and Yee et al. (2000) is revisited. The entropy splitting is a form of skew-symmetric splitting of the nonlinear Euler flux derivatives. Central differencing applied to the entropy splitting form of the Euler flux derivatives together with SBP difference operators will, hereafter, be referred to as entropy split schemes. This study is prompted by the recent growing interest in numerical methods for which a discrete entropy conservation law holds, a discrete global entropy conservation can be proved and/or the numerical method possesses a stable entropy in the framework of SBP difference operators and L2-energy norm estimate. The objective of this paper is to recast the entropy split scheme as the re- cent definition of an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for non- linear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grub- ber splitting.

Sjogreen, Bjorn↗

A Parametric Study of Spur Gear Dynamics

A parametric study of a spur gear system was performed through a numerical analysis approach. This study used the gear dynamic program DANST, a computer simulator, to determine the dynamic behavior of a spur gear system. The analytical results have taken the deflection of shafts and bearings into consideration for static analysis, and the influence of these deflections on gear dynamics was investigated. Damping in the gear system usually is an unknown quantity, but it has an important effect in resonance vibration. Typical values as reported in the literature were used in the present analysis. The dynamic response due to different damping factors was evaluated and compared. The effect of the contact ratio on spur gear dynamic load and dynamic stress was investigated through a parameter study. The contact ratio was varied over the range of 1.26 to 2.46 by adjusting the tooth addendum. Gears with contact ratio near 2.0 were found to have the most favorable dynamic performance.

Lin, Hsiang Hsi↗

Parallel spatial direct numerical simulations on the Intel iPSC/860 hypercube

The implementation and performance of a parallel spatial direct numerical simulation (PSDNS) approach on the Intel iPSC/860 hypercube is documented. The direct numerical simulation approach is used to compute spatially evolving disturbances associated with the laminar-to-turbulent transition in boundary-layer flows. The feasibility of using the PSDNS on the hypercube to perform transition studies is examined. The results indicate that the direct numerical simulation approach can effectively be parallelized on a distributed-memory parallel machine. By increasing the number of processors nearly ideal linear speedups are achieved with nonoptimized routines; slower than linear speedups are achieved with optimized (machine dependent library) routines. This slower than linear speedup results because the Fast Fourier Transform (FFT) routine dominates the computational cost and because the routine indicates less than ideal speedups. However with the machine-dependent routines the total computational cost decreases by a factor of 4 to 5 compared with standard FORTRAN routines. The computational cost increases linearly with spanwise wall-normal and streamwise grid refinements. The hypercube with 32 processors was estimated to require approximately twice the amount of Cray supercomputer single processor time to complete a comparable simulation; however it is estimated that a subgrid-scale model which reduces the required number of grid points and becomes a large-eddy simulation (PSLES) would reduce the computational cost and memory requirements by a factor of 10 over the PSDNS. This PSLES implementation would enable transition simulations on the hypercube at a reasonable computational cost.

Joslin, Ronald D.↗

A numerical experiment on the equilibrium and stability of a rotating galactic bar

A self-consistent, three-dimensional numerical experiment is performed on an N-body system whose initial state is a realization of a certain theoretical model of a rotating triaxial galaxy. The model is a stellar-dynamical counterpart of a uniformly rotating polytrope of index equal to 0.5. The aim of the experiment is to study the equilibrium of the system and, in particular, to test its stability. The experimental system behaves in the mean like a realization of the theoretical model for at least seven crossing times. The principal departure of the system from equilibrium is an oscillation which is identified as a radial pulsation. There is no indication in its behavior that the system is unstable with respect to anu mode with an e-folding time shorter than or of the order of two crossing times. Certain changes that occur in the state of the system are interpreted, with the aid of the theoretical model, as secular changes which result from a slight failure of our numerical methods to conserve the mass, energy, and angular momentum of the system; these effects are small enough that they do not vitiate the experiment on a dynamical time scale.

Miller, R. H.↗

Use It or Lose It: Skeletal Muscle Function and Performance Results from Space Shuttle

The Space Shuttle Program provided a wealth of valuable information regarding the adaptations of skeletal muscle to weightlessness. Studies conducted during the Extended Duration Orbiter Medical Project (EDOMP) represented ground breaking work on the effects of spaceflight on muscle form and function from applied human research to cellular adaptations. Results from detailed supplementary objective (DSO) 477 demonstrated that muscle strength losses could occur rapidly in response to short-duration spaceflight. The effects of spaceflight-induced unloading were primarily restricted to postural muscles such as those of the back as well as the knee extensors. DSO 606 provided evidence from MRI that the observed strength losses were partially accounted for by a reduction in the size of the individual muscles. Muscle biopsy studies conducted during DSO 475 were able to show muscle atrophy in individual muscle fibers from the quadriceps muscles. Reduced quadriceps muscle size and strength was also observed during the 17-d Life and Microgravity Spacelab mission aboard STS-78. Multiple maximal strength tests were conducted in flight on the calf muscles and it has been hypothesized that these high force contractions may have acted as a countermeasure. Muscle fiber mechanics were studied on calf muscle samples pre- and postflight. While some responses were crewmember specific, the general trend was that muscle fiber force production dropped and shortening velocity increased. The increased shortening velocity helped to maintain muscle fiber power. Numerous rodent studies performed during Shuttle missions suggest that many of the effects reported in Shuttle crewmembers could be due to lesions in the cellular signaling pathways that stimulate protein synthesis as well as an increase in the mechanisms that up-regulate protein breakdown. The results have important implications regarding the overall health and performance capabilities of future crewmembers that will venture beyond low-Earth orbit. Learning Objective: Overview of the Space Shuttle Program regarding adaptive changes in skeletal muscle function and performance, including what was learned from research and what was implemented for countermeasures.

Ryder, Jeffrey↗

Investigation of turbulent flow in highly curved ducts with application to turbomachinery components

Numerical predictions were performed using a semi-elliptic calculation procedure for the case of turbulent flow in passage through a 90 deg bend of square cross section. Two versions of the isotropic turbulent viscosity two equation k-epsilon model were used. The first, the wall function model (WFM), employs the logarithmic law-of-the-wall relation and the notion of equilibrium flow to set all the necessary boundary conditions at the first grid point adjacent to a solid wall. The second model, the Van Driest model (VDM), employs Prandtl's original mixing length formulation, in conjunction with Van Driest's semi-empirical relation for the mixing length, to calculate the turbulent viscosity in the near wall regions of the flow. In this case, boundary conditions for k and epsilon, required to calculate these quantities in the core of the flow, are obtained by matching the mixing length and Reynolds number model formulations in an overlapping region of the flow near the walls. In both cases the results obtained show an improvement over earlier calculations using an elliptic numerical procedure. This is attributed to the finer grids possible in the present work. Of the two models, the VDM formulation shows better overall conformity with the mean flow measurements. Neither model reproduces well the details of the stress distribution as a result of the implied isotropic turbulent viscosity.

Humphrey, J. A. C.↗

SSME LOX post flow analysis/fluid structure interaction. Volume 1: Flow analysis

The realization of measures to improve the performance of the Space Shuttle is, to a large extent, dependent on an improved understanding of the fluid flow phenomena occurring in the main engine. The overall arrangement of the primary components of the Space Shuttle Main Engine (SSME) are presented. The impingement of the hot gases from the transfer ducts onto the LOX posts causes them to be subjected to severe thermal and gas dynamic loads, which in the past have resulted in the occasional breakage of some elements of the outer row of posts during test firings of the engine, particularly at higher power levels. Large velocities in the gap between the LOX posts may also be a contributing factor in causing dynamic stability problems. The deforming structural response of the posts to the pressure loading may likely affect the gas flowfield by producing a moving flowfield boundary, thereby creating a dynamically coupled unsteady fluid-structure system. The objective was to investigate the three-dimensional, turbulent flow around a simplified SSME LOX post array using an existing Reynolds averaged Navier-Stokes flow solver and a suitable turbulence model to parameterize the turbulent shear stresses. Numerical computations were performed to analyze the effect on the flowfield of varying the spacing between the LOX posts, which were modeled as rigid, three-dimensional circular cylinders. The methodology used in the computations is described. Results are presented.

Burke, Roger W.↗

Comparison of ice-phase microphysical parameterization schemes using numerical simulations of tropical convection

The performance of several ice parameterizations has been evaluated through a numerical cloud model. Ice effects using different schemes are contrasted with each other and with an ice-free control by incorporating them into the cloud model and by applying them to simulations of tropical squall systems. The latter are simulated in 2D so that a large domain can be used to incorporate a complete anvil. Nonsquall-type convective lines are simulated in 3D owing to their smaller horizontal scale. It is concluded that inclusion of ice microphysics in the cloud model enhanced the agreement of the simulated convection with some features of observed convection, including the proportion of surface rainfall in the anvil region and the intensity and structure of the radar brightband near the melting level in the anvil. In the experiments with bulk microphysics, three ice categories produced much better results than two ice categories, which in turn was better than no ice. For the tropical squall-type and nonsquall-type systems the optimal mix was ice, snow, and graupel.

Mccumber, Michael↗

Effect of 3D Roughness Patch on Instability Amplification in a Supersonic Boundary Layer

Surface roughness is known to have a substantial impact on the aerothermodynamic loading of high-speed vehicles, particularly via its influence on the laminar-turbulent transition process within the boundary layer. Numerical simulations are performed to investigate the effects of a distributed region of densely packed, sinusoidal shape roughness elements on a Mach 3.5 flat plate boundary layer for flow conditions corresponding to the planned conditions of an upcoming experiment in the Mach 3.5 Supersonic Low Disturbance Tunnel at the NASA Langley Research Center. Analysis of convective instabilities in the wake of the roughness patch was reported in a previous paper and the current work extends that analysis to instability amplification across the length of the roughness patch. Quasi-parallel stability analysis of the modified boundary layer flow over the patch indicates two dominant families of unstable disturbances, namely, a group of high frequency modes that amplify in localized regions along the roughness patch and another group of low frequency modes that have smaller peak amplification rates but amplify steadily both above the roughness patch and in the wake region behind it. The results suggest that the amplification factors associated with the high-frequency modes are sufficiently low, at least for the roughness patches considered in this paper, so that these modes are unlikely to have a major influence on the transition process. The amplification of the low-frequency modes within the region of the roughness patch is further quantified via direct numerical simulations. Results confirm the strongly destabilizing influence of the roughness patch on the first mode instabilities, yielding an N-factor increment of ΔN ≈ 3.6 for a roughness patch length of eight wavelengths.

Laminar turbulent transition↗

Steady-State Cycle Deck Launcher Developed for Numerical Propulsion System Simulation

One of the objectives of NASA's High Performance Computing and Communications Program's (HPCCP) Numerical Propulsion System Simulation (NPSS) is to reduce the time and cost of generating aerothermal numerical representations of engines, called customer decks. These customer decks, which are delivered to airframe companies by various U.S. engine companies, numerically characterize an engine's performance as defined by the particular U.S. airframe manufacturer. Until recently, all numerical models were provided with a Fortran-compatible interface in compliance with the Society of Automotive Engineers (SAE) document AS681F, and data communication was performed via a standard, labeled common structure in compliance with AS681F. Recently, the SAE committee began to develop a new standard: AS681G. AS681G addresses multiple language requirements for customer decks along with alternative data communication techniques. Along with the SAE committee, the NPSS Steady-State Cycle Deck project team developed a standard Application Program Interface (API) supported by a graphical user interface. This work will result in Aerospace Recommended Practice 4868 (ARP4868). The Steady-State Cycle Deck work was validated against the Energy Efficient Engine customer deck, which is publicly available. The Energy Efficient Engine wrapper was used not only to validate ARP4868 but also to demonstrate how to wrap an existing customer deck. The graphical user interface for the Steady-State Cycle Deck facilitates the use of the new standard and makes it easier to design and analyze a customer deck. This software was developed following I. Jacobson's Object-Oriented Design methodology and is implemented in C++. The AS681G standard will establish a common generic interface for U.S. engine companies and airframe manufacturers. This will lead to more accurate cycle models, quicker model generation, and faster validation leading to specifications. The standard will facilitate cooperative work between industry and NASA. The NPSS Steady-State Cycle Deck team released a batch version of the Steady-State Cycle Deck in March 1996. Version 1.1 was released in June 1996. During fiscal 1997, NPSS accepted enhancements and modifications to the Steady-State Cycle Deck launcher. Consistent with NPSS' commercialization plan, these modifications will be done by a third party that can provide long-term software support.

VanDrei, Donald E.↗