Imposing accurate wall boundary conditions in corrective‐matrix‐based moving particle semi‐implicit method for free surface flow

Corrective matrix that is derived to restore consistency of discretization schemes can significantly enhance accuracy for the inside particles in the Moving Particle Semi‐implicit method. In this situation, the error due to free surface and wall boundaries becomes dominant. Based on the recent study on Neumann boundary condition (Matsunaga et al, CMAME, 2020), the corrective matrix schemes in MPS are generalized to straightforwardly and accurately impose Neumann boundary condition. However, the new schemes can still easily trigger instability at free surface because of the biased error caused by the incomplete/biased neighbor support. Therefore, the existing stable schemes based on virtual particles and conservative gradient models are applied to free surface and nearby particles to produce a stable transitional layer at free surface. The new corrective matrix schemes are only applied to the particles under the stable transitional layer for improving the wall boundary conditions. Three numerical examples of free surface flows demonstrate that the proposed method can help to reduce the pressure/velocity fluctuations and hence enhance accuracy further.     

Deciphering the Structure and Chemical Composition of Drug Nanocarriers: From Bulk Approaches to Individual Nanoparticle Characterization

Drug nanocarriers (NCs) with sizes usually below 200 nm are gaining increasing interest in the treatment of severe diseases such as cancer and infections. Characterization methods to investigate the morphology and physicochemical properties of multifunctional NCs are key in their optimization and in the study of their in vitro and in vivo fate. Whereas a variety of methods has been developed to characterize “bulk” NCs in suspension, the scope of this review is to describe the different approaches for the NC characterization on an individual basis, for which fewer techniques are available. The accent is put on methods devoid of labelling, which could lead to artefacts. For each characterization method, the principles and approaches to analyze the data are presented in an accessible manner. Aspects related to sample preparation to avoid artefacts are indicated, and emphasis is put on examples of applications. NC characterization on an individual basis allows gaining invaluable information in terms of quality control, on: i) NC localization and fate in biological samples; ii) NC morphology and crystallinity; iii) distribution of the NC components (drugs, shells), and iv) quantification of NCs’ chemical composition. The individual characterization approaches are expected to gain increasing interest in the near future.     

托卡马克中局域磁扰动引起的高能量离子损失研究

基于单粒子导心运动代码ORBIT，采用测试粒子模拟方法，研究了托卡马克等离子体内部不同径向位置处局域磁场扰动对高能量离子的损失的影响。研究表明，在局域磁扰动主要分布在某磁面附近、其环向具有类似纹波场形式下，可造成一些靠近等离子体中心区域的高能量离子损失，但对靠近等离子边界的离子损失影响相对不大。这些损失的高能量离子均为捕获离子，离子的投掷角越大就越容易损失。此外，造成高能量离子最大损失率的局域场径向位置与这些损失离子的初始径向位置通常存在一定的偏移，而且这个偏移与这些离子的能量密切相关。当局域场出现在某些位置时，能量较低的离子会有一定的损失，能量较高的离子反而不会损失。     

A high‐order PIC method for advection‐dominated flow with application to shallow water waves

A hybrid Eulerian‐Lagrangian particle‐in‐cell–type numerical method is developed for the solution of advection‐dominated flow problems. Particular attention is given over to the high‐order transfer of flow properties from the particles to the grid. For smooth flows, the method presented is of formal high‐order accuracy in space. The method is applied to solve the nonlinear shallow water equations resulting in a new, and novel, shock capturing shallow water solver. The approach is able to simulate complex shallow water flows, which can contain an arbitrary number of discontinuities. Both trivial and nontrivial bottom topography is considered, and it is shown that the new scheme is inherently well balanced, exactly satisfying the ‐property. The scheme is verified against several one‐dimensional benchmark shallow water problems. These include cases that involve transcritical flow regimes, shock waves, and nontrivial bathymetry. In all the test cases presented, very good results are obtained.     

Space potential particles to enhance the stability of projection-based particle methods

Particle methods have been seldom verified by a Karman vortex simulation, which is commonly performed as a typical benchmark in computational fluid dynamics. This is mainly due to a difficulty in suppression of occurrence of unphysical voids manifested usually in a strong vortex on account of definition of free surface by the Lagrangian tracking framework with inconsistency in volume conservation. This paper presents a simple and effective scheme as a free-surface boundary condition of projection-based particle methods, namely the MPS (moving particle semi-implicit) and Incompressible SPH (ISPH) methods to handle the free surface with consistency in volume conservation. The new scheme is introduced into the Poisson pressure equation (PPE) with consideration of a potential in void space as space potential particle (SPP), to reproduce physical motions of particles around free surface through a particle–void interaction. The enhancing effect of the newly proposed SPP scheme is shown by simulating a few numerical tests, including a whirling water flow, a two-phase surfacing flow, and a set of Karman vortex simulations.     

Bammann-Chiesa-Johnson粘塑性本构模型材料参数的一种识别方法

粘塑性本构模型能否成功模拟金属高应变率大应变变形过程依赖于材料参数识别结果的好坏。由于BCJ模型考虑了应变率、温度与材料硬化之间的耦合效应以及应变率、温度历史效应,同时模型中包含了多个材料参数,因此很难通过试验直接识别模型的材料参数。本文针对BCJ模型中的耦合效应和历史效应,基于对模型中材料参数物理涵义的界定,给出了一种对材料参数解耦、分离并进行估计的方法,获得了模型材料参数估计公式,估计了材料参数的取值范围。在此基础上,编制了BCJ模型应力积分径向返回算法和粒子群优化算法的计算程序,应用重新设计了BCJ模型耦合效应和历史效应的反分析方法,在参数取值范围内对材料参数进行了优化识别。以OFHC Cu为例,应用提出的识别方法对BCJ模型的材料参数进行了识别,计算结果和试验结果符合较好。     

Modified Lagrangian vortex method with improved boundary conditions for water waves past a thin bottom‐standing barrier

Herein, the modified Lagrangian vortex method (LVM), a hybrid analytical‐numerical algorithm per se, is devised to simulate the process of vortex formation and shedding from the sharp edge of a zero‐thickness vertical plate under linear water‐wave attack. Application of the Helmholtz decomposition facilitates a convenient switch between the inviscid‐ and viscous‐flow models, thereby enabling easy incorporation of vorticity effects into the potential‐flow calculations for the viscous‐dominated region. In evaluating the potential‐flow component, making good use of the quickly convergent technique with singular basis functions, correctly capturing the singular behavior in velocity fields near the tip of the plate, leads to a considerable reduction of computational burdens and to 12‐decimal‐place accuracy. The viscous correction is carried out via the meshless LVM with improved boundary conditions. Comparisons with previously published results show good agreement. Simulations of vortex generation and evolution illuminate the ability of the present method, and provide a supplement to pertinent experimental works. The hybrid scheme proposed herein allows flexibility for the former LVM and convenience in the code development. Such a compromise fits particularly well for the high‐resolution modeling of sharp‐edged vortex shedding without heavy numerical developments. Copyright © 2014 John Wiley & Sons, Ltd.     

颗粒毛细效应影响因素的离散元分析

颗粒毛细效应是指将一根细管插入填充有颗粒物质的容器中并对管施加竖直振动时颗粒在管内上升并最终达到一个稳定的高度的现象, 该现象为颗粒物料的逆重力输运提供了一种潜在的技术途径. 为探究颗粒毛细效应的影响因素, 采用离散元方法, 模拟再现了颗粒毛细效应过程,展示了不同管径下颗粒竖直方向速度演变特性, 考察了不同容器宽度和振动条件下颗粒最终毛细上升高度随管径的演变规律. 结果表明, 在容器宽度与粒径比为40、管振幅与粒径比为14.33、管振动频率为12 Hz情况下, 管径与粒径比$D/d = 3.33$时, 管内颗粒堵塞严重, 使得颗粒上升缓慢,并造成颗粒柱中断; $D/d = 8.33$时, 起初毛细上升高度增加迅速, 随后毛细上升高度的增大逐渐减缓, 管内颗粒在管径方向几乎不存在速度梯度; $D/d =15$时, 随着颗粒毛细上升高度的增大, 管内颗粒柱分离为速度截然不同的两层, 上层颗粒在管径方向几乎不存在速度梯度, 而下层颗粒存在明显的速度梯度.研究还发现, 在毛细效应能够发生的管径范围内, 存在一个对应于颗粒最终毛细上升高度最大值的临界管径, 当管径小于临界管径时, 颗粒最终毛细上升高度随管径的增大而增大, 当管径大于临界管径时, 颗粒最终毛细上升高度随管径的增大而趋于减小; 增大容器宽度,临界管径有所增大; 增大振幅、适当提高频率能够有效促进临界管径的增大.      

基于粒子滤波的运动目标光电定位仿真研究

运动目标的光电定位不能像静止目标那样简单做均值滤波，鉴于此，引入粒子滤波算法，它不仅可以应用于线性系统，而且还适用于非线性系统。结合光电定位需求，详细推导了计算公式及初值和参数选取公式，对只含测量噪声以及含有测量和运动噪声等的海面运动目标光电无源定位算法进行了仿真计算，验证了算法的有效性，讨论了噪声强度对滤波效果的影响，滤波参数选择对滤波效果的影响，目标运动方式对滤波跟随性的影响，重采样算法对滤波效果的影响等。所得结论为:粒子滤波可用于运动目标光电定位过程，可有效降低定位误差；粒子滤波算法具有较强鲁棒性，适用于噪声较大、目标运动形态变化大等情况。