132,288 research outputs found
Adjoint recovery of superconvergent functionals from PDE approximations
Motivated by applications in computational fluid dynamics, a method is presented for obtaining estimates of integral functionals, such as lift or drag, that have twice the order of accuracy of the computed flow solution on which they are based. This is achieved through error analysis that uses an adjoint PDE to relate the local errors in approximating the flow solution to the corresponding global errors in the functional of interest. Numerical evaluation of the local residual error together with an approximate solution to the adjoint equations may thus be combined to produce a correction for the computed functional value that yields the desired improvement in accuracy. Numerical results are presented for the Poisson equation in one and two dimensions and for the nonlinear quasi-one-dimensional Euler equations. The theory is equally applicable to nonlinear equations in complex multi-dimensional domains and holds great promise for use in a range of engineering disciplines in which a few integral quantities are a key output of numerical approximations
Analysis of Adjoint Error Correction for Superconvergent Functional Estimates
Earlier work introduced the notion of adjoint error correction for obtaining superconvergent estimates of functional outputs from approximate PDE solutions. This idea is based on a posteriori error analysis suggesting that the leading order error term in the functional estimate can be removed by using an adjoint PDE solution to reveal the sensitivity of the functional to the residual error in the original PDE solution. The present work provides a priori error analysis that correctly predicts the behaviour of the remaining leading order error term. Furthermore, the discussion is extended from the case of homogeneous boundary conditions and bulk functionals, to encompass the possibilities of inhomogeneous boundary conditions and boundary functionals. Numerical illustrations are provided for both linear and nonlinear problems.\ud
\ud
This research was supported by EPSRC under grant GR/K91149, and by NASA/Ames Cooperative Agreement No. NCC 2-5431
Semidefinite Relaxations for Stochastic Optimal Control Policies
Recent results in the study of the Hamilton Jacobi Bellman (HJB) equation
have led to the discovery of a formulation of the value function as a linear
Partial Differential Equation (PDE) for stochastic nonlinear systems with a
mild constraint on their disturbances. This has yielded promising directions
for research in the planning and control of nonlinear systems. This work
proposes a new method obtaining approximate solutions to these linear
stochastic optimal control (SOC) problems. A candidate polynomial with variable
coefficients is proposed as the solution to the SOC problem. A Sum of Squares
(SOS) relaxation is then taken to the partial differential constraints, leading
to a hierarchy of semidefinite relaxations with improving sub-optimality gap.
The resulting approximate solutions are shown to be guaranteed over- and
under-approximations for the optimal value function.Comment: Preprint. Accepted to American Controls Conference (ACC) 2014 in
Portland, Oregon. 7 pages, colo
Linearly Solvable Stochastic Control Lyapunov Functions
This paper presents a new method for synthesizing stochastic control Lyapunov
functions for a class of nonlinear stochastic control systems. The technique
relies on a transformation of the classical nonlinear Hamilton-Jacobi-Bellman
partial differential equation to a linear partial differential equation for a
class of problems with a particular constraint on the stochastic forcing. This
linear partial differential equation can then be relaxed to a linear
differential inclusion, allowing for relaxed solutions to be generated using
sum of squares programming. The resulting relaxed solutions are in fact
viscosity super/subsolutions, and by the maximum principle are pointwise upper
and lower bounds to the underlying value function, even for coarse polynomial
approximations. Furthermore, the pointwise upper bound is shown to be a
stochastic control Lyapunov function, yielding a method for generating
nonlinear controllers with pointwise bounded distance from the optimal cost
when using the optimal controller. These approximate solutions may be computed
with non-increasing error via a hierarchy of semidefinite optimization
problems. Finally, this paper develops a-priori bounds on trajectory
suboptimality when using these approximate value functions, as well as
demonstrates that these methods, and bounds, can be applied to a more general
class of nonlinear systems not obeying the constraint on stochastic forcing.
Simulated examples illustrate the methodology.Comment: Published in SIAM Journal of Control and Optimizatio
Accurate and efficient evaluation of the a posteriori error estimator in the reduced basis method
The reduced basis method is a model reduction technique yielding substantial
savings of computational time when a solution to a parametrized equation has to
be computed for many values of the parameter. Certification of the
approximation is possible by means of an a posteriori error bound. Under
appropriate assumptions, this error bound is computed with an algorithm of
complexity independent of the size of the full problem. In practice, the
evaluation of the error bound can become very sensitive to round-off errors. We
propose herein an explanation of this fact. A first remedy has been proposed in
[F. Casenave, Accurate \textit{a posteriori} error evaluation in the reduced
basis method. \textit{C. R. Math. Acad. Sci. Paris} \textbf{350} (2012)
539--542.]. Herein, we improve this remedy by proposing a new approximation of
the error bound using the Empirical Interpolation Method (EIM). This method
achieves higher levels of accuracy and requires potentially less
precomputations than the usual formula. A version of the EIM stabilized with
respect to round-off errors is also derived. The method is illustrated on a
simple one-dimensional diffusion problem and a three-dimensional acoustic
scattering problem solved by a boundary element method.Comment: 26 pages, 10 figures. ESAIM: Mathematical Modelling and Numerical
Analysis, 201
Tchebychev Polynomial Approximations for Order Boundary Value Problems
Higher order boundary value problems (BVPs) play an important role modeling
various scientific and engineering problems. In this article we develop an
efficient numerical scheme for linear order BVPs. First we convert the
higher order BVP to a first order BVP. Then we use Tchebychev orthogonal
polynomials to approximate the solution of the BVP as a weighted sum of
polynomials. We collocate at Tchebychev clustered grid points to generate a
system of equations to approximate the weights for the polynomials. The
excellency of the numerical scheme is illustrated through some examples.Comment: 21 pages, 10 figure
- …