In order to implement the optimal design of the indoor thermal comfort based on the numerical modeling method, the numerical calculation platform is combined seamlessly with the data-processing platform, and an intera...In order to implement the optimal design of the indoor thermal comfort based on the numerical modeling method, the numerical calculation platform is combined seamlessly with the data-processing platform, and an interactive numerical calculation platform which includes the functions of numerical simulation and optimization is established. The artificial neural network （ANN） and the greedy strategy are introduced into the hill-climbing pattern heuristic search process, and the optimizing search direction can be predicted by using small samples; when searching along the direction using the greedy strategy, the optimal values can be quickly approached. Therefore, excessive external calling of the numerical modeling process can be avoided, and the optimization time is decreased obviously. The experimental results indicate that the satisfied output parameters of air conditioning can be quickly given out based on the interactive numerical calculation platform and the improved search method, and the optimization for indoor thermal comfort can be completed.展开更多
Numerical methods commonly used for the calculating phase diagrams were listed. A new method to calculate stable phase diagrams was presented which possesse the advantages of both the New-ton Raphson method and the si...Numerical methods commonly used for the calculating phase diagrams were listed. A new method to calculate stable phase diagrams was presented which possesse the advantages of both the New-ton Raphson method and the simplex method. This method is suitable for the calculation of the stable equilibria in complicated systems. For example. calculated results in comparison with experimental results as well as the prediction of new systems were shown in the present work展开更多
In wall-bounded turbulent flow calculations, the past focus has been directed to the modelling of the Reynolds-stress gradients. Not much attention has been paid to the effects of the numerical methods used to calcula...In wall-bounded turbulent flow calculations, the past focus has been directed to the modelling of the Reynolds-stress gradients. Not much attention has been paid to the effects of the numerical methods used to calculate these terms and the modelled equations. Discrepancies between model calculations and measurements are quite often attributed to incorrect modelling, while the suitability and accuracy of the numerical methods used are seldom scrutinized. Instead, alternate near-wall and Reynolds-stress models are proposed to remedy the incorrect turbulent flow calculations. On the other hand, if care is not taken in the numerical treatment of the Reynolds-stress gradient terms, physically unrealistic results and solution instability could occur. Previous studies by the author and his collaborators on the effects of numerical methods have shown that some of the more commonly used numerical methods could enhance numerical stability in the solution procedure but would introduce considerable inaccuracy to the results. The flow cases chosen to demonstrate these inaccuracies are a backstep flow and flow in a square duct, where flow complexities are present. The current investigation attempts to show that the above-mentioned effects of numerical methods could also occur in the calculation of a developing plane channel flow, where flow complexities are absent. In addition, this study shows that the results thus obtained lead to a predicted skin friction coefficient that is influenced more by the numerical method used than by the turbulence model invoked. Together, these results show that numerical treatment of the Reynolds-stress gradients in the equations play an important role, even for a developing plane channel flow.展开更多
The model established in this paper for calculating the unsteady temperature field, in which physical parameters varies with temperatures, is simplified as compared with the classical one by defining the heat conducti...The model established in this paper for calculating the unsteady temperature field, in which physical parameters varies with temperatures, is simplified as compared with the classical one by defining the heat conductivity as function of temperature and dealing with the latent heat of phase transformation and boundary conditions. The results show that the probability of absolute error less 2℃ between the calculated and measured values in temperature field calculation reaches above 80%.展开更多
The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise w...The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise when the algorithm is used for non-homogeneous dynamic systems due to the inverse matrix calculation required. In this paper, the structural dynamic equalibrium equations are converted into a special form, the inverse matrix calculation is replaced by the Crout decomposition method to solve the dynamic equilibrium equations, and the precise integration method without the inverse matrix calculation is obtained. The new algorithm enhances the present precise integration method by improving both the computational accuracy and efficiency. Two numerical examples are given to demonstrate the validity and efficiency of the proposed algorithm.展开更多
The Eulerian？Lagrangian method（ELM） has been used by many ocean models as the solution of the advection equation,but the numerical error caused by interpolation imposes restriction on its accuracy.In the present st...The Eulerian？Lagrangian method（ELM） has been used by many ocean models as the solution of the advection equation,but the numerical error caused by interpolation imposes restriction on its accuracy.In the present study,hybrid N-order Lagrangian interpolation ELM（Li ELM） is put forward in which the N-order Lagrangian interpolation is used at first,then the lower order Lagrangian interpolation is applied in the points where the interpolation results are abnormally higher or lower.The calculation results of a step-shaped salinity advection model are analyzed,which show that higher order（N=3？8） Li ELM can reduce the mean numerical error of salinity calculation,but the numerical oscillation error is still significant.Even number order Li ELM makes larger numerical oscillation error than its adjacent odd number order Li ELM.Hybrid N-order Li ELM can remove numerical oscillation,and it significantly reduces the mean numerical error when N is even and the current is in fixed direction,while it makes less effect on mean numerical error when N is odd or the current direction changes periodically.Hybrid odd number order Li ELM makes less mean numerical error than its adjacent even number order Li ELM when the current is in the fixed direction,while the mean numerical error decreases as N increases when the current direction changes periodically,so odd number of N may be better for application.Among various types of Hybrid N-order Li ELM,the scheme reducing N-order directly to 1st-order may be the optimal for synthetic selection of accuracy and computational efficiency.展开更多
The tube inside and outside heat transfer mechanism of Passive Residual Heat Removal Heat Exchanger (PRHR HX) was analyzed. The calculation method of this special heat exchanger under natural convection condition in I...The tube inside and outside heat transfer mechanism of Passive Residual Heat Removal Heat Exchanger (PRHR HX) was analyzed. The calculation method of this special heat exchanger under natural convection condition in In-containment Refueling Water Storage Tank (IRWST) was carried out. The single-tube coupling model three-dimensional natural circulation in the IRWST was simulated numerically using Fluent. The heat transfer and flow characteristics of the fluid in IRWST were obtained. The comparison of the results between theoretical arithmetic and numerical simulation showed that the theoretical calculation method is suitable for the heat transfer calculation of PRHR HX.展开更多
This article reviews the application and progress of deep learning in efficient numerical computing methods.Deep learning,as an important branch of machine learning,provides new ideas for numerical computation by cons...This article reviews the application and progress of deep learning in efficient numerical computing methods.Deep learning,as an important branch of machine learning,provides new ideas for numerical computation by constructing multi-layer neural networks to simulate the learning process of the human brain.The article explores the application of deep learning in solving partial differential equations,optimizing problems,and data-driven modeling,and analyzes its advantages in computational efficiency,accuracy,and adaptability.At the same time,this article also points out the challenges faced by deep learning numerical computation methods in terms of computational efficiency,interpretability,and generalization ability,and proposes strategies and future development directions for integrating with traditional numerical methods.展开更多
The method of two-dimensional viscous space-time conservation element and solution element (CE/SE) can be used to calculate the gas-liquid two-phase interior flow field in pulse detonation engine (PDE). In this paper,...The method of two-dimensional viscous space-time conservation element and solution element (CE/SE) can be used to calculate the gas-liquid two-phase interior flow field in pulse detonation engine (PDE). In this paper, the evolution of the detonation wave and the distribution of its physical parameters were analyzed. The numerical results show that the change of axial velocity of gas is the same as that of detonation pressure. The larger the liquid droplet radius is, the longer the time to get stable detonation wave is. The calculated results coincide with the experimented results better.展开更多
Widespread usage of bump-type foil journal bearing(BFJB) in oil-free microturbomachinery requires accurate predictions of dynamic performance characteristics, although it remains a challenging issue because BFJB refle...Widespread usage of bump-type foil journal bearing(BFJB) in oil-free microturbomachinery requires accurate predictions of dynamic performance characteristics, although it remains a challenging issue because BFJB reflects nonlinear both structurally and aerodynamically.This paper presented a simple experimental method to semi-directly obtain the minimum film thickness and dynamic stiffness of BFJB using the journal orbit. Numerical calculations and simulations are conducted to validate the experimental method. The micro-deformation and interaction of various foils are taken into consideration to improve the model precision. The results from the numerical model regarding the BFJB dynamic characteristics are compared with the experimental results coming from a dedicated test rig, which shows that the experimental results fluctuate obviously and agree not well with the numerical results at the start stage due to the presence of dry friction at that time, nevertheless, they show fantastic agreement as soon as a gas film is gradually generated to separate the shaft from the top foil. Therefore, the proposed experimental method is effective to predict film thickness and dynamic characteristics during the period from the lift-off time to the land-off time. The dynamic characteristics, along with the journal orbits also can be used to rapidly predict the dynamics behavior of rotor-bearing systems.展开更多
As a further development of the authors' work (Huang and Qian, 1993), in this paper a new numerical method based on the time domain boundary element technique is proposed for solving fluid-solid coupling problems,...As a further development of the authors' work (Huang and Qian, 1993), in this paper a new numerical method based on the time domain boundary element technique is proposed for solving fluid-solid coupling problems, in which a rigid body impacts normally on the calm surface of a half-space fluid. A fundamental solution to the half-space potential flow problem is first derived with the method of images. Then, an equivalent boundary integral equation in the Laplace transform domain is established by means of Green's second identity. Through the inverse Laplace transform and discretization in both time and boundary of the fluid region, the numerical calculation for the problem under consideration has been carried out. Several examples demonstrate that the present method is more efficient than existing ones, from which it is also seen that the shape of the impacting body has a considerable effect on the total impact force.展开更多
Based on the study on electromagnetic field migration by Zhdanov, we have proposed an improved method for the weak points in the research. Firstly, the initial background resistivity should be determined by using 1-D ...Based on the study on electromagnetic field migration by Zhdanov, we have proposed an improved method for the weak points in the research. Firstly, the initial background resistivity should be determined by using 1-D inversion results. Then in the process of continuation, the results are corrected and calculated layer by layer by the iteration method, so that more exact resistivity can be obtained. Secondly, an improved algorithm for finite-difference equation is studied. According to the property of electromagnetic migration field, the algorithm is designed by means of grids varying with geometric progression in the longitudinal direction. Being improved by the techniques mentioned above, better results are obtained by the new method, which has been verified by both the theory model and practical data.展开更多
In this paper two classes of equivalence transform methods for solving ordinary differential equations are proposed. One class of method is the equivalence integral transform method for special differential algebraic ...In this paper two classes of equivalence transform methods for solving ordinary differential equations are proposed. One class of method is the equivalence integral transform method for special differential algebraic problems. The advantage of this class of method is such that the amount of work calculating one integration with parameters becomes that of two interpolations, when the system of nonlinear equations is solved on the right hand side function. The other class of method is the equivalence substitution method for avoiding calculating derivative on the right hand side function. In order to avoid calculation derivatives, two equivalence substitution methods are proposed here. The application instances of some special effect of the equivalence substitution methods are given.展开更多
Supernova remnant（SNR） G309.2–0.6 has a peculiar radio morphology with two bright ears to the southwest and northeast, although the main shell outside the ears is roughly circular. Based on an earlier proposal that...Supernova remnant（SNR） G309.2–0.6 has a peculiar radio morphology with two bright ears to the southwest and northeast, although the main shell outside the ears is roughly circular. Based on an earlier proposal that the supernova ejecta has a jet component with extra energy, the dynamical evolution of the remnant is solved using 3 D hydrodynamical（HD） simulation to investigate the formation of the periphery of the remnant. Assuming the ejecta with a kinetic energy of 10^51 erg and a mass of 3 M⊙evolved in a uniform ambient medium for a time of-4000 yr and the jet component has cylindrical symmetry with a half open angle of 10°, the result indicates that the energy contained in the jet is about10%–15% of the kinetic energy of the entire ejecta to reproduce the detected profile. This study supports that the remnant originated from a jet-driven core-collapse supernova.展开更多
In order to achieve highly accurate and efficient numerical calculations of structural dynamics, time collocation method is presented. For a given time interval, the numerical solution of the method is approximated by...In order to achieve highly accurate and efficient numerical calculations of structural dynamics, time collocation method is presented. For a given time interval, the numerical solution of the method is approximated by a polynomial. The polynomial coefficients are evaluated by solving algebraic equation. Once the polynomial coefficients are evaluated, the numerical solutions at any time in the interval can be easily calculated. New formulae are derived for the polynomial coefficients,which are more practical and succinct than those previously given. Two structural dynamic equations are calculated by the proposed method. The numerical solutions are compared with the traditional fourth-order Runge-Kutta method. The results show that the method proposed is highly accurate and computationally efficient. In addition, an important advantage of the method is the simplicity in software programming.展开更多
Based on the mathematical model of one dimension transient flow of the polymer foam in porous media, the numerical calculation method of the flow mentioned above by using the finite difference method is given. Through...Based on the mathematical model of one dimension transient flow of the polymer foam in porous media, the numerical calculation method of the flow mentioned above by using the finite difference method is given. Through the experiments of one dimension transient flow of HPAM (Hydrolytic Polyacrylamide) foam in the artificial sandstone core, the HPAM foam generation and coalescence coefficient of the mathematical model mentioned above are determined. The profiles of the liquid phase saturation, the pressure drop and the number density of one dimension transient flow of HPAM foam with the dimensionless time in artificial sandstone core are numerically calculated and analyzed by using the numerical calculation method.展开更多
We investigate the forming of gamma-ray burst pulses with a simple onedimensional relativistic shock model. The mechanism is that a ＂central engine＂ drives forward the nearby plasma inside the fireball to generate a...We investigate the forming of gamma-ray burst pulses with a simple onedimensional relativistic shock model. The mechanism is that a ＂central engine＂ drives forward the nearby plasma inside the fireball to generate a series of pressure waves. We give a relativistic geometric recurrence formula that connects the time when the pressure waves are produced and the time when the corresponding shocks occurred. This relation enables us to relate the pulse magnitude with the observation time. Our analysis shows that the evolution of the pressure waves leads to a fast rise and an exponential decay pulses. In determining the width of the pulses, the acceleration time is more important than that of the deceleration.展开更多
In the field of discretization-based meshfree/meshless methods,the improvements in the higher-order consistency,stability,and computational efficiency are of great concerns in computational science and numerical solut...In the field of discretization-based meshfree/meshless methods,the improvements in the higher-order consistency,stability,and computational efficiency are of great concerns in computational science and numerical solutions to partial differential equations.Various alternative numerical methods of the finite particle method(FPM)frame have been extended from mathematical theories to numerical applications separately.As a comprehensive numerical scheme,this study suggests a unified resolved program for numerically investigating their accuracy,stability,consistency,computational efficiency,and practical applicability in industrial engineering contexts.The high-order finite particle method(HFPM)and corrected methods based on the multivariate Taylor series expansion are constructed and analyzed to investigate the whole applicability in different benchmarks of computational fluid dynamics.Specifically,four benchmarks are designed purposefully from statical exact solutions to multifaceted hydrodynamic tests,which possess different numerical performances on the particle consistency,numerical discretized forms,particle distributions,and transient time evolutional stabilities.This study offers a numerical reference for the current unified resolved program.展开更多
基金Sponsored by the National Program"973"Project (2005CB623906)
文摘In order to implement the optimal design of the indoor thermal comfort based on the numerical modeling method, the numerical calculation platform is combined seamlessly with the data-processing platform, and an interactive numerical calculation platform which includes the functions of numerical simulation and optimization is established. The artificial neural network （ANN） and the greedy strategy are introduced into the hill-climbing pattern heuristic search process, and the optimizing search direction can be predicted by using small samples; when searching along the direction using the greedy strategy, the optimal values can be quickly approached. Therefore, excessive external calling of the numerical modeling process can be avoided, and the optimization time is decreased obviously. The experimental results indicate that the satisfied output parameters of air conditioning can be quickly given out based on the interactive numerical calculation platform and the improved search method, and the optimization for indoor thermal comfort can be completed.
文摘Numerical methods commonly used for the calculating phase diagrams were listed. A new method to calculate stable phase diagrams was presented which possesse the advantages of both the New-ton Raphson method and the simplex method. This method is suitable for the calculation of the stable equilibria in complicated systems. For example. calculated results in comparison with experimental results as well as the prediction of new systems were shown in the present work
文摘In wall-bounded turbulent flow calculations, the past focus has been directed to the modelling of the Reynolds-stress gradients. Not much attention has been paid to the effects of the numerical methods used to calculate these terms and the modelled equations. Discrepancies between model calculations and measurements are quite often attributed to incorrect modelling, while the suitability and accuracy of the numerical methods used are seldom scrutinized. Instead, alternate near-wall and Reynolds-stress models are proposed to remedy the incorrect turbulent flow calculations. On the other hand, if care is not taken in the numerical treatment of the Reynolds-stress gradient terms, physically unrealistic results and solution instability could occur. Previous studies by the author and his collaborators on the effects of numerical methods have shown that some of the more commonly used numerical methods could enhance numerical stability in the solution procedure but would introduce considerable inaccuracy to the results. The flow cases chosen to demonstrate these inaccuracies are a backstep flow and flow in a square duct, where flow complexities are present. The current investigation attempts to show that the above-mentioned effects of numerical methods could also occur in the calculation of a developing plane channel flow, where flow complexities are absent. In addition, this study shows that the results thus obtained lead to a predicted skin friction coefficient that is influenced more by the numerical method used than by the turbulence model invoked. Together, these results show that numerical treatment of the Reynolds-stress gradients in the equations play an important role, even for a developing plane channel flow.
文摘The model established in this paper for calculating the unsteady temperature field, in which physical parameters varies with temperatures, is simplified as compared with the classical one by defining the heat conductivity as function of temperature and dealing with the latent heat of phase transformation and boundary conditions. The results show that the probability of absolute error less 2℃ between the calculated and measured values in temperature field calculation reaches above 80%.
文摘The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise when the algorithm is used for non-homogeneous dynamic systems due to the inverse matrix calculation required. In this paper, the structural dynamic equalibrium equations are converted into a special form, the inverse matrix calculation is replaced by the Crout decomposition method to solve the dynamic equilibrium equations, and the precise integration method without the inverse matrix calculation is obtained. The new algorithm enhances the present precise integration method by improving both the computational accuracy and efficiency. Two numerical examples are given to demonstrate the validity and efficiency of the proposed algorithm.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.40906044 and 41076048)the Fundamental Research Funds for the Central Universities Project(Grant No.2011B05714)
文摘The Eulerian？Lagrangian method（ELM） has been used by many ocean models as the solution of the advection equation,but the numerical error caused by interpolation imposes restriction on its accuracy.In the present study,hybrid N-order Lagrangian interpolation ELM（Li ELM） is put forward in which the N-order Lagrangian interpolation is used at first,then the lower order Lagrangian interpolation is applied in the points where the interpolation results are abnormally higher or lower.The calculation results of a step-shaped salinity advection model are analyzed,which show that higher order（N=3？8） Li ELM can reduce the mean numerical error of salinity calculation,but the numerical oscillation error is still significant.Even number order Li ELM makes larger numerical oscillation error than its adjacent odd number order Li ELM.Hybrid N-order Li ELM can remove numerical oscillation,and it significantly reduces the mean numerical error when N is even and the current is in fixed direction,while it makes less effect on mean numerical error when N is odd or the current direction changes periodically.Hybrid odd number order Li ELM makes less mean numerical error than its adjacent even number order Li ELM when the current is in the fixed direction,while the mean numerical error decreases as N increases when the current direction changes periodically,so odd number of N may be better for application.Among various types of Hybrid N-order Li ELM,the scheme reducing N-order directly to 1st-order may be the optimal for synthetic selection of accuracy and computational efficiency.
文摘The tube inside and outside heat transfer mechanism of Passive Residual Heat Removal Heat Exchanger (PRHR HX) was analyzed. The calculation method of this special heat exchanger under natural convection condition in In-containment Refueling Water Storage Tank (IRWST) was carried out. The single-tube coupling model three-dimensional natural circulation in the IRWST was simulated numerically using Fluent. The heat transfer and flow characteristics of the fluid in IRWST were obtained. The comparison of the results between theoretical arithmetic and numerical simulation showed that the theoretical calculation method is suitable for the heat transfer calculation of PRHR HX.
文摘This article reviews the application and progress of deep learning in efficient numerical computing methods.Deep learning,as an important branch of machine learning,provides new ideas for numerical computation by constructing multi-layer neural networks to simulate the learning process of the human brain.The article explores the application of deep learning in solving partial differential equations,optimizing problems,and data-driven modeling,and analyzes its advantages in computational efficiency,accuracy,and adaptability.At the same time,this article also points out the challenges faced by deep learning numerical computation methods in terms of computational efficiency,interpretability,and generalization ability,and proposes strategies and future development directions for integrating with traditional numerical methods.
基金Sponsored by the National Natural Science Foundation of China (10672080)
文摘The method of two-dimensional viscous space-time conservation element and solution element (CE/SE) can be used to calculate the gas-liquid two-phase interior flow field in pulse detonation engine (PDE). In this paper, the evolution of the detonation wave and the distribution of its physical parameters were analyzed. The numerical results show that the change of axial velocity of gas is the same as that of detonation pressure. The larger the liquid droplet radius is, the longer the time to get stable detonation wave is. The calculated results coincide with the experimented results better.
基金funded by the Natural Science Foundation of China (No. 51775025, 51205015)China Key Research and Development Plan (No. 2017YFB0102102, 2018YFB0104100)。
文摘Widespread usage of bump-type foil journal bearing(BFJB) in oil-free microturbomachinery requires accurate predictions of dynamic performance characteristics, although it remains a challenging issue because BFJB reflects nonlinear both structurally and aerodynamically.This paper presented a simple experimental method to semi-directly obtain the minimum film thickness and dynamic stiffness of BFJB using the journal orbit. Numerical calculations and simulations are conducted to validate the experimental method. The micro-deformation and interaction of various foils are taken into consideration to improve the model precision. The results from the numerical model regarding the BFJB dynamic characteristics are compared with the experimental results coming from a dedicated test rig, which shows that the experimental results fluctuate obviously and agree not well with the numerical results at the start stage due to the presence of dry friction at that time, nevertheless, they show fantastic agreement as soon as a gas film is gradually generated to separate the shaft from the top foil. Therefore, the proposed experimental method is effective to predict film thickness and dynamic characteristics during the period from the lift-off time to the land-off time. The dynamic characteristics, along with the journal orbits also can be used to rapidly predict the dynamics behavior of rotor-bearing systems.
基金National Defence Foundation of Science & Technology of China,Grant No.J14.8.1JW05National Key Project of China(No.PD9521907)
文摘As a further development of the authors' work (Huang and Qian, 1993), in this paper a new numerical method based on the time domain boundary element technique is proposed for solving fluid-solid coupling problems, in which a rigid body impacts normally on the calm surface of a half-space fluid. A fundamental solution to the half-space potential flow problem is first derived with the method of images. Then, an equivalent boundary integral equation in the Laplace transform domain is established by means of Green's second identity. Through the inverse Laplace transform and discretization in both time and boundary of the fluid region, the numerical calculation for the problem under consideration has been carried out. Several examples demonstrate that the present method is more efficient than existing ones, from which it is also seen that the shape of the impacting body has a considerable effect on the total impact force.
文摘Based on the study on electromagnetic field migration by Zhdanov, we have proposed an improved method for the weak points in the research. Firstly, the initial background resistivity should be determined by using 1-D inversion results. Then in the process of continuation, the results are corrected and calculated layer by layer by the iteration method, so that more exact resistivity can be obtained. Secondly, an improved algorithm for finite-difference equation is studied. According to the property of electromagnetic migration field, the algorithm is designed by means of grids varying with geometric progression in the longitudinal direction. Being improved by the techniques mentioned above, better results are obtained by the new method, which has been verified by both the theory model and practical data.
基金The project was supported by the National Natural Science Faundation of China
文摘In this paper two classes of equivalence transform methods for solving ordinary differential equations are proposed. One class of method is the equivalence integral transform method for special differential algebraic problems. The advantage of this class of method is such that the amount of work calculating one integration with parameters becomes that of two interpolations, when the system of nonlinear equations is solved on the right hand side function. The other class of method is the equivalence substitution method for avoiding calculating derivative on the right hand side function. In order to avoid calculation derivatives, two equivalence substitution methods are proposed here. The application instances of some special effect of the equivalence substitution methods are given.
基金supported by the Yunnan Applied Basic Research Project (2016FD105)the Foundations of Yunnan Province (2016ZZX180 and 2016DG006)+4 种基金Kunming University (YJL15004 and XJL15015)supported by the National Natural Science Foundation of China (No. 11563009)the Yunnan Applied Basic Research Project (2016FB001)the Candidate Talents Training Fund of Yunnan Province (2017HB003)the Program for Excellent Young Talents, Yunnan University (WX069051 and 2017YDYQ01)
文摘Supernova remnant（SNR） G309.2–0.6 has a peculiar radio morphology with two bright ears to the southwest and northeast, although the main shell outside the ears is roughly circular. Based on an earlier proposal that the supernova ejecta has a jet component with extra energy, the dynamical evolution of the remnant is solved using 3 D hydrodynamical（HD） simulation to investigate the formation of the periphery of the remnant. Assuming the ejecta with a kinetic energy of 10^51 erg and a mass of 3 M⊙evolved in a uniform ambient medium for a time of-4000 yr and the jet component has cylindrical symmetry with a half open angle of 10°, the result indicates that the energy contained in the jet is about10%–15% of the kinetic energy of the entire ejecta to reproduce the detected profile. This study supports that the remnant originated from a jet-driven core-collapse supernova.
基金Supported by Liu Hui Applied Mathematics Center of Nankai University-Tianjin University( No. H10124).
文摘In order to achieve highly accurate and efficient numerical calculations of structural dynamics, time collocation method is presented. For a given time interval, the numerical solution of the method is approximated by a polynomial. The polynomial coefficients are evaluated by solving algebraic equation. Once the polynomial coefficients are evaluated, the numerical solutions at any time in the interval can be easily calculated. New formulae are derived for the polynomial coefficients,which are more practical and succinct than those previously given. Two structural dynamic equations are calculated by the proposed method. The numerical solutions are compared with the traditional fourth-order Runge-Kutta method. The results show that the method proposed is highly accurate and computationally efficient. In addition, an important advantage of the method is the simplicity in software programming.
文摘Based on the mathematical model of one dimension transient flow of the polymer foam in porous media, the numerical calculation method of the flow mentioned above by using the finite difference method is given. Through the experiments of one dimension transient flow of HPAM (Hydrolytic Polyacrylamide) foam in the artificial sandstone core, the HPAM foam generation and coalescence coefficient of the mathematical model mentioned above are determined. The profiles of the liquid phase saturation, the pressure drop and the number density of one dimension transient flow of HPAM foam with the dimensionless time in artificial sandstone core are numerically calculated and analyzed by using the numerical calculation method.
基金Supported by the National Natural Science Foundation of China.
文摘We investigate the forming of gamma-ray burst pulses with a simple onedimensional relativistic shock model. The mechanism is that a ＂central engine＂ drives forward the nearby plasma inside the fireball to generate a series of pressure waves. We give a relativistic geometric recurrence formula that connects the time when the pressure waves are produced and the time when the corresponding shocks occurred. This relation enables us to relate the pulse magnitude with the observation time. Our analysis shows that the evolution of the pressure waves leads to a fast rise and an exponential decay pulses. In determining the width of the pulses, the acceleration time is more important than that of the deceleration.
基金supported by the National Natural Science Foundation of China(No.12002290)。
文摘In the field of discretization-based meshfree/meshless methods,the improvements in the higher-order consistency,stability,and computational efficiency are of great concerns in computational science and numerical solutions to partial differential equations.Various alternative numerical methods of the finite particle method(FPM)frame have been extended from mathematical theories to numerical applications separately.As a comprehensive numerical scheme,this study suggests a unified resolved program for numerically investigating their accuracy,stability,consistency,computational efficiency,and practical applicability in industrial engineering contexts.The high-order finite particle method(HFPM)and corrected methods based on the multivariate Taylor series expansion are constructed and analyzed to investigate the whole applicability in different benchmarks of computational fluid dynamics.Specifically,four benchmarks are designed purposefully from statical exact solutions to multifaceted hydrodynamic tests,which possess different numerical performances on the particle consistency,numerical discretized forms,particle distributions,and transient time evolutional stabilities.This study offers a numerical reference for the current unified resolved program.
