SEARCH · Engineering Papers

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 991 records · Page 55

All sources NASA NTRS DOE OSTI

NASA NTRS1979

Numerical simulation of steady supersonic viscous flow

A noniterative, implicit, space-marching, finite-difference algorithm is developed for the steady thin-layer Navier-Stokes equations in conservation-law-form. The numerical algorithm is applicable to steady supersonic viscous flow over bodies of arbitrary shape. In addition, the same code can be used to compute supersonic inviscid flow or three-dimensional boundary layers. Computed results from two-dimensional and three-dimensional versions of the numerical algorithm are in good agreement with those obtained from more costly time-marching techniques.

Schiff, L. B.↗

NASA NTRS1990

A numerical simulation of the microburst phenomenon - Comparison of a discrete and a parameterized microphysical representation

Two versions of a microburst model that differ in how they represent microphysical processes are described and comparison of the numerical solutions is made. This study focuses on how the versions differ in representing the evaporation of raindrops and the differential sedimentation of drops that develops because of their fallspeeds. The two versions are then compared to demonstrate how predicted microburst characteristics are influenced by the choice of microphysical scheme.

Coen, Janice L.↗

NASA NTRS1991

Numerical simulation of swept-wing flows

Efforts of the last six months to computationally model the transition process characteristics of flow over swept wings are described. Specifically, the crossflow instability and crossflow/Tollmien-Schlichting wave interactions are analyzed through the numerical solution of the full 3D Navier-Stokes equations including unsteadiness, curvature, and sweep. This approach is chosen because of the complexity of the problem and because it appears that linear stability theory is insufficient to explain the discrepancies between different experiments and between theory and experiment. The leading edge region of a swept wing is considered in a 3D spatial simulation with random disturbances as the initial conditions.

Reed, Helen L.↗

NASA NTRS2001

High Order Numerical Simulation of Sound Generated by the Kirchhoff Vortex

An improved high order finite difference method for low Mach number computational aeroacoustics (CAA) is described. The improvements involve the conditioning of the Euler equations in perturbation form to minimize numerical cancellation error, and the use of a stable non-dissipative sixth-order central spatial differencing for the interior points and third-order at the boundary points. The spatial difference operator satisfies the summation-by-parts property to guarantee strict stability for linear hyperbolic systems. Spurious high frequency oscillations are damped by a third-order characteristic-based filter. The objective of this paper is to apply these improvements in the simulation of sound generated by the Kirchhoff vortex.

Mueller, Bernhard↗

NASA NTRS1991

Numerical simulations of thermal instabilities in stratified gases. II - Exploration of the parameter space

The temporal evolution of density perturbations in an initially hydrostatic isothermal atmosphere consisting of an optically thin radiating compressible plasma is studied. Numerical techniques are used to describe the nonlinear evolution of the perturbations, and the relative equilibrium between dynamic and thermal instabilities as governed by three independent control parameters are examined, namely, the initial density contrast of the perturbation, the ratio of the local buoyancy oscillation period to the local radiative cooling time, and the ratio of the perturbation radius to the local scaleheight. Four orders of magnitude of initial density contrasts and ratios of buoyancy and cooling times, and one order of magnitude of the bubble dimensions are explored. Well-defined oscillations were found to occur in a limited parameter range, and thermal instability to occur even within secondary condensations deriving from the bubble fragmentation.

Reale, F.↗

DOE OSTI2025

Numerical simulations of three-dimensional ion crystal dynamics in a Penning trap using the fast multipole method

We simulate the dynamics, including laser cooling, of three-dimensional (3-D) ion crystals confined in a Penning trap using a newly developed molecular dynamics-like code. The numerical integration of the ions’ equations of motion is accelerated using the fast multipole method to calculate the Coulomb interaction between ions, which allows us to efficiently study large ion crystals with thousands of ions. In particular, we show that the simulation time scales linearly with ion number, rather than with the square of the ion number. By treating the ions’ absorption of photons as a Poisson process, we simulate individual photon scattering events to study laser cooling of 3-D ellipsoidal ion crystals. Initial simulations suggest that these crystals can be efficiently cooled to ultracold temperatures, aided by the mixing of the easily cooled axial motional modes with the low frequency planar modes. In our simulations of a spherical crystal of 1000 ions, the planar kinetic energy is cooled to several millikelvin in a few milliseconds while the axial kinetic energy and total potential energy are cooled even further. This suggests that 3-D ion crystals could be well suited as platforms for future quantum science experiments.

Zaris, John (ORCID:0009000196476323)↗

NASA NTRS1994

Direct numerical simulation of sheared turbulent flow

The summer assignment to study sheared turbulent flow was divided into three phases which were: (1) literature survey, (2) computational familiarization, and (3) pilot computational studies. The governing equations of fluid dynamics or Navier-Stokes equations describe the velocity, pressure, and density as functions of position and time. In principle, when combined with conservation equations for mass, energy, and thermodynamic state of the fluid a determinate system could be obtained. In practice the Navier-Stokes equations have not been solved due to the nonlinear nature and complexity of these equations. Consequently, the importance of experiments in gaining insight for understanding the physics of the problem has been an ongoing process. Reasonable computer simulations of the problem have occured as the computational speed and storage of computers has evolved. The importance of the microstructure of the turbulence dictates the need for high resolution grids in extracting solutions which contain the physical mechanisms which are essential to a successful simulation. The recognized breakthrough occurred as a result of the pioneering work of Orzag and Patterson in which the Navier-Stokes equations were solved numerically utilizing a time saving toggling technique between physical and wave space, known as a spectral method. An equally analytically unsolvable problem, containing the same quasi-chaotic nature as turbulence, is known as the three body problem which was studied computationally as a first step this summer. This study was followed by computations of a two dimensional (2D) free shear layer.

Harris, Vascar G.↗

NASA NTRS1989

Numerical simulation of separated and vortical flows on bodies at large angles of attack

Navier-Stokes computations were carried out for high-angle-of-attack flow about bodies of revolution in supersonic and subsonic flow. Numerical and physical factors affecting flow predictions are discussed. Accurate prediction of the complex three-dimensional separated flow is found to depend on using grids sufficiently fine to resolve the details of both the viscous boundary layer and the off-surface separated flow structures, and, for turbulent flow cases, use of an eddy-viscosity turbulence model which accounts for the leeward vortical flow structures. An example of high-incidence flow computed about a more complex geometry, the F-18 fighter fuselage forebody and leading edge extension, is also presented.

Schiff, Lewis B.↗

NASA NTRS1991

Numerical simulation of spatially-evolving instability in plane channel flow

The spatial stability of plane channel flow is analyzed using a three-dimensional, time-dependent spectral/finite difference code (Danabasoglu et al., 1990) which integrates numerically the Navier-Stokes equations. The study centers on inflow disturbance amplitudes effects on the secondary instability. The resolution requirements along the spacewise direction, which become critical before the breakdown stage, are of particular interest. A direct comparison is made with the experiments of Nishioka et al. (1980).

Danabasoglu, G.↗

NASA NTRS1991

Flight test and numerical simulation of transonic flow around YAV-8B Harrier II wing

A computational fluid dynamics (CFD) method is used to study the aerodynamics of the YAV-8B Harrier II wing in the transonic region. A numerical procedure is developed to compute the flow field around the complicated wing-pylon-fairing geometry. The surface definition of the wing and pylons were obtained from direct measurement using theodolite triangulation. A thin-layer Navier-Stokes code with the Chimera technique is used to compute flow solutions. The computed pressure distributions at several span stations are compared with flight test data and show good agreement. Computed results are correlated with flight test data that show the flow is severely separated in the vicinity of the wing-pylon junction. Analysis shows that shock waves are induced by pylon swaybrace fairings, that the flow separation is much stronger at the outboard pylon and that the separation is caused mainly by the crossflow passing the geometry of wing-pylon junction.

Gea, Lie-Mine↗

NASA NTRS1988

Numerical simulation of heat transfer during the crystal growth of HgCdTe alloys

The effects of growth parameters on the thermal distribution during a Bridgman-Stockbarger-type crystal growth of HgCdTe alloys using short gradient zones are analyzed numerically. The analysis takes into account the change in the thermophysical properties upon freezing and considers both the temperature and composition dependences of the properties as well as translation rate effects. The calculated results compare favorably with previously published empirical results.

Dakhoul, Y. M.↗

NASA NTRS1988

Primitive numerical simulation of circular Couette flow - Carrousel wind tunnel nonturbulent solutions

The azimuthal-invariant, three-dimensional cylindrical, incompressible Navier-Stokes equations are solved to steady state for a finite-length, physically realistic model. The numerical method relies on an alternating-direction implicit scheme that is formally second-order accurate in space and first-order accurate in time. The equations are linearized and uncoupled by evaluating variable coefficients at the previous time iteration. Wall grid clustering is provided by a Roberts transformation in radial and axial directions. A vorticity-velocity formulation is found to be preferable to a vorticity-streamfunction approach. Subject to no-slip, Dirichlet boundary conditions, except for the inner cylinder rotation velocity (impulsive start-up) and zero-flow initial conditions, nonturbulent solutions are obtained for sub- and supercritical Reynolds numbers of 100 to 400 for a finite geometry where R(outer)/R(inner) = 1.5, H/R(inner) = 0.73, and H/Delta-R = 1.5. An axially-stretched model solution is shown to asymptotically approach the one-dimensional analytic Couette solution at the cylinder midheight. Flowfield change from laminar to Taylor-vortex flow is discussed as a function of Reynolds number. Three-dimensional velocities, vorticity, and streamfunction are presented via two-dimensional graphs and three-dimensional surface and contour plots.

Hasiuk, Jan↗

NASA NTRS1993

Three-dimensional numerical simulation of gradual opening in a wave rotor passage

The evolution of the contact interface and the propagation of compression waves inside a single wave rotor passage gradually opening to and traversing an inlet port is studied numerically using an inviscid formulation of the governing equations. Insights into the response of the interface and kinematics of the flow field to various opening times are given. Since the opening time is inversely proportional to the rotational speed of the rotor, the effects of passage rotation such as centripetal and Coriolis accelerations are intrinsically coupled to the gradual opening process. Certain three-dimensional features associated with the gradual opening process as a result of centripetal and Coriolis accelerations are illustrated. For the range of opening times or rotational speeds considered, a portion of the interface behaves like a vortex sheet that can degenerate into a complex interfacial structure. The vortices produced along the interface can serve as a stirring mechanism to promote local mixing. Coriolis and centripetal accelerations can introduce three dimensional effects such as interfacial distortions in meridional planes and spanwise migration of fluid elements.

Larosiliere, Louis M.↗

NASA NTRS1993

Three-dimensional numerical simulation of gradual opening in a wave rotor passage

Larosiliere, Louis M.↗

DOE OSTI2013

Accurate Numerical Simulations Of Chemical Phenomena Involved in Energy Production and Storage with MADNESS and MPQC: ALCF-2 Early Science Program Technical Report

In order to solve the electronic structure of large molecular systems on petascale computers using MADNESS, a numerical tool kit, are required fast and accurate implementations for linear algebra. MADNESS uses multiresolution analysis and low separation rank which translates high dimensional functions in tensor products using Legendre polynomial. The multiple tensor products make to the singular value decomposition and matrix multiplication the most intense operations in MADNESS. This work discusses the interfacing of Eigen3 as a C++ substitute of LAPACK and introduces Elemental for the diagonalization of large matrices. Furthermore, the present paper shows the performance these libraries on Blue Gene/ Q.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2001

Numerical Simulations of Precipitation Processes, Microphysics, and Microwave Radiative Properties of flood Producing Storms in Mediterranean & Adriatic Basins

A comprehensive understanding of the meteorological and microphysical nature of Mediterranean storms requires a combination of in situ data analysis, radar data analysis, and satellite data analysis, effectively integrated with numerical modeling studies at various scales. An important aspect of understanding microphysical controls of severe storms, is first understanding the meteorological controls under which a storm has evolved, and then using that information to help characterize the dominant microphysical processes. For hazardous Mediterranean storms, highlighted by the October 5-6, 1998 Friuli flood event in northern Italy, a comprehensive microphysical interpretation requires an understanding of the multiple phases of storm evolution. This involves intense convective development, Sratiform decay, orographic lifting, and sloped frontal lifting processes, as well as the associated vertical motions and thermodynamical instabilities governing physical processes that effect details of the size distributions and fall rates of the various types of hydrometeors found within the storm environment. This talk overviews the microphysical elements of a severe Mediterranean storm in such a context, investigated with the aid of TRMM satellite and other remote sensing measurements, but guided by a nonhydrostatic mesoscale model simulation of the Friuli flood event. The data analysis for this paper was conducted by my research groups at the Global Hydrology and Climate Center in Huntsville, AL and Florida State University in Tallahassee, and in collaboration with Dr. Alberto Mugnai's research group at the Institute of Atmospheric Physics in Rome. The numerical modeling was conducted by Professor Oreg Tripoli and Ms. Giulia Panegrossi at the University of Wisconsin in Madison, using Professor Tripoli's nonhydrostatic modeling system (NMS). This is a scalable, fully nested mesoscale model capable of resolving nonhydrostatic circulations from regional scale down to cloud scale and below.

Smith, Eric A.↗

NASA NTRS1994

Numerical simulation of the flow about the F-18 HARV at high angle of attack

This report summarizes research done over the past two years as part of NASA Grant NCC 2-729. This research has been aimed at validating numerical methods for computing the flow about the complete F-18 HARV at alpha = 30 deg and alpha = 45 deg. At 30 deg angle of attack, the flow about the F-18 is dominated by the formation, and subsequent breakdown, of strong vortices over the wing leading-edge extensions (LEX). As the angle of attack is increased to alpha = 45 deg, the fuselage forebody of the F-18 contains significant laminar and transitional regions which are not present at alpha = 30 deg. Further, the flow over the LEX at alpha = 45 deg is dominated by an unsteady shedding in time, rather than strong coherent vortices. This complex physics, combined with the complex geometry of a full aircraft configuration, provides a challenge for current computational fluid dynamics (CFD) techniques. The following sections present the numerical method and grid generation scheme that was used, a review of prior research done to numerically model the F-18 HARV, and a discussion of the current research. The current research is broken into two main topics: the effect of engine-inlet mass-flow rate on the F-18 vortex breakdown position, and the results using a refined F-18 computational model to compute the flow at alpha = 30 deg and alpha = 45 deg.

Murman, Scott M.↗

NASA NTRS1996

Numerical simulation of film-cooled ablative rocket nozzles

The objective of this research effort was to evaluate the impact of incorporating an additional cooling port downstream between the injector and nozzle throat in the NASA Fast Track chamber. A numerical model of the chamber was developed for the analysis. The analysis did not model ablation but instead correlated the initial ablation rate with the initial nozzle wall temperature distribution. The results of this study provide guidance in the development of a potentially lighter, second generation ablative rocket nozzle which maintains desired performance levels.

Landrum, D. B.↗