Hydrodynamics: Setting the Scene* *a selective view Prof. Rodney Eatock Taylor and Prof. Paul Taylor University of Oxford The Lloyd s Register Educational Trust (LRET) Marine & Offshore Research Workshop 16-18 February, 2010 at Engineering Auditorium, NUS
Hydrodynamics: Setting the Scene* *a selective view Rodney Eatock Taylor and Paul Taylor University of Oxford LRET Marine and Offshore Research Workshop National University of Singapore, 16-18 February 2010
Hydrodynamics and Modelling One of our aims must be to model the interaction of the most extreme waves with structures in the sea. And we must also consider hydrodynamic models motivated by design for operational conditions in waves Also Resistance, Manoeuvring, Sloshing, VIV/VIM, Hydroelasticity,.., which are not considered here Typhoon TLP from http://blog.kir.com/archives/002473.asp
Modelling using CFD From: Summary of CESOS Workshop 2007, by Marilena Greco
Issues raised Wave analysis Wave structure interaction Near-trapping by multiple structures in waves Wave-in-deck issues (waves + interaction): jackets and GBS Wave-induced sloshing in gaps Floatover A few comments on CFD
Wave analysis Time series plot of some large waves
NewWave Average of largest 50% linear profiles * compared with NewWave profile + * linearised using (average crest -average trough profile)/2 + average shape of extreme in linear random Gaussian process (Boccotti, Lindgren)
NewWave embedded in random record
Freak waves? Consider behaviour near an extreme crest: set-down will decrease crest elevation, while set-up will increase it. The famous New Year Wave measured at the Draupner platform in the North Sea. For convenience the time is set to zero under the extreme crest.
Setdown under extreme crest? T.A.A. Adcock & P. H. Taylor, 11th International Workshop on Wave Hindcasting and Forecasting, 2009.
Setup under Draupner crest The second order difference waves under the Draupner giant wave.comparison of filtered data and predicted second order difference waves based on a model of crossing seas. (Evidence of crossing seas from meteorology, measured forces, and absence of breaking in wave record from laser) T.A.A. Adcock & P. H. Taylor, 11th International Workshop on Wave Hindcasting and Forecasting, 2009.
Other evidence of set-up in crossing seas A. Toffoli et al. / European Journal of Mechanics B/Fluids 25 (2006) 649 661
Wave climate variability T.A.A. Adcock & P. H. Taylor, 11th International Workshop on Wave Hindcasting and Forecasting, 2009.
Wave-structure interaction 6 ka = 0.4 6 ka = 1 4 2 1.5 4 2 1.5 0-2 1 0-2 1-4 0.5-4 0.5-6 -5 0 5-6 -5 0 5 6 ka = 3 ray theory approx 6 ka = 5 4 2 1.5 4 2 1.5 0-2 1 0-2 1-4 0.5-4 0.5-6 -5 0 5-6 -5 0 5 Diffracted wave field near single vertical cylinder
Diffraction by arrays, near-trapping and wave-in-deck 6 ka = 2.2729, β = 0 o 6 ka = 2.2729, β = 15 o 4 4 2 4 3 2 4 3 0-2 1 2 0-2 1 2-4 -4-6 -5 0 5-6 -5 0 5 6 4 2 0-2 -4-6 ka = 2.2729, β = 30 o 4 3 1 2-5 0 5 ka = 2.2729, β = 45 o 6 4 2 4 3 0-2 1 2-4 -6-5 2 4 6 0 5
Near-trapping in square array Runup at cylinder 3, β = 45
Influence of local geometry Dimensionless wave amplitude near arrays of circular and non-circular cylinders ka= 1.831, β= 45 (simple models of a specific semisubmersible).
Application to a GBS
2 nd order behaviour simulated using Oxford program DIFFRACT near trapping frequency doubling : input waves at frequency f strong local surface response at 2f Body and free-surface grid for quadrant of the structure (using two planes of symmetry)
Incident wave End-on waves at near-trapping. β=0 o, f=0.126 Hz
Where does the high response occur? Linear amplitude response profile along the center-line of the structure at f=0.126hz. 2 nd order sum amplitude profile along the center-line of the structure at 2f =0.126Hz. Linear diffraction calculations provide useful guidance to 2 nd order behaviour
11m Predicted surface elevation between the rear legs with 11m incident focused wave group T p = 14.3 s η (1) η (2+) η (1+2) =TOTAL Practical implications for concrete GBS with closely spaced columns - water projection to high level - hitting the deck - air down to just above the caisson (slam loading on the top?) Walker et al., Ocean Engineering 2008
Waves in gap between closely spaced vessels SAFE OFFLOAD programme (EU) Noble Denton Experiments at DHI and Imperial College have confirmed importance of large free surface motions in the gap As has second order boundary element analysis at Oxford
Linear motions and elevations, beam sea 1 Floating Fixed
Relationship between first-order and second-order results for elevation(beam sea 1) 20 0.755 point (-2,0) 7 point (-2,0) 0.75 15 LINEAR 6 5 2 nd order quadratic η (1) /A 10 0.81 0.88 η (2) q 4 3 5 0.95 2 1.04 1.11 1.20 1.29 1.35 0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 15 12 9 ω (rad/s) point (-2,0) TOTAL 2 nd order 0.76 1 0 0.3 0.4 0.5 0.6 0.7 0.8 12 10 8 ω (rad/s) point (-2,0) 2 nd order potential 0.76 η (2) 6 0.44 0.405 0.475 0.52 0.56 0.60 0.64 0.68 η (2) p 6 4 0.44 0.405 0.475 0.52 0.56 0.60 0.64 0.67 0.71 3 0.38 2 0.38 0 0.3 0.4 0.5 0.6 0.7 0.8 ω (rad/s) 0 0.3 0.4 0.5 0.6 0.7 0.8 ω (rad/s) First-order 0.755 0.81 0.88 0.95 1.04 1.11 1.20 1.29 1.35 Second-order 0.38 0.405 0.440 0.475 0.52 0.56 0.60 0.64 0.67
Behaviour in transient waves Response in gap between 2 boxes with NewWavefocused group, beam seas, broad band spectrum
Experiments on GBS and Tanker (fixed) IC model scale unit: m water depth=0.82m
First-order wave elevation at centre line 5 4 T=0.780s(beam sea 2) Numerical (Oxford) Experimental (IC) 7 6 5 T=0.889s(beam sea 2) Numerical (Oxford) Experimental (IC) 3 4 η (1) /A 2 η (1) /A 3 1 2 1 0-1.2-0.8-0.4 0.0 0.4 0.8 1.2 0-1.2-0.8-0.4 0.0 0.4 0.8 1.2 8 7 6 5 y T=1.067s(beam sea 2) Numerical (Oxford) Experimental (IC) 3.0 2.5 2.0 y T=1.455s(beam sea 2) Numerical (Oxford) Experimental (IC) η (1) /A 4 3 2 1 η (1) /A 1.5 1.0 0.5 0-1.2-0.8-0.4 0.0 0.4 0.8 1.2 0.0-1.2-0.8-0.4 0.0 0.4 0.8 1.2 y y
30 20 10 T p =1.2s incident y=0 Wave group excitation Beam sea, tanker exposed ζ 0-10 -20-30 0 10 20 30 40 50 60 70 t(s) 20 15 First-order elevations at point (-0.03,0) in Beam Sea 2 0.9901 (T=1.01s) 0.9703 (T=1.03s) Numerical(Oxford) Experimental(IC) RAO for elevation at midpoint of gap η (1) /A 10 5 1.13636 (T=0.88s) 1.2987 (T=0.77s) Overall good agreement -freq of peak ~2% out - magnitude over-estimated 0 0.4 0.6 0.8 1.0 1.2 1.4 1.6 f(hz)
Floatover operations Behaviour of flexibly connected bodies Chevron Olokola Floatover Installation Study Ocean Dynamics LLC
Relative heave motions of barges based on extension of diffraction analysis for flexibly connected vessels
Link to hydroelastic analysis of flexibly connected structures From: Wang, C. M. & Tay, Z. Y. HydroelasticAnalysis and Response of Pontoon-Type Very Large Floating Structures
Areas needing further understanding As discussed here Extreme wave statistics Directional nonlinear wave kinematics (with and without breaking) Wave structure interaction effects various, including wave-in-deck Other topics Coupled tank sloshing and vessel motions (parametric resonance, ) Impact pressures (compressibility effects, ) Hydroelasticity
Areas for further CFD development Unstructured adaptive grids for VOF modelling Overlapping grids Grid size (Re) 9/4 (for representation of disparate length scales in turbulent flow) Turbulence modelling for unsteady free surface flows LES SPH Hybrid approaches: simple modelling to identify critical parameters (wave frequency, direction, load case, geometry), followed by very detailed CFD e.g. strip theory + RANSE Coupled and multi-scale formulations Benchmark results for validation Integration of codes into optimisation tools (as has been done for nonlinear wave resistance)
Opportunities for marine CFD From: http://wiki.manchester.ac.uk/spheric
Opportunities for marine CFD Courtesy:
Opportunities for marine CFD Courtesy:
Thank you Questions?