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 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS2006

Scalar Dissipation Modeling for Passive and Active Scalars: a priori Study Using Direct Numerical Simulation

Transitional databases from Direct Numerical Simulation (DNS) of three-dimensional mixing layers for single-phase flows and two-phase flows with evaporation are analyzed and used to examine the typical hypothesis that the scalar dissipation Probability Distribution Function (PDF) may be modeled as a Gaussian. The databases encompass a single-component fuel and four multicomponent fuels, two initial Reynolds numbers (Re), two mass loadings for two-phase flows and two free-stream gas temperatures. Using the DNS calculated moments of the scalar-dissipation PDF, it is shown, consistent with existing experimental information on single-phase flows, that the Gaussian is a modest approximation of the DNS-extracted PDF, particularly poor in the range of the high scalar-dissipation values, which are significant for turbulent reaction rate modeling in non-premixed flows using flamelet models. With the same DNS calculated moments of the scalar-dissipation PDF and making a change of variables, a model of this PDF is proposed in the form of the (beta)-PDF which is shown to approximate much better the DNS-extracted PDF, particularly in the regime of the high scalar-dissipation values. Several types of statistical measures are calculated over the ensemble of the fourteen databases. For each statistical measure, the proposed (beta)-PDF model is shown to be much superior to the Gaussian in approximating the DNS-extracted PDF. Additionally, the agreement between the DNS-extracted PDF and the (beta)-PDF even improves when the comparison is performed for higher initial Re layers, whereas the comparison with the Gaussian is independent of the initial Re values. For two-phase flows, the comparison between the DNS-extracted PDF and the (beta)-PDF also improves with increasing free-stream gas temperature and mass loading. The higher fidelity approximation of the DNS-extracted PDF by the (beta)-PDF with increasing Re, gas temperature and mass loading bodes well for turbulent reaction rate modeling.

scalar dissipation modeling↗

NASA NTRS2016

Numerical Simulation of HIWC Conditions with the Terminal Area Simulation System

Three-dimensional, numerical simulation of a mesoconvective system is conducted in order to better understand conditions associated with High Ice Water Content (HIWC) and its threat to aviation safety. Although peak local values of ice water content may occur early in the storm lifetime, large areas of high concentrations expand with time and persist even when the storm tops begin to warm. The storm canopy which contains HIWC, has low radar reflectivity factor and is fed by an ensemble of regenerating thermal pulses.

Proctor, Fred H.↗

NASA NTRS2001

Numerical Simulation of One- and Two-Phase Flows in Propulsion Systems

In this report, we present some results of problems investigated during joint research between the Hampton University Fluid Mechanics and Acoustics Laboratory (HU/FM&AL), NASA GRC, and the LaRC Hyper-X Program. This work is supported by joint research between the NASA GRC and the Institute of Mechanics at Moscow State University (IM/MSU) in Russia under a CRDF grant. The main areas of current scientific interest of the HU/FM&AL include an investigation of the proposed and patented advanced methods for aircraft engine thrust and noise benefits. These methods are based on nontraditional 3D corrugated and composite nozzle, inlet, propeller and screw designs such as a Bluebell and Telescope nozzle, Mobius-shaped screw, etc. This is the main subject of our other projects, of which one is presented at the current conference. Here we analyze additional methods for exhaust jet noise reduction without essential thrust loss and even with thrust augmentation. Such additional approaches are: (1) to add some solid, fluid, or gas mass at discrete locations to the main supersonic gas stream to minimize the negative influence of strong shock waves formed in propulsion systems. This mass addition may be accompanied by heat addition to the main stream as a result of the fuel combustion or by cooling of this stream as a result of the liquid mass evaporation and boiling; (2) Use of porous or permeable nozzles and additional shells at the nozzle exit for preliminary cooling of the hot jet exhaust and pressure compensation for off-design conditions (so-called continuous ejector with small mass flow rate); and (3) to propose and analyze new effective methods of fuel injection into the flow stream in air-breathing engines. The research is focused on a wide regime of problems in the propulsion field as well as in experimental testing and theoretical and numerical simulation analyses for advanced aircraft and rocket engines. The FM&AL Team uses analytical methods, numerical simulations, and experimental tests at the Hampton University campus, NASA, and IM/MSU. The main outcomes during this reporting period are: (1) Publications: The AIAA Paper #01-1893 has been accepted for the AIAA/NAL-NASDA-ISAS 10th International Space Planes and Hypersonic Systems and Technologies Conference, 24-27 April 2001, Kyoto, Japan. The AIAA Paper #01-3204 has been accepted for presentation at the 37th AIAA/ASME/SAE/ASEE Joint Propulsion Conference being held on 08-11 July, in Salt Lake City, UT; (2) Grants and proposals: The HU/FM&AL was awarded the NASA grant NAG3-2495 in October 2000 and the CRDF award was granted to the NASA GRC-HU FM&AL and IM/MSU (Russia) in July 2000. A solicited proposal was submitted for the NASA NRA-01GRC-02 competition and two unsolicited proposals to NASA are in preparation; (3) Theory and numerical simulations: Analytical theory, numerical simulation, comparison of theoretical with experimental results, and modification of theoretical approaches, models, grids, etc. Such investigations have been conducted for three main problems: (a) Combustion efficiency optimization in the half-duct combustor system; (b) Drag reduction effects for blunt bodies with solid needles; and (c) Solid particle injection from the butt-end against a supersonic flow. The NASA CFL3D, HU/FM&AL, and IM/MSU GODUNOV codes were used and modified for solution of these problems. The codes are based on full Euler and Navier-Stokes solvers with and without nonequilibrium oxygen-nitrogen and air-hydrogen chemical reactions in laminar and turbulent gas flow regimes; (4) Experimental tests: In the IM/MSU supersonic wind tunnel, experimental tests were conducted with different number, location, and geometric parameters of solid needles mounted at the front of the butt-end in supersonic flow. Optimal parameters were determined that provide minimal butt-end drag in the stationary flow regime. Experimental and numerical simulation results are in good agreement; (5) Student Research Activity: Involvement of one graduate and two undergraduate students as research assistants in the current project.

Gilinsky, Mikhail↗

NASA NTRS1982

Numerical simulation of magnetohydrodynamic shock propagation in the corona

Recent developments in the field of numerical simulation models for the study of shock wave propagation in the corona are presented. These models are based on gasdynamic (GD) and ideal (that is, dissipationless, except at shocks) magnetohydrodynamic (MHD) theories. The characteristics and physical interpretations of the results derived from these models are discussed in some detail. The most significant physical results obtained to date are provided by the two-dimensional non-planar, time-dependent, MHD numerical simulation model. In this model, the non-linear interaction among the three essential MHD waves, i.e., fast-, slow-, and Alfven waves are demonstrated. Finally, the physical relevance of these numerical simulation models in relation to observed solar activity is presented.

Wu, S. T.↗

NASA NTRS1993

Numerical simulation and physical aspects of supersonic vortex breakdown

Existing numerical simulations and physical aspects of subsonic and supersonic vortex-breakdown modes are reviewed. The solution to the problem of supersonic vortex breakdown is emphasized in this paper and carried out with the full Navier-Stokes equations for compressible flows. Numerical simulations of vortex-breakdown modes are presented in bounded and unbounded domains. The effects of different types of downstream-exit boundary conditions are studied and discussed.

Liu, C. H.↗

NASA NTRS2002

Numerical Simulation of Flow Field Within Parallel Plate Plastometer

Parallel Plate Plastometer (PPP) is a device commonly used for measuring the viscosity of high polymers at low rates of shear in the range 10 4 to 10 9 poises. This device is being validated for use in measuring the viscosity of liquid glasses at high temperatures having similar ranges for the viscosity values. PPP instrument consists of two similar parallel plates, both in the range of 1 inch in diameter with the upper plate being movable while the lower one is kept stationary. Load is applied to the upper plate by means of a beam connected to shaft attached to the upper plate. The viscosity of the fluid is deduced from measuring the variation of the plate separation, h, as a function of time when a specified fixed load is applied on the beam. Operating plate speeds measured with the PPP is usually in the range of 10 -3 cm/s or lower. The flow field within the PPP can be simulated using the equations of motion of fluid flow for this configuration. With flow speeds in the range quoted above the flow field between the two plates is certainly incompressible and laminar. Such flows can be easily simulated using numerical modeling with computational fluid dynamics (CFD) codes. We present below the mathematical model used to simulate this flow field and also the solutions obtained for the flow using a commercially available finite element CFD code.

Basil N. Antar↗

NASA NTRS2000

Recommendations for Achieving Accurate Numerical Simulation of Tip Clearance Flows in Transonic Compressor Rotors

The tip clearance flows of transonic compressor rotors are important because they have a significant impact on rotor and stage performance. While numerical simulations of these flows are quite sophisticated. they are seldom verified through rigorous comparisons of numerical and measured data because these kinds of measurements are rare in the detail necessary to be useful in high-speed machines. In this paper we compare measured tip clearance flow details (e.g. trajectory and radial extent) with corresponding data obtained from a numerical simulation. Recommendations for achieving accurate numerical simulation of tip clearance flows are presented based on this comparison. Laser Doppler Velocimeter (LDV) measurements acquired in a transonic compressor rotor, NASA Rotor 35, are used. The tip clearance flow field of this transonic rotor was simulated using a Navier-Stokes turbomachinery solver that incorporates an advanced k-epsilon turbulence model derived for flows that are not in local equilibrium. Comparison between measured and simulated results indicates that simulation accuracy is primarily dependent upon the ability of the numerical code to resolve important details of a wall-bounded shear layer formed by the relative motion between the over-tip leakage flow and the shroud wall. A simple method is presented for determining the strength of this shear layer.

VanZante, Dale E.↗

NASA NTRS2015

Numerical Simulations of Thermographic Responses in Composites

Numerical simulations of thermographic responses in composite materials have been a useful for evaluating and optimizing thermographic analysis techniques. Numerical solutions are particularly beneficial for thermographic techniques, since the fabrication of specimens with realistic flaws is difficult. Simulations are presented with different ply layups that incorporated the anisotropic thermal properties that exist in each ply. The results are compared to analytical series solutions and thermal measurements on composites with flat bottom holes and delaminations.

Winfree, William P.↗

NASA NTRS1978

The direct numerical simulations of the turbulent wakes of axisymmetric bodies

Results of direct numerical simulations of turbulence are compared with both laboratory data and self-similarity theory for the case of the turbulent wakes of towed, axisymmetric bodies. In general, the agreement of the simulation results with both the laboratory data and the self-similarity theory is good, although the comparisons are hampered by inadequate procedures for initializing the numerical simulations.

Riley, J. J.↗

NASA NTRS2003

Numerical Simulations of Non-Newtonian Convection in Ice: Application to Europa

Numerical simulations of solid state convection in Europa's ice shell have so far been limited to consideration of Newtonian flow laws, where the viscosity of ice is strongly dependent upon temperature, predicting that a stagnant lid should form at the top (10-40%) of a convecting ice shell. Such large thicknesses seem to contradict estimates of the effective elastic thickness of Europa s ice shell during its geologically active period. Recent laboratory experiments characterize the rheology of ice as the sum of contributions from several temperature and strain rate-dependent creep mechanisms. We present the results of numerical simulations of convection within Europa s ice shell using the finite-element model Citcom, applying the non-Newtonian rheology of grain boundry sliding. Our calculations suggest a shallower brittle/ductile transition and larger interior convective velocities compared to Newtonian rheology. The flow field is time-dependent, with small, localized upwellings and downwellings at the thermal boundary layers that have minimal topographic expression at the surface.

Barr, A. C.↗

NASA NTRS1977

Numerical simulation of spacecraft charging phenomena

A numerical simulation program is being constructed having the following features: (1) infinite circular cylindrical geometry with angle-dependence, (2) inclusion of incident particles, photoelectrons, secondary electrons, backscattered electrons, any gun emissions, and any internal current pathways including surface conductive layers, (3) quasistatic time-dependent iteration, in which sheath potential changes during particle transit times are ignored, (4) use of approximate, locally-dependent space charge density expressions in solving Poisson's equation for sheath potentials, with use of numerical orbit-following to determine surface currents, and (5) incident particle velocity distributions isotropic or beam-like, or some superposition of these. Rationales for each of these features are discussed.

Laframboise, J. G.↗

NASA NTRS1981

Numerical simulation of transonic flows in diffusers

Numerical simulations were made of two-dimensional transonic flows in diffusers, including flow separation induced by a shock or adverse pressure gradient. The mass-averaged, time-dependent, compressible Navier-Stokes equations, simplified by the thin-layer approximation, were solved using MacCormack's hybrid method. The eddy-viscosity formulation was described by the Wilcox-Rubesin's two-equation, k-omega model. Detailed comparison of the computed results with measurements showed good agreement in all cases, including one with massive separation induced by a strong shock. The computation correctly predicted the details of a distinct lambda shock pattern, closely duplicating the configuration observed experimentally in spark-schlieren photographs.

Liou, M.-S.↗

NASA NTRS1983

Numerical simulation of turbulent trailing edge flows

Numerical simulations of the time-dependent, Reynolds-averaged, Navier-Stokes equations, employing a two-equation turbulence model, are presented and compared with measurements from a series of trailing edge experiments at transonic Mach numbers. The test flows include an asymmetric flow with no separation, an asymmetric flow with a small region of separation and a symmetric flow with a large shock-wave induced separated zone. Comparisons are made for mean surface quantities as well as for mean and fluctuating flow-field quantities. For the trailing-edge flows with little or no separation, the solutions correctly predict all the major features of the flow field. Treatment of the viscous-inviscid interaction was found to be important for predicting these test cases. Two-equation eddy-viscosity turbulence models were found to be adequate for these flows. However, for the shock-wave induced separation case, these turbulence models were inadequate to predict this flow field. Modifications of the turbulence model to correct these deficiencies are discussed.

Horstman, C. C.↗

NASA NTRS1987

Numerical simulation of viscous transonic airfoil flows

Numerical simulations of transonic airfoil flows using the Reynolds-averaged Navier-Stokes equations and various turbulence models are presented and compared with experimental data. Three different airfoils were investigated under varying flow conditions ranging from subcritical unseparated flows to supercritical separated flows. The turbulence models investigated consisted of three zero-equation models and one two-equation model. For unseparated flows involving weak viscous-inviscid interactions, the four models were comparable in their agreement with experiment. For separated flows involving strong viscous-inviscid interactions, the nonequilibrium zero-equation model of Johnson and King gave the best overall agreement with experiment.

Coakley, Thomas J.↗

NASA NTRS1987

Numerical simulations of collisions between rotating particles

Numerical simulations of Keplerian systems consisting of 200 mutually colliding rotating particles show that friction and surface irregularity reduce the equilibrium velocity dispersion and transfer some of the energy of random velocities to rotational velocities. The results confirm the theoretical predictions of Salo (1987). Simulations are performed for identical and different particle sizes, and for a power-law distribution of sizes. Taking rotation into account is not found to alter the previously observed Rayleigh distribution of eccentricities and inclinations, while the components of the rotational velocities are found to follow a Gaussian distribution. Application to the rarefied regions of Saturn's rings suggests that friction is able to reduce the equilibrium geometric thickness by about one-half.

Salo, H.↗

NASA NTRS1987

Collisional interchange instability. I - Numerical simulations of intermediate-scale irregularities

Numerical simulations of the generalized Rayleigh-Taylor instability are presented. The model and simulations are applicable to bottomside and topside spread F, unstable barium cloud dynamics, and collisional interchange instability in general. The principal result is that the evolution of the instability tends to an anisotropic state consisting of nearly sinusoidal (quasi-periodic) variation along the effective electric field, and shocklike structures propagating perpendicular to E(eff) along the extrema of the quasi-periodic structures. The spectral properties of the nonlinear state are analyzed using one-dimensional power spectra calculated along spatial trajectories for selected angles to E(eff). In this way a direct comparison to in situ probe data can be made. The inherent anisotropy of the nonlinear state is reflected in major qualitative differences between the spectra taken parallel to and perpendicular to E(eff). The fundamental finding of the present work is that anisotropy in interchange dynamics is much greater than had been previously reported. This strong anisotropy can explain much of the spectral and spatial structural characteristics of both bottomside and topside spread F. In a companion paper a comparison of the simulation results to various in situ data sets is given.

Zargham, S.↗

NASA NTRS1989

Numerical simulations in combustion

This paper reviews numerical simulations in reacting flows in general and combustion phenomena in particular. It is shown that use of implicit schemes and/or adaptive mesh strategies can improve convergence, stability, and accuracy of the solution. Difficulties increase as turbulence and multidimensions are considered, particularly when finite-rate chemistry governs the given combustion problem. Particular attention is given to the areas of solid-propellant combustion dynamics, turbulent diffusion flames, and spray droplet vaporization.

Chung, T. J.↗

NASA NTRS1991

The use of direct numerical simulation data in turbulence modeling

Direct numerical simulations (DNS) of turbulent flows provide a complete data base to develop and to test turbulence models. In this article, the progress made in developing models for the dissipation rate equation is reviewed. New scaling arguments for the various terms in the dissipation rate equation were tested using data from DNS of homogeneous shear flows. Modifications to the epsilon-equation model that take into account near-wall effects were developed using DNS of turbulent channel flows. Testing of new models for flows under mean compression was carried out using data from DNS of isotropically compressed turbulence. In all of these studies the data from the simulations was essential in guiding the model development. The next generation of DNS will be at higher Reynolds numbers, and will undoubtedly lead to improved models for computations of flows of practical interest.

Mansour, N. N.↗