Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

The subcycled Newmark algorithm

The popular Newmark algorithm, used for implicit direct integration of structural dynamics, is extended by means of a nodal partition to permit use of different timesteps in different regions of a structural model. The algorithm developed has as a special case an explicit-explicit subcycling algorithm previously reported by Belytschko, Yen and Mullen. That algorithm has been shown, in the absence of damping or other energy dissipation, to exhibit instability over narrow timestep ranges that become narrower as the number of degrees of freedom increases, making them unlikely to be encountered in practice. The present algorithm avoids such instabilities in the case of a one to two timestep ratio (two subcycles), achieving unconditional stability in an exponential sense for a linear problem. However, with three or more subcycles, the trapezoidal rule exhibits stability that becomes conditional, falling towards that of the central difference method as the number of subcycles increases. Instabilities over narrow timestep ranges, that become narrower as the model size increases, also appear with three or more subcycles. However by moving the partition between timesteps one row of elements into the region suitable for integration with the larger timestep these the unstable timestep ranges become extremely narrow, even in simple systems with a few degrees of freedom. As well, accuracy is improved. Use of a version of the Newmark algorithm that dissipates high frequencies minimises or eliminates these narrow bands of instability. Viscous damping is also shown to remove these instabilities, at the expense of having more effect on the low frequency response.

Explicit/implicit partitioning and a new explicit form of the generalized alpha method

Most finite element packages use the Newmark algorithm for time integration of structural dynamics. Various algorithms have been proposed to better optimize the high frequency dissipation of this algorithm. Hulbert and Chung proposed both implicit and explicit forms of the generalized alpha method. The algorithms optimize high frequency dissipation effectively, and despite recent work on algorithms that possess momentum conserving/energy dissipative properties in a non-linear context, the generalized alpha method remains an efficient way to solve many problems, especially with adaptive timestep control. However, the implicit and explicit algorithms use incompatible parameter sets and cannot be used together in a spatial partition, whereas this can be done for the Newmark algorithm, as Hughes and Liu demonstrated, and for the HHT-alpha algorithm developed from it. The present paper shows that the explicit generalized alpha method can be rewritten so that it becomes compatible with the implicit form. All four algorithmic parameters can be matched between the explicit and implicit forms. An element interface between implicit and explicit partitions can then be used, analogous to that devised by Hughes and Liu to extend the Newmark method. The stability of the explicit/implicit algorithm is examined in a linear context and found to exceed that of the explicit partition. The element partition is significantly less dissipative of intermediate frequencies than one using the HHT-alpha method. The explicit algorithm can also be rewritten so that the discrete equation of motion evaluates forces from displacements and velocities found at the predicted mid-point of a cycle. Copyright (C) 2003 John Wiley Sons, Ltd.

Comparative study of power converter topologies and control strategies for the harmonic performance of variable-speed wind turbine generator systems

Power converters play a vital role in the integration of wind power into the electrical grid. Variable-speed wind turbine generator systems have a considerable interest of application for grid connection at constant frequency. In this paper, comprehensive simulation studies are carried out with three power converter topologies: matrix, two-level and multilevel. A fractional-order control strategy is studied for the variable-speed operation of wind turbine generator systems. The studies are in order to compare power converter topologies and control strategies. The studies reveal that the multilevel converter and the proposed fractional-order control strategy enable an improvement in the power quality, in comparison with the other power converters using a classical integer-order control strategy. (C) 2010 Elsevier Ltd. All rights reserved.

Ground Vibrations Induced by High-Speed Trains in Regions with Sudden Foundation Stiffness Change

In this paper analytical transient solutions of dynamic response of one-dimensional systems with sudden change of foundation stiffness are derived. In more details, cantilever dynamic response, expressed in terms of vertical displacement, is extended to account for elastic foundation and then two cantilever solutions, corresponding to beams clamped on left and right hand side, with different value of Winkler constant are connected together by continuity conditions. The internal forces, as the unknowns, can be introduced by the same values in both clamped beam solutions and solved. Assumption about time variation of internal forces at the section of discontinuity must be adopted and originally analytical solution will have to include numerical procedure.

Dynamic compact thermal models with multiple power sources: application to an ultrathin chip stacking technology

Whereas numerical modeling using finite-element methods (FEM) can provide transient temperature distribution in the component with enough accuracy, it is of the most importance the development of compact dynamic thermal models that can be used for electrothermal simulation. While in most cases single power sources are considered, here we focus on the simultaneous presence of multiple sources. The thermal model will be in the form of a thermal impedance matrix containing the thermal impedance transfer functions between two arbitrary ports. Eachindividual transfer function element ( ) is obtained from the analysis of the thermal temperature transient at node ¿ ¿ after a power step at node ¿ .¿ Different options for multiexponential transient analysis are detailed and compared. Among the options explored, small thermal models can be obtained by constrained nonlinear least squares (NLSQ) methods if the order is selected properly using validation signals. The methods are applied to the extraction of dynamic compact thermal models for a new ultrathin chip stack technology (UTCS).

Häiriöanalyysejä tehonkorotetun kiehutusvesivoimalaitoksen käytettävyyden varmistamiseksi

TVO suunnittelee reaktoritehon 10 %:n korotusta Olkiluoto 1 ja 2 -voimalaitoksille. Reaktoriteho nostetaan 2500 MW:sta 2750 MW:iin polttoaineen rikastusastetta nostamalla ja pääkiertovirtausta kasvattamalla. Samalla syöttövesivirtaus reaktoriin ja tuorehöyryvirtaus turpiineille kasvaa. Lauhteenpuhdistusjärjestelmän kapasiteettia ei voida kuitenkaan kasvattaa, joten massavirran lisäys toteutetaan ottamalla käyttöön korkeapainesivulauhteen eteenpäinpumppaus. Lauhteen esilämmityslinjojen, lauhteenpuhdistuksen ja syöttövesipumppujen massavirta säilyy siten nykyisellään. Muita merkittäviä tehonkorotukseen liittyviä laitosmuutoksia ovat pääkiertopumppujen uusinta ja korkeapaineturpiinin muutokset. Tehonkorotetun prosessin käytettävyyden varmistamiseksi tehdään häiriöanalyysejä Apros-prosessisimulointiohjelmistoa käyttäen. OL1 ja OL2 -laitoksista on olemassa validoitu 2500 MW:n laitosmalli, josta muokatulla 2750 MW:n laitosmallilla simuloinnit tehdään. Häiriöanalyysien avulla selvitetään säätöjärjestelmien kyky pitää prosessin tila hallinnassa ilman suojausautomaation laukeamista. Simuloituihin tapauksiin kuuluu pumppujen ja venttiilien vikaantumistapauksia sekä turpiini- ja reaktoripuolen pikasulku- ja osittaispikasulkutapauksia. Myös meriveden lämpötilan vaikutusta häiriötilanteisiin tarkastellaan. Analyysien perusteella voimalaitosten ohjaus- ja suojausautomaatio toimivat hyvin myös korotetulla teholla. Tehonkorotuksen jälkeiset suuremmat massavirrat aiheuttavat kuitenkin voimakkaampia reaktoripaineen ja -tehon vaihteluita varsinkin venttiilien sulkeutumistapauksissa. Simuloinnit osoittivat, että tehonkorotus 2500 MW:sta 2750 MW:iin on mahdollinen, mutta aiheuttaa pieniä muutoksia laitoksen suojausjärjestelmien laukaisurajoihin.

Investigations on broadband planar monopole and slot radiators and their suitability for UWB applications

In recent years, there is a visible trend for products/services which demand seamless integration of cellular networks, WLANs and WPANs. This is a strong indication for the inclusion of high speed short range wireless technology in future applications. In this context UWB radio has a significant role to play as an extension/complement to existing cellular/access technology. In the present work, we have investigated two major types of wide band planar antennas: Monopole and Slot. Four novel compact broadband antennas, suitable for poratble applications, are designed and characterized, namely 1. Elliptical monopole 2. Inverted cone monopole 3. Koch fractal slot 4. Wide band slot The performance of these designs have been studied using standard simulation tools used in industry/academia and they have been experimentally verified. Antenna design guidelines are also deduced by accounting the resonances in each structure. In addition to having compact sized, high efficiency and broad bandwidth antennas, one of the major criterion in the design of impulse-UWB systems have been the transmission of narrow band pulses with minimum distortion. The key challenge is not only to design a broad band antenna with constant and stable gain but to maintain a flat group delay or linear phase response in the frequency domain or excellent transient response in time domain. One of the major contributions of the thesis lies in the analysis of the frequency and time-domain response of the designed UWB antennas to confirm their suitability for portable pulsed-UWB systems. Techniques to avoid narrowband interference by engraving narrow slot resonators on the antenna is also proposed and their effect on a nano-second pulse have been investigated.

Ultra-wideband slot antenna with band-notch characteristics for wireless USB dongle applications

A compact ultra-wideband (UWB) printed slot antenna is described, suitable for integration with the printed circuit board (PCB) of a wireless, universal, serial-bus dongle. The design comprises of a near-rectangular slot fed by a coplanar waveguide (CPW) printed on a PCB of size 20 × 30 mm2. It has a large bandwidth covering the 3.1–10.6 GHz UWB band, with omnidirectional radiation patterns. Further, a notched band centered at 5.45 GHz wireless local area network bands is obtained within the wide bandwidth by inserting a narrow slot inside the tuning stub. Details of the antenna design are described, and the experimental results of the constructed prototype are presented. The time domain studies on the antenna shows a linear phase response throughout the band except at the notched frequency. The transient analysis of the antenna indicates very little pulse distortion confirming its suitability for high speed wireless connectivity.

Dedução de Funções de Alívio para Priorização de Variantes de Chaveamento Adequadas à Eliminação de Sobrecargas

Eventually, violations of voltage limits at buses or admissible loadings of transmission lines and/or power transformers may occur by the power system operation. If violations are detected in the supervision process, corrective measures may be carried out in order to eliminate them or to reduce their intensity. Loading restriction is an extreme solution and should only be adopted as the last control action. Previous researches have shown that it is possible to control constraints in electrical systems by changing the network topology, using the technique named Corrective Switching, which requires no additional costs. In previous works, the proposed calculations for verifying the ability of a switching variant in eliminating an overload in a specific branch were based on network reduction or heuristic analysis. The purpose of this work is to develop analytical derivation of linear equations to estimate current changes in a specific branch (due to switching measures) by means of few calculations. For bus-bar coupling, derivations will be based on short-circuit theory and Relief Function methodology. For bus-bar splitting, a Relief Function will be derived based on a technique of equivalent circuit. Although systems of linear equations are used to substantiate deductions, its formal solution for each variant, in real time does not become necessary. A priority list of promising variants is then assigned for final check by an exact load flow calculation and a transient analysis using ATP Alternative Transient Program. At last, results obtained by simulation in networks with different features will be presented

Avaliação da aplicação de modelos para a secagem de pastas e suspensões em leito de jorro com alimentação contínua

This work studies the development, implementation and improvement of a macroscopic model to describe the behavior of the spouted bed dryer with continuous feeding for pastes and suspensions drying. This model is based on the CST model (Freire et al., 2009) and the model of Fernandes (2005), whose theoretical foundation is based on macroscopic mass and heat balances for the three phases involved in the process: gas, liquid and solid. Because this technique is quite relevant, the studies of modeling and simulation of spouted bed drying are essential in the analysis of the process as a whole, because through them it is possible to predict and understand the behavior of the process, which contributes significantly to more efficient project and operation. The development and understanding of the phenomena involved in the drying process can be obtained by comparing the experimental data with those from computer simulations. Such knowledge is critical for choosing properly the process conditions in order to obtain a good drying efficiency. Over the past few years, researches and development of works in the field of pastes and suspensions drying in spouted bed has been gaining ground in Brazil. The Particulate Systems Laboratory at Universidade Federal do Rio Grande do Norte, has been developing several researches and generating a huge collection of experimental data concerning the drying of fruit pulps, vegetables pastes, goat milk and suspensions of agro-industrial residues. From this collection, some data of goat milk and residue from acerola (Malpighia glabra L.) drying were collected. For the first time, these data were used for the development and validation of a model that can describe the behavior of spouted bed dryer. Thus, it was possible to model the dryer and to evaluate the influence of process variables (paste feeding, temperature and flow rate of the drying air) in the drying dynamics. We also performed water evaporation experiments in order to understand and to study the behavior of the dryer wall temperature and the evaporation rate. All these analysis will contribute to future works involving the implementation of control strategies in the pastes and suspensions drying. The results obtained in transient analysis were compared with experimental data indicating that this model well represents the process

Modal representation of three-phase lines applying two transformation matrices: Evaluation of its eigenvectors

The objective of this paper is to show an alternative representation in time domain of a non-transposed three-phase transmission line decomposed in its exact modes by using two transformation matrices. The first matrix is Clarke's matrix that is real, frequency independent, easily represented in computational transient programs (EMTP) and separates the line into Quasi-modes α, β and zero. After that, Quasi-modes a and zero are decomposed into their exact modes by using a modal transformation matrix whose elements can be synthesized in time domain through standard curve-fitting techniques. The main advantage of this alternative representation is to reduce the processing time because a frequency dependent modal transformation matrix of a three-phase line has nine elements to be represented in time domain while a modal transformation matrix of a two-phase line has only four elements. This paper shows modal decomposition process and eigenvectors of a non-transposed three-phase line with a vertical symmetry plane whose nominal voltage is 440 kV and line length is 500 km. ©2006 IEEE.

An alternative procedure to decrease the dimension of the frequency dependent modal transformation matrices: Application in three-phase transmission lines with a vertical symmetry plane

The objective of this paper is to show an alternative representation in time domain of a non-transposed three-phase transmission line decomposed in its exact modes by using two transformation matrices. The first matrix is Clarke's matrix that is real, frequency independent, easily represented in computational transient programs (EMTP) and separates the line into Quasi-modes α, β and zero. After that, Quasi-modes α and zero are decomposed into their exact modes by using a modal transformation matrix whose elements can be synthesized in time domain through standard curve-fitting techniques. The main advantage of this alternative representation is to reduce the processing time because a frequency dependent modal transformation matrix of a three-phase line has nine elements to be represented in time domain while a modal transformation matrix of a two-phase line has only four elements. This paper shows modal decomposition process and eigenvectors of a non-transposed three-phase line with a vertical symmetry plane whose nominal voltage is 440 kV and line length is 500 km. © 2006 IEEE.

Análise numérica de um sistema integrado coletor solar/armazenador térmico

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Modeling, design and experimental characterization of Micro-Electro-Mechanical-Systems for gas chromatographic applications

Design parameters, process flows, electro-thermal-fluidic simulations and experimental characterizations of Micro-Electro-Mechanical-Systems (MEMS) suited for gas-chromatographic (GC) applications are presented and thoroughly described in this thesis, whose topic belongs to the research activities the Institute for Microelectronics and Microsystems (IMM)-Bologna is involved since several years, i.e. the development of micro-systems for chemical analysis, based on silicon micro-machining techniques and able to perform analysis of complex gaseous mixtures, especially in the field of environmental monitoring. In this regard, attention has been focused on the development of micro-fabricated devices to be employed in a portable mini-GC system for the analysis of aromatic Volatile Organic Compounds (VOC) like Benzene, Toluene, Ethyl-benzene and Xylene (BTEX), i.e. chemical compounds which can significantly affect environment and human health because of their demonstrated carcinogenicity (benzene) or toxicity (toluene, xylene) even at parts per billion (ppb) concentrations. The most significant results achieved through the laboratory functional characterization of the mini-GC system have been reported, together with in-field analysis results carried out in a station of the Bologna air monitoring network and compared with those provided by a commercial GC system. The development of more advanced prototypes of micro-fabricated devices specifically suited for FAST-GC have been also presented (silicon capillary columns, Ultra-Low-Power (ULP) Metal OXide (MOX) sensor, Thermal Conductivity Detector (TCD)), together with the technological processes for their fabrication. The experimentally demonstrated very high sensitivity of ULP-MOX sensors to VOCs, coupled with the extremely low power consumption, makes the developed ULP-MOX sensor the most performing metal oxide sensor reported up to now in literature, while preliminary test results proved that the developed silicon capillary columns are capable of performances comparable to those of the best fused silica capillary columns. Finally, the development and the validation of a coupled electro-thermal Finite Element Model suited for both steady-state and transient analysis of the micro-devices has been described, and subsequently implemented with a fluidic part to investigate devices behaviour in presence of a gas flowing with certain volumetric flow rates.