Optimization of pipeline transport for CO2 sequestration

Coal fired power generation will continue to provide energy to the world for the foreseeable future. However, this energy use is a significant contributor to increased atmospheric CO2 concentration and, hence, global warming. Capture and disposal Of CO2 has received increased R&D attention in the last decade as the technology promises to be the most cost effective for large scale reductions in CO2 emissions. This paper addresses CO2 transport via pipeline from capture site to disposal site, in terms of system optimization, energy efficiency and overall economics. Technically, CO2 can be transported through pipelines in the form of a gas, a supercritical. fluid or in the subcooled liquid state. Operationally, most CO2 pipelines used for enhanced oil recovery transport CO2 as a supercritical fluid. In this paper, supercritical fluid and subcooled liquid transport are examined and compared, including their impacts on energy efficiency and cost. Using a commercially available process simulator, ASPEN PLUS 10.1, the results show that subcooled liquid transport maximizes the energy efficiency and minimizes the Cost Of CO2 transport over long distances under both isothermal and adiabatic conditions. Pipeline transport of subcooled liquid CO2 can be ideally used in areas of cold climate or by burying and insulating the pipeline. In very warm climates, periodic refrigeration to cool the CO2 below its critical point of 31.1 degrees C, may prove economical. Simulations have been used to determine the maximum safe pipeline distances to subsequent booster stations as a function of inlet pressure, environmental temperature and ground level heat flux conditions. (c) 2005 Published by Elsevier Ltd.

Application of a generalized thermodynamic approach to characterize mesoporous materials

Equilibrium adsorption and desorption in mesoporous adsorbents is considered on the basis of rigorous thermodynamic analysis, in which the curvature-dependent solid-fluid potential and the compressibility of the adsorbed phase are accounted for. The compressibility of the adsorbed phase is considered for the first time in the literature in the framework of a rigorous thermodynamic approach. Our model is a further development of continuum thermodynamic approaches proposed by Derjaguin and Broekhoff and de Boer, and it is based on a reference isotherm of a non-porous material having the same chemical structure as that of the pore wall. In this improved thermodynamic model, we incorporated a prescription for transforming the solid-fluid potential exerted by the flat reference surface to the potential inside cylindrical and spherical pores. We relax the assumption that the adsorbed film density is constant and equal to that of the saturated liquid. Instead, the density of the adsorbed fluid is allowed to vary over the adsorbed film thickness and is calculated by an equation of state. As a result, the model is capable to describe the adsorption-desorption reversibility in cylindrical pores having diameter less than 2 nm. The generalized thermodynamic model may be applied to the pore size characterization of mesoporous materials instead of much more time-consuming molecular approaches. (c) 2005 Elsevier B.V. All rights reserved.

Life Cycle Assessment (LCA) Applied to the Design of an Innovative Drying Process for Sewage Sludge

Most adverse environmental impacts result from design decisions made long before manufacturing or usage. In order to prevent this situation, several authors have proposed the application of life cycle assessment (LCA) at the very first phases of the design of a process, a product or a service. The study in this paper presents an innovative thermal drying process for sewage sludge called fry-drying, in which dewatered sludge is directly contacted in the dryer with hot recycled cooking oils (RCO) as the heat medium. Considering the practical difficulties for the disposal of these two wastes, fry-drying presents a potentially convenient method for their combined elimination by incineration of the final fry-dried sludge. An analytical comparison between a conventional drying process and the new proposed fry-drying process is reported, with reference to some environmental impact categories. The results of this study, applied at the earliest stages of the design of the process, assist evaluation of the feasibility of such system compared to a current disposal process for the drying and incineration of sewage sludge.

On the chaotic rotation of a liquid-filled gyrostat via the Melnikov-Holmes-Marsden integral

The non-linear motions of a gyrostat with an axisymmetrical, fluid-filled cavity are investigated. The cavity is considered to be completely filled with an ideal incompressible liquid performing uniform rotational motion. Helmholtz theorem, Euler's angular momentum theorem and Poisson equations are used to develop the disturbed Hamiltonian equations of the motions of the liquid-filled gyrostat subjected to small perturbing moments. The equations are established in terms of a set of canonical variables comprised of Euler angles and the conjugate angular momenta in order to facilitate the application of the Melnikov-Holmes-Marsden (MHM) method to investigate homoclinic/heteroclinic transversal intersections. In such a way, a criterion for the onset of chaotic oscillations is formulated for liquid-filled gyrostats with ellipsoidal and torus-shaped cavities and the results are confirmed via numerical simulations. (c) 2006 Elsevier Ltd. All rights reserved.

Analysis of chaotic instabilities in a rotating body with internal energy dissipation

Melnikov's method is used to analytically predict the onset of chaotic instability in a rotating body with internal energy dissipation. The model has been found to exhibit chaotic instability when a harmonic disturbance torque is applied to the system for a range of forcing amplitude and frequency. Such a model may be considered to be representative of the dynamical behavior of a number of physical systems such as a spinning spacecraft. In spacecraft, disturbance torques may arise under malfunction of the control system, from an unbalanced rotor, from vibrations in appendages or from orbital variations. Chaotic instabilities arising from such disturbances could introduce uncertainties and irregularities into the motion of the multibody system and consequently could have disastrous effects on its intended operation. A comprehensive stability analysis is performed and regions of nonlinear behavior are identified. Subsequently, the closed form analytical solution for the unperturbed system is obtained in order to identify homoclinic orbits. Melnikov's method is then applied on the system once transformed into Hamiltonian form. The resulting analytical criterion for the onset of chaotic instability is obtained in terms of critical system parameters. The sufficient criterion is shown to be a useful predictor of the phenomenon via comparisons with numerical results. Finally, for the purposes of providing a complete, self-contained investigation of this fundamental system, the control of chaotic instability is demonstated using Lyapunov's method.

Suppressing chaos via Lyapunov-Krasovskii's method

An algorithm for suppressing the chaotic oscillations in non-linear dynamical systems with singular Jacobian matrices is developed using a linear feedback control law based upon the Lyapunov-Krasovskii (LK) method. It appears that the LK method can serve effectively as a generalised method for the suppression of chaotic oscillations for a wide range of systems. Based on this method, the resulting conditions for undisturbed motions to be locally or globally stable are sufficient and conservative. The generalized Lorenz system and disturbed gyrostat equations are exemplified for the validation of the proposed feedback control rule. (c) 2005 Elsevier Ltd. All rights reserved.

Prediction of chaotic instabilities in a dragline bucket swing

The occurrence of chaotic instabilities is investigated in the swing motion of a dragline bucket during operation cycles. A dragline is a large, powerful, rotating multibody system utilised in the mining industry for removal of overburden. A simplified representative model of the dragline is developed in the form of a fundamental non-linear rotating multibody system with energy dissipation. An analytical predictive criterion for the onset of chaotic instability is then obtained in terms of critical system parameters using Melnikov's method. The model is shown to exhibit chaotic instability due to a harmonic slew torque for a range of amplitudes and frequencies. These chaotic instabilities could introduce irregularities into the motion of the dragline system, rendering the system difficult to control by the operator and/or would have undesirable effects on dragline productivity and fatigue lifetime. The sufficient analytical criterion for the onset of chaotic instability is shown to be a useful predictor of the phenomenon under steady and unsteady slewing conditions via comparisons with numerical results. (c) 2005 Elsevier Ltd. All rights reserved.

Atoll lagoon flushing forced by waves

Water level and current measurements from two virtually enclosed South Pacific atolls, Manihiki and Rakahanga, support a new lagoon flushing mechanism which is driven by waves and modulated by the ocean tide for virtually enclosed atolls. This is evident because the lagoon water level remains above the ocean at all tidal phases (i.e., ruling out tidal flushing) and because the average lagoon water level rises significantly during periods with large waves. Hence, we develop a model by which the lagoons are flushed by waves pumping of ocean water into the lagoon and gravity draining water from the lagoon over the reef rim. That is, the waves on the exposed side push water into the lagoon during most of the tidal cycle while water leaves the lagoon on the protected side for most of the tidal cycle. This wave-driven through flow flushing is shown to be more efficient than alternating tidal flushing with respect to water renewal. Improved water quality should therefore be sought through enhancement of the natural wave pumping rather than by blasting deep channels which would change the system to an alternating tide-driven one.

Numerical solutions of the sediment conservation law: A review and improved formulation for coastal morphological modelling

Numerical solutions of the sediment conservation law are reviewed in terms of their application to bed update schemes in coastal morphological models. It is demonstrated that inadequately formulated numerical techniques lead to the introduction of diffusion, dispersion and the bed elevation oscillations previously reported in the literature. Four different bed update schemes are then reviewed and tested against benchmark analytical solutions. These include a first order upwind scheme, two Lax-Wendroff schemes and a non-oscillating centred scheme (NOCS) recently applied to morphological modelling by Saint-Cast [Saint-Cast, F., 2002. Modelisation de la morphodynamique des corps sableux en milieu littoral (Modelling of coastal sand banks morphodynamics), University Bordeaux 1, Bordeaux, 245 pp.]. It is shown that NOCS limits and controls numerical errors while including all the sediment flux gradients that control morphological change. Further, no post solution filtering is required, which avoids difficulties with selecting filter strength. Finally, NOCS is compared to a recent Lax-Wendroff scheme with post-solution filtering for a longer term simulation of the morphological evolution around a trained river entrance. (C) 2006 Elsevier B.V. All rights reserved.

Application of a coupled ground-surface water flow model to simulate periodic groundwater flow influenced by a sloping boundary, capillarity and vertical flows

Theoretical developments as well as field and laboratory data have shown the influence of the capillary fringe on water table fluctuations to increase with the fluctuation frequency. The numerical solution of a full, partially saturated flow equation can be computationally expensive. In this paper, the influence of the capillary fringe on water table fluctuations is simplified through its parameterisation into the storage coefficient of a fully-saturated groundwater flow model using the complex effective porosity concept [Nielsen, P., Perrochet, P., 2000. Water table dynamics under capillary fringes: experiments and modelling. Advances in Water Resources 23 (1), 503-515; Nielsen, P., Perrochet, P., 2000. ERRATA: water table dynamics under capillary fringes: experiments and modelling (Advances in Water Resources 23 (2000) 503-515). Advances in Water Resources 23, 907-908]. The model is applied to sand flume observations of periodic water table fluctuations induced by simple harmonic forcing across a sloping boundary, analogous to many beach groundwater systems. While not providing information on the moisture distribution within the aquifer, this approach can reasonably predict the water table fluctuations in response to periodic forcing across a sloping boundary. Furthermore, he coupled ground-surface water model accurately predicts the extent of the seepage face formed at the sloping boundary. (C) 2005 Elsevier Ltd. All rights reserved.

Sheet flow sediment transport under waves with acceleration skewness and boundary layer streaming

The effect of acceleration skewness on sheet flow sediment transport rates (q) over bar (s) is analysed using new data which have acceleration skewness and superimposed currents but no boundary layer streaming. Sediment mobilizing forces due to drag and to acceleration (similar to pressure gradients) are weighted by cosine and sine, respectively, of the angle phi(.)(tau)phi(tau) = 0 thus corresponds to drag dominated sediment transport, (q) over bar (s)similar to vertical bar u(infinity)vertical bar u(infinity), while phi(tau) = 90 degrees corresponds to total domination by the pressure gradients, (q) over bar similar to du(infinity)/dt. Using the optimal angle, phi = 51 degrees based on that data, good agreement is subsequently found with data that have strong influence from boundary layer streaming. Good agreement is also maintained with the large body of U-tube data simulating sine waves with superimposed currents and second-order Stokes waves, all of which have zero acceleration skewness. The recommended model can be applied to irregular waves with arbitrary shape as long as the assumption negligible time lag between forcing and sediment transport rate is valid. With respect to irregular waves, the model is much easier to apply than the competing wave-by-wave models. Issues for further model developments are identified through a comprehensive data review.

Large scale experiments on gravel and mixed beaches: Experimental procedure, data documentation and initial results

This paper provides information on the experimental set-up, data collection methods and results to date for the project Large scale modelling of coarse grained beaches, undertaken at the Large Wave Channel (GWK) of FZK in Hannover by an international group of researchers in Spring 2002. The main objective of the experiments was to provide full scale measurements of cross-shore processes on gravel and mixed beaches for the verification and further development of cross-shore numerical models of gravel and mixed sediment beaches. Identical random and regular wave tests were undertaken for a gravel beach and a mixed sand/gravel beach set up in the flume. Measurements included profile development, water surface elevation along the flume, internal pressures in the swash zone, piezometric head levels within the beach, run-up, flow velocities in the surf-zone and sediment size distributions. The purpose of the paper is to present to the scientific community the experimental procedure, a summary of the data collected, some initial results, as well as a brief outline of the on-going research being carried out with the data by different research groups. The experimental data is available to all the scientific community following submission of a statement of objectives, specification of data requirements and an agreement to abide with the GWK and EU protocols. (C) 2005 Elsevier B.V. All rights reserved.

Hydraulics of skimming flows on stepped chutes: The effects of inflow conditions?

Modern stepped spillways are typically designed for large discharge capacities corresponding to a skimming flow regime for which flow resistance is predominantly form drag. The writer demonstrates that the inflow conditions have some effect on the skimming flow properties. Boundary layer calculations show that the flow properties at inception of free-surface aeration are substantially different with pressurized intake. The re-analysis of experimental results highlights that the equivalent Darcy friction factor is f similar to 0.2 in average on uncontrolled stepped Chute and f similar to 0.1 on stepped chute with pressurized intake. A simple design chart is presented to estimate the residual flow velocity, and the agreement of the calculations with experimental results is deemed satisfactory for preliminary design.

Dam break wave of thixotropic fluid

Thixotropy is the characteristic of a fluid to form a gelled structure over time when it is not subjected to shearing, and to liquefy when agitated. Thixotropic fluids are commonly used in the construction industry (e.g., liquid concrete and drilling fluids), and related applications include some forms of mud flows and debris flows. This paper describes a basic study of dam break wave with thixotropic fluid. Theoretical considerations were developed based upon a kinematic wave approximation of the Saint-Venant equations down a prismatic sloping channel. A very simple thixotropic model, which predicts the basic theological trends of such fluids, was used. It describes the instantaneous state of fluid structure by a single parameter. The analytical solution of the basic flow motion and theology equations predicts three basic flow regimes depending upon the fluid properties and flow conditions, including the initial degree of jamming of the fluid (related to its time of restructuration at rest). These findings were successfully compared with systematic bentonite suspension experiments. The present work is the first theoretical analysis combining the basic principles of unsteady flow motion with a thixotropic fluid model and systematic laboratory experiments.

Estimation of water retention curve from mercury intrusion porosimetry and van Genuchten model

The water retention curve (WRC) is a hydraulic characteristic of concrete required for advanced modeling of water (and thus solute) transport in variably saturated, heterogeneous concrete. Unfortunately, determination by a direct experimental method (for example, measuring equilibrium moisture levels of large samples stored in constant humidity cells) is a lengthy process, taking over 2 years for large samples. A surrogate approach is presented in which the WRC is conveniently estimated from mercury intrusion porosimetry (MIP) and validated by water sorption isotherms: The well-known Barrett, Joyner and Halenda (BJH) method of estimating the pore size distribution (PSD) from the water sorption isotherm is shown to complement the PSD derived from conventional MIP. This provides a basis for predicting the complete WRC from MIP data alone. The van Genuchten equation is used to model the combined water sorption and MIP results. It is a convenient tool for describing water retention characteristics over the full moisture content range. The van Genuchten parameter estimation based solely on MIP is shown to give a satisfactory approximation to the WRC, with a simple restriction on one. of the parameters.