微平行管道内Jeffrey流体的非定常电渗流动

研究了微平行管道内线性黏弹性流体的非定常电渗流动, 其中线性黏弹性流体的本构关系是由Jeffrey流体模型来描述的. 利用Laplace变换法, 求解了线性化的Poisson-Boltzmann方程、 非定常的柯西动量方程和Jeffrey流体本构方程, 给出了黏弹性Jeffrey流体电渗速度的解析表达式, 分析了无量纲弛豫时间λ₁和滞后时间λ₂对速度剖面的影响. 发现滞后时间为零时, 弛豫时间越小, 速度剖面图越接近牛顿流体的速度剖面图; 随着弛豫时间和滞后时间的增加, 速度振幅也变得越来越大, 随着时间的增加, 速度逐渐趋于恒定. 关键词： 双电层 微平行管道 Jeffrey流体 非定常电渗流动     

金属氢超导性能的研究

 本文用强耦合超导理论研究了金属氢的一些超导性能，求出了具有超导性的金属氢的同位素效应值，热力学临界磁场，比热和T_c，H_c及Δ对有效声子谱α²F(ω)的泛函导数：dT_c/dα²F(ω)，dH_c/dα²F(ω)和d(2Δ/k_BT_c)/dα²F(ω)等与ω的关系曲线，说明超导性的一些参数与金属氢的电子-声子作用的关系。     

第Ⅰ类多模叠加态|ψ1(3)〉q中广义磁场的高次和压缩

构造了由多模相干态|{Z_j} > _q、多模真空态|{O_j} >_q和多模相干态的相反态|{-Z_j} > _q三者的线性叠加所组成的第Ⅰ类三态叠加多模叠加态光场|ψ₁⁽³⁾ > _q,利用多模压缩态理论,研究了态|ψ₁⁽³⁾ > _q中广义磁场分量的等幂次高次和压缩效应.结果发现:态|ψ₁⁽³⁾ > _q是一种典型的三态叠加多模非经典光场;当各模的初始相位和 满足一定的取值条件、并且态|ψ₁⁽³⁾〉_q中任意两态的态间初始相位差(θ_pq^(R)-θ_0q⁽⁰⁾)、(θ_nq^(R)-θ_0q⁽⁰⁾)和(θ_pq^(R)-θ_nq^(R))等分别在各自的闭区间内连续变化时,则态|ψ₁⁽³⁾ > _q的广义磁场分量(即第一正交相位分量)总可分别呈现出周期性变化的、等幂次的奇数模-偶数次、偶数模-奇数次、偶数模-偶数次或者奇数模-奇数次的高次和压缩效应.     

单口径相干合成系统激光光束的M2因子研究

单口径相干合成系统激光光束的光束质量是一个亟待解决的重要问题.基于二阶矩定义, 文中给出了单口径TEM₀₀, TEM₀₁及TEM₁₀两两相干光束M²因子的解析表达式, 并比较分析了束腰宽度、传输距离、振幅之比,以及源场位置矢量对相干光束M²因子的影响, 得到了诸如源场位置参量d₁<100λ时,各相干光束M²因子恒定,反之, 其随位置参量d₁的增大而增大等一些结论.最后,文章对两TEM₀₀模相干光束M²因子的 部分理论进行了实验验证.     

OH, OCI, HOCI(1A')的从头算与势能曲线

采用Gassian09程序包中的多种方法对OH, OCI, HOCI分子的基态结构进行优化计算, 优选出QCISD/6-311G(2df), B3P86/6-311+G(2df)方法分别对OH(X²), OCI(X²)分子进行计算, 得到平衡核间距R_OH=0.09696 nm, R_OCI=0.1569 nm, 谐振频率ω(OH)=3745.37 cm^-1, ω(OCI)=892.046 cm^-1, 与实验结果非常符合. 用Murrell-Sorbie势能函数对OH和OCI分子的扫描势能点进行拟合, 其扫描点都与四参数Murrell-Sorbie函数拟合曲线符合得很好．优选出QCISD(T)/D95(df, pd)方法对HOCI分子进行计算, 得到基态为X¹A', 键长R_OH =0.0966 nm, 键角∠HOCI=102.3°, 谐振频率ω₁(a₁)=738.69 cm^-1, ω₂(b₂)=1260.25 cm^-1, 离解能D_e=2.24eV. 通过比较发现这些结果与实验值符合得很好,并优于文献报道的结果. 随后计算出了力常数, 在此基础上,推导出HOCI分子的多体展式势能函数.报道了HOCI分子对称伸缩振动势能图中在H+OCI →HOCI反应通道上有一鞍点, H原子需要越过1.74eV的能垒才能生成HOCI的稳定结构, 在Cl+OH→HOCI通道上不存在明显势垒, 容易形成稳定的HOCI分子.     

微扩张管道内幂律流体非定常电渗流动

研究了电渗驱动下幂律流体在有限长微扩张管道内非稳态流动特性.基于Ostwald-de Wael幂律模型,采用高精度紧致差分离散二维Poisson-Nernst-Planck方程及修正的Cauchy动量方程,数值模拟了初始及稳态时刻微扩张管道内幂律流体电渗流流场分布情况,研究了管道截面改变对幂律流体无量纲剪切应变率及无量纲表观黏度的影响,以及无量纲表观黏度对拟塑性流体与胀流型流体流速分布的影响.数值模拟结果显示,当扩张角和无量纲电动宽度一定时,电场驱动下的幂律流体在近壁区域速度响应都很快;初始时刻,近壁处表观黏度的变化受到剪切应变率变化的影响,从而影响了三种幂律流体速度峰值的分布,出现拟塑性流体流速在扩张段上游及扩张段近壁处速度峰值均为幂律流体中最大、而在扩张段下游三种幂律流体速度峰值相近的现象;稳态时刻,幂律流体速度剖面呈现塞型分布,且满足连续性条件下,幂律流体流速随扩张管半径增大而减小,牛顿流体流动规律与宏观尺度下流动规律相同;初始时刻,在相同电动宽度、不同壁面电势作用下,幂律流体在扩张管近壁处剪切应变率分布的差异导致表观黏度分布的差异,并最终导致拟塑性流体与胀流型流体流速分布的差异.     

光受激发射的稳定性

本文推导出描述三能级Laser工作过程的准经典方程组,并分析了输出振动的稳定性。在阈值以上,当T₁?T₂,q^-1时,只在1/(qT₂)>1时,输出振幅是稳定的(其中T₁,T₂,q^-1分别是分子纵向、横向及谐振腔的弛豫时间)。在稳定区域,趋向平衡的时间与T₁成正比。当分子线宽小于谐振腔宽度时,输出是不稳定的,而在1/(qT₂)减小时,平衡点由稳定变到不稳定时产生一个稳定的极限环,即输出振幅逐渐开始振动。关于稳定性的结论在气体Laser中是可以检验的。本文指出,在红宝石Laser中看到的输出不稳定,可能就是谐振腔的q很大的结果。     

微通道中一类生物流体在高Zeta势下的电渗流及传热特性

在高壁面Zeta势下,研究滑移边界条件下满足牛顿流体模型的一类生物流体的电渗流动及传热特性,流体在外加电场、磁场和焦耳加热共同作用下流动.首先,在不使用Debye-Hückel线性近似条件时,利用切比雪夫谱方法给出非线性Poisson-Boltzmann方程和流函数满足的四阶微分方程及热能方程的数值解,将所得结果与利用Debye-Hückel线性近似所得结果进行比较,证明本文数值方法的有效性.其次,讨论电磁环境下壁面Zeta势、哈特曼数H、电渗参数m、滑移参数β对流动特性、泵送特性和捕获现象的影响,并探究焦耳加热参数γ和布林克曼数Br等参数对传热特性的影响.结果表明,壁面Zeta势、电渗参数m、滑移参数β的增大对流体速度有促进作用,而哈特曼数H的增大会抵抗流体流动.研究进一步表明,焦耳加热参数γ和布林克曼数Br的增大会导致温度升高.     

SiN自由基X2∑+, A2Π和B2∑+ 电子态的光谱常数研究

采用内收缩多参考组态相互作用(MRCI)方法, 结合价态范围内的最大相关一致基aug-cc-pV6Z, 计算了SiN自由基X²∑⁺, A²Π 和B²∑⁺电子态的势能曲线. 采用Davidson修正来避免由于MRCI方法本身的大小一致性缺陷产生的误差. 为了提高计算精度, 进一步考虑了相对论修正和核价相关修正对势能曲线的影响. 本文的相对论修正是利用二阶Douglas-Kroll 哈密顿近似在cc-pV5Z基组水平进行的; 同时核价相关修正是在cc-pCV5Z基组水平进行的. 对这些势能曲线进行拟合, 得到各种水平下三个电子态的光谱常数(T_e, R_e, ω_e, ω_ex_e, α_e和B_e), 并详细分析了Davidson修正、相对论修正和核价相关修正对光谱常数的影响. 与其他理论结果和实验数据进行比较, 可知本文的结果更精确、更完整.     

AsO+同位素离子X2∑+和A2∏电子态的多参考组态相互作用方法研究

采用包含Davidson修正多参考组态相互作用(MRCI)方法结合价态范围内的最大相关一致基As/aug-cc-pV5Z和O/aug-cc-pV6Z, 计算了AsO⁺ (X²∑⁺)和AsO⁺(A²∏)的势能曲线. 利用AsO⁺离子的势能曲线在同位素质量修正的基础上, 拟合出了同位素离子⁷⁵As¹⁶O⁺和⁷⁵As¹⁸O⁺的两个电子态光谱常数. 对于X²∑⁺态的主要同位素离子⁷⁵As¹⁶O⁺, 其光谱常数R_e, ω_e, ω_ex_e, B_e和α_e分别为 0.15770 nm, 1091.07 cm^-1, 5.02017 cm^-1, 0.514826 cm^-1和0.003123 cm^-1; 对于A²∏态的主要同位素离子⁷⁵As¹⁶O⁺, 其T_e, R_e, ω_e, ω_ex_e, B_e和α_e分别为5.248 eV, 0.16982 nm, 776.848 cm^-1, 6.71941 cm^-1, 0.443385 cm^-1和0.003948 cm^-1. 这些数据与已有的实验结果均符合很好. 通过求解核运动的径向薛定谔方程, 找到了J=0时AsO⁺(X²∑⁺)和AsO⁺(A²∏)的前20个振动态. 对于每一振动态, 还分别计算了它的振动能级、转动惯量及离心畸变常数, 并进行了同位素质量修正, 得到各同位素离子的分子常数. 这些结果与已有的实验值非常一致. 本文对于同位素离子⁷⁵As¹⁶O⁺(X¹∑⁺), ⁷⁵As¹⁸O⁺(X¹∑⁺), ⁷⁵As¹⁶O⁺(A¹∏)和⁷⁵As¹⁶O⁺(A¹∏)的光谱常数和分子常数属首次报导.     

Time Periodic Electroosmotic Flow of The Generalized Maxwell Fluids in a Semicircular Microchannel

Analytical solutions are presented using method of separation of variables for the time periodic electroosmotic flow (EOF) of linear viscoelastic fluids in semicircular microchannel. The linear viscoelastic fluids used here are described by the general Maxwell model. The solution involves analytically solving the linearized Poisson-Boltzmann (P -B) equation, together with the Cauchy momentum equation and the general Maxwell constitutive equation. By numerical computations, the influences of electric oscillating Reynolds number Re and Deborah number De on velocity amplitude are presented. For small Re, results show that the larger velocity amplitude is confined to the region near the charged wall when De is small. With the increase of the Deborah number De, the velocity far away the charged wall becomes larger for large Deborah number De. However, for larger Re, the oscillating characteristic of the velocity amplitude occurs and becomes significant with the increase of De, especially for larger Deborah number.     

MHD Maxwell flow modeled by fractional derivatives with chemical reaction and thermal radiation

Heat transfer in a time-dependent flow of incompressible viscoelastic Maxwell fluid induced by a stretching surface has been investigated under the effects of heat radiation and chemical reaction. The magnetic field is applied perpendicular to the direction of flow. Velocity, temperature, and concentration are functions of z and t for the modeled boundary-layer flow problem. To have a hereditary effect, the time-fractional Caputo derivative is incorporated. The pressure gradient is assumed to be zero. The governing equations are non-linear, coupled and Boussinesq approximation is assumed for the formulation of the momentum equation. To solve the derived model numerically, the spatial variables are discretized by employing the finite element method and the Caputo-time derivatives are approximated using finite difference approximations. It reveals that the fractional derivative strengthens the flow field. We also observe that the magnetic field and relaxation time suppress the velocity. The lower Reynolds number enhances the viscosity and thus motion weakens slowly. The velocity initially decreases with increasing unsteadiness parameter δ. Temperature is an increasing function of heat radiation parameter but a decreasing one for the volumetric heat absorption parameter. The increasing value of the chemical reaction parameter decreases concentration. The Prandtl and Schmidt numbers adversely affect the temperature and concentration profiles respectively. The fractional parameter changes completely the velocity profiles. The Maxwell fluids modeled by the fractional differential equations flow faster than the ordinary fluid at small values of the time t but become slower for large values of the time t.     

Metachronal propulsion of non-Newtonian viscoelastic mucus in an axisymmetric tube with ciliated walls

Cilia-induced flow of viscoelastic mucus through an idealized two-dimensional model of the human trachea is presented.The cilia motion is simulated by a metachronal wave pattern which enables the mobilization of highly viscous mucus even at nonzero Reynolds numbers.The viscoelastic mucus is analyzed with the upper convected Maxwell viscoelastic formulation which features a relaxation time and accurately captures normal stress generation in shear flows.The governing equations are transformed from fixed to wave(laboratory)frame with appropriate variables and resulting differential equations are perturbed about wave number.The trachea is treated as an axisymmetric ciliated tube.Radial and axial distributions in axial velocity are calculated via the regular perturbation method and pressure rise is computed with numerical integration using symbolic software MATHEMATICA^‘TM’.The influence of selected parameters which is cilia length,and Maxwell viscoelastic material parameter i.e.relaxation time for prescribed values of wave number are visualized graphically.Pressure rise is observed to increase considerably with elevation in both cilia length and relaxation time whereas the axial velocity is markedly decelerated.The simulations provide some insight into viscous-dominated cilia propulsion of rheological mucus and also serve as a benchmark for more advanced modeling.     

Lattice Boltzmann simulation of transverse wave travelling in Maxwell viscoelastic fluid

A nine-velocity lattice Boltzmann method for Maxwell viscoelastic fluid is proposed. Travelling of transverse wavein Maxwell viscoelastic fluid is simulated. The instantaneous oscillating velocity, transverse shear speed and decay rateagree with theoretical results very well.     

弱电离大气等离子体电子能量分布函数的理论研究

使用球谐展开的方法求解玻尔兹曼方程,得到了弱电离大气等离子体(79％氮气和21％的氧气)的电子能量分布函数(EEDF).发现当约化电场较小时(E/N＜100 Td),EEDF在2-3 eV急剧下降,在此情况下,高能尾部比麦氏分布要小;当约化电场增加,E/N＞ 400 Td,分布函数趋近于麦氏分布;当约化电场进一步增加,E/N＞ 2000 Td,EEDF的高能尾部(超过200 eV)相对于麦氏分布增加,在高频场作用下,EEDF更倾向于麦氏分布.当ω》vm时,有效电子温度只依赖于E/ω,而与碰撞频率无关;当ω《vm时,有效电子温度只依赖于E/N,与微波频率无关.与一些单原子分子等离子体中电子-电子碰撞在电离度大于10-6时就会影响EEDF不同,空气等离子体中,只有当电离度大于0.1％时,电子-电子碰撞才会对EEDF有明显影响.     

ρ-ω mixing in J/ø→VP decays

The study of ρ-ω mixing has mainly focused on vector meson decays with isospin I=1, namely the ρ(ω)→π⁺π^- process. In this paper, we present a study of ρ-ω mixing in ρ(ω)→π⁺π^-π⁰ (I=0) using a flavor parameterization model for the J/ø→ VP process. By fitting a theoretical framework to PDG data, we obtain the SU(3)-breaking effect parameters s_V=0.03±0.12, s_P=0.17±0.17 and the ρ-ω mixing polarization operator ∏_ρω=(0.006±0.011) GeV². New values are found for the branching ratios when the mixing effect is incorporated: Br(J/ψ→ ωπ⁰) = (3.64 ±0.37)×10^-4, Br(J/ψ→ ωη) = (1.48 ±0.17)×10^-3, Br(J/ψ→ ωη') = (1.55±0.56)×10^-4, these are different from the corresponding PDG2012 values by 19%, 15% and 15%, respectively.     

Sloshing waves in a heated viscoelastic fluid layer in an excited rectangular tank

In this paper, we have investigated the motion of a heated viscoelastic fluid layer in a rectangular tank that is subjected to a horizontal periodic oscillation. The mathematical model of the current problem is communicated with the linearized Navier–Stokes equation of the viscoelastic fluid and heat equation together with the boundary conditions that are solved by means of Laplace transform. Time domain solutions are consequently computed by using Durbin's numerical inverse Laplace transform scheme. Various numerical results are provided and thereby illustrated graphically to show the effects of the physical parameters on the free-surface elevation time histories and heat distribution. The numerical applications revealed that increasing the Reynolds number as well as the relaxation time parameter leads to a wider range of variation of the free-surface elevation, especially for the short time history.     

Viscoelasticity and metastability limit in supercooled liquids

A supercooled liquid is said to have a kinetic spinodal if a temperature Tsp exists below which the liquid relaxation time exceeds the crystal nucleation time. We revisit classical nucleation theory taking into account the viscoelastic response of the liquid to the formation of crystal nuclei and find that the kinetic spinodal is strongly influenced by elastic effects. We introduce a dimensionless parameter lambda, which is essentially the ratio between the infinite frequency shear modulus and the enthalpy of fusion of the crystal. In systems where lambda is larger than a critical value lambda(c) the metastability limit is totally suppressed, independently of the surface tension. On the other hand, if lambda相似文献   

Collisional damping of transverse waves in a plasma

The effect of collisions on transverse waves in a homogeneous, field free plasma is investigated by means of Gross-Krook collision model. The dispersion relation is calculated by assuming the collision frequency to be small andKλ _D ?1. The damping rate ω _I is obtained as $$\omega _I = \frac{{\nu _{ei} }}{2}\frac{{\omega _p^2 }}{{\omega _0^2 }}\left[ {1 + \frac{{3K^2 \lambda _D^2 \omega _p^2 }}{{\omega _0^2 }} - \frac{{K^2 \lambda _D^2 \omega _p^4 }}{{\omega _0^4 }}} \right] + \frac{{\nu _{ee} }}{2}\frac{{\omega _p^2 }}{{\omega _0^2 }}\left( {\frac{{K^2 \lambda _D^2 \omega _p^2 }}{{\omega _0^2 }}} \right)$$ where ω ₀ ² =c ² K ²+ω ₂ ^p , andv _ei andv _ee are electron-ion and electron-electron collision frequency respectively.