Fourth-order Problem Research Articles

Synchronous consistent integration for superconvergent isogeometric analysis of structural vibrations

A commonly used procedure to improve the frequency accuracy of isogeometric structural vibration analysis is the design of special superconvergent quadrature rules in accordance with the frequency error measure derived from the generic stencil equations. However, it is also numerically observed that such superconvergent quadrature rules may not consistently lead to the expected accuracy. The reason underlying this inconsistency is related to the failure of quadrature rules to meet the condition of variational integration consistency. Meanwhile, by means of the tensor product nature of isogeometric basis functions and quadrature rules, it is theoretically proved that in one- as well as multi-dimensional scenarios, the variational integration consistency of the isogeometric weak formulation necessitates the quadrature rules be exact through 2(p-m)th degree polynomials, where p is the complete degree of multi-dimensional basis functions or the degree of their correspondent one-dimensional basis functions, m is the order of differentiation in Galerkin weak form. A synchronous consistent superconvergent quadrature rule ensuring frequency superconvergence should simultaneously fulfill both the constraint of frequency error measure and the condition of variational integration consistency. Subsequently, the issues associated with certain quadrature rules that violate the variational integration consistency are examined for both second- and fourth-order problems, and the corresponding synchronous consistent quadrature rules are particularly elaborated.

A new spline method on graded mesh for fourth-order time-dependent PDEs: application to Kuramoto-Sivashinsky and extended Fisher-Kolmogorov equations

In this research, we introduce a two-tier non-polynomial spline approach with graded mesh discretization for addressing fourth-order time-dependent partial differential equations, which find applications in various physical scenarios like the nonlinear Kuramoto-Sivashinsky equation and extended Fisher-Kolmogorov equation. Our method involves considering three spatial points at each time step for scheme development, achieving spatial accuracy of three and temporal accuracy of two. Notably, our approach offers the advantage of directly tackling singular biharmonic problems without the need for discretizing boundary conditions or incorporating fictitious points. It demonstrates unconditional stability when tested on fourth-order problems and outperforms previous methods, as evidenced by comparative analyses. Additionally, we present 2D and 3D plots illustrating numerical solutions for different benchmark scenarios, including the depiction of the chaotic nature of the Kuramoto-Sivashinsky equation under Gaussian initial conditions at various time levels.

Existence results for a fourth order problem with functional perturbed clamped beam boundary conditions

Abstract In this paper, we study the existence of positive solutions for a fourth order boundary value problem coupled to functional perturbed clamped beam boundary conditions. Our main ingredient is the classical fixed point index. The problem investigated is an extension of other problems studied in previous papers by covering very general nonlocal perturbed conditions on the boundary.

Local discontinuous Galerkin method for a singularly perturbed fourth-order problem of convection-diffusion type

We develop a local discontinuous Galerkin (LDG) method for a fourth-order singularly perturbed problem of convection-diffusion type. The existence and uniqueness of the computed solution are verified. Using the Shishkin mesh we derive an optimal-order energy-norm error estimate which is uniformly valid in the perturbation parameter. Numerical experiments are also given to support our theoretical findings.

Analysis of the boundary conditions for the ultraweak-local discontinuous Galerkin method of time-dependent linear fourth-order problems

In this paper, we study the ultraweak-local discontinuous Galerkin (UWLDG) method for time-dependent linear fourth-order problems with four types of boundary conditions. In one dimension and two dimensions, stability and optimal error estimates of order k + 1 k+1 are derived for the UWLDG scheme with polynomials of degree at most k k ( k ≥ 1 k\ge 1 ) for solving initial-boundary value problems. The main difficulties are the design of suitable penalty terms at the boundary for numerical fluxes and the construction of projections. More precisely, in two dimensions with the Dirichlet boundary condition, an elaborate projection of the exact boundary condition is proposed as the boundary flux, which, in combination with some proper penalty terms, leads to the stability and optimal error estimates. For other three types of boundary conditions, optimal error estimates can also be proved for fluxes without any penalty terms when special projections are designed to match different boundary conditions. Numerical experiments are presented to confirm the sharpness of theoretical results.

Discrete Helmholtz Decompositions of Piecewise Constant and Piecewise Affine Vector and Tensor Fields

AbstractDiscrete Helmholtz decompositions dissect piecewise polynomial vector fields on simplicial meshes into piecewise gradients and rotations of finite element functions. This paper concisely reviews established results from the literature which all restrict to the lowest-order case of piecewise constants. Its main contribution consists of the generalization of these decompositions to 3D and of novel decompositions for piecewise affine vector fields in terms of Fortin–Soulie functions. While the classical lowest-order decompositions include one conforming and one nonconforming part, the decompositions of piecewise affine vector fields require a nonconforming enrichment in both parts. The presentation covers two and three spatial dimensions as well as generalizations to deviatoric tensor fields in the context of the Stokes equations and symmetric tensor fields for the linear elasticity and fourth-order problems. While the proofs focus on contractible domains, generalizations to multiply connected domains and domains with non-connected boundary are discussed as well.

Thermophoretic particle deposition on double-diffusive Ree-Eyring fluid flow across two deformable rotating disks with Hall current and Ion slip

This work emphasizes the characteristics of thermophoretic particle deposition on Ree-Eyring fluid flow with Hall current and Ion slip effects between two deformable spinning disks. The novelty of the envisioned model is strengthened by considering generalized Fourier and Fick laws with thermal source/sink. The proposed model is optimized by taking multiple slip boundary constraints at the surface. The equations that govern the flow are simplified with the use of the proper transformations. The numerical solution is obtained using the boundary value problem of fourth order (bvp4c) method. Graphical representations of the emerging parameters are presented. Numerical values of skin friction coefficient, local Nusselt, and Stanton number are highlighted in tabular form. The velocity field depicts a deteriorating behavior on incrementing the Reynold number. For augmented values of Hall and ion slip parameters fluid velocity elevates. An upsurge is noticed in the axial and radial velocities on amplifying the rotation parameter. On augmenting the thermophoretic parameter the particle deposition at the interface of the upper disk diminishes. Growing values of rotation parameter, Weissenberg, and Reynold numbers imply an additional resistance to the fluid flow that ultimately reinforced the surface drag force causing a decline in the fluid velocity. The rate of heat transmission accelerates on enhancing the Reynold number, however, heat flux exhibits a decay on augmenting the thermal slip parameter. On augmenting the thermophoretic parameter particle deposition deteriorates at the interface of the upper disk. An excellent concordance is noticed when the outcomes are contrasted with the previous literature. For the value of the rotation parameter β=-1 and β=-8 the percentage error of the present analysis with Imtiaz et al. [68] are 0.00001 and 0.00002.

A novel spectral method and error estimates for fourth-order problems with mixed boundary conditions in a cylindrical domain

In this paper, a novel spectral method is presented and studied for the fourth order problem with mixed boundary in a cylindrical domain. The basic idea of our approach is to reduce the original problem into a series of decoupled two-dimensional fourth-order problems first, by using the cylindrical coordinate transformation and Fourier expansion, and then adopt the standard spectral method to solve the decoupled problems. A new essential pole condition is proposed to overcome the difficulty caused by the introduction of singularity and variable coefficients in cylindrical coordinate transformation. Existence and uniqueness of the weak solution and the discrete numerical solution are proved, and error estimates of the spectral method are derived. Furthermore, the efficient implementation of our algorithm is discussed, where a set of effective basis functions are constructed to ensure the sparsity of the mass matrix and stiffness matrix. Numerical examples are presented to validate the theoretical findings and the efficiency of our algorithm.

Local discontinuous Galerkin method for a singularly perturbed fourth-order problem of reaction–diffusion type

A fourth-order singularly perturbed problem of reaction–diffusion type is solved numerically by a local discontinuous Galerkin (LDG) method. Under suitable hypotheses, we prove optimal convergence of the LDG method on a Shishkin mesh; that is, when piecewise polynomials of degree k are used on a Shishkin mesh with N mesh intervals, we obtain O((N−1lnN)k+1/2) convergence in the energy norm. The error bound is uniformly valid with respect to the singular perturbation parameter. In the error analysis, we exploit a relationship between the numerical solution of the third-order derivative with the gradient, the numerical solution and its element interface jump. We discuss also the convergence of the LDG method on two Bakhvalov-type meshes. Numerical experiments indicate that our error estimate is sharp.

Numerical shape optimization of the Canham-Helfrich-Evans bending energy

In this paper we propose a novel numerical scheme for the Canham-Helfrich-Evans bending energy based on a three-field lifting procedure of the distributional shape operator to an auxiliary mean curvature field. Together with its energetic conjugate scalar stress field as Lagrange multiplier the resulting fourth order problem is circumvented and reduced to a mixed saddle point problem involving only second order differential operators. Further, we derive its analytical first variation (also called first shape derivative), which is valid for arbitrary polynomial order, and discuss how the arising shape derivatives can be computed automatically in the finite element software NGSolve. We finish the paper with several numerical simulations showing the pertinence of the proposed scheme and method.

An Efficient Legendre–Galerkin Approximation for Fourth-Order Elliptic Problems with SSP Boundary Conditions and Variable Coefficients

Under simply supported plate (SSP) boundary conditions, a numerical method based on the higher-order Legendre polynomial approximation was studied and developed for fourth-order problems with variable coefficients. We first divide the SSP boundary conditions into two types, namely, forced boundary conditions and natural boundary conditions. According to the forced boundary conditions, an appropriate Sobolev space is defined, and a variational formulation and a discrete scheme associated with the original problem are established. Then, the existence and uniqueness of this weak solution and approximate solution are proved. By using the Céa lemma and the tensor Jacobian polynomial approximation, we further obtain the error estimation for the numerical solutions. In addition, we use the orthogonality of Legendre polynomials to construct a set of effective basis functions and derive the equivalent tensor product linear system associated with the discrete scheme, respectively. Finally, some numerical tests were carried out to validate our algorithm and theoretical analysis.

GMLS-based numerical manifold method in mechanical analysis of thin plates with complicated shape or cutouts

We proposed a generalized moving least squares (GMLS) based numerical manifold method (NMM) and applied to investigate mechanical behaviors of thin plates with various shape and cutouts. The influence domains of nodes in generalized moving least squares (GMLS) interpolation are employed to construct the mathematical cover, which is specifically fitted for the fourth-order problems with rotational degrees of freedom. Moreover, the GMLS-based NMM can naturally treat plates with arbitrary complex geometry and openings without resorting to those complicated criteria and tailored operations while simplifying the physical cover generation process. In numerical study, a suitable mass lumping scheme is formulated and advised for transverse vibration analysis. Bending, vibration and buckling behaviors of thin plates with different shapes and cutouts are studied. Numerical results manifest that GMLS-based NMM can accurately model transverse mechanical behaviors.

A fully discrete local discontinuous Galerkin method for variable-order fourth-order equation with Caputo-Fabrizio derivative based on generalized numerical fluxes

<abstract><p>In this paper, an effective numerical method for the variable-order(VO) fourth-order problem with Caputo-Fabrizio derivative will be constructed and analyzed. Based on generalized alternating numerical flux, appropriate spatial and temporal discretization, we get a fully discrete local discontinuous Galerkin(LDG) scheme. The theoretic properties of the fully discrete LDG scheme are proved in detail by mathematical induction, and the method is proved to be unconditionally stable and convergent with $ {\rm O}(\tau+{h^{k+1}}) $, where $ h $ is the spatial step, $ \tau $ is the temporal step and $ k $ is the degree of the piecewise $ P^k $ polynomial. In order to show the efficiency of our method, some numerical examples are carried out by Matlab.</p></abstract>

An efficient Fourier spectral method and error analysis for the fourth order problem with periodic boundary conditions and variable coefficients

<abstract><p>We propose in this paper an efficient algorithm based on the Fourier spectral-Galerkin approximation for the fourth-order elliptic equation with periodic boundary conditions and variable coefficients. First, by using the Lax-Milgram theorem, we prove the existence and uniqueness of weak solution and its approximate solution. Then we define a high-dimensional $ L^2 $ projection operator and prove its approximation properties. Combined with Céa lemma, we further prove the error estimate of the approximate solution. In addition, from the Fourier basis function expansion and the properties of the tensor, we establish the equivalent matrix form based on tensor product for the discrete scheme. Finally, some numerical experiments are carried out to demonstrate the efficiency of the algorithm and correctness of the theoretical analysis.</p></abstract>

A new mixed finite-element method for H2 elliptic problems

Fourth-order differential equations play an important role in many applications in science and engineering. In this paper, we present a three-field mixed finite-element formulation for fourth-order problems, with a focus on the effective treatment of the different boundary conditions that arise naturally in a variational formulation. Our formulation is based on introducing the gradient of the solution as an explicit variable, constrained using a Lagrange multiplier. The essential boundary conditions are enforced weakly, using Nitsche's method where required. As a result, the problem is rewritten as a saddle-point system, requiring analysis of the resulting finite-element discretization and the construction of optimal linear solvers. Here, we discuss the analysis of the well-posedness and accuracy of the finite-element formulation. Moreover, we develop monolithic multigrid solvers for the resulting linear systems. Two and three-dimensional numerical results are presented to demonstrate the accuracy of the discretization and efficiency of the multigrid solvers proposed.

Compact difference scheme for time-fractional nonlinear fourth-order diffusion equation with time delay

We construct a compact difference scheme for two-dimensional time-fractional nonlinear fourth-order diffusion equation with time delay. By choosing the second-order spatial derivative of the original unknown as an additional variable, the fourth-order problem is transformed into a second-order system. Then the linearized compact difference scheme is formulated by employing the L2−1σ formula and the compact difference operator to approximate temporal Caputo derivative and spatial second-order derivatives, respectively. The unique solvability of the proposed scheme is proved. Then the scheme is proved to be stable and convergent with the rate of second order in time and fourth order in space. The efficiency of the scheme is supported by some numerical experiments.

Almost-[formula omitted] splines: Biquadratic splines on unstructured quadrilateral meshes and their application to fourth order problems

Isogeometric Analysis generalizes classical finite element analysis and intends to integrate it with the field of Computer-Aided Design. A central problem in achieving this objective is the reconstruction of analysis-suitable models from Computer-Aided Design models, which is in general a non-trivial and time-consuming task. In this article, we present a novel spline construction, that enables model reconstruction as well as simulation of high-order PDEs on the reconstructed models. The proposed almost-C1splines are biquadratic splines on fully unstructured quadrilateral meshes (without restrictions on placements or number of extraordinary vertices). They are C1 smooth at all regular and extraordinary vertices. Moreover, they are C1 smooth across all edges between regular vertices and C0 smooth across all edges that are adjacent to an extraordinary vertex. The splines thus form H2-nonconforming analysis-suitable discretization spaces. This is the lowest-degree unstructured spline construction that can be used to solve fourth-order problems. The associated spline basis is non-singular and has several B-spline-like properties (e.g., partition of unity, non-negativity, local support), the almost-C1 splines are described in an explicit Bézier-extraction-based framework that can be easily implemented. Numerical tests suggest that the basis is well-conditioned and exhibits optimal approximation behaviour.

An approximate [formula omitted] multi-patch space for isogeometric analysis with a comparison to Nitsche’s method

We present an approximately C1-smooth multi-patch spline construction which can be used in isogeometric analysis (IGA). The construction extends the one presented in [1] for two-patch domains. A key property of IGA is that it is simple to achieve high order smoothness within a single patch. However, to represent more complex geometries one often uses a multi-patch construction. In this case, the global continuity for the basis functions is in general only C0. Therefore, to obtain C1-smooth isogeometric functions, a special construction for the basis is needed. Such spaces are of interest when solving numerically fourth-order problems, such as the biharmonic equation or Kirchhoff–Love plate/shell formulations, using an isogeometric Galerkin method.Isogeometric spaces that are globally C1 over multi-patch domains can be constructed as in [2–6]. The constructions require geometry parametrizations that satisfy certain constraints along the interfaces, so-called analysis-suitable G1 parametrizations. To allow C1 spaces over more general multi-patch parametrizations, one needs to increase the polynomial degree and/or to relax the C1 conditions. Thus, we define function spaces that are not exactly C1 but only approximately. We adopt the construction for two-patch domains, as developed in [1], and extend it to more general multi-patch domains.We employ the construction for a biharmonic model problem and compare the results with Nitsche’s method. We compare both methods over complex multi-patch domains with non-trivial interfaces. The numerical tests indicate that the proposed construction converges optimally under h-refinement, comparable to the solution using Nitsche’s method. In contrast to weakly imposing coupling conditions, the approximate C1 construction is explicit and no additional terms need to be introduced to stabilize the method/penalize the jump of the derivative at the interface. Thus, the new proposed method can be used more easily as no parameters need to be estimated.

A fully discrete plates complex on polygonal meshes with application to the Kirchhoff–Love problem

In this work we develop a novel fully discrete version of the plates complex, an exact Hilbert complex relevant for the mixed formulation of fourth-order problems. The derivation of the discrete complex follows the discrete de Rham paradigm, leading to an arbitrary-order construction that applies to meshes composed of general polygonal elements. The discrete plates complex is then used to derive a novel numerical scheme for Kirchhoff–Love plates, for which a full stability and convergence analysis are performed. Extensive numerical tests complete the exposition.

Multiple solutions for a fourth order problem involving Leray-Lions type operator

Multiple solutions for a fourth order problem involving Leray-Lions type operator