Long wave propagation: Difference between revisions

From CIRPwiki
Jump to navigation Jump to search
No edit summary
Line 1: Line 1:
= Analytical Solution =
= Analytical Solution =
The purpose of this verification test is to assess the model performance in simulating long wave propagation.  The case is also useful for testing the model performance and symmetry for a non-rectangular domain. The test case consists of a quarter annulus shaped domain with a linear bottom in the radial direction. The analytical
The purpose of this verification test is to assess the model performance in simulating long wave propagation.  The case is also useful for testing the model performance and symmetry for a non-rectangular domain. The test case consists of a quarter annulus with a linear bed elevation in the radial direction. The inner radius is 60.96 km (200,000 ft), and the outer radius is 152.4 km (500,000 ft). The analytical is given by Lynch and Gray (1978).


= Setup =
Table 1. General Settings for Flow over a Bump
{|border="1"
|'''Parameter''' ||'''Value'''
|-
|Deepwater Tidal Amplitude|| 0.3048 m
|-
|Tidal Period || 12.42 hrs (M2)
|-
|Bottom  friction || None
|-
|Horizontal eddy viscosity || None
|}
 
= Model Setup =
A Cartesian grid is used with two levels of refinement from 8 km spacing to 2 km spacing.
 
Table 1. General Settings for Flow over a Bump
{|border="1"
|'''Parameter''' ||'''Value'''
|-
|Time step|| 600 sec (10 min)
|-
|Ramp period  || 24 hrs
|-
|Simulation time  || 120 hrs
|-
|Grid resolution || 1 - 8 km
|-
|Number of active ocean cells || 1160
|}


[[Image:QA_Depth_8km.png|thumb|right|600px| Figure 1. Quarter Annulus]]
[[Image:QA_Depth_8km.png|thumb|right|600px| Figure 1. Quarter Annulus]]


*Forcing: 0.3048 m amplitude M2 tide (red cell string)
*Bottom friction: 0.0 sec/m^(1/3) (Manning N)
*Horizontal eddy viscosity: :0 m^2/sec
*Horizontal eddy viscosity: :0 m^2/sec
*Run duration: 5 days
*Run duration: 5 days

Revision as of 23:24, 3 March 2011

Analytical Solution

The purpose of this verification test is to assess the model performance in simulating long wave propagation. The case is also useful for testing the model performance and symmetry for a non-rectangular domain. The test case consists of a quarter annulus with a linear bed elevation in the radial direction. The inner radius is 60.96 km (200,000 ft), and the outer radius is 152.4 km (500,000 ft). The analytical is given by Lynch and Gray (1978).

Table 1. General Settings for Flow over a Bump

Parameter Value
Deepwater Tidal Amplitude 0.3048 m
Tidal Period 12.42 hrs (M2)
Bottom friction None
Horizontal eddy viscosity None

Model Setup

A Cartesian grid is used with two levels of refinement from 8 km spacing to 2 km spacing.

Table 1. General Settings for Flow over a Bump

Parameter Value
Time step 600 sec (10 min)
Ramp period 24 hrs
Simulation time 120 hrs
Grid resolution 1 - 8 km
Number of active ocean cells 1160
Figure 1. Quarter Annulus
  • Horizontal eddy viscosity: :0 m^2/sec
  • Run duration: 5 days
  • Ramp period: 2 days
  • Time step: 10 min
  • Wall Friction: Off
  • Grid resolution: 1000 x 1000 m (constant)
  • Number of cells: 23409


Results

The figure below shows a time series of water levels at the inner edge of the model domain. The goodness of fit statistics are listed in the below. The model accurately predicts the wave phase but slightly overestimates the amplitude by approximately 0.01 m.

Table 2. Goodness of fit statistics for the water elevation

Statistic Value
RMSE 0.014 m
RMAE 0.028
R^2 0.999
Bias 0.003 m

QA Water Level 8km.png


References

  • Lynch, D.R., and Gray, W.G. (1978). "Analytical solutions for computer flow model testing," J. Hydraulics Division, 104, 1409-28.

Test_Cases