1 / 38

WISE 2008, Helsinki, Finland

WISE 2008, Helsinki, Finland. On perspectives for the numerical Solution of the WAE: Experiments and Decomposition Methods Aron Roland and Jürgen Geiser Technische Universität Darmstadt Humboldt Universität Berlin. Overview of the talk.

cleary
Download Presentation

WISE 2008, Helsinki, Finland

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. WISE 2008, Helsinki, Finland On perspectives for the numerical Solution of the WAE: Experiments and Decomposition Methods Aron Roland and Jürgen Geiser Technische Universität Darmstadt Humboldt Universität Berlin

  2. Overview of the talk • Available numerical schemes for the solution of the WAE • Operator Splitting (OSM) and Iterative Direct (IDM) Methods • Splitting error demonstrated on simple cases • Some comparison between IDM (SWAN) and OSM (WWM) • Perspectives for the solution of the WAE on unstructured meshes using OSM • Splitting errors of classical methods (Splitting Errors) • Iterative Operator Splitting methods (Skip Splitting Errors) • Benefits of iterative methods (Efficiency, Simplicity) • Adaptive time-decomposition using the eigenvalues

  3. The Wave Action Equation (WAE)

  4. Numerical schemesIterative Direct Methods (IDM) • Implicit Iterative Direct Methods (SWAN) • Discretization of the whole equation as suggested by Patankar (1980) see Booij et al. (1999) for the WAE

  5. Numerical schemesOperator Splitting Methods • Operator Splitting Methods (OSM) e.g.: WAM, MIKE21, TOMAWAC, CREST … • 1st Step – Advection • 2nd Step – Source Terms

  6. Numerical schemesOperator Splitting Methods II • Operator Splitting Methods (OSM) e.g. WWIII or WWM • 1st Step – Spectral part • 2nd Step – Geographical space • 3rd Step – Integration of the source terms

  7. Numerical schemesOperator Splitting Methods III • Operator Splitting Methods and Iterative Sources (OSM-IS); WWM • 1st Step – Spectral part • 2nd Step – Geographical space and sources

  8. Two important mathematical theorems • The Godunov theorem • Linear schemes (linear combination of the unknown variables) cannot be monotone and higher order! • The Lax theorem • The numerical scheme is convergent if it can be proven that it is stable and consistent. • Consequences: • Implicit schemes must be nonlinear to achieve higher order accuracy while retaining monotonicity. This leads directly to nonlinear equation systems that have to be solved • Explicit methods can be easily defined to be nonlinear and achieve in this way higher order accuracy though retaining monotonicity. • The convergence of OSM methods is much easier to prove as for the direct methods since the chosen scheme for the solution of the sub-problems are in most cases shown to be consistent and stable by the authors itself.

  9. Numerical schemes in the WWM II • Numerical methods for the sub-problems • Geographical space • Galerkin schemes (non-monotone, conservative) • Crank-Nicolson Taylor Galerkin (implicit, 2nd order in space and time) • Euler Taylor Galerkin (implicit, 2nd order in space) • Residual distribution schemes (monotone, conservative) • CRD-N1 schemes (implicit) • CRD-N3 schemes (implicit, 2nd order in time) • CRD-N scheme (explicit) • CRD-PSI scheme (explicit, better then 1st order) • CRD-FCT scheme (explicit, 2nd order in space and time) • Source term integration • Semi-implicit (WAM) • Dynamical (WWIII) • Iteratively within the implicit advection schemes in geographical space • Spectral space • Ultimate Quickest (explicit, 3rd order in space and time)

  10. Splitting Error between advection and strong local sources CFLX=14 For Δt=1.0s Left: Significant wave height along a cross section for the unsplitted solution (blue) with Δt = 0.005 compared to the splitted solution using the explicit CRD-N scheme (red) and the implicit CRD-N1 scheme (green)

  11. Splitting Error between advection and strong local sources

  12. Splitting Error Advection in Geographical and spectral space Bathymetry (left) and Computational mesh (right). Figure 38: Computed significant wave height (left) and mean wave direction (right) using the CRD-N1 explicit scheme for the Vincent & Briggs (1989) experiment.

  13. Splitting Error Advection in Geographical and spectral space Wave height (left) and Average wave direction (right). Figure 38: Computed significant wave height (left) and mean wave direction (right) using the CRD-N1 explicit scheme for the Vincent & Briggs (1989) experiment.

  14. Splitting Error Advection in Geographical and spectral space CFLX=11 CFLθ=3.2 For Δt=1.0s Figure 38: Computed significant wave height (left) and mean wave direction (right) using the CRD-N1 explicit scheme for the Vincent & Briggs (1989) experiment.

  15. Splitting Error Advection in Geographical and spectral space CFLX=11 CFLθ=3.2 For Δt=1.0s Figure 38: Computed significant wave height (left) and mean wave direction (right) using the CRD-N1 explicit scheme for the Vincent & Briggs (1989) experiment.

  16. Splitting Error Advection in Geographical and spectral space CFLX=11 CFLθ=3.2 For Δt=1.0s Figure 38: Computed significant wave height (left) and mean wave direction (right) using the CRD-N1 explicit scheme for the Vincent & Briggs (1989) experiment.

  17. Splitting Error between advection in geographical space and spectral space L2 error norm of the CNTG and ETG scheme in combination with the CN and the EI scheme for directional space respectively (left). L2 error norm for the CNTG, the CRD-N1EXP, the CRD-FCT the CRD-N1IMP, the CRD-N2IMP and the CRD-N3IMP scheme using the UQ approach for the integration in directional space (right). Figure 38: Computed significant wave height (left) and mean wave direction (right) using the CRD-N1 explicit scheme for the Vincent & Briggs (1989) experiment.

  18. Continental Shelf refraction (SWAN)

  19. Continental Shelf refraction (WWM)

  20. Continental Shelf refraction(WWM)

  21. + - • Implicit Iterative Direct Methods (SWAN) • The whole equation is discretized at once. • Lower order diffusive linear implicit schemes (BSBT) must be used in order to maintain monotonicity • Higher order schemes (e.g. S&L) are not monotone and will lead in the vicinity of strong gradients in the solution to non physical results ( • A advantage is that the convergent solution is independed of the integration time step. • One deficiency is the strong numerical diffusion- • Another deficiency is that the convergence of the whole scheme is difficult to prove (LAX) and showed in practice to be a problem of such an approach. One solution is under relaxation (e.g. Ferzinger & Peric) applied by Zijlema & Van der Westerhuysen for the SWAN model. However, under-relaxation increases usually the amount of iterations. • The limiter acts also on the advective part! This is very bad in non-stationary situations … or somebody tell me why it should be good to limit the advective part

  22. + - • Operator Splitting Methods (OSM) • Godunov theorem can be easily applied. • Lax equivalence theorem can be easily proved in most cases. • Decoupling of the certain time scales of importance and accounting for them numerically. • Optimal numerical schemes can be used for the certain sub-problems, e.g. adaptive schemes. • Source term integration can be done adaptive (Tolman) reducing the necessity of application of the limiter to the change in time. • The global scheme converges to the analytical solution when the discretization scales approaches the infitesimal limits. • Parallelization and vectorization can be done at the same time. - Multi-scale, multi-code and multi-grid can be embedded.

  23. Development of the WWM II • The Residual Distribution Framework was successfully introduced in the WWM code and verified for different cases e.g. • U.S. East Coast • Haringvliet Estuary • Continental shelf refraction (Ardhuin & Herbers) • The WWM II was coupled with SHYFEM a SWE Finite Element Model. The coupled models SHYFEM/WWM was used for the simulation of the storm surge in the northern Adriatic and the Gulf of Mexico (Ferrarin et al., 2006) • Implementation of the source terms according to Donelan et al., Babanin et al., Young et al. … • In the cue … 3-d wave – current coupling, Ardhuin et al. • Implementation of the WWM advection schemes in the WWIII by Florent, Fabrice and me (done) • Parallelization using MPI (Florent & Fabrice) is (done) • Multi-scale validation using SHOWEX, EPEL, YANGTZE, TAIWAN, NORTH SEA with Ardhuin et al. • Multi-Grid approach using P1 (3 points, 2nd order space-time) and P2(6 points; 3rd space-time) elements • Publish something … hopefully …

  24. Questions? • How is the diffuseness of the numerical schemes included in the parameterizations of the physics and how about the conclusions from numerical results with respect to real physics? • How can we interpret physics in numerical results with no or very little knowledge about the numerical schemes and their errors when we are using them? • Possible solutions for the outlined problems of OSM are discussed in the next part of the lecture.

  25. OSM in multi-physic applications Decomposition methods are a powerful method of numerical investigation of complex (physical) time-dependent models, where the stationary part (elliptic) part consists of simpler operators, e.g.: • Transport-Reaction Processes, see [Geiser, 2006; Hundsdorfer, Verwer, 2003] (Physical splitting) • Hamiltonian Systems, see [McLachlan, 1994; Hairer, Lubich, Wanner, 2002] (Symplectic splitting) • Air pollutant models, see [Zlatev, 1995] (Operator Splitting) • Maxwell equations, see [Horvath, 2006] (Operator Splitting) • Wave Action Equation (e.g. Tolman, 1992 and many others )

  26. Attributes of OSM Available decomposition methods: • Time-decomposition methods • Spatial-decomposition methods Contribution of the decomposition of the whole equation • Decoupling the time-scales, space-scales (Reduce the stiffness in single operators) • Decoupling the multi-physics. (Reduce the unphysical behavior with best choice of discretization and solver methods, e.g. CFL-conditions) • Time-adaptivity, Space-adaptivity. (Efficiency and accuracy in computations). • Parallelization in Time and Space. (Reduction of computational time). • Results : More efficient and fast algorithms with high accuracy, simple implementable.

  27. Basic Idea • Idea : Decoupling into simpler parts with respect to operators, dimensions, times, models • Effect : More efficient computations, parallelization, adequate discretization and solver methods for each part Deficiencies: • Decoupling of the global time step – unphysical results. • Order reduction in the stiff case

  28. Splitting Error and Comutativity • The X and Y operator compute on a FDM mesh and so we have no splitting error. • This is why the splitted Ultimate Quickest scheme is better then the un-splitt counterpart maintaining the mixed terms due to the 2-d discretization • Since theta and sigma are orthogonal operators there is no splitting error when splitting the intra-spectral propagation. However, care must be taken in the stiff case.

  29. Problems of the Classical Splitting Methods Error analysis • Error is of first order global in time : • Local and global error for the decomposition and the full solution • Here is the global error and the error describes the situation if non of the terms is comuting otherwise one gets exact solution (Strang, 1968)

  30. Physical Error • In each equation, we solve a different problem due to the full equation, e.g. no influence of the intermediate time-steps to the other equations. • We only interact at the initial conditions, so the time-scales of the operators A and B are independent and therefore not interacting as in the original equation, e.g. coupled transport equations [Geiser 07].

  31. Iterative Operator Splitting Methods Benefits: • Larger time-steps are possible in each iterative step • Spectrum of the operator allow to control the stiffness of the operators • Higher order can be claimed with more iterative steps • Efficient and simple implementation of the iterative schemes • Parallel algorithms can be used

  32. Application to the Wave Action Equation We solve iteratively four equations: where i = 1,2,3…, N0 = 0, Nin = Nn

  33. Analysis: Iterative Operator Splitting Methods • We concentrate on two operators: Where N0(t) is any fixed function for each iteration. (Here, as before, un denotes the known split approximation at the time level t = tn) The split approximation at the time-level t = tn+1 is defined as Nsp,n+1 = N2m+1(tn+1). (Clearly, the functions Nk(t) (k = i−1, i, i + 1) depend on the interval [tn, tn+1], too, but, for the sake of simplicity, in our notation we omit the dependence on n)

  34. Decomposition based on operator spectrum • To detect the operators in the differential equation as stiff or non-stiff operators, we can apply the Eigen values of each operator and use them as reciprocal time scales. The operator equations are analyzed with the Eigen value problem: • where the operators A and B result form the spatial discretization. • The Eigen values are detected in the decoupled equations: • Based on the Eigen values λA, λB we can propose the time steps ΔtA=1/ λA and ΔtB=1/ λB

  35. Decomposition MethodsError for the Iterative splitting-method • Theorem: The error for the splitting methods is given as : where n is the time-step, e0 the initial error e0(t) = N(t) − N0(t) and m the number of iteration-steps, K and Km are constants, ||B|| is the maximum norm of operator B and A, B are bounded, monotone operators. The error can be controlled by the operator B, e.g. non stiff operator. Proof : Taylor-expansion and estimation of exp-functions. See the work Geiser,Farago (2005).

  36. The idea for parallelization in time are the windowing, that the processors has an amount of time-steps to compute and to share the end-result of the computation as an initial-condition for the next processor. Parallelization of the Time-Decomposition method: Windowing

  37. Outlook • Experiments for the Wave Action Equation • Theoretical Investigations in multiple iterative operator splitting methods • Parallel implementation

  38. Jaws, Hawaii Björn Dunkerbeck, E11 World Champion On perspectives for the numericalSolution of the WAE: Experiments andDecomposition Methods Mast length = 4.2m Dipl.-Ing. A. Roland Dr. Geiser Prof. Dr.-Ing. habil Prof. h.c. U.C.E. Zanke Prof. Dr. S-H. Ou Prof. Dr. T-W. Hsu Dr. J-M. Liau

More Related