The transport processes of plasmas in grid systems of krypton（Kr） ion thrusters at different acceleration voltages were simulated with a 3 D-PIC model, and the result was compared with xenon（Xe） ion thrusters. The...The transport processes of plasmas in grid systems of krypton（Kr） ion thrusters at different acceleration voltages were simulated with a 3 D-PIC model, and the result was compared with xenon（Xe） ion thrusters. The variation of the screen grid transparency, the accelerator grid current ratio and the divergence loss were explored. It is found that the screen grid transparency increases with the acceleration voltage and decreases with the beam current, while the accelerator grid current ratio and divergence loss decrease first and then increase with the beam current. This result is the same with Xe ion thrusters. Simulation results also show that Kr ion thrusters have more advantages than Xe ion thrusters, such as higher screen grid transparency, smaller accelerator grid current ratio, larger cut-off current threshold, and better divergence loss characteristic. These advantages mean that Kr ion thrusters have the ability of operating in a wide range of current. Through comprehensive analyses, it can be concluded that using Kr as propellant is very suitable for a multimode ion thruster design.展开更多
A positivity-preserving conservative semi-Lagrangian transport model by multi-moment finite volume method has been developed on the cubed-sphere grid.Two kinds of moments(i.e.,point values(PV moment) at cell interface...A positivity-preserving conservative semi-Lagrangian transport model by multi-moment finite volume method has been developed on the cubed-sphere grid.Two kinds of moments(i.e.,point values(PV moment) at cell interfaces and volume integrated average(VIA moment) value) are defined within a single cell.The PV moment is updated by a conventional semi-Lagrangian method,while the VIA moment is cast by the flux form formulation to assure the exact numerical conservation.Different from the spatial approximation used in the CSL2(conservative semi-Lagrangian scheme with second order polynomial function) scheme,a monotonic rational function which can effectively remove non-physical oscillations is reconstructed within a single cell by the PV moments and VIA moment.To achieve exactly positive-definite preserving,two kinds of corrections are made on the original conservative semi-Lagrangian with rational function(CSLR)scheme.The resulting scheme is inherently conservative,non-negative,and allows a Courant number larger than one.Moreover,the spatial reconstruction can be performed within a single cell,which is very efficient and economical for practical implementation.In addition,a dimension-splitting approach coupled with multi-moment finite volume scheme is adopted on cubed-sphere geometry,which benefitsthe implementation of the 1 D CSLR solver with large Courant number.The proposed model is evaluated by several widely used benchmark tests on cubed-sphere geometry.Numerical results show that the proposed transport model can effectively remove nonphysical oscillations and preserve the numerical nonnegativity,and it has the potential to transport the tracers accurately in a real atmospheric model.展开更多
To meet the requirements of fast and automatic computation of subsonic and transonic aerodynamics in aircraft conceptual design,a novel finite volume solver for full potential flows on adaptive Cartesian grids is deve...To meet the requirements of fast and automatic computation of subsonic and transonic aerodynamics in aircraft conceptual design,a novel finite volume solver for full potential flows on adaptive Cartesian grids is developed in this paper.Cartesian grids with geometric adaptation are firstly generated automatically with boundary cells processed by cell-cutting and cell-merging algorithms.The nonlinear full potential equation is discretized by a finite volume scheme on these Cartesian grids and iteratively solved in an implicit fashion with a generalized minimum residual（GMRES） algorithm.During computation,solution-based mesh adaptation is also applied so as to capture flow features more accurately.An improved ghost-cell method is proposed to implement the non-penetration wall boundary condition where the velocity-potential of a ghost cell is modified by an analytic method instead.According to the characteristics of the Cartesian grids,the Kutta condition is applied by specially computing the gradients on Kutta-faces without directly assigning the potential jump to cells adjacent wake faces,which can significantly improve the solution converging speed.The feasibility and accuracy of the proposed method are validated by several typical cases of sub/transonic flows around an ONERA M6 wing,a DLR-F4 wing-body,and an unconventional figuration of a blended wing body（BWB）.The validation cases demonstrate a fast convergence with fully automatic grid treatment and computation,and the results suggest its capacity in application for aircraft conceptual design.展开更多
As wave propagates into shallow water, the shoaling effect leads to increaseof wave height, and at a certain position, the wave will be breaking. The breaking wave is powerfulagents for generating turbulence, which pl...As wave propagates into shallow water, the shoaling effect leads to increaseof wave height, and at a certain position, the wave will be breaking. The breaking wave is powerfulagents for generating turbulence, which plays an important role in most of the fluid dynamicalprocesses in the surf zone, so a proper numerical model for describing the turbulent effect isneeded urgently. A numerical model is set up to simulate the wave breaking process, which consistsof a free surface model using the surface marker method and the vertical two-dimensional model thatsolves the flow equations. The turbulence is described by Large Eddy Simulation (LES) method wherethe larger turbulent features are simulated by solving the flow equations, and the small-scaleturbulence that is represented by a sub-grid model. A dynamic eddy viscosity sub-grid scale stressmodel has been used for the present simulation. The large eddy simulation model, which we presentedin this paper, can be used to study the propagation of a solitary wave in constant water depth andthe shoaling of a non-breaking solitary wave on a beach. To track free-surface movements, The TUMMACmethod is employed. By applying the model to wave breaking problem in the surf zone, we found thatthese model results compared very well with experimental data. In addition, this model is able toreproduce the complicated flow phenomena, especially the plunging breaker.展开更多
基金co-supported by the National Natural Science Foundation of China (No. 11675040)the Fundamental Research Funds for the Central Universities of China (Nos. 3102014KYJD005 and 1191329723)
文摘The transport processes of plasmas in grid systems of krypton（Kr） ion thrusters at different acceleration voltages were simulated with a 3 D-PIC model, and the result was compared with xenon（Xe） ion thrusters. The variation of the screen grid transparency, the accelerator grid current ratio and the divergence loss were explored. It is found that the screen grid transparency increases with the acceleration voltage and decreases with the beam current, while the accelerator grid current ratio and divergence loss decrease first and then increase with the beam current. This result is the same with Xe ion thrusters. Simulation results also show that Kr ion thrusters have more advantages than Xe ion thrusters, such as higher screen grid transparency, smaller accelerator grid current ratio, larger cut-off current threshold, and better divergence loss characteristic. These advantages mean that Kr ion thrusters have the ability of operating in a wide range of current. Through comprehensive analyses, it can be concluded that using Kr as propellant is very suitable for a multimode ion thruster design.
基金supported by the National Key Research and Development Program of China (Grant Nos.2017YFC1501901 and 2017YFA0603901)the Beijing Natural Science Foundation (Grant No.JQ18001)。
文摘A positivity-preserving conservative semi-Lagrangian transport model by multi-moment finite volume method has been developed on the cubed-sphere grid.Two kinds of moments(i.e.,point values(PV moment) at cell interfaces and volume integrated average(VIA moment) value) are defined within a single cell.The PV moment is updated by a conventional semi-Lagrangian method,while the VIA moment is cast by the flux form formulation to assure the exact numerical conservation.Different from the spatial approximation used in the CSL2(conservative semi-Lagrangian scheme with second order polynomial function) scheme,a monotonic rational function which can effectively remove non-physical oscillations is reconstructed within a single cell by the PV moments and VIA moment.To achieve exactly positive-definite preserving,two kinds of corrections are made on the original conservative semi-Lagrangian with rational function(CSLR)scheme.The resulting scheme is inherently conservative,non-negative,and allows a Courant number larger than one.Moreover,the spatial reconstruction can be performed within a single cell,which is very efficient and economical for practical implementation.In addition,a dimension-splitting approach coupled with multi-moment finite volume scheme is adopted on cubed-sphere geometry,which benefitsthe implementation of the 1 D CSLR solver with large Courant number.The proposed model is evaluated by several widely used benchmark tests on cubed-sphere geometry.Numerical results show that the proposed transport model can effectively remove nonphysical oscillations and preserve the numerical nonnegativity,and it has the potential to transport the tracers accurately in a real atmospheric model.
基金co-supported by the National Natural Science Foundation of China(No.11672133)the Fundamental Research Funds for the Central UniversitiesThe support from the Priority Academic Program Development(PAPD)of Jiangsu Higher Education Institutions
文摘To meet the requirements of fast and automatic computation of subsonic and transonic aerodynamics in aircraft conceptual design,a novel finite volume solver for full potential flows on adaptive Cartesian grids is developed in this paper.Cartesian grids with geometric adaptation are firstly generated automatically with boundary cells processed by cell-cutting and cell-merging algorithms.The nonlinear full potential equation is discretized by a finite volume scheme on these Cartesian grids and iteratively solved in an implicit fashion with a generalized minimum residual（GMRES） algorithm.During computation,solution-based mesh adaptation is also applied so as to capture flow features more accurately.An improved ghost-cell method is proposed to implement the non-penetration wall boundary condition where the velocity-potential of a ghost cell is modified by an analytic method instead.According to the characteristics of the Cartesian grids,the Kutta condition is applied by specially computing the gradients on Kutta-faces without directly assigning the potential jump to cells adjacent wake faces,which can significantly improve the solution converging speed.The feasibility and accuracy of the proposed method are validated by several typical cases of sub/transonic flows around an ONERA M6 wing,a DLR-F4 wing-body,and an unconventional figuration of a blended wing body（BWB）.The validation cases demonstrate a fast convergence with fully automatic grid treatment and computation,and the results suggest its capacity in application for aircraft conceptual design.
文摘As wave propagates into shallow water, the shoaling effect leads to increaseof wave height, and at a certain position, the wave will be breaking. The breaking wave is powerfulagents for generating turbulence, which plays an important role in most of the fluid dynamicalprocesses in the surf zone, so a proper numerical model for describing the turbulent effect isneeded urgently. A numerical model is set up to simulate the wave breaking process, which consistsof a free surface model using the surface marker method and the vertical two-dimensional model thatsolves the flow equations. The turbulence is described by Large Eddy Simulation (LES) method wherethe larger turbulent features are simulated by solving the flow equations, and the small-scaleturbulence that is represented by a sub-grid model. A dynamic eddy viscosity sub-grid scale stressmodel has been used for the present simulation. The large eddy simulation model, which we presentedin this paper, can be used to study the propagation of a solitary wave in constant water depth andthe shoaling of a non-breaking solitary wave on a beach. To track free-surface movements, The TUMMACmethod is employed. By applying the model to wave breaking problem in the surf zone, we found thatthese model results compared very well with experimental data. In addition, this model is able toreproduce the complicated flow phenomena, especially the plunging breaker.
