二维流场、热结构松耦合模拟研究

本文针对轴对称型飞行器前缘部位,采用松耦合方法,开展了二维高超声速流固界面的热耦合计算。满足非稳态N-S方程的外部流场和内部非稳态热传导均采用商用软件FLUENT进行计算。通过对相同时刻,不同松耦合推进时间步长引起的壁面热流密度、壁面温度、表面传热系数的差异对比,得出了松耦合推进时间对以上参数的影响规律,这对有效实施松耦合方法进行高超声速飞行器模拟具有积极意义。     

气泡与弹性膜的耦合效应数值模拟

本文在前人研究的基础上,计入浮力、表面张力、不同流体密度比等因素对弹性膜附近气泡运动的影响,结合不可压缩理想流体理论,建立气泡与弹性膜耦合动力学数值模型,采用边界元方法进行求解,计算值与Turangan等的实验结果符合良好.通过对弹性膜附近气泡运动的数值模拟,详细分析了弹性膜两侧为同种密度液体以及不同密度液体时气泡的运动,随后又分析气泡在弹性膜和浮力的共同作用下气泡的射流特性.旨在为相关气泡与弹性膜相互作用特性的研究提供参考. 关键词： 弹性膜 边界元法 气泡 浮力     

点间耦合强度对三耦合量子点系统微分电导的影响

本文利用非平衡格林函数运动方程方法,研究了与两个电极耦合在一起的三耦合量子点系统的微分电导及量子干涉的AB振荡问题.通过理论计算发现,由于量子点上的局域态密度的不同从而导致系统电导或隧穿性质的不同,而且量子点间耦合强度、量子点能级等都会对输运性质产生影响. 关键词： 量子点 非平衡格林函数 运动方程 局域态密度     

多弹簧振子耦合系统运动研究

通过理论分析,得到了多弹簧振子耦合系统的解析解,利用Matlab对其进行数值模拟,作出了系统的运动曲线及相轨迹,直观地显示了系统的运动情况,并通过对不同数量弹簧振子耦合系统运动曲线及相轨迹的比较、分析,给出了系统的运动特征.     

双间隙双耦合口输出回路冷测模拟方法

根据微波网络理论,将单耦合口输出回路的冷测模拟物理模型进行扩展,使其适用于双耦合口回路。并在此模型的基础上提出了计算双间隙双耦合口输出回路间隙阻抗频率特性的冷测模拟方法。运用该方法,结合软件仿真,计算分析了某S波段双间隙双耦合口输出腔,分别得到该回路工作在2模式下的阻抗实部与虚部的频率特性,并对计算结果的合理性作了论证。最后通过将计算结果与采用场分析法及散射参数方法得到的结果进行比较和分析,进一步验证该冷测模拟方法的可靠性。     

利用等效电路模型计算耦合腔行波管注-波互作用

通过研究Curnow等效电路模型,得到进行单腔计算的矩阵方程,进而获得线路场方程, 结合运动方程和空间电荷场方程,从而推导出耦合腔行波管一维注波互作用非线性理论模型. 该理论能计算任意的切断、跳变、渐变等结构以及多信号模拟.利用该理论编制的软件计算了 59 GHz-64 GHz耦合腔行波管AM-PM相位失真,三阶互调以及五阶互调分量. 同时模拟得到了59 GHz-64 GHz带内饱和输出功率分布,理论结果与热测结果误差在5%以内.     

棒板耦合结构的振动能量传输研究

本文以棒、板耦合结构为例,从理论和实验上证明和验证了统计能量方法对强耦合结构仍然适用;进而又对该结构的耦合损耗因子进行了较为详细的计算和讨论,并与弱耦合条件下的统计能量方法结果进行了比较,得出了该方法适用的条件以及在低模态密度情况下耦合损耗因子的表示式;最后通过大量实验对所得结论进行了验证,结果是令人满意的。     

二维拉格朗日坐标系下气粒混合双向耦合对激波流场影响的计算

喷射颗粒与气体混合是内爆压缩领域的热点和难点. 针对喷射混合中的气粒双向耦合问题, 开展了理论建模、离散算法以及颗粒反馈对激波流场的影响研究. 建立了拉格朗日计算框架下的数学模型; 给出了耦合源项的离散算法; 开展了平面及汇聚构型条件下, 气粒双向耦合的数值模拟研究; 发现了颗粒反馈导致气体激波提速现象以及气区流场物理量分布形态的改变, 初步获得了量化分析结果. 本文建立的数学模型、计算方法和获得的新的物理认识, 为深入理解喷射混合现象、解决相关工程应用问题提供了重要理论支撑.     

圆波导-矩形波导小孔耦合定向耦合器设计

 设计了圆波导-矩形波导通过小孔耦合的定向耦合器装置,应用于高功率微波器件输出功率测量。利用弱耦合波理论分析了单孔耦合,对单孔耦合进行了理论分析和数值计算,并在HFSS中模拟了中心频率34 GHz、耦合度－40 dB的TE₀₁TE₁₀单孔耦合,将数值计算结果和模拟结果进行了对比,两者最大差值为0.4 dB, 相对误差近1%,从而证实理论分析和数值计算的正确性,证实了此种单孔耦合器不具备定向耦合的条件。在单孔耦合的基础上利用相位叠加原理详细分析了多孔定向耦合阵列原理,得出多孔耦合器且满足反向相位相消的条件才具备定向耦合;为设计此种定向耦合器提供了设计理论和参考依据。在实际加工中只要给出标定的耦合度和方向性系数就可以确定耦合孔尺寸。     

紧耦合方法研究离子原子碰撞的单俘获过程

介绍了双中心原子轨道紧耦合方法计算离子原子碰撞的单俘获截面的基本原理，对于双中心 原子轨道紧耦合方法，找到一种系统化计算模型参数的方法，并根据该方法得到有关模型参 数，计算了O⁶⁺离子与H碰撞系统，所得计算结果在实验误差范围内与实验值很 好符合． 关键词： 双中心原子轨道紧耦合 单俘获截面 模型势 半经典碰撞参数模型 Numer ov算法     

用SPH-MD耦合算法模拟微通道中的液氩流动

本文构造了SPH(光滑粒子动力学)-MD(分子动力学)耦合程序,在不同宽度微通道内,对不同驱动力下的液氩流动进行了模拟计算。结果表明,在通道宽度较窄时,流体速度分布和黏性系数分布均表现出很强的微尺度效应,随着通道宽度增大,这种微尺度效应逐渐减弱,流动特征逐渐接近于宏观Poiseuille流。     

Origin of biquadratic exchange in Fe/Si/Fe

The thickness and temperature dependences of the interlayer exchange coupling in well-defined molecular beam epitaxy-grown Fe/Si/Fe sandwich structures have been studied. The biquadratic coupling shows a strong temperature dependence in contrast to the bilinear coupling. Both depend exponentially on thickness. These observations can be well understood in the framework of Slonczewski's loose spins model [J. Appl. Phys. 73, 5957 (1993)]. No bilinear contribution of the loose spins to the coupling was observed.     

考虑双向耦合效应的格子Boltzmann气固两相湍流模型

在流体粒子概率密度函数输运方程中考虑颗粒对流体的反作用力,发展了考虑双向耦合效应的LB气固两相流模型,引入Smagorinsky亚格子模型模拟高雷诺数气相流场.对经典后台阶气固两相流动进行模拟,气相和颗粒相速度分布与实验结果进行比较,发现考虑双向耦合效应的LB气固两相流模型结果明显优于单向耦合结果.进一步研究不同惯性颗粒在流场中的弥散特性,小颗粒(St~O(0.1))对流体的跟随性较好,在流场中分布较为均匀;而St~O(1)的颗粒难被流场涡卷吸进入涡内,呈现倾向性弥散现象;大颗粒(St~O(10))由于自身惯性进入流场涡,在流场中分布较为均匀.     

Strong Exciton-Plasmon Coupling and Hybridization of Organic-Inorganic Exciton-Polaritons in Plasmonic Nanocavity

We investigate strong exciton-plasmon coupling and plasmon-mediated hybridization between the Frenkel(F)and Wannier-Mott(WM)excitons of an organic-inorganic hybrid system consisting of a silver ring separated from a monolayer WS_2 by J-aggregates.The extinction spectra of the hybrid system calculated by employing the coupled oscillator model are consistent with the results simulated by the finite-difference time-domain method.The calculation results show that strong couplings among F excitons,WM excitons,and localized surface plasmon resonances(LSPRs) lead to the appearance of three plexciton branches in the extinction spectra.The weighting efficiencies of the F exciton,WM exciton and LSPR modes in three plexciton branches are used to analyze the exciton-polaritons in the system.Furthermore,the strong coupling between two different excitons and LSPRs is manipulated by tuning F or WM exciton resonances.     

XFEM modeling of hydraulic fracture in porous rocks with natural fractures

Hydraulic fracture (HF) in porous rocks is a complex multi-physics coupling process which involves fluid flow, diffusion and solid deformation. In this paper, the extended finite element method (XFEM) coupling with Biot theory is developed to study the HF in permeable rocks with natural fractures (NFs). In the recent XFEM based computational HF models, the fluid flow in fractures and interstitials of the porous media are mostly solved separately, which brings difficulties in dealing with complex fracture morphology. In our new model the fluid flow is solved in a unified framework by considering the fractures as a kind of special porous media and introducing Poiseuille-type flow inside them instead of Darcy-type flow. The most advantage is that it is very convenient to deal with fluid flow inside the complex fracture network, which is important in shale gas extraction. The weak formulation for the new coupled model is derived based on virtual work principle, which includes the XFEM formulation for multiple fractures and fractures intersection in porous media and finite element formulation for the unified fluid flow. Then the plane strain Kristianovic-Geertsma-de Klerk (KGD) model and the fluid flow inside the fracture network are simulated to validate the accuracy and applicability of this method. The numerical results show that large injection rate, low rock permeability and isotropic in-situ stresses tend to lead to a more uniform and productive fracture network.     

极性晶体中表面极化子的温度效应

有不少的极性晶体,电子与体纵光学声子的耦合弱,但与表面光学声予的耦合强.本文讨论电子和体纵光学声子耦合弱,与表面光学声子耦合强时对表面极化子的温度特性的影响,用线性组合算符法研究表面极化子的振动频率、诱生势和有效质量的温度依赖性.对AgBr晶体进行了数值计算,结果表明极化子的振动频率,诱生势和有效质量随温度的升高而减小.     

超声波振动辅助下填充填缝工艺的数值研究

本文基于数值方法研究了超声波振动辅助倒装芯片成型下填充填缝工艺。采用了一种基于紧耦合的流固耦合算法,对结构和流场之间的耦合运动进行了数值模拟,其中流场部分采用有限体积法对任意拉格朗日-欧拉方法的不可压缩N-S方程进行离散计算,而结构部分采用有限元法对拉格朗日坐标下的弹性动力学方程进行离散计算。在一个时间步内,结构与流场的计算区域的交界面上进行了多次的数值传递和插值,以保证满足耦合面边界条件。计算了超声波振动作为边界条件作用于结构场,而引起的流场中流体流速的变化规律,并通过实验验证了本方法的正确性。同时,监测了流场入口体积流量的变化规律,探究了超声波振幅、频率以及流体粘度等因素对流体流速的影响,为超声波振动辅助倒装芯片成型下填充工艺的应用提供了理论依据。     

Three-dimensional multi-relaxation-time lattice Boltzmann front-tracking method for two-phase flow

We developed a three-dimensional multi-relaxation-time lattice Boltzmann method for incompressible and immiscible two-phase flow by coupling with a front-tracking technique. The flow field was simulated by using an Eulerian grid, an adaptive unstructured triangular Lagrangian grid was applied to track explicitly the motion of the two-fluid interface, and an indicator function was introduced to update accurately the fluid properties. The surface tension was computed directly on a triangular Lagrangian grid, and then the surface tension was distributed to the background Eulerian grid. Three benchmarks of two-phase flow, including the Laplace law for a stationary drop, the oscillation of a three-dimensional ellipsoidal drop,and the drop deformation in a shear flow, were simulated to validate the present model.     

变声速弹性沉积层下压缩波与剪切波的耦合影响

针对压缩波声速n²线性分布、剪切波声速平方线性分布的变参数弹性海底沉积层中声场建模问题,建立了关于压缩波与剪切波本征函数的偏微分耦合方程组,基于微扰法,通过逐次微分给出了本征函数的近似解析解.理论分析表明,只有剪切波声速平方的梯度不为零的情况下,沉积层中的压缩波才会与剪切波产生耦合.压缩波与剪切波耦合的影响会带来本征值的变化,其对近场声传播的影响较小,但远距离声传播损失的预报在考虑耦合后与实测数据符合得更好.此外,耦合还会导致沉积层中本征函数及其导数的变化,使得耦合与非耦合情况下沉积层中质点位移场不同.通过与商用声场软件COMSOL计算结果的对比分析检验了本文方法的正确性,且本文方法计算时间远小于COMSOL软件运行时间.     

A Consistent Grid Coupling Method for Lattice-Boltzmann Schemes

A method of coupling grids of different mesh size is developed for classical Lattice-Boltzmann (LB) algorithms on uniform grids. The approach is based on an asymptotic analysis revealing suitable quantities equalized along the grid interfaces for exchanging information between the subgrids. In contrast to other couplings the method works without overlap zones. Moreover the grid velocity (Mach number) is not kept constant, as the time step depends not linearly but quadratically on the grid spacing. To illustrate the basic idea we use a simple LB algorithm solving the advection-diffusion equation. The proposed grid coupling is validated by numerical convergence studies indicating, that the coupling does not affect the second-order convergence behavior of the LB algorithm which is observed on uniform grids. In order to demonstrate its principal applicability to other LB models, the coupling is generalized to the standard D2P9 model for (Navier-)Stokes flow and tested numerically. As we use analytic tools different from the Chapman-Enskog expansion, the theoretical background material is given in two appendices. In particular, the results of numerical experiments are confirmed with a consistency analysis.