Finite volume evolution Galerkin method for hyperbolic conservation laws with spatially varying flux functions

We present a generalization of the finite volume evolution Galerkin scheme [M. Lukacova-Medvid'ova,J. Saibertov'a, G. Warnecke, Finite volume evolution Galerkin methods for nonlinear hyperbolic systems, J. Comp. Phys. (2002) 183 533-562; M. Luacova-Medvid'ova, K.W. Morton, G. Warnecke, Finite volume evolution Galerkin (FVEG) methods for hyperbolic problems, SIAM J. Sci. Comput. (2004) 26 1-30] for hyperbolic systems with spatially varying flux functions. Our goal is to develop a genuinely multi-dimensional numerical scheme for wave propagation problems in a heterogeneous media. We illustrate our methodology for acoustic waves in a heterogeneous medium but the results can be generalized to more complex systems. The finite volume evolution Galerkin (FVEG) method is a predictor-corrector method combining the finite volume corrector step with the evolutionary predictor step. In order to evolve fluxes along the cell interfaces we use multi-dimensional approximate evolution operator. The latter is constructed using the theory of bicharacteristics under the assumption of spatially dependent wave speeds. To approximate heterogeneous medium a staggered grid approach is used. Several numerical experiments for wave propagation with continuous as well as discontinuous wave speeds confirm the robustness and reliability of the new FVEG scheme.

Mixed convection in non-Newtonian fluids along a vertical plate in a porous medium

The problem of mixed convection from vertical surfaces in a porous medium saturated with a power-law type non-Newtonian fluid is investigated. The transformed conservation laws are solved numerically for the case of variable wall temperature conditions. Results for the details of the velocity and temperature fields as well as the Nusselt number have been presented. The viscosity index ranged from 0.5-2.0.

Vertical structure and horizontal gradients of aerosol extinction coefficients over coastal India inferred from airborne lidar measurements during the Integrated Campaign for Aerosol, Gases and Radiation Budget (ICARB) field campaign

Quantitative estimates of the vertical structure and the spatial gradients of aerosol extinction coefficients have been made from airborne lidar measurements across the coastline into offshore oceanic regions along the east and west coasts of India. The vertical structure revealed the presence of strong, elevated aerosol layers in the altitude region of similar to 2-4 km, well above the atmospheric boundary layer (ABL). Horizontal gradients also showed a vertical structure, being sharp with the e(-1) scaling distance (D-0H) as small as similar to 150 km in the well-mixed regions mostly under the influence of local source effects. Above the ABL, where local effects are subdued, the gradients were much shallower (similar to 600-800 km); nevertheless, they were steep compared to the value of similar to 1500-2500 km reported for columnar AOD during winter. The gradients of these elevated layers were steeper over the east coast of India than over the west coast. Near-simultaneous radio sonde (Vaisala, Inc., Finland) ascents made over the northern Bay of Bengal showed the presence of convectively unstable regions, first from surface to similar to 750-1000 m and the other extending from 1750 to 3000 m separated by a stable region in between. These can act as a conduit for the advection of aerosols and favor the transport of continental aerosols in the higher levels (> 2 km) into the oceans without entering the marine boundary layer below. Large spatial gradient in aerosol optical and hence radiative impacts between the coastal landmass and the adjacent oceans within a short distance of < 300 km (even at an altitude of 3 km) during summer and the premonsoon is of significance to the regional climate.

Using lead sulphate instead of lead oxide in flux growth of yttrium iron garnet(YIG) crystals

Crystal growth of YIG from fluxes containing lead sulphate in place of lead oxide in the usual lead oxide-lead fluoride-boron oxide flux system has been tried. Lead sulphate decomposes during crystal growth giving lead oxide and sulphur trioxide. Due to the influence of sulphur trioxide in the system the yield of crystals almost doubles. There is no change either in the morphology of the crystals or their lattice parameter. It is possible that solubility of YIG is different in the new flux and the changed solubility causes the increase in yield of crystals.

Flux balance analysis of biological systems: applications and challenges

Systems level modelling and simulations of biological processes are proving to be invaluable in obtaining a quantitative and dynamic perspective of various aspects of cellular function. In particular, constraint-based analyses of metabolic networks have gained considerable popularity for simulating cellular metabolism, of which flux balance analysis (FBA), is most widely used. Unlike mechanistic simulations that depend on accurate kinetic data, which are scarcely available, FBA is based on the principle of conservation of mass in a network, which utilizes the stoichiometric matrix and a biologically relevant objective function to identify optimal reaction flux distributions. FBA has been used to analyse genome-scale reconstructions of several organisms; it has also been used to analyse the effect of perturbations, such as gene deletions or drug inhibitions in silico. This article reviews the usefulness of FBA as a tool for gaining biological insights, advances in methodology enabling integration of regulatory information and thermodynamic constraints, and finally addresses the challenges that lie ahead. Various use scenarios and biological insights obtained from FBA, and applications in fields such metabolic engineering and drug target identification, are also discussed. Genome-scale constraint-based models have an immense potential for building and testing hypotheses, as well as to guide experimentation.

CFD Analysis of conjugate natural convection in a vertical annulus

Conjugate natural convection in a vertical annulus with a centrally located vertical heat generating rod is studied numerically. The governing equations are discretized on a staggered mesh and are solved using a pressure-correction algorithm. A parametric study is performed by varying the Grashof number, aspect ratio, and the solid-to-fluid thermal conductivity ratio over wide ranges with the Prandtl number fixed at 0.7. Results are presented for the variation of several quantities of interest such as the local Nusselt numbers on the inner and outer boundaries, the axial variation of the centerline and interface temperatures, maximum solid, average solid and average interface temperature variations with Grashof number, and the average Nusselt number variation for the inner and outer boundaries with Grashof number. The average Nusselt number from the conjugate analysis is found to be between the Nusselt numbers of the isothermal and the isoflux cases. The average Nusselt numbers on the inner and outer boundaries show an increasing trend with the Grashof number. Correlations are presented for the Nusselt number and the dimensionless temperatures of interest in terms of the parameters of the problem.

Vertical uplift capacity of two interfering horizontal anchors in sand using an upper bound limit analysis

The vertical uplift resistance of two interfering rigid rough strip anchors embedded horizontally in sand at shallow depths has been examined. The analysis is performed by using an upper bound theorem o limit analysis in combination with finite elements and linear programming. It is specified that both the anchors are loaded to failure simultaneously at the same magnitude of the failure load. For different clear spacing (S) between the anchors, the magnitude of the efficiency factor (xi(gamma)) is determined. On account of interference, the magnitude of xi(gamma) is found to reduce continuously with a decrease in the spacing between the anchors. The results from the numerical analysis were found to compare reasonably well with the available theoretical data from the literature.

On the Landau damping of sound waves in a plasma with magnetic flux tubes

The Landau damping of sound wave in a plasma consisting of an ensemble of magnetic flux tubes with reference to the work by Ryutov and Ryutova (1976) is discussed. Sound waves cannot be Landau damped in general but under certain restriction conditions on plasma parameters the possibility of absorption of these waves can exist.

Axial heat conduction effects in natural convection along a vertical cylinder

MANY TRANSPORprTo cesses occur in nature and in industrial applications in which the transfer of heat is governed by the process of natural convection. Natural convection arises in fluids when the temperature changes cause density variations leading to buoyancy forces. An excellent review of natural convection flows has been given by Ede [I]. Recently, Minkowycz and Sparrow [2, 31, Cebeci [4], and Aziz and Na [S] have studied the steady, laminar, incompressible, natural convection flow over a vertical cylinder using a local nonsimilarity method, a finite-difference scheme, and an improved perturbation method, respectively. However, they did not take into account the effect ofaxial heat conduction for small Prandtl number. It is known that the axial heat conductioneffect becomesimportant for low-Prandtl-number fluids such as a liquid metal.

Axisymmetric melting of a long cylinder due to an infinite flux

By employing a new embedding technique, a short-time analytical solution for the axisymmetric melting of a long cylinder due to an infinite flux is presented in this paper. The sufficient condition for starting the instantaneous melting of the cylinder has been derived. The melt is removed as soon as it is formed. The method of solution is simple and straightforward and consists of assuming fictitious initial temperature for some fictitious extension of the actual region.

Unsteady nonsimilar natural convection over a vertical flat plate in a thermally stratified fluid

Abstract is not available.

Free convection heat transfer to mercury in vertical annuli

Data on free convection heat transfer to water and mercury are collected using a test rig in vertical annuli of three radii ratios, the walls of which are maintained at uniform temperatures. A theoretical analysis of the boundary layer equations has been attempted using local similarity transformation and double boundary layer approach. Correlations derived from the present theoretical analysis are compared with the analysis and the experimental data available in literature for non-metallic fluids and also with the present experimental data on water and mercury. Generalised correlations are set up for expressing the ratio of heat transferred by convection to the heat transferred by pure conduction and Nusselt's number, in terms of Grashof, Rayleigh and Prandtl numbers, based on the theoretical analysis and the present data on mercury and water. The present generalised correlations agree with the reported and present data for non-metallic fluids and liquid metals with an average deviation of 9% and maximum deviation of ± 13.7%.

Water-wave scattering by two submerged plane vertical barriers-Abel integral-equation approach

The classical problem of surface water-wave scattering by two identical thin vertical barriers submerged in deep water and extending infinitely downwards from the same depth below the mean free surface, is reinvestigated here by an approach leading to the problem of solving a system of Abel integral equations. The reflection and transmission coefficients are obtained in terms of computable integrals. Known results for a single barrier are recovered as a limiting case as the separation distance between the two barriers tends to zero. The coefficients are depicted graphically in a number of figures which are identical with the corresponding figures given by Jarvis (J Inst Math Appl 7:207-215, 1971) who employed a completely different approach involving a Schwarz-Christoffel transformation of complex-variable theory to solve the problem.