Sharif Digital Repository

Free surface flows are frequently encountered in hydraulic engineering problems including water jets, weirs and around gates. An iterative solution to the incompressible two-dimensional vertical steady Navier-Stokes equations, comprising momentum and continuity equations, is used to solve for the priori unknown free surface, the velocity and the pressure fields. The entire water body is covered by a unstructured finite element grid which is locally refined. The dynamic boundary condition is imposed for the free surface where the pressure vanishes. This procedure is done continuously until the normal velocities components vanish. To overcome numerical errors and oscillations encountering in... 

In this paper, the truly Meshless Local Petrov-Galerkin (MLPG) method is extended for computation of unsteady incompressible flows, governed by the Navier-Stokes equations (NSE), in vorticity-stream function formulation. The present method is a truly meshless method based only on a number of randomly located nodes. The formulation is based on two equations including stream function Poisson equation and vorticity advection-dispersion-reaction eq. (ADRE). The meshless method is based on a local weighted residual method with the Heaviside step function and quartic spline as the test functions respectively over a local subdomain. Moving Least Square approximation (MLS) is employed in shape... 

Velocity distributions in (r,θ) directions are evaluated by solving simultaneous simplified Navier-Stokes equations (NSE) and continuity equation (CE) in polar coordinate. The analysis is based on steady state laminar flow of thin falling liquid film on a horizontal circular tube, for cases in which traction on the film surface is considered negligible. It is a common geometry for part of engineering problems such as evaporator, condenser, absorber, generator of absorption chillers and other similar units in mechanical and chemical engineering. Knowledge of the velocity profiles is usually needed for: 1- solving governing energy and species equations 2- estimating the average and film... 

A new perspective suitable for understanding the details of nonlinear pumping (formation of traveling shocks) inside a pressurized cavity is constructed in this paper. Full compressible axisymmetric three-dimensional Navier-Stokes equations are used as the starting point to cover all complexities of the problem that exceedingly increase for particular ranges of Mach, Reynolds and Prandtl numbers. Then a very high-order numerical method is introduced to preserve the user-defined order of accuracy for practical simulations. For removal of spurious waves, higher-order compact filters are derived. All equations are marched in time using the classical Runge-Kutta algorithm which is appropriate... 

The Adomian Decomposition Method (ADM) for solving the highly non-linear vorticity-stream function formulation of 2D incompressible Navier-Stokes equations has been implemented. The analysis is accompanied by numerical boundary conditions. Also, numerical simulation, using finite difference method (FDM), is performed for comparison purposes. The obtained results only for few terms of the expansion are presented. Because present software such as Mathematica/Maple can not calculate many terms (for example: up to 10 terms) of solution and then ADM approach of this problem is an open problem case  

Adomian decomposition method was employed to obtain an approximate solution to two-dimensional and axisymmetric jet impinging flows in this work. Assumptions have been made to reduce the related full Navier-Stokes equations to a non-linear ordinary differential equation. A trial and error strategy has been used to obtain the constant coefficient in the approximated solution. Velocity Profiles, Shear stresses and displacement thickness were chosen to be compared with accurate numerical data. Sensitivity of the results to the number of terms have been discussed. The results showed for open boundary problems having at least one boundary condition at infinity, the applicable range of the... 

A hydroelastic hybrid model is developed to simulate the fluid-structure interaction in water entry problems using the partitioned approach. The interactions between a flat plate and the water are modeled by a hydroelastic model using explicit and implicit couplings. Both couplings are unstable due to numerical instability associated with the fluid added mass. To overcome the instability, an extended Wagner's model is combined with the hydroelastic model, and a hybrid model is developed. The extended Wagner's model is the extension of the classical Wagner's model that is used to estimate the fluid inertial, damping, and restoring forces of a flexible plate within the potential flow theory.... 

Oil trapping behavior during the pre-flush stage is critically important to evaluate the effectiveness of matrix acidizing for the oil well stimulation. In this study, the visco-capillary behavior of the two-phase flow in the pore-scale is analyzed to investigate the influence of wetting properties for a natural rock sample. A two-dimensional model, based on Cahn-Hilliard phase-field and Navier-Stokes equations, was established and solved using the finite element method. A stability phase diagram for log capillary number (Ca)-log viscosity ratio (M) was constructed and then compared with the reported experimental works. The maximum and minimum ranges of capillary number and viscosity ratio... 

Manta ray swimming or bio-inspiration propulsion system, as a special type of marine propulsion system, is used for submersible vehicles that require high-speed maneuverability and stability, such as glider and AUV. In a manta ray swimming, the thrust force is generated by a couple of undulation and oscillation of wing, so that the direction of undulation wave and oscillation is upright and perpendicular to the direction of thrust force, respectively. It is possible to combine these two movement modes (flapping motion) on the three-dimensional model without considering the effects of wing twisting and flexibility to simplify and better understand the physical behaviors or special study of... 

This study provides a stability analysis of flexible rotating pipes taking into account the simultaneous effects of internal and external fluid loading. Using the Euler-Bernoulli beam assumptions, governing equations for flexural vibrations of rotating pipes are obtained. The internal flow characteristics and the double gyroscopic effect are taken into account when deriving the structural equations coupled with the internal flow. External fluid loading is determined by a special linearization of the Navier-Stokes equations. Considering the circular wall of the pipe as an impermeable boundary to the flow, fluid-induced forcing functions are obtained and then applied to the structural... 

Using our recently developed model, for the first time in the literature, the effect of fluid/fluid and rock/fluid interactions on the performance of Low Salinity Waterflooding (LSWF, as an Enhanced Oil Recovery process) at pore-doublet scale is investigated. The model is incorporated into OpenFOAM and both the Navier-Stokes equation for oil/water two-phase flow and the advection-diffusion equation for ion transport (at both fluid/fluid and rock/fluid interface) are solved via direct numerical simulation (DNS). The model is validated against imbibition and drainage pore-doublet experiments reported in the literature, and then applied to investigate the sole effect of wettability alteration... 

Low salinity waterflooding (LSWF) is a process in which by lowering the ionic strength and/or manipulation of the composition of the injection water, the long term equilibrium in oil/brine/rock system is disturbed to reach a new state of equilibrium through which the oil production will be enhanced due to fluid/fluid and/or rock/fluid interactions. In spite of recent advances in the simulation of the LSWF at core scale and beyond, there are very few works that have modelled and simulated this process at the pore scale specially using direct numerical simulation (DNS). As a result the effects of wettability alteration and/or Interfacial Tension (IFT) change on the distribution of the phases... 

This study provides a stability analysis of flexible rotating pipes taking into account the simultaneous effects of internal and external fluid loading. Using the Euler-Bernoulli beam assumptions, governing equations for flexural vibrations of rotating pipes are obtained. The internal flow characteristics and the double gyroscopic effect are taken into account when deriving the structural equations coupled with the internal flow. External fluid loading is determined by a special linearization of the Navier-Stokes equations. Considering the circular wall of the pipe as an impermeable boundary to the flow, fluid-induced forcing functions are obtained and then applied to the structural... 

Using our recently developed model, for the first time in the literature, the effect of fluid/fluid and rock/fluid interactions on the performance of Low Salinity Waterflooding (LSWF, as an Enhanced Oil Recovery process) at pore-doublet scale is investigated. The model is incorporated into OpenFOAM and both the Navier-Stokes equation for oil/water two-phase flow and the advection-diffusion equation for ion transport (at both fluid/fluid and rock/fluid interface) are solved via direct numerical simulation (DNS). The model is validated against imbibition and drainage pore-doublet experiments reported in the literature, and then applied to investigate the sole effect of wettability alteration... 

Low salinity waterflooding (LSWF) is a process in which by lowering the ionic strength and/or manipulation of the composition of the injection water, the long term equilibrium in oil/brine/rock system is disturbed to reach a new state of equilibrium through which the oil production will be enhanced due to fluid/fluid and/or rock/fluid interactions. In spite of recent advances in the simulation of the LSWF at core scale and beyond, there are very few works that have modelled and simulated this process at the pore scale specially using direct numerical simulation (DNS). As a result the effects of wettability alteration and/or Interfacial Tension (IFT) change on the distribution of the phases... 

A continuum hydrodynamic model with immersed solid/fluid interface is developed for simulating calcite dissolution by hydrochloric acid (HCl) at the pore scale, and is most accurate for a mass-transfer-controlled dissolution regime under laminar flow conditions. The model uses averaged Navier-Stokes equations to model momentum transfer in porous media and adopts a theoretically developed mass-transfer formulation with assumptions. The model includes no fitting parameter and is validated using experimental results. The findings of previous research and existing models are briefly discussed and their shortcomings and advantages are elucidated. The present model is used in some pore-scale... 

In puffer type Sulfur Hexafluoride (SF6) Circuit Breakers (CBs), the refill valve is the main inlet duct for the SF6 gas to enter the puffer volume. When the CB is closing and moving contacts and the nozzle are moving toward fixed contacts, the SF6 gas enters the puffer volume via the refill valve. Since changing the refill valve diameter can alter the direction of the SF6 gas motion in the closing operation of CBs, understanding its effects on CBs' function and performance may help to improve the CB design. This paper has focused on the effects of refill valve diameter on prestrike occurrence in the closing operation of puffer type SF6 CBs. In order to accurately study these effects, some... 

In this paper we will present some approaches on Reynolds stress production by vortex transport phenomena and nonlinear vorticity generation in momentum equation. First of all we represent a history of recent works to describe how fluid particle motions can be associated with Reynolds stress through either displacement or acceleration terms. In the next section we will describe how vortex stretching causes the Reynolds stress production and what is the dominant effect near and far from the boundary where viscous effects have to be considered. On the other hand, some vortex considered methodologies such as those synthesize boundary layer, as a collection of vortical objects seem to be...