Numerical Simulation and Modelling Convention of Unsteady Fluidelastic Forces of Tube Arrays

Author(s):  
Mingjie Zhang ◽  
Xu Wang ◽  
Ole Øiseth

Abstract This paper presents a numerical investigation on the unsteady fluidelastic forces of tube arrays. The key focus is on the consistency between the unsteady fluidelastic force model and the quasi-steady model for tube arrays at large reduced flow velocities, as well as comparing two well-known conventions for the unsteady model. Two-dimensional unsteady Reynolds-averaged Navier-Stokes (URANS) simulations are used to prove that the viscous damping coefficients of Tanaka's convention (Tanaka and Takahara, 1981) approach their quasi-steady values as the reduced flow velocity approaches infinity, whereas the hysteretic damping coefficients of Chen's modified convention (Chen et al., 1983) always approach zero and hence result in low-resolution data plots as the reduced flow velocity becomes large. The non-constant viscous damping coefficients of Tanaka's experimental data at high reduced flow velocities (which motivated the introduction of Chen's modified convention) might be induced by a systematic identification error in the phase of the fluidelastic force. A row of three flexible cylinders is used as a numerical example to analyse the effect of systematic phase error on the predicted stability boundary of the fluidelastic instability. Although identical fluidelastic forces are simulated by using the two conventions, Tanaka's convention is recommended due to its compatibility with the quasi-steady theory and optimal resolutions of data plots over any range of reduced flow velocities.

Author(s):  
Roger Y. Lu ◽  
David D. Steel

PWR fuel assembly damping is a key parameter in seismic/LOCA safety analysis. The damping coefficients of a fuel assembly in air, still water and flowing water are significantly different. Several researchers and engineers have published their results and methods in the past. With this paper, PWR fuel assembly damping was studied and tested in air, still water, and flowing water (including flowrate and temperature variation). The damping coefficients were obtained by the initial displacement and first response method. The coefficients are also compared with published data. Several conclusions are obtained. • The damping obtained from the tests in air gives the damping component of assembly structure damping. From the comparison of the damping in air with still water the amount of viscous damping can be determined. The viscous damping component is the effect of still water on damping. The amount of viscous damping is represented by the increase in the damping ratio from air to still water at room temperature. The results show that damping in still water is approximately two times the damping in air. • The temperature effect on damping in still water is minimal. In flowing water, the results show a very slight effect of temperature, as the damping slightly decreases with an increase in temperature. This temperature effect is much smaller than the data scatter observed in most damping measurement tests under the same test conditions. • The damping is significantly affected by flowing water. For relatively low flow velocities, compared to in-core conditions, the damping coefficient is around two times the damping in still water. For intermediate to high flow velocities, all damping coefficients are 2.5 times higher than that in still water. For high velocities and large displacement, the damping coefficient can be over 3 times higher than that in still water. The flow velocity appears to be acting on the system by suppressing the motion of the assembly. Additional damping due to flowing water is called hydraulic damping, which is generated by hydraulic force. When a fuel assembly vibrates in flowing water, the assembly is trying to change the flow direction and momentum, but the flow mass wants to retain its pure axial direction which suppresses the motion of the assembly.


Author(s):  
Marwan Hassan

Fluidelastic instability is regarded as the most complex and destructive flow excitation mechanism in heat exchanger tube arrays subjected to cross fluid flow. Several attempts have been made for modelling fluidelastic instability in tube arrays in order to predict the stability threshold. However, fretting wear prediction requires a nonlinear computation of the tube dynamics in which proper modelling of the fluid forcing function is essential. In this paper, a time domain simulation of fluidelastic instability is presented for a single flexible tube in an otherwise rigid array subjected to cross fluid flow. The model is based on the unsteady flow theory proposed by Lever and Weaver [1] and Yetisir and Weaver [2]. The developed model has been implemented in INDAP (Incremental Nonlinear Dynamic Analysis Program), an in-house finite element code. Numerical investigations were performed for two linear tube-array geometries and compared with published experimental data. A reasonable agreement between the numerical simulation and the experimental results was obtained. The fluidelastic force model was also coupled with a tube/support interaction model. The developed numerical model was utilized to study a loosely-supported cantilever tube subjected to air flow. Tube-to-support clearance, random excitation level, and flow velocity were then varied. The results indicated that the loose support has a stabilizing effect on the tube response. Both rms impact force and normal work rate increased as a result of increasing the flow velocity or the support radial clearance. Contact ratio exhibited a sharp increase at a flow velocity higher than the instability threshold of the first unsupported mode. In addition, an interesting behaviour has been observed, namely the change of tube’s equilibrium position due to fluid forces. This causes a single-sided impact. At a higher turbulence level, double-sided impact conditions were dominant. The influence of these dynamic regimes on the tube/support parameters was also addressed.


2020 ◽  
Vol 41 (Supplement_2) ◽  
Author(s):  
A Zagatina ◽  
M Novikov ◽  
N Zhuravskaya ◽  
V Balakhonov ◽  
S Efremov ◽  
...  

Abstract Background Stenosis of a coronary artery results in an increase in flow velocity in the pathologic segment. Effective grafting should decrease the stenotic native coronary velocity according to hemodynamic law. The range of decreased velocity before and after cardiac surgery can hypothetically reflect the effectiveness of a graft. The aim of the study is to determine if measuring coronary flow velocity changes during coronary artery bypass grafting (CABG) can predict intraoperative myocardial infarction. Methods One hundred sixty-six (166) consecutive patients (121 men, 64±9 years old) referred for cardiac surgery, were prospectively included in the study. A standard basic perioperative transesophageal echocardiography (TEE) examination was performed with additional scans of the left main, left anterior descending (LAD), and circumflex (LCx) arteries' proximal segments. Measurements of coronary flow velocities were performed before and after grafting in the same sites of the arteries. The maximal value of cardiac troponin I (cTnI) after CABG and the additive criteria were accounted for in the analysis as it is described in the expert consensus document for Type 5 myocardial infarction (MI) definition. Results One hundred sixty-three patients (98%) had arterial hypertension, 28 patients (17%) had diabetes mellitus, 35 patients (21%) were currently smokers. The feasibility of coronary flow assessment during cardiac operations was 95%. Before grafting, the mean velocity in the left main artery was 91±49 cm/s, in LAD 101±35 cm/s, and in LCx 117±49 cm/s. There was a significant correlation between changes in coronary flow velocities during operation and the value of cTnI (R=0.34, p<0.0001). Ten patients met the criteria for Type 5 MI. There were no differences in age, body mass index, number of coronary arteries with stenoses, frequency of prior MI, ejection fraction or coronary flow velocity before surgery in patients with and without Type 5 MI. The group of patients with Type 5 MI had an increase in native artery velocities during surgery in comparison with patients without MI, who had a significant decrease in coronary flow velocity after grafting (30±48 vs. −10±30 cm/s; p<0.0006). Increases in native coronary velocities greater than 3 cm/s predicted Type 5 MI with 81% accuracy (sensitivity 88%, specificity 70%). Conclusion Coronary flow velocity assessment during cardiac surgery could predict an elevation of cardiac troponins and Type 5 MI. Funding Acknowledgement Type of funding source: None


1997 ◽  
Vol 134 (4) ◽  
pp. 557-561
Author(s):  
KATSUHIRO NAKAYAMA

Miocene subtidal sandwave deposits in southwest Japan were influenced by periodic flow and steady flow. The sandwave deposits can be divided into five units, based on lithofacies and thickness. In order of accretion, unit 1 consists of unidirectional sand bedforms without mud drapes, unit 2 of unidirectional sand bedforms with thin, discontinuous mud drapes, unit 3 of bidirectional sand bedforms with thin continuous mud drapes, and units 4 and 5 of relatively thinner and smaller bidirectional sand bedforms with continuous mud drapes. The thickness of units 1 to 3 increase progressively to 2.6 m, and units 4 to 5 subsequently decrease from 2.0 to 1.0 m. Variations between the units are due to differing combinations of periodic and steady flow velocities. Palaeoflow velocity is estimated from grain size and unit thickness. Depth-mean velocities of steady flow components gradually decrease from 0.72 ms−1 to 0.16 ms−1 with unit accumulation.


2013 ◽  
Vol 59 (214) ◽  
pp. 315-326 ◽  
Author(s):  
A. Richter ◽  
D.V. Fedorov ◽  
M. Fritsche ◽  
S.V. Popov ◽  
V.Ya. Lipenkov ◽  
...  

AbstractRepeated Global Navigation Satellite Systems (GNSS) observations were carried out at 50 surface markers in the Vostok Subglacial Lake (East Antarctica) region between 2001 and 2011. The horizontal ice flow velocity vectors were derived with accuracies of 1 cm a−1 and 0.5°, representing the first reliable information on ice flow kinematics in the northern part of the lake. Within the lake area, ice flow velocities do not exceed 2 m a−1. The ice flow azimuth is southeast in the southern part of the lake and turns gradually to east-northeast in the northern part. In the northern part, as the ice flow enters the lake at the western shore, the velocity decreases towards the central lake axis, then increases slightly past the central axis. In the southern part, a continued acceleration is observed from the central lake axis across the downstream grounding line. Based on the observed flow velocity vectors and ice thickness data, mean surface accumulation rates are inferred for four surface segments between Ridge B and Vostok Subglacial Lake and show a steady increase towards the north.


Author(s):  
Tomomichi Nakamura

Fluidelastic vibration of tube arrays caused by cross-flow has recently been highlighted by a practical event. There have been many studies on fluidelastic instability, but almost all works have been devoted to the tube-vibration in the transverse direction to the flow. For this reason, there are few data on the fluidelastic forces for the in-flow movement of the tubes, although the measured data on the stability boundary has gradually increased. The most popular method to estimate the fluidelastic force is to measure the force acting on tubes due to the flow, combined with the movement of the tubes. However, this method does not give the physical explanation of the root-cause of fluidelastic instability. In the work reported here, the in-flow instability is assumed to be a nonlinear phenomenon with a retarded or delayed action between adjacent tubes. The fluid force acting on tubes are estimated, based on the measured data in another paper for the fixed cylinders with distributed pressure sensors on the surface of the cylinders. The fluid force acting on the downstream-cylinder is assumed in this paper to have a delayed time basically based on the distance between the separation point of the upstream-cylinder to the re-attachment point, where the fluid flows with a certain flow velocity. Two models are considered: a two-cylinder and three–cylinder models, based on the same dimensions as our experimental data to check the critical flow velocity. Both models show the same order of the critical flow velocity and a similar trend for the effect of the pitch-to-diameter ratio of the tube arrays, which indicates this analysis has a potential to explain the in-flow instability if an adequate fluid force is used.


1981 ◽  
Vol 103 (2) ◽  
pp. 130-135 ◽  
Author(s):  
S. S. Chen ◽  
J. A. Jendrzejczyk

Experiments are conducted to determine the damping for a tube in tube arrays subjected to liquid cross-flow; damping factors in the lift and drag directions are measured for in-line and staggered arrays. It is found that: 1) fluid damping is not a constant, but a function of flow velocity; 2) damping factors in the lift and drag directions are different; 3) fluid damping depends on the tube location in an array; 4) flow velocity-dependent damping is coupled with vortex shedding process and fluid-elastic instability; and 5) flow velocity-dependent damping may be negative. This study demonstrates that flow velocity-dependent damping is important. These characteristics should be properly taken into account in the mathematical modeling of tube arrays subjected to cross-flow.


2015 ◽  
Vol 52 (2) ◽  
pp. 3-12
Author(s):  
A. Kalnacs ◽  
J. Kalnacs ◽  
A. Mutule ◽  
V. Entins

Abstract In the rivers of Latvia and of many other countries the flow velocity in the places that are most suitable for installation of hydrokinetic devices is 0.4 to 0.9 m/s. In a stream or a river the hydrokinetic devices can reach full efficiency starting from about twice higher flow velocities. It is advisable to at least double this velocity thus increasing the efficiency and power output of the hydrokinetic devices installed in such places. Since Latvia has abundance of slow rivers and almost none are fast, research in this field is of high importance. Diversified technical methods are known that allow increasing substantially the efficiency of hydrokinetic devices. These methods include the use of diffusers, concentrators, different types of other channelling devices and means of flow control. Desirable effects are achieved through changing the cross-section and/or direction of a flow, its pressure, minimizing the turbulence, etc. This work substantiates the use of such devices for increasing the efficiency of hydrokinetic devices. A method is proposed for evaluation of the effects on power output gained owing to the use of channelling devices. Results show that the efficiency of hydrokinetic devices can be increased by at least 110%.


Sign in / Sign up

Export Citation Format

Share Document