Fourth-order Equation Research Articles

A mixed virtual element method for the two-dimensional Navier-Stokes equations in stream-function formulation

This work presents the formulation and analysis of a H1− conforming mixed virtual element method (VEM) for the two-dimensional stationary incompressible Navier-Stokes (NS) equations in stream-function formulation. By representing the velocity field as the curl of a stream function, we recast the second-order NS system into a fourth-order nonlinear equation for the scalar stream function, inherently satisfying the incompressibility constraint. Introducing a vorticity variable enables construction of H1− conforming VEM spaces for both stream function and vorticity and circumventing stringent C1− continuity constraints. The proposed method provides an initial exploration of stream function-vorticity discretizations on general polygonal meshes using the flexible VEM of arbitrary order. Existence and uniqueness of discrete solutions are established theoretically under a small data assumption. Optimal error estimates are then derived in the energy norm for the stream function, H1− norm for the stream function and L2− norm for the vorticity, rigorously demonstrating convergence. Numerical results validate the error analysis and illustrate the accuracy and robustness of the mixed VEM for simulation of incompressible flows on complex geometries.

Petrov-Galerkin method for small deflections in fourth-order beam equations in civil engineering

Abstract This study explores the Petrov–Galerkin method’s application in solving a linear fourth-order ordinary beam equation of the form u ″ ″ + q u = f u^{\prime\prime} ^{\prime\prime} +qu=f . The equation entails two distinct boundary conditions: pinned–pinned conditions on u u and u ′ u^{\prime} , and clamped–clamped conditions on u u and u ″ {u}^{^{\prime\prime} } . To satisfy these boundary conditions, we have built two sets of basis functions. The explicit forms of all spectral matrices were reported. The nonhomogeneous boundary conditions were easily handled using perfect transformations, ensuring the numerical solution’s accuracy. Detailed analysis of the method’s convergence was studied. Some numerical examples were presented, accompanied by comparisons with other existing methods in the literature.

Differential Equations of Fourth-Order with p-Laplacian-like Operator: Oscillation Theorems

Differential Equations of Fourth-Order with p-Laplacian-like Operator: Oscillation Theorems

On the radius of spatial analyticity for the quintic fourth-order nonlinear Schrodinger equation on ℝ2

On the radius of spatial analyticity for the quintic fourth-order nonlinear Schrodinger equation on ℝ²

Qualitative analysis of fourth-order hyperbolic equations

We investigate the qualitative properties of weak solutions to the boundary value problems for fourth-order linear hyperbolic equations with constant coefficients in a plane bounded domain convex with respect to characteristics. Our main scope is to prove some analog of the maximum principle, solvability, uniqueness and regularity results for weak solutions of initial and boundary value problems in the space L2. The main novelty of this paper is to establish some analog of the maximum principle for fourth-order hyperbolic equations. This question is very important due to natural physical interpretation and helps to establish the qualitative properties for solutions (uniqueness and existence results for weak solutions). The challenge to prove the maximum principle for weak solutions remains more complicated and at that time becomes more interesting in the case of fourth-order hyperbolic equations, especially, in the case of non-classical boundary value problems with data of weak regularity. Unlike second-order equations, qualitative analysis of solutions to fourth-order equations is not a trivial problem, since not only a solution is involved in boundary or initial conditions, but also its high- order derivatives. Other difficulty concerns the concept of weak solution of the boundary value problems with L2 – data. Such solutions do not have usual traces, thus, we have to use a special notion for traces to poss correctly the boundary value problems. This notion is traces associated with operator L or L-traces. We also derive an interesting interpretation (as periodicity of characteristic billiard or the John's mapping) of the Fredholm's property violation. Finally, we discuss some potential challenges in applying the results and proposed methods.

Multi-level physics informed deep learning for solving partial differential equations in computational structural mechanics

Physics-informed neural network has emerged as a promising approach for solving partial differential equations. However, it is still a challenge for the computation of structural mechanics problems since it involves solving higher-order partial differential equations as the governing equations are fourth-order nonlinear equations. Here we develop a multi-level physics-informed neural network framework where an aggregation model is developed by combining multiple neural networks, with each one involving only first-order or second-order partial differential equations representing different physics information such as geometrical, constitutive, and equilibrium relations of the structure. The proposed framework demonstrates a remarkable advancement over the classical neural networks in terms of the accuracy and computation time. The proposed method holds the potential to become a promising paradigm for structural mechanics computation and facilitate the intelligent computation of digital twin systems.

Generalized Linear Differential Equation using Hyers - Ulam Stability Approach

In this paper, We demonstrate the Hyers - Ulam stability of linear differential equation of fourth order. We interact with the differential equation\begin{align*}\gamma^{iv} (\omega) + \rho_1 \gamma{'''} (\omega)+ \rho_2 \gamma{''} (\omega) + \rho_3 \gamma' (\omega) + \rho_4 \gamma(\omega) = \chi(\omega),\end{align*}where $\gamma \in c^4 [\alpha,\beta], \chi \in [\alpha,\beta]$. Hyers-Ulam stability concerns the robustness of solutions of functional equations under small perturbations, ensuring that a solution approximately satisfying the equation is close to an exact solution. We extend this concept to fourth-order linear differential equations and continuous functions. Using fixed-point methods and various norms, we establish conditions under which such equations exhibit Hyers-Ulam stability. Several illustrative examples are provided to demonstrate the application of these results in specific cases, contributing to the growing understanding of stability in higher-order differential equations. Our findings have implications in both theoretical research and practical applications in physics and engineering.

A Novel and Accurate Algorithm for Solving Fractional Diffusion-Wave Equations

The main objective of this work is to apply a novel and accurate algorithm for solving the second-order and fourth-order fractional diffusion-wave equations (FDWEs). First, the desired equation is reduced to the corresponding Volterra integral equation (VIE). Then, the collocation method is applied, for which the Chebyshev cardinal functions (CCFs) have been considered as the bases. In this paper, the CCFs based on a Lobatto grid are introduced and used for the first time to solve these kinds of equations. To this end, the derivative and fractional integral operators are represented in CCFs. The main features of the method are simplicity, compliance with boundary conditions, and good accuracy. An exact analysis to show the convergence of the scheme is presented, and illustrative examples confirm our investigation.

Displacement-potential analysis of ply-stresses in a partially guided non-uniform short-column of hybrid laminated composites

Reliable and accurate analysis of local stresses in structural problems of fiber-reinforced laminated composites involving nonuniform geometry is of utmost practical importance, especially when they are subjected to local loads and supports. In this research, a new study of ply stresses in a partially guided, nonuniform, short composite column of hybrid laminates, subjected to an eccentrically distributed loading is conducted using an efficient computational scheme. The analytical modeling of the present mixed-boundary-value problem of laminated structure is realized in terms of a single scalar function of space variables (called as displacement-potential function), which is defined by the displacement components of plane elasticity, and is governed by a single fourth-order partial-differential equation of equilibrium. In this modeling approach, the original three-dimensional symmetric laminated column is reduced to a plane-stress scalar-field problem of an equivalent imaginary lamina in an attempt to obtain the corresponding strain field (global laminate strain), which is then used to determine the stress fields of individual plies through the respective transformed reduced stiffness matrices. A finite-difference based computational scheme is developed to solve the representative nonuniform equivalent lamina in terms of the displacement-potential function, which is however well capable of dealing with different ply materials as well as fiber orientations, together with the associated mixed and changeable physical conditions along the boundaries of the column. A 28-ply balanced hybrid laminate composed of boron-epoxy and glass-epoxy plies with various fiber orientation is considered as the structural material for the present column. The stress fields obtained for individual plies as well as the overall laminate are analyzed in the perspective of a number of practical issues of structural design, which include material hybridization, eccentricity of loading and attachment of partial guides to the column. The design stresses of the nonuniform laminated column are identified to be affected significantly by the material and orientation of fibers as well as the laminate hybridization. An attempt is also made to verify the soundness and accuracy of the present computational scheme through the comparison of the present potential-function solutions with those obtained by the conventional computational method.

Critical buckling load analysis of Euler-Bernoulli beam on two-parameter foundations using Galerkin method

The critical buckling load determination of Euler-Bernoulli beams on two-parameter elastic foundations (EBBo2PFs) is important to avert buckling failures. The governing equation for buckling of thin beam on two-parameter elastic foundation is a homogeneous ordinary differential equation (HODE) of fourth order and constant parameters when the beam is prismatic and homogeneous. The HODE has been solved in this work by Galerkin method for simply supported, clamped and clamped-simply supported ends. One-parameter algebraic shape function formulation was used to reduce the problem to an algebraic eigenvalue problem, which is solved to find the critical buckling load for each studied case. The critical buckling load for EBBo2PF for simply supported boundary conditions was found to be closely identical to the exact solutions. The solutions for clamped-clamped edges and clamped-simple supports were found to be accurate. The merit of the Galerkin method is the simplicity and the accuracy even when one-parameter shape function has been used.

Novel complexiton, rational wave, multi-lumps and the kink solitary wave solutions to the new (3+1)-dimensional integrable fourth-order equation for shallow water waves

In this exploration, we aim to seek a number of new exact solutions to the new (3+1)-dimensional integrable fourth-order nonlinear equation, which is widely used to describe the shallow water waves. Employing the Cole-Hopf transformation, we develop its bilinear form. Then, taking advantage of the ansatz function method, a new functional form is utilized to probe the singular complexiton solutions. Based on which, the non-singular complexiton solutions are derived by imposing the constraint conditions. In addition, we find the rational wave solutions and multi-lumps solutions wielding the rational function method and new homoclinic method respectively. At the end, we investigate the kink solitary wave solutions using the variational approach that is based on the variational principle and Ritz method. Meanwhile, the Hamiltonian of the system is also elaborated. Correspondingly, the graphic descriptions of the extracted results are presented to unfold their dynamic behaviors through Maple. As we all know, the findings of this paper are firstly reported and can enlarge the exact solutions of the considered PDE.

Existence of Three Positive Solutions for Boundary Value Problem of Fourth Order with Sign-Changing Green’s Function

In this paper, we examine a fourth-order equation that has parameter dependency and boundary conditions in three different places. We prove some of the features of the relevant asymmetric Green’s function and infer its exact form. The resulting solutions are still positive and decreasing functions on the entire interval of the Green’s function definition, and they are concave in a specific subinterval, despite the fact that the function’s sign changes on the square of its definition. The fixed point theorem of Krasnoselskii is the foundation of the existence arguments. Next, using the Leggett–Williams fixed point theorem, it is concluded that there are at least three positive solutions. Lastly, an example is provided, to highlight the primary findings of the manuscript.

Mild solutions to fourth-order parabolic equations modeling thin film growth with time fractional derivative

In this article, we study initial-boundary problems for fourth-order nonlinear parabolic equations modeling thin film growth with Caputo-type time fractional derivative. By means of the theory of abstract fractional calculus and \(L^p-L^q\) estimates, we establish the existence and uniqueness of local mild solutions in the spaces \(C([0,T]; L^{\frac{\beta N}{2-\beta}}(\Omega))\) with \(1<\beta<2\). Moreover, the local solutions can be extended globally if the initial data is sufficiently small. For more information see https://ejde.math.txstate.edu/Volumes/2024/58/abstr.html

Computation of dynamic deflection in thin elastic beam via symmetries

The deflection profiles governed by Euler Bernoulli's fourth-order equations under varied applied loads are investigated in this research. This study provides essential insights for engineers designing aircraft components, bridges, and similar structures, ensuring system safety and efficiency. The investigation emphasizes critical factors such as amplitude and frequency, load history, and material properties. Initially, conservation laws of the equations with applied loads are derived by expressing them in the Euler-Lagrange form, where the resultant conservation laws satisfy the divergence expression. The association between symmetries and conservation laws is demonstrated, followed by the application of double reduction theory, which reduces both the variables and the order of the equation. Graphical representations of the outcomes illustrate the impact of load variations on the beam's deflection profiles. These visual aids facilitate a deeper understanding of the influence of different loading conditions. A comparison between varying loads is presented, showcasing the impact of these variations on structural behavior. The findings are crucial for enhancing structural design and ensuring safety under varied loading conditions, showcasing the novelties in the analytical approach and the practical applications of the derived results.

Investigation of the Oscillatory Properties of Fourth-Order Delay Differential Equations Using a Comparison Approach with First- and Second-Order Equations

This paper investigates the oscillatory behavior of solutions to fourth-order functional differential equations (FDEs) with multiple delays and a middle term. By employing a different comparison method approach with lower-order equations, the study introduces enhanced oscillation criteria. A key strength of the proposed method is its ability to reduce the complexity of the fourth-order equation by converting it into first- and second-order forms, allowing for the application of well-established oscillation theories. This approach not only extends existing criteria to higher-order FDEs but also offers more efficient and broadly applicable results. Detailed comparisons with previous research confirm the method’s effectiveness and broader relevance while demonstrating the feasibility and significance of our results as an expansion and improvement of previous results.

Chaotic systems based on higher-order oscillatory equations

This paper discusses the design process toward new lumped chaotic systems that originates in higher-order ordinary differential equations commonly used as description of ideal oscillators. In investigated third-order case, two chaotic oscillators were constructed. These systems are dual in the sense of vector field geometry local to fixed points. The existence of robust chaos was proved by both standard routines of numerical analysis and practical measurement. For the fourth-order oscillatory equation, the concept based on interaction between superinductor and supercapacitor was examined in detail. Since both “superelements” are active, the nonlinearity essential to the evolution of chaos is fully passive. It is demonstrated that complex motion is robust and does not represent long transient behavior or numerical artefact. The existence of chaos was verified using standard quantifiers of the flow, such as the largest Lyapunov exponents, recurrence plots, approximate entropy and sensitivity calculation. A good final agreement between theoretical assumptions and practical results will be concluded, on a visual comparison basis.

Global Solution for a Fourth-Order Nonlinear Elliptic Equation

Objective: In this work, the existence of a global solution for a fourth-order nonlinear elliptic equation with boundary conditions is studied. Method: Variational methods are applied to the proposed model, naturally using the Banach Fixed Point Contraction Theorem. Results and Discussion: The conclusion of the result obtained is that for the nonlinear elliptic equation of the fourth order they have equivalent norms in the conveniently defined spaces. On the other hand, the existence of a solution exists whenever Banach's fixed-point theorem can be applied. Implications of the research: Apply in the first university years the application of global solution for a fourth-order nonlinear elliptic equation in all university careers and specialties as part of the management of the tools of Information and Communication Technology (TIC). Originality/Value: This paper also provides an alternative approach to Analysis of global solution for a fourth-order nonlinear elliptic equation.

Nodal solutions for a nonlocal fourth order equation of Kirchhoff type

We study the bifurcation behavior of nodal solutions for the Kirchhoff type beam equation u′′′′−M(∫01|u′|2dx)u′′=λf(x,u),x∈(0,1),u(0)=u(1)=u′′(0)=u′′(1)=0, where λ∈R is a parameter, M:[0,∞)→[0,∞) and f:[0,1]×R→R are smooth functions. We obtain the existence of nodal solutions under some suitable conditions. The proof of our main result is based upon bifurcation techniques.

Asymptotic behavior of solutions to a fourth-order degenerate parabolic equation

Asymptotic behavior of solutions to a fourth-order degenerate parabolic equation

Simplified conservative discretization of the Cahn-Hilliard-Navier-Stokes equations

In this paper we construct a novel discretization of the Cahn-Hilliard equation coupled with the Navier-Stokes equations. The Cahn-Hilliard equation models the separation of a binary mixture. We construct a very simple time integration scheme for simulating the Cahn-Hilliard equation, which is based on splitting the fourth-order equation into two second-order Helmholtz equations. We combine the Cahn-Hilliard equation with the Navier-Stokes equations to simulate phase separation in a two-phase fluid flow in two dimensions. The scheme conserves mass and momentum and exhibits consistency between mass and momentum, allowing it to be used with large density ratios. We introduce a novel discretization of the surface tension force from the phase-field variable that has finite support around the transition region. The model has a parameter that allows it to transition from a smoothed continuum surface force to a fully sharp interface formulation. We show that our method achieves second-order accuracy, and we compare our method to previous work in a variety of experiments.