Localized Method Of Approximate Particular Solutions For Solving Optimal Control Problems Governed By PDES

In this thesis, the method of approximate particular solutions(MAPS) and localized method of approximate particular solutions(LMAPS) with polynomial basis, and radial basis functions are proposed and applied on the optimal control problems(OCPs) governed by partial differential equations(PDEs). This study proceeds in several steps. First, polynomial basis and radial basis functions are used to globally approximate solutions for the PDEs which have been combined into a single matrix system numerically from the optimality conditions of the OCPs. Secondly, polynomial and radial basis functions are used to locally approximate particular solutions for the same matrix system numerically. We use these approaches to two types of problems, a smooth and singular problem. The first example numerically experiments on a square domain and the second example on an L-shaped disc domain. These approaches are tested and compared. The results show our proposed method for solving optimal control problems governed by partial differential equations works

Hybrid Chebyshev Polynomial Scheme for the Numerical Solution of Partial Differential Equations

In the numerical solution of partial differential equations (PDEs), it is common to find situations where the best choice is to use more than one method to arrive at an accurate solution. In this dissertation, hybrid Chebyshev polynomial scheme (HCPS) is proposed which is applied in two-step approach and one-step approach. In the two-step approach, first, Chebyshev polynomials are used to approximate a particular solution of a PDE. Chebyshev nodes which are the roots of Chebyshev polynomials are used in the polynomial interpolation due to its spectral convergence. Then, the resulting homogeneous equation is solved by boundary type methods including the method of fundamental solution (MFS) and the equilibrated collocation Trefftz method. However, this scheme can be applied to solve PDEs with constant coefficients only. So, for solving a wide variety of PDEs, one-step hybrid Chebyshev polynomial scheme is proposed. This approach combines two matrix systems of two-step approach into a single matrix system. The solution is approximated by the sum of particular solution and homogeneous solution. The Laplacian or biharmonic operator is kept on the left hand side and all the other terms are moved to the right hand side and treated as the forcing term. Various boundary value problems governed by the Poisson equation in two and three dimensions are considered for the numerical experiments. HCPS is also applied to solve an inhomogeneous Cauchy-Navier equations of elasticity in two dimensions. Numerical results show that HCPS is direct, easy to implement, and highly accurate

Moving-boundary problems solved by adaptive radial basis functions

The objective of this paper is to present an alternative approach to the conventional level set methods for solving two-dimensional moving-boundary problems known as the passive transport. Moving boundaries are associated with time-dependent problems and the position of the boundaries need to be determined as a function of time and space. The level set method has become an attractive design tool for tracking, modeling and simulating the motion of free boundaries in fluid mechanics, combustion, computer animation and image processing. Recent research on the numerical method has focused on the idea of using a meshless methodology for the numerical solution of partial differential equations. In the present approach, the moving interface is captured by the level set method at all time with the zero contour of a smooth function known as the level set function. A new approach is used to solve a convective transport equation for advancing the level set function in time. This new approach is based on the asymmetric meshless collocation method and the adaptive greedy algorithm for trial subspaces selection. Numerical simulations are performed to verify the accuracy and stability of the new numerical scheme which is then applied to simulate a bubble that is moving, stretching and circulating in an ambient flow to demonstrate the performance of the new meshless approach. (C) 2010 Elsevier Ltd. All rights reserved

Hybrid PDE solver for data-driven problems and modern branching

The numerical solution of large-scale PDEs, such as those occurring in data-driven applications, unavoidably require powerful parallel computers and tailored parallel algorithms to make the best possible use of them. In fact, considerations about the parallelization and scalability of realistic problems are often critical enough to warrant acknowledgement in the modelling phase. The purpose of this paper is to spread awareness of the Probabilistic Domain Decomposition (PDD) method, a fresh approach to the parallelization of PDEs with excellent scalability properties. The idea exploits the stochastic representation of the PDE and its approximation via Monte Carlo in combination with deterministic high-performance PDE solvers. We describe the ingredients of PDD and its applicability in the scope of data science. In particular, we highlight recent advances in stochastic representations for nonlinear PDEs using branching diffusions, which have significantly broadened the scope of PDD. We envision this work as a dictionary giving large-scale PDE practitioners references on the very latest algorithms and techniques of a non-standard, yet highly parallelizable, methodology at the interface of deterministic and probabilistic numerical methods. We close this work with an invitation to the fully nonlinear case and open research questions.Comment: 23 pages, 7 figures; Final SMUR version; To appear in the European Journal of Applied Mathematics (EJAM

Parameter estimation of partial differential equations using artificial neural network

The work presented in this paper aims at developing a novel meshless parameter estimation framework for a system of partial differential equations (PDEs) using artificial neural network (ANN) approximations. The PDE models to be treated consist of linear and nonlinear PDEs, with Dirichlet and Neumann boundary conditions, considering both regular and irregular boundaries. This paper focuses on testing the applicability of neural networks for estimating the process model parameters while simultaneously computing the model predictions of the state variables in the system of PDEs representing the process. The capability of the proposed methodology is demonstrated with five numerical problems, showing that the ANN-based approach is very efficient by providing accurate solutions in reasonable computing times

Adaptive radial basis function generated finite-difference (RBF-FD) on non-uniform nodes using pp-refinement

Radial basis functions-generated finite difference methods (RBF-FDs) have been gaining popularity recently. In particular, the RBF-FD based on polyharmonic splines (PHS) augmented with multivariate polynomials (PHS+poly) has been found significantly effective. For the approximation order of RBF-FDs' weights on scattered nodes, one can already find mathematical theories in the literature. Many practical problems in numerical analysis, however, do not require a uniform node-distribution. Instead, they would be better suited if specific areas of the domain, where complicated physics needed to be resolved, had a relatively higher node-density compared to the rest of the domain. In this work, we proposed a practical adaptive RBF-FD with a user-defined order of convergence with respect to the total number of (possibly scattered and non-uniform) data points $N$. Our algorithm outputs a sparse differentiation matrix with the desired approximation order. Numerical examples are provided to show that the proposed adaptive RBF-FD method yields the expected $N$-convergence even for highly non-uniform node-distributions. The proposed method also reduces the number of non-zero elements in the linear system without sacrificing accuracy.Comment: An updated version with seismic modeling will be included in version