Functions
	cavity (N, Ra, Pr=0.71, mu=0.05, dt=8.0, steps=3000, tol=1e-5, verbose=False)

Variables
dict	REF = {1e3: 2.243}

dict	NU_REF = {1e3: 1.118, 1e4: 2.243, 1e5: 4.519}

	r = cavity(N, Ra)

dict	e = abs(r['Nu'] - NU_REF[Ra]) / NU_REF[Ra] * 100

Detailed Description

Differentially heated square cavity (de Vahl Davis 1983) — validates the Boussinesq
field->momentum coupling (property closures + per-cell body force + scalar transport).

Left wall hot T=1, right wall cold T=0, top/bottom adiabatic, no-slip everywhere, gravity in -y
with a Boussinesq buoyancy body force. Benchmark average Nusselt number on the hot wall:

    Ra    Nu_avg   u_max*   v_max*
    1e3   1.118    3.649    3.697
    1e4   2.243   16.178   19.617
    1e5   4.519   34.73    68.59

(* velocity extrema normalised by alpha/L.) Run:  PYTHONPATH=<build> python dvd_cavity.py

Function Documentation

◆ cavity()

dvd_cavity.cavity	(	N,
		Ra,
		Pr = `0.71`,
		mu = `0.05`,
		dt = `8.0`,
		steps = `3000`,
		tol = `1e-5`,
		verbose = `False`
	)

Definition at line 22 of file dvd_cavity.py.

Variable Documentation

◆ REF

dict dvd_cavity.REF = {1e3: 2.243}

Definition at line 18 of file dvd_cavity.py.

◆ NU_REF

dict dvd_cavity.NU_REF = {1e3: 1.118, 1e4: 2.243, 1e5: 4.519}

Definition at line 19 of file dvd_cavity.py.

◆ r

dvd_cavity.r = cavity(N, Ra)

Definition at line 56 of file dvd_cavity.py.

◆ e

dict dvd_cavity.e = abs(r['Nu'] - NU_REF[Ra]) / NU_REF[Ra] * 100

Definition at line 60 of file dvd_cavity.py.

Functions

Variables