对流扩散方程的QC3无网格法

对流扩散方程作为偏微分运动方程的分支,在流体力学、气体动力学等领域有着重要应用.为解决对流扩散方程难以通过解析法得到解析解的难题,采用二阶一致3点积分(Quadratically Consistent 3-Point Integration,简称QC3)提高无网格法的计算效率,通过对积分点上形函数导数的修正,改善无网格法的精度和收敛性.本文将QC3无网格法拓展到对流扩散方程问题中,时域离散采用广义特征线Galerkin法,空间离散采用QC3法.数值结果表明,应用QC3无网格法得到的对流扩散问题数值解与解析解十分接近,验证了QC3无网格法解决对流扩散问题的可行性.     

More efficient time integration for Fourier pseudospectral DNS of incompressible turbulence

Time integration of Fourier pseudospectral DNS is usually performed using the classical fourth-order accurate Runge-Kutta method or other second- or third-order methods, with a fixed step size. We investigate the use of higher-order Runge-Kutta pairs and automatic step size control based on local error estimation. We find that the fifth-order accurate Runge-Kutta pair of Bogacki and Shampine gives much greater accuracy at a significantly reduced computational cost. Specifically, we demonstrate speedups of 2× to 10× for the same accuracy. Numerical tests (including the Taylor-Green vortex, Rayleigh-Taylor instability, and homogeneous isotropic turbulence) confirm the reliability and efficiency of the method. We also show that adaptive time stepping provides a significant computational advantage for some problems (like the development of a Rayleigh-Taylor instability) without compromising accuracy.     

Effects of van der Waals force on nonlinear vibration of electromechanical integrated electrostatic harmonic actuator

An electromechanical integrated electrostatic harmonic actuator is promising for the miniaturization of electromechanical devices. As the dimensions of the actuator decrease, the effects of the van der Waals force become obvious. In this study, by considering the nonlinearity of electrostatic and van der Waals forces, nonlinear vibration equations of the flexible ring of an electrostatic harmonic actuator are deduced. Using these equations, the nonlinear free vibration and nonlinear forced response of the actuator are investigated. The effects of the van der Waals force on the nonlinear vibration of the flexible ring are analyzed. Results show that these effects of the van der Waals force are relatively obvious under some conditions and should be considered.     

Generating Excess Protons in Microsolvated Acid Clusters under Ambient Conditions: An Issue of Configurational Entropy versus Internal Energy

Acid dissociation, and thus liberation of excess protons in small water droplets, impacts on diverse fields such as interstellar, atmospheric or environmental chemistry. At cryogenic temperatures below 1 K, it is now well established that as few as four water molecules suffice to dissociate the generic strong acid HCl, yet temperature-driven recombination sets in simply upon heating that cluster. Here, the fundamental question is posed of how many more water molecules are required to stabilize a hydrated excess proton at room temperature. Ab initio path integral simulations disclose that not five, but six water molecules are needed at 300 K to allow for HCl dissociation independently from nuclear quantum effects. In order to provide the molecular underpinnings of these observations, the classical and quantum free energy profiles were decomposed along the dissociation coordinate in terms of the corresponding internal energy and entropy profiles. What decides in the end about acid dissociation, and thus ion pair formation, in a specific microsolvated water cluster at room temperature is found to be a fierce competition between classical configurational entropy and internal energy, where the former stabilizes the undissociated state whereas the latter favors dissociation. It is expected that these are generic findings with broad implications on acid–base chemistry depending on temperature in small water assemblies.     

非齐次动力学方程的一种精细积分单步方法

针对非齐次动力学方程■,结合精细积分法和微分求积法,利用同阶的显式龙格-库塔法对计算过程中待求的v_(k+i/s)(i=1,2,…,s)进行预估,提出了一种避免状态矩阵求逆的高效精细积分单步方法。该方法采用精细积分法计算e~(Ht),而Duhamel积分项采用s级s阶的时域微分求积法,计算格式统一且易于编程,可灵活实现变阶变步长。仿真结果表明,与其他单步法及预估校正-辛时间子域法进行数值比较,该方法具有高精度、高效率及良好的稳定性,在求解大规模动力系统时间响应问题中具有较大的优势。     

Cover Image

This paper presents a framework for incorporating arbitrary implicit multistep schemes into the lattice Boltzmann method. While the temporal discretization of the lattice Boltzmann equation is usually derived using a second-order trapezoidal rule, it appears natural to augment the time discretization by using multistep methods. The effect of incorporating multistep methods into the lattice Boltzmann method is studied in terms of accuracy and stability. Numerical tests for the third-order accurate Adams-Moulton method and the second-order backward differentiation formula show that the temporal order of the method can be increased when the stability properties of multistep methods are considered in accordance with the second Dahlquist barrier.     

A discontinuous Galerkin immersed boundary solver for compressible flows: Adaptive local time stepping for artificial viscosity–based shock-capturing on cut cells

We present a higher-order cut cell immersed boundary method (IBM) for the simulation of high Mach number flows. As a novelty on a cut cell grid, we evaluate an adaptive local time stepping (LTS) scheme in combination with an artificial viscosity–based shock-capturing approach. The cut cell grid is optimized by a nonintrusive cell agglomeration strategy in order to avoid problems with small or ill-shaped cut cells. Our approach is based on a discontinuous Galerkin discretization of the compressible Euler equations, where the immersed boundary is implicitly defined by the zero isocontour of a level set function. In flow configurations with high Mach numbers, a numerical shock-capturing mechanism is crucial in order to prevent unphysical oscillations of the polynomial approximation in the vicinity of shocks. We achieve this by means of a viscous smoothing where the artificial viscosity follows from a modal decay sensor that has been adapted to the IBM. The problem of the severe time step restriction caused by the additional second-order diffusive term and small nonagglomerated cut cells is addressed by using an adaptive LTS algorithm. The robustness, stability, and accuracy of our approach are verified for several common test cases. Moreover, the results show that our approach lowers the computational costs drastically, especially for unsteady IBM problems with complex geometries.     

高中物理“教-学-评一体化”课堂教学探索——以“示波器的原理——带电粒子在电场中的偏转”为例

"教-学-评一体化"的课堂是指围绕教学目标,教师的教学、学生的学习以及教师对学生的评价组成一个有机的、整体的有效课堂.教学"有效"的唯一证据在于目标的达成,在于学生学习结果的质量,在于何以证明学生学会了什么.因此,教学中要关注对学生的评价.本文以"示波器的原理——带电粒子在电场中的偏转"为例,论述在"教-学-评一体化"的课堂中如何用评价促进学生思维发展.     

The differentiation operator in the space of uniformly convergent Dirichlet series

Continuity, compactness, the spectrum and ergodic properties of the differentiation operator are investigated, when it acts in the Fréchet space of all Dirichlet series that are uniformly convergent in all half-planes for each . The properties of the formal inverse of the differentiation are also investigated.     

Micro‐ and Nano‐Technologies for Lipid Bilayer‐Based Ion‐Channel Functional Assays

Ion channel proteins provide gated pores that allow ions to passively flow across cell membranes. Owing to their crucial roles in regulating transmembrane ion flow, ion channel proteins have attracted the attention of pharmaceutical investigators as drug targets for use in the studies of both therapeutics and side effects. In this review, we discuss the current technologies that are used in the formation of ion channel‐integrated bilayer lipid membranes (BLMs) in microfabricated devices as a potential platform for next‐generation drug screening systems. Advances in BLM fabrication methodology have allowed the preparation of BLMs in sophisticated formats, such as microfluidic, automated, and/or array systems, which can be combined with channel current recordings. A much more critical step is the integration of the target channels into BLMs. Current technologies for the functional reconstitution of ion channel proteins are presented and discussed. Finally, the remaining issues of the BLM‐based methods for recording ion channel activities and their potential applications as drug screening systems are discussed.