Exact free-surface flows for shallow-water equations .1. - the incompressible case

A comprehensive exact treatment of free surface flows governed by shallow water equations (in sigma variables) is given. Several new families of exact solutions of the governing PDEs are found and are shown to embed the well-known self-similar or traveling wave solutions which themselves are governed by reduced ODEs. The classes of solutions found here are explicit in contrast to those found earlier in an implicit form. The height of the free surface for each family of solutions is found explicitly. For the traveling or simple wave, the free surface is governed by a nonlinear wave equation, but is arbitrary otherwise. For other types of solutions, the height of the free surface is constant either on lines of constant acceleration or on lines of constant speed; in another case, the free surface is a horizontal plane while the flow underneath is a sine wave. The existence of simple waves on shear flows is analytically proved. The interaction of large amplitude progressive waves with shear flow is also studied.

Flat rotation curves: A result of the helical inverse cascade in turbulent media

We have modeled the rotation curves of 21 galaxies observed by Amram et al. (1992), by combining the effects of rigid rotation, gravity, and turbulence. The main motivation behind such modeling is to study the formation of coherent structures in turbulent media and explore its role in the large-scale structures of the universe. The values of the parameters such as mass, turbulent velocity, and angular velocity derived from the rotation curve fits are in good agreement with those derived from the prevalent models.

Experimental studies of a distorted turbulent spot in a three-dimensional flow

We report here on the results of a series of experiments carried out on a turbulent spot in a distorted duct to study the effects of a divergence with straight streamlines preceded by a short stretch of transverse streamline curvature, both in the absence of any pressure gradient. It is found that the distortion produces substantial asymmetry in the spot: the angles at which the spot cuts across the local streamlines are altered dramatically (in contradiction of a hypothesis commonly made in transition zone modelling), and the Tollmien-Schlichting waves that accompany the wing tips of the spot are much stronger on the outside of the bend than on the inside. However there is no strong effect on the internal structure of the spot and the eddies therein, or on such propagation characteristics as overall spread rate and the celerities of the leading and trailing edges. Both lateral streamline curvature and non-homogeneity of the laminar boundary layer into which the spot propagates are shown to be strong factors responsible for the observed asymmetry. It is concluded that these factors produce chiefly a geometric distortion of the coherent structure in the spot, but do not otherwise affect its dynamics in any significant way.

Exact free-surface flows for shallow-water equations .2. the compressible case

Exact free surface flows with shear in a compressible barotropic medium are found, extending the authors' earlier work for the incompressible medium. The barotropic medium is of finite extent in the vertical direction, while it is infinite in the horizontal direction. The ''shallow water'' equations for a compressible barotropic medium, subject to boundary conditions at the free surface and at the bottom, are solved in terms of double psi-series, Simple wave and time-dependent solutions are found; for the former the free surface is of arbitrary shape while for the latter it is a damping traveling wave in the horizontal direction, For other types of solutions, the height of the free surface is constant either on lines of constant acceleration or on lines of constant speed. In the case of an isothermal medium, when gamma = 1, we again find simple wave and time-dependent solutions.

Turbulent heat flux variation over the Monsoon-Trough region during MONTBLEX-90

Tower data collected during the Monsoon-Trough Boundary Layer Experiment (MONTBLEX-90) have been analysed to understand the observed structure of the surface layer over an arid region (Jodhpur) and a moist region (Kharagpur) during active and weak phases of the 1990 southwest monsoon. Turbulent heat and momentum fluxes are estimated by the eddy correlation method using sonic data. The turbulent momentum flux at both Jodhpur and Kharagpur was larger when the winds were stronger, reaching a maximum of the order of 0.5 N m(-2) on 5 and 6 August when a low pressure system was located over the region. The heat flux at Jodhpur is high during weak monsoon days, the maximum being 450 W m(-2), whereas during active days the flux never exceeds 200 W m(-2). At Kharagpur, the flux does not vary significantly between active and weak monsoon days, the maximum in either phase being 160 W m(-2) At Jodhpur, there is significant contrast in the near-surface air temperature, being higher during weak monsoon days as compared to active days. Cloud cover did not vary significantly in both the regions. The turbulent heat flux variation at both the sites appears to be correlated mainly with soil mixture, and less sensitive to cloud cover.

A finite element method for compressible viscous fluid flows

This work presents a mixed three-dimensional finite element formulation for analyzing compressible viscous flows. The formulation is based on the primitive variables velocity, density, temperature and pressure. The goal of this work is to present a `stable' numerical formulation, and, thus, the interpolation functions for the field variables are chosen so as to satisfy the inf-sup conditions. An exact tangent stiffness matrix is derived for the formulation, which ensures a quadratic rate of convergence. The good performance of the proposed strategy is shown in a number of steady-state and transient problems where compressibility effects are important such as high Mach number flows, natural convection, Riemann problems, etc., and also on problems where the fluid can be treated as almost incompressible. Copyright (C) 2010 John Wiley & Sons, Ltd.

Turbulent near-wake development behind a flat plate

An asymptotic analysis of the two-dimensional turbulent near-wake flow behind a Rat plate with sharp trailing edge has been formulated, The feature that the near-wake, which is dominated by the mixing of the oncoming turbulent boundary layers retains, to a large extent, the memory of the turbulent structure of the upstream boundary layer has been exploited to develop the analysis. This analysis leads to two regions of the near-wake flow (the inner near-wake and the outer near-wake) for which the governing equations are derived. The matching conditions among these regions lead to a logarithmic variation in the normal direction in the overlapping region surrounding the inner near-wake. These features are validated by the available experimental data. Similarity solutions for the velocity distribution (which satisfy the required matching conditions) in the inner near-wake and outer near-wake regions have been obtained by making the appropriate eddy-viscosity assumptions, Uniformly valid solutions for velocity distribution have been constructed for the near-wake. The solutions show good agreement with available experimental data. (C) Elsevier, Paris.

Assessment of conditional moment closure models of turbulent autoignition using DNS data

A systematic assessment of the submodels of conditional moment closure (CMC) formalism for the autoignition problem is carried out using direct numerical simulation (DNS) data. An initially non-premixed, n-heptane/air system, subjected to a three-dimensional, homogeneous, isotropic, and decaying turbulence, is considered. Two kinetic schemes, (1) a one-step and (2) a reduced four-step reaction mechanism, are considered for chemistry An alternative formulation is developed for closure of the mean chemical source term , based on the condition that the instantaneous fluctuation of excess temperature is small. With this model, it is shown that the CMC equations describe the autoignition process all the way up to near the equilibrium limit. The effect of second-order terms (namely, conditional variance of temperature excess sigma(2) and conditional correlations of species q(ij)) in modeling is examined. Comparison with DNS data shows that sigma(2) has little effect on the predicted conditional mean temperature evolution, if the average conditional scalar dissipation rate is properly modeled. Using DNS data, a correction factor is introduced in the modeling of nonlinear terms to include the effect of species fluctuations. Computations including such a correction factor show that the species conditional correlations q(ij) have little effect on model predictions with a one-step reaction, but those q(ij) involving intermediate species are found to be crucial when four-step reduced kinetics is considered. The "most reactive mixture fraction" is found to vary with time when a four-step kinetics is considered. First-order CMC results are found to be qualitatively wrong if the conditional mean scalar dissipation rate is not modeled properly. The autoignition delay time predicted by the CMC model compares excellently with DNS results and shows a trend similar to experimental data over a range of initial temperatures.

Mean and fluctuating pressure in boat-tail separated flows at transonic speeds

Experiments were carried out investigating the features of mean and unsteady surface pressure fluctuations in boat-tail separated flows relevant to launch vehicle configurations at transonic speeds. The tests were performed on a generic axisymmetric body in the Mach-number range of 0.7-1.2, and the important geometrical parameters, namely, the boat-tail angle and diameter ratio, were varied systematically. The measurements made included primarily the mean and unsteady surface-pressure fluctuations on nine different model configurations. Flow-visualization studies employing a surface oil flow, and schlieren techniques were carried out to infer features like boundary-layer separation, reattachment, and shock waves in the flow. The features of mean and fluctuating surface pressures are discussed in detail including aspects of similarity. It has been observed that, on a generic configuration employed in the present study, the maximum levels of surface-pressure fluctuations in the reattachment zone are appreciably lower than those found on launch vehicle configurations having a bulbous or hammerhead nose shape. A simple correlation is suggested for the maximum value of rms pressure fluctuations in the reattachment zone at different freestream Mach numbers.

Turbulent near-wake studies behind an infinitely swept wing

The three-dimensional asymmetric turbulent near-Rake behind an infinitely swept wing with GAW(2) airfoil has been investigated at low speeds. The near-wake in the present study is asymmetric because the boundary layers on the top and bottom surfaces of the model develop under different streamwise pressure gradients. Distributions of mean velocity, three turbulent normal stresses, and two important Reynolds shear stresses have been measured using hot-wire anemometry. The profiles of mean velocity and Reynolds shear stress exhibit asymmetry near the trailing edge and seem to have become symmetric within a short distance of 6 trailing edge momentum thicknesses. Results of computation using K-epsilon turbulence model with a simple scheme to predict the near-wake behind the swept wing have also been presented and compared with the experimental data. The agreement of the predicted mean How development with the experiment is fair considering the simplicity of the scheme.

Modified kinetic flux vector splitting (m-KFVS) method for compressible flows

To resolve many flow features accurately, like accurate capture of suction peak in subsonic flows and crisp shocks in flows with discontinuities, to minimise the loss in stagnation pressure in isentropic flows or even flow separation in viscous flows require an accurate and low dissipative numerical scheme. The first order kinetic flux vector splitting (KFVS) method has been found to be very robust but suffers from the problem of having much more numerical diffusion than required, resulting in inaccurate computation of the above flow features. However, numerical dissipation can be reduced by refining the grid or by using higher order kinetic schemes. In flows with strong shock waves, the higher order schemes require limiters, which reduce the local order of accuracy to first order, resulting in degradation of flow features in many cases. Further, these schemes require more points in the stencil and hence consume more computational time and memory. In this paper, we present a low dissipative modified KFVS (m-KFVS) method which leads to improved splitting of inviscid fluxes. The m-KFVS method captures the above flow features more accurately compared to first order KFVS and the results are comparable to second order accurate KFVS method, by still using the first order stencil. (C) 2011 Elsevier Ltd. All rights reserved.

Equilibrium and relaxation in turbulent wakes

In order to study the memory of the larger eddies in turbulent shear flow, experiments have been conducted on plane turbulent wakes undergoing transition from an initial (carefully prepared) equilibrium state to a different final one, as a result of a nearly impulsive pressure gradient. It is shown that under the conditions of the experiments the equations of motion possess self-preserving solutions in the sense of Townsend (1956), but the observed behaviour of the wake is appreciably different when the pressure gradient is not very small, as the flow goes through a slow relaxation process before reaching final equilibrium. Measurements of the Reynolds stresse show that the approach to a new equilibrium state is exponential, with a relaxation length of the order of 103 momentum thicknesses. It is suggested that a flow satisfying the conditions required by a self-preservation analysis will exhibit equilibrium only if the relaxation length is small compared with a characteristic streamwise length scale of the flow.

Turbulent non-equlibrium wakes

We consider here the detailed application of a model Reynolds stress equation (Narasimha 1969) to plane turbulent wakes subjected to pressure gradients. The model, which is a transport equation for the stress exhibiting relaxation and diffusion, is found to be consistent with the observed response of a wake to a nearly impulsive pressure gradient (Narasimha & Prabhu 1971). It implies in particular that a wake can be in equilibrium only if the longitudinal strain rate is appreciably less than the wake shear. We then describe a further series of experiments, undertaken to investigate the range of validity of the model. It is found that, with an appropriate convergence correction when necessary, the model provides excellent predictions of wake development under favourable, adverse and mixed pressure gradients. Furthermore, the behaviour of constant-pressure distorted wakes, as reported by Keffer (1965, 1967), is also explained very well by the model when account is taken of the effective flow convergence produced by the distortion. In all these calculations, only a simple version of the model is used, involving two non-dimensional constants both of which have been estimated from a single relaxation experiment.

The ‘bursting’ phenomenon in a turbulent boundary layer

Using a hot wire in a turbulent boundary layer in air, an experimental study has been made of the frequent periods of activity (to be called ‘bursts’) noticed in a turbulent signal that has been passed through a narrow band-pass filter. Although definitive identification of bursts presents difficulties, it is found that a reasonable characteristic value for the mean interval between such bursts is consistent, at the same Reynolds number, with the mean burst periods measured by Kline et al. (1967), using hydrogen-bubble techniques in water. However, data over the wider Reynolds number range covered here show that, even in the wall or inner layer, the mean burst period scales with outer rather than inner variables; and that the intervals are distributed according to the log normal law. It is suggested that these ‘bursts’ are to be identified with the ‘spottiness’ of Landau & Kolmogorov, and the high-frequency intermittency observed by Batchelor & Townsend. It is also concluded that the dynamics of the energy balance in a turbulent boundary layer can be understood only on the basis of a coupling between the inner and outer layers.

Transition to turbulence in the free convection boundary layers on an inclined heated plate

An experimental study has been made of transition to turbulence in the free convective flows on a heated plate. Observations have been made with the plate vertical and inclined at angles up to about 50° to the vertical, both above and below the plate. A fibre anemometer was used to survey the region of intermittent turbulence. Information has thus been obtained about the range of Grashof numbers over which transition takes place. Even when the plate is vertical the region of intermittent turbulence is long. When it is inclined, this region becomes still longer in the flow below the plate as a result of the stabilizing stratification, a Richardson number effect. It is possible to have a whole flow such that it should be described as transitional, not laminar or turbulent. It was noticed that in this flow and the vertical plate one, the velocity during the laminar periods could be either of two characteristic values, one of them close to zero. The behaviour above an inclined plate could be interpreted largely as a trend towards the behaviour described in a preceding paper.