激光输出功率随耦合率及放大率的变化规律

建立了涵盖稳定腔和非稳腔的激光有源理论模型和激光开式谐振腔的边界条件。根据定态薛定谔方程与亥姆霍兹方程的等价性,采用量子力学方法求解谐振腔的亥姆霍兹波动方程,得到满足方程、边界条件、矩阵光学、稳定条件的本征解和本征方程。根据本征方程推导出横模数目随耦合率变化的规律,进而推导出激光输出功率随耦合率变化公式,以及激光输出功率随放大率的变化规律。该理论模型能够从稳定腔自然过渡到非稳腔,在小耦合率情况下退化到与传统公式基本一致的形式,又能在大耦合率情形下与实验结果符合得很好。     

非紧凑低马赫数运动边界散射流动噪声的预测

提出了一种求解非紧凑低马赫数运动边界散射流动诱发噪声的预测方法。该方法首先基于运动坐标系下的连续性方程和动量方程推导得到该坐标系下的波动方程及其积分解,然后在该坐标系下采用边界元方法求解得到非紧凑运动边界表面的声压,最后将求解得到的壁面声压回代到静止或运动坐标系下的积分方程中实现对远场噪声的预测。推导得到的积分方程适用于流体与飞机机身、高速列车车身及旋转叶片等非紧凑结构边界作用诱发噪声的预测。     

无限深方势阱本征值和本征态的三种求解方法

一维无限深方势阱模型是量子力学理想模型,经典教材中势阱的边界一般取得比较特殊.或关于坐标原点具有对称性,或势阱左边界位于坐标原点.本文首先展示了如何利用3种方法求解一维任意边界无限深方势阱能量本征值和对应的本征态,不同方法得到的结果彼此之间等价,讨论分析了这3种方法的推导结果,然后得到关于一维任意边界无限深方势阱能量本征值和本征态的通式,从中比较容易看出这两个物理量均与阱宽有关,并且本征波函数与边界值有关,最后将一维结果拓展到二维和三维任意边界无限深方势阱情况.     

导热问题的有限单元边界积分方程方法

一、引言 边界元计算方法是七十年代迅速发展起来的一种数值计算方法,其主要优点是:将求解区域微分方程的问题转化成求解边界积分方程的问题,因而一般都把物理问题降了一维求解,使该方法计算效率和求解精度都较高.但它用于时关问题和非线性问题时,积分方程中还含有物理量的区域积分项,该方法的优点几乎全部消失.另外,在边界积分方程离散后,代数方程的系数矩阵为满阵.如果边界单元划分很多,其效率不如具     

开口/封闭薄壳体声辐射和散射的统一边界积分方程解法

推导了在二维和三维空间下开口和封闭薄壳体在任意阻抗边界条件下声辐射和散射的统一边界积分方程.相对于以前的求解方法,该方程求解声辐射和散射问题具有相同的影响矩阵,能够同时求解薄壳体气动和振动噪声的辐射和散射现象,以及分析壳体声阻抗对声波传播的影响.推导的方程可以应用于叶轮机械、管道等噪声和消声器消声性能的预测等方面.在此方程基础上,可以进一步考虑运动边界和运动介质对声辐射和散射的影响. 关键词： 薄壳体 声阻抗 积分方程 边界元方法     

抛物方程基于自然边界归化的耦合法

将冯康和余德浩提出的自然边界归化方法^[1~4]应用于求解抛物方程初边值外区域问题,提出一种自然边界元与有限元耦合算法。先将控制方程对时间进行离散化,得到关于时间步长的离散化格式,给出圆外域上的自然积分方程,基于此研究抛物方程无界区域问题的自然边界元与有限元耦合法,最后给出相应的数值例子。     

基于等几何分析方法求解任意截面波导本征问题

基于等几何分析方法具有自由度花费少、高精度、高阶连续性等特点,通过加权余量法对椭圆波导本征问题的亥姆霍兹方程等几何离散得出等几何分析方程.解决了传统数值方法的求解域与几何模型的非一致性问题,实现了将问题的分析计算构架于精确几何模型基础之上.分析任意截面波导的本征问题,对不同偏心率的椭圆波导以及三角形和五边形波导的截止波数的求解结果显示等几何分析方法求解波导本征问题的高效及高精度特性.与传统方法相比,此方法以较少的自由度消耗便会达到较高的求解精度,并且数值解的收敛率较快.     

新的对角形式快速多极边界元法求解声学Helmholtz方程

传统外部声学Helmholtz边界积分方程无法在个人计算机上求解大规模工程问题. 为了有效解决这个问题, 将快速多极方法引入到边界积分方程中, 加速系统矩阵方程组的迭代求解. 由于在边界积分方程中引入基本解的对角形式多极扩展, 新的快速多极边界元法的计算效率与传统边界元相比显著提高, 计算量和存储量减少到O(N)量级(N为问题的自由度数). 包括含有420000个自由度的大型潜艇模型数值算例验证了快速多极边界元法的准确性和高效性, 清楚表明新算法在求解大规模声学问题中的优势,     

利用IEO方法求解多原子分子的振动能谱

对于线性多原子分子体系模型,借助从海森伯方程出发的不变本征算符(IEO)方法,很方便地求解了相应哈密顿量的本征能谱与振动频率,从而给出IEO方法在分子物理中的进一步应用.     

新的对角形式快速多极边界元法求解声学Helmholtz方程

传统外部声学Helmholtz边界积分方程无法在个人计算机上求解大规模工程问题. 为了有效解决这个问题, 将快速多极方法引入到边界积分方程中, 加速系统矩阵方程组的迭代求解. 由于在边界积分方程中引入基本解的对角形式多极扩展, 新的快速多极边界元法的计算效率与传统边界元相比显著提高, 计算量和存储量减少到O(N)量级(N为问题的自由度数). 包括含有420000个自由度的大型潜艇模型数值算例验证了快速多极边界元法的准确性和高效性, 清楚表明新算法在求解大规模声学问题中的优势, 具有良好的工程应用前景.     

Heat Transfe for Power Law Non-Newtonian Fluids

We present a theoretical analysis for heat transfer in power law non-Newtonian fluid by assuming that the thermal diffusivity is a function of temperature gradient. The laminar boundary layer energy equation is considered as an example to illustrate the application. It is shown that the boundary layer energy equation subject to the corresponding boundary conditions can be transformed to a boundary value problem of a nonlinear ordinary differential equation when similarity variables are introduced. Numerical solutions of the similarity energy equation are presented.     

Thomas-Fermi equation with non-spherical boundary conditions

The confined atom Thomas-Fermi equation with non-spherical boundary conditions is considered. A 2-D finite element code for solving the Thomas-Fermi equation with general boundary conditions is demonstrated. Results for both Dirichlet and Neumann boundary conditions for ellipsoids of revolution are presented.     

Boundary element-free method for elastodynamics

The moving least-square approximation is discussed first. Sometimes the method can form an ill-conditioned equation system, and thus the solution cannot be obtained correctly. A Hilbert space is presented on which an orthogonal function system mixed a weight function is defined. Next the improved moving least-square approximation is discussed in detail. The improved method has higher computational efficiency and precision than the old method, and cannot form an ill-conditioned equation system. A boundary element-free method (BEFM) for elastodynamics problems is presented by combining the boundary integral equation method for elastodynamics and the improved moving least-square approximation. The boundary element-free method is a meshless method of boundary integral equation and is a direct numerical method compared with others, in which the basic unknowns are the real solutions of the nodal variables and the boundary conditions can be applied easily. The boundary element-free method has a higher computational efficiency and precision. In addition, the numerical procedure of the boundary element-free method for elastodynamics problems is presented in this paper. Finally, some numerical examples are given.     

Boundary element-free method for elastodynamics

1 Introduction In recent years, more and more attention has been paid to researches on the meshless (or meshfree) method, which makes it a hot direction of computational mechanics[1,2]. The meshless method is the approximation based on nodes, then the large deformation and crack growth problems can be simulated with the method without the re-meshing technique. And the meshless method has some advantages over the traditional computa- tional methods, such as finite element method (FEM) and boun…     

The Yang-Mills Heat Semigroup on Three-Manifolds with Boundary

Long time existence and uniqueness of solutions to the Yang-Mills heat equation is proven over a compact 3-manifold with smooth boundary. The initial data is taken to be a Lie algebra valued connection form in the Sobolev space H ₁. Three kinds of boundary conditions are explored, Dirichlet type, Neumann type and Marini boundary conditions. The last is a nonlinear boundary condition, specified by setting the normal component of the curvature to zero on the boundary. The Yang-Mills heat equation is a weakly parabolic nonlinear equation. We use gauge symmetry breaking to convert it to a parabolic equation and then gauge transform the solution of the parabolic equation back to a solution of the original equation. Apriori estimates are developed by first establishing a gauge invariant version of the Gaffney-Friedrichs inequality. A gauge invariant regularization procedure for solutions is also established. Uniqueness holds upon imposition of boundary conditions on only two of the three components of the connection form because of weak parabolicity. This work is motivated by possible applications to quantum field theory.     

自由界面等离子体平衡的变分原理

在自由界面情况下用求位能汛函极值的变分原理推出了等离子体平衡所需满足的平衡方程和边界条件;其泛函的Euler方程即为具有自由界面条件的磁面函数方程。也给出了便于计算的变分泛函,它相当于等值面边界问题。对于导体壁为简单几何形状的情况,用Rit法进行了数值计算。 关键词：     

Phase Transitions in Edge-Weighted Exponential Random Graphs: Near-Degeneracy and Universality

We consider the boundary value problem of the stationary transport equation in the slab domain of general dimensions. In this paper, we discuss the relation between discontinuity of the incoming boundary data and that of the solution to the stationary transport equation. We introduce two conditions posed on the boundary data so that discontinuity of the boundary data propagates along positive characteristic lines as that of the solution to the stationary transport equation. Our analysis does not depend on the celebrated velocity averaging lemma, which is different from previous works. We also introduce an example in two dimensional case which shows that piecewise continuity of the boundary data is not a sufficient condition for the main result.     

Boundary layer effects on sound in a circular duct

Boundary layer effects on an acoustic field in a unidirectional flow with transverse shear are studied. The acoustic pressure variation in the direction normal to that of the flow is governed in the boundary layer by a second order differential equation. The problem in the boundary layer is reduced from a two point boundary value problem to a one point boundary value problem by transforming the governing equation into the Riccati equation. The Riccati equation is easily integrated with standard numerical procedures. The integration process yields the effective admittance of the wall-boundary layer combination. The acoustic field in the uniform flow is then determined for this effective admittance. Further complications imposed by the boundary layer are thus eliminated. The simplicity of the technique allows calculation of the propagation and decay constants in a circular duct over a wide range of parameters and duct modes.     

Scattering by cavities of arbitrary shape in an infinite plate and associated vibration problems

This paper presents a solution for the displacement of a uniform elastic thin plate with an arbitrary cavity, modelled using the biharmonic plate equation. The problem is formulated as a system of boundary integral equations by factorizing the biharmonic equation, with the unknown boundary values expanded in terms of a Fourier series. At the edge of the cavity we consider free-edge, simply-supported and clamped boundary conditions. Methods to suppress ill-conditioning which occurs at certain frequencies are discussed, and the combined boundary integral equation method is implemented to control this problem. A connection is made between the problem of an infinite plate with an arbitrary cavity and the vibration problem of an arbitrarily shaped finite plate, using the jump discontinuity present in single-layer distributions at the boundary. The first few frequencies and modes of displacement are computed for circular and elliptic cavities, which provide a check on our numerics, and results for the displacement of an infinite plate are given for four specific cavity geometries and various boundary conditions.