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
Orthonormal eigenfunction expansions for sixth-order boundary value problems
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
- Award ID(s):
- 2029540
- PAR ID:
- 10481629
- Editor(s):
- Todorov, M D
- Publisher / Repository:
- IOP Publishing
- Date Published:
- Journal Name:
- Journal of Physics: Conference Series
- Volume:
- 2675
- ISSN:
- 1742-6588
- Page Range / eLocation ID:
- 012016
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Altenbach, Holm; Eremeyev, Victor A (Ed.)We propose an analytical approach to solving nonlocal generalizations of the Euler–Bernoulli beam. Specifically, we consider a version of the governing equation recently derived under the theory of peridynamics. We focus on the clamped–clamped case, employing the natural eigenfunctions of the fourth derivative subject to these boundary conditions. Static solutions under different loading conditions are obtained as series in these eigenfunctions. To demonstrate the utility of our proposed approach, we contrast the series solution in terms of fourth-order eigenfunctions to the previously obtained Fourier sine series solution. Our findings reveal that the series in fourth-order eigenfunctions achieve a given error tolerance (with respect to a reference solution) with ten times fewer terms than the sine series. The high level of accuracy of the fourth-order eigenfunction expansion is due to the fact that its expansion coefficients decay rapidly with the number of terms in the series, one order faster than the Fourier series in our examples.more » « less
-
In this paper, a non-trivial system governed by a continuum PT-symmetric Hamiltonian is discussed. We show that this Hamiltonian is iso-spectral to the simple harmonic oscillator. We find its eigenfunctions and the path in the complex plane along which these functions form an orthonormal set. We also find the hidden symmetry operator, C, for this system. All calculations are performed analytically and without approximation.more » « less
-
Adhesive contact between a thin elastic sheet and a substrate arises in a range of biological, physical and technological applications. By considering the dynamics of this process that naturally couples fluid flow, long-wavelength elastic deformations and microscopic adhesion, we analyse a sixth-order thin-film equation for the short-time dynamics of the onset of adhesion and the long-time dynamics of a steadily propagating adhesion front. Numerical solutions corroborate scaling laws and asymptotic analyses for the characteristic waiting time for adhesive contact and for the speed of the adhesion front. A similarity analysis of the governing partial differential equation further allows us to determine the shape of a fluid-filled blister ahead of the adhesion front. Finally, our analysis reveals a near-singular behaviour at the moving elastohydrodynamic contact line with an effective boundary condition that might be useful in other related problems.more » « less
An official website of the United States government

