1、 采用原始变量法,即以速度U、V及压力P作为直接求解的变量 2、 守恒型的差分格式,离散方程系对守恒型的控制方程通过对控制容积作积分而得出的,无论网格疏密程度如何,均满足在计算区域内守恒的条件; 3、 采用区域离散化方法B,即先定控制体界面、再定节点位置 4、 采用交叉网格,速度U、V与其他变量分别存储于三套网格系统中; 5、 不同的项在空间离散化过程中去不同的型线假设,源项采用局部线性化方法;扩散——对流项采用乘方格式(但很容易转化为中心差分、迎风差分或混合格式);街面上的扩散系数采用调和平均法,而密度与流速则用线性插值; 6、 不稳态问题采用全隐格式,以保证在任何时间步长下均可获得具有物理意义的解; 7、 边界条件采用附加源项法处理; 8、 耦合的流速与压力采用SIMPLE算法来求解; 9、 迭代式的求解方法,对非线性问题,整个求解过程具有迭代性质;对于代数方程也采用迭代法求解; 10、 采用交替方向先迭代法求解代数方程并补以块修正技术以促进收敛。
标签: 变量
上传时间: 2016-12-28
上传用户:heart520beat
1、 采用原始变量法,即以速度U、V及压力P作为直接求解的变量 2、 守恒型的差分格式,离散方程系对守恒型的控制方程通过对控制容积作积分而得出的,无论网格疏密程度如何,均满足在计算区域内守恒的条件; 3、 采用区域离散化方法B,即先定控制体界面、再定节点位置 4、 采用交叉网格,速度U、V与其他变量分别存储于三套网格系统中; 5、 不同的项在空间离散化过程中去不同的型线假设,源项采用局部线性化方法;扩散——对流项采用乘方格式(但很容易转化为中心差分、迎风差分或混合格式);街面上的扩散系数采用调和平均法,而密度与流速则用线性插值; 6、 不稳态问题采用全隐格式,以保证在任何时间步长下均可获得具有物理意义的解; 7、 边界条件采用附加源项法处理; 8、 耦合的流速与压力采用SIMPLE算法来求解; 9、 迭代式的求解方法,对非线性问题,整个求解过程具有迭代性质;对于代数方程也采用迭代法求解; 10、 采用交替方向先迭代法求解代数方程并补以块修正技术以促进收敛。
标签: 变量
上传时间: 2013-12-28
上传用户:Avoid98
利用verilog语言,从上至下层次管理的设计思想;Verilog HDL的行为描述和结构描述,实现8位频率计,4个0检测修正电路的原理说明
上传时间: 2016-12-29
上传用户:caozhizhi
基于GM算法和QR分解实现的稳健奇异值分解算法,通过该算法可以获取一个数居阵的稳健方差-协方差阵。在该稳健方差-协方差阵上可以进行诸如,稳健主成份分解、稳健聚类、稳健因子分析等操作。
上传时间: 2013-12-13
上传用户:15071087253
Euler函数: m = p1^r1 * p2^r2 * …… * pn^rn ai >= 1 , 1 <= i <= n Euler函数: 定义:phi(m) 表示小于等于m并且与m互质的正整数的个数。 phi(m) = p1^(r1-1)*(p1-1) * p2^(r2-1)*(p2-1) * …… * pn^(rn-1)*(pn-1) = m*(1 - 1/p1)*(1 - 1/p2)*……*(1 - 1/pn) = p1^(r1-1)*p2^(r2-1)* …… * pn^(rn-1)*phi(p1*p2*……*pn) 定理:若(a , m) = 1 则有 a^phi(m) = 1 (mod m) 即a^phi(m) - 1 整出m 在实际代码中可以用类似素数筛法求出 for (i = 1 i < MAXN i++) phi[i] = i for (i = 2 i < MAXN i++) if (phi[i] == i) { for (j = i j < MAXN j += i) { phi[j] /= i phi[j] *= i - 1 } } 容斥原理:定义phi(p) 为比p小的与p互素的数的个数 设n的素因子有p1, p2, p3, … pk 包含p1, p2…的个数为n/p1, n/p2… 包含p1*p2, p2*p3…的个数为n/(p1*p2)… phi(n) = n - sigm_[i = 1](n/pi) + sigm_[i!=j](n/(pi*pj)) - …… +- n/(p1*p2……pk) = n*(1 - 1/p1)*(1 - 1/p2)*……*(1 - 1/pk)
上传时间: 2014-01-10
上传用户:wkchong
//Euler 函数前n项和 /* phi(n) 为n的Euler原函数 if( (n/p) % i == 0 ) phi(n)=phi(n/p)*i else phi(n)=phi(n/p)*(i-1) 对于约数:divnum 如果i|pr[j] 那么 divnum[i*pr[j]]=divsum[i]/(e[i]+1)*(e[i]+2) //最小素因子次数加1 否则 divnum[i*pr[j]]=divnum[i]*divnum[pr[j]] //满足积性函数条件 对于素因子的幂次 e[i] 如果i|pr[j] e[i*pr[j]]=e[i]+1 //最小素因子次数加1 否则 e[i*pr[j]]=1 //pr[j]为1次 对于本题: 1. 筛素数的时候首先会判断i是否是素数。 根据定义,当 x 是素数时 phi[x] = x-1 因此这里我们可以直接写上 phi[i] = i-1 2. 接着我们会看prime[j]是否是i的约数 如果是,那么根据上述推导,我们有:phi[ i * prime[j] ] = phi[i] * prime[j] 否则 phi[ i * prime[j] ] = phi[i] * (prime[j]-1) (其实这里prime[j]-1就是phi[prime[j]],利用了欧拉函数的积性) 经过以上改良,在筛完素数后,我们就计算出了phi[]的所有值。 我们求出phi[]的前缀和 */
上传时间: 2016-12-31
上传用户:gyq
C51最常用的IDE环境,KEIL C51的最新版本Keil C51 V8.15。包含官方程序安装包和最新的注册机,已经修正了0xfd的bug,测试有效,笔者现在用的就是这个注册机破解版本。
上传时间: 2014-01-11
上传用户:jcljkh
delphi老牌的浏览器控件,现在重新开始升级,这是目前最新的版本。 基于IE,功能十分强大。 个人对此版本做了一些修正。
上传时间: 2014-12-08
上传用户:zl5712176
内点法的c源程序,实现惩罚因子,循环次数
标签: 源程序
上传时间: 2013-12-19
上传用户:moerwang
包含DVB-T标准中的内交织,内编码,卷积交织,和基4FFT 1024。所有程序为MATLAB实现,特别推荐FFT旋转因子求法。
上传时间: 2013-12-22
上传用户:洛木卓