Three modified sine-Hilbert(sH)-type equations, i.e., the modified sH equation, the modified damped sH equation, and the modified nonlinear dissipative system, are proposed, and their bilinear forms are provided.Based...Three modified sine-Hilbert(sH)-type equations, i.e., the modified sH equation, the modified damped sH equation, and the modified nonlinear dissipative system, are proposed, and their bilinear forms are provided.Based on these bilinear equations, some exact solutions to the three modified equations are derived.展开更多
The purpose of this article is to develop and analyze least-squares approximations for the incompressible magnetohydrodynamic equations. The major advantage of the least-squares finite element method is that it is not...The purpose of this article is to develop and analyze least-squares approximations for the incompressible magnetohydrodynamic equations. The major advantage of the least-squares finite element method is that it is not subjected to the so-called Ladyzhenskaya-Babuska-Brezzi （LBB） condition. The authors employ least-squares functionals which involve a discrete inner product which is related to the inner product in H^-1（Ω）.展开更多
A novel statistical second-order reduced multiscale(SSRM)approach is established for nonlinear composite materials with random distribution of grains.For these composites considered in this work,the complex microstruc...A novel statistical second-order reduced multiscale(SSRM)approach is established for nonlinear composite materials with random distribution of grains.For these composites considered in this work,the complex microstructure of grains,including their shape,orientation,size,spatial distribution,volume fraction and so on,results in changing of the macroscopic mechanical properties.The first-and second-order unit cell functions based on two-scale asymptotic expressions are constructed at first.Then,the expected homogenized parameters are defined,and the nonlinear homogenization equation on global structure is established,successively.Further,an effective reduced model format for analyzing second-order nonlinear unit cell problem with less computation cost is introduced in detail.Finally,some numerical examples for the materials with varying distribution models are evaluated and compared with the data by theoretical models and experimental results.These examples illustrate that the proposed SSRM approaches are effective for predicting the macroscopic properties of the random composite materials and supply a potential application in actual engineering computation.展开更多
This paper presents a discrete vaxiational principle and a method to build first-integrals for finite dimensional Lagrange-Maxwell mechanico-electrical systems with nonconservative forces and a dissipation function. T...This paper presents a discrete vaxiational principle and a method to build first-integrals for finite dimensional Lagrange-Maxwell mechanico-electrical systems with nonconservative forces and a dissipation function. The discrete variational principle and the corresponding Euler-Lagrange equations are derived from a discrete action associated to these systems. The first-integrals are obtained by introducing the infinitesimal transformation with respect to the generalized coordinates and electric quantities of the systems. This work also extends discrete Noether symmetries to mechanico-electrical dynamical systems. A practical example is presented to illustrate the results.展开更多
We develop a new full waveform inversion (FWI) method for slowness with the crosshole data based on the acoustic wave equation in the time domain. The method combines the total variation (TV) regularization with the c...We develop a new full waveform inversion (FWI) method for slowness with the crosshole data based on the acoustic wave equation in the time domain. The method combines the total variation (TV) regularization with the constrained optimization together which can inverse the slowness effectively. One advantage of slowness inversion is that there is no further approximation in the gradient derivation. Moreover, a new algorithm named the skip method for solving the constrained optimization problem is proposed. The TV regularization has good ability to inverse slowness at its discontinuities while the constrained optimization can keep the inversion converging in the right direction. Numerical computations both for noise free data and noisy data show the robustness and effectiveness of our method and good inversion results are yielded.展开更多
This paper presents a method to find Noether-type conserved quantities and Lie point symmetries for discrete mechanico-electrical dynamical systems,which leave invuriant the set of solutions of the corresponding diffe...This paper presents a method to find Noether-type conserved quantities and Lie point symmetries for discrete mechanico-electrical dynamical systems,which leave invuriant the set of solutions of the corresponding difference scheme. This approach makes it possible to devise techniques for solving the Lagrange Maxwell equations in differences which correspond to mechanico-electrical systems,by adapting existing differential equations.In particular,it obtains a new systematic method to determine both the one-parameter Lie groups and the discrete Noether conserved quantities of Lie point symmetries for mechanico-electrical systems.As an application,it obtains the Lie point symmetries and the conserved quantities for the difference equation of a model that represents a capacitor microphone.展开更多
In this paper,we consider 3 D tomographic reconstruction for axially symmetric objects from a single radiograph formed by cone-beam X-rays.All contemporary density reconstruction methods in high-energy X-ray radiograp...In this paper,we consider 3 D tomographic reconstruction for axially symmetric objects from a single radiograph formed by cone-beam X-rays.All contemporary density reconstruction methods in high-energy X-ray radiography are based on the assumption that the cone beam can be treated as fan beams located at parallel planes perpendicular to the symmetric axis,so that the density of the whole object can be recovered layer by layer.Considering the relationship between different layers,we undertake the cone-beam global reconstruction to solve the ambiguity effect at the material interfaces of the reconstruction results.In view of the anisotropy of classical discrete total variations,a new discretization of total variation which yields sharp edges and has better isotropy is introduced in our reconstruction model.Furthermore,considering that the object density consists of continually changing parts and jumps,a high-order regularization term is introduced.The final hybrid regularization model is solved using the alternating proximal gradient method,which was recently applied in image processing.Density reconstruction results are presented for simulated radiographs,which shows that the proposed method has led to an improvement in terms of the preservation of edge location.展开更多
We propose efficient numerical methods for nonseparable non-canonical Hamiltonian systems which are explicit,K-symplectic in the extended phase space with long time energy conservation properties. They are based on ex...We propose efficient numerical methods for nonseparable non-canonical Hamiltonian systems which are explicit,K-symplectic in the extended phase space with long time energy conservation properties. They are based on extending the original phase space to several copies of the phase space and imposing a mechanical restraint on the copies of the phase space. Explicit K-symplectic methods are constructed for two non-canonical Hamiltonian systems. Numerical tests show that the proposed methods exhibit good numerical performance in preserving the phase orbit and the energy of the system over long time, whereas higher order Runge–Kutta methods do not preserve these properties. Numerical tests also show that the K-symplectic methods exhibit better efficiency than that of the same order implicit symplectic, explicit and implicit symplectic methods for the original nonseparable non-canonical systems. On the other hand, the fourth order K-symplectic method is more efficient than the fourth order Yoshida’s method, the optimized partitioned Runge–Kutta and Runge–Kutta–Nystr ¨om explicit K-symplectic methods for the extended phase space Hamiltonians, but less efficient than the the optimized partitioned Runge–Kutta and Runge–Kutta–Nystr ¨om extended phase space symplectic-like methods with the midpoint permutation.展开更多
In this paper, a kind of second-order two-scale （SOTS） computation is developed for conductive-radiative heat trans- fer problem in periodic porous materials. First of all, by the asymptotic expansion of the tempera...In this paper, a kind of second-order two-scale （SOTS） computation is developed for conductive-radiative heat trans- fer problem in periodic porous materials. First of all, by the asymptotic expansion of the temperature field, the cell problem, homogenization problem, and second-order correctors are obtained successively. Then, the corresponding finite element al- gorithms are proposed. Finally, some numerical results are presented and compared with theoretical results. The numerical results of the proposed algorithm conform with those of the FE algorithm well, demonstrating the accuracy of the present method and its potential applications in thermal engineering of porous materials.展开更多
In this paper,a statistical second-order twoscale（SSOTS） method is developed to simulate the dynamic thcrmo-mechanical performances of the statistically inhomogeneous materials.For this kind of composite material,th...In this paper,a statistical second-order twoscale（SSOTS） method is developed to simulate the dynamic thcrmo-mechanical performances of the statistically inhomogeneous materials.For this kind of composite material,the random distribution characteristics of particles,including the shape,size,orientation,spatial location,and volume fractions,are all considered.Firstly,the repre.sentation for the microscopic configuration of the statistically inhomogeneous materials is described.Secondly,the SSOTS formulation for the dynamic thermo-mechanical coupled problem is proposed in a constructive way,including the cell problems,effective thermal and mechanical parameters,homogenized problems,and the SSOTS formulas of the temperatures,displacements,heat flux densities and stresses.And then the algorithm procedure corresponding to the SSOTS method is brought forward.The numerical results obtained by using the SSOTS algorithm are compared with those by classical methods.In addition,the thermo-mechanical coupling effect is studied by comparing the results of coupled case with those of uncoupled case.It demonstrates that the coupling effect on the temperatures,heat flux densities,displacements,and stresses is very distinct.The results show that the SSOTS method is valid to predict the dynamic thermo-mechanical coupled performances of statistically inhomogeneous materials.展开更多
This paper focuses on the dynamic thermo-mechanical coupled response of random particulate composite materials. Both the inertia term and coupling term are considered in the dynamic coupled problem. The formulation of...This paper focuses on the dynamic thermo-mechanical coupled response of random particulate composite materials. Both the inertia term and coupling term are considered in the dynamic coupled problem. The formulation of the problem by a statistical second-order two-scale （SSOTS） analysis method and the algorithm procedure based on the finite-element difference method are presented. Numerical results of coupled cases are compared with those of uncoupled cases. It shows that the coupling effects on temperature, thermal flux, displacement, and stresses are very distinct, and the micro- characteristics of particles affect the coupling effect of the random composites. Furthermore, the coupling effect causes a lag in the variations of temperature, thermal flux, displacement, and stresses.展开更多
Two general local Cm triangular interpolation schemes by rational functions from Cm data are proposed for any nonnegative integer m. The schemes can have either 2m+1 order algebraic precision if the required data are ...Two general local Cm triangular interpolation schemes by rational functions from Cm data are proposed for any nonnegative integer m. The schemes can have either 2m+1 order algebraic precision if the required data are given on vertices and edges, or m+E[m/2]+1 or m+1 order algebraic precision if the data are given only at vertices. The orders of the interpolation error are estimated. Examples that show the correctness and effectiveness of the scheme are presented.展开更多
A new variational method is proposed to investigate the dynamics of the thin film in a coating flow where a liquid is delivered through a fixed slot gap onto a moving substrate. A simplified ODE system has also been d...A new variational method is proposed to investigate the dynamics of the thin film in a coating flow where a liquid is delivered through a fixed slot gap onto a moving substrate. A simplified ODE system has also been derived for the evolution of the thin film whose thickness hf is asymptotically constant behind the coating front. We calculate the phase diagram as well as the film profiles and approximate the film thickness theoretically, and agreement with the well-known scaling law as Ca2/3 is found.展开更多
基金supported by the National Natural Science Foundation of China (Grant Nos. 11931017 and 12071447)。
文摘Three modified sine-Hilbert(sH)-type equations, i.e., the modified sH equation, the modified damped sH equation, and the modified nonlinear dissipative system, are proposed, and their bilinear forms are provided.Based on these bilinear equations, some exact solutions to the three modified equations are derived.
基金supported by the National Basic Research Program of China (2005CB321701)NSF of mathematics research special fund of Hebei Province(08M005)
文摘The purpose of this article is to develop and analyze least-squares approximations for the incompressible magnetohydrodynamic equations. The major advantage of the least-squares finite element method is that it is not subjected to the so-called Ladyzhenskaya-Babuska-Brezzi （LBB） condition. The authors employ least-squares functionals which involve a discrete inner product which is related to the inner product in H^-1（Ω）.
基金This study was funded by the National Natural Science Foundation of China(Grant 11701123)Fundamental Research Funds for the Central Universities(Grant HIT.NSRIF.2020017).
文摘A novel statistical second-order reduced multiscale(SSRM)approach is established for nonlinear composite materials with random distribution of grains.For these composites considered in this work,the complex microstructure of grains,including their shape,orientation,size,spatial distribution,volume fraction and so on,results in changing of the macroscopic mechanical properties.The first-and second-order unit cell functions based on two-scale asymptotic expressions are constructed at first.Then,the expected homogenized parameters are defined,and the nonlinear homogenization equation on global structure is established,successively.Further,an effective reduced model format for analyzing second-order nonlinear unit cell problem with less computation cost is introduced in detail.Finally,some numerical examples for the materials with varying distribution models are evaluated and compared with the data by theoretical models and experimental results.These examples illustrate that the proposed SSRM approaches are effective for predicting the macroscopic properties of the random composite materials and supply a potential application in actual engineering computation.
基金Project supported by State Key Laboratory of Scientific and Engineering Computing, Chinese Academy of Sciences and the National Natural Science Foundation of China （Grant Nos 10672143 and 10471145） and the Natural Science Foundation of Henan Province Government, China （Grant Nos 0311011400 and 0511022200）.
文摘This paper presents a discrete vaxiational principle and a method to build first-integrals for finite dimensional Lagrange-Maxwell mechanico-electrical systems with nonconservative forces and a dissipation function. The discrete variational principle and the corresponding Euler-Lagrange equations are derived from a discrete action associated to these systems. The first-integrals are obtained by introducing the infinitesimal transformation with respect to the generalized coordinates and electric quantities of the systems. This work also extends discrete Noether symmetries to mechanico-electrical dynamical systems. A practical example is presented to illustrate the results.
文摘We develop a new full waveform inversion (FWI) method for slowness with the crosshole data based on the acoustic wave equation in the time domain. The method combines the total variation (TV) regularization with the constrained optimization together which can inverse the slowness effectively. One advantage of slowness inversion is that there is no further approximation in the gradient derivation. Moreover, a new algorithm named the skip method for solving the constrained optimization problem is proposed. The TV regularization has good ability to inverse slowness at its discontinuities while the constrained optimization can keep the inversion converging in the right direction. Numerical computations both for noise free data and noisy data show the robustness and effectiveness of our method and good inversion results are yielded.
基金Project supported by the National Natural Science Foundation of China (Grants Nos 10672143 and 60575055)State Key Laboratory of Scientific and Engineering Computing,Chinese Academy of Sciences+1 种基金Tang Yi-Fa acknowledges the support under Sabbatical Program (SAB2006-0070) of the Spanish Ministry of Education and ScienceJimnez S and Vzquez L acknowledge support of the Spanish Ministry of Education and Science (Grant No MTM2005-05573)
文摘This paper presents a method to find Noether-type conserved quantities and Lie point symmetries for discrete mechanico-electrical dynamical systems,which leave invuriant the set of solutions of the corresponding difference scheme. This approach makes it possible to devise techniques for solving the Lagrange Maxwell equations in differences which correspond to mechanico-electrical systems,by adapting existing differential equations.In particular,it obtains a new systematic method to determine both the one-parameter Lie groups and the discrete Noether conserved quantities of Lie point symmetries for mechanico-electrical systems.As an application,it obtains the Lie point symmetries and the conserved quantities for the difference equation of a model that represents a capacitor microphone.
基金supported by National Postdoctoral Program for Innovative Talents(BX201700038)supported by NSFC(11571003)+1 种基金supported by NSFC(11675021)supported by Beijing Natural Science Foundation(Z180002)。
文摘In this paper,we consider 3 D tomographic reconstruction for axially symmetric objects from a single radiograph formed by cone-beam X-rays.All contemporary density reconstruction methods in high-energy X-ray radiography are based on the assumption that the cone beam can be treated as fan beams located at parallel planes perpendicular to the symmetric axis,so that the density of the whole object can be recovered layer by layer.Considering the relationship between different layers,we undertake the cone-beam global reconstruction to solve the ambiguity effect at the material interfaces of the reconstruction results.In view of the anisotropy of classical discrete total variations,a new discretization of total variation which yields sharp edges and has better isotropy is introduced in our reconstruction model.Furthermore,considering that the object density consists of continually changing parts and jumps,a high-order regularization term is introduced.The final hybrid regularization model is solved using the alternating proximal gradient method,which was recently applied in image processing.Density reconstruction results are presented for simulated radiographs,which shows that the proposed method has led to an improvement in terms of the preservation of edge location.
基金supported by the National Natural Science Foundation of China (Grant Nos. 11901564 and 12171466)。
文摘We propose efficient numerical methods for nonseparable non-canonical Hamiltonian systems which are explicit,K-symplectic in the extended phase space with long time energy conservation properties. They are based on extending the original phase space to several copies of the phase space and imposing a mechanical restraint on the copies of the phase space. Explicit K-symplectic methods are constructed for two non-canonical Hamiltonian systems. Numerical tests show that the proposed methods exhibit good numerical performance in preserving the phase orbit and the energy of the system over long time, whereas higher order Runge–Kutta methods do not preserve these properties. Numerical tests also show that the K-symplectic methods exhibit better efficiency than that of the same order implicit symplectic, explicit and implicit symplectic methods for the original nonseparable non-canonical systems. On the other hand, the fourth order K-symplectic method is more efficient than the fourth order Yoshida’s method, the optimized partitioned Runge–Kutta and Runge–Kutta–Nystr ¨om explicit K-symplectic methods for the extended phase space Hamiltonians, but less efficient than the the optimized partitioned Runge–Kutta and Runge–Kutta–Nystr ¨om extended phase space symplectic-like methods with the midpoint permutation.
基金Project supported by the National Basic Research Program of China(Grant No.2010CB832702)the National Natural Science Foundation of China(Grant No.90916027)
文摘In this paper, a kind of second-order two-scale （SOTS） computation is developed for conductive-radiative heat trans- fer problem in periodic porous materials. First of all, by the asymptotic expansion of the temperature field, the cell problem, homogenization problem, and second-order correctors are obtained successively. Then, the corresponding finite element al- gorithms are proposed. Finally, some numerical results are presented and compared with theoretical results. The numerical results of the proposed algorithm conform with those of the FE algorithm well, demonstrating the accuracy of the present method and its potential applications in thermal engineering of porous materials.
基金supported by the National Natural Science Foundation of China(Grants 11471262,11202032)the Basic Research Project of National Defense(Grant B 1520132013)supported by the State Key Laboratory of Science and Engineering Computing and Center for high performance computing of Northwestem Polytechnical University
文摘In this paper,a statistical second-order twoscale（SSOTS） method is developed to simulate the dynamic thcrmo-mechanical performances of the statistically inhomogeneous materials.For this kind of composite material,the random distribution characteristics of particles,including the shape,size,orientation,spatial location,and volume fractions,are all considered.Firstly,the repre.sentation for the microscopic configuration of the statistically inhomogeneous materials is described.Secondly,the SSOTS formulation for the dynamic thermo-mechanical coupled problem is proposed in a constructive way,including the cell problems,effective thermal and mechanical parameters,homogenized problems,and the SSOTS formulas of the temperatures,displacements,heat flux densities and stresses.And then the algorithm procedure corresponding to the SSOTS method is brought forward.The numerical results obtained by using the SSOTS algorithm are compared with those by classical methods.In addition,the thermo-mechanical coupling effect is studied by comparing the results of coupled case with those of uncoupled case.It demonstrates that the coupling effect on the temperatures,heat flux densities,displacements,and stresses is very distinct.The results show that the SSOTS method is valid to predict the dynamic thermo-mechanical coupled performances of statistically inhomogeneous materials.
基金supported by the Special Funds for the National Basic Research Program of China(Grant No.2012CB025904)the National Natural ScienceFoundation of China(Grant Nos.90916027 and 11302052)
文摘This paper focuses on the dynamic thermo-mechanical coupled response of random particulate composite materials. Both the inertia term and coupling term are considered in the dynamic coupled problem. The formulation of the problem by a statistical second-order two-scale （SSOTS） analysis method and the algorithm procedure based on the finite-element difference method are presented. Numerical results of coupled cases are compared with those of uncoupled cases. It shows that the coupling effects on temperature, thermal flux, displacement, and stresses are very distinct, and the micro- characteristics of particles affect the coupling effect of the random composites. Furthermore, the coupling effect causes a lag in the variations of temperature, thermal flux, displacement, and stresses.
基金NSFC under Project 1967108 and Croucher Foundation of Hong Kong, Supported also by FRGof Hong Kong Baptist University.
文摘Two general local Cm triangular interpolation schemes by rational functions from Cm data are proposed for any nonnegative integer m. The schemes can have either 2m+1 order algebraic precision if the required data are given on vertices and edges, or m+E[m/2]+1 or m+1 order algebraic precision if the data are given only at vertices. The orders of the interpolation error are estimated. Examples that show the correctness and effectiveness of the scheme are presented.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.91630208,91641107,and 11771437)
文摘A new variational method is proposed to investigate the dynamics of the thin film in a coating flow where a liquid is delivered through a fixed slot gap onto a moving substrate. A simplified ODE system has also been derived for the evolution of the thin film whose thickness hf is asymptotically constant behind the coating front. We calculate the phase diagram as well as the film profiles and approximate the film thickness theoretically, and agreement with the well-known scaling law as Ca2/3 is found.
