Legado · Validação térmica

Material homogêneo

Corpo de 1m×1m1\,\mathrm{m}\times 1\,\mathrm{m}, isolamento em todo o contorno. A formulação isotérmica, com a Q(teq)Q(t_{eq}) abaixo, tem de devolver a elevação adiabática analítica. A temperatura permanece uniforme.

Artigo: Validação analítica e numérica do problema termo-químico do concreto nas primeiras idades · DOI 10.70271/rirc.v1n2a008

Quadro de dados

Item Valor
Geometria quadrado 1m×1m1\,\mathrm{m}\times 1\,\mathrm{m}
Contorno q=0q=0 em todo o bordo
T0T_0 25∘C25\,^{\circ}\mathrm{C}
TrefT_{\mathrm{ref}} 25∘C25\,^{\circ}\mathrm{C}
Ea/RE_a/R 4000K4000\,\mathrm{K}
CC 300kg/m3300\,\mathrm{kg/m^3}
kk 9,36kJ/(m⋅h⋅∘C)9{,}36\,\mathrm{kJ/(m\cdot h\cdot^{\circ}C)}
ρ\rho 2300kg/m32300\,\mathrm{kg/m^3}
cec_e 0,900kJ/(kg⋅∘C)0{,}900\,\mathrm{kJ/(kg\cdot^{\circ}C)}
Elevação adiabática T(t)=25+30t2/(122+t2)T(t)=25+30\,t^{2}/(12^{2}+t^{2}), tt em horas
Horizonte da série 00 a 480h480\,\mathrm{h}, passo de 1h1\,\mathrm{h}

A geração que entra no modelo é isotérmica, dois termos de Hill, com aa em horas:

ΔTQ(teq)=16,4teq1,93(0,43⋅24)1,93+teq1,93+13,6teq1,57(1,35⋅24)1,57+teq1,57.\Delta T^{Q}(t_{eq})=16{,}4\,\frac{t_{eq}^{1{,}93}}{(0{,}43\cdot 24)^{1{,}93}+t_{eq}^{1{,}93}}+13{,}6\,\frac{t_{eq}^{1{,}57}}{(1{,}35\cdot 24)^{1{,}57}+t_{eq}^{1{,}57}}.

C⋅Q=ρceΔTQC\cdot Q=\rho c_e\,\Delta T^{Q}. A soma 16,4+13,616{,}4+13{,}6 é os 30∘C30\,^{\circ}\mathrm{C} da curva adiabática.

Resultado de referência

A série desta página é a curva adiabática analítica na coluna Tmax_C (o corpo é uniforme, então qualquer ponto serve) e a QQ isotérmica correspondente ao teqt_{eq} integrado sobre essa temperatura.

tt (h) TT (°C) ΔT\Delta T (°C)
12 40,00 15,00
24 49,00 24,00
480 54,98 29,98

O que se confere no MEF: campo uniforme e histórico sobre esta tabela. O relatório não publica um RMSE deste corpo; a comparação é a da figura do comparativo analítico-numérico.

Série t × Tmáx × teq × Q(t)

serie_homo_t_Tmax_teq_Q.csv

Tmax_C é o alvo. Q_kJ_kg é a fonte, já em função de teqt_{eq} da própria linha. Com Ea/R=4000KE_a/R=4000\,\mathrm{K} e Tref=25∘CT_{\mathrm{ref}}=25\,^{\circ}\mathrm{C}, teqt_{eq} adianta em relação a tt porque a temperatura sobe.

Figuras do relatório

Malha.

Malha do corpo homogêneo

PDF da malha

Contorno de temperatura, uniforme.

Contorno de temperatura

PDF do contorno

Histórico no MEF e comparativo com a curva adiabática.

Histórico de temperatura

PDF do histórico

Comparativo analítico-numérico

PDF do comparativo

Passo a passo da interconversão, como está no relatório.

Tabela da interconversão

PDF da tabela

Curvas adiabáticas em outras temperaturas iniciais, usadas para ler Ea/RE_a/R.

Elevação adiabática e atividade térmica