In oil and mineral exploration, gravity gradient tensor data include higher- frequency signals than gravity data, which can be used to delineate small-scale anomalies. However, full-tensor gradiometry (FTG) data are...In oil and mineral exploration, gravity gradient tensor data include higher- frequency signals than gravity data, which can be used to delineate small-scale anomalies. However, full-tensor gradiometry (FTG) data are contaminated by high-frequency random noise. The separation of noise from high-frequency signals is one of the most challenging tasks in processing of gravity gradient tensor data. We first derive the Cartesian equations of gravity gradient tensors under the constraint of the Laplace equation and the expression for the gravitational potential, and then we use the Cartesian equations to fit the measured gradient tensor data by using optimal linear inversion and remove the noise from the measured data. Based on model tests, we confirm that not only this method removes the high- frequency random noise but also enhances the weak anomaly signals masked by the noise. Compared with traditional low-pass filtering methods, this method avoids removing noise by sacrificing resolution. Finally, we apply our method to real gravity gradient tensor data acquired by Bell Geospace for the Vinton Dome at the Texas-Louisiana border.展开更多
We propose an efficient and robust algorithm to solve the steady Euler equa- tions on unstructured grids.The new algorithm is a Newton-iteration method in which each iteration step is a linear multigrid method using b...We propose an efficient and robust algorithm to solve the steady Euler equa- tions on unstructured grids.The new algorithm is a Newton-iteration method in which each iteration step is a linear multigrid method using block lower-upper symmetric Gauss-Seidel(LU-SGS)iteration as its smoother To regularize the Jacobian matrix of Newton-iteration,we adopted a local residual dependent regularization as the replace- ment of the standard time-stepping relaxation technique based on the local CFL number The proposed method can be extended to high order approximations and three spatial dimensions in a nature way.The solver was tested on a sequence of benchmark prob- lems on both quasi-uniform and local adaptive meshes.The numerical results illustrated the efficiency and robustness of our algorithm.展开更多
Tensile stress-strain curves of five metallic alloys,i.e.,SKH51,STS316L,Ti-6Al-4V,Al6061and Inconel600were analyzed to investigate the working hardening behavior.The constitutive parameters of three constitutive equat...Tensile stress-strain curves of five metallic alloys,i.e.,SKH51,STS316L,Ti-6Al-4V,Al6061and Inconel600were analyzed to investigate the working hardening behavior.The constitutive parameters of three constitutive equations,i.e.,the Hollomon,Swift and Voce equations,were compared by using different methods.A new working hardening parameter was proposed to characterize the working hardening behavior in different deformation stages.It is found that Voce equation is suitable to describe stress-strain curves in large strain region.Meanwhile,the predicting accuracy of ultimate tensile strength by Voce equation is the best.The working hardening behavior of SKH51is different from the other four metallic alloys.展开更多
A new approach is proposed to analyze the settlement behavior for single pile embedded in layered soils. Firstly, soil layers surrounding pile shaft are simulated by using distributed Voigt model, and finite soil laye...A new approach is proposed to analyze the settlement behavior for single pile embedded in layered soils. Firstly, soil layers surrounding pile shaft are simulated by using distributed Voigt model, and finite soil layers under the pile end are assumed to be virtual soil-pile whose cross-section area is the same as that of the pile shaft. Then, by means of Laplace transform and impedance function transfer method to solve the static equilibrium equation of pile, the analytical solution of the displacement impedance fimction at the pile head is derived. Furthermore, the analytical solution of the settlement at the head of single pile is theoretically derived by virtue of convolution theorem. Based on these solutions, the influences of parameters of soil-pile system on the settlement behavior for single pile are analyzed. Also, comparison of the load-settlement response for two well-instrumented field tests in multilayered soils is given to demonstrate the effectiveness and accuracy of the proposed approach. It can be noted that the presented solution can be used to calculate the settlement of single pile for the preliminary design of pile foundation.展开更多
We find that in a supersymmetric quantum mechanics (SUSY QM) system, in addition to supersymmetric algebra, an associated SU(2) algebra can be obtained by using semiunitary (SUT) operator and projection operator...We find that in a supersymmetric quantum mechanics (SUSY QM) system, in addition to supersymmetric algebra, an associated SU(2) algebra can be obtained by using semiunitary (SUT) operator and projection operator, and the relevant constants of motion can be constructed. Two typical quantum systems are investigated as examples to demonstrate the above finding. The first example is the quantum system of a nonrelativistic charged particle moving in x-y plane and coupled to a magnetic field along z-axis. The second example is provided with the Dirac particle in a magnetic field. Similarly there exists an SUτ(2) SUσ(2) symmetry in the context of the relativistic Pauli Hamilt onian squared. We show that there exists also an SU(2) symmetry associated with the supersvmmetrv of the Dirac particle.展开更多
The Dirac equations with vector and scalar potentials of the Coulomb types in two and three dimensions are solved using the supersymmetric quantum mechanics method. For the system of such potentials, the analytical ex...The Dirac equations with vector and scalar potentials of the Coulomb types in two and three dimensions are solved using the supersymmetric quantum mechanics method. For the system of such potentials, the analytical expressions of the matrix dements for both position and momentum operators are obtained.展开更多
Based on unified theory of electromagnetic interactions and gravitational interactions, the non-relativistic limit of the equation of motion of a charged Dirac particle in gravitational field is studied. From the Schr...Based on unified theory of electromagnetic interactions and gravitational interactions, the non-relativistic limit of the equation of motion of a charged Dirac particle in gravitational field is studied. From the Schroedinger equation obtained from this non-relativistic limit, we can see that the classical Newtonian gravitational potential appears as a part of the potential in the Schroedinger equation, which can explain the gravitational phase effects found in COW experiments. And because of this Newtonian gravitational potential, a quantum particle in the earth's gravitational field may form a gravitationally bound quantized state, which has already been detected in experiments. Three different kinds of phase effects related to gravitational interactions are studied in this paper, and these phase effects should be observable in some astrophysical processes. Besides, there exists direct coupling between gravitomagnetic field and quantum spin, and radiation caused by this coupling can be used to directly determine the gravitomagnetic field on the surface of a star.展开更多
Tensile impact tests of aramid (Twaron) fiber bundles were carried out under high strain rates with a wide range of 0. 01/s -1 000/s by using MTS and bar-bar tensile impact apparatus. Based on the statistical constitu...Tensile impact tests of aramid (Twaron) fiber bundles were carried out under high strain rates with a wide range of 0. 01/s -1 000/s by using MTS and bar-bar tensile impact apparatus. Based on the statistical constitutive model of fiber bundles, statistical constitutive equations of aramid fiber bundles are derived from statistical analysis of test data at different strain rates. Comparison between the theoretical predictions and experimental data indicates statistical constitutive equations fit well with the experimental data, and statistical constitutive equations of fiber bundles at different strain rates are valid.展开更多
We present a general approach to the construction of conservation laws for the nonholonomic singular Lagrange system. Firstly, the differential equations of motion of the system are written, the definition of integrat...We present a general approach to the construction of conservation laws for the nonholonomic singular Lagrange system. Firstly, the differential equations of motion of the system are written, the definition of integrating factors is given for the system. Next, the necessary conditions for the existence of the conserved quantity are studied in detail. Finally, the conservation theorem and its inverse are established for the system, an example is given to illustrate the application of the result.展开更多
Gravitational field produced by high-power laser is calculated according to the linearized Einstein field equation in weak field approximation. Gravitational Faraday effect of electromagnetic wave propagating in the a...Gravitational field produced by high-power laser is calculated according to the linearized Einstein field equation in weak field approximation. Gravitational Faraday effect of electromagnetic wave propagating in the above gravitational field is studied and the rotation angle of polarization plane of electromagnetic wave is derived. The result is discussed and estimated under the condition of present experiment facility.展开更多
In order to present a new method for analyzing the reliability of a two-link flexible robot manipulator,Lagrange dynamics differential equations of the two-link flexible robot manipulator were established by using the...In order to present a new method for analyzing the reliability of a two-link flexible robot manipulator,Lagrange dynamics differential equations of the two-link flexible robot manipulator were established by using the integrated modal method and the multi-body system dynamics method.By using the Monte Carlo method,the random sample values of the dynamic parameters were obtained and Lagrange dynamics differential equations were solved for each random sample value which revealed their displacement,speed and acceleration.On this basis,dynamic stresses and deformations were obtained.By taking the maximum values of the stresses and the deformations as output responses and the random sample values of dynamic parameters as input quantities,extremum response surface functions were established.A number of random samples were then obtained by using the Monte Carlo method and then the reliability was analyzed by using the extremum response surface method.The results show that the extremum response surface method is an efficient and fast reliability analysis method with high-accuracy for the two-link flexible robot manipulator.展开更多
In recent years,defunct satellites mitigation in the geostationary orbit(GEO) has become a hot issue in the space field.How to transfer defunct geostationary satellites to the graveyard orbit safely,economically and e...In recent years,defunct satellites mitigation in the geostationary orbit(GEO) has become a hot issue in the space field.How to transfer defunct geostationary satellites to the graveyard orbit safely,economically and efficiently presents new challenges to spacecraft dynamics and control.This paper conducts an in-depth investigation on tether-tugging de-orbit issues of defunct geostationary satellites.Firstly,a four-phase tether-tugging de-orbit scheme including acceleration,equilibrium,rotation and return is proposed.This scheme takes into consideration how to avoid the risks of tether ripping,tug-target collision,and tether twist,and how to achieve the mission objective of fuel saving.Secondly,the dynamics model of the tether combination system is established based on Lagrange equation,and the four phases of tether-tugging de-orbit scheme are simulated respectively.Simulation results indicate that the scheme is theoretically feasible and satisfies the design objectives of safety,economy and efficiency,providing a technical approach for engineering application.展开更多
Regional gravity field modeling with high-precision and high-resolution is one of the most important scientific objectives in geodesy, and can provide fundamental information for geophysics, geodynamics, seismology, a...Regional gravity field modeling with high-precision and high-resolution is one of the most important scientific objectives in geodesy, and can provide fundamental information for geophysics, geodynamics, seismology, and mineral exploration. Rectangular harmonic analysis (RHA) is proposed for regional gravity field modeling in this paper. By solving the Laplace's equation of gravitational potential in local Cartesian coordinate system, the rectangular harmonic expansions of disturbing potential, gravity anomaly, gravity disturbance, geoid undulation and deflection of the vertical are derived, and so are the formula for signal degree variance and error degree variance of the rectangular harmonic coefficients (RHC). We also present the mathematical model and detailed algorithm for the solution of RHC using RHA from gravity observations. In order to reduce the edge effects caused by periodic continuation in RHA, we propose the strategy of extending the size of computation domain. The RHA-based modeling method is validated by conducting numerical experiments based on simulated ground and airborne gravity data that are generated from geopotential model EGM2008 and contaminated by Gauss white noise with standard deviation of 2 mGal. The accuracy of the 2.5'×2.5' geoid undulations computed from ground and airborne gravity data is 1 and 1.4 cm, respectively. The standard error of the gravity disturbances that downward continued from the flight height of 4 km to the geoid is only 3.1 reGal. Numerical results confirm that RHA is able to provide a reliable and accurate regional gravity field model, which may be a new option for the representation of the fine structure of regional gravity field.展开更多
A consistent focus in theoretical mechanics has been on how to apply Lagrange's equation to continuum mechanics.This paper uses the concept of a variational derivative and its laws of operation to investigate the ...A consistent focus in theoretical mechanics has been on how to apply Lagrange's equation to continuum mechanics.This paper uses the concept of a variational derivative and its laws of operation to investigate the derivation of Lagrange's equation,which is then applied to nonlinear elasto-dynamics.In accordance with the work-energy principle and the energy conservation law,kinetic and potential energies are proposed for rigid-elastic coupling dynamics,whose governing equation is established by manipulating Lagrange's equation.In addition,case studies are used to demonstrate the application of the proposed method to spacecraft dynamics.展开更多
In this paper,we solve the Dirac equation under spin symmetry limit for attractive radial potential including a Coulomb-like tensor interaction.By using the parametric generalization of the Nikiforov-Uvarov method,the...In this paper,we solve the Dirac equation under spin symmetry limit for attractive radial potential including a Coulomb-like tensor interaction.By using the parametric generalization of the Nikiforov-Uvarov method,the energy eigenvalues equation and the corresponding wave functions have been obtained in closed forms.Some numerical results are given too.展开更多
In this paper,firstly,by solving the Riemann problem of the zero-pressure flow in gas dynamics with a flux approximation,we construct parameterized delta-shock and constant density solutions,then we show that,as the f...In this paper,firstly,by solving the Riemann problem of the zero-pressure flow in gas dynamics with a flux approximation,we construct parameterized delta-shock and constant density solutions,then we show that,as the flux perturbation vanishes,they converge to the delta-shock and vacuum state solutions of the zero-pressure flow,respectively.Secondly,we solve the Riemann problem of the Euler equations of isentropic gas dynamics with a double parameter flux approximation including pressure.Furthermore,we rigorously prove that,as the two-parameter flux perturbation vanishes,any Riemann solution containing two shock waves tends to a delta-shock solution to the zero-pressure flow;any Riemann solution containing two rarefaction waves tends to a two-contact-discontinuity solution to the zero-pressure flow and the nonvacuum intermediate state in between tends to a vacuum state.Finally,numerical results are given to present the formation processes of delta shock waves and vacuum states.展开更多
Based on an upwind compact difference scheme and the idea of monotonicity-preserving, a 5th order monotonicity-preserving upwind compact difference scheme (m-UCD5) is proposed. The new difference scheme not only ret...Based on an upwind compact difference scheme and the idea of monotonicity-preserving, a 5th order monotonicity-preserving upwind compact difference scheme (m-UCD5) is proposed. The new difference scheme not only retains the advantage of good resolution of high wave number but also avoids the Gibbs phenomenon of the original upwind compact difference scheme. Compared with the classical 5th order WENO difference scheme, the new difference scheme is simpler and small in diffusion and computation load. By employing the component-wise and characteristic-wise methods, two forms of the new difference scheme are proposed to solve the N-S/Euler equation. Through the Sod problem, the Shu-Osher problem and tbe two-dimensional Double Mach Reflection problem, numerical solutions have demonstrated this new scheme does have a good resolution of high wave number and a robust ability of capturing shock waves, leading to a conclusion that the new difference scheme may be used to simulate complex flows containing shock waves.展开更多
基金financially supported by the SinoProbe-09-01(201011078)
文摘In oil and mineral exploration, gravity gradient tensor data include higher- frequency signals than gravity data, which can be used to delineate small-scale anomalies. However, full-tensor gradiometry (FTG) data are contaminated by high-frequency random noise. The separation of noise from high-frequency signals is one of the most challenging tasks in processing of gravity gradient tensor data. We first derive the Cartesian equations of gravity gradient tensors under the constraint of the Laplace equation and the expression for the gravitational potential, and then we use the Cartesian equations to fit the measured gradient tensor data by using optimal linear inversion and remove the noise from the measured data. Based on model tests, we confirm that not only this method removes the high- frequency random noise but also enhances the weak anomaly signals masked by the noise. Compared with traditional low-pass filtering methods, this method avoids removing noise by sacrificing resolution. Finally, we apply our method to real gravity gradient tensor data acquired by Bell Geospace for the Vinton Dome at the Texas-Louisiana border.
文摘We propose an efficient and robust algorithm to solve the steady Euler equa- tions on unstructured grids.The new algorithm is a Newton-iteration method in which each iteration step is a linear multigrid method using block lower-upper symmetric Gauss-Seidel(LU-SGS)iteration as its smoother To regularize the Jacobian matrix of Newton-iteration,we adopted a local residual dependent regularization as the replace- ment of the standard time-stepping relaxation technique based on the local CFL number The proposed method can be extended to high order approximations and three spatial dimensions in a nature way.The solver was tested on a sequence of benchmark prob- lems on both quasi-uniform and local adaptive meshes.The numerical results illustrated the efficiency and robustness of our algorithm.
基金Project(51275414)supported by the National Natural Science Foundation of ChinaProject(3102015BJ(Ⅱ)ZS007)supported by the Fundamental Research Funds for the Central Universities,ChinaProject(130-QP-2015)supported by the Research Fund of the State Key Laboratory of Solidification Processing(NWPU),China
文摘Tensile stress-strain curves of five metallic alloys,i.e.,SKH51,STS316L,Ti-6Al-4V,Al6061and Inconel600were analyzed to investigate the working hardening behavior.The constitutive parameters of three constitutive equations,i.e.,the Hollomon,Swift and Voce equations,were compared by using different methods.A new working hardening parameter was proposed to characterize the working hardening behavior in different deformation stages.It is found that Voce equation is suitable to describe stress-strain curves in large strain region.Meanwhile,the predicting accuracy of ultimate tensile strength by Voce equation is the best.The working hardening behavior of SKH51is different from the other four metallic alloys.
基金Project(50879077) supported by the National Natural Science Foundation of China
文摘A new approach is proposed to analyze the settlement behavior for single pile embedded in layered soils. Firstly, soil layers surrounding pile shaft are simulated by using distributed Voigt model, and finite soil layers under the pile end are assumed to be virtual soil-pile whose cross-section area is the same as that of the pile shaft. Then, by means of Laplace transform and impedance function transfer method to solve the static equilibrium equation of pile, the analytical solution of the displacement impedance fimction at the pile head is derived. Furthermore, the analytical solution of the settlement at the head of single pile is theoretically derived by virtue of convolution theorem. Based on these solutions, the influences of parameters of soil-pile system on the settlement behavior for single pile are analyzed. Also, comparison of the load-settlement response for two well-instrumented field tests in multilayered soils is given to demonstrate the effectiveness and accuracy of the proposed approach. It can be noted that the presented solution can be used to calculate the settlement of single pile for the preliminary design of pile foundation.
基金supported by National Natural Science Foundation of China under Grant Nos.10375039 and 90503008the Doctoral Fund of Ministry of Education of Chinapartly by the Center of Theoretical Nuclear Physics of HIRFL of China
文摘We find that in a supersymmetric quantum mechanics (SUSY QM) system, in addition to supersymmetric algebra, an associated SU(2) algebra can be obtained by using semiunitary (SUT) operator and projection operator, and the relevant constants of motion can be constructed. Two typical quantum systems are investigated as examples to demonstrate the above finding. The first example is the quantum system of a nonrelativistic charged particle moving in x-y plane and coupled to a magnetic field along z-axis. The second example is provided with the Dirac particle in a magnetic field. Similarly there exists an SUτ(2) SUσ(2) symmetry in the context of the relativistic Pauli Hamilt onian squared. We show that there exists also an SU(2) symmetry associated with the supersvmmetrv of the Dirac particle.
基金National Natural Science Foundation of China under Grant Nos.10125521 and 60371013the 973 State Key Basic Research Development Project of China under Grant No.G2000077400
文摘The Dirac equations with vector and scalar potentials of the Coulomb types in two and three dimensions are solved using the supersymmetric quantum mechanics method. For the system of such potentials, the analytical expressions of the matrix dements for both position and momentum operators are obtained.
文摘Based on unified theory of electromagnetic interactions and gravitational interactions, the non-relativistic limit of the equation of motion of a charged Dirac particle in gravitational field is studied. From the Schroedinger equation obtained from this non-relativistic limit, we can see that the classical Newtonian gravitational potential appears as a part of the potential in the Schroedinger equation, which can explain the gravitational phase effects found in COW experiments. And because of this Newtonian gravitational potential, a quantum particle in the earth's gravitational field may form a gravitationally bound quantized state, which has already been detected in experiments. Three different kinds of phase effects related to gravitational interactions are studied in this paper, and these phase effects should be observable in some astrophysical processes. Besides, there exists direct coupling between gravitomagnetic field and quantum spin, and radiation caused by this coupling can be used to directly determine the gravitomagnetic field on the surface of a star.
基金The project is supported by Zhejiang Provincial Natural Science Foundaion of China(599113)Science and Technology Foundation of Ministy of Educationd of China(DF 02064)
文摘Tensile impact tests of aramid (Twaron) fiber bundles were carried out under high strain rates with a wide range of 0. 01/s -1 000/s by using MTS and bar-bar tensile impact apparatus. Based on the statistical constitutive model of fiber bundles, statistical constitutive equations of aramid fiber bundles are derived from statistical analysis of test data at different strain rates. Comparison between the theoretical predictions and experimental data indicates statistical constitutive equations fit well with the experimental data, and statistical constitutive equations of fiber bundles at different strain rates are valid.
基金The project supported by National Natural Science Foundation of China under Grant No. 10272034 and the Doctoral Program Foundation of China under Grnt No. 20030558025
文摘We present a general approach to the construction of conservation laws for the nonholonomic singular Lagrange system. Firstly, the differential equations of motion of the system are written, the definition of integrating factors is given for the system. Next, the necessary conditions for the existence of the conserved quantity are studied in detail. Finally, the conservation theorem and its inverse are established for the system, an example is given to illustrate the application of the result.
文摘Gravitational field produced by high-power laser is calculated according to the linearized Einstein field equation in weak field approximation. Gravitational Faraday effect of electromagnetic wave propagating in the above gravitational field is studied and the rotation angle of polarization plane of electromagnetic wave is derived. The result is discussed and estimated under the condition of present experiment facility.
基金Project(2006AA04Z405) supported by the National High Technology Research and Development Program of ChinaProject(3102019) supported by Beijing Municipal Natural Science Foundation,China
文摘In order to present a new method for analyzing the reliability of a two-link flexible robot manipulator,Lagrange dynamics differential equations of the two-link flexible robot manipulator were established by using the integrated modal method and the multi-body system dynamics method.By using the Monte Carlo method,the random sample values of the dynamic parameters were obtained and Lagrange dynamics differential equations were solved for each random sample value which revealed their displacement,speed and acceleration.On this basis,dynamic stresses and deformations were obtained.By taking the maximum values of the stresses and the deformations as output responses and the random sample values of dynamic parameters as input quantities,extremum response surface functions were established.A number of random samples were then obtained by using the Monte Carlo method and then the reliability was analyzed by using the extremum response surface method.The results show that the extremum response surface method is an efficient and fast reliability analysis method with high-accuracy for the two-link flexible robot manipulator.
基金supported by the National Hi-Tech Research and Development Program of China ("863" Project) (Grant No. 2011AA7044026)
文摘In recent years,defunct satellites mitigation in the geostationary orbit(GEO) has become a hot issue in the space field.How to transfer defunct geostationary satellites to the graveyard orbit safely,economically and efficiently presents new challenges to spacecraft dynamics and control.This paper conducts an in-depth investigation on tether-tugging de-orbit issues of defunct geostationary satellites.Firstly,a four-phase tether-tugging de-orbit scheme including acceleration,equilibrium,rotation and return is proposed.This scheme takes into consideration how to avoid the risks of tether ripping,tug-target collision,and tether twist,and how to achieve the mission objective of fuel saving.Secondly,the dynamics model of the tether combination system is established based on Lagrange equation,and the four phases of tether-tugging de-orbit scheme are simulated respectively.Simulation results indicate that the scheme is theoretically feasible and satisfies the design objectives of safety,economy and efficiency,providing a technical approach for engineering application.
基金jointly supported by the National Basic Research Program of China (Grant No. 2013CB733301)the National Science and Technology Support Program of China (Grant No. 2012BAB16B01)+1 种基金the National Natural Science Foundation of China (Grant No. 41204008)the Basic Research Program of National Administration of Surveying, Mapping and Geoinformation of China
文摘Regional gravity field modeling with high-precision and high-resolution is one of the most important scientific objectives in geodesy, and can provide fundamental information for geophysics, geodynamics, seismology, and mineral exploration. Rectangular harmonic analysis (RHA) is proposed for regional gravity field modeling in this paper. By solving the Laplace's equation of gravitational potential in local Cartesian coordinate system, the rectangular harmonic expansions of disturbing potential, gravity anomaly, gravity disturbance, geoid undulation and deflection of the vertical are derived, and so are the formula for signal degree variance and error degree variance of the rectangular harmonic coefficients (RHC). We also present the mathematical model and detailed algorithm for the solution of RHC using RHA from gravity observations. In order to reduce the edge effects caused by periodic continuation in RHA, we propose the strategy of extending the size of computation domain. The RHA-based modeling method is validated by conducting numerical experiments based on simulated ground and airborne gravity data that are generated from geopotential model EGM2008 and contaminated by Gauss white noise with standard deviation of 2 mGal. The accuracy of the 2.5'×2.5' geoid undulations computed from ground and airborne gravity data is 1 and 1.4 cm, respectively. The standard error of the gravity disturbances that downward continued from the flight height of 4 km to the geoid is only 3.1 reGal. Numerical results confirm that RHA is able to provide a reliable and accurate regional gravity field model, which may be a new option for the representation of the fine structure of regional gravity field.
基金supported by the National Natural Science Foundation of China(Grant No.10272034)
文摘A consistent focus in theoretical mechanics has been on how to apply Lagrange's equation to continuum mechanics.This paper uses the concept of a variational derivative and its laws of operation to investigate the derivation of Lagrange's equation,which is then applied to nonlinear elasto-dynamics.In accordance with the work-energy principle and the energy conservation law,kinetic and potential energies are proposed for rigid-elastic coupling dynamics,whose governing equation is established by manipulating Lagrange's equation.In addition,case studies are used to demonstrate the application of the proposed method to spacecraft dynamics.
文摘In this paper,we solve the Dirac equation under spin symmetry limit for attractive radial potential including a Coulomb-like tensor interaction.By using the parametric generalization of the Nikiforov-Uvarov method,the energy eigenvalues equation and the corresponding wave functions have been obtained in closed forms.Some numerical results are given too.
基金supported by National Natural Science Foundation of China(Grant No.11361073)
文摘In this paper,firstly,by solving the Riemann problem of the zero-pressure flow in gas dynamics with a flux approximation,we construct parameterized delta-shock and constant density solutions,then we show that,as the flux perturbation vanishes,they converge to the delta-shock and vacuum state solutions of the zero-pressure flow,respectively.Secondly,we solve the Riemann problem of the Euler equations of isentropic gas dynamics with a double parameter flux approximation including pressure.Furthermore,we rigorously prove that,as the two-parameter flux perturbation vanishes,any Riemann solution containing two shock waves tends to a delta-shock solution to the zero-pressure flow;any Riemann solution containing two rarefaction waves tends to a two-contact-discontinuity solution to the zero-pressure flow and the nonvacuum intermediate state in between tends to a vacuum state.Finally,numerical results are given to present the formation processes of delta shock waves and vacuum states.
基金supported by the National Natural Science Foundation of China (Grant Nos. 110632050, 10872205)the National Basic Research Program of China (Grant No. 2009CB724100)Projects of CAS INFO-115-B01
文摘Based on an upwind compact difference scheme and the idea of monotonicity-preserving, a 5th order monotonicity-preserving upwind compact difference scheme (m-UCD5) is proposed. The new difference scheme not only retains the advantage of good resolution of high wave number but also avoids the Gibbs phenomenon of the original upwind compact difference scheme. Compared with the classical 5th order WENO difference scheme, the new difference scheme is simpler and small in diffusion and computation load. By employing the component-wise and characteristic-wise methods, two forms of the new difference scheme are proposed to solve the N-S/Euler equation. Through the Sod problem, the Shu-Osher problem and tbe two-dimensional Double Mach Reflection problem, numerical solutions have demonstrated this new scheme does have a good resolution of high wave number and a robust ability of capturing shock waves, leading to a conclusion that the new difference scheme may be used to simulate complex flows containing shock waves.