Many important problems in science and engineering require solving the so-called parametric partial differential equations(PDEs),i.e.,PDEs with different physical parameters,boundary conditions,shapes of computational...Many important problems in science and engineering require solving the so-called parametric partial differential equations(PDEs),i.e.,PDEs with different physical parameters,boundary conditions,shapes of computational domains,etc.Typical reduced order modeling techniques accelerate the solution of the parametric PDEs by projecting them onto a linear trial manifold constructed in the ofline stage.These methods often need a predefined mesh as well as a series of precomputed solution snapshots,and may struggle to balance between the efficiency and accuracy due to the limitation of the linear ansatz.Utilizing the nonlinear representation of neural networks(NNs),we propose the Meta-Auto-Decoder(MAD)to construct a nonlinear trial manifold,whose best possible performance is measured theoretically by the decoder width.Based on the meta-learning concept,the trial manifold can be learned in a mesh-free and unsupervised way during the pre-training stage.Fast adaptation to new(possibly heterogeneous)PDE parameters is enabled by searching on this trial manifold,and optionally fine-tuning the trial manifold at the same time.Extensive numerical experiments show that the MAD method exhibits a faster convergence speed without losing the accuracy than other deep learning-based methods.展开更多
In this paper, we shall be concerned with the numerical solution of parabolic equations in one space variable and the time variable t. We expand Taylor series to derive a higher-order approximation for U<sub>t&l...In this paper, we shall be concerned with the numerical solution of parabolic equations in one space variable and the time variable t. We expand Taylor series to derive a higher-order approximation for U<sub>t</sub>. We begin with the simplest model problem, for heat conduction in a uniform medium. For this model problem, an explicit difference method is very straightforward in use, and the analysis of its error is easily accomplished by the use of a maximum principle. As we shall show, however, the numerical solution becomes unstable unless the time step is severely restricted, so we shall go on to consider other, more elaborate, numerical methods which can avoid such a restriction. The additional complication in the numerical calculation is more than offset by the smaller number of time steps needed. We then extend the methods to problems with more general boundary conditions, then to more general linear parabolic equations. Finally, we shall discuss the more difficult problem of the solution of nonlinear equations.展开更多
The leaderless and leader-following finite-time consensus problems for multiagent systems(MASs)described by first-order linear hyperbolic partial differential equations(PDEs)are studied.The Lyapunov theorem and the un...The leaderless and leader-following finite-time consensus problems for multiagent systems(MASs)described by first-order linear hyperbolic partial differential equations(PDEs)are studied.The Lyapunov theorem and the unique solvability result for the first-order linear hyperbolic PDE are used to obtain some sufficient conditions for ensuring the finite-time consensus of the leaderless and leader-following MASs driven by first-order linear hyperbolic PDEs.Finally,two numerical examples are provided to verify the effectiveness of the proposed methods.展开更多
In this paper, a numerical solution of nonlinear partial differential equation, Benjamin-Bona-Mahony (BBM) and Cahn-Hilliard equation is presented by using Adomain Decomposition Method (ADM) and Variational Iteration ...In this paper, a numerical solution of nonlinear partial differential equation, Benjamin-Bona-Mahony (BBM) and Cahn-Hilliard equation is presented by using Adomain Decomposition Method (ADM) and Variational Iteration Method (VIM). The results reveal that the two methods are very effective, simple and very close to the exact solution.展开更多
We considerer partial differential equations of second order, for example the Klein-Gordon equation, the Poisson equation, on a region E = (a1, b1 ) × (a2, b2 ) x (a3, b3 ). We will see that with a common p...We considerer partial differential equations of second order, for example the Klein-Gordon equation, the Poisson equation, on a region E = (a1, b1 ) × (a2, b2 ) x (a3, b3 ). We will see that with a common procedure in all cases, we can write the equation in partial derivatives as an Fredholm integral equation of first kind and will solve this latter with the techniques of inverse problem moments. We will find an approximated solution and bounds for the error of the estimated solution using the techniques on problem of moments.展开更多
Lie algorithm combined with differential form Wu's method is used to complete the symmetry classification of partial differential equations(PDEs)containing arbitrary parameter.This process can be reduced to solve ...Lie algorithm combined with differential form Wu's method is used to complete the symmetry classification of partial differential equations(PDEs)containing arbitrary parameter.This process can be reduced to solve a large system of determining equations,which seems rather difficult to solve,then the differential form Wu's method is used to decompose the determining equations into a series of equations,which are easy to solve.To illustrate the usefulness of this method,we apply it to some test problems,and the results show the performance of the present work.展开更多
In this paper, by using the optimal stopping theory, the semilinear Black-Scholes partial differential equation (PDE) was invesigated in a fixed domain for valuing two assets of American (call-max/put-min) options...In this paper, by using the optimal stopping theory, the semilinear Black-Scholes partial differential equation (PDE) was invesigated in a fixed domain for valuing two assets of American (call-max/put-min) options. From the viscosity solution of a PDE, a unique viscosity solution was obtained for the semilinear Black-Scholes PDE.展开更多
A neural network(NN) is a powerful tool for approximating bounded continuous functions in machine learning. The NN provides a framework for numerically solving ordinary differential equations(ODEs) and partial differe...A neural network(NN) is a powerful tool for approximating bounded continuous functions in machine learning. The NN provides a framework for numerically solving ordinary differential equations(ODEs) and partial differential equations(PDEs)combined with the automatic differentiation(AD) technique. In this work, we explore the use of NN for the function approximation and propose a universal solver for ODEs and PDEs. The solver is tested for initial value problems and boundary value problems of ODEs, and the results exhibit high accuracy for not only the unknown functions but also their derivatives. The same strategy can be used to construct a PDE solver based on collocation points instead of a mesh, which is tested with the Burgers equation and the heat equation(i.e., the Laplace equation).展开更多
High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of th...High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of the WENO algorithm,large amount of computational costs are required for solving multidimensional problems.In our previous work(Lu et al.in Pure Appl Math Q 14:57–86,2018;Zhu and Zhang in J Sci Comput 87:44,2021),sparse-grid techniques were applied to the classical finite difference WENO schemes in solving multidimensional hyperbolic equations,and it was shown that significant CPU times were saved,while both accuracy and stability of the classical WENO schemes were maintained for computations on sparse grids.In this technical note,we apply the approach to recently developed finite difference multi-resolution WENO scheme specifically the fifth-order scheme,which has very interesting properties such as its simplicity in linear weights’construction over a classical WENO scheme.Numerical experiments on solving high dimensional hyperbolic equations including Vlasov based kinetic problems are performed to demonstrate that the sparse-grid computations achieve large savings of CPU times,and at the same time preserve comparable accuracy and resolution with those on corresponding regular single grids.展开更多
In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier ser...In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier series. The technique described in this work is suitable for modeling initial-boundary value problems governed by the wave equation on a rectangular domain with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The new numerical scheme is based on the standard approach of decomposing the global initial-boundary value problem into a steady-state component and a time-dependent component. The steady-state component is governed by the Laplace PDE and is modeled with the CVBEM. The time-dependent component is governed by the wave PDE and is modeled using a generalized Fourier series. The approximate global solution is the sum of the CVBEM and generalized Fourier series approximations. The boundary conditions of the steady-state component are specified as the boundary conditions from the global BVP. The boundary conditions of the time-dependent component are specified to be identically zero. The initial condition of the time-dependent component is calculated as the difference between the global initial condition and the CVBEM approximation of the steady-state solution. Additionally, the generalized Fourier series approximation of the time-dependent component is fitted so as to approximately satisfy the derivative of the initial condition. It is shown that the strong formulation of the wave PDE is satisfied by the superposed approximate solutions of the time-dependent and steady-state components.展开更多
The construction of complex stratigraphic surfaces is widely employed in many fields, such as petroleum exploration, geological modeling, and geological structure analysis. It also serves as an important foundation fo...The construction of complex stratigraphic surfaces is widely employed in many fields, such as petroleum exploration, geological modeling, and geological structure analysis. It also serves as an important foundation for data visualization and visual analysis in these fields. The existing surface construction methods have several deficiencies and face various difficulties, such as the presence of multitype faults and roughness of resulting surfaces. In this paper, a surface modeling method that uses geometric partial differential equations (PDEs) is introduced for the construction of stratigraphic surfaces. It effectively solves the problem of surface roughness caused by the irregularity of stratigraphic data distribution. To cope with the presence of multitype complex faults, a two-way projection algorithm between three- dimensional space and a two-dimensional plane is proposed. Using this algorithm, a unified method based on geometric PDEs is developed for dealing with multitype faults. Moreover, the corresponding geometric PDE is derived, and an algorithm based on an evolutionary solution is developed. The algorithm proposed for constructing spatial surfaces with real data verifies its computational efficiency and its ability to handle irregular data distribution. In particular, it can reconstruct faulty surfaces, especially those with overthrust faults.展开更多
In this paper the Black Scholes differential equation is transformed into a parabolic heat equation by appropriate change in variables. The transformed equation is semi-discretized by the Method of Lines (MOL). The ev...In this paper the Black Scholes differential equation is transformed into a parabolic heat equation by appropriate change in variables. The transformed equation is semi-discretized by the Method of Lines (MOL). The evolving system of ordinary differential equations (ODEs) is integrated numerically by an L-stable trapezoidal-like integrator. Results show accuracy of relative maximum error of order 10–10.展开更多
基金supported by the National Key R&D Program of China under Grant No.2021ZD0110400.
文摘Many important problems in science and engineering require solving the so-called parametric partial differential equations(PDEs),i.e.,PDEs with different physical parameters,boundary conditions,shapes of computational domains,etc.Typical reduced order modeling techniques accelerate the solution of the parametric PDEs by projecting them onto a linear trial manifold constructed in the ofline stage.These methods often need a predefined mesh as well as a series of precomputed solution snapshots,and may struggle to balance between the efficiency and accuracy due to the limitation of the linear ansatz.Utilizing the nonlinear representation of neural networks(NNs),we propose the Meta-Auto-Decoder(MAD)to construct a nonlinear trial manifold,whose best possible performance is measured theoretically by the decoder width.Based on the meta-learning concept,the trial manifold can be learned in a mesh-free and unsupervised way during the pre-training stage.Fast adaptation to new(possibly heterogeneous)PDE parameters is enabled by searching on this trial manifold,and optionally fine-tuning the trial manifold at the same time.Extensive numerical experiments show that the MAD method exhibits a faster convergence speed without losing the accuracy than other deep learning-based methods.
文摘In this paper, we shall be concerned with the numerical solution of parabolic equations in one space variable and the time variable t. We expand Taylor series to derive a higher-order approximation for U<sub>t</sub>. We begin with the simplest model problem, for heat conduction in a uniform medium. For this model problem, an explicit difference method is very straightforward in use, and the analysis of its error is easily accomplished by the use of a maximum principle. As we shall show, however, the numerical solution becomes unstable unless the time step is severely restricted, so we shall go on to consider other, more elaborate, numerical methods which can avoid such a restriction. The additional complication in the numerical calculation is more than offset by the smaller number of time steps needed. We then extend the methods to problems with more general boundary conditions, then to more general linear parabolic equations. Finally, we shall discuss the more difficult problem of the solution of nonlinear equations.
基金the National Natural Science Foundation of China(Nos.11671282 and 12171339)。
文摘The leaderless and leader-following finite-time consensus problems for multiagent systems(MASs)described by first-order linear hyperbolic partial differential equations(PDEs)are studied.The Lyapunov theorem and the unique solvability result for the first-order linear hyperbolic PDE are used to obtain some sufficient conditions for ensuring the finite-time consensus of the leaderless and leader-following MASs driven by first-order linear hyperbolic PDEs.Finally,two numerical examples are provided to verify the effectiveness of the proposed methods.
文摘In this paper, a numerical solution of nonlinear partial differential equation, Benjamin-Bona-Mahony (BBM) and Cahn-Hilliard equation is presented by using Adomain Decomposition Method (ADM) and Variational Iteration Method (VIM). The results reveal that the two methods are very effective, simple and very close to the exact solution.
文摘We considerer partial differential equations of second order, for example the Klein-Gordon equation, the Poisson equation, on a region E = (a1, b1 ) × (a2, b2 ) x (a3, b3 ). We will see that with a common procedure in all cases, we can write the equation in partial derivatives as an Fredholm integral equation of first kind and will solve this latter with the techniques of inverse problem moments. We will find an approximated solution and bounds for the error of the estimated solution using the techniques on problem of moments.
基金National Natural Science Foundation of China(No.61862048)。
文摘Lie algorithm combined with differential form Wu's method is used to complete the symmetry classification of partial differential equations(PDEs)containing arbitrary parameter.This process can be reduced to solve a large system of determining equations,which seems rather difficult to solve,then the differential form Wu's method is used to decompose the determining equations into a series of equations,which are easy to solve.To illustrate the usefulness of this method,we apply it to some test problems,and the results show the performance of the present work.
基金Project supported by the National Natural Science Foundation of China (Grant No.10271072)
文摘In this paper, by using the optimal stopping theory, the semilinear Black-Scholes partial differential equation (PDE) was invesigated in a fixed domain for valuing two assets of American (call-max/put-min) options. From the viscosity solution of a PDE, a unique viscosity solution was obtained for the semilinear Black-Scholes PDE.
基金Project supported by the National Natural Science Foundation of China(No.11521091)
文摘A neural network(NN) is a powerful tool for approximating bounded continuous functions in machine learning. The NN provides a framework for numerically solving ordinary differential equations(ODEs) and partial differential equations(PDEs)combined with the automatic differentiation(AD) technique. In this work, we explore the use of NN for the function approximation and propose a universal solver for ODEs and PDEs. The solver is tested for initial value problems and boundary value problems of ODEs, and the results exhibit high accuracy for not only the unknown functions but also their derivatives. The same strategy can be used to construct a PDE solver based on collocation points instead of a mesh, which is tested with the Burgers equation and the heat equation(i.e., the Laplace equation).
文摘High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of the WENO algorithm,large amount of computational costs are required for solving multidimensional problems.In our previous work(Lu et al.in Pure Appl Math Q 14:57–86,2018;Zhu and Zhang in J Sci Comput 87:44,2021),sparse-grid techniques were applied to the classical finite difference WENO schemes in solving multidimensional hyperbolic equations,and it was shown that significant CPU times were saved,while both accuracy and stability of the classical WENO schemes were maintained for computations on sparse grids.In this technical note,we apply the approach to recently developed finite difference multi-resolution WENO scheme specifically the fifth-order scheme,which has very interesting properties such as its simplicity in linear weights’construction over a classical WENO scheme.Numerical experiments on solving high dimensional hyperbolic equations including Vlasov based kinetic problems are performed to demonstrate that the sparse-grid computations achieve large savings of CPU times,and at the same time preserve comparable accuracy and resolution with those on corresponding regular single grids.
文摘In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier series. The technique described in this work is suitable for modeling initial-boundary value problems governed by the wave equation on a rectangular domain with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The new numerical scheme is based on the standard approach of decomposing the global initial-boundary value problem into a steady-state component and a time-dependent component. The steady-state component is governed by the Laplace PDE and is modeled with the CVBEM. The time-dependent component is governed by the wave PDE and is modeled using a generalized Fourier series. The approximate global solution is the sum of the CVBEM and generalized Fourier series approximations. The boundary conditions of the steady-state component are specified as the boundary conditions from the global BVP. The boundary conditions of the time-dependent component are specified to be identically zero. The initial condition of the time-dependent component is calculated as the difference between the global initial condition and the CVBEM approximation of the steady-state solution. Additionally, the generalized Fourier series approximation of the time-dependent component is fitted so as to approximately satisfy the derivative of the initial condition. It is shown that the strong formulation of the wave PDE is satisfied by the superposed approximate solutions of the time-dependent and steady-state components.
基金financially supported by the National Natural Science foundation of China(No.U1562218)
文摘The construction of complex stratigraphic surfaces is widely employed in many fields, such as petroleum exploration, geological modeling, and geological structure analysis. It also serves as an important foundation for data visualization and visual analysis in these fields. The existing surface construction methods have several deficiencies and face various difficulties, such as the presence of multitype faults and roughness of resulting surfaces. In this paper, a surface modeling method that uses geometric partial differential equations (PDEs) is introduced for the construction of stratigraphic surfaces. It effectively solves the problem of surface roughness caused by the irregularity of stratigraphic data distribution. To cope with the presence of multitype complex faults, a two-way projection algorithm between three- dimensional space and a two-dimensional plane is proposed. Using this algorithm, a unified method based on geometric PDEs is developed for dealing with multitype faults. Moreover, the corresponding geometric PDE is derived, and an algorithm based on an evolutionary solution is developed. The algorithm proposed for constructing spatial surfaces with real data verifies its computational efficiency and its ability to handle irregular data distribution. In particular, it can reconstruct faulty surfaces, especially those with overthrust faults.
文摘In this paper the Black Scholes differential equation is transformed into a parabolic heat equation by appropriate change in variables. The transformed equation is semi-discretized by the Method of Lines (MOL). The evolving system of ordinary differential equations (ODEs) is integrated numerically by an L-stable trapezoidal-like integrator. Results show accuracy of relative maximum error of order 10–10.