A multi-moment finite volume formulation for shallow water equations on unstructured mesh

A novel and accurate finite volume method has been presented to solve the shallow water equations on unstructured grid in plane geometry. In addition to the volume integrated average (VIA moment) for each mesh cell, the point values (PV moment) defined on cell boundary are also treated as the model variables. The volume integrated average is updated via a finite volume formulation, and thus is numerically conserved, while the point value is computed by a point-wise Riemann solver. The cell-wise local interpolation reconstruction is built based on both the VIA and the PV moments, which results in a scheme of almost third order accuracy. Efforts have also been made to formulate the source term of the bottom topography in a way to balance the numerical flux function to satisfy the so-called C-property. The proposed numerical model is validated by numerical tests in comparison with other methods reported in the literature. (C) 2010 Elsevier Inc. All rights reserved.

A global shallow water model using high order multi-moment constrained finite volume method and icosahedral grid

A novel accurate numerical model for shallow water equations on sphere have been developed by implementing the high order multi-moment constrained finite volume (MCV) method on the icosahedral geodesic grid. High order reconstructions are conducted cell-wisely by making use of the point values as the unknowns distributed within each triangular cell element. The time evolution equations to update the unknowns are derived from a set of constrained conditions for two types of moments, i.e. the point values on the cell boundary edges and the cell-integrated average. The numerical conservation is rigorously guaranteed. in the present model, all unknowns or computational variables are point values and no numerical quadrature is involved, which particularly benefits the computational accuracy and efficiency in handling the spherical geometry, such as coordinate transformation and curved surface. Numerical formulations of third and fourth order accuracy are presented in detail. The proposed numerical model has been validated by widely used benchmark tests and competitive results are obtained. The present numerical framework provides a promising and practical base for further development of atmospheric and oceanic general circulation models. (C) 2009 Elsevier Inc. All rights reserved.

Effect of volume ratio on thermocapillary flow in liquid bridges of high-Prandtl-number fluids

In present study, the transition of thermocapillary convection from the axisymmetric stationary flow to oscillatory flow in liquid bridges of 5cst silicon oil (aspect ratio 1.0 and 1.6) is investigated in microgravity conditions by the linear instability analysis. The corresponding marginal instability boundary is closely related to the gas/liquid configuration of the liquid bridge noted as volume ratio. With the increasing volume ratio, the marginal instability boundary consists of the increasing branch and the decreasing branch. A gap region exists between the branches where the critical Marangoni number of the corresponding axisymmetric stationary flow increases drastically. Particularly, a unique axisymmetric oscillatory flow (the critical azimuthal wave number is m=0) in the gap region is reported for the liquid bridge of aspect ratio 1.6. Moreover, the energy transfer between the basic state and the disturbance fields of the thermocapillary convection is analyzed at the corresponding critical Marangoni number, which reveals different major sources of the energy transfer for the development of the disturbances in regimes of the increasing branch, the gap region and the decreasing branch, respectively.

A study of the influence of initial liquid volume on the capillary flow in an interior corner under microgravity

In this work the influence of initial liquid volume on the capillary flow in an interior corner is studied systematically by microgravity experiments using the drop tower, under three different conditions: the Concus-Finn condition is satisfied,close to and dissatisfied. The capillary flow is studied by discussing the movement of tip of the meniscus in the corner. Experimental results show that with the increase of initial liquid volume the tip location increases for a given microgravity time, the achievable maximum tip velocity increases and the flow reaches its maximum tip velocity earlier However, the results for the three different conditions show some difference. (C) 2010 Elsevier Ltd All rights reserved

Spatial structure of internal tidal waves generated on weak topographies

The generation of internal gravity waves by barotropic tidal flow passing over a two-dimensional topography is investigated. Rather than calculating the conversion of tidal energy, this study focuses on delineating the geometric characteristics of the spatial structure of the resulting internal wave fields (i.e., the configurations of the internal beams and their horizontal projections) which have usually been ignored. it is found that the various possible wave types can be demarcated by three characteristic frequencies: the tidal frequency, wo; the buoyancy frequency, N; and the vertical component of the Coriolis vector or earth's rotation.f. When different possibilities arising from the sequence of these frequencies are considered, there occur 12 kinds of wave structures in the full 3D space in contrast to the 5 kinds identified by the 2D theory. The constant wave phase lines may form as ellipses or hyperbolic lines on the horizontal plane, provided the buoyancy frequency is greater or less than the tidal frequency. The effect that stems from the consideration of the basic flow is also found, which not only serves as the reason for the occurrence of higtter harmonics but also increases the wave strength in the direction of basic flow. (C) 2009 Elsevier B.V. All rights reserved.

Determination of arsenic species by capillary zone electrophoresis with large-volume field-amplified stacking injection

Determination of arsenic species by large-volume field amplified stacking injection-capillary zone electrophoresis (LV-FASI-CZE) is reported in this paper. Whole column injection was employed. The optimum buffer pH for the separation of weak acids was discussed. It was found that the optimum buffer to analyze the stacked arsenate (As(V)), monomethylarsonate (MMA), and dimethylarsinate (DMA) was 25 mm phosphate at pH 6.5. However, the optimum buffer to analyze the concentrated arsenite (As(III)) was 20 mm phosphate - 10 mm borate at pH 9.28. The limits of detection of the method developed were 0.026 mg/L for As(III), 0.023 mg/L for As(V), 0.043 mg/L for MMA, and 0.018 mg/L for DMA. An enrichment factor of 34-100 for several arsenic species was obtained. In the end, this method was applied to determine the arsenic concentration in the environmental reference materials to show the usefulness of the method developed.

Stability shear stress and equilibrium cross-sectional geometry of sheltered tidal channels

This study relates tidal channel cross-sectional area (A) to peak spring discharge (Q) via a physical mechanism, namely the stability shear stress ( tau sub(S)) just necessary to maintain a zero gradient in net along-channel sediment transport. It is assumed that if bed shear stress ( tau ) is greater than tau sub(S), net erosion will occur, increasing A, and reducing tau similar to (Q/A) super(2) back toward tau sub(S). If tau < tau sub(S) there will be net deposition, reducing A and increasing tau toward tau sub(S). A survey of the literature allows estimates of Q and A at 242 sections in 26 separate sheltered tidal systems. Assuming a single value of tau sub(S) characterizes the entire length of a given tidal channel, it is predicted that along-channel geometry will follow the relation Ah sub(R) super(1) super(/) super(6) similar to Q. Along-channel regressions of the form Ah sub(R) super(1) super(/) super(6) similar to Q super( beta ) give a mean observed value for beta of 1.00 plus or minus 0.06, which is consistent with this concept. Results indicate that a lower bound on tau sub(S) (and an upper bound on A) for stable channels is provided by the critical shear stress ( tau sub(C)) just capable of initiating sediment motion. Observed tau sub(S) is found to vary among all systems as a function of spring tidal range (R sub(sp)) according to the relation tau sub(S) approximately 2.3 R sub(sp) super(0.79) tau sub(C). Observed deviations from uniform tau sub(S) along individual channels are associated with along-channel variation in the direction of maximum discharge (i.e., flood-versus ebb-dominance).

the ebb-tidal delta of the venice lagoon,Italy

Very high-resolution seismic lines were collected offshore the Venice Lagoon in the framework of the Co.Ri.La'New very high resolution seismic methods to study the Venice Lagoon subsoil' project. The 140 km boomer profiles led to the identification of a convex-upward,lens-shaped body just outside the Lido inlet,one of the three lagoon inlets,which has been interpreted as an ebb-tidal delta.The comparison between the scismostratigraphic setting og this deposit, as revealed by the seismic lines,and the analyses of historicalbathymetric maps highlighted the key role of human interventions in the formation and evolution of the Lido inlet cbb-tidal delta.To preserve the lagoon environment but also to ensure a navigable way, human interventions at the Lidoinlet,performed since the fourteenth centtury, caused profound variations in the inlet dynamics,ldading to a progressive increase in the sediment dispersion from the lagoon interior towards the sea. The ebb-tidal delta of the Lido nilet is thus a very recent feature compared with the formation of the Venice lagoon and formed mainly as a consequence of the construction of the two jetties that have bound the inlet form AD 1886 to the present day.

Spontaneous volume transition of polyampholyte nanocomposite hydrogels based on pure electrostatic interaction

A circular system is employed in this paper to investigate the swelling behaviors of polyampholyte hydrogels; this circular system can effectively eliminate the disturbance of various factors and keep the surrounding environment constant. It is found that there exists a spontaneous volume transition to the collapsed state of polyampholyte hydrogels, which is attributed to the overshooting effect, and the transition can occur repeatedly under certain conditions. C-13 NMR is employed to investigate the swelling behavior of polyampholyte hydrogels.

The time delay in electrochemical measurements of a finite-volume system

For a sphere electrode enclosed in finite-volume electrolyte, the measured current will deviate from the result predicted by the semi-infinite diffusion theory after some time. By random-walk simulation, we compared this time to the one needed for diffusion layer to reach electrolyte boundary, and revealed a clear signal delay of electrochemical current. Further we presented a quantitative description of this delay time. The simulation results suggested that the semi-infinite diffusion theory can even be applied when the theoretical diffusion layer grows to 1.28 electrolyte thicknesses, with an accuracy better than 0.5%. We attributed this time delay to the molecules' finite propagation velocity. Finally, we discussed how this delay can influence and facilitate the following electrochemical detection towards the nanometer and single-cell scale.

Studies of effect of phase volume ratio on transfer of ionizable species across the water/1,2-dichloroethane interface by a three-electrode setup

The electrochemical behavior of pyridine distribution at the water/1,2-dichloroethane interface with variable phase volume ratios (r=V-0/V-W) was investigated by cyclic voltammetry. The system was composed of an aqueous droplet supported on a Ag/AgCl disk electrode covered with an organic solution or an organic droplet supported on a Ag/AgTPBCl disk electrode covered with an aqueous solution. In this way, a conventional three-electrode potentiostat can be used to study an ionizable compound transfer process at a liquid/liquid interface with a wide range of phase volume ratios (from 0.0004 to 1 and from 1 to 2500). Using this special cell we designed, only very small volumes of both phase were needed for r equal to unity, which is very useful for the investigation of the distribution of ionizable species at a biphasic system when the available amount of species is limited. The ionic partition diagrams were obtained for different phase volume ratios.

Study of ionizable drugs transfer across the water/1,2-dichloroethane interface with phase volume ratio equal to unity using a three-electrode system

The electrochemical behavior of ionizable drugs (Amitriptyline, Diphenhydramine and Trihexyphenedyl) at the water/1,2-dichloroethane interface with the phase volume ratio (r = V-o/V-w) equal to 1 are investigated by cyclic voltammetry. The system is composed of an aqueous droplet supported at an Ag/AgCl disk electrode and it was covered with an organic solution. In this manner, a conventional three-electrode potentiostat can be used to study the ionizable drugs transfer process at a liquid/liquid interface. Physicochemical parameters such as the formal transfer potential, the Gibbs energy of transfer and the standard partition coefficients of the ionized forms of these drugs can be evaluated from cyclic voltammograms obtained. The obtained results have been summarized in ionic partition diagrams, which are a useful tool for predicting and interpreting the transfer mechanisms of ionizable drugs at the liquid/liquid interfaces and biological membranes.

Study of facilitated potassium ion transfer across a water vertical bar 1,2-dichloroethane interface using different phase volume ratio systems

A study of potassium ion transfer across a water \ 1,2-dichloroethane (W \ DCE) interface facilitated by dibenzo-18-crown-6 (DB18C6) with various phase volume ratio systems is presented. The key point was that a droplet of aqueous solution containing a redox couple, Fe(CN)(6)(3-)/Fe(CN)(6)(4-), with equal molar ratio, was first attached to a platinum electrode surface, and the resulting droplet electrode was then immersed into the organic solution containing a hydrophobic electrolyte to construct a platinum electrode/aqueous phase/organic phase system. The interfacial potential of the W \ DCE within the series could be externally controlled because the specific compositions in the aqueous droplet make the Pt electrode function like a reference electrode as long as the concentration ratio of Fe(CN)(6)(3-)/Fe(CN)(6)(4-) remains constant. In this way, a conventional three-electrode potentiostat can be used to study the ion transfer process at a liquid \ liquid (L \ L) interface facilitated by an ionophore with variable phase volume ratio (r = V-o/V-w). The effect of r on ion transfer and facilitated ion transfer was studied in detail experimentally. We also demonstrated that as low as 5 x 10(-8) M DB18C6 could be determined using this method due to the effect of the high phase volume ratio.

The glass transition temperatures of PS/PPO blends: Couchman volume-based equation and its verification

The glass transition temperatures (T-g) of PS/PPO blends with different compositions were studied under various pressures by means of a PVT-100 analyzer. A general relation of T-g and pressure of the PS/PPO system was deduced by fitting the experimental T-g's. Couchman volume-based equation was testified with the aid of those data. It was found that the experimental T-g's do not obey the Couchman equation of glass transition temperature based on thermodynamic theory. According to our studies, the major reason of the deviation is caused by the neglect of DeltaV(mix). (C) 2001 Published by Elsevier Science Ltd.