Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two dimensional flow”

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 1,639 records · Page 91

New Space Domain Processing Technique for Pulsed Laser Velocimetry

A new and fast data processing technique for pulsed-laser-velocimetry data is described. A simplified cross correlation algorithm was developed which uses the time history information of the recorded particle images. It is shown that, the use of the particle time history information (which can be obtained using an electronic camera to record the particle exposures) facilitates high speed processing and unambiguous velocity vector information. It was found that a complex convection vortex flow can be recorded and completely processed in less than 2 min on a personal computer, producing a two-dimensional map of the flow field containing over 300 velocity vectors.

Wernet, Mark P.↗

LIF measurements of scalar mixing in turbulent shear layers

The structure of shear layer flows at high Reynolds numbers remains a very interesting problem. Straight mixing layers have been studied and yielded information on the probability density function (pdf) of a passive scalar across the layer. Konrad and Koochesfahani & Dimotakis measured the pdf of the mixture fraction for mixing layers of moderate Reynolds numbers, each about 25,000 (Re based on velocity difference and visual thickness). Their measurements showed a 'non-marching' pdf (central hump which is invariant from edge to edge across the layer), a result which is linked to the visualizations of the spanwise Kelvin-Helmholtz (K-H) instability mode, which is the primary instability for plane shear layer flows. A secondary instability mode, the Taylor-Gortler (T-G) instability, which is associated with streamwise vortical structures, has also been observed in shear layers. Image reconstruction by Jimenez et al. and volume renderings by Karasso & Mungal at low Re numbers have demonstrated that the K-H and the T-G instability modes occur simultaneously in a non-mutually destructive way, evidence that supports the quasi two-dimensional aspect of these flows and the non-marching character of the pdf at low Reynolds numbers. At higher Re numbers though, the interaction of these two instability modes is still unclear and may affect the mixing process. In this study, we perform measurements of the concentration pdf of plane mixing layers for different operating conditions. At a speed ratio of r = U(sub 1)/U(sub 2) = 4:1, we examine three Reynolds number cases: Re = 14,000, Re = 31,000, and Re = 62,000. Some other Re number cases' results, not presented in detail, are invoked to explain the behavior of the pdf of the concentration field. A case of r = 2.6:1 at Re = 20,000 is also considered. The planar laser-induced fluorescence technique is used to yield quantitative measurements. The different Re are obtained by changing the velocity magnitudes of the two streams. The question of resolution of these measurements is addressed. In order to investigate the effects of the initial conditions on the development and the structure of the mixing layer, the boundary layer on the high-speed side of the splitter plate is tripped. The average concentration and the average mixed fluid concentration are also calculated to further understand the changes in the shear layer for the different cases examined.

Karasso, Paris S.↗

Neutral stability calculations for boundary-layer flows

An analysis is presented of the parallel neutral stability of three-dimensional incompressible, isothermal boundary-layer flows. A Taylor-series expansion of the dispersion relation is used to derive the general eigenvalues. These equations are functions of the complex group velocity. These relations are verified by numerical results obtained for two- and three-dimensional disturbances in two- and three-dimensional flows.

Nayfeh, A. H.↗

Computations of Flow over a Hump Model Using Higher Order Method with Turbulence Modeling

Turbulent separated flow over a two-dimensional hump is computed by solving the RANS equations with k - omega (SST) turbulence model for the baseline, steady suction and oscillatory blowing/suction flow control cases. The flow equations and the turbulent model equations are solved using a fifth-order accurate weighted essentially. nonoscillatory (WENO) scheme for space discretization and a third order, total variation diminishing (TVD) Runge-Kutta scheme for time integration. Qualitatively the computed pressure distributions exhibit the same behavior as those observed in the experiments. The computed separation regions are much longer than those observed experimentally. However, the percentage reduction in the separation region in the steady suction case is closer to what was measured in the experiment. The computations did not predict the expected reduction in the separation length in the oscillatory case. The predicted turbulent quantities are two to three times smaller than the measured values pointing towards the deficiencies in the existing turbulent models when they are applied to strong steady/unsteady separated flows.

Balakumar, P.↗

Investigation of particle–wall interaction in a pseudo-2D fluidized bed using CFD-DEM simulations

Here, we report on discrete element method simulations of a pseudo-two-dimensional (pseudo-2D) fluidized bed to investigate particle–wall interactions. Detailed information on macroscopic flow field variables, including solids pressure, granular temperature, and normal and tangential wall stresses are analyzed. The normal wall stress differs from the solids pressure because of the strong anisotropic flow behavior in the pseudo-2D system. A simple linear relationship exists between normal wall stress and solids pressure. In addition, an effective friction coefficient can be derived to characterize particle–wall flow interaction after evaluating the normal and tangential wall stresses. The effects of inter-particle and particle–wall friction coefficients are evaluated. Strong anisotropic flow behavior in the pseudo-2D system needs to be considered to validate the two-fluid model where the boundary condition is usually developed based on an isotropic assumption. The conclusion has been confirmed by simulation with different particle stiffnesses. Assumptions in the newly developed model for 2D simulation are further examined against the discrete element method simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Application of two-point difference schemes to the conservative Euler equations for one-dimensional flows

An implicit finite difference method is presented for obtaining steady state solutions to the time dependent, conservative Euler equations for flows containing shocks. The method used the two-point differencing approach of Keller with dissipation added at supersonic points via the retarded density concept. Application of the method to the one-dimensional nozzle flow equations for various combinations of subsonic and supersonic boundary conditions shows the method to be very efficient. Residuals are typically reduced to machine zero in approximately 35 time steps for 50 mesh points. It is shown that the scheme offers certain advantages over the more widely used three-point schemes, especially in regard to application of boundary conditions.

Wornom, S. F.↗

Inviscid analysis of jet injection between two moving streams

An analytical method is developed for determining the flow interaction when a two-dimensional jet is injected between two moving streams. The jet is flowing out of channel and is turned as it enters between the external streams. The local velocity variation resulting from the flow interaction provides a static pressure variation along the jet bounding streamlines that is a priori unknown. Hense, the flow must be obtained by coupling the three flow regions (the jet and the free stream on either side) along the jet boundaries. Both external streams have the same total pressure, which is different from that in the jet. The solution is for the condition that the total pressure in the jet does not differ from the free-stream value by a large amount compared with the free-stream dynamic head. Results are given for the shape of the jet boundaries for various injection configurations.

Goldstein, M. E.↗

Iterative solution of transonic flows over airfoils and wings, including flows at Mach 1

A new method of calculating transonic flows based on a 'rotated' difference scheme is described. It is suitable for the calculation of both two- and three-dimensional flows without restriction on the speed at infinity and is well adapted to computer use. The Murman procedure is modified to eliminate any assumptions about the direction of flow when constructing the difference scheme. The proper directional property is obtained by rotating the difference scheme to conform with the local stream direction. In the hyperbolic region retarded difference formulas are used for all contributions to the streamwise second derivative, producing a correctly oriented positive artificial viscosity. In the absence of a simple implicit scheme in the hyperbolic and elliptic regions, the concept of iterations as steps in artificial time is introduced. Computer testing of this procedure provides numerical confirmation of the existence and uniqueness of weak solutions of the potential equation when a suitable entropy inequality is enforced.

Jameson, A.↗

Flow separation in inlets at incidence angles

Wind-tunnel pressure data and flow pictures obtained for two two-dimensional inlet models have been examined to study the internal flow structure and separation at large incidence angles. The inlet models were 12-in. high (diffuser exit height) and had internal contraction ratio of 1.21 and 1.17. They were tested at low forward speeds over a wide range of throat Mach numbers (inlet mass flow rates) and angles of incidence. Characteristic features of the internal flow such as a drastic change of pressure gradient near the highlight, local separation bubbles and shock/boundary-layer interactions have been indicated and discussed. For a few specific cases, the experimental surface pressure distributions have been compared with theoretical predictions.

Jakubowski, A. K.↗

Computations of Flow Over the Hump Model Using Higher-Order Method With Turbulence Modeling

The flow over the two-dimensional hump model is computed by solving the RANS equations with kappa-omega (SST) model. The governing equations, the flow equations and the turbulent equations, are solved using the 5th order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using explicit third order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The WENO and the TVD methods and the formulas are explained in [1] and the application of ENO method to N-S equations is given in [2]. The solution method implemented in this computation is described in detail in [3].

Balakumar, P.↗

Linear stability theory and three-dimensional boundary layer transition

The viewgraphs and discussion of linear stability theory and three dimensional boundary layer transition are provided. The ability to predict, using analytical tools, the location of boundary layer transition over aircraft-type configurations is of great importance to designers interested in laminar flow control (LFC). The e(sup N) method has proven to be fairly effective in predicting, in a consistent manner, the location of the onset of transition for simple geometries in low disturbance environments. This method provides a correlation between the most amplified single normal mode and the experimental location of the onset of transition. Studies indicate that values of N between 8 and 10 correlate well with the onset of transition. For most previous calculations, the mean flows were restricted to two-dimensional or axisymmetric cases, or have employed simple three-dimensional mean flows (e.g., rotating disk, infinite swept wing, or tapered swept wing with straight isobars). Unfortunately, for flows over general wing configurations, and for nearly all flows over fuselage-type bodies at incidence, the analysis of fully three-dimensional flow fields is required. Results obtained for the linear stability of fully three-dimensional boundary layers formed over both wing and fuselage-type geometries, and for both high and low speed flows are discussed. When possible, transition estimates form the e(sup N) method are compared to experimentally determined locations. The stability calculations are made using a modified version of the linear stability code COSAL. Mean flows were computed using both Navier Stokes and boundary-layer codes.

Spall, Robert E.↗

A study in three-dimensional chaotic dynamics: Granular flow and transport in a bi-axial spherical tumbler

We study three-dimensional (3D) chaotic dynamics through an analysis of transport in a granular flow in a half-full spherical tumbler rotated sequentially about two orthogonal axes (a bi-axial “blinking” tumbler). The flow is essentially quasi-two-dimensional in any vertical slice of the sphere during rotation about a single axis, and we provide an explicit exact solution to the model in this case. Hence, the cross-sectional flow can be represented by a twist map, allowing us to express the 3D flow as a linked twist map (LTM). We prove that if the rates of rotation about each axis are equal, then (in the absence of stochasticity) particle trajectories are restricted to two-dimensional (2D) surfaces consisting of a portion of a hemispherical shell closed by a “cap''; if the rotation rates are unequal, then particles can leave the surface they start on and traverse a volume of the tumbler. The period-one structures of the governing LTM are examined in detail: analytical expressions are provided for the location of period-one curves, their extent into the bulk of the granular material, and their dependence on the protocol parameters (rates and durations of rotations). Exploiting the restriction of trajectories to 2D surfaces in the case of equal rotation rates about the axes, a method is proposed for identifying and constructing 3D Kolmogorov--Arnold--Moser (KAM) tubes around the normally elliptic period-one curves. The invariant manifold structure arising from the normally hyperbolic period-one curves is also examined. When the motion is restricted to 2D surfaces, the structure of manifolds of the hyperbolic points in the bulk differs from that corresponding to hyperbolic points in the flowing layer. Each is reminiscent of a template provided by a non-integrable perturbation to a Hamiltonian system, though the governing LTM is not. This highlights the novel 3D chaotic behaviors observed in this model dynamical system.

42 ENGINEERING↗

Computations of Flow Over a Hump Model Using Higher Order Method With Turbulence Modeling

Turbulent separated flow over a two-dimensional hump is computed by solving the RANS equations with k-omega (SST) turbulence model for the baseline, steady suction and oscillatory blowing/suction flow control cases. The flow equations and the turbulent model equations are solved using a fifth-order accurate weighted essentially nonoscillatory (WENO) scheme for space discretization and a third order, total variation diminishing (TVD) Runge-Kutta scheme for time integration. Qualitatively the computed pressure distributions exhibit the same behavior as they are observed in the experiments. The computed separation regions are much longer than that are observed. However, the percentage reduction in the separation region in the steady suction case is closer to that was measured in the experiment. The computations did not predict the expected reduction in the separation length in the oscillatory case. The predicted turbulent quantities are two to three times smaller than that are measured and it points towards the deficiencies in the existing turbulent models when they are applied to strong steady/unsteady separated flows.

Balakumar, Ponnampalam↗

Automated domain decomposition for computational fluid dynamics

Automation of flow-field zoning in two-dimensions is an important step towards easing the three-dimensional grid generation bottleneck in computational fluid dynamics. A knowledge-based approach works well, but several aspects of flow-field zoning make the use of such an approach challenging. A proposed model and language to describe the process of zoning a flow field are presented, followed by a discussion of the implementation of EZGrid, a knowledge-based two-dimensional (2-D) flow-field zoner. Results are shown for representative two-dimensional aerodynamic configurations. Finally, an approach to the evaluation of flow-field zonings is described and used to compare the performance of EZGrid with that of a human expert.

Vogel, Alison Andrews↗

Clogging and jamming transitions in periodic obstacle arrays

We numerically examine clogging transitions for bidisperse disks flowing through a two-dimensional periodic obstacle array. Here, we show that clogging is a probabilistic event that occurs through a transition from a homogeneous flowing state to a heterogeneous or phase-separated jammed state where the disks form dense connected clusters. The probability for clogging to occur during a fixed time increases with increasing particle packing and obstacle number. For driving at different angles with respect to the symmetry direction of the obstacle array, we show that certain directions have a higher clogging susceptibility. It is also possible to have a size-specific clogging transition in which one disk size becomes completely immobile while the other disk size continues to flow.

36 MATERIALS SCIENCE↗

Application of the method of characteristics to supersonic rotational flow

A system for calculating the physical properties of supersonic rotational flow with axial symmetry and supersonic rotational flow in a two-dimensional field was determined by use of the characteristics method. The system was applied to the study of external and internal flow for supersonic inlets with axial symmetry. For a circular conical inlet the shock that occurred at the lip of the inlet became stronger as it approached the axis of the inlet and became a normal shock at the axis. The region in which strong shock occurred increased with increase of the angle of internal cone at the lip of the inlet. For an inlet with a central body the method of characteristics was applied to the design of an internal channel shape that, theoretically, results in very efficient recompression in the inlet. It was shown that if an effuser is connected with the diffuser a body of revolution with very small shock-wave drag can be determined. (author)

Ferri, Antonio↗

The stability to two-dimensional wakes and shear layers at high Mach numbers

This study is concerned with the stability properties of laminar free-shear-layer flows, and in particular symmetric two-dimensional wakes, for the supersonic through the hypersonic regimes. Emphasis is given to the use of proper wake profiles that satisfy the equations of motion at high Reynolds numbers. In particular the inviscid stability of a developing two-dimensional wake is studied as it accelerates at the trailing edge of a splitter plate. The nonparallelism of the flow is a leading-order effect in the calculation of the basic state, which is obtained numerically. Neutral stability characteristics are computed and the hypersonic stability is obtained by increasing the Mach number. It is found that the stability characteristics are altered significantly as the wake develops. Multiple modes (secondary modes) are found in the near wake that are closely related to the corresponding Blasius ones, but as the wake develops mode multiplicity is delayed to higher and higher Mach numbers. At a distance of about one plate length from the trailing edge, there is only one mode in a Mach number range of 0-20. The dominant mode emerging at all wake stations, and for high enough Mach numbers, is the so-called vorticity mode that is centered around the generalized inflection point layer. The structure of the dominant mode is also obtained analytically for all streamwise wake locations and it is shown how the far-wake limit is approached. Asymptotic results for the hypersonic mixing layer given by a tanh and a Lock distribution are also given.

Papageorgiou, Demetrios T.↗

Application of two-point difference schemes to the conservative Euler equations for one-dimensional flows

An implicit finite-difference method is presented for obtaining steady-state solutions to the time-dependent, conservative Euler equations for flows containing shocks. The method uses a two-point central difference scheme with dissipation added at supersonic points via the retarded density concept. Application of the method to the one-dimensional nozzle flow equations for various combinations of subsonic and supersonic boundary conditions show the method to be very efficient. Residuals are typically reduced to machine zero in approximately 35 time steps for 50 mesh points. It is shown that the scheme offers certain advantages over the more widely-used three-point schemes, especially in regard to application of boundary conditions.

Wornom, S. F.↗