Changes between Version 13 and Version 14 of fortran
- Timestamp:
- Oct 18, 2016, 9:28:34 PM (8 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
fortran
v13 v14 177 177 178 178 == 7. Electron motion only in magnetic field == 179 The main idea in plasma physics is to understand the motion of the charge particles. For that reason first we will make a solution of the electron motion but only in magnetic filed. Solving this problem by using the equation of electron motion can be obtain the phase trajectory. Electron filed E0=0 180 Part 1. Using only dependence of B0 or initial magnetic filed 179 181 {{{ 180 182 #!fortran 181 !The main idea in plasma physics is to understand the motion of the charge particles. For182 !that reason first we will make a solution of the electron motion but only in magnetic183 !filed. Solving this problem by using the equation of electron motion can be obtain the184 !phase trajectory. Electron filed E0=0185 !Part 1. Using only dependence of B0 or initial magnetic filed186 183 integer :: n=1000 187 184 real(8) :: dt, t, vm, b, rl, om, x, y, vx, vy, B0 … … 225 222 226 223 == 8. Electron motion in magnetic filed. == 224 The introduction is like in program 7. The difference is to solve the problem how will be the electron motion if the velocities vx and vy depends of the magnetic filed Bz which is a function of the coordinate y. E0=0. Part 2 227 225 {{{ 228 226 #!fortran 229 !The introduction is like in program 7. The difference is to solve the problem how will be230 !the electron motion if the velocities vx and vy depends of the magnetic filed Bz which is a231 !function of the coordinate y. E0=0232 !Part 2233 227 integer :: n=1000 234 228 real(8) :: dt, t, vm, om, x, y, vx, vy, B0,b … … 297 291 }}} 298 292 == 9. Electron motion in electric and magnetic filed == 293 To make a conclusion how is the electric motion looks like using plus electric filed. And to get a picture of the motion in the plasma in this program we are using different functions of dependencies for vx and vy. And the electric filed is E0=3. The particle begins to drift. 299 294 {{{ 300 295 #!fortran 301 !To make a conclusion how is the electric motion looks like using plus electric filed. And to302 !get a picture of the motion in the plasma in this program we are using different functions303 !of dependencies for vx and vy. And the electric filed is E0=3. The particle begins to drift.304 305 296 integer :: n=1000 306 297 real(8) :: dt, t, vm, om, x, y, vx, vy, B0,b … … 356 347 }}} 357 348 == 10. Electron motion only in magnetic field == 349 The introduction is like in program 8, but the difference is that here we have 350 trigonometrical dependencies of the initial velocities. Part 3 358 351 {{{ 359 352 #!fortran 360 !The introduction is like in program 8, but the difference is that here we have361 !trigonometrical dependencies of the initial velocities.362 !Part 3363 353 integer :: N = 1500 364 354 real(8) :: x = 2.5, y = 0.5, B0 = 500.0, We = 4.0054e-08, b = 2.0, pi = 3.14159, m = 9.1e-28, e = -4.8e-10, c = 2.99e10 … … 397 387 398 388 == 11. Determination the proton energy and power == 389 Also for one of the main characteristics in the plasma is the energy and the power of the charge particles. In this program are obtained the energy and the power of the protons in the plasma. Using the basics equations for it. 399 390 {{{ 400 391 #!fortran 401 !Also for one of the main characteristics in the plasma is the energy and the power of the402 !charge particles. In this program are obtained the energy and the power of the protons in403 !the plasma. Using the basics equations for it.404 392 integer :: n=1000 405 393 real(8) :: dt, t, x, y, vx, vy, v, E, W … … 445 433 }}} 446 434 == 12. Runge-Kutta method for harmonic oscillations == 435 To solve a partial differential equations need to satisfy the condition y(x0)=y0 which is called Cauchy task. The most effective and the commonly used method of the solution of Cauchy task is Runge-Kuta method. It is based on a approximation of the function y, which is obtained by expanding the function in Taylor series. 436 Using method Runge-Kuta can be obtain the basic differential equation for oscillator. Part 1 437 447 438 {{{ 448 439 #!fortran 449 !To solve a partial differential equations need to satisfy the condition y(x0)=y0 which is450 !called Cauchy task. The most effective and the commonly used method of the solution of451 !Cauchy task is Runge-Kuta method. It is based on a approximation of the function y, which is452 !obtained by expanding the function in Taylor series.453 !Using method Runge-Kuta can be obtain the basic differential equation for oscillator.454 !Part 1455 456 457 440 external fct, out 458 441 real aux(8,2)