Analysis of tracer tomography using temporal moments of tracer breakthrough curves     

Junfeng Zhu  Xing Cai  Tian-Chyi Jim Yeh 《Advances in water resources》2009

Hydraulic/partitioning tracer tomography (HPTT) was recently developed by Yeh and Zhu [Yeh T-CJ, Zhu J. Hydraulic/partitioning tracer tomography for characterization of dense nonaqueous phase liquid source zones, Water Resour Res 2007;43:W06435. doi:10.1029/2006WR004877.] for estimating spatial distribution of dense nonaqueous phase liquids (DNAPLs) in the subsurface. Since discrete tracer concentration data are directly utilized for the estimation of DNAPLs, this approach solves the hyperbolic convection–dispersion equation. Solution to the convection–dispersion equation however demands fine temporal and spatial discretization, resulting in high computational cost for an HPTT analysis. In this work, we use temporal moments of tracer breakthrough curves instead of discrete concentration data to estimate DNAPL distribution. This approach solves time independent partial differential equations of the temporal moments, and therefore avoids solving the convection–dispersion equation using a time marching scheme, resulting in a dramatic reduction of computational cost. To reduce numerical oscillations associated with convection dominated transport problems such as in inter-well tracer tests, the approach uses a finite element solver adopting the streamline upwind Petrov–Galerkin method to calculate moments and sensitivities. We test the temporal moment approach through numerical simulations. Comparing the computational costs between utilizing moments and discrete concentrations, we find that temporal moments significantly reduce the computation time. We also find that tracer moment data collected through a tomographic survey alone are able to yield reasonable estimates of hydraulic conductivity, as indicated by a correlation of 0.588 between estimated and true hydraulic conductivity fields in the synthetic case study.  相似文献   

Stochastic analysis of bounded unsaturated flow in heterogeneous aquifers: Spectral/perturbation approach     

Ching-Min Chang  Hund-Der Yeh 《Advances in water resources》2009

This paper describes a stochastic analysis of steady state flow in a bounded, partially saturated heterogeneous porous medium subject to distributed infiltration. The presence of boundary conditions leads to non-uniformity in the mean unsaturated flow, which in turn causes non-stationarity in the statistics of velocity fields. Motivated by this, our aim is to investigate the impact of boundary conditions on the behavior of field-scale unsaturated flow. Within the framework of spectral theory based on Fourier–Stieltjes representations for the perturbed quantities, the general expressions for the pressure head variance, variance of log unsaturated hydraulic conductivity and variance of the specific discharge are presented in the wave number domain. Closed-form expressions are developed for the simplified case of statistical isotropy of the log hydraulic conductivity field with a constant soil pore-size distribution parameter. These expressions allow us to investigate the impact of the boundary conditions, namely the vertical infiltration from the soil surface and a prescribed pressure head at a certain depth below the soil surface. It is found that the boundary conditions are critical in predicting uncertainty in bounded unsaturated flow. Our analytical expression for the pressure head variance in a one-dimensional, heterogeneous flow domain, developed using a nonstationary spectral representation approach [Li S-G, McLaughlin D. A nonstationary spectral method for solving stochastic groundwater problems: unconditional analysis. Water Resour Res 1991;27(7):1589–605; Li S-G, McLaughlin D. Using the nonstationary spectral method to analyze flow through heterogeneous trending media. Water Resour Res 1995; 31(3):541–51], is precisely equivalent to the published result of Lu et al. [Lu Z, Zhang D. Analytical solutions to steady state unsaturated flow in layered, randomly heterogeneous soils via Kirchhoff transformation. Adv Water Resour 2004;27:775–84].  相似文献   

A Tracer Method to Determine Hydraulic Conductivity and Effective Porosity of Saturated Clays Under Low Gradients   总被引：2，自引：0，他引：2  

Yi-Jang Yeh  Cheng-Haw Lee  Shih-Tsu Chen 《Ground water》2000,38(4):522-529

  相似文献   

Large-scale laboratory experiment on erosion of sand beds by moving circular vertical jets     

Po-Hung Yeh  Kuang-An Chang  John Henriksen  Billy Edge  Peter Chang  Andrew Silver  Abel Vargas 《Ocean Engineering》2009,36(3-4):248-255

This study investigates the topographic deformation due to the erosion of a sand bed impinged by a moving submerged turbulent round jet in a large-scale laboratory. The test conditions represent the case of discharges beneath a vessel while operating in water with a limited clearance such as a shallow navigation channel. The jet moves horizontally and discharges water vertically downward towards the bed. The distance between the jet nozzle and the bed equals six times the jet diameter so the jet flow is in the potential core region. The speed of the jet horizontal motion was varied to examine its effect on the scour profile. The characteristic lengths of the scour profile in the asymptotic state were determined by modifying the empirical formulas in Aderibigbe and Rajaratnam [1996. Erosion of loose beds by submerged circular impinging vertical turbulent jets. Journal of Hydraulic Research 34(1), 19–33]. The maximum scour depth, the scour hole radius, and the ridge height were found to be a function of the ratio of the jet exit to jet translation velocities and were modeled using a hyperbolic function. Empirical equations describing the scour profile were developed and the scour profile was found to be self-similar when normalized by appropriate length scales.  相似文献   

Hydromagnetic buoyancy force in the solar atmosphere     

Tyan Yeh 《Solar physics》1985,95(1):83-97

An extraneous magnetized body, either a flux tube or a plasmoid, immersed in the solar atmosphere is subjected to a hydromagnetic buoyancy force. It results from the peripheral inhomogeneity of ambient hydromagnetic pressure, which is caused or enhanced by the presence of the extraneous body. This extra-caused force acts at various mass elements of the immersed body through its distribution as a nearly uniform force density, just like the gravitational force. Since hydromagnetic buoyancy force comprises hydrostatic buoyancy force, hydrodynamic lift force, and magnetostatic diamagnetic force, this constitutes a magnetohydrodynamic generalization of Archimedes' principle which deals with hydrostatic buoyancy force.In the solar atmosphere hydromagnetic buoyancy force has an obliquely upward direction, with a component in the direction opposite to the downward gravity. It provides an upward force to counterbalance or even to exceed the downward gravitational force. Such an upward force is the dynamic cause for the stationary equilibrium of quiescent prominences and outward motion of coronal transients.  相似文献   

Determination of the parameter pattern and values for a one-dimensional multi-zone unconfined aquifer     

Ya-Chi Chang  Hund-Der Yeh  Yen-Chen Huang 《Hydrogeology Journal》2008,16(2):205-214

In an aquifer, heterogeneity plays an important role in governing groundwater flow. Hence, aquifer characterization should involve both the pattern and values of the hydrogeological parameters. A new analytical solution describing the one-dimensional groundwater flow in a multi-zone unconfined aquifer is presented, and a methodology developed from the analytical solution and a heuristic approach for determining the pattern and values of the aquifer parameters are proposed. The analytical solution demonstrates that the hydraulic head varies spatially and is influenced by aquifer heterogeneity. Simulated annealing, a heuristic approach, is incorporated with the solution to simultaneously identify the pattern and values of the hydraulic conductivity for a horizontal multi-zone unconfined aquifer. This approach may be used to give an approximate result for a two-dimensional problem by dividing the model area into a number of transects along the transverse direction, identifying the parameter values along the longitudinal direction for each transect, and then smoothing the identified results.  相似文献   

Tabulation of Mie properties for an effective microwave radiative model     

H. -Y. M. Yeh  N. Prasad  R. F. Adler 《Meteorology and Atmospheric Physics》1990,42(2):105-112

Summary In microwave radiative transfer model simulations of atmospheric convective systems, the Mie calculations usually consume the majority of the computer time necessary for the calculations (70 to 90% for frequencies ranging from 6 to 300 GHz). For a large array of atmospheric profiles, the repeated Mie scattering calculations make the radiative transfer modeling not only expensive, but often impossible. A set of Mie tables, therefore, is developed to replace the actual Mie calculations in the microwave radiation model.The tables are divided into categories of liquid water (cloud water and rain), ice, snow, graupel and hail. The tables are free from restriction on variation of size distribution, particle density and mixing phases. Results from this study show that by using the Mie tables, the total CPU time is reduced by at least one order of magnitude, depending on the frequency. Compared to using full Mie calculation in the radiative model, the bias and root-mean-square (rms) deviation from the Mie tables are less than 1.2 K for any atmosphere with precipitation rates up to 125 mm hr^–1 at any frequency and any viewing angle.The Mie tables are particularly useful in simulating the microwave radiation field over a mesoscale meteorological system. Such a study would not be satisfactory using a Mie parameterization scheme, because one would have to trade off the accuracy and flexibility for modeling efficiency.With 7 FiguresThis research was done under contract at the Goddard Space Flight Center, Greenbelt, Maryland.  相似文献   

New data and interpretations of the shallow and deep deformation of Huangtupo No. 1 riverside sliding mass during seasonal rainfall and water level fluctuation   总被引：1，自引：0，他引：1  

Jinge Wang  Aijun Su  Wei Xiang  Hsin-Fu Yeh  Chengren Xiong  Zongxing Zou  Cheng Zhong  Qingbing Liu 《Landslides》2016,13(4):795-804

The original Badong County, Hubei, China, was mainly below the highest water level of the Three Gorges Reservoir, which is 175 m above sea level. The new downtown of Badong was rebuilt in the Huangtupo area between 1982 and 1991. After detailed geological investigation in the Huangtupo area, four independent landslides were identified, making it one of the largest and most harmful landslide group in the Three Gorges Reservoir area. Since 2003, abundant data have been obtained from the Huangtupo No. 1 sliding mass about rainfall, water level, earth surface deformation and deep deformation. The monitoring data indicate that the earth surface and deep deformation of this landslide is closely related to the seasonal rainfall and water level fluctuation of the reservoir. During increases in the water level, the earth surface deformation velocity decreases, and then increases obviously in the subsequent water level decreasing stage. Because the water level drawdown period overlaps with the rainy season in this area, the earth surface deformation is affected by both rainfall and water level. The deformation velocity of the earth surface caused by rainfall is about 5 mm/month, while that caused by water level decrease is 5–7 mm/month. On the contrary, the deformation velocity of the deep sliding mass accelerates 2 to 3 times faster than average during water level increase. The distinction of surface and deep deformation regulations indicates that the effects of seasonal rainfall and water level fluctuation on the stability of reservoir wading landslides are different. Based on all monitoring data, we also found that the Huangtupo No. 1 riverside sliding mass is creeping seasonally during the seasonal rainfall and periodic reservoir water level fluctuation. The deformation velocities of the east regions of the sliding body indicate acceleration, making these regions even more dangerous.  相似文献   

GBAS airport availability simulation tool     

Shuo-Ju Yeh  Shau-Shiun Jan 《GPS Solutions》2016,20(2):283-288

A global navigation satellite system augmentation system availability analysis tool has been developed to simulate a ground-based augmentation system (GBAS) prototype, an integrity monitor test bed, to evaluate its operational benefits at an airport of interest. The proposed availability simulation tool includes all GBAS ground facility algorithms as well as a graphical user interface that allows the user to modify simulation options and parameters. The output of the simulation tool is presented in a Stanford chart to help visualize the performance. The chart indicates the availability and integrity. The performance is evaluated primarily in the vertical position domain because of the weaker satellite geometry and more stringent required navigation performance as compared to those of the horizontal position domain. The simulation tool is implemented in Qt (http://www.qt.io/), an open-source cross-platform toolkit, allowing the tool to run on various devices. The computations are performed in the associated C++ code. The Newark Liberty International Airport (ICAO code: KEWR) is used as a simulation example to demonstrate the utility of the developed tool for investigating how reduced error models impact GBAS availability at the airport.  相似文献   

The characteristics of signal versus noise sst variability in the north pacific and the tropical pacific ocean     

Sang-Wook Yeh  Ben P. Kirtman 《Ocean Science Journal》2006,41(1):1-10

Total sea surface temperature (SST) in a coupled GCM is diagnosed by separating the variability into signal variance and noise variance. The signal and the noise is calculated from multi-decadal simulations from the COLA anomaly coupled GCM and the interactive ensemble model by assuming both simulations have a similar signal variance. The interactive ensemble model is a new coupling strategy that is designed to increase signal to noise ratio by using an ensemble of atmospheric realizations coupled to a single ocean model. The procedure for separating the signal and the noise variability presented here does not rely on any ad hoc temporal or spatial filter. Based on these simulations, we find that the signal versus the noise of SST variability in the North Pacific is significantly different from that in the equatorial Pacific. The noise SST variability explains the majority of the total variability in the North Pacific, whereas the signal dominates in the deep tropics. It is also found that the spatial characteristics of the signal and the noise are also distinct in the North Pacific and equatorial Pacific.  相似文献