A simplified approach was proposed to analyze the negative skin friction calculation of special-shaped pile considering pile-soil interaction under surcharge. Based on the concentric cylinder shearing theory, consider...A simplified approach was proposed to analyze the negative skin friction calculation of special-shaped pile considering pile-soil interaction under surcharge. Based on the concentric cylinder shearing theory, considering the changes of pile shape(such as, taper angle and diameters of pile base, etc.), the load-transfer of special-shaped pile was built. The accuracy of the developed simplified approach was verified by numerical simulation model with the same condition. Then, the influence factors, such as, taper angles, the diameter of pile base, surcharge, and pile-soil interface parameters were analyzed and discussed. The results show that the developed simplified approach can calculate NSF of special-shaped pile under surcharge effectively. A limited parametric study indicates that in many practical situations special-shaped piles(such as belled wedge pile shown in this work) offer a design option that is more economical than traditional uniform cross-section piles.展开更多
Robust numerical models that describe the complex behaviors of risers are needed because these constitute dynamically sensitive systems. This paper presents a simple and efficient algorithm for the nonlinear static an...Robust numerical models that describe the complex behaviors of risers are needed because these constitute dynamically sensitive systems. This paper presents a simple and efficient algorithm for the nonlinear static and dynamic analyses of marine risers. The proposed approach uses the vector form intrinsic finite element(VFIFE) method, which is based on vector mechanics theory and numerical calculation. In this method, the risers are described by a set of particles directly governed by Newton's second law and are connected by weightless elements that can only resist internal forces. The method does not require the integration of the stiffness matrix, nor does it need iterations to solve the governing equations. Due to these advantages, the method can easily increase or decrease the element and change the boundary conditions, thus representing an innovative concept of solving nonlinear behaviors, such as large deformation and large displacement. To prove the feasibility of the VFIFE method in the analysis of the risers, rigid and flexible risers belonging to two different categories of marine risers, which usually have differences in modeling and solving methods, are employed in the present study. In the analysis, the plane beam element is adopted in the simulation of interaction forces between the particles and the axial force, shear force, and bending moment are also considered. The results are compared with the conventional finite element method(FEM) and those reported in the related literature. The findings revealed that both the rigid and flexible risers could be modeled in a similar unified analysis model and that the VFIFE method is feasible for solving problems related to the complex behaviors of marine risers.展开更多
A new method based on the multi-wedge translation mechanism is presented to calculate the lateral force acting on the stabilizing piles. At first, there is no assumption for the shape of potential sliding surface, it ...A new method based on the multi-wedge translation mechanism is presented to calculate the lateral force acting on the stabilizing piles. At first, there is no assumption for the shape of potential sliding surface, it is just considered that the potential sliding surface is a composite of a number of straight lines. And then, the potential sliding mass is divided into a number of triangular wedges take with these straight lines as its base. The kinematic theorem of limit analysis is adopted to calculate the rate of external work and the rate of energy dissipation for each triangular wedge, respectively. Furthermore, the multivariate functions are established to calculate the lateral force acting on the stabilizing piles. The lateral force and the corresponding potential sliding surfaces can be obtained by an optimizational technique. At last, an example is taken to illustrate the method. The effect of soil strength parameters, slope angle and pile roughness on the lateral force and the corresponding potential sliding surface are analyzed.The result are compared with those obtained using other methods.展开更多
The state of deep mine rigid shaft fUrniture horizontal force calculation in China is introduced and the calculating model of horizontal force is determined. Based on the interaction between hoisting vessel and shaft ...The state of deep mine rigid shaft fUrniture horizontal force calculation in China is introduced and the calculating model of horizontal force is determined. Based on the interaction between hoisting vessel and shaft furniture, the vertical mobel test device is built and a large number of model tests and analyses are carried out. At last, the relationship between horizontal force and hoisting end load, hoisting speed or bunton level interval is gained. The three parameters horizontal force calculating equation of beep mine rigid shaft furniture is given and it can guide the design and calculation of shaft engineering.展开更多
The reaction of CH2SH radical with fluorine atom was studied at the levels of B3LYP/6-311G(d,p) and MP2(Full)/6-311G(d,p). The computational results show that the reaction has three channels and proceeds by the ...The reaction of CH2SH radical with fluorine atom was studied at the levels of B3LYP/6-311G(d,p) and MP2(Full)/6-311G(d,p). The computational results show that the reaction has three channels and proceeds by the addition of fluorine atoms on carbon or sulfur sites of CH2SH, forming initial intermediates. The calculated results show that the channel, in which fluorine attaches to the carbon atom to form CH2S and HF, is the most likely reaction pathway. Topological analysis of electron density was carried out for the three channels. The change trends of the chemical bonds on the reaction paths were discussed. The energy transition states and the structure transition regions (states) of the three channels were found. The calculated results show that the structure transition regions are broad in unobvious exothermic reactions or unobvious endothermic reactions, and are narrow in obvious exothermic reactions or obvious endothermic reactions.展开更多
A computer code based on the double-body potential flow model and the classic source panel method has been developed to study various problems of hydrodynamic interaction between ships and other objects with solid bou...A computer code based on the double-body potential flow model and the classic source panel method has been developed to study various problems of hydrodynamic interaction between ships and other objects with solid boundaries including the seabed. A peculiarity of the proposed implementation is the application of the so-called "moving-patch" method for simulating steady boundaries of large extensions. The method is based on an assumption that at any moment just the part of the boundary ("moving patch") which lies close to the interacting ship is significant for the near-field interaction. For a specific case of the fiat bottom, comparative computations were performed to determine optimal dimensions of the patch and of the constituting panels based on the trade-off between acceptable accuracy and reasonable efficiency. The method was applied to estimate the sway force on a ship hull moving obliquely across a dredged channel. The method was validated for a case of ship-to-ship interaction when tank data were available. This study also contains a description of a newly developed spline approximation algorithm necessary for creating consistent discretizations of ship hulls with various degrees of refinement.展开更多
Although several experimental and numerical studies have been carried out on the shear capacity of RC beams retrofitted by carbon or glass fibre-reinforced polymers, there has been little work on hybrid FRP sheet appl...Although several experimental and numerical studies have been carried out on the shear capacity of RC beams retrofitted by carbon or glass fibre-reinforced polymers, there has been little work on hybrid FRP sheet applications, particularly under cyclic loading. In the present research, five RC beams were constructed, and four of which were retrofitted using various schemes of FRP sheets. All beams were subjected to quasi-static cyclic loading in an attempt to represent the effect repetitive loading. The ultimate load, and deflection response at mid-span of the beams were measured and compared with predictions of a computational model based on finite element analysis. Experimental results demonstrated that hybrid applications of FRP sheets can improve the shear performance of retrofitted RC beams and increase the ultimate strain of the FRP sheets at failure. The results of the computational model were in reasonable agreement with the corresponding experimental results.展开更多
At present,the thrust of an anti-slide pile can be worked out with some calculation methods. However,the resistance in front of the pile,the distributions of resistance and thrust,and appropriate pile length cannot be...At present,the thrust of an anti-slide pile can be worked out with some calculation methods. However,the resistance in front of the pile,the distributions of resistance and thrust,and appropriate pile length cannot be easily obtained. In this paper,the authors applied the strength-reduction finite element method (FEM) to several design cases of anti-slide piles. Using this method,it is possible to take the pile-soil interactions into consideration,obtain reasonable resistance in front of pile and the distributions of thrust and resistance,and reasonable lengths of anti-slide piles. In particular,the thrust and resistance imposed on embedded anti-slide piles can be calculated and composite anti-slide pile structures such as anchored piles and braced piles can be optimized. It is proved through the calculation examples that this method is more reliable and economical in the design of anti-slide pile.展开更多
New global three dimensional potential energy surfaces for the Cl+H2 reactive system have been constructed using accurate multireference configuration interaction calculations with a large basis set. The three lowest...New global three dimensional potential energy surfaces for the Cl+H2 reactive system have been constructed using accurate multireference configuration interaction calculations with a large basis set. The three lowest adiabatic potential energy surfaces correlating asymptotically with Cl(^2p)+H2 have been transformed to adiabatic representation, which leads to a fourth coupling potential for non-linear geometries. In addition, the spin-orbit coupling surfaces have also been computed using the Breit-Pauli Hamiltonian. Properties of the new potential are described. Reaction dynamics based on the new potential agrees with the recent experimental results quite well.展开更多
In modem missile design, the operation of a missile aerodynamics with angles of attack is required to serve a demand on the maneuverability. The key aero-physics is the development of vortices and its interaction to t...In modem missile design, the operation of a missile aerodynamics with angles of attack is required to serve a demand on the maneuverability. The key aero-physics is the development of vortices and its interaction to the control surface such as wing and fins. This paper thus presents the investigation of the missile flow field at 4° and 8° degrees of angles of attack. The Mach numbers for both case were varied from 0.6 to 5.5. Here, the Steady Reynolds-Averaged Navier-Stokes (SRANS) equations with standard κ-ε turbulence model were selected. The numerical results of aerodynamics coefficients (both force and moment) were compared against semi-empirical data computed using Missile DatCOM. The results revealed the development of vortices observed and their interaction with fin at the rear part of the missile.展开更多
To make sure the intermediate products of the carbothermic reduction of Al2O3 process, such as Al4CO4, Al2CO and Al4C3,and the interaction of AlCl, AlCl2, AlCl3 with Al4CO4, Al2CO and Al4C3, respectively, thermodynami...To make sure the intermediate products of the carbothermic reduction of Al2O3 process, such as Al4CO4, Al2CO and Al4C3,and the interaction of AlCl, AlCl2, AlCl3 with Al4CO4, Al2CO and Al4C3, respectively, thermodynamic analyses were used to study the chloride reaction production of them under the vacuum situation. The stable structures and electronic properties of AlCl, AlCl2 and AlCl3 adsorbed on Al4CO4, Al2CO and Al4C3 were calculated by first-principles calculations by the CASTEP module in the Materials Studio program. The results show that the AlCl3 and AlCl2 molecules have decomposed on the plane of Al4C3(001), while there are no obvious decomposition of AlCl3 and AlCl2 on Al4CO4(001) and Al2CO(001) planes. The adsorption of AlCl on the Al4CO4(001)and Al2CO(001) planes is stronger than that on the Al4C3(001) plane. The interaction strength of AlCl3, as well as AlCl2, with Al4CO4,Al2CO and Al4C3 is in the sequence of Al4CO4Al2CO〉Al4C3.展开更多
A parallelized resolved method for the simulation of the dynamics of immersed bodies within fluids is presented. The algorithm uses a FDM (fictitious domain method) and combines the Lagrangian DEM (discrete element...A parallelized resolved method for the simulation of the dynamics of immersed bodies within fluids is presented. The algorithm uses a FDM (fictitious domain method) and combines the Lagrangian DEM (discrete element method) for tracking the bodies with a CFD (computational fluid dynamics) method for calculating the dynamics of the fluid phase. First the CFD-calculation is carried out, disregarding the solid bodies. Afterwards, the velocity information from the bodies is included and the force, the fluid imposes onto the bodies, is computed. The last step consists of a correction-operation which ensures the fulfillment of the conservation equation. Dynamic local mesh refinement is used for minimizing the number of fluid cells. The CFD-DEM coupling is realized within the Open Source framework CFDEMcoupling (www.cfdem.com), where the DEM software LIGGGHTS (www.liggghts.com) is linked against an OpenFOAM^-based CFD solver. While both LIGGGHTS and the CFD solver were already parallelized, only a recent improvement of the algorithm permits the fully parallel computation of resolved problems. This parallelization permits the treatment of large-scale problems. The enclosed validation and application examples show the dynamics of the flow around settling and rotating spheres as well as an investigation of the settling of spheres regarding the Boycott effect.展开更多
The objective of this paper is to improve the understanding of the influence of multiphase flow on the turbulent closure model, the interplay between vorticity fields and cavity dynamics around a pitching hydrofoil. T...The objective of this paper is to improve the understanding of the influence of multiphase flow on the turbulent closure model, the interplay between vorticity fields and cavity dynamics around a pitching hydrofoil. The effects of pitching rate on the sub- cavitating and cavitating response of the pitching hydrofoil are also investigated. In particular, we focus on the interactions between cavity inception, growth, and shedding and the vortex flow structures, and their impacts on the hydrofoil performance. The calculations are 2-D and performed by solving the incompressible, multiphase Unsteady Reynolds Averaged Navier Stokes (URANS) equations via the commercial CFD code CFX. The k-co SST (Shear Stress Transport) turbulence model is used along with the transport equation-based cavitation models. The density correction function is considered to reduce the eddy viscosity according to the computed local fluid mixture density. The calculation results are validated with experiments conducted by Ducoin et al. (see Computational and experimental investigation of flow over a transient pitching hydrofoil, Eur J Mech/B Fluids, 2009, 28:728-743 and An experimental analysis of fluid structure interaction of a flexible hydrofoil in vari- ous flow regimes including cavitating flow, Eur J Mech B/fluids, 2012, 36: 63-74). Results are shown for a NACA66 hydro- foil subject to slow (quasi static, t2=6~/s, &* =0.18) and fast (dynamic, &=63~/s, dr" =1.89) pitching motions from a =0~ to a =15~. Both subcavitaing (or =8.0) and cavitating (cr=3.0) flows are considered. For subcavitating flow (or=8.0), low frequency fluctuations have been observed when the leading edge vortex shedding occurs during stall, and delay of stall is ob- served with increasing pitching velocity. For cavitating flow (tr=3.0), small leading edge cavities are observed with the slow pitching case, which significantly modified the vortex dynamics at high angles of attack, leading to high frequency fluctuations of the hydrodynamic coefficients and different stall behaviors compared to the subcavitating flow at the same pitching rate. On the other hand, for the fast pitching case at or=3.0, large-scale sheet/cloud cavitation is observed, the cavity behavior is un- steady and has a strong impact on the hydrodynamic response, which leads to high amplitude fluctuations of the hydrodynamic coefficients, as well as significant changes in the stall and post-stall behavior. The numerical results also show that the local density modification helps to reduce turbulent eddy viscosity in the cavitating region, which significantly modifies the cavity lengths and shedding frequencies, particularly for the fast pitching case. In general, compared with the experimental visualiza- tions, the numerical results with local density correction have been found to agree well with experimental measurements and observations for both slow and fast transient pitching cases.展开更多
Our objective is a better understanding of the role of physical properties of real fluids in the thermodynamics of cavitation in impure water. An extension to the classical homogenous nucleation theory suitable for mi...Our objective is a better understanding of the role of physical properties of real fluids in the thermodynamics of cavitation in impure water. An extension to the classical homogenous nucleation theory suitable for mixtures is presented in attempt to address the discrepancy between the theoretical predictions and practical observations of cavitation rates in water at normal temperatures. The extension takes into account the non-equilibrium (dissipative) effects involved in nuclei formation through a substance dependent correction coefficient to be determined experimentally. The theory of thermodynamic fluctuations is applied to derive the work of formation of a bubble nucleus. The value of the correction coefficient is estimated using preliminary experimental data from a convergent-divergent nozzle. An application of the results to the numerical prediction of the cavitation zones in a radial-flow water pump is shown.展开更多
基金Projects(51278170,U1134207)supported by the National Natural Science Foundation of ChinaProject(B13024)supported by the Program of Introducing Talents of Discipline to Universities,China
文摘A simplified approach was proposed to analyze the negative skin friction calculation of special-shaped pile considering pile-soil interaction under surcharge. Based on the concentric cylinder shearing theory, considering the changes of pile shape(such as, taper angle and diameters of pile base, etc.), the load-transfer of special-shaped pile was built. The accuracy of the developed simplified approach was verified by numerical simulation model with the same condition. Then, the influence factors, such as, taper angles, the diameter of pile base, surcharge, and pile-soil interface parameters were analyzed and discussed. The results show that the developed simplified approach can calculate NSF of special-shaped pile under surcharge effectively. A limited parametric study indicates that in many practical situations special-shaped piles(such as belled wedge pile shown in this work) offer a design option that is more economical than traditional uniform cross-section piles.
基金supported by the National Key Research and Development Program (No. 2016YFC0802301)the Shandong Province Science and Technology Major Project (No. 2015ZDZX04003)the Natural Science Foundation of Shandong Province (No. ZR2016GM06)
文摘Robust numerical models that describe the complex behaviors of risers are needed because these constitute dynamically sensitive systems. This paper presents a simple and efficient algorithm for the nonlinear static and dynamic analyses of marine risers. The proposed approach uses the vector form intrinsic finite element(VFIFE) method, which is based on vector mechanics theory and numerical calculation. In this method, the risers are described by a set of particles directly governed by Newton's second law and are connected by weightless elements that can only resist internal forces. The method does not require the integration of the stiffness matrix, nor does it need iterations to solve the governing equations. Due to these advantages, the method can easily increase or decrease the element and change the boundary conditions, thus representing an innovative concept of solving nonlinear behaviors, such as large deformation and large displacement. To prove the feasibility of the VFIFE method in the analysis of the risers, rigid and flexible risers belonging to two different categories of marine risers, which usually have differences in modeling and solving methods, are employed in the present study. In the analysis, the plane beam element is adopted in the simulation of interaction forces between the particles and the axial force, shear force, and bending moment are also considered. The results are compared with the conventional finite element method(FEM) and those reported in the related literature. The findings revealed that both the rigid and flexible risers could be modeled in a similar unified analysis model and that the VFIFE method is feasible for solving problems related to the complex behaviors of marine risers.
基金Projects(SKLGP2012K024,SKLGP2013K012)supported by the Opening Fund of State Key Laboratory of Geohazard Prevention and Ceoenvironment Protection,ChinaProject(2011BAK12B03)supported by the National Technology Project,ChinaProject(41401004)supported by the National Natural Science Foundation of China
文摘A new method based on the multi-wedge translation mechanism is presented to calculate the lateral force acting on the stabilizing piles. At first, there is no assumption for the shape of potential sliding surface, it is just considered that the potential sliding surface is a composite of a number of straight lines. And then, the potential sliding mass is divided into a number of triangular wedges take with these straight lines as its base. The kinematic theorem of limit analysis is adopted to calculate the rate of external work and the rate of energy dissipation for each triangular wedge, respectively. Furthermore, the multivariate functions are established to calculate the lateral force acting on the stabilizing piles. The lateral force and the corresponding potential sliding surfaces can be obtained by an optimizational technique. At last, an example is taken to illustrate the method. The effect of soil strength parameters, slope angle and pile roughness on the lateral force and the corresponding potential sliding surface are analyzed.The result are compared with those obtained using other methods.
文摘The state of deep mine rigid shaft fUrniture horizontal force calculation in China is introduced and the calculating model of horizontal force is determined. Based on the interaction between hoisting vessel and shaft furniture, the vertical mobel test device is built and a large number of model tests and analyses are carried out. At last, the relationship between horizontal force and hoisting end load, hoisting speed or bunton level interval is gained. The three parameters horizontal force calculating equation of beep mine rigid shaft furniture is given and it can guide the design and calculation of shaft engineering.
基金This work was supported by the National Natural Science Foundation of China (No.20573032),the Natural Science Foundation of Hebei Province (No.B2006000137), the Education Department Foundation of Hebei Province (No.2007123), and the Doctoral Foundation of Hebei Normal University (No.L2005B12).
文摘The reaction of CH2SH radical with fluorine atom was studied at the levels of B3LYP/6-311G(d,p) and MP2(Full)/6-311G(d,p). The computational results show that the reaction has three channels and proceeds by the addition of fluorine atoms on carbon or sulfur sites of CH2SH, forming initial intermediates. The calculated results show that the channel, in which fluorine attaches to the carbon atom to form CH2S and HF, is the most likely reaction pathway. Topological analysis of electron density was carried out for the three channels. The change trends of the chemical bonds on the reaction paths were discussed. The energy transition states and the structure transition regions (states) of the three channels were found. The calculated results show that the structure transition regions are broad in unobvious exothermic reactions or unobvious endothermic reactions, and are narrow in obvious exothermic reactions or obvious endothermic reactions.
基金Supported by the Portuguese Foundation for Science and Technology under Grant No.PTDC/ECM/100686/2008
文摘A computer code based on the double-body potential flow model and the classic source panel method has been developed to study various problems of hydrodynamic interaction between ships and other objects with solid boundaries including the seabed. A peculiarity of the proposed implementation is the application of the so-called "moving-patch" method for simulating steady boundaries of large extensions. The method is based on an assumption that at any moment just the part of the boundary ("moving patch") which lies close to the interacting ship is significant for the near-field interaction. For a specific case of the fiat bottom, comparative computations were performed to determine optimal dimensions of the patch and of the constituting panels based on the trade-off between acceptable accuracy and reasonable efficiency. The method was applied to estimate the sway force on a ship hull moving obliquely across a dredged channel. The method was validated for a case of ship-to-ship interaction when tank data were available. This study also contains a description of a newly developed spline approximation algorithm necessary for creating consistent discretizations of ship hulls with various degrees of refinement.
文摘Although several experimental and numerical studies have been carried out on the shear capacity of RC beams retrofitted by carbon or glass fibre-reinforced polymers, there has been little work on hybrid FRP sheet applications, particularly under cyclic loading. In the present research, five RC beams were constructed, and four of which were retrofitted using various schemes of FRP sheets. All beams were subjected to quasi-static cyclic loading in an attempt to represent the effect repetitive loading. The ultimate load, and deflection response at mid-span of the beams were measured and compared with predictions of a computational model based on finite element analysis. Experimental results demonstrated that hybrid applications of FRP sheets can improve the shear performance of retrofitted RC beams and increase the ultimate strain of the FRP sheets at failure. The results of the computational model were in reasonable agreement with the corresponding experimental results.
基金This research was funded by the National Natural Science Foundation of China (Grant No.40872191)NSF of Chongqing(Grant No. CSTC2009BB6178)
文摘At present,the thrust of an anti-slide pile can be worked out with some calculation methods. However,the resistance in front of the pile,the distributions of resistance and thrust,and appropriate pile length cannot be easily obtained. In this paper,the authors applied the strength-reduction finite element method (FEM) to several design cases of anti-slide piles. Using this method,it is possible to take the pile-soil interactions into consideration,obtain reasonable resistance in front of pile and the distributions of thrust and resistance,and reasonable lengths of anti-slide piles. In particular,the thrust and resistance imposed on embedded anti-slide piles can be calculated and composite anti-slide pile structures such as anchored piles and braced piles can be optimized. It is proved through the calculation examples that this method is more reliable and economical in the design of anti-slide pile.
基金V. ACKNOWLEDGMENTS This work was supported by the National Natural Science Foundation of China (No.20725312 and No.20533060) and the Ministry of Science and Technology (No.2007CB815201).
文摘New global three dimensional potential energy surfaces for the Cl+H2 reactive system have been constructed using accurate multireference configuration interaction calculations with a large basis set. The three lowest adiabatic potential energy surfaces correlating asymptotically with Cl(^2p)+H2 have been transformed to adiabatic representation, which leads to a fourth coupling potential for non-linear geometries. In addition, the spin-orbit coupling surfaces have also been computed using the Breit-Pauli Hamiltonian. Properties of the new potential are described. Reaction dynamics based on the new potential agrees with the recent experimental results quite well.
文摘In modem missile design, the operation of a missile aerodynamics with angles of attack is required to serve a demand on the maneuverability. The key aero-physics is the development of vortices and its interaction to the control surface such as wing and fins. This paper thus presents the investigation of the missile flow field at 4° and 8° degrees of angles of attack. The Mach numbers for both case were varied from 0.6 to 5.5. Here, the Steady Reynolds-Averaged Navier-Stokes (SRANS) equations with standard κ-ε turbulence model were selected. The numerical results of aerodynamics coefficients (both force and moment) were compared against semi-empirical data computed using Missile DatCOM. The results revealed the development of vortices observed and their interaction with fin at the rear part of the missile.
基金Projects(51104078,51264023)supported by the National Natural Science Foundation of ChinaProject(2010CD022)supported by Yunnan Provincial Applied Basic Research Fund,China+2 种基金Project(IRT1250)supported by the Program for Innovative Research Team in University of Ministry of Education of ChinaProject(U1202271)supported by the National Natural Science Foundation of China-Yunnan United FundProject(KKZ3201252020)supported by Fund of Kunming University of Science and Technology for Talent Training,China
文摘To make sure the intermediate products of the carbothermic reduction of Al2O3 process, such as Al4CO4, Al2CO and Al4C3,and the interaction of AlCl, AlCl2, AlCl3 with Al4CO4, Al2CO and Al4C3, respectively, thermodynamic analyses were used to study the chloride reaction production of them under the vacuum situation. The stable structures and electronic properties of AlCl, AlCl2 and AlCl3 adsorbed on Al4CO4, Al2CO and Al4C3 were calculated by first-principles calculations by the CASTEP module in the Materials Studio program. The results show that the AlCl3 and AlCl2 molecules have decomposed on the plane of Al4C3(001), while there are no obvious decomposition of AlCl3 and AlCl2 on Al4CO4(001) and Al2CO(001) planes. The adsorption of AlCl on the Al4CO4(001)and Al2CO(001) planes is stronger than that on the Al4C3(001) plane. The interaction strength of AlCl3, as well as AlCl2, with Al4CO4,Al2CO and Al4C3 is in the sequence of Al4CO4Al2CO〉Al4C3.
文摘A parallelized resolved method for the simulation of the dynamics of immersed bodies within fluids is presented. The algorithm uses a FDM (fictitious domain method) and combines the Lagrangian DEM (discrete element method) for tracking the bodies with a CFD (computational fluid dynamics) method for calculating the dynamics of the fluid phase. First the CFD-calculation is carried out, disregarding the solid bodies. Afterwards, the velocity information from the bodies is included and the force, the fluid imposes onto the bodies, is computed. The last step consists of a correction-operation which ensures the fulfillment of the conservation equation. Dynamic local mesh refinement is used for minimizing the number of fluid cells. The CFD-DEM coupling is realized within the Open Source framework CFDEMcoupling (www.cfdem.com), where the DEM software LIGGGHTS (www.liggghts.com) is linked against an OpenFOAM^-based CFD solver. While both LIGGGHTS and the CFD solver were already parallelized, only a recent improvement of the algorithm permits the fully parallel computation of resolved problems. This parallelization permits the treatment of large-scale problems. The enclosed validation and application examples show the dynamics of the flow around settling and rotating spheres as well as an investigation of the settling of spheres regarding the Boycott effect.
基金supported by the National Natural Science Foundation of China(Grant Nos.11172040 and 51306020)
文摘The objective of this paper is to improve the understanding of the influence of multiphase flow on the turbulent closure model, the interplay between vorticity fields and cavity dynamics around a pitching hydrofoil. The effects of pitching rate on the sub- cavitating and cavitating response of the pitching hydrofoil are also investigated. In particular, we focus on the interactions between cavity inception, growth, and shedding and the vortex flow structures, and their impacts on the hydrofoil performance. The calculations are 2-D and performed by solving the incompressible, multiphase Unsteady Reynolds Averaged Navier Stokes (URANS) equations via the commercial CFD code CFX. The k-co SST (Shear Stress Transport) turbulence model is used along with the transport equation-based cavitation models. The density correction function is considered to reduce the eddy viscosity according to the computed local fluid mixture density. The calculation results are validated with experiments conducted by Ducoin et al. (see Computational and experimental investigation of flow over a transient pitching hydrofoil, Eur J Mech/B Fluids, 2009, 28:728-743 and An experimental analysis of fluid structure interaction of a flexible hydrofoil in vari- ous flow regimes including cavitating flow, Eur J Mech B/fluids, 2012, 36: 63-74). Results are shown for a NACA66 hydro- foil subject to slow (quasi static, t2=6~/s, &* =0.18) and fast (dynamic, &=63~/s, dr" =1.89) pitching motions from a =0~ to a =15~. Both subcavitaing (or =8.0) and cavitating (cr=3.0) flows are considered. For subcavitating flow (or=8.0), low frequency fluctuations have been observed when the leading edge vortex shedding occurs during stall, and delay of stall is ob- served with increasing pitching velocity. For cavitating flow (tr=3.0), small leading edge cavities are observed with the slow pitching case, which significantly modified the vortex dynamics at high angles of attack, leading to high frequency fluctuations of the hydrodynamic coefficients and different stall behaviors compared to the subcavitating flow at the same pitching rate. On the other hand, for the fast pitching case at or=3.0, large-scale sheet/cloud cavitation is observed, the cavity behavior is un- steady and has a strong impact on the hydrodynamic response, which leads to high amplitude fluctuations of the hydrodynamic coefficients, as well as significant changes in the stall and post-stall behavior. The numerical results also show that the local density modification helps to reduce turbulent eddy viscosity in the cavitating region, which significantly modifies the cavity lengths and shedding frequencies, particularly for the fast pitching case. In general, compared with the experimental visualiza- tions, the numerical results with local density correction have been found to agree well with experimental measurements and observations for both slow and fast transient pitching cases.
文摘Our objective is a better understanding of the role of physical properties of real fluids in the thermodynamics of cavitation in impure water. An extension to the classical homogenous nucleation theory suitable for mixtures is presented in attempt to address the discrepancy between the theoretical predictions and practical observations of cavitation rates in water at normal temperatures. The extension takes into account the non-equilibrium (dissipative) effects involved in nuclei formation through a substance dependent correction coefficient to be determined experimentally. The theory of thermodynamic fluctuations is applied to derive the work of formation of a bubble nucleus. The value of the correction coefficient is estimated using preliminary experimental data from a convergent-divergent nozzle. An application of the results to the numerical prediction of the cavitation zones in a radial-flow water pump is shown.