Galerkin Method Research Articles

Superconvergent scheme for a system of green Fredholm integral equations

In this study, a numerical scheme to a system of second-kind linear Fredholm integral equations featuring a Green's kernel function is proposed. This involves introducing Galerkin and iterated Galerkin (IG) methods based on piecewise polynomials to tackle the integral model. A thorough analysis of convergence and error for these proposed methods is carried out. Firstly, the existence and uniqueness of solutions for the Galerkin and iterated Galerkin methods are established. Later, the order of convergence is derived using tools from functional analysis and the boundedness property of Green's kernel. The Galerkin scheme has O(hα) order of convergence. Next, the superconvergence of the iterated Galerkin (IG) method is established. The IG method exhibits O(hα+α⁎) order of convergence. Theoretical findings are validated through extensive numerical experiments.

Augmented subspace method for semilinear parabolic problems based on symmetric interior penalty discontinuous Galerkin method

Augmented subspace method for semilinear parabolic problems based on symmetric interior penalty discontinuous Galerkin method

An exponential spectral deferred correction method for multidimensional parabolic problems

We present some efficient algorithms based on an exponential time differencing spectral deferred correction (ETDSDC) method for multidimensional second and fourth-order parabolic problems with non-periodic boundary conditions including Dirichlet, Neumann, Robin boundary conditions. Similar to the Fourier spectral method for periodic problems, the key to the efficiency of our algorithms is to construct diagonal discrete linear operators via Legendre–Galerkin methods with Fourier-like basis functions. In combination with the ETDSDC scheme, the proposed methods are spectrally accurate in space and up to 10th-order accurate in time (as shown in this work). We demonstrate the high-order of convergence and efficiency of our algorithms in solving parabolic equations through a series of two-dimensional and three-dimensional examples including Ginzburg–Landau and Allen–Cahn equations.

Structure-preserving finite element methods for computing dynamics of rotating Bose-Einstein condensate

This work is concerned with the construction and analysis of structure-preserving Galerkin methods for computing the dynamics of rotating Bose-Einstein condensate (BEC) based on the Gross-Pitaevskii equation with angular momentum rotation. Due to the presence of the rotation term, constructing finite element methods (FEMs) that preserve both mass and energy remains an unresolved issue, particularly in the context of nonconforming FEMs. Furthermore, in comparison to existing works, we provide a comprehensive convergence analysis, offering a thorough demonstration of the methods’ optimal and high-order convergence properties. Finally, extensive numerical results are presented to check the theoretical analysis of the structure-preserving numerical method for rotating BEC, and the quantized vortex lattice’s behavior is scrutinized through a series of numerical tests.

On the existence, uniqueness and regularity of strong solutions to a stochastic 2D Cahn–Hilliard-Magnetohydrodynamic model

Abstract We investigate a stochastic coupled model of the Cahn–Hilliard equations and the stochastic magnetohydrodynamic equations in a bounded domain of ℝ 2 {\mathbb{R}^{2}} . The model describes the flow of the mixture of two incompressible and immiscible fluids under the influence of an electromagnetic field with stochastic perturbations. We prove the existence, uniqueness and regularity of a probabilistic strong solution. The proof of the existence is based on the Galerkin approximation, the stopping time technique and some weak convergence principles in functional analysis.

Explicit exponential Runge–Kutta methods for semilinear time-fractional integro-differential equations

In this work, we consider and analyze explicit exponential Runge–Kutta methods for solving semilinear time-fractional integro-differential equation, which involves two nonlocal terms in time. Firstly, the temporal Runge–Kutta discretizations follow the idea of exponential integrators. Subsequently, we utilize the spectral Galerkin method to introduce a fully discrete scheme. Then, we mainly focus on discussing the one-stage and two-stage methods for solving the proposed semilinear problem. Based on special abstract settings, we perform the convergence analysis for the proposed two different stage methods. In this process, we heavily use estimates about the operator family {S̃(t)}, and in combination with Lipschitz continuous condition. Finally, some numerical experiments confirm theoretical results. Meanwhile, applying this scheme to the related linear problem yields high-order convergence, highlighting the advantages of explicit exponential Runge–Kutta methods.

Strongly quasilinear elliptic systems with perturbations through yong measures

We use the Galerkin method to obtain weak solutions for strongly quasilinear elliptic systems in the Sobolev space W1,p0(Ω; ℝm) of the form Gϖ(x) = f(x,ϖ,Dϖ) in Ω, with ϖ(x) = 0 on ∂Ω, where the term Gϖ(x) includes perturbations. Our approach involves the use of Young measures to obtain the desired results.

Chaos and resonance of nonlinear FG shallow shells reinforced by oblique stiffeners

This paper presents an analysis of the nonlinear dynamics of functionally graded (FG) shallow shells that are simply supported, reinforced with oblique stiffeners, and subjected to external excitation. The properties of these oblique stiffened functionally graded (OSFG) shallow shells vary across their thickness, following a power law distribution based on the volume fractions. Employing the first-order shear deformation theory (FSDT) coupled with the von Kármán nonlinear strain-displacement relations, the derivation of the nonlinear equations of motion is accomplished through the application of Hamilton's principle. The model of FG shallow shells reinforced with oblique two-series stiffeners is developed, allowing for the angles of the two stiffeners to be similar or different from each other. In this context, the employment of the first-order shear deformation theory (FSDT) results in the enhancement of the complexity in modeling the oblique stiffeners. Galerkin's method is applied to discretize the given equations, converting them into a nonlinear dynamical system with two degrees of freedom. This system incorporates both quadratic and cubic nonlinearities and considers the effects of external excitation. The analysis of this research particularly focuses on the resonant scenario involving principal resonance and 1:1 internal resonance. Through the asymptotic perturbation method, a four-dimensional nonlinear averaged equation is obtained. For the first time, a numerical investigation of the research reveals critical insights into the behavior of the system, such as waveforms, phase diagrams, and Poincaré maps, highlighting the effects of different angles of the stiffeners and variations of material parameters on the nonlinear dynamics of these shells.

Solving linear elasticity benchmark problems via the overset improved element-free Galerkin-finite element method

A novel approach for the solution of linear elasticity problems is introduced in this communication, which uses a hybrid chimera-type technique based on both finite element and improved element-free Galerkin methods. The proposed overset improved element-free Galerkin-finite element method (Ov-IEFG-FEM) for linear elasticity uses the finite element method (FEM) throughout the entire problem geometry, whereas a fine distribution of overlapping nodes is used to perform higher order approximations via the improved element-free Galerkin (IEFG) technique in regions demanding more computational accuracy. The method relies on keeping the FEM-based results in those regions where low order of approximation is enough to provide the required accuracy, i.e. outside the region where the solution will be enriched via the IEFG technique. The overlapping domains perform an iterative transfer of kinematics information through well-defined immersed boundaries, and a detailed explanation on this regard is also presented in this communication. The Ov-IEFG-FEM is used in a set of increasingly complex linear elasticity problems, and the outcomes demonstrate the suitability and reliability of this technique to solve such problems in an accurate and remarkably simple manner.

A hybrid interpolating element-free Galerkin method for 3D steady-state convection diffusion problems

This study investigates a hybrid interpolating element-free Galerkin (HIEFG) method for solving 3D convection diffusion problems. The HIEFG approach divides a 3D solution domain into a sequence of interconnected 2D sub-domains, and in these 2D sub-domains, the interpolating element-free Galerkin (IEFG) method is applied to form the discretized equations. The improved interpolating moving least-squares (IMLS) method is used to obtain the shape function of the IEFG method for 2D problems. The finite difference method is employed to combine the discretized equations in 2D sub-domains in the splitting direction. Then, the HIEFG method's formulas are derived for steady-state convection diffusion problems in 3D solution domain. Three numerical examples are used to discuss the impacts of the number of nodes, the number of split layers, and the scaling parameters of the influence domain on the computational precision and CPU time of the HIEFG technique. Imposing boundary conditions directly and the dimension splitting technique in this method significantly improves the computational speed greatly.

Two-route formation of soot nuclei: experimental and modeling evidence

The results of shock-tube measurements of acetylene concentration and characteristics of the ensemble of soot particles during the pyrolysis of benzene, ethylene, ethylene–methane and ethylene–propane diluted with argon are compared with the predictions of the unified kinetic model of soot formation implemented in the MACRON program based on the Galerkin method. The simulation results demonstrate good agreement with experimental data for the pyrolysis of ethylene alone and with the specified additives, thereby confirming the two route mechanism of soot nuclei formation. The agreement for benzene pyrolysis is less satisfactory.

A bound preserving cut discontinuous Galerkin method for one dimensional hyperbolic conservation laws

In this paper we present a family of high order cut finite element methods with bound preserving properties for hyperbolic conservation laws in one space dimension. The methods are based on the discontinuous Galerkin framework and use a regular background mesh, where interior boundaries are allowed to cut through the mesh arbitrarily. Our methods include ghost penalty stabilization to handle small cut elements and a new reconstruction of the approximation on macro-elements, which are local patches consisting of cut and un-cut neighboring elements that are connected by stabilization. We show that the reconstructed solution retains conservation and order of convergence. Our lowest order scheme results in a piecewise constant solution that satisfies a maximum principle for scalar hyperbolic conservation laws. When the lowest order scheme is applied to the Euler equations, the scheme is positivity preserving in the sense that positivity of pressure and density are retained. For the high order schemes, suitable bound preserving limiters are applied to the reconstructed solution on macro-elements. In the scalar case, a maximum principle limiter is applied, which ensures that the limited approximation satisfies the maximum principle. Correspondingly, we use a positivity preserving limiter for the Euler equations, and show that our scheme is positivity preserving. In the presence of shocks additional limiting is needed to avoid oscillations, hence we apply a standard TVB limiter to the reconstructed solution. The time step restrictions are of the same order as for the corresponding discontinuous Galerkin methods on the background mesh. Numerical computations illustrate accuracy, bound preservation, and shock capturing capabilities of the proposed schemes.

An ultraweak-local discontinuous Galerkin method for nonlinear biharmonic Schrödinger equations

This paper proposes and analyzes a fully discrete scheme for nonlinear biharmonic Schrödinger equations. We first write the single equation into a system of problems with second-order spatial derivatives and then discretize the space variable with an ultraweak discontinuous Galerkin scheme and the time variable with the Crank–Nicolson method. The proposed scheme proves to be computationally more efficient compared to the local discontinuous Galerkin method in terms of the number of equations needed to be solved at each single time step, and it is unconditionally stable without imposing any penalty terms. It also achieves optimal error convergence in L2 norm both in the solution and in the auxiliary variable with general nonlinear terms. We also prove several physically relevant properties of the discrete schemes, such as the conservation of mass and the Hamiltonian for the nonlinear biharmonic Schrödinger equations. Several numerical studies demonstrate and support our theoretical results.

Thermal effect on the flow induced by a single-dielectric-barrier-discharge plasma actuator under steady actuation

The thermal effect of a single-dielectric-barrier-discharge plasma actuator under steady actuation is numerically investigated. A new actuator model is proposed and validated using experimental data. A discrete Galerkin method based on high-order flux reconstruction schemes is employed to solve the flow governing equations and the actuator model equations on unstructured quadrilateral grids. By comparing the induced heated and cold flow fields of the actuator with and without a plasma thermal source, its thermal effect is revealed. The actuator generates a thermal wall jet with rich vorticity, forming a monopolar starting vortex with a high-temperature and low-density core. Over time, the starting vortex becomes unstable and transforms into a dipole. Actuator heating enhances jet velocity and width, as well as vortex stability, while slowing down vorticity generation. The relative change in density and temperature fields due to actuator heating is four orders of magnitude greater than that without actuator heating. Additionally, the actuator heating causes the background thermodynamic fields to increase approximately linearly with time. Two stages in the actuator's thermal effect are distinguished due to time accumulation. Initially, the actuator heating minimally affects the monopolar starting vortex motion, and the temperature and density fields are treated as passive variables driven by the velocity field. During this stage, the momentum and thermal effects of the actuator can be studied separately. However, after the starting vortex becomes unstable, the actuator heating significantly impacts its motion and morphology, and these two effects are coupled with each other.

Denoising graph neural network based on zero-shot learning for Gibbs phenomenon in high-order DG applications

With the availability of high-performance computing technology and the development of advanced numerical simulation methods, Computational Fluid Dynamics (CFD) is becoming more and more practical and efficient in engineering. As one of the high-precision representative algorithms, the high-order Discontinuous Galerkin Method (DGM) has not only attracted widespread attention from scholars in the CFD research community, but also received strong development. However, when DGM is extended to high-speed aerodynamic flow field calculations, non-physical numerical Gibbs oscillations near shock waves often significantly affect the numerical accuracy and even cause calculation failure. Data driven approaches based on machine learning techniques can be used to learn the characteristics of Gibbs noise, which motivates us to use it in high-speed DG applications. To achieve this goal, labeled data need to be generated in order to train the machine learning models. This paper proposes a new method for denoising modeling of Gibbs phenomenon using a machine learning technique, the zero-shot learning strategy, to eliminate acquiring large amounts of CFD data. The model adopts a graph convolutional network combined with graph attention mechanism to learn the denoising paradigm from synthetic Gibbs noise data and generalize to DGM numerical simulation data. Numerical simulation results show that the Gibbs denoising model proposed in this paper can suppress the numerical oscillation near shock waves in the high-order DGM. Our work automates the extension of DGM to high-speed aerodynamic flow field calculations with higher generalization and lower cost.

A Weak Galerkin Method for the Coupled Darcy-Stokes Problem with a Nonstandard Transmission Condition on the Stokes Boundary

A Weak Galerkin Method for the Coupled Darcy-Stokes Problem with a Nonstandard Transmission Condition on the Stokes Boundary

Cartesian mesh adaptation: Immersed boundary method based on high-order discontinuous Galerkin method

Immersed boundary method (IBM) can easily distinguish fluid and solid regions in the computational region, thereby the workload of complex grid generation can be reduced. To accurately characterize the solid geometry, a large number of cells are required near the solid surface. The h-adaptive algorithm is adopted to reduce the requirement for the number of cells. In addition, considering the inherent adaptability to the h-adaptive Cartesian grids of the discontinuous Galerkin method, a high-order discontinuous Galerkin solver with an IBM is developed. To validate the h-adaptive algorithm and the solver, three cases are tested, including the steady flow past the National Advisory Committee for Aeronautics 0012 airfoil, the steady flow past a cylinder, and the unsteady flow past a cylinder. Compared with the non-adaptive cases, the h-adaptive cases need smaller total number of cells, and the numerical accuracy is significantly improved with an increasing degree of mesh refinement.

Nonlinear free vibration analysis of the rectangular conductive elastic plate in magnetic field based on homotopy perturbation method

AbstractThis article tries to investigate the nonlinear free vibrations of the rectangular conductive elastic plate in uniform magnetic fields under the classic plate theory considering nonlinear strain‐displacement. The formulation of the governing equations integrates the electromagnetic forces arising from the motion of the plate. The nonlinear motion equations are dimensionless, which takes the effect of in‐plane inertia into account. The equations are solved by the Galerkin method and homotopy perturbation method (HPM). The effectiveness of the solution is verified. According to the solutions by HPM, the effects of the initial conditions, length‐to‐thickness ratio, and magnetic induction intensity on the nonlinear free vibrations behavior of conductive elastic plates are discussed.

General non-linear fragmentation with discontinuous Galerkin methods

ABSTRACT Dust grains play a significant role in several astrophysical processes, including gas/dust dynamics, chemical reactions, and radiative transfer. Replenishment of small-grain populations is mainly governed by fragmentation during pair-wise collisions between grains. The wide spectrum of fragmentation outcomes, from complete disruption to erosion and/or mass transfer, can be modelled by the general non-linear fragmentation equation. Efficiently solving this equation is crucial for an accurate treatment of the dust fragmentation in numerical modelling. However, similar to dust coagulation, numerical errors in current fragmentation algorithms employed in astrophysics are dominated by the numerical overdiffusion problem – particularly in three-dimensional hydrodynamic simulations where the discrete resolution of the mass-density distribution tends to be highly limited. With this in mind, we have derived the first conservative form of the general non-linear fragmentation with a mass flux highlighting the mass transfer phenomenon. Then, to address cases of limited mass density resolution, we applied a high-order discontinuous Galerkin scheme to efficiently solve the conservative fragmentation equation with a reduced number of dust bins. An accuracy of $0.1{\!-\!}1~{{\ \rm per\ cent}}$ is reached with 20 dust bins spanning a mass range of 9 orders of magnitude.

Dynamic response of a pipe conveying pulsating fluid under random vibration

This paper studies the dynamic response of a pipe conveying pulsating fluid under random vibration. Its nonlinear equation of motion is first developed using the method of Newton. The method of Galerkin is then used to discretize the equation of motion into a system of ordinary differential equations. Afterwards, the statistical moments and mean peak values of the displacement and velocity responses of the pipe are obtained using an efficient Monte Carlo simulation (MCS) method based on the explicit time-domain method (ETDM). Furthermore, the displacement-sensitive and velocity-sensitive locations of the pipe are obtained for different boundary conditions. The impacts of the fluid, random vibration, and structural parameters on the mean peak value of the dynamic response of the pipe are finally analyzed. The obtained results provide a reference for better understanding the dynamic characteristics of pulsating fluid transport pipes under random vibration and for accurately designing pipe structures.