Galerkin Approach Research Articles

Legendre spectral volume methods for Allen–Cahn equations by the direct discontinuous Galerkin formula

In this paper, we introduce novel class of Legendre spectral volume (LSV) methods for solving Allen–Cahn equations. Each spectral volume (SV) is refined with k Gauss–Legendre points to define an arbitrary order control volume (CV). Moreover, the second derivative is handled using the direct discontinuous Galerkin (DDG) approach. Furthermore, four numerical experiments are detailed including 1D and 2D Allen–Cahn equations with Neumann and periodic boundary conditions. These experiments demonstrate the stability and accuracy in capturing phase transitions of the approach. Meanwhile, we also show the LSV methods can maintain physical properties such as energy dissipation and uniform boundedness. It is worth mentioning that we observe that the LSV methods achieve both optimal convergence and superconvergence as the numerical flux parameter is carefully selected.

Nonlinear aeroelastic behavior of a two-dimensional heated panel by irregular shock reflection considering viscoelastic damping

This paper investigates the aeroelastic stability and nonlinear aeroelastic behavior of a two-dimensional heated panel in irregular shock reflection and extends prior work to include the effects of viscoelasticity. The aeroelastic model is formulated using the von Kármán large deflection plate theory and the Kelvin–Voigt damping model, accompanied by the quasi-steady thermal stress theory. The unsteady aerodynamic pressure is evaluated through the piston theory and the compressibility-corrected potential theory. The Galerkin approach is used to discretize the governing equation. The Lyapunov indirect method is applied to conduct theoretical analysis, obtaining the aeroelastic stability boundary. Also, the nonlinear aeroelastic response is numerically simulated via the fourth-order Runge–Kutta method. The proper orthogonal decomposition is applied to the panel deflection to manifest the influence of various system parameters. It is demonstrated that the shock wave aggravates the aerodynamic heating, lowering the critical buckling temperature. The viscoelastic damping restricts the impact of shock impingement location and shock strength on the stability boundary and also transforms the chaotic motions into periodic LCOs.

Radical Petrov–Galerkin Approach for the Time-Fractional KdV–Burgers’ Equation

Radical Petrov–Galerkin Approach for the Time-Fractional KdV–Burgers’ Equation

Effects of rotation and chemical reactions on an electroconvection in Maxwell dielectric nanofluids within anisotropic porous media

ABSTRACT This study explores the complex interactions between rotation, chemical reactions, and electroconvection in Maxwell dielectric nanofluids within anisotropic porous media. The research aims to elucidate how these factors influence the heat transfer efficiency and fluid dynamics under the influence of external alternating current (AC) electric fields. The normal mode method is employed to analyse the stability of the system, allowing for the examination of oscillatory modes and their growth rates. By introducing perturbations in the physical parameters, the study derives a set of ordinary linear differential equations that characterize the system's response to external influences. The coupled differential equations were solved using the Galerkin approach for first approximation while meeting the necessary boundary conditions. Analytical solution of the representative equation for stress-free boundary states yields Rayleigh number formulas for nonoscillatory and oscillatory mode occur. Oscillatory modes are seen for both bottom and top-heavy nanoparticle distributions. The electric Rayleigh number, thermal Prandtl number, and stress relaxation parameter increases, whereas the Brinkman-Darcy number delays the onset of stationary and oscillatory convection. The study finds that positive concentration Rayleigh numbers stabilize the system, while negative ones lead to instability. Additionally, it demonstrates that external alternating current electric fields influence convection behavior and provides formulas for critical Rayleigh numbers, enhancing our understanding of fluid interactions in various engineering and technological applications.

Galerkin approximation of Burgers-Huxley equation with fractional order arising in reaction mechanism and diffusion transport

Abstract Burgers-Huxley model depicts a prototype model of the interaction of convection effect, reaction mechanisms and diffusion transport, used to study the liquid crystal and nerve fibers. This study introduces Galerkin approximation for time-fractional Burgers-Huxley equation (TFBHE). The Caputo derivative is used to evaluate the temporal part using the L 1 formula. The Galerkin approach employs cubic B-spline as a shape and test function, resulting in a symmetric matrix that is easily convergent. In addition, the three-point quadrature rule is implemented to evaluate the integration of complex function . The Von Neumann analysis is used to discuss stability of the scheme. The performance and robustness of the technique is measured using various error norms The results are compared with the exact solution, demonstrating effectiveness of the proposed method.

Effect of Variable Gravity Field on Dual Component Convection in a Couple Stress Fluid Saturated Anisotropic Porous Layer With Temperature‐Dependent Heat Source

ABSTRACTThis study examines how gravity fluctuations, the couple stress parameter, anisotropic parameters, and heat source collectively influence dual‐component convection in the porous layer. The linear analysis is conducted utilizing the normal mode technique. The authors proposed three categories of gravity fluctuation, namely: (a) linear, (b) parabolic, and (c) exponential. Expressions for both stationary and oscillatory Rayleigh numbers are derived using the Galerkin approach. Neutral stability curves for both stationary and oscillatory modes are analyzed, with graphical representations to show the effects of various stability parameters, including gravity fluctuations, couple stress, anisotropy, and heat source. The results show that the mechanical anisotropy parameter and Vadasz number lead to system destabilization, while the couple stress parameter, Lewis number, gravity parameter, solute Rayleigh number, and thermal anisotropy parameter help to stabilize the system. Furthermore, the system is more stable with exponential gravity fluctuations and less stable with parabolic gravity fluctuations. This finding offers insights into thermal convective instability in porous media, impacting applications in geoscience, engineering, and environmental science.

Chaotic vibration control of an axially moving string of multidimensional nonlinear dynamic system with an improved FSMC

A new control approach based on fuzzy sliding mode control (FSMC) is proposed to regulate the chaotic vibration of an axial string. Hamilton’s principle is used to formulate the nonlinear equation of motion of the axial translation string, and the von Kármán equations are used to analyse the geometric nonlinearity. The governing equations are nondimensionalized as partial differential equations and transformed into a nonlinear 3-dimensional system via the third-order Galerkin approach. An active control technique based on the FSMC approach is suggested for the derived dynamic system. By using a recurrent neural network model, we can accurately predict and effectively apply a control strategy to suppress chaotic movements. The necessity of the suggested active control method in the regulation of the nonlinear axial translation string system is proven using different chaotic vibrations. The results show that the study of the chaotic vibrations of axially translating strings requires nonlinear multidimensional dynamic systems of axially moving strings; the validity of the proposed control strategy in controlling the chaotic vibration of axially moving strings in a multidimensional form is demonstrated.

Existence and stability of a weakly damped laminated beam with a nonlinear delay

AbstractA weakly damped laminated beam system with nonlinear time delay is studied. The existence and uniqueness are proved by Faedo–Galerkin approach. We prove that the system is stable under some specific conditions on the weight of the delay and the equal wave speeds of propagation. The general energy decay rate is established by using multiplier method and some properties of convex functions. This decay result is obtained without imposing any restrictive growth assumption on the damping term at the origin. In addition, our result improves and develops some existing results in the literature.

A Higher-Order Theory for Nonlinear Dynamic of an FG Porous Piezoelectric Microtube Exposed to a Periodic Load

This paper investigates the nonlinear dynamic deflection, natural frequency, and wave propagation in functionally graded (FG) porous piezoelectric microscale tubes under periodic load, hygrothermal conditions, and an external electric field. The piezoelectric material used to make the smart microtubes has pores that may be smoothly changed or uniformly distributed over the tube wall. Here, three types of porosity distribution are taken into consideration. The nonlinear motion equations are constructed using a novel shear deformation beam theory and the modified couple stress theory (MCST). The nonlinear motion equations are solved using the fourth-order Runge–Kutta technique and the Galerkin approach. The effects of various geometric parameters, porosity distribution type, porosity factor, periodic load amplitude and frequency, material length scale parameter, moisture, and temperature on the nonlinear dynamic deflection, natural frequency, and wave frequency of FG porous piezoelectric microtubes are explored through a number of parametric investigations.

Computational Study of Shocked V-Shaped N2/SF6 Interface across Varying Mach Numbers

The Mach number effect on the Richtmyer–Meshkov instability (RMI) evolution of the shocked V-shaped N2/SF6 interface is numerically studied in this research. Four distinct Mach numbers are taken into consideration for this purpose: Ms=1.12,1.22,1.42, and 1.62. A two-dimensional space of compressible two-component Euler equations is simulated using a high-order modal discontinuous Galerkin approach to computational simulations. The numerical results show good consistency when compared to the available experimental data. The computational results show that the RMI evolution in the shocked V-shaped N2/SF6 interface is critically dependent on the Mach number. The flow field, interface deformation, intricate wave patterns, inward jet development, and vorticity generation are all strongly impacted by the shock Mach number. As the Mach number increases, the V-shaped interface deforms differently, and the distance between the Mach stem and the triple points varies depending on the Mach number. Compared to lower Mach numbers, higher ones produce larger rolled-up vortex chains. A thorough analysis of the Mach number effect identifies the factors that propel the creation of vorticity during the interaction phase. Moreover, kinetic energy and enstrophy both dramatically rise with increasing Mach number. Lastly, a detailed analysis is carried out to determine how the Mach number affects the temporal variations in the V-shaped interface’s features.

Numerical Study of Shock Wave Interaction with V-Shaped Heavy/Light Interface

This paper investigates numerically the shock wave interaction with a V-shaped heavy/light interface. For numerical simulations, we choose six distinct vertex angles (θ=40∘,60∘,90∘,120∘,150∘, and 170∘), five distinct shock wave strengths (Ms=1.12,1.22,1.30,1.60, and 2.0), and three different Atwood numbers (At=−0.32,−0.77, and −0.87). A two-dimensional space of compressible two-component Euler equations are solved using a third-order modal discontinuous Galerkin approach for the simulations. The present findings demonstrate that the vertex angle has a crucial influence on the shock wave interaction with the V-shaped heavy/light interface. The vertex angle significantly affects the flow field, interface deformation, wave patterns, spike generation, and vorticity production. As the vertex angle decreases, the vorticity production becomes more dominant. A thorough analysis of the vertex angle effect identifies the factors that propel the creation of vorticity during the interaction phase. Notably, smaller vertex angles lead to stronger vorticity generation due to a steeper density gradient, while larger angles result in weaker, more dispersed vorticity and a less complex interaction. Moreover, kinetic energy and enstrophy both dramatically rise with decreasing vortex angles. A detailed analysis is also carried out to analyze the vertex angle effects on the temporal variations of interface features. Finally, the impacts of different Mach and Atwood numbers on the V-shaped interface are briefly presented.

Time-domain simulation of sound propagation through anisotropic porous media for room acoustic applications

Wave-based room acoustic simulations require accurate modeling of sound propagation in the presence of extended (bulk) reacting porous absorbers. The sound absorption in porous media is typically addressed by means of equivalent fluid models that provide frequency-dependent effective material parameters, such as density and bulk modulus. In the time domain, these frequency dependent properties result in the computation of convolution integrals, which is effectively handled via an auxiliary differential equations (ADE) method. In this work, we apply the time-explicit discontinuous Galerkin approach to modelling of sound propagation in anisotropic porous media (with an anisotropic equivalent fluid model) together with the ADE method. We demonstrate the importance of using the extended approach for porous media, where the sound absorption depends not only on the frequency and incident angle, but also on the travelling wave direction.

Discretization of Continuous Spaces using Barycentric Subdivision Method with Metric Space Constraints on the Nearest Neighbour Edges

The Rao-Wilton-Glisson (RWG) is a commonly used basis function in the numerical solution of the electric field integral equation (EFIE) using the Method of Moments (MoM) and Galerkin approach. This method relies on triangular patches to approximate the surface current. Traditionally, barycentric subdivision of a primary triangle into n-sub-triangles has been used with RWG basis function to solve the EFIE using MoM. This paper presents a method of approximating a surface using triangular patches by sub-dividing a primary triangle into (2n-1) subtriangles. which creates a denser mesh than the widely used nine (9) points quadrature method. The structure is approximated by small square patches, which are further sub-divided into two primary triangles. By applying barycentric sub-division, the primary triangles are decomposed into sub-triangles to create a mesh over the surface of the structure. Using graph theory, the triangular meshes are defined by the function, G (V, E), where V and E are the vertices and edges of the triangles in the mesh space, respectively. The connectivity matrix of shared edges is found by imposing a constraint on the edge length using metric spaces. This method approximates a square patch with 32 scalene triangles and shows that it can be used to reconstruct equilateral patches and doubly split double ring into mesh structures.

An exact solution of the lubrication equations for the Oldroyd-B model in a hyperbolic pipe

An exact analytical solution of the lubrication equations for the steady, isothermal, incompressible flow of a viscoelastic Oldroyd-B fluid in a hyperbolic cylindrical contracting pipe is derived. The solution is valid for values of the Deborah number, De, up to order unity (De is defined as the ratio of the longest relaxation time of the polymer to the characteristic residence time of the fluid in the pipe), all values of the ratio of the polymer viscosity to the total viscosity of the fluid, η, and typical values of the contraction ratio, Λ, encountered in experiments and practical applications. It is provided in terms of the streamfunction only and is used in the momentum balance to derive a strongly non-linear ordinary differential equation of second order with unknown a function which corresponds to a modified fluid velocity along the main flow direction. The final equation is solved semi-numerically using a fully spectral (Legendre)-Galerkin approach to resolve the unknown function almost down to machine accuracy. The exact solution for the polymer extra-stresses, which is emphasized is not the full solution of the complete lubrication equations, allows for the derivation of a variety of theoretical expressions for the average pressure-drop along the pipe. In all cases, a decrease in the pressure drop compared to the Newtonian value with increasing De, η and/or Λ is predicted. The differences between the corresponding analytical solution for the planar geometrical configuration are also identified and discussed.

Global existence and blow-up of solutions for mixed local and nonlocal hyperbolic equations

In this paper, we consider the following mixed local and nonlocal hyperbolic equation: u t t − Δ u + μ ( − Δ ) s u = | u | p − 2 u , in Ω × R + , u ( x , 0 ) = u 0 ( x ) , u t ( x , 0 ) = u 1 ( x ) , in Ω , u ( x , t ) = 0 , in ( R N ∖ Ω ) × R 0 + , where s ∈ ( 0 , 1 ), N > 2, p ∈ ( 2 , 2 s ∗ ], μ is a nonnegative real parameter, Ω ⊂ R N is a bounded domain with Lipschitz boundary ∂ Ω, Δ is the Laplace operator, ( − Δ ) s is the fractional Laplace operator. By combining the Galerkin approach with the modified potential well method, we obtain the global existence, vacuum isolating, and blow-up of solutions for the aforementioned problem, provided certain assumptions are fulfilled. Specifically, we study the existence of global solutions for the above problem in the cases of subcritical and critical initial energy levels, as well as the finite time blow-up of solutions. Then, we investigate the blow-up of solutions for the above problem in the case of supercritical initial energy level, as well as upper and lower bounds of blow-up time of solutions.

H(div)-conforming and discontinuous Galerkin approach for Herschel–Bulkley flow with density-dependent viscosity and yield stress

This paper presents a comprehensive study on Herschel–Bulkley flow, where the flow parameters are dependent on the density. The Herschel–Bulkley model is a generalized power-law model used to simulate viscoplastic fluids defined by a plasticity threshold. We consider the case where the plasticity threshold and the viscosity depend on the shear rate and fluid density. To analyze this model, we use a Huber regularization of the stress and propose an H(div)-conforming and discontinuous Galerkin (DG) numerical approximation for the coupled equations governing the flow. We discuss the stability and existence of discrete solutions and propose a semismooth Newton linearization for the numerical solution of the discretized system. Our numerical scheme is validated through several experiments that explore the behavior of Herschel–Bulkley flow under different conditions. The results demonstrate the robustness of our numerical method.

A novel multi-frame image super-resolution model based on regularized nonlinear diffusion with Caputo time fractional derivative

In this work, we introduce an innovative fractional nonlinear parabolic model using a time-fractional order derivative, specifically employing the Caputo sense for fractional differentiation. This model aims to enhance traditional super-resolution models, particularly in the context of multi-frame image super-resolution. Additionally, we incorporate a regularized Perona–Malik diffusion mechanism to control the speed and direction of diffusion at each image location. We begin our study by exploring the theoretical solvability of our proposed model. Firstly, we employ the Faedo–Galerkin approach to establish the existence and uniqueness of a weak solution for an auxiliary fractional super-resolution model. Subsequently, we use the Schauder fixed point method to demonstrate the existence and uniqueness of a weak solution for our model. To validate the effectiveness of our model in the multi-frame super-resolution (SR) context, we conduct numerical experiments on images featuring diverse characteristics, including corners and edges, while applying various warping, decimation, and blurring matrices to the low-resolution (LR) images. We start the evaluation by introducing an adaptive discrete scheme tailored to the proposed model. To prove the robustness of our approach, we subject our images to varying levels of noise. Additionally, we perform simulations on real data (videos). The obtained high-resolution (HR) results demonstrate notable efficiency and robustness against noise, outperforming competitive models both visually and quantitatively.

Finite Element Approach for Rheological Behavior in Colloidal Electrolytes in Lithium-Ion Battery Performance.

The main implication of articulating electrolyte performance is studying the energy density, charging aspects, formation of precipitates, thermal fluctuations during charging-discharging, and safety of batteries against fire or spark. One of the most significant aspects is the ability to design colloidal electrolytes that can enhance the overall performance of batteries along with dealing with all internal problems within a battery system. Through this optimization progression, the general performance and efficiency of Li-ion batteries can be improved. This work is presented in the study of the boundary value problem for rheological properties of colloidal electrolytes as a fourth grade fluid for lithium ion (Li-ion) batteries down a vertical cylinder. They have exceptional characteristics, such as low volatility and high thermal stability. The practical usage of the exact flow is restricted, as it involves very complicated integrals. The nonlinear problem that arises is solved by Galerkin's finite element approach based on the weighted-residual formulation, which is used to find the approximate solutions of the fourth-grade problem. This approach utilizes a piecewise linear approximation using linear Lagrange polynomials. Convergence of the solutions, which briefly describes the flow characteristics, includes the effects of the emerging parameters. The results obtained after implementation are not restrictive to small values of the flow parameters. Numerical studies have shown the superior accuracy and lesser computational cost of this scheme in comparison to collocation, the homotopy analysis method, and the homotopy perturbation method. The impact of the relevant parameters is examined through graphical results after implementation of a number of iterations.

Dynamics of axially moving spinning microbeams composed of tri-directional graded porous materials with axisymmetric cross-sections and piezoelectric layers in complex fields

The current article appraises the vibration and stability of tri-directional functionally graded porous microscale beams with rectangular cross-sections integrated with piezoelectric layers under spinning and axial movements in complex environments. The microbeam is surrounded by a three-parameter Winkler-Pasternak-Hetenyi medium, and its material characteristics are graded in thickness, width, and longitudinal spatial directions by considering non-uniform and uniform porosity models. Dynamic equations, vibration frequencies, and stability criteria of the system are determined with the aid of the Galerkin approach and Laplace transform. The Campbell diagram and stability maps are drawn. Frequency and stability analyses, as well as comparison and parametric analyses, are conducted. The impacts of piezoelectric voltage, magneto-hygro-thermal fields, axial and tangential distributed follower forces, substrate characteristics, scale parameter, aspect ratio, porosity factor, and material gradation on flutter and divergence instability boundaries are assessed in detail. It is deduced that instability regions are condensed, and the instability threshold is enhanced by fine-adjusting the porosity and material gradient. It is discovered that destructive environmental effects can be alleviated by regulating the piezoelectric voltage. In addition, compared with the case of a square cross-section, the divergence/flutter instability region of the microbeam with a rectangular cross-section is smaller/larger. The outcomes of the present research can be helpful in the design of next-generation bi-gyroscopic systems.

Dynamical behavior of fluid–structure interaction in ducts with rigid and flexible interfaces: Modeling and analysis

This study investigates the propagation of fluid–structure coupled waves in a duct featuring both rigid and flexible interfaces. Two cases are considered: in the first case, with a rigid interface, the truncated form of solution preserves the conservation law, while in the second setting, the power balance is not preserved initially, but is rectified as the number of modes, denoted as N, approaches infinity. The discrepancy is explained by specific equations and the presence of delta functions to enforce edge conditions on the membrane at the interface. The Galerkin approach is adopted to address this issue, expressing the membrane’s displacement using a generalized Fourier series. The investigation establishes a framework for modeling higher-order boundary conditions within a duct, employing the mode-matching (MM) method. The study also delves into the scattering components of power, which vary significantly for the two cases. These findings have implications for noise control devices, particularly in active noise control measures utilizing flexible components to dampen both noise and structural vibrations.