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 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Numerical simulation of one-dimensional heat transfer in composite bodies with phase change

A numerical simulation was developed to investigate the one dimensional heat transfer occurring in a system composed of a layered aircraft blade having an ice deposit on its surface. The finite difference representation of the heat conduction equations was done using the Crank-Nicolson implicit finite difference formulation. The simulation considers uniform or time dependent heat sources, from heaters which can be either point sources or of finite thickness. For the ice water phase change, a numerical method which approximates the latent heat effect by a large heat capacity over a small temperature interval was applied. The simulation describes the temperature profiles within the various layers of the de-icer pad, as well as the movement of the ice water interface. The simulation could also be used to predict the one dimensional temperature profiles in any composite slab having different boundary conditions.

Dewitt, K. J.↗

NASA NTRS2018

Numerical Simulation of a Simplified High-Lift CRM Configuration Embedded with Fluidic Actuators

Numerical simulations have been performed for a simplified high-lift configuration that is representative of a modern transport airplane. This configuration includes a leading-edge slat, fuselage, wing, nacelle-pylon and a simple hinged flap. The suction surface of the flap is embedded with multiple rows of fluidic actuators to reduce the extent of reversed flow regions and improve the aerodynamic performance of the configuration with flap in a deployed state. In the current paper, a Lattice Boltzmann Method based high-fidelity computational fluid dynamics (CFD) code, known as PowerFLOW® is used to simulate the entire flow field associated with this configuration, including the flow inside the actuators. A fully compressible version of the PowerFLOW® code that has been validated for high speed flows is used for the present simulations to accurately represent the transonic flow regimes that are encountered in the flow field generated by the actuators operating at higher mass flow (momentum) rates required to mitigate reverse flow regions on the suction surfaces of the main wing and the flap. The numerical solutions predict the expected trends in aerodynamic forces as the actuation levels are increased. More efficient active flow control (AFC) systems and actuator arrangement for lift augmentation are emerging based on the parametric studies conducted here prior to wind tunnel tests. These numerical solutions will be compared with experimental data, once such data becomes available.

Vatsa, Veer N.↗

NASA NTRS1997

Direct Numerical Simulation of Incompressible Pipe Flow Using a B-Spline Spectral Method

A numerical method based on b-spline polynomials was developed to study incompressible flows in cylindrical geometries. A b-spline method has the advantages of possessing spectral accuracy and the flexibility of standard finite element methods. Using this method it was possible to ensure regularity of the solution near the origin, i.e. smoothness and boundedness. Because b-splines have compact support, it is also possible to remove b-splines near the center to alleviate the constraint placed on the time step by an overly fine grid. Using the natural periodicity in the azimuthal direction and approximating the streamwise direction as periodic, so-called time evolving flow, greatly reduced the cost and complexity of the computations. A direct numerical simulation of pipe flow was carried out using the method described above at a Reynolds number of 5600 based on diameter and bulk velocity. General knowledge of pipe flow and the availability of experimental measurements make pipe flow the ideal test case with which to validate the numerical method. Results indicated that high flatness levels of the radial component of velocity in the near wall region are physical; regions of high radial velocity were detected and appear to be related to high speed streaks in the boundary layer. Budgets of Reynolds stress transport equations showed close similarity with those of channel flow. However contrary to channel flow, the log layer of pipe flow is not homogeneous for the present Reynolds number. A topological method based on a classification of the invariants of the velocity gradient tensor was used. Plotting iso-surfaces of the discriminant of the invariants proved to be a good method for identifying vortical eddies in the flow field.

Loulou, Patrick↗

NASA NTRS2016

Numerical Simulation of a High-Lift Configuration Embedded with High Momentum Fluidic Actuators

Numerical simulations have been performed for a vertical tail configuration with deflected rudder. The suction surface of the main element of this configuration, just upstream of the hinge line, is embedded with an array of 32 fluidic actuators that produce oscillating sweeping jets. Such oscillating jets have been found to be very effective for flow control applications in the past. In the current paper, a high-fidelity computational fluid dynamics (CFD) code known as the PowerFLOW R code is used to simulate the entire flow field associated with this configuration, including the flow inside the actuators. A fully compressible version of the PowerFLOW R code valid for high speed flows is used for the present simulations to accurately represent the transonic flow regimes encountered in the flow field due to the actuators operating at higher mass flow (momentum) rates required to mitigate reverse flow regions on a highly-deflected rudder surface. The computed results for the surface pressure and integrated forces compare favorably with measured data. In addition, numerical solutions predict the correct trends in forces with active flow control compared to the no control case. The effect of varying the rudder deflection angle on integrated forces and surface pressures is also presented.

Vatsa, Veer N.↗

NASA NTRS1995

Numerical simulations of small amplitude acoustic wave propagation in a converging-diverging nozzle

Numerical simulations of the small amplitude acoustic wave propagation in a converging-diverging nozzle proposed in Category 5 are performed and presented here. The quasi 1-D unsteady flow equations in conservative form are discretized by the DRP scheme with the artificial damping developed by Tam and Webb. Characteristic boundary conditions are used at both the inlet and outlet of the nozzle. The effect of different numerical implementations of the subsonic inflow boundary conditions on the convergence of the solution to the steady state is studied. In the case of a subsonic outflow in which a shock is formed in the nozzle, the interaction between acoustics and the shock is also investigated.

Dong, Thomas Z.↗

NASA NTRS1982

Numerical simulation of wall-bounded turbulent shear flows

Developments in three dimensional, time dependent numerical simulation of turbulent flows bounded by a wall are reviewed. Both direct and large eddy simulation techniques are considered within the same computational framework. The computational spatial grid requirements as dictated by the known structure of turbulent boundary layers are presented. The numerical methods currently in use are reviewed and some of the features of these algorithms, including spatial differencing and accuracy, time advancement, and data management are discussed. A selection of the results of the recent calculations of turbulent channel flow, including the effects of system rotation and transpiration on the flow are included.

Moin, P.↗

NASA NTRS1982

Numerical simulation of wall-bounded turbulent shear flows

Moin, P.↗

NASA NTRS1970

Numerical simulation of the clustering of constant-volume balloons in the global domain

Numerical simulation of motion and clustering of balloons originally randomly spaced in upper and lower troposphere

Mesinger, F.↗

NASA NTRS1977

Numerical simulation of transition to turbulence

The purpose of the present study is the numerical simulation of transition to turbulence in a boundary layer. Starting from a three-dimensional disturbance, the evolution of the Navier-Stokes solution until breakdown is followed. Qualitative comparison with available experimental observations is carried out.

Wray, A. A.↗

NASA NTRS1981

Numerical simulation of a compressible homogeneous, turbulent shear flow

A direct, low Reynolds number, numerical simulation was performed on a homogeneous turbulent shear flow. The full compressible Navier-Stokes equations were used in a simulation on the ILLIAC IV computer with a 64,000 mesh. The flow fields generated by the code are used as an experimental data base, to examine the behavior of the Reynols stresses in this simple, compressible flow. The variation of the structure of the stresses and their dynamic equations as the character of the flow changed is emphasized. The structure of the tress tensor is more heavily dependent on the shear number and less on the fluctuating Mach number. The pressure-strain correlation tensor in the dynamic uations is directly calculated in this simulation. These correlations are decomposed into several parts, as contrasted with the traditional incompressible decomposition into two parts. The performance of existing models for the conventional terms is examined, and a model is proposed for the 'mean fluctuating' part.

Feiereisen, W. J.↗

NASA NTRS1984

Numerical Simulations of Convectively Excited Gravity Waves

Magneto-convection and gravity waves are numerically simulated with a nonlinear, three-dimensional, time-dependent model of a stratified, rotating, spherical fluid shell heated from below. A Solar-like reference state is specified while global velocity, magnetic field, and thermodynamic perturbations are computed from the anelastic magnetohydrodynamic equations. Convective overshooting from the upper (superadiabatic) part of the shell excites gravity waves in the lower (subadiabatic) part. Due to differential rotation and Coriolis forces, convective cell patterns propagate eastward with a latitudinally dependent phase velocity. The structure of the excited wave motions in the stable region is more time-dependent than of the convective motions above. The magnetic field tends to be concentrated over giant-cell downdrafts in the convective zone but is affected very little by the wave motion in the stable region.

Glatzmaier, G. A.↗

NASA NTRS1987

Numerical simulations of turbulent spots in plane Poiseuille and boundary-layer flow

Direct numerical simulations of turbulent spots in plane Poiseuille and boundary-layer flows are performed. Mature, self-similar spots are obtained. The propagation velocities and spreading angles are found to compare well with corresponding experiments. The difference in shape of the two spots is also clearly discernible: the turbulent parts are contained within arrowhead regions that point in opposite directions for the two cases. The wing-tip region of the Poiseuille spot is also found to consist of a large-amplitude semiturbulent wave packet.

Henningson, Dan↗

NASA NTRS1988

Turbulence spectrum of a passive temperature field - Results of a numerical simulation

The spectrum of a passive temperature field G(k) is determined via numerical simulation using three kinds of isotropic turbulent velocity fields. For a velocity field developed concurrently with the temperature field, G(k) collapses onto the spectrum obtained using a frozen developed velocity field only for sufficiently high conductivity values. For lower conductivity values, the power law behavior of G(k) is less clear.

Chasnov, J.↗

NASA NTRS1990

Numerical simulations of the rebound shock model for solar spicules

Using time-dependent numerical simulations, the proposed rebound shock mechanism for spicules has been examined. At temperatures above a critical value, T(c), the radiation is characteristic of the conditions in the optically thin corona and near optically thin transition region. When T less than T(c), the atmosphere has a radiative cooling time, tau(rad) characteristic of chromosphere. The spicule is initiated with a quasi-impulsive force in the low chromosphere, which drives a train of upward propagating rebound shocks along the rigid magnetic flux tube. These shocks then move the transition region upward. The material below the displaced transition region has temperatures and densities similar to those of spicules when T(c) = 20,000 K or more and tau(rad) = 500 s or more, but not when T(c) = 10,000 K, and probably not when tau(rad) = 100 s. For all the cases where the cross sectional area diverges rapidly with height, the upward velocity of the transition region is less than that of spicules. Moreover, the maximum height is less than that of average spicules. Taller, higher velocity spicules result when the magnetic field cross sectional area is constant. In all cases, the rebound shock mechanism produces substantial motions and temperature and density variations in chromospheric and transition region material. It is suggested that this may be a partial explanation for the continuous dynamic state of the lower solar atmosphere.

Sterling, Alphonse C.↗

NASA NTRS1992

Numerical simulation: Its contributions to turbulence modeling

The relative advantages and disadvantages of Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) are discussed. The topics are presented in viewgraph form.

Ferziger, Joel H.↗

NASA NTRS1993

A numerical simulation of galaxy subcluster mergers

We present preliminary results of a 3-D numerical simulation of two merging subclusters of galaxies. By self-consistently modelling the intracluster gas and dark matter dynamics, we hope to gain insight as to how the dynamics of both relate to such observables as the cluster x-ray emission, radio source morphology, and velocity dispersions.

Roettiger, Kurt↗

NASA NTRS1993

Direct numerical simulation of turbulent H2-O2 combustion using reduced chemistry

Results of direct numerical simulations of hydrogen-oxygen combustion using a partial-equilibrium chemistry scheme in constant density, decaying, isotropic turbulence are reported. The simulations qualitatively reproduce many features of experimental results, such as superequilibrium radical species mole fractions, with temperature and major species mole fractions closer to chemical equilibrium. It was also observed that the peak reaction rates occur in narrow zones where the stoichiometric surface intersects regions of high scalar dissipation, as might be expected for combustion conditions close to chemical equilibrium. Another finding was that high OH mole fraction correspond more closely to the stoichiometric surface than to areas of high reaction rate for conditions of the simulations. Simulation results were compared to predictions of the Conditional Moment Closure model. This model was found to give good results for all quantities of interest when the conditionally averaged scalar dissipation was used in the prediction. When the nonconditioned average dissipation was used, the predictions compared well to the simulations for most of the species and temperature, but not for the reaction rate. The comparison would be expected to improve for higher Reynolds number flows, however.

Montgomery, Christopher J.↗

NASA NTRS1999

Numerical Simulation of Low-Density Shock-Wave Interactions

Computational Fluid Dynamics (CFD) numerical simulations of low-density shock-wave interactions for an incident shock impinging on a cylinder have been performed. Flow-field density gradient and surface pressure and heating define the type of interference pattern and corresponding perturbations. The maximum pressure and heat transfer level and location for various interaction types (i.e., shock-wave incidence with respect to the cylinder) are presented. A time-accurate solution of the Type IV interference is employed to demonstrate the establishment and the steadiness of the low-density flow interaction.

Glass, Christopher E.↗