Gas flows through micro shock tubes are widely used in many engineering applications such as micro engines, particle delivery devices etc. Recently, few studies have been carried out to explore the shock wave excursio...Gas flows through micro shock tubes are widely used in many engineering applications such as micro engines, particle delivery devices etc. Recently, few studies have been carried out to explore the shock wave excursions through micro shock tubes at very low Reynolds number and at rarefied gas condition. But these studies assumed centered shock and expansion waves, which are generally produced by instantaneous diaphragm rupture process. But in real scenario, the diaphragm ruptures with a finite rupture time and this phenomenon will significantly alter the shock wave propagation characteristics. In the present research, numerical simulations have been carried out on a two dimensional micro shock tube model to simulate the effect of finite diaphragm rupture process on the wave characteristics. The rarefaction effect was simulated using Maxwell’s slip wall equations. The results show that shock wave attenuates rapidly in micro shock tubes compared to conventional macro shock tubes. Finite diaphragm rupture causes the formation of non-centered shock wave at some distance ahead of the diaphragm. The shock propagation distance is also drastically reduced by the rupture effects.展开更多
In the present paper, random-choice method (RCM) and second-order GRP difference method, which are high resolution methods used for pure gas flows with shocks, are extended and employed to study the problem of one-dim...In the present paper, random-choice method (RCM) and second-order GRP difference method, which are high resolution methods used for pure gas flows with shocks, are extended and employed to study the problem of one-dimensional unsteady two-phase flows. The two-phase shock wave and the flow field behind it in a dusty gas shock tube are calculated and the time-dependent change of the flow parameters for the gas and particle phase are obtained. The numerical results indicate that both the two methods can give the relaxation structure of the two-phase shocks with a sharp discontinuous front and that the GRP method has the advantages of less time-consuming and higher accuracy over the RCM method.展开更多
Several kinds of explicit and implicit finite-difference schemes directly solving the discretized velocity distribution functions are designed with precision of different orders by analyzing the inner characteristics ...Several kinds of explicit and implicit finite-difference schemes directly solving the discretized velocity distribution functions are designed with precision of different orders by analyzing the inner characteristics of the gas-kinetic numerical algorithm for Boltzmann model equation. The peculiar flow phenomena and mechanism from various flow regimes are revealed in the numerical simulations of the unsteady Sod shock-tube problems and the two-dimensional channel flows with different Knudsen numbers. The numerical remainder-effects of the difference schemes are investigated aad analyzed based on the computed results. The ways of improving the computational efficiency of the gaskinetic numerical method and the computing principles of difference discretization are discussed.展开更多
Recently micro shock tubes have been widely used in many engineering and industrial fields, but the characteristics of unsteady flow are not well known to date in micro shock tubes. Compared to conventional shock tube...Recently micro shock tubes have been widely used in many engineering and industrial fields, but the characteristics of unsteady flow are not well known to date in micro shock tubes. Compared to conventional shock tubes with macro scales, flows related to shock waves in micro shock tubes are highly complicated. Stronger viscous and dissipative interactions make shock wave dynamic behaviors significantly different from theoretical predictions. In the present study, a CFD work was applied to the unsteady compressible Navier-Stokes equations which were solved using a fully implicit finite volume scheme. The diaphragm pressure ratio and shock tube diameter were varied to investigate their effects on micro shock tube flows. Different wall boundary conditions were also performed to observe shock wave and contact surface propagation with no slip and slip walls. Detailed flow characteristics at the foot of shock wave and contact surface propagation were known from the present numerical simulations.展开更多
Fine bubbles have been widely applied in many fields such as industry,medical engineering and agricultures.Therefore,many attentions have been paid to the study of fine bubble generations in order to increase the yiel...Fine bubbles have been widely applied in many fields such as industry,medical engineering and agricultures.Therefore,many attentions have been paid to the study of fine bubble generations in order to increase the yield while decrease the cost.However,the generation process of fine bubbles is a quite complicated process in which multiple hydrodynamic forces are interacted in the gas-liquid two-phase flow.Many studies focus on the techniques of the converging-diverging nozzle(venturi tube)generator,which is famous for its simple and cheap features,and generates fine bubbles by using the miniaturization phenomenon of bubbles occurring in the venturi tube.However,the impact conditions on the amount and size of bubbles such as nozzle geometry and bleed air haven’t been investigated clearly.In this work,we implement many experiments on the venturi tube fine bubble generators with different geometries and generating conditions,and evaluate different factors impacting the production components such as the volume and the bubble size.The experimental results show that the supersonic flow filed in the venturi tube promotes the miniaturization of the bubbles,and the convergent angle of the nozzle and air bleed have a great impact on the size and volume of bubbles.展开更多
In the present study, the flow visualizations were performed around the NACA 0012 models which differ in aspect ratios. We discussed the effects of the aspect ratio in the test models. Additionally the unsteady, two-d...In the present study, the flow visualizations were performed around the NACA 0012 models which differ in aspect ratios. We discussed the effects of the aspect ratio in the test models. Additionally the unsteady, two-dimensional, compressible Euler equations were solved for the NACA 0012 airfoil. Experiments were performed utilizing the conventional gas driven shock tube as the intermittent transonic wind tunnel. The aspect ratios of the models are about 0.86 and 1.5, respectively. The Mach numbers M 2 are about 0.84. The Reynolds numbers of the present experimental conditions were constant that Re based on chord length is about 4.0×10 5 . The results are as follows: in different aspect ratios, the difference of the shock wave location is confirmed though the Mach number and Reynolds number are same. It indicates the different correction Mach number by the effects of the side wall boundary layer though the nominal Mach number measured the same value. Also, on the difference of shock wave location for the effects of the aspect ratio, the tend of CFD shows the qualitative agreement with the result of an experiment.展开更多
The comparative study on seven equation models with two different six equations modelfor compressible two-phase flow analysis is proposed. The seven equations model isderived for compressible two-phase flow that is in...The comparative study on seven equation models with two different six equations modelfor compressible two-phase flow analysis is proposed. The seven equations model isderived for compressible two-phase flow that is in the nonconservation form. In thepresent work, two different six equations model are derived for two pressures, two velocities and single temperature with the derivation of the equation of state. The closingequation for one of the six equations model is energy conservation equation while anotherone is closed by entropy balance equation. The partial differential form of governingequations is hyperbolic and written in the conservative form. At this point, the set ofgoverning equations are derived based on the principle of extended thermodynamics.The method of solving single temperature from both six equation models are simple anddirect solution can be obtained. Numerical simulation has been tried using one of the sixequation models for air–water shock tube problems. Explicit fourth order Runge–Kuttascheme is used with Finite Volume Shock Capturing method for solving the governingequations numerically. The pressure, velocity and volume fraction variations are captured along the shock tube length through flow solver. Experimental work is carried outto magnify the initial stage of liquid injection into a gas. The outcome of six equationsmodel for compressible two-phase flow has revealed the multi-phase flow characteristicsthat are similar to the actual conditions.展开更多
The time-dependent behavior of non-equilibrium condensation of moist air through a Ludwieg tube with a diaphragm downstream is investigated by using a computational fluid dynamics work. The two-dimensional, compressib...The time-dependent behavior of non-equilibrium condensation of moist air through a Ludwieg tube with a diaphragm downstream is investigated by using a computational fluid dynamics work. The two-dimensional, compressible, Navier-Stokes equations, fully coupled with the condensate droplet growth equations, are numerically solved by a third-order MUSCL type TVD finite-difference scheme with a second-order fractional time step. Baldwin-Lomax turbulence model is employed to close the governing equations. The present computations represent the experimental flows well. The results obtained show that for an initial relative humidity over 40 %, the periodic excursions of the condensation shock occurs in the Ludwieg tube, and the frequency increases with the initial relative humidity. It is also found that total pressure loss due to non-equilibrium condensation in the Ludwieg tube should not be ignored even for a very low initial relative humidity. Furthermore, the variations of condensation properties are also展开更多
文摘Gas flows through micro shock tubes are widely used in many engineering applications such as micro engines, particle delivery devices etc. Recently, few studies have been carried out to explore the shock wave excursions through micro shock tubes at very low Reynolds number and at rarefied gas condition. But these studies assumed centered shock and expansion waves, which are generally produced by instantaneous diaphragm rupture process. But in real scenario, the diaphragm ruptures with a finite rupture time and this phenomenon will significantly alter the shock wave propagation characteristics. In the present research, numerical simulations have been carried out on a two dimensional micro shock tube model to simulate the effect of finite diaphragm rupture process on the wave characteristics. The rarefaction effect was simulated using Maxwell’s slip wall equations. The results show that shock wave attenuates rapidly in micro shock tubes compared to conventional macro shock tubes. Finite diaphragm rupture causes the formation of non-centered shock wave at some distance ahead of the diaphragm. The shock propagation distance is also drastically reduced by the rupture effects.
文摘In the present paper, random-choice method (RCM) and second-order GRP difference method, which are high resolution methods used for pure gas flows with shocks, are extended and employed to study the problem of one-dimensional unsteady two-phase flows. The two-phase shock wave and the flow field behind it in a dusty gas shock tube are calculated and the time-dependent change of the flow parameters for the gas and particle phase are obtained. The numerical results indicate that both the two methods can give the relaxation structure of the two-phase shocks with a sharp discontinuous front and that the GRP method has the advantages of less time-consuming and higher accuracy over the RCM method.
基金supported by the National Natural Science Foundation of China (No.10621062)the Research Fund for Next Generation of General Armament Department (No.9140A13050207KG29)
文摘Several kinds of explicit and implicit finite-difference schemes directly solving the discretized velocity distribution functions are designed with precision of different orders by analyzing the inner characteristics of the gas-kinetic numerical algorithm for Boltzmann model equation. The peculiar flow phenomena and mechanism from various flow regimes are revealed in the numerical simulations of the unsteady Sod shock-tube problems and the two-dimensional channel flows with different Knudsen numbers. The numerical remainder-effects of the difference schemes are investigated aad analyzed based on the computed results. The ways of improving the computational efficiency of the gaskinetic numerical method and the computing principles of difference discretization are discussed.
基金supported by the National Research Foundation of Korea(NRF)grant funded by the Korea government(MEST)(2011-0017506)
文摘Recently micro shock tubes have been widely used in many engineering and industrial fields, but the characteristics of unsteady flow are not well known to date in micro shock tubes. Compared to conventional shock tubes with macro scales, flows related to shock waves in micro shock tubes are highly complicated. Stronger viscous and dissipative interactions make shock wave dynamic behaviors significantly different from theoretical predictions. In the present study, a CFD work was applied to the unsteady compressible Navier-Stokes equations which were solved using a fully implicit finite volume scheme. The diaphragm pressure ratio and shock tube diameter were varied to investigate their effects on micro shock tube flows. Different wall boundary conditions were also performed to observe shock wave and contact surface propagation with no slip and slip walls. Detailed flow characteristics at the foot of shock wave and contact surface propagation were known from the present numerical simulations.
基金supported by the Mishima Laboratory in Fukuoka University。
文摘Fine bubbles have been widely applied in many fields such as industry,medical engineering and agricultures.Therefore,many attentions have been paid to the study of fine bubble generations in order to increase the yield while decrease the cost.However,the generation process of fine bubbles is a quite complicated process in which multiple hydrodynamic forces are interacted in the gas-liquid two-phase flow.Many studies focus on the techniques of the converging-diverging nozzle(venturi tube)generator,which is famous for its simple and cheap features,and generates fine bubbles by using the miniaturization phenomenon of bubbles occurring in the venturi tube.However,the impact conditions on the amount and size of bubbles such as nozzle geometry and bleed air haven’t been investigated clearly.In this work,we implement many experiments on the venturi tube fine bubble generators with different geometries and generating conditions,and evaluate different factors impacting the production components such as the volume and the bubble size.The experimental results show that the supersonic flow filed in the venturi tube promotes the miniaturization of the bubbles,and the convergent angle of the nozzle and air bleed have a great impact on the size and volume of bubbles.
文摘In the present study, the flow visualizations were performed around the NACA 0012 models which differ in aspect ratios. We discussed the effects of the aspect ratio in the test models. Additionally the unsteady, two-dimensional, compressible Euler equations were solved for the NACA 0012 airfoil. Experiments were performed utilizing the conventional gas driven shock tube as the intermittent transonic wind tunnel. The aspect ratios of the models are about 0.86 and 1.5, respectively. The Mach numbers M 2 are about 0.84. The Reynolds numbers of the present experimental conditions were constant that Re based on chord length is about 4.0×10 5 . The results are as follows: in different aspect ratios, the difference of the shock wave location is confirmed though the Mach number and Reynolds number are same. It indicates the different correction Mach number by the effects of the side wall boundary layer though the nominal Mach number measured the same value. Also, on the difference of shock wave location for the effects of the aspect ratio, the tend of CFD shows the qualitative agreement with the result of an experiment.
文摘The comparative study on seven equation models with two different six equations modelfor compressible two-phase flow analysis is proposed. The seven equations model isderived for compressible two-phase flow that is in the nonconservation form. In thepresent work, two different six equations model are derived for two pressures, two velocities and single temperature with the derivation of the equation of state. The closingequation for one of the six equations model is energy conservation equation while anotherone is closed by entropy balance equation. The partial differential form of governingequations is hyperbolic and written in the conservative form. At this point, the set ofgoverning equations are derived based on the principle of extended thermodynamics.The method of solving single temperature from both six equation models are simple anddirect solution can be obtained. Numerical simulation has been tried using one of the sixequation models for air–water shock tube problems. Explicit fourth order Runge–Kuttascheme is used with Finite Volume Shock Capturing method for solving the governingequations numerically. The pressure, velocity and volume fraction variations are captured along the shock tube length through flow solver. Experimental work is carried outto magnify the initial stage of liquid injection into a gas. The outcome of six equationsmodel for compressible two-phase flow has revealed the multi-phase flow characteristicsthat are similar to the actual conditions.
文摘The time-dependent behavior of non-equilibrium condensation of moist air through a Ludwieg tube with a diaphragm downstream is investigated by using a computational fluid dynamics work. The two-dimensional, compressible, Navier-Stokes equations, fully coupled with the condensate droplet growth equations, are numerically solved by a third-order MUSCL type TVD finite-difference scheme with a second-order fractional time step. Baldwin-Lomax turbulence model is employed to close the governing equations. The present computations represent the experimental flows well. The results obtained show that for an initial relative humidity over 40 %, the periodic excursions of the condensation shock occurs in the Ludwieg tube, and the frequency increases with the initial relative humidity. It is also found that total pressure loss due to non-equilibrium condensation in the Ludwieg tube should not be ignored even for a very low initial relative humidity. Furthermore, the variations of condensation properties are also
