No Arabic abstract
In this paper we proposed two new quasi-boundary value methods for regularizing the ill-posed backward heat conduction problems. With a standard finite difference discretization in space and time, the obtained all-at-once nonsymmetric sparse linear systems have the desired block $omega$-circulant structure, which can be utilized to design an efficient parallel-in-time (PinT) direct solver that built upon an explicit FFT-based diagonalization of the time discretization matrix. Convergence analysis is presented to justify the optimal choice of the regularization parameter. Numerical examples are reported to validate our analysis and illustrate the superior computational efficiency of our proposed PinT methods.
In this paper, we present a unified analysis of the superconvergence property for a large class of mixed discontinuous Galerkin methods. This analysis applies to both the Poisson equation and linear elasticity problems with symmetric stress formulations. Based on this result, some locally postprocess schemes are employed to improve the accuracy of displacement by order min(k+1, 2) if polynomials of degree k are employed for displacement. Some numerical experiments are carried out to validate the theoretical results.
This paper presents a steady-state and transient heat conduction analysis framework using the polygonal scaled boundary finite element method (PSBFEM) with polygon/quadtree meshes. The PSBFEM is implemented with commercial finite element code Abaqus by the User Element Sub-routine (UEL) feature. The detailed implementation of the framework, defining the UEL element, and solving the stiffness/mass matrix by the eigenvalue decomposition are presented. Several benchmark problems from heat conduction are solved to validate the proposed implementation. Results show that the PSBFEM is reliable and accurate for solving heat conduction problems. Not only can the proposed implementation help engineering practitioners analyze the heat conduction problem using polygonal mesh in Abaqus, but it also provides a reference for developing the UEL to solve other problems using the scaled boundary finite element method.
This work focuses on the construction of a new class of fourth-order accurate methods for multirate time evolution of systems of ordinary differential equations. We base our work on the Recursive Flux Splitting Multirate (RFSMR) version of the Multirate Infinitesimal Step (MIS) methods and use recent theoretical developments for Generalized Additive Runge-Kutta methods to propose our higher-order Relaxed Multirate Infinitesimal Step extensions. The resulting framework supports a range of attractive properties for multirate methods, including telescopic extensions, subcycling, embeddings for temporal error estimation, and support for changes to the fast/slow time-scale separation between steps, without requiring any sacrifices in linear stability. In addition to providing rigorous theoretical developments for these new methods, we provide numerical tests demonstrating convergence and efficiency on a suite of multirate test problems.
Discrete variational methods have shown an excellent performance in numerical simulations of different mechanical systems. In this paper, we introduce an iterative method for discrete variational methods appropriate for boundary value problems. More concretely, we explore a parallelization strategy that leverages the power of multicore CPUs and GPUs (graphics cards). We study this parallel method for first-order and second-order Lagrangians and we illustrate its excellent behavior in some interesting applications, namely Zermelos navigation problem, a fuel-optimal navigation problem, and an interpolation problem.
Parameter identification problems for partial differential equations are an important subclass of inverse problems. The parameter-to-state map, which maps the parameter of interest to the respective solution of the PDE or state of the system, plays the central role in the (usually nonlinear) forward operator. Consequently, one is interested in well-definedness and further analytic properties such as continuity and differentiability of this operator w.r.t. the parameter in order to make sure that techniques from inverse problems theory may be successfully applied to solve the inverse problem. In this work, we present a general functional analytic framework suited for the study of a huge class of parameter identification problems including a variety of elliptic boundary value problems (in divergence form) with Dirichlet, Neumann, Robin or mixed boundary conditions. In particular, we show that the corresponding parameter-to-state operators fulfil, under suitable conditions, the tangential cone condition, which is often postulated for numerical solution techniques. This framework particularly covers the inverse medium problem and an inverse problem that arises in terahertz tomography.