精确算法matlab代码
Problem
Description
非齐次的三类边值条件的双边值问题
已知精确解，p=1
,q=0,f=sinπx+cosπx
从Ritz(有限元)法出发，使用Matlab编写代码解决该问题。
M:
区间大小
参考文献请见：两点边值问题的有限元算法_付小龙.pdf
基函数没有选择好
致使效果误差较大，如下图所示:
若将该问题转为经典的一类边值条件的双边值问题，则按照传统的Rizt解法
可得到如下:
2023-04-07 00:00:37
1.29MB
1
该代码用于计算基于有限不同时域的 3D 电磁场。 在数值过程中采用一阶MUR边界条件。
2023-03-12 22:03:14
4KB
1
基于零内向流边界条件的辐射传输方程是目前被广泛采用的描述有限组织体中光子传播过程的模型.由于 生物组织体及其周围介质之间折射率的差异,零内向流边界条件无法准确模拟光子到达边界时的传播过程,从而导致了实际应用中对光学参数的重建存在明显误差.为解决这一问题,建立了考虑边界内反射效应的自然边界条件下的辐射传输方程模型,发展了一种结合立体角元离散法和有限差分法的自然边界条件稳态辐射传输方程求解算法,并通过与蒙特卡洛模拟结果的对比验证了该算法的正确性.
2023-03-08 11:17:51
1.09MB
1
入口和出口条件
描述流体流入和流出的边界条件类型:
通用的
压力入口
压力出口
不可压缩流动
速度入口
出流条件
根据不同的物理模型需要设定不同的边界参数.
指导方针:
有流体流入或流出的位置.
有利于收敛.
在边界方向避免出现过大的梯度.
表示设置错误.
减小边界上网格的斜度.
可压缩流动
质量流入口
压力远场条件
特殊条件
通风入口, 通风出口,
进气风扇, 排气风扇
2023-02-20 22:17:32
651KB
1
在回顾了SPH无滑移边界条件施加方法的基础上, 针对Ellero等人提出的无滑移边界处理方法, 结合两次导数求解粘性项, 分析了利用Poiseuille流模型施加这种方法时会遇到的不稳定问题, 包括数值误差引起的不稳定问题和边界粒子场变量更新方式引起的不稳定问题, 提出了相应的解决方案并对其方法进行了一定的改进, 获得了更加精确的计算结果。最后尝试了一种与改进边界处理方法等效的方法, 并利用Poiseuille流模型进行了验证。
2023-02-16 10:32:07
1.49MB
1