Inverse parabolic problems and discrete orthogonality.
Si considera un modello discreto (per elementi finiti) di un solido o un sistema strutturale perfettamente elastoplastico, con condizioni di snervamento «linearizzate a tratti», nell’ipotesi di olonomia assunta per processi di caricamento proporzionali. Supponendo noti su base sperimentale certi spostamenti sotto assegnate azioni esterne, si formula il problema di identificare i limiti di snervamento, ossia le resistenze locali. Si dimostra che questo problema inverso di meccanica strutturale non...
The ill-posed problem of solving linear equations in the space of vector-valued finite Radon measures with Hilbert space data is considered. Approximate solutions are obtained by minimizing the Tikhonov functional with a total variation penalty. The well-posedness of this regularization method and further regularization properties are mentioned. Furthermore, a flexible numerical minimization algorithm is proposed which converges subsequentially in the weak* sense and with rate 𝒪(n-1)...
We consider the inverse scattering of time-harmonic plane waves to reconstruct the shape of a sound-soft crack from a knowledge of the given incident field and the phaseless data, and we check the invariance of far field data with respect to translation of the crack. We present a numerical method that is based on a system of nonlinear and ill-posed integral equations, and our scheme is easy and simple to implement. The numerical implementation is described and numerical examples are presented to...
We consider the problem of determining the unknown source term in a space fractional diffusion equation from the measured data at the final time . In this way, a methodology involving minimization of the cost functional is applied and shown that this cost functional is Fréchet differentiable and its derivative can be formulated via the solution of an adjoint problem. In addition, Lipschitz continuity of the gradient is proved. These results help us to prove the monotonicity and convergence...
In this paper, we propose a new numerical method for solving elliptic equations in unbounded regions of . The method is based on the mapping of a part of the domain into a bounded region. An appropriate family of weighted spaces is used for describing the growth or the decay of functions at large distances. After exposing the main ideas of the method, we analyse carefully its convergence. Some 3D computational results are displayed to demonstrate its efficiency and its high performance.
In this paper, we propose a new numerical method for solving elliptic equations in unbounded regions of . The method is based on the mapping of a part of the domain into a bounded region. An appropriate family of weighted spaces is used for describing the growth or the decay of functions at large distances. After exposing the main ideas of the method, we analyse carefully its convergence. Some 3D computational results are displayed to demonstrate its efficiency and its high performance.
We study a microfluidic flow model where the movement of several charged species is coupled with electric field and the motion of ambient fluid. The main numerical difficulty in this model is the net charge neutrality assumption which makes the system essentially overdetermined. Hence we propose to use the involutive and the associated augmented form of the system in numerical computations. Numerical experiments on electrophoresis and stacking show that the completed system significantly improves...
We study a microfluidic flow model where the movement of several charged species is coupled with electric field and the motion of ambient fluid. The main numerical difficulty in this model is the net charge neutrality assumption which makes the system essentially overdetermined. Hence we propose to use the involutive and the associated augmented form of the system in numerical computations. Numerical experiments on electrophoresis and stacking show that the completed system significantly improves...
These past few years, new types of computational architectures based on graphics processors have emerged. These technologies provide important computational resources at low cost and low energy consumption. Lots of developments have been done around GPU and many tools and libraries are now available to implement efficiently softwares on those architectures.This article contains the two contributions of the mini-symposium about GPU organized by Loïc Gouarin (Laboratoire de Mathématiques d’Orsay),...
The article is devoted to the simulation of viscous incompressible fluid flow based on solving the Navier-Stokes equations. As a numerical model we chose isogeometrical approach. Primary goal of using isogemetric analysis is to be always geometrically exact, independently of the discretization, and to avoid a time-consuming generation of meshes of computational domains. For higher Reynolds numbers, we use stabilization techniques SUPG and PSPG. All methods mentioned in the paper are demonstrated...
Estimates for the combined effect of boundary approximation and numerical integration on the approximation of (simple) eigenvalues and eigenvectors of 4th order eigenvalue problems with variable/constant coefficients in convex domains with curved boundary by an isoparametric mixed finite element method, which, in the particular case of bending problems of aniso-/ortho-/isotropic plates with variable/constant thickness, gives a simultaneous approximation to bending moment tensor field and displacement...
Estimates for the combined effect of boundary approximation and numerical integration on the approximation of (simple) eigenvalues and eigenvectors of 4th order eigenvalue problems with variable/constant coefficients in convex domains with curved boundary by an isoparametric mixed finite element method, which, in the particular case of bending problems of aniso-/ortho-/isotropic plates with variable/constant thickness, gives a simultaneous approximation to bending moment tensor field and displacement...