The penalty and hybrid methods are being much used in dealing with the general incompatible element, With the penalty method convergence can always be assured, but comparatively speaking its accuracy is lower, and the...The penalty and hybrid methods are being much used in dealing with the general incompatible element, With the penalty method convergence can always be assured, but comparatively speaking its accuracy is lower, and the condition number and sparsity are not so good. With the hybrid method, convergence can be assured only when the rank condition is satisfied. So the construction of the element is extremely limited. This paper presents the mixed hybrid penalty element method, which combines the two methods together. And it is proved theoretically that this new method is convergent, and it has the same accuracy, condition number and sparsity as the compatible element. That is to say, they are optimal to each other.Finally, a new triangle element for plate bending with nine freedom degrees is constructed with this method (three degreesof freedom are given on each corner -- one displacement and tworotations), the calculating formula of the element stiffness matrix is almost the same as that of the old triangle element for plate bending with nine degrees of freedom But it is converged to true solution with arbitrary irregrlar triangle subdivision. If the true solution u?H3 with this method the linear and quadratic rates of convergence are obtianed for three bending moments and for the displacement and two rotations respectively.展开更多
In this paper, we prove the convergence of the nodal expansion method, a new numerical method for partial differential equations and provide the error estimates of approximation solution.
Torsional guided waves have been widely utilized to inspect the surface corrosion in pipelines due to their simple displacement behaviors and the ability of longrange transmission.Especially,the torsional mode T(0,1),...Torsional guided waves have been widely utilized to inspect the surface corrosion in pipelines due to their simple displacement behaviors and the ability of longrange transmission.Especially,the torsional mode T(0,1),which is the first order of torsional guided waves,plays the irreplaceable position and role,mainly because of its non-dispersion characteristic property.However,one of the most pressing challenges faced in modern quality inspection is to detect the surface defects in pipelines with a high level of accuracy.Taking into account this situation,a quantitative reconstruction method using the torsional guided wave T(0,1)is proposed in this paper.The methodology for defect reconstruction consists of three steps.First,the reflection coefficients of the guided wave T(0,1)scattered by different sizes of axisymmetric defects are calculated using the developed hybrid finite element method(HFEM).Then,applying the boundary integral equation(BIE)and Born approximation,the Fourier transform of the surface defect profile can be analytically derived as the correlative product of reflection coefficients of the torsional guided wave T(0,1)and the fundamental solution of the intact pipeline in the frequency domain.Finally,reconstruction of defects is precisely performed by the inverse Fourier transform of the product in the frequency domain.Numerical experiments show that the proposed approach is suitable for the detection of surface defects with arbitrary shapes.Meanwhile,the effects of the depth and width of surface defects on the accuracy of defect reconstruction are investigated.It is noted that the reconstructive error is less than 10%,providing that the defect depth is no more than one half of the pipe thickness.展开更多
An iterative solution of linear systems is studied,which arises from the discretization of a wire antennas attached with dielectric objects by the hybrid finite-element method and the method of moment (hybrid FEM-MoM...An iterative solution of linear systems is studied,which arises from the discretization of a wire antennas attached with dielectric objects by the hybrid finite-element method and the method of moment (hybrid FEM-MoM).It is efficient to model such electromagnetic problems by hybrid FEM-MoM,since it takes both the advantages of FEM's and MoM's ability.But the resulted linear systems are complicated,and it is hard to be solved by Krylov subspace methods alone,so a two-level preconditioning technique will be studied and applied to accelerate the convergence rate of the Krylov subspace methods.Numerical results show the effectiveness of the proposed two-level preconditioning technique.展开更多
In this paper,a new hybridized mixed formulation of weak Galerkin method is studied for a second order elliptic problem.This method is designed by approximate some operators with discontinuous piecewise polynomials in...In this paper,a new hybridized mixed formulation of weak Galerkin method is studied for a second order elliptic problem.This method is designed by approximate some operators with discontinuous piecewise polynomials in a shape regular finite element partition.Some discrete inequalities are presented on discontinuous spaces and optimal order error estimations are established.Some numerical results are reported to show super convergence and confirm the theory of the mixed weak Galerkin method.展开更多
Background Understanding the interaction between the mitral valve(MV)and the left ventricle(LV)is very important in assessing cardiac pump function,especially when the MV is dysfunctional.Such dysfunction is a major m...Background Understanding the interaction between the mitral valve(MV)and the left ventricle(LV)is very important in assessing cardiac pump function,especially when the MV is dysfunctional.Such dysfunction is a major medical problem owing to the essential role of the MV in cardiac pump function.Computational modelling can provide new approaches to gain insight into the functions of the MV and LV.Methods In this study,a previously developed LV-MV model was used to study cardiac dynamics of MV leaflets under normal and pathological conditions,including hypertrophic cardiomyopathy(HOCM)and calcification of the valve.The coupled LV-MV model was implemented using a hybrid immersed boundary/finite element method to enable assessment of MV haemodynamic performance.Constitutive parameters of the HOCM and calcified valves were inversely determined from published experimental data.The LV compensation mechanism was further studied in the case of the calcified MV.Results Our results showed that MV dynamics and LV pump function could be greatly affected by MV pathology.For example,the HOCM case showed bulged MV leaflets at the systole owing to low stiffness,and the calcified MV was associated with impaired diastolic filling and much-reduced stroke volume.We further demonstrated that either increasing the LV filling pressure or increasing myocardial contractility could enable a calcified valve to achieve near-normal pump function.Conclusion The modelling approach developed in this study may deepen our understanding of the interactions between the MV and the LV and help in risk stratification of heart valve disease and in silico treatment planning by exploring intrinsic compensation mechanisms.展开更多
文摘The penalty and hybrid methods are being much used in dealing with the general incompatible element, With the penalty method convergence can always be assured, but comparatively speaking its accuracy is lower, and the condition number and sparsity are not so good. With the hybrid method, convergence can be assured only when the rank condition is satisfied. So the construction of the element is extremely limited. This paper presents the mixed hybrid penalty element method, which combines the two methods together. And it is proved theoretically that this new method is convergent, and it has the same accuracy, condition number and sparsity as the compatible element. That is to say, they are optimal to each other.Finally, a new triangle element for plate bending with nine freedom degrees is constructed with this method (three degreesof freedom are given on each corner -- one displacement and tworotations), the calculating formula of the element stiffness matrix is almost the same as that of the old triangle element for plate bending with nine degrees of freedom But it is converged to true solution with arbitrary irregrlar triangle subdivision. If the true solution u?H3 with this method the linear and quadratic rates of convergence are obtianed for three bending moments and for the displacement and two rotations respectively.
基金This project is supported by the National Science Foundation of China
文摘In this paper, we prove the convergence of the nodal expansion method, a new numerical method for partial differential equations and provide the error estimates of approximation solution.
基金Project supported by the National Natural Science Foundation of China(Nos.11502108 and 1611530686)the State Key Laboratory of Mechanics and Control of Mechanical Structures at Nanjing University of Aeronautics and Astronautics(NUAA)(No.MCMS-E-0520K02)and the Key Laboratory of Impact and Safety Engineering,Ministry of Education,Ningbo University(No.CJ201904)。
文摘Torsional guided waves have been widely utilized to inspect the surface corrosion in pipelines due to their simple displacement behaviors and the ability of longrange transmission.Especially,the torsional mode T(0,1),which is the first order of torsional guided waves,plays the irreplaceable position and role,mainly because of its non-dispersion characteristic property.However,one of the most pressing challenges faced in modern quality inspection is to detect the surface defects in pipelines with a high level of accuracy.Taking into account this situation,a quantitative reconstruction method using the torsional guided wave T(0,1)is proposed in this paper.The methodology for defect reconstruction consists of three steps.First,the reflection coefficients of the guided wave T(0,1)scattered by different sizes of axisymmetric defects are calculated using the developed hybrid finite element method(HFEM).Then,applying the boundary integral equation(BIE)and Born approximation,the Fourier transform of the surface defect profile can be analytically derived as the correlative product of reflection coefficients of the torsional guided wave T(0,1)and the fundamental solution of the intact pipeline in the frequency domain.Finally,reconstruction of defects is precisely performed by the inverse Fourier transform of the product in the frequency domain.Numerical experiments show that the proposed approach is suitable for the detection of surface defects with arbitrary shapes.Meanwhile,the effects of the depth and width of surface defects on the accuracy of defect reconstruction are investigated.It is noted that the reconstructive error is less than 10%,providing that the defect depth is no more than one half of the pipe thickness.
基金supported by the National Natural Science Foundation of China under Grand No. 10926190, 60973015 the Project of National Defense Key Lab under Grand No. 9140C6902030906
文摘An iterative solution of linear systems is studied,which arises from the discretization of a wire antennas attached with dielectric objects by the hybrid finite-element method and the method of moment (hybrid FEM-MoM).It is efficient to model such electromagnetic problems by hybrid FEM-MoM,since it takes both the advantages of FEM's and MoM's ability.But the resulted linear systems are complicated,and it is hard to be solved by Krylov subspace methods alone,so a two-level preconditioning technique will be studied and applied to accelerate the convergence rate of the Krylov subspace methods.Numerical results show the effectiveness of the proposed two-level preconditioning technique.
文摘In this paper,a new hybridized mixed formulation of weak Galerkin method is studied for a second order elliptic problem.This method is designed by approximate some operators with discontinuous piecewise polynomials in a shape regular finite element partition.Some discrete inequalities are presented on discontinuous spaces and optimal order error estimations are established.Some numerical results are reported to show super convergence and confirm the theory of the mixed weak Galerkin method.
基金This work was supported by the National Natural Science Foundation of China(Grant Nos.11871399,12271440)the UK EPSRC(Grant Nos.EP/S030875,EP/S014284/1,EP/S020950/1,EP/R511705/1,and EP/T017899/1).
文摘Background Understanding the interaction between the mitral valve(MV)and the left ventricle(LV)is very important in assessing cardiac pump function,especially when the MV is dysfunctional.Such dysfunction is a major medical problem owing to the essential role of the MV in cardiac pump function.Computational modelling can provide new approaches to gain insight into the functions of the MV and LV.Methods In this study,a previously developed LV-MV model was used to study cardiac dynamics of MV leaflets under normal and pathological conditions,including hypertrophic cardiomyopathy(HOCM)and calcification of the valve.The coupled LV-MV model was implemented using a hybrid immersed boundary/finite element method to enable assessment of MV haemodynamic performance.Constitutive parameters of the HOCM and calcified valves were inversely determined from published experimental data.The LV compensation mechanism was further studied in the case of the calcified MV.Results Our results showed that MV dynamics and LV pump function could be greatly affected by MV pathology.For example,the HOCM case showed bulged MV leaflets at the systole owing to low stiffness,and the calcified MV was associated with impaired diastolic filling and much-reduced stroke volume.We further demonstrated that either increasing the LV filling pressure or increasing myocardial contractility could enable a calcified valve to achieve near-normal pump function.Conclusion The modelling approach developed in this study may deepen our understanding of the interactions between the MV and the LV and help in risk stratification of heart valve disease and in silico treatment planning by exploring intrinsic compensation mechanisms.