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.展开更多
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.展开更多
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 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.展开更多
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.展开更多
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.展开更多
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.展开更多
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 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.展开更多
Decelerating open-channel flow is a type of flow that gradually moves forward with decreasing velocity and increasing water depth.Although all flow parameters change along the streamwise direction,previous studies hav...Decelerating open-channel flow is a type of flow that gradually moves forward with decreasing velocity and increasing water depth.Although all flow parameters change along the streamwise direction,previous studies have revealed that these parameters’vertical distributions at different sections can be universally described with a single profile when being nondimensionalised by appropriate scales.This study focuses on the population trends of spanwise rotational motions at various sections along the main flow direction by particle imaging velocimetry(PIV)measurement.The wall-normal population distributions of density,radius,swirling strength,and convection velocity of the prograde and retrograde motions show similar trends in uniform open-channel flows.The dimensionless representation is invariant along the main flow direction.This study’s results indicate the self-similar characteristic of population trends of spanwise rotational motions prevails in decelerating open-channel flow.展开更多
The dispersion relation of standard geodesic acoustic modes in tokamak plasmas with anisotropic distribution and a radial equilibrium electric field is derived and analyzed. Both frequencies and damping rates increase...The dispersion relation of standard geodesic acoustic modes in tokamak plasmas with anisotropic distribution and a radial equilibrium electric field is derived and analyzed. Both frequencies and damping rates increase with respect to the poloidal Mach number which indicates the strength of the radial electric field. The strength of anisotropy is denoted by the ratio of the parallel temperature（T_‖） to the perpendicular temperature（T_⊥）. It is shown that, when the parallel temperature is lower than the perpendicular temperature, the enhanced anisotropy tends to enlarge the real frequency but reduces the damping rate, and when the parallel temperature is higher than the perpendicular temperature, the effect is opposite. The radial equilibrium electric field has stronger effect on the frequency and damping rate for the case with higher parallel temperature than the case with higher perpendicular temperature.展开更多
This study proposes a quasi-one-dimensional model to predict the chemical nonequilibrium flow along the stagnation streamline of hypersonic flow past a blunt body. The model solves reduced equations along the stagnati...This study proposes a quasi-one-dimensional model to predict the chemical nonequilibrium flow along the stagnation streamline of hypersonic flow past a blunt body. The model solves reduced equations along the stagnation streamline and predicts nearly identical results as the numerical solution of the full-field Navier-Stokes equations. The high efficiency of this model makes it useful to investigate the overall quantitative behavior of related physical-chemical phenomena. In this paper two important properties of hypersonic flow, shock stand-off distance and oxygen dissociation, are studied using the quasi-one-dimensional model with the ideal dissociating gas model. It is found that the shock stand-off distance is affected by both chemical and thermal non-equilibrium.The shock stand-off distance will increase when the flow conditions are changed from equilibrium to non-equilibrium, because the average density of the shock-compressed gas will decrease as a result of the increase in translational energy. For oxygen dissociation, the maximum value of its dissociation degree along the stagnation line varies with the flight altitude. It is increased at first and decreased thereafter with the altitude, which is due to the combination effect of the equilibrium shift and chemical non-equilibrium relaxation. The overall variation of the maximum dissociation is then plotted in the speed and altitude coordinates as a reference for engineering application.展开更多
基金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.
基金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.
文摘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.
基金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.
基金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.
文摘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.
基金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.
基金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.
文摘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.
基金the National Natural Science Foundation of China(Grant No.51679020)the Science and Technology Research Program of Chongqing Municipal Education Commission(Grant No.KJQN202100731).
文摘Decelerating open-channel flow is a type of flow that gradually moves forward with decreasing velocity and increasing water depth.Although all flow parameters change along the streamwise direction,previous studies have revealed that these parameters’vertical distributions at different sections can be universally described with a single profile when being nondimensionalised by appropriate scales.This study focuses on the population trends of spanwise rotational motions at various sections along the main flow direction by particle imaging velocimetry(PIV)measurement.The wall-normal population distributions of density,radius,swirling strength,and convection velocity of the prograde and retrograde motions show similar trends in uniform open-channel flows.The dimensionless representation is invariant along the main flow direction.This study’s results indicate the self-similar characteristic of population trends of spanwise rotational motions prevails in decelerating open-channel flow.
基金supported by National Natural Science Foundation of China(Grant No.11675222)
文摘The dispersion relation of standard geodesic acoustic modes in tokamak plasmas with anisotropic distribution and a radial equilibrium electric field is derived and analyzed. Both frequencies and damping rates increase with respect to the poloidal Mach number which indicates the strength of the radial electric field. The strength of anisotropy is denoted by the ratio of the parallel temperature（T_‖） to the perpendicular temperature（T_⊥）. It is shown that, when the parallel temperature is lower than the perpendicular temperature, the enhanced anisotropy tends to enlarge the real frequency but reduces the damping rate, and when the parallel temperature is higher than the perpendicular temperature, the effect is opposite. The radial equilibrium electric field has stronger effect on the frequency and damping rate for the case with higher parallel temperature than the case with higher perpendicular temperature.
基金supported by the National Natural Science Foundation of China (Nos. 1372325 and 91116013)
文摘This study proposes a quasi-one-dimensional model to predict the chemical nonequilibrium flow along the stagnation streamline of hypersonic flow past a blunt body. The model solves reduced equations along the stagnation streamline and predicts nearly identical results as the numerical solution of the full-field Navier-Stokes equations. The high efficiency of this model makes it useful to investigate the overall quantitative behavior of related physical-chemical phenomena. In this paper two important properties of hypersonic flow, shock stand-off distance and oxygen dissociation, are studied using the quasi-one-dimensional model with the ideal dissociating gas model. It is found that the shock stand-off distance is affected by both chemical and thermal non-equilibrium.The shock stand-off distance will increase when the flow conditions are changed from equilibrium to non-equilibrium, because the average density of the shock-compressed gas will decrease as a result of the increase in translational energy. For oxygen dissociation, the maximum value of its dissociation degree along the stagnation line varies with the flight altitude. It is increased at first and decreased thereafter with the altitude, which is due to the combination effect of the equilibrium shift and chemical non-equilibrium relaxation. The overall variation of the maximum dissociation is then plotted in the speed and altitude coordinates as a reference for engineering application.
