Advanced Programming in Engineering Saurabh Srivastava, Thomas Weinhart Martin Robinson, and Stefan Luding MSM, CTW, UTwente, NL

Size: px

Start display at page:

Download "Advanced Programming in Engineering Saurabh Srivastava, Thomas Weinhart Martin Robinson, and Stefan Luding MSM, CTW, UTwente, NL"

Myron Greene
5 years ago
Views:

Advanced Programming in Engineering 2011-19115850 Saurabh Srivastava, Thomas Weinhart Martin Robinson, and Stefan Luding MSM, CTW,

1 Advanced Programming in Engineering Saurabh Srivastava, Thomas Weinhart Martin Robinson, and Stefan Luding MSM, CTW, UTwente, NL Contents Why is this class important for me? What will I learn in this class? When is this class? Where is this class? 1

! Computations are everywhere in engineering

Many problems are resolved with the aid of

2 Introduction! Why is the class important for me?! Computations are everywhere in engineering problems. Many problems are resolved with the aid of computers and dedicated programs today. It is important to be able to implement numerical algorithms It is really important for an engineer to be familiar with computers and programming languages. Introduction! 2

Introduction! What can we do with the computer?

do we do that?! Ø Compiler languages (C, C++, Fortran,.

Ø Translate problems to algorithms Ø Practical experiences (exercises,debugging ) Introduction! Program:!

3 Introduction! What can we do with the computer? Ø Evaluation of (experimental) data Ø Solving physical/engineering problems Ø Numerical experiments How do we do that?! Ø Compiler languages (C, C++, Fortran,...) Ø Interpreter languages (MATLAB, ) What do we learn? Ø Translate problems to algorithms Ø Practical experiences (exercises,debugging ) Introduction! Program:! Algorithms to solve differential equations Finite element method Molecular dynamics with Lennard-Jones for fluids Measuring pressure, temperature and diffusion in MD! Random numbers & Fractals Monte-Carlo for solving integrals and PDEs Finite Volume SPH 3

4 Introduction! When is the class? From 14 November Monday Morning ½ Theory Monday Morning ¾ Practice Introduction! When is the class? From 14 November Monday Morning ½ Theory Monday Morning ¾ Practice Stefan Ordinary Differential Equations Thomas Debugging, object orientation, profiling, Stefan Molecular Dynamics for Solids Saurabh FEM for Solids Saurabh FEM for Solids Saurabh Nonlinear FEM Stefan Random Numbers and Applications reserve 4

Introduction! When is the class? From 14 November Monday Morning ½ 8.45-10.

5 Introduction! When is the class? From 14 November Monday Morning ½ Theory Monday Morning ¾ Practice Thomas Martin Martin Thomas Thomas MD for Fluids and Statistical Analysis SPH (Smooth Particle Hydrodynamics) 1 SPH (Smooth Particle Hydrodynamics) 2 free day FV (Free Volume) Methods 1 FV (Free Volume) Methods 2 reserve reserve Introduction! Where is the class? OH 112 5

6 Introduction! Questions? Flow with friction & rolling resistance µ = 0.5 µ = 0.5 µ = 0.2 r 6

7 Sintering Vibration test p=100 p=10 tension kt k 2 = 12 7

8 P-wave animation P-wave animation 8

9 Anisotropy 3D 3D 3-dimensional modeling of sound propagation Sound P-wave shape and speed 9

10 Preparation p=20000 p=2000p=200 p=20 p=2 Contents Introduction Examples 1&2 1st Days Goal: Solve Differential Equations (ODE) 2nd Days Goal: Debugging and Matlab Optimization 3rd Days Goal: 1D/2D Molecular Dynamics (ODE) Morning 3/4: Practical Exercises 10

11 Introduction to Advanced Programming in Engineering Examples 1 (Stefan) Saurabh Srivastava, Thomas Weinhart, Stefan Luding MSM, TS, CTW, UTwente, NL 5. int main(int argc, char *argv[]) 6. { 7. const int ipmax=20; 8. int inum; 9. cout << Type number: ; 10. cin >> inum; if(inum < pow(2.0,ipmax)) // check number 13. { 14. cout << Base 10: << inum << endl; 15. cout << Base 2: ; // perform binary check for 2^ipmax 18. for(int i=ipmax; i>=0; i--) 19. { 20. if(inum >= pow(2.0,i)) 21. { 22. inum -= pow(2.0,i); 23. cout << '1'; 24. } 25. else 26. cout << '0'; 27. } 28. cout << endl; 29. } 30. else 31. { 32. cout << ERROR: input-number << inum << > 33. << pow(2.0,ipmax) << too large! << endl; 34. } 35. system(pause); 36. return EXIT_SUCCESS; 37. } Previous Course PiE: Exercise 2 0 no INPUT Loop: i=20,,0 inum>=2 i 1 END yes inum=inum 2 i 11

12 Differential Equations Physical examples Mass-spring system Pendulum Methods Euler and Euler-Cromer Verlet and friends Runge-Kutta Predictor-Corrector etc. 1. #include<iostream> 2. #include<fstream> 3. #include<cmath> 4. using namespace std; 5. int main(int argc, char *argv[]) 6. { 7. // Define field x(t) with length double x[1000], t; // output variables 9. // initial conditions 10. double A, delta; // A=ampl, delta=phase 11. double mass; // mass=mass 12. double ksprng; // ksprng=spring-const. 13. double t_max, dt; // t_max=max-time 14. // dt=time-interval 15. // Request input of the parameters 16. cout << amplitude ; cin >> A; 17. cout << phase-angle ; cin >> delta; 18. cout << mass ; cin >> mass; 19. cout << spring-const. ; cin >> ksprng; 20. cout << max-time ; cin >> t_max; 21. cout << time-interval ; cin >> dt; Exercise double omega=sqrt(ksprng/mass); 23. // Loop from t=0 to t=t_max 24. t=0.0; 25. for( int i=0; i<=t_max/dt; i++ ) 26. { 27. // compute function 28. x[i]=a*sin(omega*t+delta); 29. t=t+dt; // Step to next time 30. } 31. ofstream outfile(plot.data); 32. // Output-start 33. t=0.0; 34. for( int i=0; i<=t_max/dt; i++ ) 35. { 36. outfile << t << 37. << x[i] << endl; 38. t=t+dt; 39. } 40. } 12

13 Differential Equations Mass-spring system Differential Equations 1 Mass spring system 13

14 Differential Equations 1 Mass-spring system Differential Eqs Mass-spring system 14

15 Differential Equations 2 Pendulum Differential Equations Physical examples Mass-spring system Pendulum Methods Euler and Euler-Cromer Verlet and friends Runge-Kutta Predictor-Corrector etc. 15

16 Differential Equations 2 Pendulum Differential Equations 16

17 Differential Eqs Poincare-cut instead of x-t-plot view: v-x-plot attractors/chaos Differential Equations - summary Physical examples Mass-spring system Pendulum Methods Euler and Euler-Cromer Verlet and friends Runge-Kutta Predictor-Corrector etc. 17

18 Questions? 08:45h 10:30h room: OH112! 10:45h 12:30h room: OH112! Advanced Programming in Engineering ! Saurabh Srivastava, Thomas Weinhart, Martin Robinson, Stefan Luding 18

Simulation in Computer Graphics. Particles. Matthias Teschner. Computer Science Department University of Freiburg

Simulation in Computer Graphics Particles Matthias Teschner Computer Science Department University of Freiburg Outline introduction particle motion finite differences system of first order ODEs second