With the ever advancing state of computer systems, it is imperative to maintain the most up-to-date and reliable safety evaluation data for nuclear power systems. Commonplace now is the practice of updating old accident simulation results with more advanced models and codes using today's faster computer systems. Though it may be quite an undertaking, the benefits of using a more advanced model and code can be significant especially if the result of the new analysis provides increased safety margin for any plant component or system. A series of parametric and sensitivity studies for the Loss of Normal Feedwater Anticipated Transient without Scram (LONF ATWS) for Southern Company's Vogtle Electric Generating Plant (VEGP) Units 1&2 located near Waynesboro, GA was performed using the best-estimate thermal-hydraulics transient analysis code RETRAN-02w. This thesis includes comparison to the results of a generic plant study published by Westinghouse Electric Corporation in 1974 using an earlier code, LOFTRAN, as well as Vogtle-specific analysis. The comparative analysis exposes and seeks to explain differences between the two codes whereas the Vogtle analysis utilizes data from the Vogtle FSAR to generate plant-specific data. The purpose of this study is to validate and update the previous analysis and gather more information about the plant actions taken in response to a LONF ATWS. As a result, now there is a new and updated evaluation of the LONF ATWS for both a generic 4-loop Westinghouse plant and VEGP using a more advanced code. Beyond the reference case analysis, a series of sensitivity and parametric studies have been performed to show how well each type of plant is designed for handling an ATWS situation. These studies cover a wide range of operating conditions to demonstrate the dependability of the model. It was found that both the generic 4-loop Westinghouse PWR system and VEGP can successfully mitigate a LONF ATWS throughout the core's operating cycle.

<p>The Discontinuous Galerkin time domain (DGTD) method is promising in modeling of realistic multiscale electromagnetic systems. This method defines the basic concept for implementing the communication between multiple domains with different scales.</p><p>Constructing a DGTD system consists of several careful choices: (a) governing equations; (b) element shape and corresponding basis functions for the spatial discretization of each subdomain; (c) numerical fluxes onto interfaces to bond all subdomains together; and (d) time stepping scheme based on properties of a discretized</p><p>system. This work present the advances in each one of these steps.</p><p> </p><p>First, a unified framework based on the theory of differential forms and the finite element method is used to analyze the discretization of the Maxwell's equations. Based on this study, field intensities (<bold>E</bold> and <bold>H</bold>) are associated to 1-forms and curl-conforming basis functions; flux densities (<bold>D</bold> and <bold>B</bold>) are associated to 2-forms and divergence-conforming basis functions; and the constitutive relations are defined by Hodge operators.</p><p>A different approach is the study of numerical dispersion. Semidiscrete analysis is the traditional method, but for high order elements modal analysis is prefered. From these analyses, we conclude that a correct discretization of fields belonging to different p-form (e.g., <bold>E</bold> and <bold>B</bold>) uses basis functions with same order of interpolation; however, different order of interpolation must be used if two fields belong to the same p-form (e.g., <bold>E</bold> and <bold>H</bold>). An alternative method to evaluate numerical dispersion based on evaluation of dispersive Hodge operators is also presented. Both dispersion analyses are equivalent and reveal same fundamental results. Eigenvalues, eigenvector and transient results are studied to verify accuracy and computational costs of different schemes. </p><p>Two different approaches are used for implementing the DG Method. The first is based on <bold>E</bold> and <bold>H</bold> fields, which use curl-conforming basis functions with different order of interpolation. In this case, the Riemman solver shows the best performance to treat interfaces between subdomains. A new spectral prismatic element, useful for modeling of layer structures, is also implemented for this approach. Furthermore, a new efficient and very accurate time integration method for sequential subdomains is implemented.</p><p>The second approach for solving multidomain cases is based on <bold>E</bold> and <bold>B</bold> fields, which use curl- and divergence-conforming basis functions, respectively, with same order of interpolation. In this way, higher accuracy and lower memory consumption are obtained with respect to the first approach based on <bold>E</bold> and <bold>H</bold> fields. The centered flux is used to treat interfaces with non-conforming meshes, and both explicit Runge-Kutta method and implicit Crank-Nicholson method are implemented for time integration. </p><p>Numerical examples and realistic cases are presented to verify that the proposed methods are non-spurious and efficient DGTD schemes.</p> / Dissertation

Submitted by Jaqueline Silva (jtas29@gmail.com) on 2016-08-29T17:39:21Z No. of bitstreams: 3 Dissertação - Bárbara Pereira Silva - 2016 - parte 1.pdf: 12258499 bytes, checksum: 8af480355ae4ffedb4ed79cbf7484361 (MD5) Dissertação - Bárbara Pereira Silva - 2016 - parte 2.pdf: 12730012 bytes, checksum: d173d49860226f3a2179be40d51ef666 (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) / Approved for entry into archive by Jaqueline Silva (jtas29@gmail.com) on 2016-08-29T17:39:44Z (GMT) No. of bitstreams: 3 Dissertação - Bárbara Pereira Silva - 2016 - parte 1.pdf: 12258499 bytes, checksum: 8af480355ae4ffedb4ed79cbf7484361 (MD5) Dissertação - Bárbara Pereira Silva - 2016 - parte 2.pdf: 12730012 bytes, checksum: d173d49860226f3a2179be40d51ef666 (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) / Made available in DSpace on 2016-08-29T17:39:44Z (GMT). No. of bitstreams: 3 Dissertação - Bárbara Pereira Silva - 2016 - parte 1.pdf: 12258499 bytes, checksum: 8af480355ae4ffedb4ed79cbf7484361 (MD5) Dissertação - Bárbara Pereira Silva - 2016 - parte 2.pdf: 12730012 bytes, checksum: d173d49860226f3a2179be40d51ef666 (MD5) license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) Previous issue date: 2016-03-29 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES / In the grounding systems analysis, the investigation of the effects caused by lightning discharge is a crucial task, since it is the main electromagnetic phenomena responsible for electricity outages which can be the cause of severe mechanical and thermal stresses on electric installations. This dissertation presents a study of transient behavior for various grounding configurations subject to impulsive currents originated from incident discharges. The model is named hybrid electromagnetic model, because it is based on electromagnetic equations and it is used concepts of electrical circuit analysis. The modeling is performed in frequency domain considering the effects of electromagnetic propagation on the metal electrode and into surrounding soil. In the numerical solution of model is used the method of moments for the discretization of the electrodes and calculus of the impedances. This methodology also considers the variation of electric soil parameters when frequency varies. The modeling comprised the effects of propagation and electromagnetic couplings between the grounding electrodes. Results were obtained for simple electrodes (wire horizontally buried in the soil or a vertical rod) and grounding grid in square format from the computational implementation of the mathematical model. This tool and the results obtained lend itself to the determination of the grounding systems response subjected to an injection of a typical current of lightning discharge. The validation of the tool which resulted from this work is made by comparing the obtained results with results of similar cases reported in the technical literature. / Na análise de sistemas de aterramento elétrico é primordial a investigação dos efeitos causados pela incidência de descarga atmosférica, uma vez que esse é o principal fenômeno eletromagnético responsável pelos desligamentos não programados de instalações elétricas causando severas solicitações mecânicas e térmicas. Nesta dissertação é apresentado um estudo do comportamento transitório para diferentes configurações de eletrodos de aterramentos frente a descargas atmosféricas. O modelo adotado é nomeado de modelo eletromagnético híbrido, pois baseia-se nas equações eletromagnéticas e utiliza-se conceitos da análise de circuitos elétricos. A modelagem é realizada no domínio da frequência, considerando-se os efeitos da propagação eletromagnética no eletrodo metálico e no solo ao redor. Na solução numérica do modelo é empregado o método dos momentos para a discretizacão dos eletrodos e obtenção das impedâncias. É também levada em conta a variação dos parâmetros elétricos do solo com a frequência. São contemplados nesta modelagem os efeitos de propagação e os acoplamentos eletromagnéticos entre os eletrodos do aterramento. A partir da implementação computacional do modelo, foram obtidos resultados para eletrodos simples (em posição horizontal ou vertical) e para malhas quadradas. Esta ferramenta e os resultados obtidos permitem avaliar a resposta do aterramento submetido a uma injeção de uma corrente impulsiva típica de descargas atmosféricas. A validação da ferramenta desenvolvida é feita por meio de comparações dos resultados obtidos com resultados de casos semelhantes publicados na literatura técnica.

Le présent document a pour objet la modélisation transitoire non linéaire du comportement vibratoire des systèmes de frein aéronautiques. Le but est de reproduire numériquement l’apparition et le niveau des vibrations au cours du temps, afin de les maîtriser et d’adapter la conception du frein. Les essais de freinage mettent en évidence deux modes de vibration que sont le whirl et le squeal. Si les niveaux de ces vibrations deviennent trop importants, la structure de la roue et du train d’atterrissage peut être endommagée. Afin d’éviter de tels dommages, la conception du frein doit être adaptée. Pour réaliser cela, Messier-Bugatti-Dowty doit disposer d’un modèle capable de prédire les niveaux de vibration du frein au cours du temps pendant la phase de freinage. Le modèle doit avoir une précision suffisante, être en lien avec la maquette numérique et ne doit pas exiger de recalage. Un premier travail vise à améliorer le modèle éléments finis existant qui se révèle être trop imprécis. Une étude portant sur les effets gyroscopiques permet d’évaluer leur impact sur la fréquence et la stabilité des modes de whirl. Une modélisation plus complète du bâti d’essai améliore la précision de la fréquence du mode de squeal. Enfin, le mode de whirl est mieux simulé grâce au développement d’un modèle de pneumatique à partir de son analyse modale. Ce modèle est ensuite réduit afin de réaliser une intégration temporelle. Une sous-structuration permet de séparer l’ensemble des disques du frein, où le frottement et la non-linéarité se situent, du reste de la structure considérée comme linéaire. Trois techniques de réduction de l’ensemble des disques sont exposées. On évalue leur représentativité par rapport au modèle non-réduit en comparant les fréquences et la stabilité des modes propres. La première méthode est une représentation nodale de l’ensemble des disques. Les équations décrivant la non-linéarité et le frottement sont analytiques. Pour la deuxième méthode, la non-linéarité est déplacée à l’extrémité de l’ensemble des disques pour la découpler du frottement. La troisième méthode, plus ambitieuse et complexe, conserve à la fois l’emplacement de la nonlinéarité aux interfaces frottantes et la géométrie des disques. Une technique de réduction modale permet d’abaisser le nombre de degrés de liberté non linéaires. Pour clore ce rapport, des simulations transitoires sont calculées à partir des modèles réduits. Des études d’influences sont réalisées. Les paramètres étudiés sont le type d’algorithme d’intégration temporelle, l’amortissement introduit, la loi non linéaire, la pression hydraulique d’entrée et le coefficient de frottement. Leurs impacts sur les niveaux et la durée d’apparition des vibrations est évalué. / This report deals with the non-linear transient simulation of the dynamic behaviour of aeronautic brake systems. The objective is to reproduce the occurrence and level of vibrations versus time in order to control and adjust design consequently. The braking tests highlight two eigenmodes, which are called whirl and squeal. If the level of these vibrations becomes too high, the structures of the wheel and the landing gear may be damaged. To avoid damage, the design has to be adjusted. To achieve this, Messier-Bugatti-Dowty requires a model that is able to predict the levels of vibrations of the brake when it is braking. This model must have an adequate accuracy, be linked to the digital mockup and not require tuning. First, the existing finite element model has to be improved because its initial accuracy is not acceptable. A study about gyroscopic effects allows to assess their impact on the frequency and the stability of whirl modes. A complete modelling of the test frame improves the squeal modes’ frequency accuracy. At last, the whirl modes are better simulated due to the development of a tyre model based on modal analysis data. Then, the finite element model is reduced in order to perform a temporal integration. A substructuring allows to separate the set of brake discs (heat sink), where friction and non-linearities are located, from the rest of the structure which is considered linear. Three heat sink reduction techniques are proposed. Their representativeness are estimated compared to the non-reduced model. The first technique is a nodal description of the heat sink. The equations of friction and non-linearity are analytical. For the second technique, the non-linearity is displaced to the extremity of the heat sink to uncouple it from friction. The third technique, more ambitious and complex, keeps the location and non-linearity in friction interfaces and discs geometry. A reduction technique enables to decrease the number of non-linear degrees of freedom. As a conclusion, transient simulations are computed from reduced models. Sensitivity studies are performed. Studied parameters are the type of integration solver, introduced damping, non-linearities, hydraulic pressure, and friction coefficient. Their impacts on level and duration of occurrence of vibrations is estimated.

Parametry nahradniho obvodu synchronniho stroje znance ovlivnuji jeho chovani jak pri statickem provozu, tak predevsim pri nahlych dynamickych jevech a poruchovych stavech. Prace je zamerena na zhodnoceni dostupnych metod pro vypocet techto parametru pomoci Metody konecnych prvku. Prvni cast je venovana teoretickemu popisu zakladnich principu Metody Konecnych Prvku a jejich aplikaci na reseni problemu elektromagnetickeho pole v elektrickych strojich. Zaroven take shrnuje zakladni usporadani nahradniho obvodu synchroniho stroje, principi jeho konstrukce a zakladni funkci. Druha cast je venovana praktickemu vypoctu reaktanci nahradniho obvodu synchronniho stroje. S pomoci MKP jsou vypocteny synchronni reactance s uvazovanim vzajmeneho magnetickeho pusobeni proudu v d a q ose. Pro vypocet transientnich a subtransientnich reaktanci jsou navrzeny ctyri odlisne metody a jsou zhodnoceny z hlediska pozadovane presnosti vypoctu a narocnosti na vypocetni cas. Zaverecna cast popisuje zakladni merici metody pro urceni parametru nahrandniho obvodu na skutecnem stroji. Kapitola take obsahuje srovnani simulace trifazoveho zkratu synchronniho stroje s realnou zkouskou provedenou laboratorne. Zaver obsahuje srovnani jednotlivych metod a navrh optimalniho postupu pro vypocet zkoumanych parametru.

This diploma thesis deals with dynamic analysis of the steam turbine blade. This blade is part of the last row of low pressure level of steam turbine. Computational analysis has been performed in first part using FEM and software ANSYS. A Transient analysis has been used to solve forced vibrations. Main goal is to obtain a behavior of main stresses and its directions as a function of loading of the blade in the crack initiation area. Second part deals contain a an experiment. Experiment has been set up to perform a modal analysis which is necessary to obtain a fundamental numbers. Then the vibration of the blade has been performed. To perform this experiment same edge conditions as which has been used during the computational analysis. Goal is the same as in the first part – obtain a behavior of main stresses and its directions as a function of loading of the blade in the crack initiation area. The comparison of results obtained during experimental analysis and computational analysis has been performed in the last part of the thesis. Also the question if any geometrical nonlinearities appeared during analyses is answered.

This Master’s thesis deals with rotordynamics analysis of asynchronous generator rigid coupling with water turbine shaft. The purpose of this thesis was to determine influence of connection of generator with water turbine by rigid coupling to reactions in bearings. This problem was solved in two ways. The first option focused on model without influence of inaccuracies whereas there were inaccuracies included in the second option. The problem was solved by computational modeling by software Ansys.

The main goal of this thesis is dynamic analysis of shaft of booster pump QHD-200. In the beginning we offer basic equations generally used in dynamics and common/possible methods for their solution. Then we focus on specific problems connected with rotor dynamics and using simple Laval rotor we show their mathematical interpretation. Based on given load as accelerogram we create response spectrum in different ways (using different methods) and compare their CPU time demand. In the next part we make three finite element models of different physical level. We perform modal analysis and comparison of eigen frequencies, mode shapes. In following chapter transient analysis was performed on all models. As load we use accelerogram. We evaluate translations in particular positions of model. After that, we decide whether or not gap between impeller and stator case will be erased. Finally we make assessment of force acting upon bearings.

[pt] Testes de formação são normalmente realizados para determinar as propriedades rochosas do reservatório e os dados obtidos costumam ser interpretados com base no pressuposto que o reservatório é homogêneo na direção vertical e descrito por um modelo uni dimensional. No entanto, muitos reservatórios são compostos por diversas camadas que possuem características diferentes. Os poços de produção nesses reservatórios podem receber óleo de mais de uma camada. Em um sistema de reservatório estratificado, o comportamento da pressão e da temperatura não é necessariamente o mesmo de um sistema em camada única e raramente revela as mesmas propriedades médias de todo o sistema. Prever as características das camadas individuais é importante para descrever adequadamente o reservatório e melhorar o gerenciamento da produção. Este trabalho apresenta um modelo numérico, transiente-térmico para um sistema acoplado poço - reservatório 2D, levando-se em consideração efeitos Joule-Thompson responsáveis pelo aquecimento / resfriamento do fluido, expansão/compressão adiabática, além de efeitos de condução e convecção para o poço e o reservatório em um escoamento monofásico. A análise bidimensional do reservatório permite que se simule zonas de estratificação e barreiras. O modelo permite fluxo através de camadas adjascentes com propriedades de rocha diferentes. Pressão e temperatura a uma certa posição no poço produtor são avaliadas ao longo do tempo. Resultados mostram que a análise do transiente de pressão (PTA) e a análise do transiente de temperatura (TTA) podem ser utilizadas para caracterizar diferentes configurações de um reservatório estratificado. / [en] Well formation tests are usually performed to determine rock properties of a reservoir and the obtained data has often been interpreted based on an assumption that the reservoir is homogeneous in the vertical direction and described by a 1-D model. However, many reservoirs are found to be composed of different number of layers that have different characteristics. Production wells in such reservoirs may receive oil from more than one layer. In stratified reservoir system, the pressure and temperature behavior are not necessarily the same as in single layered system, and rarely reveals the same average properties of the entire system. The prediction of the characteristics of the individual layers is important to describe properly the reservoir and improve production management. This work presents a numerical transient-thermal model for a coupled wellbore/2D-reservoir considering Joule-Thompson heating/cooling, adiabatic fluid expansion/compression, conduction and convection effects for both wellbore and reservoir for a single-phase fluid flow. The two-dimensional reservoir model allows the analysis of stratified zones and barriers. The model allows cross flow between the adjacent layers with different rock properties.Wellbore temperature and pressure at a certain gauge depth are evaluated along the time. Results show how pressure transient analysis (PTA) and temperature transient analysis (TTA) can be used to characterize different configuration of stratified reservoirs.

For predicting the performance characteristics of highly utilized induction machines, commonly finite element analysis (FEA) is applied. If only the stationary behavior is of interest, both time-harmonic and transient calculations are feasible. Both procedures offer benefits and drawbacks regarding precision and computing time. In this paper the stationary torque-slip-characteristic of a squirrel cage induction machine for electric vehicle application is calculated by means of time-harmonic and transient FEA. The simulation results are compared to measurements taken on the test bench. For low saturation levels the time-harmonic and transient simulation produce nearly equivalent stationary torque results. For high saturation levels, the commonly used approach for nonlinear time-harmonic calculation is to use a corrected magnetization curve. It is shown that the assumption of sinusoidal time variation of the magnetic field strength usually made in this case results in a torque error increasing with the saturation level for the time-harmonic calculation.