GEOS-DEV / GEOS

GEOS Simulation Framework
GNU Lesser General Public License v2.1
210 stars 84 forks source link

Error in time history outputs #2703

Open beckyggbbh opened 1 year ago

beckyggbbh commented 1 year ago

Describe the bug I every example I played the HDF files doesnt coincide with the integrated test you made. For example in the Problem Multiphase Flow with Wells I get some strange peaks in the graph of the consumption of the production wells and the behavior does not look like what is shown on the page Screenshot 2023-09-22 090613

jhuang2601 commented 1 year ago

@beckyggbbh Are your comparing the results from the case with exactly same input xml files? or comparing the smoke test results and the ones from benchmark test? Can you give more details?

beckyggbbh commented 1 year ago

I run the exact imput file, for this example I run "deadOilSpe10Layers84_85_benchmark.xml " And I compared it with the HDF5 Files in src/docs/sphinx/basicExamples/multiphaseFlowWithWells

Also I run the example "CO2 Plume Evolution and Leakage Through an Abandoned Well" and the results are wrong

WhatsApp Image 2023-09-07 at 1 52 30 AM

jhuang2601 commented 1 year ago

@jafranc Could you please take a look at this issue?

jafranc commented 1 year ago

@jafranc Could you please take a look at this issue?

I ll have a look

paveltomin commented 1 year ago

Same issue here https://github.com/GEOS-DEV/GEOS/discussions/2658 ?

paveltomin commented 1 year ago

Similar thing for another example: https://geosx-geosx.readthedocs-hosted.com/en/latest/docs/sphinx/basicExamples/multiphaseFlowWithWells/Example.html

Plot for old hdf5 files stored in documentation: image

Plot for hdf5 files from current run: image

Don't match and not sure which one is correct...

paveltomin commented 10 months ago

@jafranc any findings here?

jafranc commented 10 months ago

@jafranc any findings here?

I reproduced your last curves and agree with you output frequency-wise it is what we should expect. Played a bit with time step size and mesh size but couldn't reproduced back the first. Couldn't trace back when we change output frequency but that would be a good starting point to id what changed since.

I think the previous problem with TH and MPI sync has been once and for all fixed so that should be smth else.