The surface dielectric barrier discharge (SDBD) plasma actuator has shown great promise as an aerodynamic flow control device. In this paper, the encapsulated electrode width of a SDBD actuator is changed to study t...The surface dielectric barrier discharge (SDBD) plasma actuator has shown great promise as an aerodynamic flow control device. In this paper, the encapsulated electrode width of a SDBD actuator is changed to study the airflow acceleration behavior. The effects of encapsulated electrode width on the actuator performance are experimentally investigated by measuring the dielectric layer surface potential, time-averaged ionic wind velocity and thrust force. Experimental results show that the airflow velocity and thrust force increase with the encapsulated electrode width. The results can be attributed to the distinct plasma distribution at different encapsulated electrode widths.展开更多
A numerical model for shallow water flow has been developed based on the unsteady Reynolds-averaged Navier-Stokes equations with the hydrodynamic pressure instead of hydrostatic pressure assumption. The equations are ...A numerical model for shallow water flow has been developed based on the unsteady Reynolds-averaged Navier-Stokes equations with the hydrodynamic pressure instead of hydrostatic pressure assumption. The equations are transformed into the σ-coordinate system and the eddy viscosity is calculated with the standard k-ε turbulence model. The control volume method is used to discrete the equations, and the boundary conditions at the bed for shallow water models only include vertical diffusion terms expressed with wall functions. And the semi-implicit method for pressure linked equation arithmetic is adopted to solve the equations. The model is applied to the 2D vertical plane flow of a current over two steep-sided trenches for which experiment data are available for comparison and good agreement is obtained. And the model is used to predicting the flow in a channel with a steep-sided submerged breakwater at the bottom, and the streamline is drawn.展开更多
A kinetic model was proposed to predict the seawater fouling process in the seawater heat exchangers.The new model adopted an expression combining depositional and removal behaviors for seawater fouling based on the K...A kinetic model was proposed to predict the seawater fouling process in the seawater heat exchangers.The new model adopted an expression combining depositional and removal behaviors for seawater fouling based on the Kern–Seaton model.The present model parameters include the integrated kinetic rate of deposition(k d)and the integrated kinetic rate of removal(k r),which have clear physical signi ficance.A seawater-fouling monitoring device was established to validate the model.The experimental data were well fitted to the model,and the parameters were obtained in different conditions.SEM and EDX analyses were performed after the experiments,and the results show that the main components of seawater fouling are magnesium hydroxide and aluminum hydroxide.The effects of surface temperature,flow velocity and surface free energy were assessed by the model and the experimental data.The results indicate that the seawater fouling becomes aggravated as the surface temperature increased in a certain range,and the seawater fouling resistance reduced as the flow velocity of seawater increased.Furthermore,the effect of the surface free energy of metals was analyzed,showing that the lower surface free energy mitigates the seawater fouling accumulation.展开更多
A constrained interpolation profile CIP-based numerical tank is developed to simulate violent free surface flows.The numerical simulation is performed by the CIP-based Cartesian grid method,which is described in the p...A constrained interpolation profile CIP-based numerical tank is developed to simulate violent free surface flows.The numerical simulation is performed by the CIP-based Cartesian grid method,which is described in the present paper.The tangent of hyperbola for interface capturing(THINC) scheme is applied for capturing complex free surfaces.The new model is capable of simulating a flow with violently varied free surface.A series of computations are conducted to assess the developed algorithm and its versatility.These tests include the collapse of water column with and without an obstacle,sloshing in a fixed tank,the generation of regular waves in a tank,the generation of extreme waves in a tank.Excellent agreements are obtained when numerical results are compared with available analytical,experimental,and other numerical results.展开更多
Effective control of hypersonic transition is essential.In order to avoid affecting the structural proflle of the aircraft,as well as reducing power consumption and electromagnetic interference,a low-frequency surface...Effective control of hypersonic transition is essential.In order to avoid affecting the structural proflle of the aircraft,as well as reducing power consumption and electromagnetic interference,a low-frequency surface arc plasma disturbance experiment to promote hypersonic transition was carried out in theΦ0.25 m double-throat Ludwieg tube wind tunnel at Huazhong University of Science and Technology.Contacting printed circuit board sensors and non-contact focused laser differential interferometry testing technology were used in combination.Experimental results showed that the low-frequency surface arc plasma actuation had obvious stimulation effects on the second-mode unstable wave and could promote boundary layer transition by changing the spectral characteristics of the second-mode unstable wave.At the same time,the plasma actuation could promote energy exchange between the second-mode unstable wave and other unstable waves.Finally,the corresponding control mechanism is discussed.展开更多
The particle simulation method is used to solve free-surface slurry flow problems that may be encountered in several scientific and engineering fields.The main idea behind the use of the particle simulation method is ...The particle simulation method is used to solve free-surface slurry flow problems that may be encountered in several scientific and engineering fields.The main idea behind the use of the particle simulation method is to treat granular or other materials as an assembly of many particles.Compared with the continuum-mechanics-based numerical methods such as the finite element and finite volume methods,the movement of each particle is accurately described in the particle simulation method so that the free surface of a slurry flow problem can be automatically obtained.The major advantage of using the particle simulation method is that only a simple numerical algorithm is needed to solve the governing equation of a particle simulation system.For the purpose of illustrating how to use the particle simulation method to solve free-surface flow problems,three examples involving slurry flow on three different types of river beds have been considered.The related particle simulation results obtained from these three examples have demonstrated that:1) The particle simulation method is a promising and useful method for solving free-surface flow problems encountered in both the scientific and engineering fields;2) The shape and irregular roughness of a river bed can have a significant effect on the free surface morphologies of slurry flow when it passes through the river bed.展开更多
In order to accurately simulate strong three-dimensional (3-D) free surface flows and sediment transport, the fully 3- D non-hydrostatic pressure models are developed based on the incompressible Navier-Stokes equati...In order to accurately simulate strong three-dimensional (3-D) free surface flows and sediment transport, the fully 3- D non-hydrostatic pressure models are developed based on the incompressible Navier-Stokes equations and convection-diffusion equation of sediment concentration with the mixing triangle and quadrilateral grids. The governing equations are discretized with the unstructured finite volume method in order to provide conservation properties of mass and momentum, and flexibility with practical application. It is shown that it is first-order accurate on nonuniform plane two-dimensional (2-D) grids and second-order accurate on uniform plane grids. A third-order approximation of the vertical velocity at the top-layer is applied. In such a way, free surface zero stress boundary condition is satisfied maturely, and very few vertical layers are needed to give an accurate solution even for complex discontinuous flow and short wave simulation. The model is applied to four examples to simulate strong 3-D free surface flows and sediment transport where non-hydrostatic pressures have a considerable effect on the velocity field. The newly developed model is verified against analytical solutions with an excellent agreement.展开更多
This paper investigates the behaviour of traffic flow in traffic systems with a new model based on the NaSch model and cluster approximation of mean-field theory. The proposed model aims at constructing a mapping rela...This paper investigates the behaviour of traffic flow in traffic systems with a new model based on the NaSch model and cluster approximation of mean-field theory. The proposed model aims at constructing a mapping relationship between the microcosmic behaviour and the macroscopic property of traffic flow. Results demonstrate that scale-free phenomenon of the evolution network becomes obvious when the density value of traffic flow reaches at the critical point of phase transition from free flow to traffic congestion, and jamming is limited in this scale-free structure.展开更多
Based on the complex dispersion relation for the two-dimensional free-surface waves generated by a moving body in the steady Oseen flows, the effect of viscosity on wavelength and wave amplitude was investigated by me...Based on the complex dispersion relation for the two-dimensional free-surface waves generated by a moving body in the steady Oseen flows, the effect of viscosity on wavelength and wave amplitude was investigated by means of an asymptotic method and a numerical analysis. A comparison between the asymptotic and numerical analysis for the viscous decay factor demonstrates the validity of the perturbation expansions for the wave profile. The numerical result shows that the wavelength of viscous wave is slightly elongated in comparison with that of inviscid wave.展开更多
The motion of incompressible fluid of a variable fluid viscosity and variable thermal conductivity with thermal radiation, Dufour, Soret with heat and mass transfer over a linearly moving porous vertical semi-infinite...The motion of incompressible fluid of a variable fluid viscosity and variable thermal conductivity with thermal radiation, Dufour, Soret with heat and mass transfer over a linearly moving porous vertical semi-infinite plate with suction is investigated. The governing equations are transformed into a system of coupled nonlinear ordinary differential equations using similarity transformations with dimensionless variables and solved numerically using shooting method with Runge-Kutta fourth-order method and Newton-Raphson’s interpolation scheme implemented in MATLAB. The result showed that with increase in Dufour and Soret parameter, fluid velocity increases and temperature increases with increase in variation of Dufour while, temperature decreases with increase in Soret. The effects of variable fluid viscosity, variable thermal conductivity, thermal radiation, Soret, Dufour, Prandtl and Schmidt parameters on the dimensionless velocity, temperature and concentration profiles are shown graphically.展开更多
Simulating turbulent liquids with breaking waves and splashes is among the most desired features in fluid animation. Lagrangian methods such as Smoothed Particle Hydrodynamics method (SPH) are a promising way to captu...Simulating turbulent liquids with breaking waves and splashes is among the most desired features in fluid animation. Lagrangian methods such as Smoothed Particle Hydrodynamics method (SPH) are a promising way to capture such properties. However, the Particle-based liquid surface simulation has not been applied very well since its consumption is way too large. This paper derives the governing equations in SPH approaches and parallelizes the dynamics-based surface simulation with the MapReduce program models which apply the SPH approach in Cloud Computing. Compared to the serial methods, this approach obtained a 3.11 times speedup on the experimental platform.展开更多
A higher-efficient three-dimensional non-hydrostatic model is developed to simulate small amplitude free surface flows based on a staggered unstructured grid. In this model, a fractional step algorithm is adopted to s...A higher-efficient three-dimensional non-hydrostatic model is developed to simulate small amplitude free surface flows based on a staggered unstructured grid. In this model, a fractional step algorithm is adopted to solve the Navier-Stokes equations in two major steps. A top-layer pressure method is proposed to minimize the number of vertical layers and subsequently the computational cost. Three classical examples of small amplitude free surface flows are used to demonstrate the capability and efficiency of the model. The satisfactory results demonstrated the capability and efficiency of modelling a range of small amplitude free surface flows with only a small number of vertical layers.展开更多
A staggered finite-volume technique for non-hydrostatic, small amplitude free surface flow governed by the incompressible Navier-Stokes equations is presented there is a proper balance between accuracy and computing t...A staggered finite-volume technique for non-hydrostatic, small amplitude free surface flow governed by the incompressible Navier-Stokes equations is presented there is a proper balance between accuracy and computing time. The advection and horizontal diffusion terms in the momentum equation are discretized by an integral interpolation method on the orthogonal unstructured staggered mesh and, while it has the attractive property of being conservative. The pressure-correction algorithm is employed for the non-hydrostatic pressure in order to achieve second-order temporal accuracy. A conservative scalar transport algorithm is also applied to discretize k - c equations in this model. The eddy viscosity is calculated from the k-c turbulent model. The resulting model is mass and momentum conservative. The model is verified by two examples to simulate unsteady small amplitude free surface flows where non-hydrostatic pressures have a considerable effect on the velocity field, and then applied to simulate the tidal flow in the Bohai Sea.展开更多
The simplified modeling for analysis on MHD stability of free surface jet flow in a gradient magnetic fields is based on the theoretical and experimental results on channel liquid metal MHD flow, especially, the resul...The simplified modeling for analysis on MHD stability of free surface jet flow in a gradient magnetic fields is based on the theoretical and experimental results on channel liquid metal MHD flow, especially, the results of MHD flow velocity distribution in cross-section of channels (rectangular duct and circular pipe), and the expected results from the modeling are well agreed with the recent experimental data obtained. It is the first modeling which can efficiently explain the experimental results of liquid-metal free surface jet flow.展开更多
The free surface flow generated by twin-cylinders in forced motion submerged beneath the free surface is studied based on the boundary element method. Two relative locations, namely, horizontal and vertical, are exami...The free surface flow generated by twin-cylinders in forced motion submerged beneath the free surface is studied based on the boundary element method. Two relative locations, namely, horizontal and vertical, are examined for the twin cylinders. In both cases, the twin cylinders are starting from rest and ultimately moving with the same constant speed through an accelerating process. Assuming that the fluid is inviscid and incompressible and the flow to be irrotational, the integral Laplace equation can be discretized based on the boundary element method. Fully-nonlinear boundary conditions are satisfied on the unknown free surface and the moving body surface. The free surface is traced by a Lagrangian technique. Regriding and remeshing are applied, which is crucial to quality of the numerical results. Single circular cylinder and elliptical cylinder are calculated by linear method and fully nonlinear method for accuracy checking and then fully nonlinear method is conducted on the twin cylinder cases, respectively. The generated wave elevation and the resultant force are analysed to discuss the influence of the gap between the two cylinders as well as the water depth. It is found that no matter the kind of distribution, when the moving cylinders are close to each other, they suffer hydrodynamic force with large absolute value in the direction of motion. The trend of force varying with the increase of gap can be clearly seen from numerical analysis. The vertically distributed twin cylinders seem to attract with each other while the horizontally distributed twin cylinders are opposite when they are close to each other.展开更多
In this paper,we present a conservative semi-Lagrangian scheme designed for the numeri-cal solution of 3D hydrostatic free surface flows involving sediment transport on unstruc-tured Voronoi meshes.A high-order recons...In this paper,we present a conservative semi-Lagrangian scheme designed for the numeri-cal solution of 3D hydrostatic free surface flows involving sediment transport on unstruc-tured Voronoi meshes.A high-order reconstruction procedure is employed for obtaining a piecewise polynomial representation of the velocity field and sediment concentration within each control volume.This is subsequently exploited for the numerical integration of the Lagrangian trajectories needed for the discretization of the nonlinear convective and viscous terms.The presented method is fully conservative by construction,since the transported quantity or the vector field is integrated for each cell over the deformed vol-ume obtained at the foot of the characteristics that arises from all the vertexes defining the computational element.The semi-Lagrangian approach allows the numerical scheme to be unconditionally stable for what concerns the advection part of the governing equations.Furthermore,a semi-implicit discretization permits to relax the time step restriction due to the acoustic impedance,hence yielding a stability condition which depends only on the explicit discretization of the viscous terms.A decoupled approach is then employed for the hydrostatic fluid solver and the transport of suspended sediment,which is assumed to be passive.The accuracy and the robustness of the resulting conservative semi-Lagrangian scheme are assessed through a suite of test cases and compared against the analytical solu-tion whenever is known.The new numerical scheme can reach up to fourth order of accu-racy on general orthogonal meshes composed by Voronoi polygons.展开更多
We have numerically and experimentally investigated the flow modes of Taylor-Couette system consisting of coaxial two cylinders with vertical axes. The inner cylinder rotates and the outer cylinder and the bottom end ...We have numerically and experimentally investigated the flow modes of Taylor-Couette system consisting of coaxial two cylinders with vertical axes. The inner cylinder rotates and the outer cylinder and the bottom end of the cylinders remain stationary. The upper top boundary is the free surface of the working liquid between the inner and outer cylinders and it contacts with the air. While this flow appears in fluid machinery and chemical reactors and includes industrial interests, it also contains problems of fluid mechanics, which is about the behavior of the free surface in the rotating field. In this paper, we concretely show the developments of the one cell mode flow and the three cell mode flow at a small aspect ratio. We also represent the bifurcation diagram of the flow at the moderate aspect ratio about 5.5. In the numerical simulation, the flow is rest in the initial state, and the inner cylinder is linearly or suddenly accelerated to attain a flow with a prescribed Reynolds number. When the acceleration of the inner cylinder is high, an imperfect bifurcation occurs and the flows of the secondary modes emerge. At high Reynolds numbers, the flow first has many vortices and then some of the vortices collapse and the final stable flow arises. The loci of the normal five cell mode, the anomalous six cell mode and the secondary seven cell mode are determined.展开更多
A doublet integral equation is formulated for the two-dimensional dissipative potential flow around a hydrofoil submerged below a free-water surface. The free-water surface is assumed to involve energy dissipation, an...A doublet integral equation is formulated for the two-dimensional dissipative potential flow around a hydrofoil submerged below a free-water surface. The free-water surface is assumed to involve energy dissipation, and thus it is the source of damping. A doublet panel method is developed from incorporation of the dissipative Green function approach and the doublet distributions on the hydrofoil surface. Numerical computations are implemented, and the derived numerical results are in good agreement with analytic solutions and experimental measurements.展开更多
Supersonic flows past two-dimensional cavities with/without control are investigated by the direct numerical simulation (DNS). For an uncontrolled cavity, as the thickness of the boundary layer declines, transition ...Supersonic flows past two-dimensional cavities with/without control are investigated by the direct numerical simulation (DNS). For an uncontrolled cavity, as the thickness of the boundary layer declines, transition of the dominant mode from the steady mode to the Rossiter Ⅱ mode and then to the Rossiter III mode is observed due to the change of vortex-corner interactions. Meanwhile, a low frequency mode appears. However, the wake mode observed in a subsonic cavity flow is absent in the current simulation. The oscillation frequencies obtained from a global dynamic mode decomposition (DMD) approach are consistent with the local power spectral density (PSD) analysis. The dominant mode transition is clearly shown by the dynamic modes obtained from the DMD. A passive control technique of substituting the cavity trailing edge with a quarter-circle is studied. As the effective cavity length increases, the dominant mode transition from the Rossiter Ⅱ mode to the Rossiter Ⅲ mode occurs. With the control, the pressure oscillations are reduced significantly. The interaction of the shear layer and the recirculation zone is greatly weakened, combined with weaker shear layer instability, responsible for the suppression of pressure oscillations. Moreover, active control using steady subsonic mass injection upstream of a cavity leading edge can stabilize the flow.展开更多
基金supported by National Natural Science Foundation of China(No.11175037)National Natural Science Foundation for Young Scientists of China(No.11305017)Special Fund for Theoretical Physics(No.11247239)
文摘The surface dielectric barrier discharge (SDBD) plasma actuator has shown great promise as an aerodynamic flow control device. In this paper, the encapsulated electrode width of a SDBD actuator is changed to study the airflow acceleration behavior. The effects of encapsulated electrode width on the actuator performance are experimentally investigated by measuring the dielectric layer surface potential, time-averaged ionic wind velocity and thrust force. Experimental results show that the airflow velocity and thrust force increase with the encapsulated electrode width. The results can be attributed to the distinct plasma distribution at different encapsulated electrode widths.
文摘A numerical model for shallow water flow has been developed based on the unsteady Reynolds-averaged Navier-Stokes equations with the hydrodynamic pressure instead of hydrostatic pressure assumption. The equations are transformed into the σ-coordinate system and the eddy viscosity is calculated with the standard k-ε turbulence model. The control volume method is used to discrete the equations, and the boundary conditions at the bed for shallow water models only include vertical diffusion terms expressed with wall functions. And the semi-implicit method for pressure linked equation arithmetic is adopted to solve the equations. The model is applied to the 2D vertical plane flow of a current over two steep-sided trenches for which experiment data are available for comparison and good agreement is obtained. And the model is used to predicting the flow in a channel with a steep-sided submerged breakwater at the bottom, and the streamline is drawn.
基金Supported by the Leading Academic Discipline Project of Shanghai Municipal Education Commission(J50502)the Construction of Shanghai Science and Technology Commission(13DZ2260900)
文摘A kinetic model was proposed to predict the seawater fouling process in the seawater heat exchangers.The new model adopted an expression combining depositional and removal behaviors for seawater fouling based on the Kern–Seaton model.The present model parameters include the integrated kinetic rate of deposition(k d)and the integrated kinetic rate of removal(k r),which have clear physical signi ficance.A seawater-fouling monitoring device was established to validate the model.The experimental data were well fitted to the model,and the parameters were obtained in different conditions.SEM and EDX analyses were performed after the experiments,and the results show that the main components of seawater fouling are magnesium hydroxide and aluminum hydroxide.The effects of surface temperature,flow velocity and surface free energy were assessed by the model and the experimental data.The results indicate that the seawater fouling becomes aggravated as the surface temperature increased in a certain range,and the seawater fouling resistance reduced as the flow velocity of seawater increased.Furthermore,the effect of the surface free energy of metals was analyzed,showing that the lower surface free energy mitigates the seawater fouling accumulation.
基金supported by the Fundamental Research Funds for the Central Universities
文摘A constrained interpolation profile CIP-based numerical tank is developed to simulate violent free surface flows.The numerical simulation is performed by the CIP-based Cartesian grid method,which is described in the present paper.The tangent of hyperbola for interface capturing(THINC) scheme is applied for capturing complex free surfaces.The new model is capable of simulating a flow with violently varied free surface.A series of computations are conducted to assess the developed algorithm and its versatility.These tests include the collapse of water column with and without an obstacle,sloshing in a fixed tank,the generation of regular waves in a tank,the generation of extreme waves in a tank.Excellent agreements are obtained when numerical results are compared with available analytical,experimental,and other numerical results.
基金supported by National Science and Technology Major Project(No.J2019-II-0014-0035)。
文摘Effective control of hypersonic transition is essential.In order to avoid affecting the structural proflle of the aircraft,as well as reducing power consumption and electromagnetic interference,a low-frequency surface arc plasma disturbance experiment to promote hypersonic transition was carried out in theΦ0.25 m double-throat Ludwieg tube wind tunnel at Huazhong University of Science and Technology.Contacting printed circuit board sensors and non-contact focused laser differential interferometry testing technology were used in combination.Experimental results showed that the low-frequency surface arc plasma actuation had obvious stimulation effects on the second-mode unstable wave and could promote boundary layer transition by changing the spectral characteristics of the second-mode unstable wave.At the same time,the plasma actuation could promote energy exchange between the second-mode unstable wave and other unstable waves.Finally,the corresponding control mechanism is discussed.
基金Project(11272359)supported by the National Natural Science Foundation of China
文摘The particle simulation method is used to solve free-surface slurry flow problems that may be encountered in several scientific and engineering fields.The main idea behind the use of the particle simulation method is to treat granular or other materials as an assembly of many particles.Compared with the continuum-mechanics-based numerical methods such as the finite element and finite volume methods,the movement of each particle is accurately described in the particle simulation method so that the free surface of a slurry flow problem can be automatically obtained.The major advantage of using the particle simulation method is that only a simple numerical algorithm is needed to solve the governing equation of a particle simulation system.For the purpose of illustrating how to use the particle simulation method to solve free-surface flow problems,three examples involving slurry flow on three different types of river beds have been considered.The related particle simulation results obtained from these three examples have demonstrated that:1) The particle simulation method is a promising and useful method for solving free-surface flow problems encountered in both the scientific and engineering fields;2) The shape and irregular roughness of a river bed can have a significant effect on the free surface morphologies of slurry flow when it passes through the river bed.
基金financially supported by the Science and Technology Project of the Ministry of Transport(Grant No.2013328352570)
文摘In order to accurately simulate strong three-dimensional (3-D) free surface flows and sediment transport, the fully 3- D non-hydrostatic pressure models are developed based on the incompressible Navier-Stokes equations and convection-diffusion equation of sediment concentration with the mixing triangle and quadrilateral grids. The governing equations are discretized with the unstructured finite volume method in order to provide conservation properties of mass and momentum, and flexibility with practical application. It is shown that it is first-order accurate on nonuniform plane two-dimensional (2-D) grids and second-order accurate on uniform plane grids. A third-order approximation of the vertical velocity at the top-layer is applied. In such a way, free surface zero stress boundary condition is satisfied maturely, and very few vertical layers are needed to give an accurate solution even for complex discontinuous flow and short wave simulation. The model is applied to four examples to simulate strong 3-D free surface flows and sediment transport where non-hydrostatic pressures have a considerable effect on the velocity field. The newly developed model is verified against analytical solutions with an excellent agreement.
基金supported by the National Basic Research Program of China (973) (Grant No 2006CB705500)the National Natural Science Foundation of China (Grant No 70671008)
文摘This paper investigates the behaviour of traffic flow in traffic systems with a new model based on the NaSch model and cluster approximation of mean-field theory. The proposed model aims at constructing a mapping relationship between the microcosmic behaviour and the macroscopic property of traffic flow. Results demonstrate that scale-free phenomenon of the evolution network becomes obvious when the density value of traffic flow reaches at the critical point of phase transition from free flow to traffic congestion, and jamming is limited in this scale-free structure.
文摘Based on the complex dispersion relation for the two-dimensional free-surface waves generated by a moving body in the steady Oseen flows, the effect of viscosity on wavelength and wave amplitude was investigated by means of an asymptotic method and a numerical analysis. A comparison between the asymptotic and numerical analysis for the viscous decay factor demonstrates the validity of the perturbation expansions for the wave profile. The numerical result shows that the wavelength of viscous wave is slightly elongated in comparison with that of inviscid wave.
文摘The motion of incompressible fluid of a variable fluid viscosity and variable thermal conductivity with thermal radiation, Dufour, Soret with heat and mass transfer over a linearly moving porous vertical semi-infinite plate with suction is investigated. The governing equations are transformed into a system of coupled nonlinear ordinary differential equations using similarity transformations with dimensionless variables and solved numerically using shooting method with Runge-Kutta fourth-order method and Newton-Raphson’s interpolation scheme implemented in MATLAB. The result showed that with increase in Dufour and Soret parameter, fluid velocity increases and temperature increases with increase in variation of Dufour while, temperature decreases with increase in Soret. The effects of variable fluid viscosity, variable thermal conductivity, thermal radiation, Soret, Dufour, Prandtl and Schmidt parameters on the dimensionless velocity, temperature and concentration profiles are shown graphically.
基金supported by National High Technical Research and Development Program of China (863 program) under GrantNo.2009AA062801National Natural Science Foundation of China under Grant No.60973063+2 种基金Beijing Natural Science Foundation of China under Grant No.4092028China Fundamental Research Funds for the Central Universities under Grant No.FRF-TP-09-016BNew Century Personnel Plan for the Ministry of Education of China under Grant No.NCET-10-0221
文摘Simulating turbulent liquids with breaking waves and splashes is among the most desired features in fluid animation. Lagrangian methods such as Smoothed Particle Hydrodynamics method (SPH) are a promising way to capture such properties. However, the Particle-based liquid surface simulation has not been applied very well since its consumption is way too large. This paper derives the governing equations in SPH approaches and parallelizes the dynamics-based surface simulation with the MapReduce program models which apply the SPH approach in Cloud Computing. Compared to the serial methods, this approach obtained a 3.11 times speedup on the experimental platform.
基金financially supported by the Science and Technology Project of the Ministry of Transport(Grant No.2013328352570)partly supported by the National Natural Science Foundation of China(Grant No.51209112)
文摘A higher-efficient three-dimensional non-hydrostatic model is developed to simulate small amplitude free surface flows based on a staggered unstructured grid. In this model, a fractional step algorithm is adopted to solve the Navier-Stokes equations in two major steps. A top-layer pressure method is proposed to minimize the number of vertical layers and subsequently the computational cost. Three classical examples of small amplitude free surface flows are used to demonstrate the capability and efficiency of the model. The satisfactory results demonstrated the capability and efficiency of modelling a range of small amplitude free surface flows with only a small number of vertical layers.
基金financially supported by the Science and Technology Project of the Ministry of Transport (Grant No. 2011329224170)
文摘A staggered finite-volume technique for non-hydrostatic, small amplitude free surface flow governed by the incompressible Navier-Stokes equations is presented there is a proper balance between accuracy and computing time. The advection and horizontal diffusion terms in the momentum equation are discretized by an integral interpolation method on the orthogonal unstructured staggered mesh and, while it has the attractive property of being conservative. The pressure-correction algorithm is employed for the non-hydrostatic pressure in order to achieve second-order temporal accuracy. A conservative scalar transport algorithm is also applied to discretize k - c equations in this model. The eddy viscosity is calculated from the k-c turbulent model. The resulting model is mass and momentum conservative. The model is verified by two examples to simulate unsteady small amplitude free surface flows where non-hydrostatic pressures have a considerable effect on the velocity field, and then applied to simulate the tidal flow in the Bohai Sea.
基金The project partly supported by National Nature Science Foundation of China (No. 10275019)
文摘The simplified modeling for analysis on MHD stability of free surface jet flow in a gradient magnetic fields is based on the theoretical and experimental results on channel liquid metal MHD flow, especially, the results of MHD flow velocity distribution in cross-section of channels (rectangular duct and circular pipe), and the expected results from the modeling are well agreed with the recent experimental data obtained. It is the first modeling which can efficiently explain the experimental results of liquid-metal free surface jet flow.
基金Foundation item: Supported by the Lloyd's Register Foundation, the Fundamental Research Funds for the Central Universities (Gram No. HEUCF140115), the National Natural Science Foundation of China (11102048, 11302057), the Research Funds for State Key Laboratory of Ocean Engineering in Shanghai Jiao Tong University (Grant No. 1310), and the Specialized Research Fund for the Doctoral Program of Higher Education of China (20132304120028).
文摘The free surface flow generated by twin-cylinders in forced motion submerged beneath the free surface is studied based on the boundary element method. Two relative locations, namely, horizontal and vertical, are examined for the twin cylinders. In both cases, the twin cylinders are starting from rest and ultimately moving with the same constant speed through an accelerating process. Assuming that the fluid is inviscid and incompressible and the flow to be irrotational, the integral Laplace equation can be discretized based on the boundary element method. Fully-nonlinear boundary conditions are satisfied on the unknown free surface and the moving body surface. The free surface is traced by a Lagrangian technique. Regriding and remeshing are applied, which is crucial to quality of the numerical results. Single circular cylinder and elliptical cylinder are calculated by linear method and fully nonlinear method for accuracy checking and then fully nonlinear method is conducted on the twin cylinder cases, respectively. The generated wave elevation and the resultant force are analysed to discuss the influence of the gap between the two cylinders as well as the water depth. It is found that no matter the kind of distribution, when the moving cylinders are close to each other, they suffer hydrodynamic force with large absolute value in the direction of motion. The trend of force varying with the increase of gap can be clearly seen from numerical analysis. The vertically distributed twin cylinders seem to attract with each other while the horizontally distributed twin cylinders are opposite when they are close to each other.
基金support of MIUR-PRIN Project 2017,No.2017KKJP4X“Innovative numerical methods for evolutionary partial differential equations and applications”.
文摘In this paper,we present a conservative semi-Lagrangian scheme designed for the numeri-cal solution of 3D hydrostatic free surface flows involving sediment transport on unstruc-tured Voronoi meshes.A high-order reconstruction procedure is employed for obtaining a piecewise polynomial representation of the velocity field and sediment concentration within each control volume.This is subsequently exploited for the numerical integration of the Lagrangian trajectories needed for the discretization of the nonlinear convective and viscous terms.The presented method is fully conservative by construction,since the transported quantity or the vector field is integrated for each cell over the deformed vol-ume obtained at the foot of the characteristics that arises from all the vertexes defining the computational element.The semi-Lagrangian approach allows the numerical scheme to be unconditionally stable for what concerns the advection part of the governing equations.Furthermore,a semi-implicit discretization permits to relax the time step restriction due to the acoustic impedance,hence yielding a stability condition which depends only on the explicit discretization of the viscous terms.A decoupled approach is then employed for the hydrostatic fluid solver and the transport of suspended sediment,which is assumed to be passive.The accuracy and the robustness of the resulting conservative semi-Lagrangian scheme are assessed through a suite of test cases and compared against the analytical solu-tion whenever is known.The new numerical scheme can reach up to fourth order of accu-racy on general orthogonal meshes composed by Voronoi polygons.
文摘We have numerically and experimentally investigated the flow modes of Taylor-Couette system consisting of coaxial two cylinders with vertical axes. The inner cylinder rotates and the outer cylinder and the bottom end of the cylinders remain stationary. The upper top boundary is the free surface of the working liquid between the inner and outer cylinders and it contacts with the air. While this flow appears in fluid machinery and chemical reactors and includes industrial interests, it also contains problems of fluid mechanics, which is about the behavior of the free surface in the rotating field. In this paper, we concretely show the developments of the one cell mode flow and the three cell mode flow at a small aspect ratio. We also represent the bifurcation diagram of the flow at the moderate aspect ratio about 5.5. In the numerical simulation, the flow is rest in the initial state, and the inner cylinder is linearly or suddenly accelerated to attain a flow with a prescribed Reynolds number. When the acceleration of the inner cylinder is high, an imperfect bifurcation occurs and the flows of the secondary modes emerge. At high Reynolds numbers, the flow first has many vortices and then some of the vortices collapse and the final stable flow arises. The loci of the normal five cell mode, the anomalous six cell mode and the secondary seven cell mode are determined.
文摘A doublet integral equation is formulated for the two-dimensional dissipative potential flow around a hydrofoil submerged below a free-water surface. The free-water surface is assumed to involve energy dissipation, and thus it is the source of damping. A doublet panel method is developed from incorporation of the dissipative Green function approach and the doublet distributions on the hydrofoil surface. Numerical computations are implemented, and the derived numerical results are in good agreement with analytic solutions and experimental measurements.
基金supported by the National Natural Science Foundation of China(Nos.11232011 and11402262)the 111 Project of China(No.B07033)+1 种基金the China Postdoctoral Science Foundation(No.2014M561833)the Fundamental Research Funds for the Central Universities
文摘Supersonic flows past two-dimensional cavities with/without control are investigated by the direct numerical simulation (DNS). For an uncontrolled cavity, as the thickness of the boundary layer declines, transition of the dominant mode from the steady mode to the Rossiter Ⅱ mode and then to the Rossiter III mode is observed due to the change of vortex-corner interactions. Meanwhile, a low frequency mode appears. However, the wake mode observed in a subsonic cavity flow is absent in the current simulation. The oscillation frequencies obtained from a global dynamic mode decomposition (DMD) approach are consistent with the local power spectral density (PSD) analysis. The dominant mode transition is clearly shown by the dynamic modes obtained from the DMD. A passive control technique of substituting the cavity trailing edge with a quarter-circle is studied. As the effective cavity length increases, the dominant mode transition from the Rossiter Ⅱ mode to the Rossiter Ⅲ mode occurs. With the control, the pressure oscillations are reduced significantly. The interaction of the shear layer and the recirculation zone is greatly weakened, combined with weaker shear layer instability, responsible for the suppression of pressure oscillations. Moreover, active control using steady subsonic mass injection upstream of a cavity leading edge can stabilize the flow.