How do i write the code for navierstokes equation in matlab/comsol?

Question

Bartholomew Osegbe 2018-5-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/400044-how-do-i-write-the-code-for-navierstokes-equation-in-matlab-comsol

回答： Bartholomew Osegbe 2018-5-13

The idea is to write the code for the navierstokes equation with the laplacian velocity and pressure term

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

KSSV 2018-5-10

You want it in comsol? So this is not the forum. You want in MATLAB? What have you attempted?

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Precise Simulation 2018-5-10

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/400044-how-do-i-write-the-code-for-navierstokes-equation-in-matlab-comsol#answer_319563

The incompressible Navier-Stokes equations is also available as a built-in pre-defined Navier-Stokes physics mode in the FEATool FEM Matlab toolbox. In this case the equations are in 2D defined as

rho_ns*u' - miu_ns*(2*ux_x + uy_y + vx_y) + rho_ns*(u*ux_t + v*uy_t) + p_x = Fx_ns
rho_ns*v' - miu_ns*(vx_x + uy_x + 2*vy_y) + rho_ns*(u*vx_t + v*vy_t) + p_y = Fy_ns
ux_t + vy_t = 0

the Navier-Stokes equations in axisymmetry (2D cylindrical coordinates) and 3D are defined similarly.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Bartholomew Osegbe 2018-5-10

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/400044-how-do-i-write-the-code-for-navierstokes-equation-in-matlab-comsol#answer_319540

I have attempted the weak method fem.equ.weak = { ... { ... {'(-eta*ux)*test(ux)+(-eta*uy)*test(uy)-test(u)*px'} ... {'(-eta*vx)*test(vx)+(-eta*vy)*test(vy)-test(v)*py'} ... {'-test(px)*u-test(py)*v'} ... } ... }; but I have errors and when I sourced online I found a code that uses the poisson equation but just wanted a simple code with a laplacian velocity and pressure term.