Multiscale Simulation of Breaking Wave Impacts

Abstract

Formålet med dette projekt er at lave en nøjagtig, robust, geometrisk fleksibel og effektiv beregningsmodel til beregning af kræfter påstrukturer fra ikke-lineære brydende bølger. Et bedre estimat af tryk og laster påhavvindøller, tidevand - og bølgeenergi-anlæg kan bruges til forbedring af det strukturelle design.En strømningsmodel uden kompressibilitet/sammentrykkelighed og viskositet med fri overflade er præsenteret, og de styrende ligninger er formuleret i et arbitrært Lagransk- Eulersk bevægeligt koordinatsystem (ALE). Runge-Kutta-metoden er anvendt til tidsintegration, og massebevarelse er tilfredsstillet med en prædiktion-korrektionsberegning af trykket. Weighted least squares metoden (WLS) er kombineret med approksimative Riemann-løsere for at udglatte løsningen tæt ved stejle gradienter i hastigheds- og trykfelterne. Poissons ligning er løst sammen med tryk-randbetingelser ved hjælp af en generaliseret finite-pointset metode (GFPM). Dette giver en geometrisk fleksibel- og stabil løsning af væsketrykket. Den numeriske approksimation af ligningerne er udført pået ustruktureret punktsæt, og løsningerne for hastigheder og tryk er repræsenteret med multivariable polynomier. De lokale punktsæt til anvendelse i WLS- og GFPM-metoderne er fundet ved hjælp af brede først søgning i en modifiseret Delaunay graf. Grafen er opdateret i forhold til udviklingen af væskedomænet. Efter hvert Runge-Kutta skridt bliver væskedomænet tjekket for ændringer i topologi, forårsaget af væskens interaktion med en struktur eller hvis en fri overflade interagerer med andre frie overflader. Væske domænets rande bliver opdateret i forhold til topologi ændringerne i væsken. Beregningspunkternes bevægelse er næsten Lagrangsk, og det forørsager en dårlig konditionering af den generaliserede Vandermonde matrix i WLS- og GFPM-metoderne. For at undgådette gøres punktsættet finere eller grovere med en afstand til en nabo baseret på adaptivitets-metode og gendistribution med en punkt position filtrerings metode.ALE-WLS-modellen uden kompressibilitet/sammentrykkelighed og viskositet er anvendt påstandard tests: Deformation af en ellipseformet draa be, stående bålger med lille amplitude og dæmningsbrudv. Testen med deformation af en ellipseformet dråbe viser, at modellen er i stand til at beregne kinematikken og dynamikken af en strømning med fri overflade både nøjagtigt og robust. Testen med en bålgen med lille amplitude giver en lignende konklusion. En lang tidsintegration af bølger med lav amplitude er mulig, fordi bevægelsen af den fri overflade er beregnet nøjagtigt, og fordi fejlen i volumenberegningen er lille. Testen med dæmningsbrud viser, at modellen kan beregne ikke-lineære væskebevægelser, væske-struktur interaktion og brydende bølger. Udbredelseshastigheden af væskefronten og trykket påvæggen er sammenlignet med eksperimenter. Sammenligningen er rimelig god.ALE-WLS-modellen uden kompressibilitet/sammentrykkelighed og viskositet er koblet med en potential strømningsmodel, Engsig-Karup et al. [2009], for at lave multiskala beregning af brydene boelger påen lodret mole. Potential strømningsmodellen giver en nøjagtig beregning af bølgeudbredelsen fra åbent hav til tæt påmolen. Bølgebrydningen tæt ved molen, og når bølgen rammer molen, er beregnet med ALE-WLS-modellen. De største bølgelaster påden lodrette mole er ˜ 1-2 gange størrelsen af bølgelasterne beregnet ved hjælp af en kompressibel model, Bredmose et al. [2009].De vigtigste dele af dette projekt er ALE-WLS-modellen kombineret med approksimative Riemann-løsere og generaliseringen af finite-pointset-metoden til højere ordens nøjagtighed. De lokale punkt sæt til WLS-og GFPM-beregningerne er fundet med en BFS-datastruktur, som opdateres under hensyntagen til topologiændringer i væskedomænet. Denne model giver mulighed for nøjagtig multiskala-beregning af laster fra ikke-lineære brydende bølger påkonstruktioner ved hjælp af koblingen med en potential strømnings model.The purpose of this project is to make an accurate, robust, geometric flexible and efficient model for calculation of forces on structures from nonlinear ocean waves and breaking wave impacts. Accurate prediction of the extreme forces on wind turbine foundations, breakwaters and tidal or wave power devises are important for enhancement structural designs.The proposed model is based on an incompressible and inviscid flow approximation and the governing equations are applied in an arbitrary Lagrangian-Eulerian moving frame of reference (ALE). The Runge-Kutta method (RK) is used for time integration and mass conservation is satisfied through a pressure-corrector type calculation of the pressure. The weighted least squares method (WLS) is combined with approximate Riemann solvers to introduce numerical smoothing of the solution around steep gradients in the velocity and pressure fields. The Poisson equation is solved and the pressure boundary conditions are satisfied by a generalized finite pointset method (GFPM); This provides a geometrically flexible and stable solution for the fluid pressure. The numerical approximations of these equations are performed on unstructured point distributions and the solutions for velocity and pressure are represented by WLS approximation of multivariate polynomials. The stencils for the ALE-WLS and GFPM methods are found through a breadth first search (BFS) in a modified Delaunay graph. This graph is the discrete representation of the fluid domain and the connectivity between the calculation points. The graph is updated according to the evolving topology of the fluid domain caused by the fluid reaching or leaving a solid boundary or the free surface colliding with itself or another free-surface. After each time step the fluid domain is checked for any of these intersections and the topology is updated accordingly in its graph representation. The calculation points move in a Lagrangian way and this can cause ill-conditioning of the generalized Vandermonde matrix in the WLS and GFPM methods. To prevent this the point set is refined and coarsened by a fill-distance based adaptivity method and redistributed via a point position filtering method. The incompressible and inviscid ALE-WLS model is applied to the following standard validation test cases: deforming elliptical drop, small amplitude standing waves and the dam break problem. The deforming elliptical drop test show that the model can calculate the kinematics and dynamics of this free surface flow accurately and robustly. The small amplitude standing wave gives the same conclusions. Long time integration of this small amplitude periodic motion is possible due to accurate free surface evolution and small errors in the fluid volume. The dam break test case shows that the incompressible and inviscid ALE-WLS model can calculate nonlinear fluid motion, fluid structure impacts and overturning waves. The propagation speed of the wetting front and impact pressures are compared to experiments and the results compare reasonably well. The incompressible and inviscid ALE-WLS model is coupled with the potential flow model of Engsig-Karup et al. [2009], to perform multiscale calculation of breaking wave impacts on a vertical breakwater. The potential flow model provides accurate calculation of the wave transformation from offshore to the vicinity of the breakwater. The wave breaking close to the breakwater and the wave impact are calculated by the incompressible ALE-WLS model. The forces calculated with the incompressible and inviscid ALE-WLS model are ˜ 1 - 2 times the corresponding compressible calculations in Bredmose et al. [2009] for the calculations without trapped air.Among the contributions of this project are the ALE-WLS method combined with approximate Riemann solvers and the generalization of the FPM method to arbitrary order of accuracy. The WLS and GFPM stencils found using the BFS data structure, which is updated due to topology changes of the evolving fluid domain. This extension combined with ALE-WLS and approximate Riemann solvers gives a numerical model capable of calculation of forces due to breaking wave impacts. The incompressible and inviscid ALE-WLS model has been coupled with a potential flow model to provide multiscale calculation of forces from breaking wave impacts on structures