In this paper,we study the electromagnetic scattering from a two dimen- sional large rectangular open cavity embedded in an infinite ground plane,which is modelled by Helmholtz equations.By introducing nonlocal transp...In this paper,we study the electromagnetic scattering from a two dimen- sional large rectangular open cavity embedded in an infinite ground plane,which is modelled by Helmholtz equations.By introducing nonlocal transparent boundary con- ditions,the problem in the open cavity is reduced to a bounded domain problem.A hypersingular integral operator and a weakly singular integral operator are involved in the TM and TE cases,respectively.A new second-order Toeplitz type approximation and a second-order finite difference scheme are proposed for approximating the hyper- singular integral operator on the aperture and the Helmholtz in the cavity,respectively. The existence and uniqueness of the numerical solution in the TE case are established for arbitrary wavenumbers.A fast algorithm for the second-order approximation is pro- posed for solving the cavity model with layered media.Numerical results show the second-order accuracy and efficiency of the fast algorithm.More important is that the algorithm is easy to implement as a preconditioner for cavity models with more general media.展开更多
We report our recent work on a second-order Krylov subspace and the corresponding second-order Arnoldi procedure for generating its orthonormal basis. The second-order Krylov subspace is spanned by a sequence of vecto...We report our recent work on a second-order Krylov subspace and the corresponding second-order Arnoldi procedure for generating its orthonormal basis. The second-order Krylov subspace is spanned by a sequence of vectors defined via a second-order linear homogeneous recurrence relation with coefficient matrices A and B and an initial vector u. It generalizes the well-known Krylov subspace K n(A;v), which is spanned by a sequence of vectors defined via a first-order linear homogeneous recurrence relation with a single coefficient matrix A and an initial vector v. The applications are shown for the solution of quadratic eigenvalue problems and dimension reduction of second-order dynamical systems. The new approaches preserve essential structures and properties of the quadratic eigenvalue problem and second-order system, and demonstrate superior numerical results over the common approaches based on linearization of these second-order problems.展开更多
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 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.展开更多
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 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 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, 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.
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 θ.展开更多
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 Zn_(3)O_(3) was optimized by the DFT B3LYP method. With the same basis sets, dipole moments, polarizabilities and secondorder hyperpolarizabilities have be...The geometrical structure of semiconductor clusters including Zn_(3)O_(3) 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 discusses a statistical second-order two-scale（SSOTS） analysis and computation for a heat conduction problem with a radiation boundary condition in random porous materials.Firstly,the microscopic configur...This paper discusses a statistical second-order two-scale（SSOTS） analysis and computation for a heat conduction problem with a radiation boundary condition in random porous materials.Firstly,the microscopic configuration for the structure with random distribution is briefly characterized.Secondly,the SSOTS formulae for computing the heat transfer problem are derived successively by means of the construction way for each cell.Then,the statistical prediction algorithm based on the proposed two-scale model is described in detail.Finally,some numerical experiments are proposed,which show that the SSOTS method developed in this paper is effective for predicting the heat transfer performance of porous materials and demonstrating its significant applications in actual engineering computation.展开更多
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.展开更多
The Stirling engine,as a closed-cycle power machine,exhibits excellent emission characteristics and broad energy adaptability.Second-order analysis methods are extensively used during the foundational design and therm...The Stirling engine,as a closed-cycle power machine,exhibits excellent emission characteristics and broad energy adaptability.Second-order analysis methods are extensively used during the foundational design and thermodynamic examination of Stirling engines,owing to their commendable model precision and remarkable efficiency.To scrutinize the effect of Stirling engine design parameters on the cyclical work output and efficiency,this study formulates a series of differential equations for the Stirling cycle by employing second-order analysis methods,subsequently augmenting the predictive accuracy by integrating considerations of loss mechanisms.In addition,an iterative method for the convergence of the average pressure was introduced.The predictive capability of the established model was validated using GPU-3 and RE-1000 experimental data.According to the model,parameters such as the operational fluid,porosity of the regenerator,and diameter of the wire mesh and their influence on the resulting work output and cyclic efficiency of the Stirling engine were analyzed,thereby facilitating a broader understanding of the engine's functional characteristics.These findings suggest that hydrogen,owing to its lower dynamic viscosity coefficient,can provide superior output power.The loss due to flow resistance tends to increase with the rotational speed.Additionally,under conditions of elevated rotational speed,the loss from flow resistance declines in cases of increased porosity,and the enhancement of the porosity to diminish flow resistance losses can boost both the output work and the cyclic efficiency of the engine.As the porosity increased further,the hydraulic diameter and dead volume in the regenerator continued to expand,causing the pressure drop within the engine to become the dominant factor in the gradual reduction of output power.Furthermore,extending the length of the regenerator results in a decrease in the output work,although the thermal cycle efficiency initially increases before eventually decreasing.Based on these insights,this study pursues the optimal designs for Stirling engines.展开更多
基金supported in part by a grant from the Research Grants Council of the Hong Kong Special Administrative Region,China (Project No.CityU 102204).
文摘In this paper,we study the electromagnetic scattering from a two dimen- sional large rectangular open cavity embedded in an infinite ground plane,which is modelled by Helmholtz equations.By introducing nonlocal transparent boundary con- ditions,the problem in the open cavity is reduced to a bounded domain problem.A hypersingular integral operator and a weakly singular integral operator are involved in the TM and TE cases,respectively.A new second-order Toeplitz type approximation and a second-order finite difference scheme are proposed for approximating the hyper- singular integral operator on the aperture and the Helmholtz in the cavity,respectively. The existence and uniqueness of the numerical solution in the TE case are established for arbitrary wavenumbers.A fast algorithm for the second-order approximation is pro- posed for solving the cavity model with layered media.Numerical results show the second-order accuracy and efficiency of the fast algorithm.More important is that the algorithm is easy to implement as a preconditioner for cavity models with more general media.
文摘We report our recent work on a second-order Krylov subspace and the corresponding second-order Arnoldi procedure for generating its orthonormal basis. The second-order Krylov subspace is spanned by a sequence of vectors defined via a second-order linear homogeneous recurrence relation with coefficient matrices A and B and an initial vector u. It generalizes the well-known Krylov subspace K n(A;v), which is spanned by a sequence of vectors defined via a first-order linear homogeneous recurrence relation with a single coefficient matrix A and an initial vector v. The applications are shown for the solution of quadratic eigenvalue problems and dimension reduction of second-order dynamical systems. The new approaches preserve essential structures and properties of the quadratic eigenvalue problem and second-order system, and demonstrate superior numerical results over the common approaches based on linearization of these second-order problems.
基金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.
文摘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(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 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.
基金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.
文摘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.
基金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 θ.
文摘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.
基金supported by the National Natural Science Foundation of China (No. 60176003 and 60376006)
文摘The geometrical structure of semiconductor clusters including Zn_(3)O_(3) 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 China Postdoctoral Science Foundation(Grant Nos.2015M580256 and 2016T90276)
文摘This paper discusses a statistical second-order two-scale（SSOTS） analysis and computation for a heat conduction problem with a radiation boundary condition in random porous materials.Firstly,the microscopic configuration for the structure with random distribution is briefly characterized.Secondly,the SSOTS formulae for computing the heat transfer problem are derived successively by means of the construction way for each cell.Then,the statistical prediction algorithm based on the proposed two-scale model is described in detail.Finally,some numerical experiments are proposed,which show that the SSOTS method developed in this paper is effective for predicting the heat transfer performance of porous materials and demonstrating its significant applications in actual engineering computation.
基金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.
基金supported by Sichuan Science and Technology Program(No.24NSFSC4579)National Natural Science Foundation of China(No.12305193)+2 种基金Sichuan Science and Technology Program(No.23NSFSC6149)National Natural Science Foundation of China(No.12305194)Technology on Reactor System Design Technology Laboratory Stable support Funding(No.2023_JCJQ_LB_003).
文摘The Stirling engine,as a closed-cycle power machine,exhibits excellent emission characteristics and broad energy adaptability.Second-order analysis methods are extensively used during the foundational design and thermodynamic examination of Stirling engines,owing to their commendable model precision and remarkable efficiency.To scrutinize the effect of Stirling engine design parameters on the cyclical work output and efficiency,this study formulates a series of differential equations for the Stirling cycle by employing second-order analysis methods,subsequently augmenting the predictive accuracy by integrating considerations of loss mechanisms.In addition,an iterative method for the convergence of the average pressure was introduced.The predictive capability of the established model was validated using GPU-3 and RE-1000 experimental data.According to the model,parameters such as the operational fluid,porosity of the regenerator,and diameter of the wire mesh and their influence on the resulting work output and cyclic efficiency of the Stirling engine were analyzed,thereby facilitating a broader understanding of the engine's functional characteristics.These findings suggest that hydrogen,owing to its lower dynamic viscosity coefficient,can provide superior output power.The loss due to flow resistance tends to increase with the rotational speed.Additionally,under conditions of elevated rotational speed,the loss from flow resistance declines in cases of increased porosity,and the enhancement of the porosity to diminish flow resistance losses can boost both the output work and the cyclic efficiency of the engine.As the porosity increased further,the hydraulic diameter and dead volume in the regenerator continued to expand,causing the pressure drop within the engine to become the dominant factor in the gradual reduction of output power.Furthermore,extending the length of the regenerator results in a decrease in the output work,although the thermal cycle efficiency initially increases before eventually decreasing.Based on these insights,this study pursues the optimal designs for Stirling engines.
