Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Simulations”

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 181 records · Page 10

A numerical simulation of finite-length Taylor-Couette flow

Results from numerical simulations of finite-length Taylor-Couette flow are presented. Included are time-accurate and steady-state studies of the change in the nature of the symmetric two-cell/asymmetric one-cell bifurcation with varying aspect ratio and of the Reynolds number/aspect ratio locus of the two-cell/four-cell bifurcation. Preliminary results from wavy-vortex simulations at low aspect ratios are also presented.

Streett, C. L.↗

Direct numerical simulation of isotropic turbulence interacting with a weak shock wave

Direct numerical simulations are used to investigate the interaction of isotropic quasi-incompressible turbulence with a weak shock wave. A linear analysis of the interaction is conducted for comparison with the simulations. Both the simulations and the analysis show that turbulence is enhanced during the interaction. Turbulent kinetic energy and transverse vorticity components are amplified, and turbulent lengthscales are decreased. It is suggested that the amplification mechanism is primarily linear. Simulations also showed a rapid evolution of turbulent kinetic energy just downstream of the shock, a behavior not reproduced by the linear analysis. Analysis of the budget of the turbulent kinetic energy transport equation shows that this behavior can be attributed to the pressure transport term. Multiple compression peaks were found along the mean streamlines at locations where the local shock thickness had increased significantly.

Lee, Sangsan↗

Direct numerical simulation of turbulent plane Couette flow

Turbulent plane Couette flow was numerically simulated at a Reynolds number (U(sub w)h/nu) of 6000, where U(sub w) is the relative wall speed and h is half the channel-height. Unlike in Poiseuille flow, where the mean shear rate changes its sign at the centerline, the sign of mean shear rate in plane Couette flow remains the same across the whole channel. This difference is expected to yield several differences between the two flows, especially in the core region. The most significant and dramatic difference observed was the existence of large-scale structures in the core region of the plane Couette flow. The large eddies are extremely long in the flow direction and fill the entire channel (i.e., their vertical extent is 2h). The large-scale structures have the largest contribution from the wavenumber (k(sub x)h,k(sub z)h) = (0, plus or minus 1.5), corresponding to a wavelength lambda(sub z)/h is approximately equal to 4. The secondary motion associated with the k(sub x)h = 0 mode consists of the large-scale vortices. The large eddies contribute about 30 percent of turbulent kinetic energy.

Lee, Moon Joo↗

Numerical simulation of thermally induced near-surface flows over Martian terrain

Numerical simulations of the Martian near-surface wind regime using a mesoscale atmospheric model have shown that the thermally induced near-surface winds are analogous to terrestrial circulations. In particular, katabatic wind displays a striking similarity to flow observed over Antarctica. Introduction of solar radiation strongly perturbs the slope flows; anabatic conditions develop in middle to high latitudes during the daytime hours due to the solar heating of the sloping terrain. There appears to be a rapid transition from the katabatic to the anabatic flow regimes, emphasizing the primary importance of radiative exchanges at the surface in specifying the horizontal pressure gradient force.

Parish, T. R.↗

Numerical Simulations of Mass Loading in the Solar Wind Interaction with Venus

Numerical simulations are performed in the framework of nonlinear two-dimensional magnetohydrodynamics to investigate the influence of mass loading on the solar wind interaction with Venus. The principal physical features of the interaction of the solar wind with the atmosphere of Venus are presented. The formation of the bow shock, the magnetic barrier, and the magnetotail are some typical features of the interaction. The deceleration of the solar wind due to the mass loading near Venus is an additional feature. The effect of the mass loading is to push the shock farther outward from the planet. The influence of different values of the magnetic field strength on plasma evolution is considered.

Murawski, K.↗

Direct Numerical Simulation of a Coolant Jet in a Periodic Crossflow

A Direct Numerical Simulation of a coolant jet injected normally into a periodic crossflow is presented. The physical situation simulated represents a periodic module in a coolant hole array with a heated crossflow. A collocated finite difference scheme is used which is fifth-order accurate spatially and second-order accurate temporally. The scheme is based on a fractional step approach and requires the solution of a pressure-Poisson equation. The simulations are obtained for a blowing ratio of 0.25 and a channel Reynolds number of 5600. The simulations reveal the dynamics of several large scale structures including the Counter-rotating Vortex Pair (CVP), the horse-shoe vortex, the shear layer vortex, the wall vortex and the wake vortex. The origins and the interactions of these vortical structures are identified and explored. Also presented are the turbulence statistics and how they relate to the flow structures.

Sharma, Chirdeep↗

Numerical Simulation of Non-rotating and Rotating Coolant Channel Flow Fields

Numerical simulations of the three-dimensional flow and heat transfer in a rectangular duct with a bend were performed. Results are presented for Reynolds numbers of 17,000 and 37,000 and for aspect ratios of 0.5 and 1.0. A k-(omega) turbulence model with no reference to distance to a wall is used. Direct comparison between single block and multiblock grid calculations are made. Heat transfer and velocity distributions are compared to available literature with good agreement. The multi-block grid system is seen to produce more accurate results compared to a single-block grid with the same number of cells.

Rigby, David L.↗

Aspects of Supercritical Turbulence: Direct Numerical Simulation of O2/H2 and C7H16/N2 Temporal Mixing Layers

Results from Direct Numerical Simulations of temporal, supercritical mixing layers for two species systems are analyzed to elucidate species-specific turbulence aspects. The two species systems, O2/H2 and C7HG16/N2, have different thermodynamic characteristics; thus, although the simulations are performed at similar reduced pressure (ratio of the pressure to the critical pressure), the former system is dose to mixture ideality and has a relatively high solubility with respect to the latter, which exhibits strong departures from mixture ideality Due to the specified, smaller initial density stratification, the C7H16/N2 layers display higher growth and increased global molecular mixing as well as larger turbulence levels. However, smaller density gradients at the transitional state for the O2/H2 system indicate that on a local basis, the layer exhibits an enhanced mixing, this being attributed to the increased solubility and to mixture ideality. These thermodynamic features are shown to affect the irreversible entropy production (i.e. the dissipation), which is larger for the O2/H2 layer and is primarily concentrated in high density-gradient magnitude regions that are distortions of the initial density stratification boundary. In contrast, the regions of largest dissipation in the C7H16/N2 layer are located in high density-gradient magnitude regions resulting from the mixing of the two fluids.

Bellan, J.↗

Numerical Simulations of Cosmogenic Neutron Production and Transport in Planetary Surfaces

The numerical simulation code MCNPX was used to calculate the production and transport of cosmic-ray-produced neutrons in Mars and meteorites. These calculations help to understand the processes involved and the parameters that control the neutron fluxes. Results are presented here for neutrons in Mars and for the distribution of cosmic-ray neutrons and protons in a 50 cm-radius L-chondrite.

Kim, Kyeong J.↗

Numerical Simulations Studies of the Convective Instability Onset in a Supercritical Fluid

Numerical simulation studies are reported for the convection of a supercritical fluid, He-3, in a Rayleigh-Benard cell. The calculations provide the temporal profile DeltaT(t) of the temperature drop across the fluid layer. In a previous article, systematic delays in the onset of the convective instability in simulations relative to experiments were reported, as seen from the DeltaT(t) profiles. They were attributed to the smallness of the noise which is needed to start the instability. Therefore i) homogeneous temperature noise and ii) spatial lateral periodic temperature variations in the top plate were programmed into the simulations, and DeltaT(t) compared with that of an experiment with the same fluid parameters. An effective speed-up in the instability onset was obtained, with the best results obtained through the spatial temperature variations with a period of 2L, close to the wavelength of a pair of convections rolls. For a small amplitude of 0.5 micro-K, this perturbation gave a semiquantitative agreement with experimental observations. Results for various noise amplitudes are presented and discussed in relation to predictions by El Khouri and Carl es.

Furukawa, A.↗

Direct Numerical Simulation of Transition in a Swept-Wing Boundary Layer

Direct numerical simulation (DNS) is performed to examine laminar to turbulent transition due to high-frequency secondary instability of stationary crossflow vortices in a subsonic swept-wing boundary layer for a realistic natural-laminar-flow airfoil configuration. The secondary instability is introduced via inflow forcing derived from a two-dimensional, partial-differential-equation based eigenvalue computation; and the mode selected for forcing corresponds to the most amplified secondary instability mode which, in this case, derives a majority of its growth from energy production mechanisms associated with the wall-normal shear of the stationary basic state. Both the growth of the secondary instability wave and the resulting onset of laminar-turbulent transition are captured within the DNS computations. The growth of the secondary instability wave in the DNS solution compares well with linear secondary instability theory when the amplitude is small; the linear growth is followed by a region of reduced growth resulting from nonlinear effects before an explosive onset of laminar breakdown to turbulence. The peak fluctuations are concentrated near the boundary layer edge during the initial stage of transition, but rapidly propagates towards the surface during the process of laminar breakdown. Both time-averaged statistics and flow visualization based on the DNS reveal a sawtooth transition pattern that is analogous to previously documented surface flow visualizations of transition due to stationary crossflow instability. The memory of the stationary crossflow vortex is found to persist through the transition zone and well beyond the location of the maximum skin friction.

Duan, Lian↗

Direct Numerical Simulations of High-Speed Turbulent Boundary Layers over Riblets

Direct numerical simulations (DNS) of spatially developing turbulent boundary layers over riblets with a broad range of riblet spacings are conducted to investigate the effects of riblets on skin friction at high speeds. Zero-pressure gradient boundary layers under two flow conditions (Mach 2:5 with T(sub w)/T(sub r) = 1 and Mach 7:2 with T(sub w)/T(sub r) = 0:5) are considered. The DNS results show that the drag-reduction curve (delta C(sub f)/C(sub f) vs l(sup +)(sub g )) at both supersonic speeds follows the trend of low-speed data and consists of a `viscous' regime for small riblet size, a `breakdown' regime with optimal drag reduction, and a `drag-increasing' regime for larger riblet sizes. At l l(sup +)(sub g) approx. 10 (corresponding to s+ approx 20 for the current triangular riblets), drag reduction of approximately 7% is achieved at both Mach numbers, and con rms the observations of the few existing experiments under supersonic conditions. The Mach- number dependence of the drag-reduction curve occurs for riblet sizes that are larger than the optimal size, with smaller slopes of (delta C(sub f)/C(sub f) for larger freestream Mach numbers. The Reynolds analogy holds with 2(C(sub h)=C(sub f) approximately equal to that of at plates for both drag-reducing and drag-increasing configurations.

Duan, Lian↗

Numerical Simulation of a Complete Low-Speed Wind Tunnel Circuit

A numerical simulation of the complete circuit of the NASA Langley 14 x 22-ft low-speed wind tunnel is described. Inside the circuit, all turning vanes are modeled as well as the five flow control vanes downstream of the 1st corner. The fan drive system is modeled using an actuator disk for the fan blades coupled with the fan nacelle. All the surfaces are modeled as viscous walls except the turning vanes, which were modeled as inviscid surfaces. NASA Langley's TetrUSS unstructured grid software was used for grid generation and flow simulation. Two turbulence models were employed in the present study, namely, the one-equation Spalart-Allmaras model and the shear stress transport (SST) model of Menter. The paper shows the flow characteristics in the circuit and compares the results with experimental data where available.

Nayani, Sudheer N.↗

Numerical simulation of natural convection in a spherical container due to cooling at the center (idealization of the Lal/Kroes experiment)

Natural convection in a spherical container with cooling at the center was numerically simulated using a numerical fluid dynamics computer program. The numerical analysis was simplified by assuming axisymmetric flow in the spherical container, with the symmetry axis being a sphere diagonal parallel to the gravity vector. This axisymmetric spherical geometry was intended as an idealization of the proposed Lal/Kroes crystal growing experiment to be performed on Spacelab. Results were obtained for a range of Rayleigh numbers from 25 to 10,000. The computed velocities were found to be approximately proportional to the Rayleigh number over the range of Rayleigh numbers investigated.

Robertson, S. J.↗

Direct Numerical Simulation of a Temporally Evolving Incompressible Plane Wake: Effect of Initial Conditions on Evolution and Topology

Direct numerical simulations have been used to examine the effect of the initial disturbance field on the development of three-dimensionality and the transition to turbulence in the incompressible plane wake. The simulations were performed using a new numerical method for solving the time-dependent, three-dimensional, incompressible Navier-Stokes equations in flows with one infinite and two periodic directions. The method uses standard Fast Fourier Transforms and is applicable to cases where the vorticity field is compact in the infinite direction. Initial disturbances fields examined were combinations of two-dimensional waves and symmetric pairs of 60 deg oblique waves at the fundamental, subharmonic, and sub-subharmonic wavelengths. The results of these simulations indicate that the presence of 60 deg disturbances at the subharmonic streamwise wavelength results in the development of strong coherent three-dimensional structures. The resulting strong three-dimensional rate-of-strain triggers the growth of intense fine scale motions. Wakes initiated with 60 deg disturbances at the fundamental streamwise wavelength develop weak coherent streamwise structures, and do not develop significant fine scale motions, even at high Reynolds numbers. The wakes which develop strong three-dimensional structures exhibit growth rates on par with experimentally observed turbulent plane wakes. Wakes which develop only weak three-dimensional structures exhibit significantly lower late time growth rates. Preliminary studies of wakes initiated with an oblique fundamental and a two-dimensional subharmonic, which develop asymmetric coherent oblique structures at the subharmonic wavelength, indicate that significant fine scale motions only develop if the resulting oblique structures are above an angle of approximately 45 deg.

Sondergaard, R.↗

Numerical Simulation of a High Mach Number Jet Flow

The recent efforts to develop accurate numerical schemes for transition and turbulent flows are motivated, among other factors, by the need for accurate prediction of flow noise. The success of developing high speed civil transport plane (HSCT) is contingent upon our understanding and suppression of the jet exhaust noise. The radiated sound can be directly obtained by solving the full (time-dependent) compressible Navier-Stokes equations. However, this requires computational storage that is beyond currently available machines. This difficulty can be overcome by limiting the solution domain to the near field where the jet is nonlinear and then use acoustic analogy (e.g., Lighthill) to relate the far-field noise to the near-field sources. The later requires obtaining the time-dependent flow field. The other difficulty in aeroacoustics computations is that at high Reynolds numbers the turbulent flow has a large range of scales. Direct numerical simulations (DNS) cannot obtain all the scales of motion at high Reynolds number of technological interest. However, it is believed that the large scale structure is more efficient than the small-scale structure in radiating noise. Thus, one can model the small scales and calculate the acoustically active scales. The large scale structure in the noise-producing initial region of the jet can be viewed as a wavelike nature, the net radiated sound is the net cancellation after integration over space. As such, aeroacoustics computations are highly sensitive to errors in computing the sound sources. It is therefore essential to use a high-order numerical scheme to predict the flow field. The present paper presents the first step in a ongoing effort to predict jet noise. The emphasis here is in accurate prediction of the unsteady flow field. We solve the full time-dependent Navier-Stokes equations by a high order finite difference method. Time accurate spatial simulations of both plane and axisymmetric jet are presented. Jet Mach numbers of 1.5 and 2.1 are considered. Reynolds number in the simulations was about a million. Our numerical model is based on the 2-4 scheme by Gottlieb & Turkel. Bayliss et al. applied the 2-4 scheme in boundary layer computations. This scheme was also used by Ragab and Sheen to study the nonlinear development of supersonic instability waves in a mixing layer. In this study, we present two dimensional direct simulation results for both plane and axisymmetric jets. These results are compared with linear theory predictions. These computations were made for near nozzle exit region and velocity in spanwise/azimuthal direction was assumed to be zero.

Hayder, M. Ehtesham↗

Investigations of Pole Erosion Mechanisms in the 12.5 kW HERMeS Hall Thruster Using Numerical Simulations and Ion Velocity Measurements

When the Hall Effect Rocket with Magnetic Shielding (HERMeS) was tested at the NASA GRC’s VF-5 the erosion rate on the pole covers was found to be at least a factor of two higher at 300 V and 20.8 A than at other operating conditions in the 400 V to 600 V range. Simulations using our hydrodynamics code Hall2De do not predict increased erosion rates at 300 V but accurately compute erosion rates similar to the measurements at the other operating conditions. We investigate the source of the discrepancy between measurements and simulations at 300 V using a combination of numerical simulations and experimental measurements of the ion velocity fields in the acceleration region of the thruster. By examining previous simulations that predicted the measured erosion rates for other thrusters and operating conditions, we determine that sputtering of the pole surfaces by high energy ions is the most likely mechanism behind the erosion rates at 300 V. High energy ions sputter the pole surfaces when the acceleration region of the thruster is downstream of the pole surface plane. When the latter occurs, the curvature of the plasma potential contours at the edges of the channel accelerates a small fraction of the high energy ions radially. We find that neither our Hall2De simulations nor the experimental measurements produce ions of sufficiently high energy and flux to the inner pole cover to explain the measured erosion there. Only at the outer edge of the cover is where we find the ions needed to yield simulation results that are comparable to the erosion measurements. We argue that one possible mechanism for the higher erosion is high energy ions that graze the channel corner of the pole and become trapped in a sheath that develops downstream. Within the sheath, we find that the electric field is large enough to turn the ions towards the surface of the pole. A simplified calculation shows that the erosion rates produced by sputtering of trapped ions are similar to the measurements. We also propose that the presence of local oscillations and high-energy cathode ions may be sources of the enhanced erosion We conclude this article by proposing a combination of analysis, simulation, and experimental measurements that can be used to address the validity of these hypotheses.

Mikellides, Ioannis G.↗