Temperature modelling of an homogeneous medium using genetically selected RBF (LIC)

Temperature modelling of human tissue subjected to ultrasound for therapeutic use is essencial for an accurate instrumental assessment and calibration. In this paper punctual temperature modeling of a homogeneous medium, radiated by therapeutic ultrasound, is presented.

Front spreading in population dynamic models. Theory and application to the Neolithic transition

This thesis presents population dynamics models that can be applied to predict the rate of spread of the Neolithic transition (change from hunter-gathering to farming economics) across the European continent, which took place about 9000 to 5000 years ago. The first models in this thesis provide predictions at a continental scale. We develop population dynamics models with explicit kernels and apply realistic data. We also derive a new time-delayed reaction-diffusion equation which yields speeds about a 10% slower than previous models. We also deal with a regional variability: the slowdown of the Neolithic front when reaching the North of Europe. We develop simple reaction-diffusion models that can predict the measured speeds in terms of the non-homogeneous distribution of pre-Neolithic (Mesolithic) population in Europe, which were present in higher densities at the North of the continent. Such models can explain the observed speeds.

Homogeneous and heterogeneous reactions of atmospheric mercury(II) with sulfur(IV)

Atmospheric models suggest that the reduction of Hg(II) to Hg(O) by S(W) prolongs the residence time of mercury. The redox reaction was investigated both in the aqueous phase (where the reductant is sulfite) and on particulate matter (where the reductant in SO2(g)). In both cases, one of the ultimate products is HgS. A mechanism is proposed involving formation of Hg(O) followed by mercury-induced disproportionation of SO2.

Wave-activity conservation laws for the three-dimensional anelastic and Boussinesq equations with a horizontally homogeneous background flow

Wave-activity conservation laws are key to understanding wave propagation in inhomogeneous environments. Their most general formulation follows from the Hamiltonian structure of geophysical fluid dynamics. For large-scale atmospheric dynamics, the Eliassen–Palm wave activity is a well-known example and is central to theoretical analysis. On the mesoscale, while such conservation laws have been worked out in two dimensions, their application to a horizontally homogeneous background flow in three dimensions fails because of a degeneracy created by the absence of a background potential vorticity gradient. Earlier three-dimensional results based on linear WKB theory considered only Doppler-shifted gravity waves, not waves in a stratified shear flow. Consideration of a background flow depending only on altitude is motivated by the parameterization of subgrid-scales in climate models where there is an imposed separation of horizontal length and time scales, but vertical coupling within each column. Here we show how this degeneracy can be overcome and wave-activity conservation laws derived for three-dimensional disturbances to a horizontally homogeneous background flow. Explicit expressions for pseudoenergy and pseudomomentum in the anelastic and Boussinesq models are derived, and it is shown how the previously derived relations for the two-dimensional problem can be treated as a limiting case of the three-dimensional problem. The results also generalize earlier three-dimensional results in that there is no slowly varying WKB-type requirement on the background flow, and the results are extendable to finite amplitude. The relationship A E =cA P between pseudoenergy A E and pseudomomentum A P, where c is the horizontal phase speed in the direction of symmetry associated with A P, has important applications to gravity-wave parameterization and provides a generalized statement of the first Eliassen–Palm theorem.

Reducing noise associated with the Monte Carlo Independent Column Approximation for weather forecasting models

The Monte Carlo Independent Column Approximation (McICA) is a flexible method for representing subgrid-scale cloud inhomogeneity in radiative transfer schemes. It does, however, introduce conditional random errors but these have been shown to have little effect on climate simulations, where spatial and temporal scales of interest are large enough for effects of noise to be averaged out. This article considers the effect of McICA noise on a numerical weather prediction (NWP) model, where the time and spatial scales of interest are much closer to those at which the errors manifest themselves; this, as we show, means that noise is more significant. We suggest methods for efficiently reducing the magnitude of McICA noise and test these methods in a global NWP version of the UK Met Office Unified Model (MetUM). The resultant errors are put into context by comparison with errors due to the widely used assumption of maximum-random-overlap of plane-parallel homogeneous cloud. For a simple implementation of the McICA scheme, forecasts of near-surface temperature are found to be worse than those obtained using the plane-parallel, maximum-random-overlap representation of clouds. However, by applying the methods suggested in this article, we can reduce noise enough to give forecasts of near-surface temperature that are an improvement on the plane-parallel maximum-random-overlap forecasts. We conclude that the McICA scheme can be used to improve the representation of clouds in NWP models, with the provision that the associated noise is sufficiently small.

THE CO-INFLUENCE OF THE SEA BREEZE AND THE COASTAL UPWELLING AT CABO FRIO: A NUMERICAL INVESTIGATION USING COUPLED MODELS

A coupled atmospheric-oceanic model was used to investigate whether there is a positive feedback between the coastal upwelling and the sea breeze at Cabo Frio - RJ (Brazil). Two experiments were performed to ascertain the influence of the sea breeze on the coastal upwelling: the first one used the coupled model forced with synoptic NE winds of 8 m s(-1) and the sign of the sea breeze circulation was set by the atmospheric model; the second experiment used only the oceanic model with constant 8 m s(-1) NE winds. Then, to study the influence of the coastal upwelling on the sea breeze, two more experiments were performed: one using a coastal upwelling representative SST initial field and the other one using a constant and homogeneous SST field of 26 degrees C. Finally, two more experiments were conducted to verify the influence of the topography and the spatial distribution of the sea surface temperature on the previous results. The results showed that the sea breeze can intensify the coastal upwelling, but the coastal upwelling does not intensify the sea breeze circulation, suggesting that there is no positive feedback between these two phenomena at Cabo Frio.

GHASP: an H alpha kinematic survey of spiral and irregular galaxies - IX. The near-infrared, stellar and baryonic Tully-Fisher relations

We studied, for the first time, the near-infrared, stellar and baryonic Tully-Fisher relations for a sample of field galaxies taken from a homogeneous Fabry-Perot sample of galaxies [the Gassendi HAlpha survey of SPirals (GHASP) survey]. The main advantage of GHASP over other samples is that the maximum rotational velocities were estimated from 2D velocity fields, avoiding assumptions about the inclination and position angle of the galaxies. By combining these data with 2MASS photometry, optical colours, HI masses and different mass-to-light ratio estimators, we found a slope of 4.48 +/- 0.38 and 3.64 +/- 0.28 for the stellar and baryonic Tully-Fisher relation, respectively. We found that these values do not change significantly when different mass-to-light ratio recipes were used. We also point out, for the first time, that the rising rotation curves as well as asymmetric rotation curves show a larger dispersion in the Tully-Fisher relation than the flat ones or the symmetric ones. Using the baryonic mass and the optical radius of galaxies, we found that the surface baryonic mass density is almost constant for all the galaxies of this sample. In this study we also emphasize the presence of a break in the NIR Tully-Fisher relation at M(H,K) similar to -20 and we confirm that late-type galaxies present higher total-to-baryonic mass ratios than early-type spirals, suggesting that supernova feedback is actually an important issue in late-type spirals. Due to the well-defined sample selection criteria and the homogeneity of the data analysis, the Tully-Fisher relation for GHASP galaxies can be used as a reference for the study of this relation in other environments and at higher redshifts.

Bayesian analysis and constraints on kinematic models from union SNIa

The kinematic expansion history of the universe is investigated by using the 307 supernovae type Ia from the Union Compilation set. Three simple model parameterizations for the deceleration parameter ( constant, linear and abrupt transition) and two different models that are explicitly parametrized by the cosmic jerk parameter ( constant and variable) are considered. Likelihood and Bayesian analyses are employed to find best fit parameters and compare models among themselves and with the flat Lambda CDM model. Analytical expressions and estimates for the deceleration and cosmic jerk parameters today (q(0) and j(0)) and for the transition redshift (z(t)) between a past phase of cosmic deceleration to a current phase of acceleration are given. All models characterize an accelerated expansion for the universe today and largely indicate that it was decelerating in the past, having a transition redshift around 0.5. The cosmic jerk is not strongly constrained by the present supernovae data. For the most realistic kinematic models the 1 sigma confidence limits imply the following ranges of values: q(0) is an element of [-0.96, -0.46], j(0) is an element of [-3.2,-0.3] and z(t) is an element of [0.36, 0.84], which are compatible with the Lambda CDM predictions, q(0) = -0.57 +/- 0.04, j(0) = -1 and z(t) = 0.71 +/- 0.08. We find that even very simple kinematic models are equally good to describe the data compared to the concordance Lambda CDM model, and that the current observations are not powerful enough to discriminate among all of them.

Estimating the number of ozone peaks in Mexico City using a non-homogeneous Poisson model

In this paper, we consider the problem of estimating the number of times an air quality standard is exceeded in a given period of time. A non-homogeneous Poisson model is proposed to analyse this issue. The rate at which the Poisson events occur is given by a rate function lambda(t), t >= 0. This rate function also depends on some parameters that need to be estimated. Two forms of lambda(t), t >= 0 are considered. One of them is of the Weibull form and the other is of the exponentiated-Weibull form. The parameters estimation is made using a Bayesian formulation based on the Gibbs sampling algorithm. The assignation of the prior distributions for the parameters is made in two stages. In the first stage, non-informative prior distributions are considered. Using the information provided by the first stage, more informative prior distributions are used in the second one. The theoretical development is applied to data provided by the monitoring network of Mexico City. The rate function that best fit the data varies according to the region of the city and/or threshold that is considered. In some cases the best fit is the Weibull form and in other cases the best option is the exponentiated-Weibull. Copyright (C) 2007 John Wiley & Sons, Ltd.

On the Casimir energy for a massive quantum scalar field and the cosmological constant

We present a rigorous, regularization-independent local quantum field theoretic treatment of the Casimir effect for a quantum scalar field of mass mu not equal 0 which yields closed form expressions for the energy density and pressure. As an application we show that there exist special states of the quantum field in which the expectation value of the renormalized energy-momentum tensor is, for any fixed time, independent of the space coordinate and of the perfect fluid form g(mu,nu)rho with rho > 0, thus providing a concrete quantum field theoretic model of the cosmological constant. This rho represents the energy density associated to a state consisting of the vacuum and a certain number of excitations of zero momentum, i.e., the constituents correspond to lowest energy and pressure p <= 0. (C) 2009 Elsevier Inc. All rights reserved.

Stability of the curvature perturbation in dark sectors` mutual interacting models

We consider perturbations in a cosmological model with a small coupling between dark energy and dark matter. We prove that the stability of the curvature perturbation depends on the type of coupling between dark sectors. When the dark energy is of quintessence type, if the coupling is proportional to the dark matter energy density, it will drive the instability in the curvature perturbations: however if the coupling is proportional to the energy density of dark energy, there is room for the stability in the curvature perturbations. When the dark energy is of phantom type, the perturbations are always stable, no matter whether the coupling is proportional to the one or the other energy density. (C) 2008 Elsevier B.V. All rights reserved.

Physically based simulation of heterogeneous deformable models using XFEM

This paper addresses the problem of heterogeneous deformable model accuracy using the finite element methods (FEM). Classic FEM uses predefined shape functions for interpolation and does not account easily for regions of discontinuities. Extended finite element methods (XFEM) use enrichment functions to compensate for the change in an element degrees of freedom (DoFs) in deformable objects. The XFEM is an accurate and fast method as no remeshing is required. In this study we investigate the performance of XFEM and demonstrate how it may be applied to discontinuities of materials that exist in heterogeneous (piece-wise homogeneous) models. The results show realistic stress prediction compared to modeling the same objects with classic FEM.

Non-homogeneous volatility correlations in the bivariate multifractal model

In this paper, we consider an extension of the recently proposed bivariate Markov-switching multifractal model of Calvet, Fisher, and Thompson [2006. "Volatility Comovement: A Multifrequency Approach." Journal of Econometrics {131}: 179-215]. In particular, we allow correlations between volatility components to be non-homogeneous with two different parameters governing the volatility correlations at high and low frequencies. Specification tests confirm the added explanatory value of this specification. In order to explore its practical performance, we apply the model for computing value-at-risk statistics for different classes of financial assets and compare the results with the baseline, homogeneous bivariate multifractal model and the bivariate DCC-GARCH of Engle [2002. "Dynamic Conditional Correlation: A Simple Class of Multivariate Generalized Autoregressive Conditional Heteroskedasticity Models." Journal of Business & Economic Statistics 20 (3): 339-350]. As it turns out, the multifractal model with heterogeneous volatility correlations provides more reliable results than both the homogeneous benchmark and the DCC-GARCH model. © 2014 Taylor & Francis.

Efeitos das inomogeneidades da matéria em Cosmologias Aceleradas

The recent observational advances of Astronomy and a more consistent theoretical framework turned Cosmology in one of the most exciting frontiers of contemporary science. In this thesis, homogeneous and inhomogeneous Universe models containing dark matter and different kinds of dark energy are confronted with recent observational data. Initially, we analyze constraints from the existence of old high redshift objects, Supernovas type Ia and the gas mass fraction of galaxy clusters for 2 distinct classes of homogeneous and isotropic models: decaying vacuum and X(z)CDM cosmologies. By considering the quasar APM 08279+5255 at z = 3.91 with age between 2-3 Gyr, we obtain 0,2 < OM < 0,4 while for the j3 parameter which quantifies the contribution of A( t) is restricted to the intervalO, 07 < j3 < 0,32 thereby implying that the minimal age of the Universe amounts to 13.4 Gyr. A lower limit to the quasar formation redshift (zJ > 5,11) was also obtained. Our analyzes including flat, closed and hyperbolic models show that there is no an age crisis for this kind of decaying A( t) scenario. Tests from SN e Ia and gas mass fraction data were realized for flat X(z)CDM models. For an equation of state, úJ(z) = úJo + úJIZ, the best fit is úJo = -1,25, úJl = 1,3 and OM = 0,26, whereas for models with úJ(z) = úJo+úJlz/(l+z), we obtainúJo = -1,4, úJl = 2,57 and OM = 0,26. In another line of development, we have discussed the influence of the observed inhomogeneities by considering the Zeldovich-Kantowski-DyerRoeder (ZKDR) angular diameter distance. By applying the statistical X2 method to a sample of angular diameter for compact radio sources, the best fit to the cosmological parameters for XCDM models are OM = O, 26,úJ = -1,03 and a = 0,9, where úJ and a are the equation of state and the smoothness parameters, respectively. Such results are compatible with a phantom energy component (úJ < -1). The possible bidimensional spaces associated to the plane (a , OM) were restricted by using data from SNe Ia and gas mass fraction of galaxy clusters. For Supernovas the parameters are restricted to the interval 0,32 < OM < 0,5(20") and 0,32 < a < 1,0(20"), while to the gas mass fraction we find 0,18 < OM < 0,32(20") with alI alIowed values of a. For a joint analysis involving Supernovas and gas mass fraction data we obtained 0,18 < OM < 0,38(20"). In general grounds, the present study suggests that the influence of the cosmological inhomogeneities in the matter distribution need to be considered with more detail in the analyses of the observational tests. Further, the analytical treatment based on the ZKDR distance may give non-negligible corrections to the so-calIed background tests of FRW type cosmologies