Selecting methods to solve multi-well master equations

There are several competing methods commonly used to solve energy grained master equations describing gas-phase reactive systems. When it comes to selecting an appropriate method for any particular problem, there is little guidance in the literature. In this paper we directly compare several variants of spectral and numerical integration methods from the point of view of computer time required to calculate the solution and the range of temperature and pressure conditions under which the methods are successful. The test case used in the comparison is an important reaction in combustion chemistry and incorporates reversible and irreversible bimolecular reaction steps as well as isomerizations between multiple unimolecular species. While the numerical integration of the ODE with a stiff ODE integrator is not the fastest method overall, it is the fastest method applicable to all conditions.

Comportamento dinâmico da rede eléctrica da ilha de São Vicente em Cabo Verde: estabilidade de tensão

Neste trabalho é efectuado, não só o diagnóstico em regime permanente, mas também o estudo, simulação e análise do comportamento dinâmico da rede eléctrica da ilha de São Vicente em Cabo Verde. Os estudos de estabilidade transitória desempenham um importante papel, tanto no planeamento como na operação dos sistemas de potência. Tais estudos são realizados, em grande parte, através de simulação digital no domínio do tempo, utilizando integração numérica para resolver as equações não-lineares que modelam a dinâmica do sistema e dependem da existência de registos reais de perturbação (ex: osciloperturbografia). O objectivo do trabalho será também verificar a aplicabilidade dos requisitos técnicos que as unidades geradoras devem ter, no que concerne ao controlo de tensão, estabelecidos na futura regulamentação europeia desenvolvida pela ENTSO-E (European Network Transmission System Operator for Electricity). De entre os requisitos analisou-se a capacidade das máquinas existentes suportarem cavas de tensão decorrentes de curto-circuitos trifásicos simétricos, Fault Ride Through, no ponto de ligação à rede. Identificaram-se para o efeito os factores que influenciam a estabilidade desta rede, em regime perturbado nomeadamente: (i) duração do defeito, (ii) caracterização da carga, com e sem a presença do sistema de controlo de tensão (AVR) em unidades de geração síncronas. Na ausência de registos reais sobre o comportamento do sistema, conclui-se que este é sensível à elasticidade das cargas em particular do tipo potência constante, existindo risco de perda de estabilidade, neste caso, para defeitos superiores a 5ms sem AVR. A existência de AVR nesta rede afigura-se como indispensável para garantir estabilidade de tensão sendo contudo necessário proceder a uma correcta parametrização.

Sistema de multiprocessamento para simulação de N-corpos em FPGA

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia de Electrónica e Telecomunicações

Identificação modal e simulação da resposta sísmica de barragens abóbada

Trabalho Final de Mestrado elaborado no Laboratório Nacional de Engenharia Civil (LNEC) para a obtenção do grau de Mestre em Engenharia Civil pelo Instituto Superior de Engenharia de Lisboa no âmbito do protocolo de cooperação ente o ISEL e o LNEC

Domain growth in binary mixtures at low temperatures

We have studied domain growth during spinodal decomposition at low temperatures. We have performed a numerical integration of the deterministic time-dependent Ginzburg-Landau equation with a variable, concentration-dependent diffusion coefficient. The form of the pair-correlation function and the structure function are independent of temperature but the dynamics is slower at low temperature. A crossover between interfacial diffusion and bulk diffusion mechanisms is observed in the behavior of the characteristic domain size. This effect is explained theoretically in terms of an equation of motion for the interface.

Systematic weakly nonlinear analysis of radial viscous fingering

We present a weakly nonlinear analysis of the interface dynamics in a radial Hele-Shaw cell driven by both injection and rotation. We extend the systematic expansion introduced in [E. Alvarez-Lacalle et al., Phys. Rev. E 64, 016302 (2001)] to the radial geometry, and compute explicitly the first nonlinear contributions. We also find the necessary and sufficient condition for the uniform convergence of the nonlinear expansion. Within this region of convergence, the analytical predictions at low orders are compared satisfactorily to exact solutions and numerical integration of the problem. This is particularly remarkable in configurations (with no counterpart in the channel geometry) for which the interplay between injection and rotation allows that condition to be satisfied at all times. In the case of the purely centrifugal forcing we demonstrate that nonlinear couplings make the interface more unstable for lower viscosity contrast between the fluids.

Interfacial instabilities of a fluid annulus in a rotating Hele-Shaw cell

We have studied the interfacial instabilities experienced by a liquid annulus as it moves radially in a circular Hele-Shaw cell rotating with angular velocity Omega. The instability of the leading interface (oil displacing air) is driven by the density difference in the presence of centrifugal forcing, while the instability of the trailing interface (air displacing oil) is driven by the large viscosity contrast. A linear stability analysis shows that the stability of the two interfaces is coupled through the pressure field already at a linear level. We have performed experiments in a dry cell and in a cell coated with a thin fluid layer on each plate, and found that the stability depends substantially on the wetting conditions at the leading interface. Our experimental results of the number of fingers resulting from the instability compare well with the predictions obtained through a numerical integration of the coupled equations derived from a linear stability analysis. Deep in the nonlinear regime we observe the emission of liquid droplets through the formation of thin filaments at the tip of outgrowing fingers.

Nonlinear relaxation time and the detection of weak signals

Laser systems can be used to detect very weak optical signals. The physical mechanism is the dynamical process of the relaxation of a laser from an unstable state to a steady stable state. We present an analysis of this process based on the study of the nonlinear relaxation time. Our analytical results are compared with numerical integration of the stochastic differential equations that model this process.

Bistability driven by dichotomous noise

We consider mean-first-passage times and transition rates in bistable systems driven by dichotomous colored noise. We carry out an asymptotic expansion for short correlation times ¿c of the colored noise and find results that differ from those reported earlier. In particular, to retain corrections to O(¿c) we find that it is necessary to retain up to four derivatives of the potential function. We compare our asymptotic results to existing ones and also to exact ones obtained from numerical integration.

Equidistribution of Fekete Points on the Sphere

Fekete points are the points that maximize a Vandermonde-type determinant that appears in the polynomial Lagrange interpolation formula. They are well suited points for interpolation formulas and numerical integration. We prove the asymptotic equidistribution of Fekete points in the sphere. The way we proceed is by showing their connection to other arrays of points, the so-called Marcinkiewicz-Zygmund arrays and interpolating arrays, that have been studied recently.

Equidistribution of Fekete Points on the Sphere

Fekete points are the points that maximize a Vandermonde-type determinant that appears in the polynomial Lagrange interpolation formula. They are well suited points for interpolation formulas and numerical integration. We prove the asymptotic equidistribution of Fekete points in the sphere. The way we proceed is by showing their connection to other arrays of points, the so-called Marcinkiewicz-Zygmund arrays and interpolating arrays, that have been studied recently.

CR1Dmod: A Matlab program to model 1D complex resistivity effects in electrical and electromagnetic surveys

We have constructed a forward modelling code in Matlab, capable of handling several commonly used electrical and electromagnetic methods in a 1D environment. We review the implemented electromagnetic field equations for grounded wires, frequency and transient soundings and present new solutions in the case of a non-magnetic first layer. The CR1Dmod code evaluates the Hankel transforms occurring in the field equations using either the Fast Hankel Transform based on digital filter theory, or a numerical integration scheme applied between the zeros of the Bessel function. A graphical user interface allows easy construction of 1D models and control of the parameters. Modelling results are in agreement with other authors, but the time of computation is less efficient than other available codes. Nevertheless, the CR1Dmod routine handles complex resistivities and offers solutions based on the full EM-equations as well as the quasi-static approximation. Thus, modelling of effects based on changes in the magnetic permeability and the permittivity is also possible.

Sähkökäyttösimulaattori jasen suorituskyky

Erilaisten simulaatioiden tekeminen tutkimustyössä on tärkeää. Simulaatioiden avulla voidaan vähentää prototyyppitestauksen tarvetta. Diplomityössä on esitelty kehitettävää sähkökäyttösimulaattoria, jolla voidaan tarkastella erilaisten sähkökäyttöjen häviöiden muodostumista. Diplomityössä on keskitytty vertailemaan kehitettävän simulaattorin simulointituloksia todelliselta sähkökäytöltä mitattuihin suureisiin. Vertailun kohteena on taajuusmuuttajalla syötetty oikosulkumoottori, minkä virtojen ja jännitteiden vertailu on tehty aika- ja taajuustasossa.

An Assessment ofFatigue Strength of Load Carrying Welded Joints with Lack of Penetration

On yleisesti tiedossa, että väsyttävän kuormituksen alaisena olevat hitsatut rakenteet rikkoutuvat juuri hitsausliitoksista. Täyden tunkeuman hitsausliitoksia sisältävien rakenteiden asiantunteva suunnittelu janykyaikaiset valmistusmenetelmät ovat lähes eliminoineet väsymisvauriot hitsatuissa rakenteissa. Väsymislujuuden parantaminen tiukalla täyden tunkeuman vaatimuksella on kuitenkin epätaloudellinen ratkaisu. Täyden tunkeuman hitsausliitoksille asetettavien laatuvaatimuksien on määriteltävä selkeät tarkastusohjeet ja hylkäämisperusteet. Tämän diplomityön tarkoituksena oli tutkia geometristen muuttujien vaikutusta kuormaa kantavien hitsausliitosten väsymislujuuteen. Huomio kiinnitettiin pääasiassa suunnittelumuuttujiin, joilla on vaikutusta väsymisvaurioiden syntymiseen hitsauksen juuren puolella. Nykyiset määräykset ja standardit, jotka perustuvat kokeellisiin tuloksiin; antavat melko yleisiä ohjeita hitsausliitosten väsymismitoituksesta. Tämän vuoksi muodostettiin kokonaan uudet parametriset yhtälöt sallitun nimellisen jännityksen kynnysarvon vaihteluvälin, ¿¿th, laskemiseksi, jotta vältettäisiin hitsausliitosten juuren puoleiset väsymisvauriot. Lisäksi, jokaiselle liitostyypille laskettiin hitsin juuren puolen väsymisluokat (FAT), joita verrattiin olemassa olevilla mitoitusohjeilla saavutettuihin tuloksiin. Täydentäviksi referensseiksi suoritettiin useita kolmiulotteisia (3D) analyysejä. Julkaistuja kokeellisiin tuloksiin perustuvia tietoja käytettiin apuna hitsausliitosten väsymiskäyttäytymisen ymmärtämiseksi ja materiaalivakioiden määrittämiseksi. Kuormaa kantavien vajaatunkeumaisten hitsausliitosten väsymislujuus määritettiin käyttämällä elementtimenetelmää. Suurimman pääjännityksen kriteeriä hyödynnettiin murtumiskäyttäytymisen ennakoimiseksi. Valitulle hitsatulle materiaalille ja koeolosuhteille murtumiskäyttäytymistä mallinnettiin särön kasvunopeudella da/dN ja jännitysintensiteettikertoimen vaihteluvälillä, 'K. Paris:n yhtälön numeerinen integrointi suoritettiin FRANC2D/L tietokoneohjelmalla. Saatujen tulosten perusteella voidaan laskea FAT tutkittavassa tapauksessa. ¿¿th laskettiin alkusärön jännitysintensiteettikertoimen vaihteluvälin ja kynnysjännitysintensiteettikertoimen, 'Kth, perusteella. ¿Kth arvoa pienemmällä vaihteluvälillä särö ei kasva. Analyyseissäoletuksena oli hitsattu jälkikäsittelemätön liitos, jossa oli valmis alkusärö hitsin juuressa. Analyysien tulokset ovat hyödyllisiä suunnittelijoille, jotka tekevät päätöksiä koskien geometrisiä parametreja, joilla on vaikutusta hitsausliitosten väsymislujuuteen.

A computationally efficient shear deformable beam element for large deformation multibody applications

In this paper, a new two-dimensional shear deformable beam element based on the absolute nodal coordinate formulation is proposed. The nonlinear elastic forces of the beam element are obtained using a continuum mechanics approach without employing a local element coordinate system. In this study, linear polynomials are used to interpolate both the transverse and longitudinal components of the displacement. This is different from other absolute nodal-coordinate-based beam elements where cubic polynomials are used in the longitudinal direction. The accompanying defects of the phenomenon known as shear locking are avoided through the adoption of selective integration within the numerical integration method. The proposed element is verified using several numerical examples, and the results are compared to analytical solutions and the results for an existing shear deformable beam element. It is shown that by using the proposed element, accurate linear and nonlinear static deformations, as well as realistic dynamic behavior, can be achieved with a smaller computational effort than by using existing shear deformable two-dimensional beam elements.