Oscillation of Mixing Jets in a Long Horizontal Tank

Author(s):  
Zhao Yu ◽  
Kishore Kar ◽  
Tyrone Storrs ◽  
Christopher Jian ◽  
Richard Cope ◽  
...  

In contrast to mixing in vertical tanks, jet mixing in long horizontal tanks is scarcely investigated in the literature. It is known that jet mixing in long horizontal tanks is more difficult when compared to short tanks, as more liquid volume must be recirculated through the jets. In this study, computational fluid dynamics (CFD) simulations are conducted for the flow in a horizontal cylindrical tank with a length-to-diameter ratio of 3:1 and a nominal volume of 112,560 L (liquid volume of 75,708 L, i.e. 20,000 gallons). A pair of back-to-back Coldrey nozzles is placed near the center of the tank bottom, and each nozzle directs its jet towards the corresponding vessel end. An intriguing phenomenon is observed in the transient simulations, where the turbulent jets oscillate in both horizontal and vertical directions with a low frequency. In order to determine the source of such oscillation, a number of simulations are conducted to explore the effects of mesh type and size, boundary condition on the free surface, turbulence model, and time step. Oscillation persists in all cases, indicating that it is unlikely the result of some numerical instability. The oscillation also appears to be insensitive to the Reynolds number or symmetry in the nozzle or tank geometry. Another simulation with a single jet also shows the oscillatory flow behavior, and thus the oscillation is more likely to be caused by interaction between the jet and the recirculating flow in the tank, rather than interaction between the two jets. Further analysis of the jet velocity profile suggests that the secondary flow on the cross section of the jets might also contribute to the oscillation. While similar confined jet oscillations due to Coanda effect and blind cavity effect have been previously observed in small cavities by both 2D numerical simulations and laboratory scale experiments, this study shows that such oscillation also exists in industrial scale horizontal tanks. The oscillatory motion of the liquid may lead to improved mixing in the tank.

2008 ◽  
Vol 617 ◽  
pp. 231-253 ◽  
Author(s):  
DANIEL J. BODONY ◽  
SANJIVA K. LELE

An analysis of the sound radiated by three turbulent, high-speed jets is conducted using Lighthill's acoustic analogy (Proc. R. Soc. Lond. A, vol. 211, 1952, p. 564). Computed by large eddy simulation the three jets operate at different conditions: a Mach 0.9 cold jet, a Mach 2.0 cold jet and a Mach 1.0 heated jet. The last two jets have the same jet velocity and differ only by temperature. None of the jets exhibit Mach wave characteristics. For these jets the comparison between the Lighthill-predicted sound and the directly computed sound is favourable for all jets and for the two angles (30° and 90°, measured from the downstream jet axis) considered. The momentum (ρuiuj) and the so-called entropy [p − p∞ − a∞2(ρ − ρ∞)] contributions are examined in the acoustic far field. It is found that significant phase cancellation exists between the momentum and entropy components. It is observed that for high-speed jets one cannot consider ρuiuj and (p′ − a∞2ρ′)δij as independent sources. In particular the ρ′ūxūx component of ρuiuj is strongly coupled with the entropy term as a consequence of compressibility and the high jet velocity and not because of a linear sound-generation mechanism. Further, in more usefully decoupling the momentum and entropic contributions, the decomposition of Tij due to Lilley (Tech. Rep. AGARD CP-131 1974) is preferred. Connections are made between the present results and the quieting of high-speed jets with heating.


2008 ◽  
Vol 130 (12) ◽  
Author(s):  
D. Keith Walters ◽  
Davor Cokljat

An eddy-viscosity turbulence model employing three additional transport equations is presented and applied to a number of transitional flow test cases. The model is based on the k-ω framework and represents a substantial refinement to a transition-sensitive model that has been previously documented in the open literature. The third transport equation is included to predict the magnitude of low-frequency velocity fluctuations in the pretransitional boundary layer that have been identified as the precursors to transition. The closure of model terms is based on a phenomenological (i.e., physics-based) rather than a purely empirical approach and the rationale for the forms of these terms is discussed. The model has been implemented into a commercial computational fluid dynamics code and applied to a number of relevant test cases, including flat plate boundary layers with and without applied pressure gradients, as well as a variety of airfoil test cases with different geometries, Reynolds numbers, freestream turbulence conditions, and angles of attack. The test cases demonstrate the ability of the model to successfully reproduce transitional flow behavior with a reasonable degree of accuracy, particularly in comparison with commonly used models that exhibit no capability of predicting laminar-to-turbulent boundary layer development. While it is impossible to resolve all of the complex features of transitional and turbulent flows with a relatively simple Reynolds-averaged modeling approach, the results shown here demonstrate that the new model can provide a useful and practical tool for engineers addressing the simulation and prediction of transitional flow behavior in fluid systems.


2013 ◽  
Vol 117 (1188) ◽  
pp. 111-132 ◽  
Author(s):  
T. L. Grigorie ◽  
R. M. Botez

Abstract This paper presents a new adaptive algorithm for the statistical filtering of miniaturised inertial sensor noise. The algorithm uses the minimum variance method to perform a best estimate calculation of the accelerations or angular speeds on each of the three axes of an Inertial Measurement Unit (IMU) by using the information from some accelerometers and gyros arrays placed along the IMU axes. Also, the proposed algorithm allows the reduction of both components of the sensors’ noise (long term and short term) by using redundant linear configurations for the sensors dispositions. A numerical simulation is performed to illustrate how the algorithm works, using an accelerometer sensor model and a four-sensor array (unbiased and with different noise densities). Three cases of ideal input acceleration are considered: 1) a null signal; 2) a step signal with a no-null time step; and 3) a low frequency sinusoidal signal. To experimentally validate the proposed algorithm, some bench tests are performed. In this way, two sensors configurations are used: 1) one accelerometers array with four miniaturised sensors (n = 4); and 2) one accelerometers array with nine miniaturised sensors (n = 9). Each of the two configurations are tested for three cases of input accelerations: 0ms−1, 9·80655m/s2 and 9·80655m/s2.


Author(s):  
Karim Lahmer ◽  
Rachid Bessaïh ◽  
Angel Scipioni ◽  
Mohammed El Ganaoui

This paper summarizes numerical results of hydrogen absorption simulated in an axisymmetric tank geometry containing magnesium hydride heated to 300 °C and at moderate storage pressure 1 MPa. The governing equations are solved with a fully implicit finite volume numerical scheme used by a commercial software FLUENT. The effect of the different kinetic reaction equations modeling hydrogen absorption was studied by the introduction of a specific subroutine at each time step in order to consider which one will provide results close to available experimental results. Spatial and temporal profiles of temperature and concentration in hydride bed are plotted. Results show that suitable method for our two-dimensional study is a CV-2D technique because it generates the smallest error especially during the beginning of the reaction. Also, its computational time is the shortest one compared to the other methods.


Author(s):  
Shuangfeng Wang ◽  
Shigefumi Nishio

Heat transport rates of micro scale SEMOS (Self-Exciting Mode Oscillating) heat pipe with inner diameter of 1.5mm, 1.2mm and 0.9mm, were investigated by using R141b, ethanol and water as working fluids. The effects of inner diameter, liquid volume faction, and material properties of the working fluids are examined. It shows that the smaller the inner diameter, the higher the thermal transport density is. For removing high heat flux, the water is the most promising working fluid as it has the largest critical heat transfer rate and the widest operating range among the three kinds of working fluids. A one-dimensional numerical simulation is carried out to describe the heat transport characteristics and the two-phase flow behavior in the closed loop SEMOS heat pipe. The numerical prediction agrees with the experimental results fairly well, when the input heat through was not very high and the flow pattern was slug flow.   This paper was also originally published as part of the Proceedings of the ASME 2005 Pacific Rim Technical Conference and Exhibition on Integration and Packaging of MEMS, NEMS, and Electronic Systems.


Author(s):  
Takeshi Fujikawa ◽  
Etsujiro Imanishi

Abstract A method of time integration algorithm is presented for solving stiff vibration and motion problems. It is absolutely stable, numerically dissipative, and much accurate than other dissipative time integration methods. It achieves high-frequency dissipation, while minimizing unwanted low-frequency dissipation. In this method change of acceleration during time step is expressed as quadratic function including some parameters, whose appropriate values are determined through numerical investigation. Two calculation examples are demonstrated to show the usefulness of this method.


Author(s):  
Huiyu Wang ◽  
D. Keith Walters ◽  
Keisha B. Walters

Abstract This paper investigates the performance of a commercially available computational fluid dynamics (CFD) solver (Ansys FLUENT) to predict the flow and heat transfer characteristics of a two-phase closed thermosyphon (TPCT). Specifically, the study compares two different discretization schemes for the volume fraction equation with different time step sizes using three different sets of mass transfer time relaxation parameters for evaporation and condensation. The present study evaluates use of the Compressive scheme to increase the time step size compared to the Geo-Reconstruct scheme. In addition, a model is proposed to adjust the global saturation temperature of the system based on the volume of the vapor phase in order to balance the mass transfer inside TPCT and more accurately represent the realistic operating conditions of a TPCT. In this study a total of nineteen simulations are performed, and two types of boundary conditions for the condenser are investigated to determine the effect on the accuracy of the simulation results. The baseline simulation uses the Geo-Reconstruct method with a fixed saturation temperature. Other additional cases are performed using the Geo-Reconstruct method with variable saturation temperature, and the Compressive method with and without variable saturation temperature using different sets of mass transfer time relaxation parameters. Results show that the case using the Compressive method with the variable saturation temperature model has good agreement with the reference experimental data and is less computationally expensive than the Geo-Reconstruct method. The 3D CFD models implemented in this study successfully predict the phase change process and flow behavior inside the TPCTs, at least in a qualitative sense.


Author(s):  
Kang Shi ◽  
Xuhui He ◽  
Yunfeng Zou ◽  
Zhi Zheng

The dynamic interaction problem for the train–rail–bridge (TRB) systems presents a computational challenge, especially for the analysis of large-size TRB coupling systems. To address this issue, an efficient non-iterative hybrid method (NHM) is proposed. With this method, the integrated TRB system is divided into three subsystems, i.e. the train subsystem, the rail subsystem, and the bridge subsystem. Based on the individual subsystems, a multi-step[Formula: see text] technique is adopted in which a fine time step is used to analyze the high-frequency coupling vibration for the train and rail subsystems, and a coarse time step is adopted to calculate the low-frequency coupling vibration for the rail and bridge subsystem. Additionally, Zhais explicit integral method is used to predict the displacement of the wheelsets and the rail at the current time step before using the Newmark method. The proposed method incorporates the advantages of Zhai’s explicit method and the MS technique to avoid the iteration that may be required for the train–rail coupled analysis. The simulation fidelity and computational efficiency of the proposed method are demonstrated in the analysis of two examples of typical high-speed railway bridges. It was demonstrated that the proposed method can significantly enhance the computational efficiency, while maintaining a higher precision with a larger time step in comparison with other existing methods.


2015 ◽  
Vol 32 (2) ◽  
pp. 289-307 ◽  
Author(s):  
Yexiang Xiao ◽  
Zhengwei Wang ◽  
Jidi Zeng ◽  
jintai Zheng ◽  
Jiayang Lin ◽  
...  

Purpose – The purpose of this paper is to experimentally and numerically investigate the interference characteristics between two ski-jump jets on the flip bucket in a large dam spillway when two floodgates are running. Design/methodology/approach – The volume of fluid (VOF) method together with the Realizable k-ε turbulence model were used to predict the flow in two ski-jump jets and the free surface motion in a large dam spillway. The movements of the two gates were simulated using a dynamic mesh controlled by a User Defined Function (UDF). The simulations were run using the prototype dam as the field test to minimize errors due to scale effects. The simulation results are compared with field test observations. Findings – The transient flow calculations, accurately predict the two gate discharges compared to field data with the predicted ski-jump jet interference flow pattern similar to the observed shapes. The transient simulations indicate that the main reason for the deflected nappe is the larger opening difference between the two gates as the buttress side gate closes. When both gates are running, the two ski-jump jets interfere in the flip bucket and raise the jet nappe to near the buttress to form a secondary flow on this jet nappe surface. As the gate continues to close, the nappe surface continues to rise and the surface secondary flow become stronger, which deflects the nappe over the side buttress. Originality/value – A dynamic mesh is used to simulate the transient flow behavior of two prototype running gates. The transient flow simulation clarifies the hydraulics mechanism for how the two ski-jump jets interfere and deflect the nappe.


2016 ◽  
Vol 16 (01) ◽  
pp. 1640009 ◽  
Author(s):  
Chuanguo Jia ◽  
Zhou Leng ◽  
Yingmin Li ◽  
Hongliu Xia ◽  
Liping Liu

Systems of ordinary differential equations (ODEs) arising from transient structural dynamics very often exhibit high-frequency/low-frequency and linear/nonlinear behaviors of subsets of the state variables. With this in mind, the paper resorts to the use of different time integrators with different time steps for subsystems, which tailors each method and its time step to the solution behaviors of the corresponding subsystem. In detail, a partitioned integration method is introduced which imposes continuity of velocities at the interface to couple arbitrary Newmark schemes with different time steps in different subdomains. It is proved that the velocity continuity of the method is the primal factor of its reduction to first-order accuracy. To maintain second-order accuracy without increasing drift and computational cost, a novel method with the acceleration continuity is proposed whose velocity constraint is also ensured by means of the projection strategy. Both its stability and accuracy properties are examined through numerical analysis of a Single-degree-of-freedom (DoF) split mass system. Finally, numerical validations are conducted on Single- and Two-DoF split mass systems and a four-DoF nonlinear structure showing the feasibility of the proposed method.


Sign in / Sign up

Export Citation Format

Share Document