Solution for Flow Rates across the Wellbore in a Two-Zone Confined Aquifer

A closed-form solution for transient flow rates across the wellbore in a confined aquifer is derived from a two-zone radial ground-water flow equation subject to the boundary condition of keeping a constant head at the well radius. An aquifer may be considered as a two-zone system if the formation properties near the wellbore are significantly changed due to the well construction and/or well development. An efficient numerical approach is used to evaluate this newly derived solution. Values of the transient flow rate are provided in a tabular form and compared with those obtained by numerical inversion for the Laplace-domain solution. The results show that the two solutions are in good agreement. This newly derived solution can be used not only for predicting the transient flow rate across the wellbore but also for identifying the effects of a skin with a finite thickness on the estimation of transient flow rates in a ground-water system with two different formation properties.     

Analytical Solution for Two-Dimensional Solute Transport in Finite Aquifer with Time-Dependent Source Concentration

Using the Hankel Transform Technique, an analytical solution is derived for two-dimensional solute transport in a homogeneous isotropic aquifer. The aquifer is subjected to time-dependent point source contamination. The solution is derived under two conditions: (1) the flow velocity in the aquifer is a sinusoidally varying function and (2) the flow velocity is an exponentially decreasing function. Initially the aquifer is assumed solute free. The analytical solution is illustrated using an example.     

A mathematical model for marangoni flow and mass transfer in electrostatically positioned droplets

A mathematical model is developed for the free surface deformation, full three-dimensional (3-D) Marangoni convection and solute transport phenomena in electrostatically levitated droplets under microgravity. The electric field is calculated by the boundary element method and the shape deformation by the weighted residuals method. The numerical model for the transport phenomena is developed based on the Galerkin finite-element solution of the Navier-Stokes equations, the energy balance equation, and the mass transport equation. Numerical simulations are carried out for droplet deformation by electrostatic forces and 3-D Marangoni convection in droplets heated by three different heating source arrangements. Results show that the electric forces deform a droplet into an oval shape under microgravity by pulling the droplet apart at the two poles. A two-beam heating arrangement results in an axisymmetric flow and temperature distribution in the droplet. Complex 3-D Marangoni flow structure occurs when a tetrahedral or octahedral heating arrangement is applied. The thermal transport in the droplet is conduction dominant for the cases studied. In general, the convection is stronger with higher melting point melts. The internal convection has a strong effect on the concentration distribution in the droplet. For melts with high viscosities, a significant reduction in velocity can be achieved with an appropriate laser beam arrangement, thereby permitting a diffusion-controlled condition to be developed.     

Simplified Numerical and Analytical Approach for Solutes in Turbulent Flow Reacting with Smooth Pipe Walls

A large group of reactions that affect water quality in distribution networks occur on the pipe wall surface. Existing simulation models are usually based on cross-sectionally averaged variables that use mass-transfer coefficients derived for constant-concentration (Dirichlet) boundary conditions to account for cross-sectional variations. In the case of a first-order wall-demand problem, the boundary condition is however of Robin type. We derive a simple one-dimensional (1D) model for the radial concentration profile of a solute of arbitrary Schmidt number (Sc) reacting with pipe walls in a fully developed turbulent flow. A modified van Driest mixing length model was used to approximate the Reynolds-averaged velocity and eddy diffusivity. Numerical solutions of the 1D model agree well with a two-dimensional mass transport model and experimental data. An asymptotic solution for high Sc is derived, which is in excellent agreement with the 1D model for Sc>100. A comparison with the mass-transfer coefficients for constant-concentration boundary conditions shows that the differences between the two boundary conditions are small.     

Coupled Surface–Subsurface Solute Transport Model for Irrigation Borders and Basins. I. Model Development

Surface fertigation is widely practiced in irrigated crop production systems. Lack of design and management tools limits the effectiveness of surface fertigation practices. The availability of a process-based coupled surface–subsurface hydraulic and solute transport model can lead to improved surface fertigation management. This paper presents the development of a coupled surface–subsurface solute transport model. A hydraulic model described in a previous paper by the writers provided the hydrodynamic basis for the solute transport model presented here. A numerical solution of the area averaged advection–dispersion equation, based on the split-operator approach, forms the surface solute transport component of the coupled model. The subsurface transport process is simulated using HYDRUS-1D, which also solves the one-dimensional advection–dispersion equation. A driver program is used for the internal coupling of the surface and subsurface transport models. Solute fluxes calculated using the surface transport model are used as upper boundary conditions for the subsurface model. Evaluation of the model is presented in a companion paper.     

Centrifuge Modeling of Unstable Infiltration and Solute Transport

Unstable flow can result in the formation of fingers during infiltration into dry soil. Centrifuge modeling is a potentially useful tool to study the relationship between finger size, spacing, and velocity. It can also be used to investigate solute transport in such fingers. Physical properties of the fingers are obtained for three tests conducted at elevated acceleration levels. A fourth test was conducted at 1g. The physical parameters compare well with theoretical predictions. To assess solute transport in fingers, a known concentration of solute was introduced after the fingers had formed. The resulting breakthrough curves were analyzed using the two-region model as well as the advection dispersion equation with appropriate initial and boundary conditions. Although the two-region model is physically more plausible, it was found to match the extensive tailing observed in the breakthrough curves only marginally better than the advection-dispersion equation.     

Diffusion-controlled growth of disordered interphase boundaries in finite matrix

The diffusion-controlled growth of interphase boundaries in finite matrices, which is influenced by the enrichment or depletion of solute in the untransformed matrix, is investigated using the solution to the diffusion equation obtained from the general formulation of free boundary problems by Kolodner. The integro-differential equation for the boundary position and velocity derived from the flux balance of solute is numerically solved for both one-dimensional (planar growth) and three-dimensional cases (growth of a spherical interface toward its center) and compared with the quasi-stationary solution in corresponding situations. The results show that the latter solution coupled with the mass balance of solute between precipitate and matrix as previously utilized by some authors tends to overestimate the progress of transformation. The error becomes quite large at higher supersaturations where the quasi-stationary approximation obviously breaks down.     

Laplace-Domain Solutions for Radial Two-Zone Flow Equations under the Conditions of Constant-Head and Partially Penetrating Well

A mathematical model is presented for a constant-head test performed in a partially penetrating well with a finite-thickness skin. The model uses a no-flow boundary condition for the casing and a constant-head boundary condition for the screen to represent the partially penetrating well. The Laplace-domain solutions for the dimensionless flow rate at the wellbore and the hydraulic heads in the skin and formation zones are derived using the Laplace and finite Fourier cosine transforms. The solutions of hydraulic heads have been shown to satisfy the governing equations, related boundary conditions, and continuity requirements for the pressure head and flow rate at the interface of the skin zone and undisturbed formation. In addition, an efficient algorithm for evaluating those solutions is also presented. The dimensionless flow rates obtained from new solutions have been shown to be better than those of Novakowski’s solutions, especially when the penetration ratio is large.     

Extended Thermodynamics Derivation of Energy Dissipation in Unsteady Pipe Flow

Extended irreversible thermodynamics (EIT) provides a framework for deriving extensions to phenomenological equations (e.g., Newton's law of viscosity, Fick's law of mass transport, and Darcy's law for porous media flow) for problems involving high frequencies (i.e., rapid transients). In this paper, a phenomenological equation is derived for energy loss in 1D unsteady pipe flow using an EIT formalism. The resulting wall shear stress is equal to the sum of (1) the steady-state shear stress; (2) a term that is proportional to the local (i.e., temporal) acceleration; and (3) a term that is proportional to the product of the velocity and the convective (i.e., spatial) acceleration. The form of this EIT-based wall shear stress formula shows that EIT provides a physical basis for instantaneous acceleration based unsteady friction formulas. It also illustrates the limitations and underlying assumptions of these models. For example, instantaneous acceleration based unsteady friction formulas are limited to fast transients (i.e., transients in which the water hammer timescale is significantly smaller than the diffusion timescale). A characteristics solution for unsteady pipe flow is proposed in which the phenomenological equation is used to model energy dissipation. Comparison of numerical test results with measured data from upstream and downstream valve closure laboratory experiments shows excellent agreement.     

Flow Depletion Induced by Pumping Well from Stream Perpendicularly Intersecting Impermeable/Recharge Boundary

Analytical solutions for rate and volume of flow depletion induced by pumping a well from a stream that intersects an impermeable or a recharge boundary at right angles are derived using the basic flow depletion factor defined earlier by the author. A new concept of directly obtaining stream flow depletion using the method of images is proposed. The solutions are derived for five different management cases of a stream and boundary intersecting at right-angles, assuming the aquifer to be confined with semi-infinite areal extent. A computationally simple function is proposed for accurately approximating the error function. The existing analytical solution in the case of a right-angle bend of stream given by Hantush was obtained for unconfined aquifers using a linearization of the governing partial differential equation. The solution for this case obtained using the proposed method for confined aquifer is the same as obtained by Hantush for unconfined aquifers, which shows that the linearization adopted by Hantush does not actually solve this problem for unconfined aquifers.     

Coupled Large Strain Consolidation and Solute Transport. II: Model Verification and Simulation Results

The results of numerical simulations for coupled large strain consolidation and solute transport, obtained using the CST1 model, are presented. CST1 accounts for advection, longitudinal and transverse dispersion, first-order decay reactions, and linear equilibrium sorption. Verification checks of CST1 show excellent agreement with analytical solutions for one-dimensional (1D) transport in rigid porous media, including various Peclet numbers and concentration boundary conditions. Similarly excellent agreement is observed for two-dimensional advection-dispersion transport in rigid media and 1D advection-dispersion transport in compressible media undergoing large strain consolidation. CST1 is then used to investigate consolidation-induced solute transport for a single composite liner system and a confined disposal facility for dredged contaminated sediments. In both cases, solute transport was found to be strongly affected by consolidation-induced advection both during and after the consolidation period. Consolidation has a lasting effect on solute migration because transient advective flows change the distribution of solute mass, which then becomes the initial condition for subsequent transport processes.     

Analytical Solutions for Contaminant Diffusion in Double-Layered Porous Media

Analytical solutions for conservative solute diffusion in one-dimensional double-layered porous media are presented in this paper. These solutions are applicable to various combinations of fixed solute concentration and zero-flux boundary conditions (BC) applied at each end of a finite one-dimensional domain and can consider arbitrary initial solute concentration distributions throughout the media. Several analytical solutions based on several initial and BCs are presented based on typical contaminant transport problems found in geoenvironmental engineering including (1) leachate diffusion in a compacted clay liner (CCL) and an underlying stratum; (2) contaminant removal from soil layers; and (3) contaminant diffusion in a capping layer and underlying contaminated sediments. The analytical solutions are verified against numerical solutions from a finite-element method based model. Problems related to leachate transport in a CCL and an underlying stratum of a landfill and contaminant transport through a capping layer over contaminated sediments are then investigated, and the suitable definition of the average degree of diffusion is considered.     

Model for Consolidation-Induced Solute Transport with Nonlinear and Nonequilibrium Sorption

A numerical model, called CST2, is presented for coupled large strain consolidation and solute transport in saturated porous media. The consolidation and solute transport algorithms include the capabilities of a previous model, CST1, with the addition of a variable effective diffusion coefficient during consolidation and nonlinear nonequilibrium sorption. The model is based on a dual-Lagrangian framework that tracks separately the motions of fluid and solid phases. Verification checks of CST2 show excellent agreement with analytical and numerical solutions for solute transport in rigid porous media. A parametric study illustrates that, for the test cases considered, variation of effective diffusion coefficient during consolidation has an important effect on solute transport, whereas nonlinearity of the sorption isotherm has a less important effect. Additional simulations show that nonequilibrium sorption can have a strong effect on consolidation-induced solute transport and that this effect becomes more important as the rate of consolidation increases. The simulations also corroborate previous findings that consolidation can have a lasting effect on solute migration because transient advective flows change the distribution of solute mass which then becomes the initial condition for subsequent transport processes.     

Turbulent Boundary Layer Flows Above a Porous Surface Subject to Flow Injection

A new analytical expression for velocity profile in a fully developed turbulent boundary layer above a porous surface subject to flow injection is derived by solving the coupled Reynolds equations and turbulent kinetic energy equation. The advection of turbulent kinetic energy is considered during the derivation, whereas the earlier studies have neglected it. The new solution reduces to the universal logarithmic law in the case of no flow injection. For the small injection, the solution can be expanded into a series form in terms of the normalized injection velocity. The leading order terms are found to be equivalent to those in the earlier works in which the advection of turbulent kinetic energy has been neglected in the derivation. The new solution can provide more accurate prediction of bed shear stress for a wide range of flow injection rate, fluid type (e.g., from air to water), and surface roughness. On the other hand, the earlier theories may significantly underestimate bed shear stress under high injection rates.     

A statistical mechanical theory for activity coefficients of a dilute solute in a binary solvent

A statistical mechanical calculation of the activity coefficients of a dilute solute,C, in a binary solvent,A-B, was made using a model for the interactions of a coordination cluster consisting of a solute atom and its neighboring solvent atoms. The derived equations are applicable to a variety of types of ternary solutions including substitutional and interstitial alloys as well as additive and reciprocal molten salts. The theory takes into account the interactions between solute and solvent atoms (ions) as well as changes in interactions of the solvent atoms (ions) which are neighbors of solute atoms (ions). Prior theories such as those of Wagner and the quasi-chemical theories of Alcock and Richardson and Jacob and Alcock can all be shown to be special cases of the present theory. The dependence of the activity coefficients of a solute on the solvent composition is calculated from a knowledge of the activity coefficients of the solvent components, solute activity coefficients in the two pure solvent components, a coordination number, a geometric factor which depends upon the type of solution, and a term which represents the nonadditivity of pair bond interactions within the cluster of a solute atom (ion) and its neighboring solvent atoms (ions). In the model, the thermodynamic properties of the solute are related to the relative concentrations of the different coordination clusters as well as to the thermodynamic properties of the solvent.     

Scanning gel chromatography: determination of transport parameters from dynamic profiles measured at low solute concentration

Direct optical scanning of solute boundaries in large zone gel chromatography experiments provides an accurate means of determining boundary profile shapes and rates of motion. A method has been developed for correcting such boundaries to a constant time frame, eliminating the distortion which arises from finite column scanning rate. Centroids of the corrected profiles can be used to determine the partition cross section for the solute of interest. The partition cross section and flow rate determine translational motion within the column. The axial dispersion coefficient, L, which characterizes rate of boundary spreading may also be calculated from the profiles. In order to explore these procedures a study of four noninteracting solutes was conducted. Partition cross sections determined from rates of motion of boundary centroids were found to be in good agreement with those determined by the equilibrium saturation method on the same column. In order to explore the lowest concentration limits of the technique and to illustrate the boundary characteristics for a self-associating solute, a study of carboxyhemoglobin was conducted over a wide concentration range. From measurements at 220 nm the lowest concentration where useful data could be obtained was 2 micrograms per ml (0.12muM heme). These results establish validity of the procedures used in analyzing the rates of boundary transport and in studying solute transport over a wide range of conditions.     

Three-Dimensional Mathematical Model of Suspended-Sediment Transport

This paper presents the basic equations for a mathematical model of sediment-laden flow in a nonorthogonal curvilinear coordinate system. The equations were derived using a tensor analysis of two-phase flow and incorporate a natural variable-density turbulence model with nonequilibrium sediment transport. Correspondingly, a free-surface and the bottom sediment concentration are employed to provide the boundary conditions at the river surface and the riverbed. The finite analytic method is used to solve the equations of mass and momentum conservation and also the transport equation for suspended sediment. To demonstrate the method, the sediment deposition for the Three Gorges Project is considered. The mathematical model specifies the boundary conditions for the inlet and outlet using data from physical model experiments. The results for the mathematical model were tested against laboratory measurements from the physical model experiment. Good agreement and accuracy were obtained.     

Diffusive Modeling of Aggradation and Degradation in Artificial Channels

The unsteady flow and solid transport simulation problem in artificial channels is solved using a three-equation model, coupled with a local erosion law. The three equations are the water mass and momentum balance equations, as well as the total solid load balance equation. It is shown that even during severe hydrological events inertial terms can be neglected in the momentum equation without any substantial change in the solution sought. Empirical equilibrium formulas were used to estimate the solid load as a function of the flow variables. Local erosion, due to the scour generated at the jump between two channels connected at different bottom elevations, was estimated adapting a literature formulation. The double order approximation time and space marching scheme, previously proposed for the solution of the unsteady flow problem in the fixed-bed case, is applied to the solution of the new system. The model was validated with both literature and new laboratory experimental data. No parameter calibration was used to fit the computed results to the experimental ones.     

Dynamic Instability of Nanorods/Nanotubes Subjected to an End Follower Force

This paper presents an investigation on the dynamic instability of cantilevered nanorods/nanotubes subjected to an end follower force. Eringen’s nonlocal elasticity theory is employed to allow for the small length scale effect in the considered dynamic instability problem. The general solution for the governing differential equation is obtained and the dynamic instability characteristic equation is derived by applying the boundary conditions. Exact critical load factors are obtained. These nonlocal solutions are compared with the classical local solutions to assess the sensitivity of the small length scale effect on the critical load factors and flutter mode shapes. It is found that the small length scale effect decreases the critical load and the corresponding frequency parameters as well as reduces the severity of the double-curvature flutter mode shape.