Regular ArticleDispersion-Relation-Preserving Finite Difference Schemes for Computational Acoustics
Abstract
Acoustics problems are governed by the linearized Euler equations. According to wave propagation theory, the number of wave modes and their wave propagation characteristics are all encoded in the dispersion relations of the governing equations. Thus one is assured that the numerical solutions of a high order finite difference scheme will have the same number of wave modes (namely, the acoustic, vorticity, and entropy waves), the same wave propagation characteristics (namely, nondispersive, nondissipative, and isotropic) and the same wave speeds as those of the solutions of the Euler equations if both systems of equations have the same dispersion relations. Finite difference schemes which have the same dispersion relations as the original partial differential equations are referred to as dispersion-relation-preserving (DRP) schemes. A way to construct time marching DRP schemes by optimizing the finite difference approximations of the space and time derivatives in the wave number and frequency space is proposed. The stability of these schemes is analyzed and a sufficient condition for numerical stability is established. A set of radiation and outflow boundary conditions compatible with the DRP schemes is constructed. These conditions are derived from the asymptotic solutions of the governing equations. The asymptotic solutions are found by the use of Fourier-Laplace transforms and the method of stationary phase. A sequence of numerical simulations has been carried out. These simulations are designed to test the effectiveness of the DRP schemes and the radiation and outflow boundary conditions. The computed solutions agree very favorably with the exact solutions. The radiation boundary conditions perform satisfactorily causing little acoustic reflections. The outflow boundary conditions are found to be quite transparent to outgoing disturbances even when the disturbances are made up of a combination of acoustic, vorticity, and entropy waves.
References (0)
Cited by (2200)
Nonlinear acoustic radiation induced by in-plane vibration of hyperelastic rubber-like plates subject to dynamic loads
2024, Wave MotionThis work is concerned with numerical studies on nonlinear vibration and acoustic radiation behaviors of hyperelastic plates made of rubber material. Considering both the geometric and material nonlinearity of the rubber material, structural model of the hyperelastic plate is developed based on the nonlinear finite element method and the Mooney-Rivlin constitutive model. Acoustic waves in an inviscid and compressible fluid perturbed by the in-plane vibrations of the hyperelastic plate are governed by a set of first-order linearized partial differential equations. A fourth-order dispersion-relation-preserving (DRP) finite difference scheme is utilized to compute numerical solutions of the acoustic pressure responses. The structural-acoustic interface between the hyperelastic plate and exterior fluid is constructed through a robust ghost-cell sharp-interface immersed boundary method (IBM) of high numerical stability such that the compatibility conditions on the interface can be satisfied implicitly, although the structural Lagrangian meshes of the rubber plate and the fluid Eulerian grids are not matched. Several numerical examples are designed to check the effectiveness, convergence, and availability of the numerical structural-acoustic coupling model. Based on the numerical model, nonlinear vibro-acoustics response behaviors of the hyperelastic rubber plate subject to uniform dynamic loads are analyzed. Relevance between the in-plane vibrations of the plate and the spatial distribution of the acoustic pressure in exterior fluid is revealed. Effects of the excitation frequency and amplitude of the external loads on the nonlinear vibro-acoustics behaviors of the hyperelastic plate are discussed. Radiation directivity patterns of the fundamental and high-order components show significant differences, and may change dramatically against the excitation frequency and amplitude. In addition, the influences of the geometric and material nonlinearity of the hyperelastic plate on the higher-order vibro-acoustics responses are evaluated.
Machine learning optimization of compact finite volume methods on unstructured grids
2024, Journal of Computational PhysicsFourier analysis based optimization techniques have been developed for numerical schemes on structured grids and result in significant performance improvements. However, there lacks an effective optimization technique for numerical schemes on general unstructured grids, due to the mesh geometry complexity that makes Fourier transformation infeasible. To address this issue, an optimization framework based on machine learning is developed in this paper for numerical schemes on unstructured grids. The optimization of a compact high-order variational reconstruction on triangular grids is performed to illustrate the framework. An artificial neural network (ANN) is used to predict the optimal values of the derivative weights on cell interfaces that are the free parameters of the variational reconstruction, given the local geometry as input. The ANN is trained by minimizing the solution errors of a linear advection equation on a set of generated compact stencils, with a set of sampled sine waves as initial conditions. The developed optimization framework is applicable to general numerical schemes on unstructured grids as the training does not require explicit computation of dispersion or dissipation. Numerical results for inviscid flow problems show that the optimized variational finite volume scheme is remarkably more accurate and efficient than its non-optimized counterpart, demonstrating the effectiveness of the machine learning optimization.
Analysis of spectral properties of WENO schemes in smooth regions
2024, Communications in Nonlinear Science and Numerical SimulationWeighted essentially non-oscillatory (WENO) schemes are a class of numerical schemes designed for solving problems with discontinuities. The spectral properties of them are hard to analyze because of their nonlinearities. A method of analyzing the average scheme in effect is proposed, by which we found that WENO schemes with a special kind of smoothness indicators can even obtain better spectral properties than their underlying linear schemes. This facilitates the design of WENO schemes with good spectral properties. The formulas of these smoothness indicators on sub-stencils have the same form except for different subscripts, which makes the nonlinear weights slightly close to the mean value. With this weight distribution, the resultant schemes behave like optimized upwind schemes with low dissipation and optimized dispersion properties. However, these WENO schemes may suffer from anti-dissipation and accuracy order degradation. Hence, they are suggested to be combined with the mapped WENO or the Z-type WENO methodology to address these problems. Some WENO schemes are taken as examples to illustrate the proposed methodology, and the results are in good agreement with those from the approximated dispersion relation (ADR) method.
Dispersion analysis of SPH as a way to understand its order of approximation
2024, Journal of Computational and Applied MathematicsSmoothed Particle Hydrodynamics (SPH) is a numerical method to solve dynamical partial differential equations (PDE). The basis of the method is a ¡¡kernel-based¿¿ way to compute the spatial derivatives of a function whose values are given in moving irregularly located nodes (Lagrangian particles). Accuracy of the SPH is determined by independent parameters — the shape of the kernel, the kernel size , the distance between the particles . Constructing high-order SPH-schemes for different types of PDE is a state-of-the-art problem of computational mathematics.
For the classical SPH-approximation of one-dimensional hyperbolic equations (isothermal gas dynamics) we found that the order of approximation of smooth solution correlates to the dispersion properties of the method. To this end we analyzed the dispersion relation for the approximation and found analytical representation of the numerical wave phase velocity. Moreover, for the first time, the order of approximation with respect to was confirmed in computational experiments on a dynamic problem of sound wave propagation. For two kernels with 2 and 4 continuum derivatives, the second and the fourth order of approximation, respectively, was found. This finding may be generalized as follows. The solution error in the one-dimensional case for a quasi-uniformly located particles has the form , where is a parameter determined by the shape of kernel (its smoothness, i.e. the number of continuum derivatives), is a parameter that does not depend on the shape of kernel (for classical non-negative kernels ), is the wavelength.
Our results indicates that to develop high-order SPH-schemes for hyperbolic equations besides improving the order of approximation with respect to one need to ensure the order of approximation with respect to . To this end kernels of which smoothness is at least 4 are necessary.
A stable staggered-grid finite-difference scheme for acoustic modeling beyond conventional stability limit
2024, Petroleum ScienceStaggered-grid finite-difference (SGFD) schemes have been widely used in acoustic wave modeling for geophysical problems. Many improved methods are proposed to enhance the accuracy of numerical modeling. However, these methods are inevitably limited by the maximum Courant-Friedrichs-Lewy (CFL) numbers, making them unstable when modeling with large time sampling intervals or small grid spacings. To solve this problem, we extend a stable SGFD scheme by controlling SGFD dispersion relations and maximizing the maximum CFL numbers. First, to improve modeling stability, we minimize the error between the FD dispersion relation and the exact relation in the given wave-number region, and make the FD dispersion approach a given function outside the given wave-number area, thus breaking the conventional limits of the maximum CFL number. Second, to obtain high modeling accuracy, we use the SGFD scheme based on the Remez algorithm to compute the FD coefficients. In addition, the hybrid absorbing boundary condition is adopted to suppress boundary reflections and we find a suitable weighting coefficient for the proposed scheme. Theoretical derivation and numerical modeling demonstrate that the proposed scheme can maintain high accuracy in the modeling process and the value of the maximum CFL number of the proposed scheme can exceed that of the conventional SGFD scheme when adopting a small maximum effective wavenumber, indicating that the proposed scheme improves stability during the modeling.
High fidelity simulation of dense vapours with thermodynamic consistent modelling
2024, Computers and FluidsThe paper describes the implementation of thermodynamic modelling within a high-order finite difference scheme. A new method to ensure thermodynamic consistency for Navier–Stokes simulations of real-gas1 flows is outlined. The paper shows that this method is suited to high fidelity simulations of dense-vapour flows with complex equations of state, such as found in Organic Rankine Cycle turbines. The method is illustrated with two examples of siloxane MM flow in a channel and in a turbine vane. The results show that thermodynamic consistency plays a vital role in the accurate prediction of irreversibility and therefore turbine loss.