Simulation of the interaction of light with 3-D metallic nanostructures using a proper orthogonal decomposition-Galerkin reduced-order discontinuous Galerkin time-domain method

In this artice, we report on a reduced-order model (ROM) based on the proper orthogonal decomposition (POD) technique for the system of 3-D time-domain Maxwell's equations coupled to a Drude dispersion model, which is employed to describe the interaction of light with nanometer scale metallic structures. By using the singular value decomposition (SVD) method, the POD basis vectors are extracted offline from the snapshots produced by a high order discontinuous Galerkin time-domain (DGTD) solver. With a Galerkin projection and a second order leap-frog (LF₂) time discretization, a discrete ROM is constructed. The stability condition of the ROM is then analyzed. In particular, when the boundary is a perfect electric conductor condition, the global energy of the ROM is bounded, which is consistent with the characteristics of global energy in the DGTD method. It is shown that the ROM based on Galerkin projection can maintain the stability characteristics of the original high dimensional model. Numerical experiments are presented to verify the accuracy, demonstrate the capabilities of the POD-based ROM and assess its efficiency for 3-D nanophotonic problems.     

Uniform Convergence of Multigrid V-Cycle on Adaptively Refined Finite Element Meshes for Elliptic Problems with Discontinuous Coefficients

The multigrid V-cycle methods for adaptive finite element discretizations of two-dimensional elliptic problems with discontinuous coefficients are considered. Under the conditions that the coefficient is quasi-monotone up to a constant and the meshes are locally refined by using the newest vertex bisection algorithm, some uniform convergence results are proved for the standard multigrid V-cycle algorithm with Gauss-Seidel relaxations performed only on new nodes and their immediate neighbours. The multigrid V-cycle algorithm uses $\mathcal{O}(N)$ operations per iteration and is optimal.     

Flux calculation at the interface between two rock types for two-phase flow in porous media

When simulating two-phase flow in porous media, one has to consider the case where there is a discontinuity in the medium. There relative permeabilities and capillary pressure functions may change and we address the problem of calculating the convective part of the numerical flux at the interface between the two rock types. Several solutions are compared.     

STATIC STUDY OF CANTILEVER BEAM STICTION UNDER ELECTROSTATIC FORCE INFLUENCE

I. INTRODUCTION For capacitor-like microelectromechanical systems (MEMS) structure[1??6], the voltage betweenthe structure and substrate causes attractive force. The sources of the voltage can be an arti?ciallymounted device[2 , 7??9] or the temporar…     

DDA与FEM耦合法在分缝重力坝非线性分析中的应用

混凝土重力坝一般设有数条纵缝。纵缝使坝体的连续性遭到破坏，一般的有限单元法很难处理这样的不连续变形问题。本文首先阐述了DDA(Discontinuous Deformation Analysis)与FEM(Finite Element Method)耦合算法的原理，然后采用作者自行研制的DDA与FEM耦合法程序对一实际重力坝进行了详细的计算分析。重点讨论了纵缝的不同设置形式和缝宽大小对坝体变形和应力分布、坝底面接触应力分布以及坝基面应力分布的影响，为工程设计提供依据。     

一类指数型间断函数的嵌入非连续等参有限元法

在嵌入非连续有限元的基本思想下,提出一类附加位移形函数———指数型间断函数,来模拟由于非连续结构,如裂纹和节理,所导致的位移不连续规律,该附加函数是以到间断处的垂直距离为自变量,且随距离的增大而呈指数衰减的函数.指数型间断函数具有在数学上的便于积分和求导的优点,且比阶梯间断函数更能反映实际破裂后的变形情况.本文用弱解形式推导了嵌入非连续有限元格式,编制了二维4节点和三维8节点的嵌入非连续等参有限元程序,并分别给出了算例.算例表明在模拟裂纹追踪时,指数型间断函数的嵌入非连续等参有限元法可行且有效.     

Mode III fracture analysis by Trefftz boundary element method

This paper presents a hybrid Trefftz (HT) boundary element method (BEM) by using two indirect techniques for mode III fracture problems. Two Trefftz complete functions of Laplace equation for normal elements and a special purpose Trefftz function for crack elements are proposed in deriving the Galerkin and the collocation techniques of HT BEM. Then two auxiliary functions are introduced to improve the accuracy of the displacement field near the crack tips, and stress intensity factor (SIF) is evaluated by local crack elements as well. Furthermore, numerical examples are given, including comparisons of the present results with the analytical solution and the other numerical methods, to demonstrate the efficiency for different boundary conditions and to illustrate the convergence influenced by several parameters. It shows that HT BEM by using the Galerkin and the collocation techniques is effective for mode III fracture problems. The project supported by the National Natural Science Foundation of China(10472082). The English text was polished by Keren Wang.     

饱和-非饱和土壤中污染物运移过程的数值模拟

本文提出了一个模拟饱和 非饱和土壤中溶和污染物运移过程的数值模型．模拟的控制污染物运移的物理 化学现象包括：对流，机械逸散，分子弥散，吸附，蜕变，不动水效应．发展了一个修正的特征线Ｇａｌｅｒｋｉｎ方法以离散污染物运移过程的控制方程并导出了一个用于有限元方程求解的显式算法．数值例题结果表明所提出模型和算法的功能     

Regularization of Discontinuous Vector Fields on $${\mathbb{R}^3}$$ via Singular Perturbation

Singular perturbations problems in dimension three which are approximations of discontinuous vector fields are studied in this paper. The main result states that the regularization process developed by Sotomayor and Teixeira produces a singular problem for which the discontinuous set is a center manifold. Moreover, the definition of sliding vector field coincides with the reduced problem of the corresponding singular problem for a class of vector fields.      

Flows Resulting from the Incidence of a Discontinuous Wave on a Bottom Step

The solvability of the problem of the flows resulting from the incidence of a discontinuous wave on a bottom step is studied using a single-layer shallow water model. Solutions in which the total energy of the flow is conserved at the step and those in which it is lost at the step are considered. Regions of double and triple hystereses in the obtained self-similar solutions are found. An analogy is drawn with the problem of single-layer flow over a bottom obstacle. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 2, pp. 8–22, March–April, 2006.