SEARCH · Search NASA

Results for “numerical simulation”

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 757 records · Page 42

All sources NASA NTRS DOE OSTI

NASA NTRS1996

Numerical simulation of jet aerodynamics using the three-dimensional Navier-Stokes code PAB3D

This report presents a unified method for subsonic and supersonic jet analysis using the three-dimensional Navier-Stokes code PAB3D. The Navier-Stokes code was used to obtain solutions for axisymmetric jets with on-design operating conditions at Mach numbers ranging from 0.6 to 3.0, supersonic jets containing weak shocks and Mach disks, and supersonic jets with nonaxisymmetric nozzle exit geometries. This report discusses computational methods, code implementation, computed results, and comparisons with available experimental data. Very good agreement is shown between the numerical solutions and available experimental data over a wide range of operating conditions. The Navier-Stokes method using the standard Jones-Launder two-equation kappa-epsilon turbulence model can accurately predict jet flow, and such predictions are made without any modification to the published constants for the turbulence model.

Pao, S. Paul↗

NASA NTRS1992

Direct numerical simulation of laminar breakdown in high-speed, axisymmetric boundary layers

The compressible Navier-Stokes equations are solved using spectral collocation and high-order compact-difference techniques to simulate the laminar breakdown in high-speed, axisymmetric boundary-layer flow. Mach 4.5 flow along a hollow cylinder and Mach 6.8 flow along a sharp cone are considered. Data obtained replicate two previously unexplained phenomena, namely, the appearance of so-called 'rope-like waves' and 'the precursor transition effect', in which transitional flow originates near the critical layer well upstream of the transition location at the wall. The numerical data also reveal that neither of these effects can be explained, even qualitatively, by linear stability theory alone. It is shown that rope-like appearance arises from secondary instability. Certain features of the precursor transition effect also emerge from secondary instability but its nature is revealed to be fundamentally nonlinear.

Pruett, C. D.↗

NASA NTRS1999

User's Manual for LINER: FORTRAN Code for the Numerical Simulation of Plane Wave Propagation in a Lined Two-Dimensional Channel

LINER is a system of Fortran 77 codes which performs a 2D analysis of acoustic wave propagation and noise suppression in a rectangular channel with a continuous liner at the top wall. This new implementation is designed to streamline the usage of the several codes making up LINER, resulting in a useful design tool. Major input parameters are placed in two main data files, input.inc and nurn.prm. Output data appear in the form of ASCII files as well as a choice of GNUPLOT graphs. Section 2 briefly describes the physical model. Section 3 discusses the numerical methods; Section 4 gives a detailed account of program usage, including input formats and graphical options. A sample run is also provided. Finally, Section 5 briefly describes the individual program files.

Reichert, R, S.↗

NASA NTRS1973

A field investigation and numerical simulation of coastal fog

A field investigation of the microphysical and micrometeorological features of fogs occurring near Los Angeles and Vandenberg, California was conducted. Observations of wind speed and direction, temperature, dew point, vertical wind velocity, dew deposition, drop-size distribution, liquid water content, and haze and cloud nucleus concentration were obtained. These observations were initiated in late evening prior to fog formation and continued until the time of dissipation in both advection and radiation fogs. Data were also acquired in one valley fog and several dense haze situations. The behavior of these parameters prior to and during fog are discussed in detail. A two-dimensional numerical model was developed to investigate the formation and dissipation of advection fogs under the influence of horizontal variations in surface temperature. The model predicts the evolution of potential temperature, water vapor content, and liquid water content in a vertical plane as determined by vertical turbulent transfer and horizontal advection. Results are discussed from preliminary numerical experiments on the formation of warm-air advection fog and dissipation by natural and artificial heating from the surface.

Mack, E. J.↗

NASA NTRS1989

Two-dimensional numerical simulation of a Stirling engine heat exchanger

The first phase of an effort to develop multidimensional models of Stirling engine components is described; the ultimate goal is to model an entire engine working space. More specifically, parallel plate and tubular heat exchanger models with emphasis on the central part of the channel (i.e., ignoring hydrodynamic and thermal end effects) are described. The model assumes: laminar, incompressible flow with constant thermophysical properties. In addition, a constant axial temperature gradient is imposed. The governing equations, describing the model, were solved using Crank-Nicloson finite-difference scheme. Model predictions were compared with analytical solutions for oscillating/reversing flow and heat transfer in order to check numerical accuracy. Excellent agreement was obtained for the model predictions with analytical solutions available for both flow in circular tubes and between parallel plates. Also the heat transfer computational results are in good agreement with the heat transfer analytical results for parallel plates.

Ibrahim, Mounir B.↗

NASA NTRS1989

Two-dimensional numerical simulation of a Stirling engine heat exchanger

The first phase of an effort to develop multidimensional models of Stirling engine components is described. The ultimate goal is to model an entire engine working space. Parallel plate and tubular heat exchanger models are described, with emphasis on the central part of the channel (i.e., ignoring hydrodynamic and thermal end effects). The model assumes laminar, incompressible flow with constant thermophysical properties. In addition, a constant axial temperature gradient is imposed. The governing equations describing the model have been solved using the Crack-Nicloson finite-difference scheme. Model predictions are compared with analytical solutions for oscillating/reversing flow and heat transfer in order to check numerical accuracy. Excellent agreement is obtained for flow both in circular tubes and between parallel plates. The computational heat transfer results are in good agreement with the analytical heat transfer results for parallel plates.

Ibrahim, Mounir↗

NASA NTRS1990

Numerical simulation of high-incidence flow over the F-18 aircraft

Navier-Stokes solutions have been obtained using the Chimera overset grid scheme for flow over the wing, fuselage, and wing leading-edge extension (LEX) of the F-18 aircraft at high incidence. Solutions are also presented for flow over the fuselage forebody at high angles of attack. The solutions are for turbulent flows at high-Reynolds number flight-test conditions, and are compared with available qualitative and quantitative experimental data. Comparisons of predicted surface flow patterns, off-surface flow visualizations, and surface-pressure distributions are in good agreement with flight-test data. The ability of the numerical method to predict the bursting of the LEX vortex as it encounters the adverse pressure gradient field of the wing is demonstrated.

Cummings, Russell M.↗

NASA NTRS1991

The role of the sea-surface temperature distribution on numerically simulated cyclogenesis during ERICA

The goal was to quantify the extent to which a sea surface temperature (SST) front can influence cyclogenesis. The approach was to use the Drexel Limited-Area Mesoscale Prediction System (LAMPS) dynamical model to simulate cyclogenesis over various SST fields. Research during the past year focused on the development and testing of a four dimensional data assimilation (FDDA) technique within LAMPS. The technique is a continuous dynamical assimilation where forcing terms are added to the governing model equations to gradually nudge the model solution toward a gridded analysis. Here, the nudging is used as a dynamic initialization tool during a 12 hour preforecast to generate model balanced initial conditions for a subsequent 24 hour numerical prediction. Tests were performed to determine which variables to nudge and how to specify the four dimensional weighting function used to scale the nudging terms. To date, optimal results were obtained by nudging the u and v components of the wind along with the potential temperature. The weighting function ranged from 0 to 1 and varies in time as a quadratic polynomial. It was initialized at 0, reached its maximum at 9 hours into the preforecast, and fell back at 0 to 12 hours. The nudging terms are included in the model equations for all grid points except those within the model predicted oceanic boundary layer. This design attempts to confine changes imposed by the specified SST field to the oceanic boundary layer during the preforecast period.

Lapenta, William M.↗

NASA NTRS1995

Numerical simulation of heat and mass transport during space crystal growth with MEPHISTO

The MEPHISTO space experiments are collaborative United States and French investigations aimed at understanding the fundamentals of crystal growth. Microgravity experiments were conducted aboard the USMP-1 and -2 missions on STS-52 and 62 in October 1992 and March 1994 respectively. MEPHISTO is a French designed and built Bridgman type furnace which uses the Seebeck technique to monitor the solid/liquid interface temperature and Peltier pulsing to mark the location and shape of the solid/liquid interface. In this paper the Bridgman growth of Sn-Bi and Bi-Sn under terrestrial and microgravity conditions is modeled using the finite element code, FIDAP*. The numerical model considers fully coupled heat and mass transport, fluid motion and solid/liquid phase changes in the crystal growth process. The primary goals of this work are: to provide a quantitative study of the thermal buoyancy-induced convection in the melt for the two flight experiments; to compare the vertical and horizontal growth configurations and systematically evaluate the effects of various gravity levels on the solute segregation. Numerical results of the vertical and horizontal Bridgman growth configurations are presented.

Yao, Minwu↗

NASA NTRS1991

Shear-induced instability and arch filament eruption - A magnetohydrodynamic (MHD) numerical simulation

A situation wherein a bipolar magnetic field embedded in a stratified solar atmosphere undergoes symmetrical shear motion at the footpoints is investigated via a 2D (nonplanar) MHD simulation. It was found that the vertical plasma flow velocities grow exponentially, leading to a new type of global MHD instability. The growth rate increases almost linearly until it reaches the same order of magnitude as the Alfven speed. Then a nonlinear MHD instability occurs beyond this point. It was found that the central loops are pinched by opposing Lorentz forces, and the outer closed loops stretch upward with the vertically-rising mass flow. The nonlinear dynamical shearing instability is illustrated by a numerical example that is given for three different values of the plasma beta that span several orders of magnitude.

Wu, S. T.↗

NASA NTRS1998

A Numerical Simulator for Three-Dimensional Flows Through Vibrating Blade Rows

The three-dimensional, multi-stage, unsteady, turbomachinery analysis, TURBO, has been extended to predict the aeroelastic and aeroacoustic response behaviors of a single blade row operating within a cylindrical annular duct. In particular, a blade vibration capability has been incorporated so that the TURBO analysis can be applied over a solution domain that deforms with a vibratory blade motion. Also, unsteady far-field conditions have been implemented to render the computational boundaries at inlet and exit transparent to outgoing unsteady disturbances. The modified TURBO analysis is applied herein to predict unsteady subsonic and transonic flows. The intent is to partially validate this nonlinear analysis for blade flutter applications, via numerical results for benchmark unsteady flows, and to demonstrate the analysis for a realistic fan rotor. For these purposes, we have considered unsteady subsonic flows through a 3D version of the 10th Standard Cascade, and unsteady transonic flows through the first stage rotor of the NASA Lewis, Rotor 67, two-stage fan.

Chuang, H. Andrew↗

NASA NTRS1996

Numerical simulation of iodine speciation in relation to water disinfection aboard manned spacecraft I. Equilibria

Elemental iodine (I2) is currently used as the drinking water disinfectant aboard the Shuttle Orbiter and will also be incorporated into the water recovery and distribution system for the International Space Station Alpha. Controlled release of I2 is achieved using the Microbial Check Valve (MCV), a flow-through device containing an iodinated polymer which imparts a bacteriostatic residual concentration of approximately 2mg/L to the aqueous stream. During regeneration of MCV canisters, I2 concentrations of approximately 300 mg/L are used. Dissolved iodine undergoes a series of hydrolytic disproportionation and related reactions which result in the formation of an array of inorganic species including: I-, I3-, HOI, OI-, IO3-, HIO3, I2OH-, I2O(-2), and H2OI+. Numerical estimation of the steady-state distribution of inorganic iodine containing species in pure water at 25 degrees C has been achieved by simultaneous solution of the multiple equilibrium expressions as a function of pH. The results are reported herein.

NASA Discipline Environmental Health↗

NASA NTRS1994

Meso-beta scale numerical simulation studies of terrain-induced jet streak mass/momentum perturbations

An in-depth analysis of observed gravity waves and their relationship to precipitation bands over the Montana mesonetwork during the 1981 CCOPE case study indicates that there were two episodes of coherent internal gravity waves. One of the fundamental unanswered questions from this research, however, concerns the dynamical processes which generated the observed waves, all of which originated from the region encompassing the borders of Montana, Idaho, and Wyoming. While geostrophic adjustment, shearing instability, and terrain where all implicated separately or in concert as possible wave generation mechanisms, the lack of upper-air data within the wave genesis region made it difficult to rigorously define the genesis processes from observations alone. In this report we employ a mesoscale numerical model to help diagnose the intricate early wave generation mechanisms during the first observed wave episode.

Lin, Yuh-Lang↗

NASA NTRS1988

Numerical simulation of the integrated space shuttle vehicle in ascent

A simulation of the flow about the integrated space shuttle vehicle in ascent mode has been undertaken for various flight conditions using the Chimera composite grid discretization approach. Overset body-conforming grids were used to represent each geometric component, and an implicit approximately factored finite-difference procedure was used to solve the three-dimensional thin-layer Navier-Stokes equations. The computational results have been compared with both wind tunnel and flight test data. Although relatively good agreement is obtained with the experimental data, further refinement and evaluation of numerical error is under way.

Buning, P. G.↗

NASA NTRS1986

Numerical simulation of electrophoresis separation processes

A new Petrov-Galerkin finite element formulation has been proposed for transient convection-diffusion problems. Most Petrov-Galerkin formulations take into account the spatial discretization, and the weighting functions so developed give satisfactory solutions for steady state problems. Though these schemes can be used for transient problems, there is scope for improvement. The schemes proposed here, which consider temporal as well as spatial discretization, provide improved solutions. Electrophoresis, which involves the motion of charged entities under the influence of an applied electric field, is governed by equations similiar to those encountered in fluid flow problems, i.e., transient convection-diffusion equations. Test problems are solved in electrophoresis and fluid flow. The results obtained are satisfactory. It is also expected that these schemes, suitably adapted, will improve the numerical solutions of the compressible Euler and the Navier-Stokes equations.

Ganjoo, D. K.↗

NASA NTRS1994

Numerical Simulation of Shock-Induced Combustion Past Blunt Bodies Using Shock-Fitting Technique

Two-dimensional axisymmetric, reacting viscous flow over blunt projectiles is computed to study shock induced combustion at Mach 5.11 and Mach 6.46 in hydrogen-air mixture. A finite-difference, shock-fitting method is used to solve the complete set of Navier Stokes and species conservation equations. In this approach, the bow shock represents a boundary of the computational domain and is treated as a discontinuity across which Rankine-Hugoniot conditions are applied. All interior details of the flow such as compression waves, reaction front, and the wall boundary layer are captured automatically in the solution. Since shock-fitting approach reduces the amount of artificial dissipation, all the intricate details of the flow are captured much more clearly than has been possible with the shock-capturing approach. This has allowed an improved understanding of the physics of shock-induced combustion over blunt projectiles and the numerical results can now be explained more readily with one dimensional wave-interaction model than before.

Ahuja, J. K.↗

NASA NTRS1994

Numerical simulation of shock-induced combustion past blunt bodies using shock-fitting technique

Two-dimensional axisymmetric, reacting viscous flow over blunt projectiles is computed to study shock-induced combustion at Mach 5.11 and Mach 6.46 in hydrogen-air mixture. A finite-difference, shock-fitting method is used to solve the complete set of Navier-Stokes and species conservation equations. In this approach, the bow shock represents a boundary of the computational domain and is treated as a discontinuity across which Rankine-Hugoniot conditions are applied. All interior details of the flow such as compression waves, reaction front, and the wall boundary layer are captured automatically in the solution. Since shock-fitting approach reduces the amount of artificial dissipation, all the intricate details of the flow are captured much more clearly than has been possible with the shock-capturing approach. This has allowed an improved understanding of the physics of shock-induced combustion over blunt projectiles and the numerical results can now be explained more readily with one-dimensional wave-interaction model than before.

Ahuja, J. K.↗

NASA NTRS1997

Preconditioning for Numerical Simulation of Low Mach Number Three-Dimensional Viscous Turbomachinery Flows

A preconditioning scheme has been implemented into a three-dimensional viscous computational fluid dynamics code for turbomachine blade rows. The preconditioning allows the code, originally developed for simulating compressible flow fields, to be applied to nearly-incompressible, low Mach number flows. A brief description is given of the compressible Navier-Stokes equations for a rotating coordinate system, along with the preconditioning method employed. Details about the conservative formulation of artificial dissipation are provided, and different artificial dissipation schemes are discussed and compared. The preconditioned code was applied to a well-documented case involving the NASA large low-speed centrifugal compressor for which detailed experimental data are available for comparison. Performance and flow field data are compared for the near-design operating point of the compressor, with generally good agreement between computation and experiment. Further, significant differences between computational results for the different numerical implementations, revealing different levels of solution accuracy, are discussed.

Tweedt, Daniel L.↗