Direct solution of th-order IVPs by homotopy analysis method.
In this paper, a new numerical method for solving the nonlinear constrained optimal control with quadratic performance index is presented. The method is based upon B-spline functions. The properties of B-spline functions are presented. The operational matrix of derivative () and integration matrix () are introduced. These matrices are utilized to reduce the solution of nonlinear constrained quadratic optimal control to the solution of nonlinear programming one to which existing well-developed...
This paper is devoted to computational problems related to Markov chains (MC) on a finite state space. We present formulas and bounds for characteristics of MCs using directed forest expansions given by the Matrix Tree Theorem. These results are applied to analysis of direct methods for solving systems of linear equations, aggregation algorithms for nearly completely decomposable MCs and the Markov chain Monte Carlo procedures.
Although many words have been written about two recent directional (regression) quantile concepts, their applications, and the algorithms for computing associated (regression) quantile regions, their software implementation is still not widely available, which, of course, severely hinders the dissemination of both methods. Wanting to partly fill in the gap here, we provide all the codes needed for computing and plotting the multivariate (regression) quantile regions in Octave and MATLAB, describe...
Recently, the eminently popular standard quantile regression has been generalized to the multiple-output regression setup by means of directional regression quantiles in two rather interrelated ways. Unfortunately, they lead to complicated optimization problems involving parametric programming, and this may be the main obstacle standing in the way of their wide dissemination. The presented R package modQR is intended to address this issue. It originates as a quite faithful translation of the authors'...
In this paper, we study the boundary penalty method for optimal control of unsteady Navier–Stokes type system that has been proposed as an alternative for Dirichlet boundary control. Existence and uniqueness of solutions are demonstrated and existence of optimal control for a class of optimal control problems is established. The asymptotic behavior of solution, with respect to the penalty parameter ϵ, is studied. In particular, we prove convergence of solutions of penalized control problem to the...
In this paper, by the Kirchhoff transformation, a Dirichlet-Neumann (D-N) alternating algorithm which is a non-overlapping domain decomposition method based on natural boundary reduction is discussed for solving exterior anisotropic quasilinear problems with circular artificial boundary. By the principle of the natural boundary reduction, we obtain natural integral equation for the anisotropic quasilinear problems on circular artificial boundaries and construct the algorithm and analyze its convergence....
We propose a discontinuous Galerkin method for linear elasticity, based on discontinuous piecewise linear approximation of the displacements. We show optimal order a priori error estimates, uniform in the incompressible limit, and thus locking is avoided. The discontinuous Galerkin method is closely related to the non-conforming Crouzeix–Raviart (CR) element, which in fact is obtained when one of the stabilizing parameters tends to infinity. In the case of the elasticity operator, for which the...
We propose a discontinuous Galerkin method for linear elasticity, based on discontinuous piecewise linear approximation of the displacements. We show optimal order a priori error estimates, uniform in the incompressible limit, and thus locking is avoided. The discontinuous Galerkin method is closely related to the non-conforming Crouzeix–Raviart (CR) element, which in fact is obtained when one of the stabilizing parameters tends to infinity. In the case of the elasticity operator, for...
We present our work on the numerical solution of a continuum model of flocking dynamics in two spatial dimensions. The model consists of the compressible Euler equations with a nonlinear nonlocal term which requires special treatment. We use a semi-implicit discontinuous Galerkin scheme, which proves to be efficient enough to produce results in 2D in reasonable time. This work is a direct extension of the authors' previous work in 1D.
The paper is devoted to the analysis of the discontinuous Galerkin finite element method (DGFEM) applied to the space semidiscretization of a nonlinear nonstationary convection-diffusion problem with mixed Dirichlet-Neumann boundary conditions. General nonconforming meshes are used and the NIPG, IIPG and SIPG versions of the discretization of diffusion terms are considered. The main attention is paid to the impact of the Neumann boundary condition prescribed on a part of the boundary on the truncation...
In this article we study discontinuous Galerkin finite element discretizations of linear second-order elliptic partial differential equations with Dirac delta right-hand side. In particular, assuming that the underlying computational mesh is quasi-uniform, we derive an a priori bound on the error measured in terms of the L2-norm. Additionally, we develop residual-based a posteriori error estimators that can be used within an adaptive mesh refinement ...
In this article we study discontinuous Galerkin finite element discretizations of linear second-order elliptic partial differential equations with Dirac delta right-hand side. In particular, assuming that the underlying computational mesh is quasi-uniform, we derive an a priori bound on the error measured in terms of the L2-norm. Additionally, we develop residual-based a posteriori error estimators that can be used within an adaptive mesh refinement ...
We present a discontinuous Galerkin scheme for the numerical approximation of the one-dimensional periodic Vlasov-Poisson equation. The scheme is based on a Galerkin-characteristics method in which the distribution function is projected onto a space of discontinuous functions. We present comparisons with a semi-Lagrangian method to emphasize the good behavior of this scheme when applied to Vlasov-Poisson test cases.