Numerical methods and tangible interfaces for pollutant dispersion simulation

Dissertação apresentada para obtenção do Grau de Doutor em Engenharia do Ambiente, pela Universidade Nova de Lisboa, Faculdade de Ciências e Tecnologia

Two different 3D biogeochemical models for the NW mediterranean sea: validation with Meris data and intercomparison

Report for the scientific sojourn carried out at the University of New South Wales from February to June the 2007. Two different biogeochemical models are coupled to a three dimensional configuration of the Princeton Ocean Model (POM) for the Northwestern Mediterranean Sea (Ahumada and Cruzado, 2007). The first biogeochemical model (BLANES) is the three-dimensional version of the model described by Bahamon and Cruzado (2003) and computes the nitrogen fluxes through six compartments using semi-empirical descriptions of biological processes. The second biogeochemical model (BIOMEC) is the biomechanical NPZD model described in Baird et al. (2004), which uses a combination of physiological and physical descriptions to quantify the rates of planktonic interactions. Physical descriptions include, for example, the diffusion of nutrients to phytoplankton cells and the encounter rate of predators and prey. The link between physical and biogeochemical processes in both models is expressed by the advection-diffusion of the non-conservative tracers. The similarities in the mathematical formulation of the biogeochemical processes in the two models are exploited to determine the parameter set for the biomechanical model that best fits the parameter set used in the first model. Three years of integration have been carried out for each model to reach the so called perpetual year run for biogeochemical conditions. Outputs from both models are averaged monthly and then compared to remote sensing images obtained from sensor MERIS for chlorophyll.

Analytical expressions for the flame front speed in the downward combustion of thin solid fuels and comparison to experiments

We derive analytical expressions for the propagation speed of downward combustion fronts of thin solid fuels with a background flow initially at rest. The classical combustion model for thin solid fuels that consists of five coupled reaction-convection-diffusion equations is here reduced into a single equation with the gas temperature as the single variable. For doing so we apply a two-zone combustion model that divides the system into a preheating region and a pyrolyzing region. The speed of the combustion front is obtained after matching the temperature and its derivative at the location that separates both regions.We also derive a simplified version of this analytical expression expected to be valid for a wide range of cases. Flame front velocities predicted by our analyticalexpressions agree well with experimental data found in the literature for a large variety of cases and substantially improve the results obtained from a previous well-known analytical expression

Anisotropic dispersion, space competition, and the slowdown of the Neolithic transition

The front speed of the Neolithic (farmer) spread in Europe decreased as it reached Northern latitudes, where the Mesolithic (huntergatherer) population density was higher. Here, we describe a reaction diffusion model with (i) an anisotropic dispersion kernel depending on the Mesolithicpopulation density gradient and (ii) a modified population growth equation. Both effects are related to the space available for the Neolithic population. The model is able to explain the slowdown of the Neolithic front as observed from archaeological data

Speed of reaction-transport processes

We present an approach to determining the speed of wave-front solutions to reaction-transport processes. This method is more accurate than previous ones. This is explicitly shown for several cases of practical interest: (i) the anomalous diffusion reaction, (ii) reaction diffusion in an advective field, and (iii) time-delayed reaction diffusion. There is good agreement with the results of numerical simulations

A methodology for the adaptation of a PMP at the determination of a PMF

This paper presents a new type of very fine grid hydrological model based on the spatiotemporal repartition of a PMP (Probable Maximum Precipitation) and on the topography. The goal is to estimate the influence of this rain on a PMF (Probable Maximum Flood) on a catchment area in Switzerland. The spatiotemporal distribution of the PMP was realized using six clouds modeled by the advection-diffusion equation. The equation shows the movement of the clouds over the terrain and also gives the evolution of the rain intensity in time. This hydrological modeling is followed by a hydraulic modeling of the surface and subterranean flow, done considering the factors that contribute to the hydrological cycle, such as the infiltration, the resurgence and the snowmelt. These added factors make the developed model closer to reality and also offer flexibility in the initial condition that is added to the factors concerning the PMP, such as the duration of the rain, the speed and direction of the wind. All these initial conditions taken together offer a complete image of the PMF.

Time-Delayed Theory of the Neolithic Transition in Europe

The classical wave-of-advance model of the neolithic transition (i.e., the shift from hunter-gatherer to agricultural economies) is based on Fisher's reaction-diffusion equation. Here we present an extension of Einstein's approach to Fickian diffusion, incorporating reaction terms. On this basis we show that second-order terms in the reaction-diffusion equation, which have been neglected up to now, are not in fact negligible but can lead to important corrections. The resulting time-delayed model agrees quite well with observations

Accelerated tumor invasion under non-isotropic cell dispersal in glioblastomas

Glioblastomas are highly diffuse, malignant tumors that have so far evaded clinical treatment. The strongly invasive behavior of cells in these tumors makes them very resistant to treatment, and for this reason both experimental and theoretical efforts have been directed toward understanding the spatiotemporal pattern of tumor spreading. Although usual models assume a standard diffusion behavior, recent experiments with cell cultures indicate that cells tend to move in directions close to that of glioblastoma invasion, thus indicating that a biasedrandom walk model may be much more appropriate. Here we show analytically that, for realistic parameter values, the speeds predicted by biased dispersal are consistent with experimentally measured data. We also find that models beyond reaction–diffusion–advection equations are necessary to capture this substantial effect of biased dispersal on glioblastoma spread

Dispersal probability distributions and the wave-front speed problem

The speed and width of front solutions to reaction-dispersal models are analyzed both analytically and numerically. We perform our analysis for Laplace and Gaussian distribution kernels, both for delayed and nondelayed models. The results are discussed in terms of the characteristic parameters of the models

Kinetics of K release from soils of Brazilian coffee regions: effect of organic acids

Kinetic studies on soil potassium release can contribute to a better understanding of K availability to plants. This study was conducted to evaluate K release rates from the whole soil, clay, silt, and sand fractions of B-horizon samples of a basalt-derived Oxisol and a sienite-derived Ultisol, both representative soils from coffee regions of Minas Gerais State, Brazil. Potassium was extracted from each fraction after eight different shaking time periods (0-665 h) with either 0.001 mol L-1 citrate or oxalate at a 1:10 solid:solution ratio. First-order, Elovich, zero-order, and parabolic diffusion equations were used to parameterize the time dependence of K release. For the Oxisol, the first-order equation fitted best to the experimental data of K release, with similar rates for all fractions and independent of the presence of citrate or oxalate in the extractant solution. For all studied Ultisol fractions, in which K release rates increased when extractions were performed with citrate solution, the Elovich model described K release kinetics most adequately. The highest potassium release rate of the Ultisol silt fraction was probably due to the transference of "non-exchangeable" K to the extractant solution, whereas in the Oxisol exchangeable potassium represented the main K source in all studied fractions.

Regular wave propagation out of noise in chemical active media

A pacemaker, regularly emitting chemical waves, is created out of noise when an excitable photosensitive Belousov-Zhabotinsky medium, strictly unable to autonomously initiate autowaves, is forced with a spatiotemporal patterned random illumination. These experimental observations are also reproduced numerically by using a set of reaction-diffusion equations for an activator-inhibitor model, and further analytically interpreted in terms of genuine coupling effects arising from parametric fluctuations. Within the same framework we also address situations of noise-sustained propagation in subexcitable media.

Excitability transitions and wave dynamics under spatiotemporal structured noise

We present an analytic and numerical study of the effects of external fluctuations in active media. Our analytical methodology transforms the initial stochastic partial differential equations into an effective set of deterministic reaction-diffusion equations. As a result we are able to explain and make quantitative predictions on the systematic and constructive effects of the noise, for example, target patterns created out of noise and traveling or spiral waves sustained by noise. Our study includes the case of realistic noises with temporal and spatial structures.

Spatial correlations and cross sections of clusters in the A + B -> 0 reaction

We consider the distribution of cross sections of clusters and the density-density correlation functions for the A+B¿0 reaction. We solve the reaction-diffusion equations numerically for random initial distributions of reactants. When both reactant species have the same diffusion coefficients the distribution of cross sections and the correlation functions scale with the diffusion length and obey superuniversal laws (independent of dimension). For different diffusion coefficients the correlation functions still scale, but the scaling functions depend on the dimension and on the diffusion coefficients. Furthermore, we display explicitly the peculiarities of the cluster-size distribution in one dimension.

Uniformity Study of Amorphous and Microcrystalline Silicon Thin Films Deposited on 10cmx10cm Glass Substrate using Hot Wire CVD Technique

The scaling up of the Hot Wire Chemical Vapor Deposition (HW-CVD) technique to large deposition area can be done using a catalytic net of equal spaced parallel filaments. The large area deposition limit is defined as the limit whenever a further increment of the catalytic net area does not affect the properties of the deposited film. This is the case when a dense catalytic net is spread on a surface considerably larger than that of the film substrate. To study this limit, a system able to hold a net of twelve wires covering a surface of about 20 cm x 20 cm was used to deposit amorphous (a-Si:H) and microcrystalline (μc-Si:H) silicon over a substrate of 10 cm x 10 cm placed at a filament-substrate distance ranging from 1 to 2 cm. The uniformity of the film thickness d and optical constants, n(x, λ) and α(x,¯hω), was studied via transmission measurements. The thin film uniformity as a function of the filament-substrate distance was studied. The experimental thickness profile was compared with the theoretical result obtained solving the diffusion equations. The optimization of the filament-substrate distance allowed obtaining films with inhomogeneities lower than ±2.5% and deposition rates higher than 1 nm/s and 4.5 nm/s for (μc-Si:H) and (a-Si:H), respectively.

Excitability transitions and wave dynamics under spatiotemporal structured noise

We present an analytic and numerical study of the effects of external fluctuations in active media. Our analytical methodology transforms the initial stochastic partial differential equations into an effective set of deterministic reaction-diffusion equations. As a result we are able to explain and make quantitative predictions on the systematic and constructive effects of the noise, for example, target patterns created out of noise and traveling or spiral waves sustained by noise. Our study includes the case of realistic noises with temporal and spatial structures.