如何在R语言中实现输入参数T_air随时间的阶跃变化？

Hey there! To implement a step change for T_air in your heat conduction simulation, you can dynamically adjust its value based on the current simulation time within your Diffusion function. Since the ode.1D solver passes the current time t to the function at each step, we can use this to switch between temperature values seamlessly.

Modified Code with Step Change #

library(ReacTran)
library(deSolve)  # Ensure deSolve is loaded (required for ode.1D)

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 # specific heat capacity, J/kg-K (fixed your comment here!)
rho = 1100 # density, kg/m3
T_int = 57.2 # Initial temperature , degC

# Define step change parameters
T_air_pre_step = 19  # Air temp for first 1 hour
T_air_post_step = 25 # Air temp after 1 hour
step_time = 1        # Step occurs at 1 hour
h_air = 20 # Convective heat transfer coeff of air, W/m2-K

alpha.coeff <- (k*3600)/(rho*cp)
xgrid <- setup.grid.1D(x.up = 0, x.down = l, N = N)
x <- xgrid$x.mid

Diffusion <- function (t, Y, parms){
  # Switch T_air based on current time
  current_T_air <- ifelse(t < step_time, T_air_pre_step, T_air_post_step)
  
  tran <- tran.1D(C=Y, flux.down = 0, C.up = current_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 the temperature at the boundary (first grid point) over time
plot(times, out[,(N+1)], type = "l", lwd = 2, xlab = "time, hr", ylab = "Temperature")
abline(v=step_time, lty=2, col="red", lwd=1.5) # Add vertical line to mark step change

Key Changes Explained #

• Step Parameters: We defined T_air_pre_step, T_air_post_step, and step_time to clearly set up the step change conditions—you can tweak these values to match your actual requirements.
• Dynamic Temperature Switch: Inside the Diffusion function, we use ifelse() to check if the current time t is before or after the step time, then select the corresponding T_air value for the boundary condition.
• Visualization: Added a vertical dashed red line to the plot to clearly mark when the step change occurs, making it easier to verify the effect.
• Fixed Comment: Corrected your comment for cp (it's specific heat capacity, not thermal conductivity) to avoid confusion later.

Notes for Extensions #

If you need multiple step changes (e.g., switch temperatures at 1 hour and 1.5 hours), you can extend this logic using findInterval() to map time ranges to different temperature values:

temp_values <- c(19, 25, 22)
time_breaks <- c(0, 1, 1.5)
current_T_air <- temp_values[findInterval(t, time_breaks)]

内容的提问来源于stack exchange，提问作者KBH