999 research outputs found
An -Adaptive Newton-Galerkin Finite Element Procedure for Semilinear Boundary Value Problems
In this paper we develop an -adaptive procedure for the numerical
solution of general, semilinear elliptic boundary value problems in 1d, with
possible singular perturbations. Our approach combines both a prediction-type
adaptive Newton method and an -version adaptive finite element
discretization (based on a robust a posteriori residual analysis), thereby
leading to a fully -adaptive Newton-Galerkin scheme. Numerical experiments
underline the robustness and reliability of the proposed approach for various
examples.Comment: arXiv admin note: text overlap with arXiv:1408.522
Adaptive Pseudo-Transient-Continuation-Galerkin Methods for Semilinear Elliptic Partial Differential Equations
In this paper we investigate the application of pseudo-transient-continuation
(PTC) schemes for the numerical solution of semilinear elliptic partial
differential equations, with possible singular perturbations. We will outline a
residual reduction analysis within the framework of general Hilbert spaces,
and, subsequently, employ the PTC-methodology in the context of finite element
discretizations of semilinear boundary value problems. Our approach combines
both a prediction-type PTC-method (for infinite dimensional problems) and an
adaptive finite element discretization (based on a robust a posteriori residual
analysis), thereby leading to a fully adaptive PTC-Galerkin scheme. Numerical
experiments underline the robustness and reliability of the proposed approach
for different examples.Comment: arXiv admin note: text overlap with arXiv:1408.522
A Numerical Slow Manifold Approach to Model Reduction for Optimal Control of Multiple Time Scale ODE
Time scale separation is a natural property of many control systems that can
be ex- ploited, theoretically and numerically. We present a numerical scheme to
solve optimal control problems with considerable time scale separation that is
based on a model reduction approach that does not need the system to be
explicitly stated in singularly perturbed form. We present examples that
highlight the advantages and disadvantages of the method
An integrated approach to global synchronization and state estimation for nonlinear singularly perturbed complex networks
This paper aims to establish a unified framework to handle both the exponential synchronization and state estimation problems for a class of nonlinear singularly perturbed complex networks (SPCNs). Each node in the SPCN comprises both 'slow' and 'fast' dynamics that reflects the singular perturbation behavior. General sector-like nonlinear function is employed to describe the nonlinearities existing in the network. All nodes in the SPCN have the same structures and properties. By utilizing a novel Lyapunov functional and the Kronecker product, it is shown that the addressed SPCN is synchronized if certain matrix inequalities are feasible. The state estimation problem is then studied for the same complex network, where the purpose is to design a state estimator to estimate the network states through available output measurements such that dynamics (both slow and fast) of the estimation error is guaranteed to be globally asymptotically stable. Again, a matrix inequality approach is developed for the state estimation problem. Two numerical examples are presented to verify the effectiveness and merits of the proposed synchronization scheme and state estimation formulation. It is worth mentioning that our main results are still valid even if the slow subsystems within the network are unstable
Fitted non-polynomial spline method for singularly perturbed differential difference equations with integral boundary condition
The aim of this paper is to present fitted non-polynomial spline method for singularly perturbed differential-difference equations with integral boundary condition. The stability and uniform convergence of the proposed method are proved. To validate the applicability of the scheme, two model problems are considered for numerical experimentation and solved for different values of the perturbation parameter, ε and mesh size, h. The numerical results are tabulated in terms of maximum absolute errors and rate of convergence and it is observed that the present method is more accurate and uniformly convergent for h ≥ ε where the classical numerical methods fails to give good result and it also improves the results of the methods existing in the literature
- …