资源说明:自己写了一个二维Euler方程的间断有限元程序
上次发了一个三角形单元的程序 因为不是曲边单元 所以在圆柱后面容易形成涡
现在把程序改为曲边四边形单元了 没有涡出现
单元是8节点四边形单元 节点编号顺序是 1 5 2 6 3 7 4 8 也就是四个角点依次
是1 2 3 4 然后是边的中点编号 5 6 7 8.
时间推进采用 Runge-Kutta 方法
数值通量采用全局Lax-Friedrichs通量
仍然不能捕捉激波 因为没有做重构或者加人工粘性 等这个做出来了 再发一次。
程序没有进行优化 比如说内存的消耗没有优化 比如直边单元的边界积分仍然采
用了曲边的积分方法 增加了计算量 比如面积分、线积分都是采用的是Gauss-
Legendre-Lobatto积分 积分精度会比一般的Gauss-Legendre积分精度低一阶 等
等问题。 二维的 纯属交流性质 就没有考虑这些问题 ^_^
如果物面全部是直边 那么只要改变一个参数N 就可以获得不同的计算精度 且具
有谱精度 因为单元的节点是Gauss-Legendre-Lobatto积分点。 其实就是谱元法
(物面是曲边的情况我不清楚是不是也可以通过提高基函数的阶数 也就是增加N
来提高计算精度)
本源码包内暂不包含可直接显示的源代码文件,请下载源码包。
-
- gmp-5.0.1.tar.gzGMP是The GNU MP Bignum Library,是一个开源的数学运算库,它可以用于任意精度的数学 ...
- fretrt.rar所有满足下列条件的分数, 1)它大于1/6,小于1/5; 2)分子,分母都是素数; 3 ...
- check.rar有15个数按由小到大顺序放在一个数组中,输入一个数,要求用折半查找法找出该数是数组 ...
- fft.rarFast fourier transform FFT
- idct.rarIDCT Discrete cossine transform
- SignalProcess.rar基于Visual C++的地震信号处理程序源代码,包括常用的FFT和Hilbert变换代码。 ...
- study.zip随机数的一个程序,里面有组数,每个随机数为8位16进制数字,随机数保存在程序目录下 ...
- dafsdf.rarFFT算法 C++ 快速傅里叶变换 逆傅里叶变换
- Kalman1960.rarkalman滤波的应用,在电力系统中的相关应用
- FFT_1.zipfft test source in gcc
-
- pca.rarPCA算法程序设计步骤: 1、去均值 2、计算协方差矩阵及其特征值和特征向量 3、计算协 ...
- C语言精彩编程百例_源代码.rarC语言精彩编程百例_源代码
- FFT2Dsource.rar二维FFT(快速傅立叶变换)变换源代码,可编译通过并能正常运行. ...
- VC++cac.rarVisual C++常用数值算法集 本收共不数值计算中常用的Visual C++子过程近200个,内容包 ...
- 科学与工程数值计算算法配套源码vc++.zip清华大学2002年出版的《科学与工程数值计算算法VC++》配套源码,包括复数计算、矩阵 ...
- 矩阵特征值的计算.rar矩阵特征值与特征向量的计算约化对称矩阵为对称对角阵的豪斯荷尔镕变换法似对称三对角 ...
- 求解线性代数方程组.rar线性代数方程组的求解全选主元高斯消去法全选主元高斯约当消去法复系数方程组的全选主 ...
- 矩阵的运算.rar矩阵运算实矩阵相乘复矩阵相乘实矩阵求逆的全选主元高斯约当法 ...
- AR(5).rar利用AR模型进行时间序列预测的程序源代码,使用最小二乘估计法进行参数估计。拟合效果 ...
- 实数据快速傅立叶变换算法.rar傅立叶变换是信号处理时最常用的算法之一,实现时域与频域之间的转换,对人们研究各种 ...