In the paper, we extend the implicit iterative method for linear ill-posed operator equations to solve nonlinear ill-posed problems. We show that under some conditions the error sequence of solutions of the nonlinear ...In the paper, we extend the implicit iterative method for linear ill-posed operator equations to solve nonlinear ill-posed problems. We show that under some conditions the error sequence of solutions of the nonlinear implicit iterative method is monotonically decreasing and, with this monotonicity, prove convergence of the new method for both the exact and perturbed equations.展开更多
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 meth...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 Schro¨dinger (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.展开更多
We applied the projection and contraction method to nonlinear complementarity problem (NCP). Moveover,we proposed an inexact implicit method for (NCP) and proved the convergence.
In this paper, an alternating direction Galerkin finite element method is presented for solving 2D time fractional reaction sub-diffusion equation with nonlinear source term. Firstly, one order implicit-explicit metho...In this paper, an alternating direction Galerkin finite element method is presented for solving 2D time fractional reaction sub-diffusion equation with nonlinear source term. Firstly, one order implicit-explicit method is used for time discretization, then Galerkin finite element method is adopted for spatial discretization and obtain a fully discrete linear system. Secondly, Galerkin alternating direction procedure for the system is derived by adding an extra term. Finally, the stability and convergence of the method are analyzed rigorously. Numerical results confirm the accuracy and efficiency of the proposed method.展开更多
We study the propagation of N-soliton bound state in a triangular gradient refractive index waveguide with nonlocal nonlinearity. The study is based on the direct numerical solutions of the model and subsequent eigenv...We study the propagation of N-soliton bound state in a triangular gradient refractive index waveguide with nonlocal nonlinearity. The study is based on the direct numerical solutions of the model and subsequent eigenvalues evolution of the corresponding Zakharov-Shabat spectral problem. In the waveguide with local nonlinearity, the velocity of a single soliton is found to be symmetric around zero and therefore the soliton oscillates periodically inside the waveguide. If the nonlocality is presence in the medium, the periodic motion of soliton is destroyed due to the soliton experiences additional positive acceleration induced by the nonlocality. In the waveguide with the same strength of nonlocality, a higher amplitude soliton experiences higher nonlocality effects, i.e. larger acceleration. Based on this soliton behavior we predict the break up of N-soliton bound state into their single-soliton constituents. We notice that the splitting process does not affect the amplitude of each soliton component.展开更多
In this work we propose a numerical scheme for a nonlinear and degenerate parabolic problem having application in petroleum reservoir and groundwater aquifer simulation. The degeneracy of the equation includes both lo...In this work we propose a numerical scheme for a nonlinear and degenerate parabolic problem having application in petroleum reservoir and groundwater aquifer simulation. The degeneracy of the equation includes both locally fast and slow diffusion (i.e. the diffusion coefficients may explode or vanish in some point). The main difficulty is that the true solution is typically lacking in regularity. Our numerical approach includes a regularization step and a standard discretization procedure by means of C0-piecewise linear finite elements in space and backward-differences in time. Within this frame work, we analyze the accuracy of the scheme by using an integral test function and obtain several error estimates in suitable norms.展开更多
In this paper, we are going to derive four numerical methods for solving the Modified Kortweg-de Vries (MKdV) equation using fourth Pade approximation for space direction and Crank Nicolson in the time direction. Two ...In this paper, we are going to derive four numerical methods for solving the Modified Kortweg-de Vries (MKdV) equation using fourth Pade approximation for space direction and Crank Nicolson in the time direction. Two nonlinear schemes and two linearized schemes are presented. All resulting schemes will be analyzed for accuracy and stability. The exact solution and the conserved quantities are used to highlight the efficiency and the robustness of the proposed schemes. Interaction of two and three solitons will be also conducted. The numerical results show that the interaction behavior is elastic and the conserved quantities are conserved exactly, and this is a good indication of the reliability of the schemes which we derived. A comparison with some existing is presented as well.展开更多
An inverse problem for identification of the coefficient in heat-conduction equation is considered. After reducing the problem to a nonlinear ill-posed operator equation, Newton type iterative methods are considered. ...An inverse problem for identification of the coefficient in heat-conduction equation is considered. After reducing the problem to a nonlinear ill-posed operator equation, Newton type iterative methods are considered. The implicit iterative method is applied to the linearized Newton equation, and the key step in the process is that a new reasonable a posteriori stopping rule for the inner iteration is presented. Numerical experiments for the new method as well as for Tikhonov method and Bakushikskii method are given, and these results show the obvious advantages of the new method over the other ones.展开更多
A PLU-SGS method based on a time-derivative preconditioning algorithm and LU-SGS method is developed in order to calculate the Navier-Stokes equations at all speeds. The equations were discretized using AUSMPW scheme ...A PLU-SGS method based on a time-derivative preconditioning algorithm and LU-SGS method is developed in order to calculate the Navier-Stokes equations at all speeds. The equations were discretized using AUSMPW scheme in conjunction with the third-order MUSCL scheme with Van Leer limiter. The present method was applied to solve the multidimensional compressible Navier-Stokes equations in curvilinear coordinates. Characteristic boundary conditions based on the eigensystem of the preconditioned equations were employed. In order to examine the performance of present method, driven-cavity flow at various Reynolds numbers and viscous flow through a convergent-divergent nozzle at supersonic were selected to test this method. The computed results were compared with the experimental data or the other numerical results available in literature and good agreements between them are obtained. The results show that the present method is accurate, self-adaptive and stable for a wide range of flow conditions from low speed to supersonic flows.展开更多
基金supported by the Key Disciplines of Shanghai Municipality (Operations Research & Cybernetics, No. S30104)the Shanghai Leading Academic Discipline Project (No. J50101)
文摘In the paper, we extend the implicit iterative method for linear ill-posed operator equations to solve nonlinear ill-posed problems. We show that under some conditions the error sequence of solutions of the nonlinear implicit iterative method is monotonically decreasing and, with this monotonicity, prove convergence of the new method for both the exact and perturbed equations.
文摘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 Schro¨dinger (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.
基金Supported by the National Natural Science Foundation of China (No. 202001036)
文摘We applied the projection and contraction method to nonlinear complementarity problem (NCP). Moveover,we proposed an inexact implicit method for (NCP) and proved the convergence.
文摘In this paper, an alternating direction Galerkin finite element method is presented for solving 2D time fractional reaction sub-diffusion equation with nonlinear source term. Firstly, one order implicit-explicit method is used for time discretization, then Galerkin finite element method is adopted for spatial discretization and obtain a fully discrete linear system. Secondly, Galerkin alternating direction procedure for the system is derived by adding an extra term. Finally, the stability and convergence of the method are analyzed rigorously. Numerical results confirm the accuracy and efficiency of the proposed method.
文摘We study the propagation of N-soliton bound state in a triangular gradient refractive index waveguide with nonlocal nonlinearity. The study is based on the direct numerical solutions of the model and subsequent eigenvalues evolution of the corresponding Zakharov-Shabat spectral problem. In the waveguide with local nonlinearity, the velocity of a single soliton is found to be symmetric around zero and therefore the soliton oscillates periodically inside the waveguide. If the nonlocality is presence in the medium, the periodic motion of soliton is destroyed due to the soliton experiences additional positive acceleration induced by the nonlocality. In the waveguide with the same strength of nonlocality, a higher amplitude soliton experiences higher nonlocality effects, i.e. larger acceleration. Based on this soliton behavior we predict the break up of N-soliton bound state into their single-soliton constituents. We notice that the splitting process does not affect the amplitude of each soliton component.
文摘In this work we propose a numerical scheme for a nonlinear and degenerate parabolic problem having application in petroleum reservoir and groundwater aquifer simulation. The degeneracy of the equation includes both locally fast and slow diffusion (i.e. the diffusion coefficients may explode or vanish in some point). The main difficulty is that the true solution is typically lacking in regularity. Our numerical approach includes a regularization step and a standard discretization procedure by means of C0-piecewise linear finite elements in space and backward-differences in time. Within this frame work, we analyze the accuracy of the scheme by using an integral test function and obtain several error estimates in suitable norms.
文摘In this paper, we are going to derive four numerical methods for solving the Modified Kortweg-de Vries (MKdV) equation using fourth Pade approximation for space direction and Crank Nicolson in the time direction. Two nonlinear schemes and two linearized schemes are presented. All resulting schemes will be analyzed for accuracy and stability. The exact solution and the conserved quantities are used to highlight the efficiency and the robustness of the proposed schemes. Interaction of two and three solitons will be also conducted. The numerical results show that the interaction behavior is elastic and the conserved quantities are conserved exactly, and this is a good indication of the reliability of the schemes which we derived. A comparison with some existing is presented as well.
文摘An inverse problem for identification of the coefficient in heat-conduction equation is considered. After reducing the problem to a nonlinear ill-posed operator equation, Newton type iterative methods are considered. The implicit iterative method is applied to the linearized Newton equation, and the key step in the process is that a new reasonable a posteriori stopping rule for the inner iteration is presented. Numerical experiments for the new method as well as for Tikhonov method and Bakushikskii method are given, and these results show the obvious advantages of the new method over the other ones.
文摘A PLU-SGS method based on a time-derivative preconditioning algorithm and LU-SGS method is developed in order to calculate the Navier-Stokes equations at all speeds. The equations were discretized using AUSMPW scheme in conjunction with the third-order MUSCL scheme with Van Leer limiter. The present method was applied to solve the multidimensional compressible Navier-Stokes equations in curvilinear coordinates. Characteristic boundary conditions based on the eigensystem of the preconditioned equations were employed. In order to examine the performance of present method, driven-cavity flow at various Reynolds numbers and viscous flow through a convergent-divergent nozzle at supersonic were selected to test this method. The computed results were compared with the experimental data or the other numerical results available in literature and good agreements between them are obtained. The results show that the present method is accurate, self-adaptive and stable for a wide range of flow conditions from low speed to supersonic flows.