To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyur...To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyurea coating).The failure characteristics and dynamic responses of the specimens were compared through conducting explosion tests.The reliability of the numerical simulation using LS-DYNA software was verified by the test results.The effects of different scaled distances,reinforcement ratios,concrete strengths,coating thicknesses and ranges of polyurea were studied.The results show that the polyurea coating can effectively enhance the anti-explosion performance of the girder.The top plate of middle chamber in specimen G forms an elliptical penetrating hole,while that in specimen PCG only shows a very slight local dent.The peak vertical displacement and residual displacement of PCG decrease by 74.8% and 73.7%,respectively,compared with those of specimen G.For the TNT explosion with small equivalent,the polyurea coating has a more significant protective effect on reducing the size of fracture.With the increase of TNT equivalent,the protective effect of polyurea on reducing girder displacement becomes more significant.The optimal reinforcement ratio,concrete strength,thickness and range of polyurea coating were also drawn.展开更多
Numerical simulation is the most powerful computational and analysis tool for a large variety of engineering and physical problems.For a complex problem relating to multi-field,multi-process and multi-scale,different ...Numerical simulation is the most powerful computational and analysis tool for a large variety of engineering and physical problems.For a complex problem relating to multi-field,multi-process and multi-scale,different computing tools have to be developed so as to solve particular fields at different scales and for different processes.Therefore,the integration of different types of software is inevitable.However,it is difficult to perform the transfer of the meshes and simulated results among software packages because of the lack of shared data formats or encrypted data formats.An image processing based method for three-dimensional model reconstruction for numerical simulation was proposed,which presents a solution to the integration problem by a series of slice or projection images obtained by the post-processing modules of the numerical simulation software.By means of mapping image pixels to meshes of either finite difference or finite element models,the geometry contour can be extracted to export the stereolithography model.The values of results,represented by color,can be deduced and assigned to the meshes.All the models with data can be directly or indirectly integrated into other software as a continued or new numerical simulation.The three-dimensional reconstruction method has been validated in numerical simulation of castings and case studies were provided in this study.展开更多
As a western boundary current, the Kuroshio is closely related to the large scale oceanic circulation and at the same time, is greatly influenced by the local topography because of its narrow width. Numerical studies ...As a western boundary current, the Kuroshio is closely related to the large scale oceanic circulation and at the same time, is greatly influenced by the local topography because of its narrow width. Numerical studies of the Kuroshio are usually confined to portions of it in different geographical regions since the computer execution time required to run a numerical model of the Pacific using a sufficiently fine grid to resolve adequately the flow structure of the Kuroshio is enormous. In order to circumvent the problems of multiple spatial scales and consistent boundary conditions, nested models are employed in which a coarse grid model of the Pacific is used to supply the open boundary conditions for a finer grid model of the northwestern Pacific to simulate the flow and temperature fields of the Kuroshio in summer and winter. The major features of the Kuroshio have in general been successfully simulated by the nested models.展开更多
Saline intrusion into marine sewage ouffalls will greatly decrease the efficiency of sewage disposal. In order to investigate the mechanisms of this flow, in this paper, a three-dimensional numerical model based on F...Saline intrusion into marine sewage ouffalls will greatly decrease the efficiency of sewage disposal. In order to investigate the mechanisms of this flow, in this paper, a three-dimensional numerical model based on FVM (Finite Volume Method) is established, The RNG κ-ε model is selected for turbulence modeling. The time-averaged vohtme fraction equations are introduced to simulate the stratification and inteffaeial exchange of sewage and seawater in outfalls. Validity of the established three-dimensional numerical model is evaluated by comparisons of numerical results with experimental data. With this three-dimensional numerical model, the internal flow characteristics in ouffalls for different sewage discharges are simulated. The results indicate that for a low sewage discharge, saline circulates in the outfall due to intrusion and both the inflowing momentum and the inteffaeial turbulent mixing are important mechanisms to extrude the saline. For a high sewage discharge, saline intrusion could be avoided. The inflow momentum is the main mechanism to extrude the saline and the inteffacial turbulent mixing is nut important relatively. Even at a high sewage discharge, the saline wedge would be retained in the main ouffall pipe after the risers are purged. It takes a long time for this saline wedge to be extruded by interracial turbulent mixing.展开更多
In accordance to the anisotropic feature of turbulent flow, an anisotropic algebraic stress model is adopted to predict the turbulent flow field and turbulent characteristics generated by a Rushton disc turbine with t...In accordance to the anisotropic feature of turbulent flow, an anisotropic algebraic stress model is adopted to predict the turbulent flow field and turbulent characteristics generated by a Rushton disc turbine with the improved inner-outer iterative procedure. The predicted turbulent flow is compared with experimental data and the simulation by the standard κ-ε turbulence model. The anisotropic algebraic stress model is found to give better prediction than the standard κ-ε turbulence model. The predicted turbulent flow field is in accordance to experimental data and the trend of the turbulence intensity can be effectively reflected in the simulation. The distribution of turbulent shear rate in the stirred tanks was simulated with the established numerical procedure.展开更多
This study is to understand the impact of operating conditions, especially initial operation temperature (T<sub>ini</sub>) which is set in a high temperature range, on the temperature profile of the interf...This study is to understand the impact of operating conditions, especially initial operation temperature (T<sub>ini</sub>) which is set in a high temperature range, on the temperature profile of the interface between the polymer electrolyte membrane (PEM) and the catalyst layer at the cathode (i.e., the reaction surface) in a single cell of polymer electrolyte fuel cell (PEFC). A 1D multi-plate heat transfer model based on the temperature data of the separator measured using the thermograph in a power generation experiment was developed to evaluate the reaction surface temperature (T<sub>react</sub>). In addition, to validate the proposed heat transfer model, T<sub>react</sub> obtained from the model was compared with that from the 3D numerical simulation using CFD software COMSOL Multiphysics which solves the continuity equation, Brinkman equation, Maxwell-Stefan equation, Butler-Volmer equation as well as heat transfer equation. As a result, the temperature gap between the results obtained by 1D heat transfer model and those obtained by 3D numerical simulation is below approximately 0.5 K. The simulation results show the change in the molar concentration of O<sub>2</sub> and H<sub>2</sub>O from the inlet to the outlet is more even with the increase in T<sub>ini</sub> due to the lower performance of O<sub>2</sub> reduction reaction. The change in the current density from the inlet to the outlet is more even with the increase in T<sub>ini</sub> and the value of current density is smaller with the increase in T<sub>ini </sub>due to the increase in ohmic over-potential and concentration over-potential. It is revealed that the change in T<sub>react</sub> from the inlet to the outlet is more even with the increase in T<sub>ini</sub> irrespective of heat transfer model. This is because the generated heat from the power generation is lower with the increase in T<sub>ini </sub>due to the lower performance of O<sub>2</sub> reduction reaction.展开更多
Numerical simulation or evolutionary history of an oil and gas-bearing basin is to repeat geological and thermodyanomic history of basin evolution on a computer and then to quantitate petroleum generation,accumulatio...Numerical simulation or evolutionary history of an oil and gas-bearing basin is to repeat geological and thermodyanomic history of basin evolution on a computer and then to quantitate petroleum generation,accumulation and migration.The mathemat-ical model describing geological and thermodynamic history of the basin evolution ischaracterised by an initial-boundary value problem of a system of nonlinear partial dif-ferential equations. In the present paper, a numerical method for three-dimensionalproblem and the analysis of its stability are established and a numerical result for apractical model is given, which shows that the abnormal pressure and paleo-temperat-ure computed are reasonable and display physical characteristics clearly as well.展开更多
The motion of gas bubbles beneath a free surface will lead to a spike of fluid on the free surface. The distance of the bubbles to the free surface is the key factor to different phenomena. When the inception distance...The motion of gas bubbles beneath a free surface will lead to a spike of fluid on the free surface. The distance of the bubbles to the free surface is the key factor to different phenomena. When the inception distance varies in some range, crown phenomenon would happen after the impact of weak buoyancy bubbles, so this kind of spike is defined as crown spike in the present paper. Based on potential flow theory, a three-dimensional numerical model is established to simulate the motion of the free-surface spike generated by one bubble or a horizontal line of two in-phase bubbles. After the downward jet formed near the end of the collapse phase, the simulation of the free surface is performed to study the crown spike without regard to the toroidal bubble's effect. Calculations about the interaction between one bubble and free surface agree well with the experimental results conducted with a high-speed camera, and relative error is within 15%. Crown spike in both single- and two-bubble cases are simulated numerically. Different features and laws of the motion of crown spike, depending on the bubble-boundary distances and the inter-bubble distances, have been investigated.展开更多
The model equations with tbree-dimensional, time-dependent, nonlinear Navier-Stokes equations are transformed by sigma-transformation.On the basis of the process splitting technique, the fluid flow problems are divid...The model equations with tbree-dimensional, time-dependent, nonlinear Navier-Stokes equations are transformed by sigma-transformation.On the basis of the process splitting technique, the fluid flow problems are divided into two parts:the vertically-intopated equations (external mode) and the vertical structure equations(internal mace). The first set of equations being the propagation of the tidal weves and the ADI numerical scheme has ben chosen to solve them. Conerning the vertical structure equations, they are solved by means of leapfrog stepping procedure.The main features of the tide and associated tidal current in the Bohai Sea are examined with this 3-D model.To have a good reproduction of vertical structure, the column is divided into 10 layers and the M2 tidal current is computed in detail. The simulation reveal the spetial structure and some important characteristics of the tidal current of the Bohai Sea. The application of the 3-D madel to forecasting of the tidal current in the Bobal Sea has been Performed as an illustration.展开更多
Cemented paste backfill(CPB)technology is a green mining method used to control underground goaves and tailings ponds.The curing process of CPB in the stope is the product of a thermo-hydro-mechanical-chemical multi-f...Cemented paste backfill(CPB)technology is a green mining method used to control underground goaves and tailings ponds.The curing process of CPB in the stope is the product of a thermo-hydro-mechanical-chemical multi-field performance interaction.At present,research on the multi-field performance of CPB mainly includes indoor similar simulation experiments,in-situ multi-field performance monitoring experiments,multi-field performance coupling model construction of CPB,and numerical simulation of the multi-field performance of CPB.Because it is hard to study the in-situ multi-field performance of CPB in the real stope,most current research on in-situ multi-field performance adopts the numerical simulation method.By simulating the conditions of CPB in the real stope(e.g.,maintenance environment,stope geometry,drainage conditions,and barricade and backfilling rates),the multi-field performance of CPB is further studied.This paper summarizes the mathematical models employed in the numerical simulation and lists the engineering application cases of numerical simulation in the in-situ multi-field performance of CPB.Finally,it proposes that the multi-field performance of CPB needs to strengthen the theoretical study of multi-field performance,form the strength design criterion based on the multi-field performance of CPB,perform a full-range numerical simulation of the multi-field performance of CPB,develop a pre-warning technology for the CPB safety of CPB,develop automatic and wireless sensors for the multi-field performance monitoring of CPB,and realize the application and popularization of CPB monitoring technology.展开更多
A three-dimensional baroclinic shelf sea model’s numerical simulation of the South China Sea (SCS) middle and deep layer circulation structure showed that: 1. In the SCS middle and deep layer, a southward boundary cu...A three-dimensional baroclinic shelf sea model’s numerical simulation of the South China Sea (SCS) middle and deep layer circulation structure showed that: 1. In the SCS middle and deep layer, a southward boundary current exists along the east shore of the Indo-China Peninsula all year long. A cyclonic eddy (gyre) is formed by the current in the above sea areas except in the middle layer in spring, when an anticyclonic eddy exists on the eastern side of the current. In the deep layer, a large-scale anticyclonic eddy often exists in the sea areas between the Zhongsha Islands and west shore of southern Luzon Island. 2. In the middle layer in summer and autumn, and in the deep layer in autumn and winter, there is an anticyclonic eddy (gyre) in the northeastern SCS, while in the middle layer in winter and spring, and in the deep layer in spring and summer, there is a cyclonic one. 3. In the middle layer, there is a weak northeastward current in the Nansha Trough in spring and summer, while in autumn and winter it evolves into an anticyclonic eddy (gyre), which then spreads westward to the whole western Nansha Islands sea areas.展开更多
By coupling the three-dimensional hydrodynamic model with the wave model, numerical simulations of the three- dimensional wave-induced current are carried out in this study. The wave model is based on the numerical so...By coupling the three-dimensional hydrodynamic model with the wave model, numerical simulations of the three- dimensional wave-induced current are carried out in this study. The wave model is based on the numerical solution of the modified wave action equation and eikonal equation, which can describe the wave refraction and diffraction. The hydrodynamic model is driven by the wave-induced radiation stresses and affected by the wave turbulence. The numerical implementation of the module has used the finite-volume schemes on unstructured grid, which provides great flexibility for modeling the waves and currents in the complex actual nearshore, and ensures the conservation of energy propagation. The applicability of the proposed model is evaluated in calculating the cases of wave set-up, longshore currents, undertow on a sloping beach, rip currents and meandering longshore currents on a tri-cuspate beach. The results indicate that it is necessary to introduce the depth-dependent radiation stresses into the numerical simulation of wave-induced currents, and comparisons show that the present model makes better prediction on the wave procedure as well as both horizontal and vertical structures in the wave-induced current field.展开更多
A finite-rate method is used to simulate the three-dimensional combustion process in a plasma generator with CH4 as the fuel. The simulation was run with RNG k-ε model to simulate turbulence, with eddy-dissipation-co...A finite-rate method is used to simulate the three-dimensional combustion process in a plasma generator with CH4 as the fuel. The simulation was run with RNG k-ε model to simulate turbulence, with eddy-dissipation-concept (EDC) model to simulate the combustion and with discrete ordinates model to simulate radiation. The numerical results show that the flow field characteristics and the parameter distributions are under the condition of rich fuels, and these results provide valuable information when optimizing the plasma generator design and organizing its flow fields.展开更多
In the present paper, a three-dimensional (3D) Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solut...In the present paper, a three-dimensional (3D) Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solution element scheme for conservative hyperbolic governing equations with source terms is given. A modified ghost fluid method is proposed for the treatment of the boundary conditions. Numerical simulations of the Taylor bar problem and the ricochet phenomenon of a sphere impacting a plate target at an angle of 60~ are carried out. The numerical results are in good agreement with the corresponding experimental observations. It is proved that our computational technique is feasible for analyzing 3D high-velocity impact problems.展开更多
Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were ex...Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were examined to find out the control factors and formation mechanisms of the salt structures. By using the three-dimensional discrete element numerical simulation method, the formation mechanisms of typical salt structures of western Kuqa foreland thrust belt in Keshen and Dabei work areas were comprehensively analyzed. The simulation results show that the salt deformation in Keshen and Dabei work areas is of forward spread type, with deformation concentrated in the piedmont zone;the salt deformation is affected by the early uplift near the compression end, pre-existing basement faults, synsedimentary process and the initial salt depocenter;in the direction perpendicular to the compression direction, salt rocks near the compression end have strong lateral mobility with the velocity component moving towards the middle part, and the closer to the middle, the larger the velocity will be, so that salt rocks will aggregate towards the middle and deform intensely, forming complex folds and separation of salt structures from salt source, and local outcrop with thrust faults. Compared with 2 D simulation, 3 D simulation can analyze salt structures in the principal stress direction and direction perpendicular to the principal stress, give us a full view of the formation mechanisms of salt structures, and guide the exploration of oil and gas reservoirs related to salt structures.展开更多
In this study, the Reynolds-averaged Navier-Stokes (RANS) method is employed to simulate the flow within and over an intersection model with three kinds of k-ε turbulence closure schemes, namely, standard model, re...In this study, the Reynolds-averaged Navier-Stokes (RANS) method is employed to simulate the flow within and over an intersection model with three kinds of k-ε turbulence closure schemes, namely, standard model, renormalization group (RNG) model and realizable k-ε model. The comparison between the simulated and observed flow fields shows that the RANS simulation with all the three turbulence models cannot completely and accurately reproduce the observed flow field in all details. A detailed comparison between the predicted profiles of wind velocities and the measured data shows that the realizble k-ε model is the best one among the three turbulence closure models in general. However, the extent to which the improvement is achieved by the realizable k-ε model is still not enough to completely and accurately describe the turbulent flow in a relatively complex environment.展开更多
The present study focuses on the breaching process and failure of barrier dams due to overtopping. In this work, a series of centrifugal model tests is presented to examine the failure mechanisms of landslide dams. Ba...The present study focuses on the breaching process and failure of barrier dams due to overtopping. In this work, a series of centrifugal model tests is presented to examine the failure mechanisms of landslide dams. Based on the experimental results, failure process and mechanism of barrier dam due to overtopping are analyzed and further verified by simulating the experimental overtopping failure process. The results indicate that the barrier dam will develop during the entire process of overtopping in the width direction, whereas the breach will cease to develop at an early stage in the depth direction because of the large particles that accumulate on the downstream slope. Moreover, headcut erosion can be clearly observed in the first two stages of overtopping, and coarsening on the downstream slope occurs in the last stage of overtopping. Thus, the bottom part of the barrier dam can survive after dam breaching and full dam failure becomes relatively rare for a barrier dam. Furthermore, the remaining breach would be smaller than that of a homogeneous cohesive dam under the same conditions.展开更多
A composite model, which is the combination of Boussinesq equations and Volume of Fluid (VOF) method, has been developed for 2-D time-domain computations of nonlinear waves in a large region. The whole computational r...A composite model, which is the combination of Boussinesq equations and Volume of Fluid (VOF) method, has been developed for 2-D time-domain computations of nonlinear waves in a large region. The whole computational region Omega is divided into two subregions. In the near-field around a structure, Omega(2), the flow is governed by 2-D Reynolds Averaged Navier-Stokes equations with a turbulence closure model of k-epsilon equations and numerically solved by the improved VOF method; whereas in the subregion Omega(1) (Omega(1) = Omega - Omega(2)) the flow is governed by one-D Boussinesq equations and numerically solved with the predictor-corrector algorithm. The velocity and the wave surface elevation are matched on the common boundary of the two subregions. Numerical tests have been conducted for the case of wave propagation and interaction with a wave barrier. It is shown that the composite model can help perform efficient computation of nonlinear waves in a large region with the complicated flow fields near structures taken into account.展开更多
The optimization of flow control devices(FCDs) for a T-type five-strand billet caster tundish was carried out by water modeling and numerical simulation. In water modeling experiments, flow characteristics of the bare...The optimization of flow control devices(FCDs) for a T-type five-strand billet caster tundish was carried out by water modeling and numerical simulation. In water modeling experiments, flow characteristics of the bare tundish and tundish conf igurations with designed U-type baff les and a round turbulence inhibitor were analyzed using residence time distribution(RTD) curves. Mathematical models for liquid steel in the real plant tundish were established using the fluid dynamics software package Fluent. The flow field, the temperature field, and the RTD curves of liquid steel in the proposed tundish conf igurations were obtained. The results of numerical simulation and water modeling were validated with each other by the predicted and experimental RTD curves. The results of flow field and temperature field were used to ref lect the actual state of a real plant tundish and to choose the optimal FCD. Finaly, from the whole performance of the multi-strand tundish, the optimal scheme was determined by combining the results of water modeling and numerical simulation. With the optimal tundish equipped with U-type baffle with def lector holes and round turbulence inhibitor, not only was the flow characteristic of each strand improved, but also the difference of flow characteristics between multiple strands was smaller.展开更多
基金the Natural Science Foundation of Jiangsu Province(Grant No.BK20200494)China Postdoctoral Science Foundation(Grant No.2021M701725)+3 种基金Jiangsu Postdoctoral Research Funding Program(Grant No.2021K522C)Fundamental Research Funds for the Central Universities(Grant No.30919011246)National Natural Science Foundation of China(Grant No.52278188)Natural Science Foundation of Jiangsu Province(Grant No.BK20211196)。
文摘To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyurea coating).The failure characteristics and dynamic responses of the specimens were compared through conducting explosion tests.The reliability of the numerical simulation using LS-DYNA software was verified by the test results.The effects of different scaled distances,reinforcement ratios,concrete strengths,coating thicknesses and ranges of polyurea were studied.The results show that the polyurea coating can effectively enhance the anti-explosion performance of the girder.The top plate of middle chamber in specimen G forms an elliptical penetrating hole,while that in specimen PCG only shows a very slight local dent.The peak vertical displacement and residual displacement of PCG decrease by 74.8% and 73.7%,respectively,compared with those of specimen G.For the TNT explosion with small equivalent,the polyurea coating has a more significant protective effect on reducing the size of fracture.With the increase of TNT equivalent,the protective effect of polyurea on reducing girder displacement becomes more significant.The optimal reinforcement ratio,concrete strength,thickness and range of polyurea coating were also drawn.
基金funded by National Key R&D Program of China(No.2021YFB3401200)the National Natural Science Foundation of China(No.51875308)the Beijing Nature Sciences Fund-Haidian Originality Cooperation Project(L212002).
文摘Numerical simulation is the most powerful computational and analysis tool for a large variety of engineering and physical problems.For a complex problem relating to multi-field,multi-process and multi-scale,different computing tools have to be developed so as to solve particular fields at different scales and for different processes.Therefore,the integration of different types of software is inevitable.However,it is difficult to perform the transfer of the meshes and simulated results among software packages because of the lack of shared data formats or encrypted data formats.An image processing based method for three-dimensional model reconstruction for numerical simulation was proposed,which presents a solution to the integration problem by a series of slice or projection images obtained by the post-processing modules of the numerical simulation software.By means of mapping image pixels to meshes of either finite difference or finite element models,the geometry contour can be extracted to export the stereolithography model.The values of results,represented by color,can be deduced and assigned to the meshes.All the models with data can be directly or indirectly integrated into other software as a continued or new numerical simulation.The three-dimensional reconstruction method has been validated in numerical simulation of castings and case studies were provided in this study.
文摘As a western boundary current, the Kuroshio is closely related to the large scale oceanic circulation and at the same time, is greatly influenced by the local topography because of its narrow width. Numerical studies of the Kuroshio are usually confined to portions of it in different geographical regions since the computer execution time required to run a numerical model of the Pacific using a sufficiently fine grid to resolve adequately the flow structure of the Kuroshio is enormous. In order to circumvent the problems of multiple spatial scales and consistent boundary conditions, nested models are employed in which a coarse grid model of the Pacific is used to supply the open boundary conditions for a finer grid model of the northwestern Pacific to simulate the flow and temperature fields of the Kuroshio in summer and winter. The major features of the Kuroshio have in general been successfully simulated by the nested models.
文摘Saline intrusion into marine sewage ouffalls will greatly decrease the efficiency of sewage disposal. In order to investigate the mechanisms of this flow, in this paper, a three-dimensional numerical model based on FVM (Finite Volume Method) is established, The RNG κ-ε model is selected for turbulence modeling. The time-averaged vohtme fraction equations are introduced to simulate the stratification and inteffaeial exchange of sewage and seawater in outfalls. Validity of the established three-dimensional numerical model is evaluated by comparisons of numerical results with experimental data. With this three-dimensional numerical model, the internal flow characteristics in ouffalls for different sewage discharges are simulated. The results indicate that for a low sewage discharge, saline circulates in the outfall due to intrusion and both the inflowing momentum and the inteffaeial turbulent mixing are important mechanisms to extrude the saline. For a high sewage discharge, saline intrusion could be avoided. The inflow momentum is the main mechanism to extrude the saline and the inteffacial turbulent mixing is nut important relatively. Even at a high sewage discharge, the saline wedge would be retained in the main ouffall pipe after the risers are purged. It takes a long time for this saline wedge to be extruded by interracial turbulent mixing.
基金the National Natural Science Foundation of China (No. 29792074).
文摘In accordance to the anisotropic feature of turbulent flow, an anisotropic algebraic stress model is adopted to predict the turbulent flow field and turbulent characteristics generated by a Rushton disc turbine with the improved inner-outer iterative procedure. The predicted turbulent flow is compared with experimental data and the simulation by the standard κ-ε turbulence model. The anisotropic algebraic stress model is found to give better prediction than the standard κ-ε turbulence model. The predicted turbulent flow field is in accordance to experimental data and the trend of the turbulence intensity can be effectively reflected in the simulation. The distribution of turbulent shear rate in the stirred tanks was simulated with the established numerical procedure.
文摘This study is to understand the impact of operating conditions, especially initial operation temperature (T<sub>ini</sub>) which is set in a high temperature range, on the temperature profile of the interface between the polymer electrolyte membrane (PEM) and the catalyst layer at the cathode (i.e., the reaction surface) in a single cell of polymer electrolyte fuel cell (PEFC). A 1D multi-plate heat transfer model based on the temperature data of the separator measured using the thermograph in a power generation experiment was developed to evaluate the reaction surface temperature (T<sub>react</sub>). In addition, to validate the proposed heat transfer model, T<sub>react</sub> obtained from the model was compared with that from the 3D numerical simulation using CFD software COMSOL Multiphysics which solves the continuity equation, Brinkman equation, Maxwell-Stefan equation, Butler-Volmer equation as well as heat transfer equation. As a result, the temperature gap between the results obtained by 1D heat transfer model and those obtained by 3D numerical simulation is below approximately 0.5 K. The simulation results show the change in the molar concentration of O<sub>2</sub> and H<sub>2</sub>O from the inlet to the outlet is more even with the increase in T<sub>ini</sub> due to the lower performance of O<sub>2</sub> reduction reaction. The change in the current density from the inlet to the outlet is more even with the increase in T<sub>ini</sub> and the value of current density is smaller with the increase in T<sub>ini </sub>due to the increase in ohmic over-potential and concentration over-potential. It is revealed that the change in T<sub>react</sub> from the inlet to the outlet is more even with the increase in T<sub>ini</sub> irrespective of heat transfer model. This is because the generated heat from the power generation is lower with the increase in T<sub>ini </sub>due to the lower performance of O<sub>2</sub> reduction reaction.
文摘Numerical simulation or evolutionary history of an oil and gas-bearing basin is to repeat geological and thermodyanomic history of basin evolution on a computer and then to quantitate petroleum generation,accumulation and migration.The mathemat-ical model describing geological and thermodynamic history of the basin evolution ischaracterised by an initial-boundary value problem of a system of nonlinear partial dif-ferential equations. In the present paper, a numerical method for three-dimensionalproblem and the analysis of its stability are established and a numerical result for apractical model is given, which shows that the abnormal pressure and paleo-temperat-ure computed are reasonable and display physical characteristics clearly as well.
基金Project supported by the Major Basic Research Project of National Security of China(Grant No.613157)the Excellent Young Scientists Fund of China(Grant No.51222904)
文摘The motion of gas bubbles beneath a free surface will lead to a spike of fluid on the free surface. The distance of the bubbles to the free surface is the key factor to different phenomena. When the inception distance varies in some range, crown phenomenon would happen after the impact of weak buoyancy bubbles, so this kind of spike is defined as crown spike in the present paper. Based on potential flow theory, a three-dimensional numerical model is established to simulate the motion of the free-surface spike generated by one bubble or a horizontal line of two in-phase bubbles. After the downward jet formed near the end of the collapse phase, the simulation of the free surface is performed to study the crown spike without regard to the toroidal bubble's effect. Calculations about the interaction between one bubble and free surface agree well with the experimental results conducted with a high-speed camera, and relative error is within 15%. Crown spike in both single- and two-bubble cases are simulated numerically. Different features and laws of the motion of crown spike, depending on the bubble-boundary distances and the inter-bubble distances, have been investigated.
文摘The model equations with tbree-dimensional, time-dependent, nonlinear Navier-Stokes equations are transformed by sigma-transformation.On the basis of the process splitting technique, the fluid flow problems are divided into two parts:the vertically-intopated equations (external mode) and the vertical structure equations(internal mace). The first set of equations being the propagation of the tidal weves and the ADI numerical scheme has ben chosen to solve them. Conerning the vertical structure equations, they are solved by means of leapfrog stepping procedure.The main features of the tide and associated tidal current in the Bohai Sea are examined with this 3-D model.To have a good reproduction of vertical structure, the column is divided into 10 layers and the M2 tidal current is computed in detail. The simulation reveal the spetial structure and some important characteristics of the tidal current of the Bohai Sea. The application of the 3-D madel to forecasting of the tidal current in the Bobal Sea has been Performed as an illustration.
基金financially supported by the National Natural Science Foundation of China(No.52130404)the State Key Laboratory for GeoMechanics and Deep Underground Engineering China University of Mining&Technology,Beijing(No.SKLGDUEK2127)+1 种基金the Fundamental Research Funds for the Central Universities(Nos.FRF-TP-19-002C2Z,FRF-IDRY-20-031)and the Fundamental Research Funds for the Central Universities and the Youth Teacher International Exchange&Growth Program(No.QNXM20220002)。
文摘Cemented paste backfill(CPB)technology is a green mining method used to control underground goaves and tailings ponds.The curing process of CPB in the stope is the product of a thermo-hydro-mechanical-chemical multi-field performance interaction.At present,research on the multi-field performance of CPB mainly includes indoor similar simulation experiments,in-situ multi-field performance monitoring experiments,multi-field performance coupling model construction of CPB,and numerical simulation of the multi-field performance of CPB.Because it is hard to study the in-situ multi-field performance of CPB in the real stope,most current research on in-situ multi-field performance adopts the numerical simulation method.By simulating the conditions of CPB in the real stope(e.g.,maintenance environment,stope geometry,drainage conditions,and barricade and backfilling rates),the multi-field performance of CPB is further studied.This paper summarizes the mathematical models employed in the numerical simulation and lists the engineering application cases of numerical simulation in the in-situ multi-field performance of CPB.Finally,it proposes that the multi-field performance of CPB needs to strengthen the theoretical study of multi-field performance,form the strength design criterion based on the multi-field performance of CPB,perform a full-range numerical simulation of the multi-field performance of CPB,develop a pre-warning technology for the CPB safety of CPB,develop automatic and wireless sensors for the multi-field performance monitoring of CPB,and realize the application and popularization of CPB monitoring technology.
文摘A three-dimensional baroclinic shelf sea model’s numerical simulation of the South China Sea (SCS) middle and deep layer circulation structure showed that: 1. In the SCS middle and deep layer, a southward boundary current exists along the east shore of the Indo-China Peninsula all year long. A cyclonic eddy (gyre) is formed by the current in the above sea areas except in the middle layer in spring, when an anticyclonic eddy exists on the eastern side of the current. In the deep layer, a large-scale anticyclonic eddy often exists in the sea areas between the Zhongsha Islands and west shore of southern Luzon Island. 2. In the middle layer in summer and autumn, and in the deep layer in autumn and winter, there is an anticyclonic eddy (gyre) in the northeastern SCS, while in the middle layer in winter and spring, and in the deep layer in spring and summer, there is a cyclonic one. 3. In the middle layer, there is a weak northeastward current in the Nansha Trough in spring and summer, while in autumn and winter it evolves into an anticyclonic eddy (gyre), which then spreads westward to the whole western Nansha Islands sea areas.
基金financially supported by the the National Natural Science Foundation of China(Grant No.51709054)the Public Science and Technology Research Funds Projects of Ocean(Grant Nos.201405025 and 201505019)
文摘By coupling the three-dimensional hydrodynamic model with the wave model, numerical simulations of the three- dimensional wave-induced current are carried out in this study. The wave model is based on the numerical solution of the modified wave action equation and eikonal equation, which can describe the wave refraction and diffraction. The hydrodynamic model is driven by the wave-induced radiation stresses and affected by the wave turbulence. The numerical implementation of the module has used the finite-volume schemes on unstructured grid, which provides great flexibility for modeling the waves and currents in the complex actual nearshore, and ensures the conservation of energy propagation. The applicability of the proposed model is evaluated in calculating the cases of wave set-up, longshore currents, undertow on a sloping beach, rip currents and meandering longshore currents on a tri-cuspate beach. The results indicate that it is necessary to introduce the depth-dependent radiation stresses into the numerical simulation of wave-induced currents, and comparisons show that the present model makes better prediction on the wave procedure as well as both horizontal and vertical structures in the wave-induced current field.
文摘A finite-rate method is used to simulate the three-dimensional combustion process in a plasma generator with CH4 as the fuel. The simulation was run with RNG k-ε model to simulate turbulence, with eddy-dissipation-concept (EDC) model to simulate the combustion and with discrete ordinates model to simulate radiation. The numerical results show that the flow field characteristics and the parameter distributions are under the condition of rich fuels, and these results provide valuable information when optimizing the plasma generator design and organizing its flow fields.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10732010,10972010,and 11332002)
文摘In the present paper, a three-dimensional (3D) Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solution element scheme for conservative hyperbolic governing equations with source terms is given. A modified ghost fluid method is proposed for the treatment of the boundary conditions. Numerical simulations of the Taylor bar problem and the ricochet phenomenon of a sphere impacting a plate target at an angle of 60~ are carried out. The numerical results are in good agreement with the corresponding experimental observations. It is proved that our computational technique is feasible for analyzing 3D high-velocity impact problems.
基金Supported by the China National Science and Technology Major Project(2016ZX05033002,2016ZX05033001).
文摘Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were examined to find out the control factors and formation mechanisms of the salt structures. By using the three-dimensional discrete element numerical simulation method, the formation mechanisms of typical salt structures of western Kuqa foreland thrust belt in Keshen and Dabei work areas were comprehensively analyzed. The simulation results show that the salt deformation in Keshen and Dabei work areas is of forward spread type, with deformation concentrated in the piedmont zone;the salt deformation is affected by the early uplift near the compression end, pre-existing basement faults, synsedimentary process and the initial salt depocenter;in the direction perpendicular to the compression direction, salt rocks near the compression end have strong lateral mobility with the velocity component moving towards the middle part, and the closer to the middle, the larger the velocity will be, so that salt rocks will aggregate towards the middle and deform intensely, forming complex folds and separation of salt structures from salt source, and local outcrop with thrust faults. Compared with 2 D simulation, 3 D simulation can analyze salt structures in the principal stress direction and direction perpendicular to the principal stress, give us a full view of the formation mechanisms of salt structures, and guide the exploration of oil and gas reservoirs related to salt structures.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 40233030, 40405004 and 40405014) and the Special Program of the Scientific and Social Practices for Graduate Students in Chinese Academy of Sciences, China.
文摘In this study, the Reynolds-averaged Navier-Stokes (RANS) method is employed to simulate the flow within and over an intersection model with three kinds of k-ε turbulence closure schemes, namely, standard model, renormalization group (RNG) model and realizable k-ε model. The comparison between the simulated and observed flow fields shows that the RANS simulation with all the three turbulence models cannot completely and accurately reproduce the observed flow field in all details. A detailed comparison between the predicted profiles of wind velocities and the measured data shows that the realizble k-ε model is the best one among the three turbulence closure models in general. However, the extent to which the improvement is achieved by the realizable k-ε model is still not enough to completely and accurately describe the turbulent flow in a relatively complex environment.
基金financial support from the National Natural Science Foundation of China (Grant No. 51709025)the Chongqing Science and Technology Commission of China (Grant No. cstc2018jcyjAX0084, cstc2018jcyjAX0391 and cstc2016jcyjA0551)Open Research Fund of Key Laboratory of Failure Mechanism and Safety Control Techniques of Earth-Rock Dam of the Ministry of Water Resources (Grant No. YK319006)
文摘The present study focuses on the breaching process and failure of barrier dams due to overtopping. In this work, a series of centrifugal model tests is presented to examine the failure mechanisms of landslide dams. Based on the experimental results, failure process and mechanism of barrier dam due to overtopping are analyzed and further verified by simulating the experimental overtopping failure process. The results indicate that the barrier dam will develop during the entire process of overtopping in the width direction, whereas the breach will cease to develop at an early stage in the depth direction because of the large particles that accumulate on the downstream slope. Moreover, headcut erosion can be clearly observed in the first two stages of overtopping, and coarsening on the downstream slope occurs in the last stage of overtopping. Thus, the bottom part of the barrier dam can survive after dam breaching and full dam failure becomes relatively rare for a barrier dam. Furthermore, the remaining breach would be smaller than that of a homogeneous cohesive dam under the same conditions.
基金Trans-Century Training program Fund for the Talent,Ministry of Education of China
文摘A composite model, which is the combination of Boussinesq equations and Volume of Fluid (VOF) method, has been developed for 2-D time-domain computations of nonlinear waves in a large region. The whole computational region Omega is divided into two subregions. In the near-field around a structure, Omega(2), the flow is governed by 2-D Reynolds Averaged Navier-Stokes equations with a turbulence closure model of k-epsilon equations and numerically solved by the improved VOF method; whereas in the subregion Omega(1) (Omega(1) = Omega - Omega(2)) the flow is governed by one-D Boussinesq equations and numerically solved with the predictor-corrector algorithm. The velocity and the wave surface elevation are matched on the common boundary of the two subregions. Numerical tests have been conducted for the case of wave propagation and interaction with a wave barrier. It is shown that the composite model can help perform efficient computation of nonlinear waves in a large region with the complicated flow fields near structures taken into account.
基金financially supported by the National Natural Science Foundation of China(51504002)
文摘The optimization of flow control devices(FCDs) for a T-type five-strand billet caster tundish was carried out by water modeling and numerical simulation. In water modeling experiments, flow characteristics of the bare tundish and tundish conf igurations with designed U-type baff les and a round turbulence inhibitor were analyzed using residence time distribution(RTD) curves. Mathematical models for liquid steel in the real plant tundish were established using the fluid dynamics software package Fluent. The flow field, the temperature field, and the RTD curves of liquid steel in the proposed tundish conf igurations were obtained. The results of numerical simulation and water modeling were validated with each other by the predicted and experimental RTD curves. The results of flow field and temperature field were used to ref lect the actual state of a real plant tundish and to choose the optimal FCD. Finaly, from the whole performance of the multi-strand tundish, the optimal scheme was determined by combining the results of water modeling and numerical simulation. With the optimal tundish equipped with U-type baffle with def lector holes and round turbulence inhibitor, not only was the flow characteristic of each strand improved, but also the difference of flow characteristics between multiple strands was smaller.