Recursive algorithms for the computation of the potential harmonic coefficients of a constant density polyhedron

The gravitational potential of a constant density general polyhedron can be expressed both in terms of a closed analytical expression and as a series expansion involving the corresponding spherical harmonic coefficients. The latter can be obtained from two independent algorithms, which differ not only in their algorithmic architecture but in their efficiency and overall performance, especially when computing the coefficients of higher degree and order. In the present paper a comparative study of all these three approaches is carried out focusing on the numerical implementation of the recursive relations appearing in the two algorithms for the computation of the polyhedral potential harmonic coefficients. The performed numerical investigations show that the linear algorithm proposed by Jamet and Thomas (Proceedings of the second international GOCE user workshop, ‘GOCE, The Geoid and Oceanography’, ESA-ESRIN, Frascati, Italy, 8–10 March 2004, ESA SP-569, 2004), but so far not implemented, achieves a reasonable accuracy at a computational expense that opens to practical applications, for instance in the field of satellite gravimetry/gradiometry interpretation. The convergence behavior of the linear recursion algorithm is studied thoroughly and a computational procedure is proposed that enables the stable computation of potential harmonic coefficients up to degree 60 when referring to an arbitrarily shaped polyhedral body.     

Adsorption of cobalt by using inorganic components of sediment samples from water bodies

The adsorption of cobalt ions was evaluated using sediment samples from water bodies to investigate the adsorption properties of sediment and the behavior of these natural materials in the presence of nuclear and industrial waste.The two sediments(S1 and S2)were treated to eliminate humic and fulvic acids and then they were characterized by several techniques.The minerals found in both the sediments(X-ray diffraction)were quartz and albite;plus,goethite and muscovite in S1,and kaolinite and montmorillonite in S2.Point of zero charge(PZC)of S1 and S2 was 6.00 and 5.22,respectively.The specific area of S1(63.3 m²/g)is higher than S2(1.5 m²/g).Adsorption kinetics data for S1 and S2 were best fitted to the pseudo second-order model.The removal efficiency of S1 for cobalt was 96%with an adsorption capacity(qe)of 0.93 mg/g,and for S2 was 45%with a qeof 0.40 mg/g.The experimental data of the adsorption isotherms were adjusted to Langmuir and Freundlich models for S1 and S2,respectively.The thermodynamic parameters(enthalpy,entropy,and Gibb’s free energy)indicated that the adsorption processes were endothermic,spontaneous,and chemisorption mechanism.The results show that the adsorption capacities of the sediments depend on their composition.These water sediments have important adsorption properties for cobalt,and they can be used in the treatment of nuclear and industrial aqueous wastes.     

Record of enhanced late Pliocene tectono‐volcanic activity in the Lesser Antilles forearc basin: Insights into the long‐term behaviour of the Lesser Antilles subduction zone

In this study, we investigate the Lesser Antilles forearc basin, focusing on the late Pliocene to Pleistocene sedimentary archives in order to track the occurrence of extreme events triggered by enhanced subduction‐related tectono‐volcanic activity. We identify late Piacenzian deposits covering a major regional erosional surface, displaying sedimentary dykes and large marine boulders embedded in a mixed continental–marine matrix, characteristic of tsunamites. We interpret this episode of platform emersion and the successive cataclysmic deposits as resulting from enhanced tectonic activity at the interface of the subduction zone, synchronous with the initiation of the Lesser Antilles volcanic arc. We then discuss the implications in terms of the mechanical behaviour of the Lesser Antilles subduction zone.     

Numerical prediction of an intense convective system associated with the July 1987 montreal flood. part II: A trailing stratiform rainband

Abstract

In this study, the internal circulation structures of the 14 July 1987 intense mesoscale convective system (MCS) are investigated using an improved high‐resolution version of the Canadian regional finite‐element model. It is found that although the MCS is characterized by a leading convective line followed by a trailing stratiform rainband, the associated circulation structures differ substantially from those in the classical midlatitude squall system. These include the rapid propagation and separation of the leading convection from the trailing rainband, the development of a surface‐based instead of an elevated rear‐to‐front descending flow and a shallow front‐to‐rear ascending flow associated with the stratiform precipitation, the generation of low‐ and mid‐level rather than mid‐ to upper‐level stratiform cloudiness and the development of a strong anticyclonic vorticity band at the back edge of the stratiform region. It is shown that the trailing stratiform rainband is dynamically forced by frontogenetical processes, and aided by the release of conditional symmetric instability and local orographical lifting. The intense leading and trailing circulations result from latent heat released by the convective and explicit cloud schemes, respectively. Sensitivity experiments reveal that the proper coupling of these two cloud schemes is instrumental in obtaining a realistic prediction of the above‐mentioned various mesoscale components. Vorticity budget calculations show that tilting of horizontal vorticity contributes the most to the amplification of the anticyclonic vorticity band, particularly during the squall's incipient stage. The sensitivity of the simulated squall system to other model physical parameters is also examined.     

Anticipating the effects of groundwater withdrawal on seawater intrusion and soil settlement in urban coastal areas

Intensive pumping in urban coastal areas is a common threat to water resource quality due to seawater intrusion. In those areas where subsurface water resources are not usually used for human consumption or irrigation, intensive pumping is associated with other activities like the lowering of the water table necessary to support underground structures and building foundations. This activity also increases the likelihood of soil settlement that affects building stability and the corrosion of concrete structures due to groundwater salinity. Under these circumstances, the awareness of a certain municipality (Calonge, NE Spain) of the potential effects of groundwater withdrawal upon foundations has led to an integrated approach to anticipate seawater intrusion related to urban development. Geological mapping and correlation of borehole logs, electrical resistivity tomography, and hydrochemical data provide comprehensive knowledge of the geology and hydrogeology of the area and act as screening tools necessary to discern the influence of hydrological processes in coastal areas. Developing Strack's analytical solution, new comprehensive, dimensionless expressions are herein derived to determine the critical pumping rate necessary to prevent seawater intrusion, as well as to reproduce the evolution of the wedge toe and the water table stagnation point under different withdrawal rates. Furthermore, the Dupuit–Forchheimer well discharge formula allows the estimation of the effects of the water table lowering due to such critical pumping in the surrounding building foundations. Field data from the Calonge coastal plain illustrate this approach and provide assessment criteria for future urban development and planning. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical analysis of frost effects in porous media. Benefits and limits of the finite element poroelasticity formulation

The aim of this paper is to analyse the performance of a finite element formulation usable for predicting the mechanical consequences of frost effects on porous media. It considers the characteristics of porous media and how the frost action can be assessed. The problem is then separated into two parts: thermal and poromechanical calculations. The constitutive equations developed in the framework of poromechanics are presented and the implementation in a usual finite element poroelasticity formulation based on Zuber's method is adopted. An analysis of the time‐step influence on the convergence rate is given and leads us to propose a simple method in order to obtain objectivity of the finite element response and avoid over‐long calculations. Frost effect simulations are carried out on real porous media (two fired clays) as a case study. Although the experimental behaviour of the porous media subjected to frost action is in accordance with some observations, the calculated strains appear to be overestimated compared with measurements. The problem could be largely attributable to the difficulty of assessing permeability evolution during frost development. Copyright © 2011 John Wiley & Sons, Ltd.     

Non linear waves in two dimensional keplerian flows

Abstract

Accretion discs in astrophysics are fundamental for converting gravitational binding energy into observed electromagnetic radiation. We study the behavior of waves in a two dimensional supersonic Keplerian flow inside a given gravitational potential. We present the effects of shearing and rotation on short waves, and the numerical study of the dynamical stability of such flows with respect to various perturbations. We show that a large class of dynamical effects, due to pressure and associated to short time scales, may be excited.     

TITAN2D simulations of pyroclastic flows at Cerro Machín Volcano,Colombia: Hazard implications

Cerro Machín is a dacitic tuff ring located in the central part of the Colombian Andes. It lies at the southern end of the Cerro Bravo–Cerro Machín volcanic belt. This volcano has experienced at least six major explosive eruptions during the last 5000 years. These eruptions have generated pyroclastic flows associated with Plinian activity that have traveled up to 8 km from the crater, and pyroclastic flows associated with Vulcanian activity with shorter runouts of 5 km from the source. Today, some 21,000 people live within a 8 km radius of Cerro Machín. The volcano is active with fumaroles and has shown increasing seismic activity since 2004, and therefore represents a potentially increasing threat to the local population. To evaluate the possible effects of future eruptions that may generate pyroclastic density currents controlled by granular flow dynamics we performed flow simulations with the TITAN2D code. These simulations were run in all directions around the volcano, using the input parameters of the largest eruption reported. The results show that an eruption of 0.3 km³ of pyroclastic flows from a collapsing Plinian column would travel up to 9 km from the vent, emplacing a deposit thicker than 60 m within the Toche River valley. Deposits >45 m thick can be expected in the valleys of San Juan, Santa Marta, and Azufral creeks, while 30 m thick deposits could accumulate within the drainages of the Tochecito, Bermellón, and Coello Rivers. A minimum area of 56 km² could be affected directly by this kind of eruption. In comparison, Vulcanian column-collapse pyroclastic flows of 0.1 km³ would travel up to 6 km from the vent depositing >45 m thick debris inside the Toche River valley and more than 30 m inside the valleys of San Juan, Santa Marta, and Azufral creeks. The minimum area that could be affected directly by this kind of eruption is 33 km². The distribution and thickness of the deposits obtained by these simulations are consistent with the hazard map presented by INGEOMINAS (Geological Survey of Colombia) in 2002. The composite map of the simulated flow deposits suggests that after major explosive events such as these, the generation of lahars is probable.