COMPUTATION OF UNSTEADY FLOWS IN THE ALABAMA RIVER1

ABSTRACT: An application is described of the branch-network flow model, BRANCH, to the upper Alabama River system in central Alabama. The model is used to simulate one-dimensional unsteady flows and water surface elevations in approximately 60 river miles of the Alabama River system. Preliminary calibration was made using 72 hours of observed data. Simulated discharges are about 10 percent lower than observed discharges at higher discharge rates and computed flows lag observed flows by about 30 minutes.     

MEASURED AND PREDICTED VELOCITY AND LONGITUDINAL DISPERSION AT STEAI)Y AND UNSTEADY FLOW,COLORADO RIVER,GLEN CANYON DAM TO LAKE MEAD1

ABSTRACT: The effect of unsteadiness of dam releases on velocity and longitudinal dispersion of flow was evaluated by injecting a fluorescent dye into the Colorado River below Glen Canyon Dam and sampling for dye concentration at selected sites downstream. Measurements of a 26-kilometer reach of Glen Canyon, just below Glen Canyon Dam, were made at nearly steady dam releases of 139, 425, and 651 cubic meters per second. Measurements of a 380-kilometer reach of Grand Canyon were made at steady releases of 425 cubic meters per second and at unsteady releases with a daily mean of about 425 cubic meters per second. In Glen Canyon, average flow velocity through the study reach increased directly with discharge, but dispersion was greatest at the lowest of the three flows measured. In Grand Canyon, average flow velocity varied slightly from subreach to subreach at both steady and unsteady flow but was not significantly different at steady and unsteady flow over the entire study reach. Also, longitudinal dispersion was not significantly different during steady and unsteady flow. Long tails on the time-concentration curves at a site, characteristic of most rivers but not predicted by the one-dimensional theory, were not found in this study. Absence of tails on the curves shows that, at the measured flows, the eddies that are characteristic of the Grand Canyon reach do not trap water for a significant length of time. Data from the measurements were used to calibrate a one-dimensional flow model and a solute-transport model. The combined set of calibrated flow and solute-transport models was then used to predict velocity and dispersion at potential dam-release patterns.     

FLOOD ROUTING THROUGH A FLAT,COMPLEX FLOODPLAIN USING A ONE-DIMENSIONAL UNSTEADY FLOW COMPUTER PROGRAM1

ABSTRACT: The routing of flood waves through the Central Basin of the Passaic River in New Jersey is complex because of flat gradients and flow reversals. The one-dimensional unsteady flow program DWOPER, developed by the National Weather Service, was used to simulate flood wave movement through the Basin. A historical event was used for calibration and two synthetic events were simulated. Boundary conditions consisted of discharge hydrographs at inflow points to the study area, local flow hydrographs at interior points, and a stage discharge relation for flow over the crest of a diversion dam at the basin outlet. Manning's n values were adjusted based on stage and discharge data for the historical event; however, verification data were not available for events comparable in magnitude to the synthetic events. Aspects of the investigation reported include techniques for characterizing the flow system, model calibration, techniques for representing a tunnel diversion, and simulation results.     

CRITICAL EVALUATION OF CERTAIN METHODS OF UNSTEADY GROUNDWATER HYDRAULICS1

The basic theories and fundamental assumptions usually employed in the solution of unsteady groundwater flow problems are reviewed critically. The best known method of analysis for such problems is based on the Dupuit-Forchheimer approximation and leads to a nonlinear parabolic differential equation which is generally solved by linearization or numerical methods. The accuracy of the solution to this equation can be improved by use of a different approach which does not employ the Dupuit Forchheimer assumption, but rather is based on a semi-numerical solution of the Laplace equation for quasi-steady conditions. The actual unsteady process is replaced by a sequence of steady-state conditions, and it is assumed that the actual unsteady flow characteristics during a short time interval can be approximated by those associated with “average” steady state flow. The Laplace equation is solved by a semi-discretization method according to which the horizontal coordinate is divided into subintervals, while the vertical coordinate is maintained continuous. The proposed method is applied to a typical tile drainage problem, and, based on a comparison of calculated results with experimental data, the method is evaluated and practical conclusions regarding its applicability are advanced.     

EFFECTS OF TIME STEP SIZE IN IMPLICIT DYNAMIC ROUTING1

ABSTRACT The effects of the size of the Δt time step used in the integration of the implicit difference equations of unsteady open-channel flow are determined for numerous typical hydrographs with durations in the order of days or even weeks. Truncation errors related to the size of the Δt time step cause a numerical distortion (dispersion and attenuation) of the computed transient. The magnitude of the distortion is related directly to the size of the time step, the length of channel reach, and the channel resistance and inversely to the time of rise of the hydrograph. The type of finite difference expression which replaces spatial derivatives and non-derivative terms in the partial differential equations of unsteady flow has an important influence on the magnitude of the numerical distortion, as well as the numerical stability of the implicit difference equations. Time step sizes in the range of 3 to 6 hrs generally tend to minimize the combination of required computation time and numerical distortion of transients having a time of rise of the order of several days.     

基于MIKE 21的长江江山制药排污口水环境影响分析——以江山制药公司污水处理站为例

为分析长江沿岸入河排污口设置对水环境的影响,采用MIKE 21建立了长江马鞍山-高桥段二维非稳态水量水质模型,并以江山制药公司污水处理站为例,对其排污口迁址改建前后的水环境影响进行模拟.结果 表明,较于原排污口排放,拟设排污口在执行新排放标准情况下,对九圩港入江断面水质有改善作用,COD、NH3-N、TP的改善率分别为...     

基于非恒定流的污染物迁移扩散随机模拟

针对目前大部分河流都受到不同程度污染的现状,研究了基于非恒定流的污染物迁移扩散随机模拟方法。通过一维非恒定流数值计算和流场模拟,得到了水库开闸放水后下游各个断面的水位、流量、流速等数据,并以此为基础,利用蒙特卡洛法对水库开闸放水前后下游典型断面的COD和NH3-N等进行污染物迁移扩散随机模拟,并将模型应用于文峪河下游。结果表明,文峪河水库开闸放水后,COD和NH3-N的污染带长度和宽度均有不同程度的减小,污染带向排污一侧岸边集中。研究成果可为解决河流污染现状与水库调度提供参考。     

非稳态条件下人工湿地水流规律及其对水力停留时间的影响

由于降水的随机性,处理农业或城市径流的人工湿地通常在进水水量波动的非稳态条件下工作.为研究在非稳态条件下人工湿地的水流规律,在基于时间轴的多点示踪剂试验的基础上,利用RTD(水力停留时间分布)探讨了单峰和双峰两种水力冲击非稳态条件下水平潜流人工湿地中的水流规律.结果表明,非稳态过程的水量变化对不同时刻流入人工湿地的液体的RTD影响明显.用基于出水流量的非线性无量纲时间变量(φ)对原始RTD进行归一化处理,所得的RTD浓度曲线[C'(φ)]在整个非稳态过程中相对稳定,具有相似的倾斜度、峰值时间(单峰冲击下为0.83~1.00,双峰冲击下为0.87~1.08)、重心位置(单峰冲击下为1.08~1.25,双峰冲击下为1.05~1.25)及相关统计学参数,表现出一种独立于流量变化之外的相对稳定的特征水流规律.与稳态基流量下的t_mean(平均停留时间)相比,水力冲击使人工湿地t_mean呈冲击前持续缩短、冲击后又快速回升的趋势(单峰、双峰冲击下最大值与最小值相差分别约为30和32 h);而t_mean'(归一化平均停留时间)在非稳态过程中却相对稳定(单峰、双峰冲击下最大值与最小值相差分别约为0.15和0.20),适合作为非稳态下表征人工湿地水流规律和停留时间的特征参数.研究显示,归一化RTD易于模拟,可用其描述非稳态下人工湿地的水流规律以及水流规律对污染物接触反应时间的影响.      

NONSTEADY FLOW TO A WELL WITH TIME DEPENDENT DRAWDOWN1

Confined flow toward a single well of finite radius in an extensive aquifer of uniform transmissibility is studied under the assumption of time-dependent drawdown. Three particular cases are considered: (a) linear drawdown (including constant drawdown); (b) exponential drawdown; (c) periodic (sinusoidal) drawdown. The differential equation governing unsteady axial symmetric flow toward a single well in a confined aquifer is solved for the three different situations by the use of the Laplace transform method. The resulting expressions are integrated by adapting a modified Gemant scheme. General computer programs have been developed and operated for several combinations of characteristics. The results are plotted to show the effect of time dependent drawdown on the variation of the well discharge and the piezometric head distribution.     

微重条件下非稳态燃烧过程的数值模拟

在建立气固相界面耦合燃烧理论模型的基础上，对微重条件下火焰在固相材料表面的蔓延进行了三维非稳态数值模拟，模型中考虑了气相燃烧过程、固相的热解和传热过程，给出了相界面的处理方法。不同算例的计算研究了环境氧浓度、环境压力对微重火蔓延的影响，并与正常重力条件下的计算结果进行了分析对比．结果基本合理。