Sharif Digital Repository

An unstructured grid, finite-volume, and three-dimensional (3D) primitive equation ocean model has been developed to predict oceanic pollutant dispersions in depth and surface of the Caspian Sea for non-chemical reactions of dissolved constituents. The model consists of momentum, continuity, temperature, salinity, and density equations. Physical and mathematical closure has been achieved using Mellor and Yamada turbulent closure sub-models. Since determining a practical definition of salinity that enjoys acceptable accuracy is difficult; therefore, various definitions have been used in this work. A recent definition of salinity stated in the UNESCO Practical Salinity Scale of 1978, PSS78,... 

The hydrodynamics and energetics of bioinspired oscillating mechanisms have received significant attentions by engineers and biologists to develop the underwater and air vehicles. Undulating and pure heaving (or plunging) motions are two significant mechanisms which are utilized in nature to provide propulsive, maneuvering, and stabilization forces.This study aims to elucidate and compare the propulsive vortical signature and performance of these two important natural mechanisms through a systematic numerical study. Navier-Stokes equations are solved, by a pressure-based finite volume method solver, in an arbitrary Lagrangian- Eulerian (ALE) framework domain containing a 2D NACA0012 foil... 

In the present work, two cavitation modeling strategies, namely the barotropic cavitation model and the transport equation-based model are applied and assessed for the numerical simulation of inviscid cavitating flows over two-dimensional and axisymmetric geometries. The algorithm uses the preconditioned Euler equations employing the interface capturing method for both the cavitation models. A same numerical solution procedure is used herein for discretizing the governing equations resulting from these two cavitation modeling strategies for the assessment to be valid and reliable. A central difference finite-volume scheme employing the suitable dissipation terms to account for density jumps... 

In this paper a locally mass conservative finite volume method is employed to model the one-dimensional, two-phase immiscible flow in a poroelastic media. Since, an appropriate choice of primary variables is critical in simulating multiphase subsurface flow, depending on such a choice, the governing equations can be expressed in different forms. By implementing Picard iteration to a highly nonlinear system of equations, three numerical models including pressure form, mixed form and mixed form with a modified Picard linearization are developed in this study. These models have been evaluated in terms of stability, convergence and mass conservation in various one-dimensional test cases.... 

In this paper, two adaptive refinement techniques are presented for enhancing the capability of the new kind of finite volume method called meshless finite volume (MFV) method in which the moving least squares (MLS) approximation technique is employed. The proposed algorithms perform by inserting new nodes in large error regions identified using the Zienkiewicz-Zhu (Z-Z) error estimator and the T-Belytschko (TB) stress recovery scheme. In the first method referred as CV-based adaptive refinement, the new nodes are inserted at the vertices of control volumes with higher errors. The second method, referred as GA-based adaptive refinement, contains two schemes where an adaptive refinement... 

In order to be a powerful tool, finite-element and finite-volume methods must be capable of handling complex flow in complex geometries. In this work, a structured finite volume element method is suitably developed for solving incompressible flow on a collocated grid topology. The method is generally-applicable to arbitrarily shaped elements and orientations and, thus, challenges the potential to unify many of the different grid topologies into a single formulation. The correct estimation of the convec-tive and diffusive flux terms at cell faces remarkably enhances the solution accuracy of the extended formulation. It is shown that the current formulation is enough robust to treat any... 

Despite great advancement in the lattice Boltzmann method and its application in fluid flow problems, there are still major restrictions in treating either the solution domains with complex boundaries or buoyant flow problems. The past experience shows that the heat equation is a source for instabilities which jeopardizes the stable solution of the lattice Boltzmann method in solving fluid flow problems with heat transfer. The instabilities Increase with increasing buoyant force strength. In this work, we suggest a new approach to overcome the restrictions through implementing the advantages of finite volume method in LBM. In this regard, the lattice Boltzznann equation is incorporated with... 

The main objective of this work is to develop a novel moving-mesh finite-volume method capable of solving the seepage problem in domains with arbitrary geometries. One major difficulty in analysing the seepage problem is the position of phreatic boundary which is unknown at the beginning of solution. In the current algorithm, we first choose an arbitrary solution domain with a hypothetical phreatic boundary and distribute the finite volumes therein. Then, we derive the conservative statement on a curvilinear co-ordinate system for each cell and implement the known boundary conditions all over the solution domain. Defining a consistency factor, the inconsistency between the hypothesis... 

This paper presents an extension of the multiscale finite volume (MsFV) method to multi-resolution coarse grid solvers for single phase incompressible flows. To achieve this, a grid one level coarser than the coarse grids used in the MsFV method is constructed and the local problems are redefined to compute the basis and correction functions associated with this new grid. To construct the coarse-scale pressure equations, the coarse-scale transmissibility coefficients are calculated using a new multi-point flux approximation (MPFA) method. The estimated coarse-scale pressures are utilized to compute the multiscale pressure solution. Finally a reconstruction step is performed to produce a... 

In this study, the hydrodynamics and noise behaviour of marine propellers are examined through finite volume method under various operational conditions. The hydrodynamics of these propellers are studied both numerically and experimentally, and the characteristic curves are produced. Sheet cavitation inception and development conditions are also considered in order to understand the impact of different rotational speeds on propeller noise. Ffowcs Williams and Hawkings equations are used to extract the total sound pressure levels under non-cavitating and sheet cavitating conditions. Finally, the impact of cavitation on increasing the propeller noise is thoroughly explored and the results are... 

In this study, the hydrodynamics and noise behaviour of marine propellers are examined through finite volume method under various operational conditions. The hydrodynamics of these propellers are studied both numerically and experimentally, and the characteristic curves are produced. Sheet cavitation inception and development conditions are also considered in order to understand the impact of different rotational speeds on propeller noise. Ffowcs Williams and Hawkings equations are used to extract the total sound pressure levels under non-cavitating and sheet cavitating conditions. Finally, the impact of cavitation on increasing the propeller noise is thoroughly explored and the results are... 

Adomian decomposition method has been applied to evaluate the conduction-convection heat transfer through a straight fin, property distribution due to convection-diffusion, and conduction heat transfer through a slab with temperature dependent thermal conductivity. The Adomian decomposition method (ADM), provides the closed form solution of the non-linear problems without applying any non-realistic simplifications and/or approximations. The obtained analytical solutions are compared with exact and numerical solutions, using finite volume method. It is shown that the numerical simulation has some limitations and may not always produce correct results. However, the Adomian decomposition method... 

A precise thermal-hydraulics model is of great importance for developing more effective designs of High Temperature Gas Cooled Reactors (HTGR). Recently, several advancements have been made in the methods of analysis of porous media which could be of significant value in the development of more precise and robust codes. The objective of this research is to incorporate some of the most recent improvements in the development of a new 2D program for thermal-hydraulics analysis of modular high temperature reactors. The program is mainly based on the solution of a coupled set of mass, energy and momentum conservation equations for the gas flow, along with the energy conservation equation in the... 

The objective of this article is to study gas temperature estimation in a partially filled rotating cylinder. From the measured temperatures on the shell, an inverse analysis is presented for estimating the gas temperature in an arbitrary cross-section of the aforementioned system. A finite-volume method is employed to solve the direct problem. By minimizing the objective function, a hybrid effective algorithm, which contains a local optimization algorithm, is adopted to estimate the unknown parameter. The measured data are simulated by adding random errors to the exact solution. The effects of measurement errors on the accuracy of the inverse analysis are investigated. Two optimization... 

The objective of this article is to study the estimation of an overall heat transfer coefficient in a partially filled rotating cylinder. Herein is an inverse analysis for estimating the overall heat transfer coefficient in an arbitrary cross-section of the aforementioned system from the temperatures measured on the shell. The material employs the finite-volume method to solve the direct problem. The hybrid effective algorithm applied here contains the local optimization algorithm to estimate the unknown parameter by minimizing the objective function. The data measured here are simulated by adding random errors to the exact solution. An investigation is made of the impact of the measurement... 

A numerical method based on the SIMPLE algorithm has been developed for the analysis of vapor flow in a concentric annular heat pipe. The steady-state response of a concentric annular heat pipe to various heat fluxes in the evaporator and condenser sections are studied. The fluid flow and heat transfer in the annular vapor space are simulated using Navier-Stokes equations. The governing equations are solved numerically, using finite volume approach. The vapor pressure and temperature distributions along a concentric annular heat pipe are predicted for a number of symmetric test cases. The vapor flow reversal and transition to turbulence phenomena are also predicted. The results are compared... 

A numerical model is developed to simulate the borehole heat exchanger both in the short and long time. In this regard, the computational domain is divided into the inside and outside borehole regions. A two-dimensional finite volume method is implemented in a cylindrical coordinate system for modeling of the outside borehole. Also, a thermal resistance-capacity model is presented for the borehole cross section. This model is extended to take into account the fluid transport through the U-tube and the temperature variation of the borehole components with depth. The governing equations of the two regions are solved iteratively in each time step. The proposed model is verified with the... 

In resin injection/compression molding (RI/CM), a preform often comprises layers of different fiber reinforcements. Each fiber reinforcement has unique through thickness and in-plane permeabilities as well as compressibility, creating a heterogeneous porous medium in the mold cavity. In the present article, numerical simulation is utilized to investigate the filling process of RI/CM in such a heterogeneous porous medium. The filling stage is simulated in a full three-dimensional space by using control volume/finite element method and based upon an appropriate filling algorithm. The flow in the open gap which may be present in the mold cavity is modeled by Darcy's law using an equivalent... 

The objective of this study is to develop and apply an arbitrary Lagrangian-Eulerian unstructured finite-volume lattice-Boltzmann method (ALE-FVLBM) for solving two-dimensional compressible inviscid flows around moving bodies. The two-dimensional compressible form of the LB equation is considered and the resulting LB equation is formulated in the ALE framework on an unstructured body-fitted mesh to correctly model the body shape and properly incorporate the mesh movement due to the body motion. The spatial discretization of the resulting system of equations is performed by a second-order cell-centered finite-volume method on arbitrary quadrilateral meshes and an implicit dual-time stepping... 

This paper presents a new landslide-generated wave (LGW) model based on incompressible Euler equations with Savage-Hutter assumptions. A two-layer model is developed including a layer of granular-type flow beneath a layer of an inviscid fluid. Landslide is modeled as a two-phase Coulomb mixture. A well-balanced second-order finite volume formulation is applied to solve the model equations. Wet/dry transitions are treated properly using a modified non-linear method. The numerical model is validated using two sets of experimental data on subaerial and submarine LGWs. Impulsive wave characteristics and landslide deformations are estimated with a computational error less than 5 %. Then, the...