• Nuclear Engineering and Technology
• /
• v.49 no.4
• /
• pp.675-691
• /
• 2017

Jets in cross flow are of fundamental industrial importance and play an important role in validating turbulence models. Two jet configurations related to thermal fatigue phenomena are investigated: ${\bullet}$ T-junction of circular tubes where a heated jet discharges into a cold main flow and ${\bullet}$ Rectangular jet marked by a scalar discharging into a main flow in a rectangular channel. The T-junction configuration is a classical test case for thermal fatigue phenomena. The Vattenfall T-junction experiment was already subject of an OECD/NEA benchmark. A LES modelling and calculation strategy is developed and validated on this data. The rectangular-jet configuration is important for basic physical understanding and modelling and has been analyzed experimentally at CEA. The experimental work was focused on turbulent mixing between a slightly heated rectangular jet which is injected perpendicularly into the cold main flow of a rectangular channel. These experiments are analyzed for the first time with LES. The overall results show a good agreement between the experimental data and the CFD calculation. Mean values of velocity and temperature are well captured by both RANS calculation and LES. The range of critical frequencies and their amplitudes, however, are only captured by LES.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.3
• /
• pp.234-244
• /
• 1994

Longshore currents driven by monochromatic waves have been described using 2-equation $textsc{k}$-ι turbulence transport model. When using $textsc{k}$-ι closure both profiles of eddy viscosity and current velocity are found to be satisfactory. Several terms of ι equation are related to various variables concerned with turbulence mechanism. New form of turbulence frequency used in ι equation is suggested in the present approach, and non-dimensional parameters are evaluated by comparing the computational results with the laboratory measurements. Various values of a large range are applied to the non-dimensional parameters for the sensitivity test and in order to improve the predictability common values of constants are deduced, which produce similarly good computational results for the well-controlled laboratory measurements.

• Nuclear Engineering and Technology
• /
• v.38 no.2
• /
• pp.147-154
• /
• 2006

Thermal fatigue of the coolant circuits of PWR plants is a major issue for nuclear safety. The problem is especially accute in mixing zones, like T-junctions, where large differences in water temperature between the two inlets and high levels of turbulence can lead to large temperature fluctuations at the wall. Until recently, studies on the matter had been tackled at EDF using steady methods: the fluid flow was solved with a CFD code using an averaged turbulence model, which led to the knowledge of the mean temperature and temperature variance at each point of the wall. But, being based on averaged quantities, this method could not reproduce the unsteady and 3D effects of the problem, like phase lag in temperature oscillations between two points, which can generate important stresses. Benefiting from advances in computer power and turbulence modelling, a new methodology is now applied, that allows to take these effects into account. The CFD tool Code_Saturne, developped at EDF, is used to solve the fluid flow using an unsteady L.E.S. approach. It is coupled with the thermal code Syrthes, which propagates the temperature fluctuations into the wall thickness. The instantaneous temperature field inside the wall can then be extracted and used for structure mechanics computations (mainly with EDF thermomechanics tool Code_Aster). The purpose of this paper is to present the application of this methodology to the simulation of a straight T-junction mock-up, similar to the Residual Heat Remover (RHR) junction found in N4 type PWR nuclear plants, and designed to study thermal striping and cracks propagation. The results are generally in good agreement with the measurements; yet, in certain areas of the flow, progress is still needed in L.E.S. modelling and in the treatment of instantaneous heat transfer at the wall.

• Wind and Structures
• /
• v.13 no.2
• /
• pp.169-189
• /
• 2010

This paper describes the use of coupled Computational Fluid Dynamics (CFD) and Rigid Body Dynamics (RBD) in modelling the aerodynamic behaviour of wind-borne plate type objects. Unsteady 2D and 3D Reynolds Averaged Navier-Stokes (RANS) CFD models are used to simulate the unsteady and non-uniform flow field surrounding static, forced rotating, auto-rotating and free-flying plates. The auto-rotation phenomenon itself is strongly influenced by vortex shedding, and the realisable k-epsilon turbulence modelling approach is used, with a second order implicit time advancement scheme and equal or higher order advection schemes for the flow variables. Sequentially coupling the CFD code with a RBD solver allows a more detailed modelling of the Fluid-Structure Interaction (FSI) behaviour of the plate and how this influences plate motion. The results are compared against wind tunnel experiments on auto-rotating plates and an existing 3D analytical model.

• Wind and Structures
• /
• v.5 no.2_3_4
• /
• pp.177-192
• /
• 2002

Computation solutions for the flow around a cube, which were generated as part of the Computational Wind Engineering 2000 Conference Competition, are compared with full-scale measurements. The three solutions shown all use the RANS approach to predict mean flow fields. The major differences appear to be related to the use of the standard $k-{\varepsilon}$, the MMK $k-{\varepsilon}$ and the RNG $k-{\varepsilon}$ turbulence models. The inlet conditions chosen by the three modellers illustrate one of the dilemmas faced in computational wind engineering. While all modeller matched the inlet velocity profile to the full-scale profile, only one of the modellers chose to match the full-scale turbulence data. This approach led to a boundary layer that was not in equilibrium. The approach taken by the other modeller was to specify lower inlet turbulent kinetic energy level, which are more consistent with the turbulence models chosen and lead to a homogeneous boundary layer. For the $0^{\circ}$ case, wind normal to one face of the cube, it is shown that the RNG solution is closest to the full-scale data. This result appears to be associated with the RNG solution showing the correct flow separation and reattachment on the roof. The other solutions show either excessive separation (MMK) or no separation at all (K-E). For the $45^{\circ}$ case the three solutions are fairly similar. None of them correctly predicting the high suctions along the windward edges of the roof. In general the velocity components are more accurately predicted than the pressures. However in all cases the turbulence levels are poorly matched, with all of the solutions failing to match the high turbulence levels measured around the edges of separated flows. Although all of the computational solutions have deficiencies, the variability of results is shown to be similar to that which has been obtained with a similar comparative wind tunnel study. This suggests that the computational solutions are only slightly less reliable than the wind tunnel.

• Journal of the korean society of oceanography
• /
• v.32 no.2
• /
• pp.75-84
• /
• 1997

Relation between the length scale and the wall proximity function in the Mellor-Yamada level 2.5 turbulence closure model has been investigated through various experiments using a range of wall proximity functions. The model performance has been evaluated quantitatively by comparing with laboratory data for wind-driven flow (Baines and Knapp, 1965) and for open-channel flows without and with adverse wind action (Tsuruya, 1985). Comparison shows that a symmetric wall proximity function used by Blumberg and Mellor(1987) gives rise to current profiles with better accuracy than asymmetric wall proximity functions considered. It is noted that in modelling homogeneous flows the length scale 1= 0.31${\|}$z${\|}$(1+z/h) can be used with tolerable accuracy.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.3
• /
• pp.653-661
• /
• 1994

In this paper the CSCM type upwind flux difference splitting Navier-Stokes method has been applied to study the ARL-SL19 supersonic/transonic compressor cascade flow. H-type grid was chosen for its simplicity in applying cyclic tridiagonal matrix algorithm along with conventional slip/no-slip boundary conditions. The thin-layer algebraic model of Baldwin-Lomax was employed for the calculation of turbulent flows. The test case inlet Mach No. was 1.612 and inlet/exit pressure ratio($P_2/P_1$) was 2.15. The results were compared with experimental results from current method were compared well in suction surface with the experiments and other computational results; however, not well in pressure surface. It might be due to the complex flowfields such as shock/boundary layer interaction, turbulence, and flow separation, etc. In the future, a proper turbulence modelling and adaptive grid system will be studied to improve the solution quality.

• Wind and Structures
• /
• v.5 no.2_3_4
• /
• pp.89-100
• /
• 2002

Accidental gaseous losses from industrial processes can pose considerable health and environmental risks but assessing their health, safety and environmental impact is problematic. Improved understanding and simulation of the dispersion of emissions in the vicinity of storage tanks is required. The present study aims to assess the capability of the turbulence closures and meshing alternatives in a commercially available CFD code for predicting dispersion in the vicinity of cubes and circular cylindrical storage tanks. The performance of the $k-{\varepsilon}$ and Reynolds Stress turbulence models and meshing alternatives for these cases are compared to experimental data. The CFD simulations are very good qualitatively and, in many cases, quantitatively. A mesh with prismatic elements is more accurate than a tetrahedral mesh. Overall the Reynolds stress model performs slightly better than the $k-{\varepsilon}$ model.

• Wind and Structures
• /
• v.4 no.6
• /
• pp.469-480
• /
• 2001

The Leipzig Wind Profile is generally known as a typical neutral planetary boundary layer flow. But it became clear from the present research that it was not completely neutral but weakly stable. We examined whether we could simulate the Leipzig Wind Profile by using a ($k-{\varepsilon}$) turbulence model including the equation of potential temperature. By solving analytically the Second Moment Closure Model under the assumption of local equilibrium and under the condition of a stratified flow, we expressed the turbulent diffusion coefficients (both momentum and thermal) as functions of flux Richardson number. Our ($k-{\varepsilon}$) turbulence model which included the equation of potential temperature and the turbulent diffusion coefficients varying with flux Richardson number reproduced the Leipzig Wind Profile.

• International Journal of High-Rise Buildings
• /
• v.11 no.4
• /
• pp.331-346
• /
• 2022

Reproducing the horizontally homogeneous atmospheric boundary layer in computational wind engineering is essential for predicting the wind loads on structures. One of the important issues is to use fully developed inflow conditions, which will lead to the consistence problem between inflow condition and internal roughness. Thus, by analyzing the previous results of computational fluid dynamic modeling turbulent horizontally homogeneous atmospheric boundary layer, we modify the past hypotheses, detailly derive a new type of inflow condition for standard k-ε turbulence model. A group of remedial approaches including formulation for wall shear stress and fixing the values of turbulent kinetic energy and turbulent dissipation rate in first wall adjacent layer cells, are also derived to realize the consistence of inflow condition and internal roughness. By combing the approaches with four different sets of inflow conditions, the well-maintained atmospheric boundary layer flow verifies the feasibility and capability of the proposed inflow conditions and remedial approaches.