The circumferentially averaged equation of the inlet flow radial equilibrium in axial compressor was deduced. It indicates that the blade inlet radial pressure gradient is closely related to the radial component of th...The circumferentially averaged equation of the inlet flow radial equilibrium in axial compressor was deduced. It indicates that the blade inlet radial pressure gradient is closely related to the radial component of the circumferential fluctuation（CF） source item. Several simplified cascades with/without aerodynamic loading were numerically studied to investigate the effects of blade bowing on the inlet flow radial equilibrium. A data reduction program was conducted to obtain the CF source from three-dimensional（3D） simulation results. Flow parameters at the passage inlet were focused on and each term in the radial equilibrium equation was discussed quantitatively. Results indicate that the inviscid blade force is the inducement of the inlet CF due to geometrical asymmetry. Blade bowing induces variation of the inlet CF, thus changes the radial pressure gradient and leads to flow migration before leading edge（LE） in the cascades. Positive bowing drives the inlet flow to migrate from end walls to mid-span and negative bowing turns it to the reverse direction to build a new equilibrium. In addition, comparative studies indicate that the inlet Mach number and blade loading can efficiently impact the effectiveness of blade bowing on radial equilibrium in compressor design.展开更多
We derived a theoretical solution of the shock stand-off distance for a non-equilibrium flow over spheres based on Wen and Hornung＇s solution and Olivier＇s solution. Compared with previous approaches, the main advan...We derived a theoretical solution of the shock stand-off distance for a non-equilibrium flow over spheres based on Wen and Hornung＇s solution and Olivier＇s solution. Compared with previous approaches, the main advantage of the present approach is allowing an analytic solution without involving any semi-empirical parameter for the whole non-equilibrium flow regimes. The effects of some important physical quantities therefore can be fully revealed via the analytic solution. By combining the current solution with Ideal Dissociating Gas（IDG） model, we investigate the effects of free stream kinetic energy and free stream dissociation level（which can be very different between different facilities） on the shock stand-off distance.展开更多
Numerical simulations have been carried out for a supersonic three-dimensional rectangular arc nozzle, where a secondary flow toward the center of the curvature occurs due to the shape of the nozzle. It is known that ...Numerical simulations have been carried out for a supersonic three-dimensional rectangular arc nozzle, where a secondary flow toward the center of the curvature occurs due to the shape of the nozzle. It is known that secondary flow causes longitudinal vortices to form near the wall of the nozzle corner, making the nozzle outlet flow unstable and induces loss of transport energy. When the working fluid is a condensable gas with relatively large latent heat such as moist air or steam, rapid accelerated expansion in the nozzle causes non-equilibrium condensation due to supersaturation. After the release of latent heat during phase transition, nozzle flow continues expanding at an equilibrium saturation condition. In the absence of foreign particles, e.g. ions or dust particles, condensation nuclei are formed in the gas itself causing non-equilibrium homogeneous condensation. Supersonic nozzle flow properties vary considerably due to the occurrence of condensation phenomenon. The objective of this study is to investigate the effect of non-equilibrium homogeneous condensation on the longitudinal vortices which form in the range close to the corner of rectangular arc nozzle numerically.展开更多
A code named LARWM with non-ideal magnetohydrodynamic equations in cylindrical model is used to describe the instability in Tokamak plasma surrounded by a conducting wall with finite resistivity. We mainly take three ...A code named LARWM with non-ideal magnetohydrodynamic equations in cylindrical model is used to describe the instability in Tokamak plasma surrounded by a conducting wall with finite resistivity. We mainly take three factors related to the shear equilibrium plasma flow into consideration to study the stabilizing effect of the shear flow on the resistive wall modes （RWMs）. The three factors are the velocity amplitude of flow, the shear rate of flow on plasma surface, and the inertial energy of equilibrium plasma flow. In addition, a local shear plasma flow is also calculated by the LARWM code. Consequently, it is found that the inertial energy of the shear equilibrium plasma flow has an important role in the stabilization of the RWMs.展开更多
In the traffic equilibrium problem, we introduce capacity constraints of arcs, extend Beckmann’s formula to include these constraints, and give an algorithm for traffic equilibrium flows with capacity constraints on ...In the traffic equilibrium problem, we introduce capacity constraints of arcs, extend Beckmann’s formula to include these constraints, and give an algorithm for traffic equilibrium flows with capacity constraints on arcs. Using an example, we illustrate the application of the algorithm and show that Beckmann’s formula is a sufficient condition only, not a necessary condition, for traffic equilibrium with capacity constraints of arcs.展开更多
The effect of the toroidal flow on the equilibrium of tokamak plasmas is a sensitive point for high performance plasma and its precise control. In this paper the effect is studied numerically using the EFIT （Equilibr...The effect of the toroidal flow on the equilibrium of tokamak plasmas is a sensitive point for high performance plasma and its precise control. In this paper the effect is studied numerically using the EFIT （Equilibrium Fitting） code on EAST （Experimental Advanced Superconducting Tokamak）. Firstly, the numerical calculation exhibits a clear outward shift of pressure contour from the magnetic surfaces in the plasma core and the shift grows with the increase of the toroidal velocity. The peak shift of 8% is observed when the ratio between the plasma velocity and the Alfvdn speed equals to 0.15. Secondly, it is shown that the magnetic surfaces shift outwards from those without flow. With a certain plasma current the safety factor on the magnetic axis decreases as the plasma flow velocity increases. The magnetic shear increases about 10% on the plasma boundary compared with the case without flow.展开更多
Problem of the axisymmetric toroidal equilibrium with pure sheared toroidal flow isinvolved. For standard tokamak equilibrium, general approximate solutions are analytically pur-sued for arbitrary current profile and ...Problem of the axisymmetric toroidal equilibrium with pure sheared toroidal flow isinvolved. For standard tokamak equilibrium, general approximate solutions are analytically pur-sued for arbitrary current profile and non-circular cross-section. Equilibrium properties includingthe fiow-induced density asymmetry are analyzed.展开更多
The multipolar velocity field structures are investigated by 2D momentum conservation equation with 3D equilibrium sheared flows in the full toroidal system. Numerical results show that the non-existence of radial vel...The multipolar velocity field structures are investigated by 2D momentum conservation equation with 3D equilibrium sheared flows in the full toroidal system. Numerical results show that the non-existence of radial velocity field in equilibrium surfaces is suitable only for the zero-order term of our 2D simulation. The non-zero-order radial velocity field is still preserved, even when converted to conventional magnetic surface coordinates. The distribution of velocity field vectors of the order of 1, 2, and 3 are presented respectively in 2, 4, and 6 polar fields with the local vortex structure. The excitation mechanisms of these velocity vortices are the coupling effects of the magneto-fluid structure patterns and the toroidal effects. These results can help us understand the complexity of core physics, the transverse transport across magnetic field by the radial plasma flow and the formation of velocity vortices.展开更多
In this paper, a third-order exponential time differencing scheme, named ETDRK3, was investigated for large time stepping in the computation of hypersonic non-equilibrium flow. The second-order Harten-TVD scheme was u...In this paper, a third-order exponential time differencing scheme, named ETDRK3, was investigated for large time stepping in the computation of hypersonic non-equilibrium flow. The second-order Harten-TVD scheme was used for the spatial discretization. The efficient implementation of the scheme with diagonalization of Jacobin matrix was established and carried out for the semi-cylindrical around flow. Current observations showed that the numerical results were in good agreement with those obtained by the classical explicit three-stage Runge-Kutta scheme (RK3) and implicit LU scheme. Efficiency assessments promised the effectiveness of the ETDRK3 scheme. The rationality of the application of this scheme was proved by its preferable accuracy and efficiency.展开更多
The transition between regular reflection （RR） and Mach reflection （MR） of type V shock-shock interaction on a double-wedge geometry with high temperature non-equilibrium effects is investigated by extended shock-...The transition between regular reflection （RR） and Mach reflection （MR） of type V shock-shock interaction on a double-wedge geometry with high temperature non-equilibrium effects is investigated by extended shock-polar method and numerical simulation. First, the critical angles of transition from detachment criterion and yon Neumann criterion are determined by the extended shock-polar method considering the non-equilibrium effects. Then wave patterns and the transition process are numerically obtained. Results of the critical transition angles from shock-polar calculation and numerical simulation show evident disagreement, indicating transition mechanism between RR and MR of type V interaction is changed. By comparing with the frozen counterpart, it is also found that non-equilibrium effects lead to a larger critical wedge angle and a larger hysteresis interval.展开更多
基金supported by the National Natural Science Foundation of China (Nos.51236001,51006005)the National Basic Research Program of China (No. 2012CB720201)Beijing Natural Science Foundation (No. 3151002)
文摘The circumferentially averaged equation of the inlet flow radial equilibrium in axial compressor was deduced. It indicates that the blade inlet radial pressure gradient is closely related to the radial component of the circumferential fluctuation（CF） source item. Several simplified cascades with/without aerodynamic loading were numerically studied to investigate the effects of blade bowing on the inlet flow radial equilibrium. A data reduction program was conducted to obtain the CF source from three-dimensional（3D） simulation results. Flow parameters at the passage inlet were focused on and each term in the radial equilibrium equation was discussed quantitatively. Results indicate that the inviscid blade force is the inducement of the inlet CF due to geometrical asymmetry. Blade bowing induces variation of the inlet CF, thus changes the radial pressure gradient and leads to flow migration before leading edge（LE） in the cascades. Positive bowing drives the inlet flow to migrate from end walls to mid-span and negative bowing turns it to the reverse direction to build a new equilibrium. In addition, comparative studies indicate that the inlet Mach number and blade loading can efficiently impact the effectiveness of blade bowing on radial equilibrium in compressor design.
基金co-supported by the Research Grants Council of Hong Kong,China(No.C5010-14E)the National Natural Science Foundation of China(No.11372265)
文摘We derived a theoretical solution of the shock stand-off distance for a non-equilibrium flow over spheres based on Wen and Hornung＇s solution and Olivier＇s solution. Compared with previous approaches, the main advantage of the present approach is allowing an analytic solution without involving any semi-empirical parameter for the whole non-equilibrium flow regimes. The effects of some important physical quantities therefore can be fully revealed via the analytic solution. By combining the current solution with Ideal Dissociating Gas（IDG） model, we investigate the effects of free stream kinetic energy and free stream dissociation level（which can be very different between different facilities） on the shock stand-off distance.
文摘Numerical simulations have been carried out for a supersonic three-dimensional rectangular arc nozzle, where a secondary flow toward the center of the curvature occurs due to the shape of the nozzle. It is known that secondary flow causes longitudinal vortices to form near the wall of the nozzle corner, making the nozzle outlet flow unstable and induces loss of transport energy. When the working fluid is a condensable gas with relatively large latent heat such as moist air or steam, rapid accelerated expansion in the nozzle causes non-equilibrium condensation due to supersaturation. After the release of latent heat during phase transition, nozzle flow continues expanding at an equilibrium saturation condition. In the absence of foreign particles, e.g. ions or dust particles, condensation nuclei are formed in the gas itself causing non-equilibrium homogeneous condensation. Supersonic nozzle flow properties vary considerably due to the occurrence of condensation phenomenon. The objective of this study is to investigate the effect of non-equilibrium homogeneous condensation on the longitudinal vortices which form in the range close to the corner of rectangular arc nozzle numerically.
基金the National Natural Science Foundation of China(Grant Nos.11105065 and 11275041)the National Basic Research Program of China(Grant Nos.2008CB717801,2008CB787103,2009GB105004,and 2010GB106002)
文摘A code named LARWM with non-ideal magnetohydrodynamic equations in cylindrical model is used to describe the instability in Tokamak plasma surrounded by a conducting wall with finite resistivity. We mainly take three factors related to the shear equilibrium plasma flow into consideration to study the stabilizing effect of the shear flow on the resistive wall modes （RWMs）. The three factors are the velocity amplitude of flow, the shear rate of flow on plasma surface, and the inertial energy of equilibrium plasma flow. In addition, a local shear plasma flow is also calculated by the LARWM code. Consequently, it is found that the inertial energy of the shear equilibrium plasma flow has an important role in the stabilization of the RWMs.
文摘In the traffic equilibrium problem, we introduce capacity constraints of arcs, extend Beckmann’s formula to include these constraints, and give an algorithm for traffic equilibrium flows with capacity constraints on arcs. Using an example, we illustrate the application of the algorithm and show that Beckmann’s formula is a sufficient condition only, not a necessary condition, for traffic equilibrium with capacity constraints of arcs.
基金supported by National Natural Science Foundation of China (No. 10475077)
文摘The effect of the toroidal flow on the equilibrium of tokamak plasmas is a sensitive point for high performance plasma and its precise control. In this paper the effect is studied numerically using the EFIT （Equilibrium Fitting） code on EAST （Experimental Advanced Superconducting Tokamak）. Firstly, the numerical calculation exhibits a clear outward shift of pressure contour from the magnetic surfaces in the plasma core and the shift grows with the increase of the toroidal velocity. The peak shift of 8% is observed when the ratio between the plasma velocity and the Alfvdn speed equals to 0.15. Secondly, it is shown that the magnetic surfaces shift outwards from those without flow. With a certain plasma current the safety factor on the magnetic axis decreases as the plasma flow velocity increases. The magnetic shear increases about 10% on the plasma boundary compared with the case without flow.
基金National Science Foundation (19975015) and the China Nuclear Science Foundation(Y7l00C030l).
文摘Problem of the axisymmetric toroidal equilibrium with pure sheared toroidal flow isinvolved. For standard tokamak equilibrium, general approximate solutions are analytically pur-sued for arbitrary current profile and non-circular cross-section. Equilibrium properties includingthe fiow-induced density asymmetry are analyzed.
基金Supported by National Natural Science Foundation of China under Grant No.11575066Domestic ITER under Grant No.2009GB105003
文摘The multipolar velocity field structures are investigated by 2D momentum conservation equation with 3D equilibrium sheared flows in the full toroidal system. Numerical results show that the non-existence of radial velocity field in equilibrium surfaces is suitable only for the zero-order term of our 2D simulation. The non-zero-order radial velocity field is still preserved, even when converted to conventional magnetic surface coordinates. The distribution of velocity field vectors of the order of 1, 2, and 3 are presented respectively in 2, 4, and 6 polar fields with the local vortex structure. The excitation mechanisms of these velocity vortices are the coupling effects of the magneto-fluid structure patterns and the toroidal effects. These results can help us understand the complexity of core physics, the transverse transport across magnetic field by the radial plasma flow and the formation of velocity vortices.
文摘In this paper, a third-order exponential time differencing scheme, named ETDRK3, was investigated for large time stepping in the computation of hypersonic non-equilibrium flow. The second-order Harten-TVD scheme was used for the spatial discretization. The efficient implementation of the scheme with diagonalization of Jacobin matrix was established and carried out for the semi-cylindrical around flow. Current observations showed that the numerical results were in good agreement with those obtained by the classical explicit three-stage Runge-Kutta scheme (RK3) and implicit LU scheme. Efficiency assessments promised the effectiveness of the ETDRK3 scheme. The rationality of the application of this scheme was proved by its preferable accuracy and efficiency.
文摘The transition between regular reflection （RR） and Mach reflection （MR） of type V shock-shock interaction on a double-wedge geometry with high temperature non-equilibrium effects is investigated by extended shock-polar method and numerical simulation. First, the critical angles of transition from detachment criterion and yon Neumann criterion are determined by the extended shock-polar method considering the non-equilibrium effects. Then wave patterns and the transition process are numerically obtained. Results of the critical transition angles from shock-polar calculation and numerical simulation show evident disagreement, indicating transition mechanism between RR and MR of type V interaction is changed. By comparing with the frozen counterpart, it is also found that non-equilibrium effects lead to a larger critical wedge angle and a larger hysteresis interval.
