C1-Conforming Quadrilateral Spectral Element Method for Fourth-Order Equations

This paper is devoted to Professor Benyu Guo's open question on the C1-conforming quadrilateral spectral element method for fourth-order equations which has been endeavored for years. Starting with generalized Jacobi polynomials on the reference square, we construct the C1-conforming basis functions using the bilinear mapping from the reference square onto each quadrilateral element which fall into three categories-interior modes, edge modes, and vertex modes. In contrast to the triangular element, compulsively compensatory requirements on the global C1-continuity should be imposed for edge and vertex mode basis functions such that their normal derivatives on each common edge are reduced from rational functions to polynomials, which depend on only parameters of the common edge. It is amazing that the C1-conforming basis functions on each quadrilateral element contain polynomials in primitive variables, the completeness is then guaranteed and further confirmed by the numerical results on the Petrov-Galerkin spectral method for the non-homogeneous boundary value problem of fourth-order equations on an arbitrary quadrilateral. Finally, a C1-conforming quadrilateral spectral element method is proposed for the biharmonic eigenvalue problem, and numerical experiments demonstrate the effectiveness and efficiency of our spectral element method.

A Fourth-Order Unstructured NURBS-Enhanced Finite Volume WENO Scheme for Steady Euler Equations in Curved Geometries

In Li and Ren(Int.J.Numer.Methods Fluids 70:742–763,2012),a high-order k-exact WENO finite volume scheme based on secondary reconstructions was proposed to solve the two-dimensional time-dependent Euler equations in a polygonal domain,in which the high-order numerical accuracy and the oscillations-free property can be achieved.In this paper,the method is extended to solve steady state problems imposed in a curved physical domain.The numerical framework consists of a Newton type finite volume method to linearize the nonlinear governing equations,and a geometrical multigrid method to solve the derived linear system.To achieve high-order non-oscillatory numerical solutions,the classical k-exact reconstruction with k=3 and the efficient secondary reconstructions are used to perform the WENO reconstruction for the conservative variables.The non-uniform rational B-splines(NURBS)curve is used to provide an exact or a high-order representation of the curved wall boundary.Furthermore,an enlarged reconstruction patch is constructed for every element of mesh to significantly improve the convergence to steady state.A variety of numerical examples are presented to show the effectiveness and robustness of the proposed method.

Fourth-order phase-field modeling for brittle fracture in piezoelectric materials

Failure analyses of piezoelectric structures and devices are of engineering and scientific significance.In this paper,a fourth-order phase-field fracture model for piezoelectric solids is developed based on the Hamilton principle.Three typical electric boundary conditions are involved in the present model to characterize the fracture behaviors in various physical situations.A staggered algorithm is used to simulate the crack propagation.The polynomial splines over hierarchical T-meshes(PHT-splines)are adopted as the basis function,which owns the C1continuity.Systematic numerical simulations are performed to study the influence of the electric boundary conditions and the applied electric field on the fracture behaviors of piezoelectric materials.The electric boundary conditions may influence crack paths and fracture loads significantly.The present research may be helpful for the reliability evaluation of the piezoelectric structure in the future applications.

Some Modified Equations of the Sine-Hilbert Type

Three modified sine-Hilbert(sH)-type equations, i.e., the modified sH equation, the modified damped sH equation, and the modified nonlinear dissipative system, are proposed, and their bilinear forms are provided.Based on these bilinear equations, some exact solutions to the three modified equations are derived.

An Extended Numerical Method by Stancu Polynomials for Solution of Integro-Differential Equations Arising in Oscillating Magnetic Fields

In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.

Theoretical study of particle and energy balance equations in locally bounded plasmas

In this study,new particle and energy balance equations have been developed to predict the electron temperature and density in locally bounded plasmas.Classical particle and energy balance equations assume that all plasma within a reactor is completely confined only by the reactor walls.However,in industrial plasma reactors for semiconductor manufacturing,the plasma is partially confined by internal reactor structures.We predict the effect of the open boundary area(A′_(L,eff))and ion escape velocity(u_(i))on electron temperature and density by developing new particle and energy balance equations.Theoretically,we found a low ion escape velocity(u_(i)/u_(B)≈0.2)and high open boundary area(A′_(L,eff)/A_(T,eff)≈0.6)to result in an approximately 38%increase in electron density and an 8%decrease in electron temperature compared to values in a fully bounded reactor.Additionally,we suggest that the velocity of ions passing through the open boundary should exceedω_(pi)λ_(De)under the condition E^(2)_(0)?(Φ/λ_(De))^(2).

Matrix Riccati Equations in Optimal Control

In this paper, the matrix Riccati equation is considered. There is no general way for solving the matrix Riccati equation despite the many fields to which it applies. While scalar Riccati equation has been studied thoroughly, matrix Riccati equation of which scalar Riccati equations is a particular case, is much less investigated. This article proposes a change of variable that allows to find explicit solution of the Matrix Riccati equation. We then apply this solution to Optimal Control.

Analytical solutions fractional order partial differential equations arising in fluid dynamics

This article describes the solution procedure of the fractional Pade-Ⅱ equation and generalized Zakharov equation(GSEs)using the sine-cosine method.Pade-Ⅱ is an important nonlinear wave equation modeling unidirectional propagation of long-wave in dispersive media and GSEs are used to model the interaction between one-dimensional high,and low-frequency waves.Classes of trigonometric and hyperbolic function solutions in fractional calculus are discussed.Graphical simulations of the numerical solutions are flaunted by MATLAB.

THE STABILITY OF BOUSSINESQ EQUATIONS WITH PARTIAL DISSIPATION AROUND THE HYDROSTATIC BALANCE

This paper is devoted to understanding the stability of perturbations around the hydrostatic equilibrium of the Boussinesq system in order to gain insight into certain atmospheric and oceanographic phenomena.The Boussinesq system focused on here is anisotropic,and involves only horizontal dissipation and thermal damping.In the 2D case R^(2),due to the lack of vertical dissipation,the stability and large-time behavior problems have remained open in a Sobolev setting.For the spatial domain T×R,this paper solves the stability problem and gives the precise large-time behavior of the perturbation.By decomposing the velocity u and temperatureθinto the horizontal average(ū,θ)and the corresponding oscillation(ū,θ),we can derive the global stability in H~2 and the exponential decay of(ū,θ)to zero in H^(1).Moreover,we also obtain that(ū_(2),θ)decays exponentially to zero in H^(1),and thatū_(1)decays exponentially toū_(1)(∞)in H^(1)as well;this reflects a strongly stratified phenomenon of buoyancy-driven fluids.In addition,we establish the global stability in H^(3)for the 3D case R^(3).

The Maxwell-Heaviside Equations Explained by the Theory of Informatons

In the articles “Newtons Law of Universal Gravitation Explained by the Theory of Informatons” and “The Gravitational Interaction between Moving Mass Particles Explained by the Theory of Informatons” the gravitational interaction has been explained by the hypothesis that information carried by informatons is the substance of gravitational fields, i.e. the medium that the interaction in question makes possible. From the idea that “information carried by informatons” is its substance, it has been deduced that—on the macroscopic level—a gravitational field manifests itself as a dual entity, always having a field- and an induction component (Egand Bg) simultaneously created by their common sources. In this article we will mathematically deduce the Maxwell-Heaviside equations from the kinematics of the informatons. These relations describe on the macroscopic level how a gravitational field (Eg, Bg) is generated by whether or not moving masses and how spatial and temporal changes of Egand Bgare related. We show that there is no causal link between Egand Bg.

ALMOST SURE GLOBAL WELL-POSEDNESS FOR THE FOURTH-ORDER NONLINEAR SCHR?DINGER EQUATION WITH LARGE INITIAL DATA

We consider the fourth-order nonlinear Schr?dinger equation(4NLS)(i?t+εΔ+Δ2)u=c1um+c2(?u)um-1+c3(?u)2um-2,and establish the conditional almost sure global well-posedness for random initial data in Hs(Rd)for s∈(sc-1/2,sc],when d≥3 and m≥5,where sc:=d/2-2/(m-1)is the scaling critical regularity of 4NLS with the second order derivative nonlinearities.Our proof relies on the nonlinear estimates in a new M-norm and the stability theory in the probabilistic setting.Similar supercritical global well-posedness results also hold for d=2,m≥4 and d≥3,3≤m<5.

THE SMOOTHING EFFECT IN SHARP GEVREY SPACE FOR THE SPATIALLY HOMOGENEOUS NON-CUTOFF BOLTZMANN EQUATIONS WITH A HARDPOTENTIAL

In this article, we study the smoothing effect of the Cauchy problem for the spatially homogeneous non-cutoff Boltzmann equation for hard potentials. It has long been suspected that the non-cutoff Boltzmann equation enjoys similar regularity properties as to whose of the fractional heat equation. We prove that any solution with mild regularity will become smooth in Gevrey class at positive time, with a sharp Gevrey index, depending on the angular singularity. Our proof relies on the elementary L^(2) weighted estimates.

Besov Estimates for Sub-Elliptic Equations in the Heisenberg Group

In this article, we deal with weak solutions to non-degenerate sub-elliptic equations in the Heisenberg group, and study the regularities of solutions. We establish horizontal Calderón-Zygmund type estimate in Besov spaces with more general assumptions on coefficients for both homogeneous equations and non-homogeneous equations. This study of regularity estimates expands the Calderón-Zygmund theory in the Heisenberg group.

On entire solutions of some Fermat type differential-difference equations

On one hand,we study the existence of transcendental entire solutions with finite order of the Fermat type difference equations.On the other hand,we also investigate the existence and growth of solutions of nonlinear differential-difference equations.These results extend and improve some previous in[5,14].

ELLIPTIC EQUATIONS IN DIVERGENCE FORM WITH DISCONTINUOUS COEFFICIENTS IN DOMAINS WITH CORNERS

We study equations in divergence form with piecewise Cαcoefficients.The domains contain corners and the discontinuity surfaces are attached to the edges of the corners.We obtain piecewise C^(1,α) estimates across the discontinuity surfaces and provide an example to illustrate the issue regarding the regularity at the corners.

ON THE STABILITY OF PERIODIC SOLUTIONS OF PIECEWISE SMOOTH PERIODIC DIFFERENTIAL EQUATIONS

In this paper,we address the stability of periodic solutions of piecewise smooth periodic differential equations.By studying the Poincarémap,we give a sufficient condition to judge the stability of a periodic solution.We also present examples of some applications.

On Two Types of Stability of Solutions to a Pair of Damped Coupled Nonlinear Evolution Equations

The stability of a set of spatially constant plane wave solutions to a pair of damped coupled nonlinear Schrödinger evolution equations is considered. The equations could model physical phenomena arising in fluid dynamics, fibre optics or electron plasmas. The main result is that any small perturbation to the solution remains small for all time. Here small is interpreted as being both in the supremum sense and the square integrable sense.

Sparse-Grid Implementation of Fixed-Point Fast Sweeping WENO Schemes for Eikonal Equations

Fixed-point fast sweeping methods are a class of explicit iterative methods developed in the literature to efficiently solve steady-state solutions of hyperbolic partial differential equations(PDEs).As other types of fast sweeping schemes,fixed-point fast sweeping methods use the Gauss-Seidel iterations and alternating sweeping strategy to cover characteristics of hyperbolic PDEs in a certain direction simultaneously in each sweeping order.The resulting iterative schemes have a fast convergence rate to steady-state solutions.Moreover,an advantage of fixed-point fast sweeping methods over other types of fast sweeping methods is that they are explicit and do not involve the inverse operation of any nonlinear local system.Hence,they are robust and flexible,and have been combined with high-order accurate weighted essentially non-oscillatory(WENO)schemes to solve various hyperbolic PDEs in the literature.For multidimensional nonlinear problems,high-order fixed-point fast sweeping WENO methods still require quite a large amount of computational costs.In this technical note,we apply sparse-grid techniques,an effective approximation tool for multidimensional problems,to fixed-point fast sweeping WENO methods for reducing their computational costs.Here,we focus on fixed-point fast sweeping WENO schemes with third-order accuracy(Zhang et al.2006[41]),for solving Eikonal equations,an important class of static Hamilton-Jacobi(H-J)equations.Numerical experiments on solving multidimensional Eikonal equations and a more general static H-J equation are performed to show that the sparse-grid computations of the fixed-point fast sweeping WENO schemes achieve large savings of CPU times on refined meshes,and at the same time maintain comparable accuracy and resolution with those on corresponding regular single grids.

THE GLOBAL EXISTENCE OF STRONG SOLUTIONS TO THERMOMECHANICAL CUCKER-SMALE-STOKES EQUATIONS IN THE WHOLE DOMAIN

We study the global existence and uniqueness of a strong solution to the kinetic thermomechanical Cucker-Smale(for short,TCS) model coupled with Stokes equations in the whole space.The coupled system consists of the kinetic TCS equation for a particle ensemble and the Stokes equations for a fluid via a drag force.In this paper,we present a complete analysis of the existence of global-in-time strong solutions to the coupled model without any smallness restrictions on the initial data.

Multi-soliton solutions of coupled Lakshmanan-Porsezian-Daniel equations with variable coefficients under nonzero boundary conditions

This paper aims to investigate the multi-soliton solutions of the coupled Lakshmanan–Porsezian–Daniel equations with variable coefficients under nonzero boundary conditions.These equations are utilized to model the phenomenon of nonlinear waves propagating simultaneously in non-uniform optical fibers.By analyzing the Lax pair and the Riemann–Hilbert problem,we aim to provide a comprehensive understanding of the dynamics and interactions of solitons of this system.Furthermore,we study the impacts of group velocity dispersion or the fourth-order dispersion on soliton behaviors.Through appropriate parameter selections,we observe various nonlinear phenomena,including the disappearance of solitons after interaction and their transformation into breather-like solitons,as well as the propagation of breathers with variable periodicity and interactions between solitons with variable periodicities.