Search for: finite-difference-method
0.007 seconds
Total 132 records

    Assessment of characteristic boundary conditions based on the artificial compressibility method in generalized curvilinear coordinates for solution of the euler equations

    , Article Computational Methods in Applied Mathematics ; 2017 ; 16094840 (ISSN) Parseh, K ; Hejranfar, K ; Sharif University of Technology
    The characteristic boundary conditions are applied and assessed for the solution of incompressible inviscid flows. The two-dimensional incompressible Euler equations based on the artificial compressibility method are considered and then the characteristic boundary conditions are formulated in the generalized curvilinear coordinates and implemented on both the far-field and wall boundaries. A fourth-order compact finite-difference scheme is used to discretize the resulting system of equations. The solution methodology adopted is more suitable for this assessment because the Euler equations and the high-accurate numerical scheme applied are quite sensitive to the treatment of boundary... 

    Preconditioned characteristic boundary conditions based on artificial compressibility method for solution of incompressible flows

    , Article Journal of Computational Physics ; Volume 345 , 2017 , Pages 543-564 ; 00219991 (ISSN) Hejranfar, K ; Parseh, K ; Sharif University of Technology
    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... 

    Wideband maximally flat fractional-delay allpass filters

    , Article Electronics Letters ; Volume 46, Issue 10 , May , 2010 , Pages 722-723 ; 00135194 (ISSN) Jahani Yekta, M. M ; Sharif University of Technology
    The maximally flat (MF) fractional-delay (FD) allpass filter, also known as Thiran's allpass filter, is one of the most popular IIR FD systems which is typically deployed in its causal forms. It is shown that if this causality constraint is removed, MFFD allpass filters with considerably wider bandwidths can be obtained. In many applications this extra bandwidth is worth having a non-causal system  

    Influence of Darcy number on the onset of convection in a porous layer with a uniform heat source

    , Article International Journal of Thermal Sciences ; Volume 47, Issue 8 , August , 2008 , Pages 1020-1025 ; 12900729 (ISSN) Nouri Borujerdi, A ; Noghrehabadi, A. R ; Rees, D. A. S ; Sharif University of Technology
    This note considers the effect of different Darcy numbers on the onset of natural convection in a horizontal, fluid-saturated porous layer with uniform internal heating. It is assumed that the two bounding surfaces are maintained at constant but equal temperatures and that the fluid and porous matrix are in local thermal equilibrium. Linear stability theory is applied to the problem, and numerical solutions obtained using compact fourth order finite differences are presented for all Darcy numbers between Da = 0 (Darcian porous medium) and Da → ∞ (the clear fluid limit). The numerical work is supplemented by an asymptotic analysis for small values Da. © 2007 Elsevier Masson SAS. All rights... 

    Numerical Simulation of Incompressible Flows over two Dimensional Geometries by Means of Immersed Boundary Method

    , M.Sc. Thesis Sharif University of Technology Shahmardi, Armin (Author) ; Hejranfar, Kazem (Supervisor)
    Two-dimensional incompressible flow analysis is one the most important applied issues in engineering and applied science field. Numerical solution of governing equations of flow requires exact computational grid generation.In complex geometries, generation of the grid which is coincident to the body is very difficult and time consuming. Immersed boundary method is an appropriate superseded method of body conformal grid generation in flow field numerical solution. In this method a grid which is not coincidentto bodyis generated and flow field properties are modified on points adjacent to the boundary of the object (Ghost Cell Method) to satisfy boundary conditions.
    The purpose of this... 

    Numerical Solution of Hypersonic Axisymmetric Flows Including Real Gas Effects Using Compact Finite-Difference Scheme

    , M.Sc. Thesis Sharif University of Technology Khodadadi, Polin (Author) ; Hejranfar, Kazem (Supervisor)

    The numerical solution of the parabolized Navier-Stokes (PNS) equations for accurate computation of hypersonic axisymmetric flowfield with real gas effects is obtained by using the fourth-order compact finite-difference method. The PNS 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 the compact PNS scheme to obtain accurate solutions in the vicinity of the shock. To stabilize the numerical solution, numerical dissipation term and filters are used. The main advantage of the present formulation is that the basic flow variables... 

    The effects of proportional loading, plane stress, and constant thickness assumptions on hydro-mechanical deep drawing process

    , Article International Journal of Mechanical Sciences ; Volume 53, Issue 5 , 2011 , Pages 329-337 ; 00207403 (ISSN) Taghipour, E ; Assempour, A ; Sharif University of Technology
    The goal of this study is to evaluate the effects of proportional loading, plane stress, and constant thickness assumptions on hydro-mechanical deep drawing (HDD) by developing analytical models. The main model includes no simplifying assumption, and then each of the mentioned assumptions is considered in a specific model. The interrelationships between geometrical and mechanical variables are obtained in the finite difference form based on the incremental strain theory, thereby being solved by Broyden algorithm. Published experimental and FE results are used for evaluation of the results obtained in the present work. The results of models under proportional loading, plane stress, and... 

    Investigating the propagation noise in PWRs via closed-loop neutron-kinetic/thermal-hydraulic noise calculations

    , Article Annals of Nuclear Energy ; Volume 80 , 2015 , Pages 101-113 ; 03064549 (ISSN) Malmir, H ; Vosoughi, N ; Sharif University of Technology
    Elsevier Ltd  2015
    Neutron noise induced by propagating thermal-hydraulic disturbances (propagation noise for short) in pressurized water reactors is investigated in this paper. A closed-loop neutron-kinetic/thermal-hydraulic noise simulator (named NOISIM) has been developed, with the capability of modeling the propagation noise in both Western-type and VVER-type pressurized water reactors. The neutron-kinetic/thermal-hydraulic noise equations are on the basis of the first-order perturbation theory. The spatial discretization among the neutron-kinetic noise equations is based on the box-scheme finite difference method (BSFDM) for rectangular-z, triangular-z and hexagonal-z geometries. Furthermore, the finite... 

    Simulation of wellbore stability with thermo-hydro-chemo-mechanical coupling in troublesome formations: an example from Ahwaz oil field, SW Iran

    , Article Arabian Journal of Geosciences ; Volume 8, Issue 1 , 2015 , Pages 379-396 ; 18667511 (ISSN) Rafieepour, S ; Jalayeri, H ; Ghotbi, C ; Pishvaie, M. R ; Sharif University of Technology
    Wellbore stability is a main concern in drilling operation. Troublesome drilling issues are chemically active formations and/or high-pressure–high-temperature environments. These are mainly responsible for most of wellbore instabilities. Wellbore failure is mostly controlled by the interaction between active shales and drilling fluid in shale formations. The factors influencing this interaction consist of fluid pressure, temperature, composition of drilling fluid, and exposure time. In this paper, a non-linear fully coupled chemo-thermo-poroelasticity model is developed. At first, a fully implicit finite difference model is presented to analyze the problem, and then, it is verified through... 

    Evaluation of behaviors of earth and rockfill dams during construction and initial impounding using instrumentation data and numerical modeling

    , Article Journal of Rock Mechanics and Geotechnical Engineering ; Volume 9, Issue 4 , 2017 , Pages 709-725 ; 16747755 (ISSN) Rashidi, M ; Haeri, S. M ; Sharif University of Technology
    In this study, the behavior of gavoshan dam was evaluated during construction and the first impounding. A two-dimensional (2D) numerical analysis was conducted based on a finite difference method on the largest cross-section of the dam using the results of instrument measurements and back analysis. These evaluations will be completed in the case that back analysis is carried out in order to control the degree of the accuracy and the level of confidence of the measured behavior since each of the measurements could be controlled by comparing it to the result obtained from the numerical model. Following that, by comparing the results of the numerical analysis with the measured values, it is... 

    A moving-mesh finite-volume method to solve free-surface seepage problem in arbitrary geometries

    , Article International Journal for Numerical and Analytical Methods in Geomechanics ; Volume 31, Issue 14 , 2007 , Pages 1609-1629 ; 03639061 (ISSN) Darbandi, M ; Torabi, S. O ; Saadat, M ; Daghighi, Y ; Jarrahbashi, D ; Sharif University of Technology
    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 higher-order two-dimensional Boussinesq wave model

    , Article Journal of Coastal Research ; Issue SPEC. ISSUE 50 , 2007 , Pages 1183-1187 ; 07490208 (ISSN) Ataie Ashtiani, B ; Najafi Jilani, A ; Sharif University of Technology
    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... 

    On the use of high-order accurate solutions of PNS schemes as basic flows for stability analysis of hypersonic axisymmetric flows

    , Article Journal of Fluids Engineering, Transactions of the ASME ; Volume 129, Issue 10 , 2007 , Pages 1328-1338 ; 00982202 (ISSN) Heiranfar, K ; Esfahanian, V ; Mahmoodi Darian, H ; Sharif University of Technology
    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... 

    2D parallel and stable group explicit finite difference method for solution of diffusion equation

    , Article Applied Mathematics and Computation ; Volume 188, Issue 2 , 2007 , Pages 1184-1192 ; 00963003 (ISSN) Tavakoli, R ; Davami, P ; Sharif University of Technology
    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... 

    Unconditionally stable fully explicit finite difference solution of solidification problems

    , Article Metallurgical and Materials Transactions B: Process Metallurgy and Materials Processing Science ; Volume 38, Issue 1 , 2007 , Pages 121-142 ; 10735615 (ISSN) Tavakoli, R ; Davami, P ; Sharif University of Technology
    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,... 

    New stable group explicit finite difference method for solution of diffusion equation

    , Article Applied Mathematics and Computation ; Volume 181, Issue 2 , 2006 , Pages 1379-1386 ; 00963003 (ISSN) Tavakoli, R ; Davami, P ; Sharif University of Technology
    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  

    Dynamics of multi layer microplates considering nonlinear squeeze film damping

    , Article 2006 ASME International Mechanical Engineering Congress and Exposition, IMECE2006, Chicago, IL, 5 November 2006 through 10 November 2006 ; 2006 ; 1096665X (ISSN); 0791837904 (ISBN); 9780791837900 (ISBN) Ahmadian, M. T ; Moghimi Zand, M ; Borhan, H ; Sharif University of Technology
    American Society of Mechanical Engineers (ASME)  2006
    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... 

    An optimized large-stencil approach for capturing near-PI frequencies

    , Article 12th AIAA/CEAS Aeroacoustics Conference, Cambridge, MA, 8 May 2006 through 10 May 2006 ; Volume 5 , 2006 , Pages 3010-3022 ; 1563478099 (ISBN); 9781563478093 (ISBN) Ghasemi, A ; Sharif University of Technology
    American Institute of Aeronautics and Astronautics Inc  2006
    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... 

    Transient radiative heat transfer analysis of a polypropylene layer using hottel's zonal method

    , Article 8th Biennial ASME Conference on Engineering Systems Design and Analysis, ESDA2006, Torino, 4 July 2006 through 7 July 2006 ; Volume 2006 , 2006 ; 0791837793 (ISBN); 9780791837795 (ISBN) Safavisohi, B ; Sharbati, E ; Aghanajafi, C ; Khatami Firoozabadi, R ; Sharif University of Technology
    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... 

    A study on mutual interaction between atomistic and macroscopic phenomena during electrochemical processes using coupled finite difference - Kinetic Monte Carlo model: Application to potential step test in simple copper sulfate bath

    , Article Journal of Electroanalytical Chemistry ; Volume 588, Issue 2 , 2006 , Pages 267-284 ; 15726657 (ISSN) Saedi, A ; Sharif University of Technology
    Elsevier  2006
    A novel method is used to link a 2D kinetic Monte Carlo code to a 1D finite difference code to construct a more realistic and efficient tool for simulating various electrochemical processes. This multiscale model is able to simulate the long-scale mass transfer of electroactive species in bath along with electrode surface phenomena at atomic scale simultaneously. An embedded atom method (EAM) has been used to evaluate the barrier energies of diffusion and redox reactions on electrode surface. The FD code provides the ion concentration on OHP for KMC code, while the KMC code provides the surface activity and rate constants of redox reactions as an input data for FD code. The electrochemical...