Unsteady-state operation has been widely applied in chemical engineering, such as optimizing a process, increasing yield and saving energy, etc. But the knowledge of the flow characteristics in bubble column reactors(...Unsteady-state operation has been widely applied in chemical engineering, such as optimizing a process, increasing yield and saving energy, etc. But the knowledge of the flow characteristics in bubble column reactors(BCRs) under unsteady state control is far from enough. In order to study the flow structures in this operation, the volume of fluid (VOF) model and the standard k-ε model to simulate the evolution of gas-liquid flow in BCRs under the start-up state are combined. For both the symmetry and asymmetry flow, the layout of the gas-inlets, the gas-in velocity, the liquid viscosity and the aspect ratio of the BCR all have effects on the liquid velocity distribution. The simulation results could provide some information for the design and scale-up of the BCRs.展开更多
In a magnetohydrodynamic(MHD)driven fluid cell,a plane non-parallel flow in a square domain satisfying a free-slip boundary condition is examined.The energy dissipation of the flow is controlled by the viscosity and l...In a magnetohydrodynamic(MHD)driven fluid cell,a plane non-parallel flow in a square domain satisfying a free-slip boundary condition is examined.The energy dissipation of the flow is controlled by the viscosity and linear friction.The latter arises from the influence of the Hartmann bottom boundary layer in a three-dimensional(3D)MHD experiment in a square bottomed cell.The basic flow in this fluid system is a square eddy flow exhibiting a network of N~2 vortices rotating alternately in clockwise and anticlockwise directions.When N is odd,the instability of the flow gives rise to secondary steady-state flows and secondary time-periodic flows,exhibiting similar characteristics to those observed when N=3.For this reason,this study focuses on the instability of the square eddy flow of nine vortices.It is shown that there exist eight bi-critical values corresponding to the existence of eight neutral eigenfunction spaces.Especially,there exist non-real neutral eigenfunctions,which produce secondary time-periodic flows exhibiting vortices merging in an oscillatory manner.This Hopf bifurcation phenomenon has not been observed in earlier investigations.展开更多
In this paper, we design high-order Runge-Kutta discontinuous Galerkin (RKDG) methods with multi-resolution weighted essentially non-oscillatory (multi-resolution WENO) limiters to compute compressible steady-state pr...In this paper, we design high-order Runge-Kutta discontinuous Galerkin (RKDG) methods with multi-resolution weighted essentially non-oscillatory (multi-resolution WENO) limiters to compute compressible steady-state problems on triangular meshes. A troubled cell indicator extended from structured meshes to unstructured meshes is constructed to identify triangular cells in which the application of the limiting procedures is required. In such troubled cells, the multi-resolution WENO limiting methods are used to the hierarchical L^(2) projection polynomial sequence of the DG solution. Through using the RKDG methods with multi-resolution WENO limiters, the optimal high-order accuracy can be gradually reduced to first-order in the triangular troubled cells, so that the shock wave oscillations can be well suppressed. In steady-state simulations on triangular meshes, the numerical residual converges to near machine zero. The proposed spatial reconstruction methods enhance the robustness of classical DG methods on triangular meshes. The good results of these RKDG methods with multi-resolution WENO limiters are verified by a series of two-dimensional steady-state problems.展开更多
This study aims to develop a mathematical analysis for one-dimensional modeling of a radial flow through a production well drilled in a confined aquifer, in the case of steady-state flow conditions. An analytical solu...This study aims to develop a mathematical analysis for one-dimensional modeling of a radial flow through a production well drilled in a confined aquifer, in the case of steady-state flow conditions. An analytical solution has derived from that expression for estimation of drawdowns according to different flowrates. Through that process, the evaluation of static pressure, the calculation of hydraulic charge due to the waterflow through the well is evaluated, the drawdowns curves are drawn and at last, the obtained curves are analyzed. The curves obtained for the different flow rates have an asymptotic direction, the axis of the hydraulic charges. The variation of the hydraulic charge depends on the radial distance for different flow rates. The P point, is a common point of all curves obtained for different production flowrates in the well. This point is where the well production flowrate is optimum for the optimal hydraulic charge.展开更多
The current research about the flow ripple of axial piston pump mainly focuses on the effect of the structure of parts on the flow ripple. Therein, the structure of parts are usually designed and optimized at rated wo...The current research about the flow ripple of axial piston pump mainly focuses on the effect of the structure of parts on the flow ripple. Therein, the structure of parts are usually designed and optimized at rated working conditions. However, the pump usually has to work in large-scale and time-variant working conditions. Therefore, the flow ripple characteristics of pump and analysis for its test accuracy with respect to variant steady-state conditions and transient conditions in a wide range of operating parameters are focused in this paper. First, a simulation model has been constructed, which takes the kinematics of oil film within friction pairs into account for higher accuracy. Afterwards, a test bed which adopts Secondary Source Method is built to verify the model. The simulation and tests results show that the angular position of the piston, corresponding to the position where the peak flow ripple is produced, varies with the different pressure. The pulsating amplitude and pulsation rate of flow ripple increase with the rise of pressure and the variation rate of pressure. For the pump working at a constant speed, the flow pulsation rate decreases dramatically with the increasing speed when the speed is less than 27.78% of the maximum speed, subsequently presents a small decrease tendency with the speed further increasing. With the rise of the variation rate of speed, the pulsating amplitude and pulsation rate of flow ripple increase. As the swash plate angle augments, the pulsating amplitude of flow ripple increases, nevertheless the flow pulsation rate decreases. In contrast with the effect of the variation of pressure, the test accuracy of flow ripple is more sensitive to the variation of speed. It makes the test accuracy above 96.20% available for the pulsating amplitude of pressure deviating within a range of ~6% from the mean pressure. However, with a variation of speed deviating within a range of ±2% from the mean speed, the attainable test accuracy of flow ripple is above 93.07%. The model constructed in this research proposes a method to determine the flow ripple characteristics of pump and its attainable test accuracy under the large-scale and time-variant working conditions. Meanwhile, a discussion about the variation of flow ripple and its obtainable test accuracy with the conditions of the pump working in wide operating ranges is given as well.展开更多
Liquid hydrogen storage and transportation is an effective method for large-scale transportation and utilization of hydrogen energy. Revealing the flow mechanism of cryogenic working fluid is the key to optimize heat ...Liquid hydrogen storage and transportation is an effective method for large-scale transportation and utilization of hydrogen energy. Revealing the flow mechanism of cryogenic working fluid is the key to optimize heat exchanger structure and hydrogen liquefaction process(LH2). The methods of cryogenic visualization experiment, theoretical analysis and numerical simulation are conducted to study the falling film flow characteristics with the effect of co-current gas flow in LH2spiral wound heat exchanger.The results show that the flow rate of mixed refrigerant has a great influence on liquid film spreading process, falling film flow pattern and heat transfer performance. The liquid film of LH2mixed refrigerant with column flow pattern can not uniformly and completely cover the tube wall surface. As liquid flow rate increases, the falling film flow pattern evolves into sheet-column flow and sheet flow, and liquid film completely covers the surface of tube wall. With the increase of shear effect of gas-phase mixed refrigerant in the same direction, the liquid film gradually becomes unstable, and the flow pattern eventually evolves into a mist flow.展开更多
Gassolid hydrodynamic steadystate operation is the operating basis in a chemical looping dualreactor system.This study reported the experimental results on the steadystate operation characteristics of gassolid flow in...Gassolid hydrodynamic steadystate operation is the operating basis in a chemical looping dualreactor system.This study reported the experimental results on the steadystate operation characteristics of gassolid flow in a 15.5 m high dual circulating fluidized bed(CFB)cold test system.The effects of superficial gas velocity,static bed material height and solid returning modes on the steadystate operation characteristics between the two CFBs were investigated.Results suggest that the solid distributions in the dual CFB test system was mainly determined by the superficial gas velocity and larger solid inventory may help to improve the solid distributions.Besides,crossreturning mode coupled with selfreturning is good for steadystate running in the dualreactor test system.展开更多
The majority of nonlinear stochastic systems can be expressed as the quasi-Hamiltonian systems in science and engineering. Moreover, the corresponding Hamiltonian system offers two concepts of integrability and resona...The majority of nonlinear stochastic systems can be expressed as the quasi-Hamiltonian systems in science and engineering. Moreover, the corresponding Hamiltonian system offers two concepts of integrability and resonance that can fully describe the global relationship among the degrees-of-freedom(DOFs) of the system. In this work, an effective and promising approximate semi-analytical method is proposed for the steady-state response of multi-dimensional quasi-Hamiltonian systems. To be specific, the trial solution of the reduced Fokker–Plank–Kolmogorov(FPK) equation is obtained by using radial basis function(RBF) neural networks. Then, the residual generated by substituting the trial solution into the reduced FPK equation is considered, and a loss function is constructed by combining random sampling technique. The unknown weight coefficients are optimized by minimizing the loss function through the Lagrange multiplier method. Moreover, an efficient sampling strategy is employed to promote the implementation of algorithms. Finally, two numerical examples are studied in detail, and all the semi-analytical solutions are compared with Monte Carlo simulations(MCS) results. The results indicate that the complex nonlinear dynamic features of the system response can be captured through the proposed scheme accurately.展开更多
We propose an adaptive stencil construction for high-order accurate finite volume schemes a posteriori stabilized devoted to solve one-dimensional steady-state hyperbolic equations.High accuracy(up to the sixth-order ...We propose an adaptive stencil construction for high-order accurate finite volume schemes a posteriori stabilized devoted to solve one-dimensional steady-state hyperbolic equations.High accuracy(up to the sixth-order presently)is achieved,thanks to polynomial recon-structions while stability is provided with an a posteriori MOOD method which controls the cell polynomial degree for eliminating non-physical oscillations in the vicinity of dis-continuities.We supplemented this scheme with a stencil construction allowing to reduce even further the numerical dissipation.The stencil is shifted away from troubles(shocks,discontinuities,etc.)leading to less oscillating polynomial reconstructions.Experimented on linear,Burgers',and Euler equations,we demonstrate that the adaptive stencil technique manages to retrieve smooth solutions with optimal order of accuracy but also irregular ones without spurious oscillations.Moreover,we numerically show that the approach allows to reduce the dissipation still maintaining the essentially non-oscillatory behavior.展开更多
A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes...A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes due to stress sensitivity, which plays a crucial role in controlling pressure propagation and oil flow. This paper proposes a multi-scale coupled flow mathematical model of matrix nanopores, induced fractures, and hydraulic fractures. In this model, the micro-scale effects of shale oil flow in fractal nanopores, fractal induced fracture network, and stress sensitivity of multi-scale media are considered. We solved the model iteratively using Pedrosa transform, semi-analytic Segmented Bessel function, Laplace transform. The results of this model exhibit good agreement with the numerical solution and field production data, confirming the high accuracy of the model. As well, the influence of stress sensitivity on permeability, pressure and production is analyzed. It is shown that the permeability and production decrease significantly when induced fractures are weakly supported. Closed induced fractures can inhibit interporosity flow in the stimulated reservoir volume (SRV). It has been shown in sensitivity analysis that hydraulic fractures are beneficial to early production, and induced fractures in SRV are beneficial to middle production. The model can characterize multi-scale flow characteristics of shale oil, providing theoretical guidance for rapid productivity evaluation.展开更多
文摘Unsteady-state operation has been widely applied in chemical engineering, such as optimizing a process, increasing yield and saving energy, etc. But the knowledge of the flow characteristics in bubble column reactors(BCRs) under unsteady state control is far from enough. In order to study the flow structures in this operation, the volume of fluid (VOF) model and the standard k-ε model to simulate the evolution of gas-liquid flow in BCRs under the start-up state are combined. For both the symmetry and asymmetry flow, the layout of the gas-inlets, the gas-in velocity, the liquid viscosity and the aspect ratio of the BCR all have effects on the liquid velocity distribution. The simulation results could provide some information for the design and scale-up of the BCRs.
基金Project supported by the National Natural Science Foundation of China(No.11571240)the Shenzhen Natural Science Fund of China(the Stable Support Plan Program No.20220805175116001)。
文摘In a magnetohydrodynamic(MHD)driven fluid cell,a plane non-parallel flow in a square domain satisfying a free-slip boundary condition is examined.The energy dissipation of the flow is controlled by the viscosity and linear friction.The latter arises from the influence of the Hartmann bottom boundary layer in a three-dimensional(3D)MHD experiment in a square bottomed cell.The basic flow in this fluid system is a square eddy flow exhibiting a network of N~2 vortices rotating alternately in clockwise and anticlockwise directions.When N is odd,the instability of the flow gives rise to secondary steady-state flows and secondary time-periodic flows,exhibiting similar characteristics to those observed when N=3.For this reason,this study focuses on the instability of the square eddy flow of nine vortices.It is shown that there exist eight bi-critical values corresponding to the existence of eight neutral eigenfunction spaces.Especially,there exist non-real neutral eigenfunctions,which produce secondary time-periodic flows exhibiting vortices merging in an oscillatory manner.This Hopf bifurcation phenomenon has not been observed in earlier investigations.
基金supported by the NSFC Grant No.11872210 and Grant No.MCMS-I-0120G01Chi-Wang Shu:Research is supported by the AFOSR Grant FA9550-20-1-0055 and the NSF Grant DMS-2010107Jianxian Qiu:Research is supported by the NSFC Grant No.12071392.
文摘In this paper, we design high-order Runge-Kutta discontinuous Galerkin (RKDG) methods with multi-resolution weighted essentially non-oscillatory (multi-resolution WENO) limiters to compute compressible steady-state problems on triangular meshes. A troubled cell indicator extended from structured meshes to unstructured meshes is constructed to identify triangular cells in which the application of the limiting procedures is required. In such troubled cells, the multi-resolution WENO limiting methods are used to the hierarchical L^(2) projection polynomial sequence of the DG solution. Through using the RKDG methods with multi-resolution WENO limiters, the optimal high-order accuracy can be gradually reduced to first-order in the triangular troubled cells, so that the shock wave oscillations can be well suppressed. In steady-state simulations on triangular meshes, the numerical residual converges to near machine zero. The proposed spatial reconstruction methods enhance the robustness of classical DG methods on triangular meshes. The good results of these RKDG methods with multi-resolution WENO limiters are verified by a series of two-dimensional steady-state problems.
文摘This study aims to develop a mathematical analysis for one-dimensional modeling of a radial flow through a production well drilled in a confined aquifer, in the case of steady-state flow conditions. An analytical solution has derived from that expression for estimation of drawdowns according to different flowrates. Through that process, the evaluation of static pressure, the calculation of hydraulic charge due to the waterflow through the well is evaluated, the drawdowns curves are drawn and at last, the obtained curves are analyzed. The curves obtained for the different flow rates have an asymptotic direction, the axis of the hydraulic charges. The variation of the hydraulic charge depends on the radial distance for different flow rates. The P point, is a common point of all curves obtained for different production flowrates in the well. This point is where the well production flowrate is optimum for the optimal hydraulic charge.
基金Supported by National Basic Research Program of China(973 Program,Grant No.2014CB046403)National Key Technology R&D Program of the Twelfth Five-year Plan of China(Grant No.2013BAF07B01)
文摘The current research about the flow ripple of axial piston pump mainly focuses on the effect of the structure of parts on the flow ripple. Therein, the structure of parts are usually designed and optimized at rated working conditions. However, the pump usually has to work in large-scale and time-variant working conditions. Therefore, the flow ripple characteristics of pump and analysis for its test accuracy with respect to variant steady-state conditions and transient conditions in a wide range of operating parameters are focused in this paper. First, a simulation model has been constructed, which takes the kinematics of oil film within friction pairs into account for higher accuracy. Afterwards, a test bed which adopts Secondary Source Method is built to verify the model. The simulation and tests results show that the angular position of the piston, corresponding to the position where the peak flow ripple is produced, varies with the different pressure. The pulsating amplitude and pulsation rate of flow ripple increase with the rise of pressure and the variation rate of pressure. For the pump working at a constant speed, the flow pulsation rate decreases dramatically with the increasing speed when the speed is less than 27.78% of the maximum speed, subsequently presents a small decrease tendency with the speed further increasing. With the rise of the variation rate of speed, the pulsating amplitude and pulsation rate of flow ripple increase. As the swash plate angle augments, the pulsating amplitude of flow ripple increases, nevertheless the flow pulsation rate decreases. In contrast with the effect of the variation of pressure, the test accuracy of flow ripple is more sensitive to the variation of speed. It makes the test accuracy above 96.20% available for the pulsating amplitude of pressure deviating within a range of ~6% from the mean pressure. However, with a variation of speed deviating within a range of ±2% from the mean speed, the attainable test accuracy of flow ripple is above 93.07%. The model constructed in this research proposes a method to determine the flow ripple characteristics of pump and its attainable test accuracy under the large-scale and time-variant working conditions. Meanwhile, a discussion about the variation of flow ripple and its obtainable test accuracy with the conditions of the pump working in wide operating ranges is given as well.
基金supported by the National Natural Science Foundation of China(52304067,62273213)the Natural Science Foundation of Shandong Province of China(ZR2021QE073)+1 种基金the Natural Science Foundation of Shandong Province for Innovation and Development Joint Funds(ZR2022LZH001)the China Postdoctoral Science Foundation(2023M732111)。
文摘Liquid hydrogen storage and transportation is an effective method for large-scale transportation and utilization of hydrogen energy. Revealing the flow mechanism of cryogenic working fluid is the key to optimize heat exchanger structure and hydrogen liquefaction process(LH2). The methods of cryogenic visualization experiment, theoretical analysis and numerical simulation are conducted to study the falling film flow characteristics with the effect of co-current gas flow in LH2spiral wound heat exchanger.The results show that the flow rate of mixed refrigerant has a great influence on liquid film spreading process, falling film flow pattern and heat transfer performance. The liquid film of LH2mixed refrigerant with column flow pattern can not uniformly and completely cover the tube wall surface. As liquid flow rate increases, the falling film flow pattern evolves into sheet-column flow and sheet flow, and liquid film completely covers the surface of tube wall. With the increase of shear effect of gas-phase mixed refrigerant in the same direction, the liquid film gradually becomes unstable, and the flow pattern eventually evolves into a mist flow.
文摘Gassolid hydrodynamic steadystate operation is the operating basis in a chemical looping dualreactor system.This study reported the experimental results on the steadystate operation characteristics of gassolid flow in a 15.5 m high dual circulating fluidized bed(CFB)cold test system.The effects of superficial gas velocity,static bed material height and solid returning modes on the steadystate operation characteristics between the two CFBs were investigated.Results suggest that the solid distributions in the dual CFB test system was mainly determined by the superficial gas velocity and larger solid inventory may help to improve the solid distributions.Besides,crossreturning mode coupled with selfreturning is good for steadystate running in the dualreactor test system.
基金Project supported by the National Natural Science Foundation of China (Grant No. 12072118)the Natural Science Funds for Distinguished Young Scholar of the Fujian Province, China (Grant No. 2021J06024)the Project for Youth Innovation Fund of Xiamen, China (Grant No. 3502Z20206005)。
文摘The majority of nonlinear stochastic systems can be expressed as the quasi-Hamiltonian systems in science and engineering. Moreover, the corresponding Hamiltonian system offers two concepts of integrability and resonance that can fully describe the global relationship among the degrees-of-freedom(DOFs) of the system. In this work, an effective and promising approximate semi-analytical method is proposed for the steady-state response of multi-dimensional quasi-Hamiltonian systems. To be specific, the trial solution of the reduced Fokker–Plank–Kolmogorov(FPK) equation is obtained by using radial basis function(RBF) neural networks. Then, the residual generated by substituting the trial solution into the reduced FPK equation is considered, and a loss function is constructed by combining random sampling technique. The unknown weight coefficients are optimized by minimizing the loss function through the Lagrange multiplier method. Moreover, an efficient sampling strategy is employed to promote the implementation of algorithms. Finally, two numerical examples are studied in detail, and all the semi-analytical solutions are compared with Monte Carlo simulations(MCS) results. The results indicate that the complex nonlinear dynamic features of the system response can be captured through the proposed scheme accurately.
基金support by FEDER-Fundo Europeu de Desenvolvimento Regional,through COMPETE 2020-Programa Operational Fatores de Competitividade,and the National Funds through FCT-Fundacao para a Ciencia e a Tecnologia,project no.UID/FIS/04650/2019support by FEDER-Fundo Europeu de Desenvolvimento Regional,through COMPETI E 2020-Programa Operacional Fatores de Competitividade,and the National Funds through FCT-Fundacao para a Ciencia e a Tecnologia,project no.POCI-01-0145-FEDER-028118
文摘We propose an adaptive stencil construction for high-order accurate finite volume schemes a posteriori stabilized devoted to solve one-dimensional steady-state hyperbolic equations.High accuracy(up to the sixth-order presently)is achieved,thanks to polynomial recon-structions while stability is provided with an a posteriori MOOD method which controls the cell polynomial degree for eliminating non-physical oscillations in the vicinity of dis-continuities.We supplemented this scheme with a stencil construction allowing to reduce even further the numerical dissipation.The stencil is shifted away from troubles(shocks,discontinuities,etc.)leading to less oscillating polynomial reconstructions.Experimented on linear,Burgers',and Euler equations,we demonstrate that the adaptive stencil technique manages to retrieve smooth solutions with optimal order of accuracy but also irregular ones without spurious oscillations.Moreover,we numerically show that the approach allows to reduce the dissipation still maintaining the essentially non-oscillatory behavior.
基金This study was supported by the National Natural Science Foundation of China(U22B2075,52274056,51974356).
文摘A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes due to stress sensitivity, which plays a crucial role in controlling pressure propagation and oil flow. This paper proposes a multi-scale coupled flow mathematical model of matrix nanopores, induced fractures, and hydraulic fractures. In this model, the micro-scale effects of shale oil flow in fractal nanopores, fractal induced fracture network, and stress sensitivity of multi-scale media are considered. We solved the model iteratively using Pedrosa transform, semi-analytic Segmented Bessel function, Laplace transform. The results of this model exhibit good agreement with the numerical solution and field production data, confirming the high accuracy of the model. As well, the influence of stress sensitivity on permeability, pressure and production is analyzed. It is shown that the permeability and production decrease significantly when induced fractures are weakly supported. Closed induced fractures can inhibit interporosity flow in the stimulated reservoir volume (SRV). It has been shown in sensitivity analysis that hydraulic fractures are beneficial to early production, and induced fractures in SRV are beneficial to middle production. The model can characterize multi-scale flow characteristics of shale oil, providing theoretical guidance for rapid productivity evaluation.
