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.

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.

Natural convection from a vertical cone in a porous medium due to the combined effects of heat and mass diffusion with non-uniform wall temperature/concentration or heat/mass flux and suction/injection

An analysis has been carried out to study the non-Darcy natural convention flow of Newtonian fluids on a vertical cone embedded in a saturated porous medium with power-law variation of the wall temperature/concentration or heat/mass flux and suction/injection with the streamwise distance x. Both non-similar and self-similar solutions have been obtained. The effects of non-Darcy parameter, ratio of the buoyancy forces due to mass and heat diffusion, variation of wall temperature/concentration or heat/mass flux and suction/injection on the Nusselt and Sherwood numbers have been studied.

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.

Moments of flux in intermediary liability for copyright infringement

This chapter provides an overview of a recent shift in regulatory strategies to address copyright infringement toward enlisting the assistance of general purpose Internet Service Providers. In Australia, the High Court held in 2012 that iiNet, a general purpose ISP, had no legal duty to police what its subscribers did with their internet connections. We provide an overview of three recent developments in Australian copyright law since that decision that demonstrate an emerging shift in the way that obligations are imposed on ISPs to govern the actions of their users without relying on secondary liability. The first is a new privately negotiated industry code that introduces a 'graduated response' system that requires ISPs to pass on warnings to subscribers who receive allegations of infringement. The second involves a recent series of Federal Court cases where rightsholders made a partially successful application to require ISPs to hand over the identifying details of subscribers whose households are alleged to have infringed copyright. The third is a new legislative scheme that will require ISPs to block access to foreign websites that 'facilitate' infringement. We argue that these shifts represent a greater sophistication in approaches to enrolling general purpose intermediaries in the regulatory project. We also suggest that these shifts represent a potentially disturbing trend towards enforcement of copyright law in a way that does not provide strong safeguards for the legitimate constitutional due process interests of users. We conclude with a call for greater attention and research to better understand how intermediaries make decisions when governing the conduct of users, how those decisions may be influenced by both state and non-state actors, and how the rights of individuals to due process can be adequately protected.

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.

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.

Forward osmosis as a pre-treatment for treating coal seam gas associated water: Flux and fouling behaviours

In this study, a bench scale forward osmosis (FO) process was operated using two commonly available FO membranes in different orientations in order to examine the removal of foulants in the coal seam gas (CSG) associated water, the water flux and fouling behaviours of the process were also investigated. After 48 h of fouling simulation experiment, the water flux declined by approximately 55 and 35% of its initial level in the TFC-PRO and CTA-PRO modes (support layer facing the feed), respectively, while the flux decline in the TFC-FO and CTA-FO modes (active layer facing the feed) was insignificant. The flux decline in PRO modes was caused by the compounding effects of internal concentration polarisation and membrane fouling. However, the declined flux was completely recovered to its initial level following the hydraulic cleaning using deionised water. Dissolved organic carbon (DOC), adenosine tri-phosphate (ATP) and major inorganic scalants (Ca, Mg and silica) in the CSG feed were effectively removed by using the FO process. The results of this study suggest that the FO process shows promising potential to be employed as an effective pre-treatment for membrane purification of CSG associated water.

Urban energy fluxes in Helsinki and their high frequency spectral corrections

Inadvertent climate modification has led to an increase in urban temperatures compared to the surrounding rural area. The main reason for the temperature rise is the altered energy portioning of input net radiation to heat storage and sensible and latent heat fluxes in addition to the anthropogenic heat flux. The heat storage flux and anthropogenic heat flux have not yet been determined for Helsinki and they are not directly measurable. To the contrary, turbulent fluxes of sensible and latent heat in addition to net radiation can be measured, and the anthropogenic heat flux together with the heat storage flux can be solved as a residual. As a result, all inaccuracies in the determination of the energy balance components propagate to the residual term and special attention must be paid to the accurate determination of the components. One cause of error in the turbulent fluxes is the fluctuation attenuation at high frequencies which can be accounted for by high frequency spectral corrections. The aim of this study is twofold: to assess the relevance of high frequency corrections to water vapor fluxes and to assess the temporal variation of the energy fluxes. Turbulent fluxes of sensible and latent heat have been measured at SMEAR III station, Helsinki, since December 2005 using the eddy covariance technique. In addition, net radiation measurements have been ongoing since July 2007. The used calculation methods in this study consist of widely accepted eddy covariance data post processing methods in addition to Fourier and wavelet analysis. The high frequency spectral correction using the traditional transfer function method is highly dependent on relative humidity and has an 11% effect on the latent heat flux. This method is based on an assumption of spectral similarity which is shown not to be valid. A new correction method using wavelet analysis is thus initialized and it seems to account for the high frequency variation deficit. Anyhow, the resulting wavelet correction remains minimal in contrast to the traditional transfer function correction. The energy fluxes exhibit a behavior characteristic for urban environments: the energy input is channeled to sensible heat as latent heat flux is restricted by water availability. The monthly mean residual of the energy balance ranges from 30 Wm-2 in summer to -35 Wm-2 in winter meaning a heat storage to the ground during summer. Furthermore, the anthropogenic heat flux is approximated to be 50 Wm-2 during winter when residential heating is important.

The sensitivity of the northern hemisphere summer mean meridional circulation to changes in the large scale eddy forcing

A zonally averaged version of the Goddard Laboratory for Atmospheric Sciences (GLAS) climate model is used to study the sensitivity of the northern hemisphere (NH) summer mean meridional circulation to changes in the large scale eddy forcing. A standard solution is obtained by prescribing the latent heating field and climatological horizontal transports of heat and momentum by the eddies. The radiative heating and surface fluxes are calculated by model parameterizations. This standard solution is compared with the results of several sensitivity studies. When the eddy forcing is reduced to 0.5 times or increased to 1.5 times the climatological values, the strength of the Ferrel cells decrease or increase proportionally. It is also seen that such changes in the eddy forcing can influence the strength of theNH Hadley cell significantly. Possible impact of such changes in the large scale eddy forcing on the monsoon circulation via changes in the Hadley circulation is discussed. Sensitivity experiments including only one component of eddy forcing at a time show that the eddy momentum fluxes seem to be more important in maintaining the Ferrel cells than the eddy heat fluxes. In the absence of the eddy heat fluxes, the observed eddy momentum fluxes alone produce subtropical westerly jets which are weaker than those in the standard solution. On the other hand, the observed eddy heat fluxes alone produce subtropical westerly jets which are stronger than those in the standard solution.

Algorithms for 13C metabolic flux analysis

The metabolism of an organism consists of a network of biochemical reactions that transform small molecules, or metabolites, into others in order to produce energy and building blocks for essential macromolecules. The goal of metabolic flux analysis is to uncover the rates, or the fluxes, of those biochemical reactions. In a steady state, the sum of the fluxes that produce an internal metabolite is equal to the sum of the fluxes that consume the same molecule. Thus the steady state imposes linear balance constraints to the fluxes. In general, the balance constraints imposed by the steady state are not sufficient to uncover all the fluxes of a metabolic network. The fluxes through cycles and alternative pathways between the same source and target metabolites remain unknown. More information about the fluxes can be obtained from isotopic labelling experiments, where a cell population is fed with labelled nutrients, such as glucose that contains 13C atoms. Labels are then transferred by biochemical reactions to other metabolites. The relative abundances of different labelling patterns in internal metabolites depend on the fluxes of pathways producing them. Thus, the relative abundances of different labelling patterns contain information about the fluxes that cannot be uncovered from the balance constraints derived from the steady state. The field of research that estimates the fluxes utilizing the measured constraints to the relative abundances of different labelling patterns induced by 13C labelled nutrients is called 13C metabolic flux analysis. There exist two approaches of 13C metabolic flux analysis. In the optimization approach, a non-linear optimization task, where candidate fluxes are iteratively generated until they fit to the measured abundances of different labelling patterns, is constructed. In the direct approach, linear balance constraints given by the steady state are augmented with linear constraints derived from the abundances of different labelling patterns of metabolites. Thus, mathematically involved non-linear optimization methods that can get stuck to the local optima can be avoided. On the other hand, the direct approach may require more measurement data than the optimization approach to obtain the same flux information. Furthermore, the optimization framework can easily be applied regardless of the labelling measurement technology and with all network topologies. In this thesis we present a formal computational framework for direct 13C metabolic flux analysis. The aim of our study is to construct as many linear constraints to the fluxes from the 13C labelling measurements using only computational methods that avoid non-linear techniques and are independent from the type of measurement data, the labelling of external nutrients and the topology of the metabolic network. The presented framework is the first representative of the direct approach for 13C metabolic flux analysis that is free from restricting assumptions made about these parameters.In our framework, measurement data is first propagated from the measured metabolites to other metabolites. The propagation is facilitated by the flow analysis of metabolite fragments in the network. Then new linear constraints to the fluxes are derived from the propagated data by applying the techniques of linear algebra.Based on the results of the fragment flow analysis, we also present an experiment planning method that selects sets of metabolites whose relative abundances of different labelling patterns are most useful for 13C metabolic flux analysis. Furthermore, we give computational tools to process raw 13C labelling data produced by tandem mass spectrometry to a form suitable for 13C metabolic flux analysis.

Heat transfer to power-law fluids in thermal entrance region with viscous dissipation for constant heat flux conditions

An analytical solution of the heat transfer problem with viscous dissipation for non-Newtonian fluids with power-law model in the thermal entrance region of a circular pipe and two parallel plates under constant heat flux conditions is obtained using eigenvalue approach by suitably replacing one of the boundary conditions by total energy balance equation. Analytical expressions for the wall and the bulk temperatures and the local Nusselt number are presented. The results are in close agreement with those obtained by implicit finite-difference scheme. It is found that the role of viscous dissipation on heat transfer is completely different for heating and cooling conditions at the wall. The results for the case of cooling at the wall are of interest in the design of the oil pipe line.

A possible explanation of the Maunder minimum from a flux transport dynamo model

We propose that the poloidal field at the end of the last sunspot cycle before the Maunder minimum fell to a very low value due to fluctuations in the Babcock-Leighton process. With this assumption, a flux transport dynamo model is able to explain various aspects of the historical records of the Maunder minimum remarkably well by suitably choosing the parameters of the model to give the correct growth time.

Response of a turbulent boundary layer to a step change in surface heat flux

Measurements of both the velocity and the temperature field have been made in the thermal layer that grows inside a turbulent boundary layer which is subjected to a small step change in surface heat flux. Upstream of the step, the wall heat flux is zero and the velocity boundary layer is nearly self-preserving. The thermal-layer measurements are discussed in the context of a self-preserving analysis for the temperature disturbance which grows underneath a thick external turbulent boundary layer. A logarithmic mean temperature profile is established downstream of the step but the budget for the mean-square temperature fluctuations shows that, in the inner region of the thermal layer, the production and dissipation of temperature fluctuations are not quite equal at the furthest downstream measurement station. The measurements for both the mean and the fluctuating temperature field indicate that the relaxation distance for the thermal layer is quite large, of the order of 1000θ0, where θ0 is the momentum thickness of the boundary layer at the step. Statistics of the thermal-layer interface and conditionally sampled measurements with respect to this interface are presented. Measurements of the temperature intermittency factor indicate that the interface is normally distributed with respect to its mean position. Near the step, the passive heat contaminant acts as an effective marker of the organized turbulence structure that has been observed in the wall region of a boundary layer. Accordingly, conditional averages of Reynolds stresses and heat fluxes measured in the heated part of the flow are considerably larger than the conventional averages when the temperature intermittency factor is small.

Carbon balance and component CO2 fluxes in boreal Scots pine stands

The forest vegetation takes up atmospheric carbon dioxide (CO2) in photosynthesis. Part of the fixed carbon is released back into the atmosphere during plant respiration but a substantial part is stored as plant biomass, especially in the stems of trees. Carbon also accumulates in the soil as litter and via the roots. CO2 is released into the atmosphere from these carbon stocks in the decomposition of dead biomass. Carbon balance of a forest stand is the difference between the CO2 uptake and CO2 efflux. This study quantifies and analyses the dynamics of carbon balance and component CO2 fluxes in four Southern Finnish Scots pine stands that covered the typical economic rotation time of 80 years. The study was based on direct flux measurements with chambers and eddy covariance (EC), and modelling of component CO2 fluxes. The net CO2 exchange of the stand was partitioned into component fluxes: photosynthesis of trees and ground vegetation, respiration of tree foliage and stems, and CO2 efflux from the soil. The relationships between the component fluxes and the environmental factors (light, temperature, atmospheric CO2, air humidity and soil moisture) were studied with mathematical modelling. The annual CO2 balance varied from a source of about 400 g C/m2 at a recently clearcut site to net CO2 uptake of 200 300 g C/m2 in a middle-aged (40-year-old) and a mature (75-year-old) stand. A 12-year-old sapling site was at the turning point from source to a sink of CO2. In the middle-aged stand, photosynthetic production was dominated by trees. Under closed pine canopies, ground vegetation accounted for 10 20% of stand photosynthesis whereas at the open sites the proportion and also the absolute photosynthesis of ground vegetation was much higher. The aboveground respiration was dominated by tree foliage which accounted for one third of the ecosystem respiration. Rate of wood respiration was in the order of 10% of total ecosystem respiration. CO2 efflux from the soil dominated the ecosystem respiratory fluxes in all phases of stand development. Instantaneous and delayed responses to the environmental driving factors could predict well within-year variability in photosynthetic production: In the short term and during the growing season photosynthesis follows primarily light while the seasonal variation is more strongly connected to temperature. The temperature relationship of the annual cycle of photosynthesis was found to be almost equal in the southern boreal zone and at the timberline in the northern boreal zone. The respiratory fluxes showed instantaneous and seasonal temperature relationships but they could also be connected to photosynthesis at an annual timescale.