Finite element modelling of heat transfer through permeable cracks in hydrothermal systems with upward throughflow

We use the finite element method to model the heat transfer phenomenon through permeable cracks in hydrothermal systems with upward throughflow. Since the finite element method is an approximate numerical method, the method must be validated before it is used to soh,e any new, kind of problem. However, the analytical solution, which can be used to validate the finite element method and other numerical methods, is rather limited in the literature, especially, for the problem considered here. Keeping this in mind, we have derived analytical solutions for the temperature distribution along the vertical axis of a crack in a fluid-saturated porous layer. After the finite element method is validated by comparing the numerical solution with the analytical solution for the same benchmark problem, it is used to investigate the pore-fluid flow and heat transfer in layered hydrothermal systems with vertical permeable cracks. The related analytical and numerical results have demonstrated that vertical cracks are effective and efficient members to transfer heat energy from the bottom section to the top section in hydrothermal systems with upward throughflow.

A simple solution of the laminar dam break wave

The sudden release of a mass of fluid in a channel generates a highly unsteady flow motion, called dam break wave. While industrial fluids exhibit sometimes non-Newtonian behaviours, the viscous fluid flow assumption remains a useful approximation for simplified analyses. In this study, new solutions of laminar dam break wave are proposed for a semi-infinite reservoir based upon the method of characteristics. The solutions yield simple explicit expressions of the wave front location, wave front celerity and instantaneous free-surface profiles that compare favourably with experimental observations. Both horizontal and sloping channel configurations are treated. The simplicity of the equations may allow future extension to more complicated fluid flows.

Propagation of nonstationary curved and stretched premixed flames with multistep reaction mechanisms

The propagation speed of a thin premixed flame disturbed by an unsteady fluid flow of a larger scale is considered. The flame may also have a general shape but the reaction zone is assumed to be thin compared to the flame thickness. Unlike in preceding publications, the presented asymptotic analysis is performed for a general multistep reaction mechanism and, at the same time, the flame front is curved by the fluid flow. The resulting equations define the propagation speed of disturbed flames in terms of the properties of undisturbed planar flames and the flame stretch. Special attention is paid to the near-equidiffusion limit. In this case, the flame propagation speed is shown to depend on the effective Zeldovich number Z(f) , and the flame stretch. Unlike the conventional Zeldovich number, the effective Zeldovich number is not necessarily linked directly to the activation energies of the reactions. Several examples of determining the effective Zeldovich number for reduced combustion mechanisms are given while, for realistic reactions, the effective Zeldovich number is determined from experiments. Another feature of the present approach is represented by the relatively simple asymptotic technique based on the adaptive generalized curvilinear system of coordinates attached to the flame (i.e., intrinsic disturbed flame equations [IDFE]).

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.

Efeitos tridimensionais no escoamento em torno de cilindros confinados

Mestrado em Engenharia Química

Recuperação de calor dos efluentes da tinturaria

Foram utilizados como elemento principal de estudo, os efluentes líquidos dos processos industriais da ENDUTEX, Tinturaria e Acabamento de Malhas, S. A. localizada no município de Caldas de Vizela, distrito de Braga. O estudo foi realizado na empresa devido ao interesse da mesma em poder reaproveitar o calor libertado nos efluentes para aquecimento de parte da água captada no rio de Vizela. O objectivo do trabalho consiste no dimensionamento de um permutador de calor que permita satisfazer o interesse da empresa, assim como, um estudo económico relativo aos custos envolventes. Com o intuito de concretizar os objectivos propostos foram realizadas visitas semanais à empresa para se proceder ao levantamento de dados e para a realização de amostragens do efluente para posterior caracterização. Depois de efectuado o dimensionamento do permutador de placas para diferentes caudais e temperaturas dos fluidos, frio (água do rio) e quente (efluentes), concluiu-se que as condições mais rentáveis correspondiam a um caudal de fluido frio de 17 m3/h em que a temperatura de entrada e de saída no permutador seria de 14 ºC e 48 ºC, respectivamente. O caudal de fluido quente seria de 20 m3/h, sendo a temperatura de entrada e de saída no permutador de 62 ºC e 33,1 ºC, respectivamente. Como resultado do dimensionamento obteve-se um permutador de placas com 167 placas em que o coeficiente global de transferência de calor (U) é de 726,9 W/m2ºC, a área projectada de 55,7 m2 e a queda de pressão de 0,904 KPa. Foi consultada a empresa ARSOPI-THERMAL para verificação das características dos permutadores existentes no mercado. No entanto, para as mesmas condições foi sugerido um permutador com 31 placas em que o coeficiente global de transferência de calor (U) é de 6267 W/m2ºC, a área projectada de 7,39 m2 e a queda de pressão de 76 KPa. A diferença verificada nos resultados apresentados pode ter origem na utilização de diferentes expressões no cálculo do coeficiente pelicular de transferência de calor (h) e pelo facto da ARSOPI desprezar o factor de sujamento no seu dimensionamento Na análise económica do projecto é de referir que para o arranque do projecto foi feito o levantamento das necessidades de investimento, situando-se este num valor total de 9640€, sendo o investimento financiado apenas por capitais próprios. O prazo de recuperação do investimento (Pay Back Period) é de cerca de 2 meses.

Systems of Navier-Stokes equations on cantor sets

We present systems of Navier-Stokes equations on Cantor sets, which are described by the local fractional vector calculus. It is shown that the results for Navier-Stokes equations in a fractal bounded domain are efficient and accurate for describing fluid flow in fractal media.

Compact SMB chromatography for binary separation

A thesis submitted for the degree of Doctor of Philosophy

Analysis of diffusion process in fractured reservoirs using fractional derivative approach

The fractal geometry is used to model of a naturally fractured reservoir and the concept of fractional derivative is applied to the diffusion equation to incorporate the history of fluid flow in naturally fractured reservoirs. The resulting fractally fractional diffusion (FFD) equation is solved analytically in the Laplace space for three outer boundary conditions. The analytical solutions are used to analyze the response of a naturally fractured reservoir considering the anomalous behavior of oil production. Several synthetic examples are provided to illustrate the methodology proposed in this work and to explain the diffusion process in fractally fractured systems.

Cilia motility studies in zebrafish embryos

A thesis submitted in fulfilment of the requirements for the degree of Masters in Molecular Genetics and Biomedicine

Using cilia mutants to study left-right asymmetry in zebrafish

A thesis submitted in fulfillment of the requirements for the degree of the Masters in Molecular Genetics and Biomedicine

Using the GPU to design complex profile extrusion dies

In the present work the benefits of using graphics processing units (GPU) to aid the design of complex geometry profile extrusion dies, are studied. For that purpose, a3Dfinite volume based code that employs unstructured meshes to solve and couple the continuity, momentum and energy conservation equations governing the fluid flow, together with aconstitutive equation, was used. To evaluate the possibility of reducing the calculation time spent on the numerical calculations, the numerical code was parallelized in the GPU, using asimple programing approach without complex memory manipulations. For verificationpurposes, simulations were performed for three benchmark problems: Poiseuille flow, lid-driven cavity flow and flow around acylinder. Subsequently, the code was used on the design of two real life extrusion dies for the production of a medical catheter and a wood plastic composite decking profile. To evaluate the benefits, the results obtained with the GPU parallelized code were compared, in terms of speedup, with a serial implementation of the same code, that traditionally runs on the central processing unit (CPU). The results obtained show that, even with the simple parallelization approach employed, it was possible to obtain a significant reduction of the computation times.

An 40Ar/39Ar, Rb/Sr, and stable isotope study of micas in low-grade fold-and-thrust belt: An example from the Swiss Helvetic Alps

White micas in carbonate-rich tectonites and a few other rock types of large thrusts in the Swiss Helvetic fold-and-thrust belt have been analyzed by Ar-40/Ar-39 and Rb/Sr techniques to better constrain the timing of Alpine deformation for this region. Incremental Ar-40/Ar-39 heating experiments of 25 weakly metamorphosed (anchizone to low greenschist) samples yield plateau and staircase spectra. We interpret most of the staircase release spectra result from variable mixtures of syntectonic (neoformed) and detrital micas. The range in dates obtained within individual spectra depends primarily on the duration of mica nucleation and growth, and relative proportions of neoformed and detrital mica. Rb/Sr analyses of 12 samples yield dates of ca. 10-39 Ma (excluding one anomalously young sample). These dates are slightly younger than the Ar-40/Ar-39 total gas dates obtained for the same samples. The Rb/ Sr dates were calculated using initial Sr-87/Sr-86 ratios obtained from the carbonate-dominated host rocks, which are higher than normal Mesozoic carbonate values due to exchange with fluids of higher Sr-87/Sr-86 ratios (and lower O-18/O-16 ratios). Model dates calculated using Sr-87/Sr-86 values typical of Mesozoic marine carbonates more closely approximate the Ar-40/Ar-39 total gas dates for most of the samples. The similarities of Rb/Sr and Ar-40/Ar-39 total gas dates are consistent with limited amounts of detrital mica in the samples. The delta(18)O values range from 24-15%. (VSMOW) for 2-6 mum micas and 27-16parts per thousand for the carbonate host rocks. The carbonate values are significantly lower than their protolith values due to localized fluid-rock interaction and fluid flow along most thrust surfaces. Although most calcite-mica pairs are not in oxygen isotope equilibrium at temperatures of ca. 200-400 degreesC, their isotopic fractionations are indicative of either 1) partial exchange between the minerals and a common external fluid, or 2) growth or isotopic exchange of the mica with the carbonate after the carbonate had isotopically exchanged with an external fluid. The geological significance of these results is not easily or uniquely determined, and exemplifies the difficulties inherent in dating very fine-grained micas of highly deformed tectonites in low-grade metamorphic terranes. Two generalizations can be made regarding the dates obtained from the Helvetic thrusts: 1) samples from the two highest thrusts (Mt. Gond and Sublage) have all of their Ar-40/Ar-39 steps above 20 Ma, and 2) most samples from the deepest Helvetic thrusts have steps (often accounting for more than 80% of Ar-39 release) between 15 and 25 Ma. These dates are consistent with the order of thrusting in the foreland-imbricating system and increase proportions of neoformed to detrital mica in the more metamorphosed hinterland and deeply buried portions of the nappe pile. Individual thrusts accommodated the majority of their displacement during their initial incorporation into the foreland-imbricating system, and some thrusts remained active or were reactivated down to 15 Ma.

Numerical simulation of combined heat transfer phenomena (conduction, convection, and radiation). Application to he optimisation of thermal systems

Thermal systems interchanging heat and mass by conduction, convection, radiation (solar and thermal ) occur in many engineering applications like energy storage by solar collectors, window glazing in buildings, refrigeration of plastic moulds, air handling units etc. Often these thermal systems are composed of various elements for example a building with wall, windows, rooms, etc. It would be of particular interest to have a modular thermal system which is formed by connecting different modules for the elements, flexibility to use and change models for individual elements, add or remove elements without changing the entire code. A numerical approach to handle the heat transfer and fluid flow in such systems helps in saving the full scale experiment time, cost and also aids optimisation of parameters of the system. In subsequent sections are presented a short summary of the work done until now on the orientation of the thesis in the field of numerical methods for heat transfer and fluid flow applications, the work in process and the future work.

The role of a transcrustal shear zone in orogenic gold mineralization at the Aijanahalli Mine, Dharwar Craton, South India

The Ajjanahalli gold mine is spatially associated with a Late Archean craton-scale shear zone in the eastern Chitradurga greenstone belt of the Dharwar craton, India. Gold mineralization is hosted by an similar to100-m-wide antiform in a banded iron formation. Original magnetite and siderite are replaced by a peak metamorphic alteration assemblage of chlorite, stilpnomelane, minnesotaite, sericite, ankerite, arsenopyrite, pyrite, pyrrhotite, and gold at ca. 300degrees to 350degreesC. Elements enriched in the banded iron formation include Ca, Mg, C, S, An, As, Bi. Cu, Sb, Zn, Pb, Se, Ag, and Te, whereas in the wall rocks As, Cu, Zn, Bi, Ag, and An are only slightly enriched. Strontium correlates with CaO, MgO, CO2, and As, which indicates cogenetic formation of arsenopyrite and Mg-Ca carbonates. The greater extent of alteration in the Fe-rich banded iron formation layers than in the wall rock reflects the greater reactivity of the banded iron formation layers. The ore fluids, as interpreted from their isotopic composition (delta(18)O = 6.5-8.5parts per thousand; initial Sr-87/Sr-86 = 0.7068-0.7078), formed by metamorphic devolatilization of deeper levels of the Chitradurga greenstone belt. Arsenopyrite, chalcopyrite, and pyrrhotite have delta(34)S values within a narrow range between 2.1 and 2.7 per mil, consistent with a sulfur source in Chitradurga greenstone belt lithologies. Based on spatial and temporal relationships between mineralization, local structure development, and sinistral strike-slip deformation in the shear zone at the eastern contact of the Chitradurga greenstone belt, we suggest that the Ajjanahalli gold mineralization formed by fluid infiltration into a low strain area within the first-order structure. The ore fluids were transported along this shear zone into relatively shallow crustal levels during lateral terrane accretion and a change from thrust to transcurrent tectonics. Based on this model of fluid flow, exploration should focus on similar low strain areas or potentially connected higher order splays of the first-order shear zone.