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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Numerical simulation of baroclinic Jovian vortices

We examine the evolution of baroclinic vortices in a time-dependent, nonlinear numerical model of a Jovian atmosphere. The model uses a normal-mode expansion in the vertical, using the barotropic and first two baroclinic modes. Results for the stability of baroclinic vortices on an f plane in the absence of a mean zonal flow are similar to results of Earth vortex models, although the presence of a fluid interior on the Jovian planets shifts the stability boundaries to smaller length scales. The presence of a barotropic mean zonal flow in the interior stabilizes vortices against instability and significantly modifies the finite amplitude form of baroclinic instabilities. The effect of a zonal flow on a form of barotropic instability produces periodic oscillations in the latitude and longitude of the vortex as observed at the level of the cloud tops. This instability may explain some, but not all, observations of longitudinal oscillations of vortices on the outer planets. Oscillations in aspect ratio and orientation of stable vortices in a zonal shear flow are observed in this baroclinic model, as in simpler two-dimensional models. Such oscillations are also observed in the atmospheres of Jupiter and Neptune. The meridional propagation and decay of vortices on a beta plane is inhibited by the presence of a mean zonal flow. The direction of propagation of a vortex relative to the mean zonal flow depends upon the sign of the meridional potential vorticity gradient; combined with observations of vortex drift rates, this may provide a constraint on model assumption for the flow in the deep interior of the Jovian planets.

Achterberg, Richard K.↗
NASA NTRS1999

Numerical Simulation of Thermocapillary Drop Motion with Internal Circulation

The thermocapillary motion of drops in zero gravity is analyzed numerically. When convective transport is important, the internal circulation in the drop has a profound effect on the temperature distribution in its vicinity and hence on its migration speed. For sufficiently large values of the Marangoni number number Ma, for steady motion of the drop, the temperature difference on the drop surface and its scaled speed increase with Ma. This is in contrast to (1) existing computational results for liquid drops whose scaled speed decreases with Ma and (2) asymptotic results for gas bubbles whose scaled speed is independent of Ma when it is large.

Ma, Xiangjiang↗
NASA NTRS1994

Numerical simulations of particle orbits around 2060 Chiron

Scattered light from orbiting or coorbiting dust is a primary signature by which Earth-based observers study the activity and atmosphere of the unusual outer solar system object 2060 Chiron. Therefore, it is important to understand the lifetime, dynamics, and loss rates of dust in its coma. We report here dynamical simulations of particles in Chiron's collisionless coma. The orbits of 17,920 dust particles were numerically integrated under the gravitational influence of Chiron, the Sun, and solar radiation pressure. These simulations show that particles ejected from Chiron are more likely to follow suborbital trajectories, or to escape altogether, than to enter quasistable orbits. Significant orbital lifetimes can only be achieved for very specific launch conditions. These results call into question models of a long-term, bound coma generated by discrete outbursts, and instead suggest that Chiron's coma state is closely coupled to the nearly instantaneous level of Chiron's surface activity.

Stern, S. A.↗
NASA NTRS1980

Numerical simulation of stratospheric sudden warmings with a primitive equation spectral model

A 26-level primitive equation spherical harmonic spectral model allowing for wave-wave and wave-zonal flow interactions is presented for the study of stratospheric sudden warmings. The warmings are simulated by the forcing of a single planetary wave at the tropopause. Four numerical experiments were performed. Nonlinear wave-wave interactions appear to play an important role in the evolution of the flow and temperature fields in the middle to upper stratosphere. In the case involving both wave-wave and wave-zonal flow interactions, this was manifested by the split in the initial polar vortex into a quasi-wave number 2 pattern. In the cases at 60 deg N, easterlies develop first in the upper mesosphere and descend gradually. About the same time or a little later, easterlies also develop in the mid-stratosphere. The linear cases exhibit warmings which are more shallow and more intense at 30 km than the nonlinear cases.

Lordi, N. J.↗
NASA NTRS1999

Numerical Simulations of Acoustically Driven, Burning Droplets

This computational study focuses on understanding and quantifying the effects of external acoustical perturbations on droplet combustion. A one-dimensional, axisymmetric representation of the essential diffusion and reaction processes occurring in the vicinity of the droplet stagnation point is used here in order to isolate the effects of the imposed acoustic disturbance. The simulation is performed using a third order accurate, essentially non-oscillatory (ENO) numerical scheme with a full methanol-air reaction mechanism. Consistent with recent microgravity and normal gravity combustion experiments, focus is placed on conditions where the droplet is situated at a velocity antinode in order for the droplet to experience the greatest effects of fluid mechanical straining of flame structures. The effects of imposed sound pressure level and frequency are explored here, and conditions leading to maximum burning rates are identified.

Kim, H.-C.↗
NASA NTRS1993

Numerical simulation of steady and unsteady viscous flow in turbomachinery using pressure based algorithm

The objective of this research is to simulate steady and unsteady viscous flows, including rotor/stator interaction and tip clearance effects in turbomachinery. The numerical formulation for steady flow developed here includes an efficient grid generation scheme, particularly suited to computational grids for the analysis of turbulent turbomachinery flows and tip clearance flows, and a semi-implicit, pressure-based computational fluid dynamics scheme that directly includes artificial dissipation, and is applicable to both viscous and inviscid flows. The values of these artificial dissipation is optimized to achieve accuracy and convergency in the solution. The numerical model is used to investigate the structure of tip clearance flows in a turbine nozzle. The structure of leakage flow is captured accurately, including blade-to-blade variation of all three velocity components, pitch and yaw angles, losses and blade static pressures in the tip clearance region. The simulation also includes evaluation of such quantities of leakage mass flow, vortex strength, losses, dominant leakage flow regions and the spanwise extent affected by the leakage flow. It is demonstrated, through optimization of grid size and artificial dissipation, that the tip clearance flow field can be captured accurately. The above numerical formulation was modified to incorporate time accurate solutions. An inner loop iteration scheme is used at each time step to account for the non-linear effects. The computation of unsteady flow through a flat plate cascade subjected to a transverse gust reveals that the choice of grid spacing and the amount of artificial dissipation is critical for accurate prediction of unsteady phenomena. The rotor-stator interaction problem is simulated by starting the computation upstream of the stator, and the upstream rotor wake is specified from the experimental data. The results show that the stator potential effects have appreciable influence on the upstream rotor wake. The predicted unsteady wake profiles are compared with the available experimental data and the agreement is good. The numerical results are interpreted to draw conclusions on the unsteady wake transport mechanism in the blade passage.

Lakshminarayana, B.↗
NASA NTRS1986

Numerical simulation of hypersonic, axisymmetric flowfields

Flow conditions corresponding to those of three previously conducted hypersonic axisymmetric flow experiments are presently treated by numerical solutions of the axisymmetric, thin layer Navier-Stokes equations with the Chemistry and Gasdynamic Implicit computer code's Version 2. While simulation and experimental results agree for ideal gas flows, the simulation for the case of a real gas flow reveals that the particular coupling of the gasdynamics and the chemistry used in the code degrades the convergence level obtainable, by comparison with ideal gas cases. The computed stagnation point heating is 60 percent lower than for the experiment.

Green, M. J.↗
NASA NTRS1991

Numerical simulation of shock-induced combustion/detonation in a premixed H2-air mixture using Navier-Stokes equations

A numerical study was conducted to address the structural stability of the oblique detonation wave. An attempt was made to isolate the structural instability predicted by Buckmaster's (1990) linear stability analysis. For this study a detailed viscous flowfield with finite-rate chemistry was computed past a 20-deg wedge. The finite-rate chemistry was modeled by using a seven-species and seven-reaction model. Instability modes were isolated by analyzing the Fourier power spectrum of the H2O mass fraction at selected sample stations. No unstable temporal modes with dominant frequency were isolated. It was concluded that the instability which is predicted by the Buckmaster work was a high-frequency low-amplitude phenomenon. For engineering purposes, the oblique detonation wave is a stable phenomenon as long as sufficient levels of overdrive are present.

Singh, D. J.↗
NASA NTRS1993

High resolution numerical simulation of the linearized Euler equations in conservation law form

A linearized Euler solver based on a high resolution numerical scheme is presented. The approach is to linearize the flux vector as opposed to carrying through the complete linearization analysis with the dependent variable vector written as a sum of the mean and the perturbed flow. This allows the linearized equations to be maintained in conservation law form. The linearized equations are used to compute unsteady flows in turbomachinery blade rows arising due to blade vibrations. Numerical solutions are compared to theoretical results (where available) and to numerical solutions of the nonlinear Euler equations.

Sreenivas, Kidambi↗
NASA NTRS1992

Numerical simulation of satellite-ring interactions - Resonances and satellite-ring torques

The Krook kinetic equation for planetary rings is numerically solved in two spatial dimensions and in time, with (1) interparticle collisions and (2) satellite-forcing, but (3) without self-gravity, for the case of a flattened planetary ring that undergoes gravitational perturbation by a nearby satellite. It is noted that the amplitude of wakes is limited by purely kinematic effects, even in the absence of collisions. Attention is given to the results of a simulation of an inner Lindblad-resonance location, as the distribution approaches steady state; these simulations do not show an increase in velocity dispersion in the resonance zone, obviating a net torque.

Brophy, Thomas G.↗
NASA NTRS1995

Full-Scale Direct Numerical Simulation of Two- and Three-Dimensional Instabilities and Rivulet Formulation in Heated Falling Films

A thin film draining on an inclined plate has been studied numerically using finite element method. Three-dimensional governing equations of continuity, momentum and energy with a moving boundary are integrated in an arbitrary Lagrangian Eulerian frame of reference. Kinematic equation is solved to precisely update interface location. Rivulet formation based on instability mechanism has been simulated using full-scale computation. Comparisons with long-wave theory are made to validate the numerical scheme. Detailed analysis of two- and three-dimensional nonlinear wave formation and spontaneous rupture forming rivulets under the influence of combined thermocapillary and surface-wave instabilities is performed.

Krishnamoorthy, S.↗
NASA NTRS1990

Numerical simulation of an F-16A at angle of attack

The transonic flowfield around an F-16A fighter configuration at a moderate incidence angle is simulated by solving the Navier-Stokes equations on a single-block grid. The numerical solution matches experimental freestream conditions with a mach number of 0.85, 16 degrees angle of attack, and a characteristic Reynolds number of 12.75 million. MacCormack's explicit algorithm is used in conjuction with a local time step and consecutive mesh refinement procedure to accelerate numerical convergence. The Baldwin-Lomax algebraic model provides turbulent closure. Computed surface pressure distributions and the aircraft lift coefficient compare favorably with wind tunnel data. The drag coefficient in the simulation overpredicts the experimental value by 8 percent.

Huband, G. W.↗
NASA NTRS2003

Numerical Simulation of Regional Circulation in the Monterey Bay Region

The objective of this study is to produce a high-resolution numerical model of Mon- terey Bay area in which the dynamics are determined by the complex geometry of the coastline, steep bathymetry, and the in uence of the water masses that constitute the CCS. Our goal is to simulate the regional-scale ocean response with realistic dynamics (annual cycle), forcing, and domain. In particular, we focus on non-hydrostatic e ects (by comparing the results of hydrostatic and non-hydrostatic models) and the role of complex geometry, i.e. the bay and submarine canyon, on the nearshore circulation. To the best of our knowledge, the current study is the rst to simulate the regional circulation in the vicinity of Monterey Bay using a non-hydrostatic model. Section 2 introduces the high resolution Monterey Bay area regional model (MBARM). Section 3 provides the results and veri cation with mooring and satellite data. Section 4 compares the results of hydrostatic and non-hydrostatic models.

Tseng, Y. H.↗
NASA NTRS1994

Numerical simulations of upstream blocking, columnar disturbances, and bores in stably stratified shear flows over an obstacle

A two-dimensional, nonhydrostatic, elastic numerical model has been used to study the generation of gravity waves for a stably stratified shear flow over an obstacle. When a low-level wind shear is included in the simulation, we find that the predictions for noticeable upstream effects based on Froude number for a uniform flow are no longer accurate. Upstream effects are encountered in the form of upstream propagating columnar disturbances and internal bores away from the obstacle. The limited parameter space studies conducted in this study suggest that the ratio of the shear depth to the obstacle heigh (d/H), the obstacle aspect ratio (H/L), and the Froude number (U/NH) are instrumental in determining the strength and the existence of these upstream disturbances. Thus, the present theoretical and empirical understanding of the importance of the Froude number for determining the nature of upstream effects should be modified substantially to include additional nondimensional parameters when shear is present.

Chen, Chaing↗
NASA NTRS1985

Numerical simulation of vortex breakdown

The breakdown of an isolated axisymmetric vortex embedded in an unbounded uniform flow is examined by numerical integration of the complete Navier-Stokes equations for unsteady axisymmetric flow. Results show that if the vortex strength is small, the solution approaches a steady flow and the vortex is stable. If the strength is large enough, the solution remains unsteady and a recirculating zone will appear near the axis, its form and internal structure resembling those of the axisymmetric breakdown bubbles with multi-cells observed by Faler and Leibovich (1978). For apppropriate combinations of flow parameters, the flow reveals quasi-periodicity. Parallel calculations with the quasi-cylindrical approximation indicate that so far as predicting of breakdown is concerned, its results coincide quite well with the results mentioned above. Both show that the vortex breakdown has little concern with the Reynolds number or with the critical classification of the upstream flow, at least for the lower range of Reynolds numbers.

Shi, X.↗
NASA NTRS1988

A 3-dimensional numerical simulation of the atmospheric injection of aerosols by a hypothetical basaltic fissure eruption

Researchers simulated the atmospheric response to a hypothetical basaltic fissure eruption using heating rates based on the Roza flow eruption. The simulation employs the Colorado State University Regional Atmospheric Model (RAMS) with scavenging effects. The numerical model is a three-dimensional non-hydrostatic time-split compressible cloud/mesoscale model. Explicit microphysics include prediction of cloud, rain, crystal, and hail precipitation types. Nucleation and phoretic scavenging are predicted assuming that the pollutant makes an effective cloud droplet nucleus. Smoke is carried as a passive tracer. Long and short wave radiation heating tendencies, including the effects of the smoke, are parameterized. The longwave emission by the lava surface is neglected in the parameterization and included as an explicit heating term instead. A regional scale domain of 100 x 100 km in the horizontal and 22 km high is used. The horizontal grid spacing is taken to be 2 km and the vertical spacing is taken to be 0.75 km. The initial atmospheric state is taken to be horizontally homogenous and based on the standard atmospheric sounding. The fissure is assumed to be 90 km long and oriented in a zig/zag pattern.

Tripoli, Gregory J.↗
NASA NTRS1990

Numerical simulation of spatially-evolving instability

A computational study of the spatial stability of plane Poiseuille flow is presented. The numerical scheme employs a time-splitting method to integrate the full Navier-Stokes equations using spectral collocation/finite-difference discretization on a non-staggered mesh. The eigenvalue decomposition procedure is applied for the solution of the Poisson equations using the capacitance matrix technique. The buffer domain method is incorporated for the outflow boundary conditions. The input perturbation velocities are obtained by solving the Orr-Sommerfeld equation for the nonlinear eigenvalue problem employing the companion matrix method. Computational results are compared with the linear theory for two-dimensional disturbances.

Danabasoglu, G.↗