Ticket #3755: advectionCos.mo

File advectionCos.mo, 447 bytes (added by Jan Silar, 10 years ago)
Line 
1model advectionSin "advection equation"
2 parameter Real L = 1;
3 parameter DomainLineSegment1D omega(L = L, N = 100);
4 field Real u(domain = omega);
5 parameter Real c = 1;
6initial equation
7 u = if omega.x < 0.25 then cos(2*3.14159*omega.x) else 0 indomain omega;
8equation
9 der(u) + c*pder(u,x) = 0 indomain omega;
10 u = 1 indomain omega.left;
11 u = 0 indomain omega.right;
12end advectionSin;