اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدThe variational particle-mesh method for matching curves
35
0
0.0
(
0
)
تأليف
C. J. Cotter
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
Diffeomorphic matching (only one of several names for this technique) is a technique for non-rigid registration of curves and surfaces in which the curve or surface is embedded in the flow of a time-series of vector fields. One seeks the flow between two topologically-equivalent curves or surfaces which minimises some metric defined on the vector fields, emph{i.e.} the flow closest to the identity in some sense. In this paper, we describe a new particle-mesh discretisation for the evolution of the geodesic flow and the embedded shape. Particle-mesh algorithms are very natural for this problem because Lagrangian particles (particles moving with the flow) can represent the movement of the shape whereas the vector field is Eulerian and hence best represented on a static mesh. We explain the derivation of the method, and prove conservation properties: the discrete method has a set of conserved momenta corresponding to the particle-relabelling symmetry which converge to conserved quantities in the continuous problem. We also introduce a new discretisation for the geometric current matching condition of (Vaillant and Glaunes, 2005). We illustrate the method and the derived properties with numerical examples.
قيم البحث
اقرأ أيضاً
The radiative transfer equation models the interaction of radiation with scattering and absorbing media and has important applications in various fields in science and engineering. It is an integro-differential equation involving time, space and angu
lar variables and contains an integral term in angular directions while being hyperbolic in space. The challenges for its numerical solution include the needs to handle with its high dimensionality, the presence of the integral term, and the development of discontinuities and sharp layers in its solution along spatial directions. Its numerical solution is studied in this paper using an adaptive moving mesh discontinuous Galerkin method for spatial discretization together with the discrete ordinate method for angular discretization. The former employs a dynamic mesh adaptation strategy based on moving mesh partial differential equations to improve computational accuracy and efficiency. Its mesh adaptation ability, accuracy, and efficiency are demonstrated in a selection of one- and two-dimensional numerical examples.
A Lagrangian-type numerical scheme called the comoving mesh method or CMM is developed for numerically solving certain classes of moving boundary problems which include, for example, the classical Hele-Shaw flow problem and the well-known mean curvat
ure flow problem. This finite element scheme exploits the idea that the normal velocity field of the moving boundary can be extended throughout the entire domain of definition of the problem using, for instance, the Laplace operator. Then, the boundary as well as the finite element mesh of the domain are easily updated at every time step by moving the nodal points along this velocity field. The feasibility of the method, highlighting its practicality, is illustrated through various numerical experiments. Also, in order to examine the accuracy of the proposed scheme, the experimental order of convergences between the numerical and manufactured solutions for these examples are also calculated.
The Half-Space Matching (HSM) method has recently been developed as a new method for the solution of 2D scattering problems with complex backgrounds, providing an alternative to Perfectly Matched Layers (PML) or other artificial boundary conditions.
Based on half-plane representations for the solution, the scattering problem is rewritten as a system of integral equations in which the unknowns are restrictions of the solution to the boundaries of a finite number of overlapping half-planes contained in the domain: this integral equation system is coupled to a standard finite element discretisation localised around the scatterer. While satisfactory numerical results have been obtained for real wavenumbers, wellposedness and equivalence to the original scattering problem have been established only for complex wavenumbers. In the present paper, by combining the HSM framework with a complex-scaling technique, we provide a new formulation for real wavenumbers which is provably well-posed and has the attraction for computation that the complex-scaled solutions of the integral equation system decay exponentially at infinity. The analysis requires the study of double-layer potential integral operators on intersecting infinite lines, and their analytic continuations. The effectiveness of the method is validated by preliminary numerical results.
A moving mesh discontinuous Galerkin method is presented for the numerical solution of hyperbolic conservation laws. The method is a combination of the discontinuous Galerkin method and the mesh movement strategy which is based on the moving mesh par
tial differential equation approach and moves the mesh continuously in time and orderly in space. It discretizes hyperbolic conservation laws on moving meshes in the quasi-Lagrangian fashion with which the mesh movement is treated continuously and no interpolation is needed for physical variables from the old mesh to the new one. Two convection terms are induced by the mesh movement and their discretization is incorporated naturally in the DG formulation. Numerical results for a selection of one- and two-dimensional scalar and system conservation laws are presented. It is shown that the moving mesh DG method achieves the theoretically predicted order of convergence for problems with smooth solutions and is able to capture shocks and concentrate mesh points in non-smooth regions. Its advantage over uniform meshes and its insensitiveness to mesh smoothness are also demonstrated.
We investigate several important issues regarding the Random Batch Method (RBM) for second order interacting particle systems. We first show the uniform-in-time strong convergence for second order systems under suitable contraction conditions. Second
ly, we propose the application of RBM for singular interaction kernels via kernel splitting strategy, and investigate numerically the application to molecular dynamics.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
