Sharif Digital Repository

The numerical solution to the parabolized Navier-Stokes (PNS) and globally iterated PNS (IPNS) equations for accurate computation of hypersonic axisymmetric flowfields is obtained by using the fourth-order compact finite-difference method. The PNS and IPNS equations in the general curvilinear coordinates are solved by using the implicit finite-difference algorithm of Beam and Warming type with a high-order compact accuracy. A shock-fitting procedure is utilized in both compact PNS and IPNS schemes to obtain accurate solutions in the vicinity of the shock. The main advantage of the present formulation is that the basic flow variables and their first and second derivatives are simultaneously... 

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... 

A two-dimensional Boussinesq-type model is presented accurate to O(μ)6 , μ = h0/l0, in dispersion and all consequential order for non-linearity with arbitrary bottom boundary, where h0 is the water depth and l0 is the characteristic wave length. The mathematical formulation is an extension of (4,4) the Padé approximant to include varying bottom boundary in two horizontal dimensions. A higher order perturbation method is used for mathematical derivation of the presented model. A two horizontal dimension numerical model is developed based on derived equations using the Finite Difference Method in higher-order scheme for time and space. The numerical wave model is verified successfully in... 

Simulation of brittle regime machining of materials (such as ceramics) is often difficult because of the complex material removal mechanisms involved. In this study, the discrete element method is used to simulate the dynamic process for machining of slip-cast fused silica ceramics. Flat-joint contact model is exploited to model contacts between particles in synthetic discrete element method models. This contact model is suitable for modeling of brittle materials with high ratios (higher than 10) of unconfined compressive strength to tensile strength. The discrete element method has the ability to simulate initiation, propagation, and coalescence of cracks leading to chip formation in the... 

High-order accurate solutions of parabolized Navier-Stokes (PNS) schemes are used as basic flow models for stability analysis of hypersonic axisymmetric flows over blunt and sharp cones at Mach 8. Both the PNS and the globally iterated PNS (IPNS) schemes are utilized. The IPNS scheme can provide the basic flow field and stability results comparable with those of the thin-layer Navier-Stokes (TLNS) scheme. As a result, using the fourth-order compact IPNS scheme, a high-order accurate basic flow model suitable for stability analysis and transition prediction can be efficiently provided. The numerical solution of the PNS equations is based on an implicit algorithm with a shock fitting procedure... 

Recently various versions of alternating group explicit or alternating group explicit-implicit methods were proposed for solution of diffusion equation. The main benefits of these methods are: good stability, accuracy and parallelizing. But these methods were developed for 1D case and stability and accuracy were investigated for 1D case too. In the present study we extend the new group explicit method [R. Tavakoli, P. Davami, New stable group explicit finite difference method for solution of diffusion equation, Appl. Math. Comput. 181 (2006) 1379-1386] to 2D with operator splitting method. The implementation of method is discussed in details. Our numerical experiment shows that such 2D... 

A new parallel Gauss-Seidel method is presented for solution of system of linear equations related to finite difference discretization of partial differential equations. This method is based on domain decomposition method and local coupling between interfaces of neighbor sub-domains, same as alternating group explicit method. This method is convergent and number of iterations for achieving convergence criteria is near the original Gauss-Seidel method (sometimes better and sometimes worse but difference is very small). The convergence theory is discussed in details. Numerical results are given to justify the convergence and performance of the proposed iterative method. © 2006 Elsevier Inc.... 

A two-dimensional depth-integrated numerical model is developed using a fourth-order Boussinesq approximation for an arbitrary time-variable bottom boundary and is applied for submarine-landslide-generated waves. The mathematical formulation of model is an extension of (4,4) Padé approximant for moving bottom boundary. The mathematical formulations are derived based on a higher-order perturbation analysis using the expanded form of velocity components. A sixth-order multi-step finite difference method is applied for spatial discretization and a sixth-order Runge-Kutta method is applied for temporal discretization of the higher-order depth-integrated governing equations and boundary... 

An unconditionally stable fully explicit finite difference method for solution of conduction dominated phase-change problems is presented. This method is based on an asymmetric stable finite difference scheme for approximation of diffusion terms and application of the temperature recovery method as a phase-change modeling method. The computational cost of the presented method is the same as the explicit method per time-step, while it is free from time-step limitation due to stability criteria. It robustly handles isothermal and nonisothermal phase-change problems and is very efficient when the global temperature field is desirable (not accurate front position). The robustness, stability,... 

A new group explicit method for solution of diffusion equation is presented. This method is based on domain decomposition concept and using asymmetric Saul'yev schemes for internal nodes of each sub-domain and alternating group explicit method for sub-domain's interfacial nodes. This new method has several advantages such as: good parallelism, unconditional stability, fully explicit nature and better accuracy than original Saul'yev schemes. The details of implementation and proving stability are briefly discussed. Numerical experiments on stability and accuracy are also presented. © 2006 Elsevier Inc. All rights reserved  

In the present study a finite difference method has been developed to model the transient fluid flow and heat transfer. A single fluid has been selected for modeling of mold filling and The SOLA-VOF 3D technique was modified to increase the accuracy of simulation of filling phenomena for shape castings. The model was then evaluated with the experimental methods. Refereeing to the experimental and simulation results a good consistency and the accuracy of the suggested model are confirmed. © 2005 Published by Elsevier B.V  

This paper presents a model to analyze pull-in phenomenon and dynamics of multi layer microplates using coupled finite element and finite difference methods. Firstorder shear deformation theory is used to model dynamical system using finite element method, while Finite difference method is applied to solve the nonlinear Reynolds equation of squeeze film damping. Using this model, Pull-in analysis of single layer and multi layer microplates are studied. The results of pull-in analysis are in good agreement with literature. Validating our model by pull-in results, an algorithm is presented to study dynamics of microplates. These simulations have many applications in designing multi layer... 

The numerical solution of the parabolized Navier-Stokes (PNS) and globally iterated PNS (IPNS) equations for accurate computation of hypersonic axisymmetric flowfields is obtained by using the fourth-order compact finite-difference method. The PNS and IPNS equations in the general curvilinear coordinates are solved by using the implicit finite-difference algorithm of Beam and Warming type with a high-order compact accuracy. A shock fitting procedure is utilized in both the compact PNS and IPNS schemes to obtain accurate solutions in the vicinity of the shock. The main advantage of the present formulation is that the basic flow variables and their first and second derivatives are... 

Large-stencil schemes which their spectral properties are acceptable in the vicinity of ω = π are analyzed for the first time. A machine independent model for evaluating the efficiency of generalized time-marching finite-difference algorithms over periodic domains is developed. This model which is based on operation count reveals that for small values of Total Computational Cost(TCC), the previous low-order small-stencil schemes are more efficient while for moderate TCC, the efficiency of optimized large-stencil schemes abruptly increases. This important result is the motivation for developing optimized large-stencil schemes. The current schemes are successfully implemented in a full... 

This study is concerned with numerical modeling of viscous surface wave motion using boundary element method (BEM). The equations of motion for thin boundary layers at the solid surfaces are coupled with the potential flow in the bulk of the fluid, and a mixed BEM-finite difference technique is used to obtain the viscosity-related quantities such as wave damping rate, shear stress, and velocity distribution inside the boundary layer. The technique is presented for standing surface wave motion. An excellent agreement is obtained between the numerical predictions and the previous results. The extension to other free surface problems is straightforward. © 2005 Elsevier Ltd. All rights reserved  

A numerical model for simulation of liquid/gas phase flow during mold filling is presented. The incompressible Navier-Stokes equations are discretized on a fixed Cartesian mesh with finite difference method. The fractional-step scheme is employed for enforcing incompressibility constraint. The free surface effects are calculated using the volume of fluid method based on the piecewise-linear interface reconstruction and split Lagrangian advection of volume fraction field. Adding limited compressibility to the gas phase led to improvement in convergence rate of Poisson equation solver (about 2-fold). This new concept permits simulation of two-phase incompressible free surface flow during mold... 

Polymer films possess excellent optical properties, such as high transparency, and thermal characteristics, like low heat conductivity, as well as further polymer specific advantages. Consequently, polymer films have an outstanding potential for many solar applications. They are already used for encapsulation of photovoltaic (PV) cells, as convection barrier in solar collectors and as substrate or adhesive layers for glazing. In translucent polymers, energy can be transferred internally by radiation in addition to conduction. Since radiant propagation is very rapid, it can provide energy within the layer more quickly than diffusion by heat conduction. Thus, the transient thermal response of... 

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  

A multiscale FD-KMC model has been developed to simulate the cyclic voltammetry test of a copper electrode in simple copper sulfate bath. In this coupled model, the FD code provides the cupric ion concentration on OHP for KMC code, while the KMC code provides the electrochemical properties of the copper electrode (surface activity and rate constants of redox reactions) as an input data for FD code. The changes in the electrode properties due to the atomistic phenomena (deposition dissolution and surface diffusion) have been studied for the present potentiodynamic system. The results showed that the CV process consists of some distinct stages, so that the electrode exhibits a specific... 

An effective numerical technique is presented to model turbulent motion of a standing surface wave in a tank. The equations of motion for turbulent boundary layers at the solid surfaces are coupled with the potential flow in the bulk of the fluid, and a mixed BEM-finite difference technique is used to model the wave motion and the corresponding boundary layer flow. A mixing-length theory is used for turbulence modelling. The model results are in good agreement with previous physical and numerical experiments. Although the technique is presented for a standing surface wave, it can be easily applied to other free surface problems. Copyright © 2005 John Wiley & Sons, Ltd