共查询到20条相似文献，搜索用时 31 毫秒
1.
This paper presents the numerical investigation of the distribution of inertial sediment particles in the turbulent boundary layer of an open channel flow with the particle Stokes number ranging from 0.6 to 20.4. The methodology is a combination of three numerical approaches, i.e. direct numerical simulation of turbulent flow, the point-particle immersed boundary method, and the discrete particle method. By applying the Vorono? analysis, the preferential concentration characteristics of sediment particles were investigated quantitatively. It was found that the normalized area of the Vorono? cells follows a log-normal particle distribution. The inertial sediment particles distributed unevenly in the turbulent boundary layer and the unevenness, governed by the particle Stokes number, was more significant as the particle Stokes number approaches unity. The inertial sediment particles in the turbulent boundary layer accumulated preferentially in streamwise-aligned streaky structures and this pattern was less significant with increasing particle Stokes number. 相似文献
2.
《国际泥沙研究》2017,(3)
This paper presents the numerical investigation of the distribution of inertial sediment particles in the turbulent boundary layer of an open channel flow with the particle Stokes number ranging from 0.6 to 20.4. The methodology is a combination of three numerical approaches, i.e. direct numerical simulation of turbulent flow, the point-particle immersed boundary method, and the discrete particle method. By applying the Vorono analysis, the preferential concentration characteristics of sediment particles were investigated quantitatively. It was found that the normalized area of the Voronoi cells follows a lognormal particle distribution. The inertial sediment particles distributed unevenly in the turbulent boundary layer and the unevenness, governed by the particle Stokes number, was more significant as the particle Stokes number approaches unity. The inertial sediment particles in the turbulent boundary layer accumulated preferentially in streamwise-aligned streaky structures and this pattern was less significant with increasing particle Stokes number. 相似文献
3.
Renske C. Terwisscha van Scheltinga Giovanni Coco Heide Friedrich 《地球表面变化过程与地形》2019,44(10):1881-1894
This work presents measurements and analysis of sand particle velocities over a subaqueous dune with median sand diameter of 0.85 mm. Time-lapse images of the mobile bed and an automated particle image velocimetry (PIV)-based cross-correlation method are used to obtain mean velocity of sand particles. This technique is shown to be consistent with measurements obtained with manual tracing. The measurements indicate an increase in mean particle velocity over a dune slope. Three regions are distinguished over the dune slope: (1) region of fluctuating particle velocity, (2) region of increasing particle velocity, and (3) region of maximum particle velocity. The observations are aligned with experimental and numerical modelling studies, indicating fluctuations in flow velocity over a dune stoss slope. We furthermore show that the standard deviation of the mean particle velocity is affected by the slope location and decreases from the lower slope towards the upper slope. The particle velocity variability is discussed in the context of general onset and cessation of sediment transport, the effect of the reattachment zone, sweep-transport events, and the existence of superimposed bedforms. With this work we bridge the gap between measurements of bedload transport at the particle-scale and at the bedform-scale. © 2019 John Wiley & Sons, Ltd. 相似文献
4.
Simonetta Rubol Daniele Tonina Lionel Vincent Jill A. Sohm William Basham Ralph Budwig Pratixa Savalia Eva Kanso Douglas G. Capone Kenneth H. Nealson 《水文研究》2018,32(3):402-407
We describe a novel inexpensive method, utilizing particle image velocimetry (PIV) and refractive index‐matching (RIM) for visualizing and quantifying the flow field within bio‐amended porous media. To date, this technique has been limited to idealized particles, whose refractive index does not match that of fresh water and thus requires specialized and often toxic or hazardous fluids. Here, we use irregularly shaped grains made of hydrogel as the solid matrix and water as the fluid. The advantage of using water is that it provides, for the first time, the opportunity to study both hydraulic and biological processes, which typically occur in soils and streambeds. By using RIM coupled with PIV (RIM‐PIV), we measured the interstitial flow field within a cell packed with granular material consisting of hydrogel grains in a size range of 1–8 mm, both in the presence and in the absence of Sinorhizobium meliloti bacteria (strain Rm8530). We also performed experiments with fluorescent tracer (fluorescein) and fluorescent microbes (Shewanella GPF MR‐1) to test the capability of visualizing solute transport and microbial movements. Results showed that the RIM‐PIV can measure the flow field for both biofilm‐free and biofilm‐covered hydrogel grains. The fluorescent tracer injection showed the ability to visualize both physical (concave surfaces and eddies) and biological (biofilms) transient storage zones, whereas the fluorescent microbe treatment showed the ability to track microbial movements within fluids. We conclude that the proposed methodology is a promising tool to visualize and quantify biofilm attachment, growth, and detachment in a system closer to natural conditions than a 2D flow cell experiment. 相似文献
5.
We present in this note an innovative software environment, called Interactive Ground Water (IGW), for unified deterministic and stochastic ground water modeling. Based on efficient computational algorithms, IGW allows simulating three-dimensional (3D) unsteady flow and transport in saturated media subject to systematic and "random" stresses and geological and chemical heterogeneity. Adopting a new computing paradigm, IGW eliminates the fragmentation in the traditional modeling schemes and allows fully utilizing today's dramatically increased computing power. For many problems, IGW enables real-time modeling, visualization, mapping, and analysis. The software environment functions as a "numerical laboratory" in which an investigator may freely explore the following: creating visually an aquifer system of desired configurations, interactively applying stresses and boundary conditions, and then investigating and visualizing on the fly the geology and flow and transport dynamics. At any time, a researcher can pause to interact dynamically with virtually any aspects of the modeling process and then resume the integrated visual exploration; he or she can initiate, pause, or resume particle tracking, plume modeling, subscale modeling, stochastic modeling, monitoring, and budget analyses. IGW continually provides results that are dynamically processed, overlaid, and displayed. It dynamically merges modeling inputs and outputs into composite two-dimensional/3D images-integrating related data to provide a more complete view of the complex interplay among the geology, hydrology, flow system, and transport. These unique capabilities of real-time modeling, steering, analysis, and mapping expand the utility of models as tools for research, education, and professional investigations. 相似文献
6.
A Lagrangian particle tracking scheme has been extended to simulate advective transport through coupled surface and subsurface flows. This extended scheme assumes a two-dimensional flow field for the overland domain and a three-dimensional flow field for the subsurface domain. Moreover it is assumed that the flow fields are simulated by a cell centered finite difference method. Pathlines in both the subsurface and the overland domain are simulated by classical particle tracking methods. Exchange of particles between the two domains is simulated by newly-developed algorithms presented in this study. Different algorithms are used depending on the direction of the exchange across the interface between the two domains. In the subsurface domain knowledge about a particle’s pathline is enough to detect a transfer to the surface domain and the solution is straightforward. However, in the two-dimensional overland domain pathlines are parallel to the land surface. Therefore the velocity field in the overland domain alone cannot be used to detect a transfer to the subsurface. We propose a relatively simple algorithm to estimate the probability of transfer to the subsurface domain. It is shown that this algorithm can also be used to handle the transfer from the overland domain to the atmosphere domain. The algorithm to estimate the transfer probabilities is based on the mass balance of water on a streamtube aligned with the particle’s pathline. This newly developed technique ensures that transit time distributions can be simulated accurately. These new relationships are implemented in an existing particle tracking code and are verified using analytical solutions for transit times. 相似文献
7.
Modeling of suspended sediment particle movement in surface water can be achieved by stochastic particle tracking model approaches.In this paper,different mathematical forms of particle tracking models are introduced to describe particle movement under various flow conditions,i.e.,the stochastic diffusion process,stochastic jump process,and stochastic jump diffusion process.While the stochastic diffusion process can be used to represent the stochastic movement of suspended particles in turbulent flows,the stochastic jump and the stochastic jump diffusion processes can be used to describe suspended particle movement in the occurrences of a sequence of extreme flows.An extreme flow herein is defined as a hydrologic flow event or a hydrodynamic flow phenomenon with a low probability of occurrence and a high impact on its ambient flow environment.In this paper,the suspended sediment particle is assumed to immediately follow the extreme flows in the jump process（i.e.the time lag between the flow particle and the sediment particle in extreme flows is considered negligible）.In the proposed particle tracking models,a random term mainly caused by fluid eddy motions is modeled as a Wiener process,while the random occurrences of a sequence of extreme flows can be modeled as a Poisson process.The frequency of occurrence of the extreme flows in the proposed particle tracking model can be explicitly accounted for by the Poisson process when evaluating particle movement.The ensemble mean and variance of particle trajectory can be obtained from the proposed stochastic models via simulations.The ensemble mean and variance of particle velocity are verified with available data.Applicability of the proposed stochastic particle tracking models for sediment transport modeling is also discussed. 相似文献
8.
A novel technique for visualizing turbulent flow data from a gravel-bed river is presented. The time development of flow velocity and shear stress at three heights is displayed using a computer program. This can be used to observe how the fluctuations of velocity and shear stress interact both spatially and temporally. We highlight examples of flow events which are important for the understanding of flow dynamics. The visualization suggests that the turbulent flow-field is characterized not only by coherence over time at a point, but also by spatial interdependence between points. We suggest that this new visualization approach will assist further interpretation of statistical analyses of turbulent signals, as well as focusing future measurement strategies by providing a clearer spatio-temporal picture of the flow structure. 相似文献
9.
Nafiz Maden 《Pure and Applied Geophysics》2010,167(3):349-358
Curie-point depth and heat flow values of the Erciyes region are determined to identify the thermal regime of the Central Anatolia by applying the spectral analysis method to the magnetic anomaly data. To compute the spectrum of the data, the magnetic anomaly of the region is transformed into 2-D Fourier domain to attain the average Curie depth. This method is useful in determining the top boundary of magnetic anomaly sources and reveals the Curie depth as 13.7 km in the study area. The obtained results imply a high thermal gradient (42.3°C km?1) and corresponding heat flow values (88.8 mWm?2) in the research area. Using the temperature value measured at borehole drilled by the General Directorate of Mineral Research and Exploration of Turkey (MTA), the values for the thermal gradient and heat flow value were computed as 50.7°C km?1, 106.5 mWm?2. From the heat flow value, the Curie-point depth was determined as 11.4 km in this region. It is concluded from the obtained values that the region has very high geothermal potential caused by partial melting of the lower crust. 相似文献
10.
Optimization of groundwater and other subsurface resources requires analysis of multiple‐well systems. The usual modeling approach is to apply a linear flow equation (e.g., Darcy's law in confined aquifers). In such conditions, the composite response of a system of wells can be determined by summating responses of the individual wells (the principle of superposition). However, if the flow velocity increases, the nonlinear losses become important in the near‐well region and the principle of superposition is no longer valid. This article presents an alternative method for applying analytical solutions of non‐Darcy flow for a single‐ to multiple‐well systems. The method focuses on the response of the central injection well located in an array of equally spaced wells, as it is the well that exhibits the highest pressure change within the system. This critical well can be represented as a single well situated in the center of a closed square domain, the width of which is equal to the well spacing. It is hypothesized that a single well situated in a circular region of the equivalent plan area adequately represents such a system. A test case is presented and compared with a finite‐difference solution for the original problem, assuming that the flow is governed by the nonlinear Forchheimer equation. 相似文献
11.
Water Resources - This study explores the influences of flow discharge and particle size on bedload transport trajectory by applying a depth-averaged two-dimensional model to a 110°... 相似文献
12.
Gauging extreme floods on YouTube: application of LSPIV to home movies for the post‐event determination of stream discharges 
下载免费PDF全文
下载免费PDF全文 Raphaël Le Boursicaud Lionel Pénard Alexandre Hauet Fabien Thollet Jérôme Le Coz 《水文研究》2016,30(1):90-105
Movies taken by witnesses of extreme flood events are increasingly available on video sharing websites. They potentially provide highly valuable information on flow velocities and hydraulic processes that can help improve the post‐flood determination of discharges in streams and flooded areas. We investigated the troubles and potential of applying the now mature large‐scale particle image velocimetry (LSPIV) technique to such flood movies that are recorded under non‐ideal conditions. Processing was performed using user‐friendly, free software only, such as Fudaa‐LSPIV. Typical issues related to the image processing and to the hydrological analysis are illustrated using a selected example of a pulsed flash‐flood flow filmed in a mountainous torrent. Simple corrections for lens distortion (fisheye) and limited incoherent camera movement (shake) were successfully applied, and the related errors were reduced to a few percents. Testing the different image resolution levels offered by YouTube showed that the difference in time‐averaged longitudinal velocity was less than 5% compared with full resolution. A limited number of GRPs, typically 10, is required, but they must be adequately distributed around the area of interest. The indirect determination of the water level is the main source of uncertainty in the results, usually much more than errors because of the longitudinal slope and waviness of the free‐surface of the flow. The image‐based method yielded direct discharge estimates of the base flow between pulses, of the pulse waves, and of the time‐averaged flow over a movie sequence including a series of five pulses. A comparison with traditional indirect determination methods showed that the critical‐depth method may produce significantly biassed results for such a fast, unsteady flow, while the slope‐area method seems to be more robust but would overestimate the time‐averaged flow rate if applied to the high‐water marks of a pulsed flow. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
13.
Koji Matsuki Katsuya Nakatani Takashi Arai Kazuo Ohmura Ryuji Takeuchi Yasushi Arai Shinji Takeuchi 《Journal of Hydrology》2008,360(1-4):217-229
By modifying a previous method with constant elements, we developed a quadratic element method for more accurately estimating groundwater flow by the inversion of tilt data. In this method: (1) a region of groundwater flow is divided into quadratic elements in which the change in groundwater volume per unit volume of rock (Δv) and the Skempton coefficient (B) vary in a quadratic manner with the coordinates, (2) the values of Δv are set to zero at the boundaries of the region of groundwater flow and (3) the sum of the squared second derivatives of Δv is adopted as a constraining condition that is weighted and added to the sum of the squared errors in tilt. First, analyses were performed for a flow model to determine the accuracy of this method for estimating groundwater flow and also to clarify the effect of the assumed size of a region of groundwater flow. These analyses showed that the quadratic element method proposed in this study gives a much better estimation of Δv than the constant element method and that a large region of groundwater flow should be assumed, rather than a small region, since the values of Δv at points outside of the actual region of groundwater flow are estimated to be nearly zero when a large region is assumed while these values are greatly overestimated when an excessively small region is assumed. Finally, the quadratic element method was applied to the site of the Mizunami Underground Research Laboratory in the Tono area, Japan. Inverse analyses were performed for tilt data measured by four tiltmeters with a resolution of 10−9 radians during the excavation of two shafts under the assumption that the rock mass is an isotropic and homogeneous half- space. The results showed that the method proposed in this study reproduced the tilt data very accurately. Thus, the distribution of Δv was estimated without sacrificing the reproducibility of the tilt data. The contour maps of B(1 + ν)Δv (ν: Poisson’s ratio) showed that the heterogeneous flow of groundwater occurred at the site and that groundwater volume decreased mainly in the area surrounded by two faults. The latter result is consistent with the finding obtained by previous investigations that these faults have low permeability in the direction perpendicular to the strike and may act as a flow barrier. 相似文献
14.
15.
LinTRODUCTIONANDAIMSTheaimsofthisStUdyaretointroducenewtechniquesforanalysingfloodflowsfromsedimentarrangementmsapproachincorporatestheeffectsofflotvonthedevelopmentofrouglmessontheriverbed,andnot,asismostlythecase,theeffectsofbedroughnessonflow.Differentacsofflows,includingthemaalmumfloodflowarereconsmictedfromfloworientahonsofparticlesoveranentiregravCIbar.DuetOtheddricultiesinvolvedinobtaininglargeparticlesamples,studiesonpAncleOrientationsintheliteratUrearerare.WhilstmoStanalyse… 相似文献
16.
In this paper we describe the transport of pollution in groundwater in the neighbourhood of a well in a uniform background flow. We compute the rate at which contaminated particles reach the well as a function of the place of the source of pollution. The motion of a particle in a dispersive flow is seen as a random walk process. The Fokker-Planck equation for the random motion of a particle is transformed using the complex potential for the advective flow field. The resulting equation is solved asymptotically after a stretching transformation. Finally, the analytical solution is compared with results from Monte Carlo simulations with the random walk model. The method can be extended to arbitrary flow fields. Then by a numerical coordinate transformation the analytical results can still be employed. 相似文献
17.
Abstract The problem of non-steady flow of water in a soil-plant system can be described by adding a sink term to the continuity equation for soil water flow. In this paper the sink term is defined in two different ways. Firstly it is considered to be dependent on the hydraulic conductivity of the soil, on the difference in pressure head between the soil and the root-soil interface and some root effectiveness function. Secondly the sink is taken to be a prescribed function of the soil water content. The partial differential equation applying to the first problem is solved by both a finite difference (FD 1) and a finite element (FE 1) technique, that applying to the second problem by a finite difference approach (FD 2). The purpose of this paper is to verify the numerical models against field measurements, to compare the results obtained by the three numerical methods and to show how the finite element method can be applied to complex but realistic two-dimensional flow situations. Two examples are given. The first concerns one-dimensional flow and it compares numerical results with those obtained experimentally in the field from water balance studies on red cabbage (Brassica oleracea L. ‘Rode Herfst’) grown on a clay soil in the presence of a water table. The second example describes two-dimensional flow in a complex field situation in the Netherlands where flow takes place under cropped field conditions through five anisotropic layers. Water is supplied to the system by infiltration from two unlined ditches and is withdrawn from the system by evapotranspiration and by leakage to an underlying pumped aquifer. 相似文献
18.
V. V. Pak 《Izvestiya Physics of the Solid Earth》2007,43(9):783-790
A numerical compaction model of a fluid in a viscous skeleton is developed with regard for a phase transition. The temperatures of phases are different. The solution is found by the method of asymptotic expansion relative to the incompressible variant, which removes a number of computational problems related to the weak compressibility of the skeleton. For each approximation, the problem is solved by the finite element method. The process of 2-D compaction of a magmatic melt in the asthenosphere under a fault zone is examined for one-and two-temperature cases. The magmatic flow concentrates in this region due to a lower pore pressure. Higher temperature magma entering from lower levels causes a local heating of the skeleton and intense melting of its fusible component. In the two-temperature model, a magma concentration anomaly develops under the fault zone. The fundamental limitations substantially complicating the corresponding calculations within the framework of a one-temperature model are pointed out and the necessity of applying a multitemperature variant is substantiated. 相似文献
19.
The work describes experimental observations of enhancements in the electron density of the ionospheric F-region created by cusp/cleft particle precipitation at the dayside entry to the polar-cap convection flow. Measurements by meridian scanning photometer and all-sky camera of optical red-line emissions from aurora are used to identify latitudinally narrow bands of soft-particle precipitation responsible for structured enhancements in electron density determined from images obtained by radio tomography. Two examples are discussed in which the electron density features with size scales and magnitudes commensurate with those of patches are shown to be formed by precipitation at the entry region to the anti-sunward flow. In one case the spectrum of the incoming particles results in ionisation being created, for the most part below 250 km, so that the patch will persist only for minutes after convecting away from the auroral source region. However in a second example, at a time when the plasma density of the solar wind was particularly high, a substantial part of the particle-induced enhancement formed above 250 km. It is suggested that, with the reduced recombination loss in the upper F-region, this structure will retain form as a patch during passage in the anti-sunward flow across the polar cap. 相似文献
20.
This study formulates and analyzes continuous time random walk (CTRW) models in radial flow geometries for the quantification of non-local solute transport induced by heterogeneous flow distributions and by mobile–immobile mass transfer processes. To this end we derive a general CTRW framework in radial coordinates starting from the random walk equations for radial particle positions and times. The particle density, or solute concentration is governed by a non-local radial advection–dispersion equation (ADE). Unlike in CTRWs for uniform flow scenarios, particle transition times here depend on the radial particle position, which renders the CTRW non-stationary. As a consequence, the memory kernel characterizing the non-local ADE, is radially dependent. Based on this general formulation, we derive radial CTRW implementations that (i) emulate non-local radial transport due to heterogeneous advection, (ii) model multirate mass transfer (MRMT) between mobile and immobile continua, and (iii) quantify both heterogeneous advection in a mobile region and mass transfer between mobile and immobile regions. The expected solute breakthrough behavior is studied using numerical random walk particle tracking simulations. This behavior is analyzed by explicit analytical expressions for the asymptotic solute breakthrough curves. We observe clear power-law tails of the solute breakthrough for broad (power-law) distributions of particle transit times (heterogeneous advection) and particle trapping times (MRMT model). The combined model displays two distinct time regimes. An intermediate regime, in which the solute breakthrough is dominated by the particle transit times in the mobile zones, and a late time regime that is governed by the distribution of particle trapping times in immobile zones. These radial CTRW formulations allow for the identification of heterogeneous advection and mobile-immobile processes as drivers of anomalous transport, under conditions relevant for field tracer tests. 相似文献