Engineering Papers⌕ Search

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.

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

At least 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Higher-Order Compact Schemes for Numerical Simulation of Incompressible Flows

A higher order accurate numerical procedure has been developed for solving incompressible Navier-Stokes equations for 2D or 3D fluid flow problems. It is based on low-storage Runge-Kutta schemes for temporal discretization and fourth and sixth order compact finite-difference schemes for spatial discretization. The particular difficulty of satisfying the divergence-free velocity field required in incompressible fluid flow is resolved by solving a Poisson equation for pressure. It is demonstrated that for consistent global accuracy, it is necessary to employ the same order of accuracy in the discretization of the Poisson equation. Special care is also required to achieve the formal temporal accuracy of the Runge-Kutta schemes. The accuracy of the present procedure is demonstrated by application to several pertinent benchmark problems.

Wilson, Robert V.↗
NASA NTRS1999

Implicit and Multigrid Method for Ideal Multigrid Convergence: Direct Numerical Simulation of Separated Flow Around NACA 0012 Airfoil

Numerical investigation of flow separation over a NACA 0012 airfoil at large angles of attack has been carried out. The numerical calculation is performed by solving the full Navier-Stokes equations in generalized curvilinear coordinates. The second-order LU-SGS implicit scheme is applied for time integration. This scheme requires no tridiagonal inversion and is capable of being completely vectorized, provided the corresponding Jacobian matrices are properly selected. A fourth-order centered compact scheme is used for spatial derivatives. In order to reduce numerical oscillation, a sixth-order implicit filter is employed. Non-reflecting boundary conditions are imposed at the far-field and outlet boundaries to avoid possible non-physical wave reflection. Complex flow separation and vortex shedding phenomenon have been observed and discussed.

Liu, Chao-Qun↗
NASA NTRS1986

Numerical simulation of boundary-layer excitation by surface heating/cooling

The concept of active control of growing disturbances in an unstable compressible flow by using time periodic, localized surface heating is studied numerically. The simulations are calculated by a fourth-order accurate solution of the compressible, laminar Navier-Stokes equations. Fourth-order accuracy is particularly important for this problem because the solution must be computed over many wavelengths. The numerical results demonstrate the growth of an initially small fluctuation into the nonlinear regime where a local breakdown into smaller scale disturbances can be observed. It is shown that periodic surface heating over a small strip can reduce the level of the fluctuation provided that the phase of the heating current is properly chosen.

Bayliss, A.↗
DOE OSTI2013

A new model for two-dimensional numerical simulation of pseudo-2D gas-solids fluidized beds

Pseudo-two dimensional (pseudo-2D) fluidized beds, for which the thickness of the system is much smaller than the other two dimensions, is widely used to perform fundamental studies on bubble behavior, solids mixing, or clustering phenomenon in different gas-solids fluidization systems. The abundant data from such experimental systems are very useful for numerical model development and validation. However, it has been reported that two-dimensional (2D) computational fluid dynamic (CFD) simulations of pseudo-2D gas-solids fluidized beds usually predict poor quantitative agreement with the experimental data, especially for the solids velocity field. In this paper, a new model is proposed to improve the 2D numerical simulations of pseudo-2D gas-solids fluidized beds by properly accounting for the frictional effect of the front and back walls. Two previously reported pseudo-2D experimental systems were simulated with this model. Compared to the traditional 2D simulations, significant improvements in the numerical predictions have been observed and the predicted results are in better agreement with the available experimental data.

42 ENGINEERING↗
DOE OSTI2014

A numerical simulation of solar energetic particle dropouts during impulsive events

This paper investigates the conditions for producing rapid variations of solar energetic particle (SEP) intensity commonly known as 'dropouts'. In particular, we use numerical model simulations based on solving the focused transport equation in the three-dimensional Parker interplanetary magnetic field to put constraints on the properties of particle transport coefficients in both directions perpendicular and parallel to the magnetic field. Our calculations of the temporal intensity profile of 0.5 and 5 MeV protons at the Earth show that the perpendicular diffusion must be small while the parallel mean free path is long in order to reproduce the phenomenon of SEP dropouts. When the parallel mean free path is a fraction of 1 AU and the observer is located at 1 AU, the perpendicular to parallel diffusion ratio must be below 10{sup –5} if we want to see the particle flux dropping by at least several times within 3 hr. When the observer is located at a larger solar radial distance, the perpendicular to parallel diffusion ratio for reproducing the dropouts should be even lower than that in the case of 1 AU distance. A shorter parallel mean free path or a larger radial distance from the source to observer will cause the particles to arrive later, making the effects of perpendicular diffusion more prominent and SEP dropouts disappear. All of these effects require the magnetic turbulence that resonates with the particles to be low everywhere in the inner heliosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

Tip vortices of wings in subsonic and transonic flow: A numerical simulation

Thin layer Navier-Stokes and Euler equations are numerically solved using a multi-block zonal approach to simulate the formation and roll up of tip vortices of wings in subsonic and transonic flows. Several wing planforms were considered to examine the influence of tip-cap shape, planform geometry and free stream Mach number on the formation process. A good definition of the formation and qualitative roll up of tip vortices was achieved.

Srinivasan, G. R.↗
NASA NTRS1985

Three-dimensional tomographic reconstruction of the flow around a revolving helicopter rotorblade - A numerical simulation

This paper discusses tomographic data acquisition and reconstruction procedures for determining the flow around a revolving helicopter rotor blade. The data acquisition procedure is simulated using, as a testbed, a numerical representation of the flow and projections through this known field are computed. An investigation has been made of: (1) the effect of the number of projections on the accuracy of the reconstruction and (2) the fidelity of the reconstruction when the region of intersection of all projections is smaller than the spatial extent of the flow.

Hesselink, L.↗
NASA NTRS1986

Numerical simulations of comets - Predictions for Comet Giacobini-Zinner

Simulations of Comet Giacobini-Zinner's interaction with solar wind are described and results are presented. The simulations are carried out via the numerical solution of the ideal MHD equations as an initial value problem in a uniform solar wind. The calculations are performed on a Cartesian mesh centered at the comet. Results reveal that the first significant modifications of the solar wind along the ISEE/ICE trajectory will occur 100,000 km from the solar wind comet axis.

Fedder, J. A.↗
NASA NTRS1991

Numerical simulation of the helium gas spin-up channel performance of the relativity gyroscope

The dependence of the spin-up system efficiency on each geometrical parameter of the spin-up channel and the exhaust passage of the Gravity Probe-B (GPB) is individually investigated. The spin-up model is coded into a computer program which simulates the spin-up process. Numerical results reveal optimal combinations of the geometrical parameters for the ultimate spin-up performance. Comparisons are also made between the numerical results and experimental data. The experimental leakage rate can only be reached when the gap between the channel lip and the rotor surface increases beyond physical limit. The computed rotating frequency is roughly twice as high as the measured ones although the spin-up torques fairly match.

Karr, Gerald R.↗
NASA NTRS1997

Numerical Simulation of Turbulent Jets with Rectangular Cross-Section

Three-dimensional turbulent jets with rectangular cross-section are simulated with a finite-difference numerical method. The full Navier-Stokes equations are solved at low Reynolds numbers, whereas at the high Reynolds numbers filtered forms of the equations are solved along with a sub-grid scale model to approximate effects of the unresolved scales. A 2-N storage, third-order Runge-Kutta scheme is used for temporal discretization and a fourth-order compact scheme is used for spatial discretization. Computations are performed for different inlet conditions which represent different types of jet forcing. The phenomenon of axis-switching is observed, and it is confirmed that this is based on self-induction of the vorticity field. Budgets of the mean streamwise velocity show that convection is balanced by gradients of the Reynolds stresses and the pressure.

Wilson, Robert V.↗
DOE OSTI2013

Direct numerical simulation of electrokinetic instability and transition to chaotic motion

A new type of instability—electrokinetic instability—and an unusual transition to chaotic motion near a charge-selective surface (semiselective electric membrane, electrode, or system of micro-/nanochannels) was studied by the numerical integration of the Nernst-Planck-Poisson-Stokes system and a weakly nonlinear analysis near the threshold of instability. A special finite-difference method was used for the space discretization along with a semi-implicit 31/3 -step Runge-Kutta scheme for the integration in time. Two kinds of initial conditions were considered: (a) white-noise initial conditions to mimic “room disturbances” and subsequent natural evolution of the solution, and (b) an artificial monochromatic ion distribution with a fixed wave number to simulate regular wave patterns. The results were studied from the viewpoint of hydrodynamic stability and bifurcation theory. The threshold of electroconvective movement was found by the linear spectral stability theory, the results of which were confirmed by numerical simulation of the entire system. Our weakly nonlinear analysis and numerical integration of the entire system predict possibility of both kinds of bifurcations at the critical point, supercritical and subcritical, depending on the system parameters. The following regimes, which replace each other as the potential drop between the selective surfaces increases, were obtained: one-dimensional steady solution, two-dimensional steady electroconvective vortices (stationary point in a proper phase space), unsteady vortices aperiodically changing their parameters (homoclinic contour), periodic motion (limit cycle), and chaotic motion. The transition to chaotic motion does not include Hopf bifurcation. The numerical resolution of the thin concentration polarization layer showed spike-like charge profiles along the surface, which could be, depending on the regime, either steady or aperiodically coalescent. The numerical investigation confirmed the experimentally observed absence of regular (near-sinusoidal) oscillations for the overlimiting regimes. There is a qualitative agreement of the experimental and the theoretical values of the threshold of instability, the dominant size of the observed coherent structures, and the experimental and theoretical volt–current characteristics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1992

Numerical simulation of the circulation of the atmosphere of Titan

A three dimensional General Circulation Model (GCM) of Titan's atmosphere is described. Initial results obtained with an economical two dimensional (2D) axisymmetric version of the model presented a strong superrotation in the upper stratosphere. Because of this result, a more general numerical study of superrotation was started with a somewhat different version of the GCM. It appears that for a slowly rotating planet which strongly absorbs solar radiation, circulation is dominated by global equator to pole Hadley circulation and strong superrotation. The theoretical study of this superrotation is discussed. It is also shown that 2D simulations systemically lead to instabilities which make 2D models poorly adapted to numerical simulation of Titan's (or Venus) atmosphere.

Hourdin, F.↗
NASA NTRS1996

Numerical Simulation of Wake Vortices Measured During the Idaho Falls and Memphis Field Programs

A numerical large-eddy simulation model is under modification and testing for application to aircraft wake vortices. The model, having a meteorological framework, permits the interaction of wake vortices with environments characterized by crosswind shear, stratification, and humidity. As part of the validation process, model results are compared with measured field data from the 1990 Idaho Falls and the 1994-1995 Memphis field experiments. Cases are selected that represent different aircraft and a cross section of meteorological environments. Also included is one case with wake vortex generation in ground effect. The model simulations are initialized with the appropriate meteorological conditions and a post roll-up vortex system. No ambient turbulence is assumed in our initial set of experiments, although turbulence can be self generated by the interaction of the model wakes with the ground and environment.

Proctor, Fred H.↗
DOE OSTI2016

Comprehensive numerical methodology for direct numerical simulations of compressible Rayleigh-Taylor instability

A comprehensive numerical methodology has been developed that handles the challenges introduced by considering the compressive nature of Rayleigh-Taylor instability (RTI) systems, which include sharp interfacial density gradients on strongly stratified background states, acoustic wave generation and removal at computational boundaries, and stratification-dependent vorticity production. The computational framework is used to simulate two-dimensional single-mode RTI to extreme late-times for a wide range of flow compressibility and variable density effects. The results show that flow compressibility acts to reduce the growth of RTI for low Atwood numbers, as predicted from linear stability analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2012

SO2 over Central China: Measurements, Numerical Simulations and the Tropospheric Sulfur Budget

SO2 in central China was measured in situ from an aircraft and remotely using the Ozone Monitoring Instrument (OMI) from the Aura satellite; results were used to develop a numerical tool for evaluating the tropospheric sulfur budget - sources, sinks, transformation and transport. In April 2008, measured ambient SO2 concentrations decreased from approx.7 ppbv near the surface to approx. 1 ppbv at 1800 m altitude (an effective scale height of approx.800 m), but distinct SO2 plumes were observed between 1800 and 4500 m, the aircraft's ceiling. These free tropospheric plumes play a major role in the export of SO2 and in the accuracy of OMI retrievals. The mean SO2 column contents from aircraft measurements (0.73 DU, Dobson Units) and operational OMI SO2 products (0.63+/-0.26 DU) were close. The OMI retrievals were well correlated with in situ measurements (r = 0.84), but showed low bias (slope = 0.54). A new OMI retrieval algorithm was tested and showed improved agreement and bias (r = 0.87, slope = 0.86). The Community Multiscale Air Quality (CMAQ) model was used to simulate sulfur chemistry, exhibiting reasonable agreement (r = 0.62, slope = 1.33) with in situ SO2 columns. The mean CMAQ SO2 loading over central and eastern China was 54 kT, approx.30% more than the estimate from OMI SO2 products, 42 kT. These numerical simulations, constrained by observations, indicate that ",50% (35 to 61 %) of the anthropogenic sulfur emissions were transported downwind, and the overall lifetime of tropospheric SO2 was 38+/-7 h.

He, Hao↗
NASA NTRS1975

Numerical simulation of turbulence in the presence of shear

The numerical calculations are presented of the large eddy structure of turbulent flows, by use of the averaged Navier-Stokes equations, where averages are taken over spatial regions small compared to the size of the computational grid. The subgrid components of motion are modeled by a local eddy-viscosity model. A new finite-difference scheme is proposed to represent the nonlinear average advective term which has fourth-order accuracy. This scheme exhibits several advantages over existing schemes with regard to the following: (1) the scheme is compact as it extends only one point away in each direction from the point to which it is applied; (2) it gives better resolution for high wave-number waves in the solution of Poisson equation, and (3) it reduces programming complexity and computation time. Examples worked out in detail are the decay of isotropic turbulence, homogeneous turbulent shear flow, and homogeneous turbulent shear flow with system rotation.

Shaanan, S.↗