A FINITE ELEMENT SOLVER FOR NAVIER-STOKES EQUATIONS VIA VORTICITY AND VELOCITY

INTRODUCTIONThe stationary incompressible Navier- Stokesequations in the primitive variables( velocity uand pressure p)- vΔu + u . u + p =f . u =0u =g　inΩinΩonΓ( 1 )whereΩ is a set of R2 ,Γ its boundary,v thekinematic viscosity of the fluid,f a field of givenbody forces,and g a given field defined onΓ sa-tisfying the global conservation property∫Γg . n =0 ( 2 )n denotes the unit outer normal vector toΓ.The vorticity- velocity system of Eq.( 1 ) ,byintroducing the vorticityω…     

Weibull分布定时截尾寿命试验参数的近似联合置信域

研究了Weibull分布大样本定时截尾寿命试验，利用一种新的途径给出了试验参数的近似联合置信域。设产品寿命X服从Weibull分布，对n个受试样本x1,x2，…，xn进行定时（时间为t0）截尾试验，得到观察数据Si=min(x1,x0)，i=1,2,…，n。利用似然方程，分别得到参数λ，b的极大似然估计λ和b。进一步考察联合分布（λ，b）可以用二元正态分布来近似。再利用多元正态分布与X^2分布的关系，可以推导出未知参数u=（λ，b）的置信度为1－a的联合置信区域（x-u)′I0（x-y)≤xa^2。     

非稳态四阶椭圆方程的Galerkin有限元法

针对平面有界凸多角形区域上一类非稳态四阶椭圆方程 u t+Δ2 u =f (x,t) (x,t)∈Ω× [0 ,T]u = u n=0 (x,t)∈ Ω× [0 ,T]u(x,0 ) =u0 (x) x∈Ω的初边值问题进行了 Galerkin有限元方法分析。首先给出了所讨论问题的 Galerkin有限元方法的离散格式 ,其次对所讨论问题的解与其离散问题的解之间的误差进行了分析研究 ,最后利用椭圆投影算子的性质 ,获得了一定模意义下的一些误差估计     

一类非线性双曲型方程的Galerkin方法

主要讨论了平面有界凸多角形区域上的一类非线性双曲型方程utt- .( a( x,u) u) =f ( x,u)u( x,0 ) =u0 ( x)ut( x,0 ) =φ( x)u( x,t) =0　　( x,t)∈Ω× [0 ,T]x∈Ωx∈Ω( x,t)∈ Ω× [0 ,T]的 Galerkin有限元方法 ,首先给出了所讨论问题的 Galerkin有限元方法的离散格式 ,其次对所讨论问题的解与其离散问题的解之间的误差进行了分析研究 ,最后利用椭圆投影算子的性质 ,得到了 L2 模和能量模方面的一些误差估计。     

多项式系统几何性质的注记

证明：在研究多项式系统的几何性质时，多项式系统应当定义成射影空间中奇点集之余维数至少为2的线场。在这个定义中，多项式系统的次数与通常的次数，在概念上不同相同，通常的n 1次退化系统属于这里的n次系统。值得注意的是，本给出的定义与坐标系的选取无关，在这种定义下，多项式系统的某些几何性质变得非常明显。     

负二项分布可靠性增长的经典分析方法

本文用经典分析方法研究了负二项分布的可靠性增长问題。研制规划由m个阶段组成,每个阶段的失败概率(即不可靠性)满足条件:q_1≥q_2≥…≥q_m,在此条件下,我们得到了qi,i=1,…,m的约束极大似然估计以及最末阶段不可靠性q_m的经典置信上限。我们也讨论了在Bayes分析方法中,如何选取参数的验前分布问题并给出了趋势检验方法。几何分布可靠性增长的相应结果可从负二项分布可靠性增长的结果直接给出。最后用数值例说明了这些方法。     

一类二维样条插值函数的最佳误差估计

本文主要结果为下述定理。 定理:设x(uw)是矩形域上关于该矩形上均匀分割的二维双三次样条插值函数,且x(uw)满足条件(5),则x(uw)在矩形域R边界上的节点处的四阶混合偏导数有估计式: |S_(i,0)|≦|A[i,n—1]||ε_(n,0)| |B[i,n—2]||ε_(0,0)|=[0,-4,(-1)~2 4,…(-1)~i 4]/[0,-4,(-1)~2…(-1)~n 4]|ε_(n,0)| sum from h=i to n-2 (-1)~(k(k-2)-(i 1)(i-2))[0,-4,(-j)~2 4…(-1)~i 4]/[0,-4,(-1)~2 4,…(-1)~(k 1) 4][0,-4,(-1)~2 4,…,(-1)~(k 2)4] (-1)~(i(i 1)/2)/[0,-4,(-1)~2 4,…(-1)~n 4]|ε_(0,0)|其中等号成立的条件分别为: A[i,n—1] B[i,n—2] ε_(n0),ε_(00)>0 A[i,n—1] B[i,n—2] ε_(nm),ε_(0m)>0 其中 i=1,2,…,n—1. j=1,2 …,m—1.     

航空发动机的降低ZP／LTR控制

研究了航空发动机的多变量鲁棒控制，提出了代价设计对象模型的构造方法。该模型是通过对被控对象的发动机执行机构进行特定输出反馈而得到的，它与被控对象模型阶次相同，明显低于目前常用的增广设计方法设计对象模型的阶次。基于低阶设计对象模型，作者以前研究的ZP／LTR设计方法转化为低价ZP／LTR设计方法，应用该方法使设计控制器的阶次明显降低，从而更易实现。通过对某型涡扇发动机的仿真研究验证了低阶ZP／LTR设计方法的有效性。     

先验分布的拟合

本文用经验Bayes方法(简记为EB方法)中的间接法,对二项分布、负二项分布、指数分布的定数及有替换定时截尾等四种情况,进行了产品失效率、可靠性等的先验分布的拟合问题的讨论,给出了可靠性、失效率等的EB点估计及EB区间估计。还用数值说明了这些方法.     

强激波阵面的非平衡结构研究

利用测量强激波波后N2+第一负系(0,0)带和(1,2)带的辐射,对强激波后振动温度历程的测量过程进行了探索,并利用Langmuir探针技术,在低密度激波管中对强激波后电子数密度历程进行了测量.测量和计算结果进行了对比.结果表明:N2+B2∑u+态的激发比振动能的激发更快;实验测得的振动温度有明显的周期性振荡;在激波速度7.65~7.85km/s、p1=1.33Pa、实验段内径0.8m下,实验有效时间只有约6.5μs,实验中的电子数密度不能达到峰值.在约10倍波前自由程的实验有效区域内,电子数密度的测量值与计算值吻合很好.     

基于多点优化法进行可靠性数据折算

为了将非成败型试验数据折算为成败型数据,在分析现有方法局限性的基础上提出多点优化法。该方法于原始置信分布上选取多个折算控制点,并取这些位置点上的误差平方和作为目标函数,进而利用优化算法求解获得折算后的成败型数据和置信分布。该方法没有现有方法的局限性,是一种普遍适用的折算方法。算例表明,该方法较其他现有方法,如二阶矩法,点估计下限法,两点法以及近似优化两点法等,在应用范围、折算精度和计算量方面具有明显的优势,是一种简单而行之有效的折算方法。     

小孔法测量非均匀残余应力的数据处理与误差分析

讨论了逐层钻孔法测量非均匀残余应力的释放应变处理和测量误差对残余应力计算精度的影响。利用转换应力、转换应变,可方便地计算残余应力值。受小孔法本身缺陷的限制,计算所得残余应力误差随孔深增加而加大,测量深度约为应变片平均半径的一半。     

A Real-Valued 2D DOA Estimation Algorithm of Noncircular Signal via Euler Transformation and Rotational Invariance Property

The problem of two-dimensional(2 D)direction of arrival(DOA)estimation for double parallel uniform linear arrays is investigated in this paper.A real-valued DOA estimation algorithm of noncircular(NC)signal is proposed,which combines the Euler transformation and rotational invariance(RI)property between subarrays.In this work,the effective array aperture is doubled by exploiting the noncircularity of signals.The complex arithmetic is converted to real arithmetic via Euler transformation.The main contribution of this work is not only extending the NC-Euler-ESPRIT algorithm from uniform linear array to double parallel uniform linear arrays,but also constructing a new 2 Drotational invariance property between subarrays,which is more complex than that in NCEuler-ESPRIT algorithm.The proposed 2 DNC-Euler-RI algorithm has much lower computational complexity than2 DNC-ESPRIT algorithm.The proposed algorithm has better angle estimation performance than 2 DESPRIT algorithm and 2 D NC-PM algorithm for double parallel uniform linear arrays,and is very close to that of 2 D NC-ESPRIT algorithm.The elevation angles and azimuth angles can be obtained with automatically pairing.The proposed algorithm can estimate up to 2(M-1)sources,which is two times that of 2 D ESPRIT algorithm.Cramer-Rao bound(CRB)of noncircular signal is derived for the proposed algorithm.Computational complexity comparison is also analyzed.Finally,simulation results are presented to illustrate the effectiveness and usefulness of the proposed algorithm.     

多层壁瞬态温度分布计算

风挡玻璃电防冰和除冰装置、飞机机翼及直升机旋翼电除冰装量都是有内加热的多层壁。这些多层壁中温度分布计算不仅在这些装置的设计计算中是必需的，而且在考虑多层壁的热变形及热应力时也是必要的。本文用数值法对有内加热的多层壁温度分布进行了研究，根据导热微分方程及各种边界条件，用克兰克－尼可尔辛法推出了典型结构的有限差表示式，然后用高斯－赛德迭代法求解线性方程组，得到了可用于多层壁瞬态温度分布的计算软件。本文用该软件对实际的多层金属壁及多层玻璃作了计算。计算表明，计算结果与文［３—５］所得结果相符。本文所述的方法可推广到二维及三维瞬态多层壁温度分布计算，也可用于无内加热层的多层壁温度计算。     

高维数值积分的蒙特卡罗方法

给出了对于任意概率密度函数产生随机数的一种方法,同时对随机数进行均匀性及独立性检验,将产生的随机数用于计算高维数值积分的蒙特卡罗平均值方法,得到了一种计算高维数值积分的改进平均值方法,并进行复化.最后,给出了几个数值算例以验证方法的有效性.     

用于非定常粘性流动计算的无网格算法(英文)

把无网格算法发展用于求解涉及动边界的非定常粘性流动问题。在处理粘性流动的布点问题时,通过在远离物面的区域采用无粘流动计算时所用的各向同性点云,而在物面附近引入各向异性点云,从而在保证物面法线方向布点较密以准确模拟边界层的同时,有效控制了布点总量,减少了计算时间。对于动边界问题,本文在上述方法获得的初始布点的基础上,采用基于扰动衰减规律的点云移动技术,快速得到物面运动到其他位置时流场求解所需点云。在所获得的点云上,采用一种带权系数的二次极小曲面逼近方法来离散Navier-Stokes方程的空间导数,权系数的引入使得流场求解时对各向同性点云和各向异性点云无需分开考虑,而是采用统一的求解方法,从而简化了编程求解。用发展的无网格算法,结合Navier-Stokes方程求解双时间推进方法,并耦合Spalart-Allmaras湍流模型,先成功地模拟出NLR7301翼型后缘摆动诱发的非定常流,并通过与实验比较,验证了本方法.接着进行了NACA0012失速模拟,给出了动态失速过程,展示出用本文方法处理复杂动边界问题的效果。     

用MUSIC法确定声源频率和方位时误差的校正方法

为确定多个平面波窄带声源的频率和方位本文在阵列信号处理理论的基础上,提出用均匀线阵接收信号数据,并采用多信号分类法（ＭＵＳＩＣ法）对其进行频率,方位参数的分离估计,利用其正交性的谱估计来判断声源的频率和方位,在此基础地实际阵列中存在的阵元位置扰动,通道特性不一致和阵元之间的耦合这三种模型误差,提出了一种基于辅助源的校正方法,并得模型误差不影响声源的频率,只影响声源方位的判断,声学实验证实,通过校正     

对数正态分布下基于截断数据的统计分析

加速寿命试验是在短时间内获得高可靠产品的失效数据的一种有效办法,因而得到普遍应用。但在低应力下,要在短时间内得到失效数据较为困难,常常进行到一段时间被迫停止,因而得到的是截断数据。本文考虑在有截断数据情形下,对数正态分布未知参数的估计。具体地讲,试验一共有k个加速应力S_1,S_2,…,S_k,有R个产品进行试验,在应力S_i下有n_i个产品。进行加速寿命试验,但是我们得到的不全是寿命失效数据y_j~i,i=1,2,…,k,j=1,2,…,n_i,n_1+n_2+…+n_i=n,而是Z_j~i=min(y_j~i,τ_j~i是预先给定的截断时间,记δ_j~i=I(y_j~i≤τ_j~i),则我们只能观察到(Z_j~i,δ_j~i),i=1,2,…,k,j=1,…,n_i。本文给出了用(Z_j~i,δ_j~i)来估计未知参数的一种方法,证明了这种参数估计的存在性,无偏性及相合性。     

冗余系统两种应力-强度分布模型的可靠度计算

本文运用应力-强度干涉理论,说明如何确定相同设备(或元件)构成的冗余系统的可靠性。讨论了冗余系统的并联和从n中取k两种结构形式的数学模型,推导了应力和强度服从瑞利分布的并联冗余系统,以及应力和强度分别服从麦克思韦分布和瑞利分布n中取k冗余系统的数学公式。最后,给出这两个模型在应力作用下,对应于不同的强度参数,设备或元件的可靠性曲线,直观地说明应力强度参数的改变对可靠性的影响。