Pore Process

[discoleo]

Pore Process

Generate a coupled process with pores:

• based on the details of function rgrid.channel (rliniar.gen);
• generate n linear channels;
• each channel has a probability pblock to be blocked at a random position across the length of the channel;
• each wall of the channel contains a number of potential pores;
• position of pores across wall:

  npos = sample(seq(w), round(w * ppore)), where:

  • w = width of material = length of channel;
  • round(w * ppore) = number of pores;

• those positions will contain the probability of the pore beeing open or closed:

  channel = rep(-1, w);\ channel[npos] = runif(length(npos));

The process will generate 1 specific instance of the material by applying the cutoff on the probability of the pores being open or closed.