Spring, 2025
Jewett
Program #4 -- Test D
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 4.
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 D - same as test C, but with select nonlinear U, W advection terms disabled.
Do NOT DO THIS TEST unless you had problems getting your test C to work! 
     If test C was OK, there is no need to run test case D.
  • Parameters:  all same as test C:
    • 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 C results due to compiler/system changes since testC was run.
  • >> Diffusion on theta used time level (n); we will use (n-1), so these results are no longer correct for us -- expect differences <<<

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

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