A new physical structure of vortical flow, i.e., tubular limiting stream surface(TLSS), is reported. It is defined as a general mathematical structure for the physical flow field in the neighborhood of a singularity, ...A new physical structure of vortical flow, i.e., tubular limiting stream surface(TLSS), is reported. It is defined as a general mathematical structure for the physical flow field in the neighborhood of a singularity, and has a close relationship with limit cycles.The TLSS is a tornado-like structure, which separates a vortex into two regions, i.e., the inner region near the vortex axis and the outer region further away from the vortex axis.The flow particles in these two regions can approach to(or leave) the TLSS, but never could reach it.展开更多
Upon the conservation of mass, momentum and energy, volume fraction and surface penetrative rate were employed to modify the conservative equations to simulate the effect of blockages on fluid flows and heat transfer....Upon the conservation of mass, momentum and energy, volume fraction and surface penetrative rate were employed to modify the conservative equations to simulate the effect of blockages on fluid flows and heat transfer. These equations were solved numerically with the finite differential method and the primitive variable approach. This method uses staggered grid and pressure correction schemes. A computer code FASTOR3D integrated the aforementioned algorithm. The preliminary results have been compared with conventional benchmark solutions. With auxiliary software DV, the numerical results were visualized in colorful images to demonstrate the variation of flow patterns and temperature profiles during the transient process. The results of the simulation code for the fluid flows and heat transfer in the sodium pool of a fast breeder reactor are acceptable.展开更多
An improved algorithm for computing multiphase flows is presented in which the multimaterial Moment-of-Fluid(MOF)algorithm for multiphase flows,initially described by Li et al.(2015),is enhanced addressing existing MO...An improved algorithm for computing multiphase flows is presented in which the multimaterial Moment-of-Fluid(MOF)algorithm for multiphase flows,initially described by Li et al.(2015),is enhanced addressing existing MOF difficulties in computing solutions to problems in which surface tension forces are crucial for understanding salient flow mechanisms.The Continuous MOF(CMOF)method is motivated in this article.The CMOF reconstruction method inherently removes the"checkerboard instability"that persists when using the MOF method on surface tension driven multiphase(multimaterial)flows.The CMOF reconstruction algorithm is accelerated by coupling the CMOF method to the level set method and coupling the CMOF method to a decision tree machine learning(ML)algorithm.Multiphase flow examples are shown in the two-dimensional(2D),three-dimensional(3D)axisymmetric"RZ",and 3D coordinate systems.Examples include two material and three material multiphase flows:bubble formation,the impingement of a liquid jet on a gas bubble in a cryogenic fuel tank,freezing,and liquid lens dynamics.展开更多
This study is concerned with the three-dimensional(3D)stagnation-point for the mixed convection flow past a vertical surface considering the first-order and secondorder velocity slips.To the authors’knowledge,this is...This study is concerned with the three-dimensional(3D)stagnation-point for the mixed convection flow past a vertical surface considering the first-order and secondorder velocity slips.To the authors’knowledge,this is the first study presenting this very interesting analysis.Nonlinear partial differential equations for the flow problem are transformed into nonlinear ordinary differential equations(ODEs)by using appropriate similarity transformation.These ODEs with the corresponding boundary conditions are numerically solved by utilizing the bvp4c solver in MATLAB programming language.The effects of the governing parameters on the non-dimensional velocity profiles,temperature profiles,skin friction coefficients,and the local Nusselt number are presented in detail through a series of graphs and tables.Interestingly,it is reported that the reduced skin friction coefficient decreases for the assisting flow situation and increases for the opposing flow situation.The numerical computations of the present work are compared with those from other research available in specific situations,and an excellent consensus is observed.Another exciting feature for this work is the existence of dual solutions.An important remark is that the dual solutions exist for both assisting and opposing flows.A linear stability analysis is performed showing that one solution is stable and the other solution is not stable.We notice that the mixed convection and velocity slip parameters have strong effects on the flow characteristics.These effects are depicted in graphs and discussed in this paper.The obtained results show that the first-order and second-order slip parameters have a considerable effect on the flow,as well as on the heat transfer characteristics.展开更多
Cattaneo-Christov heat and mass flux models are considered rather than Fourier and Fick laws due to the presence of thermal and concentration transport hyperbolic phenomena. The generalized form of the Navier-Stokes m...Cattaneo-Christov heat and mass flux models are considered rather than Fourier and Fick laws due to the presence of thermal and concentration transport hyperbolic phenomena. The generalized form of the Navier-Stokes model is considered in hydromagnetic flow. Three-dimensional (3D) unsteady fluid motion is generated by the periodic oscillations of a rotating disk. Similarity transformations are used to obtain the normalized fluid flow model. The successive over relaxation (SOR) method with finite difference schemes are accomplished for the numerical solution of the obtained partial differential non-linear system. The flow features of the velocity, microrotation, temperature, and concentration fields are discussed in pictorial forms for various physical flow parameters. The couple stresses and heat and mass transfer rates for different physical quantities are explained via tabular forms. For better insight of the physical fluid model, 3D fluid phenomena and two-dimensional (2D) contours are also plotted. The results show that the micropolar fluids contain microstructure having non-symmetric stress tensor and are useful in lubrication theory. Moreover, the thermal and concentration waves in Cattaneo-Christov models have a significance role in the laser heating and enhancement in thermal conductivity.展开更多
A three-dimensional (3D) predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is ...A three-dimensional (3D) predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.展开更多
It is important to calibrate micro-parameters for applying partied flow code(PFC)to study mechanical characteristics and failure mechanism of rock materials.Uniform design method is firstly adopted to determine the mi...It is important to calibrate micro-parameters for applying partied flow code(PFC)to study mechanical characteristics and failure mechanism of rock materials.Uniform design method is firstly adopted to determine the microscopic parameters of parallel-bonded particle model for three-dimensional discrete element particle flow code(PFC3D).Variation ranges of microscopic of the microscopic parameters are created by analyzing the effects of microscopic parameters on macroscopic parameters(elastic modulus E,Poisson ratio v,uniaxial compressive strengthσc,and ratio of crack initial stress to uniaxial compressive strengthσci/σc)in order to obtain the actual uniform design talbe.The calculation equations of the microscopic and macroscopic parameters of rock materials can be established by the actual uniform design table and the regression analysis and thus the PFC3D microscopic parameters can be quantitatively determined.The PFC3D simulated results of the intact and pre-cracked rock specimens under uniaxial and triaxial compressions(including the macroscopic mechanical parameters,stress−strain curves and failure process)are in good agreement with experimental results,which can prove the validity of the calculation equations of microscopic and macroscopic parameters.展开更多
Streamwise vortex instability is one of the most potent mechanisms for the transition of the three-dimensional boundary layers.By using the global stability analysis methods,stability characteristics of the leeward vo...Streamwise vortex instability is one of the most potent mechanisms for the transition of the three-dimensional boundary layers.By using the global stability analysis methods,stability characteristics of the leeward vortex over a blunt cone with an angle of attack under a typical wind tunnel condition are studied and are compared to the case with a smaller wall temperature ratio(corresponding to a flight condition).The vortical structure features inward and outward vortices,similar to that in the flight con-dition.Unlike the flight condition,the outward vortices appear stronger than the inward vortices,resulting in stronger outer-mode instabilities.Although the inner mode is heavily stabilized compared to the flight condition,it can still radiate apparent acoustics.The acoustic sources are computed based on Lighthill's acoustic analogy,showing that the entropy term measuring the deviation from the isentropic relation is dominant.While Mack second mode is shown to most likely trigger the transition in the flight condition,it is absent in the wind tunnel condition,and a shear-layer mode turns out to be the most dangerous instead.Moreover,the instability frequencies and growth rates of the wind tunnel case are much smaller than those of the flight case,indicating that wall heating may stabilize the leeward vortices.展开更多
基金Project supported by the National Natural Science Foundation of China(Nos.11372340 and 11732016)
文摘A new physical structure of vortical flow, i.e., tubular limiting stream surface(TLSS), is reported. It is defined as a general mathematical structure for the physical flow field in the neighborhood of a singularity, and has a close relationship with limit cycles.The TLSS is a tornado-like structure, which separates a vortex into two regions, i.e., the inner region near the vortex axis and the outer region further away from the vortex axis.The flow particles in these two regions can approach to(or leave) the TLSS, but never could reach it.
文摘Upon the conservation of mass, momentum and energy, volume fraction and surface penetrative rate were employed to modify the conservative equations to simulate the effect of blockages on fluid flows and heat transfer. These equations were solved numerically with the finite differential method and the primitive variable approach. This method uses staggered grid and pressure correction schemes. A computer code FASTOR3D integrated the aforementioned algorithm. The preliminary results have been compared with conventional benchmark solutions. With auxiliary software DV, the numerical results were visualized in colorful images to demonstrate the variation of flow patterns and temperature profiles during the transient process. The results of the simulation code for the fluid flows and heat transfer in the sodium pool of a fast breeder reactor are acceptable.
基金supported by the National Aeronautics and Space Administration under grant number 80NSSC20K0352.
文摘An improved algorithm for computing multiphase flows is presented in which the multimaterial Moment-of-Fluid(MOF)algorithm for multiphase flows,initially described by Li et al.(2015),is enhanced addressing existing MOF difficulties in computing solutions to problems in which surface tension forces are crucial for understanding salient flow mechanisms.The Continuous MOF(CMOF)method is motivated in this article.The CMOF reconstruction method inherently removes the"checkerboard instability"that persists when using the MOF method on surface tension driven multiphase(multimaterial)flows.The CMOF reconstruction algorithm is accelerated by coupling the CMOF method to the level set method and coupling the CMOF method to a decision tree machine learning(ML)algorithm.Multiphase flow examples are shown in the two-dimensional(2D),three-dimensional(3D)axisymmetric"RZ",and 3D coordinate systems.Examples include two material and three material multiphase flows:bubble formation,the impingement of a liquid jet on a gas bubble in a cryogenic fuel tank,freezing,and liquid lens dynamics.
基金Project supported by the Executive Agency for Higher Education Research Development and Innovation Funding of Romania(No.PN-III-P4-PCE-2021-0993)。
文摘This study is concerned with the three-dimensional(3D)stagnation-point for the mixed convection flow past a vertical surface considering the first-order and secondorder velocity slips.To the authors’knowledge,this is the first study presenting this very interesting analysis.Nonlinear partial differential equations for the flow problem are transformed into nonlinear ordinary differential equations(ODEs)by using appropriate similarity transformation.These ODEs with the corresponding boundary conditions are numerically solved by utilizing the bvp4c solver in MATLAB programming language.The effects of the governing parameters on the non-dimensional velocity profiles,temperature profiles,skin friction coefficients,and the local Nusselt number are presented in detail through a series of graphs and tables.Interestingly,it is reported that the reduced skin friction coefficient decreases for the assisting flow situation and increases for the opposing flow situation.The numerical computations of the present work are compared with those from other research available in specific situations,and an excellent consensus is observed.Another exciting feature for this work is the existence of dual solutions.An important remark is that the dual solutions exist for both assisting and opposing flows.A linear stability analysis is performed showing that one solution is stable and the other solution is not stable.We notice that the mixed convection and velocity slip parameters have strong effects on the flow characteristics.These effects are depicted in graphs and discussed in this paper.The obtained results show that the first-order and second-order slip parameters have a considerable effect on the flow,as well as on the heat transfer characteristics.
文摘Cattaneo-Christov heat and mass flux models are considered rather than Fourier and Fick laws due to the presence of thermal and concentration transport hyperbolic phenomena. The generalized form of the Navier-Stokes model is considered in hydromagnetic flow. Three-dimensional (3D) unsteady fluid motion is generated by the periodic oscillations of a rotating disk. Similarity transformations are used to obtain the normalized fluid flow model. The successive over relaxation (SOR) method with finite difference schemes are accomplished for the numerical solution of the obtained partial differential non-linear system. The flow features of the velocity, microrotation, temperature, and concentration fields are discussed in pictorial forms for various physical flow parameters. The couple stresses and heat and mass transfer rates for different physical quantities are explained via tabular forms. For better insight of the physical fluid model, 3D fluid phenomena and two-dimensional (2D) contours are also plotted. The results show that the micropolar fluids contain microstructure having non-symmetric stress tensor and are useful in lubrication theory. Moreover, the thermal and concentration waves in Cattaneo-Christov models have a significance role in the laser heating and enhancement in thermal conductivity.
基金supported by the Yunnan Provincial Applied Basic Research Program of China(No. KKSY201207019)
文摘A three-dimensional (3D) predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.
基金Projects(51474251,51874351)supported by the National Natural Science Foundation,China。
文摘It is important to calibrate micro-parameters for applying partied flow code(PFC)to study mechanical characteristics and failure mechanism of rock materials.Uniform design method is firstly adopted to determine the microscopic parameters of parallel-bonded particle model for three-dimensional discrete element particle flow code(PFC3D).Variation ranges of microscopic of the microscopic parameters are created by analyzing the effects of microscopic parameters on macroscopic parameters(elastic modulus E,Poisson ratio v,uniaxial compressive strengthσc,and ratio of crack initial stress to uniaxial compressive strengthσci/σc)in order to obtain the actual uniform design talbe.The calculation equations of the microscopic and macroscopic parameters of rock materials can be established by the actual uniform design table and the regression analysis and thus the PFC3D microscopic parameters can be quantitatively determined.The PFC3D simulated results of the intact and pre-cracked rock specimens under uniaxial and triaxial compressions(including the macroscopic mechanical parameters,stress−strain curves and failure process)are in good agreement with experimental results,which can prove the validity of the calculation equations of microscopic and macroscopic parameters.
基金supported by the National Natural Science Foundation of China(Grant No.92052301).
文摘Streamwise vortex instability is one of the most potent mechanisms for the transition of the three-dimensional boundary layers.By using the global stability analysis methods,stability characteristics of the leeward vortex over a blunt cone with an angle of attack under a typical wind tunnel condition are studied and are compared to the case with a smaller wall temperature ratio(corresponding to a flight condition).The vortical structure features inward and outward vortices,similar to that in the flight con-dition.Unlike the flight condition,the outward vortices appear stronger than the inward vortices,resulting in stronger outer-mode instabilities.Although the inner mode is heavily stabilized compared to the flight condition,it can still radiate apparent acoustics.The acoustic sources are computed based on Lighthill's acoustic analogy,showing that the entropy term measuring the deviation from the isentropic relation is dominant.While Mack second mode is shown to most likely trigger the transition in the flight condition,it is absent in the wind tunnel condition,and a shear-layer mode turns out to be the most dangerous instead.Moreover,the instability frequencies and growth rates of the wind tunnel case are much smaller than those of the flight case,indicating that wall heating may stabilize the leeward vortices.