Higher-order Approximations Research Articles

A local projection stabilised HHO method for the Oseen problem

In this article, we consider a local projection stabilisation for a Hybrid High-Order (HHO) approximation of the Oseen problem. We prove an existence-uniqueness result under a stronger SUPG-like norm. We improve the stability and provide error estimation in stronger norm for convection dominated Oseen problem. We also derive an optimal order error estimate under the SUPG-like norm for equal-order polynomial discretisation of velocity and pressure spaces. Numerical experiments are performed to validate the theoretical results.

Curved-pipe flow with boundary conditions involving Bernoulli pressure

In this article, we study the steady-state flow of the incompressible viscous fluid through a thin distorted pipe with an arbitrary central curve. We prescribe the inflow and outflow boundary conditions involving the Bernoulli pressure with a given pressure drop. Using the multiscale expansion technique with respect to the pipe's thickness, we construct the higher-order asymptotic approximation of the flow given by the explicit formulae for the velocity and pressure. We also perform a detailed error analysis justifying the usage of the proposed solution and indicating its order of accuracy. For more information see https://ejde.math.txstate.edu/Volumes/2024/63/abstr.html

First Derivative Approximations and Applications

In this paper, we consider constructions of first derivative approximations using the generating function. The weights of the approximations contain the powers of a parameter whose modulus is less than one. The values of the initial weights are determined, and the convergence and order of the approximations are proved. The paper discusses applications of approximations of the first derivative for the numerical solution of ordinary and partial differential equations and proposes an algorithm for fast computation of the numerical solution. Proofs of the convergence and accuracy of the numerical solutions are presented and the performance of the numerical methods considered is compared with the Euler method. The main goal of constructing approximations for integer-order derivatives of this type is their application in deriving high-order approximations for fractional derivatives, whose weights have specific properties. The paper proposes the construction of an approximation for the fractional derivative and its application for numerically solving fractional differential equations. The theoretical results for the accuracy and order of the numerical methods are confirmed by the experimental results presented in the paper.

Topology optimization based on the improved high-order RAMP interpolation model and bending properties research for curved square honeycomb sandwich structures

Curved honeycomb sandwich structures with larger surface areas effectively reduce the number of fasteners and connectors, resulting in weight reduction, cost savings, and reliability improvement. Square honeycombs exhibit higher in-plane tensile strength, and are more compatible with mechanical components. Topology optimization can realize the variable-density design of square honeycombs, enhancing the strength and stiffness of curved sandwich structures, but the high-order Rational Approximation of Material Properties (RAMP) interpolation model with a fast convergence rate has the relatively poor clarity and stability in topology boundaries. Hence, a variable-density topology optimization method based on the improved high-order RAMP model is developed for curved square honeycomb sandwich structures. The high-order RAMP interpolation model is improved by incorporating a minimum modulus term into the material interpolation function and employing the Bi-directional Evolutionary Structural Optimization (BESO) to refine the Optimality Criteria (OC) method. A functional relationship between the wall thickness of cross-shaped cells (i.e., simplified models of four adjacent square cells arranged in a cross) and the relative density of topology units is constructed using a density mapping method, and the topology optimization with the relative density of cross-shaped cells as the design variable is performed to minimize the compliance of the core. Three-point bending tests are conducted on 3D printed non-optimized and optimized curved honeycomb sandwich structures with different core material retention rates and panel thicknesses, and the experimental results are compared with those of simulations to explore the bending properties and failure behaviors. The results indicate that the modified high-order RAMP interpolation model enhances the clarity and stability of topology structures, the variable-density topology optimization significantly improves the bending properties of curved square honeycomb sandwich structures, and the experimental and numerical results are largely consistent.

Advancing instantaneous detonation prediction of condensed energetic materials based on high-order compressible multiphase fluid dynamics

The instantaneous detonation model (IDM) is widely used in simulating underwater explosions due to its efficiency and ability to ignore the detonation reaction process. In this study, we propose a new IDM to predict the fluid structure in the detonation zone of an RS211 explosive charge. This model is based on high-order solutions provided by the detonation shock dynamics model, where the spatial term is discretized using fifth-order WENO reconstruction in characteristic space and Lax–Friedrichs’s splitting and the temporal terms are discretized using a third-order TVD Runge–Kutta scheme. The interface motion is captured using the level-set method combined with MGFM, and a programmed burn model is provided to describe the generation and propagation of the detonation wave. The self-similarity of detonation wave propagation is validated, and the quantitative calculation formula of the instantaneous detonation model is obtained by averaging or curve fitting the dimensionless results. Consequently, the IDM of the RS211 charge is established using high-order polynomial approximations of the Taylor rarefaction zone and a constant static zone for 1D planar, cylindrical, and spherical RS211 charges. The application of the IDM involves direct mapping from the radial direction to the spatial structured grid for 1D planar, 2D cylindrical, and 3D spherical charges. Numerical results demonstrate that the IDM proposed in this paper shows good accuracy and high computational efficiency.

On a parabolic Monge–Ampère type equation on compact almost Hermitian manifolds

AbstractWe investigate a parabolic Monge–Ampère type equation on compact almost Hermitian manifolds and derive a priori gradient and second‐order derivative estimates for solutions to this parabolic equation. These a priori estimates give us higher order estimates and a long‐time solution. Then, we can observe its behavior as .

Vibrational and stability analysis of planar double pendulum dynamics near resonance

AbstractThe focus of this paper is to examine the motion of a novel double pendulum (DP) system with two degrees of freedom (DOF). This system operates under specific constraints to follow a Lissajous curve, with its pivot point moving along this path in a plane. The nonlinear differential equations governing this system are derived using Lagrange's equations. Their analytical solutions (AS) are subsequently calculated using the multiple-scales method (MSM), which provides higher-order approximations. These solutions are considered new, as the traditional MSM has been applied to this novel system for the first time. Additionally, the accuracy of these solutions is validated through numerical results obtained using the fourth-order Runge–Kutta method. The solvability conditions and characteristic exponents are determined based on resonance cases. The Routh–Hurwitz criteria (RHC) are employed to assess the stability of the fixed points corresponding to the steady-state solutions. They are also used to demonstrate the frequency response curves. The nonlinear stability analysis is performed by examining the stability and instability ranges. Resonance curves and time history plots are presented to analyze the behavior of the system for specific parameter values. The investigation delves into a comprehensive analysis of bifurcation diagrams (BDs) and Lyapunov exponent spectra (LEs), aiming to uncover the various types of motion present within the system. Systematic examination of these charts reveals critical insights into transitions between stable, quasi-stable, and chaotic dynamical behaviors. This work has practical applications in various fields, such as robotics, pump compressors, rotor dynamics, and transportation devices. It can be used to study the vibrational motion of these systems.

High-order approximation of Caputo–Prabhakar derivative with applications to linear and nonlinear fractional diffusion models

Abstract In this study, we devise a high-order numerical scheme to approximate the Caputo–Prabhakar derivative of order α ∈ (0, 1), using an rth-order time stepping Lagrange interpolation polynomial, where 3 ≤ r ∈ N $3\le r\in \mathbb{N}$ . The devised scheme is a generalization of the existing schemes developed earlier. Further, we adopt the discussed scheme for solving a linear time fractional advection–diffusion equation and a nonlinear time fractional reaction–diffusion equation with Dirichlet type boundary conditions. We show that the discussed method is unconditionally stable, uniquely solvable and convergent with convergence order O(τ r+1−α , h 2), where τ and h are the temporal and spatial step sizes, respectively. Without loss of generality, applicability of the discussed method is established by illustrative examples for r = 4, 5.

Graph Neural Networks With High-Order Polynomial Spectral Filters.

General graph neural networks (GNNs) implement convolution operations on graphs based on polynomial spectral filters. Existing filters with high-order polynomial approximations can detect more structural information when reaching high-order neighborhoods but produce indistinguishable representations of nodes, which indicates their inefficiency of processing information in high-order neighborhoods, resulting in performance degradation. In this article, we theoretically identify the feasibility of avoiding this problem and attribute it to overfitting polynomial coefficients. To cope with it, the coefficients are restricted in two steps, dimensionality reduction of the coefficients' domain and sequential assignment of the forgetting factor. We transform the optimization of coefficients to the tuning of a hyperparameter and propose a flexible spectral-domain graph filter, which significantly reduces the memory demand and the adverse impacts on message transmission under large receptive fields. Utilizing our filter, the performance of GNNs is improved significantly in large receptive fields and the receptive fields of GNNs are multiplied as well. Meanwhile, the superiority of applying a high-order approximation is verified across various datasets, notably in strongly hyperbolic datasets. Codes are publicly available at: https://github.com/cengzeyuan/TNNLS-FFKSF.

Convergence and high order of approximation by Steklov sampling operators

In this paper we introduce a new class of sampling-type operators, named Steklov sampling operators. The idea is to consider a sampling series based on a kernel function that is a discrete approximate identity, and which constitutes a reconstruction process of a given signal f, based on a family of sample values which are Steklov integrals of order r evaluated at the nodes k/w, k∈Z\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$k \\in {\\mathbb {Z}}$$\\end{document}, w>0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$w>0$$\\end{document}. The convergence properties of the introduced sampling operators in continuous functions spaces and in the Lp\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L^p$$\\end{document}-setting have been studied. Moreover, the main properties of the Steklov-type functions have been exploited in order to establish results concerning the high order of approximation. Such results have been obtained in a quantitative version thanks to the use of the well-known modulus of smoothness of the approximated functions, and assuming suitable Strang-Fix type conditions, which are very typical assumptions in applications involving Fourier and Harmonic analysis. Concerning the quantitative estimates, we proposed two different approaches; the first one holds in the case of Steklov sampling operators defined with kernels with compact support, its proof is substantially based on the application of the generalized Minkowski inequality, and it is valid with respect to the p-norm, with 1≤p≤+∞\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$1 \\le p \\le +\\infty $$\\end{document}. In the second case, the restriction on the support of the kernel is removed and the corresponding estimates are valid only for 1<p≤+∞\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$1 < p\\le +\\infty $$\\end{document}. Here, the key point of the proof is the application of the well-known Hardy–Littlewood maximal inequality. Finally, a deep comparison between the proposed Steklov sampling series and the already existing sampling-type operators has been given, in order to show the effectiveness of the proposed constructive method of approximation. Examples of kernel functions satisfying the required assumptions have been provided.

On the Approximation of the Hardy Z-Function via High-Order Sections

The Z-function is the real function given by Z(t)=eiθ(t)ζ12+it, where ζ(s) is the Riemann zeta function, and θ(t) is the Riemann–Siegel theta function. The function, central to the study of the Riemann hypothesis (RH), has traditionally posed significant computational challenges. This research addresses these challenges by exploring new methods for approximating Z(t) and its zeros. The sections of Z(t) are given by ZN(t):=∑k=1Ncos(θ(t)−ln(k)t)k for any N∈N. Classically, these sections approximate the Z-function via the Hardy–Littlewood approximate functional equation (AFE) Z(t)≈2ZN˜(t)(t) for N˜(t)=t2π. While historically important, the Hardy–Littlewood AFE does not sufficiently discern the RH and requires further evaluation of the Riemann–Siegel formula. An alternative, less common, is Z(t)≈ZN(t)(t) for N(t)=t2, which is Spira’s approximation using higher-order sections. Spira conjectured, based on experimental observations, that this approximation satisfies the RH in the sense that all of its zeros are real. We present a proof of Spira’s conjecture using a new approximate equation with exponentially decaying error, recently developed by us via new techniques of acceleration of series. This establishes that higher-order approximations do not need further Riemann–Siegel type corrections, as in the classical case, enabling new theoretical methods for studying the zeros of zeta beyond numerics.

High-Order Approximation to Caputo Derivative on Graded Mesh and Time-Fractional Diffusion Equation for Nonsmooth Solutions

Abstract In this paper, a high-order approximation to Caputo-type time-fractional diffusion equations (TFDEs) involving an initial-time singularity of the solution is proposed. At first, we employ a numerical algorithm based on the Lagrange polynomial interpolation to approximate the Caputo derivative on the nonuniform mesh. The truncation error rate and the optimal grading constant of the approximation on a graded mesh are obtained as min{4−α,rα} and (4−α)/α, respectively, where α∈(0,1) is the order of fractional derivative and r≥1 is the mesh grading parameter. Using this new approximation, a difference scheme for the Caputo-type time-fractional diffusion equation on the graded temporal mesh is formulated. The scheme proves to be uniquely solvable for general r. Then, we derive the unconditional stability of the scheme on uniform mesh. The convergence of the scheme, in particular for r = 1, is analyzed for nonsmooth solutions and concluded for smooth solutions. Finally, the accuracy of the scheme is verified by analyzing the error through a few numerical examples.

Fast implicit update schemes for Cahn–Hilliard-type gradient flow in the context of Fourier-spectral methods

This work discusses a way of allowing fast implicit update schemes for the temporal discretization of phase-field models for gradient flow problems that employ Fourier-spectral methods for their spatial discretization. Through the repeated application of the Sherman–Morrison formula we provide a rule for approximations of the inverted tangent matrix of the corresponding Newton–Raphson method with a selectable order. Since the representation of this inversion is exact for a sufficiently high approximation order, the proposed scheme is shown to provide a fixed-point-type iterative solver for gradient flow problems that require the solution of linear systems in the context of an implicit time-integration. While the proposed scheme is applicable to general gradient flow phase-field models, we discuss the scheme in the context of the Cahn–Hilliard equation, the Swift–Hohenberg equation, and the phase-field crystal equation for which we demonstrate the performance of the proposed method in comparison with classical solvers.

Subharmonic response suppression of a quasi-zero stiffness system

Quasi-zero stiffness (QZS) isolators can effectively suppress low-frequency vibration due to high-static-low-dynamic stiffness (HSLDS) characteristics. The existence of nonlinearity induces intricate nonlinear phenomena in QZS isolators, particularly under conditions of small damping and large excitation. These nonlinear characteristics can seriously affect the vibration isolation performance, which has rarely been focused on in existing studies. This study examines a typical QZS isolator with nonlinear damping. The modified incremental harmonic balance (IHB) method is applied to directly analyze its dynamic characteristics and investigate the suppression of the subharmonic response. The dynamic analysis results indicate multiple subharmonic vibration intervals, including the 1/2, the 1/3, and the 1/9 subharmonic vibration. The basin of attraction at different excitation frequencies shows that subharmonic vibration is highly likely to occur. The effect of different damping on the transmissibility proves to be different. Simultaneously, rational parameter matching can guarantee the absence of subharmonic vibration while maintaining the primary resonance frequency and amplitude unchanged. The research results of this study provide theoretical support for the parameter selection of QZS isolators. In addition, the higher-order approximations of most isolators have the same form, so these results can guide other QZS isolators in suppressing subharmonic vibrations and improving their operating range.

Simulation of coupled groundwater flow and contaminant transport using quintic B-spline collocation method

PurposeThe purpose of this study is to introduce a new numerical scheme with no stability condition and high-order accuracy for the solution of two-dimensional coupled groundwater flow and transport simulation problems with regular and irregular geometries and compare the results with widely acceptable programs such as Modular Three-Dimensional Finite-Difference Ground-Water Flow Model (MODFLOW) and Modular Three-Dimensional Multispecies Transport Model (MT3DMS).Design/methodology/approachThe newly proposed numerical scheme is based on the method of lines (MOL) approach and uses high-order approximations both in space and time. Quintic B-spline (QBS) functions are used in space to transform partial differential equations, representing the relevant physical phenomena in the system of ordinary differential equations. Then this system is solved with the DOPRI5 algorithm that requires no stability condition. The obtained results are compared with the results of the MODFLOW and MT3DMS programs to verify the accuracy of the proposed scheme.FindingsThe results indicate that the proposed numerical scheme can successfully simulate the two-dimensional coupled groundwater flow and transport problems with complex geometry and parameter structures. All the results are in good agreement with the reference solutions.Originality/valueTo the best of the authors' knowledge, the QBS-DOPRI5 method is used for the first time for solving two-dimensional coupled groundwater flow and transport problems with complex geometries and can be extended to high-dimensional problems. In the future, considering the success of the proposed numerical scheme, it can be used successfully for the identification of groundwater contaminant source characteristics.

High-Order Numerical Approximation for 2D Time-Fractional Advection–Diffusion Equation under Caputo Derivative

High-Order Numerical Approximation for 2D Time-Fractional Advection–Diffusion Equation under Caputo Derivative

Holimap: an accurate and efficient method for solving stochastic gene network dynamics

Gene-gene interactions are crucial to the control of sub-cellular processes but our understanding of their stochastic dynamics is hindered by the lack of simulation methods that can accurately and efficiently predict how the distributions of gene product numbers vary across parameter space. To overcome these difficulties, here we present Holimap (high-order linear-mapping approximation), an approach that approximates the protein or mRNA number distributions of a complex gene regulatory network by the distributions of a much simpler reaction system. We demonstrate Holimap’s computational advantages over conventional methods by applying it to predict the stochastic time-dependent dynamics of various gene networks, including transcriptional networks ranging from simple autoregulatory loops to complex randomly connected networks, post-transcriptional networks, and post-translational networks. Holimap is ideally suited to study how the intricate network of gene-gene interactions results in precise coordination and control of gene expression.

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

Analytical solutions for the Boltzmann transport equation are typically limited to specific scenarios, prompting the research transport community to extensively explore numerical approaches with special attention to the discrete ordinates (SN) formulation. The spectral nodal methods are widely recognized for their capability to preserve the analytical solution of the SN mathematical model. However, these formulations exhibit notable complexity with high-order polynomial approximations for the transverse leakage terms, motivating researchers to pursue simpler hybrid approaches by partially preserving the analytical solution. Both spectral nodal and hybrid spectral nodal methods introduce nonstandard auxiliary equations not amenable to standard iterative schemes due to a highly coupling in angular directions. Following these ideas, we present a finite element – spectral nodal hybrid method referred to as Extended Diamond Difference-Constant Nodal (EDD-CN) method for two-dimensional fixed-source discrete ordinate transport problems with constant nodal approximation for the leakage terms. The proposed method seeks to preserve eigenfunctions associated with the dominant pair of eigenvalues in the transverse integrated SN equations. In addition, we propose a novel iterative scheme called Decoupled Cell Iteration (DCI) scheme which decouples the angular directions in the nonstandard auxiliary equations and follows conventional iteration transport sweeps for solving SN problems.

Increasingly high-order hybrid multi-resolution WENO schemes in multi-dimensions

In this paper, a new type of increasingly high-order hybrid multi-resolution weighted essentially non-oscillatory (HMR-WENO) schemes is presented in the finite difference framework for solving hyperbolic conservation laws in one, two, and three dimensions. Based on the reconstruction polynomials defined on the one-point, three-point, five-point, seven-point, and nine-point spatial stencils, we reconstruct one zeroth degree reconstruction polynomial, one quadratic reconstruction polynomial, one quartic reconstruction polynomial, one sextic reconstruction polynomial, one octave reconstruction polynomial, and their derivative polynomials together with a new hierarchical bisection method to design a series of new troubled cell indicators which can precisely find all extreme points of associated unequal degree reconstruction polynomials located inside the smallest interval in one dimension. The new troubled cell indicators do not introduce any manual parameters related to different problems. The new hybrid methodology is divided into two parts: if all extreme points of the reconstruction polynomials are nonexistent or outside the smallest interval, the target cell is not a troubled cell and the simple linear upwind schemes are utilized to obtain high-order approximations. Otherwise, the target cell is a troubled cell and the MR-WENO spatial reconstruction procedures with excellent shock-capture ability are adopted. Then a series of HMR-WENO schemes are proposed by using these new troubled cell indicators, which can be easily expanded to arbitrarily high-order accuracies in multi-dimensions. The main benefits of these HMR-WENO schemes are their efficiency, since they could save about 22%-75% CPU time than that of the same order MR-WENO schemes when simulating some benchmark examples in multi-dimensions.

An Unfitted Finite Element Poisson-Boltzmann Solver with Automatic Resolving of Curved Molecular Surface.

So far, the existing Poisson-Boltzmann (PB) solvers that accurately take into account the interface jump conditions need a pregenerated body-fitted mesh (molecular surface mesh). However, qualified biomolecular surface meshing and its implementation into numerical methods remains a challenging and laborious issue, which practically hinders the progress of further developments and applications of a bunch of numerical methods in this field. In addition, even with a molecular surface mesh, it is only a low-order approximation of the original curved surface. In this article, an interface-penalty finite element method (IPFEM), which is a typical unfitted finite element method, is proposed to solve the Poisson-Boltzmann equation (PBE) without requiring the user to generate a molecular surface mesh. The Gaussian molecular surface is used to represent the molecular surface and can be automatically resolved with a high-order approximation within our method. Theoretical convergence rates of the IPFEM for the linear PB equation have been provided and are well validated on a benchmark problem with an analytical solution (we also noticed from numerical examples that the IPFEM has similar convergence rates for the nonlinear PBE). Numerical results on a set of different-sized biomolecules demonstrate that the IPFEM is numerically stable and accurate in the calculation of biomolecular electrostatic solvation energy.