For real-time dynamic substructure testing(RTDST),the influence of the inertia force of fluid specimens on the stability and accuracy of the integration algorithms has never been investigated.Therefore,this study prop...For real-time dynamic substructure testing(RTDST),the influence of the inertia force of fluid specimens on the stability and accuracy of the integration algorithms has never been investigated.Therefore,this study proposes to investigate the stability and accuracy of the central difference method(CDM)for RTDST considering the specimen mass participation coefficient.First,the theory of the CDM for RTDST is presented.Next,the stability and accuracy of the CDM for RTDST considering the specimen mass participation coefficient are investigated.Finally,numerical simulations and experimental tests are conducted for verifying the effectiveness of the method.The study indicates that the stability of the algorithm is affected by the mass participation coefficient of the specimen,and the stability limit first increases and then decreases as the mass participation coefficient increases.In most cases,the mass participation coefficient will increase the stability limit of the algorithm,but in specific circumstances,the algorithm may lose its stability.The stability and accuracy of the CDM considering the mass participation coefficient are verified by numerical simulations and experimental tests on a three-story frame structure with a tuned liquid damper.展开更多
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.
A class of multidomain hybrid methods of direct discontinuous Galerkin(DDG)methods and central difference(CD)schemes for the viscous terms is pro-posed in this paper.Both conservative and nonconservative coupling mode...A class of multidomain hybrid methods of direct discontinuous Galerkin(DDG)methods and central difference(CD)schemes for the viscous terms is pro-posed in this paper.Both conservative and nonconservative coupling modes are dis-cussed.To treat the shock wave,the nonconservative coupling mode automatically switch to conservative coupling mode to preserve the conservative property when discontinuities pass through the artificial interface.To maintain the accuracy of the hybrid methods,the Lagrange interpolation polynomials and their derivatives are reconstructed to handle the coupling cells in the DDG subdomain,while the values of ghost points for the CD subdomain are calculated by the approximate polynomials from the DDG methods.The linear stabilities of these methods are demonstrated in detail through von-Neumann analysis.The multidomain hybrid DDG and CD meth-ods are then extended to one-and two-dimensional hyperbolic-parabolic equations.Numerical results validate that the multidomain hybrid methods are high-order ac-curate in the smooth regions,robust for viscous shock simulations and capable to save computational cost.展开更多
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.展开更多
Determining the venting time of a gas trunk pipeline segment provides an important basis for formulating an emergency plan in the advent of unexpected accidents.As the natural gas venting process corresponds to the tr...Determining the venting time of a gas trunk pipeline segment provides an important basis for formulating an emergency plan in the advent of unexpected accidents.As the natural gas venting process corresponds to the transient flow,it is necessary to establish a transient hydraulic-thermal simulation model in order to determine the venting time.In this paper,based on two kinds of venting scenarios in which there is only one venting point in the venting system of a gas trunk pipeline segment—namely,where the venting point is either at one of the two ends or at the junction of two gas trunk pipeline segments—transient hydraulic-thermal simulation models are established.The models consist of gas flow governing equations,the gas state equation,gas physical property equations,initial conditions,and appropriate boundary conditions.The implicit central difference method is used to discretize the gas flow partial differential equations,and the trust-region-dogleg algorithm is used to solve the equations corresponding to each time step,in order to dynamically simulate the whole venting process.The judgment condition for the end of the venting process is that the average pressure of gas trunk pipeline segment is less than 0.11 MPa(actual pressure).Comparing the simulation results of the proposed model with those of the OLGA software and real operational data,we find that the venting time error is less than 10%.On this basis,a venting valve opening control principle is proposed,which prevents the venting noise from exceeded the specified noise value(85 d B)in the venting design of domestic gas pipeline projects.The established calculation model for venting time(dynamic simulation model)for a gas trunk pipeline segment and the proposed opening control principle of venting valve provide reference for the optimal operation of gas pipeline venting systems.展开更多
In the practical problems such as nuclear waste pollution and seawater intrusion etc., many problems are reduced to solving the convection-diffusion equation, so the research of convection-diffusion equation is of gre...In the practical problems such as nuclear waste pollution and seawater intrusion etc., many problems are reduced to solving the convection-diffusion equation, so the research of convection-diffusion equation is of great value. In this work, a spectral method is presented for solving one and two dimensional convection-diffusion equation with source term. The finite difference method is also used to solve the convection diffusion equation. The numerical experiments show that the spectral method is more efficient than other methods for solving the convection-diffusion equation.展开更多
A numerical treatment for the prediction of cavitating flows is presented and assessed. The algorithm uses the preconditioned multiphase Euler equations with appropriate mass transfer terms. A central difference finit...A numerical treatment for the prediction of cavitating flows is presented and assessed. The algorithm uses the preconditioned multiphase Euler equations with appropriate mass transfer terms. A central difference finite volume scheme with suitable dissipation terms to account for density jumps across the cavity interface is shown to yield an effective method for solving the multiphase Euler equations. The Euler equations are utilized herein for the cavitation modeling, because some certain characteristics of cavitating flows can be obtained using the solution of this system of equations with relative low computational effort. In addition, the Euler equations are appropriate for the assessment of the numerical method used, because of the sensitivity of the solution to the numerical instabilities. For this reason, a sensitivity study is conducted to evaluate the effects of various parameters, such as numerical dissipation coefficients and grid size, on the accuracy and performance of the solution. The computations are performed for steady cavitating flows around the NACA 0012 and NACA 66 （MOD） hydrofoils and also an axisymmetric hemispherical fore-body under different conditions and the results are compared with the available numerical and experimental data. The solution procedure presented is shown to be accurate and efficient for predicting steady sheet- and super-cavitation for 2D/axisymmetric geometries.展开更多
基金National Natural Science Foundation of China under Grant Nos.51978213 and 51778190the National Key Research and Development Program of China under Grant Nos.2017YFC0703605 and 2016YFC0701106。
文摘For real-time dynamic substructure testing(RTDST),the influence of the inertia force of fluid specimens on the stability and accuracy of the integration algorithms has never been investigated.Therefore,this study proposes to investigate the stability and accuracy of the central difference method(CDM)for RTDST considering the specimen mass participation coefficient.First,the theory of the CDM for RTDST is presented.Next,the stability and accuracy of the CDM for RTDST considering the specimen mass participation coefficient are investigated.Finally,numerical simulations and experimental tests are conducted for verifying the effectiveness of the method.The study indicates that the stability of the algorithm is affected by the mass participation coefficient of the specimen,and the stability limit first increases and then decreases as the mass participation coefficient increases.In most cases,the mass participation coefficient will increase the stability limit of the algorithm,but in specific circumstances,the algorithm may lose its stability.The stability and accuracy of the CDM considering the mass participation coefficient are verified by numerical simulations and experimental tests on a three-story frame structure with a tuned liquid damper.
文摘In this paper, we construct a uniform second-order difference scheme for a class of boundary value problems of fourth-order ordinary differential equations. Finally, a numerical example is given.
基金supported by the National Natural Science Foundation of China(Grant No.12001031).
文摘A class of multidomain hybrid methods of direct discontinuous Galerkin(DDG)methods and central difference(CD)schemes for the viscous terms is pro-posed in this paper.Both conservative and nonconservative coupling modes are dis-cussed.To treat the shock wave,the nonconservative coupling mode automatically switch to conservative coupling mode to preserve the conservative property when discontinuities pass through the artificial interface.To maintain the accuracy of the hybrid methods,the Lagrange interpolation polynomials and their derivatives are reconstructed to handle the coupling cells in the DDG subdomain,while the values of ghost points for the CD subdomain are calculated by the approximate polynomials from the DDG methods.The linear stabilities of these methods are demonstrated in detail through von-Neumann analysis.The multidomain hybrid DDG and CD meth-ods are then extended to one-and two-dimensional hyperbolic-parabolic equations.Numerical results validate that the multidomain hybrid methods are high-order ac-curate in the smooth regions,robust for viscous shock simulations and capable to save computational cost.
文摘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(Grant No.52174064)
文摘Determining the venting time of a gas trunk pipeline segment provides an important basis for formulating an emergency plan in the advent of unexpected accidents.As the natural gas venting process corresponds to the transient flow,it is necessary to establish a transient hydraulic-thermal simulation model in order to determine the venting time.In this paper,based on two kinds of venting scenarios in which there is only one venting point in the venting system of a gas trunk pipeline segment—namely,where the venting point is either at one of the two ends or at the junction of two gas trunk pipeline segments—transient hydraulic-thermal simulation models are established.The models consist of gas flow governing equations,the gas state equation,gas physical property equations,initial conditions,and appropriate boundary conditions.The implicit central difference method is used to discretize the gas flow partial differential equations,and the trust-region-dogleg algorithm is used to solve the equations corresponding to each time step,in order to dynamically simulate the whole venting process.The judgment condition for the end of the venting process is that the average pressure of gas trunk pipeline segment is less than 0.11 MPa(actual pressure).Comparing the simulation results of the proposed model with those of the OLGA software and real operational data,we find that the venting time error is less than 10%.On this basis,a venting valve opening control principle is proposed,which prevents the venting noise from exceeded the specified noise value(85 d B)in the venting design of domestic gas pipeline projects.The established calculation model for venting time(dynamic simulation model)for a gas trunk pipeline segment and the proposed opening control principle of venting valve provide reference for the optimal operation of gas pipeline venting systems.
文摘In the practical problems such as nuclear waste pollution and seawater intrusion etc., many problems are reduced to solving the convection-diffusion equation, so the research of convection-diffusion equation is of great value. In this work, a spectral method is presented for solving one and two dimensional convection-diffusion equation with source term. The finite difference method is also used to solve the convection diffusion equation. The numerical experiments show that the spectral method is more efficient than other methods for solving the convection-diffusion equation.
基金Sharif University of Technology for financial support of this research
文摘A numerical treatment for the prediction of cavitating flows is presented and assessed. The algorithm uses the preconditioned multiphase Euler equations with appropriate mass transfer terms. A central difference finite volume scheme with suitable dissipation terms to account for density jumps across the cavity interface is shown to yield an effective method for solving the multiphase Euler equations. The Euler equations are utilized herein for the cavitation modeling, because some certain characteristics of cavitating flows can be obtained using the solution of this system of equations with relative low computational effort. In addition, the Euler equations are appropriate for the assessment of the numerical method used, because of the sensitivity of the solution to the numerical instabilities. For this reason, a sensitivity study is conducted to evaluate the effects of various parameters, such as numerical dissipation coefficients and grid size, on the accuracy and performance of the solution. The computations are performed for steady cavitating flows around the NACA 0012 and NACA 66 （MOD） hydrofoils and also an axisymmetric hemispherical fore-body under different conditions and the results are compared with the available numerical and experimental data. The solution procedure presented is shown to be accurate and efficient for predicting steady sheet- and super-cavitation for 2D/axisymmetric geometries.
