skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


This content will become publicly available on November 1, 2026

Title: Non-self-adjoint sixth-order eigenvalue problems arising from clamped elastic thin films on closed domains
Sixth-order boundary value problems (BVPs) arise in thin-film flows with a surface that has elastic bending resistance. We consider the case in which the elastic interface is clamped at the lateral walls of a closed trough and thus encloses a finite amount of fluid. For a slender film undergoing infinitesimal deformations, the displacement of the elastic surface from its initial equilibrium position obeys a sixth-order (in space) initial boundary value problem (IBVP). To solve this IBVP, we construct a set of odd and even eigenfunctions that intrinsically satisfy the boundary conditions (BCs) of the original IBVP. These eigenfunctions are the solutions of a non-self-adjoint sixth-order eigenvalue problem (EVP). To use the eigenfunctions for series expansions, we also construct and solve the adjoint EVP, leading to another set of even and odd eigenfunctions, which are orthogonal to the original set (biorthogonal). The eigenvalues of the adjoint EVP are the same as those of the original EVP, and we find accurate asymptotic formulas for them. Next, employing the biorthogonal sets of eigenfunctions, a Petrov–Galerkin spectral method for sixth-order problems is proposed, which can also handle lower-order terms in the IBVP. The proposed method is tested on two model sixth-order BVPs, which admit exact solutions. We explicitly derive all the necessary formulas for expanding the quantities that appear in the model problems into the set(s) of eigenfunctions. For both model problems, we find that the approximate Petrov–Galerkin spectral solution is in excellent agreement with the exact solution. The convergence rate of the spectral series is rapid, exceeding the expected sixth-order algebraic rate.  more » « less
Award ID(s):
2029540
PAR ID:
10650481
Author(s) / Creator(s):
;
Publisher / Repository:
IOP Publishing
Date Published:
Journal Name:
Journal of Physics: Conference Series
Volume:
3145
ISSN:
1742-6588
Page Range / eLocation ID:
012005
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; (, Journal of Physics: Conference Series)
    Todorov, M D (Ed.)
    Sixth-order boundary value problems (BVPs) arise in thin-film flows with a surface that has elastic bending resistance. To solve such problems, we first derive a complete set of odd and even orthonormal eigenfunctions — resembling trigonometric sines and cosines, as well as the so-called “beam” functions. These functions intrinsically satisfy boundary conditions (BCs) of relevance to thin-film flows, since they are the solutions of a self-adjoint sixth-order Sturm–Liouville BVP with the same BCs. Next, we propose a Galerkin spectral approach for sixth-order problems; namely the sought function as well as all its derivatives and terms appearing in the differential equation are expanded into an infinite series with respect to the derived complete orthonormal (CON) set of eigenfunctions. The unknown coefficients in the series expansion are determined by solving the algebraic system derived by taking successive inner products with each member of the CON set of eigenfunctions. The proposed method and its convergence are demonstrated by solving two model sixth-order BVPs. 
    more » « less
  2. ; (, Journal of Engineering Mathematics)
    A linear sixth-order partial differential equation (PDE) of “parabolic” type describes the dynamics of thin liquid films beneath surfaces with elastic bending resistance when deflections from the equilibrium film height are small. On a finite domain, the associated sixth-order eigenvalue problem is self-adjoint for the boundary conditions corresponding to a thin film in a closed trough, and the eigenfunctions form a complete orthonormal set. Using these eigenfunctions, we derive the Green’s function for the governing sixth-order PDE on a finite interval and compare it to the known infinite-line solution. Further, we propose a Galerkin spectral method based on the constructed sixth-order eigenfunctions and their derivative expansions. The system of ordinary differential equations for the time-dependent expansion coefficients is solved by standard numerical methods. The numerical approach is applied to versions of the governing PDE with a second-order spatial derivative (in addition to the sixth-order one), which arises from gravity acting on the film. In the absence of gravity, we demonstrate the self-similar intermediate asymptotics of initially localized disturbances on the film surface, at least until the disturbances “feel” the finite boundaries, and show that the derived Green’s function is an attractor for such solutions. In the presence of gravity, we use the proposed Galerkin numerical method to demonstrate that self-similar behavior persists, albeit for shortened intervals of time, even for large values of the gravity-to-bending ratio. 
    more » « less
  3. ; (, Journal of Physics: Conference Series)
    The flow of a thin viscous liquid layer under an elastic film arises in natural processes, such as magmatic intrusions between rock strata, and industrial applications, such as coating surfaces with cured polymeric films. We study the linear dynamics of small perturbations to the equilibrium state of the film in a closed trough (finite film). Specifically, we are interested in the spreading (early-time) and leveling (late-time) dynamics as the film adjusts to equilibrium, starting from different initial perturbations. We consider both smooth and non-smooth spatially symmetric and localized initial conditions (perturbations). We find the exact series solutions for the film height, using the sixth-order complete orthonormal eigenfunctions associated with the posed initial-boundary-value problem derived in our previous work [Papanicolaou N C and Christov I C 2023J. Phys.: Conf. Ser.2675012016]. We show that the evolution of the perturbations begins with the spreading of the localized perturbations, followed by their mutual interaction as they spread, and finally, interactions with the confining lateral boundaries of the domain as the perturbations level. In particular, we highlight how the leading eigenvalues of the problem determine the scalings of certain figures of merit with time. 
    more » « less
  4. ; (, Computational methods in applied mathematics)
    Numerical integration of the stiffness matrix in higher-order finite element (FE) methods is recognized as one of the heaviest computational tasks in an FE solver. The problem becomes even more relevant when computing the Gram matrix in the algorithm of the Discontinuous Petrov Galerkin (DPG) FE methodology. Making use of 3D tensor-product shape functions, and the concept of sum factorization, known from standard high-order FE and spectral methods, here we take advantage of this idea for the entire exact sequence of FE spaces defined on the hexahedron. The key piece to the presented algorithms is the exact sequence for the one-dimensional element, and use of hierarchical shape functions. Consistent with existing results, the presented algorithms for the integration of H1, H(curl), H(div), and L2 inner products, have the O(p7) computational complexity in contrast to the O(p9) cost of conventional integration routines. Use of Legendre polynomials for shape functions is critical in this implementation. Three boundary value problems under different variational formulations, requiring combinations of H1, H(div) and H(curl) test shape functions, were chosen to experimentally assess the computation time for constructing DPG element matrices, showing good correspondence with the expected rates. 
    more » « less
  5. ; (, Acta Numerica)
    The discontinuous Petrov–Galerkin (DPG) method is a Petrov–Galerkin finite element method with test functions designed for obtaining stability. These test functions are computable locally, element by element, and are motivated by optimal test functions which attain the supremum in an inf-sup condition. A profound consequence of the use of nearly optimal test functions is that the DPG method can inherit the stability of the (undiscretized) variational formulation, be it coercive or not. This paper combines a presentation of the fundamentals of the DPG ideas with a review of the ongoing research on theory and applications of the DPG methodology. The scope of the presented theory is restricted to linear problems on Hilbert spaces, but pointers to extensions are provided. Multiple viewpoints to the basic theory are provided. They show that the DPG method is equivalent to a method which minimizes a residual in a dual norm, as well as to a mixed method where one solution component is an approximate error representation function. Being a residual minimization method, the DPG method yields Hermitian positive definite stiffness matrix systems even for non-self-adjoint boundary value problems. Having a built-in error representation, the method has the out-of-the-box feature that it can immediately be used in automatic adaptive algorithms. Contrary to standard Galerkin methods, which are uninformed about test and trial norms, the DPG method must be equipped with a concrete test norm which enters the computations. Of particular interest are variational formulations in which one can tailor the norm to obtain robust stability. Key techniques to rigorously prove convergence of DPG schemes, including construction of Fortin operators, which in the DPG case can be done element by element, are discussed in detail. Pointers to open frontiers are presented. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email