声学超表面的非对称声分束特性研究

非对称声分束超表面是由人工微单元结构按照特定序列构建的二维平面结构,可将垂直入射的声波分成两束传播方向和分束比自由调控的透射波,在声功能器件设计及声通信领域具有广泛的应用前景。本文系统研究了一种实现非对称声分束的设计理论和实现方法,基于局域声功率守恒条件研究了声分束器的设计理论、阻抗矩阵分布、法向声强分布、声压场分布等。利用遗传算法对四串联共振腔结构进行参数优化实现了声分束器所需的阻抗矩阵分布,声压场分布表明声波入射到声分束器后在入射侧激发出两列传播方向相反且幅值和衰减系数均相同的表面波,实现了入射侧与透射侧的局域声功率相互匹配。声波经过声分束器后被分为两束透射波,两束透射波的折射角和透射系数与理论值十分吻合,证明了设计理论及实现方法的正确性和可行性。本文的研究工作可以为新型非对称声分束结构设计提供理论参考、设计方法和技术支持,并促进其在工程领域的实际应用。     

Effect of wall slip on the viscoelastic particle ordering in a microfluidic channel

The formation of a line of equally spaced particles at the centerline of a microchannel, referred as “particle ordering,” is desired in several microfluidic applications. Recent experiments and simulations highlighted the capability of viscoelastic fluids to form a row of particles characterized by a preferential spacing. When dealing with non-Newtonian fluids in microfluidics, the adherence condition of the liquid at the channel wall may be violated and the liquid can slip over the surface, possibly affecting the ordering efficiency. In this work, we investigate the effect of wall slip on the ordering of particles suspended in a viscoelastic liquid by numerical simulations. The dynamics of a triplet of particles in an infinite cylindrical channel is first addressed by solving the fluid and particle governing equations. The relative velocities computed for the three-particle system are used to predict the dynamics of a train of particles flowing in a long microchannel. The distributions of the interparticle spacing evaluated at different slip coefficients, linear particle concentrations, and distances from the channel inlet show that wall slip slows down the self-assembly mechanism. For strong slipping surfaces, no significant change of the initial microstructure is observed at low particle concentrations, whereas strings of particles in contact form at higher concentrations. The detrimental effect of wall slip on viscoelastic ordering suggests care when designing microdevices, especially in case of hydrophobic surfaces that may enhance the slipping phenomenon.     

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.     

Regulating Phase Junction and Oxygen Vacancies of TiO2 Nanoarrays for Boosted Photoelectrochemical Water Oxidation

In this study, we have provided a facile solution to synthesize well-aligned titanium dioxide nanorods by using hydrothermal reaction. By calcining the materials under different atmospheres and temperatures, a batch of titanium dioxides with excellent oxygen evolution reaction(OER) catalytic efficiency were obtained. This new structured TiO₂ photoanode material yields a high photocurrent density of 5.69 mA/cm² at 1.23 V vs. reversible hydrogen electrode(RHE) under simulated solar light(100 mW/cm²). Surface photovoltage techniques and other measurements were carried out to confirm that the enhanced photoelectrochemical performances were attributed to the synergistic effect of the phase junction and a certain content of surface states, which accelerate the separation and transmission of the photogenerated charges. This material with phase junction and surface states promises a potential application in the field of photoelectric catalysis under solar light.     

A modified SPH method to model the coalescence of colliding non-Newtonian liquid droplets

This paper develops a modified smoothed particle hydrodynamics (SPH) method to model the coalescence of colliding non-Newtonian liquid droplets. In the present SPH, a van der Waals (vdW) equation of state is particularly used to represent the gas-to-liquid phase transition similar to that of a real fluid. To remove the unphysical behavior of the particle clustering, also known as tensile instability, an optimized particle shifting technique is implemented in the simulations. To validate the numerical method, the formation of a Newtonian vdW droplet is first tested, and it clearly demonstrates that the tensile instability can be effectively removed. The method is then extended to simulate the head-on binary collision of vdW liquid droplets. Both Newtonian and non-Newtonian fluid flows are considered. The effect of Reynolds number on the coalescence process of droplets is analyzed. It is observed that the time up to the completion of the first oscillation period does not always increase as the Reynolds number increases. Results for the off-center binary collision of non-Newtonian vdW liquid droplets are lastly presented. All the results enrich the simulations of the droplet dynamics and deepen understandings of flow physics. Also, the present SPH is able to model the coalescence of colliding non-Newtonian liquid droplets without tensile instability.     

Monte Carlo simulations of magnetic particle suspensions with a simple assessment method for the particle overlap between magnetic spheroids

We have developed a simple assessment method for the overlap between spheroidal particles, which neither requires the complex manipulation of vectors and matrices that is indispensable in the ordinary methods, nor is based on a model potential. Moreover, we have developed an evaluation method for the interaction energy arising from the overlap of the steric layer coating spheroidal particles. This is based on a sphere-connected particle model, but some modifications are introduced in order to express an appropriate repulsive interaction energy at the deepest overlapping position. We have investigated the phase change in a magnetic spheroidal particle suspension for a two-dimensional system by means of Monte Carlo simulations. In the case of no external magnetic field, if the magnetic particle-particle interaction is sufficiently strong to favour cluster formation, long raft-like clusters tend to be formed in a dilute situation. With decreasing values of area fraction, a chain-like structure in a dense situation transforms into a raft-like structure within a narrow range of the particle area fraction. Similarly, the raft-like clusters are preferred in a weak applied magnetic field, but an increase in the field strength induces a phase change from a raft-like into a chain-like structure.

Highlights of the present paper:

1. A simple assessment method has been proposed for the overlap between two spheroidal particles.

2. The particle overlap assessment is free from a complex mathematical manipulation regarding vectors and matrices.

3. A modified sphere-connected model has been proposed in order to more accurately evaluate a repulsive interaction due to the overlap of the steric layers coating spheroidal particles.

4. 2D Monte Carlo simulations have been performed to elucidate the phenomenon of a phase change by magnetic spheroidal particles on a material plane surface.

5. A phase change between a raft-like and a chain-like aggregate structure is able to be controlled by the area fraction of particles and an external magnetic field.

     

利用电磁法研究HMX与TATB混合钝感炸药的冲击起爆特性

为研究含有少量奥克托金(HMX)且以三氨基三硝基苯(TATB)为基的高能钝感炸药PBX-3的冲击起爆反应增长规律,采用火炮驱动蓝宝石飞片的方法和铝基组合式电磁粒子速度计技术进行了一维平面冲击实验。通过实验测量撞击表面及内部不同深度处的冲击波后粒子速度,得到PBX-3炸药的Hugoniot关系。根据冲击波示踪器所测数据绘制了炸药到爆轰的时间-距离(x-t)图,获得了反映炸药冲击起爆性能的Pop关系。将入射压力为12.964 GPa时达到爆轰的6条速度曲线修整成相同零点,通过读取6条曲线的分离点即反应区末端的C-J点,计算出化学反应区时间和宽度。     

运用三维离散元技术模拟落锤撞击下奥克托今颗粒的点火燃烧过程

炸药颗粒的点火燃烧过程一直是人们关注的热点问题。近年来，三维离散元技术在中尺度观测颗粒材料的动力学过程中拥有显著优势。炸药燃烧属于颗粒材料的反应动力学，运用三维离散元技术(DM3)可以有效地观测炸药燃烧传播的过程。以奥克托今(HMX)颗粒为例，本文成功模拟并观测到了HMX颗粒的燃烧反应程度，确定了颗粒开始燃烧反应的时间，以及燃烧反应传播的时间。同时，结合落锤冲击颗粒的三维图像以及其表观压强和放热功率，得到了HMX颗粒燃烧反应、燃烧传播的整个反应动力学过程，包括颗粒在冲击加载下碎化塑性变形的过程，颗粒燃烧反应放热的过程，落锤回弹颗粒喷射的过程等。同时，进一步说明了尖顶颗粒更利于颗粒点火，平顶颗粒有抑制颗粒点火的能力。     

Finite-difference calculation of donor energy levels in a spherical quantum dot subject to a magnetic field

The ground and excited states of a donor impurity at the center of a spherical quantum dot subject to a magnetic field are calculated within the effective-mass approximation. The barriers are infinitely high and the differential equation is solved by combining the finite-difference method with the Richardson extrapolation. The binding and transition energies are more accurate than the available variational values, and excellent agreement is found with the hydrogen atom. The transition energies for a medium-size quantum dot are given.