Transient Heat equation: comparison between CA methods and FE methods

This pages shows the resolution of the 2-D transient heat equation: Ft+Kx*Fxx+Ky*Fyy=0, this The boundary conditions are chosen to be two sinoidal functions in the upper and right edges, the initial condition is chosen to be an exponential function, they are shown in the following graphics:


Initial Condition

Conditions for the problem

The following applet shows the solution found by both models.

Theory pages:
Main page
FEM (i)
FEM (ii)
FDM (i)
Example pages:
1-d Heat Equation
2-d steady state Heat Equation
2-d Heat Equation
1-d non diffusive transport equation
1-d diffusive transport equation
2-d Non diffusive transport equation
Mesh generation with OOCSMP
Moving grids
Application pages:
Heating of two beams
Heating of two moving beams
Solving the equation Ut+Uxx+Uxy+Uyx=0
Solving the equation Ut+Uxx+Uxy+Uyx=0 using MGEN
Heat 1d using several outputs
Solving the Heat equation with a CA
Comparing a CA with the FEM
Gordon's sine equation

Other courses
Other pages

Last modified 16/02/2000 by Juan de Lara (, need help for using this courses?.

The OOCSMP code - The SODA code