Unsteady convection-diffusion-reaction problem

版本 3.3 (18.9 KB) 作者: Andrea La Spina
1D scalar equation of a convection-diffusion-reaction problem with piecewise linear approximation
919.0 次下载
更新时间 2021/10/13

查看许可证

The "UNSTEADY_CONVECTION_DIFFUSION_REACTION_1D" script solves the 1D scalar equation of a convection-diffusion-reaction problem with piecewise linear approximation.
For the space discretization the user can choose the standard Galerkin or the Full Upwind approach.
For the time integration the theta-method has been implemented. According to the value of theta these schemes are obtained:
0 -> Forward Euler
1/2 -> Crank-Nicolson
3/4 -> Galerkin
1 -> Backward Euler
The FEM parameters such as the number of finite elements and the number of Gauss integration points can be easily chosen.
The the boundary conditions can be expressed a function of time and a nonlinear source term s(x) can be taken into account.
The functions and the examples are developed according with Chapter 5 "Unsteady convection-diffusion
problems" of the book "Finite Element Methods for flow problems" of Jean Donea and Antonio Huerta.
If you like the file, please give feedbacks.

引用格式

Andrea La Spina (2024). Unsteady convection-diffusion-reaction problem (https://www.mathworks.com/matlabcentral/fileexchange/60623-unsteady-convection-diffusion-reaction-problem), MATLAB Central File Exchange. 检索来源 .

MATLAB 版本兼容性
创建方式 R2014a
兼容任何版本
平台兼容性
Windows macOS Linux

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!
版本 已发布 发行说明
3.3

.

3.2

.

3.1.0.0

.

3.0.0.0

A nice animation has been included.

2.0.0.0

A nice animation has been included in the script.

1.0.0.0