Three‐dimensional vortex simulation of unsteady flow in hydraulic turbines

On the basis of the Helmholtz decomposition, a grid‐free numerical scheme is provided for the solution of unsteady flow in hydraulic turbines. The Lagrangian vortex method is utilized to evaluate the convection and stretch of the vorticity, and the BEM is used to solve the Neumann problem to define the potential flow. The no‐slip boundary condition is satisfied by generating vortex sticks at the solid surface. A semi‐analytical regularization technique is applied to evaluate the singular boundary surface integrals of the potential velocity and its gradients accurately. The fast multipole method was extended to evaluate the velocity and velocity gradients induced by the discretized vortex blobs in the Lagrangian vortex method. The successful simulation for the unsteady flow through a hydraulic turbine's runner has manifested the effectiveness of the proposed method. Copyright © 2011 John Wiley & Sons, Ltd.     

Topological textures and their bifurcation processes in 2D ferromagnetic thin films

In this paper, by the use of the topological current theory, the topological structures and the dynamic processes in thin-film ferromagnetic systems are investigated directly from the viewpoint of topology. It is found that the topological charge of a thin-film ferromagnetic system can be changed by annihilation or creation processes of opposite polarized vortex–antivortex pairs taking place at space–time singularities of the normalized magnetization vector field of the system, the variation of the topological charge is integer and can further be expressed in terms of the Hopf indices and Brouwer degrees of the magnetization vector field around the singularities. Moreover, the change of the topological charge of the system is crucial to vortex core reversal processes in ferromagnetic thin films. With the help of the topological current theory and implicit function theorem, the processes of vortex merging, splitting as well as vortex core reversal are discussed in detail.     

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.     

Multiferroicity in La1/2Nd1/2FeO3 nanoparticles

Nano-sized La_1/2Nd_1/2FeO₃ (LNF) powder is synthesized by the sol–gel citrate method. The Rietveld refinement of the X-ray diffraction profile of the sample at room temperature (303 K) shows the orthorhombic phase with Pbnm symmetry. The particle size is obtained by transmission electron microscope. The antiferromagnetic nature of the sample is explained using zero field cooled and field cooled magnetisation and the corresponding hysteresis loop. A signature of weak ferromagnetic phase is observed in LNF at low temperature which is explained on the basis of spin glass like behaviour of surface spins. The dielectric relaxation of the sample has been investigated using impedance spectroscopy in the frequency range from 42 Hz to 1 MHz and in the temperature range from 303 K to 513 K. The Cole–Cole model is used to analyse the dielectric relaxation of LNF. The frequency dependent conductivity spectra follow the power law. The magneto capacitance measurement of the sample confirms its multiferroic behaviour.     

Structural,magnetic and optical properties of two concomitant molecular crystals

A new 1D complex has been prepared and characterized. X-ray single crystal structure confirms that the Cu(II) ions assemble in alternating chains with Cu … Cu distances of 2.5685(4) and 3.1760(4) Å. The temperature dependence of the magnetic susceptibility reveals an antiferromagnetic interaction between the paddle-wheel copper centers with an exchange of −300 cm⁻¹. The exchange integral was also determined by quantum chemical ab-initio calculations, using polarised and unpolarised basis sets reproducing well the experimental value. The second harmonic generation efficiency of a concomitantly crystallized material was evaluated and was found to be comparable to urea.     

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₂群理论来进行类型的划分. 此外, 自旋波模式的频率随着方位角指标的变化而线性增加. 模拟结果显示样品的平均交换能量密度明显的高于平均静磁能量密度; 局域交换能量密度主要集中在涡核初始位置, 而局域静磁能量密度主要分布在长轴附近. 交换作用对受限于铁磁薄膜椭圆盘中的单个涡旋态的能量要起主导作用, 从而导致方位角自旋波模式频率随着方位角指标的增加而增加.