Mathematical and computational modeling for the generation and propagation of waves in marine and coastal environments
|
|
- Clementine Reynolds
- 6 years ago
- Views:
Transcription
1 Mathematical and computational modeling for the generation and propagation of waves in marine and coastal environments Maria Kazolea School of Environmental Engineering Technical University of Crete September 27, 213 MINISTRY OF EDUCATION & RELIGIOUS AFFAIRS, CULTURE & SPORTS
2 Motivation Discretize depth-integrated equations that model free surface flows, using mass and momentum conservation, by (unstructured) FV schemes. Most popular (applied): Nonlinear Shallow Water Equations (SWE)
3 Motivation Discretize depth-integrated equations that model free surface flows, using mass and momentum conservation, by (unstructured) FV schemes. Most popular (applied): Nonlinear Shallow Water Equations (SWE) Limitation: Not applicable for wave propagation in intermediate / deeper waters (dispersion has an effect on free surface flow) Popular Boussinesq-type (BT) models for intermediate water depths: Peregrine s standard equations (1967) but for h L 1 5
4 Motivation Discretize depth-integrated equations that model free surface flows, using mass and momentum conservation, by (unstructured) FV schemes. Most popular (applied): Nonlinear Shallow Water Equations (SWE) Limitation: Not applicable for wave propagation in intermediate / deeper waters (dispersion has an effect on free surface flow) Popular Boussinesq-type (BT) models for intermediate water depths: Peregrine s standard equations (1967) but for h 1 L 5 Madsen and Sörensen s (MS) equations (1992) Nowgu s equations (1993) for Beji and Nadaoka (BN) equations (1996) h L 1 2
5 Motivation Discretize depth-integrated equations that model free surface flows, using mass and momentum conservation, by (unstructured) FV schemes. Most popular (applied): Nonlinear Shallow Water Equations (SWE) Limitation: Not applicable for wave propagation in intermediate / deeper waters (dispersion has an effect on free surface flow) Popular Boussinesq-type (BT) models for intermediate water depths: Peregrine s standard equations (1967) but for h 1 L 5 Madsen and Sörensen s (MS) equations (1992) h Nowgu s equations (1993) for 1 L 2 Beji and Nadaoka (BN) equations (1996) Gobbi, Kirby and Wei BT model (2) Variety of BT models that include higher-order nonlinear and dispersive terms: P.A. Madsen et al. (22-29), Lynett et al. (24-21).
6 (Some) Numerical Works Until recently, Finite-Differences (FD) was the predominant method based on the work of Wei & Kirby (1995), for 1D & 2D computations
7 (Some) Numerical Works Until recently, Finite-Differences (FD) was the predominant method based on the work of Wei & Kirby (1995), for 1D & 2D computations More recently hybrid FV/FD schemes in 1D & 2D: Nwogu s, MS, BN, Serre Green-Green Naghdi S-GN equations using: Riemann solvers (Roe s and HLL-type), MUSCL-type reconstructions. Erduran et al. (25 & 27), Cienfuegos et al (26 & 27), Tonelli & Petti (29), Shiach & Mingham (29), Roeber et al. (21), Bonneton et al. (21), Kazolea & Delis (211), Dutykh et al. (211).
8 (Some) Numerical Works Until recently, Finite-Differences (FD) was the predominant method based on the work of Wei & Kirby (1995), for 1D & 2D computations More recently hybrid FV/FD schemes in 1D & 2D: Nwogu s, MS, BN, Serre Green-Green Naghdi S-GN equations using: Riemann solvers (Roe s and HLL-type), MUSCL-type reconstructions. Erduran et al. (25 & 27), Cienfuegos et al (26 & 27), Tonelli & Petti (29), Shiach & Mingham (29), Roeber et al. (21), Bonneton et al. (21), Kazolea & Delis (211), Dutykh et al. (211). MS equations in 2D Tonelli & Petti (29 & 21), two-layer BT equations Lynnet et al., (26-21), TVD Boussinesq solver Shi, Kirby et. al. (212).
9 (Some) Numerical Works Until recently, Finite-Differences (FD) was the predominant method based on the work of Wei & Kirby (1995), for 1D & 2D computations More recently hybrid FV/FD schemes in 1D & 2D: Nwogu s, MS, BN, Serre Green-Green Naghdi S-GN equations using: Riemann solvers (Roe s and HLL-type), MUSCL-type reconstructions. Erduran et al. (25 & 27), Cienfuegos et al (26 & 27), Tonelli & Petti (29), Shiach & Mingham (29), Roeber et al. (21), Bonneton et al. (21), Kazolea & Delis (211), Dutykh et al. (211). MS equations in 2D Tonelli & Petti (29 & 21), two-layer BT equations Lynnet et al., (26-21), TVD Boussinesq solver Shi, Kirby et. al. (212). 2D Finite Element (FE) on unstructured meshes: Walkey & Berzins (22), Sorensen et al. (24), Escilsson & Sherwin (26) and Engsig-Karup et al. (28).
10 (Some) Numerical Works Until recently, Finite-Differences (FD) was the predominant method based on the work of Wei & Kirby (1995), for 1D & 2D computations More recently hybrid FV/FD schemes in 1D & 2D: Nwogu s, MS, BN, Serre Green-Green Naghdi S-GN equations using: Riemann solvers (Roe s and HLL-type), MUSCL-type reconstructions. Erduran et al. (25 & 27), Cienfuegos et al (26 & 27), Tonelli & Petti (29), Shiach & Mingham (29), Roeber et al. (21), Bonneton et al. (21), Kazolea & Delis (211), Dutykh et al. (211). MS equations in 2D Tonelli & Petti (29 & 21), two-layer BT equations Lynnet et al., (26-21), TVD Boussinesq solver Shi, Kirby et. al. (212). 2D Finite Element (FE) on unstructured meshes: Walkey & Berzins (22), Sorensen et al. (24), Escilsson & Sherwin (26) and Engsig-Karup et al. (28). FV for unstructured meshes: Only one work by Asmar and Nwogu (26) using a low-order staggered scheme
11 Outline Physical problem set up
12 Outline Physical problem set up Discretization of Nwogu s and Madsen and Sørensen s extened BT models in 1D Numerical results in 1D
13 Outline Physical problem set up Discretization of Nwogu s and Madsen and Sørensen s extened BT models in 1D Numerical results in 1D Comparison of two 2D FV schemes on triangles
14 Outline Physical problem set up Discretization of Nwogu s and Madsen and Sørensen s extened BT models in 1D Numerical results in 1D Comparison of two 2D FV schemes on triangles A new unstructured FV scheme for BT equations Numerical results in 2D
15 Outline Physical problem set up Discretization of Nwogu s and Madsen and Sørensen s extened BT models in 1D Numerical results in 1D Comparison of two 2D FV schemes on triangles A new unstructured FV scheme for BT equations Numerical results in 2D Conclusions
16 Physical problem setup η: free surface elevation; h: steel water level; H = η + h: total water depth; b: bottom topography; L: wave length; A: wave amplitude Intermediate water: Deep water: h L > < h L 1 2 Shallow water: h L 1 2
17 Physical problem setup η: free surface elevation; h: steel water level; H = η + h: total water depth; b: bottom topography; L: wave length; A: wave amplitude Intermediate water: Deep water: h L > < h L 1 2 Shallow water: h L 1 2
18 Mathematical models: Nwogu s equations Vector conservative form (for both models): U t + F(U ) x = S(U), [ ] [ H Hu U = P, F(U) = Hu gh2 ] [, U = H Hu ]. (1) Using z a =.53753h as optimal reference depth (Roeber et al., 21).
19 Mathematical models: Nwogu s equations Vector conservative form (for both models): U t + F(U ) x = S(U), [ ] [ H Hu U = P, F(U) = Hu gh2 ] [, U = H Hu ]. (1) Using z a =.53753h as optimal reference depth (Roeber et al., 21). P = Hu + Hz a ( za 2 u xx + (hu) xx ) [ Velocity function]
20 Mathematical models: Nwogu s equations Vector conservative form (for both models): U t + F(U ) x = S(U), [ ] [ H Hu U = P, F(U) = Hu gh2 ] [, U = H Hu ]. (1) Using z a =.53753h as optimal reference depth (Roeber et al., 21). P = Hu + Hz a ( za 2 u xx + (hu) xx ) [ Velocity function] S(U) = S b + S f + S d [Source terms]
21 Mathematical models: Nwogu s equations Vector conservative form (for both models): U t + F(U ) x = S(U), [ ] [ H Hu U = P, F(U) = Hu gh2 ] [, U = H Hu ]. (1) Using z a =.53753h as optimal reference depth (Roeber et al., 21). P = Hu + Hz a ( za 2 u xx + (hu) xx ) [ Velocity function] S(U) = S b + S f + S d [Source terms] S b = [ ghb x ] T, S f = [ ghs f ], S d = [ ψ C uψ C + ψ M R b ] u u S f = nm 2 Friction force, n H 4/3 m = Manning coeff.
22 Mathematical models: Nwogu s equations Vector conservative form (for both models): U t + F(U ) x = S(U), [ ] [ H Hu U = P, F(U) = Hu gh2 ] [, U = H Hu ]. (1) Using z a =.53753h as optimal reference depth (Roeber et al., 21). P = Hu + Hz a ( za 2 u xx + (hu) xx ) [ Velocity function] S(U) = S b + S f + S d [Source terms] S b = [ ghb x ] T, S f = [ ghs f ], S d = [ ψ C uψ C + ψ M R b ] u u S f = nm 2 Friction force, n H 4/3 m = Manning coeff. ( ψ M = H t z za ) ( ) ( ) a 2 u z xx + (hu) xx, a ψc = 2 h2 hu 2 6 xx + z a + h h(hu) 2 xx x R b parametrization of wave breaking characteristics
23 Mathematical models: Madsen & Sørensen s equations P = Hu (B )h2 (Hu) xx 1 3 h x(hu) x [ Velocity function]
24 Mathematical models: Madsen & Sørensen s equations P = Hu (B )h2 (Hu) xx 1 3 h x(hu) x [ Velocity function] S(U) = S b + S f + S d [Source term] where now S d = [ ψ R b ] and
25 Mathematical models: Madsen & Sørensen s equations P = Hu (B )h2 (Hu) xx 1 3 h x(hu) x [ Velocity function] S(U) = S b + S f + S d [Source term] where now S d = [ ψ R b ] and ψ = Bgh 3 η xxx 2h 2 h x Bgη xx B = 1 determines the dispersion properties of the system. 15
26 Numerical Model Advective part and topography source: Well-balanced FV formulation Dispersive terms: Finite differences.
27 Numerical Model Advective part and topography source: Well-balanced FV formulation Dispersive terms: Finite differences. Roe s Riemann solver is used (Roe, 1981).
28 Numerical Model Advective part and topography source: Well-balanced FV formulation Dispersive terms: Finite differences. Roe s Riemann solver is used (Roe, 1981). Upwinding of the topography source term (Bermudez et al., 1994, Delis et al., 28).
29 Numerical Model Advective part and topography source: Well-balanced FV formulation Dispersive terms: Finite differences. Roe s Riemann solver is used (Roe, 1981). Upwinding of the topography source term (Bermudez et al., 1994, Delis et al., 28). High-order spatial accuracy: fourth order MUSCL-type scheme (Yamamoto et al., 1998).
30 Numerical Model Advective part and topography source: Well-balanced FV formulation Dispersive terms: Finite differences. Roe s Riemann solver is used (Roe, 1981). Upwinding of the topography source term (Bermudez et al., 1994, Delis et al., 28). High-order spatial accuracy: fourth order MUSCL-type scheme (Yamamoto et al., 1998). For dispersive terms: fourth order FD of first-order spatial derivatives, second and third-order FD for second and third-order derivatives. Satisfy the C-property (flow at rest) to higher spatial order: Addition of an extra term to bed upwinding (Hubbard and Garcia-Navarro, 2 and Delis and Nikolos, 29)
31 Numerical Model (cont.) Special treatment wet/dry fronts: Identify dry cells: through a tolerance parameter Consistent depth reconstruction: satisfy h = b to high-order on x x wet/dry fronts Satisfy an extended C-property: Redefinition of the bed slope, numerical fluxes are computed assuming temporarily zero velocity at wet/dry faces (Brufau et al., 24) Time Integration (should at least match the order of truncation errors from dispersion terms): Third order Adams-Basforth predictor and fourth-order Adams-Moulton corrector stage.
32 Numerical Model (cont.) Special treatment wet/dry fronts: Identify dry cells: through a tolerance parameter Consistent depth reconstruction: satisfy h = b to high-order on x x wet/dry fronts Satisfy an extended C-property: Redefinition of the bed slope, numerical fluxes are computed assuming temporarily zero velocity at wet/dry faces (Brufau et al., 24) Time Integration (should at least match the order of truncation errors from dispersion terms): Third order Adams-Basforth predictor and fourth-order Adams-Moulton corrector stage. Extract depth averaged velocities, u, from the velocities functions P by solving a tridiagonal system,
33 Numerical Model (cont.) Special treatment wet/dry fronts: Identify dry cells: through a tolerance parameter Consistent depth reconstruction: satisfy h = b to high-order on x x wet/dry fronts Satisfy an extended C-property: Redefinition of the bed slope, numerical fluxes are computed assuming temporarily zero velocity at wet/dry faces (Brufau et al., 24) Time Integration (should at least match the order of truncation errors from dispersion terms): Third order Adams-Basforth predictor and fourth-order Adams-Moulton corrector stage. Extract depth averaged velocities, u, from the velocities functions P by solving a tridiagonal system, Thomas algorithm.
34 Numerical Test & Results in 1D: (Solitary wave propagation) Two cases: A/ =.1 Dx =.5 C r =.4 (MS left, Nwogu s right) η/h η/h x(m) Analytical t=5s t=15s x(m) Analytical t=5s t=15s η/h x Ct η/d x Ct η/a Nwogu MS
35 Solitary wave propagation Two cases: A/h =.2 Dx =.5 C r =.4 (MS left, Nwogu s right) η/h η/h x(m) x(m).2 Analytical t=5s t=15s.2 Analytical t=5s t=15s η/h.1 η/h x Ct x Ct η/a Nwogu MS
36 Head on collision of two solitary waves A/h =.3 Dx =.1m C r =.4, x [, 3].3 η/h.2.1 η/h Analytical.6 Nwogu MS.4.2 η/h Analytical Nwogu MS x/d
37 Regular wave propagation over a submerged bar T = 2.2s, H =.2, h/l =.11, Dx =.4m.1 Depth(m) Bar 5 1 x(m)
38 Regular wave propagation over a submerged bar T = 2.2s, H =.2, h/l =.11, Dx =.4m η(m).4.2 WG 4 MS Nwogu Experimental η(m).4.2 WG t(s) t(s)
39 Regular wave propagation over a submerged bar T = 2.2s, H =.2, h/l =.11, Dx =.4m η(m).4.2 WG 4 MS Nwogu Experimental η(m).4.2 WG t(s) t(s).4.4 WG 6 WG 7 η(m).2 η(m) t(s) t(s)
40 Regular wave propagation over a submerged bar T = 2.2s, H =.2, h/l =.11, Dx =.4m η(m).4.2 WG 4 MS Nwogu Experimental η(m).4.2 WG t(s) t(s).4.4 WG 6 WG 7 η(m).2 η(m) t(s) t(s).4.4 WG 8 WG 1 η(m).2 η(m) t(s) t(s)
41 Mathematical Model: The NSWE U t U = + H( U) = L( U) on Ω [, t] R 2 R +, H Hu Hv, H(U) = [F, G] = Hu Hu gh2 Huv Hv Huv Hv gh2, R 1 = [ L(U) = [R 1 + R 2 + S b ] gh b(x,y) x ] T and R 2 = [ gh b(x,y) y ] T. S = [ ghs f x ghs f y ] T with S f x = n2 m u u H 4 3 and Sy f = n2 m v u, H 4 3
42 NCFV approach C P U t dω + HdΩ = SdΩ C P C P UdΩ + H ñdl = SdΩ t C P C P C P
43 NCFV approach C P U t dω + HdΩ = SdΩ C P C P Introducing the flux vectors ) Φ PQ = (Fñ x + Gñ y dl and Φ P,Γ = C PQ Hence, FV scheme reads U P = 1 Φ PQ 1 t C P C Q K P Φ P,Γ + 1 C P P UdΩ + H ñdl = SdΩ t C P C P C P C P Γ C P. (Fñ x + Gñ y ) dl (S b + S d + S f )dω
44 NCFV approach C P U t dω + HdΩ = SdΩ C P C P Introducing the flux vectors ) Φ PQ = (Fñ x + Gñ y dl and Φ P,Γ = C PQ Hence, FV scheme reads U P = 1 Φ PQ 1 t C P C Q K P Φ P,Γ + 1 C P P Φ PQ Numerical flux UdΩ + H ñdl = SdΩ t C P C P C P C P Γ C P. (Fñ x + Gñ y ) dl (S b + S d + S f )dω
45 NCFV approach: Gradient & divergence formulas Green -Gauss linear reconstruction on Ω P C P = 1 3 Ω P wda = Ω P ( w) P = wñdl Ω P 1 ) 1 (w P + w Q n PQ C P 2 Q K P
46 NCFV approach: Gradient & divergence formulas Green -Gauss linear reconstruction on Ω P C P = 1 3 Ω P wda = Ω P Ω P udω = ( u) P = ( w) P = 3 u ñdl = Ω P 2 (u P + u Q ) n PQ Q K P 1 (u P + u Q ) n PQ 2 C P Q K P wñdl Ω P 1 ) 1 (w P + w Q n PQ C P 2 Q K P
47 NCFV approach: Gradient & divergence formulas Green -Gauss linear reconstruction on Ω P C P = 1 3 Ω P wda = Ω P ( w) P = wñdl Ω P 1 ) 1 (w P + w Q n PQ C P 2 Q K P
48 NCFV approach: Gradient & divergence formulas Green -Gauss linear reconstruction on Ω P C P = 1 3 Ω P wda = Ω P ( w) P = wñdl Ω P 1 ) 1 (w P + w Q n PQ C P 2 Q K P w L i,pq = w i,p LIM ( ( w i ) upw P r PQ, ( w i ) cent r PQ ) w R i,pq = w i,q 1 2 LIM ( ( w i ) upw Q r PQ, ( w i ) cent r PQ ) ( w i ), cent r PQ = w i,q w i,p, ( w i ) P, upw r PQ = 2( w i ) P ( w i ) cent
49 CCFV approach U p T p = t q K(p) Φ q + SdΩ. T p Φ q Numerical flux. Linear reconstruction Naive calculation (at point D) (w i,p ) L D = w i,p + r pd r pq LIM ( ( w i ) upw p r pq, ( w i ) cent r pq ) ; (w i,q ) R D = w i,q r Dq r pq LIM ( ( w i ) upw q r pq, ( w i ) cent r pq ),
50 CCFV approach U p T p = t q K(p) Φ q + SdΩ. T p Φ q Numerical flux. Linear reconstruction Naive calculation (at point D) (w i,p ) L D = w i,p + r pd r pq LIM ( ( w i ) upw p r pq, ( w i ) cent r pq ) ; (w i,q ) R D = w i,q r Dq r pq LIM ( ( w i ) upw q r pq, ( w i ) cent r pq ), Corrected calculation (at point M) (Delis et al., 211) w L i,p = (w i,p ) L D + r DM ( w i ) p, w R i,q = (w i,q ) R D + r DM ( w i ) q.
51 CCFV approach:gradient formulas Three element (compact stencil) gradient w i,p = 1 C c p q,r K(p) r q 1 2 (w i,q + w i,r )n q,r Extended element (wide stencil) gradient w i,p = 1 ) 1 (w Cp w i,l + w i,r n lr 2 l,r K (p) r l
52 CCFV approach:gradient formulas Three element (compact stencil) gradient w i,p = 1 C c p q,r K(p) r q 1 2 (w i,q + w i,r )n q,r Extended element (wide stencil) gradient w i,p = 1 ) 1 (w Cp w i,l + w i,r n lr 2 l,r K (p) r l Boundary conditions: use the theory of characteristics for weak formulation for the NCFV and ghost cells for CCFV scheme
53 Topography source discretization Wel-ballanced schemes introduce topography source flux vectors S b : S b (U )dxdy = S bq (CCFV) T p q K(p) S b (U )dxdy = (NCFV) C P Q K P S bpq S b = ( 1 2 P I ) Λ Λ 1 P 1 Sb where (for 1st order scheme): S b q = g HL + H R ( b R b L) n qx 2 S b PQ = g HL + H R ( b R b L) n PQx 2 g HL + H R ( b R b L) n qy g HL + H R ( b R b L) n PQy 2 q 2 PQ
54 Topography source discretization Wel-ballanced schemes introduce topography source flux vectors S b : S b (U )dxdy = S bq (CCFV) T p q K(p) S b (U )dxdy = (NCFV) C P Q K P S bpq S b = ( 1 2 P I ) Λ Λ 1 P 1 Sb + S /2nd order scheme, correction b term S b q = g HL + H p ( ) b L b p nqx 2 S b PQ = g HL + H P ( ) b L b P npqx 2 g HL + H p ( ) b L b p nqy g HL + H P ( ) b L b P npqy 2 2.
55 Topography source discretization (wet/dry) Extended C-property, (Castro et al, 25) In the MUSCL scheme for hydrostatic conditions we must have, at i-cell b L b i = (h L h i ) ( B) i = ( h) i If in the gradient calculation, of a wet cell, a dry node is involved we correct the h L and/or h R by imposing h L = h i (b L b i ) and/or h R = h j (b R b j ) Redefine the bed value in the dry node for emerging bed situations in Sb to maintain hydrostatic conditions, (Brufau et al, 22). For water in motion over emerging slopes: If h L > ɛ wd and h R ɛ wd and h L < (b B b L ), set temporarily for the wet i-cell u L = v L =
56 Numerical Results Grids used: Consistently refined grids, for N =degrees of freedom, characteristic length h N = L x L y /N Equivalent meshes Scheme NCFV CCFVc1 CCFVc2 CCFVw1 CCFVw2 Description Node-Centered FV Scheme Cell-Centered FV compact (naive) reconstruction stencil Cell-Centered FV compact reconstruction stencil (corrected) Cell-Centered FV wide (naive) reconstruction stencil Cell-Centered FV wide reconstruction stencil (corrected)
57 Traveling vortex: Convergence results The traveling vortex solution
58 Comparison 1. Convergence behavior to second order: NCFV is not grid dependent/ CCFV depents on the grid used.
59 Comparison 1. Convergence behavior to second order: NCFV is not grid dependent/ CCFV depents on the grid used. 2. Edge based limiting procedure in the MUSCL reconstruction: inadequate for the CCFV schemes if the proposed correction not used.
60 Comparison 1. Convergence behavior to second order: NCFV is not grid dependent/ CCFV depents on the grid used. 2. Edge based limiting procedure in the MUSCL reconstruction: inadequate for the CCFV schemes if the proposed correction not used. 3. CCFVw2 has an almost identical behavior with the NCFV scheme is achieved, but with the extra computational cost introduced.
61 Comparison 1. Convergence behavior to second order: NCFV is not grid dependent/ CCFV depents on the grid used. 2. Edge based limiting procedure in the MUSCL reconstruction: inadequate for the CCFV schemes if the proposed correction not used. 3. CCFVw2 has an almost identical behavior with the NCFV scheme is achieved, but with the extra computational cost introduced. 4. Ghost cells for the boundary treatment in CCFV formulations can lead to an order reduction.
62 Comparison 1. Convergence behavior to second order: NCFV is not grid dependent/ CCFV depents on the grid used. 2. Edge based limiting procedure in the MUSCL reconstruction: inadequate for the CCFV schemes if the proposed correction not used. 3. CCFVw2 has an almost identical behavior with the NCFV scheme is achieved, but with the extra computational cost introduced. 4. Ghost cells for the boundary treatment in CCFV formulations can lead to an order reduction. 5. Wet/dry treatment: accurate for both FV approaches on all grid types.
63 Comparison 1. Convergence behavior to second order: NCFV is not grid dependent/ CCFV depents on the grid used. 2. Edge based limiting procedure in the MUSCL reconstruction: inadequate for the CCFV schemes if the proposed correction not used. 3. CCFVw2 has an almost identical behavior with the NCFV scheme is achieved, but with the extra computational cost introduced. 4. Ghost cells for the boundary treatment in CCFV formulations can lead to an order reduction. 5. Wet/dry treatment: accurate for both FV approaches on all grid types.
64 An unstructured FV scheme for BT Equations Vector conservative form for Nwogu s equations: U t + H( U ) = L( U ) on Ω [, t] R 2 R +, U vector of the new variables, U = [H, Hu, Hv] T and H = [F, G] H U = P 1, L(U) = [S b + S d + S f ] P 2 [ ] [ ] P1 z 2 with P = = H a P ( u) + 2 z a ( hu) + u and 2 ψ c S b = ghb x, S d = uψ c + ψ Mx, S f = Sx f + R bx ghb y vψ c + ψ My Sy f + R by
65 Vector conservative form (cont.) [( ) ( z 2 ψ c = a 2 h2 h ( u) + z a + h ) ] h ( hu) 6 2 ψ M = [ ψmx ψ My ] = H t z 2 a 2 ( u) + H tz a ( hu)
66 Vector conservative form (cont.) [( ) ( z 2 ψ c = a 2 h2 h ( u) + z a + h ) ] h ( hu) 6 2 ψ M = [ ψmx ψ My ] = H t z 2 a 2 ( u) + H tz a ( hu) R b = [R bx, R by ] T = parametrization of wave breaking characteristics where: [ ν ] T R bx = R bx, where R bx = ν(hu) x 2 ((Hu) y + (Hv) x ) and [ ν ] T R by = R by, where R by = 2 ((Hu) y + (Hv) x ) ν(hv) y. where ν = Bδ 2 b Hη t is the eddy viscosity coefficient with < B < 1 and δ b is a mixing length coefficient.
67 Numerical Model: Spatial discretization Advective (nonlinear) part and topography source term: Well-balanced FV formulation. Roe s approximate Riemann solver is used (Roe, 1981). Upwinding of the topography source term (Bermudez et al., 1994). High-order spatial accuracy: third-order MUSCL-type scheme (Barth, 1993). Satisfy the C-property (flow at rest) to higher spatial order: Addition of an extra term to bed upwinding (Hubbard and Garcia-Navarro, 2 and Delis and Nikolos, 29). Special treatment of wet/dry fronts. Dispersion terms: consistent FV approximations based on gradient and divergence computations (Kazolea et al., 212).
68 Numerical Model: Spatial discretization Advective (nonlinear) part and topography source term: Well-balanced FV formulation. Roe s approximate Riemann solver is used (Roe, 1981). Upwinding of the topography source term (Bermudez et al., 1994). High-order spatial accuracy: third-order MUSCL-type scheme (Barth, 1993). Satisfy the C-property (flow at rest) to higher spatial order: Addition of an extra term to bed upwinding (Hubbard and Garcia-Navarro, 2 and Delis and Nikolos, 29). Special treatment of wet/dry fronts. Dispersion terms: consistent FV approximations based on gradient and divergence computations (Kazolea et al., 212).
69 Numerical Model: Spatial discretization Advective (nonlinear) part and topography source term: Well-balanced FV formulation. Roe s approximate Riemann solver is used (Roe, 1981). Upwinding of the topography source term (Bermudez et al., 1994). High-order spatial accuracy: third-order MUSCL-type scheme (Barth, 1993). Satisfy the C-property (flow at rest) to higher spatial order: Addition of an extra term to bed upwinding (Hubbard and Garcia-Navarro, 2 and Delis and Nikolos, 29). Special treatment of wet/dry fronts. Dispersion terms: consistent FV approximations based on gradient and divergence computations (Kazolea et al., 212).
70 Numerical Model: Spatial discretization Advective (nonlinear) part and topography source term: Well-balanced FV formulation. Roe s approximate Riemann solver is used (Roe, 1981). Upwinding of the topography source term (Bermudez et al., 1994). High-order spatial accuracy: third-order MUSCL-type scheme (Barth, 1993). Satisfy the C-property (flow at rest) to higher spatial order: Addition of an extra term to bed upwinding (Hubbard and Garcia-Navarro, 2 and Delis and Nikolos, 29). Special treatment of wet/dry fronts. Dispersion terms: consistent FV approximations based on gradient and divergence computations (Kazolea et al., 212).
71 Numerical Model: Spatial discretization Advective (nonlinear) part and topography source term: Well-balanced FV formulation. Roe s approximate Riemann solver is used (Roe, 1981). Upwinding of the topography source term (Bermudez et al., 1994). High-order spatial accuracy: third-order MUSCL-type scheme (Barth, 1993). Satisfy the C-property (flow at rest) to higher spatial order: Addition of an extra term to bed upwinding (Hubbard and Garcia-Navarro, 2 and Delis and Nikolos, 29). Special treatment of wet/dry fronts. Dispersion terms: consistent FV approximations based on gradient and divergence computations (Kazolea et al., 212).
72 Numerical Model: Spatial discretization Advective (nonlinear) part and topography source term: Well-balanced FV formulation. Roe s approximate Riemann solver is used (Roe, 1981). Upwinding of the topography source term (Bermudez et al., 1994). High-order spatial accuracy: third-order MUSCL-type scheme (Barth, 1993). Satisfy the C-property (flow at rest) to higher spatial order: Addition of an extra term to bed upwinding (Hubbard and Garcia-Navarro, 2 and Delis and Nikolos, 29). Special treatment of wet/dry fronts. Dispersion terms: consistent FV approximations based on gradient and divergence computations (Kazolea et al., 212).
73 Numerical Model: Spatial discretization Advective (nonlinear) part and topography source term: Well-balanced FV formulation. Roe s approximate Riemann solver is used (Roe, 1981). Upwinding of the topography source term (Bermudez et al., 1994). High-order spatial accuracy: third-order MUSCL-type scheme (Barth, 1993). Satisfy the C-property (flow at rest) to higher spatial order: Addition of an extra term to bed upwinding (Hubbard and Garcia-Navarro, 2 and Delis and Nikolos, 29). Special treatment of wet/dry fronts. Dispersion terms: consistent FV approximations based on gradient and divergence computations (Kazolea et al., 212).
74 Discretization of the dispersive terms (mass equation) Integral averaging: (ψ c ) P = 1 [( z 2 a C P C P 2 h2 6 = 1 { [( z 2 a C P 2 h2 6 Q K P C PQ ) ( h ( u) + z a + h ) ] h ( hu) dω 2 ) ] [( h ( u) ñdl + z a + h ) ] } h ( hu) ñdl C PQ 2 C PQ C PQ ( ) [( ) ] z 2 a 2 h2 z 2 h ( u) ñdl a 6 2 h2 h [ ( u) n PQ ] 6 M, M ( z a + h 2 ) [( h ( hu) ñdl z a + h ) ] h [ ( hu) n PQ ] 2 M M K PQ := {R N R is a vertex of M PQ and RQ M PQ } ( w) M = 1 M PQ R,Q K PQ R Q 1 2 (w R + w Q ) n RQ
75 Discretization of the dispersive terms (momentum equations) 1 ( uψ c + ψ M ) dω = u P ψ c dω + 1 ψ M dω. C P C P C P C P C P C P The ψ c is discretized as before and the second term takes the discrete form: (ψ M ) P = 1 C P = 1 C P [ za 2 H t 2 ψ M dω = 1 C P C P z H t CP 2 ] P CP H t z 2 a a 1 ( u)dω + 2 C P [ ( u)] P + [H t z a ] P [ ( hu)] P, 2 ( u) + H tz a ( hu)dω H t z a ( hu)dω C P
76 Time integration Consider the semi-discrete scheme: U P = L (U) t Time Integration (match the order of truncation errors from dispersion terms): Use 3rd order explicit Strong Stability-Preserving Runge-Kutta (SSP-RK): U (1) P = U (n) P + t n L ( U (n)) ; U (2) = 3 P 4 U(n) + 1 P 4 U(1) + P t n 1 4 L ( U (1)) ; U (n+1) = 1 P 3 U(n) + 2 P 3 U(2) + P t n 2 3 L ( U (2)) Time step t n estimated by a CFL stability condition as t n R P = CFL min P ( u2 + v 2 + ) n c P
77 Velocity field recovery From new solution variables P = [P 1, P 2 ] T At each step in the RK scheme a linear system MV = C with M R 2N 2N and C = [P 1 P 2 P N ] T, has to be solved to obtain the velocities V = [u 1 u 2 u N ] T. Each two rows of the system read as [ ] H (i) z 2 (i) a P 2 ( u) + z a ( hu) + u P = P (i), i = 1, 2, n + 1. P
78 Velocity field recovery From new solution variables P = [P 1, P 2 ] T At each step in the RK scheme a linear system MV = C with M R 2N 2N and C = [P 1 P 2 P N ] T, has to be solved to obtain the velocities V = [u 1 u 2 u N ] T. Each two rows of the system read as [ ] H (i) z 2 (i) a P 2 ( u) + z a ( hu) + u P = P (i), i = 1, 2, n + 1. P Important to (a) keep the unknown information needed at the minimum possible level (i.e neighboring nodes) and (b) exploit already computed geometrical information. (z H a) P 2 P 2 1 C P ( u) M n PQ + (z a) P C Q K P P ( hu) M n PQ + u P = P P Q K P
79 Velocity field recovery From new solution variables P = [P 1, P 2 ] T At each step in the RK scheme a linear system MV = C with M R 2N 2N and C = [P 1 P 2 P N ] T, has to be solved to obtain the velocities V = [u 1 u 2 u N ] T. Each two rows of the system read as [ ] H (i) z 2 (i) a P 2 ( u) + z a ( hu) + u P = P (i), i = 1, 2, n + 1. P Important to (a) keep the unknown information needed at the minimum possible level (i.e neighboring nodes) and (b) exploit already computed geometrical information. (z 2 a) P 2 C P (z H a) P 2 P 2 1 C P A Q u Q + A P u P + (z a) P C Q K P P ( u) M n PQ + (z a) P C Q K P P ( hu) M n PQ + u P = P P Q K P B Q u Q + B P u P + Iu P = 1 P P, P = 1... N H Q K P P
80 Solution of the linear system The matrix M is sparse and structurally symmetric but is also mesh dependent.
81 Solution of the linear system The matrix M is sparse and structurally symmetric but is also mesh dependent. Matrix M is stored in the compressed sparse row (CSR) format.
82 Solution of the linear system The matrix M is sparse and structurally symmetric but is also mesh dependent. Matrix M is stored in the compressed sparse row (CSR) format. The ILUT preconditioner from SPARSKIT package is used.
83 Solution of the linear system The matrix M is sparse and structurally symmetric but is also mesh dependent. Matrix M is stored in the compressed sparse row (CSR) format. The ILUT preconditioner from SPARSKIT package is used. The reverse Cuthill McKee (RCM) algorithm is also employed to reorder the matrix elements as to minimize the matrix bandwidth.
84 Solution of the linear system The matrix M is sparse and structurally symmetric but is also mesh dependent. Matrix M is stored in the compressed sparse row (CSR) format. The ILUT preconditioner from SPARSKIT package is used. The reverse Cuthill McKee (RCM) algorithm is also employed to reorder the matrix elements as to minimize the matrix bandwidth. System is solved using Bi-Conjugate Gradient Stabilized method (BiCGStab), with tolerance 5 1 6
85 Solution of the linear system The matrix M is sparse and structurally symmetric but is also mesh dependent. Matrix M is stored in the compressed sparse row (CSR) format. The ILUT preconditioner from SPARSKIT package is used. The reverse Cuthill McKee (RCM) algorithm is also employed to reorder the matrix elements as to minimize the matrix bandwidth. System is solved using Bi-Conjugate Gradient Stabilized method (BiCGStab), with tolerance Convergence to the solution was obtained in one or two steps with the numerical solution for the velocities at the previous time step given as initial guess.
86 Boundary conditions and the internal source function Wall (reflective) boundary condition: u ñ = for x Ω By conservation of mass (no loss or gain through the wall) [ ( ) ( z 2 HdΩ + Hu + a t Ω Ω 2 h2 h ( u) + z a + h ) ] h ( hu) ñdl = 6 2 Define the normal boundary advective flux in weak form, 1 Φ P,Γ = 2 g(h ) 2 n P,1x by the method of characteristics 1 2 g(h ) 2 n P,1y
87 Boundary conditions and the internal source function Wall (reflective) boundary condition: u ñ = for x Ω By conservation of mass (no loss or gain through the wall) [ ( ) ( z 2 HdΩ + Hu + a t Ω Ω 2 h2 h ( u) + z a + h ) ] h ( hu) ñdl = 6 2 Define the normal boundary advective flux in weak form, 1 Φ P,Γ = 2 g(h ) 2 n P,1x by the method of characteristics 1 2 g(h ) 2 n P,1y Absorbing boundaries: should dissipate the energy of incoming waves ( ) x d(x) 2 Sponge layer is defined: m(x) = 1, L L s 1.5L, L s
88 Boundary conditions and the internal source function Wall (reflective) boundary condition: u ñ = for x Ω By conservation of mass (no loss or gain through the wall) [ ( ) ( z 2 HdΩ + Hu + a t Ω Ω 2 h2 h ( u) + z a + h ) ] h ( hu) ñdl = 6 2 Define the normal boundary advective flux in weak form, 1 Φ P,Γ = 2 g(h ) 2 n P,1x by the method of characteristics 1 2 g(h ) 2 n P,1y Absorbing boundaries: should dissipate the energy of incoming waves ( ) x d(x) 2 Sponge layer is defined: m(x) = 1, L L s 1.5L, Internal source function for regular waves (Wei et al., 1993) added to the mass equation S(x, t) = D exp ( γ(x x s ) 2) sin(λy ωt) L s
89 Wave breaking models Eddy viscosity wave breaking model (R b ) P = = 1 C P 1 C P R b dω = 1 C P C P C P [ ] Rbx ñ dl 1 R by ñ C P C PQ [ R by ] dω R bx [ ] Rbx n PQ R by n PQ M
90 Wave breaking models Eddy viscosity wave breaking model (R b ) P = Hybrid models = 1 C P 1 C P R b dω = 1 C P C P C P [ ] Rbx ñ dl 1 R by ñ C P C PQ [ R by ] dω R bx [ ] Rbx n PQ R by n PQ Hybrid(ɛ) BT degenerate into NSWE as dispersive terms become negligible. Criterion: ɛ = A h.8 In post breaking region ɛ <.4 in order to switch NSWE/BT. M
91 Wave breaking models Eddy viscosity wave breaking model (R b ) P = Hybrid models = 1 C P 1 C P R b dω = 1 C P C P C P [ ] Rbx ñ dl 1 R by ñ C P C PQ [ R by ] dω R bx [ ] Rbx n PQ R by n PQ Hybrid(ɛ) BT degenerate into NSWE as dispersive terms become negligible. Criterion: ɛ = A h.8 In post breaking region ɛ <.4 in order to switch NSWE/BT. New Hybrid criteria: η t γ gh, γ [.35,.65], η x tan(φ c ). Distinguish the different breaking waves. Find non-breaking undular bores checking the Froude number. Extend the computational region of the NSWE. M
92 Suppression of the dispersive terms Not clear how the switching between the two models is implemented
93 Suppression of the dispersive terms Not clear how the switching between the two models is implemented Discontinuity at the switching point BT/NSWE, causing spurious oscillations.
94 Suppression of the dispersive terms Not clear how the switching between the two models is implemented Discontinuity at the switching point BT/NSWE, causing spurious oscillations. Matrix M can not be change.
95 Suppression of the dispersive terms Not clear how the switching between the two models is implemented Discontinuity at the switching point BT/NSWE, causing spurious oscillations. Matrix M can not be change. Need of a clever implementation.
96 Suppression of the dispersive terms Not clear how the switching between the two models is implemented Discontinuity at the switching point BT/NSWE, causing spurious oscillations. Matrix M can not be change. Need of a clever implementation. Methodology. Starting with the solution vector U n P, P = 1,..., N, at time t n, 1. H is computed C P using the BT model (named from now on H n+1 ). BT 1.1 If breaking is on for N br < N cells additional solution vector H n+1 ψ BT c named H n+1. BT/SW 2. P n+1 BT is computed C P / t H n+1 Hn+1 BT Hn for the ψ t n+1 M. 2.1 N br cells P n+1 BT/SW Pn+1 BT 3. MV = C,with C = [P n+1 1, P n+1 2,, P n+1 ψ c ψ M [(Hu), (Hv)] n+1,t N ]T BT un+1 BT 4. Final solution: H n+1 BT/SW, Pn+1 BT/SW, un+1 BT/SW un+1 BT vector with its values at the breaking nodes replaced by those of u n+1 SW.
97 Numerical Results:Spatial accuracy Solitary wave: A/h =.1, (x, y) [, 3] [, 5m] Reference solution of N = 232, 849 nodes
98 Solitary wave interaction with a vertical circular cylinder Area: (x, y) = [ 4, 1m] [,.55m], A/h =.25, N = 1, 69
99 Solitary wave interaction with a vertical circular cylinder Wave Gauges:
100 Run-up of a solitary wave on a conical island (Briggs et al. 1995) Area: (x, y) = [ 5, 28m] [, 3m], A/h =.18, N = 52, 191, CFL =.8 Using mesh h enrichment
101 Run-up of a solitary wave on a conical island (Briggs et al. 1995) Area: (x, y) = [ 5, 28m] [, 3m], A/h =.18, N = 52, 191, CFL =.8 Using mesh h enrichment
102 Run-up of a solitary wave on a conical island (Briggs et al. 1995) Area: (x, y) = [ 5, 28m] [, 3m], A/h =.18, N = 52, 191, CFL =.8 Using mesh h enrichment
103 Run-up of a solitary wave on a conical island (Briggs et al. 1995) Area: (x, y) = [ 5, 28m] [, 3m], A/h =.18, N = 52, 191, CFL =.8 Using mesh h enrichment
104 Run-up of a solitary wave on a conical island (Briggs et al. 1995) Area: (x, y) = [ 5, 28m] [, 3m], A/h =.18, N = 52, 191, CFL =.8 Using mesh h enrichment
105 Run-up of a solitary wave on a conical island (Briggs et al. 1995) Area: (x, y) = [ 5, 28m] [, 3m], A/h =.18, N = 52, 191, CFL =.8 Using mesh h enrichment
106 Run-up of a solitary wave on a conical island (Briggs et al. 1995) Area: (x, y) = [ 5, 28m] [, 3m], A/h =.18, N = 52, 191, CFL =.8 Using mesh h enrichment
107 Run-up of a solitary wave on a conical island (cont) Time series of surface elevation at wave gauges around the island:
108 Run-up of a solitary wave on a conical island (cont) Time series of surface elevation at wave gauges around the island: Experimental measurements and numerical runup around the island: Simulation 28min on a single 2.4GHz Intel Core 2 Quad Q66 processor
109 Wave propagation over a semicircular shoal (Whalin 1971) T = 2.s, h/l =.117, A/h =.165, kh =.735 and S = 1.198
110 Wave propagation over a semicircular shoal (Whalin 1971) T = 2.s, h/l =.117, A/h =.165, kh =.735 and S = Free surface and spatial evolution of the 1st, 2nd and 3rd harmonic
111 Solitary wave run-up on a plane beach (Synolakis, 1987) Area: (x, y) = [ 2, 6m] [, 1m], A/h =.28, N = 8, 816, CFL =.4, n m =.1 H(m) Bed New Hybrid Eddy Visc. Experimental t g/h =
112 Solitary wave run-up on a plane beach (Synolakis, 1987) Area: (x, y) = [ 2, 6m] [, 1m], A/h =.28, N = 8, 816, CFL =.4, n m =.1 H(m) Bed New Hybrid Eddy Visc. Experimental H(m) t g/h = t g/h =
113 Solitary wave run-up on a plane beach (Synolakis, 1987) Area: (x, y) = [ 2, 6m] [, 1m], A/h =.28, N = 8, 816, CFL =.4, n m =.1 H(m) Bed New Hybrid Eddy Visc. Experimental H(m) H(m) t g/h = t g/h = t g/h =
114 Solitary wave run-up on a plane beach (Synolakis, 1987) Area: (x, y) = [ 2, 6m] [, 1m], A/h =.28, N = 8, 816, CFL =.4, n m =.1 H(m) Bed New Hybrid Eddy Visc. Experimental H(m) H(m) H(m) t g/h = t g/h = t g/h = t g/h =
115 Solitary wave run-up on a plane beach (Synolakis, 1987) Area: (x, y) = [ 2, 6m] [, 1m], A/h =.28, N = 8, 816, CFL =.4, n m =.1 H(m) Bed New Hybrid Eddy Visc. Experimental H(m) H(m) H(m) t g/h = t g/h = H(m) t g/h = t g/h = t g/h =
116 Solitary wave run-up on a plane beach (Synolakis, 1987) Area: (x, y) = [ 2, 6m] [, 1m], A/h =.28, N = 8, 816, CFL =.4, n m =.1 H(m) Bed New Hybrid Eddy Visc. Experimental H(m) H(m) H(m) t g/h = t g/h = H(m) t g/h = t g/h = H(m) t g/h = t g/h =7 1 5 x(m)
117 Wave over a bar (Beji and Battjes, 1993) (x, y) = [ 1, 3m] [,.8m], H =.2m, T = 2.2s, N = 4, Depth(m) Bar 5 1 x(m)
118 Wave over a bar (Beji and Battjes, 1993) (x, y) = [ 1, 3m] [,.8m], H =.2m, T = 2.2s, N = 4, 364 Depth(m) Bar 5 1 x(m) H(m).5.4 Free surface.3 l NSW Bar
119 Wave over a bar (Beji and Battjes, 1993) (x, y) = [ 1, 3m] [,.8m], H =.2m, T = 2.2s, N = 4, 364 Depth(m) H(m) Bar 5 1 x(m) H(m) H(m).5.4 Free surface.3 l NSW Bar
120 Wave over a bar (Beji and Battjes, 1993) (x, y) = [ 1, 3m] [,.8m], H =.2m, T = 2.2s, N = 4, 364 Depth(m) H(m) Bar 5 1 x(m) H(m) H(m).5.4 Free surface.3 l NSW Bar H(m).4.3 H(m) x(m) 15 2
121 Wave over a bar (Beji and Battjes, 1993) Wave gauges: η(m) η(m) η(m) η(m) Hybrid (ε) New Hybrid Experimental data t(s)
122 A two-dimensional reef(roeber, 29) η/h Area: (x, y) = [, 83m] [, 1m], A/h =.3, N = 1, 9, n m = t =31.5s Bed NSW BT Experimental Shoaling
123 A two-dimensional reef(roeber, 29) η/h Area: (x, y) = [, 83m] [, 1m], A/h =.3, N = 1, 9, n m = t =31.5s Bed NSW BT Experimental η/h t =34.5s Shoaling During breaking
124 A two-dimensional reef(roeber, 29) η/h Area: (x, y) = [, 83m] [, 1m], A/h =.3, N = 1, 9, n m = t =31.5s Bed NSW BT Experimental η/h η/h t =34.5s t =35.5s Shoaling During breaking Wave jet hits still water
125 .4.2 A two-dimensional reef(roeber, 29) η/h η/h η/h Area: (x, y) = [, 83m] [, 1m], A/h =.3, N = 1, 9, n m = t =31.5s Bed NSW BT Experimental t =34.5s t =35.5s t =38.5s Shoaling During breaking Wave jet hits still water Bore propagation
126 A two-dimensional reef(roeber, 29) η/h.4.2 t =59.5s Hydraulic jump
127 A two-dimensional reef(roeber, 29) η/h.4.2 t =59.5s Hydraulic jump Time series of surface elevation at wave gauges: η/h η/h η/h Gauge at (x =44.25m, y =.5m) Gauge at (x =5.37m, y =.5m) η/h Gauge at (x =48.233m, y =.5m) Gauge at (x =8.m, y =.5m) t(s) 6 8 1
128 A two-dimensional reef(roeber, 29) Spatial snapshots along the centerline: η/h.4 t =34.5s.2 Bed l NSW Free surface η/h.4 t =36.5s.2 η/h η/h t =38.5s t =49.5s.2 η/h η/h t =4.5s t =56.5s
129 A three-dimensional reef(lynett et al., 29) Area: (x, y) = [, 45m] [ 13m, 13m], A/h =.5, N = 87, 961
130 A three-dimensional reef(lynett et al., 29) Area: (x, y) = [, 45m] [ 13m, 13m], A/h =.5, N = 87, 961
131 A three-dimensional reef(lynett et al., 29) Time series of surface elevation at wave gauges: η WG 2
132 A three-dimensional reef(lynett et al., 29) Time series of surface elevation at wave gauges: η η WG 2 WG 3
133 A three-dimensional reef(lynett et al., 29) Time series of surface elevation at wave gauges: η η η WG 2 WG 3 WG 6
134 A three-dimensional reef(lynett et al., 29) Time series of surface elevation at wave gauges: η η η η.2 WG 2 WG 3 WG 6 WG t(s)
135 A three-dimensional reef(lynett et al., 29) Time series of velocities at wave gauges: u(m/s) 2 ADV 1 Experimental New Hybrid
136 A three-dimensional reef(lynett et al., 29) Time series of velocities at wave gauges: u(m/s) 2 ADV 1 Experimental New Hybrid ADV 2 2 u(m/s)
137 A three-dimensional reef(lynett et al., 29) Time series of velocities at wave gauges: u(m/s) 2 ADV 1 Experimental New Hybrid ADV 2 2 u(m/s) u(m/s) ADV
138 A three-dimensional reef(lynett et al., 29) Time series of velocities at wave gauges: u(m/s) 2 ADV 1 Experimental New Hybrid ADV 2 2 u(m/s) u(m/s) v(m/s) ADV ADV t(s)
139 Conclusions A 2D unstructured FV scheme numerical model has been developed for solving Nwogu s extended BT equations formulated as to have identical flux terms as to the NSWE.
140 Conclusions A 2D unstructured FV scheme numerical model has been developed for solving Nwogu s extended BT equations formulated as to have identical flux terms as to the NSWE. Conservative formulation and higher-order FV scheme enhance the applicability of the model without altering its dispersion characteristics.
141 Conclusions A 2D unstructured FV scheme numerical model has been developed for solving Nwogu s extended BT equations formulated as to have identical flux terms as to the NSWE. Conservative formulation and higher-order FV scheme enhance the applicability of the model without altering its dispersion characteristics. The well-balanced topography and wet/dry front discretizations provided accurate conservative and stable wave propagation shoaling and run-up.
142 Conclusions A 2D unstructured FV scheme numerical model has been developed for solving Nwogu s extended BT equations formulated as to have identical flux terms as to the NSWE. Conservative formulation and higher-order FV scheme enhance the applicability of the model without altering its dispersion characteristics. The well-balanced topography and wet/dry front discretizations provided accurate conservative and stable wave propagation shoaling and run-up. The edge-based structure adopted can provide computational efficiency, since most of the geometric quantities needed can be calculated in a pre-processing stage.
143 Conclusions A 2D unstructured FV scheme numerical model has been developed for solving Nwogu s extended BT equations formulated as to have identical flux terms as to the NSWE. Conservative formulation and higher-order FV scheme enhance the applicability of the model without altering its dispersion characteristics. The well-balanced topography and wet/dry front discretizations provided accurate conservative and stable wave propagation shoaling and run-up. The edge-based structure adopted can provide computational efficiency, since most of the geometric quantities needed can be calculated in a pre-processing stage. Different type of wave breaking can be implemented. The New Hybrid model proved more stable and accurate than the others applied in this work.
144 Conclusions A 2D unstructured FV scheme numerical model has been developed for solving Nwogu s extended BT equations formulated as to have identical flux terms as to the NSWE. Conservative formulation and higher-order FV scheme enhance the applicability of the model without altering its dispersion characteristics. The well-balanced topography and wet/dry front discretizations provided accurate conservative and stable wave propagation shoaling and run-up. The edge-based structure adopted can provide computational efficiency, since most of the geometric quantities needed can be calculated in a pre-processing stage. Different type of wave breaking can be implemented. The New Hybrid model proved more stable and accurate than the others applied in this work. Relatively straight forward to extend existing NSWE codes that use (unstructured) FV schemes as to include dispersion characteristics for deeper water simulations.
145 Conclusions A 2D unstructured FV scheme numerical model has been developed for solving Nwogu s extended BT equations formulated as to have identical flux terms as to the NSWE. Conservative formulation and higher-order FV scheme enhance the applicability of the model without altering its dispersion characteristics. The well-balanced topography and wet/dry front discretizations provided accurate conservative and stable wave propagation shoaling and run-up. The edge-based structure adopted can provide computational efficiency, since most of the geometric quantities needed can be calculated in a pre-processing stage. Different type of wave breaking can be implemented. The New Hybrid model proved more stable and accurate than the others applied in this work. Relatively straight forward to extend existing NSWE codes that use (unstructured) FV schemes as to include dispersion characteristics for deeper water simulations.
146 References 1. A. I. Delis, M. Kazolea, and N. A Kampanis. A robust high resolution finite volume scheme for the simulation of long waves over complex domain. Int. J. Num. Meth. Fluids, 56:419, A. I. Delis, I. A Nikolos, and M. Kazolea. Performance and comparison of cellcentered and node-centered unstructured finite volume discretizations for shallow water free surface flows. Arch. Comput. Methods Eng., 18:57, M. Kazolea and A. I. Delis. A well-balanced shock-capturing hybrid finite volume-finite difference numerical scheme for extended 1D Boussinesq models. Applied Numerical Mathematics, 67:167186, M. Kazolea, A. I. Delis, I. A Nikolos, and C. E. Synolakis. An unstructured finite volume numerical scheme for extended 2D Boussinesq-type equations.coast. Eng., 69:4266, 212.
147 Thank you for your attention!!
148 Solitary wave run-up on a plane beach Area: x [ 2, 6m] A/h =.4, Dx =.5, CFL =.4, n m =.1.1 H(m).5 t g/h= x(m)
149 Solitary wave run-up on a plane beach Area: x [ 2, 6m] A/h =.4, Dx =.5, CFL =.4, n m = H(m).5 H(m).5 t g/h= x(m) t g/h= x(m)
150 Solitary wave run-up on a plane beach Area: x [ 2, 6m] A/h =.4, Dx =.5, CFL =.4, n m = H(m).5 H(m).5 H(m).5 t g/h= x(m) t g/h= x(m) t g/h= x(m)
151 Solitary wave run-up on a plane beach Area: x [ 2, 6m] A/h =.4, Dx =.5, CFL =.4, n m = H(m).5 H(m).5 H(m).5 t g/h= x(m) t g/h= x(m) t g/h= x(m) H(m).5 t g/h= x(m)
152 Solitary wave run-up on a plane beach Area: x [ 2, 6m] A/h =.4, Dx =.5, CFL =.4, n m = H(m).5 H(m).5 H(m).5 t g/h= x(m) t g/h= x(m) t g/h= x(m) H(m).5 H(m).5 t g/h= x(m) t g/h= x(m)
Well-balanced shock-capturing hybrid finite volume-finite difference schemes for Boussinesq-type models
NUMAN 2010 Well-balanced shock-capturing hybrid finite volume-finite difference schemes for Boussinesq-type models Maria Kazolea 1 Argiris I. Delis 2 1 Environmental Engineering Department, TUC, Greece
More informationComparison of cell-centered and node-centered formulations of a high-resolution well-balanced finite volume scheme: application to shallow water flows
Comparison of cell-centered and node-centered formulations of a high-resolution well-balanced finite volume scheme: application to shallow water flows Dr Argiris I. Delis Dr Ioannis K. Nikolos (TUC) Maria
More informationAdvances and perspectives in numerical modelling using Serre-Green Naghdi equations. Philippe Bonneton
Long wave & run-up workshop Santander 2012 Advances and perspectives in numerical modelling using Serre-Green Naghdi equations Philippe Bonneton EPOC, METHYS team, Bordeaux Univ., CNRS d0 µ = λ 0 2 small
More informationB O S Z. - Boussinesq Ocean & Surf Zone model - International Research Institute of Disaster Science (IRIDeS), Tohoku University, JAPAN
B O S Z - Boussinesq Ocean & Surf Zone model - Volker Roeber 1 Troy W. Heitmann 2, Kwok Fai Cheung 2, Gabriel C. David 3, Jeremy D. Bricker 1 1 International Research Institute of Disaster Science (IRIDeS),
More informationNon-Hydrostatic Pressure Shallow Flows: GPU Implementation Using Finite-Volume and Finite-Difference Scheme.
arxiv:176.4551v1 [math.na] 14 Jun 17 Non-Hydrostatic Pressure Shallow Flows: GPU Implementation Using Finite-Volume and Finite-Difference Scheme. C. Escalante 1 T. Morales de Luna and M.J. Castro 1 1 Departamento
More informationNon-Hydrostatic Pressure Shallow Flows: GPU Implementation Using Finite Volume and Finite Difference Scheme.
arxiv:1706.04551v [math.na] 1 Jul 018 Non-Hydrostatic Pressure Shallow Flows: GPU Implementation Using Finite Volume and Finite Difference Scheme. C. Escalante 1 T. Morales de Luna and M.J. Castro 1 1
More informationWave breaking for Boussinesq-type models using a turbulence kinetic energy model
Wave breaking for Boussinesq-type models using a turbulence kinetic energy model Maria Kazolea, Mario Ricchiuto To cite this version: Maria Kazolea, Mario Ricchiuto. Wave breaking for Boussinesq-type models
More information2017 年環境流體力學短期講座 Short Course on Environmental Flows
2017 年環境流體力學短期講座 Short Course on Environmental Flows 數學 海浪 與沿海動態過程 Mathematics, ocean waves and coastal dynamic processes Philip L-F. Liu National University of Singapore Cornell University September 2017
More information3 where g is gravity. S o and S f are bed slope and friction slope at x and y direction, respectively. Here, friction slope is calculated based on Man
東北地域災害科学研究第 48 巻 () 3 D FORCE MUSCL SCHEME FOR SIMULATING BREAKING SOLITARY WAVE RUNUP Mohammad Bagus Adityawan * Hitoshi Tanaka ABSTRACT Breaking wave simulation using depth averaged based model, i.e.
More informationBOUSSINESQ-TYPE EQUATIONS WITH VARIABLE COEFFICIENTS FOR NARROW-BANDED WAVE PROPAGATION FROM ARBITRARY DEPTHS TO SHALLOW WATERS
BOUSSINESQ-TYPE EQUATIONS WITH VARIABLE COEFFICIENTS FOR NARROW-BANDED WAVE PROPAGATION FROM ARBITRARY DEPTHS TO SHALLOW WATERS Gonzalo Simarro 1, Alvaro Galan, Alejandro Orfila 3 A fully nonlinear Boussinessq-type
More informationStable Semi-Discrete Schemes for the 2D Incompressible Euler Equations
Stable Semi-Discrete Schemes for the D Incompressible Euler Equations Doron Levy Department of Mathematics Stanford University dlevy@math.stanford.edu http://math.stanford.edu/ dlevy Maryland, May 004
More informationNumerical simulation of dispersive waves
Numerical simulation of dispersive waves DENYS DUTYKH 1 Chargé de Recherche CNRS 1 LAMA, Université de Savoie 73376 Le Bourget-du-Lac, France Colloque EDP-Normandie DENYS DUTYKH (CNRS LAMA) Dispersive
More informationEmpirical versus Direct sensitivity computation
Empirical versus Direct sensitivity computation Application to the shallow water equations C. Delenne P. Finaud Guyot V. Guinot B. Cappelaere HydroSciences Montpellier UMR 5569 CNRS IRD UM1 UM2 Simhydro
More informationOn the linear stability of one- and two-layer Boussinesq-type Equations for wave propagation over uneven beds
On the linear stability of one- and two-layer Boussinesq-type Equations for wave propagation over uneven beds Gonzalo Simarro Marine Sciences Institute (ICM, CSIC), 83 Barcelona, Spain Alejandro Orfila
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationUpwind Stabilized Finite Element Modelling of Non-hydrostatic Wave Breaking and Run-up
Upwind Stabilized Finite Element Modelling of Non-hydrostatic Wave Breaking and Run-up P. Bacigaluppi, M. Ricchiuto, P. Bonneton RESEARCH REPORT N 8536 May 014 Project-Team BACCHUS ISSN 049-6399 ISRN INRIA/RR--8536--FR+ENG
More informationA two-layer approach to wave modelling
1.198/rspa.4.135 A two-layer approach to wave modelling By Patrick Lynett and Philip L.-F. Liu School of Civil and Environmental Engineering, Cornell University, Ithaca, NY 14853, USA (plynett@civil.tamu.edu)
More informationarxiv: v1 [math.na] 27 Jun 2017
Behaviour of the Serre Equations in the Presence of Steep Gradients Revisited J.P.A. Pitt a,, C. Zoppou a, S.G. Roberts a arxiv:706.08637v [math.na] 27 Jun 207 a Mathematical Sciences Institute, Australian
More informationImproved Performance in Boussinesq-type Equations
Improved Performance in Boussinesq-type Equations Andrew B. Kennedy, James T. Kirby 1 & Mauricio F. Gobbi 2 Abstract In this paper, simple but effective techniques are used to improve the performance of
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationProceedings of the 7 th International Conference on HydroScience and Engineering Philadelphia, USA September 10-13, 2006 (ICHE 2006) ISBN:
Proceedings of the 7 th International Conference on HydroScience and Engineering Philadelphia, USA September 1-13, 6 (ICHE 6) ISBN: 97744745 Drexel University College of Engineering Drexel E-Repository
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Mathematics, Vol.4, No.3, 6, 39 5. ANTI-DIFFUSIVE FINITE DIFFERENCE WENO METHODS FOR SHALLOW WATER WITH TRANSPORT OF POLLUTANT ) Zhengfu Xu (Department of Mathematics, Pennsylvania
More informationDepth-integrated, non-hydrostatic model for wave breaking and run-up
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids (8) Published online in Wiley InterScience (www.interscience.wiley.com). DOI: 1.1/fld.19 Depth-integrated, non-hydrostatic
More informationModified Serre Green Naghdi equations with improved or without dispersion
Modified Serre Green Naghdi equations with improved or without dispersion DIDIER CLAMOND Université Côte d Azur Laboratoire J. A. Dieudonné Parc Valrose, 06108 Nice cedex 2, France didier.clamond@gmail.com
More informationNUMERICAL SIMULATION OF LONG WAVE RUNUP ON A SLOPING BEACH*
NUMERICAL SIMULATION OF LONG WAVE RUNUP ON A SLOPING BEACH* * presented at Long Waves Symposium (in parallel with the XXX IAHR Congress) August 5-7, 003, AUTh, Thessaloniki, Greece. by HAKAN I. TARMAN
More informationSome Benchmark Simulations for Flash Flood Modelling
Some Benchmark Simulations for Flash Flood Modelling Ekkehard Holzbecher, Ahmed Hadidi German Univ. of Technology in Oman (GUtech) E-mail: ekkehard.holzbecher@gutech.edu.om Flash Floods Rapid flooding
More informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More informationc Copyright 2014 Jihwan Kim
c Copyright 2014 Jihwan Kim Finite Volume Methods for Tsunamis Generated by Submarine Landslides Jihwan Kim A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationWhat is a flux? The Things We Does Know
What is a flux? Finite Volume methods (and others) (are based on ensuring conservation by computing the flux through the surfaces of a polyhedral box. Either the normal component of the flux is evaluated
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationTransformation of irregular waves in the inner surf zone
Transformation of irregular waves in the inner surf zone Philippe Bonneton and Hélène Dupuis 1 Abstract A numerical model based on a high order non-oscillatory MacCormack TVD scheme is presented for the
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
More informationTIME DOMAIN COMPARISONS OF MEASURED AND SPECTRALLY SIMULATED BREAKING WAVES
TIME DOMAIN COMPARISONS OF MEASRED AND SPECTRAY SIMATED BREAKING WAVES Mustafa Kemal Özalp 1 and Serdar Beji 1 For realistic wave simulations in the nearshore zone besides nonlinear interactions the dissipative
More informationA Strategy for the Development of Coupled Ocean- Atmosphere Discontinuous Galerkin Models
A Strategy for the Development of Coupled Ocean- Atmosphere Discontinuous Galerkin Models Frank Giraldo Department of Applied Math, Naval Postgraduate School, Monterey CA 93943 Collaborators: Jim Kelly
More informationAn efficient numerical method for hydraulic transient computations M. Ciccotelli," S. Sello," P. Molmaro& " CISE Innovative Technology, Segrate, Italy
An efficient numerical method for hydraulic transient computations M. Ciccotelli," S. Sello," P. Molmaro& " CISE Innovative Technology, Segrate, Italy Abstract The aim of this paper is to present a new
More informationA Flexible 2D Nonlinear Approach for Nonlinear Wave Propagation, Breaking and Run up
A Flexible 2D Nonlinear Approac for Nonlinear Wave Propagation, Breaking and Run up Andrea Gilberto Filippini, Maria Kazolea, Mario Ricciuto To cite tis version: Andrea Gilberto Filippini, Maria Kazolea,
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationThe RAMSES code and related techniques 2- MHD solvers
The RAMSES code and related techniques 2- MHD solvers Outline - The ideal MHD equations - Godunov method for 1D MHD equations - Ideal MHD in multiple dimensions - Cell-centered variables: divergence B
More informationYanlin Shao 1 Odd M. Faltinsen 2
Yanlin Shao 1 Odd M. Faltinsen 1 Ship Hydrodynamics & Stability, Det Norsk Veritas, Norway Centre for Ships and Ocean Structures (CeSOS), NTNU, Norway 1 The state-of-the-art potential flow analysis: Boundary
More informationIntroduction to Finite Volume projection methods. On Interfaces with non-zero mass flux
Introduction to Finite Volume projection methods On Interfaces with non-zero mass flux Rupert Klein Mathematik & Informatik, Freie Universität Berlin Summerschool SPP 1506 Darmstadt, July 09, 2010 Introduction
More informationLecture 12: Transcritical flow over an obstacle
Lecture 12: Transcritical flow over an obstacle Lecturer: Roger Grimshaw. Write-up: Erinna Chen June 22, 2009 1 Introduction The flow of a fluid over an obstacle is a classical and fundamental problem
More informationDesign of optimal Runge-Kutta methods
Design of optimal Runge-Kutta methods David I. Ketcheson King Abdullah University of Science & Technology (KAUST) D. Ketcheson (KAUST) 1 / 36 Acknowledgments Some parts of this are joint work with: Aron
More informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationA NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS
A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS HASEENA AHMED AND HAILIANG LIU Abstract. High resolution finite difference methods
More informationA class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, 1549 1560 www.scichina.com math.scichina.com www.springerlink.com A class of the fourth order finite volume Hermite weighted essentially
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationWell-balanced central finite volume methods for the Ripa system
Well-balanced central finite volume methods for the Ripa system R. Touma a C. Klingenberg b a Lebanese American University, Beirut, Lebanon b Würzburg University, Würzburg, Germany This paper is dedicated
More informationAvailable online at Eng. Math. Lett. 2014, 2014:17 ISSN: WAVE ATTENUATION OVER A SUBMERGED POROUS MEDIA I.
Available online at http://scik.org Eng. Math. Lett. 04, 04:7 ISSN: 049-9337 WAVE ATTENUATION OVER A SUBMERGED POROUS MEDIA I. MAGDALENA Industrial and Financial Mathematics Research Group, Faculty of
More informationGeneration of undular bores and solitary wave trains in fully nonlinear shallow water theory
Generation of undular bores and solitary wave trains in fully nonlinear shallow water theory Gennady El 1, Roger Grimshaw 1 and Noel Smyth 2 1 Loughborough University, UK, 2 University of Edinburgh, UK
More informationModeling and simulation of bedload transport with viscous effects
Introduction Modeling and simulation of bedload transport with viscous effects E. Audusse, L. Boittin, M. Parisot, J. Sainte-Marie Project-team ANGE, Inria; CEREMA; LJLL, UPMC Université Paris VI; UMR
More informationReceived 6 August 2005; Accepted (in revised version) 22 September 2005
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol., No., pp. -34 Commun. Comput. Phys. February 6 A New Approach of High OrderWell-Balanced Finite Volume WENO Schemes and Discontinuous Galerkin Methods for a
More informationA well-balanced scheme for the shallow-water equations with topography or Manning friction.
A well-balanced scheme for the shallow-water equations with topography or Manning friction. Victor Michel-Dansac a,, Christophe Berthon a, Stéphane Clain b, Françoise Foucher a,c a Laboratoire de Mathématiques
More informationStabilized Spectral Element Approximation of the Saint Venant System Using the Entropy Viscosity Technique
Stabilized Spectral Element Approximation of the Saint Venant System Using the Entropy Viscosity Technique R. Pasquetti, J.L. Guermond, and B. Popov Abstract We consider the Saint Venant system (shallow
More informationHIGH ORDER NUMERICAL METHODS FOR TIME DEPENDENT HAMILTON-JACOBI EQUATIONS
June 6, 7 :7 WSPC/Lecture Notes Series: 9in x 6in chapter HIGH ORDER NUMERICAL METHODS FOR TIME DEPENDENT HAMILTON-JACOBI EQUATIONS Chi-Wang Shu Division of Applied Mathematics, Brown University Providence,
More informationFinite volume method on unstructured grids
Finite volume method on unstructured grids Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationIntroduction to Heat and Mass Transfer. Week 9
Introduction to Heat and Mass Transfer Week 9 補充! Multidimensional Effects Transient problems with heat transfer in two or three dimensions can be considered using the solutions obtained for one dimensional
More informationThe shallow water equations Lecture 8. (photo due to Clark Little /SWNS)
The shallow water equations Lecture 8 (photo due to Clark Little /SWNS) The shallow water equations This lecture: 1) Derive the shallow water equations 2) Their mathematical structure 3) Some consequences
More informationA parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows.
A parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows Tao Xiong Jing-ei Qiu Zhengfu Xu 3 Abstract In Xu [] a class of
More informationRankine-Hugoniot-Riemann solver for steady multidimensional conservation laws with source terms
Rankine-Hugoniot-Riemann solver for steady multidimensional conservation laws with source terms Halvor Lund a,b,c,, Florian Müller a, Bernhard Müller b, Patrick Jenny a a Institute of Fluid Dynamics, ETH
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationEntropy stable schemes for compressible flows on unstructured meshes
Entropy stable schemes for compressible flows on unstructured meshes Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore deep@math.tifrbng.res.in http://math.tifrbng.res.in/
More informationHigh-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
More informationFLUX DIFFERENCE SPLITTING AND THE BALANCING OF SOURCE TERMS AND FLUX GRADIENTS
FLUX DIFFERENCE SPLITTING AND THE BALANCING OF SOURCE TERMS AND FLUX GRADIENTS M.E.Hubbard The University of Cambridge, DAMTP, Silver Street, Cambridge, CB3 9EW, United Kingdom and P.Garcia-Navarro Area
More informationAnti-diffusive finite difference WENO methods for shallow water with. transport of pollutant
Anti-diffusive finite difference WENO methods for shallow water with transport of pollutant Zhengfu Xu 1 and Chi-Wang Shu 2 Dedicated to Professor Qun Lin on the occasion of his 70th birthday Abstract
More informationarxiv: v4 [physics.flu-dyn] 23 Mar 2019
High-Order Localized Dissipation Weighted Compact Nonlinear Scheme for Shock- and Interface-Capturing in Compressible Flows Man Long Wong a, Sanjiva K. Lele a,b arxiv:7.895v4 [physics.flu-dyn] 3 Mar 9
More informationEntropy stable schemes for degenerate convection-diffusion equations
Entropy stable schemes for degenerate convection-diffusion equations Silvia Jerez 1 Carlos Parés 2 ModCompShock, Paris 6-8 Decmber 216 1 CIMAT, Guanajuato, Mexico. 2 University of Malaga, Málaga, Spain.
More informationA Fully Coupled Model of Non-linear Wave in a Harbor
Copyright 2013 Tech Science Press CMES, vol.91, no.4, pp.289-312, 2013 A Fully Coupled Model of Non-linear Wave in a Harbor Daguo Wang 1 Abstract: A 2-D time-domain numerical coupled model for non-linear
More informationA Comparative Study of Divergence-Cleaning Techniques for Multi-Dimensional MHD Schemes )
A Comparative Study of Divergence-Cleaning Techniques for Multi-Dimensional MHD Schemes ) Takahiro MIYOSHI and Kanya KUSANO 1) Hiroshima University, Higashi-Hiroshima 739-856, Japan 1) Nagoya University,
More informationEarth System Modeling Domain decomposition
Earth System Modeling Domain decomposition Graziano Giuliani International Centre for Theorethical Physics Earth System Physics Section Advanced School on Regional Climate Modeling over South America February
More informationGODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS
GODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS ARTURO S. LEON Dept. of Civil and Envir. Engng., Univ. of Illinois at Urbana-Champaign, 2519 Hydrosystems Lab., MC-250. 205 North Mathews Av., Urbana,
More informationx n+1 = x n f(x n) f (x n ), n 0.
1. Nonlinear Equations Given scalar equation, f(x) = 0, (a) Describe I) Newtons Method, II) Secant Method for approximating the solution. (b) State sufficient conditions for Newton and Secant to converge.
More informationTsunami modeling. Philip L-F. Liu Class of 1912 Professor School of Civil and Environmental Engineering Cornell University Ithaca, NY USA
Tsunami modeling Philip L-F. Liu Class of 1912 Professor School of Civil and Environmental Engineering Cornell University Ithaca, NY USA PASI 2013: Tsunamis and storm surges Valparaiso, Chile January 2-13,
More informationScalable Non-Linear Compact Schemes
Scalable Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory International Conference on Spectral and High Order Methods
More informationA Non-Hydrostatic Shallow Water Model on Triangular Meshes in Sam(oa) 2 Guided Research Project
A Non-Hydrostatic Shallow Water Model on Triangular Meshes in Sam(oa 2 Guided Research Project Philipp Samfaß Fakulät für Informatik Technische Universität München Supervisor: Univ.-Prof. Dr. Michael Bader
More informationHydraulic Modelling for Drilling Automation
Hydraulic Modelling for Drilling Automation CASA Day Harshit Bansal April 19, 2017 Where innovation starts Team: Supervisors at TU/e : W.H.A. Schilders, N. van de Wouw, B. Koren, L. Iapichino Collaborators:
More informationSpectral element schemes for the. Korteweg-de Vries and Saint-Venant equations
Spectral element schemes for the Korteweg-de Vries and Saint-Venant equations R. PASQUETTI a a. Université Côte d Azur, CNRS, Inria, LJAD, France, richard.pasquetti@unice.fr... Résumé : Les sytèmes hyperboliques
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
More informationGreen-Naghdi type solutions to the Pressure Poisson equation with Boussinesq Scaling ADCIRC Workshop 2013
Green-Naghdi type solutions to the Pressure Poisson equation with Boussinesq Scaling ADCIRC Workshop 2013 Aaron S. Donahue*, Joannes J. Westerink, Andrew B. Kennedy Environmental Fluid Dynamics Group Department
More informationForce analysis of underwater object with supercavitation evolution
Indian Journal of Geo-Marine Sciences Vol. 42(8), December 2013, pp. 957-963 Force analysis of underwater object with supercavitation evolution B C Khoo 1,2,3* & J G Zheng 1,3 1 Department of Mechanical
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationA Scalable, Parallel Implementation of Weighted, Non-Linear Compact Schemes
A Scalable, Parallel Implementation of Weighted, Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory SIAM Annual Meeting
More informationCranfield ^91. College of Aeronautics Report No.9007 March The Dry-Bed Problem in Shallow-Water Flows. E F Toro
Cranfield ^91 College of Aeronautics Report No.9007 March 1990 The Dry-Bed Problem in Shallow-Water Flows E F Toro College of Aeronautics Cranfield Institute of Technology Cranfield. Bedford MK43 OAL.
More informationMooring Model for Barge Tows in Lock Chamber
Mooring Model for Barge Tows in Lock Chamber by Richard L. Stockstill BACKGROUND: Extensive research has been conducted in the area of modeling mooring systems in sea environments where the forcing function
More informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationA Well-Balanced Gas-Kinetic Scheme for the Shallow-Water Equations with Source Terms
Journal of Computational Physics 178, 533 562 (2002) doi:10.1006/jcph.2002.7040, available online at http://www.idealibrary.com on A Well-Balanced Gas-Kinetic Scheme for the Shallow-Water Equations with
More informationTwo-layer shallow water system and its applications
Proceedings of Symposia in Applied Mathematics Two-layer shallow water system and its applications Jihwan Kim and Randall J. LeVeque Abstract. The multi-layer shallow water system is derived by depth averaging
More informationIntroduction to BASEMENT Basic Simulation Environment for Computation of Environmental Flow and Natural Hazard Simulation
Introduction to BASEMENT Basic Simulation Environment for Computation of Environmental Flow and Natural Hazard Simulation Numerical Hydraulics Autumn semester 2016 Prof. Dr. Markus Holzner Author: Pascal
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationA primal-dual mixed finite element method. for accurate and efficient atmospheric. modelling on massively parallel computers
A primal-dual mixed finite element method for accurate and efficient atmospheric modelling on massively parallel computers John Thuburn (University of Exeter, UK) Colin Cotter (Imperial College, UK) AMMW03,
More informationRotation & nonlinear effects in shallow water: the Rossby adjustment problem
Rotation & nonlinear effects in shallow water: the Rossby adjustment problem Ross Tulloch April 29, 2004 1 1 Shallow water equations Consider, as in [6], a homogeneous incompressible 2D fluid in a channel
More informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
More informationPhysical Diffusion Cures the Carbuncle Phenomenon
Physical Diffusion Cures the Carbuncle Phenomenon J. M. Powers 1, J. Bruns 1, A. Jemcov 1 1 Department of Aerospace and Mechanical Engineering University of Notre Dame, USA Fifty-Third AIAA Aerospace Sciences
More informationAlessandro Valiani & Valerio Caleffi
An engineering approach to the source terms treatment in shallow water equations Alessandro Valiani & Valerio Caleffi Dipartimento di Ingegneria Università degli Studi di Ferrara - ITALY Roscoff, September
More informationarxiv: v2 [math.na] 4 Mar 2014
Finite Difference Weighted Essentially Non-Oscillatory Schemes with Constrained Transport for Ideal Magnetohydrodynamics arxiv:1309.3344v [math.na] 4 Mar 014 Andrew J. Christlieb a, James A. Rossmanith
More informationNUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE
NUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE Kevin G. Wang (1), Patrick Lea (2), and Charbel Farhat (3) (1) Department of Aerospace, California Institute of Technology
More information