A class of finite difference schemes for interface problems with an HOC approach

In this paper, we propose a new methodology for numerically solving elliptic and parabolic equations with discontinuous coefficients and singular source terms. This new scheme is obtained by clubbing a recently developed higher‐order compact methodology with special interface treatment for the points just next to the points of discontinuity. The overall order of accuracy of the scheme is at least second. We first formulate the scheme for one‐dimensional (1D) problems, and then extend it directly to two‐dimensional (2D) problems in polar coordinates. In the process, we also perform convergence and related analysis for both the cases. Finally, we show a new direction of implementing the methodology to 2D problems in cartesian coordinates. We then conduct numerous numerical studies on a number of problems, both for 1D and 2D cases, including the flow past circular cylinder governed by the incompressible Navier–Stokes equations. We compare our results with existing numerical and experimental results. In all the cases, our formulation is found to produce better results on coarser grids. For the circular cylinder problem, the scheme used is seen to capture all the flow characteristics including the famous von Kármán vortex street. Copyright © 2016 John Wiley & Sons, Ltd.     

Generation and Propagation of Optical Superoscillatory Vortex Arrays

Superoscillation is an intriguing wave phenomenon which enables subwavelength features propagating into far field and hence has potential applications in super‐resolution microscopy as well as particle trapping and manipulation. While previous demonstrations mostly concentrate on designing complicated nanostructures for generating uncontrollable superoscillatory functions, here a new technique which allows for creating polynomially shaped superoscillatory functions that contain phase singularity arrays is demonstrated both theoretically and experimentally. Such a technique is implemented in optical experiments for the first time and controllable superoscillatory lobes with feature much below the diffraction limit is achieved. More importantly, a general theoretical framework, which, to our knowledge, has not been reported before, is developed to show how the created superoscillations propagate to a distance of many Rayleigh ranges and eventually disappear when the distance is sufficiently larger. The validity of the model is confirmed by the experiments. The results may trigger further studies in light field shaping and manipulations in subwavelength scale.     

Accuracy of high‐order density‐based compressible methods in low Mach vortical flows

At low Mach numbers, Godunov‐type approaches, based on the method of lines, suffer from an accuracy problem. This paper shows the importance of using the low Mach number correction in Godunov‐type methods for simulations involving low Mach numbers by utilising a new, well‐posed, two‐dimensional, two‐mode Kelvin–Helmholtz test case. Four independent codes have been used, enabling the examination of several numerical schemes. The second‐order and fifth‐order accurate Godunov‐type methods show that the vortex‐pairing process can be captured on a low resolution with the low Mach number correction applied down to 0.002. The results are compared without the low Mach number correction and also three other methods, a Lagrange‐remap method, a fifth‐order accurate in space and time finite difference type method based on the wave propagation algorithm, and fifth‐order spatial and third‐order temporal accurate finite volume Monotone Upwind Scheme for Conservation Laws (MUSCL) approach based on the Godunov method and Simple Low Dissipation Advection Upstream Splitting Method (SLAU) numerical flux with low Mach capture property. The ability of the compressible flow solver of the commercial software, ANSYS FLUENT , in solving low Mach flows is also demonstrated for the two time‐stepping methods provided in the compressible flow solver, implicit and explicit. Results demonstrate clearly that a low Mach correction is required for all algorithms except the Lagrange‐remap approach, where dissipation is independent of Mach number. © 2013 Crown copyright. International Journal for Numerical Methods in Fluids. © 2013 John Wiley & Sons, Ltd.     

Hydrodynamic Limit of the Gross-Pitaevskii Equation

We study dynamics of vortices in solutions of the Gross-Pitaevskii equation i? _t u = Δu + ?^?2 u(1 ? |u|²) on ?² with nonzero degree at infinity. We prove that vortices move according to the classical Kirchhoff-Onsager ODE for a small but finite coupling parameter ?. By carefully tracking errors we allow for asymptotically large numbers of vortices, and this lets us connect the Gross-Pitaevskii equation on the plane to two dimensional incompressible Euler equations through the work of Schochet [19 Schochet , S. ( 1996 ). The point vortex method for periodic weak solutions of the 2D Euler equations . Comm. Pure Appl. Math. 49 : 911 – 965 .[Crossref], [Web of Science ®] , [Google Scholar]].     

椭圆纳米盘中磁涡旋结构的方位角自旋波模式

本文针对坡莫合金椭圆形盘中的磁涡旋结构, 采用微磁学模拟与傅里叶分析相结合的技术研究了磁涡旋自旋波的本征激发模式. 通过沿样品短轴方向施加一面内方向的脉冲磁场, 观察到一系列方位角自旋波模式. 观察到的自旋波模式具有两重对称性, 可以通过C₂群理论来进行类型的划分. 此外, 自旋波模式的频率随着方位角指标的变化而线性增加. 模拟结果显示样品的平均交换能量密度明显的高于平均静磁能量密度; 局域交换能量密度主要集中在涡核初始位置, 而局域静磁能量密度主要分布在长轴附近. 交换作用对受限于铁磁薄膜椭圆盘中的单个涡旋态的能量要起主导作用, 从而导致方位角自旋波模式频率随着方位角指标的增加而增加.     

基于光楔阵列产生光学涡旋阵列的研究

在研究光楔衍射法产生单涡旋的基础上,基于长条形光楔阵列,提出了利用光束阵列衍射产生涡旋阵列的方法.该方法要求光束阵列在平行于光楔边缘方向上的光束间距等于光束直径的整数倍.利用超精密机床采用一体化加工法加工了光楔阵列元件,验证了该方法的可行性.利用空间光调制器快速灵活调整光束阵列的优点,搭建了借助空间光调制器加载达曼光栅衍射产生所需光束阵列的实验光学系统.针对光束阵列与光楔阵列的匹配问题,研究了达曼光栅掩模图基本单元对光束阵列的调控,获得了可调结构的光束阵列.实验产生了拓扑荷一致的光学涡旋阵列,与仿真结果相一致,证明所提方法的有效性.     

纳米磁流体动态散斑干涉场的奇异场分布

通过对由纳米磁流体运动引起的双扫描激光散斑干涉光场及其变化做拉盖尔-高斯滤波下的傅里叶变换,获得动态散斑干涉图对应的光学涡旋分布及变化特征。分析认为,光学涡旋分布及变化对应着由纳米磁微粒及其团族的运动所引起的动态散斑变化。当纳米磁微粒聚集到分散的过程中,动态激光散斑光场的奇异场分布发生相应变化,说明了磁流体运动过程对应涡旋密度有先大后小,再由小变大的两个变化;并且光学涡旋密度高,对应较小颗粒的散斑场,磁流体处于稳态的状况;光学涡旋密度低,对应较大颗粒的散斑场,对应着磁流体激烈的运动。研究结果体现了奇异场分布变化和纳米磁流体动后趋稳的过程存在对应关系。     

高温超导体中磁通涡旋薄饼运动的实验证明

为了验证高温超导体中涡旋薄饼运动的存在,我们设计了电流和电位引线分处样品不同表面的实验。实验结果表明,在这种电流电位引线配置下的E～j关系不同于刚性涡旋线运动所预期的E～j关系,由此支持了涡旋薄饼运动的存在。     

旋流泵全流道三维定常流场的数值模拟

本文在计算模型采用分块非结构化六面体网格的基础上,利用雷诺平均N-S方程,标准k-ε紊流模型,结合 SIMPLEC算法,数值模拟了旋流泵的内部流动,得到了旋流泵压力面、吸力面上的压力分布,同时也估算了叶轮的效率, 比较了旋流泵试验性能与数值模拟性能,有助于了解旋流泵的内部流动特性,指导旋流泵的水力设计。