In this paper Zadunaisky's technique is used to estimate the global error propagated in the numerical solution of the system of retarded differential equations by Euler's method. Some numerical examples are given.
In this article we transform a large class of parabolic inverse problems into a nonclassical parabolic equation whose coefficients consist of trace type functionals of the solution and its derivatives subject to some initial and boundary conditions. For this nonclassical problem, we study finite element methods and present an immediate analysis for global superconvergence for these problems, on basis of which we obtain a posteriori error estimators.
We describe both the classical lagrangian and the Eulerian methods for first order Hamilton–Jacobi equations of geometric optic type. We then explain the basic structure of the software and how new solvers/models can be added to it. A selection of numerical examples are presented.
We describe both the classical Lagrangian and the Eulerian methods for first order Hamilton–Jacobi equations of geometric optic type. We then explain the basic structure of the software and how new solvers/models can be added to it. A selection of numerical examples are presented.
This paper is concerned with goal-oriented a posteriori error estimates for discontinous Galerkin discretizations of linear elliptic boundary value problems. Our approach combines the Dual Weighted Residual method (DWR) with local weighted least-squares reconstruction of the discrete solution. This technique is used not only for controlling the discretization error, but also to track the influence of the algebraic errors. We illustrate the performance of the proposed method by numerical experiments....
We deal with a posteriori error control of discontinuous Galerkin approximations for linear boundary value problems. The computational error is estimated in the framework of the Dual Weighted Residual method (DWR) for goal-oriented error estimation which requires to solve an additional (adjoint) problem. We focus on the control of the algebraic errors arising from iterative solutions of algebraic systems corresponding to both the primal and adjoint problems. Moreover, we present two different reconstruction...
This paper is concerned with the numerical approximation of Cauchy problems for one-dimensional nonconservative hyperbolic systems. The theory developed by Dal Maso et al. [J. Math. Pures Appl.74 (1995) 483–548] is used in order to define the weak solutions of the system: an interpretation of the nonconservative products as Borel measures is given, based on the choice of a family of paths drawn in the phase space. Even if the family of paths can be chosen arbitrarily, it is natural to require this...
We describe a numerical technique for the solution of macroscopic traffic flow models on networks of roads. On individual roads, we consider the standard Lighthill-Whitham-Richards model which is discretized using the discontinuous Galerkin method along with suitable limiters. In order to solve traffic flows on networks, we construct suitable numerical fluxes at junctions based on preferences of the drivers. Numerical experiment comparing different approaches is presented.
Test procedures are constructed for testing the goodness-of-fit in parametric regression models. The test statistic is in the form of an L2 distance between the empirical characteristic function of the residuals in a parametric regression fit and the corresponding empirical characteristic function of the residuals in a non-parametric regression fit. The asymptotic null distribution as well as the behavior of the test statistic under contiguous alternatives is investigated. Theoretical results are...
The integration to steady state of many initial value ODEs and PDEs using the forward Euler method can alternatively be considered as gradient descent for an associated minimization problem. Greedy algorithms such as steepest descent for determining the step size are as slow to reach steady state as is forward Euler integration with the best uniform step size. But other, much faster methods using bolder step size selection exist. Various alternatives are investigated from both theoretical and practical...
The purpose of our work is to develop an automatic shape optimization tool for runner wheel blades in reaction water turbines, especially in Kaplan turbines. The fluid flow is simulated using an in-house incompressible turbulent flow solver based on recently introduced isogeometric analysis (see e.g. J. A. Cotrell et al.: Isogeometric Analysis: Toward Integration of CAD and FEA, Wiley, 2009). The proposed automatic shape optimization approach is based on a so-called hybrid optimization which combines...