We present the results of a numerical model which has been developed for estimating the contribution to the methane slip from different sources in a four-stroke dual-fuel marine engine running on natural gas. The model is a thermodynamic three-zone zero-dimensional full engine cycle model and considers methane slip contributions from short-circuiting, crevices and wall quenching. The model is applied to analyze the methane slip from a four-stroke dual-fuel medium speed marine engine using natural gas as primary fuel. At low loads, wall quenching is found to be the dominant contribution to the methane slip. At full load, the wall quenching contribution is comparable to the level of the short-circuiting and crevice contributions which only vary relatively little with load. At 75% load, the contribution from short-circuiting is highest. In addition, we found that in-cylinder post-oxidation of unburned fuel remaining after the main combustion is negligible.
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 numerical model (MOODY) for the study of the dynamics of cables is presented in Palm et al. (2013), which was developed for the design of mooring systems for floating wave energy converters. But how does it behave when it is employed together with the tools used to model floating bodies? To answer this question, MOODY was coupled to a linear potential theory code and to a computational fluid dynamics code (OpenFOAM), to model small scale experiments with a moored buoy in linear waves. The experiments are well reproduced in the simulations, with the exception of second order effects when linear potential theory is used and of the small overestimation of the surge drift when computational fluid dynamics is used. The results suggest that MOODY can be used to successfully model moored floating wave energy converters.
This thesis presents the numerical study of combustion under marine engine like condi- tions. The thesis is divided into two main parts. In the first part, combustion is studied in a large two-stroke marine engine with conventional diesel fuel. In the second part, two different dual-fuel combustion modes of diesel-methane i.e. non-premixed combus- tion and premixed combustion are studied. All numerical models are validated with the experimental data.
First, in a simplified geometry of the marine engine, conjugate heat transfer (CHT) calculations are applied to simultaneously solve the in-cylinder gas phase dynamics and the temperature field within the liner of the engine. The effects of different initial tem- peratures across the liner and the effects of the amount of water vapor in the air on the sulfuric acid formation and condensation in a large two-stroke marine engine are studied. An initial temperature is calculated based on heat transfer modeling and it is observed that the sulfuric acid vapor formation is more sensitive to the variation of the water vapor amount than the sulfuric acid condensation. In the next step, the effects of the turbulence modelling on the simulation of the full cycle of the engine including scavenging process, combustion, and emission formation is studied in a real geometry of a large two-stroke marine research engine. The Unsteady Reynolds Averaged Navier- Stokes (URANS) and Large Eddy Simulation (LES) turbulence models are utilized for modeling of in-cylinder turbulent flow. The accuracy of the tangential velocity and swirl flow in the top of the cylinder where the fuel is injected is crucial to predict the air-fuel mixing correctly. It is found that URANS predicts a solid body rotation for the tangential velocity in this region. However, LES predicts a tangential velocity that is uniformly distributed in the radial direction that is consistent with experimental results. Furthermore, during the scavenging process, LES is able to predict the Burgers vortex upstream of the cylinder near the scavenging ports. Also, LES predicts a higher angu- lar momentum inside the cylinder in comparison with URANS. During the combustion process, the LES model shows a moderately better performance in capturing the experi- mental pressure and heat release rate profiles than URANS. However, the predicted gas temperature at the liner wall is approximately 45 % higher for URANS than LES during the expansion stroke, which is attributed to a higher predicted turbulent viscosity in the URANS case. A higher temperature of gas beside the liner wall may decrease the sulfuric acid formation and increase the heat transfer. The higher predicted swirl by LES than that in URANS leads to an earlier and stronger interaction between the flame and the spray, decreasing the oxidation of the emissions. The second cycle LES simula- tion shows that the solutions after the scavenging process are in-sensitive to the initial conditions and the main governing parameters are boundary conditions and injection characteristics.
Next, two different dual-fuel combustion modes of non-premixed combustion and premixed combustion are studied. The non-premixed combustion is simulated and val- idated with the experimental data of a large two-stroke marine research engine under low and high engine loads. Based on the results, further methane jets penetration in the low load case leads to better air-fuel mixing and a higher combustion intensity than that in the high load. Effects of the pilot diesel fuel injection timing on combustion and emission formation and the governing mechanisms are also investigated in detail. Results indicate that the intense combustion of the accumulated methane expands the methane flame towards the piston when the pilot injection timing is retarded. The NO formation is lower in the high load case due to the lower combustion intensity. Also, retarding the pilot injection timing decreases the NO formation. Furthermore, the effect of the direction of pilot diesel injection is investigated which shows a significant effect on the methane start of combustion and intensity as well as flame propagation direction which leads to different heat transfer trends from the combustion chamber walls.
Premixed combustion is analysed in a constant volume combustion chamber (CVCC) and validated with experimental data. Results show that by simulation of methane-air mixing, the numerical model is able to capture the ignition delay time (IDT) within a maximum relative difference of 7 % to the measurements. A higher relative difference of 38% is obtained when methane gas injection is omitted and the methane-air and temperature are assumed homogeneous. Therefore, it is concluded that the simulation of methane-air mixing process is crucial in this type of combustion due to the presence of inhomogeneities in both methane fuel and temperature distribution after mixing. Creating the idealized inhomogeneities for separately investigation of methane and tem- perature inhomogeneities shows that the inhomogeneity in the temperature has a more profound influence on the IDT than the inhomogeneity in the methane distribution. Furthermore, the effects of the number of pilot fuel nozzle holes on the auto-ignition are studied. The auto-ignition process in two cases with 4 nozzle holes is investigated and compared with the base case with 8 nozzle holes. Considering the same amount of pilot fuel, the injection rate is assumed to be double in one of the cases, while in the other case, the injection duration is doubled. Results show that a reduction of the nozzle hole numbers can improve the pilot diesel ignition in the case with 4 nozzle holes and double injection duration compared to the base case with 8 nozzle holes. However, combustion in the case with 4 nozzle holes and a double injection rate is incomplete due to flame impingement on the walls.
Ice-breaking cones are commonly used in the design of marine structures in cold regions. This study investigates the effects of higher-harmonic wave loads and wave runup on a 5-MW offshore wind turbine with and without ice-breaking cones under extreme wave conditions on the Liaodong Peninsula in China. Two ice-breaking cones (upward-downward and inverted types) are considered. The numerical model adopts a two-phase flow by solving unsteady Reynolds-averaged Navier-Stokes (URANS) equations using the volume of fluid (VOF) method. A phase decomposition method through a ‘Stokes-like’ formulation was adopted to obtain the parameters for each harmonics. The presence of the conical part is seen to increase the second-harmonic wave loads by up to 40%, but it has only limited influence on the fourth and fifth harmonics. The upward-downward-type ice-breaking cone increases the third harmonic, while the inverted-type ice-breaking cone decreases the third harmonic. Due to the phase difference between the first-harmonic and higher harmonics, the largest wave runup occurs at 0°, and 135° is the location with the smallest wave runup. This is because at the 135-degree location, the linear component is positive but the other nonlinear components are negative. For the 0-degree location, all harmonics are positive. By contrast, the inverted type has little effect. The high harmonic wave runup of the minimum point is backwards compared with that of the monopile, and most nonlinear wave runups are different upstream of the monopile.
Traditionally, most ship hulls are optimized for ideal conditions, where the ships are sailing in calm water with full speed in full load. In the last decade, some ships have been designed for a range of draughts and speeds in calm water. However, there is still a large gap between the ideal conditions the ships are designed for and conditions (waves, wind, currents, hull roughness ets.) the ships will operate in. The target for the thesis is to develop accurate numerical models that can help ship designers narrow a part of this gap.
The main body of this thesis is three papers. The first papers compares the speed/power performance of full-scale CFD simulations, towing tank predictions, and high quality speed trial measurements from six sister vessels. Much research have been conducted comparing model- scale CFD with towing tank results. However, very few studies have compared full-scale CFD with speed trial measurements. The study includes both a ro-ro vessel and a general cargo vessel. The present study finds that including the hull and propeller roughness directly into the CFD simulations by modifying the wall-functions provides more accurate results than the traditional approach of estimating the effect of roughness using an empirical formula.
Today, most ships are designed for sailing in calm water. However, very few ships sail entirely in calm water. Before numerical simulations can be used to predict added resistance in waves and seakeeping responses, a systematic verification and validation is required to ensure the accuracy. The second paper presents such a systematic verification and validation for the KCS container ship in oblique waves. Five wave headings and six wavelengths are studied. The estimated spatial and temporal discretization errors are found by an extensive verification study to be less than 5 %. Results from the verified CFD model are compared with existing potential flow and CFD results from the literature, as well as up to three experimental data sets. The comparison shows that the present CFD results in general show significantly better agreement with the experiments than previously published CFD results.
This CFD set-up is used in the third paper to study how sailing in oblique regular waves influences the nominal wake field of the KCS ship. Five different headings are studied and the waves have a steepness of 1/60 and a wave length equal to the ship length. The present study finds that the studied incident waves make the nominal wake field highly transient. Especially the transient bilge vortex and shadow from the skeg have a significant influence on the nominal wake field. The results show that the nominal wake fraction fluctuates up to 39 % of the mean nominal wake fraction for the studied waves. The mean nominal wake fraction is higher than in calm water for all headings besides head sea waves. It is found that the stern quartering sea waves has the maximum mean nominal wake fraction, with a 16 % higher mean nominal wake fraction than in calm water. Finally the study finds that the modified advance angle on the r/R = 0.7 circle in the propeller plane varies 3.5 degrees more in stern quartering than in calm water. This increases the risk of cavitation leading to potential vibrations and loss of propulsive efficiency.
The three papers show that CFD simulations can deliver highly accuracy results, when the CFD simulations are set-up very carefully and systematic verification and validation are conducted. The results from the three papers shows that numerical simulations have a massive potential as useful tools when designing ships for the conditions, the ship will operate in.
Numerical tests are performed to investigate wave transformations of nonlinear nonbreaking regular waves with normal incidence to the shore in decreasing and increasing water depth. The wave height transformation (shoaling) of nonlinear waves can, just as for linear waves, be described by conservation of the mechanical energy flux. The numerical tests show that the mechanical energy flux for nonlinear waves on sloping foreshores is well described by stream function wave theory for horizontal foreshore. Thus, this theory can be used to estimate the shoaled wave height. Furthermore, the amplitude and the celerity of the wave components of nonlinear waves on mildly sloping foreshores can also be predicted with the stream function wave theory. The tests also show that waves propagating to deeper water (de-shoaling) on a very gentle foreshore with a slope of cot(β) = 1200 can be described in the same way as shoaling waves. For de-shoaling on steeper foreshores, free waves are released leading to waves that are not of constant form and thus cannot be modelled by the proposed approach.
The power output from many wave energy converters (WECs) is limited by a finite stroke length in the power take-off (PTO) mechanism. As the PTO approaches its maximum stroke length, an end-stop system needs to be engaged to avoid damage to the machinery. Still the on-set of the end-stop is a nonlinear trigger force, a stiff point in the system. In this respect it is similar to how snap loads in the mooring cables affect the system after a period of cable slack. This paper presents a detailed study into the dynamics of end-stop events and snap loads for a WEC. The WEC is a bottom-mounted linear generator connected to a surface buoy via a steel wire. By comparing a linear spring model with three dynamic mooring line models we conclude that large differences are observed in the low-tension and slack regions of the cable during moderate wave loads, while minor differences are seen in the estimated peak tension. By further varying end-stop parameters we observe that the peak tension in the line changes mildly with the axial stiffness for moderate wave heights. The peak tension is surprisingly unaffected by the introduction of a critical damping level to the end-stop system, despite the significant increase in end-stop force which causes the translator to come to a sudden stop. We discuss how the connection between maximum line force and end-stop parameters is highly dependent on the buoy position in the wave at the instant of end-stop onset.
The present paper describes the work carried out in the project ’Mooring Solutions for Large Wave Energy Converters’, which is a Danish research project carried out in a period of three years from September 2014, with the aim of reducing cost of the moorings for four wave energy converters and improving the applied design procedure. The paper presents the initial layouts and costs and illustrates which solutions could potentially reduce cost. Different methods for analysis of the systems were applied, ranging from simple quasi-static analysis to full dynamic analysis and experimental work. The numerical methods were compared to the experimental data, and results showed significant underestimation of tensions in the quasi-static model while reasonable overestimation was found in the dynamic analysis even without major tuning of the model. The dynamic analysis has then been implemented in a meta-model based optimization process with the aim of optimizing the mooring layout for each WEC according to cost of the systems.
Computational fluid dynamics (CFD) is becoming an increasingly popular tool in the wave energy sector, and over the last five years we have seen many studies using CFD. While the focus of the CFD studies have been on the validation phase, comparing numerically obtained results against experimental tests, the uncertainties associated with the numerical solution has so far been more or less overlooked. There is a need to increase the reliability of the numerical solutions in order to perform simulation based optimization at early stages of development. In this paper we introduce a well-established verification and validation (V&V) technique. We focus on the solution verification stage and how to estimate spatial discretization errors for simulations where no exact solutions are available. The technique is applied to the cases of a 2D heaving box and a 3D moored cylinder. The uncertainties are typically acceptable with a few percent for the 2D box, while the 3D cylinder case show double digit uncertainties. The uncertainties are discussed with regard to physical features of the flow and numerical techniques.