Extended Matlab Set

Advanced Soil Surface Energy Models

?????Cuc02.m?oil02.m?uc03.m?oil03.m嚗誑??solar?abs?WS?OHI?olrad 蝑??獢?

The page groups the air-temperature boundary model, the latent-heat model, the solar-radiation waveform, and the HO/HI reference curves under one bilingual structure.

Matlab sources are preserved in matlab-source/cuc03 for reference.

CUC02 Scenario A

Air boundary, higher radiation

Matches the first half of Cuc02.m with RP=2000, EPSA=0.75, and uniform initial soil temperatures.

RP solar radiation amplitude (kJ/m2/hr)

EPSA air emissivity

dt numerical step (hr)

tfinal simulation time (hr)

initial TF (簞C)

initial T1 (簞C)

initial T2 (簞C)

initial T3 (簞C)

initial T4 (簞C)

CUC02 Scenario B

Air boundary, lower radiation

Matches the second half of Cuc02.m with RP=100, EPSA=0.71, and non-uniform initial temperatures.

RP solar radiation amplitude (kJ/m2/hr)

EPSA air emissivity

dt numerical step (hr)

tfinal simulation time (hr)

initial TF (簞C)

initial T1 (簞C)

initial T2 (簞C)

initial T3 (簞C)

initial T4 (簞C)

CUC03

Latent heat included

Derived from soil03.m, where EPSA depends on dew point, WO is evaluated from outdoor air temperature TO, and the floor energy balance includes latent heat for evaporation.

outside dew point TD (簞C)

RP solar radiation amplitude (kJ/m2/hr)

dt numerical step (hr)

tfinal simulation time (hr)

initial TF (簞C)

initial T1 (簞C)

initial T2 (簞C)

initial T3 (簞C)

initial T4 (簞C)

Layer	Final Temp (簞C)	Minimum (簞C)	Maximum (簞C)	Peak Time (hr)

solrad Reference

Solar radiation waveforms

Recreates the ideas in solrad.m: full sine, half sine, and positive sine radiation curves.

HOHI Reference

Convective coefficients

Plots the empirical HO and HI curves from HOHI.m for wind speed and temperature difference.

Fixed Constants

Values preserved from Matlab

For soil02.m and soil03.m, the browser keeps KS=5.5, CS=2000, HS=25.2, SIG=20.4.
Layer depths remain Z0=0.01, Z1=0.05, Z2=0.1, Z3=0.2, Z4=0.5.
The outdoor temperature driver is still TO = T0 + TU * sin(ω(clk - 8)) and solar input follows RP * sin(ω(clk - 6)) clipped at zero.
KS is soil thermal conductivity, CS is volumetric heat capacity, HS is the ground heat-transfer coefficient, and SIG is the long-wave radiation constant used in the original examples.
T0, TU, TO, and TD are temperatures in °C; clk is time in hr; ω = 2π / 24 is in hr^-1; and RP is solar-radiation amplitude in kJ/m²/hr.

Formula Notes

Key relations and symbols

HO = 18*v^0.576 and HO = 6*v^0.8 are representative outdoor convection curves, where v is wind speed in m/s and HO is in W m^-2 °C^-1.
HI = 4.6*DT^1/3 and HI = 4.36*DT^0.25 are representative indoor convection curves, where DT is temperature difference in °C and HI is in W m^-2 °C^-1.
RP is solar-radiation amplitude (kJ/m²/hr), TD is outside dew point (°C), EPSA is air emissivity (dimensionless), and WO is outdoor humidity ratio (kg/kg).