Lecture of FreeFEM Update:2015/7/29
|
|
- Brittany Sharyl Powell
- 5 years ago
- Views:
Transcription
1 Lecture of FreeFEM Update:2015/7/29
2 About this lecture The goal of this lecture is: At the first, 3~4 participants are grouped, They understand how to use throughout total 10 problems, They should teach each other in the same group They can search anything from internet.
3 Schedule 9:00-10:20 Running Sample files and Mesh Generating 10:30-11:50 Poisson equations 12:00-13:30 13:30-14:50 Lunch Time Convection-Diffusion equations 15:00-16:20 Navier-Stokes equations 16:30-17:30 Free Time One lecture is 80 minutes in which, 30 minutes is used for explanations from me, 50 minutes is used for exercise.
4 Grouping and notations Grouping 3~4 participants are grouped. Notations About evaluations from A. Prof. J. Masamune
5 Schedule 9:00-10:20 Running Sample files and Mesh Generating 10:30-11:50 Poisson equations 12:00-13:30 Lunch Time 13:30-14:50 Convection-Diffusion equations 15:00-16:20 Navier-Stokes equations 16:30-17:30 Free Time
6 Running Sample files Laplace.edp diffusion.edp convection.edp tunnel.edp LapComplexEigenValue.edp Mesh_square.edp Mesh_circle.edp Mesh_circle_in_square.edp I show you how to generate finite element meshes
7 Mesh generation(mesh_square.edp) border a0(t=1,0){ x=0; y=t; label=1;} border a1(t=0,1){ x=t; y=0; label=2;} border a2(t=0,1){ x=1; y=t; label=3;} border a3(t=1,0){ x=t; y=1; label=4;} int n=5; Mesh Th =buildmesh( a0(10*n)+a1(10*n)+a2(10*n)+a3(10*n));
8 Mesh generation(mesh_square.edp) border a0(t=1,0){ x=0; y=t; label=1;} border a1(t=0,1){ x=t; y=0; label=2;} border a2(t=0,1){ x=1; y=t; label=3;} border a3(t=1,0){ x=t; y=1; label=4;} int n=5; Mesh Th =buildmesh( a0(10*n)+a1(10*n)+a2(10*n)+a3(10*n));
9 Command for definitions of borders border a0(t=1,0){ x=0; y=t; label=1;}
10 border a0(t=1,0){ x=0; y=t; label=1;} The name of border. You can use any words which you want.
11 border a0(t=1,0){ x=0; y=t; label=1;} A parameter range from 1 to 0. You can use any number which you want.
12 positions of borders (Like position vector?) border a0(t=1,0){ x=0; y=t; label=1;} y (0,1) o (0,0) x
13 positions of borders (Like position vector?) border a0(t=0,1){ x=0; y=t; label=1;} y (0,1) o (0,0) x
14 y (0,1) border a3(t=1,0){ x=t; y=1; label=4;} (1,1) border a0(t=1,0){ x=0; y=t; label=1;} border a2(t=0,1){ x=1; y=t; label=3;} border a1(t=0,1){ x=t; y=0; label=2;} (0,0) (1,0) x
15 y Notatins: definitions of borders in the Anti-Crockwise direction should be done! (0,1) border a3(t=1,0){ x=t; y=1; label=4;} (1,1) border a0(t=1,0){ x=0; y=t; label=1;} border a2(t=0,1){ x=1; y=t; label=3;} border a1(t=0,1){ x=t; y=0; label=2;} (0,0) (1,0) x
16 Mesh generation(mesh_square.edp) border a0(t=1,0){ x=0; y=t; label=1;} border a1(t=0,1){ x=t; y=0; label=2;} border a2(t=0,1){ x=1; y=t; label=3;} border a3(t=1,0){ x=t; y=1; label=4;} int n=5; Mesh Th =buildmesh( a0(10*n)+a1(10*n)+a2(10*n)+a3(10*n));
17 Commands for definitions of finite element meshes Mesh Th =buildmesh( a0(10*n)+a1(10*n)+a2(10*n)+a3(10*n)); Commands for definitions of finite element meshes
18 The name of Mesh You can use any words which you want. Mesh Th =buildmesh( a0(10*n)+a1(10*n)+a2(10*n)+a3(10*n));
19 Fineness of the mesh: should be defined by integer Mesh Th =buildmesh( a0(10*n)+a1(10*n)+a2(10*n)+a3(10*n));
20 Mesh generation(mesh_square.edp) border a0(t=1,0){ x=0; y=t; label=1;} border a1(t=0,1){ x=t; y=0; label=2;} border a2(t=0,1){ x=1; y=t; label=3;} border a3(t=1,0){ x=t; y=1; label=4;} int n=5; Mesh Th Making boundaries =buildmesh( a0(10*n)+a1(10*n)+a2(10*n)+a3(10*n)); Generating finite elements(triangles)
21 To define inside boundary Mesh_circle_in_square.edp border a0(t=1,0){ x=0; y=t; label=1;} border a1(t=0,1){ x=t; y=0; label=2;} border a2(t=0,1){ x=1; y=t; label=3;} border a3(t=1,0){ x=t; y=1; label=4;} border a4(t=0,-2*pi){ x= *cos(t); y= *sin(t); label=1;} int n=5; mesh Th= buildmesh( a0(10*n)+a1(10*n)+a2(10*n)+a3(10*n)+a4(10*n));
22 y Notatins: definitions of outside borders in the Anti-Crockwise direction should be done! (0,1) border a3(t=1,0){ x=t; y=1; label=4;} (1,1) border a0(t=1,0){ x=0; y=t; label=1;} border a2(t=0,1){ x=1; y=t; label=3;} border a1(t=0,1){ x=t; y=0; label=2;} (0,0) (1,0) x
23 y Notatins: definitions of inside borders in the Crockwise direction should be done! (0,1) (1,1) border a4(t=0,-2*pi){ x= *cos(t); y= *sin(t); label=1;} (0,0) (1,0) x
24 Problem1 y o x 5 10
25 Problem2 y o 1 2 x -1
26 Schedule 9:00-10:20 Running Sample files and Mesh Generating 10:30-11:50 Poisson equations 12:00-13:30 13:30-14:50 Lunch Time Convection-Diffusion equations 15:00-16:20 Navier-Stokes equations 16:30-17:30 Free Time
27 Strong form of Poisson eq. Γ R 2 Δu f = 0 in u = 0 on Γ
28 Weak form of Poisson eq. = lhs = Δu f v d = 0 u vd + x x + y y fvd v d + fvd = x x v + y y v d + x = v + x y v d + v x x + v y y = v + Γ x y v ndγ + v x x + v y y = Γ n vdγ + v x x + v y y v x + v y y d + fvd d + fvd d + fvd d + fvd
29 Strong form and Weak form of Poisson eq. Γ R 2 Δu f = 0 in u = 0 on Γ x v x + y v y d Γ n vdγ + fvd = 0
30 Boundary Conditions x v x + y v y d Γ n vdγ + fvd = 0 In the case of Dirichlet condition Let the trial function v be 0 on the boundary, we have v x x + v y y d + fvd = 0 But commands for the Dirichlet condition should be written in the source code in FreeFEM++. In the case of Neumann condition Substitute n v x x + y = g into boundary integration term, we have v y d g vdγ + Γ fvd = 0
31 1 Γ 1 Poisson equation: Laplace.edp o y Γ 4 R 2 Γ 3 x Δu = 0 in u = 1 on Γ 1 u = 0 on Γ 2 n = 0 on Γ 3 n = 0 on Γ 4 Γ 2 1
32 Poisson equation: Laplace.edp fespace Vh(Th,P2); Vh uh,vh; problem laplace(uh,vh,solver=lu) = int2d(th)( dx(uh)*dx(vh) + dy(uh)*dy(vh) ) + on(1,uh=1) + on(2,uh=0) ; laplace; plot(uh); Definition of variable Running Numerical calculations Visualizations Definition of finite element space Weak forms and numerical scheme
33 Poisson equation: Laplace.edp Mesh Th=buildmesh( a0(10)+a1(10)+a2(10)+a3(10)); Commands for definitions of finite element space fespace Vh(Th,P2); Vh uh,vh;
34 Poisson equation: Laplace.edp Mesh Th=buildmesh( a0(10)+a1(10)+a2(10)+a3(10)); The name of finite element space You can use any words which you want. fespace Vh(Th,P2); Vh uh,vh;
35 Poisson equation: Laplace.edp Mesh Th=buildmesh( a0(10)+a1(10)+a2(10)+a3(10)); fespace Vh(Th,P2); Vh uh,vh;
36 Poisson equation: Laplace.edp Mesh Th=buildmesh( a0(10)+a1(10)+a2(10)+a3(10)); fespace Vh(Th,P1); Variables are calculated on 3 points on one triangle Vh uh,vh; P1 element
37 Poisson equation: Laplace.edp Mesh Th=buildmesh( a0(10)+a1(10)+a2(10)+a3(10)); fespace Vh(Th,P2); Variables are calculated on 6 points on one triangle Vh uh,vh; P2 element
38 Poisson equation: Laplace.edp Mesh Th=buildmesh( a0(10)+a1(10)+a2(10)+a3(10)); fespace Vh(Th,P1 or P2); P1: low cost and low accuracy P2: high cost and high accuracy P1 element Vh uh,vh; P2 element
39 Poisson equation: Laplace.edp Commands for definition of numerical calculations problem laplace(uh,vh,solver=lu) = int2d(th)( dx(uh)*dx(vh) + dy(uh)*dy(vh) ) + on(1,uh=1) + on(2,uh=0) ; laplace;
40 Poisson equation: Laplace.edp The name of the problem You can use any words which you want. problem laplace(uh,vh,solver=lu) = int2d(th)( dx(uh)*dx(vh) + dy(uh)*dy(vh) ) + on(1,uh=1) + on(2,uh=0) ; laplace;
41 Poisson equation: Laplace.edp Variables that you need to calculate problem laplace(uh,vh,solver=lu) = int2d(th)( dx(uh)*dx(vh) + dy(uh)*dy(vh) ) + on(1,uh=1) + on(2,uh=0) ; laplace;
42 Poisson equation: Laplace.edp Numerical Scheme LU decompositions problem laplace(uh,vh,solver=lu) = int2d(th)( dx(uh)*dx(vh) + dy(uh)*dy(vh) ) + on(1,uh=1) + on(2,uh=0) ; laplace;
43 Poisson equation: Laplace.edp Integral on finite elements Th If you want to Integrate in 2 dimension, you should int2d. int2d(th)( dx(uh)*dx(vh) + dy(uh)*dy(vh) ) Weak form of Poisson equation dx() and dy() mean derivatives with respect to x and y. x v x + y v y d
44 Poisson equation: Laplace.edp problem laplace(uh,vh,solver=lu) = int2d(th)( dx(uh)*dx(vh) + dy(uh)*dy(vh) ) + on(1,uh=1) + on(2,uh=0) ; 1 Γ 4 (Dirichlet conditions) Dirichlet condition u = 1 on Γ 1 u = 0 on Γ 2 border a0 Γ 1 Γ 3 border a0(t=1,0){ x=0; y=t; label=1;} border a1(t=0,1){ x=t; y=0; label=2;} Γ 2 border a1 1
45 Poisson equation: Laplace.edp problem laplace(uh,vh,solver=lu) = int2d(th)( dx(uh)*dx(vh) + dy(uh)*dy(vh) ) + on(1,uh=1) + on(2,uh=0) ; Neumann conditions x v x + y v y d Γ n vdγ + fvd = 0 n = 0 on Γ 3 n = 0 on Γ 4 Substitute n = 0 into the second term, You should not type anything.
46 If you want to use Neumann conditions x v x + y v y d g vdγ + Γ fvd = 0 Please type boundary integration as follows; ex) +int1d(th, (label number))(- g *vh) on the boundary border a1, for g =2, you should type +int1d(th, 2)(-2*vh) border a1(t=0,1){ x=t; y=0; label=2;}
47 Poisson equation: Laplace.edp Integral on finite elements Th If you want to Integrate in 2 dimension, you should int2d. int2d(th)( dx(uh)*dx(vh) + dy(uh)*dy(vh) ) Weak form of Poisson equation dx() and dy() mean derivatives with respect to x and y. x v x + y v y d
48 If you want to use Neumann conditions x v x + y v y d g vdγ + Γ fvd = 0 Please type boundary integration as follows; +int1d(th, label number)(- g *vh) ex) on the boundary border a1, for g =2, you should type +int1d(th, 2)(-2*vh) border a1(t=0,1){ x=t; y=0; label=2;}
49 Poisson equation: Laplace.edp problem laplace(uh,vh,solver=lu) = int2d(th)( dx(uh)*dx(vh) + dy(uh)*dy(vh) ) + on(1,uh=1) + on(2,uh=0) ; Dirichlet condition laplace; Numerical Calculation
50 Output with Paraview Laplace_output.edp Please add the following commands in the line where you want output data. load "iovtk"; string vtkoutputfile="./output.vtk";savevtk(vtkoutputf ile, Th, uh, dataname="poisson");
51 Problem3 1 Γ 1 y Γ 4 R 2 Γ 3 Δu = u in u = x x x on Γ 1 u = 0 on Γ 2 n = 0 on Γ 3 n = 1 on Γ 4 o x Γ 2 1
52 Answer of Problem3
53 y Problem4 (The same domain as Problem1) Γ top Γ in Γ circle Γ out Γ bottom x
54 Problem4 (The same domain as Problem1) Δu = sin x cos(y) in n = 0 on Γ top Γ bottom u = 1 on Γ in Γ out n = 1 on Γ circle
55 Answer of Problem4
56 Γ 1 z Problem 5 In the axi-symmetric Cylindrical coordinate 1 system: Γ 4 R 2 o 1 Γ 2 Γ 3 r Δu = 1 r Space integration In this system: 2 u r r r + 2 u z 2 = 0 in u = 1 on Γ 2 u = 0 on Γ 3 n = 1 on Γ 4 r r r = 2 u r r r + 2 u z 2 Δu vrd = 0
57 Answer of Problem 5
58 Schedule 9:00-10:20 Running Sample files and Mesh Generating 10:30-11:50 Poisson equations 12:00-13:30 13:30-14:50 Lunch Time Convection-Diffusion equations 15:00-16:20 Navier-Stokes equations 16:30-17:30 Free Time
59 Diffusion equation: diffusion.edp Γ 1 Re t 1 Re Δu = 0 un+1 u n 1 Δt Re Δun+1 = 0 u n+1 u n 1 Δt Re Δu wd = 0 u n+1 Δt n+1 un w Δt w + 1 Re un+1 w n w = 0 d
60 Convection equation:convection.edp + u u = 0 t un+1 u n + u = 0 Δt u n+1 u n + u u n wd = 0 Δt u n+1 Δt u n+1 w un Δt w + u un w d = 0 Δt w w Δt convect u, Δt, un d = 0 u n+1 = convect u, Δt, u n = 0
61 u n+1 Δt u n+1 Convection diffusion equation: convection_diffusion.edp t + u u 1 Re Δu = 0 un+1 u n + u u n 1 Δt Re Δun+1 = 0 u n+1 u n + u u n 1 Δt Re Δun+1 wd = 0 w un Δt w + un u n w + 1 Re un+1 w d Γ Δt w w Δt convect u, Δt, un + 1 Re un+1 w d Γ Re 1 Re 1 n+1 n w = 0 n+1 n w = 0
62 Problem6 (The same domain as Problem1) Δu = 0 in t Initial condition: u = exp 10 x y Boundary conditions: u = 0 on Γ top [0, t] u = 0 on Γ bottom [0, t] n = 0 onγ in [0, t] n = 0 on Γ out [0, t] u = 0 on Γ circle [0, t]
63 Problem7 (The same domain as Problem1) + f u Δu = 0 in t where f = ( y y 5, 0) Initial condition: u = exp 10 x y Boundary conditions: u = 0 on Γ top [0, t] u = 0 on Γ bottom [0, t] n = 0 onγ in [0, t] n = 0 on Γ out [0, t] u = 0 on Γ circle [0, t]
64 Schedule 9:00-10:20 Running Sample files and Mesh Generating 10:30-11:50 Poisson equations 12:00-13:30 13:30-14:50 Lunch Time Convection-Diffusion equations 15:00-16:20 Navier-Stokes equations 16:30-17:30 Free Time
65 Navier-Stokes equation u u t u 0 Smac method 1st step 2nd step 3rd step n u* u t u* q t p n 1 p n u p 1 Re u 1 Re n n n u u p u* q u ( u, v) 4th step u n 1 u* q t
66 1 st step on Smac method: Only x-direction component u t u u n Δt u u n w un Δt w + Δt Δt u Δt w w Δt convect un, Δt, u n + u u + p x 1 Re Δu = 0 + u n u n + pn x 1 Re Δu = 0 + u n u n + pn x 1 Re Δu wd = 0 un u n w + pn x w + 1 Re u w d Γ + pn x w + 1 Re u w d Γ 1 Re n w = 0 1 Re n w = 0
67 Boundary condition Γ wall Γ in Γ out Γ wall q u y( 1 y), v 0, n u 0, v 0, q 0 x q u 0, v 0, 0 n 0 on Γ in on Γ out on Γ wall
68 Problem8 (The same domain as Problem1) t + u u = p + 1 Δu in Re u = 0 in Parameters: Re = 10, Δt = 0.05 Initial condition: u = y y 5, v = 0 Boundary conditions: u = 0, v = 0, q n = 0 on Γ top [0, t] u = 0, v = 0, q n = 0on Γ bottom [0, t] u = y y 5, v = 0, q n = 0onΓ in [0, t] n = 0, v = 0, p = 0 on Γ out [0, t] u = 0, v = 0, q n = 0on Γ circle [0, t]
69 Problem9 (The same domain as Problem1) t + u u = p + 1 Δu in Re u = 0 in f + u f = Δf in t Parameters: Re = 10, Δt = 0.05 Initial condition: u = y y 5, v = 0, f = exp 10 x y Boundary conditions: u = 0, v = 0, q n = 0, f = 0 on Γ top [0, t] u = 0, v = 0, q n = 0, f = 0 on Γ bottom [0, t] u = y y 5, v = 0, q f = 0, n n = 0 onγ in [0, t] f = 0, v = 0, p = 0, n n = 0 on Γ out [0, t] u = 0, v = 0, q n = 0, f = 0 on Γ circle [0, t]
70 t + u u = p + 1 Δu in Re u = 0 in Parameters: Re = 10 5, Δt = Initial condition: u = 1, v = 0.4 Boundary conditions: n Problem10 (The same domain as Problem1) = 0, v n = 0, q = 0 on Γ top [0, t] u = 1, v = 0.4, q n = 0 on Γ bottom [0, t] u = 1, v = 0.4, q n = 0 onγ in [0, t] v = 0, n n = 0, q = 0 on Γ out [0, t] u = 0, v = 0, q n = 0 on Γ wing [0, t]
71 Borders for Problem 10 border a41(t=0.5*pi,1.5*pi){ x=1+2*cos(t); y=2*sin(t); label=1;} border a42(t=-0.5*pi,0.5*pi){ x=1+2*cos(t); y=2*sin(t); label=2;} border Splus(t=0,1){ x = t; y = *sqrt(t) *t *(t^2) *(t^3) *(t^4); label=5;} border Sminus(t=1,0){ x =t; y= -( *sqrt(t) *t *(t^2) *(t^3) *(t^4)); label=5;} int n=5; mesh Th= buildmesh( a41(30*n)+a42(30*n)+splus(60*n)+sminus(60*n));plot(th);
72 Output with Paraview Please add the following commands at the top of code load "iovtk"; Please add the following commands in the line where you want output data. String vtkoutputfile="./problem10_out"+n+".vtk"; savevtk(vtkoutputfile, Th, [u,v,0],p,dataname="velocity pressure");
An (incomplete) Introduction to FreeFem++
An (incomplete) Introduction to FreeFem++ Nathaniel Mays Wheeling Jesuit University November 12, 2011 Nathaniel Mays (Wheeling Jesuit University)An (incomplete) Introduction to FreeFem++ November 12, 2011
More informationA short introduction to FreeFem++
A short introduction to FreeFem++ M. Ersoy Basque Center for Applied Mathematics 12 November 2010 Outline of the talk Outline of the talk 1 Introduction 2 HOW TO solve steady pdes : e.g. Laplace equation
More informationNumerical simulations and mesh adaptation : exercise sessions Charles Dapogny 1, Pascal Frey 2, Yannick Privat 2
Numerical simulations and mesh adaptation : exercise sessions Charles Dapogny Pascal Frey 2 Yannick Privat 2 The purpose of this exercise session is to give an overview of several aspects of scientific
More informationFinite Elements for Fluids Prof. Sheu
Finite Elements for Fluids Prof. Sheu Yannick Deleuze http://homepage.ntu.edu.tw/~ydeleuze/ff26/ Taipei, April 28, 26 Homework from April 7, 26 Use the code from the exercices. Write a small report with
More informationElements of Matlab and Simulink Lecture 7
Elements of Matlab and Simulink Lecture 7 Emanuele Ruffaldi 12th May 2009 Copyright 2009,Emanuele Ruffaldi. This work is licensed under the Creative Commons Attribution-ShareAlike License. PARTIAL DIFFERENTIAL
More informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
More informationMATH 251 Final Examination December 16, 2015 FORM A. Name: Student Number: Section:
MATH 5 Final Examination December 6, 5 FORM A Name: Student Number: Section: This exam has 7 questions for a total of 5 points. In order to obtain full credit for partial credit problems, all work must
More informationMath Homework 3 Solutions. (1 y sin x) dx + (cos x) dy = 0. = sin x =
2.6 #10: Determine if the equation is exact. If so, solve it. Math 315-01 Homework 3 Solutions (1 y sin x) dx + (cos x) dy = 0 Solution: Let P (x, y) = 1 y sin x and Q(x, y) = cos x. Note P = sin x = Q
More informationMATH 251 Final Examination May 3, 2017 FORM A. Name: Student Number: Section:
MATH 5 Final Examination May 3, 07 FORM A Name: Student Number: Section: This exam has 6 questions for a total of 50 points. In order to obtain full credit for partial credit problems, all work must be
More informationphenomena during acupuncture Yannick Deleuze Yannick Deleuze
Modeling FreeFem++, and a PDEs simulation solver : Basic on Tutorial transport phenomena during acupuncture Yannick Deleuze Yannick Deleuze Laboratoire Jacques-Louis Lions, Sorbonne Universités, UPMC Univ
More informationM.Sc. in Meteorology. Numerical Weather Prediction
M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Centre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. In this section
More informationPH2130. Week 1. Week 2. Week 3. Questions for contemplation. Why differential equations? Why usually linear diff eq n s? Why usually 2 nd order?
PH130 week 3, page 1 PH130 Questions for contemplation Week 1 Why differential equations? Week Why usually linear diff eq n s? Week 3 Why usually nd order? PH130 week 3, page Aims of Wk 3 Lect 1 Recognise
More informationFinite element programming by FreeFem++ advanced course
1 / 37 Finite element programming by FreeFem++ advanced course Atsushi Suzuki 1 1 Cybermedia Center, Osaka University atsushi.suzuki@cas.cmc.osaka-u.ac.jp Japan SIAM tutorial 4-5, June 2016 2 / 37 Outline
More informationSTAT 430/510 Probability
STAT 430/510 Probability Hui Nie Lecture 16 June 24th, 2009 Review Sum of Independent Normal Random Variables Sum of Independent Poisson Random Variables Sum of Independent Binomial Random Variables Conditional
More informationUNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering Department of Mathematics and Statistics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4006 SEMESTER: Spring 2011 MODULE TITLE:
More informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationHomework 4 in 5C1212; Part A: Incompressible Navier- Stokes, Finite Volume Methods
Homework 4 in 5C11; Part A: Incompressible Navier- Stokes, Finite Volume Methods Consider the incompressible Navier Stokes in two dimensions u x + v y = 0 u t + (u ) x + (uv) y + p x = 1 Re u + f (1) v
More informationGETTING STARTED WITH FENICS
GETTING STARTED WITH FENICS DOUGLAS N. ARNOLD 1. A first program in FEniCS 1.1. The boundary value problem. As a first problem we consider the Neumann boundary value problem: (1) u + u = f in, u = 0 on,
More informationIntroduction to Differential Equations
Chapter 1 Introduction to Differential Equations 1.1 Basic Terminology Most of the phenomena studied in the sciences and engineering involve processes that change with time. For example, it is well known
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 13
REVIEW Lecture 12: Spring 2015 Lecture 13 Grid-Refinement and Error estimation Estimation of the order of convergence and of the discretization error Richardson s extrapolation and Iterative improvements
More information- - - - - - - - - - - - - - - - - - DISCLAIMER - - - - - - - - - - - - - - - - - - General Information: This midterm is a sample midterm. This means: The sample midterm contains problems that are of similar,
More informationElectrodynamics PHY712. Lecture 3 Electrostatic potentials and fields. Reference: Chap. 1 in J. D. Jackson s textbook.
Electrodynamics PHY712 Lecture 3 Electrostatic potentials and fields Reference: Chap. 1 in J. D. Jackson s textbook. 1. Poisson and Laplace Equations 2. Green s Theorem 3. One-dimensional examples 1 Poisson
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationMathematical Methods - Lecture 9
Mathematical Methods - Lecture 9 Yuliya Tarabalka Inria Sophia-Antipolis Méditerranée, Titane team, http://www-sop.inria.fr/members/yuliya.tarabalka/ Tel.: +33 (0)4 92 38 77 09 email: yuliya.tarabalka@inria.fr
More informationFinite Difference Methods for
CE 601: Numerical Methods Lecture 33 Finite Difference Methods for PDEs Course Coordinator: Course Coordinator: Dr. Suresh A. Kartha, Associate Professor, Department of Civil Engineering, IIT Guwahati.
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More informationThe Virtual Element Method: an introduction with focus on fluid flows
The Virtual Element Method: an introduction with focus on fluid flows L. Beirão da Veiga Dipartimento di Matematica e Applicazioni Università di Milano-Bicocca VEM people (or group representatives): P.
More informationModeling a Catalytic Converter in Comsol Multiphysics
Modeling a Catalytic Converter in Comsol Multiphysics By Jacob Harding December 10 th, 2007 Chem E 499 Problem The goal of this project was to develop a model of a catalytic converter in Comsol Multiphysics.
More informationStaple or bind all pages together. DO NOT dog ear pages as a method to bind.
Math 3337 Homework Instructions: Staple or bind all pages together. DO NOT dog ear pages as a method to bind. Hand-drawn sketches should be neat, clear, of reasonable size, with axis and tick marks appropriately
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More information10.3 Parametric Equations. 1 Math 1432 Dr. Almus
Math 1432 DAY 39 Dr. Melahat Almus almus@math.uh.edu OFFICE HOURS (212 PGH) MW12-1:30pm, F:12-1pm. If you email me, please mention the course (1432) in the subject line. Check your CASA account for Quiz
More informationVideo 15: PDE Classification: Elliptic, Parabolic and Hyperbolic March Equations 11, / 20
Video 15: : Elliptic, Parabolic and Hyperbolic Equations March 11, 2015 Video 15: : Elliptic, Parabolic and Hyperbolic March Equations 11, 2015 1 / 20 Table of contents 1 Video 15: : Elliptic, Parabolic
More information0, such that. all. for all. 0, there exists. Name: Continuity. Limits and. calculus. the definitionn. satisfying. limit. However, is the limit of its
L Marizzaa A Bailey Multivariable andd Vector Calculus Name: Limits and Continuity Limits and Continuity We have previously defined limit in for single variable functions, but how do we generalize this
More informationHomework for Math , Fall 2016
Homework for Math 5440 1, Fall 2016 A. Treibergs, Instructor November 22, 2016 Our text is by Walter A. Strauss, Introduction to Partial Differential Equations 2nd ed., Wiley, 2007. Please read the relevant
More informationAn Introduction to FreeFem++ UPMC, December 2011 O. Pantz. CMAP, Ecole Polytechnique.
An Introduction to FreeFem++ UPMC, December 2011 O. Pantz CMAP, Ecole Polytechnique. FreeFem++ What is it? FreeFEM++ is a Free software to solve PDE using the Finite Element Method. It runs on Windows,
More informationMultiple Choice Answers. MA 114 Calculus II Spring 2013 Final Exam 1 May Question
MA 114 Calculus II Spring 2013 Final Exam 1 May 2013 Name: Section: Last 4 digits of student ID #: This exam has six multiple choice questions (six points each) and five free response questions with points
More informationMATH 251 Final Examination August 14, 2015 FORM A. Name: Student Number: Section:
MATH 251 Final Examination August 14, 2015 FORM A Name: Student Number: Section: This exam has 11 questions for a total of 150 points. Show all your work! In order to obtain full credit for partial credit
More informationAugust 7, 2007 NUMERICAL SOLUTION OF LAPLACE'S EQUATION
August 7, 007 NUMERICAL SOLUTION OF LAPLACE'S EQUATION PURPOSE: This experiment illustrates the numerical solution of Laplace's Equation using a relaxation method. The results of the relaxation method
More informationMTF071 Computational Fluid Dynamics of Turbulent
CHALMERS TEKNISKA HÖGSKOLA Termo- och Fluiddynamik 42 96 Göteborg MTF07 Computational Fluid Dynamics of Turbulent Flow http://wwwtfdchalmersse/gr-kurs/mtf07 Task K2 Lars Davidson 2003-02-03 In Task K you
More informationCalculus and Parametric Equations
Calculus and Parametric Equations MATH 211, Calculus II J. Robert Buchanan Department of Mathematics Spring 2018 Introduction Given a pair a parametric equations x = f (t) y = g(t) for a t b we know how
More informationLibMesh Experience and Usage
LibMesh Experience and Usage John W. Peterson peterson@cfdlab.ae.utexas.edu and Roy H. Stogner roystgnr@cfdlab.ae.utexas.edu Univ. of Texas at Austin September 9, 2008 1 Introduction 2 Weighted Residuals
More informationMATH20411 PDEs and Vector Calculus B
MATH2411 PDEs and Vector Calculus B Dr Stefan Güttel Acknowledgement The lecture notes and other course materials are based on notes provided by Dr Catherine Powell. SECTION 1: Introctory Material MATH2411
More informationFundamental Theorems of Vector
Chapter 17 Analysis Fundamental Theorems of Vector Useful Tip: If you are reading the electronic version of this publication formatted as a Mathematica Notebook, then it is possible to view 3-D plots generated
More informationSolving Differential Equations Using MATLAB
Solving Differential Equations Using MATLAB Abraham Asfaw aasfaw.student@manhattan.edu November 28, 2011 1 Introduction In this lecture, we will follow up on lecture 2 with a discussion of solutions to
More informationMath 3150 Problems Chapter 3
Name Math 15 Problems Chapter Due date: See the internet due date. Problems are collected once a week. Records are locked when the stack is returned. Records are only corrected, never appended. Submitted
More informationMATH 251 Final Examination May 4, 2015 FORM A. Name: Student Number: Section:
MATH 251 Final Examination May 4, 2015 FORM A Name: Student Number: Section: This exam has 16 questions for a total of 150 points. In order to obtain full credit for partial credit problems, all work must
More informationMATH 251 Examination II April 7, 2014 FORM A. Name: Student Number: Section:
MATH 251 Examination II April 7, 2014 FORM A Name: Student Number: Section: This exam has 12 questions for a total of 100 points. In order to obtain full credit for partial credit problems, all work must
More informationConvergence and Error Bound Analysis for the Space-Time CESE Method
Convergence and Error Bound Analysis for the Space-Time CESE Method Daoqi Yang, 1 Shengtao Yu, Jennifer Zhao 3 1 Department of Mathematics Wayne State University Detroit, MI 480 Department of Mechanics
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) Key Concepts Current Semester 1 / 25 Introduction The purpose of this section is to define
More informationIntroduction to PARTIAL DIFFERENTIAL EQUATIONS THIRD EDITION
Introduction to PARTIAL DIFFERENTIAL EQUATIONS THIRD EDITION K. SANKARA RAO Formerly Professor Department of Mathematics Anna University, Chennai New Delhi-110001 2011 INTRODUCTION TO PARTIAL DIFFERENTIAL
More informationPart I. Discrete Models. Part I: Discrete Models. Scientific Computing I. Motivation: Heat Transfer. A Wiremesh Model (2) A Wiremesh Model
Part I: iscrete Models Scientific Computing I Module 5: Heat Transfer iscrete and Continuous Models Tobias Neckel Winter 04/05 Motivation: Heat Transfer Wiremesh Model A Finite Volume Model Time ependent
More informationEMAFF: MAgnetic Equations with FreeFem++
EMAFF: MAgnetic Equations with FreeFem++ The Grad-Shafranov equation & The Current Hole Erwan DERIAZ Bruno DESPRÉS Gloria FACCANONI š Kirill Pichon GOSTAF Lise-Marie IMBERT-GÉRARD Georges SADAKA Remy
More informationJMBC Computational Fluid Dynamics I Exercises by A.E.P. Veldman
JMBC Computational Fluid Dynamics I Exercises by A.E.P. Veldman The exercises will be carried out on PC s in the practicum rooms. Several (Matlab and Fortran) files are required. How these can be obtained
More informationSeparation of variables
Separation of variables Idea: Transform a PDE of 2 variables into a pair of ODEs Example : Find the general solution of u x u y = 0 Step. Assume that u(x,y) = G(x)H(y), i.e., u can be written as the product
More informationWithout fully opening the exam, check that you have pages 1 through 11.
MTH 33 Solutions to Final Exam May, 8 Name: Section: Recitation Instructor: INSTRUCTIONS Fill in your name, etc. on this first page. Without fully opening the exam, check that you have pages through. Show
More informationE2.5 Signals & Linear Systems. Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & 2)
E.5 Signals & Linear Systems Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & ) 1. Sketch each of the following continuous-time signals, specify if the signal is periodic/non-periodic,
More informationLecture 13: Solution to Poission Equation, Numerical Integration, and Wave Equation 1. REVIEW: Poisson s Equation Solution
Lecture 13: Solution to Poission Equation, Numerical Integration, and Wave Equation 1 Poisson s Equation REVIEW: Poisson s Equation Solution Poisson s equation relates the potential function V (x, y, z)
More informationa Write down the coordinates of the point on the curve where t = 2. b Find the value of t at the point on the curve with coordinates ( 5 4, 8).
Worksheet A 1 A curve is given by the parametric equations x = t + 1, y = 4 t. a Write down the coordinates of the point on the curve where t =. b Find the value of t at the point on the curve with coordinates
More information11.3 MATLAB for Partial Differential Equations
276 3. Generate the shape functions N (i) j = a (i) j where i =1, 2, 3,..., m and j =1, 2, 3,..., m. + b (i) j x + c(i) j y (11.2.17) 4. Compute the integrals for matrix elements α ij and vector elements
More informationName: SOLUTIONS Date: 11/9/2017. M20550 Calculus III Tutorial Worksheet 8
Name: SOLUTIONS Date: /9/7 M55 alculus III Tutorial Worksheet 8. ompute R da where R is the region bounded by x + xy + y 8 using the change of variables given by x u + v and y v. Solution: We know R is
More informationStudy of Forced and Free convection in Lid driven cavity problem
MIT Study of Forced and Free convection in Lid driven cavity problem 18.086 Project report Divya Panchanathan 5-11-2014 Aim To solve the Navier-stokes momentum equations for a lid driven cavity problem
More informationParallelism in FreeFem++.
Parallelism in FreeFem++. Guy Atenekeng 1 Frederic Hecht 2 Laura Grigori 1 Jacques Morice 2 Frederic Nataf 2 1 INRIA, Saclay 2 University of Paris 6 Workshop on FreeFem++, 2009 Outline 1 Introduction Motivation
More informationENGI Parametric Vector Functions Page 5-01
ENGI 3425 5. Parametric Vector Functions Page 5-01 5. Parametric Vector Functions Contents: 5.1 Arc Length (Cartesian parametric and plane polar) 5.2 Surfaces of Revolution 5.3 Area under a Parametric
More informationSpecialist Mathematics 2017 Sample paper
South Australian Certificate of Education The external assessment requirements of this subject are listed on page 20. Question Booklet 2 Specialist Mathematics 2017 Sample paper 2 Part 2 (Questions 10
More informationModule 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 1: Finite Difference Method
file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture1/1_1.htm 1 of 1 6/19/2012 4:29 PM The Lecture deals with: Classification of Partial Differential Equations Boundary and Initial Conditions Finite Differences
More informationAn Introduction to the Discontinuous Galerkin Method
An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationMATH 251 Examination II April 4, 2016 FORM A. Name: Student Number: Section:
MATH 251 Examination II April 4, 2016 FORM A Name: Student Number: Section: This exam has 12 questions for a total of 100 points. In order to obtain full credit for partial credit problems, all work must
More informationTMA4220: Programming project - part 1
TMA4220: Programming project - part 1 TMA4220 - Numerical solution of partial differential equations using the finite element method The programming project will be split into two parts. This is the first
More informationNumerical Optimization Algorithms
Numerical Optimization Algorithms 1. Overview. Calculus of Variations 3. Linearized Supersonic Flow 4. Steepest Descent 5. Smoothed Steepest Descent Overview 1 Two Main Categories of Optimization Algorithms
More informationAP Calculus (BC) Chapter 10 Test No Calculator Section. Name: Date: Period:
AP Calculus (BC) Chapter 10 Test No Calculator Section Name: Date: Period: Part I. Multiple-Choice Questions (5 points each; please circle the correct answer.) 1. The graph in the xy-plane represented
More informationPractice Problems for Exam 3 (Solutions) 1. Let F(x, y) = xyi+(y 3x)j, and let C be the curve r(t) = ti+(3t t 2 )j for 0 t 2. Compute F dr.
1. Let F(x, y) xyi+(y 3x)j, and let be the curve r(t) ti+(3t t 2 )j for t 2. ompute F dr. Solution. F dr b a 2 2 F(r(t)) r (t) dt t(3t t 2 ), 3t t 2 3t 1, 3 2t dt t 3 dt 1 2 4 t4 4. 2. Evaluate the line
More informationMATH 124. Midterm 2 Topics
MATH 124 Midterm 2 Topics Anything you ve learned in class (from lecture and homework) so far is fair game, but here s a list of some main topics since the first midterm that you should be familiar with:
More informationLecture 2: Fundamentals 2/3: Navier-Stokes Equation + Basics of Discretization
Lecture 2: Fundamentals 2/3: Navier-Stokes Equation Basics of Discretization V.Vuorinen Aalto University School of Engineering CFD Course, Spring 2018 January 15th 2018, Otaniemi ville.vuorinen@aalto.fi
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for solving higher order ODEs and
More informationFEniCS Course. Lecture 0: Introduction to FEM. Contributors Anders Logg, Kent-Andre Mardal
FEniCS Course Lecture 0: Introduction to FEM Contributors Anders Logg, Kent-Andre Mardal 1 / 46 What is FEM? The finite element method is a framework and a recipe for discretization of mathematical problems
More informationConduction and Radiation Prof. C. Balaji Department of Mechanical Engineering Indian Institute of Technology, Madras
Conduction and Radiation Prof. C. Balaji Department of Mechanical Engineering Indian Institute of Technology, Madras Lecture No. # 11 Emissivity Contd So, we will continue the discussion on definition
More informationECE 301 Division 1 Exam 1 Solutions, 10/6/2011, 8-9:45pm in ME 1061.
ECE 301 Division 1 Exam 1 Solutions, 10/6/011, 8-9:45pm in ME 1061. Your ID will be checked during the exam. Please bring a No. pencil to fill out the answer sheet. This is a closed-book exam. No calculators
More informationLecture 16: Relaxation methods
Lecture 16: Relaxation methods Clever technique which begins with a first guess of the trajectory across the entire interval Break the interval into M small steps: x 1 =0, x 2,..x M =L Form a grid of points,
More informationNAME: ht () 1 2π. Hj0 ( ) dω Find the value of BW for the system having the following impulse response.
University of California at Berkeley Department of Electrical Engineering and Computer Sciences Professor J. M. Kahn, EECS 120, Fall 1998 Final Examination, Wednesday, December 16, 1998, 5-8 pm NAME: 1.
More informationIntroduction into Implementation of Optimization problems with PDEs: Sheet 3
Technische Universität München Center for Mathematical Sciences, M17 Lucas Bonifacius, Korbinian Singhammer www-m17.ma.tum.de/lehrstuhl/lehresose16labcourseoptpde Summer 216 Introduction into Implementation
More informationMA3D9. Geometry of curves and surfaces. T (s) = κ(s)n(s),
MA3D9. Geometry of 2. Planar curves. Let : I R 2 be a curve parameterised by arc-length. Given s I, let T(s) = (s) be the unit tangent. Let N(s) be the unit normal obtained by rotating T(s) through π/2
More informationMath 333 Exam 1. Name: On my honor, I have neither given nor received any unauthorized aid on this examination. Signature: Math 333: Diff Eq 1 Exam 1
Math 333 Exam 1 You have approximately one week to work on this exam. The exam is due at 5:00 pm on Thursday, February 28. No late exams will be accepted. During the exam, you are permitted to use your
More informationSolving Differential Equations on 2-D Geometries with Matlab
Solving Differential Equations on 2-D Geometries with Matlab Joshua Wall Drexel University Philadelphia, PA 19104 (Dated: April 28, 2014) I. INTRODUCTION Here we introduce the reader to solving partial
More informationChapter 2 Boundary and Initial Data
Chapter 2 Boundary and Initial Data Abstract This chapter introduces the notions of boundary and initial value problems. Some operator notation is developed in order to represent boundary and initial value
More informationMath 106 Answers to Exam 3a Fall 2015
Math 6 Answers to Exam 3a Fall 5.. Consider the curve given parametrically by x(t) = cos(t), y(t) = (t 3 ) 3, for t from π to π. (a) (6 points) Find all the points (x, y) where the graph has either a vertical
More informationNATIONAL UNIVERSITY OF SINGAPORE FACULTY OF SCIENCE SEMESTER 1 EXAMINATION Nov/Dec 2009 Time allowed :
NATIONAL UNIVERSITY OF SINGAPORE FACULTY OF SCIENCE SEMESTER 1 EXAMINATION 2009-2010 Mathematics II Nov/Dec 2009 Time allowed : 2 hours INSTRUCTIONS TO CANDIDATES 1 This examination paper contains SIX
More informationLecture 10: Finite Differences for ODEs & Nonlinear Equations
Lecture 10: Finite Differences for ODEs & Nonlinear Equations J.K. Ryan@tudelft.nl WI3097TU Delft Institute of Applied Mathematics Delft University of Technology 21 November 2012 () Finite Differences
More informationMath 116 Practice for Exam 2
Math 116 Practice for Exam 2 Generated October 12, 215 Name: SOLUTIONS Instructor: Section Number: 1. This exam has 8 questions. Note that the problems are not of equal difficulty, so you may want to skip
More information1 Comparison of Kansa s method versus Method of Fundamental Solution (MFS) and Dual Reciprocity Method of Fundamental Solution (MFS- DRM)
1 Comparison of Kansa s method versus Method of Fundamental Solution (MFS) and Dual Reciprocity Method of Fundamental Solution (MFS- DRM) 1.1 Introduction In this work, performances of two most widely
More informationThis exam is closed book and closed notes. (You will have access to a copy of the Table of Common Distributions given in the back of the text.
TEST #3 STA 536 December, 00 Name: Please read the following directions. DO NOT TURN THE PAGE UNTIL INSTRUCTED TO DO SO Directions This exam is closed book and closed notes. You will have access to a copy
More informationUpdated 2013 (Mathematica Version) M1.1. Lab M1: The Simple Pendulum
Updated 2013 (Mathematica Version) M1.1 Introduction. Lab M1: The Simple Pendulum The simple pendulum is a favorite introductory exercise because Galileo's experiments on pendulums in the early 1600s are
More informationLecture 18 Finite Element Methods (FEM): Functional Spaces and Splines. Songting Luo. Department of Mathematics Iowa State University
Lecture 18 Finite Element Methods (FEM): Functional Spaces and Splines Songting Luo Department of Mathematics Iowa State University MATH 481 Numerical Methods for Differential Equations Draft Songting
More informationMA424, S13 HW #6: Homework Problems 1. Answer the following, showing all work clearly and neatly. ONLY EXACT VALUES WILL BE ACCEPTED.
MA424, S13 HW #6: 44-47 Homework Problems 1 Answer the following, showing all work clearly and neatly. ONLY EXACT VALUES WILL BE ACCEPTED. NOTATION: Recall that C r (z) is the positively oriented circle
More informationMath 308 Week 8 Solutions
Math 38 Week 8 Solutions There is a solution manual to Chapter 4 online: www.pearsoncustom.com/tamu math/. This online solutions manual contains solutions to some of the suggested problems. Here are solutions
More informationVolumes of Solids of Revolution Lecture #6 a
Volumes of Solids of Revolution Lecture #6 a Sphereoid Parabaloid Hyperboloid Whateveroid Volumes Calculating 3-D Space an Object Occupies Take a cross-sectional slice. Compute the area of the slice. Multiply
More information