A one point shock capturing kinetic scheme for hyperbolic conservation laws

We have presented a new low dissipative kinetic scheme based on a modified Courant Splitting of the molecular velocity through a parameter φ. Conditions for the split fluxes derived based on equilibrium determine φ for a one point shock. It turns out that φ is a function of the Left and Right states to the shock and that these states should satisfy the Rankine-Hugoniot Jump condition. Hence φ is utilized in regions where the gradients are sufficiently high, and is switched to unity in smooth regions. Numerical results confirm a discrete shock structure with a single interior point when the shock is aligned with the grid.

Application of Artificial Viscosity in Establishing Supercritical Solutions to the Transonic Integra

The nonlinear singular integral equation of transonic flow is examined in the free-stream Mach number range where only solutions with shocks are known to exist. It is shown that, by the addition of an artificial viscosity term to the integral equation, even the direct iterative scheme, with the linear solution as the initial iterate, leads to convergence. Detailed tables indicating how the solution varies with changes in the parameters of the artificial viscosity term are also given. In the best cases (when the artificial viscosity is smallest), the solutions compare well with known results, their characteristic feature being the representation of the shock by steep gradients rather than by abrupt discontinuities. However, 'sharp-shock solutions' have also been obtained by the implementation of a quadratic iterative scheme with the 'artificial viscosity solution' as the initial iterate; the converged solution with a sharp shock is obtained with only a few more iterates. Finally, a review is given of various shock-capturing and shock-fitting schemes for the transonic flow equations in general, and for the transonic integral equation in particular, frequent comparisons being made with the approach of this paper.

Interface capturing simulations of acoustically driven bubble dynamics in and near tissue

Tissue injury during therapeutic ultrasound or lithotripsy is thought, in cases, to be due to the action of cavitation bubbles. Assessing this and mitigating it is challenging since bubble dynamics in the complex confinement of tissues or in small blood vessels are challenging to predict. Simulations tools require specialized algorithms to simultaneously represent strong acoustic waves and shocks, topologically complex liquid‐vapor phase boundaries, and the complex viscoelastic material dynamics of tissue. We discuss advances in a simulation tool for such situations. A single‐mesh Eulerian solver is used to solve the governing equations. Special sharpening terms maintain the liquid‐vapor interface in face of the finite numerical dissipation included in the scheme to accurately capture shocks. A recent enhancement to this formulation has significantly improved this interface capturing procedure, which is demonstrated for simulation of the Rayleigh collapse of a bubble. The solver also transports elastic stresses and can thus be used to assess the effects of elastic properties on bubble dynamics. A shock‐induced bubble collapse adjacent to a model elastic tissue is used to demonstrate this and draw some conclusions regarding the injury suppressing role that tissue elasticity might play.

Nonlinear preconditioning for efficient and accurate interface capturing in simulation of multicomponent compressible flows

Single fluid schemes that rely on an interface function for phase identification in multicomponent compressible flows are widely used to study hydrodynamic flow phenomena in several diverse applications. Simulations based on standard numerical implementation of these schemes suffer from an artificial increase in the width of the interface function owing to the numerical dissipation introduced by an upwind discretization of the governing equations. In addition, monotonicity requirements which ensure that the sharp interface function remains bounded at all times necessitate use of low-order accurate discretization strategies. This results in a significant reduction in accuracy along with a loss of intricate flow features. In this paper we develop a nonlinear transformation based interface capturing method which achieves superior accuracy without compromising the simplicity, computational efficiency and robustness of the original flow solver. A nonlinear map from the signed distance function to the sigmoid type interface function is used to effectively couple a standard single fluid shock and interface capturing scheme with a high-order accurate constrained level set reinitialization method in a way that allows for oscillation-free transport of the sharp material interface. Imposition of a maximum principle, which ensures that the multidimensional preconditioned interface capturing method does not produce new maxima or minima even in the extreme events of interface merger or breakup, allows for an explicit determination of the interface thickness in terms of the grid spacing. A narrow band method is formulated in order to localize computations pertinent to the preconditioned interface capturing method. Numerical tests in one dimension reveal a significant improvement in accuracy and convergence; in stark contrast to the conventional scheme, the proposed method retains its accuracy and convergence characteristics in a shifted reference frame. Results from the test cases in two dimensions show that the nonlinear transformation based interface capturing method outperforms both the conventional method and an interface capturing method without nonlinear transformation in resolving intricate flow features such as sheet jetting in the shock-induced cavity collapse. The ability of the proposed method in accounting for the gravitational and surface tension forces besides compressibility is demonstrated through a model fully three-dimensional problem concerning droplet splash and formation of a crownlike feature. (C) 2014 Elsevier Inc. All rights reserved.

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.

Heat shock protein 70 on Neuro2a cells is a putative receptor for Japanese encephalitis virus

Japanese encephalitis virus (JEV) envelope (E) protein has been shown to play a critical role in attachment to cells. However, the receptor interacting with envelope protein has not been conclusively identified. Using mouse neuroblastoma (Neuro2a) cells and purified JEV-E protein in `Virus Overlay Protein Binding Assay' followed by MALDI-TOF analysis, we identified `heat shock protein 70' (Hsp70) as a possible receptor for JEV. Indirect immunofluorescence and flow-cytometry analysis demonstrated localization of Hsp70 on Neuro2a cell surface. Co-immunoprecipitation followed by Western blot analysis reconfirmed the interaction between Hsp70 and JEV-E protein. Further, anti-Hsp70 polyclonal-antibodies were able to block JEV entry into Neuro2a cells. Additionally, using the bioinformatic tool - FTDOCK, clocking between the proteins was performed. Amongst six interacting structural poses studied one pose involving RGD motif on JEV-E and leucine(539) on Hsp70 displayed stable interaction. These observations indicate that Hsp70 serves as putative receptor for JEV in Neuro2A cells.

Analysis of an internally mountable accelerometer balance system for use with non-isotropic models in shock tunnels

Knowledge of drag force is an important design parameter in aerodynamics. Measurement of aerodynamic forces at hypersonic speed is a challenge and usually ground test facilities like shock tunnels are used to carry out such tests. Accelerometer based force balances are commonly employed for measuring aerodynamic drag around bodies in hypersonic shock tunnels. In this study, we present an analysis of the effect of model material on the performance of an accelerometer balance used for measurement of drag in impulse facilities. From the experimental studies performed on models constructed out of Bakelite HYLEM and Aluminum, it is clear that the rigid body assumption does not hold good during the short testing duration available in shock tunnels. This is notwithstanding the fact that the rubber bush used for supporting the model allows unconstrained motion of the model during the short testing time available in the shock tunnel. The vibrations induced in the model on impact loading in the shock tunnel are damped out in metallic model, resulting in a smooth acceleration signal, while the signal become noisy and non-linear when we use non-isotropic materials like Bakelite HYLEM. This also implies that careful analysis and proper data reduction methodologies are necessary for measuring aerodynamic drag for non-metallic models in shock tunnels. The results from the drag measurements carried out using a 60 degrees half angle blunt cone is given in the present analysis.

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.

On the propagation of a multi-dimensional shock of arbitrary strength

In this paper the kinematics of a curved shock of arbitrary strength has been discussed using the theory of generalised functions. This is the extension of Moslov’s work where he has considered isentropic flow even across the shock. The condition for a nontrivial jump in the flow variables gives the shock manifold equation (sme). An equation for the rate of change of shock strength along the shock rays (defined as the characteristics of the sme) has been obtained. This exact result is then compared with the approximate result of shock dynamics derived by Whitham. The comparison shows that the approximate equations of shock dynamics deviate considerably from the exact equations derived here. In the last section we have derived the conservation form of our shock dynamic equations. These conservation forms would be very useful in numerical computations as it would allow us to derive difference schemes for which it would not be necessary to fit the shock-shock explicitly.

Structure of shock waves at re-entry speeds

The unified structure of steady, one-dimensional shock waves in argon, in the absence of an external electric or magnetic field, is investigated. The analysis is based on a two-temperature, three-fluid continuum approach, using the Navier—Stokes equations as a model and including non-equilibrium collisional as well as radiative ionization phenomena. Quasi charge neutrality and zero velocity slip are assumed. The integral nature of the radiative terms is reduced to analytical forms through suitable spectral and directional approximations. The analysis is based on the method of matched asymptotic expansions. With respect to a suitably chosen small parameter, which is the ratio of atom-atom elastic collisional mean free-path to photon mean free-path, the following shock morphology emerges: within the radiation and electron thermal conduction dominated outer layer occurs an optically transparent discontinuity which consists of a chemically frozen heavy particle (atoms and ions) shock and a collisional ionization relaxation layer. Solutions are obtained for the first order with respect to the small parameter of the problem for two cases: (i) including electron thermal conduction and (ii) neglecting it in the analysis of the outer layer. It has been found that the influence of electron thermal conduction on the shock structure is substantial. Results for various free-stream conditions are presented in the form of tables and figures.

Application of shock waves in pencil manufacturing industry

A new shock wave generator has been designed, fabricated and tested for preservative impregnation studies into wood slats used for manufacturing pencils in the Shock Waves Laboratory, IISc, Bangalore. Series of experiments have been carried out in the laboratory to achieve satisfactory preservative impregnation into VATTA wood slats. The experiments have shown that it is indeed possible to impregnate preservatives into VATTA wood slats using shock waves and the depth of penetration and the retention of preservatives by wood slats is as good as the conventional methods. This method is expected to result in substantial reduction in the treatment process time compared to conventional methods that are currently being used by the pencil manufacturing industry.

Material processing and surface reaction studies in free piston driven shock tube

Recently established moderate size free piston driven hypersonic shock tunnel HST3 along with its calibration is described here. The extreme thermodynamic conditions prevalent behind the reflected shock wave have been utilized to study the catalytic and non-catalytic reactions of shock heated test gases like Ar, N2 or O2 with different material like C60 carbon, zirconia and ceria substituted zirconia. The exposed test samples are investigated using different experimental methods. These studies show the formation of carbon nitride due to the non-catalytic interaction of shock heated nitrogen gas with C60 carbon film. On the other hand, the ZrO2 undergoes only phase transformation from cubic to monoclinic structure and Ce0.5Zr0.5O2 in fluorite cubic phase changes to pyrochlore (Ce2Zr2O7±δ) phase by releasing oxygen from the lattice due to heterogeneous catalytic surface reaction.

Experimental investigation of catalytic and non-catalytic surface reactions on SiO2 thin films with shock heated oxygen gas

The present study is to investigate the interaction of strong shock heated oxygen on the surface of SiO2 thin film. The thermally excited oxygen undergoes a three-body recombination reaction on the surface of silicon dioxide film. The different oxidation states of silicon species on the surface of the shock-exposed SiO2 film are discussed based on X-ray Photoelectron Spectroscopy (XPS) results. The surface morphology of the shock wave induced damage at the cross section of SiO2 film and structure modification of these materials are analyzed using scanning electron microscopy and ion microscopy. Whether the surface reaction of oxygen on SiO2 film is catalytic or non-catalytic is discussed in this paper.

Experimental investigation of interaction of strong shock heated oxygen gas on the surface of ZrO2 - a novel method to understand re-combination heating

Interaction of shock heated test gas in the free piston driven shock tube with bulk and thin film of cubic zirconium dioxide (ZrO2) prepared by combustion method is investigated. The test samples before and after exposure to the shock wave are analyzed by X-ray diffraction (XRD), X-Ray Photoelectron Spectroscopy (XPS) and Scanning Electron Microscope (SEM). The study shows transformation of metastable cubic ZrO2 to stable monoclinic ZrO2 phase after interacting with shock heated oxygen gas due to the heterogeneous catalytic recombination surface reaction.