A Class of Diagonally Dominant Implicit Schemes with Arbitrary Numerical Dissipation

1980 ◽  
pp. 186-193
Author(s):  
Nimai Kumar Mitra ◽  
Martin Fiebig
Keyword(s):  
Numerical Dissipation ◽  
Implicit Schemes ◽  
Diagonally Dominant

scholarly journals Polyhedral approximations of the semidefinite cone and their application

2021 ◽  
Author(s):  
Yuzhu Wang ◽  
Akihiro Tanaka ◽  
Akiko Yoshise
Keyword(s):  
Initial Approximation ◽  
Semidefinite Relaxation ◽  
Stable Set ◽  
Cutting Plane Methods ◽  
Diagonally Dominant Matrices ◽  
Maximum Stable Set ◽  
Simple Expansion ◽  
Diagonally Dominant ◽  
Polyhedral Approximations ◽  
Semidefinite Cone

AbstractWe develop techniques to construct a series of sparse polyhedral approximations of the semidefinite cone. Motivated by the semidefinite (SD) bases proposed by Tanaka and Yoshise (Ann Oper Res 265:155–182, 2018), we propose a simple expansion of SD bases so as to keep the sparsity of the matrices composing it. We prove that the polyhedral approximation using our expanded SD bases contains the set of all diagonally dominant matrices and is contained in the set of all scaled diagonally dominant matrices. We also prove that the set of all scaled diagonally dominant matrices can be expressed using an infinite number of expanded SD bases. We use our approximations as the initial approximation in cutting plane methods for solving a semidefinite relaxation of the maximum stable set problem. It is found that the proposed methods with expanded SD bases are significantly more efficient than methods using other existing approximations or solving semidefinite relaxation problems directly.

scholarly journals An Improved Time Integration Method Based on Galerkin Weak Form with Controllable Numerical Dissipation

2021 ◽  
Vol 11 (4) ◽  
pp. 1932
Author(s):  
Weixuan Wang ◽  
Qinyan Xing ◽  
Qinghao Yang
Keyword(s):  
High Frequency ◽  
Integration Method ◽  
Time Integration ◽  
Low Frequency ◽  
Weak Form ◽  
High Accuracy ◽  
Numerical Dissipation ◽  
Frequency Range ◽  
Time Integration Method ◽  
Numerical Damping

Based on the newly proposed generalized Galerkin weak form (GGW) method, a two-step time integration method with controllable numerical dissipation is presented. In the first sub-step, the GGW method is used, and in the second sub-step, a new parameter is introduced by using the idea of a trapezoidal integral. According to the numerical analysis, it can be concluded that this method is unconditionally stable and its numerical damping is controllable with the change in introduced parameters. Compared with the GGW method, this two-step scheme avoids the fast numerical dissipation in a low-frequency range. To highlight the performance of the proposed method, some numerical problems are presented and illustrated which show that this method possesses superior accuracy, stability and efficiency compared with conventional trapezoidal rule, the Wilson method, and the Bathe method. High accuracy in a low-frequency range and controllable numerical dissipation in a high-frequency range are both the merits of the method.

scholarly journals A parameter-free total Lagrangian smooth particle hydrodynamics algorithm applied to problems with free surfaces

2021 ◽  
Author(s):  
Kenny W. Q. Low ◽  
Chun Hean Lee ◽  
Antonio J. Gil ◽  
Jibran Haider ◽  
Javier Bonet
Keyword(s):  
Ad Hoc ◽  
Wide Spectrum ◽  
Free Surface Flow ◽  
Surface Flow ◽  
Point Of View ◽  
Numerical Dissipation ◽  
Lagrangian Description ◽  
Smooth Particle Hydrodynamics ◽  
Particle Hydrodynamics ◽  
Sph Algorithm

AbstractThis paper presents a new Smooth Particle Hydrodynamics computational framework for the solution of inviscid free surface flow problems. The formulation is based on the Total Lagrangian description of a system of first-order conservation laws written in terms of the linear momentum and the Jacobian of the deformation. One of the aims of this paper is to explore the use of Total Lagrangian description in the case of large deformations but without topological changes. In this case, the evaluation of spatial integrals is carried out with respect to the initial undeformed configuration, yielding an extremely efficient formulation where the need for continuous particle neighbouring search is completely circumvented. To guarantee stability from the SPH discretisation point of view, consistently derived Riemann-based numerical dissipation is suitably introduced where global numerical entropy production is demonstrated via a novel technique in terms of the time rate of the Hamiltonian of the system. Since the kernel derivatives presented in this work are fixed in the reference configuration, the non-physical clumping mechanism is completely removed. To fulfil conservation of the global angular momentum, a posteriori (least-squares) projection procedure is introduced. Finally, a wide spectrum of dedicated prototype problems is thoroughly examined. Through these tests, the SPH methodology overcomes by construction a number of persistent numerical drawbacks (e.g. hour-glassing, pressure instability, global conservation and/or completeness issues) commonly found in SPH literature, without resorting to the use of any ad-hoc user-defined artificial stabilisation parameters. Crucially, the overall SPH algorithm yields equal second order of convergence for both velocities and pressure.

scholarly journals Visual Simulation of Detailed Turbulent Water by Preserving the Thin Sheets of Fluid

Symmetry ◽  
2018 ◽  
Vol 10 (10) ◽  
pp. 502 ◽  
Author(s):  
Jong-Hyun Kim ◽  
Wook Kim ◽  
Young Kim ◽  
Jung Lee
Keyword(s):  
Turbulent Flows ◽  
Visual Simulation ◽  
Numerical Dissipation ◽  
Small Scale ◽  
Thin Sheets ◽  
Global Flow ◽  
Low Mass ◽  
The Difference ◽  
Water Simulation ◽  
Turbulent Water

When we perform particle-based water simulation, water particles are often increased dramatically because of particle splitting around breaking holes to maintain the thin fluid sheets. Because most of the existing approaches do not consider the volume of the water particles, the water particles must have a very low mass to satisfy the law of the conservation of mass. This phenomenon smears the motion of the water, which would otherwise result in splashing, thereby resulting in artifacts such as numerical dissipation. Thus, we propose a new fluid-implicit, particle-based framework for maintaining and representing the thin sheets and turbulent flows of water. After splitting the water particles, the proposed method uses the ghost density and ghost mass to redistribute the difference in mass based on the volume of the water particles. Next, small-scale turbulent flows are formed in local regions and transferred in a smooth manner to the global flow field. Our results show us the turbulence details as well as the thin sheets of water, thereby obtaining an aesthetically pleasing improvement compared with existing methods.

scholarly journals Goal-Oriented Error Estimation for the Fractional Step Theta Scheme

2014 ◽  
Vol 14 (2) ◽  
pp. 203-230 ◽  
Author(s):  
Dominik Meidner ◽  
Thomas Richter
Keyword(s):  
Time Discretization ◽  
Posteriori Error ◽  
Numerical Dissipation ◽  
Error Estimator ◽  
Fractional Step ◽  
Weighted Residual Method ◽  
Time Stepping ◽  
Backward Euler ◽  
Galerkin Formulation ◽  
Theta Method

Abstract. In this work, we derive a goal-oriented a posteriori error estimator for the error due to time-discretization of nonlinear parabolic partial differential equations by the fractional step theta method. This time-stepping scheme is assembled by three steps of the general theta method, that also unifies simple schemes like forward and backward Euler as well as the Crank–Nicolson method. Further, by combining three substeps of the theta time-stepping scheme, the fractional step theta time-stepping scheme is derived. It possesses highly desired stability and numerical dissipation properties and is second order accurate. The derived error estimator is based on a Petrov–Galerkin formulation that is up to a numerical quadrature error equivalent to the theta time-stepping scheme. The error estimator is assembled as one weighted residual term given by the dual weighted residual method and one additional residual estimating the Galerkin error between time-stepping scheme and Petrov–Galerkin formulation.

Efficient semi-implicit schemes for stiff systems

2006 ◽  
Vol 214 (2) ◽  
pp. 521-537 ◽  
Author(s):  
Qing Nie ◽  
Yong-Tao Zhang ◽  
Rui Zhao
Keyword(s):  
Stiff Systems ◽  
Implicit Schemes

scholarly journals A Finite-Volume Icosahedral Shallow-Water Model on a Local Coordinate

2009 ◽  
Vol 137 (4) ◽  
pp. 1422-1437 ◽  
Author(s):  
Jin-Luen Lee ◽  
Alexander E. MacDonald
Keyword(s):  
Shallow Water ◽  
Finite Volume ◽  
Spherical Surface ◽  
Second Order ◽  
Water Model ◽  
Numerical Dissipation ◽  
Shallow Water Model ◽  
Finite Volume Model ◽  
Volume Model ◽  
Cartesian Coordinate

Abstract An icosahedral-hexagonal shallow-water model (SWM) on the sphere is formulated on a local Cartesian coordinate based on the general stereographic projection plane. It is discretized with the third-order Adam–Bashforth time-differencing scheme and the second-order finite-volume operators for spatial derivative terms. The finite-volume operators are applied to the model variables defined on the nonstaggered grid with the edge variables interpolated using polynomial interpolation. The projected local coordinate reduces the solution space from the three-dimensional, curved, spherical surface to the two-dimensional plane and thus reduces the number of complete sets of basis functions in the Vandermonde matrix, which is the essential component of the interpolation. The use of a local Cartesian coordinate also greatly simplifies the mathematic formulation of the finite-volume operators and leads to the finite-volume integration along straight lines on the plane, rather than along curved lines on the spherical surface. The SWM is evaluated with the standard test cases of Williamson et al. Numerical results show that the icosahedral SWM is free from Pole problems. The SWM is a second-order finite-volume model as shown by the truncation error convergence test. The lee-wave numerical solutions are compared and found to be very similar to the solutions shown in other SWMs. The SWM is stably integrated for several weeks without numerical dissipation using the wavenumber 4 Rossby–Haurwitz solution as an initial condition. It is also shown that the icosahedral SWM achieves mass conservation within round-off errors as one would expect from a finite-volume model.

Dissipation Improvement of MUSCL Scheme for Computational Aeroacoustics

2001 ◽  
Vol 17 (1) ◽  
pp. 39-47
Author(s):  
San-Yin Lin ◽  
Sheng-Chang Shih ◽  
Jen-Jiun Hu
Keyword(s):  
Wave Interaction ◽  
Sine Wave ◽  
Numerical Dissipation ◽  
Finite Volume Scheme ◽  
Linear Wave ◽  
Riemann Solver ◽  
Two Dimensional ◽  
Order Accuracy ◽  
Upwind Schemes ◽  
Continuous Waves

ABSTRACTAn upwind finite-volume scheme is studied for solving the solutions of two dimensional Euler equations. It based on the MUSCL (Monotone Upstream Scheme for Conservation Laws) approach with the Roe approximate Riemann solver for the numerical flux evaluation. First, dissipation and dispersion relation, and group velocity of the scheme are derived to analyze the capability of the proposed scheme for capturing physical waves, such as acoustic, entropy, and vorticity waves. Then the scheme is greatly enhanced through a strategy on the numerical dissipation to effectively handle aeroacoustic computations. The numerical results indicate that the numerical dissipation strategy allows that the scheme simulates the continuous waves, such as sound and sine waves, at fourth-order accuracy and captures the discontinuous waves, such a shock wave, sharply as well as most of upwind schemes do. The tested problems include linear wave convection, propagation of a sine-wave packet, propagation of discontinuous and sine waves, shock and sine wave interaction, propagation of acoustic, vorticity, and density pulses in an uniform freestream, and two-dimensional traveling vortex in a low-speed freestream.

Sign in / Sign up

Export Citation Format

Share Document