Sharif Digital Repository

The non-linear bubble dynamics equations in a compressible liquid have been modified by considering the effects of compressibility of both the liquid and the gas at the bubble interface. A new bubble boundary equation has been derived, which includes a new term resulting from the liquid bulk viscosity effect. The influence of this term has been numerically investigated by considering the effects of water vapour and chemical reactions inside the bubble. The results clearly indicate that the new term has an important damping role at the collapse, so that its consideration dramatically decreases the amplitude of the bubble rebounds after the collapse. This damping feature is more remarkable for... 

Supercavitation is known as the way of viscous drag reduction for the projectiles, moving in the liquid phase. In recent works, there is distinct investigation between cavitation flow and momentum transfer far away from the cavity surface. However, it seems that there is strong connection between overall flow and what takes place in the sheet cavity where a constant pressure distribution is assumed. Furthermore as we'll see, pressure distribution on cavity surface caused due to overall conditions, induct nonaxisymetric forces and they may need to be investigated. Primarily we describe how pressure distribution into the cavity can cause separation of the aero boundary layer. Then we present... 

This study presents the application of the Proposed Modified Reduced-Order Aerodynamics Modelling approach for aeroelastic analysis based on the boundary element method (BEM) as a novel approach. The used BEM has the capability to capture the thickness effect and geometric complexity of a general three-dimensional model. In this approach the reduced-order aerodynamic model is defined through the eigenvalue problem of unsteady flow based on the unknown wake singularities. Based on the used aerodynamic model an explicit algebraic form of the aeroelastic equations is derived that reduces computational efforts and complexity. This special feature enables us to determine the aeroelastic... 

In this study the developing turbulent swirling pipe flow is investigated both numerically and analytically. Governing equations are derived accompanying the boundary layer assumptions. Uniform and solid body rotation (SBR) distributions are taken into account for the axial and tangential velocities at the inlet of the pipe, respectively. Beyond the boundary layers, the flow pattern is considered to be the potential flow. Making use of the fourth-order Runge-Kutta scheme, the numerical solution of the differential equations is obtained. Further more, by simplifying the governing equations for large Rossby number, the analytical solution is performed. The results of numerical and analytical... 

Convection heat transfer below a horizontal, hot, and isothermal strip of infinite length and width of 2L embedded in fluids with different Prandtl number (Pr) and Nusselt number (Nu) is analyzed with the aid of integral method. A new concept is utilized to determine the boundary layer thickness at the strip's edge, which is based on matching the flow rate of the boundary layer below the strip at its edge and the flow rate of the plume, which forms after the heated fluid detaches from the strip's edge. In addition to these novelties, a numerical model is developed to verify the analytical framework, and an excellent agreement is observed between the analytical and numerical models  

The characteristics of a cavitation water tunnel test setup and the experiments of cavitation around different models are given in this paper. Models of test are formed from combination of head and variable length afterbody. The experiments are performed initially with the smooth head and then with the roughened head. Cavitation initiates both in the wake and after the head. If the nose of models is roughened then cavitation initiates at the lower cavitation numbers. The cavitation inception is a function of the model geometry, nose condition and upstream flow characteristics. During the measuring process of drag force from non-cavitating flow regions to supercavitation case, only the tunnel... 

In this study, factors affecting the noise generation by instability waves in a subsonic jet with acoustic Mach number of 0.5 are investigated using linear stability analysis. The base flow required for instability analysis is obtained by modeling the jet stream based on the k-ϵ turbulence model and using the empirical coefficients suggested by Thies and Tam [1]. The resulting base flow profiles are used to solve the linear instability equation, which governs the pressure perturbation for obtaining the eigenvalues and eigenfunctions. The results of linear instability analysis for phase and amplitude of pressure fluctuations are compared against the existing experimental data, which... 

In this study, factors affecting the noise generation by instability waves in a subsonic jet with acoustic Mach number of 0.5 are investigated using linear stability analysis. The base flow required for instability analysis is obtained by modeling the jet stream based on the k-ϵ turbulence model and using the empirical coefficients suggested by Thies and Tam [1]. The resulting base flow profiles are used to solve the linear instability equation, which governs the pressure perturbation for obtaining the eigenvalues and eigenfunctions. The results of linear instability analysis for phase and amplitude of pressure fluctuations are compared against the existing experimental data, which... 

In the present study, by starting from the reduced form of elasticity displacement field for a long flat laminate, an analytical method is developed in order to accurately calculate the interlaminar stresses near the free edges of generally laminated composite plates under extension. The constant parameter appearing in the reduced displacement field, which describes the global rotational deformation of a laminate, is appropriately obtained by employing an improved first-order shear deformation theory. The accuracy and effectiveness of the proposed first-order theory are verified by means of comparison with the results of Reddy's layerwise theory as a three-dimensional benchmark. Reddy's... 

A CFD code has been developed to describe the salt solution density current, which propagates three-dimensionally in deep ambient water. The height and width of the dense layer are two dominated length scales in a 3-D structure of the density current. In experimental efforts, it is common to measure the height and width of this current via its brightness. Although there are analytical relations to calculate the current height in a two-dimensional flow, these relations cannot be used to identify the width and height of a 3-D density current, due to the existence of two unknown parameters. In the present model, the height and width of the dense layer are obtained by using the boundary layer... 

When velocities in the piping systems change rapidly, spectacular accidents occur, due to tranient-state pressures where the elastic properties of the pipe and liquid must be considered. This hydraulic transient is commonly known as water hammer. A conventional widely-used technique for analyzing this phenomenon is the Method Of Characteristic (MOC), in which, by introducing the characteristic lines, two ordinary differential equations, in lieu of the governing partial differential equations, are produced. In the undisturbed form of the equations, the energy dissipation is evaluated by the steady or quasi-steady approximation. But, there is experimental and theoretical evidence which shows... 

In this paper, heat and mass transportation flow of swimming gyrotactic microorganisms (microbes) and solid nanoparticles under the viscous dissipation effect is investigated. The flow model PDEs are renovated with ordinary ones using suitable boundary layer approximations. The system governing the flow model dimensionless equations as well as boundary conditions is numerically treated with the SOR (successive over relaxation) technique. The flow, heat, and mass transport characteristics are examined against the prime parameters. A comparison is examined to be in a good agreement with the earlier results. It is found here that flow and thermal characteristics of the problem are substantially... 

This study presents a developed successive Boundary Element Method to determine the symmetric and antisymmetric sloshing natural frequencies and mode shapes for multi baffled axisymmetric containers with arbitrary geometries. The developed fluid model is based on the Laplace equation and Green's theorem. The governing equations of fluid dynamic and free surface boundary condition are also applied to proposed model. A zoning method is presented to model arbitrary arrangement of baffles in multi baffled axisymmetric tanks. The influence of each zone on neighboring zones is applied by introducing interface influence matrix which correlates the velocity potential of interfaces to their flux. By... 

A weakly compressible smoothed particle hydrodynamics (WCSPH) method is used along with a new no-slip boundary condition to simulate movement of rigid bodies in incompressible Newtonian fluid flows. It is shown that the new boundary treatment method helps to efficiently calculate the hydrodynamic interaction forces acting on moving bodies. To compensate the effect of truncated compact support near solid boundaries, the method needs specific consistent renormalized schemes for the first and second-order spatial derivatives. In order to resolve the problem of spurious pressure oscillations in the WCSPH method, a modification to the continuity equation is used which improves the stability of... 

This work simulates the turbulent boundary layer of an incompressible viscous swirling flow through a conical chamber. To model the pressure gradient normal to the wall, the radial and tangential velocity components across the boundary layer have been calculated by both the integral and numerical methods. The numerical solution is accomplished by finite difference, based on the finite volume method. The results show that the radial and tangential boundary layer thicknesses depend on the velocity ratios, Reynolds number and nozzle angle. The peak of radial and tangential boundary layer thicknesses are located at zL≈0.2 and zL≈0.8 from the nozzle inlet, respectively. Due to the short length of... 

This is an experimental study on the boundary layer over an airfoil under steady and unsteady conditions. It specifically deals with the effect of plunging oscillation on the laminar/turbulent characteristics of the boundary layer. The wind tunnel measurements involved surface-mounted hot-film sensors and boundary-layer rake. The experiments were conducted at Reynolds numbers of 0.42×10 6 to 0.84 × 10 6 and the reduced frequency was varied from 0.01 to 0.11. The results of the quasi-wall-shear stress as well as the boundary layer velocity profiles provided important information about the state of the boundary layer over the suction surface of the airfoil in both static and dynamic cases. For... 

In this research, the developing turbulent swirling flow in the entrance region of a pipe is investigated analytically by using the boundary layer integral method. The governing equations are integrated through the boundary layer and obtained differential equations are solved with forth-order Adams predictor-corrector method. The general tangential velocity is applied at the inlet region to consider both free and forced vortex velocity profiles. The comparison between present model and available experimental data demonstrates the capability of the model in predicting boundary layer parameters (e.g. boundary layer growth, shear rate and swirl intensity decay rate). Analytical results showed... 

In this study, an arbitrary Lagrangian-Eulerian (ALE) approach is incorporated with a mixed finite-volume-element (FVE) method to establish a novel moving boundary method for simulating unsteady incompressible flow on non-stationary meshes. The method collects the advantages of both finite-volume and finite-element (FE) methods as well as the ALE approach in a unified algorithm. In this regard, the convection terms are treated at the cell faces using a physical-influence upwinding scheme, while the diffusion terms are treated using bilinear FE shape functions. On the other hand, the performance of ALE approach is improved by using the Laplace method to improve the hybrid grids, involving... 

Shock control bump (SCB), suction and blowing are three flow control methods used to control the shock wave/boundary layer interaction to reduce the resulting wave drag in transonic flows. An SCB uses a small local surface deformation to reduce the shock wave strength, while the suction decreases the boundary layer thickness and the blowing delays the flow separation. Here, we will use a multi-point continuous adjoint optimisation scheme to find the optimum design of suction and blowing separately or together, or with the SCB, on two supercritical airfoils, i.e. RAE-5225 and RAE-2822, for a wide range of off-design transonic Mach numbers. The RANS flow equations are solved using the Roe’s... 

Fluid flow and heat transfer at microscale have attracted an important research interest in recent years due to the rapid development of microelectromechanical systems (MEMS). Fluid flow in microdevices has some characteristics which one of them is rarefaction effect related with gas flow. In this research, hydrodynamically and thermally fully developed laminar rarefied gas flow in annular microducts is studied using slip flow boundary conditions. Two different cases of the thermal boundary conditions are considered, namely: uniform temperature at the outer wall and adiabatic inner wall (Case A) and uniform temperature at the inner wall and adiabatic outer wall (Case B). Using the previously...