Sharif Digital Repository

A finite element solution procedure is presented for the simulation of transient incompressible fluid flows using triangular meshes. The algorithm is based on the artificial compressibility technique in connection with a dual time-stepping approach. A second-order discretization is employed to achieve the required accuracy in real-time while an explicit multistage Runge-Kutta scheme is used to march in the pseudo-time domain. A standard Galerkin finite element method, stabilized by using an artificial dissipation technique, is used for the spatial discretization. The performance of the proposed algorithm is demonstrated by solving a set of internal and external problems including flows with... 

The well-known Rayleigh-Plesset (RP) equation is the base of nearly all hydrodynamical descriptions of the sonoluminescence phenomenon. A major deficiency of this equation is that it accounts for viscosity of an incompressible liquid and compressibility, separately. By removing this approximation, we have modified the RP equation, considering effects of liquid second viscosity. This modification exhibits its importance at the end of an intense collapse, so that the new model predicts the appearance of a new picosecond bouncing during high-amplitude sonoluminescence radiation. This new bouncing produces very sharp (sub-picosecond) peaks on the top of the sonoluminescence pulse. These new... 

The steady-state laminar and incompressible vapor flow in four partially-beated concentric annular beat pipe (CAHP) is studied. The governing equations are solved numerically, using finite volume approach based on collocated grids. The first order upwind scheme and the QUICK scheme are used in the numerical solution. The vapor pressure distributions and velocity profiles along the annular vapor space are predicted for a number of test cases in the range of low to moderate radial Reynolds numbers. The results show that in a partially-heated annular beat pipe, as the radial Reynolds number increases, a number of recirculation zones may be created at both ends of the evaporator and condenser... 

This paper uses a fourth-order compact implicit operator scheme for solving 2D/3D steady incompressible flows using the artificial compressibility method. To stabilize the numerical solution, numerical dissipation terms and/or filters are used. Results obtained for test cases are in good agreement with the available numerical and experimental results. A sensitivity study is also conducted to evaluate the effects of grid resolution and pseudocompressibility parameter on accuracy and convergence rate of the solution. The effects of filtering and numerical dissipation on the solution are also investigated  

The preconditioned characteristic boundary conditions based on the artificial compressibility (AC) method are implemented at artificial boundaries for the solution of two- and three-dimensional incompressible viscous flows in the generalized curvilinear coordinates. The compatibility equations and the corresponding characteristic variables (or the Riemann invariants) are mathematically derived and then applied as suitable boundary conditions in a high-order accurate incompressible flow solver. The spatial discretization of the resulting system of equations is carried out by the fourth-order compact finite-difference (FD) scheme. In the preconditioning applied here, the value of AC parameter... 

The use of classical Bossiness approximation is a straightforward strategy to take into account the buoyancy effect in incompressible solvers. This strategy is highly effective if the density variation is low. However, ignoring the importance of density variation in high thermo buoyant flows can cause considerable deviation in predicting the correct fluid flow behavior and heat transfer phenomenon. Indeed, there are many technological and environmental problems where the Bossiness approximation is not valid. In this study, an incompressible algorithm is suitably extended in order to solve compressible flow problems with natural-convection heat transfer. In this regard, the density field is... 

In this paper, a finite point method is employed to solve the incompressible laminar Stokes flow. A moving least-squares approximation, using linear and quadratic basis functions, in conjunction with a point collocation method, has been utilized to discretize the governing equations. Two examples, including the driven cavity and the fully developed channel flow, are solved showing the accuracy and applicability of the method. In summary, the solutions for the linear basis case exhibit a large sensitivity to the size of the domain of influence of the weighting function, in contrast to the quadratic basis case  

The unsteady preconditioned characteristic boundary conditions (UPCBCs) based on the artificial compressibility (AC) method are formulated and applied at artificial boundaries for the direct numerical simulation (DNS) of incompressible flows. The compatibility equations including the unsteady terms are mathematically derived in the generalized curvilinear coordinates and then incorporated as boundary conditions (BCs) in a high-order accurate incompressible flowsolver. The spatial derivative terms of the systemof equations are discretized using the fourth-order compact finite difference (FD) scheme, consistent with the high-order accuracy required for the DNS. The time integration is carried... 

The use of the classical Boussinesq approximation is a straightforward strategy for taking into account the buoyancy effect in incompressible solvers. This strategy is highly effective if density variation is low. However, ignoring the importance of density variation in highly thermobuoyant flow fields can cause considerable deviation from the correct prediction of fluid flow behavior and the accurate estimation of heat transfer rate. In this study, an incompressible algorithm is suitably extended to solve high-density-variation fields caused by strong natural-convection influence. The key point in this research is the way that an ordinary incompressible algorithm is extended to... 

The objective of this research is to study the ability of a meshless method, called finite point method, in solving incompressible fluid flow problems using two stabilization schemes. The main goal of meshless methods is to reduce or remove the cost of grid generation. This issue is implemented using the satisfaction of governing differential equations on a regular or irregular set of nodes by interpolation functions, based on special least-squares approximations. In this research, the finite point method is used to solve the Stokes and the Navier-Stokes equations by employing two different stabilization schemes. In addition, the effects of least-squares approximations are studied. © Sharif... 

In this work, basic design features of Rayan are documented. One of the new design features presented in this work is the way Rayan handles polyhedral grids. Grid definition is combined with the definition of the structure of the sparse coefficient matrix, thereby releasing a considerable part of the memory used by the grid to store otherwise required faces belonging to the cell part of the connectivity description. The key idea is to use a uniform way for creating the structure of the coefficient matrix from the grid connectivity description and to access that data when computing the elements of the coefficient matrix. This saving requires many modifications to the computational algorithm... 

In this paper, the aeroelastic behavior and flutter instability of aircraft wings in subsonic incompressible flight speed regime are investigated. Quasi-steady and unsteady aerodynamic models are used for aerodynamic modeling and the obtained aeroelastic predictions are compared to those available in the specialized literature. Based on a number of test cases, it is shown that the quasi-steady aerodynamic models are inadequate for the determination of aeroelastic behavior and flutter boundary of aircraft wings in the incompressible flight speed range. © 2006 Elsevier Ltd. All rights reserved  

Hyperelastic materials have high deformability and nonlinearity in load-deformation behavior. Based on a phenomenological approach, these materials are treated as a continuum, and a strain energy density is considered to describe their hyperelastic behavior. In this paper, the mechanical behavior characterization of these materials is studied from the continuum viewpoint. For this purpose, the strain energy density is expressed as sum of independent functions of the mutual multiple of principal stretches. These functions are determined by applying the governing postulates on the form of the strain energy density. It is observed that a consistent strain energy density is expressible in terms... 

In this work, we present a more realistic inlet boundary condition to simulate compressible and incompressible flows through micro and nano channels considering consistent momentum and heat transfer specifications there. At solid walls, a constant wall temperature with suitable jump is applied as the wall thermal boundary condition; however, two types of thermal inlet boundary conditions are investigated at the inlet. We firstly examine the classical inlet boundary condition, which specifies a uniform temperature distribution right at the real inlet. Alternatively, we apply the same boundary condition but at a fictitious place far upstream of the real channel inlet. To validate our results,... 

In this work, we extend a numerical tool capable of solving compressible and incompressible gas flows to study the momentum and heat transfer rates in micro/nano channels with high aspect ratio (L/H = 8000), where the compressibility effect is dominant. The constant heat flux thermal boundary condition is firstly applied at the wall. Next, the flow regime is extended to the early transition regime employing a high order slip velocity boundary condition based on the kinetic theory assumptions. The accuracy of the present results in the slip flow regimes is evaluated against other available theoretical and experimental results. The thermal and compressibility effects on the pressure and... 

In this work, we extend a finite-volume pressure-based incompressible algorithm to solve three-dimensional compressible and incompressible turbulent flow regimes. To achieve a hybrid algorithm capable of solving either compressible or incompressible flows, the mass flux components instead of the primitive velocity components are chosen as the primary dependent variables in a SIMPLE-based algorithm. This choice warrants to reduce the nonlinearities arose in treating the system of conservative equations. The use of a new Favre-averaging like technique plays a key role to render this benefit. The developed formulations indicate that there is less demand to interpolate the fluxes at the cell... 

An incompressible-smoothed particle hydrodynamics (I-SPH) formulation is presented to simulate impulsive waves generated by landslides. The governing equations, Navier-Stokes equations, are solved in a Lagrangian form using a two-step fractional method. Landslides in this paper are simulated by a submerged mass sliding along an inclined plane. During sliding, both rigid and deformable landslides mass are considered. The present numerical method is examined for a rigid wedge sliding into water along an inclined plane. In addition solitary wave generated by a heavy box falling inside water, known as Scott Russell wave generator, which is an example for simulating falling rock avalanche into...