For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geomet...For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.展开更多
In this paper the conservation theorems of the constrained Birkhoffian systems are studied by using the method of integrating factors. The differential equations of motion of the system are written. The definition of ...In this paper the conservation theorems of the constrained Birkhoffian systems are studied by using the method of integrating factors. The differential equations of motion of the system are written. The definition of integrating factors is given for the system. The necessary conditions for the existence of the conserved quantity for the system are studied. The conservation theorem and its inverse for the system are established. Finally, an example is given to illustrate the application of the results.展开更多
The compact implicit integration factor (cIIF) method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF me...The compact implicit integration factor (cIIF) method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF method to some complex-valued nonlinear evolutionary equations such as the nonlinear SchrSdinger (NLS) equation and the complex Ginzburg-Landau (GL) equation. Detailed algorithm formulation and practical implementation of cIIF method are performed. The numerical results indicate that this method is very accurate and efficient.展开更多
Let (P, Q) be a C1 vector field defined in an open subset U IR2. We call inverse integrating factor a C1 solution V(x, y) of the equation . In previous works it has been shown that this function plays an important ro...Let (P, Q) be a C1 vector field defined in an open subset U IR2. We call inverse integrating factor a C1 solution V(x, y) of the equation . In previous works it has been shown that this function plays an important role in the problem of the center and in the determination of limit cycles. In this paper we obtain necessary conditions for a polynomial vector field (P, Q) to have a polynomial inverse integrating factor.展开更多
State functions play important roles in thermodynamics.Different from the process function,such as the exchanged heatδQ and the applied workδW,the change of the state function can be expressed as an exact differenti...State functions play important roles in thermodynamics.Different from the process function,such as the exchanged heatδQ and the applied workδW,the change of the state function can be expressed as an exact differential.We prove here that,for a generic thermodynamic system,only the inverse of the temperature,namely 1/T,can serve as the integration factor for the exchanged heatδQ.The uniqueness of the integration factor invalidates any attempt to define other state functions associated with the exchanged heat,and in turn,reveals the incorrectness of defining the entransy E_(vh)=CVT^(2)/2 as a state function by treating T as an integration factor.We further show the errors in the derivation of entransy by treating the heat capacity C_(V)as a temperature-independent constant.展开更多
Current Situation and Problems of the Treatment in Advanced Prostate Cancer In recent years,the incidence of prostate cancer shows a rising trend in China with an increase of 70%and has been the first place in the gro...Current Situation and Problems of the Treatment in Advanced Prostate Cancer In recent years,the incidence of prostate cancer shows a rising trend in China with an increase of 70%and has been the first place in the growth rate of malignant tumor in the male reproductive system. Prostate cancer has become a serious threat to male senior’s health.Because of the application of展开更多
A high-order accuracy time discretization method is developed in this paper to solve the one-dimensional nonlinear Dirac(NLD)equation.Based on the implicit integration factor(IIF)method,two schemes are proposed.Centra...A high-order accuracy time discretization method is developed in this paper to solve the one-dimensional nonlinear Dirac(NLD)equation.Based on the implicit integration factor(IIF)method,two schemes are proposed.Central differences are applied to the spatial discretization.The semi-discrete scheme keeps the conservation of the charge and energy.For the temporal discretization,second-order IIF method and fourth-order IIF method are applied respectively to the nonlinear system arising from the spatial discretization.Numerical experiments are given to validate the accuracy of these schemes and to discuss the interaction dynamics of the NLD solitary waves.展开更多
Weighted interior penalty discontinuous Galerkin method is developed to solve the two-dimensional non-equilibrium radiation diffusion equation on unstructured mesh.There are three weights including the arithmetic,the ...Weighted interior penalty discontinuous Galerkin method is developed to solve the two-dimensional non-equilibrium radiation diffusion equation on unstructured mesh.There are three weights including the arithmetic,the harmonic,and the geometric weight in the weighted discontinuous Galerkin scheme.For the time discretization,we treat the nonlinear diffusion coefficients explicitly,and apply the semiimplicit integration factormethod to the nonlinear ordinary differential equations arising from discontinuous Galerkin spatial discretization.The semi-implicit integration factor method can not only avoid severe timestep limits,but also takes advantage of the local property of DG methods by which small sized nonlinear algebraic systems are solved element by element with the exact Newton iteration method.Numerical results are presented to demonstrate the validity of discontinuous Galerkin method for high nonlinear and tightly coupled radiation diffusion equation.展开更多
Objective To investigate the integral dissolution model based on biological potency in order to evaluate the dissolution of Compound Chinese materia medica(CCMM) in vitro. Methods The contents of paeoniflorin, phill...Objective To investigate the integral dissolution model based on biological potency in order to evaluate the dissolution of Compound Chinese materia medica(CCMM) in vitro. Methods The contents of paeoniflorin, phillyrin, ginsenoside Rg1, and adenosine of ten batches of Compound Biejia Ruangan Tablet(CBRT) were determined at different times. The self-defined weighting coefficient based on the contents has been created to establish the integral dissolution model. In addition, the biological potency of CBRT was measured by MTT assay. Then, the f2 similar factor was used to evaluate the similarity of the batches. Results Compared with batch a, some batches’ f2 values of paeoniflorin and adenosine were less than 50, while f2 values of ginsenoside Rg1, phillyrin, and integral component were more than 50. Likewise, ginsenoside Rg1, phillyrin, and integral component were all in good correlation with biological dissolution. Conclusion The results of the integral dissolution based on biological test of CBRT demonstrate that the bioassay method may be a promising supplement for its quality evaluation.展开更多
For the planar Z2-equivariant cubic systems having twoelementary focuses, the characterization of a bi-center problem and shortened expressions of the first six Liapunov constants are completely discussed. The necessa...For the planar Z2-equivariant cubic systems having twoelementary focuses, the characterization of a bi-center problem and shortened expressions of the first six Liapunov constants are completely discussed. The necessary and sufficient conditions for the existence of the bi-center are obtained. All possible first integrals are given. Under small Z2-equivariant cubic perturbations, the conclusion that there exist at most 12 small-amplitude limit cycles with the scheme (6 II 6) is proved.展开更多
In this study, we present a conservative local discontinuous Galerkin(LDG) method for numerically solving the two-dimensional nonlinear Schrdinger(NLS) equation. The NLS equation is rewritten as a firstorder system an...In this study, we present a conservative local discontinuous Galerkin(LDG) method for numerically solving the two-dimensional nonlinear Schrdinger(NLS) equation. The NLS equation is rewritten as a firstorder system and then we construct the LDG formulation with appropriate numerical flux. The mass and energy conserving laws for the semi-discrete formulation can be proved based on different choices of numerical fluxes such as the central, alternative and upwind-based flux. We will propose two kinds of time discretization methods for the semi-discrete formulation. One is based on Crank-Nicolson method and can be proved to preserve the discrete mass and energy conservation. The other one is Krylov implicit integration factor(IIF) method which demands much less computational effort. Various numerical experiments are presented to demonstrate the conservation law of mass and energy, the optimal rates of convergence, and the blow-up phenomenon.展开更多
Private capital is one of the main driving forces in China's initiatives towards stimulating the market economy. The development of private economy in China has always been based on integrating industrial and corp...Private capital is one of the main driving forces in China's initiatives towards stimulating the market economy. The development of private economy in China has always been based on integrating industrial and corporate structures with product composition and market structures. This paper explores the development of the private economy and how it integrates different industries with specific markets by analyzing the leading private sector in Zhejiang province. It also examines the trends of industrial cluster, the formation of the agglomerative economy and their effects on private economy development. Finally, the paper explains why Zhejiang people have profited much from the Wenzhou model and discusses some existing problems and future possibilities for development of the Wenzhou model.展开更多
The effects of natural andgeochemical factors depending on heavy metal contamination in nuisance dust particles were evaluated. The nuisance dust particles were sampled using passive deposit gauge method for one year ...The effects of natural andgeochemical factors depending on heavy metal contamination in nuisance dust particles were evaluated. The nuisance dust particles were sampled using passive deposit gauge method for one year from April2010 to March2011 and the obtained samples were measured for the total contents and the contamination levels of Fe, Mn, Cu and As usinggeo-accumulation index (Igeo ), enrichment factor (EF) and the integrated pollution index (IPI). The results showed that, the contamination levels of Fe and Mn based on Igeo values, were uncontaminated (Igeo 〈 0) (variations of the Igeo index was from -3.11 to -1.751 for Fe, from -0.630 to -1.925 for Mn), while the values of Cu and As were demonstrated to have moderate contamination based on Igeo values (variations of Igeo index was from -1.125 to 0.848 for Cu, and from -2.002 to 1.249 for As). The analysis of EF also revealed minor to moderate enrichment for Mn (1.215-4.214), minor to moderately severe enrichment for Cu (2.791-6.484), and As (1.370-8.462), respectively. The variation of the IPI index also showed low to moderate level of heavy metal pollution in nuisance dust particulates (0.511-1.829). The analysis of the results also approved that the natural processes andgeochemical variables (the changing meteorological parameters) can significantly affect the availability of heavy metals in nuisance dust particles in Western Iran.展开更多
We investigate and concentrate on new infinitesimal generators of Lie symmetries for an extended(2+1)-dimensional Calogero-Bogoyavlenskii-Schif(eCBS)equation using the commutator table which results in a system of non...We investigate and concentrate on new infinitesimal generators of Lie symmetries for an extended(2+1)-dimensional Calogero-Bogoyavlenskii-Schif(eCBS)equation using the commutator table which results in a system of nonlinear ordinary differential equations(ODEs)which can be manually solved.Through two stages of Lie symmetry reductions,the eCBS equation is reduced to non-solvable nonlinear ODEs using different combinations of optimal Lie vectors.Using the integration method and the Riccati and Bernoulli equation methods,we investigate new analytical solutions to those ODEs.Back substituting to the original variables generates new solutions to the eCBS equation.These results are simulated through three-and two-dimensional plots.展开更多
This paper proposes and analyzes an efficient finite difference scheme for the two-dimensional nonlinear Schr?dinger(NLS) equation involving fractional Laplacian. The scheme is based on a weighted and shifted Grü...This paper proposes and analyzes an efficient finite difference scheme for the two-dimensional nonlinear Schr?dinger(NLS) equation involving fractional Laplacian. The scheme is based on a weighted and shifted Grünwald-Letnikov difference(WSGD) operator for the spatial fractional Laplacian. We prove that the proposed method preserves the mass and energy conservation laws in semi-discrete formulations. By introducing the differentiation matrices, the semi-discrete fractional nonlinear Schr?dinger(FNLS) equation can be rewritten as a system of nonlinear ordinary differential equations(ODEs) in matrix formulations. Two kinds of time discretization methods are proposed for the semi-discrete formulation. One is based on the Crank-Nicolson(CN) method which can be proved to preserve the fully discrete mass and energy conservation. The other one is the compact implicit integration factor(c IIF) method which demands much less computational effort. It can be shown that the cIIF scheme can approximate CN scheme with the error O(τ~2). Finally numerical results are presented to demonstrate the method’s conservation, accuracy, efficiency and the capability of capturing blow-up.展开更多
We characterize the complex differential equations of the form dy/dx=a_(n)(x)y^)n_+a_(n-1)(x)y^(n-1)+…+a_(1)(x)y+a_(0)(x) where a_(j)(x) are meromorphic functions in the variable x for j = 0,..., n that admit either ...We characterize the complex differential equations of the form dy/dx=a_(n)(x)y^)n_+a_(n-1)(x)y^(n-1)+…+a_(1)(x)y+a_(0)(x) where a_(j)(x) are meromorphic functions in the variable x for j = 0,..., n that admit either a Weierstrass first integral or a Weierstrass inverse integrating factor.展开更多
In this paper, we study the behavior of the solution at the crack edges for a nearly circular crack with developing cusps subject to shear loading. The problem of finding the resulting force can be written in the form...In this paper, we study the behavior of the solution at the crack edges for a nearly circular crack with developing cusps subject to shear loading. The problem of finding the resulting force can be written in the form of a hypersingular integral equation. The equation is then trans-formed into a similar equation over a circular region using conformal mapping. The equation is solved numerically for the unknown coefficients, which will later be used in finding the stress intensity factors. The sliding and tearing mode stress intensity factors are evaluated for cracks and displayed graphically. Our results seem to agree with the existing asymptotic solution.展开更多
文摘For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.
基金Project supported by the Heilongjiang Natural Science Foundation of China (Grant No 9507)
文摘In this paper the conservation theorems of the constrained Birkhoffian systems are studied by using the method of integrating factors. The differential equations of motion of the system are written. The definition of integrating factors is given for the system. The necessary conditions for the existence of the conserved quantity for the system are studied. The conservation theorem and its inverse for the system are established. Finally, an example is given to illustrate the application of the results.
文摘The compact implicit integration factor (cIIF) method is an efficient time discretization scheme for stiff nonlinear diffusion equations in two and three spatial dimensions. In the current work, we apply the cIIF method to some complex-valued nonlinear evolutionary equations such as the nonlinear SchrSdinger (NLS) equation and the complex Ginzburg-Landau (GL) equation. Detailed algorithm formulation and practical implementation of cIIF method are performed. The numerical results indicate that this method is very accurate and efficient.
基金the DGICYT grant, number PB96-1153 The third author is partially supported by the University of Lleida Project P98-207.
文摘Let (P, Q) be a C1 vector field defined in an open subset U IR2. We call inverse integrating factor a C1 solution V(x, y) of the equation . In previous works it has been shown that this function plays an important role in the problem of the center and in the determination of limit cycles. In this paper we obtain necessary conditions for a polynomial vector field (P, Q) to have a polynomial inverse integrating factor.
基金This work was supported by the National Natural Science Foundation of China(NSFC)(Grants No.11534002,No.12088101,No.U1530402,No.U1930403,No.11775001,No.11534002,and No.11825001)the National Basic Research Program of China(Grants No.2016YFA0301201).
文摘State functions play important roles in thermodynamics.Different from the process function,such as the exchanged heatδQ and the applied workδW,the change of the state function can be expressed as an exact differential.We prove here that,for a generic thermodynamic system,only the inverse of the temperature,namely 1/T,can serve as the integration factor for the exchanged heatδQ.The uniqueness of the integration factor invalidates any attempt to define other state functions associated with the exchanged heat,and in turn,reveals the incorrectness of defining the entransy E_(vh)=CVT^(2)/2 as a state function by treating T as an integration factor.We further show the errors in the derivation of entransy by treating the heat capacity C_(V)as a temperature-independent constant.
基金Supported by the National Natural Science Foundation of China(No.30873268)
文摘Current Situation and Problems of the Treatment in Advanced Prostate Cancer In recent years,the incidence of prostate cancer shows a rising trend in China with an increase of 70%and has been the first place in the growth rate of malignant tumor in the male reproductive system. Prostate cancer has become a serious threat to male senior’s health.Because of the application of
基金the National Natural Science Foundation of China(No.11671044)the Science Challenge Project(No.TZ2016001)the Beijing Municipal Education Commission(No.PXM2017014224000020).
文摘A high-order accuracy time discretization method is developed in this paper to solve the one-dimensional nonlinear Dirac(NLD)equation.Based on the implicit integration factor(IIF)method,two schemes are proposed.Central differences are applied to the spatial discretization.The semi-discrete scheme keeps the conservation of the charge and energy.For the temporal discretization,second-order IIF method and fourth-order IIF method are applied respectively to the nonlinear system arising from the spatial discretization.Numerical experiments are given to validate the accuracy of these schemes and to discuss the interaction dynamics of the NLD solitary waves.
基金the National Nature Science Foundation of China(11171038)R.Zhang’s work was also supported by Brazilian Young Talent Attraction Program via National Council for Scientific and Technological Development(CNPq).J.Zhu and A.Loula’s works were partially supported by CNPq.X.Cui’s work was partially supported by the National Natural Science Foundation of China(11271054)+1 种基金the Science Foundation of CAEP(2010A0202010,2012B0202026)the Defense Industrial Technology Development Program(B1520110011).
文摘Weighted interior penalty discontinuous Galerkin method is developed to solve the two-dimensional non-equilibrium radiation diffusion equation on unstructured mesh.There are three weights including the arithmetic,the harmonic,and the geometric weight in the weighted discontinuous Galerkin scheme.For the time discretization,we treat the nonlinear diffusion coefficients explicitly,and apply the semiimplicit integration factormethod to the nonlinear ordinary differential equations arising from discontinuous Galerkin spatial discretization.The semi-implicit integration factor method can not only avoid severe timestep limits,but also takes advantage of the local property of DG methods by which small sized nonlinear algebraic systems are solved element by element with the exact Newton iteration method.Numerical results are presented to demonstrate the validity of discontinuous Galerkin method for high nonlinear and tightly coupled radiation diffusion equation.
基金Major Scientific and Technological Specialized Project for Significant New.Formulation of New Drugs(No.2011ZX09201-201-14)National Natural Science Foundation of China(No.81073069)
文摘Objective To investigate the integral dissolution model based on biological potency in order to evaluate the dissolution of Compound Chinese materia medica(CCMM) in vitro. Methods The contents of paeoniflorin, phillyrin, ginsenoside Rg1, and adenosine of ten batches of Compound Biejia Ruangan Tablet(CBRT) were determined at different times. The self-defined weighting coefficient based on the contents has been created to establish the integral dissolution model. In addition, the biological potency of CBRT was measured by MTT assay. Then, the f2 similar factor was used to evaluate the similarity of the batches. Results Compared with batch a, some batches’ f2 values of paeoniflorin and adenosine were less than 50, while f2 values of ginsenoside Rg1, phillyrin, and integral component were more than 50. Likewise, ginsenoside Rg1, phillyrin, and integral component were all in good correlation with biological dissolution. Conclusion The results of the integral dissolution based on biological test of CBRT demonstrate that the bioassay method may be a promising supplement for its quality evaluation.
基金Supported by National Natural Science Foundation of China (Grant Nos. 10831003 and 10771196)
文摘For the planar Z2-equivariant cubic systems having twoelementary focuses, the characterization of a bi-center problem and shortened expressions of the first six Liapunov constants are completely discussed. The necessary and sufficient conditions for the existence of the bi-center are obtained. All possible first integrals are given. Under small Z2-equivariant cubic perturbations, the conclusion that there exist at most 12 small-amplitude limit cycles with the scheme (6 II 6) is proved.
基金supported by the Foundation of Liaoning Educational Committee (Grant No. L201604)China Scholarship Council, National Natural Science Foundation of China (Grant Nos. 11571002, 11171281 and 11671044)+1 种基金the Science Foundation of China Academy of Engineering Physics (Grant No. 2015B0101021)the Defense Industrial Technology Development Program (Grant No. B1520133015)
文摘In this study, we present a conservative local discontinuous Galerkin(LDG) method for numerically solving the two-dimensional nonlinear Schrdinger(NLS) equation. The NLS equation is rewritten as a firstorder system and then we construct the LDG formulation with appropriate numerical flux. The mass and energy conserving laws for the semi-discrete formulation can be proved based on different choices of numerical fluxes such as the central, alternative and upwind-based flux. We will propose two kinds of time discretization methods for the semi-discrete formulation. One is based on Crank-Nicolson method and can be proved to preserve the discrete mass and energy conservation. The other one is Krylov implicit integration factor(IIF) method which demands much less computational effort. Various numerical experiments are presented to demonstrate the conservation law of mass and energy, the optimal rates of convergence, and the blow-up phenomenon.
文摘Private capital is one of the main driving forces in China's initiatives towards stimulating the market economy. The development of private economy in China has always been based on integrating industrial and corporate structures with product composition and market structures. This paper explores the development of the private economy and how it integrates different industries with specific markets by analyzing the leading private sector in Zhejiang province. It also examines the trends of industrial cluster, the formation of the agglomerative economy and their effects on private economy development. Finally, the paper explains why Zhejiang people have profited much from the Wenzhou model and discusses some existing problems and future possibilities for development of the Wenzhou model.
文摘The effects of natural andgeochemical factors depending on heavy metal contamination in nuisance dust particles were evaluated. The nuisance dust particles were sampled using passive deposit gauge method for one year from April2010 to March2011 and the obtained samples were measured for the total contents and the contamination levels of Fe, Mn, Cu and As usinggeo-accumulation index (Igeo ), enrichment factor (EF) and the integrated pollution index (IPI). The results showed that, the contamination levels of Fe and Mn based on Igeo values, were uncontaminated (Igeo 〈 0) (variations of the Igeo index was from -3.11 to -1.751 for Fe, from -0.630 to -1.925 for Mn), while the values of Cu and As were demonstrated to have moderate contamination based on Igeo values (variations of Igeo index was from -1.125 to 0.848 for Cu, and from -2.002 to 1.249 for As). The analysis of EF also revealed minor to moderate enrichment for Mn (1.215-4.214), minor to moderately severe enrichment for Cu (2.791-6.484), and As (1.370-8.462), respectively. The variation of the IPI index also showed low to moderate level of heavy metal pollution in nuisance dust particulates (0.511-1.829). The analysis of the results also approved that the natural processes andgeochemical variables (the changing meteorological parameters) can significantly affect the availability of heavy metals in nuisance dust particles in Western Iran.
文摘We investigate and concentrate on new infinitesimal generators of Lie symmetries for an extended(2+1)-dimensional Calogero-Bogoyavlenskii-Schif(eCBS)equation using the commutator table which results in a system of nonlinear ordinary differential equations(ODEs)which can be manually solved.Through two stages of Lie symmetry reductions,the eCBS equation is reduced to non-solvable nonlinear ODEs using different combinations of optimal Lie vectors.Using the integration method and the Riccati and Bernoulli equation methods,we investigate new analytical solutions to those ODEs.Back substituting to the original variables generates new solutions to the eCBS equation.These results are simulated through three-and two-dimensional plots.
基金supported by National Natural Science Foundation of China(Grant Nos.61573008 and 61703290)Laboratory of Computational Physics(Grant No.6142A0502020717)National Science Foundation of USA(Grant No.DMS-1620108)
文摘This paper proposes and analyzes an efficient finite difference scheme for the two-dimensional nonlinear Schr?dinger(NLS) equation involving fractional Laplacian. The scheme is based on a weighted and shifted Grünwald-Letnikov difference(WSGD) operator for the spatial fractional Laplacian. We prove that the proposed method preserves the mass and energy conservation laws in semi-discrete formulations. By introducing the differentiation matrices, the semi-discrete fractional nonlinear Schr?dinger(FNLS) equation can be rewritten as a system of nonlinear ordinary differential equations(ODEs) in matrix formulations. Two kinds of time discretization methods are proposed for the semi-discrete formulation. One is based on the Crank-Nicolson(CN) method which can be proved to preserve the fully discrete mass and energy conservation. The other one is the compact implicit integration factor(c IIF) method which demands much less computational effort. It can be shown that the cIIF scheme can approximate CN scheme with the error O(τ~2). Finally numerical results are presented to demonstrate the method’s conservation, accuracy, efficiency and the capability of capturing blow-up.
基金partially supported by the Ministerio de Economia,Industria y Competitividad,Agencia Estatal de Investigacion grant MTM2016-77278-P (FEDER)the Agència de Gestio d’Ajuts Universitaris i de Recerca grant 2017SGR1617+1 种基金the H2020 European Research Council grant MSCA-RISE-2017-777911partially supported by FCT/Portugal through the pro ject UID/MAT/04459/2013。
文摘We characterize the complex differential equations of the form dy/dx=a_(n)(x)y^)n_+a_(n-1)(x)y^(n-1)+…+a_(1)(x)y+a_(0)(x) where a_(j)(x) are meromorphic functions in the variable x for j = 0,..., n that admit either a Weierstrass first integral or a Weierstrass inverse integrating factor.
基金supported by the Ministry Of Higher Education Malaysia for the Fundamental Research Grant scheme,project No. 01-04-10-897FRthe NSF scholarship
文摘In this paper, we study the behavior of the solution at the crack edges for a nearly circular crack with developing cusps subject to shear loading. The problem of finding the resulting force can be written in the form of a hypersingular integral equation. The equation is then trans-formed into a similar equation over a circular region using conformal mapping. The equation is solved numerically for the unknown coefficients, which will later be used in finding the stress intensity factors. The sliding and tearing mode stress intensity factors are evaluated for cracks and displayed graphically. Our results seem to agree with the existing asymptotic solution.