"Generalized mobility"is used to realize the unification of tube flow and seepage in form and the unification of commonly used linear and nonlinear flow laws in form,which makes it possible to use the same f..."Generalized mobility"is used to realize the unification of tube flow and seepage in form and the unification of commonly used linear and nonlinear flow laws in form,which makes it possible to use the same form of motion equations to construct unified governing equations for reservoirs of different scales in different regions.Firstly,by defining the generalized mobility under different flow conditions,the basic equation governing fluid flow in reservoir coupling generalized tube flow and seepage is established.Secondly,two typical well test analysis models for coupling tube flow and seepage flow are given,namely,pipe-shaped composite reservoir model and partially open cylindrical reservoir model.The log-log pressure draw-down type-curve of composite pipe-shaped reservoir model can show characteristics of two sets of linear flow.The log-log pressure drawdown plot of partially opened cylindrical reservoir model can show the characteristics of spherical flow and linear flow,as well as spherical flow and radial flow.The pressure build-up derivative curves of the two models basically coincide with their respective pressure drawdown derivative curves in the early stage,pulling down features in the late stage,and the shorter the production time is,the earlier the pulling down feature appears.Finally,the practicability and reliability of the models presented in this paper are verified by three application examples.展开更多
A number of critical problems of topology optimization concerning the thermostructural coupling field axe studied at length. The governing equations and topology optimization model for the thermal-structural coupling ...A number of critical problems of topology optimization concerning the thermostructural coupling field axe studied at length. The governing equations and topology optimization model for the thermal-structural coupling field axe derived, with an adjoint method for sensitivity analysis of the thermo-structural coupling field proposed. The optimization algorithm for coupling field topology optimization is investigated and a flowchart of coupling field topology optimization presented. The theory and algorithms axe implemented and verified by two numerical examples.展开更多
The coupling behavior of the imbedded weapon store occurring between the local unsteady flow field round the store and the structure response on the processing of opening its bay-door is simulated by using numerical m...The coupling behavior of the imbedded weapon store occurring between the local unsteady flow field round the store and the structure response on the processing of opening its bay-door is simulated by using numerical method based on computational fluid mechanics(CFD).The transient aerodynamic behaviors when opening door under various flight altitudes and the corresponding structure deformation evolution in the unsteady flow fields are analyzed respectively and presented.The rules of aircraft attitude parameters′impacting to the responses of structure and the bay-door′s opening process are obtained by comparing with the analysis results.These rules can be applied to the structure design of bay-door and route specification of missile when disengaged and launched from within store.展开更多
In this paper,we first establish a new fractional magnetohydrodynamic(MHD)coupled flow and heat transfer model for a generalized second-grade fluid.This coupled model consists of a fractional momentum equation and a h...In this paper,we first establish a new fractional magnetohydrodynamic(MHD)coupled flow and heat transfer model for a generalized second-grade fluid.This coupled model consists of a fractional momentum equation and a heat conduction equation with a generalized form of Fourier law.The second-order fractional backward difference formula is applied to the temporal discretization and the Legendre spectral method is used for the spatial discretization.The fully discrete scheme is proved to be stable and convergent with an accuracy of O(τ^(2)+N-r),whereτis the time step-size and N is the polynomial degree.To reduce the memory requirements and computational cost,a fast method is developed,which is based on a globally uniform approximation of the trapezoidal rule for integrals on the real line.The strict convergence of the numerical scheme with this fast method is proved.We present the results of several numerical experiments to verify the effectiveness of the proposed method.Finally,we simulate the unsteady fractional MHD flow and heat transfer of the generalized second-grade fluid through a porous medium.The effects of the relevant parameters on the velocity and temperature are presented and analyzed in detail.展开更多
Two full 3D steady mathematical models are developed by finite element method (FEM) to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the...Two full 3D steady mathematical models are developed by finite element method (FEM) to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the acquired electromagnetic force as source body forces in Navier-Stokes equations. Effects caused by the ferromagnetic shell, busbar system around, and open boundary problem as well as inside induced current were considered in terms of the magnetic field. Furthermore, a new modeling method is found to set up solid models and then mesh them entirely with so-called structuralized grids, namely hex-mesh. Examples of 75kA prebaked cell with two kinds of busbar arrangements are presented. Results agree with those disclosed in the literature and confirm that the coupled simulation is valid. It is also concluded that the usage of these models facilitates the consistent analysis of the electric field to magnetic field and then flow motion to the greater extent, local distributions of current density and magnetic flux density are very much dependent on the cell structure, the steel shell is a shield to reduce the magnetic field and flow pattern is two dimensional in the main body of the metal pad.展开更多
In this study, rapid drawdown scenarios were analyzed by means of numerical examples as well as modeling of real cases with in situ measurements. The aim of the study was to evaluate different approaches available for...In this study, rapid drawdown scenarios were analyzed by means of numerical examples as well as modeling of real cases with in situ measurements. The aim of the study was to evaluate different approaches available for calculating pore water pressure distributions during and after a drawdown. To do that, a single slope subjected to a drawdown was first analyzed under different calculation alternatives, and numerical results were discussed. Simple methods, such as undrained analysis and pure flow analysis, implicitly assuming a rigid soil skeleton, lead to significant errors in pore water pressure distributions when compared with coupled flow-deformation analysis. A similar analysis was performed for the upstream slope of the Glen Shira Dam, Scotland, and numerical results were compared with field measurements during a controlled drawdown. Field records indicate that classical undrained calculations are conservative but unrealistic. Then, a recent case of a major landslide triggered by a rapid drawdown in a reservoir was interpreted. A key aspect of the case was the correct characterization of permeability of a representative soil profile. This was achieved by combining laboratory test results and a back analysis of pore water pressure time records during a period of reservoir water level fluctuations. The results highlight the difficulty of predicting whether the pore water pressure is overestimated or underestimated when using simplified approaches, and it is concluded that predicting the pore water pressure distribution in a slope after a rapid drawdown requires a coupled flow-deformation analysis in saturated and unsaturated porous media.展开更多
Based on the new viewpoint of solid and gas interaction mechanics, gas leakage in a double deformable coal seam can be understood. That is, under the action of geophysical fields, the methane flow in a double deformab...Based on the new viewpoint of solid and gas interaction mechanics, gas leakage in a double deformable coal seam can be understood. That is, under the action of geophysical fields, the methane flow in a double deformable coal seam can be essentially considered to be compressible with time dependent and mixed permeation and diffusion through a pore cleat deformable heterogeneous and anisotropy medium. Based on this new viewpoint, a coupled mathematical model for coal seam deformation and gas leakage in a double coal seam was formulated and numerical simulations for gas emission from the coal seam are presented. It is found that coupled models might be closer to reality.展开更多
碎屑流是我国山区最危险的地质灾害之一,山区桥墩常受到碎屑流冲击而开裂、倾斜甚至倒塌,给山区桥梁建设、运营带来严重的安全隐患。采用离散元方法(discrete element method,DEM)和有限元方法(finite element method,FEM)耦合的三维数...碎屑流是我国山区最危险的地质灾害之一,山区桥墩常受到碎屑流冲击而开裂、倾斜甚至倒塌,给山区桥梁建设、运营带来严重的安全隐患。采用离散元方法(discrete element method,DEM)和有限元方法(finite element method,FEM)耦合的三维数值模拟方法模拟了碎屑流对双柱式桥墩的冲击效应,并结合斜槽试验,验证了耦合方法的准确性,进一步分析了碎屑流冲击坡度、距离和体积密度对桥墩冲击力的影响规律。结果表明,最大冲击力与碎屑流冲击坡度、距离和体积密度分别呈幂函数(指数大于1)、幂函数(指数小于1)和线性正相关。冲击坡度、距离和体积密度对最大冲击力的敏感度值分别为3.012、0.202、0.804,在桥梁碎屑流灾害防治时需重视冲击坡度和体积密度的影响。将冲击力的数值模拟值与流体动力学模型预测值对比分析表明,流体动力学模型理论公式能较好地预测桥墩所受的最大冲击力,最大预测误差低于23.6%。相关研究结果可为山区桥梁碎屑流灾害防治与设计提供一定的参考依据。展开更多
Heat exchanger systems(HXSs)or heat recovery steam generators(HRSGs)are commonly used in 100 kW to 50 MW combined cooling,heating,and power(CCHP)systems.Power flow coupling(PFC)is found in HXSs and is complex for rese...Heat exchanger systems(HXSs)or heat recovery steam generators(HRSGs)are commonly used in 100 kW to 50 MW combined cooling,heating,and power(CCHP)systems.Power flow coupling(PFC)is found in HXSs and is complex for researchers to quantify.This could possibly mislead the dispatch schedule and result in the inaccurate dispatch.PFC is caused by the inlet and outlet temperatures of each component,gas flow pressure variation,conductive medium flow rate,and atmosphere condition variation.In this paper,the expression of PFC is built by using quadratic functions to fit the non!inearit>of thermal dynamics.While fitting the model,the environmental condition needs prediction,which is calculated using phase space reconstruction(PSR)Kalman filter.In order to solve the complex quadratic dispatch model,a hybrid following electricity load(FEL)and following thermal load(FTL)mode for reducing the dimension of dispatch model,and a feasible zone analysis(FZA)method are proposed.As a result,the PFC problem of CCHP system is solved,and the dispatch cost,investment cost,and the maximum power requirements are optimized.In this paper,a case in Jinan,China is studied.The PFC model is proven to be more precise and accurate compared with traditional models.展开更多
In the design of rock sheds for the mitigation of risk due to rapid and long landslides, a crucial role is played by the evaluation of the impact force exerted by the flowing mass on the rock sheds. This paper is focu...In the design of rock sheds for the mitigation of risk due to rapid and long landslides, a crucial role is played by the evaluation of the impact force exerted by the flowing mass on the rock sheds. This paper is focused on the influencing factors of the impact force of dry granular flow onto rock shed and in particular on the evaluation of the maximum impact force. The coupled DEM-FEM model calibrated with small-scale physical experiment is used to simulate the movement of dry granular flow coupled with impact forces on the rock-shed. Based on the numerical results, three key stages were identified of impact process, namely startup streams slippery, impact and pile-up. The maximum impact force increases linearly with bulk density, and the maximum impact force exhibits a power law dependence on the impact height and slop angle respectively. The sensitivities of bulk density, impact height, and slope angle on the maximum impact force are: 1.0, 0.496, and 2.32 respectively in the benchmark model. The parameters with high sensitivity should be given priority in the design of the rock shed. The results obtained from this study are useful for facilitating design of shed against dry granular flow.展开更多
基金Supported by the Scientific Research Project of Key Laboratory of Shaanxi Provincial Department of Education(13JS090)。
文摘"Generalized mobility"is used to realize the unification of tube flow and seepage in form and the unification of commonly used linear and nonlinear flow laws in form,which makes it possible to use the same form of motion equations to construct unified governing equations for reservoirs of different scales in different regions.Firstly,by defining the generalized mobility under different flow conditions,the basic equation governing fluid flow in reservoir coupling generalized tube flow and seepage is established.Secondly,two typical well test analysis models for coupling tube flow and seepage flow are given,namely,pipe-shaped composite reservoir model and partially open cylindrical reservoir model.The log-log pressure draw-down type-curve of composite pipe-shaped reservoir model can show characteristics of two sets of linear flow.The log-log pressure drawdown plot of partially opened cylindrical reservoir model can show the characteristics of spherical flow and linear flow,as well as spherical flow and radial flow.The pressure build-up derivative curves of the two models basically coincide with their respective pressure drawdown derivative curves in the early stage,pulling down features in the late stage,and the shorter the production time is,the earlier the pulling down feature appears.Finally,the practicability and reliability of the models presented in this paper are verified by three application examples.
基金Project supported by the National Natural Postdoctor Scientific Foundation (No.2005037347)973 Project of China (No.2003CB716207).
文摘A number of critical problems of topology optimization concerning the thermostructural coupling field axe studied at length. The governing equations and topology optimization model for the thermal-structural coupling field axe derived, with an adjoint method for sensitivity analysis of the thermo-structural coupling field proposed. The optimization algorithm for coupling field topology optimization is investigated and a flowchart of coupling field topology optimization presented. The theory and algorithms axe implemented and verified by two numerical examples.
文摘The coupling behavior of the imbedded weapon store occurring between the local unsteady flow field round the store and the structure response on the processing of opening its bay-door is simulated by using numerical method based on computational fluid mechanics(CFD).The transient aerodynamic behaviors when opening door under various flight altitudes and the corresponding structure deformation evolution in the unsteady flow fields are analyzed respectively and presented.The rules of aircraft attitude parameters′impacting to the responses of structure and the bay-door′s opening process are obtained by comparing with the analysis results.These rules can be applied to the structure design of bay-door and route specification of missile when disengaged and launched from within store.
基金supported by the Project of the National Key R&D Program(Grant No.2021YFA1000202)National Natural Science Foundation of China(Grant Nos.12120101001,12001326 and 12171283)+2 种基金Natural Science Foundation of Shandong Province(Grant Nos.ZR2021ZD03,ZR2020QA032 and ZR2019ZD42)China Postdoctoral Science Foundation(Grant Nos.BX20190191 and 2020M672038)the Startup Fund from Shandong University(Grant No.11140082063130)。
文摘In this paper,we first establish a new fractional magnetohydrodynamic(MHD)coupled flow and heat transfer model for a generalized second-grade fluid.This coupled model consists of a fractional momentum equation and a heat conduction equation with a generalized form of Fourier law.The second-order fractional backward difference formula is applied to the temporal discretization and the Legendre spectral method is used for the spatial discretization.The fully discrete scheme is proved to be stable and convergent with an accuracy of O(τ^(2)+N-r),whereτis the time step-size and N is the polynomial degree.To reduce the memory requirements and computational cost,a fast method is developed,which is based on a globally uniform approximation of the trapezoidal rule for integrals on the real line.The strict convergence of the numerical scheme with this fast method is proved.We present the results of several numerical experiments to verify the effectiveness of the proposed method.Finally,we simulate the unsteady fractional MHD flow and heat transfer of the generalized second-grade fluid through a porous medium.The effects of the relevant parameters on the velocity and temperature are presented and analyzed in detail.
基金the National High Technical Reasearch and Development Programme of China (No. 2003AA327140) the National Natural Science Foundation of China (No. 50374081).
文摘Two full 3D steady mathematical models are developed by finite element method (FEM) to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the acquired electromagnetic force as source body forces in Navier-Stokes equations. Effects caused by the ferromagnetic shell, busbar system around, and open boundary problem as well as inside induced current were considered in terms of the magnetic field. Furthermore, a new modeling method is found to set up solid models and then mesh them entirely with so-called structuralized grids, namely hex-mesh. Examples of 75kA prebaked cell with two kinds of busbar arrangements are presented. Results agree with those disclosed in the literature and confirm that the coupled simulation is valid. It is also concluded that the usage of these models facilitates the consistent analysis of the electric field to magnetic field and then flow motion to the greater extent, local distributions of current density and magnetic flux density are very much dependent on the cell structure, the steel shell is a shield to reduce the magnetic field and flow pattern is two dimensional in the main body of the metal pad.
文摘In this study, rapid drawdown scenarios were analyzed by means of numerical examples as well as modeling of real cases with in situ measurements. The aim of the study was to evaluate different approaches available for calculating pore water pressure distributions during and after a drawdown. To do that, a single slope subjected to a drawdown was first analyzed under different calculation alternatives, and numerical results were discussed. Simple methods, such as undrained analysis and pure flow analysis, implicitly assuming a rigid soil skeleton, lead to significant errors in pore water pressure distributions when compared with coupled flow-deformation analysis. A similar analysis was performed for the upstream slope of the Glen Shira Dam, Scotland, and numerical results were compared with field measurements during a controlled drawdown. Field records indicate that classical undrained calculations are conservative but unrealistic. Then, a recent case of a major landslide triggered by a rapid drawdown in a reservoir was interpreted. A key aspect of the case was the correct characterization of permeability of a representative soil profile. This was achieved by combining laboratory test results and a back analysis of pore water pressure time records during a period of reservoir water level fluctuations. The results highlight the difficulty of predicting whether the pore water pressure is overestimated or underestimated when using simplified approaches, and it is concluded that predicting the pore water pressure distribution in a slope after a rapid drawdown requires a coupled flow-deformation analysis in saturated and unsaturated porous media.
文摘Based on the new viewpoint of solid and gas interaction mechanics, gas leakage in a double deformable coal seam can be understood. That is, under the action of geophysical fields, the methane flow in a double deformable coal seam can be essentially considered to be compressible with time dependent and mixed permeation and diffusion through a pore cleat deformable heterogeneous and anisotropy medium. Based on this new viewpoint, a coupled mathematical model for coal seam deformation and gas leakage in a double coal seam was formulated and numerical simulations for gas emission from the coal seam are presented. It is found that coupled models might be closer to reality.
文摘碎屑流是我国山区最危险的地质灾害之一,山区桥墩常受到碎屑流冲击而开裂、倾斜甚至倒塌,给山区桥梁建设、运营带来严重的安全隐患。采用离散元方法(discrete element method,DEM)和有限元方法(finite element method,FEM)耦合的三维数值模拟方法模拟了碎屑流对双柱式桥墩的冲击效应,并结合斜槽试验,验证了耦合方法的准确性,进一步分析了碎屑流冲击坡度、距离和体积密度对桥墩冲击力的影响规律。结果表明,最大冲击力与碎屑流冲击坡度、距离和体积密度分别呈幂函数(指数大于1)、幂函数(指数小于1)和线性正相关。冲击坡度、距离和体积密度对最大冲击力的敏感度值分别为3.012、0.202、0.804,在桥梁碎屑流灾害防治时需重视冲击坡度和体积密度的影响。将冲击力的数值模拟值与流体动力学模型预测值对比分析表明,流体动力学模型理论公式能较好地预测桥墩所受的最大冲击力,最大预测误差低于23.6%。相关研究结果可为山区桥梁碎屑流灾害防治与设计提供一定的参考依据。
基金the National Natural Science Foundation of China(No.61733010).
文摘Heat exchanger systems(HXSs)or heat recovery steam generators(HRSGs)are commonly used in 100 kW to 50 MW combined cooling,heating,and power(CCHP)systems.Power flow coupling(PFC)is found in HXSs and is complex for researchers to quantify.This could possibly mislead the dispatch schedule and result in the inaccurate dispatch.PFC is caused by the inlet and outlet temperatures of each component,gas flow pressure variation,conductive medium flow rate,and atmosphere condition variation.In this paper,the expression of PFC is built by using quadratic functions to fit the non!inearit>of thermal dynamics.While fitting the model,the environmental condition needs prediction,which is calculated using phase space reconstruction(PSR)Kalman filter.In order to solve the complex quadratic dispatch model,a hybrid following electricity load(FEL)and following thermal load(FTL)mode for reducing the dimension of dispatch model,and a feasible zone analysis(FZA)method are proposed.As a result,the PFC problem of CCHP system is solved,and the dispatch cost,investment cost,and the maximum power requirements are optimized.In this paper,a case in Jinan,China is studied.The PFC model is proven to be more precise and accurate compared with traditional models.
文摘In the design of rock sheds for the mitigation of risk due to rapid and long landslides, a crucial role is played by the evaluation of the impact force exerted by the flowing mass on the rock sheds. This paper is focused on the influencing factors of the impact force of dry granular flow onto rock shed and in particular on the evaluation of the maximum impact force. The coupled DEM-FEM model calibrated with small-scale physical experiment is used to simulate the movement of dry granular flow coupled with impact forces on the rock-shed. Based on the numerical results, three key stages were identified of impact process, namely startup streams slippery, impact and pile-up. The maximum impact force increases linearly with bulk density, and the maximum impact force exhibits a power law dependence on the impact height and slop angle respectively. The sensitivities of bulk density, impact height, and slope angle on the maximum impact force are: 1.0, 0.496, and 2.32 respectively in the benchmark model. The parameters with high sensitivity should be given priority in the design of the rock shed. The results obtained from this study are useful for facilitating design of shed against dry granular flow.