Search for: finite-difference-method
0.008 seconds
Total 158 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... 

    Reflection analysis of the end-facet dielectric slab waveguide by FDTD method

    , Article ICCEA 2004 - 2004 3rd International Conference on Computational Electromagnetics and its Applications, Beijing, 1 November 2004 through 4 November 2004 ; 2004 , Pages 453-456 ; 0780385624 (ISBN) Vahidpour, M ; Shishegar, A. A ; Sharif University of Technology
    The Finite Difference Time Domain (FDTD) method has been applied to the analysis of abruptly-ended dielectric waveguides. In these waveguides, incident propagating wave undergoes reflection in an interaction with the end-facet. As a result of the discontinuity, all possible propagating modes may be excited. The constituent propagating modes are extracted from the reflected wave by the least square method. Thus, we present a good estimation of the amplitudes of the reflected modes. This full wave analysis technique is also capable of analyzing any arbitrarily shaped facet. © 2004 IEEE  

    Transient and stability analysis in single-phase natural circulation

    , Article Annals of Nuclear Energy ; Volume 31, Issue 10 , 2004 , Pages 1177-1198 ; 03064549 (ISSN) Mousavian, S. K ; Misale, M ; D'Auria, F ; Salehi, M. A ; Sharif University of Technology
    This paper presents the mathematical modeling of single-phase natural circulation of the University of Genoa's rectangular loop (LOOP#1) by a computer program and using RELAP5 system code. The mass, momentum and energy conservation equations in transient form were solved numerically using the finite difference method. One-dimensional linear stability analysis was performed for the single-phase natural circulation loop and the numerical perturbation technique was used in this analysis. The Nyquist criterion was employed to find the stability map of the LOOP#1. The obtained transient results using the first order upwind scheme of the fluid temperatures in various sectors of the LOOP#1 are... 

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

    Mode identification of high-amplitude pressure waves in liquid rocket engines

    , Article Journal of Sound and Vibration ; Volume 229, Issue 4 , 2000 , Pages 973-991 ; 0022460X (ISSN) Ebrahimi, R ; Mazaheri, K ; Ghafourian, A ; Sharif University of Technology
    Identification of existing instability modes from experimental pressure measurements of rocket engines is difficult, specially when steep waves are present. Actual pressure waves are often non-linear and include steep shocks followed by gradual expansions. It is generally believed that interaction of these non-linear waves is difficult to analyze. A method of mode identification is introduced. After presumption of constituent modes, they are superposed by using a standard finite difference scheme for solution of the classical wave equation. Waves are numerically produced at each end of the combustion tube with different wavelengths, amplitudes, and phases with respect to each other. Pressure... 

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

    Numerical Simulation Cavitating Flows Using Compact Finite-difference Scheme

    , M.Sc. Thesis Sharif University of Technology Shokri, Maryam (Author) ; Hejranfar, Kazem (Supervisor)
    In the study, the simulation of two-dimensional cavitating flows is performed by applying a high-order accurate numerical method to the preconditioned, homogenous, multiphase Navier-Stokes equations. The baseline differential equations system is comprised of the mixture volume, mixture momentum and constituent volume fraction equations. A coordinate transformation is applied and the resulting system of governing equations in curvilinear coordinates is discretized using a fourth-order compact finite-difference scheme. The high-order accurate numerical scheme employing the suitable linear and nonlinear filters to account for density jumps across the cavity interface is shown to yield an... 

    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