Mesh-free approximations via the error reproducing kernel method and applications to nonlinear systems developing shocks

Error estimates for the error reproducing kernel method (ERKM) are provided. The ERKM is a mesh-free functional approximation scheme [A. Shaw, D. Roy, A NURBS-based error reproducing kernel method with applications in solid mechanics, Computational Mechanics (2006), to appear (available online)], wherein a targeted function and its derivatives are first approximated via non-uniform rational B-splines (NURBS) basis function. Errors in the NURBS approximation are then reproduced via a family of non-NURBS basis functions, constructed using a polynomial reproduction condition, and added to the NURBS approximation of the function obtained in the first step. In addition to the derivation of error estimates, convergence studies are undertaken for a couple of test boundary value problems with known exact solutions. The ERKM is next applied to a one-dimensional Burgers equation where, time evolution leads to a breakdown of the continuous solution and the appearance of a shock. Many available mesh-free schemes appear to be unable to capture this shock without numerical instability. However, given that any desired order of continuity is achievable through NURBS approximations, the ERKM can even accurately approximate functions with discontinuous derivatives. Moreover, due to the variation diminishing property of NURBS, it has advantages in representing sharp changes in gradients. This paper is focused on demonstrating this ability of ERKM via some numerical examples. Comparisons of some of the results with those via the standard form of the reproducing kernel particle method (RKPM) demonstrate the relative numerical advantages and accuracy of the ERKM.

Pressure shocks in relativistic viscous heat-conducting fluids

In this paper we have studied the propagation of pressure shocks in viscous, heat-conducting, relativistic fluids. Velocities of wave fronts and growth equations for the strength of the waves are obtained in the case of low and high temperatures with variable transport coefficients. On the basis of numerical integrations the growth equation results have been discussed. In the case of constant transport coefficients and for all admissible values of ratio of specific heats of the fluid, an analytical solution for the velocity of the wave as a function of distance along the normal trajectory to the wave front, has been obtained.

Some Exact-Solutions Describing Unsteady Plane Gas-Flows With Shocks

A new class of exact solutions of plane gasdynamic equations is found which describes piston-driven shocks into non-uniform media. The governing equations of these flows are taken in the coordinate system used earlier by Ustinov, and their similarity form is determined by the method of infinitesimal transformations. The solutions give shocks with velocities which either decay or grown in a finite or infinite time depending on the density distribution in the ambient medium, although their strength remains constant. The results of the present study are related to earlier investigations describing the propagation of shocks of constant strength into non-uniform media.

A model of the polarization position-angle swings in BL Lacertae objects

The polarization position-angle swings that have been measured in a number of BL Lacertae objects and highly variable quasars are interpreted in terms of shock waves which illuminate (by enhanced synchrotron radiation) successive transverse cross sections of a magnetized, relativistic jet. The jet is assumed to have a nonaxisymmetric magnetic field configuration of the type discussed in the companion paper on the equilibria of force-free jets. For a jet that is viewed at a small angle to the axis, the passage of a shock will give rise to an apparent rotation of the polarization position angle whose amplitude can be substantially larger than 180 deg. The effects of freely propagating shocks are compared with those of bow shocks which form in front of dense obstacles in the jet, and specific applications to 0727 - 115 and BL Lacertae are considered. In the case of 0727 - 115, it is pointed out that the nonuniformity of the swing rate and the apparent oscillations of the degree of polarization could be a consequence of relativistic aberration.

Absorption of Electromagnetic Waves in Astrophysical Plasmas

We study parametric Decay Instabilities (PDI) using the kinetic description, in the homogeneous and unmagnetic plasmas. These instabilities cause anomalous absorption of the incident electromagnetic (e.m) radiations. The maximum plasma temperatures reached are functionas of luminocity of the non-thermal radio radiation and the plasma parameters.

Radio wave heating of astrophysical plasmas

The coherent plasma process such as parametric decay instability (PDI) has been applied to a homogeneous and unmagnetized plasma. These instabilities cause anomalous absorption of strong electromagnetic radiation under specific conditions of energy and momentum conservation and thus cause anomalous heating of the plasma. The maximum plasma temperatures reached are functions of luminosity of the radio radiation and plasma parameters. We believe that these processes may be taking place in many astrophysical objects. Here, the conditions in the sources 3C 273, 3C 48 and Crab Nebula are shown to be conducive to the excitation of PDI. These processes also contribute towards the absorption of 21cm radiation

Can the viscosity in astrophysical black hole accretion disks be close to its string theory bound?

String theory and gauge/gravity duality suggest the lower bound of shear viscosity (eta) to entropy density (s) for any matter to be mu h/4 pi k(B), when h and k(B) are reduced Planck and Boltzmann constants respectively and mu <= 1. Motivated by this, we explore eta/s in black hole accretion flows, in order to understand if such exotic flows could be a natural site for the lowest eta/s. Accretion flow plays an important role in black hole physics in identifying the existence of the underlying black hole. This is a rotating shear flow with insignificant molecular viscosity, which could however have a significant turbulent viscosity, generating transport, heat and hence entropy in the flow. However, in presence of strong magnetic field, magnetic stresses can help in transporting matter independent of viscosity, via celebrated Blandford-Payne mechanism. In such cases, energy and then entropy produces via Ohmic dissipation. In,addition, certain optically thin, hot, accretion flows, of temperature greater than or similar to 10(9) K, may be favourable for nuclear burning which could generate/absorb huge energy, much higher than that in a star. We find that eta/s in accretion flows appears to be close to the lower bound suggested by theory, if they are embedded by strong magnetic field or producing nuclear energy, when the source of energy is not viscous effects. A lower bound on eta/s also leads to an upper bound on the Reynolds number of the flow.

Transonic flutter, shocks, and parameter effects

There is a drop in the flutter boundary of an aeroelastic system placed in a transonic flow due to compressibility effects and is known as the transonic dip. Viscous effects can shift the lo-cation of the shock and depending on the shock strength the boundary layer may separate leading to changes in the flutter speed. An unsteady Euler flow solver coupled with the structural dynamic equations is used to understand the effect of shock on the transonic dip. The effect of various system parameters such as mass ratio, location of the center of mass, position of the elastic axis, ratio of uncoupled natural frequencies in heave and pitch are also studied. Steady turbulent flow results are presented to demonstrate the effect of viscosity on the location and strength of the shock.

Riemann shock tube: 1D normal shocks in air, simulations and experiments

This paper presents numerical simulation of the evolution of one-dimensional normal shocks, their propagation, reflection and interaction in air using a single diaphragm Riemann shock tube and validate them using experimental results. Mathematical model is derived for one-dimensional compressible flow of viscous and conducting medium. Dimensionless form of the mathematical model is used to construct space-time finite element processes based on minimization of the space-time residual functional. The space-time local approximation functions for space-time p-version hierarchical finite elements are considered in higher order GRAPHICS] spaces that permit desired order of global differentiability of local approximations in space and time. The resulting algebraic systems from this approach yield unconditionally positive-definite coefficient matrices, hence ensure unique numerical solution. The evolution is computed for a space-time strip corresponding to a time increment Delta t and then time march to obtain the evolution up to any desired value of time. Numerical studies are designed using recently invented hand-driven shock tube (Reddy tube) parameters, high/low side density and pressure values, high- and low-pressure side shock tube lengths, so that numerically computed results can be compared with actual experimental measurements.

A stabilization approach for mesh-free simulations of systems developing shocks or extreme strain localizations

A new stabilization scheme, based on a stochastic representation of the discretized field variables, is proposed with a view to reduce or even eliminate unphysical oscillations in the mesh-free numerical simulations of systems developing shocks or exhibiting localized bands of extreme deformation in the response. The origin of the stabilization scheme may be traced to nonlinear stochastic filtering and, consistent with a class of such filters, gain-based additive correction terms are applied to the simulated solution of the system, herein achieved through the element-free Galerkin method, in order to impose a set of constraints that help arresting the spurious oscillations. The method is numerically illustrated through its Applications to inviscid Burgers' equations, wherein shocks may develop as a result of intersections of the characteristics, and to a gradient plasticity model whose response is often characterized by a developing shear band as the external load is gradually increased. The potential of the method in stabilized yet accurate numerical simulations of such systems involving extreme gradient variations in the response is thus brought forth. (C) 2014 Elsevier Ltd. 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.

Higher-Order Nonlinearity In Accretion Disks: Quasi-Periodic Oscillations Of Black Hole And Neutron Star Sources And Their Spin

We propose a unified model to explain Quasi-Periodic Oscillation (QPO), particularly of high frequency, observed from black hole and neutron star systems globally. We consider accreting systems to be damped harmonic oscillators exhibiting epicyclic oscillations with higher-order nonlinear resonance to explain QPO. The resonance is expected to be driven by the disturbance from the compact object at its spin frequency. The model explains various properties parallelly for both types of the compact object. It describes QPOs successfully for ten different compact sources. Based on this, we predict the spin frequency of the neutron star Sco X-1 and specific angular momentum of black holes GRO J1655–40, XTE J1550–564, H1743–322, and GRS 1915+105.

Electrical Characterization of Airborne Vehicle Exhaust Plume

he induced current and voltage on the skin of an airborne vehicle due to the coupling of external electromagnetic field could be altered in the presence of ionized exhaust plume. So in the present work, a theoretical analysis is done to estimate the electrical parameters such as electrical conductivity and permittivity and their distribution in the axial and radial directions of the exhaust plume of an airborne vehicle. The electrical conductivity depends on the distribution of the major ionic species produced from the propellant combustion. In addition it also depends on temperature and pressure distribution of the exhaust plume as well as the generated shock wave. The chemically reactive rocket exhaust flow is modeled in two stages. The first part is simulated from the combustion chamber to the throat of the supersonic nozzle by using NASA Chemical Equilibrium with Application (CEA) package and the second part is simulated from the nozzle throat to the downstream of the plume by using a commercial Computational Fluid Dynamics (CFD) solver. The contour plots of the exhaust parameters are presented. Eight barrel shocks which influence the distribution of the vehicle exhaust parameters are obtained in this simulation. The computed peak value of the electrical conductivity of the plume is 0.123 S/m and the relative permittivity varies from 0.89 to 0.99. The attenuation of the microwave when it is passing through the conducting exhaust plume has also been presented.

Massive black-hole binary inspirals: results from the LISA parameter estimation taskforce

The LISA Parameter Estimation Taskforce was formed in September 2007 to provide the LISA Project with vetted codes, source distribution models and results related to parameter estimation. The Taskforce's goal is to be able to quickly calculate the impact of any mission design changes on LISA's science capabilities, based on reasonable estimates of the distribution of astrophysical sources in the universe. This paper describes our Taskforce's work on massive black-hole binaries (MBHBs). Given present uncertainties in the formation history of MBHBs, we adopt four different population models, based on (i) whether the initial black-hole seeds are small or large and (ii) whether accretion is efficient or inefficient at spinning up the holes. We compare four largely independent codes for calculating LISA's parameter-estimation capabilities. All codes are based on the Fisher-matrix approximation, but in the past they used somewhat different signal models, source parametrizations and noise curves. We show that once these differences are removed, the four codes give results in extremely close agreement with each other. Using a code that includes both spin precession and higher harmonics in the gravitational-wave signal, we carry out Monte Carlo simulations and determine the number of events that can be detected and accurately localized in our four population models.