Physalis method for heterogeneous mixtures of dielectrics and conductors: Accurately simulating one million particles using a PC

Prosperetti’s seminal Physalis method, an Immersed Boundary/spectral method, had been used extensively to investigate fluid flows with suspended solid particles. Its underlying idea of creating a cage and using a spectral general analytical solution around a discontinuity in a surrounding field as a computational mechanism to enable the accommodation of physical and geometric discontinuities is a general concept, and can be applied to other problems of importance to physics, mechanics, and chemistry. In this paper we provide a foundation for the application of this approach to the determination of the distribution of electric charge in heterogeneous mixtures of dielectrics and conductors. The proposed Physalis method is remarkably accurate and efficient. In the method, a spectral analytical solution is used to tackle the discontinuity and thus the discontinuous boundary conditions at the interface of two media are satisfied exactly. Owing to the hybrid finite difference and spectral schemes, the method is spectrally accurate if the modes are not sufficiently resolved, while higher than second-order accurate if the modes are sufficiently resolved, for the solved potential field. Because of the features of the analytical solutions, the derivative quantities of importance, such as electric field, charge distribution, and force, have the same order of accuracy as the solved potential field during postprocessing. This is an important advantage of the Physalis method over other numerical methods involving interpolation, differentiation, and integration during postprocessing, which may significantly degrade the accuracy of the derivative quantities of importance. The analytical solutions enable the user to use relatively few mesh points to accurately represent the regions of discontinuity. In addition, the spectral convergence and a linear relationship between the cost of computer memory/computation and particle numbers results in a very efficient method. In the present paper, the accuracy of the method is numerically investigated by example computations using one dielectric particle, one isolated conductor particle, one conductor particle connected to an external source with imposed voltage, and two conductor/dielectric particles with strong interactions. The efficiency of the method is demonstrated with one million particles, which suggests that the method can be used for many important engineering applications of broad interest.     

Observations of the atmospheric boundary layer across the land-sea transition zone using a scanning Mie lidar

Optical properties of the atmospheric boundary layer (ABL) above the land-sea transition interface were measured using a scanning Mie lidar located 30 km away from the Adriatic coast. Based on the two-dimensional range-height-indicator scans, detailed information on the ABL was obtained, including parameters such as atmospheric optical depth, aerosol extinction coefficient and the height of the ABL. The presented case study indicates that the height of the ABL in the land-sea transition zone and the adjacent mountainous region was changing rapidly due to highly variable atmospheric conditions.     

Spontaneous 2D modulation instability in second harmonic generation process

Two-dimensional modulation instability (2DMI) is experimentally demonstrated in a classical second harmonic generation setup. The spatial spectrum is measured and reveals typical 2DMI bands, in agreement with the analytical MI model. These observations are confirmed by (2 + 1)D numerical simulations.     

一阶速度-应力Biot双相各向同性介质弹性波

提出一种等价的一阶双曲型速度一应力Biot双相各向同性介质弹性波波动方程,以实现双相介质混合波场中纯快慢纵波和纯横波波场分离的问题.应用散度和旋度理论证明双相介质等价方程波场分离的可行性,采用高阶交错网格有限差分法构建高精度正演算子,推导其PML吸收边界条件和稳定性条件,并对均匀双相介质和层状非均匀双相介质模型进行数值...     

透明陶瓷透光性能的影响因素

 研究了影响激光透明陶瓷透光性能的主要因素,讨论了陶瓷内部气孔和杂质颗粒等散射粒子、晶界结构中晶界折射率与晶粒折射率的差异以及晶界表面粗糙度等因素对陶瓷透光性能的影响,并定量分析了激光陶瓷透过率随气孔尺寸、气孔率、晶粒相对晶界折射率以及晶界表面粗糙度的变化关系。结果表明： 陶瓷的透过率随着气孔率的减小而增大,但透过率随气孔尺寸的增大而呈现出周期性振荡,且当气孔尺寸与入射光波长可比拟时,陶瓷的透过率会明显降低;在晶界结构中,晶界的折射率与晶粒的折射率相差越小,陶瓷的透过率就越高;晶界表面粗糙度越大,透过率越低。然而,晶界折射率不同于晶粒折射率,这使得其陶瓷透过率降低的程度比对晶界表面粗糙度的影响明显得多。在陶瓷制备过程中,需要重点排除尺寸与入射波长可比拟的气孔, 以抑制晶界结构中第二相的产生。     

透明激光陶瓷散射损耗分析

 为了定量分析激光陶瓷中散射损耗对其透过率的影响,通过建立气孔尺寸分布模型,引入第二相体积比概念,并结合Mie散射、瑞利散射和全散射积分等理论,讨论了激光陶瓷中气孔、晶界第二相和表面粗糙度等引起的散射损耗对激光陶瓷透光性能的影响。研究结果表明：气孔率的大小将明显影响陶瓷透过率,且气孔尺寸分布决定了透过率包络的变化趋势;晶界和表面散射对透过率的影响主要集中在短波长处;在气孔率较低情况下,晶界第二相的存在是导致短波长处透过率急剧降低的主要因素。     

Modelling of multi-bodies in close proximity under water waves——Fluid resonance in narrow gaps

Viscous fluid model and potential flow model with and without artificial damping force(f=-μV,μ the damping coefficient and V the local averaging flow velocity) are employed in this work to investigate the phenomenon of fluid resonance in narrow gaps between multi-bodies in close proximity under water waves.The numerical results are compared with experimental data available in the literature.The comparison demonstrates that both the viscous fluid model and the potential flow model are able to predict the res...     

Solution to the MHD flow over a non-linear stretching sheet by homotopy perturbation method

In this study,by means of homotopy perturbation method(HPM) an approximate solution of the magnetohydrodynamic(MHD) boundary layer flow is obtained.The main feature of the HPM is that it deforms a difficult problem into a set of problems which are easier to solve.HPM produces analytical expressions for the solution to nonlinear differential equations.The obtained analytic solution is in the form of an infinite power series.In this work,the analytical solution obtained by using only two terms from HPM soluti...     

量子系统保结构计算新进展

本文主要介绍量子系统保结构计算最新进展情况,分以下几部分内容:哈密顿系统的辛算法、适合于量子系统的哈密顿量显含时间的辛算法、A2B模型分子和双原子分子系统的经典轨迹辛算法计算、双原子分子CO在激光场中的经典轨迹的辛算法计算及其振动和解离、定态Schr dinger方程的辛形式及求解定态Schr dinger方程本征值问题的辛 打靶法、含时Schr dinger方程的保结构算法及其在激光原子物理中的应用、伪分立态模型、强激光与原子相互作用的渐近边界条件、"非齐线性正则方程"的辛算法及其在计算强激光场中一维原子的多光子电离和高次谐波发射中的应用以及Heisenberg方程的保结构计算等等。     

在时谐电磁场中用矢量势A直接表达电场强度E的一个显式

根据麦克斯韦方程组和洛伦兹条件,对时谐变场e导出用矢量势A直接表达电场强度E的一个显式,并对它进行了讨论.