Lines Matching refs:E
45 PetscScalar E; /* radiation energy */ member
121 dEm->E = 0; in RDMaterialEnergy_Saha()
132 dEm->E = 0; in RDMaterialEnergy_Reduced()
141 dsigma_R->E = 0; in RDSigma_R()
152 denom = 3. * rd->rho * sigma_R + (int)limit * PetscAbsScalar(nx->E) / n->E; in RDDiffusionCoefficient()
153 ddenom.E = -(int)limit * PetscAbsScalar(nx->E) / PetscSqr(n->E); in RDDiffusionCoefficient()
155 dxdenom.E = (int)limit * (PetscRealPart(nx->E) < 0 ? -1. : 1.) / n->E; in RDDiffusionCoefficient()
159 dD_R->E = -rd->c / PetscSqr(denom) * ddenom.E; in RDDiffusionCoefficient()
163 dxD_R->E = -rd->c / PetscSqr(denom) * dxdenom.E; in RDDiffusionCoefficient()
182 …,%10.2G)\n", i, (double)PetscRealPart(x[i].E), (double)PetscRealPart(x[i].T), (double)PetscRealPar… in RDStateView()
183 … (double)PetscRealPart(f[i].E), (double)PetscRealPart(f[i].T))); in RDStateView()
194 …ma_p / n->T, tmp = 4. * rd->sigma_b * PetscSqr(PetscSqr(n->T)) / rd->c - n->E, tmp_E = -1., tmp_T … in RDRadiation()
196 dn->E = rad_E; in RDRadiation()
208 n_L.E = 0.5 * (x[i - 1].E + x[i].E); in RDDiffusion()
210 nx_L.E = (x[i].E - x[i - 1].E) / hx; in RDDiffusion()
213 fluxL = D_L * nx_L.E; in RDDiffusion()
214 dfluxL[0].E = -ihx * D_L + (0.5 * dD_L.E - ihx * dxD_L.E) * nx_L.E; in RDDiffusion()
215 dfluxL[1].E = +ihx * D_L + (0.5 * dD_L.E + ihx * dxD_L.E) * nx_L.E; 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()
219 n_R.E = 0.5 * (x[i].E + x[i + 1].E); in RDDiffusion()
221 nx_R.E = (x[i + 1].E - x[i].E) / hx; in RDDiffusion()
224 fluxR = D_R * nx_R.E; in RDDiffusion()
225 dfluxR[0].E = -ihx * D_R + (0.5 * dD_R.E - ihx * dxD_R.E) * nx_R.E; in RDDiffusion()
226 dfluxR[1].E = +ihx * D_R + (0.5 * dD_R.E + ihx * dxD_R.E) * nx_R.E; 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()
231 d[0].E = -ihx * dfluxL[0].E; in RDDiffusion()
233 d[1].E = ihx * (dfluxR[0].E - dfluxL[1].E); in RDDiffusion()
235 d[2].E = ihx * dfluxR[1].E; in RDDiffusion()
341 Em_t = dEm.E * xdot[i].E + dEm.T * xdot[i].T; in RDIFunction_FD()
351 n.E = (1. - bcTheta) * x0[0].E + bcTheta * x[0].E; in RDIFunction_FD()
353 nx.E = ((1. - bcTheta) * (x0[1].E - x0[0].E) + bcTheta * (x[1].E - x[0].E)) / hx; in RDIFunction_FD()
358 f[0].E = hx * (n.E - 2. * D_R * nx.E - rd->Eapplied); in RDIFunction_FD()
361 f[0].E = x[1].E - x[0].E; in RDIFunction_FD()
367 f[i].E = x[i].E - x[i - 1].E; /* Homogeneous Neumann */ in RDIFunction_FD()
370 f[i].E = hx * (xdot[i].E - diff - rad); in RDIFunction_FD()
409 dEm_t.E = dEm1.E / (Theta * dt); in RDIJacobian_FD()
417 n1.E = x[i].E; in RDIJacobian_FD()
425 dEm_t.E = dEm.E * a; in RDIJacobian_FD()
436 n.E = (1. - bcTheta) * x0[0].E + bcTheta * x[0].E; in RDIJacobian_FD()
438 nx.E = ((1. - bcTheta) * (x0[1].E - x0[0].E) + bcTheta * (x[1].E - x[0].E)) / hx; in RDIJacobian_FD()
443 … = (bcTheta / Theta) * hx * (1. - 2. * D * (-1. / hx) - 2. * nx.E * dD.E + 2. * nx.E * dxD.E / hx); in RDIJacobian_FD()
444 K[0][1 * 2 + 1] = (bcTheta / Theta) * hx * (-2. * nx.E * dD.T); in RDIJacobian_FD()
445 … + 0] = (bcTheta / Theta) * hx * (-2. * D * (1. / hx) - 2. * nx.E * dD.E - 2. * nx.E * dxD.E / hx); in RDIJacobian_FD()
461 K[0][0 * 2 + 0] = -hx * ddiff[0].E; in RDIJacobian_FD()
463 K[0][1 * 2 + 0] = hx * (a - ddiff[1].E - drad.E); in RDIJacobian_FD()
465 K[0][2 * 2 + 0] = -hx * ddiff[2].E; in RDIJacobian_FD()
469 K[1][1 * 2 + 0] = hx * (rho * dEm_t.E + drad.E); in RDIJacobian_FD()
491 n->E = 0; in RDEvaluate()
493 nx->E = 0; in RDEvaluate()
496 n->E += interp[q][j] * x[i + j].E; in RDEvaluate()
498 nx->E += deriv[q][j] * x[i + j].E; in RDEvaluate()
666 Em_t = dEm.E * nt.E + dEm.T * nt.T; in RDIFunction_FE()
671 …f[i + j].E += (deriv[q][j] * weight[q] * ((1. - Theta) * D0_R * n0x.E + Theta * D_R * nx.E) + inte… in RDIFunction_FE()
683 n.E = (1 - bcTheta) * x0[0].E + bcTheta * x[0].E; in RDIFunction_FE()
685 nx.E = (x[1].E - x[0].E) / hx; in RDIFunction_FE()
692 f[0].E += -ratio * 0.5 * (rd->Eapplied - n.E); in RDIFunction_FE()
746 …j] * weight[q] * (a - drad.E) * interp[q][k] + deriv[q][j] * weight[q] * ((D_R + dxD_R.E * nx.E) *… in RDIJacobian_FE()
747 … interp[q][k]) + deriv[q][j] * weight[q] * (dxD_R.T * deriv[q][k] + dD_R.T * interp[q][k]) * nx.E); in RDIJacobian_FE()
748 K[j * 2 + 1][k * 2 + 0] += interp[q][j] * weight[q] * drad.E * interp[q][k]; in RDIJacobian_FE()
762 n.E = (1 - Theta) * x0[0].E + Theta * x[0].E; in RDIJacobian_FE()
764 nx.E = (x[1].E - x[0].E) / hx; in RDIJacobian_FE()
790 static PetscScalar RDRadiationTemperature(RD rd, PetscScalar E) in RDRadiationTemperature() argument
792 return PetscPowScalar(E * rd->c / (4. * rd->sigma_b), 0.25); in RDRadiationTemperature()
808 x[i].E = 0.001; in RDInitialState()
809 x[i].T = RDRadiationTemperature(rd, x[i].E); in RDInitialState()
812 x[i].E = 0.001 + 100. * PetscExpReal(-PetscSqr(coord / 0.1)); in RDInitialState()
813 x[i].T = RDRadiationTemperature(rd, x[i].E); in RDInitialState()
816 x[i].E = 7.56e-2 * rd->unit.Joule / PetscPowScalarInt(rd->unit.meter, 3); in RDInitialState()
817 x[i].T = RDRadiationTemperature(rd, x[i].E); in RDInitialState()
857 for (i = 0; i < m; i++) y[i] = RDRadiationTemperature(rd, x[i].E); in RDView()
879 n.E = 1.; in RDTestDifferentiation()
882 n.E = 1. + epsilon; in RDTestDifferentiation()
885 fdEm.E = (Em1 - Em0) / epsilon; in RDTestDifferentiation()
886 n.E = 1.; in RDTestDifferentiation()
890 …lPart(dEm.E), (double)PetscRealPart(dEm.T), (double)PetscRealPart(fdEm.E), (double)PetscRealPart(f… in RDTestDifferentiation()
896 n.E = 1.; in RDTestDifferentiation()
898 nx.E = 1.; in RDTestDifferentiation()
901 n.E = 1. + epsilon; in RDTestDifferentiation()
903 nx.E = 1.; in RDTestDifferentiation()
906 fdD.E = (D - D0) / epsilon; in RDTestDifferentiation()
907 n.E = 1; in RDTestDifferentiation()
909 nx.E = 1.; in RDTestDifferentiation()
913 n.E = 1; in RDTestDifferentiation()
915 nx.E = 1. + epsilon; in RDTestDifferentiation()
918 fdxD.E = (D - D0) / epsilon; in RDTestDifferentiation()
919 n.E = 1; in RDTestDifferentiation()
921 nx.E = 1.; in RDTestDifferentiation()
925 …cRealPart(dD.E), (double)PetscRealPart(dD.T), (double)PetscRealPart(fdD.E), (double)PetscRealPart(… in RDTestDifferentiation()
927 …lPart(dxD.E), (double)PetscRealPart(dxD.T), (double)PetscRealPart(fdxD.E), (double)PetscRealPart(f… in RDTestDifferentiation()
936 n0[0].E = 1.; in RDTestDifferentiation()
938 n0[1].E = 5.; in RDTestDifferentiation()
940 n0[2].E = 4.; in RDTestDifferentiation()
945 n1[i].E += epsilon; in RDTestDifferentiation()
946 fd[i].E = (RDDiffusion(rd, hx, n1, 1, 0) - a0) / epsilon; in RDTestDifferentiation()
950 …f {%g,%g}\n", i, (double)PetscRealPart(d[i].E), (double)PetscRealPart(d[i].T), (double)PetscRealPa… in RDTestDifferentiation()
951 … (double)PetscRealPart(d[i].E - fd[i].E), (double)PetscRealPart(d[i].T - fd[i].T))); in RDTestDifferentiation()
957 n.E = 1.; in RDTestDifferentiation()
960 n.E = 1. + epsilon; in RDTestDifferentiation()
963 fdrad.E = (rad - rad0) / epsilon; in RDTestDifferentiation()
964 n.E = 1.; in RDTestDifferentiation()
968 …rt(drad.E), (double)PetscRealPart(drad.T), (double)PetscRealPart(fdrad.E), (double)PetscRealPart(f… in RDTestDifferentiation()