光纤模式分析的伽辽金有限元模型

提出了采用二阶吸收边界条件的全矢量平面伽辽金有限元模型,用于分析任意横截面形状和各向异性折射率分布的光纤的传导模式和泄漏模式,能精确求出各模式传输常量的实部和虚部以及模场分布,既不出现伪解,又不漏解.推导了各向异性介质全矢量耦合波动方程及其变分形式,给出了基于结点的二阶三角形单元的离散公式和单元矩阵,成功的将二阶吸收边界条件加入外边界二次线性单元的离散公式.计算表明采用该模型分析光子晶体光纤模式有效折射率与采用多极方法和基于离散函数展开的有限差分法所得结果吻合很好,采用二阶吸收边界条件计算限制损耗比一阶吸收边界条件结果精确.     

Influence of the Spatial Model of Lithospheric Plates on an Initial Earthquake

The boundary-value problem for a block structure consisting of an elastic half-space that occupies the lower part of a Cartesian coordinate system and interacts with two block elements in the form of unbounded wedges with right angles at the vertices is considered. The horizontal boundaries of the wedges are rigidly connected with the half-space boundary, and vertical ones are parallel planes. We consider the antiplane problem for cases where the distance between the vertical boundaries of block elements is present or absent. By the example of a three-dimensional problem degenerating into a two-dimensional one, it is shown that an initial earthquake is possible, like in the case of using Kirchhoff plates instead of three-dimensional block elements.     

Atomistic finite element model for axial buckling and vibration analysis of single-layered graphene sheets

In this article, an atomistic model is developed to study the buckling and vibration characteristics of single-layered graphene sheets (SLGSs). By treating SLGSs as space-frame structures, in which the discrete nature of graphene sheets is preserved, they are modeled using three-dimensional elastic beam elements for the bonds. The elastic moduli of the beam elements are determined via a linkage between molecular mechanics and structural mechanics. Based on this model, the critical compressive forces and fundamental natural frequencies of single-layered graphene sheets with different boundary conditions and geometries are obtained and then compared. It is indicated that the compressive buckling force decreases when the graphene sheet aspect ratio increases. At low aspect ratios, the increase of aspect ratios will result in a significant decrease in the critical buckling load. It is also indicated that increasing aspect ratio at a given side length results in the convergence of buckling envelops associated with armchair and zigzag graphene sheets. The influence of boundary conditions will be studied for different geometries. It will be shown that the influence of boundary conditions is not significant for sufficiently large SLGSs.     

The Gaussian Free Field and SLE<Subscript>4</Subscript> on Doubly Connected Domains

The level lines of the Gaussian free field are known to be related to SLE₄. It is shown how this relation allows to define chordal SLE₄ processes on doubly connected domains, describing traces that are anchored on one of the two boundary components. The precise nature of the processes depends on the conformally invariant boundary conditions imposed on the second boundary component. Extensions of Schramm’s formula to doubly connected domains are given for the standard Dirichlet and Neumann conditions and a relation to first-exit problems for Brownian bridges is established. For the free field compactified at the self-dual radius, the extended symmetry leads to a class of conformally invariant boundary conditions parametrised by elements of SU(2). It is shown how to extend SLE₄ to this setting. This allows for a derivation of new passage probabilities à la Schramm that interpolate continuously from Dirichlet to Neumann conditions.     

REFLECTION OF PLANE WAVES BY THE FREE BOUNDARY OF A POROUS ELASTIC HALF-SPACE

In the present paper, we study reflection of inclined incident plane waves from a free boundary of the half-space in which the material is described by constitutive equations valid for elastic solids with voids. Both the cases of the transverse and longitudinal incident waves are considered, and it is shown that only the transverse one can propagate in the solid without attenuation, after having been reflected from the free boundary surface. The reflection coefficient and the amplitude of the surface oscillations are expressed in explicit form. The general results are demonstrated for several hypothetical porous materials, and it is shown that the reflection coefficient and the vibration amplitude are typically less than in classical media without voids. However, for relatively large transverse wave speed and high porosity, free boundary oscillation can exceed the classical one.     

Boundary form factors in finite volume

We describe the volume dependence of matrix elements of local boundary fields to all orders in inverse powers of the volume. Using the scaling boundary Lee–Yang model as testing ground, we compare the matrix elements extracted from boundary truncated conformal space approach to exact form factors obtained using the bootstrap method. We obtain solid confirmation for the boundary form factor bootstrap, which is different from all previously available tests in that it is a non-perturbative and direct comparison of exact form factors to multi-particle matrix elements of local operators, computed from the Hamiltonian formulation of the quantum field theory.     

Studying the buckling and vibration characteristics of single-walled zinc oxide nanotubes using a nanoscale finite element model

The free vibration and axial buckling of achiral zinc oxide nanotubes (ZnONTs) are studied in this paper based on a three-dimensional finite-element model in which bonds are modeled using beam elements and mass elements are placed at the joints of beams instead of atoms. To determine the mechanical properties of the nanotubes, a linkage is established between molecular mechanics and density functional theory. The fundamental frequency and critical buckling load of ZnONTs with different geometries, chiralities and boundary conditions are calculated. It is shown that zigzag nanotubes are more stable than armchair ones. Investigating the effect of aspect ratio on the critical force shows that longer nanotubes are less stable. Also, it is indicated that increasing the length of the nanotubes will result in decreasing the frequency. Moreover, as the aspect ratio increases, the effect of end conditions diminishes.     

Sum rules for the vertex function from local commutativity (II) Relations for the causal boundary values in mometum space

The sum rules derived in an earlier paper [1] are applied to the momentum-space boundary values of the vertex function and it is shown how to use the formulas to derive relations for the dynamical on-the-mass-shell quantities, i.e. the matrix elements of locally commuting field operators.     

On the vertical large eddy breakup device capability to decrease the turbulent drag

The possibility of decreasing the turbulent friction with the use of streamwise-aligned vertical large eddy breakup devices installed normal to the surface of a flat plate in an incompressible equilibrium turbulent boundary layer with a nominally gradientless flow past this plate is studied experimentally. The Reynolds number based on the boundary-layer momentum thickness is 1099 in the cross section where these vertical devices are mounted. It is shown that elements of this geometry are effective tools for modification of the turbulent boundary layer from the viewpoint of both the gain in friction drag and the suppression of turbulent fluctuations of velocity and, hence, can be considered as one of the most important factors of controlling the structure of the near-wall turbulence in the flow past a flat surface.     

Three-dimensional finite-element simulation of Zener pinning dynamics

The Zener pinning dynamics of a moving boundary interacting with one or more particles is described by a three-dimensional (3D) finite-element model. The model, based upon a variational formulation for boundary motion by viscous drag, is solved by a finite-element method to obtain the velocity at each node of triangular linear elements on the grain boundary. It is first applied to relatively simple and validated cases, for which analytical and numerical results are available. These cases correspond to an axisymmetrical geometry, in which the grain boundary interacts with a centred particle. A simple analytical pinning criterion is derived from these simulations. The model is then applied to general 3D cases, in which the grain boundary interacts with arbitrarily localized and sized particles. The aim of these 3D simulations is to quantify the influence of the position and the number of particles on the average grain-boundary velocity. It is shown, for example, that the drag effect is enhanced when the particle, or the cluster of particles, is off-centre and that pinning is less efficient with several particles than with a single particle producing the same Zener force.     

A finite element approach for predicting nozzle admittances

A finite element method is used to predict the admittances of axisymmetric nozzles. It is assumed that the flow in the nozzle is isentropic and irrotational, and the disturbances are small so that linear analyses apply. An approximate, two dimensional compressible model is used to describe the steady flow in the nozzle. The propagation of acoustic disturbances is governed by the complete linear wave equation. The differential form of the acoustic equation is transformed to an integral equation by using Galerkin's method, and Green's theorem is applied so that the acoustic boundary conditions can be introduced through the boundary residuals. The boundary conditions are described for both straight and curved sonic lines. A two dimensional FEM with linear elements is used to solve the acoustic equation. A one dimensional FEM is also used to solve the reduced equation of Crocco, and the solution verifies the sufficiency of the boundary residual formulation. Comparison between computed admittances and experimental data is shown to be quite good.     

Nested T-Duality

We identify the obstructions for Kiritsis–Obers T-duality of boundary WZW models. The open string duality pattern is much richer than in the closed strings case since it depends substantially on the geometry of branes. In particular, the duality obstructions disappear for certain brane configurations associated to non-regular elements of the Cartan torus. It is shown in this case that the boundary WZW model is “nested” in the twisted boundary WZW model as the dynamical subsystem of the latter     

A numerical comparison between two unconstrained variational formulations

In an effort to relieve the often cumbersome burden of meeting the requirements of the end conditions and to unify the solution formulation for boundary- and initial-value problems, unconstrained variational statements have been introduced in conjunction with some approximate methods. In the case of a boundary value problem, it is shown in this paper that two different variational statements can be established: one is arrived at by the use of the Lagrange multipliers, the other by energy considerations. The numerical convergence of the solutions associated with finite element schemes involving use of one of these two different variational statements is compared with that of the other. In the case of an initial value problem, both formulations can again be established when the adjoint field variable and the adjoint variational statement are introduced. The numerical data presented here indicate that while both methods generate excellent convergent results for the boundary problem, the method of stiff springs yields results which show much better convergence for the initial value problem than those achieved by Lagrange multipliers.     

On the connection between the solutions to the Dirac and Weyl equations and the corresponding electromagnetic four-potentials

In this work we study in detail the connection between the solutions to the Dirac and Weyl equations and the associated electromagnetic four-potentials.First,it is proven that all solutions to the Weyl equation are degenerate,in the sense that they correspond to an infinite number of electromagnetic four-potentials.As far as the solutions to the Dirac equation are concerned,it is shown that they can be classified into two classes.The elements of the first class correspond to one and only one four-potential,and are called non-degenerate Dirac solutions.On the other hand,the elements of the second class correspond to an infinite number of four-potentials,and are called degenerate Dirac solutions.Further,it is proven that at least two of these fourpotentials are gauge-inequivalent,corresponding to different electromagnetic fields.In order to illustrate this particularly important result we have studied the degenerate solutions to the forcefree Dirac equation and shown that they correspond to massless particles.We have also provided explicit examples regarding solutions to the force-free Weyl equation and the Weyl equation for a constant magnetic field.In all cases we have calculated the infinite number of different electromagnetic fields corresponding to these solutions.Finally,we have discussed potential applications of our results in cosmology,materials science and nanoelectronics.     

A low numerical dissipation immersed interface method for the compressible Navier–Stokes equations

A numerical method to solve the compressible Navier–Stokes equations around objects of arbitrary shape using Cartesian grids is described. The approach considered here uses an embedded geometry representation of the objects and approximate the governing equations with a low numerical dissipation centered finite-difference discretization. The method is suitable for compressible flows without shocks and can be classified as an immersed interface method. The objects are sharply captured by the Cartesian mesh by appropriately adapting the discretization stencils around the irregular grid nodes, located around the boundary. In contrast with available methods, no jump conditions are used or explicitly derived from the boundary conditions, although a number of elements are adopted from previous immersed interface approaches. A new element in the present approach is the use of the summation-by-parts formalism to develop stable non-stiff first-order derivative approximations at the irregular grid points. Second-order derivative approximations, as those appearing in the transport terms, can be stiff when irregular grid points are located too close to the boundary. This is addressed using a semi-implicit time integration method. Moreover, it is shown that the resulting implicit equations can be solved explicitly in the case of constant transport properties. Convergence studies are performed for a rotating cylinder and vortex shedding behind objects of varying shapes at different Mach and Reynolds numbers.     

Cr,Cu掺杂Ag电子结构和光学特性的第一性原理研究

基于密度泛函理论的第一性原理平面波超软赝势方法,分别计算了Cr、Cu掺杂Ag晶胞的电子结构和光学特性,计算结果表明,掺杂后Cr、Cu外层电子分别与Ag外层电子发生作用,价带会出现一定展宽;掺杂体系的反射谱和吸收谱在不同波段范围内会发生不同变化,体系的光谱特性发生改变.研究结果对于光电子器件制备及薄膜界面特性研究具有指导意义.     

Cr,Cu掺杂Ag电子结构和光学特性的第一性原理研究

基于密度泛函理论的第一性原理平面波超软赝势方法,分别计算了Cr、Cu掺杂Ag晶胞的电子结构和光学特性,计算结果表明,掺杂后Cr、Cu外层电子分别与Ag外层电子发生作用,价带会出现一定展宽;掺杂体系的反射谱和吸收谱在不同波段范围内会发生不同变化,体系的光谱特性发生改变. 研究结果对于光电子器件制备及薄膜界面特性研究具有指导意义.     

Boundary element method for bandgap computation of photonic crystals

A boundary element method for computing bandgap structures of two-dimensional photonic crystals is developed. For photonic crystals composed of a square or triangular lattice of parallel cylinders with arbitrarily shaped cross-sections, the boundary integral equations are formulated for a unit cell. Constant boundary elements are adopted to discretize the boundaries. Applying the periodic boundary conditions and the interface conditions, we obtain a linear eigenvalue equation with relatively small matrices. The solution of the eigenvalue equation yields the Bloch wave vectors for given frequencies. The convergence of the method for the desired accuracy and efficiency is examined by some typical numerical examples. It is shown that the present method is efficient and accurate and thus provides a flexible treatment of electromagnetic waves in periodic structures with inclusions of arbitrary shape.     

Boundary Quantum Entanglement of the XXZ Spin Chain with Boundary Impurities

The boundary quantum entanglement for the s = 1/2 X X Z spin chain with boundary impurities is studied via the density matrix renormalization group （DMRG） method. It is shown that the entanglement entropy of the boundary bond （the impurity and the chain spin next to it） behaves differently in different phases. The relationship between the singular points of the boundary entropy and boundary quantum critical points is discussed.     

Proposal of structures possessing high exciton concentration

The possibility of achievement of high exciton concentrations is analyzed. It is shown that high concentrations can be achieved in a three-layer thin molecular film due to the autoreduction processes taking place in it. Shortly, the appearance of high concentrations is the consequence of boundary conditions in film and of the magnitude of matrix elements of dipol-dipol interactions. The autoreduction takes place in the cases when matrix elements characterizing exciton transfer are less than statistical matrix elements. Based on numerical analysis, it was found that optical quanta concentrations of a three-layer film can achieve values of about 5×10⁻². The structures possessing so high concentration do not exist in nature, thus they have to be synthesised. For the current state of nanotechnology, it is not a problem. Fortunately achieving high concentrations requires only certain ratios of relevant characteristics of the film with a two-level exciton scheme, but not their single values.