Based on computational fluid dynamics （CFD）/computational eleetromagnetics method （CEM） coupling method and surrogate model optimization techniques, an integration design method about aerodynamic/stealth character...Based on computational fluid dynamics （CFD）/computational eleetromagnetics method （CEM） coupling method and surrogate model optimization techniques, an integration design method about aerodynamic/stealth characteristics of airfoil is established. The O-type body-fitted and orthogonal grid around airfoil is first generated by using the Poisson equations, in which the points per wave and the normal range satisfy the aerodynamic and electromagnetic calculation accuracy requirement. Then the aerodynamic performance of airfoil is calculated by sol- ving the Navier-Stokes （N-S） equations with Baldwin-Lomax （B-L） turbulence model. The stealth characteristics of airfoil are simulated by using finite volume time domain （FVTD） method based on the Maxwell＇s equations, Steger-Warming flux splitting and the third-order MUSCL scheme. In addition, based upon the surrogate model optimization technique with full factorial design （FFD） and radial basis function （RBF）, an integration design about aerodynamic/stealth characteristics of rotor airfoil is conducted by employing the CFD/CEM coupling meth- od. The aerodynamic/stealth characteristics of NACA series airfoils with different maximum thickness and camber combinations are discussed. Finally, by choosing suitable lift-to-drag ratio and radar cross section （RCS） ampli- tudes of rotor airfoil in four important scattering regions as the objective function and constraint, the compromised airfoil with high lift-to-drag ratio and low scattering characteristics is designed via systemic and comprehensive ana- lyses.展开更多
A coupling fluid-structure method with a combination of viscous wake model（VWM）,computational fluid dynamics（CFD） and comprehensive structural dynamics（CSD） modules is developed in this paper for rotor unsteady ...A coupling fluid-structure method with a combination of viscous wake model（VWM）,computational fluid dynamics（CFD） and comprehensive structural dynamics（CSD） modules is developed in this paper for rotor unsteady airload prediction. The hybrid VWM/CFD solver is employed to model the nonlinear aerodynamic phenomena and complicated rotor wake dynamics;the moderate deflection beam theory is implemented to predict the blade structural deformation; the loose coupling strategy based on the ‘delt method＇ is used to couple the fluid and structure solvers.Several cases of Helishape 7A rotor are performed first to investigate the effect of elastic deformation on airloads. Then, two challenging forward flight conditions of UH-60 A helicopter rotor are investigated, and the simulated results of wake geometry, chordwise pressure distribution and sectional normal force show excellent agreement with available test data; a comparison with traditional CFD/CSD method is also presented to illustrate the efficiency of the developed method.展开更多
In this study a mathematical model of a small scale single pellet for the oxidative coupling of methane（OCM）over titanite pervoskite is developed.The method is based on a computational fluid dynamics（CFD）code whic...In this study a mathematical model of a small scale single pellet for the oxidative coupling of methane（OCM）over titanite pervoskite is developed.The method is based on a computational fluid dynamics（CFD）code which known as Fluent may be adopted to model the reactions that take place inside the porous catalyst pellet.The steady state single pellet model is coupled with a kinetic model and the intra-pellet concentration profiles of species are provided.Subsequent to achieving this goal,a nonlinear reaction network consisting of nine catalytic reactions and one gas phase reaction as an external program is successfully implemented to CFD-code as a reaction term in solving the equations.This study is based on the experimental design which is conducted in a differential reactor with a Sn/BaTiO3 catalyst（7-8 mesh） at atmospheric pressure,GHSV of 12000 h-1,ratio of methane to oxygen of 2,and three different temperatures of 1023,1048 and 1073 K.The modeling results such as selectivity and conversion at the pellet exit are in good agreement with the experimental data.Therefore,it is suggested that to achieve high yield in OCM process the modeling of the single pellet should be considered as the heart of catalytic fixed bed reactor.展开更多
In this paper,a high-fidelity computational fluid dynamics(CFD)and rigid body dynamics(RBD)coupled platform for virtual flight simulation is developed to investigate the flight stability of fixed canard dual-spin proj...In this paper,a high-fidelity computational fluid dynamics(CFD)and rigid body dynamics(RBD)coupled platform for virtual flight simulation is developed to investigate the flight stability of fixed canard dual-spin projectile.The platform's reliability is validated by reproducing the characteristic resonance instability of such projectiles.By coupling the solution of the Unsteady Reynolds-Averaged Navier-Stokes equations and the seven-degree-of-freedom RBD equations,the virtual flight simulations of fixed canard dual-spin projectiles at various curvature trajectories are achieved,and the dynamic mechanism of the trajectory following process is analyzed.The instability mechanism of the dynamic instability during trajectory following process of the fixed canard dual-spin projectile is elucidated by simulating the rolling/coning coupled forced motion,and subsequently validated through virtual flight simulations.The findings suggest that an appropriate yaw moment can drive the projectile axis to precession in the tangential direction of the trajectory,thereby enhancing the trajectory following stability.However,the damping of the projectile attains its minimum value when the forward body equilibrium rotational speed(-128 rad/s)is equal to the negative of the fast mode frequency of the projectile.Insufficient damping leads to the fixed canard dual-spin projectile exiting the dynamic stability domain during the trajectory following,resulting in weakly damped instability.Keeping the forward body not rotating or increasing the spin rates to-192 rad/s can enhance the projectile's damping,thereby improving its dynamic stability.展开更多
For accurate aeroelastic analysis,the unsteady rotor flowfield is solved by computational fluid dynamics(CFD)module based on RANS/Euler equations and moving-embedded grid system,while computational structural dynamics...For accurate aeroelastic analysis,the unsteady rotor flowfield is solved by computational fluid dynamics(CFD)module based on RANS/Euler equations and moving-embedded grid system,while computational structural dynamics(CSD)module is introduced to handle blade flexibility.In CFD module,dual time-stepping algorithm is employed in temporal discretization,Jameson two-order central difference(JST)scheme is adopted in spatial discretization and B-L turbulent model is used to illustrate the viscous effect.The CSD module is developed based on Hamilton′s variational principles and moderate deflection beam theory.Grid deformation is implemented using algebraic method through coordinate transformations to achieve deflections with high quality and efficiency.A CFD/CSD loose coupling strategy is developed to transfer information between rotor flowfield and blade structure.The CFD and the CSD modules are verified seperately.Then the CFD/CSD loose coupling is adopted in airloads prediction of UH-60A rotor under high speed forward flight condition.The calculated results agree well with test data.Finally,effects of torsional stiffness properties on airloads of rotors with different tip swept angles(from 10° forward to 30° backward)are investigated.The results are evaluated through pressure distribution and airloads variation,and some meaningful conclusions are drawn the moderated shock wave strength and pressure gradient caused by varied tip swept angle and structural properties.展开更多
The thermal elastic hydro dynamic （TEHD） lubrication analysis for the thrust bearing is usually conducted by combining Reynolds equation with finite element analysis （FEA）. But it is still a problem to conduct the...The thermal elastic hydro dynamic （TEHD） lubrication analysis for the thrust bearing is usually conducted by combining Reynolds equation with finite element analysis （FEA）. But it is still a problem to conduct the computation by combining computational fluid dynamics （CFD） and FEA which can simulate the TEHD more accurately. In this paper, by using both direct and separate coupled solutions together, steady TEHD lubrication considering the viscosity-temperature effect for a bidirectional thrust bearing in a pump-turbine unit is simulated combining a 3D CFD model for the oil film with a 3D FEA model for the pad and mirror plate. Cyclic symmetry condition is used in the oil film flow as more reasonable boundary conditions which avoids the oil temperature assumption at the leading and trailing edge. Deformations of the pad and mirror plate are predicted and discussed as well as the distributions of oil film thickness, pressure, temperature. The predicted temperature shows good agreement with measurements, while the pressure shows a reasonable distribution comparing with previous studies. Further analysis of the three-coupled-field reveals the reason of the high pressure and high temperature generated in the film. Finally, the influence of rotational speed of the mirror plate on the lubrication characteristics is illustrated which shows the thrust load should be balanced against the oil film temperature and pressure in optimized designs. This research proposes a thrust bearing computation method by combining CFD and FEA which can do the TEHD analysis more accurately.展开更多
Flow characteristics, such as flow pattern, gas holdup, and bubble size distribution, in an internal loop reactor with external liquid circulation, are simulated to investigate the influence of reactor internals by us...Flow characteristics, such as flow pattern, gas holdup, and bubble size distribution, in an internal loop reactor with external liquid circulation, are simulated to investigate the influence of reactor internals by using the computational fluid dynamics （CFD）-population balance equations （PBE） coupled model. Numerical results reveal that introducing a downcomer tube and a draft tube can help to improve the mass and heat transfer of the reactor through enhanced liquid circulation, increased gas holdup and reduced bubble diameter. The hydrodynamic behavior in the internal loop reactor with external liquid circulation can be managed effectively by adjusting the diameter and axial position of the draft tube.展开更多
The Lagrangian DDM （discrete droplet model） is state-of-the-art for CFD （computational fluid dynamics） simulations of mixture formation and combustion in industrial engines. A commonly known drawback of the DDM ap...The Lagrangian DDM （discrete droplet model） is state-of-the-art for CFD （computational fluid dynamics） simulations of mixture formation and combustion in industrial engines. A commonly known drawback of the DDM approach is the attenuated validity in the dense spray, where the bulk liquid disintegrates into droplets. There the assumption of single droplets surrounded by a homogenous gas field is not reasonable. In this region, the Eulerian-Eulerian multi-phase approach performs better because instead of parcels the spray is represented by the volume fractions of one bulk liquid and several droplet size class phases. A further drawback of the DDM approach is that increasing the spatial resolution of the computational grid leads to a reduced statistical convergence, since the number of spray parcels per computational cell becomes smaller. It is desirable to combine the benefits of both spray approaches in coupled CFD simulations. Therefore, the dense spray region is simulated separately with the Eulerian spray approach on a highly resolved mesh covering only the region close to the nozzle orifice. The entire engine domain with combustion and emission models is simulated with the Eulerian-Lagrangian spray approach for the dilute spray region. The two simulations are coupled through exchange of boundary conditions and model source terms. An on-line coupling interface manages the data transfer between the two simulation clients, i.e., Eulerian spray and engine client. The aim of this work is to extend the coupled spray approach in terms of exchanging combustion related heat and species sources, and consequently creating the link between Eulerian spray and combustion models. The results show mixture formation and combustion in real-case engine simulations, and demonstrate the feasibility of spray model combination in engineering applications.展开更多
Greenhouse is an important place for crop growth, and it is necessary to control the temperature of growing environment in winter. In addition, the root temperature underground also plays a decisive role for plants gr...Greenhouse is an important place for crop growth, and it is necessary to control the temperature of growing environment in winter. In addition, the root temperature underground also plays a decisive role for plants growth. Adopting underground heating to increase the temperature can effectively improve the yield of crops. The objective of our study was to model the heat transfer of greenhouse underfloor heating which is analyzed and simplified based on the FLUENT software by changing the several important factors that affect the temperature distribution: pipe diameter, pipe spacing, laying depth, supplied water temperature and flow rate, as boundary conditions to simulate the changes of the soil temperature field around the winter night environment. Researching the temperature distribution of the greenhouse, the soil surface and the plant root layer under the different parameters and the basic rules of the heating system are summarized. The results show that the water supply temperature, pipe spacing and diameter of the pipe has a greater impact on the ground and room temperature, and the laying depth has greater impact on the temperature uniformity of the ground, the velocity of water in pipe has little impact on the uniformity of ground temperature.展开更多
对于高超声速飞行器在临近空间形成的连续与稀薄混合流场,DSMC(direct simulation of Monte Carlo)方法需要消耗巨大的计算资源,CFD(computational fluid dynamics)方法,无法对稀薄效应进行准确模拟。在对连续/稀薄耦合数值方法深入研...对于高超声速飞行器在临近空间形成的连续与稀薄混合流场,DSMC(direct simulation of Monte Carlo)方法需要消耗巨大的计算资源,CFD(computational fluid dynamics)方法,无法对稀薄效应进行准确模拟。在对连续/稀薄耦合数值方法深入研究的基础上,提出了一套非结构网格的CFD/DSMC耦合方法。该方法具有通用强、适应性良好的特点,进行耦合计算时对不规则复杂分界面无需进行光滑处理。对超声速圆柱绕流和钝锥体流场进行了模拟,数值结果表明:该方法对不规则分界面和复杂外形具有高度适用性,通过与重叠网格结果进行比较验证了该方法的有效性和计算效率,相较于传统的DSMC方法,计算效率分别提高了2.3倍和3.16倍,具有高效性。展开更多
文摘Based on computational fluid dynamics （CFD）/computational eleetromagnetics method （CEM） coupling method and surrogate model optimization techniques, an integration design method about aerodynamic/stealth characteristics of airfoil is established. The O-type body-fitted and orthogonal grid around airfoil is first generated by using the Poisson equations, in which the points per wave and the normal range satisfy the aerodynamic and electromagnetic calculation accuracy requirement. Then the aerodynamic performance of airfoil is calculated by sol- ving the Navier-Stokes （N-S） equations with Baldwin-Lomax （B-L） turbulence model. The stealth characteristics of airfoil are simulated by using finite volume time domain （FVTD） method based on the Maxwell＇s equations, Steger-Warming flux splitting and the third-order MUSCL scheme. In addition, based upon the surrogate model optimization technique with full factorial design （FFD） and radial basis function （RBF）, an integration design about aerodynamic/stealth characteristics of rotor airfoil is conducted by employing the CFD/CEM coupling meth- od. The aerodynamic/stealth characteristics of NACA series airfoils with different maximum thickness and camber combinations are discussed. Finally, by choosing suitable lift-to-drag ratio and radar cross section （RCS） ampli- tudes of rotor airfoil in four important scattering regions as the objective function and constraint, the compromised airfoil with high lift-to-drag ratio and low scattering characteristics is designed via systemic and comprehensive ana- lyses.
基金supported by the National Natural Science Foundation of China (No. 11302103)
文摘A coupling fluid-structure method with a combination of viscous wake model（VWM）,computational fluid dynamics（CFD） and comprehensive structural dynamics（CSD） modules is developed in this paper for rotor unsteady airload prediction. The hybrid VWM/CFD solver is employed to model the nonlinear aerodynamic phenomena and complicated rotor wake dynamics;the moderate deflection beam theory is implemented to predict the blade structural deformation; the loose coupling strategy based on the ‘delt method＇ is used to couple the fluid and structure solvers.Several cases of Helishape 7A rotor are performed first to investigate the effect of elastic deformation on airloads. Then, two challenging forward flight conditions of UH-60 A helicopter rotor are investigated, and the simulated results of wake geometry, chordwise pressure distribution and sectional normal force show excellent agreement with available test data; a comparison with traditional CFD/CSD method is also presented to illustrate the efficiency of the developed method.
文摘In this study a mathematical model of a small scale single pellet for the oxidative coupling of methane（OCM）over titanite pervoskite is developed.The method is based on a computational fluid dynamics（CFD）code which known as Fluent may be adopted to model the reactions that take place inside the porous catalyst pellet.The steady state single pellet model is coupled with a kinetic model and the intra-pellet concentration profiles of species are provided.Subsequent to achieving this goal,a nonlinear reaction network consisting of nine catalytic reactions and one gas phase reaction as an external program is successfully implemented to CFD-code as a reaction term in solving the equations.This study is based on the experimental design which is conducted in a differential reactor with a Sn/BaTiO3 catalyst（7-8 mesh） at atmospheric pressure,GHSV of 12000 h-1,ratio of methane to oxygen of 2,and three different temperatures of 1023,1048 and 1073 K.The modeling results such as selectivity and conversion at the pellet exit are in good agreement with the experimental data.Therefore,it is suggested that to achieve high yield in OCM process the modeling of the single pellet should be considered as the heart of catalytic fixed bed reactor.
基金supported by the National Natural Science Foundation of China(Grant Nos.U2141254 and U23B6009)。
文摘In this paper,a high-fidelity computational fluid dynamics(CFD)and rigid body dynamics(RBD)coupled platform for virtual flight simulation is developed to investigate the flight stability of fixed canard dual-spin projectile.The platform's reliability is validated by reproducing the characteristic resonance instability of such projectiles.By coupling the solution of the Unsteady Reynolds-Averaged Navier-Stokes equations and the seven-degree-of-freedom RBD equations,the virtual flight simulations of fixed canard dual-spin projectiles at various curvature trajectories are achieved,and the dynamic mechanism of the trajectory following process is analyzed.The instability mechanism of the dynamic instability during trajectory following process of the fixed canard dual-spin projectile is elucidated by simulating the rolling/coning coupled forced motion,and subsequently validated through virtual flight simulations.The findings suggest that an appropriate yaw moment can drive the projectile axis to precession in the tangential direction of the trajectory,thereby enhancing the trajectory following stability.However,the damping of the projectile attains its minimum value when the forward body equilibrium rotational speed(-128 rad/s)is equal to the negative of the fast mode frequency of the projectile.Insufficient damping leads to the fixed canard dual-spin projectile exiting the dynamic stability domain during the trajectory following,resulting in weakly damped instability.Keeping the forward body not rotating or increasing the spin rates to-192 rad/s can enhance the projectile's damping,thereby improving its dynamic stability.
文摘For accurate aeroelastic analysis,the unsteady rotor flowfield is solved by computational fluid dynamics(CFD)module based on RANS/Euler equations and moving-embedded grid system,while computational structural dynamics(CSD)module is introduced to handle blade flexibility.In CFD module,dual time-stepping algorithm is employed in temporal discretization,Jameson two-order central difference(JST)scheme is adopted in spatial discretization and B-L turbulent model is used to illustrate the viscous effect.The CSD module is developed based on Hamilton′s variational principles and moderate deflection beam theory.Grid deformation is implemented using algebraic method through coordinate transformations to achieve deflections with high quality and efficiency.A CFD/CSD loose coupling strategy is developed to transfer information between rotor flowfield and blade structure.The CFD and the CSD modules are verified seperately.Then the CFD/CSD loose coupling is adopted in airloads prediction of UH-60A rotor under high speed forward flight condition.The calculated results agree well with test data.Finally,effects of torsional stiffness properties on airloads of rotors with different tip swept angles(from 10° forward to 30° backward)are investigated.The results are evaluated through pressure distribution and airloads variation,and some meaningful conclusions are drawn the moderated shock wave strength and pressure gradient caused by varied tip swept angle and structural properties.
基金Supported by National Natural Science Foundation of China(Grant No.51439002)Specialized Research Fund for the Doctoral Program of Higher Education of China(Grant Nos.20120002110011,20130002110072)Special Funds for Marine Renewable Energy Projects(Grant no.GHME2012GC02)
文摘The thermal elastic hydro dynamic （TEHD） lubrication analysis for the thrust bearing is usually conducted by combining Reynolds equation with finite element analysis （FEA）. But it is still a problem to conduct the computation by combining computational fluid dynamics （CFD） and FEA which can simulate the TEHD more accurately. In this paper, by using both direct and separate coupled solutions together, steady TEHD lubrication considering the viscosity-temperature effect for a bidirectional thrust bearing in a pump-turbine unit is simulated combining a 3D CFD model for the oil film with a 3D FEA model for the pad and mirror plate. Cyclic symmetry condition is used in the oil film flow as more reasonable boundary conditions which avoids the oil temperature assumption at the leading and trailing edge. Deformations of the pad and mirror plate are predicted and discussed as well as the distributions of oil film thickness, pressure, temperature. The predicted temperature shows good agreement with measurements, while the pressure shows a reasonable distribution comparing with previous studies. Further analysis of the three-coupled-field reveals the reason of the high pressure and high temperature generated in the film. Finally, the influence of rotational speed of the mirror plate on the lubrication characteristics is illustrated which shows the thrust load should be balanced against the oil film temperature and pressure in optimized designs. This research proposes a thrust bearing computation method by combining CFD and FEA which can do the TEHD analysis more accurately.
基金Financial support from the Central Universities (12QN02)National Natural Science Foundation of China (51025624and 51076043)111 Project (B12034)
文摘Flow characteristics, such as flow pattern, gas holdup, and bubble size distribution, in an internal loop reactor with external liquid circulation, are simulated to investigate the influence of reactor internals by using the computational fluid dynamics （CFD）-population balance equations （PBE） coupled model. Numerical results reveal that introducing a downcomer tube and a draft tube can help to improve the mass and heat transfer of the reactor through enhanced liquid circulation, increased gas holdup and reduced bubble diameter. The hydrodynamic behavior in the internal loop reactor with external liquid circulation can be managed effectively by adjusting the diameter and axial position of the draft tube.
文摘The Lagrangian DDM （discrete droplet model） is state-of-the-art for CFD （computational fluid dynamics） simulations of mixture formation and combustion in industrial engines. A commonly known drawback of the DDM approach is the attenuated validity in the dense spray, where the bulk liquid disintegrates into droplets. There the assumption of single droplets surrounded by a homogenous gas field is not reasonable. In this region, the Eulerian-Eulerian multi-phase approach performs better because instead of parcels the spray is represented by the volume fractions of one bulk liquid and several droplet size class phases. A further drawback of the DDM approach is that increasing the spatial resolution of the computational grid leads to a reduced statistical convergence, since the number of spray parcels per computational cell becomes smaller. It is desirable to combine the benefits of both spray approaches in coupled CFD simulations. Therefore, the dense spray region is simulated separately with the Eulerian spray approach on a highly resolved mesh covering only the region close to the nozzle orifice. The entire engine domain with combustion and emission models is simulated with the Eulerian-Lagrangian spray approach for the dilute spray region. The two simulations are coupled through exchange of boundary conditions and model source terms. An on-line coupling interface manages the data transfer between the two simulation clients, i.e., Eulerian spray and engine client. The aim of this work is to extend the coupled spray approach in terms of exchanging combustion related heat and species sources, and consequently creating the link between Eulerian spray and combustion models. The results show mixture formation and combustion in real-case engine simulations, and demonstrate the feasibility of spray model combination in engineering applications.
文摘Greenhouse is an important place for crop growth, and it is necessary to control the temperature of growing environment in winter. In addition, the root temperature underground also plays a decisive role for plants growth. Adopting underground heating to increase the temperature can effectively improve the yield of crops. The objective of our study was to model the heat transfer of greenhouse underfloor heating which is analyzed and simplified based on the FLUENT software by changing the several important factors that affect the temperature distribution: pipe diameter, pipe spacing, laying depth, supplied water temperature and flow rate, as boundary conditions to simulate the changes of the soil temperature field around the winter night environment. Researching the temperature distribution of the greenhouse, the soil surface and the plant root layer under the different parameters and the basic rules of the heating system are summarized. The results show that the water supply temperature, pipe spacing and diameter of the pipe has a greater impact on the ground and room temperature, and the laying depth has greater impact on the temperature uniformity of the ground, the velocity of water in pipe has little impact on the uniformity of ground temperature.
文摘对于高超声速飞行器在临近空间形成的连续与稀薄混合流场,DSMC(direct simulation of Monte Carlo)方法需要消耗巨大的计算资源,CFD(computational fluid dynamics)方法,无法对稀薄效应进行准确模拟。在对连续/稀薄耦合数值方法深入研究的基础上,提出了一套非结构网格的CFD/DSMC耦合方法。该方法具有通用强、适应性良好的特点,进行耦合计算时对不规则复杂分界面无需进行光滑处理。对超声速圆柱绕流和钝锥体流场进行了模拟,数值结果表明:该方法对不规则分界面和复杂外形具有高度适用性,通过与重叠网格结果进行比较验证了该方法的有效性和计算效率,相较于传统的DSMC方法,计算效率分别提高了2.3倍和3.16倍,具有高效性。
