堪为物理学第三分支的计算物理

本文是作者1982年在中国访问时所作的讲演,论述物理学各领域的前沿工作,用许多事例说明计算物理对推动物理理论的重大意义。文中还指出计算物理在其性质、方法及需要等方面不同于和独立于解析的理论物理和实验物理,而成为物理学的第三分支。     

Vacuum birefringence in strong magnetic fields: (I) Photon polarization tensor with all the Landau levels

Photons propagating in strong magnetic fields are subject to a phenomenon called the “vacuum birefringence” where refractive indices of two physical modes both deviate from unity and are different from each other. We compute the vacuum polarization tensor of a photon in a static and homogeneous magnetic field by utilizing Schwinger’s proper-time method, and obtain a series representation as a result of double integrals analytically performed with respect to proper-time variables. The outcome is expressed in terms of an infinite sum of known functions which is plausibly interpreted as summation over all the Landau levels of fermions. Each contribution from infinitely many Landau levels yields a kinematical condition above which the contribution has an imaginary part. This indicates decay of a sufficiently energetic photon into a fermion–antifermion pair with corresponding Landau level indices. Since we do not resort to any approximation, our result is applicable to the calculation of refractive indices in the whole kinematical region of a photon momentum and in any magnitude of the external magnetic field.     

A hybrid level set method for modelling detonation and combustion problems in complex geometries

We present an accurate and fast wave tracking method that uses parametric representations of tracked fronts, combined with modifications of level set methods that use narrow bands. Our strategy generates accurate computations of the front curvature and other geometric properties of the front. We introduce data structures that can store discrete representations of the location of the moving fronts and boundaries, as well as the corresponding level set fields, that are designed to reduce computational overhead and memory storage. We present an algorithm we call stack sweeping to efficiently sort and store data that is used to represent orientable fronts. Our implementation features two reciprocal procedures, a forward ‘front parameterization’ that constructs a parameterization of a front given a level set field and a backward ‘field construction’ that constructs an approximation of the signed normal distance to the front, given a parameterized representation of the front. These reciprocal procedures are used to achieve and maintain high spatial accuracy. Close to the front, precise computation of the normal distance is carried out by requiring that displacement vectors from grid points to the front be along a normal direction. For front curves in two dimensions, a cubic interpolation scheme is used, and G ¹ surface parameterization based on triangular patches is used for the three-dimensional implementation to compute the distances from grid points near the front. To demonstrate this new, high accuracy method we present validations and show examples of combustion-like applications that include detonation shock dynamics, material interface motions in a compressible multi-material simulation and the Stephan problem associated with dendrite solidification.     

Issue Information

An integrated finite element method (FEM) is proposed to simulate incompressible two‐phase flows with surface tension effects, and three different surface tension models are applied to the FEM to investigate spurious currents and temporal stability. A Q2Q1 element is adopted to solve the continuity and Navier–Stokes equations and a Q2‐iso‐Q1 to solve the level set equation. The integrated FEM solves pressure and velocity simultaneously in a strongly coupled manner; the level set function is reinitialized by adopting a direct approach using interfacial geometry information instead of solving a conventional hyperbolic‐type equation. In addition, a consistent continuum surface force (consistent CSF) model is utilized by employing the same basis function for both surface tension and pressure variables to damp out spurious currents and to estimate the accurate pressure distribution. The model is further represented as a semi‐implicit manner to improve temporal stability with an increased time step. In order to verify the accuracy and robustness of the code, the present method is applied to a few benchmark problems of the static bubble and rising bubble with large density and viscosity ratios. The Q2Q1‐integrated FEM coupled with the semi‐implicit consistent CSF demonstrates the significantly reduced spurious currents and improved temporal stability. The numerical results are in good qualitative and quantitative agreements with those of the existing studies. Copyright © 2015 John Wiley & Sons, Ltd.     

Evolution of downsized crescent-shaped dune in wind tunnel experiment

The migration of a downsized crescent-shaped dune was investigated in a wind tunnel experiment.Quantified upwind influx and vertical oscillation of the sand bed were introduced to modulate the saturation level of the sand flux above the dune surface to affect dune evolution.The evolution was recorded by top-view photography and then abstracted as the evolution of self-defined characteristic quantities using a digital image processing algorithm.The results showed that,in contrast to the case for spanwise quantities,the evolution of streamwise quantities corresponds to a linear increase in the modulation magnitude more positively and in a monotonic and convergent manner.In contrast with quantities on the windward face,the changes in quantities with respect to the horns were nonmonotonic with time and almost uncorrelated with the variation in modulation strength,which reveals the distinctiveness of leeside evolution.     

大同市居民对室内空气质量认知程度的调查

为了解大同市居民对室内空气质量的认知情况，并分析认知程度与教育的关系，于2016年3月至5月，在大同市城区的9个社区随机对居民进行问卷调查。问卷内容包括居民对室内空气质量安全知识的了解，检测室内空气质量的意愿和治理意愿3个方面。将居民所处地域、居民的年龄与学历对室内空气质量安全认知程度的影响分别进行统计学分析。结果表明，大同市居民室内空气质量安全认知程度总体得分15.8，得分率为52.7%；新城区的居民对室内空气质量安全认知水平高于老城区居民（p<0.05）；学历越高且越年轻的居民对室内空气质量安全的认知水平越高，寻求有关机构进行室内空气检测和治理的意愿更强烈，且对检测和治理的费用承受度较高。大同市居民对室内空气质量安全认知程度总体不高，还应加强对市民的室内空气质量安全方面的教育。     

Maximum likelihood estimation for second level fMRI data analysis with expectation trust region algorithm

The trust region method which originated from the Levenberg–Marquardt (LM) algorithm for mixed effect model estimation are considered in the context of second level functional magnetic resonance imaging (fMRI) data analysis. We first present the mathematical and optimization details of the method for the mixed effect model analysis, then we compare the proposed methods with the conventional expectation-maximization (EM) algorithm based on a series of datasets (synthetic and real human fMRI datasets). From simulation studies, we found a higher damping factor for the LM algorithm is better than lower damping factor for the fMRI data analysis. More importantly, in most cases, the expectation trust region algorithm is superior to the EM algorithm in terms of accuracy if the random effect variance is large. We also compare these algorithms on real human datasets which comprise repeated measures of fMRI in phased-encoded and random block experiment designs. We observed that the proposed method is faster in computation and robust to Gaussian noise for the fMRI analysis. The advantages and limitations of the suggested methods are discussed.     

Temperature dependence of volume and surface symmetry energy coefficients of nuclei

The thermal evolution of the energies and free energies of a set of spherical and near-spherical nuclei spanning the whole periodic table are calculated in the subtracted finite-temperature Thomas–Fermi framework with the zero-range Skyrme-type KDE0 and the finite-range modified Seyler–Blanchard interaction. The calculated energies are subjected to a global fit in the spirit of the liquid-drop model. The extracted parameters in this model reflect the temperature dependence of the volume symmetry and surface symmetry coefficients of finite nuclei, in addition to that of the volume and surface energy coefficients. The temperature dependence of the surface symmetry energy is found to be very substantial whereas that of the volume symmetry energy turns out to be comparatively mild.     

A new level‐dependent coarse grid correction scheme for indefinite Helmholtz problems

In this paper, we construct and analyze a level‐dependent coarse grid correction scheme for indefinite Helmholtz problems. This adapted multigrid (MG) method is capable of solving the Helmholtz equation on the finest grid using a series of MG cycles with a grid‐dependent complex shift, leading to a stable correction scheme on all levels. It is rigorously shown that the adaptation of the complex shift throughout the MG cycle maintains the functionality of the two‐grid correction scheme, as no smooth modes are amplified in or added to the error. In addition, a sufficiently smoothing relaxation scheme should be applied to ensure damping of the oscillatory error components. Numerical experiments on various benchmark problems show the method to be competitive with or even outperform the current state‐of‐the‐art MG‐preconditioned Krylov methods, for example, complex shifted Laplacian preconditioned flexible GMRES. Copyright © 2013 John Wiley & Sons, Ltd.     

A preliminary approach of dynamic identification of slender buildings by neuronal networks

The study of the dynamic behavior of slender masonry structures is usually related to the preservation of the historic heritage. This study, for bell towers and industrial masonry chimneys, is particularly relevant in areas with an important seismic hazard. The analysis of the dynamic behavior of masonry structures is particularly complex due to the multiple effects that can affect the variation of its main frequencies along the seasons of the year: temperature and humidity. Moreover, these dynamic properties also vary considerably in structures built in areas where land subsidence due to the variation of the phreatic level along the year is particularly evident: the stiffness of the soil–structure interaction also varies. This paper presents a study to evaluate the possibility of detecting the variation of groundwater level based on the readings obtained using accelerometers in different positions on the structure. To do this a general case study was considered: a 3D numerical model of a bellower. The variation of the phreatic level was evaluated between 0 and −20 m, and 81 cases studies were developed modifying the rigidity of the soil–structure interaction associated to a position of the phreatic level. To simulate the dispositions of accelerometers on a real construction, 16 points of the numerical model were selected along the structure to obtain modal displacements in two orthogonal directions. Through an adjustment by using neural networks, a good correlation has been observed between the predicted position of the water table and acceleration readings obtained from the numerical model. It is possible to conclude that with a discrete register of accelerations on the tower it is possible to predict the water table depth.