Short pitch corrugation has been a problem for railways worldwide over one century.In this paper,a parametric investigation of fastenings is conducted to understand the corrugation formation mechanism and gain insight...Short pitch corrugation has been a problem for railways worldwide over one century.In this paper,a parametric investigation of fastenings is conducted to understand the corrugation formation mechanism and gain insights into corrugation mitigation.A three-dimensional finite element vehicle-track dynamic interaction model is employed,which considers the coupling between the structural dynamics and the contact mechanics,while the damage mechanism is assumed to be differential wear.Various fastening models with different configurations,boundary conditions,and parameters of stiffness and damping are built up and analysed.These models may represent different service stages of fastenings in the field.Besides,the effect of train speeds on corrugation features is studied.The results indicate:(1)Fastening parameters and modelling play an important role in corrugation formation.(2)The fastening longitudinal constraint to the rail is the major factor that determines the corrugation formation.The fastening vertical and lateral constraints influence corrugation features in terms of spatial distribution and wavelength components.(3)The strengthening of fastening constraints in the longitudinal dimension helps to mitigate corrugation.Meanwhile,the inner fastening constraint in the lateral direction is necessary for corrugation alleviation.(4)The increase in fastening longitudinal stiffness and damping can reduce the vibration amplitudes of longitudinal compression modes and thus reduce the track corrugation propensity.The simulation in this work can well explain the field corrugation in terms of the occurrence possibility and major wavelength components.It can also explain the field data with respect to the small variation between the corrugation wavelength and train speed,which is caused by frequency selection and jump between rail longitudinal compression modes.展开更多
The model of dislocations was used to construct the model of grain boundary (GB) with pure rare earths, and rare earth elements and impurities. The influence of the interaction between rare earth elements and impuriti...The model of dislocations was used to construct the model of grain boundary (GB) with pure rare earths, and rare earth elements and impurities. The influence of the interaction between rare earth elements and impurities on the cohesive properties of 5.3° low angle GB of Fe was investigated by the recursion method. The calculated results of environment sensitive embeding energy( E ESE ) show that the preferential segregation of rare earth elements towards GBs exists. Calculations of bond order integrals (BOI) show that rare earth elements increase the cohesive strength of low angle GB, and impurities such as S, P weaken the intergranular cohesion of the GB. So rare earth element of proper quantity added in steel not only cleanses other harmful impurities off the GBs, but also enhances the intergranular cohesion. This elucidates the action mechanism of rare earth elements in steel from electronic level and offers theoretical evidence for applications of rare earth elements in steels.展开更多
This paper uses a mathematical method to develop an analytical solution to the local buckling behaviour of long rectangular plates resting on tensionless elastic Winkler foundations and under combined uniform longitud...This paper uses a mathematical method to develop an analytical solution to the local buckling behaviour of long rectangular plates resting on tensionless elastic Winkler foundations and under combined uniform longitudinal uniaxial compressive and uniform in-plane shear loads. Fitted formulas are derived for plates with clamped edges and simplified supported edges. Two examples are given to demonstrate the application of the current method: one is a plate on tensionless spring foundations and the other is the contact between the steel sheet and elastic solid foundation. Finite element (FE) analysis is also conducted to validate the analytical results. Good agreement is obtained between the current method and FE analysis.展开更多
The percolation fields constructed around the elements of a cluster system in the phase spaces of properties are studied.It is shown that such neighborhoods significantly increase the number of structure parameters of...The percolation fields constructed around the elements of a cluster system in the phase spaces of properties are studied.It is shown that such neighborhoods significantly increase the number of structure parameters of the system under study,expanding the possibilities of analytical description.To study the structure and properties of such systems in the proposed model,a three-dimensional continuum percolation problem with interacting elements is solved.The dependences of the structure and properties of clusters on the parameters of the generation processes of the cluster system are studied,and analytical dependences are obtained.展开更多
The main goal of this study was to investigate the effects of selected ship collision parameter values on the characteristics of the absorbed energy in several ship collision scenarios. Non-linear simulations were per...The main goal of this study was to investigate the effects of selected ship collision parameter values on the characteristics of the absorbed energy in several ship collision scenarios. Non-linear simulations were performed using a finite element method (FEM) to obtain virtual experiment data. In the present research, the size of the side damage from a collision phenomenon were measured and used to verify the numerical configuration together with the calculation results using an empirical equation. Parameters in the external dynamics of a ship collision such as the location of the contact point and velocity of the striking ship were taken into consideration. The internal energy and deformation size on the side structure were discussed further in a comparative study. The effects of the selected parameters on several structural behaviors, namely energy, force, and damage extent were also observed and evaluated in this section. Stiffener on side hull was found to contribute significantly into resistance capability of the target ship against penetration of the striking bow. Remarkable force during penetration was observed to occur when inner shell was crushed as certain velocity was applied in the striking bow.展开更多
Deep coal mines in northern Anhui province, China, provide opportunities for tracing the distribution and fractionation of rare earth elements CREEs) in deep seated environments. Major ions, as well as REE concentrat...Deep coal mines in northern Anhui province, China, provide opportunities for tracing the distribution and fractionation of rare earth elements CREEs) in deep seated environments. Major ions, as well as REE concentrations were measured in groundwater from a sandstone aquifer located between -400 and -280 m. Our results indicate that this groundwater consists of CI.HCOH-Na or CI.CO3-Na water types with warm temperature (30.1-31.4~C), circumneutral pH (7.27-8.61) and high levels of total dissolved solids (TDS ~ 1306--2165 mg/L). Concentrations of REEs in groundwater are high as expressed by their Nd con- centrations (0.0086-0.018μg/L). Except for weak heavy REEs (HREE) enrichment relative to light REEs (LREE), the similarity of REE distribution patterns between groundwater and aquifer rock indicate that enrichment of REEs is considered to be controlled by aquifer rock, as well as by their minerals, whereas the fractionation of REEs is controlled by HREE enriched minerals and, to a lesser extent, by inorganic REE complexes. Ce anomalies normalized to Post Archean Average Shale (PAAS) and aquifer rock are weak, which probably reflect the contribution of reduced conditions in combination with pH, rather than a sig- nature of aouifer rock.展开更多
Dynamics and vibration of control valves under flow-induced vibration are analyzed. Hydrodynamic load characteristics and structural response under flow-induced vibration are mainly influenced by inertia, damping, ela...Dynamics and vibration of control valves under flow-induced vibration are analyzed. Hydrodynamic load characteristics and structural response under flow-induced vibration are mainly influenced by inertia, damping, elastic, geometric characteristics and hydraulic parameters. The purpose of this work is to investigate the dynamic behavior of control valves in the response to self-excited fluid flow. An analytical and numerical method is developed to simulate the dynamic and vibrational behavior of sliding dam valves, in response to flow excitation. In order to demonstrate the effectiveness of proposed model, the simulation results are validated with experimental ones. Finally, to achieve the optimal valve geometry, numerical results for various shapes of valves are compared. Rounded valve with the least amount of flow turbulence obtains lower fluctuations and vibration amplitude compared with the flat and steep valves. Simulation results demonstrate that with the optimal design requirements of valves, vibration amplitude can be reduced by an average to 30%.展开更多
The Dalnyaya kimberlite pipe(Yakutia,Russia) contains mantle peridotite xenoliths(mostly Iherzolites and harzburgites) that show both sheared porphyroclastic(deformed) and coarse granular textures,together with ...The Dalnyaya kimberlite pipe(Yakutia,Russia) contains mantle peridotite xenoliths(mostly Iherzolites and harzburgites) that show both sheared porphyroclastic(deformed) and coarse granular textures,together with ilmenite and clinopyroxene megacrysts.Deformed peridotites contain high-temperature Fe-rich clinopyroxenes,sometimes associated with picroilmenites,which are products of interaction of the lithospheric mantle with protokimberlite related melts.The orthopyroxene-derived geotherm for the lithospheric mantle beneath Dalnyaya is stepped similar to that beneath the Udachnaya pipe.Coarse granular xenoliths fall on a geotherm of 35 mWm-2 whereas deformed varieties yield a 45 mWm-2)geotherm in the 2-7.5 GPa pressure interval.The chemistry of the constituent minerals including garnet,olivine and clinopyroxene shows trends of increasing Fe~#(=Fe/(Fe+Mg))with decreasing pressure.This may suggest that the interaction with fractionating protokimberlite melts occurred at different levels.Two major mantle lithologies are distinguished by the trace element patterns of their constituent minerals,determined by LA-ICP-MS.Orthopyroxenes,some clinopyroxenes and rare garnets are depleted in Ba,Sr,HFSE and MREE and represent relic lithospheric mantle.Re-fertilized garnet and clinopyroxene are more enriched.The distribution of trace elements between garnet and clinopyroxene shows that the garnets dissolved primary orthopyroxene and clinopyroxene.Later high temperature clinopyroxenes related to the protokimberlite melts partially dissolved these garnets.Olivines show decreases in Ni and increases in Al,Ca and Ti from Mg-rich varieties to the more Fe-rich,deformed and refertilized ones.Minerals showing higher Fe~#(0.11-0.15) are found within intergrowths of low-Cr ilmenite-clinopyroxene-garnet related to the crystallization of protokimberlite melts in feeder channels.In P-f(O_2) diagrams,garnets and Cr-rich clinopyroxenes indicate reduced conditions at the base of the lithosphere at-5 log units below a FMQ buffer.However,Cr-poor clinopyroxenes,together with ilmenite and some Fe-Ca-rich garnets,demonstrate a more oxidized trend in the lower part of lithosphere at-2 to 0 log units relative to FMQ.Clinopyroxenes from xenoliths in most cases show conditions transitional between those determined for garnets and megacrystalline Cr-poor suite.The relatively low diamond grade of Dalnyaya kimberlites is explained by a high degree of interaction with the oxidized protokimberlite melts,which is greater at the base of the lithosphere.展开更多
In the case of composite girders, an effective cooperation of both parts of the section is influenced by deformability of connectors. Limited flexural stiffness of welded studs, used commonly in bridge structures, doe...In the case of composite girders, an effective cooperation of both parts of the section is influenced by deformability of connectors. Limited flexural stiffness of welded studs, used commonly in bridge structures, does not provide full interaction of a steel beam and a concrete slab. This changes strain distribution in cross-sections of a composite girder and results in redistribution of internal forces in steel and concrete element. In the paper partial interaction index defined on the basis of a neutral axis position, which can be used for verification of steel-concrete interaction in real bridge structures rather than in specimens is proposed. The range of the index value changes, obtained during load testing of a typical steel-concrete composite beam bridge, is presented. The investigation was carried out on a motorway viaduct, consisting of two parallel structures. During the testing values of strains in girders under static and quasi-static loads were measured. The readings from the gauges were used to determine the index, characterizing composite action of the girders. Results of bridge testing under movable load, changing position along the bridge span is presented and obtained in-situ influence functions of strains and index values are commented in the paper.展开更多
We report on discrete element method simulations of a pseudo-two-dimensional (pseudo-2D) fluidized bed to investigate particle-wall interactions. Detailed information on macroscopic flow field variables, including s...We report on discrete element method simulations of a pseudo-two-dimensional (pseudo-2D) fluidized bed to investigate particle-wall interactions. Detailed information on macroscopic flow field variables, including solids pressure, granular temperature, and normal and tangential wall stresses are analyzed. The normal wall stress differs from the solids pressure because of the strong anisotropic flow behavior in the pseudo-2D system. A simple linear relationship exists between normal wall stress and solids pressure. In addition, an effective friction coefficient can be derived to characterize particle-wall flow interaction after evaluating the normal and tangential wall stresses. The effects of inter-particle and particle-wall friction coefficients are evaluated. Strong anisotropic flow behavior in the pseudo-2D system needs to be considered to validate the two-fluid model where the boundary condition is usually developed based on an isotropic assumption. The conclusion has been confirmed by simulation with different particle stiffnesses. Assumptions in the newly developed model for 2D simulation are further examined against the discrete element method simulation.展开更多
This paper presents the simulation of tsunamis due to rigid and deformable landslides with consideration of submerged conditions by using particle methods. The smoothed particle hydrodynamics(SPH), as a particle bas...This paper presents the simulation of tsunamis due to rigid and deformable landslides with consideration of submerged conditions by using particle methods. The smoothed particle hydrodynamics(SPH), as a particle based method, is for solving problems of fast moving boundaries in the field of continuum mechanics. Other particle based methods, like the discrete element method(DEM), are suitable for modeling the displacement and the collision related to the rigid landslides. In the present work, we use the SPH and the DEM to simulate tsunamis generated by rigid and deformable landslides with consideration of submerged conditions. The viscous free-surface flows are solved by a weakly compressible SPH and the displacement and the rotation of the rigid body slides are calculated using a multi-sphere DEM allowing for modeling solids of arbitrarily complex shapes. The fluid-solid interactions are simulated by coupling the SPH and the DEM. A rheology model combining the Papanastasiou and the Herschel-Bulkley models is applied to represent the viscoplastic behavior of the non-Newtonian flow in the submarine deformable landslide cases. Submarine landslide tsunamis due to rigid and deformable landslides are both simulated as typical landslide cases in this investigation. Our simulated results and the previous experimental results in the literatures are in good agreement, which shows that the proposed particle based methods are capable of modeling the submarine landslide tsunamis.展开更多
A two-dimensional coupled lattice Boltzmann immersed boundary discrete element method is introduced for the simulation of polygonal particles moving in incompressible viscous fluids. A collision model of polygonal par...A two-dimensional coupled lattice Boltzmann immersed boundary discrete element method is introduced for the simulation of polygonal particles moving in incompressible viscous fluids. A collision model of polygonal particles is used in the discrete element method. Instead of a collision model of circular particles, the collision model used in our method can deal with particles of more complex shape and efficiently simulate the effects of shape on particle–particle and particle–wall interactions. For two particles falling under gravity, because of the edges and corners, different collision patterns for circular and polygonal particles are found in our simulations. The complex vortexes generated near the corners of polygonal particles affect the flow field and lead to a difference in particle motions between circular and polygonal particles. For multiple particles falling under gravity, the polygonal particles easily become stuck owing to their corners and edges, while circular particles slip along contact areas. The present method provides an efficient approach for understanding the effects of particle shape on the dynamics of non-circular particles in fluids.展开更多
The boundary element method(BEM) is a main method for analyzing the interactions between the waves and the marine structures. As with the BEM, a set of linear equations are generated with a full matrix, the required...The boundary element method(BEM) is a main method for analyzing the interactions between the waves and the marine structures. As with the BEM, a set of linear equations are generated with a full matrix, the required calculations and storage increase rapidly with the increase of the structure scale. Thus, an accelerated method with a low storage is desirable for the wave interaction with a very large structure. A systematic review is given in this paper for the BEM for solving the problem of the wave interaction with a large scale structure. Various integral equations are derived based on different Green functions, the advantages and disadvantages of different discretization schemes of the integral equations by the constant panels, the higher order elements, and the spline functions are discussed. For the higher order element discretization method, the special concerns are given to the numerical calculations of the single-layer potential, the double layer potential and the solid angle coefficients. For a large scale computation problem such as the wave interaction with a very large structure or a large number of bodies, the BEMs with the FMM and p FFT accelerations are discussed, respectively, including the principles of the FMM and the p FFT, and their implementations in various integral equations with different Green functions. Finally, some potential applications of the acceleration methods for problems with large scale computations in the ocean and coastal engineering are introduced.展开更多
文摘Short pitch corrugation has been a problem for railways worldwide over one century.In this paper,a parametric investigation of fastenings is conducted to understand the corrugation formation mechanism and gain insights into corrugation mitigation.A three-dimensional finite element vehicle-track dynamic interaction model is employed,which considers the coupling between the structural dynamics and the contact mechanics,while the damage mechanism is assumed to be differential wear.Various fastening models with different configurations,boundary conditions,and parameters of stiffness and damping are built up and analysed.These models may represent different service stages of fastenings in the field.Besides,the effect of train speeds on corrugation features is studied.The results indicate:(1)Fastening parameters and modelling play an important role in corrugation formation.(2)The fastening longitudinal constraint to the rail is the major factor that determines the corrugation formation.The fastening vertical and lateral constraints influence corrugation features in terms of spatial distribution and wavelength components.(3)The strengthening of fastening constraints in the longitudinal dimension helps to mitigate corrugation.Meanwhile,the inner fastening constraint in the lateral direction is necessary for corrugation alleviation.(4)The increase in fastening longitudinal stiffness and damping can reduce the vibration amplitudes of longitudinal compression modes and thus reduce the track corrugation propensity.The simulation in this work can well explain the field corrugation in terms of the occurrence possibility and major wavelength components.It can also explain the field data with respect to the small variation between the corrugation wavelength and train speed,which is caused by frequency selection and jump between rail longitudinal compression modes.
文摘The model of dislocations was used to construct the model of grain boundary (GB) with pure rare earths, and rare earth elements and impurities. The influence of the interaction between rare earth elements and impurities on the cohesive properties of 5.3° low angle GB of Fe was investigated by the recursion method. The calculated results of environment sensitive embeding energy( E ESE ) show that the preferential segregation of rare earth elements towards GBs exists. Calculations of bond order integrals (BOI) show that rare earth elements increase the cohesive strength of low angle GB, and impurities such as S, P weaken the intergranular cohesion of the GB. So rare earth element of proper quantity added in steel not only cleanses other harmful impurities off the GBs, but also enhances the intergranular cohesion. This elucidates the action mechanism of rare earth elements in steel from electronic level and offers theoretical evidence for applications of rare earth elements in steels.
文摘This paper uses a mathematical method to develop an analytical solution to the local buckling behaviour of long rectangular plates resting on tensionless elastic Winkler foundations and under combined uniform longitudinal uniaxial compressive and uniform in-plane shear loads. Fitted formulas are derived for plates with clamped edges and simplified supported edges. Two examples are given to demonstrate the application of the current method: one is a plate on tensionless spring foundations and the other is the contact between the steel sheet and elastic solid foundation. Finite element (FE) analysis is also conducted to validate the analytical results. Good agreement is obtained between the current method and FE analysis.
文摘The percolation fields constructed around the elements of a cluster system in the phase spaces of properties are studied.It is shown that such neighborhoods significantly increase the number of structure parameters of the system under study,expanding the possibilities of analytical description.To study the structure and properties of such systems in the proposed model,a three-dimensional continuum percolation problem with interacting elements is solved.The dependences of the structure and properties of clusters on the parameters of the generation processes of the cluster system are studied,and analytical dependences are obtained.
文摘The main goal of this study was to investigate the effects of selected ship collision parameter values on the characteristics of the absorbed energy in several ship collision scenarios. Non-linear simulations were performed using a finite element method (FEM) to obtain virtual experiment data. In the present research, the size of the side damage from a collision phenomenon were measured and used to verify the numerical configuration together with the calculation results using an empirical equation. Parameters in the external dynamics of a ship collision such as the location of the contact point and velocity of the striking ship were taken into consideration. The internal energy and deformation size on the side structure were discussed further in a comparative study. The effects of the selected parameters on several structural behaviors, namely energy, force, and damage extent were also observed and evaluated in this section. Stiffener on side hull was found to contribute significantly into resistance capability of the target ship against penetration of the striking bow. Remarkable force during penetration was observed to occur when inner shell was crushed as certain velocity was applied in the striking bow.
基金financially supported by the National Natural Science Foundation of China (No. 40873015)the Eleventh Five-Year Science and Technology Project of Anhui Province,China(No. 08010302062)
文摘Deep coal mines in northern Anhui province, China, provide opportunities for tracing the distribution and fractionation of rare earth elements CREEs) in deep seated environments. Major ions, as well as REE concentrations were measured in groundwater from a sandstone aquifer located between -400 and -280 m. Our results indicate that this groundwater consists of CI.HCOH-Na or CI.CO3-Na water types with warm temperature (30.1-31.4~C), circumneutral pH (7.27-8.61) and high levels of total dissolved solids (TDS ~ 1306--2165 mg/L). Concentrations of REEs in groundwater are high as expressed by their Nd con- centrations (0.0086-0.018μg/L). Except for weak heavy REEs (HREE) enrichment relative to light REEs (LREE), the similarity of REE distribution patterns between groundwater and aquifer rock indicate that enrichment of REEs is considered to be controlled by aquifer rock, as well as by their minerals, whereas the fractionation of REEs is controlled by HREE enriched minerals and, to a lesser extent, by inorganic REE complexes. Ce anomalies normalized to Post Archean Average Shale (PAAS) and aquifer rock are weak, which probably reflect the contribution of reduced conditions in combination with pH, rather than a sig- nature of aouifer rock.
文摘Dynamics and vibration of control valves under flow-induced vibration are analyzed. Hydrodynamic load characteristics and structural response under flow-induced vibration are mainly influenced by inertia, damping, elastic, geometric characteristics and hydraulic parameters. The purpose of this work is to investigate the dynamic behavior of control valves in the response to self-excited fluid flow. An analytical and numerical method is developed to simulate the dynamic and vibrational behavior of sliding dam valves, in response to flow excitation. In order to demonstrate the effectiveness of proposed model, the simulation results are validated with experimental ones. Finally, to achieve the optimal valve geometry, numerical results for various shapes of valves are compared. Rounded valve with the least amount of flow turbulence obtains lower fluctuations and vibration amplitude compared with the flat and steep valves. Simulation results demonstrate that with the optimal design requirements of valves, vibration amplitude can be reduced by an average to 30%.
基金supported by RBRF grants:05-05-64718,11-0500060,11-05-91060-PICS,16-05-00860the projects 77-2,65-03,02-05 UIGGM SD RAS and ALROSA Stock Company
文摘The Dalnyaya kimberlite pipe(Yakutia,Russia) contains mantle peridotite xenoliths(mostly Iherzolites and harzburgites) that show both sheared porphyroclastic(deformed) and coarse granular textures,together with ilmenite and clinopyroxene megacrysts.Deformed peridotites contain high-temperature Fe-rich clinopyroxenes,sometimes associated with picroilmenites,which are products of interaction of the lithospheric mantle with protokimberlite related melts.The orthopyroxene-derived geotherm for the lithospheric mantle beneath Dalnyaya is stepped similar to that beneath the Udachnaya pipe.Coarse granular xenoliths fall on a geotherm of 35 mWm-2 whereas deformed varieties yield a 45 mWm-2)geotherm in the 2-7.5 GPa pressure interval.The chemistry of the constituent minerals including garnet,olivine and clinopyroxene shows trends of increasing Fe~#(=Fe/(Fe+Mg))with decreasing pressure.This may suggest that the interaction with fractionating protokimberlite melts occurred at different levels.Two major mantle lithologies are distinguished by the trace element patterns of their constituent minerals,determined by LA-ICP-MS.Orthopyroxenes,some clinopyroxenes and rare garnets are depleted in Ba,Sr,HFSE and MREE and represent relic lithospheric mantle.Re-fertilized garnet and clinopyroxene are more enriched.The distribution of trace elements between garnet and clinopyroxene shows that the garnets dissolved primary orthopyroxene and clinopyroxene.Later high temperature clinopyroxenes related to the protokimberlite melts partially dissolved these garnets.Olivines show decreases in Ni and increases in Al,Ca and Ti from Mg-rich varieties to the more Fe-rich,deformed and refertilized ones.Minerals showing higher Fe~#(0.11-0.15) are found within intergrowths of low-Cr ilmenite-clinopyroxene-garnet related to the crystallization of protokimberlite melts in feeder channels.In P-f(O_2) diagrams,garnets and Cr-rich clinopyroxenes indicate reduced conditions at the base of the lithosphere at-5 log units below a FMQ buffer.However,Cr-poor clinopyroxenes,together with ilmenite and some Fe-Ca-rich garnets,demonstrate a more oxidized trend in the lower part of lithosphere at-2 to 0 log units relative to FMQ.Clinopyroxenes from xenoliths in most cases show conditions transitional between those determined for garnets and megacrystalline Cr-poor suite.The relatively low diamond grade of Dalnyaya kimberlites is explained by a high degree of interaction with the oxidized protokimberlite melts,which is greater at the base of the lithosphere.
文摘In the case of composite girders, an effective cooperation of both parts of the section is influenced by deformability of connectors. Limited flexural stiffness of welded studs, used commonly in bridge structures, does not provide full interaction of a steel beam and a concrete slab. This changes strain distribution in cross-sections of a composite girder and results in redistribution of internal forces in steel and concrete element. In the paper partial interaction index defined on the basis of a neutral axis position, which can be used for verification of steel-concrete interaction in real bridge structures rather than in specimens is proposed. The range of the index value changes, obtained during load testing of a typical steel-concrete composite beam bridge, is presented. The investigation was carried out on a motorway viaduct, consisting of two parallel structures. During the testing values of strains in girders under static and quasi-static loads were measured. The readings from the gauges were used to determine the index, characterizing composite action of the girders. Results of bridge testing under movable load, changing position along the bridge span is presented and obtained in-situ influence functions of strains and index values are commented in the paper.
文摘We report on discrete element method simulations of a pseudo-two-dimensional (pseudo-2D) fluidized bed to investigate particle-wall interactions. Detailed information on macroscopic flow field variables, including solids pressure, granular temperature, and normal and tangential wall stresses are analyzed. The normal wall stress differs from the solids pressure because of the strong anisotropic flow behavior in the pseudo-2D system. A simple linear relationship exists between normal wall stress and solids pressure. In addition, an effective friction coefficient can be derived to characterize particle-wall flow interaction after evaluating the normal and tangential wall stresses. The effects of inter-particle and particle-wall friction coefficients are evaluated. Strong anisotropic flow behavior in the pseudo-2D system needs to be considered to validate the two-fluid model where the boundary condition is usually developed based on an isotropic assumption. The conclusion has been confirmed by simulation with different particle stiffnesses. Assumptions in the newly developed model for 2D simulation are further examined against the discrete element method simulation.
基金Project supported by the National Natural Science Foun-dation of China(Grant Nos.11172321,51509248)the Scientific Research and Experiment of Regulation Engineering for the Songhua River Mainstream in Heilongjiang Province(Grant No.SGZL/KY-12)
文摘This paper presents the simulation of tsunamis due to rigid and deformable landslides with consideration of submerged conditions by using particle methods. The smoothed particle hydrodynamics(SPH), as a particle based method, is for solving problems of fast moving boundaries in the field of continuum mechanics. Other particle based methods, like the discrete element method(DEM), are suitable for modeling the displacement and the collision related to the rigid landslides. In the present work, we use the SPH and the DEM to simulate tsunamis generated by rigid and deformable landslides with consideration of submerged conditions. The viscous free-surface flows are solved by a weakly compressible SPH and the displacement and the rotation of the rigid body slides are calculated using a multi-sphere DEM allowing for modeling solids of arbitrarily complex shapes. The fluid-solid interactions are simulated by coupling the SPH and the DEM. A rheology model combining the Papanastasiou and the Herschel-Bulkley models is applied to represent the viscoplastic behavior of the non-Newtonian flow in the submarine deformable landslide cases. Submarine landslide tsunamis due to rigid and deformable landslides are both simulated as typical landslide cases in this investigation. Our simulated results and the previous experimental results in the literatures are in good agreement, which shows that the proposed particle based methods are capable of modeling the submarine landslide tsunamis.
基金This study was funded by the National Science Foundation of China (Grant No. 11272176).
文摘A two-dimensional coupled lattice Boltzmann immersed boundary discrete element method is introduced for the simulation of polygonal particles moving in incompressible viscous fluids. A collision model of polygonal particles is used in the discrete element method. Instead of a collision model of circular particles, the collision model used in our method can deal with particles of more complex shape and efficiently simulate the effects of shape on particle–particle and particle–wall interactions. For two particles falling under gravity, because of the edges and corners, different collision patterns for circular and polygonal particles are found in our simulations. The complex vortexes generated near the corners of polygonal particles affect the flow field and lead to a difference in particle motions between circular and polygonal particles. For multiple particles falling under gravity, the polygonal particles easily become stuck owing to their corners and edges, while circular particles slip along contact areas. The present method provides an efficient approach for understanding the effects of particle shape on the dynamics of non-circular particles in fluids.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.51379032,51490672 and 51479026)
文摘The boundary element method(BEM) is a main method for analyzing the interactions between the waves and the marine structures. As with the BEM, a set of linear equations are generated with a full matrix, the required calculations and storage increase rapidly with the increase of the structure scale. Thus, an accelerated method with a low storage is desirable for the wave interaction with a very large structure. A systematic review is given in this paper for the BEM for solving the problem of the wave interaction with a large scale structure. Various integral equations are derived based on different Green functions, the advantages and disadvantages of different discretization schemes of the integral equations by the constant panels, the higher order elements, and the spline functions are discussed. For the higher order element discretization method, the special concerns are given to the numerical calculations of the single-layer potential, the double layer potential and the solid angle coefficients. For a large scale computation problem such as the wave interaction with a very large structure or a large number of bodies, the BEMs with the FMM and p FFT accelerations are discussed, respectively, including the principles of the FMM and the p FFT, and their implementations in various integral equations with different Green functions. Finally, some potential applications of the acceleration methods for problems with large scale computations in the ocean and coastal engineering are introduced.