Thermal models, lithosphere thickness and heat flow in South Portugal. Some comments about the subject

Reinterpretation of old heat flow data or use of new data and new techniques of detection of the temperature under the surface have conducted to new heat flow density values in some regions of the globe. The problem of ice melting in Greenland and Antarctica caught the public's attention to the importance of knowledge on heat flow values and thermal structure of the globe. In the last years, several models were presented trying to obtain lithosphere and Moho thickness of the Iberia Peninsula. The work we intend to present is related with the SW part of the Iberia Peninsula ( south of the Ossa Morena zone, South Portuguese Zone and Algarve). The results obtained show a decrease in the thickness of the crust and the lithosphere in this region. Density anomalies in the crust are also referred. I intend to make the connection between the results of these models and the heat flow thermal conductivity, heat production and geological data available for the region, trying to explain the results of heat flow density data obtained.

Adaptivity and a posteriori error control for bifurcation problems II: Incompressible fluid flow in open systems with Z_2 symmetry

In this article we consider the a posteriori error estimation and adaptive mesh refinement of discontinuous Galerkin finite element approximations of the bifurcation problem associated with the steady incompressible Navier-Stokes equations. Particular attention is given to the reliable error estimation of the critical Reynolds number at which a steady pitchfork or Hopf bifurcation occurs when the underlying physical system possesses reflectional or Z_2 symmetry. Here, computable a posteriori error bounds are derived based on employing the generalization of the standard Dual-Weighted-Residual approach, originally developed for the estimation of target functionals of the solution, to bifurcation problems. Numerical experiments highlighting the practical performance of the proposed a posteriori error indicator on adaptively refined computational meshes are presented.

Numerical simulation of fresh SCC flow in wall and beam elements using flow dynamics models

Abstract : Recently, there is a great interest to study the flow characteristics of suspensions in different environmental and industrial applications, such as snow avalanches, debris flows, hydrotransport systems, and material casting processes. Regarding rheological aspects, the majority of these suspensions, such as fresh concrete, behave mostly as non-Newtonian fluids. Concrete is the most widely used construction material in the world. Due to the limitations that exist in terms of workability and formwork filling abilities of normal concrete, a new class of concrete that is able to flow under its own weight, especially through narrow gaps in the congested areas of the formwork was developed. Accordingly, self-consolidating concrete (SCC) is a novel construction material that is gaining market acceptance in various applications. Higher fluidity characteristics of SCC enable it to be used in a number of special applications, such as densely reinforced sections. However, higher flowability of SCC makes it more sensitive to segregation of coarse particles during flow (i.e., dynamic segregation) and thereafter at rest (i.e., static segregation). Dynamic segregation can increase when SCC flows over a long distance or in the presence of obstacles. Therefore, there is always a need to establish a trade-off between the flowability, passing ability, and stability properties of SCC suspensions. This should be taken into consideration to design the casting process and the mixture proportioning of SCC. This is called “workability design” of SCC. An efficient and non-expensive workability design approach consists of the prediction and optimization of the workability of the concrete mixtures for the selected construction processes, such as transportation, pumping, casting, compaction, and finishing. Indeed, the mixture proportioning of SCC should ensure the construction quality demands, such as demanded levels of flowability, passing ability, filling ability, and stability (dynamic and static). This is necessary to develop some theoretical tools to assess under what conditions the construction quality demands are satisfied. Accordingly, this thesis is dedicated to carry out analytical and numerical simulations to predict flow performance of SCC under different casting processes, such as pumping and tremie applications, or casting using buckets. The L-Box and T-Box set-ups can evaluate flow performance properties of SCC (e.g., flowability, passing ability, filling ability, shear-induced and gravitational dynamic segregation) in casting process of wall and beam elements. The specific objective of the study consists of relating numerical results of flow simulation of SCC in L-Box and T-Box test set-ups, reported in this thesis, to the flow performance properties of SCC during casting. Accordingly, the SCC is modeled as a heterogeneous material. Furthermore, an analytical model is proposed to predict flow performance of SCC in L-Box set-up using the Dam Break Theory. On the other hand, results of the numerical simulation of SCC casting in a reinforced beam are verified by experimental free surface profiles. The results of numerical simulations of SCC casting (modeled as a single homogeneous fluid), are used to determine the critical zones corresponding to the higher risks of segregation and blocking. The effects of rheological parameters, density, particle contents, distribution of reinforcing bars, and particle-bar interactions on flow performance of SCC are evaluated using CFD simulations of SCC flow in L-Box and T-box test set-ups (modeled as a heterogeneous material). Two new approaches are proposed to classify the SCC mixtures based on filling ability and performability properties, as a contribution of flowability, passing ability, and dynamic stability of SCC.

Simulação computacional do escoamento em bombas de cavidades progressivas

The use of Progressing Cavity Pumps (PCPs) in artificial lift applications in low deep wells is becoming more common in the oil industry, mainly, due to its ability to pump heavy oils, produce oil with large concentrations of sand, besides present high efficiency when compared to other artificial lift methods. Although this system has been widely used as an oil lift method, few investigations about its hydrodynamic behavior are presented, either experimental or numeric. Therefore, in order to increase the knowledge about the BCP operational behavior, this work presents a novel computational model for the 3-D transient flow in progressing cavity pumps, which includes the relative motion between rotor and stator, using an element based finite volume method. The model developed is able to accurately predict the volumetric efficiency and viscous looses as well as to provide detailed information of pressure and velocity fields inside the pump. In order to predict PCP performance for low viscosity fluids, advanced turbulence models were used to treat, accurately, the turbulent effects on the flow, which allowed for obtaining results consistent with experimental values encountered in literature. In addition to the 3D computational model, a simplified model was developed, based on mass balance within cavities and on simplification on the momentum equations for fully developed flow along the seal region between cavities. This simplified model, based on previous approaches encountered in literature, has the ability to predict flow rate for a given differential pressure, presenting exactness and low CPU requirements, becoming an engineering tool for quick calculations and providing adequate results, almost real-time time. The results presented in this work consider a rigid stator PCP and the models developed were validated against experimental results from open literature. The results for the 3-D model showed to be sensitive to the mesh size, such that a numerical mesh refinement study is also presented. Regarding to the simplified model, some improvements were introduced in the calculation of the friction factor, allowing the application fo the model for low viscosity fluids, which was unsuccessful in models using similar approaches, presented in previous works

Compressible modes in a square lid-driven cavity

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

FLOW REGIME DRIVEN THERMAL ENHANCEMENT IN INTERNALLY-GROOVED TUBES

Internally-grooved refrigeration tubes maximize tube-side evaporative heat transfer rates and have been identified as a most promising technology for integration into compact cold plates. Unfortunately, the absence of phenomenological insights and physical models hinders the extrapolation of grooved-tube performance to new applications. The success of regime-based heat transfer correlations for smooth tubes has motivated the current effort to explore the relationship between flow regimes and enhanced heat transfer in internally-grooved tubes. In this thesis, a detailed analysis of smooth and internally-grooved tube data reveals that performance improvement in internally-grooved tubes at low-to-intermediate mass flux is a result of early flow regime transition. Based on this analysis, a new flow regime map and corresponding heat transfer coefficient correlation, which account for the increased wetted angle, turbulence, and Gregorig effects unique to internally-grooved tubes, were developed. A two-phase test facility was designed and fabricated to validate the newly-developed flow regime map and regime-based heat transfer coefficient correlation. As part of this setup, a non-intrusive optical technique was developed to study the dynamic nature of two-phase flows. It was found that different flow regimes result in unique temporally varying film thickness profiles. Using these profiles, quantitative flow regime identification measures were developed, including the ability to explain and quantify the more subtle transitions that exist between dominant flow regimes. Flow regime data, based on the newly-developed method, and heat transfer coefficient data, using infrared thermography, were collected for two-phase HFE-7100 flow in horizontal 2.62mm - 8.84mm diameter smooth and internally-grooved tubes with mass fluxes from 25-300 kg/m²s, heat fluxes from 4-56 kW/m², and vapor qualities approaching 1. In total, over 6500 combined data points for the adiabatic and diabatic smooth and internally-grooved tubes were acquired. Based on results from the experiments and a reinterpretation of data from independent researchers, it was established that heat transfer enhancement in internally-grooved tubes at low-to-intermediate mass flux is primarily due to early flow regime transition to Annular flow. The regime-based heat transfer coefficient outperformed empirical correlations from the literature, with mean and absolute deviations of 4.0% and 32% for the full range of data collected.

Effects of hot intrusions on pore-fluid flow and heat transfer in fluid-saturated rocks

We use the finite element method to solve coupled problems between pore-fluid flow and heat transfer in fluid-saturated porous rocks. In particular, we investigate the effects of both the hot pluton intrusion and topographically driven horizontal flow on the distributions of the pore-flow velocity and temperature in large-scale hydrothermal systems. Since general mineralization patterns are strongly dependent on distributions of both the pore-fluid velocity and temperature fields, the modern mineralization theory has been used to predict the general mineralization patterns in several realistic hydrothermal systems. The related numerical results have demonstrated that: (1) The existence of a hot intrusion can cause an increase in the maximum value of the pore-fluid velocity in the hydrothermal system. (2) The permeability of an intruded pluton is one of the sensitive parameters to control the pore-fluid flow, heat transfer and ore body formation in hydrothermal systems. (3) The maximum value of the pore-fluid velocity increases when the bottom temperature of the hydrothermal system is increased. (4) The topographically driven flow has significant effects on the pore-fluid flow, temperature distribution and precipitation pattern of minerals in hydrothermal systems. (5) The size of the computational domain may have some effects on the pore-fluid flow and heat transfer, indicating that the size of a hydrothermal system may affect the pore-fluid flow and heat transfer within the system. (C) 2003 Elsevier Science B.V. All rights reserved.

Simulation of mixing in unsteady flow through a periodically obstructed channel

Fluid mixing in steady and unsteady Bow through a channel containing periodic square obstructions has been studied using a finite-difference simulation to determine fluid velocities, followed by the use of passive marker particle advection to look at fluid transport out of the cavities formed between each of the obstructions. The geometry and Bow conditions were chosen from the work by Perkins (1989, M.S. Thesis, Lehigh University; 1992, Ph.D. Thesis, Lehigh University); who investigated heat transfer enhancement due to unsteady flow through such an obstructed channel. Particle advection shows that Bow regimes which are predicted to give good mixing based on snapshots of instantaneous streamline contour plots were not necessarily able to efficiently mix fluid which started in the cavity regions throughout the channel. The use of Poincare sections shows regular regions existing under these conditions which inhibit efficient fluid transport. These regular regions are found to disappear when the unsteady Bow velocity is increased. (C) 1997 Elsevier Science Ltd.

Spontaneous Emergence of Chirality in the Limited Enantioselectivity Model: Autocatalytic Cycle Driven by an External Reagent

The emergence of chirality in enantioselective autocatalysis for compounds unable to transform according to the Frank-like reaction network is discussed with respect to the controversial limited enantioselectivity (LES) model composed of coupled enantioselective and non-enantioselective autocatalyses. The LES model cannot lead to spontaneous mirror symmetry breaking (SMSB) either in closed systems with a homogeneous temperature distribution or in closed systems with a stationary non-uniform temperature distribution. However, simulations of chemical kinetics in a two-compartment model demonstrate that SMSB may occur if both autocatalytic reactions are spatially separated at different temperatures in different compartments but coupled under the action of a continuous internal flow. In such conditions, the system can evolve, for certain reaction and system parameters, toward a chiral stationary state; that is, the system is able to reach a bifurcation point leading to SMSB. Numerical simulations in which reasonable chemical parameters have been used suggest that an adequate scenario for such a SMSB would be that of abyssal hydrothermal vents, by virtue of the typical temperature gradients found there and the role of inorganic solids mediating chemical reactions in an enzyme-like role.

Internal material logistics development in customer-driven production

Työn tavoitteena oli synkronoida yleinen materiaalivirta tuotannon kanssa. Pyrkimyksenä oli virtaviivaistaa materiaalivirta, jotta materiaalin määrä, sekä turhat työvaiheet tuotantoalueella vähenisivät. Päämääränä oli vähentää materiaalilavojen määrää lattialla, sekä parantaa materiaalitoimitusten OTD (On-Time-Delivery) prosenttia.Teoreettinen osa käsittelee nykypäivän toimitusketjun tärkeimpiä elementtejä. Keskeisenä asiana perehdytään materiaalivirtaan ja sen tehokkaaseen hallintaan. Työ esittelee myös synkronisen johtamismallin periaatteet, sekä materiaalivirran synkronoimisen tuotannon kanssa.Empiirinen osuus kuvaa yrityksen materiaalihallinnan nykypäivän tilanteen, sen keskeiset ongelmat, sekä uuden toimintamalliehdotuksen. Työ esittelee kaksi pilottia, joiden tulokset varmistivat uuden virtaviivaistetun materiaalivirran ja täydennysmallin toimivuuden.Työn tulokset osoittavat, miten uusi materiaalin täydennysmalli vähentää materiaalilavojen, sekä jalostamattoman työn määrää tuotantoalueella.

A pulsatile flow model for in vitro quantitative evaluation of prosthetic valve regurgitation

A pulsatile pressure-flow model was developed for in vitro quantitative color Doppler flow mapping studies of valvular regurgitation. The flow through the system was generated by a piston which was driven by stepper motors controlled by a computer. The piston was connected to acrylic chambers designed to simulate "ventricular" and "atrial" heart chambers. Inside the "ventricular" chamber, a prosthetic heart valve was placed at the inflow connection with the "atrial" chamber while another prosthetic valve was positioned at the outflow connection with flexible tubes, elastic balloons and a reservoir arranged to mimic the peripheral circulation. The flow model was filled with a 0.25% corn starch/water suspension to improve Doppler imaging. A continuous flow pump transferred the liquid from the peripheral reservoir to another one connected to the "atrial" chamber. The dimensions of the flow model were designed to permit adequate imaging by Doppler echocardiography. Acoustic windows allowed placement of transducers distal and perpendicular to the valves, so that the ultrasound beam could be positioned parallel to the valvular flow. Strain-gauge and electromagnetic transducers were used for measurements of pressure and flow in different segments of the system. The flow model was also designed to fit different sizes and types of prosthetic valves. This pulsatile flow model was able to generate pressure and flow in the physiological human range, with independent adjustment of pulse duration and rate as well as of stroke volume. This model mimics flow profiles observed in patients with regurgitant prosthetic valves.

Numerical Simulation of Electroosmotic Flow with Step Change in Zeta Potential

Electroosmotic flow is a convenient mechanism for transporting polar fluid in a microfluidic device. The flow is generated through the application of an external electric field that acts on the free charges that exists in a thin Debye layer at the channel walls. The charge on the wall is due to the chemistry of the solid-fluid interface, and it can vary along the channel, e.g. due to modification of the wall. This investigation focuses on the simulation of the electroosmotic flow (EOF) profile in a cylindrical microchannel with step change in zeta potential. The modified Navier-Stoke equation governing the velocity field and a non-linear two-dimensional Poisson-Boltzmann equation governing the electrical double-layer (EDL) field distribution are solved numerically using finite control-volume method. Continuities of flow rate and electric current are enforced resulting in a non-uniform electrical field and pressure gradient distribution along the channel. The resulting parabolic velocity distribution at the junction of the step change in zeta potential, which is more typical of a pressure-driven velocity flow profile, is obtained.

Passive urban ventilation by combined buoyancy-driven slope flow and wall flow: Parametric CFD studies on idealized city models

This paper reports the results of a parametric CFD study on idealized city models to investigate the potential of slope flow in ventilating a city located in a mountainous region when the background synoptic wind is absent. Examples of such a city include Tokyo in Japan, Los Angeles and Phoenix in the US, and Hong Kong. Two types of buoyancy-driven flow are considered, i.e., slope flow from the mountain slope (katabatic wind at night and anabatic wind in the daytime), and wall flow due to heated/cooled urban surfaces. The combined buoyancy-driven flow system can serve the purpose of dispersing the accumulated urban air pollutants when the background wind is weak or absent. The microscopic picture of ventilation performance within the urban structures was evaluated in terms of air change rate (ACH) and age of air. The simulation results reveal that the slope flow plays an important role in ventilating the urban area, especially in calm conditions. Katabatic flow at night is conducive to mitigating the nocturnal urban heat island. In the present parametric study, the mountain slope angle and mountain height are assumed to be constant, and the changing variables are heating/cooling intensity and building height. For a typical mountain of 500 m inclined at an angle of 20° to the horizontal level, the interactive structure is very much dependent on the ratio of heating/cooling intensity as well as building height. When the building is lower than 60 m, the slope wind dominates. When the building is as high as 100 m, the contribution from the urban wall flow cannot be ignored. It is found that katabatic wind can be very beneficial to the thermal environment as well as air quality at the pedestrian level. The air change rate for the pedestrian volume can be as high as 300 ACH.

Predicting natural ventilation in a two-zone building driven by combined forces.

Natural ventilation relies on less controllable natural forces so that it needs more artificial control, and thus its prediction, design and analysis become more important. This paper presents both theoretical and numerical simulations for predicting the natural ventilation flow in a two-zone building with multiple openings which is subjected to the combined natural forces. To our knowledge, this is the first analytical solutions obtained so far for a building with more than one zones and in each zone with possibly more than 2 openings. The analytical solution offers a possibility for validating a multi-zone airflow program. A computer program MIX is employed to conduct the numerical simulation. Good agreement is achieved. Different airflow modes are identified and some design recommendations are also provided.

Effects of water flow on ablation rate and morphological changes in human enamel and dentin after Er:YAG laser irradiation

Purpose: To investigate the laboratory effect of Er:YAG laser on ablation rate and morphological changes in human enamel and dentin with varying water flow. Methods: 23 human third molars were sectioned in mesio-distal and buccal-lingual directions. The slabs were flattened and weighted on an analytical laboratory balance (control). A 4-mm(2) area was demarcated and the samples were randomly assigned into three groups according to water flow employed during the laser irradiation (1.0, 1.5, and 2.0 mL/minute). An Er:YAG laser was used to ablate enamel (80.22-J/cm(2), 300 mJ/4Hz) and dentin (96.26-J/cm(2), 250 mJ/4Hz). After irradiation, the samples were immersed in distilled water for 1 hour and then weighted again. The final mass was obtained and laser-irradiated substrate mass loss was calculated by the difference between the initial and final mass. Afterwards, specimens were prepared for SEM. Results: Data were submitted to ANOVA and Tukey's test (P< 0.05). It was observed that the 2.0 mL/minute resulted in a higher mass loss, 1.0 mL/minute showed a lower mass loss, and 1.5 mL/minute demonstrated intermediate results (P< 0.05). The increase of water flow promoted less melting areas and cracks. Furthermore, dentin was more ablated than enamel. It may be concluded that the water flow of Er:YAG laser and the substrates affected the ablation rate. Among the tested parameters, 2.0 mL/minute improved the ability of ablation in enamel and dentin, with less morphologic surface alteration. (Am J Dent 20 12;25:332-336).