Capturing watershed water balance with a physically-based two-hydrological-variable model: Application to the Little Washita basin

Author(s):  
Fanny Picourlat ◽  
Emmanuel Mouche ◽  
Claude Mugler

<p>Hydrological processes import across scales is known to constitute a key challenge to improve their representation in large-scale land surface models. Since these models describe continental hydrology with vertical one dimensional infiltration and evapotranspiration, the challenge mainly resides in the dimensionality reduction of the processes. Departing from the catchment three-dimensional scale, previous work has shown that an equivalent two-dimensional hillslope model is able to simulate long term watershed water balance with good accuracy. This work has been done on the Little Washita basin (Ok, USA) using the integrated code HydroGeoSphere. Following this framework, we show that hillslope hydrology can be described by using realistic simplifying assumptions, such as linear water table profile. These assumptions allow the writing of an analytical model relying on two hydrological variables: the seepage face extension, which describe the intersection length between the water table and the land surface, and the water table slope. The last step of the work will be to use these key variables and this simplified description of the driving processes for importing small-scale hydrological processes into large-scale models.</p>

2008 ◽  
Vol 615 ◽  
pp. 371-399 ◽  
Author(s):  
S. DONG

We report three-dimensional direct numerical simulations of the turbulent flow between counter-rotating concentric cylinders with a radius ratio 0.5. The inner- and outer-cylinder Reynolds numbers have the same magnitude, which ranges from 500 to 4000 in the simulations. We show that with the increase of Reynolds number, the prevailing structures in the flow are azimuthal vortices with scales much smaller than the cylinder gap. At high Reynolds numbers, while the instantaneous small-scale vortices permeate the entire domain, the large-scale Taylor vortex motions manifested by the time-averaged field do not penetrate a layer of fluid near the outer cylinder. Comparisons between the standard Taylor–Couette system (rotating inner cylinder, fixed outer cylinder) and the counter-rotating system demonstrate the profound effects of the Coriolis force on the mean flow and other statistical quantities. The dynamical and statistical features of the flow have been investigated in detail.


1999 ◽  
Vol 382 ◽  
pp. 307-329 ◽  
Author(s):  
JUDITH K. FOSS ◽  
K. B. M. Q. ZAMAN

The large- and small-scale vortical motions produced by ‘delta tabs’ in a two-stream shear layer have been studied experimentally. An increase in mixing was observed when the base of the triangular shaped tab was affixed to the trailing edge of the splitter plate and the apex was pitched at some angle with respect to the flow axis. Such an arrangement produced a pair of counter-rotating streamwise vortices. Hot-wire measurements detailed the velocity, time-averaged vorticity (Ωx) and small-scale turbulence features in the three-dimensional space downstream of the tabs. The small-scale structures, whose scale corresponds to that of the peak in the dissipation spectrum, were identified and counted using the peak-valley-counting technique. The optimal pitch angle, θ, for a single tab and the optimal spanwise spacing, S, for a multiple tab array were identified. Since the goal was to increase mixing, the optimal tab configuration was determined from two properties of the flow field: (i) the large-scale motions with the maximum Ωx, and (ii) the largest number of small-scale motions in a given time period. The peak streamwise vorticity magnitude [mid ]Ωx−max[mid ] was found to have a unique relationship with the tab pitch angle. Furthermore, for all cases examined, the overall small-scale population was found to correlate directly with [mid ]Ωx−max[mid ]. Both quantities peaked at θ≈±45°. It is interesting to note that the peak magnitude of the corresponding circulation in the cross-sectional plane occurred for θ≈±90°. For an array of tabs, the two quantities also depended on the tab spacing. An array of contiguous tabs acted as a solid deflector producing the weakest streamwise vortices and the least small-scale population. For the measurement range covered, the optimal spacing was found to be S≈1.5 tab widths.


Author(s):  
Feng Jie Zheng ◽  
Fu Zheng Qu ◽  
Xue Guan Song

Reservoir-pipe-valve (RPV) systems are widely used in many industrial process. The pressure in an RPV system plays an important role in the safe operation of the system, especially during the sudden operation such as rapid valve opening/closing. To investigate the pressure especially the pressure fluctuation in an RPV system, a multidimensional and multiscale model combining the method of characteristics (MOC) and computational fluid dynamics (CFD) method is proposed. In the model, the reservoir is modeled by a zero-dimensional virtual point, the pipe is modeled by a one-dimensional MOC, and the valve is modeled by a three-dimensional CFD model. An interface model is used to connect the multidimensional and multiscale model. Based on the model, a transient simulation of the turbulent flow in an RPV system is conducted, in which not only the pressure fluctuation in the pipe but also the detailed pressure distribution in the valve are obtained. The results show that the proposed model is in good agreement with the full CFD model in both large-scale and small-scale spaces. Moreover, the proposed model is more computationally efficient than the CFD model, which provides a feasibility in the analysis of complex RPV system within an affordable computational time.


2019 ◽  
Vol 141 (5) ◽  
Author(s):  
Feng Jie Zheng ◽  
Chao Yong Zong ◽  
William Dempster ◽  
Fu Zheng Qu ◽  
Xue Guan Song

Reservoir-pipe-valve (RPV) systems are widely used in many industrial processes. The pressure in an RPV system plays an important role in the safe operation of the system, especially during the sudden operations such as rapid valve opening or closing. To investigate the pressure response, with particular interest in the pressure fluctuations in an RPV system, a multidimensional and multiscale model combining the method of characteristics (MOC) and computational fluid dynamics (CFD) method is proposed. In the model, the reservoir is modeled as a zero-dimensional virtual point, the pipe is modeled as a one-dimensional system using the MOC, and the valve is modeled using a three-dimensional CFD model. An interface model is used to connect the multidimensional and multiscale model. Based on the model, a transient simulation of the turbulent flow in an RPV system is conducted in which not only the pressure fluctuation in the pipe but also the detailed pressure distribution in the valve is obtained. The results show that the proposed model is in good agreement when compared with a high fidelity CFD model used to represent both large-scale and small-scale spaces. As expected, the proposed model is significantly more computationally efficient than the CFD model. This demonstrates the feasibility of analyzing complex RPV systems within an affordable computational time.


Author(s):  
Roger W. Ainsworth ◽  
John L. Allen ◽  
J. Julian M. Batt

The advent of a new generation of transient rotating turbine simulation facilities, where engine values of Reynolds and Mach number are matched simultaneously together with the relevant rotational parameters for dimensional similitude (Dunn et al [1988], Epstein et al [1984]. Ainsworth et al [1988]), has provided the stimulus for developing improved instrumentation for investigating the aerodynamic flows in these stages. Much useful work has been conducted in the past using hot-wire and laser anemometers. However, hot-wire anemometers are prone to breakage in the high pressure flows required for correct Reynolds numbers, Furthermore some laser techniques require a longer runtime than these transient facilites permit, and generally yield velocity information only, giving no data on loss production. Advances in semiconductor aerodynamic probes are beginning to fulfil this perceived need. This paper describes advances made in the design, construction, and testing of two and three dimensional fast response aerodynamic probes, where semiconductor pressure sensors are mounted directly on the surface of the probes, using techniques which have previously been successfully used on the surface of rotor blades (Ainsworth, Dietz and Nunn [1991]). These are to be used to measure Mach number and flow direction in compressible unsteady flow regimes. In the first section, a brief review is made of the sensor and associated technology which has been developed to permit a flexible design of fast response aerodynamic probe. Following this, an extensive programme of testing large scale aerodynamic models of candidate geometries for suitable semiconductor scale probes is described, and the results of these discussed. The conclusions of these experiments, conducted for turbine representative mean and unsteady flows, yielded new information for optimising the design of the small scale semiconductor probes, in terms of probe geometry, sensor placement, and aerodynamic performance. Details are given of a range of wedge and pyramid semiconductor probes constructed, and the procedures used in calibrating and making measurements with them. Differences in performance are discussed, allowing the experimenter to choose an appropriate probe for the particular measurement required. Finally, the application of prototype semiconductor probes in a transient rotor experiment at HP turbine representative conditions is described, and the data so obtained is compared with (PD solutions of the unsteady viscous flow-field.


1962 ◽  
Vol 13 (1) ◽  
pp. 77-81 ◽  
Author(s):  
A. M. Oboukhov

The spectrum of atmospheric turbulence is very broad by comparison with spectra in wind tunnels. We introduce the notion of small-scale and large-scale turbulence. Small-scale turbulence consists of a set of disturbances, the scales of which do not exceed the distance to the wall and for which the hypothesis of three-dimensional isotropy is valid in a certain rough approximation. Large-scale turbulence is essentially anisotropic; the horizontal scale in the atmosphere is much larger than the vertical one, the latter being confined to a certain characteristic height H. The horizontal scale varies widely according to the external conditions and characteristics of the medium.


2010 ◽  
Vol 14 (10) ◽  
pp. 1989-2001 ◽  
Author(s):  
H. Murakami ◽  
X. Chen ◽  
M. S. Hahn ◽  
Y. Liu ◽  
M. L. Rockhold ◽  
...  

Abstract. This study presents a stochastic, three-dimensional characterization of a heterogeneous hydraulic conductivity field within the Hanford 300 Area, Washington, USA, by assimilating large-scale, constant-rate injection test data with small-scale, three-dimensional electromagnetic borehole flowmeter (EBF) measurement data. We first inverted the injection test data to estimate the transmissivity field, using zeroth-order temporal moments of pressure buildup curves. We applied a newly developed Bayesian geostatistical inversion framework, the method of anchored distributions (MAD), to obtain a joint posterior distribution of geostatistical parameters and local log-transmissivities at multiple locations. The unique aspects of MAD that make it suitable for this purpose are its ability to integrate multi-scale, multi-type data within a Bayesian framework and to compute a nonparametric posterior distribution. After we combined the distribution of transmissivities with depth-discrete relative-conductivity profile from the EBF data, we inferred the three-dimensional geostatistical parameters of the log-conductivity field, using the Bayesian model-based geostatistics. Such consistent use of the Bayesian approach throughout the procedure enabled us to systematically incorporate data uncertainty into the final posterior distribution. The method was tested in a synthetic study and validated using the actual data that was not part of the estimation. Results showed broader and skewed posterior distributions of geostatistical parameters except for the mean, which suggests the importance of inferring the entire distribution to quantify the parameter uncertainty.


2011 ◽  
Vol 673 ◽  
pp. 255-285 ◽  
Author(s):  
N. HUTCHINS ◽  
J. P. MONTY ◽  
B. GANAPATHISUBRAMANI ◽  
H. C. H. NG ◽  
I. MARUSIC

An array of surface hot-film shear-stress sensors together with a traversing hot-wire probe is used to identify the conditional structure associated with a large-scale skin-friction event in a high-Reynolds-number turbulent boundary layer. It is found that the large-scale skin-friction events convect at a velocity that is much faster than the local mean in the near-wall region (the convection velocity for large-scale skin-friction fluctuations is found to be close to the local mean at the midpoint of the logarithmic region). Instantaneous shear-stress data indicate the presence of large-scale structures at the wall that are comparable in scale and arrangement to the superstructure events that have been previously observed to populate the logarithmic regions of turbulent boundary layers. Conditional averages of streamwise velocity computed based on a low skin-friction footprint at the wall offer a wider three-dimensional view of the average superstructure event. These events consist of highly elongated forward-leaning low-speed structures, flanked on either side by high-speed events of similar general form. An analysis of small-scale energy associated with these large-scale events reveals that the small-scale velocity fluctuations are attenuated near the wall and upstream of a low skin-friction event, while downstream and above the low skin-friction event, the fluctuations are significantly amplified. In general, it is observed that the attenuation and amplification of the small-scale energy seems to approximately align with large-scale regions of streamwise acceleration and deceleration, respectively. Further conditional averaging based on streamwise skin-friction gradients confirms this observation. A conditioning scheme to detect the presence of meandering large-scale structures is also proposed. The large-scale meandering events are shown to be a possible source of the strong streamwise velocity gradients, and as such play a significant role in modulating the small-scale motions.


2020 ◽  
Author(s):  
Giulia Mazzotti ◽  
Richard Essery ◽  
Johanna Malle ◽  
Clare Webster ◽  
Tobias Jonas

<p>Forest canopies strongly affect snowpack energetics during wintertime. In discontinuous forest stands, spatio-temporal variations in radiative and turbulent fluxes create complex snow distribution and melt patterns, with further impacts on the hydrological regimes and on the land surface properties of seasonally snow-covered forested environments.</p><p>As increasingly detailed canopy structure datasets are becoming available, canopy-induced energy exchange processes can be explicitly represented in high-resolution snow models. We applied the modelling framework FSM2 to obtain spatially distributed simulations of the forest snowpack in subalpine and boreal forest stands at high spatial (2m) and temporal (10min) resolution. Modelled sub-canopy radiative and turbulent fluxes were compared to detailed meteorological data of incoming irradiances, air and snow surface temperatures. These were acquired with novel observational systems, including 1) a motorized cable car setup recording spatially and temporally resolved data along a transect and 2) a handheld setup designed to capture temporal snapshots of 2D spatial distributions across forest discontinuities.</p><p>The combination of high-resolution modelling and multi-dimensional datasets allowed us to assess model performance at the level of individual energy balance components, under various meteorological conditions and across canopy density gradients. We showed which canopy representation strategies within FSM2 best succeeded in reproducing snowpack energy transfer dynamics in discontinuous forests, and derived implications for implementing forest snow processes in coarser-resolution models.</p>


1976 ◽  
Vol 77 (2) ◽  
pp. 321-354 ◽  
Author(s):  
A. Pouquet ◽  
U. Frisch ◽  
J. Léorat

To understand the turbulent generation of large-scale magnetic fields and to advance beyond purely kinematic approaches to the dynamo effect like that introduced by Steenbeck, Krause & Radler (1966)’ a new nonlinear theory is developed for three-dimensional, homogeneous, isotropic, incompressible MHD turbulence with helicity, i.e. not statistically invariant under plane reflexions. For this, techniques introduced for ordinary turbulence in recent years by Kraichnan (1971 a)’ Orszag (1970, 1976) and others are generalized to MHD; in particular we make use of the eddy-damped quasi-normal Markovian approximation. The resulting closed equations for the evolution of the kinetic and magnetic energy and helicity spectra are studied both theoretically and numerically in situations with high Reynolds number and unit magnetic Prandtl number.Interactions between widely separated scales are much more important than for non-magnetic turbulence. Large-scale magnetic energy brings to equipartition small-scale kinetic and magnetic excitation (energy or helicity) by the ‘Alfvén effect’; the small-scale ‘residual’ helicity, which is the difference between a purely kinetic and a purely magnetic helical term, induces growth of large-scale magnetic energy and helicity by the ‘helicity effect’. In the absence of helicity an inertial range occurs with a cascade of energy to small scales; to lowest order it is a −3/2 power law with equipartition of kinetic and magnetic energy spectra as in Kraichnan (1965) but there are −2 corrections (and possibly higher ones) leading to a slight excess of magnetic energy. When kinetic energy is continuously injected, an initial seed of magnetic field will grow to approximate equipartition, at least in the small scales. If in addition kinetic helicity is injected, an inverse cascade of magnetic helicity is obtained leading to the appearance of magnetic energy and helicity in ever-increasing scales (in fact, limited by the size of the system). This inverse cascade, predicted by Frischet al.(1975), results from a competition between the helicity and Alféh effects and yields an inertial range with approximately — 1 and — 2 power laws for magnetic energy and helicity. When kinetic helicity is injected at the scale linjand the rate$\tilde{\epsilon}^V$(per unit mass), the time of build-up of magnetic energy with scaleL[Gt ] linjis$t \approx L(|\tilde{\epsilon}^V|l^2_{\rm inj})^{-1/3}.$


Sign in / Sign up

Export Citation Format

Share Document