Based on the Dyson-Schwinger equations of QCD in the "rainbow" approximation, the fully dressed quarkpropagator Sf(p) is investigated, and then an algebraic parametrization form of the propagator is obtained...Based on the Dyson-Schwinger equations of QCD in the "rainbow" approximation, the fully dressed quarkpropagator Sf(p) is investigated, and then an algebraic parametrization form of the propagator is obtained as a solutionof the equations. The dressed quark amplitudes Af and Bf built up the fully dressed quark propagator and the dynamicalrunning masses Mf defined by Af and Bf for light quarks u, d and s are calculated, respectively. Using the predictedrunning masses Mf, quark condensates <0|q(0)q(0)|0> = -(0.255 GeV)a for u, d quarks, and <0|s s|0> = 0.8<0|q(0)q(0)]0)for s quark, and experimental pion decay constant fπ = 0.093 GeV, the masses of Goldstone bosons K, π, and η are alsoevaluated. The numerical results show that the masses of quarks are dependent on their momentum p2. The fully dressedquark amplitudes Af and Bf have correct behaviors which can be used for many purposes in our future researches onnonperturbative QCD.展开更多
We propose a novel inverse-free neurodynamic approach (NIFNA) for solving absolute value equations (AVE). The NIFNA guarantees global convergence and notably improves convergence speed by achieving fixed-time converge...We propose a novel inverse-free neurodynamic approach (NIFNA) for solving absolute value equations (AVE). The NIFNA guarantees global convergence and notably improves convergence speed by achieving fixed-time convergence. To validate the theoretical findings, numerical simulations are conducted, demonstrating the effectiveness and efficiency of the proposed NIFNA.展开更多
In this paper,we consider the high order method for solving the linear transport equations under diffusive scaling and with random inputs.To tackle the randomness in the problem,the stochastic Galerkin method of the g...In this paper,we consider the high order method for solving the linear transport equations under diffusive scaling and with random inputs.To tackle the randomness in the problem,the stochastic Galerkin method of the generalized polynomial chaos approach has been employed.Besides,the high order implicit-explicit scheme under the micro-macro decomposition framework and the discontinuous Galerkin method have been employed.We provide several numerical experiments to validate the accuracy and the stochastic asymptotic-preserving property.展开更多
Based on the Dyson-Schwinger equations of quark propagator in rainbow truncation with an effective gluonpropagator,the ten unknown Gasser-Leutwyler coefficients of the chiral Lagrangian for pseudoscalar Goldstone boso...Based on the Dyson-Schwinger equations of quark propagator in rainbow truncation with an effective gluonpropagator,the ten unknown Gasser-Leutwyler coefficients of the chiral Lagrangian for pseudoscalar Goldstone bosonsare predicted.The predicted values of L_i with i=1,2,...,10 are in a reasonable agreement with empirical values usedwidely in literature,and the values predicted by many other theoretical models with QCD characteristics.展开更多
The pion and tensor vacuum susceptibilities are calculated in the framework of the renormalizable DysonSchwinger equations. A comparison with the results of other nonperturbative QCD approaches is given.
We discuss the chiral phase transition of quantum chromodynamics (QCD) with a chiral chemical potential μ5 as an additional scale. Within the framework of Dyson-Schwinger equations, we focus particularly on the beh...We discuss the chiral phase transition of quantum chromodynamics (QCD) with a chiral chemical potential μ5 as an additional scale. Within the framework of Dyson-Schwinger equations, we focus particularly on the behavior of the widely accepted as well as interesting critical end point (CEP), using a separable gluon propagator and a Gaussian gluon propagator. We find that there may be no CEP5 in the T-μ5 plane, and the phase transition in the T μ5 plane might be totally crossover. Our results have apparent consistency with the Lattice QCD calculation. On the other hand, our study may also provide some useful hints to some other studies related to μ5.展开更多
In view of the properties of mesons in hot strongly interacting matter, the properties of the solutions of the truncated Dyson-Schwinger equation for the quark propagator at finite temperatures within the rainbow-ladd...In view of the properties of mesons in hot strongly interacting matter, the properties of the solutions of the truncated Dyson-Schwinger equation for the quark propagator at finite temperatures within the rainbow-ladder approximation are analysed in some detail. In Euclidean space within the Matsubara imaginary time formalism, the quark propagator is not longer a O(4) symmetric function and possesses a discrete spectrum of the fourth component of the momentum. This makes the treatment of the Dyson-Schwinger and Bethe-Salpeter equations conceptually different from the vacuum and technically much more involved. The question whether the interaction kernel known from vacuum calculations can be applied at finite temperatures remains still open. We find that, at low temperatures, the model interaction with vacuum parameters provides a reasonable description of the quark propagator, while at temperatures higher than a certain critical value T<sub>c </sub>the interaction requires stringent modifications. The general properties of the quark propagator at finite temperatures can be inferred from lattice QCD (LQCD) calculations. We argue that, to achieve a reasonable agreement of the model calculations with that from LQCD, the kernel is to be modified in such a way as to screen the infra-red part of the interaction at temperatures larger than T<sub>c </sub>. For this, we analyse the solutions of the truncated Dyson-Schwinger equation with existing interaction kernels in a large temperature range with particular attention on high temperatures in order to find hints to an adequate temperature dependence of the interaction kernel to be further implemented in the Bethe-Salpeter equation for mesons. This will allow investigating the possible in medium modifications of the meson properties as well as the conditions of quark deconfinement in hot matter.展开更多
The coupled modified nonlinear Schrodinger equations are under investigation in this work. Starting from analyzing the spectral problem of the Lax pair, a Riemann-Hilbert problem for the coupled modified nonlinear Sch...The coupled modified nonlinear Schrodinger equations are under investigation in this work. Starting from analyzing the spectral problem of the Lax pair, a Riemann-Hilbert problem for the coupled modified nonlinear Schrodinger equations is formulated. And then, through solving the obtained Riemann-Hilbert problem under the conditions of irregularity and reflectionless case, N-soliton solutions for the equations are presented. Furthermore, the localized structures and dynamic behaviors of the one-soliton solution are shown graphically.展开更多
Recently, the Clarkson and Kruskal direct method has been modified to find new similarity reductions (conditional similarity reductions) of nonlinear systems and the results obtained by the modified direct method cann...Recently, the Clarkson and Kruskal direct method has been modified to find new similarity reductions (conditional similarity reductions) of nonlinear systems and the results obtained by the modified direct method cannot be obtained by the current classical and/or non-classical Lie group approach. In this paper, we show that the conditional similarity reductions of the Jimbo-Miwa equation can be reobtained by adding an additional constraint equation to the original model to form a conditional equation system first and then solving the model system by means of the classical Lie group approach.展开更多
In this article, we are concerned with the existence of solutions of a quasilinear elliptic equation in R^N which includes the so-called modified nonlinear Schrodinger equation as a special case. Combining the dual ap...In this article, we are concerned with the existence of solutions of a quasilinear elliptic equation in R^N which includes the so-called modified nonlinear Schrodinger equation as a special case. Combining the dual approach and the nonsmooth critical point theory, we obtain the existence of a nontrivial solution.展开更多
A finite-difference approach is used to develop a time-dependent mild-slope equation incorporating the effects of bottom dissipation and nonlinearity. The Enler predietor-corrector method and the three-point finite-di...A finite-difference approach is used to develop a time-dependent mild-slope equation incorporating the effects of bottom dissipation and nonlinearity. The Enler predietor-corrector method and the three-point finite-difference method with varying spatial steps are adopted to discretize the time derivatives and the two-dimensional horizontal ones, respectively, thus leading both the time and spatial derivatives to the second-order accuracy. The boundary conditions for the present model are treated on the basis of the general conditions for open and fixed boundaries with an arbitrary reflection coefficient and phase shift. Both the linear and nonlinear versions of the numerical model are applied to the wave propagation and transformation over an elliptic shoal on a sloping beach, respectively, and the linear version is applied to the simulation of wave propagation in a fully open rectangular harbor. From comparison of numerical results with theoretical or experimental ones, it is found that they are in reasonable agreement.展开更多
In this paper we propose a numerical approach to solve the relativistic Dirac equation suitable for computational calculations of one-electron systems. A variational procedure is carried out similar to the well-known ...In this paper we propose a numerical approach to solve the relativistic Dirac equation suitable for computational calculations of one-electron systems. A variational procedure is carried out similar to the well-known Hylleraas computational method. An application of the method to hydrogen isoelectronic atoms is presented, showing its consistency and high accuracy, relative to the exact analytical eigenvalues.展开更多
One existence integral condition was obtained for the adapted solution of the general backward stochastic differential equations(BSDEs). Then by solving the integral constraint condition, and using a limit procedure, ...One existence integral condition was obtained for the adapted solution of the general backward stochastic differential equations(BSDEs). Then by solving the integral constraint condition, and using a limit procedure, a new approach method is proposed and the existence of the solution was proved for the BSDEs if the diffusion coefficients satisfy the locally Lipschitz condition. In the special case the solution was a Brownian bridge. The uniqueness is also considered in the meaning of "F0-integrable equivalent class" . The new approach method would give us an efficient way to control the main object instead of the "noise".展开更多
After considering the variable coefficient of a nonlinear equation as a new dependent variable, some special types of variable-coefficient equation can be solved from the corresponding constant-coefficient equations b...After considering the variable coefficient of a nonlinear equation as a new dependent variable, some special types of variable-coefficient equation can be solved from the corresponding constant-coefficient equations by using the general classical Lie approach. Taking the nonlinear Schr?dinger equation as a concrete example, the method is recommended in detail.展开更多
Employing the Pekeris-type approximation to deal with the pseudo-centrifugal term,we analytically study the pseudospin symmetry of a Dirac nucleon subjected to equal scalar and vector modified Rosen-Morse potential in...Employing the Pekeris-type approximation to deal with the pseudo-centrifugal term,we analytically study the pseudospin symmetry of a Dirac nucleon subjected to equal scalar and vector modified Rosen-Morse potential including the spin-orbit coupling term by using the Nikiforov-Uvarov method and supersymmetric quantum mechanics approach.The complex eigenvalue equation and the total normalized wave functions expressed in terms of Jacobi polynomial with arbitrary spin-orbit coupling quantum number k are presented under the condition of pseudospin symmetry.The eigenvalue equations for both methods reproduce the same result to affirm the mathematical accuracy of analytical calculations.The numerical solutions obtained for different adjustable parameters produce degeneracies for some quantum number.展开更多
The extended symmetry approach is used to study the general Korteweg-de Vries-type (KdV-type) equation. Several variable-coefficient equations are obtained. The solutions of these resulting equations can be construc...The extended symmetry approach is used to study the general Korteweg-de Vries-type (KdV-type) equation. Several variable-coefficient equations are obtained. The solutions of these resulting equations can be constructed by the solutions of original models if their solutions are well known, such as the standard constant coefficient KdV equation and the standard compound KdV--Burgers equation, and so on. Then any one of these variable-coefficient equations can be considered as an original model to obtain new variable-coefficient equations whose solutions can also be known by means of transformation relations between solutions of the resulting new variable-coefficient equations and the original equation.展开更多
With a new approach,the general current expressions of two typical second order catalytic reactions at microelectrodes are obtained.This allows the study of fast chemical reactions and systems where the reactants are ...With a new approach,the general current expressions of two typical second order catalytic reactions at microelectrodes are obtained.This allows the study of fast chemical reactions and systems where the reactants are present in similar concentrations.展开更多
Many applications in fluid mechanics require the numerical solution of sequences of linear systems typically issued from finite element discretization of the Navier-Stokes equations. The resulting matrices then exhibi...Many applications in fluid mechanics require the numerical solution of sequences of linear systems typically issued from finite element discretization of the Navier-Stokes equations. The resulting matrices then exhibit a saddle point structure. To achieve this task, a Newton-based root-finding algorithm is usually employed which in turn necessitates to solve a saddle point system at every Newton iteration. The involved linear systems being large scale and ill-conditioned, effective linear solvers must be implemented. Here, we develop and test several methods for solving the saddle point systems, considering in particular the LU factorization, as direct approach, and the preconditioned generalized minimal residual (ΡGMRES) solver, an iterative approach. We apply the various solvers within the root-finding algorithm for Flow over backward facing step systems. The particularity of Flow over backward facing step system is an interesting case for studying the performance and solution strategy of a turbulence model. In this case, the flow is subjected to a sudden increase of cross-sectional area, resulting in a separation of flow starting at the point of expansion, making the system of differential equations particularly stiff. We assess the performance of the direct and iterative solvers in terms of computational time, numbers of Newton iterations and time steps.展开更多
The (3+1)-dimensional Jimbo-Miwa (JM) equation is solved approximately by using the conformal invariant asymptotic expansion approach presented by Ruan. By solving the new (3+1)-dimensional integrable models, ...The (3+1)-dimensional Jimbo-Miwa (JM) equation is solved approximately by using the conformal invariant asymptotic expansion approach presented by Ruan. By solving the new (3+1)-dimensional integrable models, which are conformal invariant and possess Painlevé property, the approximate solutions are obtained for the JM equation, containing not only one-soliton solutions but also periodic solutions and multi-soliton solutions. Some approximate solutions happen to be exact and some approximate solutions can become exact by choosing relations between the parameters properly.展开更多
The Gravity Recovery and Climate Experiment(GRACE) mission can significantly improve our knowledge of the temporal variability of the Earth's gravity field.We obtained monthly gravity field solutions based on varia...The Gravity Recovery and Climate Experiment(GRACE) mission can significantly improve our knowledge of the temporal variability of the Earth's gravity field.We obtained monthly gravity field solutions based on variational equations approach from GPS-derived positions of GRACE satellites and K-band range-rate measurements.The impact of different fixed data weighting ratios in temporal gravity field recovery while combining the two types of data was investigated for the purpose of deriving the best combined solution.The monthly gravity field solution obtained through above procedures was named as the Institute of Geodesy and Geophysics(IGG) temporal gravity field models.IGG temporal gravity field models were compared with GRACE Release05(RL05) products in following aspects:(i) the trend of the mass anomaly in China and its nearby regions within 2005-2010; (ii) the root mean squares of the global mass anomaly during 2005-2010; (iii) time-series changes in the mean water storage in the region of the Amazon Basin and the Sahara Desert between 2005 and 2010.The results showed that IGG solutions were almost consistent with GRACE RL05 products in above aspects(i)-(iii).Changes in the annual amplitude of mean water storage in the Amazon Basin were 14.7 ± 1.2 cm for IGG,17.1 ± 1.3 cm for the Centre for Space Research(CSR),16.4 ± 0.9 cm for the GeoForschungsZentrum(GFZ) and 16.9 ± 1.2 cm for the Jet Propulsion Laboratory(JPL) in terms of equivalent water height(EWH),respectively.The root mean squares of the mean mass anomaly in Sahara were 1.2 cm,0.9 cm,0.9 cm and 1.2 cm for temporal gravity field models of IGG,CSR,GFZ and JPL,respectively.Comparison suggested that IGG temporal gravity field solutions were at the same accuracy level with the latest temporal gravity field solutions published by CSR,GFZ and JPL.展开更多
基金The project supported in part by National Natural Science Foundation of China under Grant Nos. 19975053, 19835010, 100750811007505, and the CAS Knowledge Innovation Pro jet No. KJCX2-SW-No2
文摘Based on the Dyson-Schwinger equations of QCD in the "rainbow" approximation, the fully dressed quarkpropagator Sf(p) is investigated, and then an algebraic parametrization form of the propagator is obtained as a solutionof the equations. The dressed quark amplitudes Af and Bf built up the fully dressed quark propagator and the dynamicalrunning masses Mf defined by Af and Bf for light quarks u, d and s are calculated, respectively. Using the predictedrunning masses Mf, quark condensates <0|q(0)q(0)|0> = -(0.255 GeV)a for u, d quarks, and <0|s s|0> = 0.8<0|q(0)q(0)]0)for s quark, and experimental pion decay constant fπ = 0.093 GeV, the masses of Goldstone bosons K, π, and η are alsoevaluated. The numerical results show that the masses of quarks are dependent on their momentum p2. The fully dressedquark amplitudes Af and Bf have correct behaviors which can be used for many purposes in our future researches onnonperturbative QCD.
文摘We propose a novel inverse-free neurodynamic approach (NIFNA) for solving absolute value equations (AVE). The NIFNA guarantees global convergence and notably improves convergence speed by achieving fixed-time convergence. To validate the theoretical findings, numerical simulations are conducted, demonstrating the effectiveness and efficiency of the proposed NIFNA.
基金supported by the Simons Foundation:Collaboration Grantssupported by the AFOSR grant FA9550-18-1-0383.
文摘In this paper,we consider the high order method for solving the linear transport equations under diffusive scaling and with random inputs.To tackle the randomness in the problem,the stochastic Galerkin method of the generalized polynomial chaos approach has been employed.Besides,the high order implicit-explicit scheme under the micro-macro decomposition framework and the discontinuous Galerkin method have been employed.We provide several numerical experiments to validate the accuracy and the stochastic asymptotic-preserving property.
基金supported in part by National Natural Science Foundation of China under Grant Nos.10647002 and 10565001the Natural Science Foundation of Guangxi under Grant Nos.0542042,0481030,and 0575020
文摘Based on the Dyson-Schwinger equations of quark propagator in rainbow truncation with an effective gluonpropagator,the ten unknown Gasser-Leutwyler coefficients of the chiral Lagrangian for pseudoscalar Goldstone bosonsare predicted.The predicted values of L_i with i=1,2,...,10 are in a reasonable agreement with empirical values usedwidely in literature,and the values predicted by many other theoretical models with QCD characteristics.
文摘The pion and tensor vacuum susceptibilities are calculated in the framework of the renormalizable DysonSchwinger equations. A comparison with the results of other nonperturbative QCD approaches is given.
基金Supported by the National Natural Science Foundation of China under Grant Nos 11275097,11475085,11265017,and 11247219the Jiangsu Planned Projects for Postdoctoral Research Funds under Grant No 1402006C+1 种基金the Natural Science Foundation of Jiangsu Province under Grant No BK20130078the Guizhou Province Outstanding Youth Science and Technology Talent Cultivation Object Special Funds under Grant No QKHRZ(2013)28
文摘We discuss the chiral phase transition of quantum chromodynamics (QCD) with a chiral chemical potential μ5 as an additional scale. Within the framework of Dyson-Schwinger equations, we focus particularly on the behavior of the widely accepted as well as interesting critical end point (CEP), using a separable gluon propagator and a Gaussian gluon propagator. We find that there may be no CEP5 in the T-μ5 plane, and the phase transition in the T μ5 plane might be totally crossover. Our results have apparent consistency with the Lattice QCD calculation. On the other hand, our study may also provide some useful hints to some other studies related to μ5.
文摘In view of the properties of mesons in hot strongly interacting matter, the properties of the solutions of the truncated Dyson-Schwinger equation for the quark propagator at finite temperatures within the rainbow-ladder approximation are analysed in some detail. In Euclidean space within the Matsubara imaginary time formalism, the quark propagator is not longer a O(4) symmetric function and possesses a discrete spectrum of the fourth component of the momentum. This makes the treatment of the Dyson-Schwinger and Bethe-Salpeter equations conceptually different from the vacuum and technically much more involved. The question whether the interaction kernel known from vacuum calculations can be applied at finite temperatures remains still open. We find that, at low temperatures, the model interaction with vacuum parameters provides a reasonable description of the quark propagator, while at temperatures higher than a certain critical value T<sub>c </sub>the interaction requires stringent modifications. The general properties of the quark propagator at finite temperatures can be inferred from lattice QCD (LQCD) calculations. We argue that, to achieve a reasonable agreement of the model calculations with that from LQCD, the kernel is to be modified in such a way as to screen the infra-red part of the interaction at temperatures larger than T<sub>c </sub>. For this, we analyse the solutions of the truncated Dyson-Schwinger equation with existing interaction kernels in a large temperature range with particular attention on high temperatures in order to find hints to an adequate temperature dependence of the interaction kernel to be further implemented in the Bethe-Salpeter equation for mesons. This will allow investigating the possible in medium modifications of the meson properties as well as the conditions of quark deconfinement in hot matter.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.61072147 and 11271008)
文摘The coupled modified nonlinear Schrodinger equations are under investigation in this work. Starting from analyzing the spectral problem of the Lax pair, a Riemann-Hilbert problem for the coupled modified nonlinear Schrodinger equations is formulated. And then, through solving the obtained Riemann-Hilbert problem under the conditions of irregularity and reflectionless case, N-soliton solutions for the equations are presented. Furthermore, the localized structures and dynamic behaviors of the one-soliton solution are shown graphically.
基金国家杰出青年科学基金,the Research Fund for the Doctoral Program of Higher Education of China
文摘Recently, the Clarkson and Kruskal direct method has been modified to find new similarity reductions (conditional similarity reductions) of nonlinear systems and the results obtained by the modified direct method cannot be obtained by the current classical and/or non-classical Lie group approach. In this paper, we show that the conditional similarity reductions of the Jimbo-Miwa equation can be reobtained by adding an additional constraint equation to the original model to form a conditional equation system first and then solving the model system by means of the classical Lie group approach.
基金supported partially by National Natural Science Foundation of China(11771385,11661083)the Youth Foundation of Yunnan Minzu University(2017QNo3)
文摘In this article, we are concerned with the existence of solutions of a quasilinear elliptic equation in R^N which includes the so-called modified nonlinear Schrodinger equation as a special case. Combining the dual approach and the nonsmooth critical point theory, we obtain the existence of a nontrivial solution.
基金This work wasjointlysupported by the National Natural Science Foundation of China(Grant No.40106008) the National Natural Science Fundfor Distinguished Young Scholars(Grant No.40225014)
文摘A finite-difference approach is used to develop a time-dependent mild-slope equation incorporating the effects of bottom dissipation and nonlinearity. The Enler predietor-corrector method and the three-point finite-difference method with varying spatial steps are adopted to discretize the time derivatives and the two-dimensional horizontal ones, respectively, thus leading both the time and spatial derivatives to the second-order accuracy. The boundary conditions for the present model are treated on the basis of the general conditions for open and fixed boundaries with an arbitrary reflection coefficient and phase shift. Both the linear and nonlinear versions of the numerical model are applied to the wave propagation and transformation over an elliptic shoal on a sloping beach, respectively, and the linear version is applied to the simulation of wave propagation in a fully open rectangular harbor. From comparison of numerical results with theoretical or experimental ones, it is found that they are in reasonable agreement.
文摘In this paper we propose a numerical approach to solve the relativistic Dirac equation suitable for computational calculations of one-electron systems. A variational procedure is carried out similar to the well-known Hylleraas computational method. An application of the method to hydrogen isoelectronic atoms is presented, showing its consistency and high accuracy, relative to the exact analytical eigenvalues.
基金National Natural Science Foundation of China ( No. 11171062 ) Natural Science Foundation for the Youth,China ( No.11101077) Innovation Program of Shanghai Municipal Education Commission,China ( No. 12ZZ063)
文摘One existence integral condition was obtained for the adapted solution of the general backward stochastic differential equations(BSDEs). Then by solving the integral constraint condition, and using a limit procedure, a new approach method is proposed and the existence of the solution was proved for the BSDEs if the diffusion coefficients satisfy the locally Lipschitz condition. In the special case the solution was a Brownian bridge. The uniqueness is also considered in the meaning of "F0-integrable equivalent class" . The new approach method would give us an efficient way to control the main object instead of the "noise".
基金国家自然科学基金,浙江省自然科学基金,Foundation of State Key Laboratory of Oil/Gas Reservoir Geology and Exploitation (PLN 0104),the Foundation of Educational Commission,浙江省宁波市博士基金
文摘After considering the variable coefficient of a nonlinear equation as a new dependent variable, some special types of variable-coefficient equation can be solved from the corresponding constant-coefficient equations by using the general classical Lie approach. Taking the nonlinear Schr?dinger equation as a concrete example, the method is recommended in detail.
文摘Employing the Pekeris-type approximation to deal with the pseudo-centrifugal term,we analytically study the pseudospin symmetry of a Dirac nucleon subjected to equal scalar and vector modified Rosen-Morse potential including the spin-orbit coupling term by using the Nikiforov-Uvarov method and supersymmetric quantum mechanics approach.The complex eigenvalue equation and the total normalized wave functions expressed in terms of Jacobi polynomial with arbitrary spin-orbit coupling quantum number k are presented under the condition of pseudospin symmetry.The eigenvalue equations for both methods reproduce the same result to affirm the mathematical accuracy of analytical calculations.The numerical solutions obtained for different adjustable parameters produce degeneracies for some quantum number.
基金Project supported by the National Natural Science Foundation of China (Grant No. 10675065)the Scientific Research Fundof the Education Department of Zhejiang Province of China (Grant No. 20070979)
文摘The extended symmetry approach is used to study the general Korteweg-de Vries-type (KdV-type) equation. Several variable-coefficient equations are obtained. The solutions of these resulting equations can be constructed by the solutions of original models if their solutions are well known, such as the standard constant coefficient KdV equation and the standard compound KdV--Burgers equation, and so on. Then any one of these variable-coefficient equations can be considered as an original model to obtain new variable-coefficient equations whose solutions can also be known by means of transformation relations between solutions of the resulting new variable-coefficient equations and the original equation.
文摘With a new approach,the general current expressions of two typical second order catalytic reactions at microelectrodes are obtained.This allows the study of fast chemical reactions and systems where the reactants are present in similar concentrations.
文摘Many applications in fluid mechanics require the numerical solution of sequences of linear systems typically issued from finite element discretization of the Navier-Stokes equations. The resulting matrices then exhibit a saddle point structure. To achieve this task, a Newton-based root-finding algorithm is usually employed which in turn necessitates to solve a saddle point system at every Newton iteration. The involved linear systems being large scale and ill-conditioned, effective linear solvers must be implemented. Here, we develop and test several methods for solving the saddle point systems, considering in particular the LU factorization, as direct approach, and the preconditioned generalized minimal residual (ΡGMRES) solver, an iterative approach. We apply the various solvers within the root-finding algorithm for Flow over backward facing step systems. The particularity of Flow over backward facing step system is an interesting case for studying the performance and solution strategy of a turbulence model. In this case, the flow is subjected to a sudden increase of cross-sectional area, resulting in a separation of flow starting at the point of expansion, making the system of differential equations particularly stiff. We assess the performance of the direct and iterative solvers in terms of computational time, numbers of Newton iterations and time steps.
基金The project supported by the Natural Science Foundation of Zhejiang Province of China under Grant No. Y604036 and State Key Laboratory of 0il/Gas Reservoir Geology and Exploitation "PLN0402" The authors would like to thank Prof. Sen-Yue Lou for his help and discussion.
文摘The (3+1)-dimensional Jimbo-Miwa (JM) equation is solved approximately by using the conformal invariant asymptotic expansion approach presented by Ruan. By solving the new (3+1)-dimensional integrable models, which are conformal invariant and possess Painlevé property, the approximate solutions are obtained for the JM equation, containing not only one-soliton solutions but also periodic solutions and multi-soliton solutions. Some approximate solutions happen to be exact and some approximate solutions can become exact by choosing relations between the parameters properly.
基金funded by the Major National Scientific Research Plan(2013CB733305,2012CB957703)the National Natural Science Foundation of China(41174066,41131067,41374087,41431070)
文摘The Gravity Recovery and Climate Experiment(GRACE) mission can significantly improve our knowledge of the temporal variability of the Earth's gravity field.We obtained monthly gravity field solutions based on variational equations approach from GPS-derived positions of GRACE satellites and K-band range-rate measurements.The impact of different fixed data weighting ratios in temporal gravity field recovery while combining the two types of data was investigated for the purpose of deriving the best combined solution.The monthly gravity field solution obtained through above procedures was named as the Institute of Geodesy and Geophysics(IGG) temporal gravity field models.IGG temporal gravity field models were compared with GRACE Release05(RL05) products in following aspects:(i) the trend of the mass anomaly in China and its nearby regions within 2005-2010; (ii) the root mean squares of the global mass anomaly during 2005-2010; (iii) time-series changes in the mean water storage in the region of the Amazon Basin and the Sahara Desert between 2005 and 2010.The results showed that IGG solutions were almost consistent with GRACE RL05 products in above aspects(i)-(iii).Changes in the annual amplitude of mean water storage in the Amazon Basin were 14.7 ± 1.2 cm for IGG,17.1 ± 1.3 cm for the Centre for Space Research(CSR),16.4 ± 0.9 cm for the GeoForschungsZentrum(GFZ) and 16.9 ± 1.2 cm for the Jet Propulsion Laboratory(JPL) in terms of equivalent water height(EWH),respectively.The root mean squares of the mean mass anomaly in Sahara were 1.2 cm,0.9 cm,0.9 cm and 1.2 cm for temporal gravity field models of IGG,CSR,GFZ and JPL,respectively.Comparison suggested that IGG temporal gravity field solutions were at the same accuracy level with the latest temporal gravity field solutions published by CSR,GFZ and JPL.