10,019 research outputs found
Finite-Element Discretization of Static Hamilton-Jacobi Equations Based on a Local Variational Principle
We propose a linear finite-element discretization of Dirichlet problems for
static Hamilton-Jacobi equations on unstructured triangulations. The
discretization is based on simplified localized Dirichlet problems that are
solved by a local variational principle. It generalizes several approaches
known in the literature and allows for a simple and transparent convergence
theory. In this paper the resulting system of nonlinear equations is solved by
an adaptive Gauss-Seidel iteration that is easily implemented and quite
effective as a couple of numerical experiments show.Comment: 19 page
An approximation algorithm for the solution of the nonlinear Lane-Emden type equations arising in astrophysics using Hermite functions collocation method
In this paper we propose a collocation method for solving some well-known
classes of Lane-Emden type equations which are nonlinear ordinary differential
equations on the semi-infinite domain. They are categorized as singular initial
value problems. The proposed approach is based on a Hermite function
collocation (HFC) method. To illustrate the reliability of the method, some
special cases of the equations are solved as test examples. The new method
reduces the solution of a problem to the solution of a system of algebraic
equations. Hermite functions have prefect properties that make them useful to
achieve this goal. We compare the present work with some well-known results and
show that the new method is efficient and applicable.Comment: 34 pages, 13 figures, Published in "Computer Physics Communications
Numerical integration of asymptotic solutions of ordinary differential equations
Classical asymptotic analysis of ordinary differential equations derives approximate solutions that are numerically stable. However, the analysis also leads to tedious expansions in powers of the relevant parameter for a particular problem. The expansions are replaced with integrals that can be evaluated by numerical integration. The resulting numerical solutions retain the linear independence that is the main advantage of asymptotic solutions. Examples, including the Falkner-Skan equation from laminar boundary layer theory, illustrate the method of asymptotic analysis with numerical integration
Modal interaction in postbuckled plates. Theory
Plates can have more than one buckled solution for a fixed set of boundary conditions. The theory for the identification and the computation of multiple solutions in buckled plates is examined. The theory predicts modal interaction (which is also called change in buckle pattern or secondary buckling) in experiments on certain plates with multiple theoretical solutions. A set of coordinate functions is defined for Galerkin's method so that the von Karman plate equations are reduced to a coupled set of cubic equations in generalized coordinates that are uncoupled in the linear terms. An iterative procedure for solving modal interaction problems is suggested based on this cubic form
A combined first and second order variational approach for image reconstruction
In this paper we study a variational problem in the space of functions of
bounded Hessian. Our model constitutes a straightforward higher-order extension
of the well known ROF functional (total variation minimisation) to which we add
a non-smooth second order regulariser. It combines convex functions of the
total variation and the total variation of the first derivatives. In what
follows, we prove existence and uniqueness of minimisers of the combined model
and present the numerical solution of the corresponding discretised problem by
employing the split Bregman method. The paper is furnished with applications of
our model to image denoising, deblurring as well as image inpainting. The
obtained numerical results are compared with results obtained from total
generalised variation (TGV), infimal convolution and Euler's elastica, three
other state of the art higher-order models. The numerical discussion confirms
that the proposed higher-order model competes with models of its kind in
avoiding the creation of undesirable artifacts and blocky-like structures in
the reconstructed images -- a known disadvantage of the ROF model -- while
being simple and efficiently numerically solvable.Comment: 34 pages, 89 figure
- …