In this paper,optimize-then-discretize,variational discretization and the finite volume method are applied to solve the distributed optimal control problems governed by a second order hyperbolic equation.A semi-discre...In this paper,optimize-then-discretize,variational discretization and the finite volume method are applied to solve the distributed optimal control problems governed by a second order hyperbolic equation.A semi-discrete optimal system is obtained.We prove the existence and uniqueness of the solution to the semidiscrete optimal system and obtain the optimal order error estimates in L ∞(J;L 2)-and L ∞(J;H 1)-norm.Numerical experiments are presented to test these theoretical results.展开更多
The maximum principle is a basic qualitative property of the solution of second-order elliptic boundary value problems.The preservation of the qualitative characteristics,such as the maximum principle,in discrete mode...The maximum principle is a basic qualitative property of the solution of second-order elliptic boundary value problems.The preservation of the qualitative characteristics,such as the maximum principle,in discrete model is one of the key requirements.It is well known that standard linear finite element solution does not satisfy maximum principle on general triangular meshes in 2D.In this paper we consider how to enforce discrete maximum principle for linear finite element solutions for the linear second-order self-adjoint elliptic equation.First approach is based on repair technique,which is a posteriori correction of the discrete solution.Second method is based on constrained optimization.Numerical tests that include anisotropic cases demonstrate how our method works for problems for which the standard finite element methods produce numerical solutions that violate the discrete maximum principle.展开更多
A two-dimensional nonlinear sloshing problem is analyzed by means of the fully nonlinear theory and time domain second order theory of water waves. Liquid sloshing in a rectangular container Subjected to a horizontal ...A two-dimensional nonlinear sloshing problem is analyzed by means of the fully nonlinear theory and time domain second order theory of water waves. Liquid sloshing in a rectangular container Subjected to a horizontal excitation is simulated by the finite element method. Comparisons between the two theories are made based on their numerical results. It is found that good agreement is obtained for the case of small amplitude oscillation and obvious differences occur for large amplitude excitation. Even though, the second order solution can still exhibit typical nonlinear features of nonlinear wave and can be used instead of the fully nonlinear theory.展开更多
Let(W,S) be a Coxeter group with S = I■J such that J consists of all universal elements of S and that I generates a finite parabolic subgroup W_I of W with w_0 the longest element of W_I. We describe all the left cel...Let(W,S) be a Coxeter group with S = I■J such that J consists of all universal elements of S and that I generates a finite parabolic subgroup W_I of W with w_0 the longest element of W_I. We describe all the left cells and two-sided cells of the weighted Coxeter group(W,S,L) that have non-empty intersection with W_J,where the weight function L of(W, S) is in one of the following cases:(i) max{L(s) | s ∈J} < min{L(t)|t∈I};(ii) min{L(s)|s ∈J} ≥L(w_0);(iii) there exists some t ∈ I satisfying L(t) < L(s) for any s ∈I-{t} and L takes a constant value L_J on J with L_J in some subintervals of [1, L(w_0)-1]. The results in the case(iii) are obtained under a certain assumption on(W, W_I).展开更多
In this paper,a residual type of a posteriori error estimator for the general second order elliptic eigenpair approximation by the mixed finite element method is derived and analyzed,based on a type of superconvergenc...In this paper,a residual type of a posteriori error estimator for the general second order elliptic eigenpair approximation by the mixed finite element method is derived and analyzed,based on a type of superconvergence result of the eigenfunction approximation.Its efficiency and reliability are proved by both theoretical analysis and numerical experiments.展开更多
he out of pile corrosion resistance of new Zirconium alloys in 500℃/10 3MPa steam has been investigated and the effect of alloying elements, Sn,Nb,Fe,Cr on the property has been analyzed. The results show that the...he out of pile corrosion resistance of new Zirconium alloys in 500℃/10 3MPa steam has been investigated and the effect of alloying elements, Sn,Nb,Fe,Cr on the property has been analyzed. The results show that the new alloys have better corrosion resistance than Zircaloy 4. That no nodular corrosion is found during test cycles shows the nodular corrosion resistance can be dramatically improved by addition of alloying elements Nb. The Fe/Cr ratio should be properly controlled if there is Cr addition in the alloys when developing new alloys.展开更多
In this paper,the micromorphic theory and the second gradient theory are proposedwhere the micromorphic model can be reduced to the second gradient model with the vanishing relative deformation between macrodeformatio...In this paper,the micromorphic theory and the second gradient theory are proposedwhere the micromorphic model can be reduced to the second gradient model with the vanishing relative deformation between macrodeformation gradient and microdeformation.Analytical solutions for the simple shear problem in the case of a general small strain isotropic elasticity micromorphic model and the second gradient model are presented,respectively.Besides,uniaxial tension of a constrained layer with two different boundary conditions is also analytically solved.Finally,the micromorphic theory is implemented numerically within a two-dimensional plane strain finite element framework by developing two isoparametric elements.展开更多
A novel multiscale algorithm based on the higher-order continuum at both micro-and macrostructural level is proposed for the consideration of the quasi-brittle damage response of heterogeneous materials.Herein,the mic...A novel multiscale algorithm based on the higher-order continuum at both micro-and macrostructural level is proposed for the consideration of the quasi-brittle damage response of heterogeneous materials.Herein,the microlevel damage is modelled by the degradation of the homogenized stress and tangent stiffness tensors,which are then upscaled to govern the localization at the macrolevel.The C^1 continuity finite element employing a modified case of Mindlin’s form II strain energy density is derived for the softening analysis.To the authors’knowledge,the finite element discretization based on the strain gradient theory is applied for the modeling of damage evolution at the microstructural level for heterogeneous materials for the first time.The advantage of the novel C1 finite element formulation in comparison with the standard finite element discretization in terms of the regularization efficiency as well as the objectivity has been shown.An isotropic damage law is used for the reduction of the constitutive and nonlocal material behaviour,which is necessary for the physically correct description of the localization formation in quasi-brittle materials.The capabilities of the derived finite element to capture the fully developed localization zones are tested on a random representative volume element(RVE)for several different loading cases.By employing the conventional second-order computational homogenization,the microstructural material constitutive response is averaged over the whole RVE area.In order to model the loss of structural integrity when sharp localization is formed across RVE,the specific conditions which detect a completely formed localization zone are developed.A new failure criterion at the microstructural level has been proposed.The derived finite element formulation,as well as the multiscale damage algorithm,are implemented into the finite element program ABAQUS.The capabilities of the presented multiscale scheme to capture the effects of the deformation localization are demonstrated by few ben展开更多
Ground-penetrating radar(GPR)is a highly efficient,fast and non-destructive exploration method for shallow surfaces.High-precision numerical simulation method is employed to improve the interpretation precision of det...Ground-penetrating radar(GPR)is a highly efficient,fast and non-destructive exploration method for shallow surfaces.High-precision numerical simulation method is employed to improve the interpretation precision of detection.Second-generation wavelet finite element is introduced into the forward modeling of the GPR.As the finite element basis function,the second-generation wavelet scaling function constructed by the scheme is characterized as having multiple scales and resolutions.The function can change the analytical scale arbitrarily according to actual needs.We can adopt a small analysis scale at a large gradient to improve the precision of analysis while adopting a large analytical scale at a small gradient to improve the efficiency of analysis.This approach is beneficial to capture the local mutation characteristics of the solution and improve the resolution without changing mesh subdivision to realize the efficient solution of the forward GPR problem.The algorithm is applied to the numerical simulation of line current radiation source and tunnel non-dense lining model with analytical solutions.Result show that the solution results of the secondgeneration wavelet finite element are in agreement with the analytical solutions and the conventional finite element solutions,thereby verifying the accuracy of the second-generation wavelet finite element algorithm.Furthermore,the second-generation wavelet finite element algorithm can change the analysis scale arbitrarily according to the actual problem without subdividing grids again.The adaptive algorithm is superior to traditional scheme in grid refinement and basis function order increase,which makes this algorithm suitable for solving complex GPR forward-modeling problems with large gradient and singularity.展开更多
Helmert’s second method of condensation is an effective method for terrain reduction in the geoid and quasi-geoid determinations. Condensing the masses outside the geoid to a surface layer on the geoid produces sever...Helmert’s second method of condensation is an effective method for terrain reduction in the geoid and quasi-geoid determinations. Condensing the masses outside the geoid to a surface layer on the geoid produces several forms of topographic effects: direct effect on gravity, secondary indirect effect on gravity and indirect effects on the (quasi-) geoid, respectively. To strike a balance between computation accuracy and numerical efficiency, the global integration region of topographic effects is usually divided into near zone and far zone. We focus on the computation of near-zone topographic effects, which are functions of actual topographic masses and condensed masses. Since there have already been mature formulas for gravitational attraction and potential of actual topographic masses using rectangular prism model, we put forward surface element model for condensed masses. Afterwards, the formulas for near-zone direct and indirect effects are obtained easily by combining the rectangular prism model and surface element model. To overcome the planar approximation errors involved with the new formulas for near-zone topographic effects, the Earth’s curvature can be taken into account. It is recommended to apply the formulas based on the rectangular prism and surface element considering the Earth’s curvature to calculate near-zone topographic effects for high-accuracy demand to determine geoid and quasi-geoid.展开更多
The paper deals with the numerical modelling of ductile damage responses in heterogeneous materials using the classical second-order homogenization approach.The scale transition methodology in the multiscale framework...The paper deals with the numerical modelling of ductile damage responses in heterogeneous materials using the classical second-order homogenization approach.The scale transition methodology in the multiscale framework is described.The structure at the macrolevel is discretized by the triangular C^(1) finite elements obeying nonlocal continuum theory,while the discretization of microstructural volume element at the microscale is conducted by means of the mixed type quadrilateral finite element with the nonlocal equivalent plastic strain as an additional nodal variable.The ductile damage evolution at the microlevel is modelled by using the gradient enhanced elastoplasticity.The macrolevel softening is governed by two criterions expressed by the increase in homogenized damage variable and the threshold of the local equivalent strain.The softening at each material point at the macrolevel is detected by the critical value of the homogenized damage,where homogenization of the damage variable is performed onlywithin softening area.Due to the nonlocal continuumtheory applied,a realistic softening behaviour is demonstrated after the damage initiation,compared to the widely used first-order homogenization approach.All algorithms derived have been embedded into the finite element code ABAQUS by means of the user subroutines and verified on the standard benchmark problems.The damage evolution at both microlevel and macrolevel has been demonstrated.展开更多
We perform a comparison in terms of accuracy and CPU time between second order BDF semi-Lagrangian and Lagrange-Galerkin schemes in combination with high order finite element method.The numerical results show that for...We perform a comparison in terms of accuracy and CPU time between second order BDF semi-Lagrangian and Lagrange-Galerkin schemes in combination with high order finite element method.The numerical results show that for polynomials of degree 2 semi-Lagrangian schemes are faster than Lagrange-Galerkin schemes for the same number of degrees of freedom,however,for the same level of accuracy both methods are about the same in terms of CPU time.For polynomials of degree larger than 2,Lagrange-Galerkin schemes behave better than semi-Lagrangian schemes in terms of both accuracy and CPU time;specially,for polynomials of degree 8 or larger.Also,we have performed tests on the parallelization of these schemes and the speedup obtained is quasi-optimal even with more than 100 processors.展开更多
A computational study on the enhancement of the second harmonic generation(SHG)in one-dimensional(1D)photonic crystals is presented.The mathematical model is derived from a nonlinear system of Maxwell’s equations,whi...A computational study on the enhancement of the second harmonic generation(SHG)in one-dimensional(1D)photonic crystals is presented.The mathematical model is derived from a nonlinear system of Maxwell’s equations,which partly overcomes the shortcoming of some existing models based on the undepleted pump approximation.We designed an iterative scheme coupled with the finite element method which can be applied to simulate the SHG in one dimensional nonlinear photonic band gap structures in our previous work.For the case that the nonlinearity is strong which is desirable to enhance the conversion efficiency,a continuation method is introduced to ensure the convergence of the iterative procedure.The convergence of our method is fast.Numerical experiments also indicate the conversion efficiency of SHG can be significantly enhanced when the frequencies of the fundamental and the second harmonic wave are tuned at the photonic band edges.The maximum total conversion efficiency available reaches more than 50%in all the cases studied.展开更多
We consider acoustic scattering of time-harmonic waves at objects composed of several homogeneous parts.Some of those may be impenetrable,giving rise to Dirichlet boundary conditions on their surfaces.We start from th...We consider acoustic scattering of time-harmonic waves at objects composed of several homogeneous parts.Some of those may be impenetrable,giving rise to Dirichlet boundary conditions on their surfaces.We start from the recent secondkind boundary integral approach of[X.Claeys,and R.Hiptmair,and E.Spindler.A second-kind Galerkin boundary element method for scattering at composite objects.BIT Numerical Mathematics,55(1):33-57,2015]for pure transmission problems and extend it to settings with essential boundary conditions.Based on so-called global multipotentials,we derive variational second-kind boundary integral equations posed in L^(2)(S),where S denotes the union of material interfaces.To suppress spurious resonances,we introduce a combined-field version(CFIE)of our new method.Thorough numerical tests highlight the low andmesh-independent condition numbers of Galerkin matrices obtained with discontinuous piecewise polynomial boundary element spaces.They also confirm competitive accuracy of the numerical solution in comparison with the widely used first-kind single-trace approach.展开更多
Dynamic equations of motional flexible beam elements were derived considering second-order effect. Non-linear finite element method and three-node Euler-Bernoulli beam elements were used. Because accuracy is higher in...Dynamic equations of motional flexible beam elements were derived considering second-order effect. Non-linear finite element method and three-node Euler-Bernoulli beam elements were used. Because accuracy is higher in non-linear structural analysis,three-node beam elements are used to deduce shape functions and stiffness matrices in dynamic equations of flexible elements. Static condensation method was used to obtain the finial dynamic equations of three-node beam elements. According to geometrical relations of nodal displacements in concomitant and global coordinate system,dynamic equations of elements can be transformed to global coordinate system by concomitant coordinate method in order to build the global dynamic equations. Analyzed amplitude condition of flexible arm support of a port crane,the results show that second-order effect should be considered in kinetic-elastic analysis for heavy load machinery of big flexibility.展开更多
基金supported by National Natural Science Foundation of China(Grant Nos.11261011,11271145 and 11031006)Foundation of Guizhou Science and Technology Department(Grant No.[2011]2098)+2 种基金Foundation for Talent Introduction of Guangdong Provincial UniversitySpecialized Research Fund for the Doctoral Program of Higher Education(Grant No. 20114407110009)the Project of Department of Education of Guangdong Province(Grant No. 2012KJCX0036)
文摘In this paper,optimize-then-discretize,variational discretization and the finite volume method are applied to solve the distributed optimal control problems governed by a second order hyperbolic equation.A semi-discrete optimal system is obtained.We prove the existence and uniqueness of the solution to the semidiscrete optimal system and obtain the optimal order error estimates in L ∞(J;L 2)-and L ∞(J;H 1)-norm.Numerical experiments are presented to test these theoretical results.
基金the National Nuclear Security Administration of the U.S.Department of Energy at Los Alamos National Laboratory under Contract No.DE-AC52-06NA25396the DOE Office of Science Advanced Scientific Computing Research(ASCR)Program in Applied Mathematics Research.The first author has been supported in part by the Czech Ministry of Education projects MSM 6840770022 and LC06052(Necas Center for Mathematical Modeling).
文摘The maximum principle is a basic qualitative property of the solution of second-order elliptic boundary value problems.The preservation of the qualitative characteristics,such as the maximum principle,in discrete model is one of the key requirements.It is well known that standard linear finite element solution does not satisfy maximum principle on general triangular meshes in 2D.In this paper we consider how to enforce discrete maximum principle for linear finite element solutions for the linear second-order self-adjoint elliptic equation.First approach is based on repair technique,which is a posteriori correction of the discrete solution.Second method is based on constrained optimization.Numerical tests that include anisotropic cases demonstrate how our method works for problems for which the standard finite element methods produce numerical solutions that violate the discrete maximum principle.
文摘A two-dimensional nonlinear sloshing problem is analyzed by means of the fully nonlinear theory and time domain second order theory of water waves. Liquid sloshing in a rectangular container Subjected to a horizontal excitation is simulated by the finite element method. Comparisons between the two theories are made based on their numerical results. It is found that good agreement is obtained for the case of small amplitude oscillation and obvious differences occur for large amplitude excitation. Even though, the second order solution can still exhibit typical nonlinear features of nonlinear wave and can be used instead of the fully nonlinear theory.
基金supported by National Natural Science Foundation of China (Grant Nos. 11131001 and 11471115)Shanghai Key Laboratory of Pure Mathematics and Mathematical PracticeScience and Technology Commission of Shanghai Municipality (Grant No.13dz2260400)
文摘Let(W,S) be a Coxeter group with S = I■J such that J consists of all universal elements of S and that I generates a finite parabolic subgroup W_I of W with w_0 the longest element of W_I. We describe all the left cells and two-sided cells of the weighted Coxeter group(W,S,L) that have non-empty intersection with W_J,where the weight function L of(W, S) is in one of the following cases:(i) max{L(s) | s ∈J} < min{L(t)|t∈I};(ii) min{L(s)|s ∈J} ≥L(w_0);(iii) there exists some t ∈ I satisfying L(t) < L(s) for any s ∈I-{t} and L takes a constant value L_J on J with L_J in some subintervals of [1, L(w_0)-1]. The results in the case(iii) are obtained under a certain assumption on(W, W_I).
基金supported by National Natural Science Foundation of China(Grant Nos.11001259,11031006,11071265,11201501 and 91230110)National Basic Research Program of China(973 Project)(Grant No. 2011CB309703)+3 种基金International S&T Cooperation Program of China(Grant No. 2010DFR00700)Croucher Foundation of Hong Kong Baptist Universitythe National Center for Mathematics and Interdisciplinary Science,CAS,the President Foundation of AMSS-CASthe Fundamental Research Funds for the CentralUniversities(Grant No. 2012121003)
文摘In this paper,a residual type of a posteriori error estimator for the general second order elliptic eigenpair approximation by the mixed finite element method is derived and analyzed,based on a type of superconvergence result of the eigenfunction approximation.Its efficiency and reliability are proved by both theoretical analysis and numerical experiments.
文摘he out of pile corrosion resistance of new Zirconium alloys in 500℃/10 3MPa steam has been investigated and the effect of alloying elements, Sn,Nb,Fe,Cr on the property has been analyzed. The results show that the new alloys have better corrosion resistance than Zircaloy 4. That no nodular corrosion is found during test cycles shows the nodular corrosion resistance can be dramatically improved by addition of alloying elements Nb. The Fe/Cr ratio should be properly controlled if there is Cr addition in the alloys when developing new alloys.
基金supported by the National Natural Science Foundation of China (No. 10772096)
文摘In this paper,the micromorphic theory and the second gradient theory are proposedwhere the micromorphic model can be reduced to the second gradient model with the vanishing relative deformation between macrodeformation gradient and microdeformation.Analytical solutions for the simple shear problem in the case of a general small strain isotropic elasticity micromorphic model and the second gradient model are presented,respectively.Besides,uniaxial tension of a constrained layer with two different boundary conditions is also analytically solved.Finally,the micromorphic theory is implemented numerically within a two-dimensional plane strain finite element framework by developing two isoparametric elements.
基金This work has been fully supported by Croatian Science Foundation under the project“Multiscale Numerical Modelling of Material Deformation Responses from Macro-to Nanolevel”(2516).
文摘A novel multiscale algorithm based on the higher-order continuum at both micro-and macrostructural level is proposed for the consideration of the quasi-brittle damage response of heterogeneous materials.Herein,the microlevel damage is modelled by the degradation of the homogenized stress and tangent stiffness tensors,which are then upscaled to govern the localization at the macrolevel.The C^1 continuity finite element employing a modified case of Mindlin’s form II strain energy density is derived for the softening analysis.To the authors’knowledge,the finite element discretization based on the strain gradient theory is applied for the modeling of damage evolution at the microstructural level for heterogeneous materials for the first time.The advantage of the novel C1 finite element formulation in comparison with the standard finite element discretization in terms of the regularization efficiency as well as the objectivity has been shown.An isotropic damage law is used for the reduction of the constitutive and nonlocal material behaviour,which is necessary for the physically correct description of the localization formation in quasi-brittle materials.The capabilities of the derived finite element to capture the fully developed localization zones are tested on a random representative volume element(RVE)for several different loading cases.By employing the conventional second-order computational homogenization,the microstructural material constitutive response is averaged over the whole RVE area.In order to model the loss of structural integrity when sharp localization is formed across RVE,the specific conditions which detect a completely formed localization zone are developed.A new failure criterion at the microstructural level has been proposed.The derived finite element formulation,as well as the multiscale damage algorithm,are implemented into the finite element program ABAQUS.The capabilities of the presented multiscale scheme to capture the effects of the deformation localization are demonstrated by few ben
基金supported by the National Natural Science Foundation of China(Nos.41574116 and 41774132)Hunan Provincial Innovation Foundation for Postgraduate(Grant Nos.CX2017B052)the Fundamental Research Funds for the Central Universities of Central South University(Nos.2018zzts693)。
文摘Ground-penetrating radar(GPR)is a highly efficient,fast and non-destructive exploration method for shallow surfaces.High-precision numerical simulation method is employed to improve the interpretation precision of detection.Second-generation wavelet finite element is introduced into the forward modeling of the GPR.As the finite element basis function,the second-generation wavelet scaling function constructed by the scheme is characterized as having multiple scales and resolutions.The function can change the analytical scale arbitrarily according to actual needs.We can adopt a small analysis scale at a large gradient to improve the precision of analysis while adopting a large analytical scale at a small gradient to improve the efficiency of analysis.This approach is beneficial to capture the local mutation characteristics of the solution and improve the resolution without changing mesh subdivision to realize the efficient solution of the forward GPR problem.The algorithm is applied to the numerical simulation of line current radiation source and tunnel non-dense lining model with analytical solutions.Result show that the solution results of the secondgeneration wavelet finite element are in agreement with the analytical solutions and the conventional finite element solutions,thereby verifying the accuracy of the second-generation wavelet finite element algorithm.Furthermore,the second-generation wavelet finite element algorithm can change the analysis scale arbitrarily according to the actual problem without subdividing grids again.The adaptive algorithm is superior to traditional scheme in grid refinement and basis function order increase,which makes this algorithm suitable for solving complex GPR forward-modeling problems with large gradient and singularity.
基金The National Natural Science Foundation of China (41674025,41674082)The Independent Research Foundation of State Key Laboratory of Geo-information Engineering (SKLGIE2018-ZZ-10).
文摘Helmert’s second method of condensation is an effective method for terrain reduction in the geoid and quasi-geoid determinations. Condensing the masses outside the geoid to a surface layer on the geoid produces several forms of topographic effects: direct effect on gravity, secondary indirect effect on gravity and indirect effects on the (quasi-) geoid, respectively. To strike a balance between computation accuracy and numerical efficiency, the global integration region of topographic effects is usually divided into near zone and far zone. We focus on the computation of near-zone topographic effects, which are functions of actual topographic masses and condensed masses. Since there have already been mature formulas for gravitational attraction and potential of actual topographic masses using rectangular prism model, we put forward surface element model for condensed masses. Afterwards, the formulas for near-zone direct and indirect effects are obtained easily by combining the rectangular prism model and surface element model. To overcome the planar approximation errors involved with the new formulas for near-zone topographic effects, the Earth’s curvature can be taken into account. It is recommended to apply the formulas based on the rectangular prism and surface element considering the Earth’s curvature to calculate near-zone topographic effects for high-accuracy demand to determine geoid and quasi-geoid.
文摘The paper deals with the numerical modelling of ductile damage responses in heterogeneous materials using the classical second-order homogenization approach.The scale transition methodology in the multiscale framework is described.The structure at the macrolevel is discretized by the triangular C^(1) finite elements obeying nonlocal continuum theory,while the discretization of microstructural volume element at the microscale is conducted by means of the mixed type quadrilateral finite element with the nonlocal equivalent plastic strain as an additional nodal variable.The ductile damage evolution at the microlevel is modelled by using the gradient enhanced elastoplasticity.The macrolevel softening is governed by two criterions expressed by the increase in homogenized damage variable and the threshold of the local equivalent strain.The softening at each material point at the macrolevel is detected by the critical value of the homogenized damage,where homogenization of the damage variable is performed onlywithin softening area.Due to the nonlocal continuumtheory applied,a realistic softening behaviour is demonstrated after the damage initiation,compared to the widely used first-order homogenization approach.All algorithms derived have been embedded into the finite element code ABAQUS by means of the user subroutines and verified on the standard benchmark problems.The damage evolution at both microlevel and macrolevel has been demonstrated.
基金funded by grant CGL2007-66440-C04-01 from Ministerio de Educacion y Ciencia de Espana.
文摘We perform a comparison in terms of accuracy and CPU time between second order BDF semi-Lagrangian and Lagrange-Galerkin schemes in combination with high order finite element method.The numerical results show that for polynomials of degree 2 semi-Lagrangian schemes are faster than Lagrange-Galerkin schemes for the same number of degrees of freedom,however,for the same level of accuracy both methods are about the same in terms of CPU time.For polynomials of degree larger than 2,Lagrange-Galerkin schemes behave better than semi-Lagrangian schemes in terms of both accuracy and CPU time;specially,for polynomials of degree 8 or larger.Also,we have performed tests on the parallelization of these schemes and the speedup obtained is quasi-optimal even with more than 100 processors.
基金supported by the fundamental research funds for the central universities(BUPT2009RC0706)the Project 11001030 supported by National Natural Science Foundation of China+3 种基金the open fund of key laboratory of Information Photonics and Optical Communications(Beijing University of Posts and Telecommunications)of Ministry of Educationsupported in part by the NSF grants DMS-0604790,DMS-0908325,CCF-0830161,EAR-0724527,and DMS-0968360the ONR grant N00014-09-1-0384a special research grant from Zhejiang University.
文摘A computational study on the enhancement of the second harmonic generation(SHG)in one-dimensional(1D)photonic crystals is presented.The mathematical model is derived from a nonlinear system of Maxwell’s equations,which partly overcomes the shortcoming of some existing models based on the undepleted pump approximation.We designed an iterative scheme coupled with the finite element method which can be applied to simulate the SHG in one dimensional nonlinear photonic band gap structures in our previous work.For the case that the nonlinearity is strong which is desirable to enhance the conversion efficiency,a continuation method is introduced to ensure the convergence of the iterative procedure.The convergence of our method is fast.Numerical experiments also indicate the conversion efficiency of SHG can be significantly enhanced when the frequencies of the fundamental and the second harmonic wave are tuned at the photonic band edges.The maximum total conversion efficiency available reaches more than 50%in all the cases studied.
基金The authors would like to thank L.Kielhorn for his great support during the development of the code for the first-and second-kind formulation in BETL2[25]The work of E.Spindler was partially supported by SNF under grant 20021137873/1X.Claeys received support from the ANR Research Grant ANR-15-CE23-0017-01.
文摘We consider acoustic scattering of time-harmonic waves at objects composed of several homogeneous parts.Some of those may be impenetrable,giving rise to Dirichlet boundary conditions on their surfaces.We start from the recent secondkind boundary integral approach of[X.Claeys,and R.Hiptmair,and E.Spindler.A second-kind Galerkin boundary element method for scattering at composite objects.BIT Numerical Mathematics,55(1):33-57,2015]for pure transmission problems and extend it to settings with essential boundary conditions.Based on so-called global multipotentials,we derive variational second-kind boundary integral equations posed in L^(2)(S),where S denotes the union of material interfaces.To suppress spurious resonances,we introduce a combined-field version(CFIE)of our new method.Thorough numerical tests highlight the low andmesh-independent condition numbers of Galerkin matrices obtained with discontinuous piecewise polynomial boundary element spaces.They also confirm competitive accuracy of the numerical solution in comparison with the widely used first-kind single-trace approach.
文摘Dynamic equations of motional flexible beam elements were derived considering second-order effect. Non-linear finite element method and three-node Euler-Bernoulli beam elements were used. Because accuracy is higher in non-linear structural analysis,three-node beam elements are used to deduce shape functions and stiffness matrices in dynamic equations of flexible elements. Static condensation method was used to obtain the finial dynamic equations of three-node beam elements. According to geometrical relations of nodal displacements in concomitant and global coordinate system,dynamic equations of elements can be transformed to global coordinate system by concomitant coordinate method in order to build the global dynamic equations. Analyzed amplitude condition of flexible arm support of a port crane,the results show that second-order effect should be considered in kinetic-elastic analysis for heavy load machinery of big flexibility.
