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

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

Finite element computations of viscoelastic two-phase flows using local projection stabilization

A three-field local projection stabilized (LPS) finite element method is developed for computations of a three-dimensional axisymmetric buoyancy driven liquid drop rising in a liquid column where one of the liquid is viscoelastic. The two-phase flow is described by the time-dependent incompressible Navier-Stokes equations, whereas the viscoelasticity is modeled by the Giesekus constitutive equation in a time-dependent domain. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the governing equations in the time-dependent domain. Interface-resolved moving meshes in ALE allows to incorporate the interfacial tension force and jumps in the material parameters accurately. A one-level LPS based on an enriched approximation space and a discontinuous projection space is used to stabilize the numerical scheme. A comprehensive numerical investigation is performed for a Newtonian drop rising in a viscoelastic fluid column and a viscoelastic drop rising in a Newtonian fluid column. The influence of the viscosity ratio, Newtonian solvent ratio, Giesekus mobility factor, and the Eötvös number on the drop dynamics are analyzed. The numerical study shows that beyond a critical Capillary number, a Newtonian drop rising in a viscoelastic fluid column experiences an extended trailing edge with a cusp-like shape and also exhibits a negative wake phenomena. However, a viscoelastic drop rising in a Newtonian fluid column develops an indentation around the rear stagnation point with a dimpled shape.     

Asymptotic results for weighted means of linear combinations of independent Poisson random variables

ABSTRACT

In this paper we prove the large deviation principle for a class of weighted means of linear combinations of independent Poisson distributed random variables, which converge weakly to a normal distribution. The interest in these linear combinations is motivated by the diffusion approximation in Lansky [On approximations of Stein's neuronal model, J. Theoret. Biol. 107 (1984), pp. 631–647] of the Stein's neuronal model (see Stein [A theoretical analysis of neuronal variability, Biophys. J. 5 (1965), pp. 173–194]). We also prove an analogue result for sequences of multivariate random variables based on the diffusion approximation in Tamborrino, Sacerdote, and Jacobsen [Weak convergence of marked point processes generated by crossings of multivariate jump processes. Applications to neural network modeling, Phys. D 288 (2014), pp. 45–52]. The weighted means studied in this paper generalize the logarithmic means. We also investigate moderate deviations.     

An inequality connecting entropy distance,Fisher Information and large deviations

In this paper we introduce a new generalisation of the relative Fisher Information for Markov jump processes on a finite or countable state space, and prove an inequality which connects this object with the relative entropy and a large deviation rate functional. In addition to possessing various favourable properties, we show that this generalised Fisher Information converges to the classical Fisher Information in an appropriate limit. We then use this generalised Fisher Information and the aforementioned inequality to qualitatively study coarse-graining problems for jump processes on discrete spaces.     

Properties of bounded stochastic processes employed in biophysics

Abstract

Realistic stochastic modeling is increasingly requiring the use of bounded noises. In this work, properties and relationships of commonly employed bounded stochastic processes are investigated within a solid mathematical ground. Four families are object of investigation: the Sine-Wiener (SW), the Doering–Cai–Lin (DCL), the Tsallis–Stariolo–Borland (TSB), and the Kessler–Sørensen (KS) families. We address mathematical questions on existence and uniqueness of the processes defined through Stochastic Differential Equations, which often conceal non-obvious behavior, and we explore the behavior of the solutions near the boundaries of the state space. The expression of the time-dependent probability density of the Sine-Wiener noise is provided in closed form, and a close connection with the Doering–Cai–Lin noise is shown. Further relationships among the different families are explored, pathwise and in distribution. Finally, we illustrate an analogy between the Kessler–Sørensen family and Bessel processes, which allows to relate the respective local times at the boundaries.     

基于改进格林元的压裂水平井动态分析方法

以建立高效的动态分析方法为出发点，以边单元作为求解点，改进传统的格林元方法，减少未知数和求解矩阵维度；并提出基于改进格林元的加密网格加密方法，保证考虑复杂裂缝网络的压裂水平井动态模拟的早期精度.退化模型与半解析解、数值模拟结果进行对比，验证本文基于加密网格的改进格林元方法的准确性和动态分析的高效性.最后进行动态响应的敏感性分析，结果表明：①格林元方法是一种高精度的动态模拟方法，将求解节点设置在网格的边上可以提高压裂水平井动态模拟的速度；②改进格林元方法的加密基于叠加原理，不需要通过插值近似，其求解精度高.在相同加密网格条件下，基于本文改进格林元方法的加密效果比有限差分加密效果更佳；③复杂裂缝导流能力、改造区渗透率提高倍数、改造区大小等参数对压裂水平井动态特征影响较大，在动态分析和参数反演时，应着重考虑这些因素的影响.     

基于六面体网格的向量式有限元分析及应用

基于向量式有限元基本原理，给出了八节点六面体等参实体单元的基本公式，通过投影方式将空间曲面六面体转换为投影六面体，采用参考面的逆向运动求解节点纯变形，通过单元形函的虚功方程计算节点内力；针对坐标模式和内力积分模式等关键问题提出了有效的处理方案。编制了六面体实体单元的数值计算程序，并进行工程结构算例分析。结果表明，所编制程序可有效模拟实体结构的静力、动力及大变形大位移行为分析，验证了本文理论和程序的正确性和实用性。     

Opposed type double stage cell for Mbar pressure experiment with large sample volume

ABSTRACT

A new opposed type double-stage large volume cell has been developed to compress large volume samples to more than 100?GPa (Mbar) pressure. A pair of second-stage diamond anvils is introduced into the first-stage Paris–Edinburgh press. The double-stage large volume cell allows the generation of ultrahigh pressures using a large culet diameter of the second-stage diamond anvils (diameters of 0.5–1.2?mm). Pressure generation up to 131?GPa has been achieved by using the culet diameter of 0.5?mm. Sample volume of the double-stage large volume cell can be more than ～100 times larger than that of conventional Mbar experiment using a diamond anvil cell. The double-stage large volume cell has a large opening in the horizontal plane for X-ray measurements, which is particularly suited for the multi-angle energy dispersive X-ray diffraction measurement, thus opening a new way of in situ structural determinations of amorphous materials at Mbar pressures.     

Meshless acoustic analysis using a weakly singular Burton-Miller boundary integral formulation

The Burton-Miller boundary integral formulation is solved by a complex variable boundary element-free method (CVBEFM) for the boundary-only meshless analysis of acoustic problems with arbitrary wavenumbers. To regularize both strongly singular and hypersingular integrals and to avoid the computation of the solid angle and its normal derivative, a weakly singular Burton-Miller formulation is derived by considering the normal derivative of the solid angle and adopting the singularity subtraction procedures. To facilitate the implementation of the CVBEFM and the approximation of gradients of the boundary variables, a stabilized complex variable moving least-square approximation is selected in the meshless discretization procedure. The results show the accuracy and efficiency of the present CVBEFM and reveal that the method can produce satisfactory results for all wavenumbers, even for extremely large wavenumbers such as k = 10 000.     

The coupled-cluster formalism – a mathematical perspective

ABSTRACT

The Coupled-Cluster (CC) theory is one of the most successful high precision methods used to solve the stationary Schrödinger equation. In this article, we address the mathematical foundation of this theory with focus on the advances made in the past decade. Rather than solely relying on spectral gap assumptions (non-degeneracy of the ground state), we highlight the importance of coercivity assumptions – Gårding type inequalities – for the local uniqueness of the CC solution. Based on local strong monotonicity, different sufficient conditions for a local unique solution are suggested. One of the criteria assumes the relative smallness of the total cluster amplitudes (after possibly removing the single amplitudes) compared to the Gårding constants. In the extended CC theory the Lagrange multipliers are wave function parameters and, by means of the bivariational principle, we here derive a connection between the exact cluster amplitudes and the Lagrange multipliers. This relation might prove useful when determining the quality of a CC solution. Furthermore, the use of an Aubin–Nitsche duality type method in different CC approaches is discussed and contrasted with the bivariational principle.