Applying Computational Fluid Dynamics Methodologies to Hydro-Dynamic Loads in the Suppression Pool of a BWR

Author(s):  
David Mauritzson ◽  
Sven Perzon ◽  
Thomas Probert

The direct and indirect hydro-dynamic loads in the suppression pool of the BWR fleet were developed well before 2000. These loads were based on scaled experiments and numerical solutions using one dimensional models. The analysis was cumbersome and conservatisms were added in multiple steps. As these loads spread throughout the containment, as global vibrations, they are generally a part of all structural verification inside the containment. This fact has made it hard and unpractical to challenge and revise these loads, as a change could lead to significant re-work. A consequence of this is that loads in the suppression pools have seldom been revisited, regardless if they cause local or global vibrations. This is problematic when new equipment is needed. The design of this equipment suffers from the very conservative loads. Experience has shown that these loads can be challenged and refined using updated techniques, leading to significantly lower loads. This was realized during the modernization of Swedish plant Oskarshamn Unit 2 with Mark II containment, where loads following pool swell proved to be particularly challenging and it was decided to investigate possibilities to reduce the conservative loads. Due to the large scale of the condensation pool coupled with the transient and small scaled condensation of steam at the drywell vent pipe nozzle full CFD resolution is not feasible. Instead lumped models in GOTHIC was used to increase resolution from the 1D approach of normal containment analysis, to a resolution that can account for the features of the condensation pool. This showed that the pool swell was less uniform than initially thought, leading to fewer objects affected as well as lower loads on objects that suffered from these loads. A full CFD analysis was then used to resolve phenomena working on even shorter time scales leading to a complete rework of all local loads. The loads addressed using updated codes and modelling techniques was pool swell impact (PSI), pool swell drag (PSD) and local drag loads due to pressure relief valve opening, LV/SRV. The current work shows that using updated modelling techniques and aligning results with previous analysis and documentation, it is possible to reduce loads for some events in the suppression pool without violating safety for the power plant. The results from the GOTHIC model shows that the cylindrical shape of the pool will create an uneven velocity distribution radially at pool swell resulting in much smaller loads at the outer boundary. CFD analysis of the LV/SRV event shows that the loads are reduced in comparison with previous methodology and this is mainly because of shadow effects. The calculations also verifies assumptions used in the previous methodology.

Acta Numerica ◽  
2003 ◽  
Vol 12 ◽  
pp. 267-319 ◽  
Author(s):  
Roland W. Freund

In recent years, reduced-order modelling techniques based on Krylov-subspace iterations, especially the Lanczos algorithm and the Arnoldi process, have become popular tools for tackling the large-scale time-invariant linear dynamical systems that arise in the simulation of electronic circuits. This paper reviews the main ideas of reduced-order modelling techniques based on Krylov subspaces and describes some applications of reduced-order modelling in circuit simulation.


Author(s):  
W. Dempster ◽  
C. K. Lee ◽  
J. Deans

The design of safety relief valves depends on knowledge of the expected force-lift and flow-lift characteristics at the desired operating conditions of the valve. During valve opening the flow conditions change from seal-leakage type flows to combinations of sub-sonic and supersonic flows It is these highly compressible flow conditions that control the force and flow lift characteristics. This paper reports the use of computational fluid dynamics techniques to investigate the valve characteristics for a conventional spring operated 1/4” safety relief valve designed for gases operating between 10 and 30 bar. The force and flow magnitudes are highly dependent on the lift and geometry of the valve and these characteristics are explained with the aid of the detailed information available from the CFD analysis. Experimental determination of the force and flow lift conditions has also been carried out and a comparison indicates good correspondence between the predictions and the experiment. However, attention requires to be paid to specific aspects of the geometry modeling including corner radii and edge chamfers to ensure satisfactory prediction.


2005 ◽  
Vol 129 (2) ◽  
pp. 227-233 ◽  
Author(s):  
Simon Li ◽  
Li Chen

We have developed a decomposition-based rapid redesign methodology for large and complex computational redesign problems. While the overall methodology consists of two general steps: diagnosis and repair, in this paper we focus on the repair step in which decomposition patterns are utilized for redesign planning. Resulting from design diagnosis, a typical decomposition pattern solution to a given redesign problem indicates the portions of the design model necessary for recomputation as well as the interaction part within the model accountable for design change propagation. Following this, in this paper we suggest repair actions with an approach derived from an input pattern solution, to generate a redesign road map allowing for taking a shortcut in the redesign solution process. To do so, a two-stage redesign planning approach from recomputation strategy selection to redesign road map generation is proposed. An example problem concerning the redesign of a relief valve is used for illustration and validation.


2012 ◽  
Vol 446-449 ◽  
pp. 837-840
Author(s):  
Yu Zhao ◽  
Shu Fang Yuan ◽  
Jian Wei Zhang

The underwater structure of power house is major structure under the dynamic loads of unit. The vibration problem is very common in operation. So the structures should have sufficient stiffness to resist dynamic loads of unit. This paper establishes three-dimensional finite element models with finite element analysis software—ANSYS. Dynamic characteristics of the power house and dynamic responses of structure under earthquake are analyzed. The results of the computation show that fluid-solid coupling may be ignored when studying dynamic characteristics of structures of the underground power house.


Author(s):  
C. Guardino ◽  
J. W. Chew ◽  
N. J. Hills

The effects of surface roughness on air-riding seals are investigated here using the Rayleigh-pad as an example. Both incompressible and compressible flows are considered using both CFD analysis and analytical/numerical solutions of the Reynolds equation for various 2D or 3D roughness patterns on the stationary wall. A ‘unit-based’ approach for incompressible flows has also been employed and is shown to be computationally much less expensive than the full-geometry solution. Results are presented showing the effect of surface roughness on the net lift force. The effects of varying the Reynolds number are demonstrated, as well as comparative results for static stiffness.


2008 ◽  
Vol 74 (2) ◽  
pp. 215-232 ◽  
Author(s):  
I. F. SHAIKHISLAMOV

AbstractReconnection physics at micro-scales is investigated in an electron magnetohydrodynamics frame. A new process of collapse of the neutral current sheet is demonstrated by means of analytical and numerical solutions. It shows how at scales smaller than ion inertia length a compression of the sheet triggers an explosive evolution of current perturbation. Collapse results in the formation of a intense sub-sheet and then an X-point structure embedded into the equilibrium sheet. Hall currents associated with this structure support high reconnection rates. Nonlinear static solution at scales of the electron skin reveals that electron inertia and small viscosity provide an efficient mechanism of field lines breaking. The reconnection rate does not depend on the actual value of viscosity, while the maximum current is found to be restricted even for space plasmas with extremely rare collisions. The results obtained are verified by a two-fluid large-scale numerical simulation.


2002 ◽  
Vol 465 ◽  
pp. 99-130 ◽  
Author(s):  
A. V. OBABKO ◽  
K. W. CASSEL

Numerical solutions of the unsteady Navier–Stokes equations are considered for the flow induced by a thick-core vortex convecting along a surface in a two-dimensional incompressible flow. The presence of the vortex induces an adverse streamwise pressure gradient along the surface that leads to the formation of a secondary recirculation region followed by a narrow eruption of near-wall fluid in solutions of the unsteady boundary-layer equations. The locally thickening boundary layer in the vicinity of the eruption provokes an interaction between the viscous boundary layer and the outer inviscid flow. Numerical solutions of the Navier–Stokes equations show that the interaction occurs on two distinct streamwise length scales depending upon which of three Reynolds-number regimes is being considered. At high Reynolds numbers, the spike leads to a small-scale interaction; at moderate Reynolds numbers, the flow experiences a large-scale interaction followed by the small-scale interaction due to the spike; at low Reynolds numbers, large-scale interaction occurs, but there is no spike or subsequent small-scale interaction. The large-scale interaction is found to play an essential role in determining the overall evolution of unsteady separation in the moderate-Reynolds-number regime; it accelerates the spike formation process and leads to formation of secondary recirculation regions, splitting of the primary recirculation region into multiple corotating eddies and ejections of near-wall vorticity. These eddies later merge prior to being lifted away from the surface and causing detachment of the thick-core vortex.


Author(s):  
Lawrence J. De Chant

PurposeAlthough most physical problems in fluid mechanics and heat transfer are governed by nonlinear differential equations, it is less common to be confronted with a “so – called” implicit differential equation, i.e. a differential equation where the highest order derivative cannot be isolated. The purpose of this paper is to derive and analyze an implicit differential equation that arises from a simple model for radiation dominated heat transfer based upon an unsteady lumped capacitance approach.Design/methodology/approachHere we discuss an implicit differential equation that arises from a simple model for radiation dominated heat transfer based upon an unsteady lumped capacitance approach. Due to the implicit nature of this problem, standard integration schemes, e.g. Runge‐Kutta, are not conveniently applied to this problem. Moreover, numerical solutions do not provide the insight afforded by an analytical solution.FindingsA predictor predictor‐corrector scheme with secant iteration is presented which readily integrates differential equations where the derivative cannot be explicitly obtained. These solutions are compared to numerical integration of the equations and show good agreement.Originality/valueThe paper emphasizes that although large‐scale, multi‐dimensional time‐dependent heat transfer simulation tools are routinely available, there are instances where unsteady, engineering models such as the one discussed here are both adequate and appropriate.


Author(s):  
Jianxin Zhu ◽  
Xuedong Chen ◽  
YunRong Lu ◽  
Zhibin Ai ◽  
Weihe Guan

The shutdown of charge gas compressor in large-scale ethylene-cracking plant always involves emergency pressure relief of charge gas through multiple safety valves. The emergency relief capacity plays an important role on the safety of the overall plant. In this paper, by studying the difference between the configuration of the pressure relief system of two 1000 KTA ethylene-cracking plants (the inner diameters of the charge gas pipeline in both plants are 2m, while the number of same-sized relief valves are 28 and 19, respectively), the relief capacity of multiple relief valves is studied and compared with empirical calculation and numerical analysis. It is found that, due to the interruption of fluid flow when compressor is emergency shutdown, the upstream pressure of each relief valve increase steadily with the continuously make-up of the charge gas, but the difference between the inlet pressure of all relief valves can be neglected. With the increase of the upstream pressure, the opening of relief valves is determined mainly by the set pressure. In multiple valves pressure relief scenario, normally the downstream valves have greater relief capacity than those upstream valves if both relief valves have the same back pressure. Also, by analysis it is noted that the pressure relief capacities of multiple relief valves in both plants are sufficient. The minimum number of relief valves required for process safety is obtained. The maximum achievable Safety Integrity Level (SIL) of pressure relief system is determined by calculation of the reliability of the redundant relief valves. The analysis is used for determination of the SIL of the pressure relief system. The finding is also significant for determination of the required capacity of multiple relief valves.


Sign in / Sign up

Export Citation Format

Share Document