Lines Matching refs:T
46 PetscScalar T; /* material temperature */ member
115 …PetscScalar Em, alpha, alpha_t, T = n->T, T_t = 1., chi = rd->I_H / (rd->k * T), chi_t = -chi / T … in RDMaterialEnergy_Saha() local
118 Em = rd->k * T / rd->m_p * (1.5 * (1. + alpha) + alpha * chi); /* Eq 6 */ in RDMaterialEnergy_Saha()
122 dEm->T = Em / T * T_t + rd->k * T / rd->m_p * (1.5 * alpha_t + alpha_t * chi + alpha * chi_t); in RDMaterialEnergy_Saha()
128 …PetscScalar alpha, alpha_t, T = n->T, T_t = 1., chi = -0.3 / T, chi_t = -chi / T * T_t, a = 1., a_… in RDMaterialEnergy_Reduced() local
130 if (Em) *Em = (1. + alpha) * T + 0.3 * alpha; in RDMaterialEnergy_Reduced()
133 dEm->T = alpha_t * T + (1. + alpha) * T_t + 0.3 * alpha_t; in RDMaterialEnergy_Reduced()
140 *sigma_R = rd->K_R * rd->rho * PetscPowScalar(n->T, -rd->gamma); in RDSigma_R()
142 dsigma_R->T = -rd->gamma * (*sigma_R) / n->T; in RDSigma_R()
154 ddenom.T = 3. * rd->rho * dsigma_R.T; in RDDiffusionCoefficient()
156 dxdenom.T = 0; in RDDiffusionCoefficient()
160 dD_R->T = -rd->c / PetscSqr(denom) * ddenom.T; in RDDiffusionCoefficient()
164 dxD_R->T = -rd->c / PetscSqr(denom) * dxdenom.T; in RDDiffusionCoefficient()
182 …lPart(x[i].E), (double)PetscRealPart(x[i].T), (double)PetscRealPart(xdot[i].E), (double)PetscRealP… in RDStateView()
183 … (double)PetscRealPart(f[i].E), (double)PetscRealPart(f[i].T))); in RDStateView()
194 …T, -rd->beta), sigma_p_T = -rd->beta * sigma_p / n->T, tmp = 4. * rd->sigma_b * PetscSqr(PetscSqr(… in RDRadiation()
197 dn->T = rad_T; in RDRadiation()
209 n_L.T = 0.5 * (x[i - 1].T + x[i].T); in RDDiffusion()
211 nx_L.T = (x[i].T - x[i - 1].T) / hx; in RDDiffusion()
216 dfluxL[0].T = (0.5 * dD_L.T - ihx * dxD_L.T) * nx_L.E; in RDDiffusion()
217 dfluxL[1].T = (0.5 * dD_L.T + ihx * dxD_L.T) * nx_L.E; in RDDiffusion()
220 n_R.T = 0.5 * (x[i].T + x[i + 1].T); in RDDiffusion()
222 nx_R.T = (x[i + 1].T - x[i].T) / hx; in RDDiffusion()
227 dfluxR[0].T = (0.5 * dD_R.T - ihx * dxD_R.T) * nx_R.E; in RDDiffusion()
228 dfluxR[1].T = (0.5 * dD_R.T + ihx * dxD_R.T) * nx_R.E; in RDDiffusion()
232 d[0].T = -ihx * dfluxL[0].T; in RDDiffusion()
234 d[1].T = ihx * (dfluxR[0].T - dfluxL[1].T); in RDDiffusion()
236 d[2].T = ihx * dfluxR[1].T; in RDDiffusion()
341 Em_t = dEm.E * xdot[i].E + dEm.T * xdot[i].T; in RDIFunction_FD()
345 f[i].T = hx * (rho * Em_t + rad); in RDIFunction_FD()
352 n.T = (1. - bcTheta) * x0[0].T + bcTheta * x[0].T; in RDIFunction_FD()
354 nx.T = ((1. - bcTheta) * (x0[1].T - x0[0].T) + bcTheta * (x[1].T - x[0].T)) / hx; in RDIFunction_FD()
410 dEm_t.T = dEm1.T / (Theta * dt); in RDIJacobian_FD()
412 const PetscScalar epsilon = x[i].T * PETSC_SQRT_MACHINE_EPSILON; in RDIJacobian_FD()
418 n1.T = x[i].T + epsilon; in RDIJacobian_FD()
423 Em_TT = (dEm1.T - dEm.T) / epsilon; in RDIJacobian_FD()
426 dEm_t.T = dEm.T * a + Em_TT * xdot[i].T; in RDIJacobian_FD()
437 n.T = (1. - bcTheta) * x0[0].T + bcTheta * x[0].T; in RDIJacobian_FD()
439 nx.T = ((1. - bcTheta) * (x0[1].T - x0[0].T) + bcTheta * (x[1].T - x[0].T)) / hx; in RDIJacobian_FD()
444 K[0][1 * 2 + 1] = (bcTheta / Theta) * hx * (-2. * nx.E * dD.T); in RDIJacobian_FD()
462 K[0][0 * 2 + 1] = -hx * ddiff[0].T; in RDIJacobian_FD()
464 K[0][1 * 2 + 1] = hx * (-ddiff[1].T - drad.T); in RDIJacobian_FD()
466 K[0][2 * 2 + 1] = -hx * ddiff[2].T; in RDIJacobian_FD()
470 K[1][1 * 2 + 1] = hx * (rho * dEm_t.T + drad.T); in RDIJacobian_FD()
492 n->T = 0; in RDEvaluate()
494 nx->T = 0; in RDEvaluate()
497 n->T += interp[q][j] * x[i + j].T; in RDEvaluate()
499 nx->T += deriv[q][j] * x[i + j].T; in RDEvaluate()
666 Em_t = dEm.E * nt.E + dEm.T * nt.T; in RDIFunction_FE()
672 f[i + j].T += interp[q][j] * weight[q] * (rho * Em_t + rad); in RDIFunction_FE()
684 n.T = (1 - bcTheta) * x0[0].T + bcTheta * x[0].T; in RDIFunction_FE()
686 nx.T = (x[1].T - x[0].T) / hx; in RDIFunction_FE()
747 …= (+interp[q][j] * weight[q] * (-drad.T * interp[q][k]) + deriv[q][j] * weight[q] * (dxD_R.T * der… in RDIJacobian_FE()
749 …K[j * 2 + 1][k * 2 + 1] += interp[q][j] * weight[q] * (a * rd->rho * dEm.T + drad.T) * interp[q][k… in RDIJacobian_FE()
763 n.T = (1 - Theta) * x0[0].T + Theta * x[0].T; in RDIJacobian_FE()
765 nx.T = (x[1].T - x[0].T) / hx; in RDIJacobian_FE()
809 x[i].T = RDRadiationTemperature(rd, x[i].E); in RDInitialState()
813 x[i].T = RDRadiationTemperature(rd, x[i].E); in RDInitialState()
817 x[i].T = RDRadiationTemperature(rd, x[i].E); in RDInitialState()
880 n.T = T0; in RDTestDifferentiation()
883 n.T = T0; in RDTestDifferentiation()
887 n.T = T1; in RDTestDifferentiation()
889 fdEm.T = (Em1 - Em0) / (T0 * epsilon); in RDTestDifferentiation()
890 …scRealPart(dEm.E), (double)PetscRealPart(dEm.T), (double)PetscRealPart(fdEm.E), (double)PetscRealP… in RDTestDifferentiation()
891 (double)PetscRealPart(dEm.T - fdEm.T))); in RDTestDifferentiation()
897 n.T = 1.; in RDTestDifferentiation()
899 n.T = 1.; in RDTestDifferentiation()
902 n.T = 1.; in RDTestDifferentiation()
904 n.T = 1.; in RDTestDifferentiation()
908 n.T = 1. + epsilon; in RDTestDifferentiation()
910 n.T = 1.; in RDTestDifferentiation()
912 fdD.T = (D - D0) / epsilon; in RDTestDifferentiation()
914 n.T = 1.; in RDTestDifferentiation()
916 n.T = 1.; in RDTestDifferentiation()
920 n.T = 1.; in RDTestDifferentiation()
922 n.T = 1. + epsilon; in RDTestDifferentiation()
924 fdxD.T = (D - D0) / epsilon; in RDTestDifferentiation()
925 …PetscRealPart(dD.E), (double)PetscRealPart(dD.T), (double)PetscRealPart(fdD.E), (double)PetscRealP… in RDTestDifferentiation()
926 (double)PetscRealPart(dD.T - fdD.T))); in RDTestDifferentiation()
927 …scRealPart(dxD.E), (double)PetscRealPart(dxD.T), (double)PetscRealPart(fdxD.E), (double)PetscRealP… in RDTestDifferentiation()
928 (double)PetscRealPart(dxD.T - fdxD.T))); in RDTestDifferentiation()
937 n0[0].T = 1.; in RDTestDifferentiation()
939 n0[1].T = 3.; in RDTestDifferentiation()
941 n0[2].T = 2.; in RDTestDifferentiation()
948 n1[i].T += epsilon; in RDTestDifferentiation()
949 fd[i].T = (RDDiffusion(rd, hx, n1, 1, 0) - a0) / epsilon; in RDTestDifferentiation()
950 …ealPart(d[i].E), (double)PetscRealPart(d[i].T), (double)PetscRealPart(fd[i].E), (double)PetscRealP… in RDTestDifferentiation()
951 … (double)PetscRealPart(d[i].E - fd[i].E), (double)PetscRealPart(d[i].T - fd[i].T))); in RDTestDifferentiation()
958 n.T = 1.; in RDTestDifferentiation()
961 n.T = 1.; in RDTestDifferentiation()
965 n.T = 1. + epsilon; in RDTestDifferentiation()
967 fdrad.T = (rad - rad0) / epsilon; in RDTestDifferentiation()
968 …ealPart(drad.E), (double)PetscRealPart(drad.T), (double)PetscRealPart(fdrad.E), (double)PetscRealP… in RDTestDifferentiation()
969 (double)PetscRealPart(drad.T - fdrad.T))); in RDTestDifferentiation()