Search NASA⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

Direct numerical simulations of on-demand vortex generators: Mathematical formulation

The objective of the present research is the development and application of efficient adaptive numerical algorithms for the study, via direct numerical simulations, of active vortex generators. We are using innovative computational schemes to investigate flows past complex configurations undergoing arbitrary motions. Some of the questions we try to answer are: Can and how may we control the dynamics of the wake? What is the importance of body shape and motion in the active control of the flow? What is the effect of three-dimensionality in laboratory experiments? We are interested not only in coupling our results to ongoing, related experimental work, but furthermore to develop an extensive database relating the above mechanisms to the vortical wake structures with the long-range objective of developing feedback control mechanisms. This technology is very important to aircraft, ship, automotive, and other industries that require predictive capability for fluid mechanical problems. The results would have an impact in high angle of attack aerodynamics and help design ways to improve the efficiency of ships and submarines (maneuverability, vortex induced vibration, and noise).

Koumoutsakos, Petros↗

Interactive visualization of numerical simulation results: A tool for mission planning and data analysis

We report on the development of an interactive system for visualizing and analyzing numerical simulation results. This system is based on visualization modules which use the Application Visualization System (AVS) and the NCAR graphics packages. Examples from recent simulations are presented to illustrate how these modules can be used for displaying and manipulating simulation results to facilitate their comparison with phenomenological model results and observations.

Berchem, J.↗

Numerical Simulation of 3-D Supersonic Viscous Flow in an Experimental MHD Channel

The 3-D supersonic viscous flow in an experimental MHD channel has been numerically simulated. The experimental MHD channel is currently in operation at NASA Ames Research Center. The channel contains a nozzle section, a center section, and an accelerator section where magnetic and electric fields can be imposed on the flow. In recent tests, velocity increases of up to 40% have been achieved in the accelerator section. The flow in the channel is numerically computed using a new 3-D parabolized Navier-Stokes (PNS) algorithm that has been developed to efficiently compute MHD flows in the low magnetic Reynolds number regime. The MHD effects are modeled by introducing source terms into the PNS equations which can then be solved in a very e5uent manner. To account for upstream (elliptic) effects, the flowfield can be computed using multiple streamwise sweeps with an iterated PNS algorithm. The new algorithm has been used to compute two test cases that match the experimental conditions. In both cases, magnetic and electric fields are applied to the flow. The computed results are in good agreement with the available experimental data.

Kato, Hiromasa↗

Direct numerical simulation of the transitional zone

Properties of the transitional zone in channel and boundary-layer flow are determined from the results of direct numerical simulations of forced transition. Three channel flow cases produce similar behavior, particularly in the 'transition function' relevant for some algebraic eddy viscosity models of the transitional zone. Moreover, the transitional zone from the computed forced boundary-layer transition is qualitatively similar to experimental data for natural transition and to the computed channel flow transition.

Zang, Thomas A.↗

Numerical simulations of space processing - The reality, the myth, and the future

The present capabilities offered by the combination of advanced special and general purpose numerical codes and supercomputers are discussed. Emphasis is placed on unique demands placed by the space environment on the computation of realistic space experiments and recommendations for improving the scientific return on the space investment. It is concluded that numerical simulations can provide timely and relevant answers which can guide and improve future space experiments.

Chait, A.↗

Direct Numerical Simulation of Turbulent Condensation in Clouds

In this brief, we investigate the turbulent condensation of a population of droplets by means of a direct numerical simulation. To that end, a coupled Navier-Stokes/Lagrangian solver is used where each particle is tracked and its growth by water vapor condensation is monitored exactly. The main goals of the study are to find out whether turbulence broadens the droplet size distribution, as observed in in situ measurements. The second issue is to understand if and for how long a correlation between the droplet radius and the local supersaturation exists for the purpose of modeling sub-grid scale microphysics in cloud-resolving codes. This brief is organized as follows. In Section 2 the governing equations are presented, including the droplet condensation model. The implementation of the forcing procedure is described in Section 3. The simulation results are presented in Section 4 together with a sketch of a simple stochastic model for turbulent condensation. Conclusions and the main outcomes of the study are given in Section 5.

Shariff, K.↗

Time-Accurate Numerical Simulations of Synthetic Jet Quiescent Air

The unsteady evolution of three-dimensional synthetic jet into quiescent air is studied by time-accurate numerical simulations using a second-order accurate mixed explicit-implicit fractional step scheme on Cartesian grids. Both two-dimensional and three-dimensional calculations of synthetic jet are carried out at a Reynolds number (based on average velocity during the discharge phase of the cycle V(sub j), and jet width d) of 750 and Stokes number of 17.02. The results obtained are assessed against PIV and hotwire measurements provided for the NASA LaRC workshop on CFD validation of synthetic jets.

Rupesh, K-A. B.↗

Numerical Simulation Of Silicon-Ribbon Growth

Mathematical model includes nonlinear effects. In development simulates growth of silicon ribbon from melt. Takes account of entire temperature and stress history of ribbon. Numerical simulations performed with new model helps in search for temperature distribution, pulling speed, and other conditions favoring growth of wide, flat, relatively defect-free silicon ribbons for solar photovoltaic cells at economically attractive, high production rates. Also applicable to materials other than silicon.

Woda, Ben K.↗

Numerical simulation of shock wave diffraction by TVD schemes

An upwind total variation diminishing (TVD) scheme and a predictor-corrector symmetric TVD scheme were used to numerically simulate the blast wave diffraction on a stationary object. The objective is to help design an optimum configuration so that lateral motion is minimized and at the same time vortex shedding and flow separation are reduced during a blast wave encounter. Results are presented for a generic configuration for both a coarse grid and a fine grid to illustrate the global and local diffraction flow fields. Numerical experiments for the shock wave reflection on a wedge are also included to validate the current approach. Numerical study indicated that these TVD schemes are more stable and produced higher shock resolution than classical shock capturing methods such as the explicit MacCormack scheme.

Young, Victor Y. C.↗

Numerical Simulation of Flow Through an Artificial Heart

A solution procedure was developed that solves the unsteady, incompressible Navier-Stokes equations, and was used to numerically simulate viscous incompressible flow through a model of the Pennsylvania State artificial heart. The solution algorithm is based on the artificial compressibility method, and uses flux-difference splitting to upwind the convective terms; a line-relaxation scheme is used to solve the equations. The time-accuracy of the method is obtained by iteratively solving the equations at each physical time step. The artificial heart geometry involves a piston-type action with a moving solid wall. A single H-grid is fit inside the heart chamber. The grid is continuously compressed and expanded with a constant number of grid points to accommodate the moving piston. The computational domain ends at the valve openings where nonreflective boundary conditions based on the method of characteristics are applied. Although a number of simplifing assumptions were made regarding the geometry, the computational results agreed reasonably well with an experimental picture. The computer time requirements for this flow simulation, however, are quite extensive. Computational study of this type of geometry would benefit greatly from improvements in computer hardware speed and algorithm efficiency enhancements.

Rogers, Stuart E.↗

A time-accurate adaptive grid method and the numerical simulation of a shock-vortex interaction

A time accurate, general purpose, adaptive grid method is developed that is suitable for multidimensional steady and unsteady numerical simulations. The grid point movement is performed in a manner that generates smooth grids which resolve the severe solution gradients and the sharp transitions in the solution gradients. The temporal coupling of the adaptive grid and the PDE solver is performed with a grid prediction correction method that is simple to implement and ensures the time accuracy of the grid. Time accurate solutions of the 2-D Euler equations for an unsteady shock vortex interaction demonstrate the ability of the adaptive method to accurately adapt the grid to multiple solution features.

Bockelie, Michael J.↗

Numerical simulation of plasma processes driven by transverse ion heating

The plasma processes driven by transverse ion heating in a diverging flux tube are investigated with numerical simulation. The heating is found to drive a host of plasma processes, in addition to the well-known phenomenon of ion conics. The downward electric field near the reverse shock generates a doublestreaming situation consisting of two upflowing ion populations with different average flow velocities. The electric field in the reverse shock region is modulated by the ion-ion instability driven by the multistreaming ions. The oscillating fields in this region have the possibility of heating electrons. These results from the simulations are compared with results from a previous study based on a hydrodynamical model. Effects of spatial resolutions provided by simulations on the evolution of the plasma are discussed.

Singh, Nagendra↗

Numerical simulations of auroral plasma processes. I - Ion beams and conics. II - Electric fields

The characteristics of ion beams and conics as seen in two-dimensional numerical simulations in which the plasmas are driven by current sheets of a finite thickness are described. It is shown that the most energetic ions in the simulations have pitch angles near 90 deg, implying a large perpendicular acceleration of the ions. Large perpendicular electric fields, similar to those measured in electrostatic shocks, are mainly confined near the edges of the current sheets, where E-perpendicular is much larger than E-parallel. Double layers with E-perpendicular about equal to, and spatially separated from, E-parallel form in the interior of wide current sheets. Such double layers have upward electric fields in the region of the upward current. Downward electric fields develop outside the current sheet.

Singh, Nagendra↗

A numerical simulation of the magnetospheric gate model for the X-ray bursters

A Lagrangian, stellar interior computer code which assumes radiation transport by diffusion is used to numerically simulate an X-ray burst occurring on a rapidly rotating, 1 solar mass neutron star accreting material from a companion in noncontinuous, episodic fashion. The blob of gas from the companion extends from the surface of the neutron star to a radius of 185 km, is optically thick, has mass of 10 to the -11th solar masses, and is initially at rest with respect to the surface. The gas initially has a temperature of about 10 million K and a surface luminosity of 0.1 that of the sun. The infall results in a burst which lasts about 0.1 second and reaches a peak luminosity and effective temperature of 240,000 suns and nine million K. The burst is followed by a phase of oscillations with a period of 0.2 second. It is found that the burst is too cool and rapid to resemble normal bursts, although the time scale is in agreement with some observations.

Starrfield, S.↗

Numerical simulations of loops heated to solar flare temperatures. I - Gasdynamics. II - X-ray and UV spectroscopy

The NRL's Dynamic Flux Tube Model is used to numerically simulate the dynamic response of a coronal magnetic loop to an energy input of the order encountered in solar flares. The coronal plasma is heated by the deposition of flare energy at the top of the loop to more than 10 million K, yielding a conduction front that moves toward the chromosphere, where the plasma is heated by the large downward conductive flux and ablates upward to the coronal part of the loop at velocities of a few hundred km/sec. The conduction front simultaneously produces chromospheric ablation and compresses the material ahead of it. With the aid of compressional instabilities, the compressed plasma grows throughout the flare heating phase, presenting a possible source of the flare optical continuum emission which is correlated with soft X-ray radiation. The observational consequences of rapidly heated loop gas dynamic processes are discussed. In the second part of this presentation, the dynamical calculation results previously obtained are used to predict the spectral line intensities, profiles and wavelengths of several X-ray lines and the UV line of Fe XXI at 1354.1 A. Three different viewing orientations of the loop are considered.

Cheng, C.-C.↗

Numerical Simulation of Forced and Free-to-Roll Delta-Wing Motions

The three-dimensional, Reynolds-averaged, Navier-Stokes (RANS) equations are used to numerically simulate nonsteady vortical flow about a 65-deg sweep delta wing at 30-deg angle of attack. Two large-amplitude, high-rate, forced-roll motions, and a damped free-to-roll motion are presented. The free-to-roll motion is computed by coupling the time-dependent RANS equations to the flight dynamic equation of motion. The computed results are in good agreement with the forces, moments, and roll-angle time histories. Vortex breakdown is present in each case. Significant time lags in the vortex breakdown motions relative to the body motions strongly influence the dynamic forces and moments.

Chaderjian, Neal M.↗

Direct numerical simulation of hot jets

The ultimate motivation of this work is to investigate the stability of two dimensional heated jets and its implications for aerodynamic sound generation from data obtained with direct numerical simulations (DNS). As pointed out in our last report, these flows undergo two types of instabilities, convective or absolute, depending on their temperature. We also described the limits of earlier experimental and theoretical studies and explained why a numerical investigation could give us new insight into the physics of these instabilities. The aeroacoustical interest of these flows was also underlined. In order to reach this goal, we first need to succeed in the DNS of heated jets. Our past efforts have been focused on this issue which encountered several difficulties. Our numerical difficulties are directly related to the physical problem we want to investigate since these absolutely or almost absolutely unstable flows are by definition very sensitive to the smallest disturbances and are very likely to reach nonlinear saturation through a numerical feedback mechanism. As a result, it is very difficult to compute a steady laminar solution using a spatial DNS. A steady state was reached only for strongly co-flowed jets, but these flows are almost equivalent to two independent mixing layers. Thus they are far from absolute instability and have much lower growth rates.

Jacob, Marc C.↗

Analysis of subgrid-scale eddy viscosity with use of results from direct numerical simulations

Without resort to any modeling, subgrid-scale eddy viscosity is computed from the results of high-resolution (64 cubed and 128 cubed grid points) direct numerical simulations of three-dimensional homogeneous isotropic decaying turbulence. In these simulations the eddy viscosity peaks sharply at the cutoff wavenumber, in rough agreement with the results of Kraichnan (1976). In addition, in the low-wavenumber range the eddy viscosity may be negative, contrary to the generally accepted concept of subgrid-scale eddy viscosity. Some possible explanations of this behavior are discussed.

Domaradzki, J. Andrzej↗