Anomalous water absorption in porous materials

The absorption of fluid by unsaturated, rigid porous materials may be characterized by the sorptivity. This is a simple parameter to determine and is increasingly being used as a measure of a material's resistance to exposure to fluids (especially moisture and reactive solutes) in aggressive environments. The complete isothermal absorption process is described by a nonlinear diffusion equation, with the hydraulic diffusivity being a strongly nonlinear function of the degree of saturation of the material. This diffusivity can be estimated from the sorptivity test. In a typical test the cumulative absorption is proportional to the square root of time. However, a number of researchers have observed deviation from this behaviour when the infiltrating fluid is water and there is some potential for chemo-mechanical interaction with the material. In that case the current interpretation of the test and estimation of the hydraulic diffusivity is no longer appropriate. Kuntz and Lavallee (2001) discuss the anomalous behaviour and propose a non-Darcian model as a more appropriate physical description. We present an alternative Darcian explanation and theory that retrieves the earlier advantages of the simple sorptivity test in providing parametric information about the material's hydraulic properties and allowing simple predictive formulae for the wetting profile to be generated.

Wall mediated transport in confined spaces: Exact theory for low density

We present a theory for the transport of molecules adsorbed in slit and cylindrical nanopores at low density, considering the axial momentum gain of molecules oscillating between diffuse wall reflections. Good agreement with molecular dynamics simulations is obtained over a wide range of pore sizes, including the regime of single-file diffusion where fluid-fluid interactions are shown to have a negligible effect on the collective transport coefficient. We show that dispersive fluid-wall interactions considerably attenuate transport compared to classical hard sphere theory.

Molecular transport in nanopores

Simulation of the transport of methane in cylindrical silica mesopores have been performed using equilibrium and nonequilibrium molecular dynamics (NEMD) as well as dual control volume grand canonical molecular dynamics methods. It is demonstrated that all three techniques yield the same transport coefficient even in the presence of viscous flow. A modified locally averaged density model for viscous flow, combined with consideration of wall slip through a frictional condition, gives a convincing interpretation of the variation of the transport coefficient over a wide range of densities, and for various pore sizes and temperatures. Wall friction coefficients extracted from NEMD simulations are found to be consistent with momentum transfer arguments, and the approach is shown to be more meaningful than the classical slip length concept. (C) 2003 American Institute of Physics.

Herbicidal activity of volatiles from coriander, winter savory, cotton lavender, and thyme isolated by hydrodistillation and supercritical fluid extraction

The volatiles from Coriandrum sativum L., Satureja montana L., Santolina chamaecyparissus L., and Thymus vulgaris L. were isolated by hydrodistillation (essential oil) and supercritical fluid extraction (volatile oil). Their effect on seed germination and root and shoot growth of the surviving seedlings of four crops (Zea mays L., Triticum durum L., Pisum sativum L., and Lactuca sativa L.) and two weeds (Portulaca oleracea L. and Vicia sativa L.) was investigated and compared with those of two synthetic herbicides, Agrocide and Prowl. The volatile oils of thyme and cotton lavender seemed to be promising alternatives to the synthetic herbicides because they were the least injurious to the crop species. The essential oil of winter savory, on the other hand, affected both crop and weeds and can be appropriate for uncultivated fields.

Composition and antioxidant activity of thymus vulgaris volatiles: comparison between supercritical fluid extraction and hydrodistillation

Supercritical fluid extraction (SEE) of the volatile oil from Thymus vulgaris L. aerial flowering parts was performed under different conditions of pressure, temperature, mean particle size and CO2 flow rate and the correspondent yield and composition were compared with those of the essential oil isolated by hydrodistillation (HD). Both the oils were analyzed by GC and GC-MS and 52 components were identified. The main volatile components obtained were p-cymene (10.0-42.6% for SFE and 28.9-34.8% for HD), gamma-terpinene (0.8-6.9% for SFE and 5.1-7.0% for HD), linalool (2.3-5.3% for SFE and 2.8-3.1% for HD), thymol (19.5-40.8% for SFE and 35.4-41.6% for HD), and carvacrol (1.4-3.1% for SFE and 2.6-3.1% for HD). The main difference was found to be the relative percentage of thymoquinone (not found in the essential oil) and carvacryl methyl ether (1.0-1.2% for HD versus t-0.4 for SFE) which can explain the higher antioxidant activity, assessed by Rancimat test, of the SFE volatiles when compared with HD. Thymoquinone is considered a strong antioxidant compound.

The condensation and ordering of models of empty liquids

We consider a simple model consisting of particles with four bonding sites ("patches"), two of type A and two of type B, on the square lattice, and investigate its global phase behavior by simulations and theory. We set the interaction between B patches to zero and calculate the phase diagram as the ratio between the AB and the AA interactions, epsilon(AB)*, varies. In line with previous work, on three-dimensional off-lattice models, we show that the liquid-vapor phase diagram exhibits a re-entrant or "pinched" shape for the same range of epsilon(AB)*, suggesting that the ratio of the energy scales - and the corresponding empty fluid regime - is independent of the dimensionality of the system and of the lattice structure. In addition, the model exhibits an order-disorder transition that is ferromagnetic in the re-entrant regime. The use of low-dimensional lattice models allows the simulation of sufficiently large systems to establish the nature of the liquid-vapor critical points and to describe the structure of the liquid phase in the empty fluid regime, where the size of the "voids" increases as the temperature decreases. We have found that the liquid-vapor critical point is in the 2D Ising universality class, with a scaling region that decreases rapidly as the temperature decreases. The results of simulations and theoretical analysis suggest that the line of order-disorder transitions intersects the condensation line at a multi-critical point at zero temperature and density, for patchy particle models with a re-entrant, empty fluid, regime. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3657406]

New physics contributions to A(FB)(t(t)over-bar) at the Tevatron

The Tevatron has measured a discrepancy relative to the standard model prediction in the forward-backward asymmetry in top quark pair production. This asymmetry grows with the rapidity difference of the two top quarks. It also increases with the invariant mass of the t (t) over bar pair, reaching, for high invariant masses, 3.4 standard deviations above the next-to-leading order prediction for the charge asymmetry of QCD. However, perfect agreement between experiment and the standard model was found in both total and differential cross section of top quark pair production. As this result could be a sign of new physics we have parametrized this new physics in terms of a complete set of dimension six operators involving the top quark. We have then used a Markov chain Monte Carlo approach in order to find the best set of parameters that fits the data, using all available data regarding top quark pair production at the Tevatron. We have found that just a very small number of operators are able to fit the data better than the standard model.

Phase diagrams of particles with dissimilar patches: X-junctions and Y-junctions

We use Wertheim's first-order perturbation theory to investigate the phase behaviour and the structure of coexisting fluid phases for a model of patchy particles with dissimilar patches (two patches of type A and f(B) patches of type B). A patch of type alpha = {A, B} can bond to a patch of type beta = {A, B} in a volume nu(alpha beta), thereby decreasing the internal energy by epsilon(alpha beta). We analyse the range of model parameters where AB bonds, or Y-junctions, are energetically disfavoured (epsilon(AB) < epsilon(AA)/2) but entropically favoured (nu(AB) >> nu(alpha alpha)), and BB bonds, or X-junctions, are energetically favoured (epsilon(BB) > 0). We show that, for low values of epsilon(BB)/epsilon(AA), the phase diagram has three different regions: (i) close to the critical temperature a low-density liquid composed of long chains and rich in Y-junctions coexists with a vapour of chains; (ii) at intermediate temperatures there is coexistence between a vapour of short chains and a liquid of very long chains with X-and Y-junctions; (iii) at low temperatures an ideal gas coexists with a high-density liquid with all possible AA and BB bonds formed. It is also shown that in region (i) the liquid binodal is reentrant (its density decreases with decreasing temperature) for the lower values of epsilon(BB)/epsilon(AA). The existence of these three regions is a consequence of the competition between the formation of X- and Y-junctions: X-junctions are energetically favoured and thus dominate at low temperatures, whereas Y-junctions are entropically favoured and dominate at higher temperatures.

Patching up dipoles: can dipolar particles be viewed as patchy colloids?

We investigate whether the liquid-vapour phase transition of strongly dipolar fluids can be understood using a model of patchy colloids. These consist of hard spherical particles with three short-ranged attractive sites (patches) on their surfaces. Two of the patches are of type A and one is of type B. Patches A on a particle may bond either to a patch A or to a patch B on another particle. Formation of an AA (AB) bond lowers the energy by epsilon AA (epsilon AB). In the limit [image omitted], this patchy model exhibits condensation driven by AB-bonds (Y-junctions). Y-junctions are also present in low-density, strongly dipolar fluids, and have been conjectured to play a key role in determining their critical behaviour. We map the dipolar Yukawa hard-sphere (DYHS) fluid onto this 2A + 1B patchy model by requiring that the latter reproduce the correct DYHS critical point as a function of the isotropic interaction strength epsilon Y. This is achieved for sensible values of epsilon AB and the bond volumes. Results for the internal energy and the particle coordination number are in qualitative agreement with simulations of DYHSs. Finally, by taking the limit [image omitted], we arrive at a new estimate for the critical point of the dipolar hard-sphere fluid, which agrees with extrapolations from simulation.

Re-entrant phase behaviour of network fluids: A patchy particle model with temperature-dependent valence

We study a model consisting of particles with dissimilar bonding sites ("patches"), which exhibits self-assembly into chains connected by Y-junctions, and investigate its phase behaviour by both simulations and theory. We show that, as the energy cost epsilon(j) of forming Y-junctions increases, the extent of the liquid-vapour coexistence region at lower temperatures and densities is reduced. The phase diagram thus acquires a characteristic "pinched" shape in which the liquid branch density decreases as the temperature is lowered. To our knowledge, this is the first model in which the predicted topological phase transition between a fluid composed of short chains and a fluid rich in Y-junctions is actually observed. Above a certain threshold for epsilon(j), condensation ceases to exist because the entropy gain of forming Y-junctions can no longer offset their energy cost. We also show that the properties of these phase diagrams can be understood in terms of a temperature-dependent effective valence of the patchy particles. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3605703]

Study of the one dimensional and transient bioheat transfer equation: Multi-layer solution development and applications

: In this work we derive an analytical solution given by Bessel series to the transient and one-dimensional (1D) bioheat transfer equation in a multi-layer region with spatially dependent heat sources. Each region represents an independent biological tissue characterized by temperature-invariant physiological parameters and a linearly temperature dependent metabolic heat generation. Moreover, 1D Cartesian, cylindrical or spherical coordinates are used to define the geometry and temperature boundary conditions of first, second and third kinds are assumed at the inner and outer surfaces. We present two examples of clinical applications for the developed solution. In the first one, we investigate two different heat source terms to simulate the heating in a tumor and its surrounding tissue, induced during a magnetic fluid hyperthermia technique used for cancer treatment. To obtain an accurate analytical solution, we determine the error associated with the truncated Bessel series that defines the transient solution. In the second application, we explore the potential of this model to study the effect of different environmental conditions in a multi-layered human head model (brain, bone and scalp). The convective heat transfer effect of a large blood vessel located inside the brain is also investigated. The results are further compared with a numerical solution obtained by the Finite Element Method and computed with COMSOL Multi-physics v4.1 (c). (c) 2013 Elsevier Ltd. All rights reserved.

Effective isotropic potential for dipolar hard spheres

A new effective isotropic potential is proposed for the dipolar hard-sphere fluid, on the basis of recent results by others for its angle-averaged radial distribution function. The new effective potential is shown to exhibit oscillations even for moderately high densities and moderately strong dipole moments, which are absent from earlier effective isotropic potentials. The validity and significance of this result are briefly discussed.

Coletores de distribuição e sua otimização hidráulica, aplicados a uma instalação de produção de água arrefecida em sistemas de AVAC ( Fórum Picoas)

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Mecânica

Biaxial Nematic Order in the Hard-boomerang Fluid

We consider a fluid of hard boomerangs, each composed of two hard spherocylinders joined at their ends at an angle Psi. The resulting particle is nonconvex and biaxial. The occurence of nematic order in such a system has been investigated using Straley's theory, which is a simplificaton of Onsager's second-virial treatment of long hard rods, and by bifurcation analysis. The excluded volume of two hard boomerangs has been approximated by the sum of excluded volumes of pairs of constituent spherocylinders, and the angle-dependent second-virial coefficient has been replaced by a low-order interpolating function. At the so-called Landau point, Psi(Landau)approximate to 107.4 degrees, the fluid undergoes a continuous transition from the isotropic to a biaxial nematic (B) phase. For Psi not equal Psi(Landau) ordering is via a first-order transition into a rod-like uniaxial nematic phase (N(+)) if Psi > Psi(Landau), or a plate-like uniaxial nematic (N(-)) phase if Psi < Psi(Landau). The B phase is separated from the N(+) and N(-) phases by two lines of continuous transitions meeting at the Landau point. This topology of the phase diagram is in agreement with previous studies of spheroplatelets and biaxial ellipsoids. We have checked the accuracy of our theory by performing numerical calculations of the angle-dependent second virial coefficient, which yields Psi(Landau)approximate to 110 degrees for very long rods, and Psi(Landau)approximate to 90 degrees for short rods. In the latter case, the I-N transitions occur at unphysically high packing fractions, reflecting the inappropriateness of the second-virial approximation in this limit.

Modellus: Learning Physics with Mathematical Modelling

Tese de doutoramento em Ciências da Educação, área de Teoria Curricular e Ensino das Ciências