Semi-analytical finite element method applied for characterizing micropolar fibrous composites

A semi-analytical finite element method(SAFEM),based on the two-scale asymptotic homogenization method(AHM)and the finite element method(FEM),is implemented to obtain the effective properties of two-phase fiber-reinforced composites(FRCs).The fibers are periodically distributed and unidirectionally aligned in a homogeneous matrix.This framework addresses the static linear elastic micropolar problem through partial differential equations,subject to boundary conditions and perfect interface contact conditions.The mathematical formulation of the local problems and the effective coefficients are presented by the AHM.The local problems obtained from the AHM are solved by the FEM,which is denoted as the SAFEM.The numerical results are provided,and the accuracy of the solutions is analyzed,indicating that the formulas and results obtained with the SAFEM may serve as the reference points for validating the outcomes of experimental and numerical computations.

Asymptotic Behavior of the Finite Difference and the Finite Element Methods for Parabolic Equations

The asymptotic convergence of the solution of the parabolic equation is proved. By the eigenvalues estimation, we obtain that the approximate solutions by the finite difference method and the finite element method are asymptotically convergent. Both methods are considered in continnous time.

ASYMPTOTICAL STABILITY OF NEUTRAL REACTION-DIFFUSION EQUATIONS WITH PCAS AND THEIR FINITE ELEMENT METHODS

This paper focuses on the analytical and numerical asymptotical stability of neutral reaction-diffusion equations with piecewise continuous arguments.First,for the analytical solutions of the equations,we derive their expressions and asymptotical stability criteria.Second,for the semi-discrete and one-parameter fully-discrete finite element methods solving the above equations,we work out the sufficient conditions for assuring that the finite element solutions are asymptotically stable.Finally,with a typical example with numerical experiments,we illustrate the applicability of the obtained theoretical results.

FINITE ELEMENT ANALYSIS OF TRASINENT TWO DIMENSIONAL CHANNEL FLOW WITH A SUDDEN EXPANSION

The finite element method is employed to simulate incompressible viscous unsteady flow in two dimensional channel with a sudden expansion. The streamline patterns and the velocity vector plots are presented at different expansion ratios and at different moment. The results obtained have certain significance to analyze the formation of eddies and energy loss in a sudden expansion channel flow and other complex channel flow.

Highly efficient H^1-Galerkin mixed finite element method (MFEM) for parabolic integro-differential equation

A highly efficient H1-Galerkin mixed finite element method （MFEM） is presented with linear triangular element for the parabolic integro-differential equation. Firstly, some new results about the integral estimation and asymptotic expansions are studied. Then, the superconvergence of order O（h^2） for both the original variable u in H1 （Ω） norm and the flux p = u in H（div, Ω） norm is derived through the interpolation post processing technique. Furthermore, with the help of the asymptotic expansions and a suitable auxiliary problem, the extrapolation solutions with accuracy O（h^3） are obtained for the above two variables. Finally, some numerical results are provided to confirm validity of the theoretical analysis and excellent performance of the proposed method.

Asymptotic expansions of finite element solutions to Robin problems in H^3 and their application in extrapolation cascadic multigrid method

For the Poisson equation with Robin boundary conditions,by using a few techniques such as orthogonal expansion(M-type),separation of the main part and the finite element projection,we prove for the first time that the asymptotic error expansions of bilinear finite element have the accuracy of O(h3)for u∈H3.Based on the obtained asymptotic error expansions for linear finite elements,extrapolation cascadic multigrid method(EXCMG)can be used to solve Robin problems effectively.Furthermore,by virtue of Richardson not only the accuracy of the approximation is improved,but also a posteriori error estimation is obtained.Finally,some numerical experiments that confirm the theoretical analysis are presented.

Finite Element Analysis on Plane Stress Crack Growth in a Power-Law Hardening Material

Aim To determine numerically the field characteristics in the vied at the tip of a place crack growing steadily in a power-law hardening material. Meteods. Methods on the Euler mode and small-scale yield assumption, the numerical results were given by nonlinear finite element analysis. Results The numerical results of the shape of the active plastic sone, the angular distribution of stresseses and Clack tip opening displacement (CTOD) in the vicinity at the hp of the steadily groWing CraCk are determined. Conclusion The comparison between the numerical results given by the present wort and those given by analytic asymptotic analysis shows that the present work reached a very high accuracy.

The Numerical Accuracy Analysis of Asymptotic Homogenization Method and Multiscale Finite Element Method for Periodic Composite Materials

In this paper,we discuss the numerical accuracy of asymptotic homogenization method(AHM)and multiscale finite element method(MsFEM)for periodic composite materials.Through numerical calculation of the model problems for four kinds of typical periodic composite materials,the main factors to determine the accuracy of first-order AHM and second-order AHM are found,and the physical interpretation of these factors is given.Furthermore,the way to recover multiscale solutions of first-order AHM and MsFEM is theoretically analyzed,and it is found that first-order AHM and MsFEM provide similar multiscale solutions under some assumptions.Finally,numerical experiments verify that MsFEM is essentially a first-order multiscale method for periodic composite materials.

NEW ERROR EXPANSION FOR ONE-DIMENSIONAL FINITE ELEMENTS AND ULTRACONVERGENCE

Based on an improved orthogonal expansion in an element, a new error expression of n-degree finite element approximation uh to two-point boundary value problem is derived, and then superconvergence of two order for both function and derivatives are obtained.

Neumann stochastic finite element method for calculating temperature field of frozen soil based on random field theory

To study the effect of uncertain factors on the temperature field of frozen soil, we propose a method to calculate the spatial average variance from just the point variance based on the local average theory of random fields. We model the heat transfer coefficient and specific heat capacity as spatially random fields instead of traditional random variables. An analysis for calculating the random temperature field of seasonal frozen soil is suggested by the Neumann stochastic finite element method, and here we provide the computational formulae of mathematical expectation, variance and variable coefficient. As shown in the calculation flow chart, the stochastic finite element calculation program for solving the random temperature field, as compiled by Matrix Laboratory （MATLAB） sottware, can directly output the statistical results of the temperature field of frozen soil. An example is presented to demonstrate the random effects from random field parameters, and the feasibility of the proposed approach is proven by compar- ing these results with the results derived when the random parameters are only modeled as random variables. The results show that the Neumann stochastic finite element method can efficiently solve the problem of random temperature fields of frozen soil based on random field theory, and it can reduce the variability of calculation results when the random parameters are modeled as spatial- ly random fields.

AN EFFICIENT GPU ACCELERATION FORMAT FOR FINITE ELEMENT ANALYSIS

This paper proposes a new Graphics Processing Unit(GPU)-accelerated storage format to speed up Sparse Matrix Vector Products(SMVPs) for Finite Element Method(FEM) analysis of electromagnetic problems.A new format called Modified Compile Time Optimization(MCTO) format is used to reduce much execution time and design for hastening the iterative solution of FEM equations especially when rows have uneven lengths.The MCTO-applied FEM is about 10 times faster than conventional FEM on a CPU,and faster than other row-major ordering formats on a GPU.Numerical results show that the proposed GPU-accelerated storage format turns out to be an excellent accelerator.

High-Order Spectral Stochastic Finite Element Analysis of Stochastic Elliptical Partial Differential Equations

This study presents an experiment of improving the performance of spectral stochastic finite element method using high-order elements. This experiment is implemented through a two-dimensional spectral stochastic finite element formulation of an elliptic partial differential equation having stochastic coefficients. Deriving this spectral stochastic finite element formulation couples a two-dimensional deterministic finite element formulation of an elliptic partial differential equation with generalized polynomial chaos expansions of stochastic coefficients. Further inspection of the performance of resulting spectral stochastic finite element formulation with adopting linear and quadratic (9-node or 8-node) quadrilateral elements finds that more accurate standard deviations of unknowns are surprisingly predicted using quadratic quadrilateral elements, especially under high autocorrelation function values of stochastic coefficients. In addition, creating spectral stochastic finite element results using quadratic quadrilateral elements is not unacceptably time-consuming. Therefore, this study concludes that adopting high-order elements can be a lower-cost method to improve the performance of spectral stochastic finite element method.

Elastic Stress Predictor for Stochastic Finite Element Problems

The paper presents a new algorithm of elastic stress predictor in non linear stochastic finite element method using the Generalized Polynomial Chaos. The statistical moments of strains calculated based on the displacement Polynomial Chaos expansion. To descretise the stochastic process of material the Karhunen-Loeve Expansion was used and it is presented. Using the strains and the material Karhunen-Loeve Expansion the stress components are calculated. A numerical example of shallow foundation was carried out and the results of stress and strain of the new algorithm were compared with those raised from Monte Carlo method which is treated as the exact solution. A great accuracy was presented.

Superconvergence and Asymptotic Expansions for Bilinear Finite Volume Element Approximations

Aiming at the isoparametric bilinear finite volume element scheme,we initially derive an asymptotic expansion and a high accuracy combination formula of the derivatives in the sense of pointwise by employing the energy-embedded method on uniform grids.Furthermore,we prove that the approximate derivatives are convergent of order two.Finally,numerical examples verify the theoretical results.

Finite element analysis of expansion-matched submounts for high-power laser diodes packaging

In order to improve the output power and increase the lifetime of laser diodes,expansion-matched submounts were investigated by finite element analysis.The submount was designed as sandwiched structure.By varying the vertical structure and material of the middle layer,the thermal expansion behavior on the mounting surface was simulated to obtain the expansion-matched design.In addition,the thermal performance of laser diodes packaged by different submounts was compared.The numerical results showed that,changing the thickness ratio of surface copper to middle layer will lead the stress and junction temperature to the opposite direction.Thus compromise needs to be made in the design of the vertical structure.In addition,the silicon carbide（SiC） is the most promising material candidate for the middle layer among the materials discussed in this paper.The simulated results were aimed at providing guidance for the optimal design of sandwich-structure submounts.

A SECOND ORDER MODIFIED CHARACTERISTICS VARIATIONAL MULTISCALE FINITE ELEMENT METHOD FOR TIME-DEPENDENT NAVIER-STOKES PROBLEMS

In this paper, by combining the second order characteristics time discretization with the variational multiscale finite element method in space we get a second order modified characteristics variational multiscale finite element method for the time dependent Navier- Stokes problem. The theoretical analysis shows that the proposed method has a good convergence property. To show the efficiency of the proposed finite element method, we first present some numerical results for analytical solution problems. We then give some numerical results for the lid-driven cavity flow with Re = 5000, 7500 and 10000. We present the numerical results as the time are sufficient long, so that the steady state numerical solutions can be obtained.

A Robust Modified Weak Galerkin Finite Element Method for Reaction-Diffusion Equations

In this paper,a robust modified weak Galerkin(MWG)finite element method for reaction-diffusion equations is proposed and investigated.An advantage of this method is that it can deal with the singularly perturbed reaction-diffusion equations.Another advantage of this method is that it produces fewer degrees of freedom than the traditional WG method by eliminating the element boundaries freedom.It is worth pointing out that,in our method,the test functions space is the same as the finite element space,which is helpful for the error analysis.Optimalorder error estimates are established for the corresponding numerical approximation in various norms.Some numerical results are reported to confirm the theory.

A MODIFIED WEAK GALERKIN FINITE ELEMENTMETHOD FOR SINGULARLY PERTURBED PARABOLIC CONVECTION-DIFFUSION-REACTION PROBLEMS

In this work,a modified weak Galerkin finite element method is proposed for solving second order linear parabolic singularly perturbed convection-diffusion equations.The key feature of the proposed method is to replace the classical gradient and divergence operators by the modified weak gradient and modified divergence operators,respectively.We apply the backward finite difference method in time and the modified weak Galerkin finite element method in space on uniform mesh.The stability analyses are presented for both semi-discrete and fully-discrete modified weak Galerkin finite element methods.Optimal order of convergences are obtained in suitable norms.We have achieved the same accuracy with the weak Galerkin method while the degrees of freedom are reduced in our method.Various numerical examples are presented to support the theoretical results.It is theoretically and numerically shown that the method is quite stable.

A FINITE ELEMENT/BOUNDARY ELEMENT——MODIFIED MODAL DECOMPOSITION METHOD FOR VIBRATION AND SOUND RADIATION FROM SUBMERGED SHELL OF REVOLUTION

A finite element / boundary element-modified modal decomposition method (FBMMD) is presented for predicting the vibration and sound radiation from submerged shell of revolution. Improvement has been made to accelerate the convergence to FBMD method by means of introducing the residual modes which take into accaunt the quasi -state contributiort of all neglected modes. As an example, the vibration and sound radiation of a submerged spherical shell excited by axisymmetric force are studied in cases of ka=l,2,3 and 4. From the calculated results we see that the FBMMD method shows a significant improvement to the accuracy of surface sound pressure, normal displacement and directivity patterns of radiating sound, especially to the directivity patterns.

FINITE ELEMENT EXPANSION FOR VARIABLE COEFFICIENT ELLIPTIC PROBLEMS

The Ritz projections of the elliptic problem（?）(aij（?）ju)+c0u=f in Ω,u=0 on （?）Ω,and the quasilinear eiliptic problem（?）（ai（x,Du））+a0（x,Du）=0 in Ω,u=0 on （?）Ωwith linear finite elements admit an asymptotic error expansion,respectively,for certain classes of“uniform” meshes.This provides the theoretical justification for the use of Richardson extrapolationfor increasing the accuracy of the scheme from O（h2） to O（h4 |ln h|）.