To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based sim...To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.展开更多
In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step ...In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.展开更多
Outwash deposit is a unique type of geological materials, and its features such as heterogeneity, discontinuity and nonlinearity determine the complexity of mechanical characteristics and failure mechanism. In this wo...Outwash deposit is a unique type of geological materials, and its features such as heterogeneity, discontinuity and nonlinearity determine the complexity of mechanical characteristics and failure mechanism. In this work, random meso-structure of outwash deposits was constructed by the technique of computer random simulation based on characteristics of its meso-structure in the statistical sense and some simplifications, and a series of large direct shear tests on numerical samples of outwash deposits with stone contents of 15%, 30%, 45% and 60% were conducted using the discrete element method to further investigate its mechanical characteristics and failure mechanism under external load. The results show that the deformation characteristics and shear strength of outwash deposits are to some extent improved with the increase of stone content, and the shear stress–shear displacement curves of outwash deposits show great differences at the post-peak stage due to the random spatial distribution and content of stones. From the mesoscopic view, normal directions of contacts between "soil" and "stone" particles undergo apparent deflection as the shear displacement continues during the shearing process, accompanying redistribution of the magnitude of contact forces during the shearing process. For outwash deposits, the shear zone formed after shear failure is an irregular stripe due to the movements of stones near the shear zone, and it expands gradually with the increase of stone content. In addition, there is an approximately linear relation between the mean increment of internal friction angle and the stone content lying between 30% and 60%, and a concave nonlinear relation between the mean increment of cohesion and stone content, which are in good agreement with the existing research results.展开更多
Fetr6 is an underground mine using the stope-and-pillar mining method. As there was some evidence regarding pillar failure in this mine, improving works such as roof support and replacing existing pillars with concret...Fetr6 is an underground mine using the stope-and-pillar mining method. As there was some evidence regarding pillar failure in this mine, improving works such as roof support and replacing existing pillars with concrete pillars (CP) were carried out. During the construction of the second CP, in the space between the remaining pillars, one of the pillars failed leading to the progressive failure of other pillars until 4 000 m 2 of mine had collapsed within a few minutes. In this work, this phenomenon is described by applying both numerical and empirical methods and the respective results are compared. The results of numerical modelling are found to be closer to the actual condition than those of the empirical method. Also, a width-to-height (W/H) ratio less than 1, an inadequate support system and the absence of a detailed program for pillar recovery are shown to be the most important causes of the Domino failure in this mine.展开更多
The behavior of single bubble rising in quiescent shear-thinning tlmds was lnvestlgateO numerically by level set metnoa. number of bubbles in a large range of Reynolds number and Eotvos number were investigated includ...The behavior of single bubble rising in quiescent shear-thinning tlmds was lnvestlgateO numerically by level set metnoa. number of bubbles in a large range of Reynolds number and Eotvos number were investigated including spherical, oblate and spherical. The bubble shape and drag coefficient were compared with experimental results. It is observed that the simulated results show good conformity to experimental results over a wide range of Reynolds number. In addition, the detailed flow field based on the reference coordinate system moving with the bubble is obtained, and the relationship among flow field, bubble shape and velocity is discussed.展开更多
Columnar jointed rock mass with unique geometric and geological properties is one spectacular example of geometrical order in nature.Columnar joints are generally accepted to be formed by spatially uniform volume cont...Columnar jointed rock mass with unique geometric and geological properties is one spectacular example of geometrical order in nature.Columnar joints are generally accepted to be formed by spatially uniform volume contraction during cooling.In this paper,substantial field work was performed to study the geological characteristics of irregular columnar jointed basalt on the left bank dam foundation in the Baihetan Hydropower Station,where the columnar jointed rock mass is extensively exposed due to excavation.The quantitative measurements of the sizing of polygonal crack pattern of columnar joints and assessment of their degree of irregularity were summarized.Considering the irregularity of polygonal crack pattern,a modified Voronoi polygon(MVP)method was developed to model the special polygonal crack pattern of columnar joints.The new polygonal pattern obtained by the MVP method consists of a large number of irregular polygons,of which the degree of irregularity is consistent with the field measurement results.This method can reproduce the rapid evolution from an initial ideal regular hexagonal pattern to a final actual irregular polygonal pattern as the degree of irregularity increases.The compression tests of columnar jointed rock mass with different irregularity show that the geometric irregularity has a great influence on its mechanical properties.This numerical construction method provides a reliable way to reconstruct columnar joint structure with specific polygonal crack pattern,which is consistent with onsite columnar jointed basalt.展开更多
Several parameter identification methods of thermal response test were evaluated through numerical and experimental study.A three-dimensional finite-volume numerical model was established under the assumption that the...Several parameter identification methods of thermal response test were evaluated through numerical and experimental study.A three-dimensional finite-volume numerical model was established under the assumption that the soil thermal conductivity had been known in the simulation of thermal response test.The thermal response curve was firstly obtained through numerical calculation.Then,the accuracy of the numerical model was verified with measured data obtained through a thermal response test.Based on the numerical and experimental thermal response curves,the thermal conductivity of the soil was calculated by different parameter identification methods.The calculated results were compared with the assumed value and then the accuracy of these methods was evaluated.Furthermore,the effects of test time,variable data quality,borehole radius,initial ground temperature,and heat injection rate were analyzed.The results show that the method based on cylinder-source model has a low precision and the identified thermal conductivity decreases with an increase in borehole radius.For parameter estimation,the measuring accuracy of the initial temperature of the deep ground soil has greater effect on identified thermal conductivity.展开更多
Partition of unity based numerical manifold method can solve continuous and discontinuous problems in a unified framework with a two-cover system,i.e.,the mathematical cover and physical cover.However,renewal of the t...Partition of unity based numerical manifold method can solve continuous and discontinuous problems in a unified framework with a two-cover system,i.e.,the mathematical cover and physical cover.However,renewal of the topology of the two-cover system poses a challenge for multiple crack propagation problems and there are few references.In this study,a robust and efficient strategy is proposed to update the cover system of the numerical manifold method in simulation of multiple crack propagation problems.The proposed algorithm updates the cover system with a bottom-up process:1)identification of fractured manifold elements according to the previous and latest crack tip position;and 2)local topological update of the manifold elements,physical patches,block boundary loops,and non-persistent joint loops according to the scenario classification of the propagating crack.The proposed crack tracking strategy and classification of the renewal cases promote a robust and efficient cover renewal algorithm for multiple crack propagation analysis.Three crack propagation examples show that the proposed algorithm performs well in updating the cover system.This cover renewal methodology can be extended for numerical manifold method with polygonal mathematical covers.展开更多
Combined with the kinetic model of liquid film spreading, a new numerical method of solid-liquid-gas three-phase flow was developed for the moving of contact line, which was a hybrid method of computational fluid dyna...Combined with the kinetic model of liquid film spreading, a new numerical method of solid-liquid-gas three-phase flow was developed for the moving of contact line, which was a hybrid method of computational fluid dynamics and lattice Boltzmalm method （LBM）. By taking the effect of molecule force in droplet and the wall surface on liquid film into account, the changing law of contact angle with different surface tensions was analyzed on glass and aluminum foil surfaces. Compared with experimental results, the standard deviation by using LBM is less than 0.5°, which validates the feasibility of LBM simulation on the dynamic process of liquid film spreading. In addition, oscillations are discovered both at the initial and end phases. The phenomenon of retraction is also found and the maximum retraction angle is 7.58°. The obtained result shows that the retraction is proved to be correlative with precursor film by tracking the volume change of liquid film contour. Furthermore, non-dimensional coefficient 2 is introduced to measure the liquid film retraction capacity.展开更多
It is well known that the incorrect results will be given using either the electric or magnetic field integral equation to calculate the radar cross section (RCS) of a closed body at the interior resonance. In this pa...It is well known that the incorrect results will be given using either the electric or magnetic field integral equation to calculate the radar cross section (RCS) of a closed body at the interior resonance. In this paper, an effective iterative technique is used to correct the calculated surface current density from the electric field integral equation. The radar cross section is computed for an infinite conducting circular cylinder at the interior resonance, and the obtained results are in good agreement with the analytical results. The backscattering cross section of an infinite triangular cylinder in the vicinity of a resonant frequency is also calculated. It is shown that the presence method is efficient and accurate.展开更多
A kind of international rapid field measurement methods of hydraulic conductivity and it's applications in Sanjiang Plain have been introduced in the paper.
A class of modified parallel combined methods of real-time numerical simulation are presented for a stiff dynamic system. By combining the parallelism across the system with the parallelism across the method, and rela...A class of modified parallel combined methods of real-time numerical simulation are presented for a stiff dynamic system. By combining the parallelism across the system with the parallelism across the method, and relaxing the dependence of stage value computation on sampling time of input function, a class of modified real-time parallel combined methods are constructed. Stiff and nonstiff subsystems are solved in parallel on a parallel computer by a parallel Rosen-brock method and a parallel RK method, respectively. Their order conditions and convergences are discussed. The numerical simulation experiments show that this class of modified algorithms can get high speed and efficiency.展开更多
This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solutio...This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solution of test equations u'(t) = a 11 u(t) + a12v(t) + b11 u(t - τ1) + b12v(t-τ2,v'(t) = a21 u(t) + a22 v(t) + b21 u(t -τ1,) + b22 v(t -τ2), t>0,with initial conditionsu(t)=u0(t),v(t) =v0(t), t≤0.where aij, bij∈C, τj >0, i,j = 1,2,, and u0(t), v0(t)are continuous and complex valued. Sufficient conditions for the asymptotic stability of test equation are derived. Furthermore, with respect to an appropriate definition of stability for the numerical method, it is proved that the linear θ-method is stable if and only if 1/2≤θ≤1 and the one-leg θ-method is stable if and only if θ= 1.展开更多
In this paper, a novel calibration integral equation is derived for resolving double-sided, two-probe inverse heat conduction problem of surface heat flux estimation. In contrast to the conventional inverse heat condu...In this paper, a novel calibration integral equation is derived for resolving double-sided, two-probe inverse heat conduction problem of surface heat flux estimation. In contrast to the conventional inverse heat conduction techniques, this calibration approach does not require explicit input of the probe locations, thermophysical properties of the host material and temperature sensor parameters related to thermal contact resistance, sensor capacitance and conductive lead losses. All those parameters and properties are inherently contained in the calibration framework in terms of Volterra integral equation of the first kind. The Laplace transform technique is applied and the frequency domain manipulations of the heat equation are performed for deriving the calibration integral equation. Due to the ill-posed nature, regularization is required for the inverse heat conduction problem, a future-time method or singular value decomposition (SVD) can be used for stabilizing the ill-posed Volterra integral equation of the first kind.展开更多
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52271317 and 52071149)the Fundamental Research Funds for the Central Universities(HUST:2019kfy XJJS007)。
文摘To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.
基金Project(41630642)supported by the Key Project of National Natural Science Foundation of ChinaProject(51974360)supported by the National Natural Science Foundation of ChinaProject(2018JJ3656)supported by the Natural Science Foundation of Hunan Province,China。
文摘In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.
基金Project(2011CB013504) supported by the National Basic Research Program(973 Program)of ChinaProject(2013BAB06B01) supported by the National Science&Technology Pillar Program during the Twelfth Five-year Plan Period+2 种基金Projects(11772118,51479049,51709282) supported by the National Natural Science Foundation of ChinaProject(2017M620838) supported by the Postdoctoral Science Foundation of ChinaProject(487237) supported by the Natural Sciences and Engineering Research Council of Canada
文摘Outwash deposit is a unique type of geological materials, and its features such as heterogeneity, discontinuity and nonlinearity determine the complexity of mechanical characteristics and failure mechanism. In this work, random meso-structure of outwash deposits was constructed by the technique of computer random simulation based on characteristics of its meso-structure in the statistical sense and some simplifications, and a series of large direct shear tests on numerical samples of outwash deposits with stone contents of 15%, 30%, 45% and 60% were conducted using the discrete element method to further investigate its mechanical characteristics and failure mechanism under external load. The results show that the deformation characteristics and shear strength of outwash deposits are to some extent improved with the increase of stone content, and the shear stress–shear displacement curves of outwash deposits show great differences at the post-peak stage due to the random spatial distribution and content of stones. From the mesoscopic view, normal directions of contacts between "soil" and "stone" particles undergo apparent deflection as the shear displacement continues during the shearing process, accompanying redistribution of the magnitude of contact forces during the shearing process. For outwash deposits, the shear zone formed after shear failure is an irregular stripe due to the movements of stones near the shear zone, and it expands gradually with the increase of stone content. In addition, there is an approximately linear relation between the mean increment of internal friction angle and the stone content lying between 30% and 60%, and a concave nonlinear relation between the mean increment of cohesion and stone content, which are in good agreement with the existing research results.
文摘Fetr6 is an underground mine using the stope-and-pillar mining method. As there was some evidence regarding pillar failure in this mine, improving works such as roof support and replacing existing pillars with concrete pillars (CP) were carried out. During the construction of the second CP, in the space between the remaining pillars, one of the pillars failed leading to the progressive failure of other pillars until 4 000 m 2 of mine had collapsed within a few minutes. In this work, this phenomenon is described by applying both numerical and empirical methods and the respective results are compared. The results of numerical modelling are found to be closer to the actual condition than those of the empirical method. Also, a width-to-height (W/H) ratio less than 1, an inadequate support system and the absence of a detailed program for pillar recovery are shown to be the most important causes of the Domino failure in this mine.
基金Project(21406141)supported by the National Natural Science Foundation of ChinaProject(20141078)supported by the Scientific Research Starting Foundation for Doctors of Liaoning Province,China+1 种基金Project(L2014060)supported by the Foundation of Department of Education of Liaoning Province,ChinaProject(157B21)supported by the Scientific Research Starting Foundation for Doctors of Shenyang Aerospace University,China
文摘The behavior of single bubble rising in quiescent shear-thinning tlmds was lnvestlgateO numerically by level set metnoa. number of bubbles in a large range of Reynolds number and Eotvos number were investigated including spherical, oblate and spherical. The bubble shape and drag coefficient were compared with experimental results. It is observed that the simulated results show good conformity to experimental results over a wide range of Reynolds number. In addition, the detailed flow field based on the reference coordinate system moving with the bubble is obtained, and the relationship among flow field, bubble shape and velocity is discussed.
基金Projects(51621006,51779251)supported by the National Natural Science Foundation of China。
文摘Columnar jointed rock mass with unique geometric and geological properties is one spectacular example of geometrical order in nature.Columnar joints are generally accepted to be formed by spatially uniform volume contraction during cooling.In this paper,substantial field work was performed to study the geological characteristics of irregular columnar jointed basalt on the left bank dam foundation in the Baihetan Hydropower Station,where the columnar jointed rock mass is extensively exposed due to excavation.The quantitative measurements of the sizing of polygonal crack pattern of columnar joints and assessment of their degree of irregularity were summarized.Considering the irregularity of polygonal crack pattern,a modified Voronoi polygon(MVP)method was developed to model the special polygonal crack pattern of columnar joints.The new polygonal pattern obtained by the MVP method consists of a large number of irregular polygons,of which the degree of irregularity is consistent with the field measurement results.This method can reproduce the rapid evolution from an initial ideal regular hexagonal pattern to a final actual irregular polygonal pattern as the degree of irregularity increases.The compression tests of columnar jointed rock mass with different irregularity show that the geometric irregularity has a great influence on its mechanical properties.This numerical construction method provides a reliable way to reconstruct columnar joint structure with specific polygonal crack pattern,which is consistent with onsite columnar jointed basalt.
基金Project(xjj20100078) supported by the Fundamental Research Funds for the Central Universities in China
文摘Several parameter identification methods of thermal response test were evaluated through numerical and experimental study.A three-dimensional finite-volume numerical model was established under the assumption that the soil thermal conductivity had been known in the simulation of thermal response test.The thermal response curve was firstly obtained through numerical calculation.Then,the accuracy of the numerical model was verified with measured data obtained through a thermal response test.Based on the numerical and experimental thermal response curves,the thermal conductivity of the soil was calculated by different parameter identification methods.The calculated results were compared with the assumed value and then the accuracy of these methods was evaluated.Furthermore,the effects of test time,variable data quality,borehole radius,initial ground temperature,and heat injection rate were analyzed.The results show that the method based on cylinder-source model has a low precision and the identified thermal conductivity decreases with an increase in borehole radius.For parameter estimation,the measuring accuracy of the initial temperature of the deep ground soil has greater effect on identified thermal conductivity.
基金Project(51321065,51479191,11672360)supported by the National Natural Science Foundation of China。
文摘Partition of unity based numerical manifold method can solve continuous and discontinuous problems in a unified framework with a two-cover system,i.e.,the mathematical cover and physical cover.However,renewal of the topology of the two-cover system poses a challenge for multiple crack propagation problems and there are few references.In this study,a robust and efficient strategy is proposed to update the cover system of the numerical manifold method in simulation of multiple crack propagation problems.The proposed algorithm updates the cover system with a bottom-up process:1)identification of fractured manifold elements according to the previous and latest crack tip position;and 2)local topological update of the manifold elements,physical patches,block boundary loops,and non-persistent joint loops according to the scenario classification of the propagating crack.The proposed crack tracking strategy and classification of the renewal cases promote a robust and efficient cover renewal algorithm for multiple crack propagation analysis.Three crack propagation examples show that the proposed algorithm performs well in updating the cover system.This cover renewal methodology can be extended for numerical manifold method with polygonal mathematical covers.
基金Project(U1261107)supported by the National Natural Science Foundation of China
文摘Combined with the kinetic model of liquid film spreading, a new numerical method of solid-liquid-gas three-phase flow was developed for the moving of contact line, which was a hybrid method of computational fluid dynamics and lattice Boltzmalm method （LBM）. By taking the effect of molecule force in droplet and the wall surface on liquid film into account, the changing law of contact angle with different surface tensions was analyzed on glass and aluminum foil surfaces. Compared with experimental results, the standard deviation by using LBM is less than 0.5°, which validates the feasibility of LBM simulation on the dynamic process of liquid film spreading. In addition, oscillations are discovered both at the initial and end phases. The phenomenon of retraction is also found and the maximum retraction angle is 7.58°. The obtained result shows that the retraction is proved to be correlative with precursor film by tracking the volume change of liquid film contour. Furthermore, non-dimensional coefficient 2 is introduced to measure the liquid film retraction capacity.
基金This project was supported by the Foundation of MOE of China (No. 00179).
文摘It is well known that the incorrect results will be given using either the electric or magnetic field integral equation to calculate the radar cross section (RCS) of a closed body at the interior resonance. In this paper, an effective iterative technique is used to correct the calculated surface current density from the electric field integral equation. The radar cross section is computed for an infinite conducting circular cylinder at the interior resonance, and the obtained results are in good agreement with the analytical results. The backscattering cross section of an infinite triangular cylinder in the vicinity of a resonant frequency is also calculated. It is shown that the presence method is efficient and accurate.
文摘A kind of international rapid field measurement methods of hydraulic conductivity and it's applications in Sanjiang Plain have been introduced in the paper.
基金This project was supported by the National Natural Science Foundation of China (19871080).
文摘A class of modified parallel combined methods of real-time numerical simulation are presented for a stiff dynamic system. By combining the parallelism across the system with the parallelism across the method, and relaxing the dependence of stage value computation on sampling time of input function, a class of modified real-time parallel combined methods are constructed. Stiff and nonstiff subsystems are solved in parallel on a parallel computer by a parallel Rosen-brock method and a parallel RK method, respectively. Their order conditions and convergences are discussed. The numerical simulation experiments show that this class of modified algorithms can get high speed and efficiency.
文摘This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solution of test equations u'(t) = a 11 u(t) + a12v(t) + b11 u(t - τ1) + b12v(t-τ2,v'(t) = a21 u(t) + a22 v(t) + b21 u(t -τ1,) + b22 v(t -τ2), t>0,with initial conditionsu(t)=u0(t),v(t) =v0(t), t≤0.where aij, bij∈C, τj >0, i,j = 1,2,, and u0(t), v0(t)are continuous and complex valued. Sufficient conditions for the asymptotic stability of test equation are derived. Furthermore, with respect to an appropriate definition of stability for the numerical method, it is proved that the linear θ-method is stable if and only if 1/2≤θ≤1 and the one-leg θ-method is stable if and only if θ= 1.
文摘In this paper, a novel calibration integral equation is derived for resolving double-sided, two-probe inverse heat conduction problem of surface heat flux estimation. In contrast to the conventional inverse heat conduction techniques, this calibration approach does not require explicit input of the probe locations, thermophysical properties of the host material and temperature sensor parameters related to thermal contact resistance, sensor capacitance and conductive lead losses. All those parameters and properties are inherently contained in the calibration framework in terms of Volterra integral equation of the first kind. The Laplace transform technique is applied and the frequency domain manipulations of the heat equation are performed for deriving the calibration integral equation. Due to the ill-posed nature, regularization is required for the inverse heat conduction problem, a future-time method or singular value decomposition (SVD) can be used for stabilizing the ill-posed Volterra integral equation of the first kind.
