(Table 2) Bottom-hole temperature measurements at DSDP Sites 68-501, 69-504, and 69-505

Drilling durin Deep Sea Drilling Project Legs 68, 69, and 70 on the southern limb of the Costa Rica Rift was used to study geothermal processes in the ocean crust. Two areas were drilled. One was a geothermally hot site on 6.2-m.y.-old crust, where topography is smooth, heat flow is close to that predicted by conductive cooling of the lithosphere (200 mWm**-2), and hydrothermal circulation may be sealed within the crust. The other was on 3.9-m.y.-old crust, where rough topography is associated with low heat flow (15 to 50 mWm**-2) and possible open convection of sea water. At both sites, about 250 m of siliceous-calcareous sediments overlies igneous basement. In the hot area, it blankets the topography, whereas in the cold area, basement outcrops still occur. Operations included numerous down-hole experiments in both areas and hydraulic piston coring of a 230-m sediment section in the hot area. Diagenesis of the sediments appears closely related to temperature. At the hot site, chert was found near basement, and the chemistry of pore fluids, sampled from both sediments and basement, is strongly influenced by reactions within the basement. Strong lateral gradients in the composition of pore fluids occur in the sediments. At the cold site, no chert was found, and bacterial processes within the sediment dominated the chemistry of the pore fluids. Basaltic basement in both areas consists mainly of pillow lavas and thin flows, with occasional more massive units. The basalt is relatively magnesian. The degree of alteration is very small in the cold area, but much more extensive in the hot area. Ease of drilling also shows a strong contrast. Basement penetration reached 562 m in the hot area and was halted because of lack of time; at the cold site, 43 m of basement was cored only with difficulty. The most intensive in-hole experiments were conducted in the hot area. Successful runs with the borehole televiewer allowed basement lithology to be determined and showed the presence of more and less fractured zones. Pulse tests using a single borehole packer gave values of basement permeability of about 2 to 40 millidarcies. Numerous temperature logs established a broadly conductive in situ temperature gradient, with temperatures reaching 120°C at 562 m into the basement. However, anomalously low temperatures in the upper part of the hole, which persisted after drilling disturbance had decayed away, showed that cold ocean water was flowing down the hole and into the basement at about 90 m below the base of the sediments, at rates of about 80 to 100 m/hr. The packer records indicate a pressure at this depth of 10 bars below hydrostatic.

Shallow water marine sediment bacterial community shifts along a natural CO2 gradient in the Mediterranean Sea Off vulcano, Italy

The effects of increasing atmospheric CO(2) on ocean ecosystems are a major environmental concern, as rapid shoaling of the carbonate saturation horizon is exposing vast areas of marine sediments to corrosive waters worldwide. Natural CO(2) gradients off Vulcano, Italy, have revealed profound ecosystem changes along rocky shore habitats as carbonate saturation levels decrease, but no investigations have yet been made of the sedimentary habitat. Here, we sampled the upper 2 cm of volcanic sand in three zones, ambient (median pCO(2) 419 µatm, minimum Omega (arag) 3.77), moderately CO(2)-enriched (median pCO(2) 592 µatm, minimum Omega (arag) 2.96), and highly CO(2)-enriched (median pCO(2) 1611 µatm, minimum Omega (arag) 0.35). We tested the hypothesis that increasing levels of seawater pCO(2) would cause significant shifts in sediment bacterial community composition, as shown recently in epilithic biofilms at the study site. In this study, 454 pyrosequencing of the V1 to V3 region of the 16S rRNA gene revealed a shift in community composition with increasing pCO(2). The relative abundances of most of the dominant genera were unaffected by the pCO(2) gradient, although there were significant differences for some 5 % of the genera present (viz. Georgenia, Lutibacter, Photobacterium, Acinetobacter, and Paenibacillus), and Shannon Diversity was greatest in sediments subject to long-term acidification (>100 years). Overall, this supports the view that globally increased ocean pCO(2) will be associated with changes in sediment bacterial community composition but that most of these organisms are resilient. However, further work is required to assess whether these results apply to other types of coastal sediments and whether the changes in relative abundance of bacterial taxa that we observed can significantly alter the biogeochemical functions of marine sediments.

Rice flows across regions in Madagascar

Although spatial disintegration of rice markets in Madagascar has been well documented, little is known about actual rice flows across regions. Using weekly collected unique data from rice markets in 22 regional capital cities for one year, this study explores the physical distribution of rice in Madagascar and reveals that rice flows from sufficient regions to deficit regions, along with geographical proximity, have positive effects on rice flow. In contrast, season factors, such as harvest/non-harvest periods and weather conditions have a negligible effect on rice flow. These findings suggest that rice flows generally follow a rational pattern despite the fact that market indicators indicate that rice trade across regions is underdeveloped.

Uso de DES y RANS no estacionario para la simulación del flujo alrededor de una carena Wigley - Use of Unsteady RANS and DES to simulate flows around a Wigley Hull)

Computer Fluid Dynamics tools have already become a valuable instrument for Naval Architects during the ship design process, thanks to their accuracy and the available computer power. Unfortunately, the development of RANSE codes, generally used when viscous effects play a major role in the flow, has not reached a mature stage, being the accuracy of the turbulence models and the free surface representation the most important sources of uncertainty. Another level of uncertainty is added when the simulations are carried out for unsteady flows, as those generally studied in seakeeping and maneuvering analysis and URANS equations solvers are used. Present work shows the applicability and the benefits derived from the use of new approaches for the turbulence modeling (Detached Eddy Simulation) and the free surface representation (Level Set) on the URANS equations solver CFDSHIP-Iowa. Compared to URANS, DES is expected to predict much broader frequency contents and behave better in flows where boundary layer separation plays a major role. Level Set methods are able to capture very complex free surface geometries, including breaking and overturning waves. The performance of these improvements is tested in set of fairly complex flows, generated by a Wigley hull at pure drift motion, with drift angle ranging from 10 to 60 degrees and at several Froude numbers to study the impact of its variation. Quantitative verification and validation are performed with the obtained results to guarantee their accuracy. The results show the capability of the CFDSHIP-Iowa code to carry out time-accurate simulations of complex flows of extreme unsteady ship maneuvers. The Level Set method is able to capture very complex geometries of the free surface and the use of DES in unsteady simulations highly improves the results obtained. Vortical structures and instabilities as a function of the drift angle and Fr are qualitatively identified. Overall analysis of the flow pattern shows a strong correlation between the vortical structures and free surface wave pattern. Karman-like vortex shedding is identified and the scaled St agrees well with the universal St value. Tip vortices are identified and the associated helical instabilities are analyzed. St using the hull length decreases with the increase of the distance along the vortex core (x), which is similar to results from other simulations. However, St scaled using distance along the vortex cores shows strong oscillations compared to almost constants for those previous simulations. The difference may be caused by the effect of the free-surface, grid resolution, and interaction between the tip vortex and other vortical structures, which needs further investigations. This study is exploratory in the sense that finer grids are desirable and experimental data is lacking for large α, especially for the local flow. More recently, high performance computational capability of CFDSHIP-Iowa V4 has been improved such that large scale computations are possible. DES for DTMB 5415 with bilge keels at α = 20º were conducted using three grids with 10M, 48M and 250M points. DES analysis for flows around KVLCC2 at α = 30º is analyzed using a 13M grid and compared with the results of DES on the 1.6M grid by. Both studies are consistent with what was concluded on grid resolution herein since dominant frequencies for shear-layer, Karman-like, horse-shoe and helical instabilities only show marginal variation on grid refinement. The penalties of using coarse grids are smaller frequency amplitude and less resolved TKE. Therefore finer grids should be used to improve V&V for resolving most of the active turbulent scales for all different Fr and α, which hopefully can be compared with additional EFD data for large α when it becomes available.

Simulation of plasma flows in divergent magnetic nozzles

Images from the simulation code DIMAGNO illustrate the roles of pressure, electric, and magnetic forces in the 2-D plasma expansion in a magnetic nozzle and the generation of thrust

WCSPH viscosity diffusion processes in vortex flows

The aim of this paper is to clarify the role played by the most commonly used viscous terms in simulating viscous laminar flows using the weakly compressible approach in the context of smooth particle hydrodynamics (WCSPH). To achieve this, Takeda et al. (Prog. Theor. Phys. 1994; 92(5):939–960), Morris et al. (J. Comput. Phys. 1997; 136:214–226) and Monaghan–Cleary–Gingold's (Appl. Math. Model. 1998; 22(12):981–993; Monthly Notices of the Royal Astronomical Society 2005; 365:199–213) viscous terms will be analysed, discussing their origins, structures and conservation properties. Their performance will be monitored with canonical flows of which related viscosity phenomena are well understood, and in which boundary effects are not relevant. Following the validation process of three previously published examples, two vortex flows of engineering importance have been studied. First, an isolated Lamb–Oseen vortex evolution where viscous effects are dominant and second, a pair of co-rotating vortices in which viscous effects are combined with transport phenomena. The corresponding SPH solutions have been compared to finite-element numerical solutions. The SPH viscosity model's behaviour in modelling the viscosity related effects for these canonical flows is adequate

Theoretical analysis and numerical verification of the consistency of viscous smoothed-particle-hydrodynamics formulations in simulating free-surface flows

The theoretical formulation of the smoothed particle hydrodynamics (SPH) method deserves great care because of some inconsistencies occurring when considering free-surface inviscid flows. Actually, in SPH formulations one usually assumes that (i) surface integral terms on the boundary of the interpolation kernel support are neglected, (ii) free-surface conditions are implicitly verified. These assumptions are studied in detail in the present work for free-surface Newtonian viscous flow. The consistency of classical viscous weakly compressible SPH formulations is investigated. In particular, the principle of virtual work is used to study the verification of the free-surface boundary conditions in a weak sense. The latter can be related to the global energy dissipation induced by the viscous term formulations and their consistency. Numerical verification of this theoretical analysis is provided on three free-surface test cases including a standing wave, with the three viscous term formulations investigated.

Latin American Agricultural Trade: The Role of the WTO in Sustainable Virtual Water Flows, European Association of Agricultural Economists

International agricultural trade has been growing significantly during the last decade. Many countries rely on imports to ensure adequate food supplies to the people. A few are becoming food baskets of the world. This process raises issues about the food security in depending countries and potentially unsustainable land and water use in exporting countries. In this paper, we analyse the impacts of amplified farm trade on natural resources, especially water. Farm exports and imports of five Latin America countries (Brazil, Argentina, Mexico, Peru and Chile) are examined carefully. A preliminary analysis indicates that virtual water imports can save valuable water resources in water-short countries, such as Mexico and Chile. Major exporting countries, including Brazil and Argentina, have become big exporters due to abundant natural resource endowments. The opportunity costs of agricultural production in those countries are identified as being low, because of the predominant green water use. It is concluded that virtual water trade can be a powerful tool to alleviate water stress in semi-arid countries. However, for exporting nations a sustainable water use can only be guaranteed if environmental production costs are fully reflected in the commodity prices. There is no basis for erecting environmental trade tariffs on exporters though. Setting up legal foundations for them in full compliance with WTOs processes would be a daunting task.

Effect of Aspect Ratio on the Three-Dimensional Global Instability Analysis of Incompressible Open Cavity Flows

The stability analysis of open cavity flows is a problem of great interest in the aeronautical industry. This type of flow can appear, for example, in landing gears or auxiliary power unit configurations. Open cavity flows is very sensitive to any change in the configuration, either physical (incoming boundary layer, Reynolds or Mach numbers) or geometrical (length to depth and length to width ratio). In this work, we have focused on the effect of geometry and of the Reynolds number on the stability properties of a threedimensional spanwise periodic cavity flow in the incompressible limit. To that end, BiGlobal analysis is used to investigate the instabilities in this configuration. The basic flow is obtained by the numerical integration of the Navier-Stokes equations with laminar boundary layers imposed upstream. The 3D perturbation, assumed to be periodic in the spanwise direction, is obtained as the solution of the global eigenvalue problem. A parametric study has been performed, analyzing the stability of the flow under variation of the Reynolds number, the L/D ratio of the cavity, and the spanwise wavenumber β. For consistency, multidomain high order numerical schemes have been used in all the computations, either basic flow or eigenvalue problems. The results allow to define the neutral curves in the range of L/D = 1 to L/D = 3. A scaling relating the frequency of the eigenmodes and the length to depth ratio is provided, based on the analysis results.

The PSE-3D instability analysis methodology for flows depending strongly on two and weakly on the third spatial dimension

The present contribution discusses the development of a PSE-3D instability analysis algorithm, in which a matrix forming and storing approach is followed. Alternatively to the typically used in stability calculations spectral methods, new stable high-order finitedifference-based numerical schemes for spatial discretization 1 are employed. Attention is paid to the issue of efficiency, which is critical for the success of the overall algorithm. To this end, use is made of a parallelizable sparse matrix linear algebra package which takes advantage of the sparsity offered by the finite-difference scheme and, as expected, is shown to perform substantially more efficiently than when spectral collocation methods are used. The building blocks of the algorithm have been implemented and extensively validated, focusing on classic PSE analysis of instability on the flow-plate boundary layer, temporal and spatial BiGlobal EVP solutions (the latter necessary for the initialization of the PSE-3D), as well as standard PSE in a cylindrical coordinates using the nonparallel Batchelor vortex basic flow model, such that comparisons between PSE and PSE-3D be possible; excellent agreement is shown in all aforementioned comparisons. Finally, the linear PSE-3D instability analysis is applied to a fully three-dimensional flow composed of a counter-rotating pair of nonparallel Batchelor vortices.