Multi-Dimensional Numerical Free Surface VOF Modeling With Moonpool Experiments

Author(s):  
Salman Sadiq ◽  
Xiong-Liang Yao

The homogeneous multiphase incompressible flow past a moonpool was investigated to determine the shape of vortex, flow field and pressure distributions. In this approach, a homogeneous flow model together with the VOF method for interface capturing is used to compute the entire flow field within the moonpool. The turbulence is represented via fully homogeneous buoyant model with air and water as continuous fluids. Numerical results are verified by conducting towing tank experiments. Simulation ship moonpool applications are verification of capabilities in ANSYS CFX multi-physics code by two and three dimensional circular and square shaped moonpool subjected to flow in a channel. The numerical results indicate that cavity location with a unit factor of 1.667(L/L1), where ‘L’ is length of cavity and ‘L1’ is the distance from inflow edge; was proved to be appropriate for such cavities to have optimum performance related to moonpool hydrodynamics. Distance ‘L2’ had a negligible effect on cavity. Free surface height ‘H’ inside the cavity was maintained at 0.4 m both in numerical calculations then in experiments. Numerical and experimental comparison of results reveals that due to multiphase modeling; results diverge from actual value near phase coupling. Shape factor like internal curvature, has a great influence on vortex shape and hydrodynamic forces inside the moonpool. By analyzing these numerical results a better understanding is established for VOF models and moonpool piston phenomenon. Square shaped moonpool results and numerical results agree very well till half of the moonpool depth from free surface height. These results can be very well used for the basis of designing complex shaped moonpool having free surface. The agreement with the experimental data is within the accuracy of other simulations. Further studies and the development of Best Practice Procedures are required as a next step to reduce/avoid numerical and set-up errors and establish CFD as an industrial tool.

Author(s):  
Fatemah Hoseini Dadmarzi ◽  
Hassan Ghassemi ◽  
Parviz Ghadimi ◽  
Babak Ommani

Fin stabilizers are very important device for controlling the ship roll motion against the external moments due to wave. This paper presents numerical results for flow field simulation and the hydrodynamic performance of fin stabilizer attached to a ship hull with free surface effects. Combination of CFD and RANS method has been used for this study. The fin is non-rectangular NACA0015 profile section with a finite aspect ratio. The numerical results include pressure distributions and flow field around the fin which are used to calculate lift coefficients and free surface elevation as the main interest. Some results are compared with available experimental and numerical data in literature and they show good agreement.


e-Polymers ◽  
2020 ◽  
Vol 20 (1) ◽  
pp. 724-732
Author(s):  
Changchun Ji ◽  
Yudong Wang

AbstractTo investigate the distribution characteristics of the three-dimensional flow field under the slot die, an online measurement of the airflow velocity was performed using a hot wire anemometer. The experimental results show that the air-slot end faces have a great influence on the airflow distribution in its vicinity. Compared with the air velocity in the center area, the velocity below the slot end face is much lower. The distribution characteristics of the three-dimensional flow field under the slot die would cause the fibers at different positions to bear inconsistent air force. The air velocity of the spinning centerline is higher than that around it, which is more conducive to fiber diameter attenuation. The violent fluctuation of the instantaneous velocity of the airflow could easily cause the meltblowing fiber to whip in the area close to the die.


Author(s):  
Huimin Tang ◽  
Shuaiqiang Liu ◽  
Hualing Luo

Profiled endwall is an effective method to improve aerodynamic performance of turbine. This approach has been widely studied in the past decade on many engines. When automatic design optimisation is considered, most of the researches are usually based on the assumption of a simplified simulation model without considering cooling and rim seal flows. However, many researchers find out that some of the benefits achieved by optimization procedure are lost when applying the high-fidelity geometry configuration. Previously, an optimization procedure has been implemented by integrating the in-house geometry manipulator, a commercial three-dimensional CFD flow solver and the optimization driver, IsightTM. This optimization procedure has been executed [12] to design profiled endwalls for a turbine cascade and a one-and-half stage axial turbine. Improvements of the turbine performance have been achieved. As the profiled endwall is applied to a high pressure turbine, the problems of cooling and rim seal flows should be addressed. In this work, the effects of rim seal flow and cooling on the flow field of two-stage high pressure turbine have been presented. Three optimization runs are performed to design the profiled endwall of Rotor-One with different optimization model to consider the effects of rim flow and cooling separately. It is found that the rim seal flow has a significant impact on the flow field. The cooling is able to change the operation condition greatly, but barely affects the secondary flow in the turbine. The influences of the profiled endwalls on the flow field in turbine and cavities have been analyzed in detail. A significant reduction of secondary flows and corresponding increase of performance are achieved when taking account of the rim flows into the optimization. The traditional optimization mechanism of profiled endwall is to reduce the cross passage gradient, which has great influence on the strength of the secondary flow. However, with considering the rim seal flows, the profiled endwall improves the turbine performance mainly by controlling the path of rim seal flow. Then the optimization procedure with consideration of rim seal flow has also been applied to the design of the profiled endwall for Stator Two.


1986 ◽  
Vol 108 (2) ◽  
pp. 285-292 ◽  
Author(s):  
W. Bra¨unling ◽  
F. Lehthaus

In a test facility for rotating annular cascades with three conical test sections of different taper angles (0, 30, 45 deg), experiments are conducted for two geometrically different turbine cascade configurations, a hub section cascade with high deflection and a tip section cascade with low deflection. The evaluation of time-averaged data derived from conventional probe measurements upstream and downstream of the test wheel in the machine-fixed absolute system is based on the assumption of axisymmetric stream surfaces. The cascade characteristics, i.e., mass flow, deflection, and losses, for a wide range of inlet flow angles and outlet Mach numbers are provided in the blade-fixed relative system with respect to the influence of annulus taper. Some of the results are compared with simple theoretical calculations. To obtain some information about the spatial structure of the flow within the cascade passages, surface pressure distributions on the profiles of the rotating test wheels are measured at three different radial blade sections. For some examples those distributions are compared with numerical results on plane cascades of the same sweep and dihedral angles and the same aspect ratios. The computer code used is based on a three-dimensional time-marching finite-volume method solving the Euler equations. Both experimental and numerical results show a fairly good qualitative agreement in the three-dimensional blade surface pressure distributions. This work will be continued with detailed investigations on the spatial flow structure.


1999 ◽  
Vol 396 ◽  
pp. 37-71 ◽  
Author(s):  
LEONID BREVDO ◽  
PATRICE LAURE ◽  
FREDERIC DIAS ◽  
THOMAS J. BRIDGES

The film flow down an inclined plane has several features that make it an interesting prototype for studying transition in a shear flow: the basic parallel state is an exact explicit solution of the Navier–Stokes equations; the experimentally observed transition of this flow shows many properties in common with boundary-layer transition; and it has a free surface, leading to more than one class of modes. In this paper, unstable wavepackets – associated with the full Navier–Stokes equations with viscous free-surface boundary conditions – are analysed by using the formalism of absolute and convective instabilities based on the exact Briggs collision criterion for multiple k-roots of D(k, ω) = 0; where k is a wavenumber, ω is a frequency and D(k, ω) is the dispersion relation function.The main results of this paper are threefold. First, we work with the full Navier–Stokes equations with viscous free-surface boundary conditions, rather than a model partial differential equation, and, guided by experiments, explore a large region of the parameter space to see if absolute instability – as predicted by some model equations – is possible. Secondly, our numerical results find only convective instability, in complete agreement with experiments. Thirdly, we find a curious saddle-point bifurcation which affects dramatically the interpretation of the convective instability. This is the first finding of this type of bifurcation in a fluids problem and it may have implications for the analysis of wavepackets in other flows, in particular for three-dimensional instabilities. The numerical results of the wavepacket analysis compare well with the available experimental data, confirming the importance of convective instability for this problem.The numerical results on the position of a dominant saddle point obtained by using the exact collision criterion are also compared to the results based on a steepest-descent method coupled with a continuation procedure for tracking convective instability that until now was considered as reliable. While for two-dimensional instabilities a numerical implementation of the collision criterion is readily available, the only existing numerical procedure for studying three-dimensional wavepackets is based on the tracking technique. For the present flow, the comparison shows a failure of the tracking treatment to recover a subinterval of the interval of unstable ray velocities V whose length constitutes 29% of the length of the entire unstable interval of V. The failure occurs due to a bifurcation of the saddle point, where V is a bifurcation parameter. We argue that this bifurcation of unstable ray velocities should be observable in experiments because of the abrupt increase by a factor of about 5.3 of the wavelength across the wavepacket associated with the appearance of the bifurcating branch. Further implications for experiments including the effect on spatial amplification rate are also discussed.


2012 ◽  
Vol 588-589 ◽  
pp. 1355-1358
Author(s):  
Xiao Xing ◽  
Guo Ming Ye

During the splicing process of pneumatic splicer, the principle of yarn splicing is closely related to the flow field inside the splicing chamber. This paper presents a numerical simulation of the flow char-acteristics inside the splicing chamber of the pneumatic splicer. A three-dimensional grid and the realizable tur¬bulence model are used in this simulation. The numerical results of veloc¬ity vectors distribution inside the chamber are shown. Streamlines starting from the two air injectors are also acquired. Based on the simulation, the principle of yarn splicing of the pneumatic splicer is discussed. The airflow in the splicing chamber can be divided into three regions. In addition, the simulation results have well sup¬ported the principle of yarn splicing of pneumatic splicer claimed by the splicing chamber makers.


2018 ◽  
Vol 140 (3) ◽  
Author(s):  
James Schock ◽  
Jason Dahl

Two methods are investigated to simultaneously obtain both three-dimensional (3D) velocity field and free surface elevations (FSEs) measurements near a surface piercing foil, while limiting the equipment. The combined velocity field and FSE measurements are obtained specifically for the validation of numerical methods requiring simultaneous field data and free surface measurements for a slender body shape. Both methods use stereo particle image velocimetry (SPIV) to measure three component velocities in the flow field and both methods use an off the shelf digital camera with a laser intersection line to measure FSEs. The first method is performed using a vertical laser sheet oriented parallel to the foil chord line. Through repetition of experiments with repositioning of the laser, a statistical representation of the three-dimensional flow field and surface elevations is obtained. The second method orients the vertical laser sheet such that the foil chord line is orthogonal to the laser sheet. A single experiment is performed with this method to measure the three-dimensional three component (3D3C) flow field and free surface, assuming steady flow conditions, such that the time dimension is used to expand the flow field in 3D space. The two methods are compared using dynamic mode decomposition and found to be comparable in the primary mode. Utilizing these methods produces results that are acceptable for use in numerical methods verification, at a fraction of the capital and computing cost associated with two plane or tomographic particle image velocimetry (PIV).


2016 ◽  
Vol 51 (16) ◽  
pp. 2363-2375 ◽  
Author(s):  
Annegret Mallach ◽  
Frank Härtel ◽  
Frieder Heieck ◽  
Jan-Philipp Fuhr ◽  
Peter Middendorf ◽  
...  

Scope of the presented work is a detailed comparison of a macroscopic draping model with real fibre architecture on a complex non-crimp-fabric preform using a new robot-based optical measurement system. By means of a preliminary analytical process design approach, a preforming test centre is set up to manufacture dry non-crimp-fabric preforms. A variable blank holder setup is used to investigate the effect of different process parameters on the fibre architecture. The real fibre architecture of those preforms is captured by the optical measurement system, which generates a three-dimensional model containing information about the fibre orientation along the entire surface of the preform. The measured and calculated fiber orientations are then compared with the simulation results in a three-dimensional overlay file. The results show that the analytical approach is able to predict local hot spots with high shear angles on the preform. Macroscopic simulations show a higher sensitivity towards changes in blank holder pressure than reality and limit the approach to precisely predict fibre architecture parameters on complex geometries.


Author(s):  
M Eslami ◽  
M M Tavakol ◽  
E Goshtasbirad

The problem of flow field around multiple bluff bodies mounted on a surface is of great significance in different fields of engineering. In this study, a computational fluid dynamics (CFD) code is developed to calculate three-dimensional (3D) steady state laminar fluid flow around two cuboids of arbitrary size and configuration mounted on a surface in free stream conditions. This study presents the results for two cubes of the same size mounted on a surface in both inline and staggered arrangements. Streamlines are plotted for various combinations of the distance between the two cubes and Reynolds number. Moreover, the effects of different parameters on vortical structures, separation, and reattachment points are discussed. Also, velocity and pressure distributions are plotted in the wake region behind the two cubes. It is clearly shown that how the presence of the second cube changes the flow field and the vortical structures in comparison with the case of a single cube.


Author(s):  
Weihua Mo ◽  
Philip L.-F. Liu

AbstractIn this paper we validate a numerical model for-structure interaction by comparing numerical results with laboratory data. The numerical model is based on the Navier-Stokes(N-S) equations for an incompressible fluid. The N-S equations are solved by two-step projection finite volume scheme and the free surface displacements are tracked by the slender vertical piles. Numerical results are compared with the laboratory data and very good agreement is observed for the time history of free surface displacement, fluid particle velocity and force. The agreement for dynamic pressure on the cylinder is less satisfactory, which is primarily caused by instrument errors.


Sign in / Sign up

Export Citation Format

Share Document