Composite scheme using localized relaxation with non-standard finite difference method for hyperbolic conservation laws

Non-standard finite difference methods (NSFDM) introduced by Mickens [Non-standard Finite Difference Models of Differential Equations, World Scientific, Singapore, 1994] are interesting alternatives to the traditional finite difference and finite volume methods. When applied to linear hyperbolic conservation laws, these methods reproduce exact solutions. In this paper, the NSFDM is first extended to hyperbolic systems of conservation laws, by a novel utilization of the decoupled equations using characteristic variables. In the second part of this paper, the NSFDM is studied for its efficacy in application to nonlinear scalar hyperbolic conservation laws. The original NSFDMs introduced by Mickens (1994) were not in conservation form, which is an important feature in capturing discontinuities at the right locations. Mickens [Construction and analysis of a non-standard finite difference scheme for the Burgers–Fisher equations, Journal of Sound and Vibration 257 (4) (2002) 791–797] recently introduced a NSFDM in conservative form. This method captures the shock waves exactly, without any numerical dissipation. In this paper, this algorithm is tested for the case of expansion waves with sonic points and is found to generate unphysical expansion shocks. As a remedy to this defect, we use the strategy of composite schemes [R. Liska, B. Wendroff, Composite schemes for conservation laws, SIAM Journal of Numerical Analysis 35 (6) (1998) 2250–2271] in which the accurate NSFDM is used as the basic scheme and localized relaxation NSFDM is used as the supporting scheme which acts like a filter. Relaxation schemes introduced by Jin and Xin [The relaxation schemes for systems of conservation laws in arbitrary space dimensions, Communications in Pure and Applied Mathematics 48 (1995) 235–276] are based on relaxation systems which replace the nonlinear hyperbolic conservation laws by a semi-linear system with a stiff relaxation term. The relaxation parameter (λ) is chosen locally on the three point stencil of grid which makes the proposed method more efficient. This composite scheme overcomes the problem of unphysical expansion shocks and captures the shock waves with an accuracy better than the upwind relaxation scheme, as demonstrated by the test cases, together with comparisons with popular numerical methods like Roe scheme and ENO schemes.

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.

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.

Hypersonic wave drag reduction in re-entry capsules using concentrated energy deposition

Experiments are carried out in a shock tunnel at a nominal Mach number of 5.75 in order to study the effect of concentrated energy deposition on the drag force experienced by a 120° blunt cone. Electrical energy was deposited along the stagnation streamline of the model using a high voltage DC discharge circuit (1.5 – 3.5KW) and the drag force was measured by a single component accelerometer balance. Numerical simulations were also carried complimenting the experiments. These simulations showed a substantial drag reduction (20% ~ 65%) whereas the experiments show no appreciable reduction in drag

Effect of gas injection on drag and surface heat transfer rates for a 30 degrees semi-apex angle blunt body flying at Mach 5.75

Effect of coolant gas injection in the stagnation region on the surface heat transfer rates and aerodynamic drag for a large angle blunt body flying at hypersonic Mach number is reported for two stagnation enthalpies. A 60° apex-angle blunt cone model is employed for this purpose with air injection at the nose through a hole of 2mm diameter. The convective surface heating rates and aerodynamic drag are measured simultaneously using surface mounted platinum thin film sensors and internally mounted accelerometer balance system, respectively. About 35–40% reduction in surface heating rates is observed in the vicinity of stagnation region whereas 15–25% reduction in surface heating rates is felt beyond the stagnation region at stagnation enthalpy of 1.6MJ/kg. The aerodynamic drag expressed in terms of drag coefficient is found to increase by 20% due to the air injection.

A new 6-component accelerometer force balance for short duration ground testing facilities

A new six-component accelerometer force balance is developed and used in the HST2 shock tunnel of Indian Institute of Science. Aerodynamic forces and moments for a hypersonic slender body measured using this balance system at a free stream Mach number of 5.75 and Reynolds number of 1.5 million and stagnation enthalpy of 1.5 and 2 MJ/kg are presented. These measured values compare well with the theoretical values estimated using modified Newtonian theory.

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.

On the propagation of a weak shock front: Theory and application

In this paper the kinematics of a weak shock front governed by a hyperbolic system of conservation laws is studied. This is used to develop a method for solving problems, involving the propagation of nonlinear unimodal waves. It consists of first solving the nonlinear wave problem by moving along the bicharacteristics of the system and then fitting the shock into this solution field, so that it satisfies the necessary jump conditions. The kinematics of the shock leads in a natural way to the definition of ldquoshock-raysrdquo, which play the same role as the ldquoraysrdquo in a continuous flow. A special case of a circular cylinder introduced suddenly in a constant streaming flow is studied in detail. The shock fitted in the upstream region propagates with a velocity which is the mean of the velocities of the linear and the nonlinear wave fronts. In the downstream the solution is given by an expansion wave.

Singular perturbation analysis of spherical and cylindrical N waves

Using a singular perturbation analysis the nonplanar Burgers' equation is solved to yield the shock wave-displacement due to diffusion for spherical and cylindrical N waves, thus supplementing the earlier results of Lighthill for the plane N waves. Physics of Fluids is copyrighted by The American Institute of Physics.

Scattering of electromagnetic waves from chirally coated cylinders

The problem of electromagnetic scattering from an isotropic homogeneous chirally coated conducting cylinder is analysed. The cylinder is assumed to be illuminated by either a transverse magnetic or a transverse electric wave. Mie's analysis is used to evaluate the scattering characteristics. The computed results include the evaluation of the normalized scattering width and the absorption efficiency. The results show that there is a significant reduction in the normalized scattering width as compared to a RAM coated cylinder. This reduction has been attributed to increased absorption.

Shock propagation in gas dynamics: explicit form of higher order compatibility conditions

With the use of tensor analysis and the method of singular surfaces, an infinite system of equations can be derived to study the propagation of curved shocks of arbitrary strength in gas dynamics. The first three of these have been explicitly given here. This system is further reduced to one involving scalars only. The choice of dependent variables in the infinite system is quite important, it leads to coefficients free from singularities for all values of the shock strength.

Thermal shock characteristics of plasma sprayed mullite coatings

Commercially available mullite (3Al(2)O(3). 2SiO(2)) powders containing oxides of calcium and iron as impurities, have been made suitable for plasma spraying by using an organic binder. Stainless steel substrates covered with Ni-22Cr-10Al-1.0Y bond coat were spray coated with mullite, The 425 mu m thick coatings were subjected to thermal shock cycling under burner rig conditions between 1000 and 1200 degrees C and less than 200 degrees C with holding times of 1, 5, and 30 min. While the coatings withstood as high as 1000 shock cycles without failure between 1000 and 200 degrees C, spallation occurred early at 120 cycles when shocked from 1200 degrees C, The coatings appeared to go through a process of self erosion at high temperatures resulting in loss of material. Also observed were changes attributable to melting of the silicate grains, which smooth down the surface. Oxidation of the bond coat did not appear to influence the failure, These observations were supported by detailed scanning electron microscopy and quantitative chemical composition analysis, differential thermal analysis, and surface roughness measurements.

Diffraction Of Elastic Waves By Two Parallel Rigid Strips Embedded In An Infinite Orthotropic Medium

The elastodynamic response of a pair of parallel rigid strips embedded in an infinite orthotropic medium due to elastic waves incident normally on the strips has been investigated. The mixed boundary value problem has been solved by the Integral Equation method. The normal stress and the vertical displacement have been derived in closed form. Numerical values of stress intensity factors at inner and outer edges of the strips and vertical displacement at points in the plane of the strips for several orthotropic materials have been calculated and plotted graphically to show the effect of material orthotropy.