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 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Numerical simulation of viscous flows over transonic aircraft configurations

Numerical simulation of compressible viscous flow fields is performed for a transonic transport configuration. A single structured grid system is constructed using analytical transformations such as conformal mapping, shearing/twisting/rotating/clustering/stretching transformations. The Reynolds-averaged, thin-layer Navier-Stokes equations are solved on a supercomputer, FACOM VP-400, using the LU-ADI factorization method.

Takanashi, Susumu↗

NASA NTRS1989

Resolution requirements for numerical simulations of transition

The resolution requirements for direct numerical simulations of transition to turbulence are investigated. A reliable resolution criterion is determined from the results of several detailed simulations of channel and boundary-layer transition.

Zang, Thomas A.↗

NASA NTRS1988

Ion-cyclotron wave heating of heavy ions in the equatorial magnetosphere - A numerical simulation theory

A 1-2/2 dimensional hybrid numerical simulation code is used to study the heating of cold H(+) ions and heavy ions by electromagnetic ion-cyclotron waves (ICWs) in the ring current region of the equatorial magnetosphere. Consideration is given to a plasma consisting of electrons, hot H(+) ions, and cold heavy ions in which the ICWs are driven by the temperature anisotropy of the hot protons. For large-amplitude ICWs, it is found that the cold H(+) ions are preferentially heated over the heavy ions although the cold H(+) ions are heated by a three-step process.

Chen, M. W.↗

NASA NTRS1989

Direct numerical simulations of transition in a compressible wake

Direct numerical simulations of a compressible wake are performed to determine the sensitivity of large-scale motions to initial phase lags between the most unstable two-dimensional wave and a pair of oblique waves, and to determine the effect of Mach number on the development of three-dimensional motions. The vorticity dynamics of the flow are studied, and it is determined that, depending upon the initial phase lag, the interaction of the streamwise vorticity with the spanwise vortices may result in the formation of closed vortex loops, or in the removal of streamwise vorticity from the saddle region between the spanwise vortices. Due to lower growth rates experienced at high Mach numbers, the roll-up process is delayed, and the resultant nonlinear motions are much more diffuse.

Chen, Jacqueline H.↗

NASA NTRS1989

Resolution requirements for numerical simulations of transition

Zang, Thomas A.↗

NASA NTRS1991

Numerical simulation of fluid flow around a scramaccelerator projectile

Numerical simulations of the fluid motion and temperature distribution around a 'scramaccelerator' projectile are obtained for Mach numbers in the 5-10 range. A finite element method is used to solve the equations of motion for inviscid and viscous two-dimensional or axisymmetric compressible flow. The time-dependent equations are solved explicitly, using bilinear isoparametric quadrilateral elements, mass lumping, and a shock-capturing Petrov-Galerkin formulation. Computed results indicate that maintaining on-design performance for controlling and stabilizing oblique detonation waves is critically dependent on projectile shape and Mach number.

Pepper, Darrell W.↗

NASA NTRS2014

Direct Numerical Simulation of a Weakly Stratified Turbulent Wake

Direct numerical simulation (DNS) is used to investigate a time-dependent turbulent wake evolving in a stably stratified background. A large initial Froude number is chosen to allow the wake to become fully turbulent and axisymmetric before stratification affects the spreading rate of the mean defect. The uncertainty introduced by the finite sample size associated with gathering statistics from a simulation of a time-dependent flow is reduced, compared to earlier simulations of this flow. The DNS reveals the buoyancy-induced changes to the turbulence structure, as well as to the mean-defect history and the terms in the mean-momentum and turbulence-kinetic-energy budgets, that characterize the various states of this flow - namely the three-dimensional (essentially unstratified), non-equilibrium (or 'wake-collapse') and quasi-two-dimensional (or 'two-component') regimes observed elsewhere for wakes embedded in both weakly and strongly stratified backgrounds. The wake-collapse regime is not accompanied by transfer (or 'reconversion') of the potential energy of the turbulence to the kinetic energy of the turbulence, implying that this is not an essential feature of stratified-wake dynamics. The dependence upon Reynolds number of the duration of the wake-collapse period is demonstrated, and the effect of the details of the initial/near-field conditions of the wake on its subsequent development is examined.

Redford, J. A.↗

NASA NTRS2000

Accuracy of Numerical Simulations of Tip Clearance Flow in Transonic Compressor Rotors Improved Dramatically

The tip clearance flows of transonic compressor rotors have a significant impact on rotor and stage performance. Although numerical simulations of these flows are quite sophisticated, they are seldom verified through rigorous comparisons of numerical and measured data because, in high-speed machines, measurements acquired in sufficient detail to be useful are rare. Researchers at the NASA Glenn Research Center at Lewis Field compared measured tip clearance flow details (e.g., trajectory and radial extent) of the NASA Rotor 35 with results obtained from a numerical simulation. Previous investigations had focused on capturing the detailed development of the jetlike flow leaking through the clearance gap between the rotating blade tip and the stationary compressor shroud. However, we discovered that the simulation accuracy depends primarily on capturing the detailed development of a wall-bounded shear layer formed by the relative motion between the leakage jet and the shroud.

VanZante, Dale E.↗

NASA NTRS1993

Multigrid direct numerical simulation of the whole process of flow transition in 3-D boundary layers

A new technology was developed in this study which provides a successful numerical simulation of the whole process of flow transition in 3-D boundary layers, including linear growth, secondary instability, breakdown, and transition at relatively low CPU cost. Most other spatial numerical simulations require high CPU cost and blow up at the stage of flow breakdown. A fourth-order finite difference scheme on stretched and staggered grids, a fully implicit time marching technique, a semi-coarsening multigrid based on the so-called approximate line-box relaxation, and a buffer domain for the outflow boundary conditions were all used for high-order accuracy, good stability, and fast convergence. A new fine-coarse-fine grid mapping technique was developed to keep the code running after the laminar flow breaks down. The computational results are in good agreement with linear stability theory, secondary instability theory, and some experiments. The cost for a typical case with 162 x 34 x 34 grid is around 2 CRAY-YMP CPU hours for 10 T-S periods.

Liu, Chaoqun↗

DOE OSTI2016

Systematic parameter study of a nonlinear electromagnetic energy harvester with matched magnetic orientation: Numerical simulation and experimental investigation

Here, this paper reports the systematic parameter study of a tristable nonlinear electromagnetic energy harvester for ambient low-frequency vibration. Numerical simulations and experimental investigations are performed on the harvester which consists of a cantilever beam, a tip coil, two tip magnets and two external side magnets. The external side magnets are deployed symmetrically along a concave surface parallel to the trajectory of the cantilever tip with a controllable distance so that the magnetic orientation of the tip magnets are matched with that of the side magnets. Therefore, instead of the ternary position parameters (d, h, α), a binary parameters pair (d 0 , d) is used to characterize the position of the side magnets and the performance of the energy harvester. The magnetic force and magnetic field on the cantilever tip therefore depend on the relative distance in the tip displacement direction between the tip magnets and side magnets, but is independent of the position of the side magnets on the concave surface. The magnetic force (field)-distance relationship is measured experimentally and curve fitted to obtain explicit expressions, in order to characterize the magnetic force (field) when the side magnets are placed at varied positions along the concave surface. Numerical simulation is, then, performed to predict the electromagnetic voltage output and the bandwidth of the energy harvester. The simulation results coincided with the measured data. Significant broadband response is obtained experimentally and the maximum RMS power output is 40.2 mW at 0.45 g of excitation. Furthermore, the proposed structure showcasing the matched magnetic orientation is characterized by the binary parameters pair (d 0 , d) and the systematic parametric approach could contribute to the design and study of nonlinear broadband energy harvesters.

30 DIRECT ENERGY CONVERSION↗

NASA NTRS1986

Numerical simulation of unsteady flow in an axisymmetric shear layer

Numerical simulation of unsteady flow in an axisymmetric subsonic shear layer is accomplished by solving the time-dependent compressible Navier-Stokes equations. The objective of the effort is to investigate by numerical means the influence of various flow parameters on the shear layer behavior. The parameters investigated include the velocity ratio of two streams, total temperature, and nozzle lip thickness. The computations are performed on a CRAY-IS computer using McCormack's explicit finite difference scheme. The computed results generally show qualitative agreement with experimental data.

Scott, J. N.↗

NASA NTRS1998

Spurious Numerical Oscillations in Numerical Simulation of Supersonic Flows Using Shock Capturing Schemes

The numerical simulation of transitional and turbulent processes in hypersonic boundary layers often involves a physical process of a shock-disturbance wave interaction in complex two-dimensional and three-dimensional flow fields. For such simulations, it is required that there be a high order of accuracy in capturing both the shock and the small disturbances. The purpose of this paper is to evaluate the viability of using high order shock capturing schemes to track small disturbances in a multi-dimensional steady hypersonic flow. The numerical methods that are to be studied are the Total Variation Diminishing (TVD) scheme, and Essentially Non-Oscillatory (ENO) scheme. This paper shows that the presence of numerical oscillations in the flow field solution may drastically hinder any attempt at tracking the propagation of any physical disturbances. It has been found that the numerical oscillations that exist for shock capturing methods may be significant enough to pollute a flow field containing small physical disturbances. The effects of the refinement of the grid do not reduce the oscillations, but rather they decrease the wavelength of the oscillations. It is shown that by aligning the shock with the grid, the amplitude of these spurious oscillations may be greatly reduced.

Lee, Theodore K.↗

NASA NTRS1994

Numerical simulation of supersonic boundary layer transition

The present contribution reviews some of the recent progress obtained at our group in the direct numerical simulation (DNS) of compressible boundary layer transition. Elements of the different simulation approaches and numerical techniques employed are surveyed. Temporal and spatial simulations, as well as comparisons with results obtained from Parabolized Stability Equations, are discussed. DNS results are given for flat plate boundary layers in the Mach number range 1.6 to 4.5. A temporal DNS at Mach 4.5 has been continued through breakdown all the way to the turbulent stage. In addition results obtained with a recently developed extended temporal DNS approach are presented, which takes into account some nonparallel effects of a growing boundary layer. Results from this approach are quite close to those of spatial DNS, while preserving the efficiency of the temporal DNS.

Guo, Y.↗

NASA NTRS1994

Direct numerical simulations of turbulent non-premixed methane-air flames modeled with reduced kinetics

Turbulent non-premixed stoichiometric methane-air flames modeled with reduced kinetics have been studied using the direct numerical simulation approach. The simulations include realistic chemical kinetics, and the molecular transport is modeled with constant Lewis numbers for individual species. The effect of turbulence on the internal flame structure and extinction characteristics of methane-air flames is evaluated. Consistent with earlier DNS with simple one-step chemistry, the flame is wrinkled and in some regions extinguished by the turbulence, while the turbulence is weakened in the vicinity of the flame due to a combination of dilatation and an increase in kinematic viscosity. Unlike previous results, reignition is observed in the present simulations. Lewis number effects are important in determining the local stoichiometry of the flame. The results presented in this work are preliminary but demonstrate the feasibility of incorporating reduced kinetics for the oxidation of methane with direct numerical simulations of homogeneous turbulence to evaluate the limitations of various levels of reduction in the kinetics and to address the formation of thermal and prompt NO(x).

Card, J. M.↗

DOE OSTI2022

1D model of tokamak scrape-off layer and divertor plasmas: Dimensionless parameters and the results of numerical simulations

Solutions of a 1D model of the scrape-off layer and a divertor plasma in a tokamak are analyzed on the basis of dimensionless parameters governing the physics of both plasma and neutral transport and atomic physics processes. Numerical simulations based on fluid plasma and Monte Carlo neutral modeling confirm earlier results obtained with simplified analytic considerations and numerical simulations performed with fluid neutral description and demonstrate bifurcation of plasma parameters. Further, these simulations also confirm that the major physical reasons for the rollover of the plasma flux to the target, which is the key feature of the phenomenon of divertor plasma detachment, are the impurity radiation loss and plasma recombination. It is suggested that simulations of the edge plasma based on dimensionless parameters provide a powerful tool for checking the consistency of the physics incorporated in existing complex 2D codes used for edge plasma modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2002

Entropy Splitting for High Order Numerical Simulation of Compressible Turbulence

A stable high-order numerical scheme for direct numerical simulation (DNS) of shock-free compressible turbulence is presented. The method is applicable to general geometries. It contains no upwinding, artificial dissipation, or filtering. Instead the method relies on the stabilizing mechanisms of an appropriate conditioning of the governing equations and the use of compatible spatial difference operators for the interior points (interior scheme) as well as the boundary points (boundary scheme). An entropy-splitting approach splits the inviscid flux derivatives into conservative and nonconservative portions. The spatial difference operators satisfy a summation-by-parts condition, leading to a stable scheme (combined interior and boundary schemes) for the initial boundary value problem using a generalized energy estimate. A Laplacian formulation of the viscous and heat conduction terms on the right hand side of the Navier–Stokes equations is used to ensure that any tendency to odd–even decoupling associated with central schemes can be countered by the fluid viscosity. The resulting methods are able to minimize the spurious high-frequency oscillations associated with pure central schemes, especially for long time integration applications such as DNS. For validation purposes, the methods are tested in a DNS of compressible turbulent plane channel flow at low values of friction Mach number, where reference turbulence data bases exist. It is demonstrated that the methods are robust in terms of grid resolution, and in good agreement with published channel data. Accurate turbulence statistics can be obtained with moderate grid sizes. Stability limits on the range of the splitting parameter are determined from numerical tests.

N.D. Sandham↗

DOE OSTI2020

Numerical simulations of onset and growth of Rayleigh–Taylor instability involving solids in converging geometry

Numerical simulation of Rayleigh–Taylor instabilities involving solids has been carried out in a spherically converging geometry. A solid shell is driven toward the center by high-pressure gas outside. Low-pressure gas inside the shell is then compressed, causing the solid shell to slow down and eventually bounce back. Here, simulation results show the absence of the onset condition unlike in the planar geometry equivalent. Results also show that disturbance growth seems to be strongly associated with the deformation direction instead of the acceleration direction. For example, when the interface is compressed in the cross-stream direction, the disturbance at the gas–solid interface grows even under the Rayleigh–Taylor stable situation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1992

Recent advances in numerical simulation and control of asymmetric flows around slender bodies

The problems of asymmetric flow around slender bodies and its control are formulated using the unsteady, compressible, thin-layer or full Navier-Stokes equations which are solved using an implicit, flux-difference splitting, finite-volume scheme. The problem is numerically simulated for both locally-conical and three-dimensional flows. The numerical applications include studies of the effects of relative incidence, Mach number and Reynolds number on the flow asymmetry. For the control of flow asymmetry, the numerical simulation cover passive and active control methods. For the passive control, the effectiveness of vertical fins placed in the leeward plane of geometric symmetry and side strakes with different orientations is studied. For the active control, the effectiveness of normal and tangential flow injection and surface heating and a combination of these methods is studied.

Kandil, Osama A.↗