-
practice6_17
用矩阵法求定积分1/(1+x^2)dx,区间是(-4,4)(Definite Integral Matrix Method)
- 2010-05-17 01:34:47下载
- 积分:1
-
Euler_DG_Quadrilateral_2D
自己写了一个二维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
来提高计算精度)(Wrote a two-dimensional Euler equations with discontinuous finite element program
Last made a triangular element of the program, not curved edge unit is so easy to form a vortex in the cylinder behind the
Program to curved edge quadrilateral element vortices appear
Unit is the order of 8-node quadrilateral element node number is 15,263,748 which is the four corners of the points in turn
Is 1234 and then the side of the midpoint of the number 5678.
Time promote the use of Runge-Kutta method
Numerical flux of the overall situation of Lax-Friedrichs, flux
Still can not capture the shock wave did not do the reconstruction or artificial viscosity do it Zaifayici.
The program is not optimized for example, memory consumption is not optimized such as straight-edge boundary integral of the unit is still mining
Integral method to increase the amount of computation such as surface integral with a curved edge, the line integral using the Gauss-
The Legendr)
- 2021-01-29 12:58:40下载
- 积分:1
-
gain saturation in mode locked laser
用于计算锁模激光腔中增益饱和强度,用于锁模激光(For calculating the gain saturation intensity in a mode-locked laser cavity for a mode locke)
- 2017-11-30 10:42:02下载
- 积分:1
-
SHAKE91
国外的地震有限元程序,美国伯克利大学的,非常经典(Earthquake abroad finite element program, the United States Berkeley, very classic)
- 2021-03-10 01:09:27下载
- 积分:1
-
The-UDF-for-phase-change-of-water
Fluent 中短语模拟水的受热蒸发相变的UDF程序代码(the UDF for phase change of water)
- 2014-01-06 19:38:21下载
- 积分:1
-
energy-leakage--dual-tree
首先根据高斯白噪声频率充满整个频带的特性,通过双树复小波包变换对高斯白噪声进行分解,利用频带能量泄漏的定量分析方法,验证了双树复小波包变换具有较低的频带能量泄漏特性;其次利用双树复小波包变换逐层分解信号,对每层分解所得分量求其FFT谱的峭度,得到基于双树复小波包变换的谱峭度图,根据图中峭度最大的原则,可以自动准确的选择信号分解最佳层数和最佳分量;最后将基于双树复小波包变换的谱峭度图的故障诊断方法应用于实际工程中,对齿轮故障振动信号进行分析,选择最佳分解层数和分量后利用希尔伯特包络解调,有效准确地提取了故障特征信息,验证了方法的可行性和有效性(The parameters of a filter were determined by experience, and that has a great influence on the
results of signal processing. The discrete wavelet packet transform has a larger energy leakage of frequency band,
which obviously affected the results of the envelope demodulation. It is necessary to have a method with a lower
energy leakage of the frequency band before envelope demodulation. The dual tree complex wavelet packet
transform (DT-CWPT) was a new signal processing method that had many good qualities. Because the energy
leakage of the frequency band was smaller when the signal was decomposed by a dual tree complex wavelet
packet transform, the dual tree complex wavelet packet transform was used to extract the fault feature information
in the field of fault diagnosis. In this paper, first, according to the characteristics of Gaussian white noise, whose
frequency was full of the whole frequency band, Gaussian white noise was decomposed by a dual-tree complex
wavel)
- 2014-05-28 17:13:41下载
- 积分:1
-
CSAMT_1D_forword
一维CSAMT正演程序,对学地球物理的有帮助(CSAMT)
- 2021-04-16 17:58:54下载
- 积分:1
-
Linear Algebra Package
说明: 特征值分解,HouseHolder分解等(matlab Linear Algebra Package)
- 2019-11-06 22:20:41下载
- 积分:1
-
multiplelinearregression
说明: 多元线性回归的C++源码,在VC6.0上可直接运行!(Multiple Linear Regression C++ source code, can be directly run in VC6.0!)
- 2009-08-25 11:35:52下载
- 积分:1
-
mul
实现有限域中乘法,输入二个普通二级制数,输出在本原多项式的乘法结果(Achieve limited multiplication field, enter the number of two-tier system of two ordinary output in primitive polynomial multiplication results)
- 2014-01-12 22:52:38下载
- 积分:1