In this paper, we derive a general formula of the flow equtions for the Harry-Dym hierarchy. And three applications in n-reduction, (2 + 1)-dimensional generalization, and Kupershmidt reduction, are considered.
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.展开更多
In this paper, the convergence turbulent flow equations are considered. By rates of solutions to the three-dimensional combining the LP-Lq estimate for the linearized equations and an elaborate energy method, the conv...In this paper, the convergence turbulent flow equations are considered. By rates of solutions to the three-dimensional combining the LP-Lq estimate for the linearized equations and an elaborate energy method, the convergence rates are obtained in various norms for the solution to the equilibrium state in the whole space when the initial perturbation of the equilibrium state is small in the H3-framework. More precisely, the optimal convergence rates of the solutions and their first-order derivatives in the L2-norm are obtained when the LP-norm of the perturbation is bounded for some p ε [1, 6).展开更多
hi this paper. the author uses the theory of fluid mechanics. dynamics of fluids in Porous media. gas seepage flow in coal seams and combines the tests in the laboratory with the actual coal infusion to have an inves...hi this paper. the author uses the theory of fluid mechanics. dynamics of fluids in Porous media. gas seepage flow in coal seams and combines the tests in the laboratory with the actual coal infusion to have an investigating and study from the theory to the mechanism of coal infusion to wet coal seams. through the analysis to the process of coal infusion the author builds up the mathematical models and has a detailed discussion to the boundary conditions of coal infusion. Because the equation sets to describe coal infusion are non-linear. we have made a simplification to them to use the dimension analysis theory by leading into the non-dimensions of water pressure of coal infusion, seepage flow rate. increment of coal seam moisture and so on Besides the analytic and approximate solutions have also been discussed. At last. we use the scientific research item of the actual coal infusion to illustrate the effects and importance of the theory to direct actual coal infusion and its designs.展开更多
We develop a numerical solution algorithm of the nonlinear potential flow equations with the nonlinear free surface boundary condition.A finite difference method with a predictor-corrector method is applied to solve t...We develop a numerical solution algorithm of the nonlinear potential flow equations with the nonlinear free surface boundary condition.A finite difference method with a predictor-corrector method is applied to solve the nonlinear potential flow equations in a two-dimensional (2D) tank.The irregular tank is mapped onto a fixed square domain with rectangular cells through a proper mapping function.A staggered mesh system is adopted in a 2D tank to capture the wave elevation of the transient fluid.The finite difference method with a predictor-corrector scheme is applied to discretize the nonlinear dynamic boundary condition and nonlinear kinematic boundary condition.We present the numerical results of wave elevations from small to large amplitude waves with free oscillation motion,and the numerical solutions of wave elevation with horizontal excited motion.The beating period and the nonlinear phenomenon are very clear.The numerical solutions agree well with the analytical solutions and previously published results.展开更多
Three dimensional Euler equations are solved in the finite volume form with van Leer's flux vector splitting technique. Block matrix is inverted by Gauss-Seidel iteration in two dimensional plane while strongly im...Three dimensional Euler equations are solved in the finite volume form with van Leer's flux vector splitting technique. Block matrix is inverted by Gauss-Seidel iteration in two dimensional plane while strongly implicit alternating sweeping is implemented in the direction of the third dimension. Very rapid convergence rate is obtained with CFL number reaching the order of 100. The memory resources can be greatly saved too. It is verified that the reflection boundary condition can not be used with flux vector splitting since it will produce too large numerical dissipation. The computed flow fields agree well with experimental results. Only one or two grid points are there within the shock transition zone.展开更多
A complete boundary integral formulation for steady compressible inviscid flows governed by nonlinear equations is established by using ρV as variable. Thus, the dimensionality of the problem to be solved is reduced ...A complete boundary integral formulation for steady compressible inviscid flows governed by nonlinear equations is established by using ρV as variable. Thus, the dimensionality of the problem to be solved is reduced by one and the computational mesh to be generated is needed only on the boundary of the domain.展开更多
AbstractThis paper gives H-O type grid generation and applies finite-volume method forwing-body flow.A 3-D Euler code for wing-body has been developed. We improve in thiscode the distribution of the number of grid poi...AbstractThis paper gives H-O type grid generation and applies finite-volume method forwing-body flow.A 3-D Euler code for wing-body has been developed. We improve in thiscode the distribution of the number of grid points over wing and body surfaces,and keepthe training edge of the wing being one of the grid lines even when the edge has sweepbackor sweep forward.The code developed can provide not only the spanwise pressure distribu-tion as usual ou cross-flow planes,but also the chordwise Pressure distributions.The com-putation for NASA TN D-712 wing-body model with the present code shows that thecomputed pressure distributions are in very good agreement with the experiment.展开更多
This paper presents analytieal solutions to the partial differential equations for unsteady flow of the second-order fluid and Maxwell fluid in tube by using the integral transform method. It can be used to analyse th...This paper presents analytieal solutions to the partial differential equations for unsteady flow of the second-order fluid and Maxwell fluid in tube by using the integral transform method. It can be used to analyse the behaviour of axial velocity and shear stress for unsteady flow of nun-Newtonian visco-elastie fluids in tube, and to provide a theoretical base for the projection of pipe-line engineering.展开更多
Inthis paper, each of the two phases in dense two-phase flow is considered as continuous medium and the fundamental equations for two-phase flow arc described in Eulerian form. The generalized constitutive relation of...Inthis paper, each of the two phases in dense two-phase flow is considered as continuous medium and the fundamental equations for two-phase flow arc described in Eulerian form. The generalized constitutive relation of the Bingham fluid is applied to the dispersed phase with the analysis oj physical mechanism of dense two-phase flow. The shearing stress of dispersed phase at a wall is used to give a boundary condition. Then a mathematical model for dense two-phase flow is obtained. In addition, the expressions of shearing stress of dispersed phase at a wall is derived according to the fundamental model of the friclional collision between dispersed-plutse particles and the wall.展开更多
Fluid flow at nanoscale is closely related to many areas in nature and technology(e.g.,unconventional hydrocarbon recovery,carbon dioxide geo-storage,underground hydrocarbon storage,fuel cells,ocean desalination,and b...Fluid flow at nanoscale is closely related to many areas in nature and technology(e.g.,unconventional hydrocarbon recovery,carbon dioxide geo-storage,underground hydrocarbon storage,fuel cells,ocean desalination,and biomedicine).At nanoscale,interfacial forces dominate over bulk forces,and nonlinear effects are important,which significantly deviate from conventional theory.During the past decades,a series of experiments,theories,and simulations have been performed to investigate fluid flow at nanoscale,which has advanced our fundamental knowledge of this topic.However,a critical review is still lacking,which has seriously limited the basic understanding of this area.Therefore herein,we systematically review experimental,theoretical,and simulation works on single-and multi-phases fluid flow at nanoscale.We also clearly point out the current research gaps and future outlook.These insights will promote the significant development of nonlinear flow physics at nanoscale and will provide crucial guidance on the relevant areas.展开更多
This paper gives an overview on nonlinear porous flow in low permeability porous media, reveals the microscopic mechanisms of flows, and clarifies properties of porous flow fluids. It shows that, deviating from Darcy...This paper gives an overview on nonlinear porous flow in low permeability porous media, reveals the microscopic mechanisms of flows, and clarifies properties of porous flow fluids. It shows that, deviating from Darcy's linear law, the porous flow characteristics obey a nonlinear law in a low-permeability porous medium, and the viscosity of the porous flow fluid and the permeability values of water and oil are not constants. Based on these characters, a new porous flow model, which can better describe low permeability reservoir~ is established. This model can describe various patterns of porous flow, as Darcy's linear law does. All the parameters involved in the model, having definite physical meanings, can be obtained directly from the experiments.展开更多
A series of simple axisymmetric compression tests were carried out on the computer servo controlled Gleeble 1 500 machine when strain rates ranged between 0.05 25 s -1 and deformation temperature 300 500 ℃. The resul...A series of simple axisymmetric compression tests were carried out on the computer servo controlled Gleeble 1 500 machine when strain rates ranged between 0.05 25 s -1 and deformation temperature 300 500 ℃. The results show that flow stress is related to the Zener Hollonom parameter Z and strain, as well as the static recrystallization fraction between passes during multipass hot deformation of 5182 aluminum alloy. Hence, a modified exponential flow stress equation was presented by considering the values of ln A and β as functions of strain, and by using the uniform softening method and incorporating the static recrystallization fraction between passes to consider the effects of residual strain during multipass hot rolling of 5182 aluminum alloy. The validity of the equation was examined by a typical non isothermal multipass deformation test.展开更多
Underwater gliders are recent innovative types of autonomous underwater vehicles (AUVs) used in ocean exploration and observation. They adjust their buoyancy to dive and to return to the ocean surface. During the ch...Underwater gliders are recent innovative types of autonomous underwater vehicles (AUVs) used in ocean exploration and observation. They adjust their buoyancy to dive and to return to the ocean surface. During the change of altitude, they use the hydrodynamic forces developed by their wings to move forward. Their flights are controlled by changing the position of their centers of gravity and their buoyancy to adjust their trim and heel angles. For better flight control, the understanding of the hydrodynamic behavior and the flight mechanics of the underwater glider is necessary. A 6-DOF motion simulator is coupled with an unsteady potential flow model for this purpose. In some specific cases, the numerical study demonstrates that an inappropriate stabilizer dimension can cause counter-steering behavior. The simulator can be used to improve the automatic flight control. It can also be used for the hydrodynamic design optimization of the devices.展开更多
In this paper, the fluid flow differential equation based on the homogenous reservoirs model is first reviewed. Then a theorem about the formal similarity of solutions in the Laplace space with outer boundary conditio...In this paper, the fluid flow differential equation based on the homogenous reservoirs model is first reviewed. Then a theorem about the formal similarity of solutions in the Laplace space with outer boundary conditions and inner boundary condition is presented and proved. Lastly, a corollary of our theorem is given particularly on inner boundary. The obtained results are very helpful for understanding inherent laws of relevant engineering science and designing practical analysis software.展开更多
This paper focuses on a two-dimensional bidirectional pedestrian flow model which involves the next-nearest-neighbor effect. The stability condition and the Korteweg-de Vries (KdV) equation are derived to describe t...This paper focuses on a two-dimensional bidirectional pedestrian flow model which involves the next-nearest-neighbor effect. The stability condition and the Korteweg-de Vries (KdV) equation are derived to describe the density wave of pedestrian congestion by linear stability and nonlinear analysis. Through theoretical analysis, the soliton solution is obtained.展开更多
This paper defines new kinds of functions——the conjugate axisymmetric poteptial functions. With the aid of them, we can prove the completeness of the solutions of the generalized axisymmetric Stokes flow equation wi...This paper defines new kinds of functions——the conjugate axisymmetric poteptial functions. With the aid of them, we can prove the completeness of the solutions of the generalized axisymmetric Stokes flow equation without the condition on the convexity of the domain.展开更多
Debris flow simulations are useful for predicting the sediment supplied to watersheds from upstream areas. However, the topographic conditions upstream are more complicated than those downstream and the relationship b...Debris flow simulations are useful for predicting the sediment supplied to watersheds from upstream areas. However, the topographic conditions upstream are more complicated than those downstream and the relationship between the topographic conditions and debris flow initiation is not well understood. This study compared the use of several entrainment rate equations in numerical simulations of debris flows to examine the effect of topographic conditions on the flow. One-dimensional numerical simulations were performed based on the shallow water equations and three entrainment rate equations were tested. These entrainment rate equations were based on the same idea that erosion and the deposition of debris flows occur via the difference between the equilibrium and current conditions of debris flows, while they differed in the expression of the concentration, channel angle, and sediment amount. The comparison was performed using a straight channel with various channel angles and a channel with a periodically undulating surface. The three entrainment rate equations gave different amounts of channel bed degradation and hydrographs for a straight channel with a channel angle greater than 21° when water was supplied from upstream at a steady rate. The difference was caused by the expression of the entrainment rate equations. For channels with little undulation, the numerical simulations gave results almost identical to those for straight channels with the same channel angle. However, for channels with large undulations, the hydrographs differed from those for straight channels with the same channel angle when the channel angle was less than 21°. Rapid erosion occurred and the hydrograph showed a significant peak, especially in cases using the entrainment equation expressed by channel angle. This was caused by the effects of the steep undulating sections, since the effect increased with the magnitude of the undulation, suggesting that a debris flow in an upstream area develops differently according to the topographic conditions. These results also inferred that numerical simulations of debris flow can differ depending on the spatial resolution of the simulation domain, as the resolution determines the reproducibility of the undulations.展开更多
Through transformations, the time-dependent boundary condition on the airfoil contour and the boundary condition at infinity are brought fixed to the boundaries of a finite domain. The boundary conditions can thus be ...Through transformations, the time-dependent boundary condition on the airfoil contour and the boundary condition at infinity are brought fixed to the boundaries of a finite domain. The boundary conditions can thus be satisfied exactly without increasing the computational time. The novel scheme is useful for computing transonic, strong disturbance, unsteady flows with high reduced frequencies. The scheme makes use of curvefitted orthogonal meshes and the lattice control technique to obtain the optimal grid distribution. The numerical results are satisfactory.展开更多
The finite difference method (FDM) is applied in the present paper to solve the unsteady NHS equations for incompressible fluids. ADI and SLOR methods are served for the vorticity equation and the Poisson equation for...The finite difference method (FDM) is applied in the present paper to solve the unsteady NHS equations for incompressible fluids. ADI and SLOR methods are served for the vorticity equation and the Poisson equation for ψ respectively. The upwind scheme is used for the convective terms. The moving boundary conditions are specially treated, and the effects of outlet conditions on the flow field are abo examined. Numerical results obtained show that the spoiler's oscillation induces forming, growing and shedding of the vortices. The shedding frequency of vortices is equal to that of the spoiler's oscillation. The forced unsteady separated flows under the present investigation depend mainly on the reduced frequency. At low reduced frequency, the vortices shed from the spoiler interact weakly with each other, and move downstream at an almost uniform speed of 038 V∞. At high reduced frequency, the interaction between the adjacent vortices strengthens. They close up to and rotate around each other, and eventually, merge into one vortex.展开更多
基金Supported by the Natural Science Foundation of China under Grant Nos.10671187 and 10971 109Supported by Program for New Century Excellent Talents in Universities under Grant No.NECT-08-0515
文摘In this paper, we derive a general formula of the flow equtions for the Harry-Dym hierarchy. And three applications in n-reduction, (2 + 1)-dimensional generalization, and Kupershmidt reduction, are considered.
基金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.
基金supported by the National Natural Science Foundation of China(Nos.11071057 and 11271052)the Special Fund Project of Mathematical Tian Yuan Fund(No.11226029)
文摘In this paper, the convergence turbulent flow equations are considered. By rates of solutions to the three-dimensional combining the LP-Lq estimate for the linearized equations and an elaborate energy method, the convergence rates are obtained in various norms for the solution to the equilibrium state in the whole space when the initial perturbation of the equilibrium state is small in the H3-framework. More precisely, the optimal convergence rates of the solutions and their first-order derivatives in the L2-norm are obtained when the LP-norm of the perturbation is bounded for some p ε [1, 6).
文摘hi this paper. the author uses the theory of fluid mechanics. dynamics of fluids in Porous media. gas seepage flow in coal seams and combines the tests in the laboratory with the actual coal infusion to have an investigating and study from the theory to the mechanism of coal infusion to wet coal seams. through the analysis to the process of coal infusion the author builds up the mathematical models and has a detailed discussion to the boundary conditions of coal infusion. Because the equation sets to describe coal infusion are non-linear. we have made a simplification to them to use the dimension analysis theory by leading into the non-dimensions of water pressure of coal infusion, seepage flow rate. increment of coal seam moisture and so on Besides the analytic and approximate solutions have also been discussed. At last. we use the scientific research item of the actual coal infusion to illustrate the effects and importance of the theory to direct actual coal infusion and its designs.
文摘We develop a numerical solution algorithm of the nonlinear potential flow equations with the nonlinear free surface boundary condition.A finite difference method with a predictor-corrector method is applied to solve the nonlinear potential flow equations in a two-dimensional (2D) tank.The irregular tank is mapped onto a fixed square domain with rectangular cells through a proper mapping function.A staggered mesh system is adopted in a 2D tank to capture the wave elevation of the transient fluid.The finite difference method with a predictor-corrector scheme is applied to discretize the nonlinear dynamic boundary condition and nonlinear kinematic boundary condition.We present the numerical results of wave elevations from small to large amplitude waves with free oscillation motion,and the numerical solutions of wave elevation with horizontal excited motion.The beating period and the nonlinear phenomenon are very clear.The numerical solutions agree well with the analytical solutions and previously published results.
文摘Three dimensional Euler equations are solved in the finite volume form with van Leer's flux vector splitting technique. Block matrix is inverted by Gauss-Seidel iteration in two dimensional plane while strongly implicit alternating sweeping is implemented in the direction of the third dimension. Very rapid convergence rate is obtained with CFL number reaching the order of 100. The memory resources can be greatly saved too. It is verified that the reflection boundary condition can not be used with flux vector splitting since it will produce too large numerical dissipation. The computed flow fields agree well with experimental results. Only one or two grid points are there within the shock transition zone.
文摘A complete boundary integral formulation for steady compressible inviscid flows governed by nonlinear equations is established by using ρV as variable. Thus, the dimensionality of the problem to be solved is reduced by one and the computational mesh to be generated is needed only on the boundary of the domain.
文摘AbstractThis paper gives H-O type grid generation and applies finite-volume method forwing-body flow.A 3-D Euler code for wing-body has been developed. We improve in thiscode the distribution of the number of grid points over wing and body surfaces,and keepthe training edge of the wing being one of the grid lines even when the edge has sweepbackor sweep forward.The code developed can provide not only the spanwise pressure distribu-tion as usual ou cross-flow planes,but also the chordwise Pressure distributions.The com-putation for NASA TN D-712 wing-body model with the present code shows that thecomputed pressure distributions are in very good agreement with the experiment.
文摘This paper presents analytieal solutions to the partial differential equations for unsteady flow of the second-order fluid and Maxwell fluid in tube by using the integral transform method. It can be used to analyse the behaviour of axial velocity and shear stress for unsteady flow of nun-Newtonian visco-elastie fluids in tube, and to provide a theoretical base for the projection of pipe-line engineering.
文摘Inthis paper, each of the two phases in dense two-phase flow is considered as continuous medium and the fundamental equations for two-phase flow arc described in Eulerian form. The generalized constitutive relation of the Bingham fluid is applied to the dispersed phase with the analysis oj physical mechanism of dense two-phase flow. The shearing stress of dispersed phase at a wall is used to give a boundary condition. Then a mathematical model for dense two-phase flow is obtained. In addition, the expressions of shearing stress of dispersed phase at a wall is derived according to the fundamental model of the friclional collision between dispersed-plutse particles and the wall.
基金the funding support from the National Natural Science Foundation of China(51974013 and 11372033)the Open Research Foundation(NEPU-EOR-2019-003)the initiative funding from the University of Science and Technology Beijing.
文摘Fluid flow at nanoscale is closely related to many areas in nature and technology(e.g.,unconventional hydrocarbon recovery,carbon dioxide geo-storage,underground hydrocarbon storage,fuel cells,ocean desalination,and biomedicine).At nanoscale,interfacial forces dominate over bulk forces,and nonlinear effects are important,which significantly deviate from conventional theory.During the past decades,a series of experiments,theories,and simulations have been performed to investigate fluid flow at nanoscale,which has advanced our fundamental knowledge of this topic.However,a critical review is still lacking,which has seriously limited the basic understanding of this area.Therefore herein,we systematically review experimental,theoretical,and simulation works on single-and multi-phases fluid flow at nanoscale.We also clearly point out the current research gaps and future outlook.These insights will promote the significant development of nonlinear flow physics at nanoscale and will provide crucial guidance on the relevant areas.
文摘This paper gives an overview on nonlinear porous flow in low permeability porous media, reveals the microscopic mechanisms of flows, and clarifies properties of porous flow fluids. It shows that, deviating from Darcy's linear law, the porous flow characteristics obey a nonlinear law in a low-permeability porous medium, and the viscosity of the porous flow fluid and the permeability values of water and oil are not constants. Based on these characters, a new porous flow model, which can better describe low permeability reservoir~ is established. This model can describe various patterns of porous flow, as Darcy's linear law does. All the parameters involved in the model, having definite physical meanings, can be obtained directly from the experiments.
文摘A series of simple axisymmetric compression tests were carried out on the computer servo controlled Gleeble 1 500 machine when strain rates ranged between 0.05 25 s -1 and deformation temperature 300 500 ℃. The results show that flow stress is related to the Zener Hollonom parameter Z and strain, as well as the static recrystallization fraction between passes during multipass hot deformation of 5182 aluminum alloy. Hence, a modified exponential flow stress equation was presented by considering the values of ln A and β as functions of strain, and by using the uniform softening method and incorporating the static recrystallization fraction between passes to consider the effects of residual strain during multipass hot rolling of 5182 aluminum alloy. The validity of the equation was examined by a typical non isothermal multipass deformation test.
文摘Underwater gliders are recent innovative types of autonomous underwater vehicles (AUVs) used in ocean exploration and observation. They adjust their buoyancy to dive and to return to the ocean surface. During the change of altitude, they use the hydrodynamic forces developed by their wings to move forward. Their flights are controlled by changing the position of their centers of gravity and their buoyancy to adjust their trim and heel angles. For better flight control, the understanding of the hydrodynamic behavior and the flight mechanics of the underwater glider is necessary. A 6-DOF motion simulator is coupled with an unsteady potential flow model for this purpose. In some specific cases, the numerical study demonstrates that an inappropriate stabilizer dimension can cause counter-steering behavior. The simulator can be used to improve the automatic flight control. It can also be used for the hydrodynamic design optimization of the devices.
文摘In this paper, the fluid flow differential equation based on the homogenous reservoirs model is first reviewed. Then a theorem about the formal similarity of solutions in the Laplace space with outer boundary conditions and inner boundary condition is presented and proved. Lastly, a corollary of our theorem is given particularly on inner boundary. The obtained results are very helpful for understanding inherent laws of relevant engineering science and designing practical analysis software.
基金Project supported by the National Natural Science Foundation of China(Grant No.11072117)the Scientific Research Fund of Zhejiang Province,China(Grant No.LY13A010005)+4 种基金the Disciplinary Project of Ningbo City,China(Grant No.SZXL1067)the Scientific Research Fund of Education Department of Zhejiang Province,China(Grant No.Z201119278)the Natural Science Foundation of Ningbo City,China(Grant Nos.2012A610152 and 2012A610038)the K.C.Wong Magna Fund in Ningbo University,Chinathe Research Grant Council,Government of the Hong Kong Administrative Region,China(Grant No.CityU119011)
文摘This paper focuses on a two-dimensional bidirectional pedestrian flow model which involves the next-nearest-neighbor effect. The stability condition and the Korteweg-de Vries (KdV) equation are derived to describe the density wave of pedestrian congestion by linear stability and nonlinear analysis. Through theoretical analysis, the soliton solution is obtained.
文摘This paper defines new kinds of functions——the conjugate axisymmetric poteptial functions. With the aid of them, we can prove the completeness of the solutions of the generalized axisymmetric Stokes flow equation without the condition on the convexity of the domain.
基金partially supported by Grant-in-Aid for Scientific Research 26292077, 2014, from the Ministry of Education, Science, Sports, and Culture of Japanby the River Fund in charge of the River Foundation, Japan
文摘Debris flow simulations are useful for predicting the sediment supplied to watersheds from upstream areas. However, the topographic conditions upstream are more complicated than those downstream and the relationship between the topographic conditions and debris flow initiation is not well understood. This study compared the use of several entrainment rate equations in numerical simulations of debris flows to examine the effect of topographic conditions on the flow. One-dimensional numerical simulations were performed based on the shallow water equations and three entrainment rate equations were tested. These entrainment rate equations were based on the same idea that erosion and the deposition of debris flows occur via the difference between the equilibrium and current conditions of debris flows, while they differed in the expression of the concentration, channel angle, and sediment amount. The comparison was performed using a straight channel with various channel angles and a channel with a periodically undulating surface. The three entrainment rate equations gave different amounts of channel bed degradation and hydrographs for a straight channel with a channel angle greater than 21° when water was supplied from upstream at a steady rate. The difference was caused by the expression of the entrainment rate equations. For channels with little undulation, the numerical simulations gave results almost identical to those for straight channels with the same channel angle. However, for channels with large undulations, the hydrographs differed from those for straight channels with the same channel angle when the channel angle was less than 21°. Rapid erosion occurred and the hydrograph showed a significant peak, especially in cases using the entrainment equation expressed by channel angle. This was caused by the effects of the steep undulating sections, since the effect increased with the magnitude of the undulation, suggesting that a debris flow in an upstream area develops differently according to the topographic conditions. These results also inferred that numerical simulations of debris flow can differ depending on the spatial resolution of the simulation domain, as the resolution determines the reproducibility of the undulations.
文摘Through transformations, the time-dependent boundary condition on the airfoil contour and the boundary condition at infinity are brought fixed to the boundaries of a finite domain. The boundary conditions can thus be satisfied exactly without increasing the computational time. The novel scheme is useful for computing transonic, strong disturbance, unsteady flows with high reduced frequencies. The scheme makes use of curvefitted orthogonal meshes and the lattice control technique to obtain the optimal grid distribution. The numerical results are satisfactory.
基金The project is supported by the National Nature Science Foundation of China(NNSFC)
文摘The finite difference method (FDM) is applied in the present paper to solve the unsteady NHS equations for incompressible fluids. ADI and SLOR methods are served for the vorticity equation and the Poisson equation for ψ respectively. The upwind scheme is used for the convective terms. The moving boundary conditions are specially treated, and the effects of outlet conditions on the flow field are abo examined. Numerical results obtained show that the spoiler's oscillation induces forming, growing and shedding of the vortices. The shedding frequency of vortices is equal to that of the spoiler's oscillation. The forced unsteady separated flows under the present investigation depend mainly on the reduced frequency. At low reduced frequency, the vortices shed from the spoiler interact weakly with each other, and move downstream at an almost uniform speed of 038 V∞. At high reduced frequency, the interaction between the adjacent vortices strengthens. They close up to and rotate around each other, and eventually, merge into one vortex.