使用非正交曲线座标与速度分量S_1流面正问题流场矩阵解

本文基于吴仲华提出的使用非正交曲线座标与相应的非正交速度分量的叶轮机械三元流动气动热力学基本方程组,引入流函数,得出求解的主要方程:流函数的二阶拟线性偏微分方程.除了与密度有关的项以外,流函数的各阶导数都置在方程的左端.这样加快了收敛的速度.用中心九点差分格式,将微分方程离散化后,所得的线性代数方程组用矩阵[L][u]分解直接求解.这种解法收敛速度较快.系数矩阵为对角线带状稀疏矩阵.采用了:(1)非零元素按对角线编号;(2)增设虚点两项办法.大量减少了计算机内存量.由流函数求密度时采用了内存密度函数表插值方法.简单地讨论了松弛因子的选取.用此程序对一些压气机和透平的叶型进行了计算,同实验结果及理论解析解进行了比较,相互是一致的.     

基于混合场积分方程的半空间上方金属目标电磁散射特性高效分析

研究一种可以高效求解半空间金属目标电磁散射积分方程方法,电场积分方程适用于任意结构电磁问题分析,但是生成的矩阵条件数大,迭代求解收敛性差;而磁场积分方程生成的矩阵条件数小,迭代收敛性好,但是仅能分析闭合结构问题,本文采用了混合场积分方程方法,同时具备电场积分方程的普适性与磁场积分方程的收敛性.由于混合场积分方程中涉及格林函数的梯度项,为了进一步加快计算效率,本文引入了一种针对半空间格林函数的高效四维空间插值方法,对组成半空间格林函数的索末菲积分进行列表和Lagrange插值,以实现高效的迭代求解,效率在传统混合场积分方程的基础上提高12.6倍.数值结果表明,该方法在保证精度的同时,可以显著降低求解问题的时间.     

跨声速任意迴转面叶栅流分区计算

本文进一步发展了文献[1]的方法:通道激波前超声区用特征线法,栅前外伸弓形波采用自动伸展;通道波后亚声或跨声区用文献[2]给出的弱守恒流函数方程、引入人工密度、用中心九点差分格式离散主方程、矩阵分解法解出{ψ}场,由文献[2]的办法决定.{ρ}场;迥转面上几个典型算例表明了本方法的工程实用性.文中讨论了τ、Δs、Ω等对流场的影响;捕获了通道波后仍是跨声流时的激波,气流穿过这道波时参数的突跃已明显反映在计算结果中.     

S_1流面跨声流场流函数矩阵解

跨声速叶栅流的计算,可采用时间相关法求解Euler方程,或用松弛方法求解势函数方程和流函数方程。一般说来,时间相关法耗费机时较多,势函数方法仅对无旋流适用。流函数方法适用于二元有旋流的计算,并且边界条件也较为简单,可方便地进行S_1和S_2两类流面迭代得到三元解。流函数方法的跨声计算最大的困难是密度双值问题     

使用非正交曲线坐标与速度分量S_2流面跨音速流函数解

基于叶轮机械两类流面迭代计算理论,在非正交曲线坐标上建立了S_2流面上弱守恒型流函数方程.使用人工密度修正方法求解S_2流面跨音流动正问题,用速度积分方法避免了密度双值问题,并编制了相应的计算机程序.     

用拟流函数求解叶轮机械三维可压流场

本文采用拟流函数方法求解叶轮机械三维可压流场,在跨音情况下,由于存在密度双值问题,采用从第三个动量方程求解密度;针对本文算例的特点,对密度方程提出了适当的边界条件,简化了计算。     

轨道角动量算符的本征值和本征函数与阶梯算符

量子力学的一个重要课题是求解力学量算符的本征值方程,以求得其本征值和本征函数.但是,大部分的本征值方程是二阶交系数微分方程.用通常使用的级数解法求解这类方程比较费时,也不好学.如果采用阶梯算符,可以变立阶变系数微分方程为两个一阶微分方程。这样一来,求解就比较简便,物理意义也很的确,通用性也较强,因此,阶梯算符法是值得推广的. 这篇文章中,我们阐述如何利用阶梯算符法求出轨道角动量平方算符的本征值和本征函数一球谐函数;讨论球谐函数与勒让特多项式及缔合勒让特函数间的关系并求出它们的正规性关系.由于有心力场问题哈密顿角…     

跨声流函数方程的多层网格-强隐式解法

本文是文献[1]工作的继续,是研究在多层网格上采用强隐式过程(即Strogly Implicit Procedure)求解跨声速流函数方程的问题。使用多层网格技术,可消除不同频率或波长的误差分量,有助于残差下降;在每层求解中采用强隐式解法将五对角系数矩阵变为七对角阵以便快速因式分解获取流函数场。文中以四层网格为例,详细研究了在每层求解流函数方程时采用强隐式解和在任意两层间使用FAS法(即Full APProximation Storage)时层间转换等一些重要技术细节,计算了一组典型跨声速双圆弧平面叶栅,所得结果与实验较吻合,并且残差收敛速度要比单层网格快得多,表明多层网格法十分有效。     

用Laplace变换求Hartmann势的精确解及波函数的递推关系

本文求解了在球坐标下Hartmann势的Schrdinger方程,得到了能量方程和归一化的波函数.用Laplace变换使径向的二阶微分方程退化为一阶微分方程,直接积分后用级数展开,应用Laplace逆变换得出本征函数.讨论了径向本征函数的像函数的递推关系,从而得出径向波函数的递推关系.     

非正交座标S_1流面流函数跨音松弛计算

本文在吴仲华教授的叶轮机械三元流动理论的基础上,推导了非正交曲线座标系下的叶轮机械流函数方程及有限差分方程的通用形式.这些方程可用于平面、任意迥转面及任意翘曲的S_1或S_2流面的跨音及亚音流场计算。数值求解中采用了混合差分格式线松弛计算方法。采用了密度预测法由流函数值唯一确定了密度值,解决了流函数方程求解跨音流场的困难,用此方法编制了计算机程序并作了计算,所得结果与实验结果比较一致。     

Separation Transformation and New Exact Solutions of the (N+1)-dimensional Dispersive Double sine-Gordon Equation

In this paper,the separation transformation approach is extended to the(N + 1)-dimensional dispersive double sine-Gordon equation arising in many physical systems such as the spin dynamics in the B phase of 3 He superfluid.This equation is first reduced to a set of partial differential equations and a nonlinear ordinary differential equation.Then the general solutions of the set of partial differential equations are obtained and the nonlinear ordinary differential equation is solved by F-expansion method.Finally,many new exact solutions of the(N + 1)-dimensional dispersive double sine-Gordon equation are constructed explicitly via the separation transformation.For the case of N 2,there is an arbitrary function in the exact solutions,which may reveal more novel nonlinear structures in the high-dimensional dispersive double sine-Gordon equation.     

The MAST FV/FE scheme for the simulation of two-dimensional thermohaline processes in variable-density saturated porous media

A novel methodology for the simulation of 2D thermohaline double diffusive processes, driven by heterogeneous temperature and concentration fields in variable-density saturated porous media, is presented. The stream function is used to describe the flow field and it is defined in terms of mass flux. The partial differential equations governing system is given by the mass conservation equation of the fluid phase written in terms of the mass-based stream function, as well as by the advection–diffusion transport equations of the contaminant concentration and of the heat. The unknown variables are the stream function, the contaminant concentration and the temperature. The governing equations system is solved using a fractional time step procedure, splitting the convective components from the diffusive ones. In the case of existing scalar potential of the flow field, the convective components are solved using a finite volume marching in space and time (MAST) procedure; this solves a sequence of small systems of ordinary differential equations, one for each computational cell, according to the decreasing value of the scalar potential. In the case of variable-density groundwater transport problem, where a scalar potential of the flow field does not exist, a second MAST procedure has to be applied to solve again the ODEs according to the increasing value of a new function, called approximated potential. The diffusive components are solved using a standard Galerkin finite element method. The numerical scheme is validated using literature tests.     

Separation Transformation and New Exact Solutions of the (N+1)-dimensional Dispersive Double sine-Gordon Equation

In this paper, the separation transformation approach is extended to the (N+1)-dimensional dispersive double sine-Gordon equation arising in many physical systems such as the spin dynamics in the B phase of ³He superfluid. This equation is first reduced to a set of partial differential equations and a nonlinear ordinary differential equation. Then the general solutions of the set of partial differential equations are obtained and the nonlinear ordinary differential equation is solved by F-expansion method. Finally, many new exact solutions of the (N+1)-dimensional dispersive double sine-Gordon equation are constructed explicitly via the separation transformation. For the case of N>2, there is an arbitrary function in the exact solutions, which may reveal more novel nonlinear structures in the high-dimensional dispersive double sine-Gordon equation.     

非定常可压等熵流非线性方程显式解析解的推导

本文对作者以前凭试凑、灵感、运气与经验得出的一系列非定常可压流动显式解析解,寻找线索,总结出其可能的推导途径,并以非定常可压等熵一维流为例,具体给出了四种新的求解方法。这些方法会对今后寻找工程热物理领域的非线性主控方程的解析解有所帮助。本文同时还给出了两个新的解析解。     

Representations of Lie algebras and the problem of noncommutative integrability of linear differential equations

An algorithm is proposed for integrating linear partial differential equations with the help of a special set of noncommuting linear differential operators — an analogue of the method of noncommutative integration of finite-dimensional Hamiltonian systems. The algorithm allows one to construct a parametric family of solutions of an equation satisfying the requirement of completeness. The case is considered when the noncommutative set of operators form a Lie algebra. An essential element of the algorithm is the representation of this algebra by linear differential operators in the space of parameters. A connection is indicated of the given method with the method of separation of variables, and also with problems of the theory of representations of Lie algebras. Let us emphasize that on the whole the proposed algorithm differs from the method of separation of variables, in which sets of commuting symmetry operators are used.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 4, pp. 95–100, April, 1991.     

对称分类在非线性偏微分方程组边值问题中的应用

研究了微分方程对称分类在非线性偏微分方程组边值问题中的应用.首先,利用偏微分方程(组)完全对称分类微分特征列集算法确定了给定非线性偏微分方程组边值问题的完全对称分类;其次,利用一个扩充对称将非线性偏微分方程组边值问题约化为常微分方程组初值问题;最后,利用龙格-库塔法求解了常微分方程组初值问题的数值解.     

Extended Determinant Solution of a （3 ＋ 1）-Dimensional KP Equation

In this paper, new extended Grammian determinant solutions to a （3 ＋ 1）-dimensional KP equation are presented by using Hirora＇s bilinear method, and a broad set of suftlcient conditions of systems of linear partial differential equations is given. Moreover, some special solutions of the representative systems are obtained through a systematic analysis.     

Dynamics of combined soliton solutions of unstable nonlinear Schrodinger equation with new version of the trial equation method

In this article, a new version of the trial equation method is suggested. With this method, it is possible to find the new exact solutions of the nonlinear partial differential equations. The developed method is applied to unstable nonlinear Schrödinger equation. New exact solutions are expressed with Jacobi elliptic function solutions, 1-soliton solutions and rational function solutions. When the obtained results are examined, we can say the unstable nonlinear Schrödinger equation shows different dynamic behaviors. In addition, the physical behaviors of these new exact solution are given with two and three dimensional graphs.     

A Novel Low-Dimensional Method for Analytically Solving Partial Differential Equations

This paper is concerned with a low-dimensional dynamical system model for analytically solving partial differential equations (PDEs). The model proposed is based on a posterior optimal truncated weighted residue (POT-WR) method, by which an infinite dimensional PDE is optimally truncated and analytically solved in required condition of accuracy. To end that, a POT-WR condition for PDE under consideration is used as a dynamically optimal control criterion with the solving process. A set of bases needs to be constructed without any reference database in order to establish a space to describe low-dimensional dynamical system that is required. The Lagrangian multiplier is introduced to release the constraints due to the Galerkin projection, and a penalty function is also employed to remove the orthogonal constraints. According to the extreme principle, a set of ordinary differential equations is thus obtained by taking the variational operation of the generalized optimal function. A conjugate gradient algorithm by FORTRAN code is developed to solve the ordinary differential equations. The two examples of one-dimensional heat transfer equation and nonlinear Burgers’ equation show that the analytical results on the method proposed are good agreement with the numerical simulations and analytical solutions in references, and the dominant characteristics of the dynamics are well captured in case of few bases used only.     

New Generalized Transformation Method and Its Application in Higher-Dimensional Soliton Equation

A new generalized transformation method is presented to find more exact solutions of nonlinear partial differential equation. As an application of the method, we choose the (3+1)-dimensional breaking soliton equation to illustrate the method. As a result many types of explicit and exact traveling wave solutions, which contain solitary wave solutions, trigonometric function solutions, Jacobian elliptic function solutions, and rational solutions, are obtained. The new method can be extended to other nonlinear partial differential equations in mathematical physics.