一类滞后相对转动动力学方程的分岔特性及其解析近似解

建立了一类含Davidenkov滞后环的非线性相对转动动力学方程.分别分析了该非线性相对转动自治方程和微外扰下非自治方程的分岔特性,并采用KBM法求解了滞后环指数n=2时该非线性相对转动方程在周期激励下的解析近似解.通过数值仿真,得到了几种分岔结构及外扰下全局分岔图,同时将数值解与本文KBM法求解结果进行比较,证明本文求解结果有较高的精度,为研究这一类滞后相对转动系统提供了理论参考依据. 关键词： 相对转动 滞后环 分岔 KBM法     

基于IB-LBFS和绝对节点坐标法的降落伞柔性结构流固耦合数值模拟

在空气动力学、水动力学和生物流体力学领域中,大变形柔性结构的流固耦合现象是一个重要的非线性力学问题.对该系统的数值模拟是分析这一问题的有效手段.将近年提出的一种Descartes流场求解器,即浸润边界-格子Boltzmann通量求解器(immersed boundary-lattice Boltzmann flux solver,IB-LBFS)作为流场求解方法,并引入绝对节点坐标法(absolute nodal coordinate formulation,ANCF)作为大变形结构分析手段,构建了流固耦合求解器以模拟三维流场中的大变形柔性体运动.使用三维来流中的旗帜摆动算例对该流固耦合求解器进行了验证计算.基于该流固耦合求解器对三维不可压流场中的矩形降落伞和十字形降落伞的展开过程进行了非定常流固耦合数值模拟.     

具有非线性控制的Chua电路的混沌同步

Chua电路是一个非光滑系统.本文通过广义哈密顿系统和观测器方法,将具有非线性控制的Chua电路的混沌同步问题转化成研究具有非线性控制的光滑系统的零解稳定性;进而利用滑模控制对该光滑系统的零解稳定性进行了研究,从而使得Chua电路达到了混沌同步.最后,将上述方法应用到具体系统,数值结果也表明其正确性.     

Birkhoff系统的离散最优控制及其在航天器交会对接中的应用

由于非线性,最优控制问题通常依赖于数值求解,即通过离散目标泛函和受控运动方程转化为一有限维的非线性最优化问题.最优控制问题中的受控运动方程在表示为受控Birkhoff方程的形式之后,可以利用受控Birkhoff方程的离散变分差分格式进行离散.与按照传统差分格式近似受控运动方程相比,此途径可以诱导更加真实可靠的非线性最优化问题,进而也会诱导更加精确有效的离散最优控制.应用于航天器交会对接问题,该种数值求解最优控制问题的方法在较大时间步长的情况下仍然求得了一个有效实现交会对接的离散最优控制.模拟结果验证了该方法的有效性.     

快速数值差分递推算法在偏振模色散作用下光脉冲传输中的应用研究

利用快速数值差分递推算法对非线性薛定谔方程进行推导分析,得到一个同时考虑群速色散效应、非线性效应以及偏振模色散效应的时域快速数值差分递推关系,并将该方法应用于研究偏振模色散作用下的光脉冲传输,得到的计算结果分别与解析结果和传统的分步傅里叶算法所得结果进行对比,表明这种快速数值差分递推算法拥有很高的计算准确度,并且物理模型更加合理,完全适合分析计算偏振模色散作用下的光脉冲传输问题.利用该方法研究了光脉冲的畸变以及脉冲展宽问题,获得的计算结果为光纤传输系统的设计提供了参考.     

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

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

基于遗传算法的非线性映射预测硝基苯化合物的生物毒性

摘要: 非线性映射是一种降维处理技术,通过Sammon算法可以将多维分析问题转化为二维问题分析。应用该法的关键在于寻求二维空间的最佳投影点,这可以转化为一个复杂的非线性优化问题来进行解决。遗传算法是解决各类函数优化和非线性映射问题的一种有效算法。以硝基苯化合物的生物毒性的模式识别为例,编制基于误差函数的目标函数和遗传算法约束参数,应用基于MATLAB的遗传算法工具箱进行优化求解。结果表明:该工具箱在求解此类非线性优化问题上是有效的, 基于遗传算法的非线性映射的预测结果比其它方法要好些。     

基于 MAT LA B的强迫振动达芬 方程的非线性幅频响应分析

利用MAT LA B研究了强迫振动达芬( D u f f i n g) 方程的非线性幅频响应特性, 分析了达芬方程非线性幅 频响应近似解析求解和数值求解的方法和步骤, 给出了相应的 MAT LA B求解程序, 并将解析解与数值解结果进行 了比较. 仿真程序和结果能够加深学生对非线性振动相关知识的理解, 提高大学物理及相关力学课程的课堂教学效 果     

一类非线性微分方程的近似解析解

从理论上研究了一类广义扩散方程的求解问题. 利用相似变换和解析拆分技巧给出了求解该类非线性微分方程近似解的一种有效方法, 方程的解可以表示为一个收敛的幂级数. 近似解结果和数值结果非常符合，证明了所提出的方法的准确性和可靠性, 该方法可以用于解决其他科学和工程技术问题. 关键词： 广义扩散方程 非线性边界值问题 解析拆分 近似解析解     

一类非线性相对转动系统的组合共振分岔与混沌

研究一类具有异宿轨道的非线性相对转动系统的分岔与混沌运动. 应用耗散系统的拉格朗日方程建立一类组合谐波激励作用下非线性相对转动系统的动力学方程. 利用多尺度法求解相对转动系统发生组合共振时满足的分岔响应方程并进行奇异性分析, 得到了系统稳态响应的转迁集. 根据相对转动系统异宿轨道参数方程, 求解了异宿轨道的Melnikov函数, 并给出了系统发生Smale马蹄变换意义下混沌的临界条件. 最后采用数值方法, 通过分岔图, 最大Lyapunov指数图, 相轨迹图和庞加莱截面图研究系统参数对混沌运动的影响.     

An approximation for the boundary optimal control problem of a heat equation defined in a variable domain

In this paper, we consider a numerical approximation for the boundary optimal control problem with the control constraint governed by a heat equation defined in a variable domain. For this variable domain problem, the boundary of the domain is moving and the shape of theboundary is defined by a known time-dependent function. By making use of the Galerkin finite element method, we first project the original optimal control problem into a semi-discrete optimal control problem governed by a system of ordinary differential equations. Then, based on the aforementioned semi-discrete problem, we apply the control parameterization method to obtain an optimal parameter selection problem governed by a lumped parameter system, which can be solved as a nonlinear optimization problem by a Sequential Quadratic Programming （SQP） algorithm. The numerical simulation is given to illustrate the effectiveness of our numerical approximation for the variable domain problem with the finite element method and the control parameterization method.     

基于Monte Carlo的热脉冲法数据分析

热脉冲法是测量聚合物介质薄膜空间电荷分布的有效方法之一, 其数据的分析涉及第一类Fredholm积分方程, 只能采用合适的数值计算方法进行求解, 而Monte Carlo法是近年来提出的数值求解该方程的方法之一. 本文尝试使用Monte Carlo法在频域内实现热脉冲数据的分析, 通过一系列模拟计算讨论Monte Carlo法的分析效果. 计算结果表明: Monte Carlo法可实现对热脉冲法实验数据的有效分析, 提取被测薄膜内的电场分布, 而且计算的电场分布在整个样品厚度上都与真实分布较好地符合, 可有效地弥补尺度变换法只在样品表面附近获得较高准确度的缺陷. 该方法的局限性在于计算结果存在一定的振荡, 且在噪声和数据误差的影响下, 其准确性很大程度上依赖于奇异值分解过程中容差的选择, 在应用的方便程度方面还有待进一步提升.     

Near-field performance analysis of locally-conformal perfectly matched absorbers via Monte Carlo simulations

In the numerical solution of some boundary value problems by the finite element method (FEM), the unbounded domain must be truncated by an artificial absorbing boundary or layer to have a bounded computational domain. The perfectly matched layer (PML) approach is based on the truncation of the computational domain by a reflectionless artificial layer which absorbs outgoing waves regardless of their frequency and angle of incidence. In this paper, we present the near-field numerical performance analysis of our new PML approach, which we call as locally-conformal PML, using Monte Carlo simulations. The locally-conformal PML method is an easily implementable conformal PML implementation, to the problem of mesh truncation in the FEM. The most distinguished feature of the method is its simplicity and flexibility to design conformal PMLs over challenging geometries, especially those with curvature discontinuities, in a straightforward way without using artificial absorbers. The method is based on a special complex coordinate transformation which is ‘locally-defined’ for each point inside the PML region. The method can be implemented in an existing FEM software by just replacing the nodal coordinates inside the PML region by their complex counterparts obtained via complex coordinate transformation. We first introduce the analytical derivation of the locally-conformal PML method for the FEM solution of the two-dimensional scalar Helmholtz equation arising in the mathematical modeling of various steady-state (or, time-harmonic) wave phenomena. Then, we carry out its numerical performance analysis by means of some Monte Carlo simulations which consider both the problem of constructing the two-dimensional Green’s function, and some specific cases of electromagnetic scattering.     

Three-dimensional computer modeling of electromagnetic fields: A global lookback lattice truncation scheme

A new lattice truncation scheme for the finite difference time domain approach to the solution of Maxwell's equations has been developed. The problem space is truncated near the sources and the field components on its boundary are generated from those field values known at retarded times on an interior surface one cell from it with an integral representation of the electromagnetic field. The numerical implementation of this global lookback scheme is discussed. Examples which have been used to determine its characteristics and its validity are given.     

Element-free Galerkin method for a kind of KdV equation

The present paper deals with the numerical solution of the third-order nonlinear KdV equation using the elementfree Galerkin (EFG) method which is based on the moving least-squares approximation.A variational method is used to obtain discrete equations,and the essential boundary conditions are enforced by the penalty method.Compared with numerical methods based on mesh,the EFG method for KdV equations needs only scattered nodes instead of meshing the domain of the problem.It does not require any element connectivity and does not suffer much degradation in accuracy when nodal arrangements are very irregular.The effectiveness of the EFG method for the KdV equation is investigated by two numerical examples in this paper.     

Near-field performance analysis of locally-conformal perfectly matched absorbers via Monte Carlo simulations

In the numerical solution of some boundary value problems by the finite element method (FEM), the unbounded domain must be truncated by an artificial absorbing boundary or layer to have a bounded computational domain. The perfectly matched layer (PML) approach is based on the truncation of the computational domain by a reflectionless artificial layer which absorbs outgoing waves regardless of their frequency and angle of incidence. In this paper, we present the near-field numerical performance analysis of our new PML approach, which we call as locally-conformal PML, using Monte Carlo simulations. The locally-conformal PML method is an easily implementable conformal PML implementation, to the problem of mesh truncation in the FEM. The most distinguished feature of the method is its simplicity and flexibility to design conformal PMLs over challenging geometries, especially those with curvature discontinuities, in a straightforward way without using artificial absorbers. The method is based on a special complex coordinate transformation which is ‘locally-defined’ for each point inside the PML region. The method can be implemented in an existing FEM software by just replacing the nodal coordinates inside the PML region by their complex counterparts obtained via complex coordinate transformation. We first introduce the analytical derivation of the locally-conformal PML method for the FEM solution of the two-dimensional scalar Helmholtz equation arising in the mathematical modeling of various steady-state (or, time-harmonic) wave phenomena. Then, we carry out its numerical performance analysis by means of some Monte Carlo simulations which consider both the problem of constructing the two-dimensional Green’s function, and some specific cases of electromagnetic scattering.     

Iterative regularization methods for atmospheric remote sensing

In this paper we present different inversion algorithms for nonlinear ill-posed problems arising in atmosphere remote sensing. The proposed methods are Landweber's method (LwM), the iteratively regularized Gauss-Newton method, and the conventional and regularizing Levenberg-Marquardt method. In addition, some accelerated LwMs and a technique for smoothing the Levenberg-Marquardt solution are proposed. The numerical performance of the methods is studied by means of simulations. Results are presented for an inverse problem in atmospheric remote sensing, i.e., temperature sounding with an airborne uplooking high-resolution far-infrared spectrometer.     

基于Newmark算法的任意磁化方向铁氧体电磁散射时域有限差分分析

利用坐标系转换矩阵给出实验室坐标系中饱和磁化铁氧体的频域磁导系数张量,再通过频域到时域的转换关系jω→?/?t得到一个二阶微分方程形式的时域本构关系.然后采用Newmark方法求解时域本构关系从而给出一种适用于处理任意磁化方向铁氧体电磁问题的Newmark时域有限差分算法.利用此算法计算了饱和磁化铁氧体层的反(透)射系数和饱和磁化铁氧体球的后向雷达散射截面,所获得的结果验证了此算法的正确有效性.     

Nonlinear waves in media with fifth order dispersion

Nonlinear waves described of the fifth order dispersive nonlinear evolution equation are numerically investigated. The numerical method for boundary value problem for this equation is proposed. Exact solutions to nonlinear evolution equation of the fifth order are given. The numerical method was tested using some exact solutions. The influence of the fifth order dispersion on the propagation of nonlinear waves and formation of the periodic structures is studied.     

非线性波动方程的新数值迭代方法

提出了一种新的求解非线性波动方程的数值迭代法,它是一种半解析的方法.与完全的数值计算方法扰法相比,它能够考虑各阶谐波的相互作用,且能够满足能量守恒定律.用它研究了非线性声波在液体中的传播性质,结果表明,在微扰法适用的声强范围内迭代法也适用,在微扰法不适用的一个较宽的声强范围内迭代法依然适用.