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 Hamilt...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.展开更多
This paper is concerned with the following fourth-order three-point boundary value problem , where , we discuss the existence of positive solutions to the above problem by applying to the fixed point theory in cones a...This paper is concerned with the following fourth-order three-point boundary value problem , where , we discuss the existence of positive solutions to the above problem by applying to the fixed point theory in cones and iterative technique.展开更多
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-...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.展开更多
Additive Runge-Kutta methods designed for preserving highly accurate solutions in mixed-precision computation were previously proposed and analyzed.These specially designed methods use reduced precision for the implic...Additive Runge-Kutta methods designed for preserving highly accurate solutions in mixed-precision computation were previously proposed and analyzed.These specially designed methods use reduced precision for the implicit computations and full precision for the explicit computations.In this work,we analyze the stability properties of these methods and their sensitivity to the low-precision rounding errors,and demonstrate their performance in terms of accuracy and efficiency.We develop codes in FORTRAN and Julia to solve nonlinear systems of ODEs and PDEs using the mixed-precision additive Runge-Kutta(MP-ARK)methods.The convergence,accuracy,and runtime of these methods are explored.We show that for a given level of accuracy,suitably chosen MP-ARK methods may provide significant reductions in runtime.展开更多
This work presents a comprehensive fourth-order predictive modeling (PM) methodology that uses the MaxEnt principle to incorporate fourth-order moments (means, covariances, skewness, kurtosis) of model parameters, com...This work presents a comprehensive fourth-order predictive modeling (PM) methodology that uses the MaxEnt principle to incorporate fourth-order moments (means, covariances, skewness, kurtosis) of model parameters, computed and measured model responses, as well as fourth (and higher) order sensitivities of computed model responses to model parameters. This new methodology is designated by the acronym 4<sup>th</sup>-BERRU-PM, which stands for “fourth-order best-estimate results with reduced uncertainties.” The results predicted by the 4<sup>th</sup>-BERRU-PM incorporates, as particular cases, the results previously predicted by the second-order predictive modeling methodology 2<sup>nd</sup>-BERRU-PM, and vastly generalizes the results produced by extant data assimilation and data adjustment procedures.展开更多
This work (in two parts) will present a novel predictive modeling methodology aimed at obtaining “best-estimate results with reduced uncertainties” for the first four moments (mean values, covariance, skewness and k...This work (in two parts) will present a novel predictive modeling methodology aimed at obtaining “best-estimate results with reduced uncertainties” for the first four moments (mean values, covariance, skewness and kurtosis) of the optimally predicted distribution of model results and calibrated model parameters, by combining fourth-order experimental and computational information, including fourth (and higher) order sensitivities of computed model responses to model parameters. Underlying the construction of this fourth-order predictive modeling methodology is the “maximum entropy principle” which is initially used to obtain a novel closed-form expression of the (moments-constrained) fourth-order Maximum Entropy (MaxEnt) probability distribution constructed from the first four moments (means, covariances, skewness, kurtosis), which are assumed to be known, of an otherwise unknown distribution of a high-dimensional multivariate uncertain quantity of interest. This fourth-order MaxEnt distribution provides optimal compatibility of the available information while simultaneously ensuring minimal spurious information content, yielding an estimate of a probability density with the highest uncertainty among all densities satisfying the known moment constraints. Since this novel generic fourth-order MaxEnt distribution is of interest in its own right for applications in addition to predictive modeling, its construction is presented separately, in this first part of a two-part work. The fourth-order predictive modeling methodology that will be constructed by particularizing this generic fourth-order MaxEnt distribution will be presented in the accompanying work (Part-2).展开更多
This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either...This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.展开更多
Several existence theorems were established for a nonlinear fourth-order two-point boundary value problem with second derivative by using Leray-Schauder fixed point theorem, equivalent norm and technique on system of ...Several existence theorems were established for a nonlinear fourth-order two-point boundary value problem with second derivative by using Leray-Schauder fixed point theorem, equivalent norm and technique on system of integral equations. The main conditions of our results are local. In other words, the existence of the solution can be determined by considering the height of the nonlinear term on a bounded set. This class of problems usually describes the equilibrium state of an elastic beam which is simply supported at both ends.展开更多
In this work, we present a computational method for solving eigenvalue problems of fourth-order ordinary differential equations which based on the use of Chebychev method. The efficiency of the method is demonstrated ...In this work, we present a computational method for solving eigenvalue problems of fourth-order ordinary differential equations which based on the use of Chebychev method. The efficiency of the method is demonstrated by three numerical examples. Comparison results with others will be presented.展开更多
Under suitable conditions on h(x) and f(u), the authors show that the following boundary value problem has at least one positive solution. Moreover, the authors also establish several existence theorems of multiple po...Under suitable conditions on h(x) and f(u), the authors show that the following boundary value problem has at least one positive solution. Moreover, the authors also establish several existence theorems of multiple positive solutions.展开更多
This paper deals with superlinear fourth-order elliptic problem under Navier boundary condition. By using the mountain pass theorem and suitable truncation, a multiplicity result is established for all λ〉 0 and some...This paper deals with superlinear fourth-order elliptic problem under Navier boundary condition. By using the mountain pass theorem and suitable truncation, a multiplicity result is established for all λ〉 0 and some previous result is extended.展开更多
To consider the effects of the interactions and interplay among microstructures, gradient-dependent models of second- and fourth-order are included in the widely used phenomenological Johnson-Cook model where the effe...To consider the effects of the interactions and interplay among microstructures, gradient-dependent models of second- and fourth-order are included in the widely used phenomenological Johnson-Cook model where the effects of strain-hardening, strain rate sensitivity, and thermal-softening are successfully described. The various parameters for 1006 steel, 4340 steel and S-7 tool steel are assigned. The distributions and evolutions of the local plastic shear strain and deformation in adiabatic shear band (ASB) are predicted. The calculated results of the second- and fourth- order gradient plasticity models are compared. S-7 tool steel possesses the steepest profile of local plastic shear strain in ASB, whereas 1006 steel has the least profile. The peak local plastic shear strain in ASB for S-7 tool steel is slightly higher than that for 4340 steel and is higher than that for 1006 steel. The extent of the nonlinear distribution of the local plastic shear deformation in ASB is more apparent for the S-7 tool steel, whereas it is the least apparent for 1006 steel. In fourth-order gradient plasticity model, the profile of the local plastic shear strain in the middle of ASB has a pronounced plateau whose width decreases with increasing average plastic shear strain, leading to a shrink of the portion of linear distribution of the profile of the local plastic shear deformation. When compared with the sec- ond-order gradient plasticity model, the fourth-order gradient plasticity model shows a lower peak local plastic shear strain in ASB and a higher magnitude of plastic shear deformation at the top or base of ASB, which is due to wider ASB. The present numerical results of the second- and fourth-order gradient plasticity models are consistent with the previous numerical and experimental results at least qualitatively.展开更多
A fourth-order variational inequality of the second kind arising in a plate frictional bending problem is considered. By using regularization method, the original problem can be formulated as a differentiable variatio...A fourth-order variational inequality of the second kind arising in a plate frictional bending problem is considered. By using regularization method, the original problem can be formulated as a differentiable variational equation, and the corresponding discrete FEM variational equation is presented afterwards. Abstract error estimates and error estimates of the approximation are derived in terms of energy norm and L^2-norm.展开更多
In this paper, existence and uniqueness of the generalized global solution and the classical global solution to the initial value problem for a class of fourth-order nonlinear wave equations are studied in the fractio...In this paper, existence and uniqueness of the generalized global solution and the classical global solution to the initial value problem for a class of fourth-order nonlinear wave equations are studied in the fractional order Sobolev space using the contraction mapping principle and the extension theorem. The sufficient conditions for the blow up of the solution to the initial value problem are given.展开更多
In this paper we consider the singularly perturbed boundary value problem for the fourth-order elliptic differential equation, establish the energy estimates of the solutionand its derivatives and construct the formal...In this paper we consider the singularly perturbed boundary value problem for the fourth-order elliptic differential equation, establish the energy estimates of the solutionand its derivatives and construct the formal asymptotic solution by Lyuternik- Vishik 's method. Finally, by means of the energy estimates we obtain the bound of the remainder of the asymptotic solution.展开更多
In this work, by extending the method of Hockney into three dimensions, the Poisson’s equation in cylindrical coordinates system with the Dirichlet’s boundary conditions in a portion of a cylinder for is solved dire...In this work, by extending the method of Hockney into three dimensions, the Poisson’s equation in cylindrical coordinates system with the Dirichlet’s boundary conditions in a portion of a cylinder for is solved directly. The Poisson equation is approximated by fourth-order finite differences and the resulting large algebraic system of linear equations is treated systematically in order to get a block tri-diagonal system. The accuracy of this method is tested for some Poisson’s equations with known analytical solutions and the numerical results obtained show that the method produces accurate results.展开更多
The finite data estimates of the complex fourth-order moments of a signal consisting of random harmonics are analyzed. Conditions for the fourth-order stationarity and ergodicity are obtained. Explicit formulas for th...The finite data estimates of the complex fourth-order moments of a signal consisting of random harmonics are analyzed. Conditions for the fourth-order stationarity and ergodicity are obtained. Explicit formulas for the estimation error and its variance, as well as their limiting large sample values are derived. Finally, a special case relevant to cubic phase coupling is considered, and these results are stated for this case, the variance is shown to comprise an ergodic and a nonergodic part.展开更多
基金Project supported by the National Natural Science Foundation of China(Nos.12072297 and12202370)the Natural Science Foundation of Sichuan Province of China(No.24NSFSC4777)。
文摘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.
文摘This paper is concerned with the following fourth-order three-point boundary value problem , where , we discuss the existence of positive solutions to the above problem by applying to the fixed point theory in cones and iterative technique.
基金supported by the NationalNatural Science Foundation of China(12001236)the Natural Science Foundation of Guangdong Province(2020A1515110494)。
文摘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.
基金supported by ONR UMass Dartmouth Marine and UnderSea Technology(MUST)grant N00014-20-1-2849 under the project S31320000049160by DOE grant DE-SC0023164 sub-award RC114586-UMD+2 种基金by AFOSR grants FA9550-18-1-0383 and FA9550-23-1-0037supported by Michigan State University,by AFOSR grants FA9550-19-1-0281 and FA9550-18-1-0383by DOE grant DE-SC0023164.
文摘Additive Runge-Kutta methods designed for preserving highly accurate solutions in mixed-precision computation were previously proposed and analyzed.These specially designed methods use reduced precision for the implicit computations and full precision for the explicit computations.In this work,we analyze the stability properties of these methods and their sensitivity to the low-precision rounding errors,and demonstrate their performance in terms of accuracy and efficiency.We develop codes in FORTRAN and Julia to solve nonlinear systems of ODEs and PDEs using the mixed-precision additive Runge-Kutta(MP-ARK)methods.The convergence,accuracy,and runtime of these methods are explored.We show that for a given level of accuracy,suitably chosen MP-ARK methods may provide significant reductions in runtime.
文摘This work presents a comprehensive fourth-order predictive modeling (PM) methodology that uses the MaxEnt principle to incorporate fourth-order moments (means, covariances, skewness, kurtosis) of model parameters, computed and measured model responses, as well as fourth (and higher) order sensitivities of computed model responses to model parameters. This new methodology is designated by the acronym 4<sup>th</sup>-BERRU-PM, which stands for “fourth-order best-estimate results with reduced uncertainties.” The results predicted by the 4<sup>th</sup>-BERRU-PM incorporates, as particular cases, the results previously predicted by the second-order predictive modeling methodology 2<sup>nd</sup>-BERRU-PM, and vastly generalizes the results produced by extant data assimilation and data adjustment procedures.
文摘This work (in two parts) will present a novel predictive modeling methodology aimed at obtaining “best-estimate results with reduced uncertainties” for the first four moments (mean values, covariance, skewness and kurtosis) of the optimally predicted distribution of model results and calibrated model parameters, by combining fourth-order experimental and computational information, including fourth (and higher) order sensitivities of computed model responses to model parameters. Underlying the construction of this fourth-order predictive modeling methodology is the “maximum entropy principle” which is initially used to obtain a novel closed-form expression of the (moments-constrained) fourth-order Maximum Entropy (MaxEnt) probability distribution constructed from the first four moments (means, covariances, skewness, kurtosis), which are assumed to be known, of an otherwise unknown distribution of a high-dimensional multivariate uncertain quantity of interest. This fourth-order MaxEnt distribution provides optimal compatibility of the available information while simultaneously ensuring minimal spurious information content, yielding an estimate of a probability density with the highest uncertainty among all densities satisfying the known moment constraints. Since this novel generic fourth-order MaxEnt distribution is of interest in its own right for applications in addition to predictive modeling, its construction is presented separately, in this first part of a two-part work. The fourth-order predictive modeling methodology that will be constructed by particularizing this generic fourth-order MaxEnt distribution will be presented in the accompanying work (Part-2).
基金supported by the NSF under Grant DMS-2208391sponsored by the NSF under Grant DMS-1753581.
文摘This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.
文摘Several existence theorems were established for a nonlinear fourth-order two-point boundary value problem with second derivative by using Leray-Schauder fixed point theorem, equivalent norm and technique on system of integral equations. The main conditions of our results are local. In other words, the existence of the solution can be determined by considering the height of the nonlinear term on a bounded set. This class of problems usually describes the equilibrium state of an elastic beam which is simply supported at both ends.
文摘In this work, we present a computational method for solving eigenvalue problems of fourth-order ordinary differential equations which based on the use of Chebychev method. The efficiency of the method is demonstrated by three numerical examples. Comparison results with others will be presented.
文摘Under suitable conditions on h(x) and f(u), the authors show that the following boundary value problem has at least one positive solution. Moreover, the authors also establish several existence theorems of multiple positive solutions.
基金The 985 Program of Jilin Universitythe Science Research Foundation for Excellent Young Teachers of College of Mathematics at Jilin University
文摘This paper deals with superlinear fourth-order elliptic problem under Navier boundary condition. By using the mountain pass theorem and suitable truncation, a multiplicity result is established for all λ〉 0 and some previous result is extended.
基金Item Sponsored by Educational Department of Liaoning Province of China (2004F052)
文摘To consider the effects of the interactions and interplay among microstructures, gradient-dependent models of second- and fourth-order are included in the widely used phenomenological Johnson-Cook model where the effects of strain-hardening, strain rate sensitivity, and thermal-softening are successfully described. The various parameters for 1006 steel, 4340 steel and S-7 tool steel are assigned. The distributions and evolutions of the local plastic shear strain and deformation in adiabatic shear band (ASB) are predicted. The calculated results of the second- and fourth- order gradient plasticity models are compared. S-7 tool steel possesses the steepest profile of local plastic shear strain in ASB, whereas 1006 steel has the least profile. The peak local plastic shear strain in ASB for S-7 tool steel is slightly higher than that for 4340 steel and is higher than that for 1006 steel. The extent of the nonlinear distribution of the local plastic shear deformation in ASB is more apparent for the S-7 tool steel, whereas it is the least apparent for 1006 steel. In fourth-order gradient plasticity model, the profile of the local plastic shear strain in the middle of ASB has a pronounced plateau whose width decreases with increasing average plastic shear strain, leading to a shrink of the portion of linear distribution of the profile of the local plastic shear deformation. When compared with the sec- ond-order gradient plasticity model, the fourth-order gradient plasticity model shows a lower peak local plastic shear strain in ASB and a higher magnitude of plastic shear deformation at the top or base of ASB, which is due to wider ASB. The present numerical results of the second- and fourth-order gradient plasticity models are consistent with the previous numerical and experimental results at least qualitatively.
基金Supported by the National Natural Science Foundation of China(10201026,10672111)
文摘A fourth-order variational inequality of the second kind arising in a plate frictional bending problem is considered. By using regularization method, the original problem can be formulated as a differentiable variational equation, and the corresponding discrete FEM variational equation is presented afterwards. Abstract error estimates and error estimates of the approximation are derived in terms of energy norm and L^2-norm.
基金supported by the National Natural Science Foundation of China (No. 10671182)
文摘In this paper, existence and uniqueness of the generalized global solution and the classical global solution to the initial value problem for a class of fourth-order nonlinear wave equations are studied in the fractional order Sobolev space using the contraction mapping principle and the extension theorem. The sufficient conditions for the blow up of the solution to the initial value problem are given.
文摘In this paper we consider the singularly perturbed boundary value problem for the fourth-order elliptic differential equation, establish the energy estimates of the solutionand its derivatives and construct the formal asymptotic solution by Lyuternik- Vishik 's method. Finally, by means of the energy estimates we obtain the bound of the remainder of the asymptotic solution.
文摘In this work, by extending the method of Hockney into three dimensions, the Poisson’s equation in cylindrical coordinates system with the Dirichlet’s boundary conditions in a portion of a cylinder for is solved directly. The Poisson equation is approximated by fourth-order finite differences and the resulting large algebraic system of linear equations is treated systematically in order to get a block tri-diagonal system. The accuracy of this method is tested for some Poisson’s equations with known analytical solutions and the numerical results obtained show that the method produces accurate results.
文摘The finite data estimates of the complex fourth-order moments of a signal consisting of random harmonics are analyzed. Conditions for the fourth-order stationarity and ergodicity are obtained. Explicit formulas for the estimation error and its variance, as well as their limiting large sample values are derived. Finally, a special case relevant to cubic phase coupling is considered, and these results are stated for this case, the variance is shown to comprise an ergodic and a nonergodic part.