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 study, the flow of a fiber filled viscoelastic matrix through planar contractions is investigated. It was found that by adding fiber to the matrix vortex, the intensity increases. Fiber orientation along "x" and "y" axes was studied too. It was found that fiber orientation could be used for determining the flow regime through the contraction geometry. The rigidity condition of fibers, which needs the trace of the orientation tensor to be unity everywhere in the domain, is correct except near walls and the reentrant corner, which is slightly less than one. In these regions, the stress magnitude is higher, which results in more numerical errors, and which further leads to some error in... 

In this paper, the truly Meshless Local Petrov-Galerkin (MLPG) method is extended for computation of steady incompressible flows, governed by the Navier-Stokes equations (NSE), in vorticity-stream function formulation. The present method is a truly meshless method based on only a number of randomly located nodes. The formulation is based on two equations including stream function Poisson equation and vorticity advection-dispersion-reaction equation (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. Radial basis functions (RBF) interpolation is employed in shape... 

The accuracy of the large-eddy simulation (LES) of turbulent flows can be increased by using high-order numerical schemes in space and time, due to a decrease in numerical errors. This work investigates a high-order compact finite-volume scheme suitable for LES. The explicit fourth-order Runge-Kutta (RK) scheme for time marching and fourth-order compact schemes for spatial derivatives using a cell-averaged approach are implemented. Different subgrid-scale models and the effect of explicit filtering in a fully turbulent channel flow are studied. In this flow, the fourth-order compact finite-volume method in space, and fourth-order RK in time in conjunction with the dynamic Smagorinsky model... 

The numerical errors associated with explicit upstream finite difference solutions of two-dimensional advection - Dispersion equation with linear sorption are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation. The numerical truncation errors are defined using Peclet and Courant numbers in the X and Y direction, a sink/source dimensionless number and new Peclet and Courant numbers in the XY plane. The effects of these truncation errors on the explicit solution of a two-dimensional advection-dispersion equation with a first-order reaction or degradation are demonstrated by comparison with an analytical solution in... 

In this article, the flow behavior around a NACA 0012 airfoil which is oscillating with different Reynolds numbers and in various amplitudes has been investigated numerically. Numerical simulations have been performed with ANSYS software. First, the 2-D geometry has been studied in different Reynolds numbers and angles of attack with various numerical methods in its static condition. This analysis was to choose the best turbulent model and comparing the grids to have the optimum one for dynamic simulations. Because the analysis was to study the blades of wind turbines, the Reynolds numbers were not arbitrary. They were in the range of 9.71e5 to 22.65e5. The angle of attack was in the range...