Secondary Buckling Analysis of Thin Rectangular Plates Based on the Wavelet Galerkin Method北大核心CSCD

Application of the wavelet Galerkin method (WGM) to numerical solution of nonlinear buckling problems was studied with classical elastic thin rectangular plates. First, the discretized scheme of the von Kármán equation were introduced, then a simple calculation approach to the Jacobian and Hessian matrices based on the WGM was proposed, and the wavelet discretized scheme-based eigenvalue equation method, the extended equation method and the pseudo arc-length method for nonlinear buckling analysis were discussed. Second, the secondary post-buckling equilibrium paths of elastic thin rectangular plates and the effects of aspect ratios, boundary conditions and bi-directional compression on the mode jumping behaviors, were discussed in detail. Numerical results show that, the WGM possesses good convergence for solving buckling loads on rectangular plates, and the obtained equilibrium paths are in good agreement with those from the stability experiments, the 2-step perturbation method and the nonlinear finite element method. Given the feasibility of combination with different bifurcation computation methods, the WGM makes an efficient spatial discretization method for complex nonlinear stability problems of typical plates and shells. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Numerical Solution for the Helmholtz Equation with Mixed Boundary Condition

We consider the numerical solution for the Helmholtz equation in R~2 with mixed boundary conditions.The solvability of this mixed boundary value problem is estab- lished by the boundary integral equation method.Based on the Green formula,we express the solution in terms of the boundary data.The key to the numerical real- ization of this method is the computation of weakly singular integrals.Numerical performances show the validity and feasibility of our method.The numerical schemes proposed in this paper have been applied in the realization of probe method for inverse scattering problems.     

On the coupled NBEM and FEM for a class of nonlinear exterior Dirichlet problem in R~(2*)

In this paper, based on the natural boundary reduction advanced by Feng and Yu, we couple the finite element approach with the natural boundary element method to study the weak solvability and Galerkin approximation of a class of nonlinear exterior boundary value problems. The analysis is mainly based on the variational formulation with constraints. We prove the error estimate of the finite element solution and obtain     

THE UNCONDITIONAL CONVERGENT DIFFERENCE METHODS WITH INTRINSIC PARALLELISM FOR QUASILINEAR PARABOLIC SYSTEMS WITH TWO DIMENSIONS

In the present work we are going to solve the boundary value problem for the quasilinear parabolic systems of partial differential equations with two space dimensions by the finite difference method with intrinsic parallelism.Some fundamental behaviors of general finite difference schemes with intrinsic parallelism for the mentioned problems are studied.By the method of a priori estimation of the discrete solutions of the nonlinear difference systems,and the interpolation formulas of the various norms of the discrete functions and the fixed-point technique in finite dimensional Euclidean space,the existennce of the discrete vector solutions of the nonliear difference system with intrinsic parallelism are proved .Moreover the convergence of the discrete vector solutions of these difference schemes to the unique generalizd solution of the original quasilinear parabolic problem is proved.     

Weak Galerkin finite element method for valuation of American options

We introduce a weak Galerkin finite element method for the valuation of American options governed by the Black-Scholes equation. In order to implement, we need to solve the optimal exercise boundary and then introduce an artificial boundary to make the computational domain bounded. For the optimal exercise boundary, which satisfies a nonlinear Volterra integral equation, it is resolved by a higher-order collocation method based on graded meshes. With the computed optimal exercise boundary, the front-fixing technique is employed to transform the free boundary problem to a one- dimensional parabolic problem in a half infinite area. For the other spatial domain boundary, a perfectly matched layer is used to truncate the unbounded domain and carry out the computation. Finally, the resulting initial-boundary value problems are solved by weak Galerkin finite element method, and numerical examples are provided to illustrate the efficiency of the method.     

AN EMBEDDED BOUNDARY METHOD FOR ELLIPTIC AND PARABOLIC PROBLEMS WITH INTERFACES AND APPLICATION TO MULTI-MATERIAL SYSTEMS WITH PHASE TRANSITIONS

The embedded boundary method for solving elliptic and parabolic problems in geometrically complex domains using Cartesian meshes by Johansen and Colella （1998, J. Comput. Phys. 147, 60） has been extended for elliptic and parabolic problems with interior boundaries or interfaces of discontinuities of material properties or solutions. Second order accuracy is achieved in space and time for both stationary and moving interface problems. The method is conservative for elliptic and parabolic problems with fixed interfaces. Based on this method, a front tracking algorithm for the Stefan problem has been developed. The accuracy of the method is measured through comparison with exact solution to a two-dimensional Stefan problem. The algorithm has been used for the study of melting and solidification problems.     

解一类半线性外边值问题的自然边界元与有限元耦合法

In this paper, we combine a finite element approach with the natural boundary element method to stduy the weak solvability and Galerkin approximations of a class of semilinear exterior boundary value problems. Our analysis is mainly based on the variational formulation with constraints. We discuss the error estimate of the finite element solution and obtain the asymptotic rate of convergence O（h^n) Finally, we also give two numerical examples.     

凹角型区域椭圆边值问题的自然边界归化

In this paper, the natural boundary reduction for some elliptic boundary value problems with concave angle domains and their natural boundary methods are investigated. The natural integral equations and the Poisson integral formulae are given. The finite element methods of the natural integral equations are discussed in details. The convergences of the approximate solutions and their error estimates are obtained. Finally, some numerical examples are presented to show that our methods are effective.     

AN EFFECTIVE CONTINUOUS ALGORITHM FOR APPROXIMATE SOLUTIONS OF LARGE SCALE MAX-CUT PROBLEMS

An effective continuous algorithm is proposed to find approximate solutions of NP-hardmax-cut problems.The algorithm relaxes the max-cut problem into a continuous nonlinearprogramming problem by replacing n discrete constraints in the original problem with onesingle continuous constraint.A feasible direction method is designed to solve the resultingnonlinear programming problem.The method employs only the gradient evaluations ofthe objective function,and no any matrix calculations and no line searches are required.This greatly reduces the calculation cost of the method,and is suitable for the solutionof large size max-cut problems.The convergence properties of the proposed method toKKT points of the nonlinear programming are analyzed.If the solution obtained by theproposed method is a global solution of the nonlinear programming problem,the solutionwill provide an upper bound on the max-cut value.Then an approximate solution to themax-cut problem is generated from the solution of the nonlinear programming and providesa lower bound on the max-cut value.Numerical experiments and comparisons on somemax-cut test problems(small and large size)show that the proposed algorithm is efficientto get the exact solutions for all small test problems and well satisfied solutions for mostof the large size test problems with less calculation costs.     

The Riemann-Hilbert problem for mixed complex equations of first order with degenerate rank 0

This article deals with the Riemann-Hilbert boundary value problem for quasilinear mixed(elliptichyperbolic) complex equations of first order with degenerate rank 0. Firstly, we give the representation theorem and prove the uniqueness of solutions for the boundary value problem. Afterwards, by using the method of successive iteration, the existence and estimates of solutions for the boundary value problem are verified. The above problem possesses the important applications to the Tricomi problem of mixed type equations of second order. In this article, the proof of H¨older continuity of a singular double integer is very difficult and interesting as in this Section 4 below.     

AN IMPROVED SEARCH-EXTENTION METHOD FOR SOLVING SEMILINEAR PDES

This article will combine the finite element method, the interpolated coefftcient finite element method, the eigenfunction expansion method, and the search-extension method to obtain the multiple solutions for semilinear elliptic equations. This strategy not only grently reduces the expensive computation, but also is successfully implemented to obtain multiple solutions for a class of semilinear elliptic boundary value problems with non-odd nonlinearity on some convex or nonconvex domains. Numerical solutions illustrated by their graphics for visualization will show the efficiency of the approach.     

The mixed problem for degenerate elliptic equations of second order

The present paper deals with the mixed boundary value problem for elliptic equations with degenerate rank 0. We first give the formulation of the problem and estimates of solutions of the problem, and then prove the existence of solutions of the above problem for elliptic equations by the above estimates and the method of parameter extension. We use the complex method, namely first discuss the corresponding problem for degenerate elliptic complex equations of first order, afterwards discuss the above problem for degenerate elliptic equations of second order.     

奇摄动问题中阶梯状解的缝接法

In view of singularly perturbed problems with complex inner layer phenomenon,including contrast structures(step-step solution and spike-type solution),corner layer behavior and right-hand side discontinuity,we carry out the process with sewing connection.The presented method of sewing connection for singularly perturbed equations is based on the two points singularly perturbed simple boundary problems.By means of sewing orbit smoothness,we get the uniformly valid solution in the whole interval.It is easy to prove the existence of solutions and deal with the high dimensional singularly perturbed problems.     

GLOBALLY DEFINED CLASSICAL SOLUTIONS TO FREE BOUNDARY PROBLEMS WITH CHARACTERISTIC BOUNDARY FOR QUASILINEAR HYPERBOLIC SYSTEMS

In this paper, the authors prove the global existence and uniqueness of classical solutions to some free boundary problems with characteristic boundary for the reducible quasilinear hyperbolic system.     

EXISTENCE OF EXTREME SOLUTION TO FIRST-ORDER IMPULSIVE DIFFERENTIAL EQUATIONS WITH NONLINEAR BOUNDARY CONDITIONS

This paper is concerned with the existence of extreme solutions to three-point boundary value problems with nonlinear boundary conditions for a class of first order impulsive differential equations. We obtain suficient conditions for the existence of extreme solutions by the upper and lower solutions method coupled with a monotone iterative technique.     

RIEMANN-HILBERT PROBLEMS AND SOLITON SOLUTIONS OF NONLOCAL REVERSE-TIME NLS HIERARCHIES

The paper aims at establishing Riemann-Hilbert problems and presenting soliton solutions for nonlocal reverse-time nonlinear Schrodinger(NLS) hierarchies associated with higher-order matrix spectral problems.The Sokhotski-Plemelj formula is used to transform the Riemann-Hilbert problems into Gelfand-Levitan-Marchenko type integral equations.A new formulation of solutions to special Riemann-Hilbert problems with the identity jump matrix,corresponding to the reflectionless inverse scattering transforms,is proposed and applied to construction of soliton solutions to each system in the considered nonlocal reversetime NLS hierarchies.     

RIEMANN-HILBERT PROBLEMS AND SOLITON SOLUTIONS OF NONLOCAL REVERSE-TIME NLS HIERARCHIES

The paper aims at establishing Riemann-Hilbert problems and presenting soliton solutions for nonlocal reverse-time nonlinear Schrodinger(NLS) hierarchies associated with higher-order matrix spectral problems.The Sokhotski-Plemelj formula is used to transform the Riemann-Hilbert problems into Gelfand-Levitan-Marchenko type integral equations.A new formulation of solutions to special Riemann-Hilbert problems with the identity jump matrix,corresponding to the reflectionless inverse scattering transforms,is proposed and applied to construction of soliton solutions to each system in the considered nonlocal reversetime NLS hierarchies.     

The Error Estimates of Semi-discrete Galerkin Methods for some Non-linear Parabolic Equations

1. Introduction Concerning the finite element analysis for linear and nonlinear parabolic equa tions there are a lot of papers, however, only a few of them devoted to the parabolic problems with mixed boundary conditions. In[1],[2] and [3], the finite element methods for some non-linear parabolic problems are systematically considered, but they are mainly restriced to the cases in which the boundary conditions are of Dirichlet     

ON A CLASS OF SINGULAR INTEGRAL EQUATIONS WITH T ANSLATIONS

In this paper,singular integral equations with upper and lower translations arediscussed and the equivalence for solving them in and are shown.Equations with asingle pair of upper and lower translations are studied in detail.Their solutions as wellas the conditions of solvability are obtained.The method used consists of transferringthem to boundary value problems of analytic functions in upper and lower half-planesand then solving the latter.     

A multiscale Galerkin method for the hypersingular integral equation reduced by the harmonic equation

The aim of this paper is to investigate the numerical solution of the hypersingular integral equation reduced by the harmonic equation. First, we transform the hypersingular integral equation into 2π-periodic hypersingular integral equation with the map x=cot(θ/2). Second, we initiate the study of the multiscale Galerkin method for the 2π-periodic hypersingular integral equation. The trigonometric wavelets are used as trial functions. Consequently, the 2j+1 × 2j+1 stiffness matrix Kj can be partitioned j×j block matrices. Furthermore, these block matrices are zeros except main diagonal block matrices. These main diagonal block matrices are symmetrical and circulant matrices, and hence the solution of the associated linear algebraic system can be solved with the fast Fourier transform and the inverse fast Fourier transform instead of the inverse matrix. Finally, we provide several numerical examples to demonstrate our method has good accuracy even though the exact solutions are multi-peak and almost singular.