Calculation of transverse hydrodynamic coefficients using computational fluid dynamic approach

Computational Fluid Dynamic (CFD) based on Reynolds Averaged Navier–Stokes equation is used for determining the transverse hydrodynamic damping force and moment coefficients that are needed in the maneuverability study of marine vehicles. Computations are performed for two geometrical shapes representing typical AUVs presently in use. Results are compared with available data on similar geometries and from some of the available semi-empirical relations. It is found that the CFD predictions compares reasonable well with these results. In particular, the CFD predictions of forces and moments are found to be nonlinear with respect to the transverse velocity, and therefore both linear and nonlinear coefficients can be derived. A discussion on the sources of the component forces reveal that the total force and moment variations should in fact be nonlinear.     

Beneath the surface: Characteristics of oceanic ecosystems under weak mixing conditions – A theoretical investigation

This study considers an important biome in aquatic environments, the subsurface ecosystem that evolves under low mixing conditions, from a theoretical point of view. Employing a conceptual model that involves phytoplankton, a limiting nutrient and sinking detritus, we use a set of key characteristics (thickness, depth, biomass amplitude/productivity) to qualitatively and quantitatively describe subsurface biomass maximum layers (SBMLs) of phytoplankton. These SBMLs are defined by the existence of two community compensation depths in the water column, which confine the layer of net community production; their depth coincides with the upper nutricline. Analysing the results of a large ensemble of simulations with a one-dimensional numerical model, we explore the parameter dependencies to obtain fundamental steady-state relationships that connect primary production, mortality and grazing, remineralization, vertical diffusion and detrital sinking. As a main result, we find that we can distinguish between factors that determine the vertically integrated primary production and others that affect only depth and shape (thickness and biomass amplitude) of this subsurface production layer. A simple relationship is derived analytically, which can be used to estimate the steady-state primary productivity in the subsurface oligotrophic ocean. The fundamental nature of the results provides further insight into the dynamics of these “hidden” ecosystems and their role in marine nutrient cycling.     

Computational study of sail performance in upwind condition

Airflow around yacht sails with imposed final geometry is simulated using a CFD code, reproducing experimental tests carried in a wind tunnel. Two configurations are considered: one, mast–main sail and the other, mast–jib–main sail. Both configurations were studied in the same flow conditions of air velocity and incidence angle. The grid is structured-like next to the sails and unstructured in the rest of the domain. The turbulence model used is Shear Stress Transport. The results are compared against experimental and numerical results.     

A virtual voyage through 3D structures generated by chaotic mixing of magmas and numerical simulations: a new approach for understanding spatial and temporal complexity of magma dynamics

In this contribution, we present a virtual voyage through 3D structures generated by chaotic mixing of magmas and numerical simulations with the aim to highlight the power of 3D representations in the understanding of this geological phenomenon. In particular, samples of mixed juveniles from Salina island (Southern Italy) are reconstructed in 3D by serial lapping and digital montage and numerical simulations are performed by using a 3D chaotic dynamical system. Natural and simulated magma mixing structures are visualized by using several multimedia tools including animations and “virtual reality” models. It is shown that magma interaction processes can generate large spatial and temporal compositional heterogeneities in magmatic systems. The same topological structures are observed in both 3D reconstructed rock samples and chaotic numerical simulations, indicating that the mixing of magmas is governed by chaotic dynamics. The use of 3D multimedia models gives the opportunity to penetrate into magma mixing structures and to understand their significance in the context of magma dynamics. Such an approach is very powerful since multimedia tools can strongly capture the attention of the reader bringing him/her into an interactive and memorable geological experience. Electronic supplementary material  enclosed:     

Probabilistic modeling of seismically triggered landslides using Monte Carlo simulations

The 2004 Chuetsu earthquakes of Niigata (Japan) triggered numerous landslides, and the most widespread types of landslides were highly disrupted, relatively shallow slides and soil (debris) flows. This paper presented a method to evaluate slope instability using Newmark displacement on a pixel-by-pixel basis in a given area. The proposed method was able to integrate Newmark displacement modeling and Monte Carlo simulations within geographical information systems. In the modeling, an empirical attenuation relationship was utilized to calculate Arias intensity over this study area, and the variability of geotechnical parameters was taken into account to calculate coseismic landslide displacement. Before deriving the displacement from related inputs, the Monte Carlo simulations ran 1,500 times and generated 1,500 displacement values for each grid cell, and then means and standard deviations of displacement were calculated and probabilistic distributions can be obtained. Finally, given 10 cm as a threshold value of displacement, estimated probabilities of displacement exceeding 10 cm were shown as a map of seismic landslide hazards. The resulting hazard map was classified into four categories from very low to high level.     

Efficient generation of conditional simulations by chebyshev matrix polynomial approximations to the symmetric square root of the covariance matrix

Consider the problem of generating a realization y₁ of a Gaussian random field on a dense grid of points ₁ conditioned on field observations y₂ collected on a sparse grid of points ₂. An approach to this is to generate first an unconditional realization y over the grid =_{1 2}, and then to produce y₁ by conditioning y on the data y₂. As standard methods for generating y, such as the turning bands, spectral or Cholesky approaches can have various limitations, it has been proposed by M. W. Davis to generate realizations from a matrix polynomial approximations to the square root of the covariance matrix. In this paper we describe how to generate a direct approximation to the conditional realization y₁, on ₁ using a variant of Davis' approach based on approximation by Chebyshev polynomials. The resulting algorithm is simple to implement, numerically stable, and bounds on the approximation error are readily available. Furthermore we show that the conditional realization y₁ can be generated directly with a lower order polynomial than the unconditional realization y, and that further reductions can be achieved by exploiting a nugget effect if one is present. A pseudocode version of the algorithm is provided that can be implemented using the fast Fourier transform if the field is stationary and the grid ₁ is rectangular. Finally, numerical illustrations are given of the algorithm's performance in generating various 2-D realizations of conditional processes on large sampling grids.     

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.     

Computer simulation studies of Crystallization under confinement conditions

Crystallization under confinement conditions is a very important process in geochemistry and geophysics. Computer simulations of fluids in nanometer scale pore spaces can provide a unique microscopic insight into the structure, dynamics and forces arising from the crystallization process. We discuss in this paper molecular dynamics computer simulations of crystallization in pores of nanometer dimensions. The crystallization pressure due to the freezing of a model of Argon in a nanopore is computed using molecular dynamics simulations. We also investigate the influence of pore geometry in determining the dynamics of confined fluids, as well as mass separation in binary mixtures. It turns out that the pore geometry reveals itself as an important variable, leading to 1) new mechanisms for fast diffusion in confined spaces, and 2) accumulation of solute in specific regions inside the pore.     

From stable dipolar towards reversing numerical dynamos

We are using a three-dimensional convection-driven numerical dynamo model without hyperdiffusivity to study the characteristic structure and time variability of the magnetic field in dependence of the Rayleigh number (Ra) for values up to 40 times supercritical. We also compare a variety of ways to drive the convection and basically find two dynamo regimes. At low Ra, the magnetic field at the surface of the model is dominated by the non-reversing axial dipole component. At high Ra, the dipole part becomes small in comparison to higher multipole components. At transitional values of Ra, the dynamo vacillates between the dipole-dominated and the multipolar regime, which includes excursions and reversals of the dipole axis. We discuss, in particular, one model of chemically driven convection, where for a suitable value of Ra, the mean dipole moment and the temporal evolution of the magnetic field resemble the known properties of the Earth’s field from paleomagnetic data.     

Simulation of individual particle movement in a gravel streambed

A two‐dimensional simulation model of travel distances of individual particles in a gravel‐bed river is presented. The model is based on a number of rules, which include particle size, entrainment, trajectory, distance of movement and entrapment. Particle interactions are controlled by resistance fields defined about each obstacle and critical elevation defined in the model. Resistance fields, particle dropping and critical elevation rules control particle interactions. The interaction rules cause the particles to develop pebble clusters, stone cells and transverse structures (transverse ribs). The simulated travel distances of individual particles are consistent with reported field results. Individual particle travel distances were simulated using two different models; one without interactions between the individual particles and the stationary bed and one with interactions. The case without interactions demonstrates the random nature of sediment transport, and narrow ranges of travel distances. Wider ranges of travel distances, similar to those for natural situations, were obtained for the cases with interactions. The more intense the interaction between the mobile stones and the stationary ones, the wider the range of distances of travel for a given particle size. Modelling the mean travel distance yielded a result similar to that published previously, which was based on empirical data. Well developed bed‐surface structures were obtained for relatively poorly sorted sediment with intense interactions between particles. Transverse structures developed when relatively large particles were allowed to move. Copyright © 2002 John Wiley & Sons, Ltd.