Convergence of FFT‐based homogenization for strongly heterogeneous media

The FFT‐based homogenization method of Moulinec–Suquet has recently attracted attention because of its wide range of applicability and short computational time. In this article, we deduce an optimal a priori error estimate for the homogenization method of Moulinec–Suquet, which can be interpreted as a spectral collocation method. Such methods are well‐known to converge for sufficiently smooth coefficients. We extend this result to rough coefficients. More precisely, we prove convergence of the fields involved for Riemann‐integrable coercive coefficients without the need for an a priori regularization. We show that our L² estimates are optimal and extend to mildly nonlinear situations and L^p estimates for p in the vicinity of 2. The results carry over to the case of scalar elliptic and curl ? curl‐type equations, encountered, for instance, in stationary electromagnetism. Copyright © 2014 John Wiley & Sons, Ltd.     

Uniform stabilization to equilibrium of a nonlinear fluid–structure interaction model

We consider uniform stability to a nontrivial equilibrium of a nonlinear fluid–structure interaction (FSI) defined on a two or three dimensional bounded domain. Stabilization is achieved via boundary and/or interior feedback controls implemented on both the fluid and the structure. The interior damping on the fluid combining with the viscosity effect stabilizes the dynamics of fluid. However, this dissipation propagated from the fluid alone is not sufficient to drive uniformly to equilibrium the entire coupled system. Therefore, additional interior damping on the wave component or boundary porous like damping on the interface is considered. A geometric condition on the interface is needed if only boundary damping on the wave is active. The main technical difficulty is the mismatch of regularity of hyperbolic and parabolic component of the coupled system. This is overcome by considering special multipliers constructed from Stokes solvers. The uniform stabilization result obtained in this article is global for the fully coupled FSI model.     

Acceleration of self‐consistent field convergence in ab initio molecular dynamics simulation with multiconfigurational wave function

The Lagrange interpolation of molecular orbital (LIMO) method, which reduces the number of self‐consistent field iterations in ab initio molecular dynamics simulations with the Hartree–Fock method and the Kohn–Sham density functional theories, is extended to the theory of multiconfigurational wave functions. We examine two types of treatments for the active orbitals that are partially occupied. The first treatment, as denoted by LIMO(C), is a simple application of the conventional LIMO method to the union of the inactive core and the active orbitals. The second, as denoted by LIMO(S), separately treats the inactive core and the active orbitals. Numerical tests to compare the two treatments clarify that LIMO(S) is superior to LIMO(C). Further applications of LIMO(S) to various systems demonstrate its effectiveness and robustness. © 2014 Wiley Periodicals, Inc.     

Wavelet multiresolution interpolation Galerkin method for nonlinear boundary value problems with localized steep gradients

The wavelet multiresolution interpolation for continuous functions defined on a finite interval is developed in this study by using a simple alternative of transformation matrix. The wavelet multiresolution interpolation Galerkin method that applies this interpolation to represent the unknown function and nonlinear terms independently is proposed to solve the boundary value problems with the mixed Dirichlet-Robin boundary conditions and various nonlinearities, including transcendental ones, in which the discretization process is as simple as that in solving linear problems, and only common two-term connection coefficients are needed. All matrices are independent of unknown node values and lead to high efficiency in the calculation of the residual and Jacobian matrices needed in Newton’s method, which does not require numerical integration in the resulting nonlinear discrete system. The validity of the proposed method is examined through several nonlinear problems with interior or boundary layers. The results demonstrate that the proposed wavelet method shows excellent accuracy and stability against nonuniform grids, and high resolution of localized steep gradients can be achieved by using local refined multiresolution grids. In addition, Newton’s method converges rapidly in solving the nonlinear discrete system created by the proposed wavelet method, including the initial guess far from real solutions.

     

闭区间[0,1]上的分形插值函数的分数阶积分

介绍了分形插值函数和迭代函数系统以及v阶黎曼-刘准尔分数阶积分的概念及相关定理.利用这些概念及定理讨论了分形插值函数的分数阶积分在[0,1]上连续性及判定[0,1]上的分形插值函数的分数阶积分也是[0,1]上的分形插值函数,并给予了证明.     

一种基于谱线幅值归零的加窗插值FFT谐波检测方法

提出了一种基于谱线幅值归零的加窗插值FFT谐波检测方法用于分析复杂电网信号中的各次谐波.该检测方法将每一次在离散频谱中找到的最大谱线及其左右两侧的数根谱线幅值归零,直到各次谐波分析完毕.通过仿真实验分析了在基波频率稳定时对模拟信号的测量精度,研究了基波频率波动对谐波分析的影响.通过对EAST实验数据的分析,验证了该谐波...     

构造低次有理插值函数的一种方法

关于有理插值的算法已有很多[1,4,5],受二元多项式插值迭加算法[6]的启发,我们给出一种简便的求低次有理插值函数的方法,同时给出有理插值函数存在的充分条件,便于检验.所给方法具有可操作性和实际应用价值,且具有较好的灵活性.     

优化端点条件的平面二次均匀B样条插值曲线

在利用反求法构造B样条插值曲线时,往往需要选取端点条件。 因此,可对端点条件进行优化选取,使得构造的B样条插值曲线满足特定要求。提出了一种利用曲线内能极小选取平面二次均匀B样条插值曲线端点条件的算法。首先给出了二次均匀B样条插值曲线分控制顶点与首个控制顶点（即端点条件）的递推关系式;然后给出了利用曲线内能极小优化选取首个控制顶点的算法,证明了利用该算法构造的C¹连续二次均匀B样条插值曲线为保形插值,并通过数值算例证明了算法的有效性;最后,为便于实际应用,基于MATLAB平台设计了算法所对应的图形用户界面,用户通过简单的操作即可获得光顺的C¹连续二次均匀B样条保形插值曲线。     

有耗介质层上多导体传输线的电磁耦合时域分析方法

目前,针对空间电磁场作用有耗介质层上传输线的电磁耦合,仍缺乏有效的数值分析方法.因此,本文提出一种高效的时域混合算法,很好地解决了有耗介质层上传输线电磁耦合建模难的问题.首先,对经典传输线方程进行改进,推导了适用于有耗介质层上多导体传输线电磁耦合分析的修正传输线方程.然后,结合时域有限差分方法和相应插值技术,求解修正传输线方程,获得多导线及其端接负载上的电压和电流响应,并实现空间电磁场辐射与多导线瞬态响应的同步计算.最后,通过相应计算实例的数值模拟,与CST软件的仿真结果进行对比,验证了时域混合算法的正确性和高效性.     

不可压粘性流体浸入边界直接力法及软件实现

浸入边界法通过在N-S方程中施加体积力模拟不可滑移固壁边界及动边界,避免生成复杂贴体网格及动网格,极大地节省了网格建模时间及动网格计算消耗。本文提出一种新型附加体积力简化计算方法,将简化附加体积力以源项形式嵌入动量方程迭代中,通过用户自定义函数对CFD软件FLUENT二次开发,实现了浸入边界法和通用流体力学求解器的耦合计算。通过静止圆柱和动圆柱绕流数值模拟进行了验证,并探讨了插值函数对计算精度的影响。研究表明,通过引入浸入边界模型,能够提高计算效率,并实现结构网格背景下复杂边界和动边界的高效建模。