This paper considers the newly introduced generalized scalar auxiliary variable approaches to construct high-efficiency energy-preserving schemes for the sine-Gordon equation with Neumann boundary conditions. The equa...
详细信息
This paper considers the newly introduced generalized scalar auxiliary variable approaches to construct high-efficiency energy-preserving schemes for the sine-Gordon equation with Neumann boundary conditions. The equation is first reformulated into an equivalent system by defining a new auxiliary variable that is not limited to square root. Then, the cosine pseudo-spectral method is applied to the system and derive a semi-discrete conservative scheme. Subsequently, we combine the auxiliary variable with the nonlinear term and use an explicit technique discretization in time to derive a fully-discrete energy-preserving scheme. Furthermore, a fast algorithm based on the discrete cosine transform technique reduces the computational complexity in practical computation. Finally, various numerical experiments are displayed to verify the accuracy, efficiency and conservation of the proposed schemes.
The proper orthogonal decomposition reduced-order model (POD-ROM) has been widely used as a computationally efficient surrogate model in large-scale numerical simulations of complex systems. However, when it is applie...
详细信息
The proper orthogonal decomposition reduced-order model (POD-ROM) has been widely used as a computationally efficient surrogate model in large-scale numerical simulations of complex systems. However, when it is applied to a Hamiltonian system, a naive application of the POD method can destroy the Hamiltonian structure in the reduced-order modelin this paper, we develop a new reduced-order modeling approach for Hamiltonian systems, which modifies the Galerkin projection-based POD -ROM so that the appropriate Hamiltonian structure is preserved. Since the POD truncation can degrade the approximation of the Hamiltonian function, we propose to use a POD basis from shifted snapshots to improve the approximation to the Hamiltonian function. We further derive a rigorous a priori error estimate for the structure-preserving ROM and demonstrate its effectiveness in several numerical examples. This approach can be readily extended to dissipative Hamiltonian systems, port-Hamiltonian systems, etc. Published by Elsevier B.V.
In this paper, we focus on the numerical solutions of Maxwell's equations with Dirichlet boundary conditions in rectangular coordinate. A class of explicit methods is derived by using an effective solver for a sys...
详细信息
In this paper, we focus on the numerical solutions of Maxwell's equations with Dirichlet boundary conditions in rectangular coordinate. A class of explicit methods is derived by using an effective solver for a system of ordinary differential equations which is obtained by approximating on spatial fields. A significant advantage of this class of methods is their simplicity and their ease of implementation. The error estimates presented in this paper show that the numerical solutions obtained by this class of methods is of high-order. The main advantage of this class of methods is that it is divergence-free.
In this paper the authors develop a new algorithm to solve the standard discrete-time algebraic Riccati equation by using a skew-Hamiltonian transformation and the square-root method. The algorithm is structure-preser...
详细信息
In this paper the authors develop a new algorithm to solve the standard discrete-time algebraic Riccati equation by using a skew-Hamiltonian transformation and the square-root method. The algorithm is structure-preserving and efficient because the Hamiltonian structure is fully exploited and only orthogonal transformations are used. The efficiency and stability of the algorithm are analyzed. Numerical examples are included.
Maxwell's equations describe the propagation of electromagnetic waves and are therefore fundamental to understanding many problems encountered in the study of antennas and electromagnetics. The aim of this paper i...
详细信息
Maxwell's equations describe the propagation of electromagnetic waves and are therefore fundamental to understanding many problems encountered in the study of antennas and electromagnetics. The aim of this paper is to propose and analyse an efficient fully discrete scheme for solving three-dimensional Maxwell's equations. This is accomplished by combining Fourier pseudospectral methods in space and exact formulation in time. Fast computation is efficiently implemented in the scheme by using the matrix diagonalisation method and fast Fourier transform algorithm which are well known in scientific computations. An optimal error estimate which is not encumbered by the CFL condition is established and the resulting scheme is proved to be of spectral accuracy in space and exact in time. Furthermore, the scheme is shown to have multiple conservation laws including discrete energy, helicity, momentum, symplecticity, and divergence-free field conservations. All the theoretical results of the accuracy and conservations are numerically illustrated by two numerical tests.
Nonlinear diffusion equations are now widely used to restore and enhance images. They allow to eliminate noise and artifacts while preserving large global features, such as object contours. In this context, we propose...
详细信息
Nonlinear diffusion equations are now widely used to restore and enhance images. They allow to eliminate noise and artifacts while preserving large global features, such as object contours. In this context, we propose a differential-geometric framework to define PDEs acting on some manifold constrained datasets. We consider the case of images taking value into matrix manifolds defined by orthogonal and spectral constraints. We directly incorporate the geometry and natural metric of the underlying configuration space ( viewed as a Lie group or a homogeneous space) in the design of the corresponding flows. Our numerical implementation relies on structure-preserving integrators that respect intrinsically the constraints geometry. The efficiency and versatility of this approach are illustrated through the anisotropic smoothing of diffusion tensor volumes in medical imaging.
Nonlinear diffusion equations are now widely used to restore and enhance images. They allow to eliminate noise and artifacts while preserving large global features, such as object contours. In this context, we propose...
详细信息
Nonlinear diffusion equations are now widely used to restore and enhance images. They allow to eliminate noise and artifacts while preserving large global features, such as object contours. In this context, we propose a differential-geometric framework to define PDEs acting on some manifold constrained datasets. We consider the case of images taking value into matrix manifolds defined by orthogonal and spectral constraints. We directly incorporate the geometry and natural metric of the underlying configuration space ( viewed as a Lie group or a homogeneous space) in the design of the corresponding flows. Our numerical implementation relies on structure-preserving integrators that respect intrinsically the constraints geometry. The efficiency and versatility of this approach are illustrated through the anisotropic smoothing of diffusion tensor volumes in medical imaging.
Generally, large-scale dynamical systems pose tremendous computational difficulties when applied in numerical simulations. In order to overcome these challenges we use several model reduction techniques. For stable li...
详细信息
Generally, large-scale dynamical systems pose tremendous computational difficulties when applied in numerical simulations. In order to overcome these challenges we use several model reduction techniques. For stable linear models these techniques work very well and provide good approximations for the full model. However, large-scale unstable systems arise in many applications. Many of the known model reduction methods are not very robust, or in some cases, may not even work if we are dealing with unstable systems. When approximating an unstable sytem by a reduced order model, accuracy is not the only concern. We also need to consider the structure of the reduced order model. Often, it is important that the number of unstable poles in the reduced system is the same as the number of unstable poles in the original system. The Iterative Rational Krylov Algorithm (IRKA) is a robust model reduction technique which is used to locally reduce stable linear dynamical systems optimally in the $mathcal{H}_2$-norm. While we cannot guarantee that IRKA reduces an unstable model optimally, there are no numerical obstacles to the reduction of an unstable model via IRKA. In this thesis, we investigate IRKA’s behavior when it is used to reduce unstable models. We also consider systems for which we cannot obtain a first order realization of the transfer function. We can use Realization-independent IRKA to obtain a reduced order model which does not preserve the structure of the original model. In this paper, we
implement a structurepreserving algorithm for systems with nonlinear frequency dependency.
In this paper, we propose a class of stochastic exponential discrete gradient schemes for SDEs with linear and gradient components in the coefficients. The root mean-square errors of the schemes are analyzed, and the ...
详细信息
In this paper, we propose a class of stochastic exponential discrete gradient schemes for SDEs with linear and gradient components in the coefficients. The root mean-square errors of the schemes are analyzed, and the structure-preserving properties of the schemes for SDEs with special structures are investigated. Numerical tests are performed to verify the theoretical results and illustrate the numerical behavior of the proposed methods. (C) 2021 Elsevier B.V. All rights reserved.
暂无评论