FreeFEM++: Solve a parabolic equation
- We seek the temperature distribution in a plate of rectangular cross section Ω = ( − 1 , 1 ) × ( − 1 , 1 ) \Omega=(-1,1) \times(-1,1) Ω=(−1,1)×(−1,1);
- We must solve the temperature equation in Ω \Omega Ω in a time interval ( 0 , T ) (0, \mathrm{~T}) (0, T).
∂ t u − Δ u = 0 in Ω × ( 0 , T ) u ( x , y , 0 ) = u 0 ∂ u ∂ n = g n on Γ 2 , Γ 4 × ( 0 , T ) u = g d on Γ 1 , Γ 3 × ( 0 , T ) \begin{aligned} \partial_{t} u-\Delta u &=0 & & \text { in } \Omega \times(0, T) \\ u(x, y, 0) &=u_{0} & & \\ \frac{\partial u}{\partial n}&=g_n & & \text { on } {\Gamma_2,\Gamma_4} \times(0, T) \\ u&=g_d & & \text { on } {\Gamma_1,\Gamma_3} \times(0, T) \end{aligned} ∂tu−Δuu(x,y,0)∂n∂uu=0=u0=gn=gd in Ω×(0,T) on Γ2,</