scholarly journals WAVE PROPAGATION ON A FLUME: NUMERICAL SIMULATION

2015 ◽  
Vol 14 (1) ◽  
pp. 95
Author(s):  
J. M. P. Conde ◽  
P. T. S. Roberto ◽  
C. J. E. M. Fortes

This paper presents the numerical simulations done by using the waves2Foam, an OpenFOAM® library, to simulate the propagation of regular waves without breaking in a three-dimensional flume. The numerical code solves the unsteady Navier-Stokes equations and uses a Volume-of-Fluid (VoF) method to identify the free-surface. A regular incident wave with a 1.5s period and 0.1m wave height was considered. This is one of the conditions, from the wide range of wave flume tests conducted at the National Laboratory for Civil Engineering (LNEC), whose objective was to analyze the hydrodynamics of wave transformation and wave breaking for different incident conditions over a variable bathymetry. Comparisons are made between the numerical and the experimental results. These comparisons include time-series of wave-gauges records at several locations along the flume and the corresponding amplitude spectra; significant wave height and average period evolution along the flume; time-series of the velocity components at one section of the flume, measured at the middle of the water column; and hodograph representation of the velocity components, in the middle of the water column, in the xy, xz, and yz planes, along the flume. It was found that the numerical results obtained are close to the experimental data. The observed differences are attributable to numerical inaccuracies as well as the differences between the wave generation method in the numerical and experimental tests.

2014 ◽  
Vol 16 (5) ◽  
pp. 901-918 ◽  

<div> <p>Three-dimensional calculations were performed to simulate the flow around a cylindrical vegetation element using the Scale Adaptive Simulation (SAS) model; commonly, this is the first step of the modeling of the flow through multiple vegetation elements. SAS solves the Reynolds Averaged Navier-Stokes equations in stable flow regions, while in regions with unstable flow it goes unsteady producing a resolved turbulent spectrum after reducing eddy viscosity according to the locally resolved vortex size represented by the von Karman length scale. A finite volume numerical code was used for the spatial discretisation of the rectangular computational domain with stream-wise, cross-flow and vertical dimensions equal to 30D, 11D and 1D, respectively, which was resolved with unstructured grids. Calculations were compared with experiments and Large Eddy Simulations (LES). Predicted overall flow parameters and mean flow velocities exhibited a very satisfactory agreement with experiments and LES, while the agreement of predicted turbulent stresses was satisfactory. Calculations showed that SAS is an efficient and relatively fast turbulence modeling approach, especially in relevant practical problems, in which the very high accuracy that can be achieved by LES at the expense of large computational times is not required.</p> </div> <p>&nbsp;</p>


Author(s):  
Azita Soleymani ◽  
Eveliina Takasuo ◽  
Piroz Zamankhan ◽  
William Polashenski

Results are presented from a numerical study examining the flow of a viscous, incompressible fluid through random packing of nonoverlapping spheres at moderate Reynolds numbers (based on pore permeability and interstitial fluid velocity), spanning a wide range of flow conditions for porous media. By using a laminar model including inertial terms and assuming rough walls, numerical solutions of the Navier-Stokes equations in three-dimensional porous packed beds resulted in dimensionless pressure drops in excellent agreement with those reported in a previous study (Fand et al., 1987). This observation suggests that no transition to turbulence could occur in the range of Reynolds number studied. For flows in the Forchheimer regime, numerical results are presented of the lateral dispersivity of solute continuously injected into a three-dimensional bounded granular bed at moderate Peclet numbers. Lateral fluid dispersion coefficients are calculated by comparing the concentration profiles obtained from numerical and analytical methods. Comparing the present numerical results with data available in the literature, no evidence has been found to support the speculations by others for a transition from laminar to turbulent regimes in porous media at a critical Reynolds number.


Author(s):  
Jack Paterson ◽  
Philipp R Thies ◽  
Roman Sueur ◽  
Jérôme Lonchampt ◽  
Federico D’Amico

This article presents a metocean modelling methodology using a Markov-switching autoregressive model to produce stochastic wind speed and wave height time series, for inclusion in marine risk planning software tools. By generating a large number of stochastic weather series that resemble the variability in key metocean parameters, probabilistic outcomes can be obtained to predict the occurrence of weather windows, delays and subsequent operational durations for specific tasks or offshore construction phases. To cope with the variation in the offshore weather conditions at each project, it is vital that a stochastic weather model is adaptable to seasonal and inter-monthly fluctuations at each site, generating realistic time series to support weather risk assessments. A model selection process is presented for both weather parameters across three locations, and a personnel transfer task is used to contextualise a realistic weather window analysis. Summarising plots demonstrate the validity of the presented methodology and that a small extension improves the adaptability of the approach for sites with strong correlations between wind speed and wave height. It is concluded that the overall methodology can produce suitable wind speed and wave time series for the assessment of marine operations, yet it is recommended that the methodology is applied to other sites and operations, to determine the method’s adaptability to a wide range of offshore locations.


Author(s):  
Dominique Legendre ◽  
Catherine Colin ◽  
Typhaine Coquard

The three-dimensional flow around a hemispherical bubble sliding and growing on a wall in a viscous linear shear flow is studied numerically by solving the full Navier–Stokes equations in a boundary-fitted domain. The main goal of the present study is to provide a complete description of the forces experienced by the bubble (drag, lift and added mass) over a wide range of sliding and shear Reynolds numbers (0.01≤ Re b , Re α ≤2000) and shear rate (0≤ Sr ≤5). The drag and lift forces are computed successively for the following situations: an immobile bubble in a linear shear flow; a bubble sliding on the wall in a fluid at rest; and a bubble sliding in a linear shear flow. The added-mass force is studied by considering an unsteady motion relative to the wall or a time-dependent radius.


Author(s):  
Md. Readul Mahmud

The fluids inside passive micromixers are laminar in nature and mixing depends primarily on diffusion. Hence mixing efficiency is generally low, and requires a long channel length and longtime compare to active mixers. Various designs of complex channel structures with/without obstacles and three-dimensional geometries have been investigated in the past to obtain an efficient mixing in passive mixers. This work presents a design of a modified T mixer. To enhance the mixing performance, circular and hexagonal obstacles are introduced inside the modified T mixer. Numerical investigation on mixing and flow characteristics in microchannels is carried out using the computational fluid dynamics (CFD) software ANSYS 15. Mixing in the channels has been analyzed by using Navier–Stokes equations with water-water for a wide range of the Reynolds numbers from 1 to 500. The results show that the modified T mixer with circular obstacles has far better mixing performance than the modified T mixer without obstacles. The reason is that fluids' path length becomes longer due to the presence of obstacles which gives fluids more time to diffuse. For all cases, the modified T mixer with circular obstacle yields the best mixing efficiency (more than 60%) at all examined Reynolds numbers. It is also clear that efficiency increase with axial length. Efficiency can be simply improved by adding extra mixing units to provide adequate mixing. The value of the pressure drop is the lowest for the modified T mixer because there is no obstacle inside the channel. Modified T mixer and modified T mixer with circular obstacle have the lowest and highest mixing cost, respectively. Therefore, the current design of modified T with circular obstacles can act as an effective and simple passive mixing device for various micromixing applications.


2021 ◽  
Vol 15 ◽  
pp. 16-21
Author(s):  
D. G. Baek ◽  
J. H. Jung ◽  
H. S. Yoon

This study numerically carried out the propeller open water test (POW) by solving Navier-Stokes equations governing the three-dimensional unsteady incompressible viscous flow with the turbulence closure model of the Κ-ω SST model. Numerical simulations were performed at wide range of advance ratios. A great difference of velocity magnitude between the inner region and the outer region of the slipstream tube forms the thick and large velocity gradient which originates from the propeller tip and develops along the downstream. Eventually, the strong shear layer appears and plays the role of the slipstream boundary. As the advance ratio increases, the vortical structures originated from the propeller tips quickly decay. The contraction of the vortices trace is considerable with decreasing the advance ratio.


Author(s):  
Kuki Junichi ◽  
Kazuyuki Toda ◽  
Makoto Yamamoto

This paper presents a numerical procedure to predict a three-dimensional sand erosion phenomenon and the interaction between the flow field and the eroded surface. To simulate this phenomenon, the turbulent flow field, the particle trajectory and the amount of erosion on the eroded wall are calculated repeatedly. In computations of the flow field, compressible Navier-Stokes equations and low-Reynolds-number type k–ε turbulence model are adopted. Assuming that the concentration of suspended particle is dilute, particle-particle collision and the influence of particle motions on the flow field are neglected. The Neilson-Gilchrist erosion model is used to estimate the weight loss due to erosion. To verify the developed code, two types of 90-degree bends are computed. The results show that the present procedure can reasonably reproduce the sand erosion process and the temporal change of both the flow field and the wall surface qualitatively.


Author(s):  
F. Giannetti ◽  
V. Narcisi ◽  
A. Subioli ◽  
A. Del Nevo

Before the definitive shutdown of Phénix fast reactor, a pool-type sodium-cooled prototype, occurred in the year 2009, the French Commissariat a l’Energie Atomique et aux Energies Alternatives (CEA) decided to carry out a final set of experimental tests, in order to gather data and additional knowledge on relevant Sodium Fast Reactors (SFR) operation and safety aspects. One of these experiments was the dissymmetrical configuration test, which has been selected as a benchmark transient on H2020 SESAME project. ENEA and Sapienza University of Rome are participating, in collaboration, in the benchmark using RELAP5-3D© code, developed by Idaho National Laboratory (INL). The analysis of the thermal hydraulic characteristics and the plant behavior focuses on the accurate core cooling prediction in natural circulation, by assessing the flow and the temperatures in the various reactor components. Aiming to investigate the asymmetric thermal hydraulic behavior inside the pool, the nodalization scheme has been developed using a detailed parallel pipe with cross-junctions for the hot and cold pool and the core has been simulated assembly per assembly, including blanket and reflector, in the first seven rings, ready for a proficient coupling with a three-dimensional neutron kinetic analysis code. The scope of the present study is to verify if a detailed thermal–hydraulics pool model could add useful data to show the differences which occur in each loop and into the pools, caused by the azimuthal asymmetry of the boundary conditions. Blind calculations results are here presented and discussed, highlighting the details added with the multidimensional nodalization scheme, compared to a mono-dimensional model. The paper will be a first step on the comparison of a Phénix mono-dimensional and three-dimensional thermal-hydraulic nodalization and on the RELAP5-3D© code assessment against the dissymmetrical experimental results.


2020 ◽  
Author(s):  
HaKun Jang ◽  
Celalettin Emre Ozdemir ◽  
Mayank Tyagi ◽  
Jun-Hong Liang

Abstract The purpose of this study is to numerically investigate the bed shear stress and near-bed mixing due to coherent vortex structures in the vicinity of a vertically wall-mounted circular cylinder subject to an imposed finite-depth oscillatory sinusoidal flow. Previous studies reveal that the Keulegan–Carpenter (KC) number influences the formation of lee-side wake vortex structures as well as the horseshoe vortex in front of a cylinder. Therefore, parametric studies in a moderately wide range of KC from 5 to 20 are numerically performed. In the present study, Direct Numerical Simulation (DNS) is conducted using the open-source software, OpenFOAM, that solves the three-dimensional unsteady incompressible Navier-Stokes equations using finite volume method. Nondimensional parameters used in the simulations are carefully chosen to represent the real physics. The numerical solutions are first validated using an analytical solution for the oscillating Stokes flow and the results are then systematically and quantitatively compared with the experimental measurements. The results show that the lee-side wake is significantly influenced by KC, and distinctive types of the lee-side wake are generated and classified based on KC. It is also found that both KC and the ratio of the thickness of the Stokes boundary layer to the water depth are heavily associated with the stability of the lee-side wake. In addition, the simulated size and lifespan of the horseshoe vortex agree well with the experimental data.


2011 ◽  
Vol 133 (4) ◽  
Author(s):  
C. I. Papadopoulos ◽  
E. E. Efstathiou ◽  
P. G. Nikolakopoulos ◽  
L. Kaiktsis

This paper presents an optimization study of the geometry of three-dimensional micro-thrust bearings in a wide range of convergence ratios. The optimization goal is the maximization of the bearing load carrying capacity. The bearings are modeled as micro-channels, consisting of a smooth moving wall (rotor), and a stationary wall (stator) with partial periodic rectangular texturing. The flow field is calculated from the numerical solution of the Navier-Stokes equations for incompressible isothermal flow; processing of the results yields the bearing load capacity and friction coefficient. The geometry of the textured channel is defined parametrically for several width-to-length ratios. Optimal texturing geometries are obtained by utilizing an optimization tool based on genetic algorithms, which is coupled to the CFD code. Here, the design variables define the bearing geometry and convergence ratio. To minimize the computational cost, a multi-objective approach is proposed, consisting in the simultaneous maximization of the load carrying capacity and minimization of the bearing convergence ratio. The optimal solutions, identified based on the concept of Pareto dominance, are equivalent to those of single-objective optimization problems for different convergence ratio values. The present results demonstrate that the characteristics of the optimal texturing patterns depend strongly on both the convergence ratio and the width-to-length ratio. Further, the optimal load carrying capacity increases at increasing convergence ratio, up to an optimal value, identified by the optimization procedure. Finally, proper surface texturing provides substantial load carrying capacity even for parallel or slightly diverging bearings. Based on the present results, we propose simple formulas for the design of textured micro-thrust bearings.


Sign in / Sign up

Export Citation Format

Share Document