Two-scale sparse finite element approximations

To reduce computational cost,we study some two-scale finite element approximations on sparse grids for elliptic partial differential equations of second order in a general setting.Over any tensor product domain ?R~d with d = 2,3,we construct the two-scale finite element approximations for both boundary value and eigenvalue problems by using a Boolean sum of some existing finite element approximations on a coarse grid and some univariate fine grids and hence they are cheaper approximations.As applications,we obtain some new efficient finite element discretizations for the two classes of problem:The new two-scale finite element approximation on a sparse grid not only has the less degrees of freedom but also achieves a good accuracy of approximation.     

Numerical Simulation for Convection-Diffusion Problem with Periodic Micro-Structure

In this article, the convection dominated convection-diffusion problems with the periodic micro-structure are discussed. A two-scale finite element scheme based on the homogenization technique for this kind of problems is provided. The error estimates between the exact solution and the approximation solution, of the homogenized equation or the two-scale finite element scheme are analyzed. It is shown that the scheme provided in this article is convergent for any fixed diffusion coefficient 5, and it may be convergent independent of δ under some conditions. The numerical results demonstrating the theoretical results are presented in this article.     

周期结构区域中压电问题的双尺度有限元方法

近年来纤维压电复合材料的力电性能预测已发展为一个重要的研究领域．对力电耦合周期结构的复合材料问题,通过引入匹配的边界层得到了电势与位移解的新型双尺度有限元计算方法,建立了电势与位移的双尺度耦合关系,分析了双尺度有限元解的误差.数值算例验证了方法的有效性．     

具有小周期孔隙复合材料弹性结构的双尺度有限元分析

对于具有小周期孔隙复合材料弹性结构,在双尺度渐近分析理论结果的基础上提出了双尺度有限元计算格式,并给出了严格的误差估计.     

整周期复合材料弹性结构的有限元计算

1．引言周期性复合材料与周期结构的弹性力学问题，由于其材料特征剧烈振荡，且周期很小，问题相当复杂，除了一些特殊的和简单的问题可以用解析法求解外，多数问题很难或不可能用解析法求解，需要采用数值方法计算，有限元法无疑是最有效的方法之一．用细观力学方法研究复合材料的力学问题时，需要涉及纤维的排列情况，纤维和基体的性能，界面的分布情况，以及细观的几何参数和物理参数等．由于复合材料细观构造的不均匀性和不规则性，损伤和缺陷的存在，以及许多难以精确测定的因素，使得复合材料的细观力学问题十分复杂，不作出一些简化…     

The heterogeneous multiscale method as a framework for coupling the finite element method and molecular dynamics

Hierarchical two-scale methods are computationally very powerful as there is no direct coupling between the macro- and microscale. Such schemes develop first a microscale model under macroscopic constraints, then the macroscopic constitutive laws are found by averaging over the microscale. The heterogeneous multiscale method (HMM) is a general top-down approach for the design of multiscale algorithms. While this method is mainly used for concurrent coupling schemes in the literature, the proposed methodology also applies to a hierarchical coupling. This contribution discusses a hierarchical two-scale setting based on the heterogeneous multi-scale method for quasi-static problems: the macroscale is treated by continuum mechanics and the finite element method and the microscale is treated by statistical mechanics and molecular dynamics. Our investigation focuses on an optimised coupling of solvers on the macro- and microscale which yields a significant decrease in computational time with no associated loss in accuracy. In particular, the number of time steps used for the molecular dynamics simulation is adjusted at each iteration of the macroscopic solver. A numerical example demonstrates the performance of the model. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

具有周期振荡系数二阶椭圆型方程的双尺度高精度算法

本文考虑了一类具有周期振荡系数二阶椭圆型方程边值问题,提出了基于双尺度渐近展开式的高精度有限元算法,并给出严格的证明。     

Two-Scale Simulation of a Disc With a Crack by X-FEM and Cell Model

The two-scale simulation of a linear-elastic orthotropic disc with a central crack under mode-I loading may be used to verify the extended finite element method implementation of orthotropic enrichment functions into finite element codes such as FEAP. The stress distribution on the finer scale is simultaneously resolved by the high fidelity generalized method of cells called at each integration point of the macro elements. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

ERROR ESTIMATES FOR TWO-SCALE COMPOSITE FINITE ELEMENT APPROXIMATIONS OF NONLINEAR PARABOLIC EQUATIONS

We study spatially semidiscrete and fully discrete two-scale composite finite element method for approximations of the nonlinear parabolic equations with homogeneous Dirich-let boundary conditions in a convex polygonal domain in the plane.This new class of finite elements,which is called composite finite elements,was first introduced by Hackbusch and Sauter[Numer.Math.,75(1997),pp.447-472]for the approximation of partial differential equations on domains with complicated geometry.The aim of this paper is to introduce an efficient numerical method which gives a lower dimensional approach for solving par-tial differential equations by domain discretization method.The composite finite element method introduces two-scale grid for discretization of the domain,the coarse-scale and the fine-scale grid with the degrees of freedom lies on the coarse-scale grid only.While the fine-scale grid is used to resolve the Dirichlet boundary condition,the dimension of the finite element space depends only on the coarse-scale grid.As a consequence,the resulting linear system will have a fewer number of unknowns.A continuous,piecewise linear composite finite element space is employed for the space discretization whereas the time discretization is based on both the backward Euler and the Crank-Nicolson methods.We have derived the error estimates in the L∞(L2)-norm for both semidiscrete and fully discrete schemes.Moreover,numerical simulations show that the proposed method is an efficient method to provide a good approximate solution.     

Multiscale computational method for thermoelastic problems of composite materials with orthogonal periodic configurations

This study develops a novel multiscale computational method for thermoelastic problems of composite materials with orthogonal periodic configurations. Firstly, the multiscale asymptotic analysis for these multiscale problems is given successfully, and the formal second-order two-scale approximate solutions for these multiscale problems are constructed based on the above-mentioned analysis. Then, the error estimates for the second-order two-scale (SOTS) solutions are obtained. Furthermore, the corresponding SOTS numerical algorithm based on finite element method (FEM) is brought forward in details. Finally, some numerical examples are presented to verify the feasibility and effectiveness of our multiscale computational method. Moreover, our multiscale computational method can accurately capture the local thermoelastic responses in composite block structure, plate, cylindrical and doubly-curved shallow shells.     

基于二重网格的定常Navier-Stokes方程的局部和并行有限元算法

对二维定常的不可压缩的Navier-Stokes方程的局部和并行算法进行了研究．给出的算法是多重网格和区域分解相结合的算法,它是基于两个有限元空间:粗网格上的函数空间和子区域的细网格上的函数空间．局部算法是在粗网格上求一个非线性问题,然后在细网格上求一个线性问题,并舍掉内部边界附近的误差相对较大的解．最后,基于局部算法,通过有重叠的区域分解而构造了并行算法,并且做了算法的误差分析,得到了比标准有限元方法更好的误差估计,也对算法做了数值试验,数值结果通过比较验证了本算法的高效性和合理性．     

TWO-SCALE FEM FOR ELLIPTIC MIXED BOUNDARY VALUE PROBLEMS WITH SMALL PERIODIC COEFFICIENTS

1. IntroductionComPosite materials have been widely used in high tedrilogy ewineering as well as or-dinary industrial products since they have mW elegat qllallties, such as high strength, bostffeess, high temperature resistance, corrosion resistance, and fatgh resistance. Most of thecomPosite materiaIs have small periodic condgurations. Thu8 the static analysis of the struc-tures of comPosite materials usually leads to the boundary valu Problezns of eiliPtic patialdtherelitial eqllations wi…     

三维高速冲击动力有限元滑移面算法

本文给出三维高速冲击动力有限元滑移面算法.该算法不但能保证单元结点的动量守恒、动量矩守恒,而且由实例计算表明该算法在处理高速穿、破甲过程中是稳定和有效的.     

Two-level multiscale finite element methods for the steady navier-stokes problem

In this article, on the basis of two-level discretizations and multiscale finite element method, two kinds of finite element algorithms for steady Navier-Stokes problem are presented and discussed. The main technique is first to use a standard finite element discretization on a coarse mesh to approximate low frequencies, then to apply the simple and Newton scheme to linearize discretizations on a fine grid. At this process, multiscale finite element method as a stabilized method deals with the lowest equal-order finite element pairs not satisfying the inf-sup condition. Under the uniqueness condition, error analyses for both algorithms are given. Numerical results are reported to demonstrate the effectiveness of the simple and Newton scheme.     

TWO-LEVEL MULTISCALE FINITE ELEMENT METHODS FOR THE STEADY NAVIER-STOKES PROBLEM

In this article, on the basis of two-level discretizations and multiscale finite element method, two kinds of finite element algorithms for steady Navier-Stokes problem are presented and discussed. The main technique is first to use a standard finite element discretization on a coarse mesh to approximate low frequencies, then to apply the simple and Newton scheme to linearize discretizations on a fine grid. At this process, multiscale finite element method as a stabilized method deals with the lowest equal-order finite element pairs not satisfying the inf-sup condition. Under the uniqueness condition, error analyses for both algorithms are given. Numerical results are reported to demonstrate the effectiveness of the simple and Newton scheme.     

A finite element model based on statistical two-scale analysis for equivalent heat transfer parameters of composite material with random grains

In this paper, a new finite element model based on statistical two-scale analysis for predicting the equivalent heat transfer parameters of the composite material with random grains is presented and its convergence, its error result and the symmetry, positive property of equivalent heat transfer parameters matrix are also proved. Firstly, some definitions of the probability space and the composite material with random grains are described and the STSA formulation predicting the equivalent heat transfer parameters of the composite material are briefly reviewed. Next, a finite element formulation and its corresponding procedure for the composite material with random grains is described. Then, the convergence, the error estimate and the symmetry, positive property of the equivalent heat transfer parameters matrix computed by FE based on STSA are proved. The numerical result shows the validity of the FE model based on STSA and the convergence and the symmetry, positive property of the equivalent heat transfer parameter matrix of the composite material with random grains by the FE model.     

Local and parallel finite element algorithms based on two-grid discretizations for the transient Stokes equations

Based on two-grid discretizations, some local and parallel finite element algorithms for the d-dimensional (d = 2,3) transient Stokes equations are proposed and analyzed. Both semi- and fully discrete schemes are considered. With backward Euler scheme for the temporal discretization, the basic idea of the fully discrete finite element algorithms is to approximate the generalized Stokes equations using a coarse grid on the entire domain, then correct the resulted residue using a finer grid on overlapped subdomains by some local and parallel procedures at each time step. By the technical tool of local a priori estimate for the fully discrete finite element solution, errors of the corresponding solutions from these algorithms are estimated. Some numerical results are also given which show that the algorithms are highly efficient.     

The heterogeneous multiscale finite element method for elliptic homogenization problems in perforated domains

In this contribution we analyze a generalization of the heterogeneous multiscale finite element method for elliptic homogenization problems in perforated domains. The method was originally introduced by E and Engquist (Commun Math Sci 1(1):87–132, 2003) for homogenization problems in fixed domains. It is based on a standard finite element approach on the macroscale, where the stiffness matrix is computed by solving local cell problems on the microscale. A-posteriori error estimates are derived in L ²(Ω) by reformulating the problem into a discrete two-scale formulation (see also, Ohlberger in Multiscale Model Simul 4(1):88–114, 2005) and using duality methods afterwards. Numerical experiments are given in order to numerically evaluate the efficiency of the error estimate.     

Numerical Simulation of the Influence of Rough Bone-Callus Interface on the Healing of Fractured Bone

The process of healing of fractured bone is known to be influenced by the mechanical environment and the loads exerted by physical activity of the patient or otherwise. We compute mechanical fields in the soft connective tissue of the healing fracture using Biot's poroelasticity model and a finite element (FE) method for low-frequency loading. A two-scale FE framework is used to model effects of the rough bone-callus contact surface. We look at the difference the interface roughness makes with respect to different possible mechanostimulating agents.     

Error analysis of finite element approximations of the stochastic Stokes equations

Numerical solutions of the stochastic Stokes equations driven by white noise perturbed forcing terms using finite element methods are considered. The discretization of the white noise and finite element approximation algorithms are studied. The rate of convergence of the finite element approximations is proved to be almost first order (h|ln h|) in two dimensions and one half order ( h^\frac12h^{\frac{1}{2}}) in three dimensions. Numerical results using the algorithms developed are also presented.