Search results for: numerical groundwater flow modeling

Authors: Murat Avci, Fatih Kosar, Ismail Yilmaz

Abstract:

In this study, structure-flow interaction for water brake mechanism is studied with Abaqus CEL approach. The water brake mechanism is used for dynamic systems such as sled system on rail. For the achievement of these system tests, structure-flow interaction should be investigated in detail. This study is about a sled test of an aircraft subsystem which rises to supersonic speeds thanks to rocket engines. To decrease or to stop the thrusting rocket sleds, water brake mechanisms are used. Water brake mechanism provides the deceleration of the structures that have supersonic speeds. Therefore, structure-flow interaction may cause damage to the water brake mechanism. To verify all design revisions with system tests are so costly so that some decisions are taken in accordance with numerical methods. In this study, structure-flow interaction that belongs to water brake mechanism is solved with Abaqus CEL approach. Fluid and deformation on the structure behaviors are modeled at the same time thanks to CEL approach. Provided analysis results are corrected with the dynamic tests. Deformation zones seen in numerical analysis are also observed in dynamic tests. Finally, Johnson-Cook material model parameters used for this analysis are proven, and it is understood that these parameters can be used for dynamic analysis like water brake mechanism.

Keywords: aircraft, rocket, structure-flow, supersonic

Procedia PDF Downloads 135

Authors: Kiranmay Sarma

Abstract:

Recent decades have witnessed rapid increase in urbanization, for which, rural-urban migration is stated to be the principal reason. Urban growth throughout the world has already outstripped the capacities of most of the cities to provide basic amenities to the citizens, including clean drinking water and consequently, they are struggling to get fresh and clean water to meet water demands. Delhi, the capital of India, is one of the rapid fast growing metropolitan cities of the country. As a result, there has been large influx of population during the last few decades and pressure exerted to the limited available water resources, mainly on groundwater. Considering this important aspect, the present research has been designed to study the effects of urbanization on the green cover and groundwater and their relationships of Delhi. For the purpose, four different land uses of the study area have been considered, viz., protected forest area, trees outside forest, maintained park and settlement area. Samples for groundwater and vegetation were collected seasonally in post-monsoon (October), winter (February) and summer (June) at each study site for two years during 2012 and 2014. The results were integrated into GIS platform. The spatial distribution of groundwater showed that the concentration of most of the ions is decreasing from northern to southern parts of Delhi, thus groundwater shows an improving trend from north to south. The depth was found to be improving from south to north Delhi, i.e., opposite to the water quality. The study concludes the groundwater properties in Delhi vary spatially with depending on the types of land cover.

Keywords: groundwater, urbanization, GIS, green cover, Delhi

Procedia PDF Downloads 269

Authors: Ridha Salhi, Mondher Labidi, Fethi Choubani

Abstract:

Metamaterials have attracted much attention in recent years because of their electromagnetic exquisite proprieties. We will present, in this paper, the modeling of three metamaterial structures by equivalent circuit model. We begin by modeling the SRR (Split Ring Resonator), then we model the HIS (High Impedance Surfaces), and finally, we present the model of the CPW (Coplanar Wave Guide). In order to validate models, we compare the results obtained by an equivalent circuit models with numerical simulation.

Keywords: metamaterials, SRR, HIS, CPW, IDC

Procedia PDF Downloads 406

Authors: Shidvash Vakilipour, Scott Ormiston, Masoud Mohammadi, Rouzbeh Riazi, Kimia Amiri, Sahar Barati

Abstract:

Two-phase and multi-phase flows are common flow types in fluid mechanics engineering. Among the basic and applied problems of these flow types, two-phase parallel flow is the one that two immiscible fluids flow in the vicinity of each other. In this type of flow, fluid properties (e.g. density, viscosity, and temperature) are different at the two sides of the interface of the two fluids. The most challenging part of the numerical simulation of two-phase flow is to determine the location of interface accurately. In the present work, a coupled interface tracking algorithm is developed based on Arbitrary Lagrangian-Eulerian (ALE) approach using a cell-centered, pressure-based, coupled solver. To validate this algorithm, an analytical solution for fully developed two-phase flow in presence of gravity is derived, and then, the results of the numerical simulation of this flow are compared with analytical solution at various flow conditions. The results of the simulations show good accuracy of the algorithm despite using a nearly coarse and uniform grid. Temporal variations of interface profile toward the steady-state solution show that a greater difference between fluids properties (especially dynamic viscosity) will result in larger traveling waves. Gravity effect studies also show that favorable gravity will result in a reduction of heavier fluid thickness and adverse gravity leads to increasing it with respect to the zero gravity condition. However, the magnitude of variation in favorable gravity is much more than adverse gravity.

Keywords: coupled solver, gravitational force, interface tracking, Reynolds number to Froude number, two-phase flow

Procedia PDF Downloads 294

Authors: Y. Lasbet, A. L. Boukhalkhal, K. Loubar

Abstract:

The study of mixed convection is, usually, focused on the straight channels in which the onset of the mixed convection is well defined as function of the ratio between Grashof number and Reynolds number, Gr/Re. This is not the case for a complex channel wherein the mixed convection is not sufficiently examined in the literature. Our paper focuses on the study of the mixed convection in a complex geometry in which our main contribution reveals that the critical value of the ratio Gr/Re for the onset of the mixed convection increases highly in the type of geometry contrary to the straight channel. Furthermore, the accentuated secondary flow in this geometry prevents the thermal stratification in the flow and consequently the buoyancy driven becomes negligible. To perform these objectives, a numerical study in complex geometry for several values of the ratio Gr/Re with prescribed wall heat flux (H2), was realized by using the CFD code.

Keywords: complex geometry, heat transfer, laminar flow, mixed convection, Nusselt number

Procedia PDF Downloads 478

Authors: Yagya Murti Aryal, Bipin Adhikari, Pradeep Gyawali

Abstract:

The Lesser Himalaya zone of Nepal consists of thrusting and folding belts, which play an important role in the sustainable management of groundwater in the Himalayan regions. The study area is located in the Dolakha and Ramechhap Districts of Bagmati Province, Nepal. Geologically, these districts are situated in the Lesser Himalayas and partly encompass the Higher Himalayan rock sequence, which includes low-grade to high-grade metamorphic rocks. Following the Gorkha Earthquake in 2015, numerous springs dried up, and many others are currently experiencing depletion due to the distortion of the natural groundwater flow. The primary objective of this study is to identify potential groundwater areas and determine suitable sites for artificial groundwater recharge. Two distinct statistical approaches were used to develop models: The Frequency Ratio (FR) and Shannon Entropy (SE) methods. The study utilized both primary and secondary datasets and incorporated significant role and controlling factors derived from field works and literature reviews. Field data collection involved spring inventory, soil analysis, lithology assessment, and hydro-geomorphology study. Additionally, slope, aspect, drainage density, and lineament density were extracted from a Digital Elevation Model (DEM) using GIS and transformed into thematic layers. For training and validation, 114 springs were divided into a 70/30 ratio, with an equal number of non-spring pixels. After assigning weights to each class based on the two proposed models, a groundwater potential map was generated using GIS, classifying the area into five levels: very low, low, moderate, high, and very high. The model's outcome reveals that over 41% of the area falls into the low and very low potential categories, while only 30% of the area demonstrates a high probability of groundwater potential. To evaluate model performance, accuracy was assessed using the Area under the Curve (AUC). The success rate AUC values for the FR and SE methods were determined to be 78.73% and 77.09%, respectively. Additionally, the prediction rate AUC values for the FR and SE methods were calculated as 76.31% and 74.08%. The results indicate that the FR model exhibits greater prediction capability compared to the SE model in this case study.

Keywords: groundwater potential mapping, frequency ratio, Shannon’s Entropy, Lesser Himalaya Zone, sustainable groundwater management

Procedia PDF Downloads 55

Authors: Y. Galerkin, K. Soldatova, A. Drozdov

Abstract:

The 6th version of Universal modeling method for centrifugal compressor stage calculation is described. Identification of the new mathematical model was made. As a result of identification the uniform set of empirical coefficients is received. The efficiency definition error is 0,86 % at a design point. The efficiency definition error at five flow rate points (except a point of the maximum flow rate) is 1,22 %. Several variants of the stage with 3D impellers designed by 6th version program and quasi three-dimensional calculation programs were compared by their gas dynamic performances CFD (NUMECA FINE TURBO). Performance comparison demonstrated general principles of design validity and leads to some design recommendations.

Keywords: compressor design, loss model, performance prediction, test data, model stages, flow rate coefficient, work coefficient

Procedia PDF Downloads 400

Authors: R. Chouatah, E. G. Filali, B. Zouzou

Abstract:

It has been well established that there are no differences between microscale and macroscale flows of incompressible liquids. However, surface roughness has been known to impact the transport phenomena. The effect of structured roughness on the dynamics and heat transfer of water flowing through minichannel was numerically investigated in this study. Our study consists in characterizing the dynamic field and heat transfer aspect of a flow in circular minichannel equipped with structured roughness using CFD software, CFX. The study is performed to understand the effect of various roughness elements (rectangular, triangular), roughness height and roughness pitch on the friction factor and heat transfer coefficient. Our work focuses on a water flow inside a circular mini-channel of 1 mm in diameter and 10 cm in length. The speed entry into the mini-channel varies from 0.1 m/s to 25 m/s. The wall of the mini-channel is submitted to a constant heat flux; q=100,000 W/m². The simulations results are compared to those obtained with smooth minichannel and the existing experimental and numerical results in the literature.

Keywords: heat transfer, laminar and turbulent flow, minichannel, structured roughness

Procedia PDF Downloads 324

Authors: Dmitry V. Fomichev, Vladimir V. Solonin

Abstract:

This paper presents the findings from a numerical simulation of the flow in 37-rod fuel assembly models spaced by a double-wire trapezoidal wrapping as applied to the BREST-OD-300 experimental nuclear reactor. Data on a high static pressure distribution within the models, and equations for determining the fuel bundle flow friction factors have been obtained. Recommendations are provided on using the closing turbulence models available in the ANSYS Fluent. A comparative analysis has been performed against the existing empirical equations for determining the flow friction factors. The calculated and experimental data fit has been shown. An analysis into the experimental data and results of the numerical simulation of the BREST-OD-300 fuel rod assembly hydrodynamic performance are presented.

Keywords: BREST-OD-300, ware-spaces, fuel assembly, computation fluid dynamics

Procedia PDF Downloads 361

Authors: Nader Pourmahmoud, Amir Hassanzadeh

Abstract:

This article investigates the effects of inlet pressure in a newly introduced vortex tube which has been equipped with an additional vortex chamber. A 3-D compressible turbulent flow computation has been carried out toward analysis of complex flow field in this apparatus. Numerical results of flows are derived by utilizing the standard k-ε turbulence model for analyzing high rotating complex flow field. The present research has focused on cooling effect and given a characteristics curve for minimum cool temperature. In addition, the effect of inlet pressure for both chambers has been studied in details. To be presented numerical results show that the effect of inlet pressure in second chamber has more important role in improving the performance of the vortex tube than first one. By increasing the pressure in the second chamber, cold outlet temperature reaches a higher decrease. When both chambers are fed with high pressure fluid, best operation condition of vortex tube occurs. However, it is not possible to feed both chambers with high pressure due to the conditions of working environment.

Keywords: energy separation, inlet pressure, numerical simulation, vortex chamber, vortex tube

Procedia PDF Downloads 352

Authors: Arash Mir Abdolah Lavasani, M. Ebrahimisabet

Abstract:

In this paper, the 2-D unsteady viscous flow around two cam shaped cylinders in tandem arrangement is numerically simulated in order to study the characteristics of the flow in turbulent regimes. The investigation covers the effects of high subcritical and supercritical Reynolds numbers and L/D ratio on total drag coefficient. The equivalent diameter of cylinders is 27.6 mm The space between center to center of two cam shaped cylinders is define as longitudinal pitch ratio and it varies in range of 1.5 < L/D < 6. Reynolds number base on equivalent circular cylinder varies in range of 27×103 < Re < 166×103 Results show that drag coefficient of both cylinders depends on pitch ratio. However drag coefficient of downstream cylinder is more dependent on the pitch ratio.

Keywords: cam shaped, tandem, numerical, drag coefficient, turbulent

Procedia PDF Downloads 447

Authors: S. M. Moustafa, A. Ezzat, Y. S. Taha, G. H. Hassib, S. Hamada

Abstract:

After the occurrence of 14, Nov, 1981 earthquake (M = 5.3), on Kalabska fault, south of Egypt, seismic stations distributed in and around the Kalabsha area, in order to monitoring, recording and studying the seismic activity in the area. In addition of that, from 1989 a number of piezometer wells drilled in the same area, distribed on at the both side of the active faults area and in different water bearing formations, in order to measuring the groundwater parameters (level, temperature, ph, and conductivity) to monitoring the relationship between those parameters and the seismic activity at Kalabsha area. The behavior of groundwater due to seismic activity over the world studied by several scientists i.e. H. Wakita (1979) on Izu-Oshima earthquake (M= 7.0) at Japan, M. E. Contadakis & G.asteriadis (1972), and Evans (1966), they found an anomalies on groundwater measurements prior, co, and post the occurrence of bigger earthquakes, referring to the probability of precursory evidence of impending earthquakes. In Kalabsha area south of Egypt, this study has been done using recorded seismic data, and the measurements of underground water parameters. same phenomena of anomalies founded on groundwater measurements pre, co. and post the occurrence of earthquakes with magnitude bigger than 3, and no systematic regularity exists for epicenter distance, duration of anomalies or time lag between anomalies appear and occurrence of events. Also the results found present strong relation between the groundwater in the upper unconfined aquifer Nubian Sandstone formation, and Kalabsha seismic activity, otherwise no relation between the seismic activities in the area with the deep groundwater in the lower confined aquifer Sandstone.

Keywords: seismicity, groundwater, Aswan, Egypt

Procedia PDF Downloads 364

Authors: Ainul Haque, Ameeye Kumar Nayak

Abstract:

Flow enhancement and species transport in a slit hydrophobic microchannel is studied for non-Newtonian fluids with the externally imposed electric field and pressure gradient. The incompressible Poisson-Nernst-Plank equations and the Navier-Stokes equations are approximated by lubrication theory to quantify the flow structure due to hydrophobic and hydrophilic surfaces. The analytical quantification of velocity and pressure of electroosmotic flow (EOF) is made with the numerical results due to the staggered grid based finite volume method for flow governing equations. The resistance force due to fluid friction and shear force along the surface are decreased by the hydrophobicity, enables the faster movement of fluid particles. The resulting flow enhancement factor Ef is increased with the low viscous fluid and provides maximum species transport. Also, the analytical comparison of EOF with pressure driven EOF justifies the flow enhancement due to hydrophobicity and shear impact on flow variation.

Keywords: electroosmotic flow, hydrophobic surface, power-law fluid, shear effect

Procedia PDF Downloads 350

Authors: R. Faiez, M. Mashhoudi, F. Najafi

Abstract:

Implicit in most large-scale numerical analyses of the crystal growth from the melt is the assumption that the shape and position of the phase boundary are determined by the transport phenomena coupled strongly to the melt hydrodynamics. In the present numerical study, the interface shape-effect on the convective interactions in a Czochralski oxide melt is described. It was demonstrated that thermos-capillary flow affects inversely the phase boundaries of distinct shapes. The in homogenity of heat flux and the location of the stagnation point at the crystallization front were investigated. The forced convection effect on the point displacement at the boundary found to be much stronger for the flat plate interface compared to the cone-shaped one with and without the Marangoni flow.

Keywords: computer simulation, fluid flow, interface shape, thermos-capillary effect

Procedia PDF Downloads 226

Authors: Tae Gwan Kim, Hyun Kyu Cho, Young Hoon Lee, Young Chul Park

Abstract:

Foam monitors are devices that are installed on cargo tank decks to suppress cargo area fires in oil tankers or hazardous chemical ship cargo ships. In general, the main design parameter of the foam monitor is the distance of the projection through the foam monitor. In this study, the relationship between flow characteristics and projection distance, depending on the shape was examined. Numerical techniques for fluid analysis of foam monitors have been developed for prediction. The flow pattern of the fluid varies depending on the shape of the flow path of the foam monitor, as the flow losses affecting projection distance were calculated through numerical analysis. The basic shape of the foam monitor was an L shape designed by N Company. The modified model increased the length of the flow path and used the S shape model. The calculation result shows that the L shape, which is the basic shape, has a problem that the force is directed to one side and the vibration and noise are generated there. In order to solve the problem, S-shaped model, which is a change model, was used. As a result, the problem is solved, and the projection distance from the nozzle is improved.

Keywords: CFD, foam monitor, projection distance, moment

Procedia PDF Downloads 321

Authors: Sumei Cai, Hong Li

Abstract:

Either during water or gas injection into reservoir, in order to understand the areal flow pressure distribution underground, associated bounding deformation is prevalently monitored by ground or downhole tiltmeters. In this paper, an inverse solution to elastic response of far field displacements induced by reservoir pressure change due to flow injection was studied. Furthermore, the fundamental theory on inverse solution to elastic problem as well as its spatial smoothing approach is presented. Taking advantage of source code development based on Boundary Element Method, numerical analysis on the monitoring data of ground surface displacements to further understand the behavior of reservoir process was developed. Numerical examples were also conducted to verify the effectiveness.

Keywords: remote displacement, inverse problem, boundary element method, BEM, reservoir process

Procedia PDF Downloads 100

Authors: Kumaresh Selvakumar, Man Young Kim

Abstract:

In selective catalytic reduction (SCR) unit, the injection system is provided with unique dosing pump to govern the urea injection phenomenon. The urea based operating liquid from the AdBlue tank links up directly with the dosing pump unit to furnish appropriate high pressure for examining the flow characteristics inside the liquid pump. This work aims in demonstrating the importance of external gear pump to provide pertinent high pressure and respective mass flow rate for each rotation. Numerical simulations are conducted using immersed solid method technique for better understanding of unsteady flow characteristics within the pump. Parametric analyses have been carried out for the gear speed and mass flow rate to find the behavior of pressure fluctuations. In the simulation results, the outlet pressure achieves maximum magnitude with the increase in rotational speed and the fluctuations grow higher.

Keywords: AdBlue tank, external gear pump, immersed solid method, selective catalytic reduction

Procedia PDF Downloads 253

Authors: D. Srinivasacharya, G. Madhava Rao

Abstract:

In this article, the pulsatile flow of blood flow in bifurcated artery with mild stenosis is investigated. Blood is treated to be a micropolar fluid with constant density. The arteries forming bifurcation are assumed to be symmetric about its axes and straight cylinders of restricted length. As the geometry of the stenosed bifurcated artery is irregular, it is changed to regular geometry utilizing the appropriate transformations. The numerical solutions, using the finite difference method, are computed for the flow rate, the shear stress, and the impedance. The influence of time, coupling number, half of the bifurcated angle and Womersley number on shear stress, flow rate and impedance (resistance to the flow) on both sides of the flow divider is shown graphically. It has been observed that the shear stress and flow rate are increasing with increase in the values of Womersley number and bifurcation angle on both sides of the apex. The shear stress is increasing along the inner wall and decreasing along the outer wall of the daughter artery with an increase in the value of coupling number. Further, it has been noticed that the shear stress, flow rate, and impedance are perturbed largely near to the apex in the parent artery due to the presence of backflow near the apex.

Keywords: micropolar fluid, bifurcated artery, stenosis, back flow, secondary flow, pulsatile flow, Womersley number

Procedia PDF Downloads 179

Authors: Ogan Karabas, Suleyman Yigit

Abstract:

Number of blade and blade type of centrifugal fans are the most decisive factor on the field of application, noise level, suction pressure and outlet mass flow rate. Nowadays, in order to determine these effects on centrifugal fans, numerical studies are carried out in addition to experimental studies. In this study, it is aimed to numerically investigate the changes of suction pressure and outlet mass flow rate values of a centrifugal fan according to the number of blade and blade type. Centrifugal fans of the same size with forward, backward and straight blade type were analyzed by using a simulation program and compared with each other. This analysis was carried out under steady state condition by selecting k-Ɛ turbulence model and air is assumed incompressible. Then, 16, 32 and 48 blade centrifugal fans were again analyzed by using same simulation program, and the optimum number of blades was determined for the suction pressure and the outlet mass flow rate. According to the results of the analysis, it was obtained that the suction pressure in the 32 blade fan was twice the value obtained in the 16 blade fan. In addition, the outlet mass flow rate increased by 45% with the increase in the number of blade from 16 to 32. There is no significant change observed on the suction pressure and outlet mass flow rate when the number of blades increased from 32 to 48. In the light of the analysis results, the optimum blade number was determined as 32.

Keywords: blade type, centrifugal fan, cfd, outlet mass flow rate, suction pressure

Procedia PDF Downloads 385

Authors: Hamed K. Arzani, Hamid K. Arzani, S.N. Kazi, A. Badarudin

Abstract:

Numerical investigation into convective heat transfer of CuO-Water based nanofluid in a pipe with return bend under laminar flow conditions has been done. The impacts of Reynolds number and the volume concentration of nanoparticles on the flow and the convective heat transfer behaviour are investigated. The results indicate that the increase in Reynolds number leads to the enhancement of average Nusselt number, and the increase in specific heat in the presence of the nanofluid results in improvement in heat transfer. Also, the presence of the secondary flow in the curve plays a key role in increasing the average Nusselt number and it appears higher than the inlet and outlet tubes. However, the pressure drop curve increases significantly in the tubes with the increase in nanoparticles concentration.

Keywords: laminar forced convection, curve pipe, return bend, nanufluid, CFD

Procedia PDF Downloads 283

Authors: Seyed Abolhasan Naeini, Eisa Aliagahei

Abstract:

Today, the use of geosynthetic materials in geotechnical activities is increasing significantly. One of the main uses of these materials is to increase the compressive strength of clay reinforced by geotextile layers. In the present study, the effect of clay reinforcement by geotextile layers in increasing the compressive strength of clay has been investigated using modeling in ABAQUS 6.11.3 software. For this purpose, the modified Drager Prager model has been chosen to simulate the stress-strain behavior of soil layers and the linear elastic model for the geotextile layer. Unreinforced samples and reinforced samples are modeled by geotextile layers (1, 2 and 3 geotextile layers) by software. In order to validate the results, an article in the same field was used and the numerical modeling results were calibrated with the laboratory results. Based on the obtained results, the software has a suitable capability for modeling and the results of the numerical model overlap with the laboratory results to a very acceptable extent, by increasing the number of geotextile layers, the error between the results of the laboratory sample and the software model increases. The highest amount of error is related to the sample reinforced with three layers of geotextile and is 7.3%.

Keywords: Abaqus, cap model, clay, geotextile layer, reinforced soil

Procedia PDF Downloads 71

Authors: R. Haoui

Abstract:

The purpose of this work is to simulate the flow at the exit of Vulcan 1 engine of European launcher Ariane 5. The geometry of the propellant nozzle is already determined using the characteristics method. The pressure in the outlet section of the nozzle is less than atmospheric pressure on the ground, causing the existence of oblique and normal shock waves at the exit. During the rise of the launcher, the atmospheric pressure decreases and the shock wave disappears. The code allows the capture of shock wave at exit of nozzle. The numerical technique uses the Flux Vector Splitting method of Van Leer to ensure convergence and avoid the calculation instabilities. The Courant, Friedrichs and Lewy coefficient (CFL) and mesh size level are selected to ensure the numerical convergence. The nonlinear partial derivative equations system which governs this flow is solved by an explicit unsteady numerical scheme by the finite volume method. The accuracy of the solution depends on the size of the mesh and also the step of time used in the discretized equations. We have chosen in this study the mesh that gives us a stationary solution with good accuracy.

Keywords: finite volume, lunchers, nozzles, shock wave

Procedia PDF Downloads 274

Authors: Omar Imine, Mohammed Aounallah, Mustapha Belkadi

Abstract:

Steady three-dimensional and two free surface waves generated by moving bodies are presented, the flow problem to be simulated is rich in complexity and poses many modeling challenges because of the existence of breaking waves around the ship hull, and because of the interaction of the two-phase flow with the turbulent boundary layer. The results of several simulations are reported. The first study was performed for NACA0012 of hydrofoil with different meshes, this section is analyzed at h/c= 1, 0345 for 2D. In the second simulation, a mathematically defined Wigley hull form is used to investigate the application of a commercial CFD code in prediction of the total resistance and its components from tangential and normal forces on the hull wetted surface. The computed resistance and wave profiles are used to estimate the coefficient of the total resistance for Wigley hull advancing in calm water under steady conditions. The commercial CFD software FLUENT version 12 is used for the computations in the present study. The calculated grid is established using the code computer GAMBIT 2.3.26. The shear stress k-ωSST model is used for turbulence modeling and the volume of the fluid technique is employed to simulate the free-surface motion. The second order upwind scheme is used for discretizing the convection terms in the momentum transport equations, the Modified HRICscheme for VOF discretization. The results obtained compare well with the experimental data.

Keywords: free surface flows, breaking waves, boundary layer, Wigley hull, volume of fluid

Procedia PDF Downloads 354

Authors: Nadim Zgheib, Sivaramakrishnan Balachandar

Abstract:

We present results on the initial formation of ripples from an initially flattened erodible bed. We use direct numerical simulations (DNS) of turbulent open channel flow over a fixed sinusoidal bed coupled with hydrodynamic stability analysis. We use the direct forcing immersed boundary method to account for the presence of the sediment bed. The resolved flow provides the bed shear stress and consequently the sediment transport rate, which is needed in the stability analysis of the Exner equation. The approach is different from traditional linear stability analysis in the sense that the phase lag between the bed topology, and the sediment flux is obtained from the DNS. We ran 11 simulations at a fixed shear Reynolds number of 180, but for different sediment bed wavelengths. The analysis allows us to sweep a large range of physical and modelling parameters to predict their effects on linear growth. The Froude number appears to be the critical controlling parameter in the early linear development of ripples, in contrast with the dominant role of particle Reynolds number during the equilibrium stage.

Keywords: direct numerical simulation, immersed boundary method, sediment-bed interactions, turbulent multiphase flow, linear stability analysis

Procedia PDF Downloads 161

Authors: Tuqa Abdulrazzaq, Hussein Togun, M. K. A. Ariffin, S. N. Kazi, A. Badarudin, N. M. Adam, S. Masuri

Abstract:

Numerical study of heat transfer and fluid flow over vertical double forward facing step were presented. The k-w model with finite volume method was employed to solve continuity, momentum, and energy equations. Different step heights were adopted for range of Reynolds number varied from 10000 to 40000, and range of temperature varied from 310K to 340 K. The straight side of duct is insulated while the side of double forward facing step is heated. The result shows augmentation of heat transfer due to the recirculation region created after and before steps. Effect of step length and Reynolds number observed on increase of local Nusselt number particularly at recirculation regions. Contour of streamline velocity is plotted to show recirculation regions after and before steps. Numerical simulation in this paper done by used ANSYS Fluent 14.

Keywords: turbulent flow, double forward, heat transfer, separation flow

Procedia PDF Downloads 449

Authors: DRIS Mohammed El-Amine, BOUNIF Abdelhamid

Abstract:

This paper considers an experimental and numerical investigation of variable density in axisymmetric turbulent free jets. Special attention is paid to the study of the scalar dissipation rate. In this case, dynamic field equations are coupled to scalar field equations by the density which can vary by the thermal effect (jet heating). The numerical investigation is based on the first and second order turbulence models. For the discretization of the equations system characterizing the flow, the finite volume method described by Patankar (1980) was used. The experimental study was conducted in order to evaluate dynamical characteristics of a heated axisymmetric air flow using the Laser Doppler Anemometer (LDA) which is a very accurate optical measurement method. Experimental and numerical results are compared and discussed. This comparison do not show large difference and the results obtained are in general satisfactory.

Keywords: Scalar dissipation rate, thermal effects, turbulent axisymmetric jets, second order modelling, Velocimetry Laser Doppler.

Procedia PDF Downloads 435

Authors: Marzieh Zarei

Abstract:

Well instability is one of the most fundamental challenges faced by the oil and gas industry. Well wall stability analysis is a gap to be filled in the oil industry. The collection of static data such as well logging leads to the construction of a geomechanical numerical model, which will help in assessing the probable risks in future drilling. In this paper, geomechanical model was designed, and mechanical properties of the rock was determined at all points of the model. It was found the safe mud window was determined and the minimum and maximum mud pressures were determined in the ranges of 70-60 MPa and 110-100 MPa, respectively.

Keywords: geomechanics, numerical model, well stability, in-situ stress, underbalanced drilling

Procedia PDF Downloads 103

Authors: Samuel David Dunstan

Abstract:

An elementary observation of a laminar cylindrical Poiseulle-Couette flow profile reveals no distinction in the parabolic streamwise profile from one without a cross-stream flow in whatever reference frame the observation is made. This is because the laminar flow is in solid-body rotation, and there is no intrinsic fluid rotation. Hence the main streamwise Poiseuille flow is unaffected. However, in turbulent (unsteady) cylindrical Poiseuille-Couette flow, the rotational reference frame must be considered, and any observation from an external inertial reference frame can give outright incorrect results. A common misconception in the study of fluid mechanics is the position of the observer does not matter. In this DNS (direct numerical simulation) study, firstly, turbulent flow in a pipe with axial rotation is established. Then in turbulent flow in the concentric pipe, with inner wall rotation, it is shown how the wall streak direction is oriented by the rotational reference frame. The Coriolis force here is not so fictitious after all!

Keywords: concentric pipe, rotational and inertial frames, frame invariance, wall streaks, flow orientation

Procedia PDF Downloads 68

Authors: Mohammad Moezzibadi, Isabelle Charpentier, Adrien Wanko, Robert Mosé

Abstract:

The calibration of hydrodynamic parameters for subsurface constructed wetlands (CWs) is a sensitive process since highly non-linear equations are involved in unsaturated flow modeling. CW systems are engineered systems designed to favour natural treatment processes involving wetland vegetation, soil, and their microbial flora. Their significant efficiency at reducing the ecological impact of urban runoff has been recently proved in the field. Numerical flow modeling in a vertical variably saturated CW is here carried out by implementing the Richards model by means of a mixed hybrid finite element method (MHFEM), particularly well adapted to the simulation of heterogeneous media, and the van Genuchten-Mualem parametrization. For validation purposes, MHFEM results were compared to those of HYDRUS (a software based on a finite element discretization). As van Genuchten-Mualem soil hydrodynamic parameters depend on water content, their estimation is subject to considerable experimental and numerical studies. In particular, the sensitivity analysis performed with respect to the van Genuchten-Mualem parameters reveals a predominant influence of the shape parameters α, n and the saturated conductivity of the filter on the piezometric heads, during saturation and desaturation. Modeling issues arise when the soil reaches oven-dry conditions. A particular attention should also be brought to boundary condition modeling (surface ponding or evaporation) to be able to tackle different sequences of rainfall-runoff events. For proper parameter identification, large field datasets would be needed. As these are usually not available, notably due to the randomness of the storm events, we thus propose a simple, robust and low-cost numerical method for the inverse modeling of the soil hydrodynamic properties. Among the methods, the variational data assimilation technique introduced by Le Dimet and Talagrand is applied. To that end, a variational data assimilation technique is implemented by applying automatic differentiation (AD) to augment computer codes with derivative computations. Note that very little effort is needed to obtain the differentiated code using the on-line Tapenade AD engine. Field data are collected for a three-layered CW located in Strasbourg (Alsace, France) at the water edge of the urban water stream Ostwaldergraben, during several months. Identification experiments are conducted by comparing measured and computed piezometric head by means of the least square objective function. The temporal variability of hydrodynamic parameter is then assessed and analyzed.

Keywords: automatic differentiation, constructed wetland, inverse method, mixed hybrid FEM, sensitivity analysis

Procedia PDF Downloads 137

Authors: Kruti Jaruriya

Abstract:

Distribution of arsenic (As) and its compound and related toxicology are serious concerns. This is particularly so since millions worldwide are suffering from toxicity due to drinking of As-contaminated groundwater. The Bengal delta plain, formed by the Ganga– Padma–Meghna–Brahmaputra river basin, covering several districts of West Bengal, India and Bangladesh is considered as the worst As affected alluvial basin. However, some equally affected, if not more, areas are emerging in upper Brahmaputra plains. The present study was carried out to examine As contamination trends in the worst affected part of Assam, India. Arsenic (As) mobilization to the groundwater of Brahmaputra floodplains was investigated in Titabor, Jorhat District, located in the North Eastern part of India. The groundwater and the aquifer geochemistry were characterized. The groundwater is characterized by high dissolved Fe, Mn, and HCO-3 and low concentrations of NO-3 and SO2-4 indicating anoxic conditions prevailing in the groundwater. Fifty groundwater samples collected from shallow and deep tubewells of Titabor, Jorhat district (Assam) were examined. Along with total As, examination of concentration levels of other key parameters, viz., pH, EC, Fe, Mn , Mg2+, Ca2+, Na+, K+, PO43- , HCO-3 , NO3- ,Cl - and SO42- was also carried out. In respect to the permissible guideline of World Health Organization (WHO: As 0.01 ppm, Fe 1.0 ppm, and Mn 0.3 ppm for potable water), the range of As concentration in the groundwater varied from 0.014 to 0.604 mg/L with mean concentration 0.184 mg/L. The present study showed that out of the 50 groundwater samples,100%, 54%, and 42% were found contaminated with higher metal contents (for total As, Fe, and Mn, respectively). The results of hydrogeochemical study revealed that the reductive dissolution of MnOOH and FeOOH represents an important mechanism of arsenic release in the study area along with major cations playing an important role in leaching of As into the groundwater. Arsenic released by oxidation of pyrite, as water levels are drawn down and air enters the aquifer, contributes negligibly to the problem of As pollution. Identification of the mechanism of As release to groundwater helps to provide a framework to guide the placement of new water wells so that they will have acceptable concentrations of As.

Keywords: arsenic, assam, brahmaputra floodplain, groundwater, hydrogeochemistry

Procedia PDF Downloads 292