2,899 research outputs found
A variational Bayesian method for inverse problems with impulsive noise
We propose a novel numerical method for solving inverse problems subject to
impulsive noises which possibly contain a large number of outliers. The
approach is of Bayesian type, and it exploits a heavy-tailed t distribution for
data noise to achieve robustness with respect to outliers. A hierarchical model
with all hyper-parameters automatically determined from the given data is
described. An algorithm of variational type by minimizing the Kullback-Leibler
divergence between the true posteriori distribution and a separable
approximation is developed. The numerical method is illustrated on several one-
and two-dimensional linear and nonlinear inverse problems arising from heat
conduction, including estimating boundary temperature, heat flux and heat
transfer coefficient. The results show its robustness to outliers and the fast
and steady convergence of the algorithm.Comment: 20 pages, to appear in J. Comput. Phy
Numerical resolution of an exact heat conduction model with a delay term
In this paper we analyze, from the numerical point of view, a dynamic thermoelastic problem. Here, the so-called exact heat conduction model with a delay term is used to obtain the heat evolution. Thus, the thermomechanical problem is written as a coupled system of partial differential equations, and its variational formulation leads to a system written in terms of the velocity and the temperature fields. An existence and uniqueness result is recalled. Then, fully discrete approximations are introduced by using the classical finite element method to approximate the spatial variable and the implicit Euler scheme to discretize the time derivatives. A priori error estimates are proved, from which the linear convergence of the algorithm could be derived under suitable additional regularity conditions. Finally, a two-dimensional numerical example is solved to show the accuracy of the approximation and the decay of the discrete energy.Peer ReviewedPostprint (published version
Recovery of a space-dependent vector source in thermoelastic systems
In this contribution, an inverse problem of determining a space-dependent vector source in a thermoelastic system of type-I, type-II and type-III is studied using information from a supplementary measurement at a fixed time. These thermoelastic systems consist of two equations that are coupled: a parabolic equation for the temperature
[GRAPHICS]
and a vectorial hyperbolic equation for the displacement
[GRAPHICS]
. In this latter one, the source is unknown, but solely space dependent. A spacewise-dependent additional measurement at the final time ensures that the inverse problem corresponding with each type of thermoelasticity has a unique solution when a damping term
[GRAPHICS]
(with
[GRAPHICS]
componentwise strictly monotone increasing) is present in the hyperbolic equation. Despite the ill-posed nature of these inverse problems, a stable iterative algorithm is proposed to recover the unknown source in the case that
[GRAPHICS]
is also linear. This method is based on a sequence of well-posed direct problems, which are numerically solved at each iteration, step by step, using the finite element method. The instability of the inverse source problem is overcome by stopping the iterations at the first iteration for which the discrepancy principle is satisfied. Numerical results support the theoretically obtained results
General treatment of essential boundary conditions in reduced order models for non-linear problems
Inhomogeneous essential boundary conditions must be carefully treated in the formulation of Reduced Order Models (ROMs) for non-linear problems. In order to investigate this issue, two methods are analysed: one in which the boundary conditions are imposed in an strong way, and a second one in which a weak imposition of boundary conditions is made. The ideas presented in this work apply to the big realm of a posteriori ROMs. Nevertheless, an a posteriori hyper-reduction method is specifically considered in order to deal with the cost associated to the non-linearity of the problems. Applications to nonlinear transient heat conduction problems with temperature dependent thermophysical properties and time dependent essential boundary conditions are studied. However, the strategies introduced in this work are of general application.Fil: Cosimo, Alejandro. Consejo Nacional de Investigaciones CientÃficas y Técnicas. Centro CientÃfico Tecnológico Conicet - Santa Fe. Centro de Investigaciones en Métodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones en Métodos Computacionales; ArgentinaFil: Cardona, Alberto. Consejo Nacional de Investigaciones CientÃficas y Técnicas. Centro CientÃfico Tecnológico Conicet - Santa Fe. Centro de Investigaciones en Métodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones en Métodos Computacionales; ArgentinaFil: Idelsohn, Sergio Rodolfo. Consejo Nacional de Investigaciones CientÃficas y Técnicas. Centro CientÃfico Tecnológico Conicet - Santa Fe. Centro de Investigaciones en Métodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones en Métodos Computacionales; Argentina. Institució Catalana de Recerca i Estudis Avancats; España. International Center for Numerical Methods in Engineering; Españ
Topology and grid adaption for high-speed flow computations
This study investigates the effects of grid topology and grid adaptation on numerical solutions of the Navier-Stokes equations. In the first part of this study, a general procedure is presented for computation of high-speed flow over complex three-dimensional configurations. The flow field is simulated on the surface of a Butler wing in a uniform stream. Results are presented for Mach number 3.5 and a Reynolds number of 2,000,000. The O-type and H-type grids have been used for this study, and the results are compared together and with other theoretical and experimental results. The results demonstrate that while the H-type grid is suitable for the leading and trailing edges, a more accurate solution can be obtained for the middle part of the wing with an O-type grid. In the second part of this study, methods of grid adaption are reviewed and a method is developed with the capability of adapting to several variables. This method is based on a variational approach and is an algebraic method. Also, the method has been formulated in such a way that there is no need for any matrix inversion. This method is used in conjunction with the calculation of hypersonic flow over a blunt-nose body. A movie has been produced which shows simultaneously the transient behavior of the solution and the grid adaption
Numerical analysis of a thermoelastic problem with dual-phase-lag heat conduction
In this paper we study, from the numerical point of view, a thermoelastic problem with dual-phase-lag heat conduction. The variational formulation is written as a coupled system of hyperbolic linear variational equations. An existence and uniqueness result is recalled. Then, fully discrete approximations are introduced by using the finite element method and the implicit Euler scheme. A discrete stability result is proved and a priori error estimates are obtained, from which the linear convergence of the algorithm is deduced under suitable additional regularity conditions. Finally, some two-dimensional numerical simulations are presented to demonstrate the accuracy of the approximation and the behaviour of the solution.Peer ReviewedPostprint (author's final draft
Generalized Forward-Backward Splitting
This paper introduces the generalized forward-backward splitting algorithm
for minimizing convex functions of the form , where
has a Lipschitz-continuous gradient and the 's are simple in the sense
that their Moreau proximity operators are easy to compute. While the
forward-backward algorithm cannot deal with more than non-smooth
function, our method generalizes it to the case of arbitrary . Our method
makes an explicit use of the regularity of in the forward step, and the
proximity operators of the 's are applied in parallel in the backward
step. This allows the generalized forward backward to efficiently address an
important class of convex problems. We prove its convergence in infinite
dimension, and its robustness to errors on the computation of the proximity
operators and of the gradient of . Examples on inverse problems in imaging
demonstrate the advantage of the proposed methods in comparison to other
splitting algorithms.Comment: 24 pages, 4 figure
- …