Numerical prediction of nozzle flow separation: Issue of turbulence modeling

Numerical simulation of separated flows in rocket nozzles is challenging because existing turbulence models are unable to predict it correctly. This paper addresses this issue with the Spalart-Allmaras and Shear Stress Transport (SST) eddy-viscosity models, which predict flow separation with moderate success. Their performances have been compared against experimental data for a conical and two contoured subscale nozzles. It is found that they fail to predict the separation location correctly, exhibiting sensitivity to the nozzle pressure ratio (NPR) and nozzle type. A careful assessment indicated how the model had to be tuned for better, consistent prediction. It is learnt that SST model's failure is caused by limiting of the shear stress inside boundary layer according to Bradshaw's assumption, and by over prediction of jet spreading rate. Accordingly, SST's coefficients were empirically modified to match the experimental wall pressure data. Results confirm that accurate RANS prediction of separation depends on the correct capture of the jet spreading rate, and that it is feasible over a wide range of NPRs by modified values of the diffusion coefficients in the turbulence model. (C) 2015 Elsevier Masson SAS. All rights reserved.

Flow of micropolar fluid through a tube with stenosis

mathematical model for the steady flow of non-Newtonian fluid through a stenotic region is presented. The results indicate that the general shape and size of the stenosis together with rheological properties of blood are important in understanding the flow characteristics and the presence of flow separation.

Surface pressure and heat transfer measurements in the unsteady separated hypersonic flow field over double cones

Reliable bench mark experimental database in the separated hypersonic flow regime is necessary to validate high resolution CFD codes. In this paper we report the surface pressure and heat transfer measurements carried out on double cones (first cone semi-apex angle = 15, 25 deg.; second cone semi-apex angle= 35, 68 deg.) at hypersonic speeds that will be useful for CFD code validation studies. The surface pressure measurements are carried out at nominal Mach number of 8.35 in the IISc hypersonic wind tunnel. On the other hand the surface heat transfer measurements are carried out at a nominal Mach number of 5.75 in the IISc hypersonic shock tunnel. The flow separation point on the first cone, flow reattachment on the second cone and the wild fluctuation of the transmitted shock on the second cone surface (25/68 deg. double cone) in the presence of severe adverse pressure gradient are some of the flow features captured in the measurements. The results from the CFD studies indicate good agreement with experiments in the attached flow regime while considerable differences are noticeable in the separated flow regime.

Shock tunnel study of spiked aerodynamic bodies flying at hypersonic Mach numbers

A three-component accelerometer balance system is used to study the drag reduction effect of an aerodisc on large angle blunt cones flying at hypersonic Mach numbers. Measurements in a hypersonic shock tunnel at a freestream Mach number of 5.75 indicate more than 50% reduction in the drag coefficient for a 120degrees apex angle blunt cone with a forward facing aerospike having a flat faced aerodisc at moderate angles of attack. Enhancement of drag has been observed for higher angles of attack due to the impingement of the flow separation shock on the windward side of the cone. The flowfields around the large angle blunt cone with aerospike assembly flying at hypersonic Mach numbers are also simulated numerically using a commercial CFD code. The pressure and density levels on the model surface, which is under the aerodynamic shadow of the flat disc tipped spike, are found very low and a drag reduction of 64.34% has been deduced numerically.

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.

Numerical and experimental investigation and optimization of an open wheeled race car cooling air duct

In this study the cooling performance due to air flow and aerodynamics of the Formula Student open wheeled race car has been investigated and optimized with the help of CFD simulations and experimental validation. The race car in context previously suffered from overheating problems. Flow analysis was carried out based on the detailed race car 3D model (NITK Racing 2012 formula student race car). Wind tunnel experiments were carried out on the same. The results obtained from the computer simulations are compared with experimental results obtained from wind tunnel testing of the full car. Through this study it was possible to locate the problem areas and hence choose the best configuration for the cooling duct. The CFD analysis helped in calculating the mass flow rate, pressure and velocity distribution for different velocities of the car which is then used to determine the heat dissipated by the radiator. Area of flow separation could be visualized and made sure smooth airflow into the radiator core area. This significantly increased the cooling performance of the car with reduction in drag.

Morphing fixed wing MAV modeling using VAM

The design and implementation of a morphing Micro Air Vehicle (MAV) wing using a smart composite is attempted in this research work. Control surfaces actuated by traditional servos are difficult to instrument and fabricate on thin composite-wings of MAVs. Piezoelectric Fiber Reinforced Composites (PFRCs) are the chosen smart structural materials in the current work for incorporation onto fixed-wing MAVs to simultaneously perform the dual functions of structural load-bearing and actuation of flexure, torsion and/or extension for morphing. Further, PFRC use can be extended towards shape control of a “fixed” wing MAV to meet changing performance requirements. Wings that can warp into desired shapes and/or have variable camber are well-known to exhibit improved efficiency in aerodynamic control. During an entire flight cycle, there are multiple optimal configurations, each of which suits a particular phase of the flight regime. Widely proposed methods of wing morphing include changes in camber, twist, sweep and span. However, camber change during flight is already established, in terms of its potential, as a major factor in improving the aerofoil efficiency and flow separation behavior. Hence, for this work, morphing with camber change is adopted with the goal to better tailor aerodynamic properties.

Enhanced Gas Adsorption on Graphitic Substrates via Defects and Local Curvature: A Density Functional Theory Study

Using van-der-Waals-corrected density functional theory calculations, we explore the possibility of engineering the local structure and morphology of high-surface-area graphene-derived materials to improve the uptake of methane and carbon dioxide for gas storage and sensing. We test the sensitivity of the gas adsorption energy to the introduction of native point defects, curvature, and the application of strain. The binding energy at topological point defect sites is inversely correlated with the number of missing carbon atoms, causing Stone-Wales defects to show the largest enhancement with respect to pristine graphene (similar to 20%). Improvements of similar magnitude are observed at concavely curved surfaces in buckled graphene sheets under compressive strain, whereas tensile strain tends to weaken gas binding. Trends for CO2 and CH4 are, similar, although CO2 binding is generally stronger by similar to 4 to 5 kJ mol(-1). However, the differential between the adsorption of CO2 and CH4 is much higher on folded graphene sheets and at concave curvatures; this could possibly be leveraged for CH4/CO2 flow separation and gasselective sensors.

A Ca2+-dependent early functional heterogeneity in amoebae of Dictyostelium discoideum, revealed by flow cytometry

When freshly starved amoebae of Dictyostelium discoideum are loaded with the Ca2+-specific dye indo-1/AM and analyzed in a fluorescence-activated cell sorter, they exhibit a quasi-bimodal distribution of fluorescence. This permits a separation of the population into two classes: H, or ''high Ca2+-indo-1 fluorescence,'' and L, or ''low Ca2+-indo-1 fluorescence.'' Simultaneous monitoring of Ca2+-indo-1 and Ca2+-chlortetracycline fluorescence shows that by and large the same cells tend to have high (or low) levels of both cytoplasmic and sequestered Ca2+. Next we label H cells with tetramethylrhodamine isothiocyanate (TRITC) and mix them in a 1:4 ratio with L cells, In the slugs that result, TRITC fluorescence is confined mainly to the anterior prestalk region. This implies that amoebae with relatively high Ca2+ at the vegetative stage tend to develop into prestalk cells and those with low Ca2+ into prespores. Polysphondylium violaceum, a cellular slime mold that does not possess prestalk and prespore cells, also does not display a Ca2+-dependent heterogeneity at the vegetative stage or in slugs. Finally, confirming earlier findings with the fluorophore fura-2 (Azhar ef al., Curr. Sci. 68, 337-342 (1995)), a prestalk-prespore difference in cellular Ca2+ is present in the cells of the slug in vivo. These findings are discussed in light of the possible roles of Ca2+ for cell differentiation in D. discoideum.

Effect of a magnetic field on the flow and blood oxygenation in channels of variable cross-section

The effect of a magnetic field on the flow and oxygenation of an incompressible Newtonian conducting fluid in channels with irregular boundaries has been investigated. The geometric parameter δ, which is a ratio of the mean half width of the channel d to the characteristic length λ along the channel over which the significant changes in the flow quantities occur, has been used for perturbing the governing equations. Closed form solutions of the various order equations are presented for the stream function. The equations for oxygen partial pressure remain nonlinear even after perturbation, therefore a numerical solution is presented. The expressions for shear stress at a wall and pressure distributions are derived. Here the separation in the flow occurs at a higher Reynolds number than the corresponding non-magnetic case. It is found that the magnetic field has an effect on local oxygen concentration but has a little effect on the saturation length.

Real-fluid effects in flow cavitation

The possible role of reaj fluid effects in two aspects of flow cavitation namely inception and separation is discussed. This is primarily qualitative in the case of inception whereas some quantitative results are presented in the case of separation. Existing evidence clearly indicates that in particular viscous effects can play a significant role in determining the conditions for cavitation inception and in determining the location of cavitation separation from smooth bodies.

Pulsatile flow in tubes of varying cross-sections

The pulsatile flow of an incompressible viscous fluid in a cylindrical tube of varying cross section is investigated for small Reynolds numbers. The solutions consist of a stedy and an oscillatory part. The shear stress distribution on the wall is evaluated and discussed in detail for special geometries like tapered tubes, locally constricted tubes and peristaltic tubes. The existence of separation in the flow field is noticed.

Localised instability in unsteady separation bubbles

We propose in experimental method to study the instability of thin unsteady separation bubbles, i.e. of unsteady boundary layers with reverse flow. The unsteady boundary layer is created by controlled temporal and spatial variations of the velocity external to the boundary layer. We present results of the evolution of instability in different temporally varying flows in a shallow angle diffuser. Depending on the extent of reverse flow in the boundary we observe that instability can be spatially localised.

Deformation and crystallization of Zr-based amorphous alloys in homogeneous flow regime

The purpose of this study is to experimentally investigate the interaction of inelastic deformation and microstructural changes of two Zr-based bulk metallic glasses (BMGs): Zr41.25Ti13.75Cu12.5Ni10Be22.5 (commercially designated as Vitreloy 1 or Vit1) and Zr46.75Ti8.25Cu7.5Ni10Be27.5 (Vitreloy 4, Vit4). High-temperature uniaxial compression tests were performed on the two Zr alloys at various strain rates, followed by structural characterization using differential scanning calorimetry (DSC) and transmission electron microscopy (TEM). Two distinct modes of mechanically induced atomic disordering in the two alloys were observed, with Vit1 featuring clear phase separation and crystallization after deformation as observed with TEM, while Vit4 showing only structural relaxation with no crystallization. The influence of the structural changes on the mechanical behaviors of the two materials was further investigated by jump-in-strain-rate tests, and flow softening was observed in Vit4. A free volume theory was applied to explain the deformation behaviors, and the activation volumes were calculated for both alloys.

Effect of natural ventilation on the boundary layer separation and near-wake vortex shedding characteristics of a sphere

Experiments were conducted in water and wind tunnels on spheres in the Reynolds number range 6 x 10(3) to 6.5 x 10(5) to study the effect of natural ventilation on the boundary layer separation and near-wake Vortex shedding characteristics. In the subcritical range of Re (<2 x 10(5)), ventilation caused a marginal downstream shift in the location of laminar boundary layer separation; there was only a small change in the vortex shedding frequency. In the supercritical range (Re > 4 x 10(5)), ventilation caused a downstream shift in the mean locations of boundary layer separation and reattachment; these lines showed significant axisymmetry in the presence of venting. No distinct vortex shedding frequency was found. Instead, a dramatic reduction occurred in the wake unsteadiness at all frequencies. The reduction of wake unsteadiness is consistent with the reduction in total drag already reported. Based on the present results and those reported earlier, the effects of natural ventilation on the flow past a sphere can be categorized in two broad regimes, viz., weak and strong interaction regimes. In the weak interaction regime (subcritical Re), the broad features of the basic sphere are largely unaltered despite the large addition of mass in the near wake. Strong interaction is promoted by the closer proximity of the inner and outer shear layers at supercritical Re. This results in a modified and steady near-wake flow, characterized by reduced unsteadiness and small drag.