基于Kronecker型替换的黑盒多项式稀疏插值

本文基于新的Kronecker型替换,给出两个由黑盒表示的稀疏多项式的新确定性插值算法.令f∈R[x1,……,xn]是一个稀疏黑盒多项式,其次数上界为D.当R是C或者是有限域时,相对于已有算法,新算法具有更好的计算复杂度或者关于D的复杂度更低.特别地,对于一般黑盒模型,D是复杂度中的主要因素,而在所有的确定性算法中,本文的第二个算法的复杂度关于D是最低的.     

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

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

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

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

Arbitrage‐free call option surface construction using regression splines

In this work, we suggest a novel quadratic programming‐based algorithm to generate an arbitrage‐free call option surface. The empirical performance of the proposed method is evaluated using S&P 500 Index call options. Our results indicate that the proposed method provides a more precise fit to observed option prices than other alternative methodologies. Copyright © 2014 John Wiley & Sons, Ltd.     

A coupled immersed boundary‐lattice Boltzmann method with smoothed point interpolation method for fluid‐structure interaction problems

The immersed boundary‐lattice Boltzmann method has been verified to be an effective tool for fluid‐structure interaction simulation associated with thin and flexible bodies. The newly developed smoothed point interpolation method (S‐PIM) can handle the largely deformable solids owing to its softened model stiffness and insensitivity to mesh distortion. In this work, a novel coupled method has been proposed by combining the immersed boundary‐lattice Boltzmann method with the S‐PIM for fluid‐structure interaction problems with large‐displacement solids. The proposed method preserves the simplicity of the lattice Boltzmann method for fluid solvers, utilizes the S‐PIM to establish the realistic constitutive laws for nonlinear solids, and avoids mesh regeneration based on the frame of the immersed boundary method. Both two‐ and three‐dimensional numerical examples have been carried out to validate the accuracy, convergence, and stability of the proposed method in consideration of comparative results with referenced solutions.     

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

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

第二类端点奇异Fredholm积分方程的分数阶退化核方法

本文针对第二类端点奇异Fredholm积分方程构造基于分数阶Taylor展开的退化核方法，设计了两种计算格式，一是在全区间上使用分数阶Taylor展开式近似核函数，二是在包含奇点的小区间上采用分数阶插值，在剩余区间上采用分段二次多项式插值逼近核函数.讨论了两种退化核方法收敛的条件，并给出了混合插值法的收敛阶估计.数值算例表明对于非光滑核函数分数阶退化核方法有着良好的计算效果，且混合二次插值法比全区间上的分数阶退化核方法有着更广泛的适用范围.     

Design of sextuple-mode triple-ring HTS UWB filter using two-round interpolation

A single-stage ring resonator capable of introducing six modes within the ultra-wideband (UWB) passband is presented. The sextuple-mode resonator consists of three rings and three sets of stepped-impedance open stubs. Based on this sextuple-mode resonator, a UWB filter fed by the interdigital-coupling line (ICL) is designed. And we propose a two-round interpolation method to obtain the filter's initial dimensions. The designed filter is fabricated on a double-sided YBCO/MgO/YBCO high-temperature superconducting (HTS) thin film for demonstration. The experimental results show that this UWB filter produces eight resonances in the passband eventually, which effectively improves the in-band reflection and the band-edge steepness. Moreover, the upper stopband performance is enhanced due to the transmission zeros (TZs) generated by the stepped-impedance open stubs and the ICL structure. The measured good performance verifies the practicability of the two-round interpolation approach, which can also be extended to other odd-even-mode filter designs.     

High-resolution semi-discrete Hermite central-upwind scheme for multidimensional Hamilton–Jacobi equations

We introduce a high resolution fifth-order semi-discrete Hermite central-upwind scheme for multidimensional Hamilton–Jacobi equations. The numerical fluxes of the scheme are constructed by Hermite polynomials which can be obtained by using the short-time assignment of the first derivatives. The extensions of the proposed semi-discrete Hermite central-upwind scheme to multidimensional cases are straightforward. The accuracy, efficiency and stability properties of our schemes are finally demonstrated via a variety of numerical examples.     

Kinematics of layered reinforced-concrete planar beam finite elements with embedded transversal cracking

In this work crack formation and development is addressed and implemented in a planar layered reinforced-concrete beam element. The crack initiation and growth is described using the strength criterion in conjunction with exact kinematics of the interlayer connection. In this way a novel embedded-discontinuity beam finite element is derived in which the tensile stresses in concrete at the crack position reaching the tensile strength will trigger a crack to open. Since the element is multi-layered, in this way the crack is allowed to propagate through the depth of the beam. The cracked layer(s) will involve discontinuity in the cross-sectional rotation equal to the crack-profile angle, as well as a discontinuity in the position vector of the layer’s reference line. A bond–slip relationship is superimposed onto this model in a kinematically consistent manner with reinforcement being treated as an additional layer of zero thickness with its own material parameters and a constitutive law implemented in the multi-layered beam element.