Hydrologic-hydraulic modeling of sediment transport along the main stem of a watershed: role of tributaries and channel geometry

ABSTRACT

Current estimations of sediment transport at the watershed scale are limited by the difficulty of accurately simulating the sediment transfer along the main stem. The typical approach to simulating watershed sediment transport involves the adoption of hydrologic sediment routing schemes that do not fully capture the contribution and timing of side tributaries, and the inclusion of a simplified channel geometry that does not include its hydraulic feedback. In this paper, we present the results of a coupled hydrologic-hydraulic model of sediment transport applied to a small watershed of Iowa. The model was developed to simulate both the hydrologic network and non-equilibrium sediment transport that occur during a flood. The model results highlight the importance of including side tributaries in order to capture a realistic duration of shear stress that ultimately affects sediment transport. Comparisons with bank erosion measurements indicate that the presented approach is also promising to estimate sediment sources along the main stem.     

Probabilistic evaluation of flood hazard in urban areas using Monte Carlo simulation

The goal of the presented research was the derivation of flood hazard maps, using Monte Carlo simulation of flood propagation at an urban site in the UK, specifically an urban area of the city of Glasgow. A hydrodynamic model describing the propagation of flood waves, based on the De Saint Venant equations in two‐dimensional form capable of accounting for the topographic complexity of the area (preferential outflow paths, buildings, manholes, etc.) and for the characteristics of prevailing imperviousness typical of the urban areas, has been used to derive the hydrodynamic characteristics of flood events (i.e. water depths and flow velocities). The knowledge of the water depth distribution and of the current velocities derived from the propagation model along with the knowledge of the topographic characteristics of the urban area from digital map data allowed for the production of hazard maps based on properly defined hazard indexes. These indexes are evaluated in a probabilistic framework to overcome the classical problem of single deterministic prediction of flood extent for the design event and to introduce the concept of the likelihood of flooding at a given point as the sum of data uncertainty, model structural error and parameterization uncertainty. Copyright © 2011 John Wiley & Sons, Ltd.     

Application of a Monte‐Carlo simulation approach for the probabilistic assessment of seismic hazard for geographically distributed portfolio

The conventional integral approach is very well established in probabilistic seismic hazard assessment (PSHA). However, Monte‐Carlo (MC) simulations can become an efficient and flexible alternative against conventional PSHA when more complicated factors (e.g. spatial correlation of ground shaking) are involved. This study aims at showing the implementation of MC simulation techniques for computing the annual exceedance rates of dynamic ground‐motion intensity measures (GMIMs) (e.g. peak ground acceleration and spectral acceleration). We use multi‐scale random field technique to incorporate spatial correlation and near‐fault directivity while generating MC simulations to assess the probabilistic seismic hazard of dynamic GMIMs. Our approach is capable of producing conditional hazard curves as well. We show various examples to illustrate the potential use of the proposed procedures in the hazard and risk assessment of geographically distributed structural systems. Copyright © 2015 John Wiley & Sons, Ltd.     

Hyperslow drainage in a porous medium: Influence of retention curve

We investigate the process of very slow drainage from pore network simulations in relation with the study of the repository of radioactive waste materials in deep geological formation. The process is characterized by a quite low desaturation of the medium. Considering the special case of a constant stabilizing gradient, we study the sensitivity of the solution given by the classical continuum two-phase flow model to the shape of the retention curve. The solution is found to be very sensitive to the evolution of the retention curve in the range of the high wetting fluid saturations. We show how to construct the retention curve so as to minimize the prediction error with the classical two-phase flow model.     

Application of FLUENT on fine-scale simulation of wind field over complex terrain

The state-of-art Computational Fluid Dynamics （CFD） codes FLUENT is applied in a fine-scale simulation of the wind field over a complex terrain. Several numerical tests are performed to validate the capability of FLUENT on describing the wind field details over a complex terrain. The results of the numerical tests show that FLUENT can simulate the wind field over extremely complex terrain, which cannot be simulated by mesoscale models. The reason why FLUENT can cope with extremely complex terrain, which can not be coped with by mesoscale models, relies on some particular techniques adopted by FLUENT, such as computer-aided design （CAD） technique, unstructured grid technique and finite volume method. Compared with mesoscale models, FLUENT can describe terrain in much more accurate details and can provide wind simulation results with higher resolution and more accuracy.     

汶川 Ms8.0地震同震倾斜应变变化分析

对汶川地震震中附近450 km范围内倾斜仪、应变仪记录到的同震变化做了筛选统计,并根据弹性半无限空间矩形位错模型,分别采用远场地震波资料、近场大地测量资料反演的同震滑移分布模型,正演计算震中附近台站的同震倾斜、应变变化.通过对实际观测与模拟计算结果的比较,初步表明:实际观测的同震倾斜变化方向与模拟计算结果存在75%的一致性,实际观测的同震倾斜变化比模拟值大约1～2个量级;实际观测的同震应变方向与模拟计算结果仅有33.3%的一致性,实际观测的同震应变变化与模拟计算值相当或仅大1个数量级.     

圆筒型FPSO阻尼结构的阻尼性能研究

针对圆筒型海上储油装置FPSO,研究其不同类型阻尼结构的阻尼性能,以便有效减小此类结构的垂荡运动。提出U型和倒U型两种不同的圆筒型FPSO阻尼结构型式,采用CFD方法模拟浮体的强迫运动,计算分析阻尼结构的阻尼特性,对比了传统垂荡板结构、U型阻尼结构、倒U型阻尼结构的阻尼性能,详细分析阻尼结构的不同参数包括高度、筒体间距、开孔情况等三个因素对于阻尼特性的影响,基于CFD方法,分析了不同阻尼结构阻尼形成的机理。计算分析表明,倒U型阻尼性能最好,U型和倒U型阻尼结构与筒体间距增大后阻尼增大,阻尼结构高度增加或开孔增加后阻尼降低。     

A multi-model study of the restratification phase in an idealized convection basin

The representation of baroclinic instability in numerical models depends strongly upon the model physics and significant differences may be found depending on the vertical discretization of the governing dynamical equations. This dependency is explored in the context of the restratification of an idealized convective basin with no external forcing. A comparison is made between an isopycnic model including a mixed layer (the Miami Isopycnic Coordinate Ocean Model, MICOM), its adiabatic version (MICOM-ADIAB) in which the mixed layer physics are removed and the convective layer is described by a deep adiabatic layer outcropping at the surface instead of a thick dense mixed layer, and a z-coordinate model (OPA model).In the absence of a buoyancy source at the surface, the mixed layer geometry in MICOM prevents almost any retreat of this layer. As a result, lateral heat exchanges in the upper layers are limited while mass transfers across the outer boundary of the deep convective mixed layer result in an unrealistic outward spreading of this layer. Such a widespread deep mixed layer maintains a low level of baroclinic instability, and therefore limits lateral heat exchanges in the upper layers over most of the model domain. The behavior of the adiabatic isopycnic model and z-coordinate model is by far more satisfactory although contrasted features can be observed between the two simulations. In MICOM-ADIAB, the more baroclinic dynamics introduce a stronger contrast between the surface and the dense waters in the eddy kinetic energy and heat flux distributions. Better preservation of the density contrasts around the dense water patch maintains more persistent baroclinic instability, essentially associated with the process of dense water spreading. The OPA simulation shows a faster growth of the eddy kinetic energy in the early stages of the restratification which is attributed to more efficient baroclinic instability and leads to the most rapid buoyancy restoring in the convective area among the three simulations. Dense water spreading and warm surface capping occur on fairly similar time scales in MICOM-ADIAB although the former is more persistent that the latter. In this model, heat is mainly transported by anticyclonic eddies in the dense layer while both cyclonic and anticyclonic eddies are involved in the upper layers. In OPA, heat is mainly brought into the convective zone through the export of cold water trapped in cyclonic eddies with a strong barotropic structure. Probably the most interesting difference between the z-coordinate and the adiabatic isopycnic model is found in the temperature distribution ultimately produced by the restratification process. OPA generates a spurious volume of intermediate water which is not seen in MICOM-ADIAB where the volume of the dense water is preserved.     

Discrepancies in Simulated Ocean Net Surface Heat Fluxes over the North Atlantic

The change in ocean net surface heat flux plays an important role in the climate system. It is closely related to the ocean heat content change and ocean heat transport, particularly over the North Atlantic, where the ocean loses heat to the atmosphere, affecting the AMOC(Atlantic Meridional Overturning Circulation) variability and hence the global climate.However, the difference between simulated surface heat fluxes is still large due to poorly represented dynamical processes involving multisca...     

The Effects of Canopy Leaf Area Index on Airflow Across Forest Edges: Large-eddy Simulation and Analytical Results

The structure of turbulent flows along a transition between tall-forested canopies and forest clearings continues to be an active research topic in canopy turbulence. The difficulties in describing the turbulent flow along these transitions stem from the fact that the vertical structure of the canopy and its leaf area distribution cannot be ignored or represented by an effective roughness length. Large-eddy simulation (LES) runs were performed to explore the effect of a homogeneous variation in the forest leaf area index (LAI) on the turbulent flow across forest edges. A nested grid numerical method was used to ensure the development of a deep boundary layer above the forest while maintaining a sufficiently high resolution in the region close to the ground. It was demonstrated that the LES here predicted first-order and second-order mean velocity statistics within the canopy that agree with reported Reynolds-Averaged Navier–Stokes (RANS) model results, field and laboratory experiments. In the simulations reported here, the LAI was varied between 2 and 8 spanning a broad range of observed LAI in terrestrial ecosystems. By increasing the forest LAI, the mean flow properties both within the forest and in the clearing near the forest edge were altered in two fundamental ways: near the forest edge and into the clearing, the flow statistical properties resembled the so-called back-facing step (BFS) flow with a mean recirculation zone near the edge. Another recirculation zone sets up downstream of the clearing as the flow enters the tall forest canopy. The genesis of this within-forest recirculation zone can be primarily described using the interplay between the mean pressure gradients (forcing the flow) and the drag force (opposing the flow). Using the LES results, a simplified analytical model was also proposed to explain the location of the recirculation zone inside the canopy and its dependence on the forest LAI. Furthermore, a simplified scaling argument that decomposes the mean velocity at the outflow edge into a superposition of ‘exit flow’ and BFS-like flow with their relative importance determined by LAI was explored.