In this paper a singular perturbation of boundary value problem for elliptic partial differential equations of higher order is considered by using the differential inequalities. The uniformly valid asymptotic expansio...In this paper a singular perturbation of boundary value problem for elliptic partial differential equations of higher order is considered by using the differential inequalities. The uniformly valid asymptotic expansion in entire region is obtained.展开更多
A boundary integral method with radial basis function approximation is proposed for numerically solving an important class of boundary value problems governed by a system of thermoelastostatic equations with variable ...A boundary integral method with radial basis function approximation is proposed for numerically solving an important class of boundary value problems governed by a system of thermoelastostatic equations with variable coe?cients. The equations describe the thermoelastic behaviors of nonhomogeneous anisotropic materials with properties that vary smoothly from point to point in space. No restriction is imposed on the spatial variations of the thermoelastic coe?cients as long as all the requirements of the laws of physics are satis?ed. To check the validity and accuracy of the proposed numerical method, some speci?c test problems with known solutions are solved.展开更多
The paper studies the well-posedness and optimal error estimates of spectral finite element approximations for the boundary value problems of semi-linear elliptic SPDEs driven by white or colored Gaussian noises.The n...The paper studies the well-posedness and optimal error estimates of spectral finite element approximations for the boundary value problems of semi-linear elliptic SPDEs driven by white or colored Gaussian noises.The noise term is approximated through the spectral projection of the covariance operator,which is not required to be commutative with the Laplacian operator.Through the convergence analysis of SPDEs with the noise terms replaced by the projected noises,the well-posedness of the SPDE is established under certain covariance operator-dependent conditions.These SPDEs with projected noises are then numerically approximated with the finite element method.A general error estimate framework is established for the finite element approximations.Based on this framework,optimal error estimates of finite element approximations for elliptic SPDEs driven by power-law noises are obtained.It is shown that with the proposed approach,convergence order of white noise driven SPDEs is improved by half for one-dimensional problems,and by an infinitesimal factor for higher-dimensional problems.展开更多
A complete process of grid generation for complex practical aircraft is described with a twin tail fighter as an example. Euler equations are discretized in the generated multiblock grid by a finite volume method and...A complete process of grid generation for complex practical aircraft is described with a twin tail fighter as an example. Euler equations are discretized in the generated multiblock grid by a finite volume method and solved by a three stage explicit time stepping scheme in each block with some extra treatments of interface at each step. The predicted aerodynamic coefficients and vortical flow field are reasonable.展开更多
The paper deals with growth and approximation of solutions (not necessarily entire) of certain elliptic partial differential equations. These solutions are called Generalized Bi-axially Symmetric Potentials (GBSP'...The paper deals with growth and approximation of solutions (not necessarily entire) of certain elliptic partial differential equations. These solutions are called Generalized Bi-axially Symmetric Potentials (GBSP's). The GBSP's are taken to be regular in a finite hyperball and influence of the growth of their maximum moduli on the rate of decay of their approximation errors in sup norm is studied. The authors obtain the characterizations of the q-type and lower q-type of a GBSP H ∈ HP,0 < R < ∞, in terms of rate of decay of approximation error E.(H,R0), 0 < R0<R <∞.展开更多
We demonstrate a new nonuniform mesh finite difference method to obtain accurate solutions for the elliptic partial differential equations in two dimensions with nonlinear first-order partial derivative terms.The meth...We demonstrate a new nonuniform mesh finite difference method to obtain accurate solutions for the elliptic partial differential equations in two dimensions with nonlinear first-order partial derivative terms.The method will be based on a geometric grid network area and included among the most stable differencing scheme in which the nine-point spatial finite differences are implemented,thus arriving at a compact formulation.In general,a third order of accuracy has been achieved and a fourth-order truncation error in the solution values will follow as a particular case.The efficiency of using geometric mesh ratio parameter has been shown with the help of illustrations.The convergence of the scheme has been established using the matrix analysis,and irreducibility is proved with the help of strongly connected characteristics of the iteration matrix.The difference scheme has been applied to test convection diffusion equation,steady state Burger’s equation,ocean model and a semi-linear elliptic equation.The computational results confirm the theoretical order and accuracy of the method.展开更多
This work proposes a generalized boundary integral method for variable coefficients elliptic partial differential equations(PDEs),including both boundary value and interface problems.The method is kernel-free in the s...This work proposes a generalized boundary integral method for variable coefficients elliptic partial differential equations(PDEs),including both boundary value and interface problems.The method is kernel-free in the sense that there is no need to know analytical expressions for kernels of the boundary and volume integrals in the solution of boundary integral equations.Evaluation of a boundary or volume integral is replaced with interpolation of a Cartesian grid based solution,which satisfies an equivalent discrete interface problem,while the interface problem is solved by a fast solver in the Cartesian grid.The computational work involved with the generalized boundary integral method is essentially linearly proportional to the number of grid nodes in the domain.This paper gives implementation details for a secondorder version of the kernel-free boundary integral method in two space dimensions and presents numerical experiments to demonstrate the efficiency and accuracy of the method for both boundary value and interface problems.The interface problems demonstrated include those with piecewise constant and large-ratio coefficients and the heterogeneous interface problem,where the elliptic PDEs on two sides of the interface are of different types.展开更多
Atmospheric variables(temperature, velocity, etc.) are often decomposed into balanced and unbalanced components that represent low-frequency and high-frequency waves, respectively. Such decompositions can be defined, ...Atmospheric variables(temperature, velocity, etc.) are often decomposed into balanced and unbalanced components that represent low-frequency and high-frequency waves, respectively. Such decompositions can be defined, for instance, in terms of eigenmodes of a linear operator. Traditionally these decompositions ignore phase changes of water since phase changes create a piecewise-linear operator that differs in different phases(cloudy versus non-cloudy). Here we investigate the following question: How can a balanced–unbalanced decomposition be performed in the presence of phase changes? A method is described here motivated by the case of small Froude and Rossby numbers,in which case the asymptotic limit yields precipitating quasi-geostrophic equations with phase changes. Facilitated by its zero-frequency eigenvalue, the balanced component can be found by potential vorticity(PV) inversion, by solving an elliptic partial differential equation(PDE), which includes Heaviside discontinuities due to phase changes. The method is also compared with two simpler methods: one which neglects phase changes, and one which simply treats the raw pressure data as a streamfunction. Tests are shown for both synthetic, idealized data and data from Weather Research and Forecasting(WRF) model simulations. In comparisons, the phase-change method and no-phase-change method produce substantial differences within cloudy regions, of approximately 5K in potential temperature, due to the presence of clouds and phase changes in the data. A theoretical justification is also derived in the form of a elliptic PDE for the differences in the two streamfunctions.展开更多
We study a functional modelling the progressive lens design,which is a combination of Willmore functional and total Gauss curvature.First,we prove the existence for the minimizers of this class of functionals among th...We study a functional modelling the progressive lens design,which is a combination of Willmore functional and total Gauss curvature.First,we prove the existence for the minimizers of this class of functionals among the class of revolution surfaces rotated by the curves y=f(x)about the x-axis.Then,choosing such a minimiser as background surfaces to approximate the functional by a quadratic functional,we prove the existence and uniqueness of the solution to the Euler-Lagrange equation for the quadratic functionals.Our results not only provide a strictly mathematical proof for numerical methods,but also give a more reasonable and more extensive choice for the background surfaces.展开更多
文摘In this paper a singular perturbation of boundary value problem for elliptic partial differential equations of higher order is considered by using the differential inequalities. The uniformly valid asymptotic expansion in entire region is obtained.
文摘A boundary integral method with radial basis function approximation is proposed for numerically solving an important class of boundary value problems governed by a system of thermoelastostatic equations with variable coe?cients. The equations describe the thermoelastic behaviors of nonhomogeneous anisotropic materials with properties that vary smoothly from point to point in space. No restriction is imposed on the spatial variations of the thermoelastic coe?cients as long as all the requirements of the laws of physics are satis?ed. To check the validity and accuracy of the proposed numerical method, some speci?c test problems with known solutions are solved.
基金partially supported by U.S.National Science Foundation,No.DMS1620150U.S.Army ARDEC,No.W911SR-14-2-0001+2 种基金partially supported by National Natural Science Foundation of China,No.91130003,No.11021101,and No.11290142partially supported by Hong Kong RGC General Research Fund,No.16307319the UGC–Research Infrastructure Grant,No.IRS20SC39。
文摘The paper studies the well-posedness and optimal error estimates of spectral finite element approximations for the boundary value problems of semi-linear elliptic SPDEs driven by white or colored Gaussian noises.The noise term is approximated through the spectral projection of the covariance operator,which is not required to be commutative with the Laplacian operator.Through the convergence analysis of SPDEs with the noise terms replaced by the projected noises,the well-posedness of the SPDE is established under certain covariance operator-dependent conditions.These SPDEs with projected noises are then numerically approximated with the finite element method.A general error estimate framework is established for the finite element approximations.Based on this framework,optimal error estimates of finite element approximations for elliptic SPDEs driven by power-law noises are obtained.It is shown that with the proposed approach,convergence order of white noise driven SPDEs is improved by half for one-dimensional problems,and by an infinitesimal factor for higher-dimensional problems.
文摘A complete process of grid generation for complex practical aircraft is described with a twin tail fighter as an example. Euler equations are discretized in the generated multiblock grid by a finite volume method and solved by a three stage explicit time stepping scheme in each block with some extra treatments of interface at each step. The predicted aerodynamic coefficients and vortical flow field are reasonable.
文摘The paper deals with growth and approximation of solutions (not necessarily entire) of certain elliptic partial differential equations. These solutions are called Generalized Bi-axially Symmetric Potentials (GBSP's). The GBSP's are taken to be regular in a finite hyperball and influence of the growth of their maximum moduli on the rate of decay of their approximation errors in sup norm is studied. The authors obtain the characterizations of the q-type and lower q-type of a GBSP H ∈ HP,0 < R < ∞, in terms of rate of decay of approximation error E.(H,R0), 0 < R0<R <∞.
文摘We demonstrate a new nonuniform mesh finite difference method to obtain accurate solutions for the elliptic partial differential equations in two dimensions with nonlinear first-order partial derivative terms.The method will be based on a geometric grid network area and included among the most stable differencing scheme in which the nine-point spatial finite differences are implemented,thus arriving at a compact formulation.In general,a third order of accuracy has been achieved and a fourth-order truncation error in the solution values will follow as a particular case.The efficiency of using geometric mesh ratio parameter has been shown with the help of illustrations.The convergence of the scheme has been established using the matrix analysis,and irreducibility is proved with the help of strongly connected characteristics of the iteration matrix.The difference scheme has been applied to test convection diffusion equation,steady state Burger’s equation,ocean model and a semi-linear elliptic equation.The computational results confirm the theoretical order and accuracy of the method.
基金supported in part by the National Science Foundation of the USA under Grant DMS-0915023is supported by the National Natural Science Foundation of China under Grants DMS-11101278 and DMS-91130012+2 种基金supported by the Young Thousand Talents Program of Chinasupported in part by National Science Committee of Taiwan under Grant 99-2115-M-007-002-MY2supported in part by National Center for Theoretical Sciences of Taiwan,too.
文摘This work proposes a generalized boundary integral method for variable coefficients elliptic partial differential equations(PDEs),including both boundary value and interface problems.The method is kernel-free in the sense that there is no need to know analytical expressions for kernels of the boundary and volume integrals in the solution of boundary integral equations.Evaluation of a boundary or volume integral is replaced with interpolation of a Cartesian grid based solution,which satisfies an equivalent discrete interface problem,while the interface problem is solved by a fast solver in the Cartesian grid.The computational work involved with the generalized boundary integral method is essentially linearly proportional to the number of grid nodes in the domain.This paper gives implementation details for a secondorder version of the kernel-free boundary integral method in two space dimensions and presents numerical experiments to demonstrate the efficiency and accuracy of the method for both boundary value and interface problems.The interface problems demonstrated include those with piecewise constant and large-ratio coefficients and the heterogeneous interface problem,where the elliptic PDEs on two sides of the interface are of different types.
基金supported by the National Science Foundation through grant AGS–1443325 and DMS-1907667the University of Wisconsin–Madison Office of the Vice Chancellor for Research and Graduate Education with funding from the Wisconsin Alumni Research Foundation
文摘Atmospheric variables(temperature, velocity, etc.) are often decomposed into balanced and unbalanced components that represent low-frequency and high-frequency waves, respectively. Such decompositions can be defined, for instance, in terms of eigenmodes of a linear operator. Traditionally these decompositions ignore phase changes of water since phase changes create a piecewise-linear operator that differs in different phases(cloudy versus non-cloudy). Here we investigate the following question: How can a balanced–unbalanced decomposition be performed in the presence of phase changes? A method is described here motivated by the case of small Froude and Rossby numbers,in which case the asymptotic limit yields precipitating quasi-geostrophic equations with phase changes. Facilitated by its zero-frequency eigenvalue, the balanced component can be found by potential vorticity(PV) inversion, by solving an elliptic partial differential equation(PDE), which includes Heaviside discontinuities due to phase changes. The method is also compared with two simpler methods: one which neglects phase changes, and one which simply treats the raw pressure data as a streamfunction. Tests are shown for both synthetic, idealized data and data from Weather Research and Forecasting(WRF) model simulations. In comparisons, the phase-change method and no-phase-change method produce substantial differences within cloudy regions, of approximately 5K in potential temperature, due to the presence of clouds and phase changes in the data. A theoretical justification is also derived in the form of a elliptic PDE for the differences in the two streamfunctions.
基金This work was supported by the National Natural Science Foundation of China(Grant No.11771237).
文摘We study a functional modelling the progressive lens design,which is a combination of Willmore functional and total Gauss curvature.First,we prove the existence for the minimizers of this class of functionals among the class of revolution surfaces rotated by the curves y=f(x)about the x-axis.Then,choosing such a minimiser as background surfaces to approximate the functional by a quadratic functional,we prove the existence and uniqueness of the solution to the Euler-Lagrange equation for the quadratic functionals.Our results not only provide a strictly mathematical proof for numerical methods,but also give a more reasonable and more extensive choice for the background surfaces.