Re: Front-fixing method for moving boundary problem

Hello everyone,

Jonathan, I think I have found a way to do what I want to do. I tested 
it and it works. It probably falls into the "ugly hack" category, but I 
like it because it's really simple and seems to work fine.

So, what I do is I include the line of code where I define the equation 
inside the solution loop along with my self-made and very basic 
Runge-Kutta-Fehlberg ODE solver. Which results in something like:

while res > tolRes:
     grad = (temp.getHarmonicFaceValue()[-1] - 

             temp.getHarmonicFaceValue()[-2]) / dx
     dSdt = intRecRate(elapsed, Snew, grad)
     eq = TransientTerm()
          == ImplicitDiffusionTerm(coeff=K/(c*rho*Snew**2)) \
          - ExponentialConvectionTerm(xi/Snew * dSdt) \
          - temp / Snew * dSdt
     res = eq.sweep(var=temp,
                    boundaryConditions=BCs,
                    dt=dt,
                    solver=solver)
     grad = (temp.getHarmonicFaceValue()[-1]
             - temp.getHarmonicFaceValue()[-2]) / dx
     resODE = RKF45(lambda t, S: intRecRate(t, S, grad), 0, S, dt)

This probably won't beat any records of speed, but I'm not after that 
for the moment.

Also, the use of french in your answers leads me to believe 2 things. 
First, you probably guessed my francophone background from my name 
(sadly, maybe also my non-native english). That is very clever! Second, 
I like to think that you want to learn some french. I would like to 
support this endeavour by correcting a tiny little mistake. The word 
"dommage" is masculin. Therefore we say "Quel dommage!"

Meilleures salutations!

definition of the eu
> 
> 
> On Jul 20, 2008, at 7:35 AM, Etienne Rivard wrote:
> 
>> Right. I remember once considering swithcing to another method like 
>> the phase field method to solve my problem. But I have to admit I was 
>> somewhat demotivated by the apparent complexity of the method. For the 
>> moment, I don't think I'll get deeper into that.
> 
> chacun à son goût
> 
> 
>> Thanks for the offer, I appreciate it but. I have time to work on this 
>> and this actually a part of my PhD work so I need to show a deep 
>> understanding of everything I use, so I would rather shoot for a more 
>> elegant solution.
> 
> Understood. When you get it done, if you wanted to contribute said 
> elegant solution to FiPy, we'd be delighted.
> 
> 
>