沿海近地层湍流特征分析

通过对杭州湾南北两侧共四个观测点近地层湍流特征分析发现:①湍流强度随风速增大而减小。低风速时,几乎与风速呈反比,高风速时基本保持为常数。可描述为:U≤U_0时,I_k=mk/U;U>U_0时,I_k=n_K。其中n_k=a_k/ln(z/z_0),mk=a_kU_0/ln(z/z_0)。a_x=1.0,ay=0.8,a_z=0.5。海岸附近,U_0=4m/s,深入内陆1.5千米以后,U_0=3m/s。②u_*/U=k/ψ(z/L)的 Businger 模式在海岸附近与观测结果吻合很好,在内陆,-1.0≤z/L≤0.3时二者吻合较好,在此范围以外,偏离较大。③不稳定条件下,水平风速标准差可用(?)描述,垂直风速标准差可用(?)描述。稳定条件下,可用(?)(?)描述。此外,稳定条件下,有σ_w/σ_u=0.48,σ_w/σ_v=0.59;不稳定条件下有σ_w/σ_u=0.53,σ_w/σ_v=O.67。     

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

本文主要结果为下述定理。 定理:设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.     

非齐次Poisson过程的若干新性质

本文在一些文献的基础上,进一步讨论非齐次Poisson过程的某些性质,给出了若干新性质。 设{N(t),t≥0}是累积强度为的非齐次Poisson过程。迄今的文献(如[1][2]等)指出,(?)n≥1,前n个到达时刻τ_1,τ_2,…,τ_n的联合概率密度为 本文定理1指出,上式不仅是必要的,而且是充分的,并给出了充分性的证明。从而,得到了描写过程统计规律的一个刻画。然后,在N(t_i)=n_i(i=1,…,k,0相似文献   

三维平壁进气道的特性研究

介绍三维进气道和超声速燃烧冲压式喷气发动机 (SCRAMJET)模型一项实验研究的结果。Scramjet模型由进气道和燃烧室组成 ,此模型组件的设计是为了研究流动结构以获得进气道的特性 ,以及研究燃烧室和进气道、燃料 (氢及碳氢燃料 )点火和燃烧的相互影响。这些试验是在下吹式风洞(M =2 ～ 6,Reunit=( 8～ 54) × 1 0 6)和热射式风洞 (M =6及 7 2 ,Reunit =( 1 0 ～ 32 )× 1 0 6,Tt =1 50 0 ～ 2 50 0K)中进行的。     

常用的振动公式

正弦、位移、速度、频率 如果知道泣侈(D)’l{J峰一峰或涪隔胜,频率川赫竺,(11寸/秒)一」’川: v=汀fD 例:i艾D二一叫·,f二二一0112,则v=3 1 .11!、」/秒。方程l变换,如:则计算0一峰矢徽或峰速伙v(1)也按止个变以,},已知的_几个以D二v/汀fla)f二 V汀D比如D用毫米,t 11。) 如果用国际单位制或米制,也一样。前例中V=0.8米/秒.把结果除l,000,v就是米l秒。 正弦位移、加速度、频率 设D用时,频率用赫芝,现在要计算。一峰矢墩或峰加速度,加速度要)lj重力单位,则计算公式为: A=0 .05llf2D(2) 例:设D=1,f二10Hz,则人=5.泣19。方程2也可变成其…     

导电薄圆筒上的电荷分布——矩量法的改进处理

设薄圆筒以X轴为轴,二端面各为x=-1,x=1。当静电平衡时面电荷密度σ(x)满足积分方程: integral from -1 to 1 G(x∣x′)σ(x′)dx′=常数 (1)设:U(x):=integral from 0 to x σ(x)dx (2)并令:g(U,U′)=G(x∣x′) (3)(1)可表为:ingegral g(U,U′)dU′=常数 (4)对于二维(即圆筒半径为无穷大)情形,(4)的解为 U=2/πsin~(-1)x (5)现以此作为一般情形的尝试解: (ⅰ)把这U区间(-1,1)作2n等分,在与U=-1,-1 2/n,……,1-2/n,1相应的n 1圆环上分布线电荷,其密度各为q_1,q_1 q_2,……,q_(n-1) q_n,q_n,使它们在与U=-1 1/n,……,1-1/n相应的n个圆环上产生相同的电位,对应于(4)可列出n阶线性方程组。 (ⅱ)解出q_1,q_2,…q_n。对于二维情形,可证: q_1=q_2=……q_(n-1) (6)对于一般情形并不如此,但可由此构成新的x-U曲线。 反复(ⅰ),(ⅱ),直到(6)近似满足而使x-U曲线稳定为止。 本法对粗圆筒特别适用,沿圆筒长度不取等分点,而是电荷越密集,取点越密,因而节省计算量,但仍提高了精密度.     

对冲击谱定义的一个讨论

国际标准组织在冲击、振动部分专业术语(ISO 2041)中有一条关于位移、速度、加速度冲击谱的解释(条款3.29),我们认为有不妥之处,在此提出,希与大家讨论。该条内容为:位移、速度、加速度冲击响应谱分别定义为S_d=X≈V/ω≈A/ω~2 (1)S_v=ωX≈V≈A/ω (2)S_a=ω~2X≈ωV≈A (3)共中X、V、A分别代表一组单自由度系统对给定的冲击激励的最大(相对)位移响应、最大相对速度响应、最大(绝对)加速度响应。上述定义可用语言解释为:位移冲击谱是最大相对位移响应谱,它近似等于相对速度     

关于Davidson—Lanczos方法的收敛率

本文对文[1]提出的求解大型对称矩阵A的极端(几个最大或最小)特征值及相应特征向量的Davidson—Lanczos方法,用Rayleigh—Ritz逼近理论,研究了该方法的收敛率。证明了由该方法产生的规范正交向量{v_i}_i~m=1是Krylov子空间K_m≡Span(v_1,Av_1,…,A~(m-1)v_1)的一组基。设A的k个最大特征值为又,λ_1>λ_2>…>λ_k,相应的近似特征值为λ_i~(m)(i=1,…,k),得到 这里γ_i(γ_i>1),W_i和W_i~(m)是常数。     

对某些奇异积分的逐次分半加速法

对于f(x)∈C~(2p 1)[0,1],尤拉-麦克洛林求和公式可写为:(h_k=1/2~k) (?)(1)对此,在逐次减半加速法中计算序列: T_(mi)=(γ_mT_(m-1,i 1)-T_(m-1),:)/(γ_m-1)(m=1,2,3,……) (2)其中γ_m=4~m,形成γ序列(4,4~2,4~3,……)。 对于f(x)∈C~(2p 1)(0,1],我们建立公式: (?)(3)其中(?) (4)对于f(x)=x~βln~nx(β>-1,n=0,1,2,…)则有(?) (5)(3),(5)右边的前半部是由于f(x)在x=0上的奇异性而引起的误差。如将α=2~(β 1)以(n 1)次加入γ序列,就可怕除此项误差。选出十个计算实例作为说明。     

关于Grothendieck代数的特征标

本文用群表示理论和复向量理论证明了Grothendieck代数的一些性质,特别是关于Grothendieck代数基的描述,即设B为K_G(X)的基,Θ={(θ,ρ)|θ是G在X上的轨道,ρ是G_x的一个不可约表示,x∈θ},其中G_x={g∈G|gx=x},则存在1—1对应f:B→Θ使得对于任意V∈B,f(v)=(θ,ρ):V_y≠0<==>y∈θ,且ρ是G_x在V_x上的一个不可约表示,x∈θ,利用这一性质,本文给出了一个求Grothendieck代数的特征标的方法,从而改进了由Luszrig,G.在文[1]中提出的方法,并且给出二面体群D_n关于其一些子群H的Grotheodieck代数的特征标表。     

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ω…     

CARDC 2.4m引射式跨声速风洞设计与运行调试

中国空气动力研究与发展中心研制了可更换喷嘴的中压气体引射器 ,利用现有中压气源驱动 ,建成一座增压回流引射式跨声速风洞。试验段截面尺寸 2 .4m×2 .4m ,M =0 .3～ 1 .2。稳定段最高工作压力为 0 .45MPa ,最高模型试验雷诺数Rec=1 5× 1 0 6(M =0 .90 ,C =0 .2 4m) ,稳定吹风时间≥ 1 5s。风洞气动回路上分别配置有多喷管引射器、栅指扩散段、跨声速试验段驻室抽气系统及特殊的主排气系统等装置。采用智能自适应解耦控制技术 ,实现总压和M数独立、快速、精确地控制。该气动布局与部段配置及其功能设计 ,在国内跨声速风洞中均是首次采用。     

用消去法求逆矩阵的标准程序设计

在本文中,我们将给出一个求逆矩阵的方法和它的程序设计。如果 Γ_(nin)…Γ_(lil)AΓ_(lil)P_1…Γ_(njn)P_nQ=E 则我们有 A~(-1)=Γ_(ljl)…I_(njn)Γ_(njn)…Γ_(ljl)EΓ_(ljl)P_1…Γ_(nfn)P_nQΓ_(nin)…Γ_(lil) 此处 detΓ_(if)≠0,detP_m≠0,(m=1,2,…,n). 这样,我们得到了一个求逆阵的消去法。     

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。     

二维电磁衍射场的静场近似

设通过保角变换: ζ=x+jy=Aζ_1+A_1ζ_1~(-1)+A_2ζ_2~(-2)+……使无限长导体柱的正截面外部变成ζ_1平面上的单位园外部。由二维的Helmholtz公式出发,求得当波长远较柱截面尺寸为大的平面电磁波以垂直于柱轴的方向投射时: (1)E_1平行于柱轴,E_1=exp[jb(ycosα-xsinα)],则远区衍射场 (2)H_1平行于柱轴,H_1=exp[jk(ycosα-xsina)],则远区衍射场 其中:S=截面积, u=cosθ+jsinθ=(x+jy/γ), p=2π(∈μ)(1/2)A(e~(jα)A_1-e~(jα)A),p=P_α+jP_y所相应的矢量P=i_xP_x+i_yP_y就是导体柱在入射波的电场下所感应的等效电矩。 在椭柱(长短半径各为a,b)的情形中: A=(a+b)/2,A_1=(a-b)/2,S=πaba=b就是圆柱的情形;b=0就是薄片的情形,利用Babinet原理,可推得平面上无限长开槽的情形——此二情形都已有准确解,与本文结果相比较,当ka→0时,只差高阶无限小。     

延性金属的断裂强度

由Griffith脆性断裂基础理论引伸,导出了延性断裂理论,求得含有穿透裂纹或表面裂纹非加劲平板结构断裂强度新的表达式。与常用的线弹性断裂力学使用一个材料参数不同,在表达式中使用两个材料参数。本理论独特之处在于两个参数可以由单向拉伸的应力一应变曲线求出;并且,对常用的结构金属,在很宽的裂纹尺寸范围内,应力超过或者低于金属屈服应力下,理论结果和试验数据相当符合。 A—半椭园表面裂纹临界面积,(πac)/2,in~2。(吋~2) Au—在σ=σ_U下半椭园表面裂纹临界面积,in~2。(吋~2) A—埃,0.394×10~(-8)in。(吋) a—半椭园表面裂纹的深度,in。(吋) a_U—在σ=σ_U下半椭园表面裂纹的深度,in。(吋) 2C—穿透裂纹或表面裂纹的长度,in。(吋) 2C_U—在σ=σ_U下穿透裂纹或表面裂纹的长度,in。(吋) 2C_L—在σ=σ_L下穿透裂纹或表面裂纹的长度,in。(吋) E—拉伸时的杨氏模量,Psi(磅/吋~2) h—滑移带的有效高度,in。(吋) h_F—裂纹前缘变形区城的有效高度,in,(吋) h_U—裂纹前缘附近变形区域的有效高度,in。(吋) K_O—线弹性平面应力或混合型的断裂韧性,Psi in~(1/2)。(磅/吋~(3/2)) K_(1C)—线弹性平面应变断裂韧性,Psi in~(1/2)。(磅/吋~(3/2)) K_(TC)—具有中心穿透裂纹的薄板或平板的断裂靱性,Psi(in)~(1/(2 ω)(磅/吋~((3 2ω)/(2 ω)) K_(pC)—具有中心表面裂纹的薄板或平板的断裂靱性,Psi(in.)~(1/(2 ω)(磅/吋~((3 2ω)/(2 ω))) K—厚度参数 L_G—单向拉伸试验中所用的应变片长度,in。(吋) n—ε_(TP)之Ramberg—Osgood关系的指数 P—单位厚度塑性能吸收率,L bs/in。(磅/吋) T—产生单位面积新裂纹表面所消耗的能量,Lbs/in。(磅/吋) t—断裂试件厚度,in。(吋) t—单向拉伸试件厚度,in。(吋) t_o—平面应力断裂的最大厚度,in。(吋) U_E—可用于产生新裂纹表面的单位厚度弹性能,Lbs(磅) U_S—产生新裂纹表面时单位厚度所消耗的能量,Lbs(磅) U_P—塑性变形时单位厚度所消耗的能量,Lbs(磅) U_F—裂纹前缘塑性变形时单位厚度所消耗的能量,Lbs(磅) U_(F1)—在σ=σ_U下,裂纹前缘塑性变形时单位厚度所消耗的能量,Lbs(磅) U_(F2)—在σ=σ_L下,裂纹前缘塑性变形时单位厚度所消耗的能量,Lbs(磅) U_U—裂纹前缘附近塑性变形时单位厚度所消耗的能量,Lbs(磅) U_(U1)—在σ=σ_U下,裂纹前缘附近塑性变形时单位厚度所消耗的能量,Lbs(磅) U_(U2)—在σ=σ_L下,裂纹前缘附近塑性变形时单位厚度所消耗的能量,Lbs(磅) W—试件宽度,in。(吋) W_F—在应力—应变曲线下面,从颈缩开始时的应变到σ_F的应变之间的塑性能密度, Psi(磅/吋~2) W_U—在应力—应变曲线下面,从σ_L的应变到颈缩开始时的应变之同的塑性能密度, Psi(磅/吋~2) β—厚度参数ε_L—在σ=σ_L下的单向拉伸应变ε_N—修正后的颈缩单向拉伸应变ε_U—颈缩开始(σ=0.995σ_U)时的单向拉伸应变ε_F—在σ=σ_F下的修正后的单向拉伸应变ε_F—在σ=σ_F下的平均单向拉伸应变(应变片长度内平均) ε_Y—在σ=σ_Y下的单向拉伸应变ε_(PL)—在σ=σ_L下的单向塑性应变ε_(PU)—在颈缩开始时的应力下的单向塑性应变ε_(PF)—断裂应力下的单向塑性应变ε_(TL)—在σ=σ_L下的单向真正拉伸应变ε_(TY)—在σ=σ_Y下的单向真正拉伸应变ε__(TU)—颈缩开始时的单向真正拉伸应变ε_(TF)—在σ=σ_F下的单向真正拉伸应变ε_(TP)—单向真正塑性拉伸应变ε_(TPU)—在σ=σ_L下的单向真正塑性拉伸应变ε_(TPY)—在σ=σ_Y下的单向真正塑性拉伸应变ε_(TPU)—颈缩开始时的单向真正塑性拉伸应变ε_(TPF)—在σ=σ_F下的单向真正塑性拉伸应变λ—裂纹形状因子μ—厚度参数ν—波松比σ—垂直于裂纹平面的总(毛)面积应力(单向拉伸应力),Psi(磅/吋~2) σ_L—相当于0.0005单向塑性应变的弹性极限拉仲应力,Psi(磅/吋~2) σ_Y—单向屈服拉伸应力,Psi(磅/吋~2) σ_U—单向极限拉伸应力,Psi(磅/吋~2) σ_(UF)—从σ_U至σ_F的平均单向拉伸应力,Psi(磅/吋~2) σ_F—单向断裂拉伸应力,Psi(磅/吋~2) σ_T—单向真正拉伸应力,Psi(磅/吋~2) σ_(TY)—单向真正屈服拉伸应力,Psi(磅/吋~2) σ_(TU)—单向真正极限拉伸应力,Psi(磅/吋~2) σ_(TUF)—从σ_(T_U)至σ(TF)的平均真正单向拉伸应力,Psi(磅/吋~2) σ_(TL)—单向真正极限拉伸应力,Psi(磅/吋~2) σ_(TF)—单向真正断裂拉伸应力,Psi(磅/吋~2) φ—裂纹形状参数ω—断裂靱性参数     

深凹露天矿通风风洞研制与模拟实验

本文扼要阐述了开展露天矿通风研究、尤其是进行实验室风洞物理模拟研究的必要性。对露天矿通风模拟风洞的研制作了简介,指出了本风洞所具有的特点。文中介绍了最近完成的二项风洞模拟实验:对大冶露天矿深部开采风流结构及其污染物扩散净化能力进行了模拟研究,发现该矿具有比较独特的三维风流结构;对大爆破后靠自然风力稀释烟尘蘑菇云团所需的停工通风时间,分别完成了以风速 U_0和炸药量 A 为单因素变量的两组试验,揭示了通风时间与地表风速 U_0、炸药量 A_0的本质联系,推导出了通风时间的计算式。     

COMPUTING REACTIVE FLOW ON A COARSE GRID

INTRODUCTIONConsider the Euler equations supplementedby an additional reactive equation which consistsof a scalar balance law for the mass fraction ofunburntgasρρuρEρY t+ρuρu2 + pρu E + upρu Y x=κ000-ρΘ( T) Ywhereρ is the density,ρu the momentum,ρEthe total specific energy,ρY the mass fraction ofunburnt gas( 0≤ Y≤ 1 ) ,κ a large number( thereaction rate) and T temperature.Θ( T) =01 　if T≤ Tignif T>Tign,Tign is the ignition temperature,E=pγ- 1 + u22 + q0 Y,p is …     

Lyapunov矩阵方程的一种数值解法及部分并行处理

本文设计了求解Lyapunov矩阵方程的一种新方法。所考虑的矩阵方程是 AX—XB=C(1)其中A,B,C分别是m×m,n×n和m×n的已知矩阵。 该方法首先是将系数矩阵A,B初等相似约化为三对角矩阵,即存在可逆矩阵U,V,使U~(-1)AU=A,V~(-1)BV=B,其中A,B为三对角矩阵。然后设计了矩阵方程AY—YB=C的公式解法,分三步: 1)求f(λ)=det(λI—A)的λ各次幂的系数a_0,…,a_m; 2)计算sum from i=1 to m (A_(m-i)-CB~(m-i)),f(B); 3)求解Y。解方程AY—YB=C的方法称为THR算法。 最后经逆变换获得原矩阵方程(1)的解X。 求解矩阵方程(1)的方法称为R—THR算法。该方法的计算量约为m~3+4/3n~3+7m~2n+5nm~2+m~2。 本文给出了R—THR的串行计算的数值例子,并给出了THR算法的并行计算格式。最后通过几种数值方法的比较,表明该方法是可行的,也是有效的。     

与坐标系无关的双圆弧样条

设通过n 2个节点P_0,P_1……P_n,P_(n 1)的拟合曲线是由2(n 1)段光滑连接的圆弧(P_0,Q_0,P_1,Q_1……,P_n,Q_n,P_(n 1)构成,且在边缘节点P_0,P_(n 1)上给定切向。则由判据 integral from P_0 to P_(n 1) k~2ds=min (k:曲率;s:弧长)可以近似得到与解三次样条相似的线性方程组。由此确定各节点上的切向。选定相邻圆弧的其他切点Q_0,Q_1,……Q_n使(?)(i=0,1,……n),因而决定每段圆弧的半径、圆心。这样一组圆弧,称为“双圆弧样条”,与坐标系无关,而且有效地降低了曲率起伏,使整个拟合曲线的曲率变化得均匀而缓和。