Abstract:In the implementation of spectral difference (SD) method, the conserved variables at the flux points are calculated from the solution points using extrapolation or interpolation schemes. The errors incurred in using extrapolation and interpolation would result in instability. On the other hand, the difference between the left and right conserved variables at the edge interface will introduce dissipation to the SD method when applying a Riemann solver to compute the flux at the element interface. In this paper, an optimization of the extrapolation and interpolation schemes for the fourth order SD method on quadrilateral element is carried out in the wavenumber space through minimizing their dispersion error over a selected band of wavenumbers. The optimized coefficients of the extrapolation and interpolation are presented. And the dispersion error of the original and optimized schemes is plotted and compared. An improvement of the dispersion error over the resolvable wavenumber range of SD method is obtained. The stability of the optimized fourth order SD scheme is analyzed. It is found that the stability of the 4th order scheme with Chebyshev–Gauss–Lobatto flux points, which is originally weakly unstable, has been improved through the optimization. The weak instability is eliminated completely if an additional second order filter is applied on selected flux points. One and two dimensional linear wave propagation analyses are carried out for the optimized scheme. It is found that in the resolvable wavenumber range the new SD scheme is less dispersive and less dissipative than the original scheme, and the new scheme is less anisotropic for 2D wave propagation. The optimized SD solver is validated with four computational aeroacoustics (CAA) workshop benchmark problems. The numerical results with optimized schemes agree much better with the analytical data than those with the original schemes.

Stabilizing the Consistent Quasidiffusion Method with Linear Prolongation

A "Consistent" Quasidiffusion Method for Solving Particle Transport Problems

An Optimized Dispersion-Relation-Preserving Combined Compact Difference Scheme to Solve Advection Equations

A Stable Condition and Adaptive Diffusion Coefficients for the Coarse-Mesh Finite Difference Method

A Weighted Least-Squares Transport Equation Compatible with Source Iteration and Voids

Numerical Study of Lattice Boltzmann Methods for a Convection–diffusion Equation Coupled with Navier–Stokes Equations

Convergence Analysis for the CMFD Accelerated Linear Axial Expansion Transport Method Based on Fourier Analysis

A Discontinuous Galerkin Finite Element Method Based Axial SN for the 2D/1D Transport Method

Nonlinear Diffusion Acceleration of the Least-Squares Transport Equation in Geometries with Voids

An Optimized Spectral Difference Scheme for CAA Problems

The Extended Diamond Difference - Constant Nodal Method with Decoupled Cell Iteration Scheme in Two-Dimensional Discrete Ordinate Transport Problems

Numerical analysis of a spherical harmonic discontinuous Galerkin method for scaled radiative transfer equations with isotropic scattering

The Improved Odcmfd Method Accelerated for Quasi-3D Method Based on Fourier Analysis

On the Accuracy of Transport Spatial Discretization Methods for Diffusive Regimes

Additive Angular Dependent Rebalance Acceleration Arithmetic for Neutron Transport Equation in Optically Thick Diffusive Region

The Characteristic Finite Difference Streamline Diffusion Method for Convection-Dominated Diffusion Problems

A Stabilized Spectral Element Method for Reaction-Convection-Diffusion Equation

A diffusion synthetic acceleration method for steady discrete unified gas kinetic scheme in radiative heat transfer

Numerical Recovery of the Diffusion Coefficient in Diffusion Equations from Terminal Measurement

Globally Conservative, Hybrid Self-Adjoint Angular Flux and Least-Squares Method Compatible with Void

An improved method for solving quasilinear convection diffusion problems on a coarse mesh