Hexagonal Grid Computation of the Derivatives of the Solution to the Heat Equation by Using Fourth-Order Accurate Two-Stage Implicit Methods

We give fourth-order accurate implicit methods for the computation of the first-order spatial derivatives and second-order mixed derivatives involving the time derivative of the solution of first type boundary value problem of two dimensional heat equation. The methods are constructed based on two stages: At the first stage of the methods, the solution and its derivative with respect to time variable are approximated by using the implicit scheme in Buranay and Arshad in 2020. Therefore, Oh4+τ of convergence on constructed hexagonal grids is obtained that the step sizes in the space variables x1, x2 and in time variable are indicated by h, 32h and τ, respectively. Special difference boundary value problems on hexagonal grids are constructed at the second stages to approximate the first order spatial derivatives and the second order mixed derivatives of the solution. Further, Oh4+τ order of uniform convergence of these schemes are shown for r=ωτh2≥116,ω>0. Additionally, the methods are applied on two sample problems.

Reconstructed Interpolating Differential Operator Method with Arbitrary Order of Accuracy for the Hyperbolic Equation

We propose a family of multi-moment methods with arbitrary orders of accuracy for the hyperbolic equation via the reconstructed interpolating differential operator (RDO) approach. Reconstruction up to arbitrary order can be achieved on a single cell from properly allocated model variables including spatial derivatives of varying orders. Then we calculate the temporal derivatives of coefficients of the reconstructed polynomial and transform them into the temporal derivatives of the model variables. Unlike the conventional multi-moment methods which evolve different types of moments by deriving different equations, RDO can update all derivatives uniformly via a simple linear transform more efficiently. Based on difference in introducing interaction from adjacent cells, the central RDO and the upwind RDO are proposed. Both schemes enjoy high-order accuracy which is verified by Fourier analysis and numerical experiments.

A higher-order numerical analysis to study the flow physics and to optimize the design of a short-dwell blade coaters for higher efficiency

Blade coaters are most commonly used for coating of paper and paperboard with higher efficiency. The efficiency of short-dwell blades coaters depends on many factors such as the properties of the coating material, design of the coating reservoir, the types of flow behaviour taking place inside the reservoir, etc. In this work, we have proposed an optimal design of the reservoir to improve the efficiency of short-dwell coaters. The reservoir has been modeled as flow inside a two-dimensional rectangular cavity. Incompressible Navier-Stokes equations in primitive variable formulation have been solved to obtain the flow fields inside the cavity. Spatial derivatives present in the momentum, and continuity equations are evaluated using a sixth-order accurate compact scheme whereas the temporal derivatives are calculated using the fourth-order Runge-Kutta method. The actual rate of convergence of the numerical scheme has been discussed in detail. In addition, the accuracy and stability of the used numerical method are also analysed in the spectral plane with the help of amplification factor and group velocity contour plot. The obtained numerical solutions have been validated with the existing literature. Four different aspect ratio cases (L/H = 3/4,4/3,4/5 and 5/4) have been considered for the simulations including the case of square cavity. It has been observed that L/H = 5/4 case provides best results among all others.

Non-polynomial B-spline and shifted Jacobi spectral collocation techniques to solve time-fractional nonlinear coupled Burgers’ equations numerically

This paper proposes two numerical approaches for solving the coupled nonlinear time-fractional Burgers’ equations with initial or boundary conditions on the interval $[0, L]$ [ 0 , L ] . The first method is the non-polynomial B-spline method based on L1-approximation and the finite difference approximations for spatial derivatives. The method has been shown to be unconditionally stable by using the Von-Neumann technique. The second method is the shifted Jacobi spectral collocation method based on an operational matrix of fractional derivatives. The proposed algorithms’ main feature is that when solving the original problem it is converted into a nonlinear system of algebraic equations. The efficiency of these methods is demonstrated by applying several examples in time-fractional coupled Burgers equations. The error norms and figures show the effectiveness and reasonable accuracy of the proposed methods.

A FOURTH ORDER ENERGY STABLE FINITE DIFFERENCE SCHEME FOR A TIMOSHENKO BEAM EQUATIONS WITH LOCALLY DISTRIBUTED FEEDBACK

This work deals with the numerical solution of a control problem governed by the Timoshenko beam equations with locally distributed feedback. We apply a fourth-order Compact Finite Difference (CFD) approximation for the discretizing spatial derivatives and a Forward second order method for the resulting linear system of ordinary differential equations. Using the energy method, we derive energy relation for the continuous model, and design numerical scheme that preserve a discrete analogue of the energy relation. Numerical results show that the CFD approximation of fourth order give an efficient method for solving the Timoshenko beam equations.

Acoustic wave propagation with new spatial implicit and temporal high-order staggered-grid finite-difference schemes

The implicit staggered-grid (SG) finite-difference (FD) method can obtain significant improvement in spatial accuracy for performing numerical simulations of wave equations. Normally, the second-order central grid FD formulas are used to approximate the temporal derivatives, and a relatively fine time step has to be used to reduce the temporal dispersion. To obtain high accuracy both in space and time, we propose a new spatial implicit and temporal high-order SG FD stencil in the time–space domain by incorporating some additional grid points to the conventional implicit FD one. Instead of attaining the implicit FD coefficients by approximating spatial derivatives only, we calculate the coefficients by approximating the temporal and spatial derivatives simultaneously through matching the dispersion formula of the seismic wave equation and compute the FD coefficients of our new stencil by two schemes. The first one is adopting a variable substitution-based Taylor-series expansion (TE) to derive the FD coefficients, which can attain (2M + 2)th-order spatial accuracy and (2N)th-order temporal accuracy. Note that the dispersion formula of our new stencil is non-linear with respect to the axial and off-axial FD coefficients, it is complicated to obtain the optimal spatial and temporal FD coefficients simultaneously. To tackle the issue, we further develop a linear optimisation strategy by minimising the L2-norm errors of the dispersion formula to further improve the accuracy. Dispersion analysis, stability analysis and modelling examples demonstrate the accuracy, stability and efficiency advantages of our two new schemes.

Optical Solitons and Vortices in Fractional Media: A Mini-Review of Recent Results

The article produces a brief review of some recent results which predict stable propagation of solitons and solitary vortices in models based on the nonlinear Schrödinger equation (NLSE) including fractional one-dimensional or two-dimensional diffraction and cubic or cubic-quintic nonlinear terms, as well as linear potentials. The fractional diffraction is represented by fractional-order spatial derivatives of the Riesz type, defined in terms of the direct and inverse Fourier transform. In this form, it can be realized by spatial-domain light propagation in optical setups with a specially devised combination of mirrors, lenses, and phase masks. The results presented in the article were chiefly obtained in a numerical form. Some analytical findings are included too, in particular, for fast moving solitons and the results produced by the variational approximation. Moreover, dissipative solitons are briefly considered, which are governed by the fractional complex Ginzburg–Landau equation.

Part 2: Spatially Dispersive Metasurfaces - IE-GSTC-SD Field Solver with Extended GSTCs

<div>An Integral Equation (IE) based field solver to compute the scattered fields from spatially dispersive metasurfaces is proposed and numerically confirmed using various examples involving physical unit cells. The work is a continuation of Part-</div><div>1 [1], which proposed the basic methodology of representing spatially dispersive metasurface structure in the spatial frequency domain, k. By representing the angular dependence of the surface susceptibilities in k as a ratio of two polynomials, the standard Generalized Sheet Transition Conditions (GSTCs) have been extended to include the spatial derivatives of both the difference and average fields around the metasurface. These extended boundary conditions are successfully integrated here into a standard IE-GSTC solver, which leads to the new IEGSTC-SD simulation framework presented here. The proposed IE-GSTC-SD platform is applied to various uniform metasurfaces, including a practical short conducting wire unit cell, as a representative practical example, for various cases of finite-sized flat and curvilinear surfaces. In all cases, computed field distributions are successfully validated, either against the semi-analytical Fourier decomposition method or the brute-force full-wave simulation of volumetric metasurfaces in the commercial Ansys FEM-HFSS simulator.</div>