Free Essay

# Finite Different Method

Submitted By vspinftu
Words 769
Pages 4
Finite Difference Methods

1

Finite Difference Methods
A finite difference method obtains a price for a derivative by solving the partial differential equation numerically
Example: • An American put option on a stock that pays a continuous dividend yield q.
ƒ ƒ  2ƒ 2 2  (r  q) S  ½ 2  S  rf t  S S

• Finite difference methods aim to represent the differential equation in the form of a difference equation • We form a grid by considering equally spaced time values and stock price values • Define ƒi,j as the value of ƒ at time iDt when the stock price is jDS
3

Implicit Finite Difference Method
In

ƒ  ƒ 1 2 2  2ƒ  (r  q) S   S  rƒ 2 t S 2 S

we set at the node (i,j):

 ƒ ƒ i ,j 1  ƒ i ,j 1  S 2DS and:  2 ƒ  ƒ i ,j 1  ƒ i ,j ƒ i ,j  ƒ i ,j 1     DS 2 DS DS S    2 ƒ ƒ i ,j 1  ƒ i ,j 1  2ƒ i ,j  2 S D S2

or

If we also set

 ƒ ƒ i 1,j  ƒi ,j  t Dt

we obtain the implicit finite difference method.

Re-arranging we get:

a j ƒi ,j 1  b j ƒi ,j  c j ƒi ,j 1  ƒi 1,j where: The Boundary Conditions
• Example: American put fN,j=0 fN,j=0 N,j

a j ƒi,j 1  b j ƒi,j  c j ƒi ,j 1  ƒi 1,j
• N-1 equations • N-1 unkowns • Exactly 1 solution fN,j=max(K-jΔS,0) f0,j=K

6

Explicit Finite Difference Method
If f S and  2 f S 2 are assumed to be the same at the (i  1,j ) point as they are at the (i,j ) point we obtain the explicit finite difference method f i 1, j 1  f i 1, j 1 f  S 2 DS and : f i 1, j 1  f i 1, j 1  2 f i 1, j  f  2 S DS 2
2

With the explicit method we do not need to solve any equations. We can immediatel y calculate the value for f i , j from the values one time - step ahead : ƒi , j  a* ƒ i 1, j 1  b* ƒ i 1, j  c* ƒ i 1, j 1 j j j
With:

8

Implicit vs Explicit Finite Difference Methods ƒi , j +1 ƒi , j ƒi , j –1
Implicit Method

ƒi +1, j +1 ƒi +1, j ƒi , j ƒi +1, j ƒi +1, j –1
Explicit Method

Implicit

10

Explicit

11

Finite Difference Methods versus Trees
• The explicit finite difference method is equivalent to the trinomial tree approach
Remember:

12

Explicit FD method as a Tree pu pm

pd

We need: pu>0,pm>0 and pu>0

This not always true!

13

• One can show that the implicit finite difference method is equivalent to a multinomial tree approach

14

Improving Efficiency
• It is computationally more to use a FD method with a grid for lnS than S itself. • For Z=lnS we have:

• For the implicit FD method we get: where: 15

• The step-size is now proportional with S • For the implicit method we get: with: 16

• The explicit equation now becomes:

• with:

• If we set DZ   3Dt , then the three probabilities are positive pu pm pd

• The explicit FD method is now identical to the trinomial tree.

17

Crank-Nicolson method
• Rolling back with the implicit method: ƒ i , j  a j ƒ i 1, j 1  b j ƒ i 1, j  c j ƒ i 1, j 1 • Rolling back with the implicit method: ƒ i 1, j  a * ƒ i , j 1  b* ƒ i , j  c * ƒ i , j 1 j j j

• If we sum the two, we get: ƒ i , j  ƒ i , j 1  a j ƒ i 1, j 1  b j ƒ i 1, j  c j ƒ i 1, j 1  a * ƒ i , j 1  b* ƒ i , j  c* ƒ i , j 1 j j j

• Defining gi,j as, g i , j  f i , j  a * ƒ i , j 1  b* ƒ i , j  c * ƒ i , j 1 j j j

• we obtain:

a j ƒ i 1, j 1  1  b j ƒ i 1, j  c j ƒ i 1, j 1  g i , j

• This is the Crank-Nicolson scheme

• Why do we care? • The Crank-Nicolson scheme converges faster then either the implicit or the explicit method.

19

### Similar Documents

#### A Novel Cogging Torque Simulation Method for Permanent-Magnet Synchronous Machines

...com/journal/energies Article A Novel Cogging Torque Simulation Method for Permanent-Magnet Synchronous Machines Chun-Yu Hsiao *, Sheng-Nian Yeh and Jonq-Chin Hwang Department of Electrical Engineering, National Taiwan University of Science and Technology, Taipei 10607, Taiwan; E-Mails: snyeh@mail.ntust.edu.tw (S.-N.Y.); jchwang@ee.ntust.edu.tw (J.-C.H.) * Author to whom correspondence should be addressed; E-Mail: D9607101@mail.ntust.edu.tw; Tel.: +886-921135781; Fax: +886-2-27376699. Received: 16 September 2011; in revised form: 22 November 2011 / Accepted: 30 November 2011 / Published: 6 December 2011 Abstract: Cogging torque exists between rotor mounted permanent magnets and stator teeth due to magnetic attraction and this is an undesired phenomenon which produces output ripple, vibration and noise in machines. The purpose of this paper is to study the existence and effects of cogging torque, and to present a novel, rapid, half magnet pole pair technique for forecasting and evaluating cogging torque. The technique uses the finite element method as well as Matlab research and development oriented software tools to reduce numerous computing jobs and simulation time. An example of a rotor-skewed structure used to reduce cogging torque of permanent magnet synchronous machines is evaluated and compared with a conventional analysis method for the same motor to verify the effectiveness of the proposed approach. The novel method is proved valuable and suitable for large-capacity......

Words: 4883 - Pages: 20

#### Resume

...competitive and supportive work environment. EDUCATIONAL HISTORY   Bachelors in Industrial & Manufacturing Engineering University of Engineering & Technology, Lahore, Pakistan F.sc (Pre-Engineering) Army Public Degree College, Sargodha Cantt. 2009 – 2013 2007 – 2009 CGPA: 3.7 85% MAJORS AND SPECIALIZATION       Management: Production and Operation Management, Operations Research, Total Quality Management, Work Study and Methods Engineering, Engineering Ergonomics Manufacturing: Machining Process, Manufacturing Process Design: Machine Design, Production Tooling Design, Engineering Design and Graphics Mechanics: Mechanics of Materials, Fluid Mechanics, Mechanics of Machines Mechanical: Thermodynamics, Plant Engineering Control Engineering: Computer Integrated Manufacturing, Finite Element Analysis, Automation and Robotics, Instrumentation and Control PROJECTS Senior Year Project  Application and Study of TIG and MIG Welding’s on mechanical and metallurgical properties of low carbon steels and Al alloys under different welding parameters (Current & Voltage) Course Projects  Formulated a business plan for setting up a bullet proof cars factory in Pakistan With complete details in terms of Marketing, Organizational and Financial plans  Analyzed & Identified misalignments and vibrations in gearing system by using vibration analyzer  Investigated & formulated techniques to increase the demand and supply of hamburgers  Studied and Improved the......

Words: 387 - Pages: 2

Free Essay

#### Indian Streams Research Journals

...ISSN:-2230-7850 “Dynamic Load and Stress Analysis of a Crankshaft” Dr. S. V. Deshmukh Prof.& Head, Bapurao Deshmukh College of Engineering, Sevagram, Wardha, India Ram.R.Wayzode Lecturer, Suryodaya College of Engineering & Technology, Nagpur, India Prof. N. G. Alvi Prof,Suryodaya College of Engineering & Technology, Nagpur,India ABSTRACT In this study a dynamic simulation was conducted on a crankshaft from a single cylinder four stroke engine. Finite element analysis was performed to obtain the variation of stress magnitude at critical locations. The pressure-volume diagram was used to calculate the load boundary condition in dynamic simulation model, and other simulation inputs were taken from the engine specification chart. The dynamic analysis was done analytically and was verified by simulation in ADAMS which resulted in the load spectrum applied to crank pin bearing. This load was applied to the FE model in ABAQUS, and boundary conditions were applied according to the engine mounting conditions. The analysis was done for different engine speeds and as a result critical engine speed and critical region on the crankshaft were obtained. Stress variation over the engine cycle and the effect of torsional load in the analysis were investigated. Results from FE analysis were verified by strain gages attached to several locations on the crankshaft. Results achieved from aforementioned analysis can be used in fatigue life calculation and optimization of this......

Words: 3588 - Pages: 15

Free Essay

#### Design and Analysis of a Chassis

...DESIGN AND ANALYSIS OF A CHASSIS FOR A 60 TON PAY LOAD Abstract: As the chassis frame forms the backbone of a heavy vehicle, its principal function is to safely carry the maximum load for all designed operating conditions. To achieve a satisfactory performance, the construction of a chassis is the result of careful design and rigorous testing. Various manufacturers have individual design concepts and different methods of achieving the desired performance standards for the complete chassis, not all chassis components are interchangeable between various makes and models of vehicles. So, there is no standard design for chassis frame. Even though start with the chassis frame design start with selection of the section for side rails and cross members. In this paper we have designed a chassis for storage cum resting fixture. Storage cum resting fixture is a structure of length 11000mm & width 2300mm, used to store cylindrical specimens of various sizes and weights in horizontal configuration. The storage cum resting fixture is used to carry the propellant stored in the cylindrical specimens. The trolley fixture should be designed for a maximum pay load of 60 tons. The Objective of my project is to design a chassis for a pay load of 60 Tons. The design process involves manual design calculations, 3d modeling using UNIGRAPHICS software and analysis to validate the design. Ansys package has been implemented to perform the structural analysis. ...

Words: 13647 - Pages: 55

#### Assignment 1: Finite Element Analysis

...[pic] ENGD3016 Solid Mechanics Assignment 1: Finite Element Analysis Name: Wei Zhang ID: P14021978 Date: Dec17th 2015 Abstract 1.0 Introduction 2.0 Objectives 3.0 Matlab 4.0 Solidworks 4.1 Model of truss 1 4.2 Model of truss 2 5.0 Comparison of the two trusses 6.0 Comparison between MATLAB and SOLIDWORKS 6.1 Comparison of results 6.2 Advantages and disadvantages 7.0 Conclusion Appendix Abstract The purpose of this report is to based on two different 2D pictures of trusses to finite element analysis from two ways. First carried out the selection of materials, and choosing the size of the pipe. Two methods were used two different software for finite element analysis, solidworks and MATLAB. By comparing the simulation results with two software to make sure which is the better way. It was found that the results of MATLAB are more accurate, solidworks results more intuitive. According to the calculation results of the two methods, truss 2 is more reliable. 1.0 Introduction Finite element analysis method using a mathematical approximation of the real physical system to simulation. Also use simple and interacting elements, namely unit, you can use a limited number of unknown quantities to approximate the infinite unknowns real system.The engineering is relying on......

Words: 2550 - Pages: 11

#### A Finite Element Analysis Model for the Behaviour of Cold-Formed Steel Channel Columns Under Axial Thrust

...A FINITE ELEMENT ANALYSIS MODEL FOR THE BEHAVIOUR OF COLD-FORMED STEEL CHANNEL COLUMNS UNDER AXIAL THRUST MASTURA BINTI RAFEK Bachelor of Engineering (Hons) Civil UNIVERSITI TEKNOLOGI MARA 2010 A FINITE ELEMENT ANALYSIS MODEL FOR THE BEHAVIOUR OF COLD-FORMED STEEL CHANNEL COLUMNS UNDER AXIAL THRUST By MASTURA BINTI RAFEK This report is submitted as a requirement for the degree of Bachelor of Engineering (Hons) Civil UNIVERSITI TEKNOLOGI MARA MARCH 2010 DECLARATION BY THE CANDIDATE I (Mastura Binti Rafek, 2007121693) confirm that the work in this report is my own work and the appropriate credit has been given where references have been made to the work of other researchers. (………………………………………….) Student Name : Mastura Binti Rafek Student ID : 2007121693 Date : 31st MARCH 2010 ACKNOWLEDGEMENT First and foremost, I would like to express my gratitude to ALLAH S.W.T for giving me the guidance and strength in making this research project proposal a success. I would like to express my sincerest gratitude to all those who helped me navigate this long and fruitful journey. Sincerest thanks go to my supervisor Mr. Mohd. Raizamzamani Bin Md. Zain for his guidance, knowledge, patience, supervision and encouragement towards me in preparation of this study. His are ever dynamic and also his dedication in encourage of young researcher. Not forget, great appreciation go to my friends for their co- operation in completion this project. Finally, to my......

Words: 19871 - Pages: 80

#### Determination Soil Strength

...is often confusion between the terms Modulus of Elasticity EYoung and Modulus of Deformation EDef for homogenous, isotropic masses, which, sometimes, plainly referred as Modulus E, without any further explanation. In this paper a distinction between the EYoung and the EDef is being done and the interest is concentrated in the Plate Bearing Test, which is commonly used in highway earthworks. Moreover, a soil classification diagram is given, in which the soils are classified by their shear strength parameters and their behavior, as totally elastic or elastoplastic, under specific loading applied by rigid circular plate. Finally, a method for the estimation of the soil parameters is proposed, using only the in-situ Plate Bearing Test. The study and the documentation were supported by a presentation of theoretical examples, based on the Finite Element Method. 1. INTRODUCTION In highway earthworks the knowledge of the Modulus of Elasticity of the soils is very important, whether the projects are in the stage of design, construction or compaction check and this is required to be obtained by an easy manner. Due to the complexity and difficulty in determining the Modulus of Elasticity in the laboratory by the Triaxial Test, which is usually not available on the worksite, the common practice for the designer is to determine the Modulus of Deformation with the in situ and easier Plate Bearing Test. It must be mentioned that, though the soils under load...

Words: 1818 - Pages: 8

Free Essay

#### Fundamentals of Design

Words: 3308 - Pages: 14

#### Bridge Truss

...we not only learned the use of the load resistance factor design but also applied that knowledge — I designed a four-story building. The professor was a practicing engineer, and he always related the subject to real life steel structures he had engineered, for example, the SB Medical Center, an all steel building with a base isolated campus. This is the kind of project on which I would like to work, designing the structure and considering how the building will respond to ground motion. After two quarters of structural analysis, I had come as close as possible to analyzing real world structures. Looking back I realize, I had learned great tools for structural analysis, but my "tool box" was still inadequate. I lacked a very important tool: finite element analysis....

Words: 893 - Pages: 4

Free Essay

#### Stress Concentration Factors Assignment

...DEPARTMENT OF MECHANICAL ENGINEERING KCEC3107 (FINITE ELEMENT ANALYSIS) LAB TITLE: STRESS CONCENTRATION FACTORS ASSIGNMENT LECTURER: MR JUDHA PURBOLAKSONO NAME: DEWI RATIH BINTI MOHAMED RAZZ ROZZMAN MATRIC NO: KEC110007 Assignment Question Produce stress concentration factors for the following geometry by using FEM (ANSYS). The stress concentration factor will have three different ratios of D/d; each D/d will have at least four different values of r/d. Please compare the findings with those in available charts of stress concentration factors. Discuss your finding. The submission will be in the form of a brief report in PDF or MS Word. Assignment Answer Table1.1 Ratio, r/d | Value of Kt | | Experimental D/d, 2.0 | Experimental D/d, 1.5 | Experimental D/d, 1.2 | Ideal D/d, 2.0 | Ideal D/d, 1.5 | Ideal D/d, 1.2 | 0.05 | 2.06 | 2.10 | 2.01 | 2.87 | 2.66 | 2.37 | 0.10 | 2.03 | 1.99 | 2.00 | 2.28 | 2.1 | 1.96 | 0.15 | 1.84 | 1.89 | 1.81 | 2 | 1.87 | 1.76 | 0.20 | 1.73 | 1.83 | 1.76 | 1.84 | 1.73 | 1.63 | Figure 1.1 D/d=2 | Ratio, r/d | | 0.05 | | 0.10 | | 0.15 | | 0.20 | D/d=1.5 | Ratio, r/d | | 0.05 | | 0.10 | | 0.15 | | 0.20 | D/d=1.3 | Ratio, r/d | | 0.05 | | 0.10 | | 0.15 | | 0.20 | In design under fatigue loading, stress concentration factor is used in modifying the values of endurance limit while in design under static loading it simply acts as stress modifier. This means Actual stress=...

Words: 378 - Pages: 2

Free Essay

#### Subsea Pipeline Buckling

...*Corresponding author: jaswar,koto@gmail.com and jaswar@mail.fkm.utm.my Paper History Received: 25-December-2014 Received in revised form: 10-January-2015 Accepted: 18-January-2015 ABSTRACT Oil and gas transportation in subsea operation continues to the extreme depth. Harsh environment in deep water lead to a challenge for especially pipeline design. The pipelines are operated at high pressure and high temperature in order to be able to transport the crude oil from the well to the end termination of loading. Such condition, the pipelines are subjected to axial compressive forces which will cause the pipelines to expand, consequently the pipelines tend to buckle for certain size and distance from the initial of pipeline. The sleeper is one of method to control the pipeline expansion by insertion of bar underneath the pipeline. The sleeper results initial imperfection for pipeline which forms a curvature. The magnitude of curvature is designed comply with DNV OS F101 where the design load will accommodate the combination load works on pipeline and the curvature configuration will validate by using ANSYS 14. ∆ Temperature Difference Axial coefficient of friction Lateral coefficient of friction Submerge weight of the pipeline Cross-sectional area of steel Young’s modulus Second moment of area Constants are dependent on the particular mode of buckling, taken from Hobbs and Liang Anchor length of moving portion of pipeline Hoop stress = ⁄2 1. INTRODUCTION Crude oil production...

Words: 2012 - Pages: 9

Free Essay

#### Introduction to Finite Element Method

...Introduction to Finite Element Method Mathematic Model Finite Element Method Historical Background Analytical Process of FEM Applications of FEM Computer Programs for FEM 1. Mathematical Model (1) Modeling Physical Problems Mathematica l Model Solution Identify control variables Assumptions (empirical law) (2) Types of solution Sol. Eq. Exact Sol. Approx. Sol. Exact Eq. Approx. Eq. ◎ ◎ ◎ ◎ (3) Methods of Solution (3) Method of Solution A. Classical methods They offer a high degree of insight, but the problems are difficult or impossible to solve for anything but simple geometries and loadings. B. Numerical methods (I) Energy: Minimize an expression for the potential energy of the structure over the whole domain. (II) Boundary element: Approximates functions satisfying the governing differential equations not the boundary conditions. (III) Finite difference: Replaces governing differential equations and boundary conditions with algebraic finite difference equations. (IV) Finite element: Approximates the behavior of an irregular, continuous structure under general loadings and constraints with an assembly of discrete elements. 2. Finite Element Method (1) Definition FEM is a numerical method for solving a system of governing equations over the domain of a continuous physical system, which is discretized into simple geometric shapes called finite element. Continuous system Time-independent......

Words: 784 - Pages: 4

Free Essay

#### Representation of Riveted Joints in Fea

...Representation of riveted joints in finite element modelling Modelling of riveted joints in Finite Element Analysis (FEA) programs can be a complicated task. Previously beam elements have been used to represent fasteners where the shear stiffness of the beam has been approximated to that of the shear stiffness of the rivet. Errors result in this method since the beam can also support bending. Test results of rivets joining sheet metal components have shown that beams are not ideal at simulating a rivet in FEA. An acceptable method for simulating riveted joints in FEA is using a spring element. The spring element has the benefit of providing the appropriate shear stiffness whilst not being able to support any bending. Riveted joints that attached thin sheet metal components are capable of transferring shear loads but due to the close proximity of the bending planes, they do not transfer any, or very little, bending loads. The spring element connects two nodes between the components as shown in Figure 1 below. DOUBLER / SPLICE PLATE  KSHEAR  KAXIAL  SKIN  Figure 1 – Spring Element Connection A rivet provides an axial stiffness which is a function of the material properties and the geometry of the rivet. The axial stiffness can be calculated as follows: K AXIAL  AE L Where A = Rivet Cross Sectional Area E = Rivet Modulus of Elasticity L = Rivet Grip Length The shear stiffness of a riveted joint is calculated by the following equation for rivets installed in...

Words: 470 - Pages: 2

#### Engineering Management

...Keeping in view that my undergraduate course would be completed in a few months, I am at the crucial stage wherein I have to decide what I would be doing to improve my career. I have a great compulsion towards practical industrial work environment and therefore I wish to pursue graduation study. Also, the sweeping changes in Global Industrialization and Design Engineering have made the role of mechanical engineer a really indispensable one. It has been my deepest desire to be a part of this rapidly burgeoning community. With due reverence to the Department and the University, I pen down this Statement of purpose. To reach the peaks of perfection, I consider studying in a suitable environment and a good university with a rewarding research program, where the research facilities and the faculty at the university are excellent is of utmost necessity. Right from my schooling I have shown keen interest in mathematics and physics that enhanced quantitative and analytical skills in me. This strength enabled me to secure XX% in XX(10th grade). I passed my senior secondary school with XX% in mathematics and XX% in physics respectively and I gained an overall percentage of XX% marks in Mathematics, Physics and Chemistry, which were my main subjects. The competitive spirit within me acted as a catalyst in all my academic successes. Pursuing engineering as my future was a natural consequence because of my love towards science and technology. To study in this field, I had to......

Words: 808 - Pages: 4

Free Essay

#### Modflow

...Introduction of MODFLOW Yangxiao Zhou PhD Associate Professor in Hydrogeology UNESCO-IHE Institute for Water Education General Information • Developed by McDonald & Harbaugh of the USGS, 1983 • Public Domain available, 1988 • Block-Centered, 3D, modular structural, finite difference groundwater flow model • Most widely used groundwater model • Steady state or transient saturated flow • Currently MODFLOW 2000, 2005 available Mathematical model 3D groundwater flow in the saturated heterogeneous and anisotropic porous media   h    h    h  h  K xx    K yy    K zz   W  Ss x  x  y  y  z  z  t where: Kxx, Kyy, Kzz = values of hydraulic conductivity along xyz axes (LT-1) h = total head (L) W = Sources and sinks (T-1) Ss = Specific storage (L-1) t = time (T) Finite Difference Model Model grid: • Rows, columns, layers • Count from the upper left corner Finite Difference Model Model grid: • Cell-Centered • Horizontal grid: rectangular cell varies with size Finite Difference Model Model grid: • Cell-Centered • Vertical layers: thickness varies MODFLOW Packages • Basic • Flow Package – BCF – LPF – HUF • River • Drain • Well • Recharge • Evapotranspiration • Changing Head Boundary • General Head Boundary • Horizontal Flow Barrier • Stream-Aquifer Interaction • Solvers • Output Control MODFLOW Basic Package (BAS) Functions: • Discretization of model domain; • Selection of major options and the......

Words: 1870 - Pages: 8