PVU and wave-particle splitting schemes for Euler equations of gas dynamics

A new way of flux-splitting, termed as the wave-particle splitting is presented for developing upwind methods for solving Euler equations of gas dynamics. Based on this splitting, two new upwind methods termed as Acoustic Flux Vector Splitting (AFVS) and Acoustic Flux Difference Splitting (AFDS) methods are developed. A new Boltzmann scheme, which closely resembles the wave-particle splitting, is developed using the kinetic theory of gases. This method, termed as Peculiar Velocity based Upwind (PVU) method, uses the concept of peculiar velocity for upwinding. A special feature of all these methods that the unidirectional and multidirectional parts of the flux vector are treated separately. Extensive computations done using these schemes demonstrate the soundness of the ideas.

A characteristic based numerical method with tracking for nonlinear wave equations

Many physical problems can be modeled by scalar, first-order, nonlinear, hyperbolic, partial differential equations (PDEs). The solutions to these PDEs often contain shock and rarefaction waves, where the solution becomes discontinuous or has a discontinuous derivative. One can encounter difficulties using traditional finite difference methods to solve these equations. In this paper, we introduce a numerical method for solving first-order scalar wave equations. The method involves solving ordinary differential equations (ODEs) to advance the solution along the characteristics and to propagate the characteristics in time. Shocks are created when characteristics cross, and the shocks are then propagated by applying analytical jump conditions. New characteristics are inserted in spreading rarefaction fans. New characteristics are also inserted when values on adjacent characteristics lie on opposite sides of an inflection point of a nonconvex flux function, Solutions along characteristics are propagated using a standard fourth-order Runge-Kutta ODE solver. Shocks waves are kept perfectly sharp. In addition, shock locations and velocities are determined without analyzing smeared profiles or taking numerical derivatives. In order to test the numerical method, we study analytically a particular class of nonlinear hyperbolic PDEs, deriving closed form solutions for certain special initial data. We also find bounded, smooth, self-similar solutions using group theoretic methods. The numerical method is validated against these analytical results. In addition, we compare the errors in our method with those using the Lax-Wendroff method for both convex and nonconvex flux functions. Finally, we apply the method to solve a PDE with a convex flux function describing the development of a thin liquid film on a horizontally rotating disk and a PDE with a nonconvex flux function, arising in a problem concerning flow in an underground reservoir.

Applicability of Butler's equation in interpreting the thermodynamic behavior of surfaces and adsorption in Fe-S-O melts

Expressions for various second-order derivatives of surface tension with respect to composition at infinite dilution in terms of the interaction parameters of the surface and those of the bulk phases of dilute ternary melts have been presented. A method of deducing the parameters, which consists of repeated differentiation of Butler's equations with subsequent application of the appropriate boundary conditions, has been developed. The present investigation calculates the surface tension and adsorption functions of the Fe-S-O melts at 1873 and 1923 K using the modified form of Butler's equations and the derived values for the surface interaction parameters of the system. The calculated values are found to be in good agreement with those of the experimental data of the system. The present analysis indicates that the energetics of the surface phase are considerably different from those of the bulk phase. The present research investigates a critical compositional range beyond which the surface tension increases with temperature. The observed increase in adsorption of sulfur with consequent desorption of oxygen as a function of temperature above the critical compositional range has been ascribed to the increase of activity ratios of oxygen to sulfur in the surface relative to those in the bulk phase of the system.

Conjectures about phase turbulence in the complex Ginzburg-Landau equation

In the complex Ginzburg-Landau equation, we consider possible ''phase turbulent'' regimes, where asymptotic correlations are controlled by phase fluctuations rather than by topological defects. Conjecturing that the decay of such correlations is governed by the Kardar-Parisi-Zhang (KPZ) model of growing interfaces, we derive the following results: (1) A scaling ansatz implies that equal-time spatial correlations in 1d, 2d, and 3d decay like e(-Ax2 zeta), where A is a nonuniversal constant, and zeta=1/2 in 1d. (2) Temporal correlations decay as exp(-t(2 beta)h(t/L(z))), with the scaling law <(beta)over bar> = <(zeta)over bar>/z, where z = 3/2, 1.58..., and 1.66..., for d = 1,2, and 3 respectively. The scaling function h(y) approaches a constant as y --> 0, and behaves like y(2(beta-<(beta)over bar>)), for large y. If in 3d the associated KPZ model turns out to be in its weak-coupling (''smooth'') phase, then, instead of the above behavior, the CGLE exhibits rotating long-range order whose connected correlations decay like 1/x in space or 1/t(1/2) in time. (3) For system sizes, L, and times t respectively less than a crossover length, L(c), and time, t(c), correlations are governed by the free-field or Edwards-Wilkinson (EW) equation, rather than the KPZ model. In 1d, we find that L(c) is large: L(c) similar to 35,000; for L < L(c) we show numerical evidence for stretched exponential decay of temporal correlations with an exponent consistent with the EW value beta(EW)= 1/4.

Effect of open joint on stress wave propagation

Stress wave characteristics are drastically altered by joints and other inhomogenities. This paper addresses the effect of an open joint on stress wave transmission. An elastodynamic analysis is developed to supplement and explain some recent observations by Fourney and Dick(1995) on open as well as filled joints. The analytical model developed here assuming spherical symmetry can be extended to filled joints between dissimilar media, but results are presented only for open joints separating identical materials. As a special case, stress wave transmission across a joint with no gap is also addressed.

Scroll-Wave Dynamics in Human Cardiac Tissue: Lessons from a Mathematical Model with Inhomogeneities and Fiber Architecture

Cardiac arrhythmias, such as ventricular tachycardia (VT) and ventricular fibrillation (VF), are among the leading causes of death in the industrialized world. These are associated with the formation of spiral and scroll waves of electrical activation in cardiac tissue; single spiral and scroll waves are believed to be associated with VT whereas their turbulent analogs are associated with VF. Thus, the study of these waves is an important biophysical problem. We present a systematic study of the combined effects of muscle-fiber rotation and inhomogeneities on scroll-wave dynamics in the TNNP (ten Tusscher Noble Noble Panfilov) model for human cardiac tissue. In particular, we use the three-dimensional TNNP model with fiber rotation and consider both conduction and ionic inhomogeneities. We find that, in addition to displaying a sensitive dependence on the positions, sizes, and types of inhomogeneities, scroll-wave dynamics also depends delicately upon the degree of fiber rotation. We find that the tendency of scroll waves to anchor to cylindrical conduction inhomogeneities increases with the radius of the inhomogeneity. Furthermore, the filament of the scroll wave can exhibit drift or meandering, transmural bending, twisting, and break-up. If the scroll-wave filament exhibits weak meandering, then there is a fine balance between the anchoring of this wave at the inhomogeneity and a disruption of wave-pinning by fiber rotation. If this filament displays strong meandering, then again the anchoring is suppressed by fiber rotation; also, the scroll wave can be eliminated from most of the layers only to be regenerated by a seed wave. Ionic inhomogeneities can also lead to an anchoring of the scroll wave; scroll waves can now enter the region inside an ionic inhomogeneity and can display a coexistence of spatiotemporal chaos and quasi-periodic behavior in different parts of the simulation domain. We discuss the experimental implications of our study.

Performance analysis of surface acoustic wave matched filters

SAW matched filter is commonly used in spread spectrum communication receivers in order to maximize the SNR prior to detection, At times the receiver would be a mobile one while the signal is processed at the IF level, In that case frequency deviations due to Doppler shift or temperature dependence of the acoustic medium used for SAW device would, severely effect it's performance, The impact of these errors on the receiver performance is analyzed on a generalised basis.

Design and realization of high-performance wave-pipelined 8×8 b multiplier in CMOS technology

Wave pipelining is a design technique for increasing the throughput of a digital circuit or system without introducing pipelining registers between adjacent combinational logic blocks in the circuit/system. However, this requires balancing of the delays along all the paths from the input to the output which comes the way of its implementation. Static CMOS is inherently susceptible to delay variation with input data, and hence, receives a low priority for wave pipelined digital design. On the other hand, ECL and CML, which are amenable to wave pipelining, lack the compactness and low power attributes of CMOS. In this paper we attempt to exploit wave pipelining in CMOS technology. We use a single generic building block in Normal Process Complementary Pass Transistor Logic (NPCPL), modeled after CPL, to achieve equal delay along all the propagation paths in the logic structure. An 8×8 b multiplier is designed using this logic in a 0.8 ?m technology. The carry-save multiplier architecture is modified suitably to support wave pipelining, viz., the logic depth of all the paths are made identical. The 1 mm×0.6 mm multiplier core supports a throughput of 400 MHz and dissipates a total power of 0.6 W. We develop simple enhancements to the NPCPL building blocks that allow the multiplier to sustain throughputs in excess of 600 MHz. The methodology can be extended to introduce wave pipelining in other circuits as well

Diffusing wave spectroscopy of dense colloids: Liquid, crystal and glassy states

Using intensity autocorrelation of multiply scattered light, we show that the increase in interparticle interaction in dense, binary colloidal fluid mixtures of particle diameters 0.115µm and 0.089µm results in freezing into a crystalline phase at volume fraction? of 0.1 and into a glassy state at?=0.2. The functional form of the field autocorrelation functiong (1)(t) for the binary fluid phase is fitted to exp[??(6k 0 2 D eff t)1/2] wherek 0 is the magnitude of the incident light wavevector and? is a parameter inversely proportional to the photon transport mean free pathl*. TheD eff is thel* weighted average of the individual diffusion coefficients of the pure species. Thel* used in calculatingD eff was computed using the Mie theory. In the solid (crystal or glass) phase, theg (1)(t) is fitted (only with a moderate success) to exp[??(6k 0 2 W(t))1/2] where the mean-squared displacementW(t) is evaluated for a harmonically bound overdamped Brownian oscillator. It is found that the fitted parameter? for both the binary and monodisperse suspensions decreases significantly with the increase of interparticle interactions. This has been justified by showing that the calculated values ofl* in a monodisperse suspension using Mie theory increase very significantly with the interactions incorporated inl* via the static structure factor.

Prediction of nonlocal scaling parameter for armchair and zigzag single-walled carbon nanotubes based on molecular structural mechanics, nonlocal elasticity and wave propagation

Atomic vibration in the Carbon Nanotubes (CNTs) gives rise to non-local interactions. In this paper, an expression for the non-local scaling parameter is derived as a function of the geometric and electronic properties of the rolled graphene sheet in single-walled CNTs. A self-consistent method is developed for the linearization of the problem of ultrasonic wave propagation in CNTs. We show that (i) the general three-dimensional elastic problem leads to a single non-local scaling parameter (e(0)), (ii) e(0) is almost constant irrespective of chirality of CNT in the case of longitudinal wave propagation, (iii) e(0) is a linear function of diameter of CNT for the case of torsional mode of wave propagation, (iv) e(0) in the case of coupled longitudinal-torsional modes of wave propagation, is a function which exponentially converges to that of axial mode at large diameters and to torsional mode at smaller diameters. These results are valid in the long-wavelength limit. (C) 2011 Elsevier Ltd. All rights reserved.

Estimation of surface temperature from MONTBLEX data

It is observed that the daily mean temperature of the soil is linear with depth and the variation of the temperature is sinusoidal with a period of a day. Based on these observations the one-dimensional heat conduction equation for the soil can be solved which gives the amplitude and phase variation of the temperature wave with depth. Given the temperature data at three levels below the surface, the amplitude and phase variation and hence the surface temperature variation over the day are estimated. The daily mean temperature of the surface is estimated from linear extrapolation of the daily means at the three levels below the surface. Estimated values of soil thermal diffusivity show a subtantial change after sudden and heavy rains.

Optical phase conjugation by two-wave mixing in iron-doped lithium niobate

A novel technique to generate forward phase conjugate wave by two-wave mixing (TWM) in photorefractive iron-doped lithium niobate crystal has been demonstrated. An optical beam from a positive transparency was forward phase conjugated by TWM technique. The experimental scheme was then extended to a specific interferometric application.

The role of Compton and Raman scattering in the quasar continuum

There are three ways in which an electromagnetic wave can undergo scattering in a plasma: (i) when the scattering of radiation occurs by a single electron, it is called Compton Scattering (CS); (ii) if it occurs by a longitudinal electron plasma mode, it is called Stimulated Raman Scattering (SRS), and (iii) if it occurs by a highly damped electron plasma mode, it is called Stimulated Compton Scattering (SCS). The non-thermal continuum of quasars is believed to be produced through the combined action of synchrotron and inverse Compton processes, which are essentially single-particle processes. Here, we investigate the role of SRS and SCS in the generation of continuum radiation from these compact objects. It is shown as an example that the complete spectrum of 3C 273 can be reproduced by suitably combining SCS and SRS. The differential contributions of SCS and SRS under different values of the plasma parameters are also calculated.

Ultrasonic wave characteristics of a monolayer graphene on silicon substrate

The present work deals with an ultrasonic type of wave propagation characteristics of monolayer graphene on silicon (Si) substrate. An atomistic model of a hybrid lattice involving a hexagonal lattice of graphene and surface atoms of diamond lattice of Si is developed to identify the carbon-silicon bond stiffness. Properties of this hybrid lattice model is then mapped into a nonlocal continuum framework. Equivalent force constant due to Si substrate is obtained by minimizing the total potential energy of the system. For this equilibrium configuration, the nonlocal governing equations are derived to analyze the ultrasonic wave dispersion based on spectral analysis. From the present analysis we show that the silicon substrate affects only the flexural wave mode. The frequency band gap of flexural mode is also significantly affected by this substrate. The results also show that, the silicon substrate adds cushioning effect to the graphene and it makes the graphene more stable. The analysis also show that the frequency bang gap relations of in-plane (longitudinal and lateral) and out-of-plane (flexural) wave modes depends not only on the y-direction wavenumber but also on nonlocal scaling parameter. In the nonlocal analysis, at higher values of the y-directional wavenumber, a decrease in the frequency band gap is observed for all the three fundamental wave modes in the graphene-silicon system. The atoms movement in the graphene due to the wave propagation are also captured for all the tree fundamental wave modes. The results presented in this work are qualitatively different from those obtained based on the local analysis and thus, are important for the development of graphene based nanodevices such as strain sensor, mass and pressure sensors, atomic dust detectors and enhancer of surface image resolution that make use of the ultrasonic wave dispersion properties of graphene. (C) 2011 Elsevier Ltd. All rights reserved.