A series of large pilot constructed wetland （CW） systems were constructed near the confluence of an urban stream to a larger fiver in Xi＇an, a northwestern megacity in China, for treating polluted stream water befo...A series of large pilot constructed wetland （CW） systems were constructed near the confluence of an urban stream to a larger fiver in Xi＇an, a northwestern megacity in China, for treating polluted stream water before it entered the receiving water body. Each CW system is a combination of surface- and subsurface-flow cells with local gravel, sand or slag as substrates and Phragmites australis and Typha orientalis as plants. During a one-year operation with an average surface loading of 0.053 m3/（m2.day）, the overall COD, BOD, NH3-N, total nitrogen （TN） and total phosphorus （TP） removals were 72.7% ~ 4.5%, 93.4% ＋ 2.1%, 54.0% ＋ 6.3%, 53.9% ~ 6.0% and 69.4% ：t： 4.6%, respectively, which brought about an effective improvement of the fiver water quality. Surface-flow cells showed better NH3-N removal than their TN removal while subsurface-flow cells showed better TN removal than their NH3-N removal. Using local slag as the substrate, the organic and phosphorus removal could be much improved. Seasonal variation was also found in the removal of all the pollutants and autumn seemed to be the best season for pollutant removal due to the moderate water temperature and well grown plants in the CWs.展开更多
Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved met...Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded interface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.展开更多
Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved met...Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded in- terface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.展开更多
The Singular Integral Operators Method (S.I.O.M.) is applied to the determination of the free-surface profile of an un-steady flow over a spillway, which defines a classical hydraulics problem in open channel flow. Th...The Singular Integral Operators Method (S.I.O.M.) is applied to the determination of the free-surface profile of an un-steady flow over a spillway, which defines a classical hydraulics problem in open channel flow. Thus, with a known flow rate Q, then the velocities and the elevations are computed on the free surface of the spillway flow. For the numerical evaluation of the singular integral equations both constant and linear elements are used. An application is finally given to the determination of the free-surface profile of a special spillway and comparing the numerical results with corresponding results by the Boundary Integral Equation Method (B.I.E.M.) and by using experiments.展开更多
Due to global warming,changes in the rainfall intensity profile(i.e.,the temporal intensity distribution within a rainfall event)increase the difficulty of accurate erosion prediction and control.Surface cover has bee...Due to global warming,changes in the rainfall intensity profile(i.e.,the temporal intensity distribution within a rainfall event)increase the difficulty of accurate erosion prediction and control.Surface cover has been widely used as a critical measure to control soil erosion worldwide.However,the effects of the rainfall intensity profile(RIP)on soil erosion under different surface covers are not fully understood.In this study,long-term in situ field observations of the rain hyetograph,surface runoff coefficient(SRC),subsurface flow rate(SFR),and soil loss rate(SLR)from bare land,litter cover and grass cover were conducted over 11 consecutive years in the red soil hilly region of southern China.According to the occurrence time of the most intense rainfall,226 erosive events were classified into four RIP patterns:advanced,intermediate,delayed,and uniform patterns.The results indicated that the advanced pattern with short durationehigh intensity and the delayed pattern with long durationehigh depth contributed to 73.45% of the total erosive events.For bare land,advanced events were the dominant pattern producing surface runoff and soil erosion,accounting for 57.24%and 75.17%,respectively,of the total surface runoff and erosion.The average SRC and SLR from the advanced pattern were 1.29-2.42 times and 2.52-39.78 times greater than those from the other patterns,respectively.The delayed pattern contributed to subsurface flow,and the average SFR was 1.27-2.17 times greater than that of the other patterns.Furthermore,surface cover significantly reduced surface runoff and erosion and increased subsurface flow,especially under the advanced pattern.Both surface cover measures were equally effective in controlling surface runoff and erosion,but the increase in subsurface flow caused by litter cover was 1.38 -2.67 times greater than that caused by grass cover.Advanced pattern events increase the erosion risk on red soil slopes,and surface cover effectively weakens the effect of variation in the RIP pattern on soil erosion.Moreover,surface cover significantly alters the surface‒subsurface flow distribution pattern for all the RIP patterns.This study highlights the crucial importance of rain intensity profiles on water erosion and provides a basis for optimizing measures to effectively control soil and water loss under climate change.展开更多
This paper presents a physically-based integrated hydrologic model that can simulate the rain-fall-induced 2D surface water flow, 3D variably saturated subsurface flow, upland soil erosion and transport, and contamina...This paper presents a physically-based integrated hydrologic model that can simulate the rain-fall-induced 2D surface water flow, 3D variably saturated subsurface flow, upland soil erosion and transport, and contaminant transport in the surface-subsurface system of a watershed. The model couples surface and subsurface flows based on the assumption of continuity conditions of pressure head and exchange flux at the ground, considering infiltration and evapotranspiration. The upland rill/interrill soil erosion and transport are simulated using a non-equilibrium transport model. Contaminant transport in the integrated surface and subsurface domains is simulated using advection-diffusion equations with mass changes due to sediment sorption and desorption and exchanges between two domains due to infiltration, diffusion, and bed change. The model requires no special treatments at the interface of upland areas and streams and is suitable for wetland areas and agricultural watersheds with shallow streams.展开更多
The flow is assumed to be potential, and a boundary integral method is used to solve the Laplace equation for the velocity potential to investigate the shape and the position of the bubble. A 3D code to study the bubb...The flow is assumed to be potential, and a boundary integral method is used to solve the Laplace equation for the velocity potential to investigate the shape and the position of the bubble. A 3D code to study the bubble dynamics is developed, and the calculation results agree well with the experimental data. Numerical analyses are carried out for the interaction between multiple bubbles near the free surface including in-phase and out-of-phase bubbles. The calculation result shows that the bubble period increases with the decrease of the distance between bubble centres because of the depression effect between multiple bubbles. The depression has no relationship with the free surface and it is more apparent for out-of-phase bubbles. There are great differences in dynamic behaviour between the in-phase bubbles and the out-of-phase bubbles due to the depression effect. Furthermore, the interaction among eight bubbles is simulated with a three-dlmensional model, and the evolving process and the relevant physical phenomena are presented. These phenomena can give a reference to the future work on the power of bubbles induced by multiple charges exploding simultaneously or continuously.展开更多
This research is concerned with new developments and practical applications of a physically-based numerical model that incorporates new approaches for a finite elements solution to the steady/transient problems of the...This research is concerned with new developments and practical applications of a physically-based numerical model that incorporates new approaches for a finite elements solution to the steady/transient problems of the joint ground/surface water flows. Python scripts are implemented in Geographic Information System (GIS) to store, represent and take decisions on the simulated conditions related to the water resources management at the scale of the watershed. The proposed surface-subsurface model considers surface and groundwater interactions to be 2-D horizontally distributed and depth-averaged through a diffusive wave approach for surface flood routing. Infiltration rates, overland flows and evapotranspiration processes are considered by a diffuse discharge from surface water, non-saturated subsoil and groundwater table. Recent developments also allow for the management of surface water flow control through the capacity of diversion on river beds, spillways and outflow operations of floodgates in weirs and dams of reservoirs. Practical application regards the actual hydrology of the Mero River watershed, with two important water bodies mainly concerned with the water resources management at the Cecebre Reservoir and the present flooding of a deep coal mining excavation. The MELEF model (Modèle d’éLéments Fluides, in French) was adapted and calibrated during a period of five years (2008/ 2012) with the help of hydrological parameters, registered flow rates, water levels and registered precipitation, water uses and water management operations in surface and groundwater bodies. The results predict the likely evolution of the Cecebre Reservoir, the flow rates in rivers, the flooding of the Meirama open pit and the local water balances for different hydrological components.展开更多
基金supported by the National Natural Science Foundation of China(No.50838005,51021140002)the Program for Innovative Research Team in Shaanxi(No.2013KCT-13)
文摘A series of large pilot constructed wetland （CW） systems were constructed near the confluence of an urban stream to a larger fiver in Xi＇an, a northwestern megacity in China, for treating polluted stream water before it entered the receiving water body. Each CW system is a combination of surface- and subsurface-flow cells with local gravel, sand or slag as substrates and Phragmites australis and Typha orientalis as plants. During a one-year operation with an average surface loading of 0.053 m3/（m2.day）, the overall COD, BOD, NH3-N, total nitrogen （TN） and total phosphorus （TP） removals were 72.7% ~ 4.5%, 93.4% ＋ 2.1%, 54.0% ＋ 6.3%, 53.9% ~ 6.0% and 69.4% ：t： 4.6%, respectively, which brought about an effective improvement of the fiver water quality. Surface-flow cells showed better NH3-N removal than their TN removal while subsurface-flow cells showed better TN removal than their NH3-N removal. Using local slag as the substrate, the organic and phosphorus removal could be much improved. Seasonal variation was also found in the removal of all the pollutants and autumn seemed to be the best season for pollutant removal due to the moderate water temperature and well grown plants in the CWs.
基金Supported by the National Natural Science Foundation of China (No.20490206) and the Special Funds for Major State Basic Research Program of China (973 Program, 2004CB217604).
文摘Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded interface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.
基金the National Natural Science Foundation of China (No.20490206) the Special Funds for Major State BasicResearch Program of China (973 Program, 2004CB217604).
文摘Parasitic flows may occur in the numerical simulation of incompressible multiphase flow due to errors in the calculation of surface tension terms, specifically for the curvature and unit normal vector. An improved method for calculating the surface tension based on the level set approach is proposed, in which the contribution of not only the center node but also the rest area of a control volume to the calculation of surface tension is considered in a balanced manner. The weighted integration method (WIM) is more consistent with the concept of a banded in- terface in the level set method. It is applied to the temporal evolution of a two-dimensional neutrally buoyant liquid drop and a buoyancy driven deformable bubble in an immiscible fluid for the validation of WIM. The results show that the parasitic flows are evidently suppressed by the weighted integration method. The weight factors for WIM in 3-D cases are also suggested.
文摘The Singular Integral Operators Method (S.I.O.M.) is applied to the determination of the free-surface profile of an un-steady flow over a spillway, which defines a classical hydraulics problem in open channel flow. Thus, with a known flow rate Q, then the velocities and the elevations are computed on the free surface of the spillway flow. For the numerical evaluation of the singular integral equations both constant and linear elements are used. An application is finally given to the determination of the free-surface profile of a special spillway and comparing the numerical results with corresponding results by the Boundary Integral Equation Method (B.I.E.M.) and by using experiments.
基金supported by the National Natural Science Foundation of China(42107378)the Human Resources-Related Program of Jiangxi Province(20232BCJ23032)+1 种基金the Natural Science Foundation of Jiangxi Province(20232BAB203081)the Science and Technology Project of Jiangxi Province(2022KSG01001,202426ZDKT06).
文摘Due to global warming,changes in the rainfall intensity profile(i.e.,the temporal intensity distribution within a rainfall event)increase the difficulty of accurate erosion prediction and control.Surface cover has been widely used as a critical measure to control soil erosion worldwide.However,the effects of the rainfall intensity profile(RIP)on soil erosion under different surface covers are not fully understood.In this study,long-term in situ field observations of the rain hyetograph,surface runoff coefficient(SRC),subsurface flow rate(SFR),and soil loss rate(SLR)from bare land,litter cover and grass cover were conducted over 11 consecutive years in the red soil hilly region of southern China.According to the occurrence time of the most intense rainfall,226 erosive events were classified into four RIP patterns:advanced,intermediate,delayed,and uniform patterns.The results indicated that the advanced pattern with short durationehigh intensity and the delayed pattern with long durationehigh depth contributed to 73.45% of the total erosive events.For bare land,advanced events were the dominant pattern producing surface runoff and soil erosion,accounting for 57.24%and 75.17%,respectively,of the total surface runoff and erosion.The average SRC and SLR from the advanced pattern were 1.29-2.42 times and 2.52-39.78 times greater than those from the other patterns,respectively.The delayed pattern contributed to subsurface flow,and the average SFR was 1.27-2.17 times greater than that of the other patterns.Furthermore,surface cover significantly reduced surface runoff and erosion and increased subsurface flow,especially under the advanced pattern.Both surface cover measures were equally effective in controlling surface runoff and erosion,but the increase in subsurface flow caused by litter cover was 1.38 -2.67 times greater than that caused by grass cover.Advanced pattern events increase the erosion risk on red soil slopes,and surface cover effectively weakens the effect of variation in the RIP pattern on soil erosion.Moreover,surface cover significantly alters the surface‒subsurface flow distribution pattern for all the RIP patterns.This study highlights the crucial importance of rain intensity profiles on water erosion and provides a basis for optimizing measures to effectively control soil and water loss under climate change.
基金Supported by the University of Mississippi and the USDA Agricultural Research Service
文摘This paper presents a physically-based integrated hydrologic model that can simulate the rain-fall-induced 2D surface water flow, 3D variably saturated subsurface flow, upland soil erosion and transport, and contaminant transport in the surface-subsurface system of a watershed. The model couples surface and subsurface flows based on the assumption of continuity conditions of pressure head and exchange flux at the ground, considering infiltration and evapotranspiration. The upland rill/interrill soil erosion and transport are simulated using a non-equilibrium transport model. Contaminant transport in the integrated surface and subsurface domains is simulated using advection-diffusion equations with mass changes due to sediment sorption and desorption and exchanges between two domains due to infiltration, diffusion, and bed change. The model requires no special treatments at the interface of upland areas and streams and is suitable for wetland areas and agricultural watersheds with shallow streams.
基金Project supported by the National Natural Science Foundation of China(Grant No50779007)
文摘The flow is assumed to be potential, and a boundary integral method is used to solve the Laplace equation for the velocity potential to investigate the shape and the position of the bubble. A 3D code to study the bubble dynamics is developed, and the calculation results agree well with the experimental data. Numerical analyses are carried out for the interaction between multiple bubbles near the free surface including in-phase and out-of-phase bubbles. The calculation result shows that the bubble period increases with the decrease of the distance between bubble centres because of the depression effect between multiple bubbles. The depression has no relationship with the free surface and it is more apparent for out-of-phase bubbles. There are great differences in dynamic behaviour between the in-phase bubbles and the out-of-phase bubbles due to the depression effect. Furthermore, the interaction among eight bubbles is simulated with a three-dlmensional model, and the evolving process and the relevant physical phenomena are presented. These phenomena can give a reference to the future work on the power of bubbles induced by multiple charges exploding simultaneously or continuously.
文摘This research is concerned with new developments and practical applications of a physically-based numerical model that incorporates new approaches for a finite elements solution to the steady/transient problems of the joint ground/surface water flows. Python scripts are implemented in Geographic Information System (GIS) to store, represent and take decisions on the simulated conditions related to the water resources management at the scale of the watershed. The proposed surface-subsurface model considers surface and groundwater interactions to be 2-D horizontally distributed and depth-averaged through a diffusive wave approach for surface flood routing. Infiltration rates, overland flows and evapotranspiration processes are considered by a diffuse discharge from surface water, non-saturated subsoil and groundwater table. Recent developments also allow for the management of surface water flow control through the capacity of diversion on river beds, spillways and outflow operations of floodgates in weirs and dams of reservoirs. Practical application regards the actual hydrology of the Mero River watershed, with two important water bodies mainly concerned with the water resources management at the Cecebre Reservoir and the present flooding of a deep coal mining excavation. The MELEF model (Modèle d’éLéments Fluides, in French) was adapted and calibrated during a period of five years (2008/ 2012) with the help of hydrological parameters, registered flow rates, water levels and registered precipitation, water uses and water management operations in surface and groundwater bodies. The results predict the likely evolution of the Cecebre Reservoir, the flow rates in rivers, the flooding of the Meirama open pit and the local water balances for different hydrological components.
