Discussion Forum

Hello,

I am currently running a boundary value simulation, using the general form PDE, which utilizes a moving mesh to account for a shrinking (or swelling, depending on Initial conditions and Boundary conditions) of circular domain. I have run a swelling case, and wish to use the final conditions of this case as initial conditions on another simulation (i.e. a shrinking case). I am already using a separate solver to extract the solution from the first solver and use the final time states, as initial states for my dependent variables. However, when I initially run the swelling case, I must impose a thin boundary layer at the outer radius, for my initial node conditions so I don't have too steep of a gradient with the environment. In the second solver (shrinking case), I redefine my initial conditions to be the values at the last time (of the initial swell case) and do not include the decaying function to the boundary. But when I run the simulation the same type of thin boundary layer appears at the initial time. Also, I find that when I set up the second solver as above, but do not include a time dependent study in it (i.e. just let it have the initial values, the initial shrinking circle does have the exact distribution as the final swelling circle. Is there any insight to what I may try to do to remedy this?

Thank you for any help