Transient‐state groundwater flow in various geometries with wells at arbitrary positions: analytical methods

We have developed a method for analytically solving the porous medium flow equation in many different geometries for horizontal (two‐dimensional), homogeneous and isotropic aquifers containing impermeable boundaries and any number of pumping or injection wells located at arbitrary positions within the system. Solutions and results are presented for rectangular and circular aquifers but the method presented here is easily extendible to many geometries. Results are also presented for systems where constant head boundary conditions can be emulated internal to the aquifer boundary. Recommendations for extensions of the present work are briefly discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

Multiple bipolar structure of normal magnetic field at the magnetopause

We describe how a local plasma structure can be changed by a transverse shear flow using numerical simulation to investigate the disturbance process near the magnetopause. The results show that magnetic field lines are bent by transverse shear flow disturbance near the current sheet region. There are multiple bipolar structures of the normal magnetic field in the numerical simulation. We term this new feature as K-point magnetic reconnection, realistic for discussing space observations.     

Circulation in the Gulf of Aqaba (Red Sea) during winter—spring

This study is devoted to oceanographic features of the semi-enclosed Gulf of Aqaba, Red Sea. The data were recorded in winter—spring 1999 on the R/V Meteor cruise leg 44/2. Temperature and salinity profiles were measured at six positions (I—VI). The shipboard NarrowBand Acoustic Doppler Current Profiler (NB ADCP) 150 kHz continuously recorded current profiles down to 350 m en route. The research revealed that the current near the Strait of Tiran front (position VI) represents a semidiurnal signal of an internal tide wave (∼12 h period; 0.2 ms⁻¹ amplitude) that might be generated by the barotropic tide at the sill of the Strait. A sequence of cyclonic and anticyclonic eddy pairs is found along the axis of the Gulf of Aqaba during winter-spring seasons. These sub-mesoscale signals are dominant above the main thermocline and might be caused by wind forcing and the narrowness of the Gulf; it might remain in other seasons with different dimensions in relation to the depth of thermocline. The total diameter of each pair was twice the baroclinic Rossby radius (R ≈ 10 km). A single anti-cyclonic eddy was observed in the upper 300 m in the northern tip of the Gulf with a diameter of about 5–8 km.     

CO<Subscript>2</Subscript> Flux Measurements over the Sea Surface by Eddy Correlation and Aerodynamic Techniques

The measurements of the vertical transport of CO₂ were carried out over the Sea of Japan using the specially designed pier of Kyoto University on September 20 to 22, 2000. CO₂ fluxes were measured by the eddy correlation and aerodynamic techniques. Both techniques showed comparable CO₂ fluxes during sea breeze conditions: −0.001 to −0.08 mg m⁻²s⁻¹ with the mean of −0.05 mg m⁻²s⁻¹. This means that the measuring site satisfies the fetch requirement for meteorological observations under sea breeze conditions. Moreover, the eddy diffusivity coefficient used in the aerodynamic technique is found to be consistent with the coefficient used in the eddy correlation technique. The present result leads us to conclude that the aerodynamic technique may be applicable to underway CO₂ flux measurements over the ocean and may be used in place of the bulk technique. The important point is the need to maintain a measuring accuracy of CO₂ concentration difference of the order of 0.1 ppmv on the research vessels or the buoys.     

浅海流体动力学一种含变湍粘性系数的流速分解模型(Ⅱ)——模型在超浅海风暴潮零阶问题中的应用

本文将流速分解模型应用于作为超浅海风暴潮的渤海风潮,并讨论了变湍粘性系数的确定。作为一个初步的,但较为成功的数值试验例子,描述了实际风场作用下的渤海风潮,比较了变湍粘性系数模型与常湍粘性系数模型的计算结果间的差异。     

A model of the formation and development of ocean shear wave

Ａｍｏｄｅｌｏｆｔｈｅｆｏｒｍａｔｉｏｎａｎｄｄｅｖｅｌｏｐｍｅｎｔｏｆｏｃｅａｎｓｈｅａｒｗａｖｅ￥ＱｉａｏＦａｎｇｌｉａｎｄＺｈａｎｇＱｉｎｇｈｕａ（ＦｉｒｓｔＩｎｓｔｉｔｕｔｅｏｆＯｃｅａｎｏｇｒａｐｈｙ，ｓｔａｔｅｏｃｅａｎｉｃＡｄｍｉｎｉｓ...     

Monitoring the Kuroshio Extension with Dynamically Constrained Synthesis of the Acoustic Tomography, Satellite Altimeter and in situ Data

A finite-difference quasigeostrophic (QG) model of an open ocean region has been employed to produce a dynamically constrained synthesis of acoustic tomography and satellite altimetry data with in situ observations. The assimilation algorithm is based upon the 4D variational data interpolation scheme controlled by the model's initial and boundary conditions. The data sets analyzed include direct and differential travel times measured at the array of five acoustic transceivers deployed by JAMSTEC in the region of the Kuroshio Extension in 1997, Topex/Poseidon altimetry, CTD soundings, and ADCP velocity profiles. The region monitored is located within the area 27.5°–36.5°N, 143°–155°. The results of assimilation show that mesoscale variability can be effectively reconstructed by five transceivers measuring direct and reciprocal travel times supported by relatively sparse in situ measurements. The misfits between model and data lie within the observational error bars for all the data types used in assimilation. We have compared the results of assimilation with the statistical inversion of travel time data and analyzed energy balances of the optimized model solution. Energy exchange between the depth-averaged and shear components of the observed currents reveals a weak decay of the barotropic mode at the rate of 0.2 ± 0.7⋅10⁻⁵ cm²/s³ due to topographic interaction. Mean currents in the region are unstable with an estimate of the available potential energy flux from the mean current to the eddies of 4.7 ± 2.3⋅10⁻⁵ cm²/s³. Kinetic energy transition has the same sign and is estimated as 2.8 ± 2.5⋅10⁻⁵ cm²/s³. Potential enstrophy is transferred to the mesoscale at a rate of 5.5 ± 2.7⋅10⁻¹⁸ s⁻³. These figures provide observational evidence of the properties of free geostrophic turbulence which were predicted by theory and observed in numerical experiments. This revised version was published online in July 2006 with corrections to the Cover Date.     

Numerical modeling of long waves in shallow water using Incremental Differential Quadrature Method

Incremental Differential Quadrature Method (IDQM) as a rapid and accurate method for numerical simulation of Nonlinear Shallow Water (NLSW) waves is employed. To the best of authors’ knowledge, this is the first endeavor to exploit DQM in coastal hydraulics. The one-dimensional NLSW equations and related boundary conditions are discretized in space and temporal directions by DQM rules and the resulting system of equations are used to compute the state variables in the entire computational domain. It was found that the splitting of total simulation time into a number of smaller time increments, could significantly enhance the performance of the proposed method. Furthermore, results of this study show two main advantages for IDQM compared with other conventional methods, namely; unconditional stability and minimal computational effort. Indeed, using IDQM, one can use a few grid points (in spatial or time direction) without imposing any stability condition on the time step to obtain an accurate convergent solution.     

Large eddy simulation of hydrocyclone—prediction of air-core diameter and shape

Hydrocyclones are widely used in the mining and chemical industries. An attempt has been made in this study, to develop a CFD (computational fluid dynamics) model, which is capable of predicting the flow patterns inside the hydrocyclone, including accurate prediction of flow split as well as the size of the air-core. The flow velocities and air-core diameters are predicted by DRSM (differential Reynolds stress model) and LES (large eddy simulations) models were compared to experimental results. The predicted water splits and air-core diameter with LES and RSM turbulence models along with VOF (volume of fluid) model for the air phase, through the outlets for various inlet pressures were also analyzed. The LES turbulence model led to an improved turbulence field prediction and thereby to more accurate prediction of pressure and velocity fields. This improvement was distinctive for the axial profile of pressure, indicating that air-core development is principally a transport effect rather than a pressure effect.     

Laplace-transform analytic element solution of transient flow in porous media

A Laplace-transform analytic element method (LT-AEM) is described for the solution of transient flow problems in porous media. Following Laplace transformation of the original flow problem, the analytic element method (AEM) is used to solve the resultant time-independent modified Helmholtz equation, and the solution is inverted numerically back into the time domain. The solution is entirely general, retaining the mathematical elegance and computational efficiency of the AEM while being amenable to parallel computation. It is especially well suited for problems in which a solution is required at a limited number of points in space–time, and for problems involving materials with sharply contrasting hydraulic properties. We illustrate the LT-AEM on transient flow through a uniform confined aquifer with a circular inclusion of contrasting hydraulic conductivity and specific storage. Our results compare well with published analytical solutions in the special case of radial flow.