In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a gene...In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.展开更多
The two-dimensional steady flow of an incompressible second-order viscoelastic fluid between two parallel plates was studied in terms of vorticity, the stream function and temperature equations. The governing equation...The two-dimensional steady flow of an incompressible second-order viscoelastic fluid between two parallel plates was studied in terms of vorticity, the stream function and temperature equations. The governing equations were expanded with respect to a snmll parameter to get the zeroth- and first-order approximate equations. By using the differenl2al quadrature method with only a few grid points, the high-accurate numerical results were obtained.展开更多
In this paper,the static output feedback stabilization for large-scale unstable second-order singular systems is investigated.First,the upper bound of all unstable eigenvalues of second-order singular systems is deriv...In this paper,the static output feedback stabilization for large-scale unstable second-order singular systems is investigated.First,the upper bound of all unstable eigenvalues of second-order singular systems is derived.Then,by using the argument principle,a computable stability criterion is proposed to check the stability of secondorder singular systems.Furthermore,by applying model reduction methods to original systems,a static output feedback design algorithm for stabilizing second-order singular systems is presented.A simulation example is provided to illustrate the effectiveness of the design algorithm.展开更多
The correspondence principle is an important mathematical technique to compute the non-ageing linear viscoelastic problem as it allows to take advantage of the computational methods originally developed for the elasti...The correspondence principle is an important mathematical technique to compute the non-ageing linear viscoelastic problem as it allows to take advantage of the computational methods originally developed for the elastic case. However, the correspon- dence principle becomes invalid when the materials exhibit ageing. To deal with this problem, a second-order two-scale (SOTS) computational method in the time domain is presented to predict the ageing linear viscoelastic performance of composite materials with a periodic structure. First, in the time domain, the SOTS formulation for calcu- lating the effective relaxation modulus and displacement approximate solutions of the ageing viscoelastic problem is formally derived. Error estimates of the displacement ap- proximate solutions for SOTS method are then given. Numerical results obtained by the SOTS method are shown and compared with those by the finite element method in a very fine mesh. Both the analytical and numerical results show that the SOTS computational method is feasible and efficient to predict the ageing linear viscoelastic performance of composite materials with a periodic structure.展开更多
In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this...In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations: the first-order (velocity-stress) and the second- order (displacement-stress) separate elastic wave equa- tions, with the first-order (velocity-stress) and the second- order (displacement-stress) full (or mixed) elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.展开更多
In this paper,a stochastic second-order two-scale(SSOTS)method is proposed for predicting the non-deterministic mechanical properties of composites with random interpenetrating phase.Firstly,based on random morphology...In this paper,a stochastic second-order two-scale(SSOTS)method is proposed for predicting the non-deterministic mechanical properties of composites with random interpenetrating phase.Firstly,based on random morphology description functions(RMDF),the randomness of the material properties of the constituents as well as the correlation among these random properties are fully characterized through the topologies of the constituents.Then,by virtue of multiscale asymptotic analysis,the random effective quantities such as stiffness parameters and strength parameters along with their numerical computation formulae are derived by a SSOTS strategy combined with the Monte-Carlo method.Finally,the SSOTS method developed in this paper shows an excellent computational accuracy,and therefore present an important advance towards computationally efficient multiscale modeling frameworks considering microstructure uncertainties.展开更多
In this paper, we construct a uniform second-order difference scheme for a class of boundary value problems of fourth-order ordinary differential equations. Finally, a numerical example is given.
In[20],a semi-implicit spectral deferred correction(SDC)method was proposed,which is efficient for highly nonlinear partial differential equations(PDEs).The semi-implicit SDC method in[20]is based on first-order time ...In[20],a semi-implicit spectral deferred correction(SDC)method was proposed,which is efficient for highly nonlinear partial differential equations(PDEs).The semi-implicit SDC method in[20]is based on first-order time integration methods,which are corrected iteratively,with the order of accuracy increased by one for each additional iteration.In this paper,we will develop a class of semi-implicit SDC methods,which are based on second-order time integration methods and the order of accuracy are increased by two for each additional iteration.For spatial discretization,we employ the local discontinuous Galerkin(LDG)method to arrive at fully-discrete schemes,which are high-order accurate in both space and time.Numerical experiments are presented to demonstrate the accuracy,efficiency and robustness of the proposed semi-implicit SDC methods for solving complex nonlinear PDEs.展开更多
We propose and analyze a single-interval Legendre-Gauss-Radau(LGR)spectral collocation method for nonlinear second-order initial value problems of ordinary differential equations.We design an efficient iterative algor...We propose and analyze a single-interval Legendre-Gauss-Radau(LGR)spectral collocation method for nonlinear second-order initial value problems of ordinary differential equations.We design an efficient iterative algorithm and prove spectral convergence for the single-interval LGR collocation method.For more effective implementation,we propose a multi-interval LGR spectral collocation scheme,which provides us great flexibility with respect to the local time steps and local approximation degrees.Moreover,we combine the multi-interval LGR collocation method in time with the Legendre-Gauss-Lobatto collocation method in space to obtain a space-time spectral collocation approximation for nonlinear second-order evolution equations.Numerical results show that the proposed methods have high accuracy and excellent long-time stability.Numerical comparison between our methods and several commonly used methods are also provided.展开更多
The heat transfer of a magnetohydrodynamics nanofluid inside an annulus considering the second-order slip condition and nanoparticle migration is theoret-ically investigated. A second-order slip condition, which appro...The heat transfer of a magnetohydrodynamics nanofluid inside an annulus considering the second-order slip condition and nanoparticle migration is theoret-ically investigated. A second-order slip condition, which appropriately represents the non-equilibrium region near the interface, is prescribed rather than the no-slip condition and the linear Navier slip condition. To impose different temperature gradients, the outer wall is subjected to q2, the inner wall is subjected to q1, and q1 〉 q2. A modified two-component four-equation non-homogeneous equilibrium model is employed for the nanofiuid, which have been reduced to two-point ordinary boundary value differential equations in the consideration of the thermally and hydrodynamically fully developed flow. The homotopy analysis method (HAM) is employed to solve the equations, and the h-curves are plotted to verify the accuracy and efficiency of the solutions. Moreover, the effects of the physical factors on the flow and heat transfer are discussed in detail, and the semi-analytical relation between NUB and NBT is obtained.展开更多
The normal viscous force of squeeze flow between two arbitrary rigid spheres with an interstitial second-order fluid was studied for modeling wet granular materials using the discrete element method. Based on the Reyn...The normal viscous force of squeeze flow between two arbitrary rigid spheres with an interstitial second-order fluid was studied for modeling wet granular materials using the discrete element method. Based on the Reynolds' lubrication theory, the small parameter method was introduced to approximately analyze velocity field and stress distribution between the two disks. Then a similar procedure was carried out for analyzing the normal interaction between two nearly touching, arbitrary rigid spheres to obtain the pressure distribution and the resulting squeeze force. It has been proved that the solutions can be reduced to the case of a Newtonian fluid when the non-Newtonian terms are neglected.展开更多
In this paper, the extremum of second-order directional derivatives, i.e. the gradient of first-order derivatives is discussed. Given second-order directional derivatives in three nonparallel directions, or given seco...In this paper, the extremum of second-order directional derivatives, i.e. the gradient of first-order derivatives is discussed. Given second-order directional derivatives in three nonparallel directions, or given second-order directional derivatives and mixed directional derivatives in two nonparallel directions, the formulae for the extremum of second-order directional derivatives are derived, and the directions corresponding to maximum and minimum are perpendicular to each other.展开更多
In this article,some high-order local discontinuous Galerkin(LDG)schemes based on some second-order θ approximation formulas in time are presented to solve a two-dimen-sional nonlinear fractional diffusion equation.T...In this article,some high-order local discontinuous Galerkin(LDG)schemes based on some second-order θ approximation formulas in time are presented to solve a two-dimen-sional nonlinear fractional diffusion equation.The unconditional stability of the LDG scheme is proved,and an a priori error estimate with O(h^(k+1)+At^(2))is derived,where k≥0 denotes the index of the basis function.Extensive numerical results with Q^(k)(k=0,1,2,3)elements are provided to confirm our theoretical results,which also show that the second-order convergence rate in time is not impacted by the changed parameter θ.展开更多
Complementary and alternative medicine(CAM)is a set of different diagnostic and therapeutic procedures,as well as the use of natural products for the treatment of patients,derived from previously known traditional met...Complementary and alternative medicine(CAM)is a set of different diagnostic and therapeutic procedures,as well as the use of natural products for the treatment of patients,derived from previously known traditional methods and enriched with modern scientific knowledge.The present article reviews the available data regarding the use of CAM and the legislation behind it in European countries.The use of CAM is recorded in Europe as a whole and varies between 1070%of the population of individual European countries.At least 300,000 registered CAM providers have been identified in the European Union(EU),of which slightly more than half includes non-medical practitioners.The most practiced discipline is acupuncture,followed by homeopathy.CAM regulation and legislation in Europe is not precisely defined and is constantly striving to find a common approach.Since legal frameworks for CAM are not defined,each European country has its own regulations and legislation.In order to define universal legislation for CAM,the EU created the CAMbrella project,a project of the EU designed to find a unique system that would include the treatment of CAM in Europe.According to the data from CAMbrella,from 39 countries in the EU,17 have general CAM legislations.The status of CAM in Europe is characterized by enormous heterogeneity in all aspects,including terminology,methods,prevalence and ultimately,legal status,regulations and legislation.展开更多
Numerical simulations have been performed in time-developing plane mixing layers of the viscoelastic second-order fluids with pseudo-spectral method. Roll-up, pairing and merging of large eddies were examined at high ...Numerical simulations have been performed in time-developing plane mixing layers of the viscoelastic second-order fluids with pseudo-spectral method. Roll-up, pairing and merging of large eddies were examined at high Reynolds numbers and low Deborah numbers. The effect of viscoelastics on the evolution of the large coherent structure was shown by making a comparison between the second-order and Newtonian fluids at the same Reynolds numbers.展开更多
In this paper, we are concerned with the numerical solution of second-order partial differential equations. We analyse the use of the Sine Transform precondilioners for the solution of linear systems arising from the ...In this paper, we are concerned with the numerical solution of second-order partial differential equations. We analyse the use of the Sine Transform precondilioners for the solution of linear systems arising from the discretization of p.d.e. via the preconditioned conjugate gradient method. For the second-order partial differential equations with Dirichlel boundary conditions, we prove that the condition number of the preconditioned system is O(1) while the condition number of the original system is O(m 2) Here m is the number of interior gridpoints in each direction. Such condition number produces a linear convergence rale.展开更多
The geometrical structure of semiconductor clusters including Zn3O3 was optimized by the DFT B3LYP method. With the same basis sets, dipole moments, polarizabilities and secondorder hyperpolarizabilities have been cal...The geometrical structure of semiconductor clusters including Zn3O3 was optimized by the DFT B3LYP method. With the same basis sets, dipole moments, polarizabilities and secondorder hyperpolarizabilities have been calculated and compared with the results obtained by TDDFT B3LYP method combined with sum-over-state (SOS) formula. The calculation results indicate that the dipole moments of the ground state depend on the atom radius and electronegative differences between elements and are their balance point as well. The polarizabilities of the clusters accord with the rule of the corresponding energy transformation from ground to excited state. The results predict an increase of second-order hyperpolarizabilities with increasing the distances between atoms in the clusters as well as a decrease of the polarizabilities and second-order hyperpolarizabilities in the same serial of semiconductor clusters with increasing the dipole moments of the ground states. The changes of dipole moments in ground states are inconsistent with transition moments. Spatial structure, charge transfer and other factors play an important role in composing the transition moments.展开更多
We study the hyperbolic–parabolic equations with rapidly oscillating coefficients. The formal second-order two-scale asymptotic expansion solutions are constructed by the multiscale asymptotic analysis. In addition, ...We study the hyperbolic–parabolic equations with rapidly oscillating coefficients. The formal second-order two-scale asymptotic expansion solutions are constructed by the multiscale asymptotic analysis. In addition, we theoretically explain the importance of the second-order two-scale solution by the error analysis in the pointwise sense. The associated explicit convergence rates are also obtained. Then a second-order two-scale numerical method based on the Newmark scheme is presented to solve the equations. Finally, some numerical examples are used to verify the effectiveness and efficiency of the multiscale numerical algorithm we proposed.展开更多
The superconvergence in the finite element method is a phenomenon in which the fi-nite element approximation converges to the exact solution at a rate higher than the optimal order error estimate. Wang proposed and an...The superconvergence in the finite element method is a phenomenon in which the fi-nite element approximation converges to the exact solution at a rate higher than the optimal order error estimate. Wang proposed and analyzed superconvergence of the conforming finite element method by L2-projections. However, since the conforming finite element method (CFEM) requires a strong continuity, it is not easy to construct such finite elements for the complex partial differential equations. Thus, the nonconforming finite element method (NCFEM) is more appealing computationally due to better stability and flexibility properties compared to CFEM. The objective of this paper is to establish a general superconvergence result for the nonconforming finite element approximations for second-order elliptic problems by L2-projection methods by applying the idea presented in Wang. MATLAB codes are published at https://github.com/annaleeharris/Superconvergence-NCFEM for anyone to use and to study. The results of numerical experiments show great promise for the robustness, reliability, flexibility and accuracy of superconvergence in NCFEM by L2- projections.展开更多
This paper presents methods for computing a second-order sensitivity matrix and the Hessian matrix of eigenvalues and eigenvectors of multiple parameter structures. Second-order perturbations of eigenvalues and eigenv...This paper presents methods for computing a second-order sensitivity matrix and the Hessian matrix of eigenvalues and eigenvectors of multiple parameter structures. Second-order perturbations of eigenvalues and eigenvectors are transformed into multiple parameter forms,and the second-order perturbation sensitivity matrices of eigenvalues and eigenvectors are developed.With these formulations,the efficient methods based on the second-order Taylor expansion and second-order perturbation are obtained to estimate changes of eigenvalues and eigenvectors when the design parameters are changed. The presented method avoids direct differential operation,and thus reduces difficulty for computing the second-order sensitivity matrices of eigenpairs.A numerical example is given to demonstrate application and accuracy of the proposed method.展开更多
基金supported by the Swiss National Science Foundation(Grant No.189882)the National Natural Science Foundation of China(Grant No.41961134032)support provided by the New Investigator Award grant from the UK Engineering and Physical Sciences Research Council(Grant No.EP/V012169/1).
文摘In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.
文摘The two-dimensional steady flow of an incompressible second-order viscoelastic fluid between two parallel plates was studied in terms of vorticity, the stream function and temperature equations. The governing equations were expanded with respect to a snmll parameter to get the zeroth- and first-order approximate equations. By using the differenl2al quadrature method with only a few grid points, the high-accurate numerical results were obtained.
基金Project supported by the National Natural Science Foundation of China(Nos.11971303 and 11871330)。
文摘In this paper,the static output feedback stabilization for large-scale unstable second-order singular systems is investigated.First,the upper bound of all unstable eigenvalues of second-order singular systems is derived.Then,by using the argument principle,a computable stability criterion is proposed to check the stability of secondorder singular systems.Furthermore,by applying model reduction methods to original systems,a static output feedback design algorithm for stabilizing second-order singular systems is presented.A simulation example is provided to illustrate the effectiveness of the design algorithm.
基金Project supported by the National Natural Science Foundation of China(No.11471262)
文摘The correspondence principle is an important mathematical technique to compute the non-ageing linear viscoelastic problem as it allows to take advantage of the computational methods originally developed for the elastic case. However, the correspon- dence principle becomes invalid when the materials exhibit ageing. To deal with this problem, a second-order two-scale (SOTS) computational method in the time domain is presented to predict the ageing linear viscoelastic performance of composite materials with a periodic structure. First, in the time domain, the SOTS formulation for calcu- lating the effective relaxation modulus and displacement approximate solutions of the ageing viscoelastic problem is formally derived. Error estimates of the displacement ap- proximate solutions for SOTS method are then given. Numerical results obtained by the SOTS method are shown and compared with those by the finite element method in a very fine mesh. Both the analytical and numerical results show that the SOTS computational method is feasible and efficient to predict the ageing linear viscoelastic performance of composite materials with a periodic structure.
基金partially supported by China National Major Science and Technology Project (Subproject No:2011ZX05024-001-03)
文摘In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations: the first-order (velocity-stress) and the second- order (displacement-stress) separate elastic wave equa- tions, with the first-order (velocity-stress) and the second- order (displacement-stress) full (or mixed) elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.
基金partially supported by China Postdoctoral Science Foundation(2018M643573)National Natural Science Foundation of Shaanxi Province(2019JQ-048)+2 种基金National Natural Science Foundation of China(51739007,61971328,11301392 and 11961009)of ChinaShanghai Peak Discipline Program for Higher Education Institutions(ClassⅠ)–Civil EngineeringFundamental Research Funds for the Central Universities(No.22120180529)。
文摘In this paper,a stochastic second-order two-scale(SSOTS)method is proposed for predicting the non-deterministic mechanical properties of composites with random interpenetrating phase.Firstly,based on random morphology description functions(RMDF),the randomness of the material properties of the constituents as well as the correlation among these random properties are fully characterized through the topologies of the constituents.Then,by virtue of multiscale asymptotic analysis,the random effective quantities such as stiffness parameters and strength parameters along with their numerical computation formulae are derived by a SSOTS strategy combined with the Monte-Carlo method.Finally,the SSOTS method developed in this paper shows an excellent computational accuracy,and therefore present an important advance towards computationally efficient multiscale modeling frameworks considering microstructure uncertainties.
文摘In this paper, we construct a uniform second-order difference scheme for a class of boundary value problems of fourth-order ordinary differential equations. Finally, a numerical example is given.
基金supported by NSFC(Grant No.11601490).Research of Y.Xu is supported by NSFC(Grant No.12071455).
文摘In[20],a semi-implicit spectral deferred correction(SDC)method was proposed,which is efficient for highly nonlinear partial differential equations(PDEs).The semi-implicit SDC method in[20]is based on first-order time integration methods,which are corrected iteratively,with the order of accuracy increased by one for each additional iteration.In this paper,we will develop a class of semi-implicit SDC methods,which are based on second-order time integration methods and the order of accuracy are increased by two for each additional iteration.For spatial discretization,we employ the local discontinuous Galerkin(LDG)method to arrive at fully-discrete schemes,which are high-order accurate in both space and time.Numerical experiments are presented to demonstrate the accuracy,efficiency and robustness of the proposed semi-implicit SDC methods for solving complex nonlinear PDEs.
基金supported in part by the National Natural Science Foundation of China(Grant Nos.12171322,11771298 and 11871043)the Natural Science Foundation of Shanghai(Grant Nos.21ZR1447200,20ZR1441200 and 22ZR1445500)the Science and Technology Innovation Plan of Shanghai(Grant No.20JC1414200).
文摘We propose and analyze a single-interval Legendre-Gauss-Radau(LGR)spectral collocation method for nonlinear second-order initial value problems of ordinary differential equations.We design an efficient iterative algorithm and prove spectral convergence for the single-interval LGR collocation method.For more effective implementation,we propose a multi-interval LGR spectral collocation scheme,which provides us great flexibility with respect to the local time steps and local approximation degrees.Moreover,we combine the multi-interval LGR collocation method in time with the Legendre-Gauss-Lobatto collocation method in space to obtain a space-time spectral collocation approximation for nonlinear second-order evolution equations.Numerical results show that the proposed methods have high accuracy and excellent long-time stability.Numerical comparison between our methods and several commonly used methods are also provided.
基金Project supported by the National Natural Science Foundation of China(Nos.51476191 and51406008)
文摘The heat transfer of a magnetohydrodynamics nanofluid inside an annulus considering the second-order slip condition and nanoparticle migration is theoret-ically investigated. A second-order slip condition, which appropriately represents the non-equilibrium region near the interface, is prescribed rather than the no-slip condition and the linear Navier slip condition. To impose different temperature gradients, the outer wall is subjected to q2, the inner wall is subjected to q1, and q1 〉 q2. A modified two-component four-equation non-homogeneous equilibrium model is employed for the nanofiuid, which have been reduced to two-point ordinary boundary value differential equations in the consideration of the thermally and hydrodynamically fully developed flow. The homotopy analysis method (HAM) is employed to solve the equations, and the h-curves are plotted to verify the accuracy and efficiency of the solutions. Moreover, the effects of the physical factors on the flow and heat transfer are discussed in detail, and the semi-analytical relation between NUB and NBT is obtained.
文摘The normal viscous force of squeeze flow between two arbitrary rigid spheres with an interstitial second-order fluid was studied for modeling wet granular materials using the discrete element method. Based on the Reynolds' lubrication theory, the small parameter method was introduced to approximately analyze velocity field and stress distribution between the two disks. Then a similar procedure was carried out for analyzing the normal interaction between two nearly touching, arbitrary rigid spheres to obtain the pressure distribution and the resulting squeeze force. It has been proved that the solutions can be reduced to the case of a Newtonian fluid when the non-Newtonian terms are neglected.
基金Supported by the National Natural Science Foundation of China (10871029,11071025)the Foundation of CAEP (2010A0202010)the Foundation of National Key Laboratory of Science and Technology on Computational Physics
文摘In this paper, the extremum of second-order directional derivatives, i.e. the gradient of first-order derivatives is discussed. Given second-order directional derivatives in three nonparallel directions, or given second-order directional derivatives and mixed directional derivatives in two nonparallel directions, the formulae for the extremum of second-order directional derivatives are derived, and the directions corresponding to maximum and minimum are perpendicular to each other.
基金This work is supported by the National Natural Science Foundation of China(11661058,11761053)the Natural Science Foundation of Inner Mongolia(2017MS0107)the Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region(NJYT-17-A07).
文摘In this article,some high-order local discontinuous Galerkin(LDG)schemes based on some second-order θ approximation formulas in time are presented to solve a two-dimen-sional nonlinear fractional diffusion equation.The unconditional stability of the LDG scheme is proved,and an a priori error estimate with O(h^(k+1)+At^(2))is derived,where k≥0 denotes the index of the basis function.Extensive numerical results with Q^(k)(k=0,1,2,3)elements are provided to confirm our theoretical results,which also show that the second-order convergence rate in time is not impacted by the changed parameter θ.
文摘Complementary and alternative medicine(CAM)is a set of different diagnostic and therapeutic procedures,as well as the use of natural products for the treatment of patients,derived from previously known traditional methods and enriched with modern scientific knowledge.The present article reviews the available data regarding the use of CAM and the legislation behind it in European countries.The use of CAM is recorded in Europe as a whole and varies between 1070%of the population of individual European countries.At least 300,000 registered CAM providers have been identified in the European Union(EU),of which slightly more than half includes non-medical practitioners.The most practiced discipline is acupuncture,followed by homeopathy.CAM regulation and legislation in Europe is not precisely defined and is constantly striving to find a common approach.Since legal frameworks for CAM are not defined,each European country has its own regulations and legislation.In order to define universal legislation for CAM,the EU created the CAMbrella project,a project of the EU designed to find a unique system that would include the treatment of CAM in Europe.According to the data from CAMbrella,from 39 countries in the EU,17 have general CAM legislations.The status of CAM in Europe is characterized by enormous heterogeneity in all aspects,including terminology,methods,prevalence and ultimately,legal status,regulations and legislation.
文摘Numerical simulations have been performed in time-developing plane mixing layers of the viscoelastic second-order fluids with pseudo-spectral method. Roll-up, pairing and merging of large eddies were examined at high Reynolds numbers and low Deborah numbers. The effect of viscoelastics on the evolution of the large coherent structure was shown by making a comparison between the second-order and Newtonian fluids at the same Reynolds numbers.
文摘In this paper, we are concerned with the numerical solution of second-order partial differential equations. We analyse the use of the Sine Transform precondilioners for the solution of linear systems arising from the discretization of p.d.e. via the preconditioned conjugate gradient method. For the second-order partial differential equations with Dirichlel boundary conditions, we prove that the condition number of the preconditioned system is O(1) while the condition number of the original system is O(m 2) Here m is the number of interior gridpoints in each direction. Such condition number produces a linear convergence rale.
基金Project supported by the National Natural Science Foundation of China (No. 60176003 and 60376006)
文摘The geometrical structure of semiconductor clusters including Zn3O3 was optimized by the DFT B3LYP method. With the same basis sets, dipole moments, polarizabilities and secondorder hyperpolarizabilities have been calculated and compared with the results obtained by TDDFT B3LYP method combined with sum-over-state (SOS) formula. The calculation results indicate that the dipole moments of the ground state depend on the atom radius and electronegative differences between elements and are their balance point as well. The polarizabilities of the clusters accord with the rule of the corresponding energy transformation from ground to excited state. The results predict an increase of second-order hyperpolarizabilities with increasing the distances between atoms in the clusters as well as a decrease of the polarizabilities and second-order hyperpolarizabilities in the same serial of semiconductor clusters with increasing the dipole moments of the ground states. The changes of dipole moments in ground states are inconsistent with transition moments. Spatial structure, charge transfer and other factors play an important role in composing the transition moments.
基金Project supported by the National Natural Science Foundation of China(Grant No.11471262)the National Basic Research Program of China(Grant No.2012CB025904)the State Key Laboratory of Science and Engineering Computing and the Center for High Performance Computing of Northwestern Polytechnical University,China
文摘We study the hyperbolic–parabolic equations with rapidly oscillating coefficients. The formal second-order two-scale asymptotic expansion solutions are constructed by the multiscale asymptotic analysis. In addition, we theoretically explain the importance of the second-order two-scale solution by the error analysis in the pointwise sense. The associated explicit convergence rates are also obtained. Then a second-order two-scale numerical method based on the Newmark scheme is presented to solve the equations. Finally, some numerical examples are used to verify the effectiveness and efficiency of the multiscale numerical algorithm we proposed.
文摘The superconvergence in the finite element method is a phenomenon in which the fi-nite element approximation converges to the exact solution at a rate higher than the optimal order error estimate. Wang proposed and analyzed superconvergence of the conforming finite element method by L2-projections. However, since the conforming finite element method (CFEM) requires a strong continuity, it is not easy to construct such finite elements for the complex partial differential equations. Thus, the nonconforming finite element method (NCFEM) is more appealing computationally due to better stability and flexibility properties compared to CFEM. The objective of this paper is to establish a general superconvergence result for the nonconforming finite element approximations for second-order elliptic problems by L2-projection methods by applying the idea presented in Wang. MATLAB codes are published at https://github.com/annaleeharris/Superconvergence-NCFEM for anyone to use and to study. The results of numerical experiments show great promise for the robustness, reliability, flexibility and accuracy of superconvergence in NCFEM by L2- projections.
基金Project supported by the 985-Engineering Innovation of Graduate Students of Jilin Universitythe Science and Technology Development Foundation of Jilin Province(20070541)
文摘This paper presents methods for computing a second-order sensitivity matrix and the Hessian matrix of eigenvalues and eigenvectors of multiple parameter structures. Second-order perturbations of eigenvalues and eigenvectors are transformed into multiple parameter forms,and the second-order perturbation sensitivity matrices of eigenvalues and eigenvectors are developed.With these formulations,the efficient methods based on the second-order Taylor expansion and second-order perturbation are obtained to estimate changes of eigenvalues and eigenvectors when the design parameters are changed. The presented method avoids direct differential operation,and thus reduces difficulty for computing the second-order sensitivity matrices of eigenpairs.A numerical example is given to demonstrate application and accuracy of the proposed method.