Finite element analysis of flow patterns near geological lenses in hydrodynamic and hydrothermal systems

We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method. Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately. Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.

Analytical solution to the zero-inertia problem for surge flow phenomena in nonprismatic channels

Surge flow phenomena. e.g.. as a consequence of a dam failure or a flash flood, represent free boundary problems. ne extending computational domain together with the discontinuities involved renders their numerical solution a cumbersome procedure. This contribution proposes an analytical solution to the problem, It is based on the slightly modified zero-inertia (ZI) differential equations for nonprismatic channels and uses exclusively physical parameters. Employing the concept of a momentum-representative cross section of the moving water body together with a specific relationship for describing the cross sectional geometry leads, after considerable mathematical calculus. to the analytical solution. The hydrodynamic analytical model is free of numerical troubles, easy to run, computationally efficient. and fully satisfies the law of volume conservation. In a first test series, the hydrodynamic analytical ZI model compares very favorably with a full hydrodynamic numerical model in respect to published results of surge flow simulations in different types of prismatic channels. In order to extend these considerations to natural rivers, the accuracy of the analytical model in describing an irregular cross section is investigated and tested successfully. A sensitivity and error analysis reveals the important impact of the hydraulic radius on the velocity of the surge, and this underlines the importance of an adequate description of the topography, The new approach is finally applied to simulate a surge propagating down the irregularly shaped Isar Valley in the Bavarian Alps after a hypothetical dam failure. The straightforward and fully stable computation of the flood hydrograph along the Isar Valley clearly reflects the impact of the strongly varying topographic characteristics on the How phenomenon. Apart from treating surge flow phenomena as a whole, the analytical solution also offers a rigorous alternative to both (a) the approximate Whitham solution, for generating initial values, and (b) the rough volume balance techniques used to model the wave tip in numerical surge flow computations.

Effect of material anisotropy on the onset of convective flow in three-dimensional fluid-saturated faults

In order to investigate the effect of material anisotropy on convective instability of three-dimensional fluid-saturated faults, an exact analytical solution for the critical Rayleigh number of three-dimensional convective flow has been obtained. Using this critical Rayleigh number, effects of different permeability ratios and thermal conductivity ratios on convective instability of a vertically oriented three-dimensional fault have been examined in detail. It has been recognized that (1) if the fault material is isotropic in the horizontal direction, the horizontal to vertical permeability ratio has a significant effect on the critical Rayleigh number of the three-dimensional fault system, but the horizontal to vertical thermal conductivity ratio has little influence on the convective instability of the system, and (2) if the fault material is isotropic in the fault plane, the thermal conductivity ratio of the fault normal to plane has a considerable effect on the critical Rayleigh number of the three-dimensional fault system, but the effect of the permeability ratio of the fault normal to plane on the critical Rayleigh number of three-dimensional convective flow is negligible.

Numerical simulation of double-diffusion driven convective flow and rock alteration in three-dimensional fluid-saturated geological fault zones

Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. (c) 2005 Elsevier B.V. All rights reserved.

Hydraulic tracer studies in a pilot scale subsurface flow constructed wetland

Current design procedures for Subsurface Flow (SSF) Wetlands are based on the simplifying assumptions of plug flow and first order decay of pollutants. These design procedures do yield functional wetlands but result in over-design and inadequate descriptions of the pollutant removal mechanisms which occur within them. Even though these deficiencies are often noted, few authors have attempted to improve modelling of either flow or pollutant removal in such systems. Consequently the Oxley Creek Wetland, a pilot scale SSF wetland designed to enable rigorous monitoring, has recently been constructed in Brisbane, Australia. Tracer studies have been carried out in order to determine the hydraulics of this wetland prior to commissioning it with sealed sewage. The tracer studies will continue during the wetland's commissioning and operational phases. These studies will improve our understanding of the hydraulics of newly built SSF wetlands and the changes brought on by operational factors such as biological films and wetland plant root structures. Results to date indicate that the flow through the gravel beds is not uniform and cannot be adequately modelled by a single parameter, plug flow with dispersion, model. We have developed a multiparameter model, incorporating four plug flow reactors, which provides a better approximation of our experimental data. With further development this model will allow improvements to current SSF wetland design procedures and operational strategies, and will underpin investigations into the pollutant removal mechanisms at the Oxley Creek Wetland. (C) 1997 IAWQ. Published by Elsevier Science Ltd.

Vertical fluxes of sediment in oscillatory sheet flow

Time series of vertical sediment fluxes are derived from concentration time series in sheet flow under waves. While the concentrations C(z,t) vary very little with time for \z\ < 10d(50), the measured vertical sediment fluxes Q(zs)(z,t) vary strongly with time in this vertical band and their time variation follows, to some extent, the variation of the grain roughness Shields parameter 02,5(t). Thus, sediment distribution models based on the pickup function boundary condition are in some qualitative agreement with the measurements. However, the pickup function models are only able to model the upward bursts of sediment during the accelerating phases of the flow. They are, so far, unable to model the following strong downward sediment fluxes, which are observed during the periods of flow deceleration. Classical pickup functions, which essentially depend on the Shields parameter, are also incapable of modelling the secondary entrainment fluxes, which sometimes occur at free stream velocity reversal. The measured vertical fluxes indicate that the effective sediment settling velocity in the high [(0.3 < C(z,t) < 0.4] concentration area is typically only a few percent of the clear water settling velocity, while the measurements of Richardson and Jeronimo [Chem. Eng. Sci. 34 (1979) 1419], from a different physical setting, lead to estimates of the order 20%. The data does not support gradient diffusion as a model for sediment entrainment from the bed. That is, detailed modelling of the observed near-bed fluxes would require diffusivities that go negative during periods of flow deceleration. An observed general trend for concentration variability to increase with elevation close to the bed is also irreconcilable with diffusion models driven by a bottom boundary condition. (C) 2002 Published by Elsevier Science B.V.

Investigation of the flow field of a highly heated jet of air

Measurements of mean and fluctuating velocity and temperature and their self- and cross-products to the third-order are presented for a heated axisymmetric air jet. Froude numbers in the range of 3500 13,190, Reynolds numbers in the range of 3470-8500 and non-dimensional streamwise distances. X*, from 0.27 to 1.98 are covered by the data. It was found that turbulence intensity decreases for the heated jet in the region between the inertia dominated and the buoyancy dominated regions which is contrary to findings with helium jets mixing with ambient air to produce density fluctuations. The effects of heating on the turbulent kinetic energy budget and the temperature variance budget show small differences for the inertia dominated region and the intermediate region which help to explain the transition process to the far field plume region. Constants are evaluated for the isotropic eddy diffusivity and generalised gradient hypothesis models as well as the scalar variance model. No significant effect of heating on the dissipation time-scale ratio was found. A novel wire array with an inclined cold wire was used. Measurements obtained with this probe are found to lead to asymmetries in some of the higher-order products. Further investigation suggested that the asymmetries are attributable to an as yet unreported interference effect produced by the leading prong of the inclined temperature wire, The effect may also have implications for inclined velocity wires which contain a temperature component when used in heated flows. (C) 2002 Elsevier Science Inc. All rights reserved.

Fast left ventricle tracking in CMR images using localized anatomical affine optical flow

In daily cardiology practice, assessment of left ventricular (LV) global function using non-invasive imaging remains central for the diagnosis and follow-up of patients with cardiovascular diseases. Despite the different methodologies currently accessible for LV segmentation in cardiac magnetic resonance (CMR) images, a fast and complete LV delineation is still limitedly available for routine use. In this study, a localized anatomically constrained affine optical flow method is proposed for fast and automatic LV tracking throughout the full cardiac cycle in short-axis CMR images. Starting from an automatically delineated LV in the end-diastolic frame, the endocardial and epicardial boundaries are propagated by estimating the motion between adjacent cardiac phases using optical flow. In order to reduce the computational burden, the motion is only estimated in an anatomical region of interest around the tracked boundaries and subsequently integrated into a local affine motion model. Such localized estimation enables to capture complex motion patterns, while still being spatially consistent. The method was validated on 45 CMR datasets taken from the 2009 MICCAI LV segmentation challenge. The proposed approach proved to be robust and efficient, with an average distance error of 2.1 mm and a correlation with reference ejection fraction of 0.98 (1.9 ± 4.5%). Moreover, it showed to be fast, taking 5 seconds for the tracking of a full 4D dataset (30 ms per image). Overall, a novel fast, robust and accurate LV tracking methodology was proposed, enabling accurate assessment of relevant global function cardiac indices, such as volumes and ejection fraction.

Simulated annealing metaheuristic to solve the optimal power flow

The optimal power flow problem has been widely studied in order to improve power systems operation and planning. For real power systems, the problem is formulated as a non-linear and as a large combinatorial problem. The first approaches used to solve this problem were based on mathematical methods which required huge computational efforts. Lately, artificial intelligence techniques, such as metaheuristics based on biological processes, were adopted. Metaheuristics require lower computational resources, which is a clear advantage for addressing the problem in large power systems. This paper proposes a methodology to solve optimal power flow on economic dispatch context using a Simulated Annealing algorithm inspired on the cooling temperature process seen in metallurgy. The main contribution of the proposed method is the specific neighborhood generation according to the optimal power flow problem characteristics. The proposed methodology has been tested with IEEE 6 bus and 30 bus networks. The obtained results are compared with other wellknown methodologies presented in the literature, showing the effectiveness of the proposed method.

Linear and sliding-mode control design for matrix converter-based unified power flow controllers

This paper presents the design and compares the performance of linear, decoupled and direct power controllers (DPC) for three-phase matrix converters operating as unified power flow controllers (UPFC). A simplified steady-state model of the matrix converter-based UPFC fitted with a modified Venturini high-frequency pulse width modulator is first used to design the linear controllers for the transmission line active (P) and reactive (Q) powers. In order to minimize the resulting cross coupling between P and Q power controllers, decoupled linear controllers (DLC) are synthesized using inverse dynamics linearization. DPC are then developed using sliding-mode control techniques, in order to guarantee both robustness and decoupled control. The designed P and Q power controllers are compared using simulations and experimental results. Linear controllers show acceptable steady-state behaviour but still exhibit coupling between P and Q powers in transient operation. DLC are free from cross coupling but are parameter sensitive. Results obtained by DPC show decoupled power control with zero error tracking and faster responses with no overshoot and no steady-state error. All the designed controllers were implemented using the same digital signal processing hardware.

Evidence for non-coaxiality of ferrimagnetic and paramagnetic fabrics, developed during magma flow and cooling in a thick mafic dyke

A detailed analysis of fabrics of the chilled margin of a thick dolerite dyke (Foum Zguid dyke, Southern Morocco) was performed in order to better understand the development of sub-fabrics during dyke emplacement and cooling. AMS data were complemented with measurements of paramagnetic and ferrimagnetic fabrics (measured with high field torque magnetometer), neutron texture and microstructural analyses. The ferrimagnetic and AMS fabrics are similar, indicating that the ferrimagnetic minerals dominate the AMS signal. The paramagnetic fabric is different from the previous ones. Based on the crystallization timing of the different mineralogical phases, the paramagnetic fabric appears related to the upward flow, while the ferrimagnetic fabric rather reflects the late-stage of dyke emplacement and cooling stresses. (C) 2014 Elsevier B.V. All rights reserved.

Matrix converter as unified power flow controller: design and implementation of decoupled direct power controllers

This paper presents the design and implementation of direct power controllers for three-phase matrix converters (MC) operating as Unified Power Flow Controllers (UPFC). Theoretical principles of the decoupled linear power controllers of the MC-UPFC to minimize the cross-coupling between active and reactive power control are established. From the matrix converter based UPFC model with a modified Venturini high frequency PWM modulator, decoupled controllers for the transmission line active (P) and reactive (Q) power direct control are synthesized. Simulation results, obtained from Matlab/Simulink, are presented in order to confirm the proposed approach. Results obtained show decoupled power control, zero error tracking, and fast responses with no overshoot and no steady-state error.

Information flow analysis using data-dependent logical propositions

Dissertação para obtenção do Grau de Mestre em Engenharia Informática

A method to correct the flow distortion of offshore wind data using CFD simulation and experimental wind tunnel tests

The assessment of wind energy resource for the development of deep offshore wind plants requires the use of every possible source of data and, in many cases, includes data gathered at meteorological stations installed at islands, islets or even oil platforms—all structures that interfere with, and change, the flow characteristics. This work aims to contribute to the evaluation of such changes in the flow by developing a correction methodology and applying it to the case of Berlenga island, Portugal. The study is performed using computational fluid dynamic simulations (CFD) validated by wind tunnel tests. In order to simulate the incoming offshore flow with CFD models a wind profile, unknown a priori, was established using observations from two coastal wind stations and a power law wind profile was fitted to the existing data (a=0.165). The results show that the resulting horizontal wind speed at 80 m above sea level is 16% lower than the wind speed at 80 m above the island for the dominant wind direction sector.

Is intra-operative blood flow predictive for early failure of radiocephalic arteriovenous fistula?

BACKGROUND: For over 50 years, radiocephalic wrist arteriovenous fistulae (RCAVF) have been the primary and best vascular access for haemodialysis. Nevertheless, early failure due to thrombosis or non-maturation is a major complication resulting in their abandonment. This prospective study was designed to investigate the predictive value of intra-operative blood flow on early failure of primary RCAVF before the first effective dialysis. METHODS: We enrolled patients undergoing creation of primary RCAVF for haemodialysis based on the pre-operative ultrasound vascular mapping discussed in a multidisciplinary approach. Intra-operative blood flow measurement was systematically performed once the anastomosis had been completed using a transit-time ultrasonic flowmeter. During the follow-up, blood flow was estimated by colour flow ultrasound at various intervals. Any events related to the RCAVF were recorded. RESULTS: Autogenous RCAVFs (n = 58) in 58 patients were constructed and followed up for an average of 30 days. Thrombosis and non-maturation occurred in eight (14%) and four (7%) patients, respectively. The intra-operative blood flow in functioning RCAVFs was significantly higher compared to non-functioning RCAVFs (230 vs 98 mL/min; P = 0.007), as well as 1 week (753 vs 228 mL/min; P = 0.0008) and 4 weeks (915 vs 245 mL/min, P < 0.0001) later. Blood flow volume measurements with a cut-off value of 120 mL/min had a sensitivity of 67%, specificity of 75% and positive predictive value of 91%. CONCLUSIONS: Blood flow <120 mL has a good predictive value for early failure in RCAVF. During the procedure, this cut-off value may be used to select appropriately which RCAVF should be investigated in the operation theatre in order to correct in real time any abnormality.