Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer experiments”

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 1,657 records · Page 92

Monthly mean forecast experiments with the GISS model

The GISS general circulation model was used to compute global monthly mean forecasts for January 1973, 1974, and 1975 from initial conditions on the first day of each month and constant sea surface temperatures. Forecasts were evaluated in terms of global and hemispheric energetics, zonally averaged meridional and vertical profiles, forecast error statistics, and monthly mean synoptic fields. Although it generated a realistic mean meridional structure, the model did not adequately reproduce the observed interannual variations in the large scale monthly mean energetics and zonally averaged circulation. The monthly mean sea level pressure field was not predicted satisfactorily, but annual changes in the Icelandic low were simulated. The impact of temporal sea surface temperature variations on the forecasts was investigated by comparing two parallel forecasts for January 1974, one using climatological ocean temperatures and the other observed daily ocean temperatures. The use of daily updated sea surface temperatures produced no discernible beneficial effect.

Spar, J.↗

Monthly mean forecast experiments with the GISS model - Correction

One of the three monthly mean forecasts computed previously by Spar et al., (1976) with the GISS (Goddard Institute for Space Studies) general circulation model was started from incorrect initial conditions. For the January 1975 prediction, a climatological mean January data set was accidentally substituted for the required January 1, 1975 initial data. In the present paper, the forecast for January 1975 is computed with the proper initial conditions, and the corrected results are presented. It is shown that the use of correct initial conditions improves the gross forecast mean energetics. In general, the use of correct initial conditions has resulted in a modest improvement in the monthly mean forecast for January 1975, but has not altered any of the previous conclusions by Spar et al.

Spar, J.↗

Prospects for computational aerodynamics

The integral representations approach, for the solution of the Navier-Stokes equations is discussed as well as experience in its development and in applying available finite-difference and finite-element techniques to the treatment of three-dimensional problems, and the computation of turbulent flow. The magnitude of efforts required to develop turbulence models and three-dimensional algorithms indicates that the computational fluid dynamics research must have a broad base. Broader access to modern computing facilities that are in existence within NASA should be promoted for active researchers not directly affiliated with that agency.

Wu, J. C.↗

Computation of transonic separated wing flows using an Euler/Navier-Stokes zonal approach

A computer program called Transonic Navier Stokes (TNS) has been developed which solves the Euler/Navier-Stokes equations around wings using a zonal grid approach. In the present zonal scheme, the physical domain of interest is divided into several subdomains called zones and the governing equations are solved interactively. The advantages of the Zonal Grid approach are as follows: (1) the grid for any subdomain can be generated easily; (2) grids can be, in a sense, adapted to the solution; (3) different equation sets can be used in different zones; and, (4) this approach allows for a convenient data base organization scheme. Using this code, separated flows on a NACA 0012 section wing and on the NASA Ames WING C have been computed. First, the effects of turbulence and artificial dissipation models incorporated into the code are assessed by comparing the TNS results with other CFD codes and experiments. Then a series of flow cases is described where data are available. The computed results, including cases with shock-induced separation, are in good agreement with experimental data. Finally, some futuristic cases are presented to demonstrate the abilities of the code for massively separated cases which do not have experimental data.

Kaynak, Uenver↗

Hypersonic rarefied flow about a compression corner - DSMC simulation and experiment

The phenomena of shock/boundary-layer interactions and flow separation are investigated using both computational and experimental methods for low-density hypersonic flow about two-dimensional compression corners. The numerical calculations are made with the direct simulation Monte Carlo (DSMC) method. Experimental measurements provide information concerning the flowfield structure and surface flow patterns by means of gas glow discharge and oil flow pictures, respectively. Comparison of the two data sets provides a qualitative basis for assessing the ability of the DSMC method to describe such flows.

Moss, James N.↗

CFD validation experiments at UTRC

Wind tunnel tests are performed in order to validate a computational fluid dynamics code. A large scale, two dimensional separation bubble is created on a flat plate, and low speed, turbulent flow is used. Extensive data sets are obtained with a nonintrusive laser velocimeter in addition to wall static pressure, total pressure, and hot-film measurements. All data and results are presented in a viewgraph format.

Patrick, William P.↗

Airfoil Ice-Accretion Aerodynamics Simulation

NASA Glenn Research Center, ONERA, and the University of Illinois are conducting a major research program whose goal is to improve our understanding of the aerodynamic scaling of ice accretions on airfoils. The program when it is completed will result in validated scaled simulation methods that produce the essential aerodynamic features of the full-scale iced-airfoil. This research will provide some of the first, high-fidelity, full-scale, iced-airfoil aerodynamic data. An initial study classified ice accretions based on their aerodynamics into four types: roughness, streamwise ice, horn ice, and spanwise-ridge ice. Subscale testing using a NACA 23012 airfoil was performed in the NASA IRT and University of Illinois wind tunnel to better understand the aerodynamics of these ice types and to test various levels of ice simulation fidelity. These studies are briefly reviewed here and have been presented in more detail in other papers. Based on these results, full-scale testing at the ONERA F1 tunnel using cast ice shapes obtained from molds taken in the IRT will provide full-scale iced airfoil data from full-scale ice accretions. Using these data as a baseline, the final step is to validate the simulation methods in scale in the Illinois wind tunnel. Computational ice accretion methods including LEWICE and ONICE have been used to guide the experiments and are briefly described and results shown. When full-scale and simulation aerodynamic results are available, these data will be used to further develop computational tools. Thus the purpose of the paper is to present an overview of the program and key results to date.

Bragg, Michael B.↗

Macrosegregation Resulting from Directional Solidification Through an Abrupt Change in Cross-Sections

Simulations of the directional solidification of two hypoeutectic alloys (Al-7Si alloy and Al-19Cu) and resulting macrosegregation patterns are presented. The casting geometries include abrupt changes in cross-section from a larger width of 9.5 mm to a narrower 3.2 mm width then through an expansion back to a width of 9.5 mm. The alloys were chosen as model alloys because they have similar solidification shrinkages, but the effect of Cu on changing the density of the liquid alloy is about an order of magnitude greater than that of Si. The simulations compare well with experimental castings that were directionally solidified in a graphite mold in a Bridgman furnace. In addition to the simulations of the directional solidification in graphite molds, some simulations were effected for solidification in an alumina mold. This study showed that the mold must be included in numerical simulations of directional solidification because of its effect on the temperature field and solidification. For the model alloys used for the study, the simulations clearly show the interaction of the convection field with the solidifying alloys to produce a macrosegregation pattern known as "steepling" in sections with a uniform width. Details of the complex convection- and segregation-patterns at both the contraction and expansion of the cross-sectional area are revealed by the computer simulations. The convection and solidification through the expansions suggest a possible mechanism for the formation of stray grains. The computer simulations and the experimental castings have been part of on-going ground-based research with the goal of providing necessary background for eventual experiments aboard the ISS. For casting practitioners, the results of the simulations demonstrate that computer simulations should be applied to reveal interactions between alloy solidification properties, solidification conditions, and mold geometries on macrosegregation. The simulations also presents the possibility of engineering the mold-material to avoid, or mitigate, the effects of thermosolutal convection and macrosegregation by selecting a mold material with suitable thermal properties, especially its thermal conductivity.

Lauer, M.↗

Computing Rotational, Rovibrational, and Vibrational Spectra for Astronomical Observations: High Accuracy Line Lists for High Temperatures, Limited Line Lists for Biosignature Molecules, and PAH Emission Spectra

Over the last several years, our group has been involved in developing approaches to compute highly accurate spectroscopic constants and vibrational frequencies for small transient molecules that may be used in the interpretation and assignment of high-resolution laboratory experiments as well as high-resolution astronomical spectra. Additionally, we have used the computed spectroscopic constants to simulate purely rotational and rovibrational spectra so that these may be compared directly with high- resolution astronomical observations, and we have worked on developing approaches that can be applied to much larger molecules, such as PAHs, where we can explicitly determine an harmonic corrections to vibrational modes as well as take into account intensity sharing due to resonances. Another part of our work in spectroscopic signatures involves computing highly accurate line lists for common molecules, such as CO2, SO2, and NH3, which occur in many astrophysical environments, including the atmospheres of exoplanets, and often need to have their lines identified in high- resolution observations in order to determine which lines are due to other molecules. In order to characterize the atmospheres of hot exoplanets, these line lists need to be very accurate and extend to very high energies. I will discuss our latest work in these areas of astrochemical spectroscopy research.

Lee, Timothy J.↗

Computing rovibrational, vibrational, and cascade emission spectra for comparison to astronomical observations

Over the last several years, our group has been involved in developing approaches to compute highly accurate spectroscopic constants and vibrational frequencies for small transient molecules that may be used in the interpretation and assignment of high-resolution laboratory experiments as well as high-resolution astronomical spectra. Additionally, we have used the computed spectroscopic constants to simulate purely rotational and rovibrational spectra so that these may be compared directly with high-resolution astronomical observations, and we have worked on developing approaches that can be applied to much larger molecules, such as polycyclic aromatic hydrocarbon (PAH) molecules, where we can explicitly determine anharmonic corrections to vibrational frequencies as well as take into account intensity sharing due to resonances. In addition, for PAH molecules we have used the anharmonic data to compute a library of temperature-dependent anharmonic vibrational spectra, with proper inclusion of polyad resonances, which can then be used to model the cascade emission spectra of PAH molecules – the type of spectra directly observed by astronomers. Another part of our work in spectroscopic signatures involves computing highly accurate line lists for common molecules, such as CO2, SO2, and NH3, which occur in many astrophysical environments, including the atmospheres of exoplanents, and often need to have their lines identified in high-resolution observations in order to determine which lines are due to other molecules. In order to characterize the atmospheres of hot exoplanets, these line lists need to be very accurate and extend to very high energies. I will discuss our latest work in these areas of astrochemical spectroscopy research.

Timothy J Lee↗

Prevention of Over-Pressurization During Combustion in a Sealed Chamber

The combustion of flammable material in a sealed chamber invariably leads to an initial pressure rise in the volume. The pressure rise is due to the increase in the total number of gaseous moles (condensed fuel plus chamber oxygen combining to form gaseous carbon dioxide and water vapor) and, most importantly, the temperature rise of the gas in the chamber. Though the rise in temperature and pressure would reduce with time after flame extinguishment due to the absorption of heat by the walls and contents of the sealed spacecraft, the initial pressure rise from a fire, if large enough, could lead to a vehicle over-pressure and the release of gas through the pressure relief valve. This paper presents a simple lumped-parameter model of the pressure rise in a sealed chamber resulting from the heat release during combustion. The transient model considers the increase in gaseous moles due to combustion, and heat transfer to the chamber walls by convection and radiation and to the fuel-sample holder by conduction, as a function of the burning rate of the material. The results of the model are compared to the pressure rise in an experimental chamber during flame spread tests as well as to the pressure falloff after flame extinguishment. The experiments involve flame spread over thin solid fuel samples. Estimates of the heat release rate profiles for input to the model come from the assumed stoichiometric burning of the fuel along with the observed flame spread behavior. The sensitivity of the model to predict maximum chamber pressure is determined with respect to the uncertainties in input parameters. Model predictions are also presented for the pressure profile anticipated in the Fire Safety-1 experiment, a material flammability and fire safety experiment proposed for the European Space Agency (ESA) Automated Transfer Vehicle (ATV). Computations are done for a range of scenarios including various initial pressures and sample sizes. Based on these results, various mitigation approaches are suggested to prevent vehicle over-pressurization and help guide the definition of the space experiment.

Gokoglu, Suleyman A.↗

Turbulence Model Validation through Joint Experimental / Computational Studies of Separated Flow Over a Three-Dimensional Tapered Bump: Part II - Computational Investigation

This report serves as documentation of activities related to Computational Fluid Dynamics (CFD) analysis applied to the Smooth Body Separation Experiment designed to produce validation-quality experimental data from wind-tunnel tests performed at the University of Notre Dame. All CFD analysis utilized the Boeing CFD (BCFD) flow solver. BCFD was used to analyze flow over the Boeing Speed Bump with second-order accuracy in space and time. Inviscid fluxes were computed using the HLLE scheme for steady-state simulations, and the Low-Dissipation Roe scheme for time-accurate simulations. The working fluid was modeled as calorically perfect air using Sutherland's law for viscosity. Thermal diffusivity was determined by assuming a constant molecular Prandtl number of 0.72, and a turbulent Prandtl number of 0.9. Reynolds Averaged Navier-Stokes (RANS) simulations were conducted using the one-equation Spalart-Allmaras turbulence model with rotation correction3 and the quadratic constitutive relation (SA-RC-QCR), as well as Menter's two-equation, Shear-Stress Transport (SST) model. Time-accurate computations were conducted using Delayed Detached eddy Simulation (DDES) with the one-equation RANS model, and Large Eddy Stress balance (LESb) with the two-equation RANS model. Implicit dual-time stepping was used to achieve second-order temporal accuracy. The experimental investigation conducted during the course of this study is described in detail in an accompanying report.

Flow separation↗

Computer simulation of supersonic rarefied gas flow in the transition region, about a spherical probe; a Monte Carlo approach with application to rocket-borne ion probe experiments

This report describes a Monte Carlo simulation of transition flow around a sphere. Conditions for the simulation correspond to neutral monatomic molecules at two altitudes (70 and 75 km) in the D region of the ionosphere. Results are presented in the form of density contours, velocity vector plots and density, velocity and temperature profiles for the two altitudes. Contours and density profiles are related to independent Monte Carlo and experimental studies, and drag coefficients are calculated and compared with available experimental data. The small computer used is a PDP-15 with 16 K of core, and a typical run for 75 km requires five iterations, each taking five hours. The results are recorded on DECTAPE to be printed when required, and the program provides error estimates for any flow field parameter.

Horton, B. E.↗

Exploiting the Temperature Dependence of Magnetic Susceptibility to Control Convection in Fundamental Studies of Solidification Phenomena

It is well known that convection is a dominant mass transport mechanism when materials are solidified on Earth's surface. This convection is caused by gradients in density (and therefore gravitational force) that are brought about by gradients in temperature, composition or both. Diffusion of solute is therefore dwarfed by convection and the study of fundamental parameters, such as dendrite tip shape and growth velocity in the absence of convection is nearly impossible. Significant experimental work has therefore been carried out in orbiting laboratories with the intent of minimizing convection by minimizing gravity. One of the best known experiments of this kind is the Isothermal Dendritic Growth Experiment (IDGE), supported by NASA. Naturally such experiments are costly and one objective of the present investigation is to develop an experimental method whereby convection can be- halted, in solidification and other experiments, on the surface. A second objective is to use the method to minimize convection resulting from the residual accelerations suffered by experiments in microgravity. The method to be used to minimize convection relies on the dependence of the magnetic susceptibility of a fluid on temperature or composition (whichever is driving convection). All materials experience a force when placed in a magnetic field gradient. The direction and magnitude of that force depend on the magnetic susceptibility of the material. Consequently the force will vary if the susceptibility varies with temperature or composition. With a magnetic field gradient in the right direction (typically upward) and of the right magnitude, this variation in the magnetic force can be made to exactly cancel the variation in the gravitational force. Expressed another way, normal buoyancy is exactly countered by a "magnetic buoyancy". To demonstrate the principle, a solution of MnC12 in water has been used. First the variation of the susceptibility of this paramagnetic solution with temperature and concentration was measured. Then a "cell", containing this solution and 50mm long by 15mm high by 155mm wide, was placed in a superconducting magnet at Marshall Space Flight Center. The magnetic field was measured at various positions within the bore of the magnet using a Hall effect probe. In this way, a position was found where the magnetic field gradient was predominantly upward; the magnitude of the gradient could then be adjusted by adjusting the current of the magnet. The ends of the cell consisted of machined copper blocks maintained at controlled temperatures by circulating water from constant temperature baths. The walls of the cell were of rectangular section glass tubing so that the cell contents could be seen. Velocities arising from thermal gradients within the cell were measured by particle image velocimetry (PIV). Particles used for this purpose were silver-coated hollow glass spheres of micrometers diameter and nearly the same density as the solution. A central vertical plane of the cell was illuminated by a laser beam passing through a cylindrical lens. Digital images of the particles were captured on a CCD camera and fed to a computer so that frame-to-frame movements of particles traveling with the fluid were captured. These images were employed to compute velocity maps using commercial PIV software. In a typical experiment the cold end of the cell was maintained at 10C and the warm end at 30 C. With no current in the magnet, i.e.- with natural -convection allowed to occur, the fluid was observed to circulate with an average speed of approximately 0.3 millimeters per second. It was visually apparent that this circulation was diminished as the current was increased. At currents of approximately 20A the flow was halted, to within the precision of the PIV measurements. At yet higher currents the convection was reversed with the hotter solution sinking and the cooler solution rising. At 40A this reversed convection had speeds averaging 0.43 millimeters per second. The measurements of susceptibility and density allow an estimate of the field gradient necessary to halt convection in the experiment. That estimate was 7.8T (squared) per meter and the convection was observed to halt in the magnet at a current giving 7.21T (squared) per meter from the magnetic field measurements. Calculations of the flow have been carried out using the computational fluid dynamics software FLUENT and show good agreement with the measurements.

Seybert, C. D.↗

Analysis of an advanced ducted propeller subsonic inlet

A time marching Navier-Stokes code called PARC (PARC2D for 2-D/axisymmetric and PARC3D for 3-D flow simulations) was validated for an advanced ducted propeller (ADP) subsonic inlet. The code validation for an advanced ducted propeller (ADP) subsonic inlet. The code validation was implemented for a non-separated flow condition associated with the inlet operating at angles-of-attack of 0 and 25 degrees. The inlet test data were obtained in the 9 x 15 ft Low Speed Wind Tunnel at NASA Lewis Research Center as part of a cooperative study with Pratt and Whitney. The experimental study focused on the ADP inlet performance for take-off and approach conditions. The inlet was tested at a free stream Mach number of 0.2, at angles-of-attack between O and 35 degrees, and at a maximum propeller speed of 12,000 RPM which induced a corrected air flow rate of about 46 lb/sec based on standard day conditions. The computational grid and flow boundary conditions (BC) were based on the actual inlet geometry and the funnel flow conditions. At the propeller face, two types of BC's were applied: a mass flow BC and a fixed flow properties BC. The fixed flow properties BC was based on a combination of data obtained from the experiment and calculations using a potential flow code. Comparison of the computational results with the test data indicates that the PARC code with the propeller face fixed flow properties BC provided a better prediction of the inlet surface static pressures than the predictions when the mass flow BC was used. For an angle-of-attack of 0 degrees, the PARC2D code with the propeller face mass flow BC provided a good prediction of inlet static pressures except in the region of high pressure gradient. With the propeller face fixed flow properties BC, the PARC2D code provided a good prediction of the inlet static pressures. For an angle-of-attack of 25 degrees with the mass flow BC, the PARC3D code predicted statis pressures which deviated significantly from the test data; however, with the fixed flow properties BC, a good comparison with the test data was obtained.

Iek, Chanthy↗

Space station Simulation Computer System (SCS) study for NASA/MSFC. Volume 1: Overview and summary

NASA's Space Station Freedom Program (SSFP) planning efforts have identified a need for a payload training simulator system to serve as both a training facility and as a demonstrator to validate operational concepts. The envisioned Marshall Space Flight Center (MSFC) Payload Training Complex (PTC) required to meet this need will train the space station payload scientists, station scientists, and ground controllers to operate the wide variety of experiments that will be onboard the Space Station Freedom. The Simulation Computer System (SCS) is the computer hardware, software, and workstations that will support the Payload Training Complex at MSFC. The purpose of this SCS study is to investigate issues related to the SCS, alternative requirements, simulator approaches, and state-of-the-art technologies to develop candidate concepts and designs. This study was performed August 1988 to October 1989. Thus, the results are based on the SSFP August 1989 baseline, i.e., pre-Langley configuration/budget review (C/BR) baseline. Some terms, e.g., combined trainer, are being redefined. An overview of the study activities and a summary of study results are given here.

Source record↗

Progress in Computational Aeroelasticity Using High Fidelity Flow and Structural Equations on Parallel Computers

Aeroelasticity which involves strong coupling of fluids, structures and controls is an important element in designing an aircraft. Computational aeroelasticity using low fidelity methods such as the linear aerodynamic flow equations coupled with the modal structural equations are well advanced. Though these low fidelity approaches are computationally less intensive, they are not adequate for the analysis of modern aircraft such as High Speed Civil Transport (HSCT) and Advanced Subsonic Transport (AST) which can experience complex flow/structure interactions. HSCT can experience vortex induced aeroelastic oscillations whereas AST can experience transonic buffet associated structural oscillations. Both aircraft may experience a dip in the flutter speed at the transonic regime. For accurate aeroelastic computations at these complex fluid/structure interaction situations, high fidelity equations such as the Navier-Stokes for fluids and the finite-elements for structures are needed. Computations using these high fidelity equations require large computational resources both in memory and speed. Current conventional supercomputers have reached their limitations both in memory and speed. As a result, parallel computers have evolved to overcome the limitations of conventional computers. This paper will address the transition that is taking place in computational aeroelasticity from conventional computers to parallel computers. The paper will address special techniques needed to take advantage of the architecture of new parallel computers. Results will be illustrated from computations made on iPSC/860 and IBM SP2 computer by using ENASERO code that directly couples the Euler/Navier-Stokes flow equations with high resolution finite-element structural equations.

Guruswamy, Guru P.↗