Hello,
I am working on developing a time dependent simulation of a system of pde's that I have already set up. The system of pde's has 32 variables to be solved for. I am having trouble defining the initial conditions for it. The initial condition that I would like to impose on the simulation is the stationary state of my pde system plus a small perturbation on only one of the variables (this variable is meant to activate other variables in the system). The perturbation is supposed to be a small deviation, say of around 1% of the original value, on only a specific patch of the domain in which the variable is defined. How could I go about this? I've looked at some threads that mentioned using the "withsol" operator but I have not found this function anywhere.