Spring, 2019
Jewett
Program #5 -- Test E
ATMS 502 / CSE 566
Numerical Fluid Dynamics

Description

This test case is provided for you to check your code before running the assigned problem for Problem 5.
Please note the settings below, and remember to change all settings back to the official problem(s) before running and handing in results.

This is test case E - same as test D, but with nonlinear U, W advection enabled or disabled.
Do NOT DO THIS TEST unless you had problems getting your test D to work! 
     If test D was OK, there is no need to run test case E.
  • Parameters:  all same as test D:
    • Grid layout: 201x41; dx=dz=100 meters
    • Quasi-compressible sound speed Cs = 100 meters/sec
    • Integration: time step ?t = 0.25 sec; run for 2400 steps (T=600s)
    • This test includes full physics, with the initial condition U = W = P' = 0, and T' = a cool bubble in the middle of the domain.
    • Initial bubble location: center (10050, 2050) meters, which is grid point (101,21).
    • Initial bubble specification: center magnitude -15 degrees K, "radii" in X: 4000 m; Z: 1000 m.
    • Diffusion: done for all k=[1,...,nz] for U and Theta (changed from original handout). 
    • Diffusion coefficient = 50 (for U, W) for all cases; diffusion for theta = 40.0 (for all, Lax-W and PL)
  • Note: plots have U and W averaged to the theta/p locations at the center of each grid cell.
  • Results for 'control' differ slightly from test D results due to compiler/system changes since testD was run.

Results: (click to see all times (400-1600 seconds) for each case)

Control (=test D) No U/W advection U but no W advection W but no U advection