Corrosive wear of cylinder liners in large two-stroke marine diesel engines that burn heavy fuel oil containing sulfur is coupled to the formation of gaseous sulfur trioxide (SO3) and subsequent combined condensation of sulfuric acid (H2SO4) and water (H2O) vapor. The present work seeks to address how fuel sulfur content, charge air humidity and liner temperature variations affects the deposition of water and sulfuric acid at low load operation. A phenomenological engine model is applied to simulate the formation of cylinder/bulk gas combustion products and dew points comply with H2O–H2SO4 vapor liquid equilibrium. By assuming homogenous cylinder gas mixtures condensation is modeled using a convective heat and mass transfer analogy combined with realistic liner temperature profiles. Condensation of water is significantly altered by the liner temperature and charge air humidity while sulfuric acid condensation (the order is a few mg per cylinder every cycle) is proportional to the fuel sulfur content. Condensation takes place primarily in the upper part of the cylinder liner where a reduction of the surface temperature or saturated charge air provides that the deposited acid can be highly diluted with water.
This article gives a review of techniques applied to make sea state estimation on the basis of measured responses on a ship. The general concept of the procedures is similar to that of a classical wave buoy, which exploits a linear assumption between waves and the associated motions. In the frequency domain, this assumption yields the mathematical relation between the measured motion spectra and the directional wave spectrum. The analogy between a buoy and a ship is clear, and the author has worked on this wave buoy analogy for about fifteen years. In the article, available techniques for shipboard sea state estimation are addressed, but with a focus on only the wave buoy analogy. Most of the existing work is based on methods established in the frequency domain but, to counteract disadvantages of the frequency-domain procedures, newer studies are working also on procedures formulated directly in the time domain. Sample results from several studies are included, and the main findings from these are mentioned.
In this work, three-dimensional computational fluid dynamics (CFD) studies of sulphur oxides (SOx) and sulphuric acid (H2SO4) formation processes in a large, low speed two-stroke marine diesel engine are carried out. The current numerical study aims to investigate the conversion of sulphuric dioxide (SO2) to sulphuric trioxide (SO3) and the possibility of H2SO4 condensation which are the prerequisites to better understand the corrosion-induced wear phenomenon. This is achieved with the aid of the implementation of a multicomponent surrogate model, which comprises a skeletal n-heptane mechanism and a reduced sulphur subset mechanism. In the present work, performance of the coupled CFD-chemical kinetic model is evaluated using both qualitative and quantitative methods. The modelling results show that the temporal and spatial evolutions of SOx predicted by the skeletal model are similar to those by the base mechanism. Predictions of the variations of SOx and the associated SO2 to SO3 conversion in response to the change of fuel sulphur content, swirl velocity, start of injection, scavenge pressure and humidity qualitatively agree with numerical and experimental results from the literature. The model is further evaluated using the measured SO2 to SO3 conversion levels in a low load, low scavenge pressure case and a low load, high scavenge pressure case. The absolute values of simulated and measured conversion levels are close, although the former appear to be higher. The current results show that the flame impinges at the cylinder liner near top dead centre. The gas is cooled rapidly by the wall temperature and H2SO4 is produced in the region where the local temperature is less than 600 K. Based on the flue gas correlation, the acid dew point temperature is higher than the wall temperature, suggesting that acid condensation may begin early at the top part of the cylinder liner. The predicted distribution corresponds well with the distribution of corroded parts observed in service engines. The model is expected to serve as an important tool to simulate the rates of SO2 absorption into lubricating oil film and H2SO4 condensation in this combustion system.
In this reported work, multi-dimensional computational fluid dynamics studies of diesel combustion and soot formation processes in a constant volume combustion chamber and a marine diesel engine are carried out. The key interest here is firstly to validate the coupling of a newly developed skeletal n-heptane mechanism and a revised multi-step soot model using laser extinction measurements of diesel soot obtained at different ambient pressure levels in an optical accessible, constant volume chamber experiment. It is revealed that ignition delay times and liftoff lengths generated using the new skeletal model are close to those produced by the larger and more comprehensive chemical mechanisms, apart from those at the low pressure condition. The current study also demonstrates that the variation of averaged soot volume fraction with respect to the change of combustion chamber pressure captured using the revised soot model agrees reasonably well with the measurements in terms of peak values. The numerical model is subsequently applied to investigate the flame development, soot/nitrogen monoxide formation and heat transfer in a two-stroke, low-speed uniflow-scavenged marine diesel engine operating at full load condition, where optical measurements are not available. Comparisons to the experimental data show that the simulated pressure rise starts 1.0 crank angle degree in advance and the calculated peak pressure is 1.7% lower. The associated flame liftoff length is negligible, yielding higher local equivalence ratio and soot volume fraction values as compared to those under similar test condition in the constant volume chamber. With the use of the revised model, the total heat transfer to the walls calculated when soot radiative heat loss is taken into account is approximately 30% higher compared to that when only convective heat loss is considered. The averaged nitrogen monoxide concentration is 7.7% lower when both convective and soot radiative heat losses are accounted for but the net soot mass production is less sensitive to soot radiation. A sensitivity study reveals that neither increasing nor decreasing the soot absorption coefficient by 30% from the baseline setup is influential to nitrogen monoxide formation, soot mass production and heat transfer. The findings here aid to gain insights and provide a better understanding of the combustion and soot processes in large, uniflow-scavenged marine engines. The numerical model developed in this work can also be applied to explore different phenomena in this combustion system.
A major challenge in next-generation industrial applications is to improve numerical analysis by quantifying uncertainties in predictions. In this work we present a formulation of a fully nonlinear and dispersive potential flow water wave model with random inputs for the probabilistic description of the evolution of waves. The model is analyzed using random sampling techniques and nonintrusive methods based on generalized polynomial chaos (PC). These methods allow us to accurately and efficiently estimate the probability distribution of the solution and require only the computation of the solution at different points in the parameter space, allowing for the reuse of existing simulation software. The choice of the applied methods is driven by the number of uncertain input parameters and by the fact that finding the solution of the considered model is computationally intensive. We revisit experimental benchmarks often used for validation of deterministic water wave models. Based on numerical experiments and assumed uncertainties in boundary data, our analysis reveals that some of the known discrepancies from deterministic simulation in comparison with experimental measurements could be partially explained by the variability in the model input. Finally, we present a synthetic experiment studying the variance-based sensitivity of the wave load on an offshore structure to a number of input uncertainties. In the numerical examples presented the PC methods exhibit fast convergence, suggesting that the problem is amenable to analysis using such methods.
We present an arbitrary-order spectral element method for general-purpose simulation of non-overturning water waves, described by fully nonlinear potential theory. The method can be viewed as a high-order extension of the classical finite element method proposed by Cai et al. (1998)[5], although the numerical implementation differs greatly. Features of the proposed spectral element method include: nodal Lagrange basis functions, a general quadrature-free approach and gradient recovery using global L2projections. The quartic nonlinear terms present in the Zakharov form of the free surface conditions can cause severe aliasing problems and consequently numerical instability for marginally resolved or very steep waves. We show how the scheme can be stabilised through a combination of over-integration of the Galerkin projections and a mild spectral filtering on a per element basis. This effectively removes any aliasing driven instabilities while retaining the high-order accuracy of the numerical scheme. The additional computational cost of the over-integration is found insignificant compared to the cost of solving the Laplace problem. The model is applied to several benchmark cases in two dimensions. The results confirm the high order accuracy of the model (exponential convergence), and demonstrate the potential for accuracy and speedup. The results of numerical experiments are in excellent agreement with both analytical and experimental results for strongly nonlinear and irregular dispersive wave propagation. The benefit of using a high-order – possibly adapted – spatial discretisation for accurate water wave propagation over long times and distances is particularly attractive for marine hydrodynamics applications.
In a number of experiments and field tests of point absorbers, snap loads have been identified to cause damage on the mooring cables. Snap loads are basically propagating shock waves, which require special care in the numerical modeling of the mooring cable dynamics. In this paper we present a mooring cable model based on a conservative formulation, discretized using the Runge-Kutta discontinuous Galerkin method. The numerical model is thus well suited for correctly capturing snap loads. The numerical model is verified and validated using analytic and experimental data and the computed results are satisfactory.
Mooring systems are required to keep floating wave energy converters (WECs) on station. The mooring concept might impact the performance of the WEC, its cost and its integrity. With the aim of clarifying the pros and cons of different mooring designs, we present the results from physical model experiments of three different mooring concepts in regular and irregular waves, including operational and survival conditions. The parameters investigated are the tension in the cables, the motions of the device in the different degrees of freedom and the seabed footprint in each case. We can see that the mooring system affects the performance of the wave energy converter, but the magnitude of the impact depends on the parameter analysed, on the mode of motion studied and on the conditions of the sea. Moreover, different configurations have similar performances in some situations and the choice of one over another might come down to factors such as the type of soil of the seabed, the spacing desired between devices, or environmental impacts. The results of our experiments provide information for a better selection of the mooring system for a wave energy converter when several constraints are taken into account (power production, maximum displacements, extreme tensions, etc).
Floating wave energy converters (WECs) operating in the resonance region are strongly affected by non-linearities arising from the interaction between the waves, the WEC motion and the mooring restraints. To compute the restrained WEC motion thus requires a method which readily accounts for these effects. This paper presents a method for coupled mooring analysis using a two-phase Navier-Stokes (VOF-RANS) model and a high-order finite element model of mooring cables. The method is validated against experimental measurements of a cylindrical buoy in regular waves, slack-moored with three catenary mooring cables. There is overall a good agreement between experimental and computational results with respect to buoy motions and mooring forces. Most importantly, the coupled numerical model accurately recreates the strong wave height dependence of the response amplitude operators seen in the experiments.