R中输入参数随时间的步进变化

时间:2018-05-13 21:15:57

标签: r numerical-methods pde

如果有人可以帮助我如何在输入参数中加入List<Topic> topicList = Ebean.find(Topic.class).where().eq("columnName", "NametoMatch").findList(); topicList.get(0) 的时间。请参阅以下代码:

step

我希望library(ReacTran) N <- 10 # No of grids L = 0.10 # thickness, m l = L/2 # Half of thickness, m k= 0.412 # thermal conductivity, W/m-K cp = 3530 # thermal conductivity, J/kg-K rho = 1100 # density, kg/m3 T_int = 57.2 # Initial temperature , degC T_air = 19 # air temperature, degC h_air = 20 # Convective heat transfer coeff of air, W/m2-K xgrid <- setup.grid.1D(x.up = 0, x.down = l, N = N) x <- xgrid$x.mid alpha.coeff <- (k*3600)/(rho*cp) Diffusion <- function (t, Y, parms){ tran <- tran.1D(C=Y, flux.down = 0, C.up = T_air, a.bl.up = h_air, D = alpha.coeff, dx = xgrid) list(dY = tran$dC, flux.up = tran$flux.up, flux.down = tran$flux.down) } # Initial condition Yini <- rep(T_int, N) times <- seq(from = 0, to = 2, by = 0.2) print(system.time( out <- ode.1D(y = Yini, times = times, func = Diffusion, parms = NULL, dimens = N))) plot(times, out[,(N+1)], type = "l", lwd = 2, xlab = "time, hr", ylab = "Temperature") 在第1小时内保持不变,并且它会更改为剩余T_air的其他值。这将是参数中的1 hr。我该怎么做? 任何帮助,将不胜感激。 谢谢,

0 个答案:

没有答案