Problem loading initial temperature conditions from file

[Vale]

Hello everyone!

I have a simple rectangular model with EOS1.

In the initial conditions panel, I need to import temperatures from a file. That is because I have not a linear gradient, but a polynomial one (third order).

So my idea was to import a .TXT file with depth versus temperature as initial conditions, since petrasim does not allow polynomial functions. I've tried to insert a linear function but the temperature does not really fix with the real ones.

I wrote a simple .txt file that looks like this: (depth vs T)

1500 15

500 25

-300 100

-6000 400

When I try to load the file, I get an error "Unknown Error: No Data Loaded".

I've also tried to invert the columns, T vs depth but the result is the same.

I used space, tab, comma to separate the columns, but nothing has worked.

Which is the right way to write the .txt file to import in the initial conditions?

I used absolute values for the elevations (meters a.s.l) and °C for the temperature.

Thank you

Best regards

[Vale]

Hello everyone, it's me again.

I'm still experiencing this problem.

I don't really know how to import a temperature file to load into the initial conditions.

I've tried different ways but without success.

Any help would be really appreciated,

Best regards

Vale

[Alison]

We hadn’t responded to this one because there is really no easy way to do this. However, because we've gotten more requests here, I thought that I would go ahead and post the following suggestion.I was able to establish the desired initial conditions using the following instructions. However, it would be important for you the experiment with adding boundary conditions such as fixed state cells, heat sources, etc. to be sure that these initial conditions are realistic steady state conditions.

1. Create the model and set the top layer as a fixed state layer. Run the model for several years to establish a hydrostatic pressure gradient.

2. Save the SIM file with a new name (and in a new folder), and then load the SAVE file created during step 1 as Initial conditions

3. Select the layers for which you have desired temperatures and make them fixed state, with a constant pressure and temperature. You can select a cell from the layer to determine the pressure that you should assign to the layer. Its important for you to assign a constant pressure to each model layer.

4. You will need to lower the permeability of the material in the model to a very low value, otherwise, fluid movement caused by convection will created convection cells.

5. Run the model for a very long time (for example, 1m years), and you should end up with results that mimic the desired temperature gradient (with linear gradients in-between the fixed state layers.

I hope that this helps. Again, this isn’t really the best way to do this, but if it might help speed up the modeling process. I think that the ideal way to establish these types of conditions would be for you to build a model that reproduces the process that created the conditions in the first place (for example, the instruction of heat sources or fixed T cells at the bottom of the model, combined with fixed state cells at the edges of the geothermal feature being modeled. This is the recommendation of the TOUGH2 authors.

Best Regards,

Alison