Discussion Forum

Hello everybody,

I want to simulate diffusion of pressure for a simple 1D model. The right boundary I set to symmetry. On the left I defined a Dirichelt-type boundary condition simulating constant pressure up to a certain time t0 which is then penetrating into my 1D medium using the diffusion equation.

After the time t0 I stop the constant pressure injection. As a results the pressure drops to zero at the left side of the boundary. But also within the medium the pressure starts to decrease. However due to diffusion the shape of the individual solutions with time are characterized by a certain maximum.

Here comes my problem:
I know want to define a medium diffusivity D depending on these particular maxima of pressure. Until the pressure increases locally at a certain location up to the corresponding pressure maximum the diffusivity should be D1 = 10. After the maximum of pressure has reached and pressure starts to drop I want to freeze the medium diffusivity to another value D2 = 3.

Thus the medium diffusivity should look like:
D = 10*(t<=t0) + ( D1*(t>t0 && p(x,t) <= PMax(x,t) + ...)

But I can't figure out the pressure values of the maxima during evaluation of my problem.

Has onyone a hint for me?