Skip to content

Commit

Permalink
[T/THM] Water freezing volumetric expansion
Browse files Browse the repository at this point in the history
Benchmark to check the volumetric expansion due to water-to-ice phase change.
  • Loading branch information
Tymofiy Gerasimov authored and endJunction committed Aug 31, 2023
1 parent 0e87889 commit 9d6f6f4
Show file tree
Hide file tree
Showing 14 changed files with 834 additions and 0 deletions.
1 change: 1 addition & 0 deletions ProcessLib/ThermoHydroMechanics/Tests.cmake
Original file line number Diff line number Diff line change
Expand Up @@ -6,6 +6,7 @@ if (NOT OGS_USE_MPI)
OgsTest(PROJECTFILE ThermoHydroMechanics/1D_freezing_column_Stefan/Stefan_problem.prj RUNTIME 15)
OgsTest(PROJECTFILE ThermoHydroMechanics/ColumnDeformationFreezing/TM.prj RUNTIME 13)
OgsTest(PROJECTFILE ThermoHydroMechanics/HeatingHomogeneousDomain/hex_THM.prj RUNTIME 30)
OgsTest(PROJECTFILE ThermoHydroMechanics/9percentWaterFreezingExpansion/UnitSquare.prj RUNTIME 1)
endif()

AddTest(
Expand Down
Original file line number Diff line number Diff line change
@@ -0,0 +1,46 @@
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
<UnstructuredGrid>
<FieldData>
<DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="166" format="appended" RangeMin="34" RangeMax="125" offset="0" />
<DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="26" format="appended" RangeMin="45" RangeMax="121" offset="176" />
<DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0" RangeMax="0" offset="268" />
<DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0" RangeMax="0" offset="328" />
</FieldData>
<Piece NumberOfPoints="4" NumberOfCells="1" >
<PointData>
<DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="388" />
<DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="448" />
<DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="508" />
<DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="3" offset="568" />
<DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="640" />
<DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="700" />
<DataArray type="Float64" Name="ice_volume_fraction" format="appended" RangeMin="6.3169798575e-36" RangeMax="6.3169798575e-36" offset="760" />
<DataArray type="Float64" Name="pressure" format="appended" RangeMin="0" RangeMax="0" offset="832" />
<DataArray type="Float64" Name="pressure_interpolated" format="appended" RangeMin="0" RangeMax="0" offset="892" />
<DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="952" />
<DataArray type="Float64" Name="sigma_ice" NumberOfComponents="4" format="appended" RangeMin="1.5630471519e-26" RangeMax="1.5630471519e-26" offset="1012" />
<DataArray type="Float64" Name="temperature" format="appended" RangeMin="277.15" RangeMax="277.15" offset="1104" />
<DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="277.15" RangeMax="277.15" offset="1172" />
</PointData>
<CellData>
<DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="1240" />
<DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="0" offset="1300" />
<DataArray type="Float64" Name="fluid_density_avg" format="appended" RangeMin="1000" RangeMax="1000" offset="1360" />
<DataArray type="Float64" Name="stress_avg" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="1424" />
<DataArray type="Float64" Name="viscosity_avg" format="appended" RangeMin="0.001" RangeMax="0.001" offset="1484" />
</CellData>
<Points>
<DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="0.14142135624" offset="1552" />
</Points>
<Cells>
<DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="1628" />
<DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="1700" />
<DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="1760" />
</Cells>
</Piece>
</UnstructuredGrid>
<AppendedData encoding="base64">
_AQAAAAAAAAAAgAAAAAAAAKYAAAAAAAAAYwAAAAAAAAA=eF6FzDEKgDAMheG7ZO4kTr2KSIgaJWCTktZBxLvb1UXH9374LhCtvDlVMcVsbSG501kgDtcrmi/sELsASokhQpEtEUqG9hxpYkdbcbaUTVlrA/o7/BCci+wtfCLj/QAlSj0gAQAAAAAAAAAAgAAAAAAAABoAAAAAAAAAIgAAAAAAAAA=eF4z0zPRM9E1NDQ21003NU9KSTFJMtVLySwqqQQAUwMHQw==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAADAAAAAAAAAA=eF5jYBhYAAAAgAABAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAADAAAAAAAAAA=eF5jYBhYAAAAgAABAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAACwAAAAAAAAA=eF5jYMAPAAAgAAE=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAACwAAAAAAAAA=eF5jYMAPAAAgAAE=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAADAAAAAAAAAA=eF5jYKAMAAAAQAABAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF5jYIAARijNBKWZoTQAAHAABw==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAADAAAAAAAAAA=eF5jYKAMAAAAQAABAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAADAAAAAAAAAA=eF5jYBhYAAAAgAABAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFQAAAAAAAAA=eF7TaiqWkT29wEIDSuug8QHAtAvrAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAACwAAAAAAAAA=eF5jYMAPAAAgAAE=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAACwAAAAAAAAA=eF5jYMAPAAAgAAE=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAADAAAAAAAAAA=eF5jYBhYAAAAgAABAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAIwAAAAAAAAA=eF4TCF6Wfja8bZcAGs0ABTxQPjoNkxeG8tFpYvUDACnwLeI=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCoEKHNBw0AMiADAU=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCoEKHNBw0AMiADAU=AQAAAAAAAAAAgAAAAAAAAAQAAAAAAAAADAAAAAAAAAA=eF5jYGBgAAAABAABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYIAAAAAIAAE=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADgAAAAAAAAA=eF5jYAACh34HAAImARA=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAACwAAAAAAAAA=eF5jYMAPAAAgAAE=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF77s/LjJd+kAHsAGYAEpw==AQAAAAAAAAAAgAAAAAAAAGAAAAAAAAAAFwAAAAAAAAA=eF5jYMAOZs0EgZ32uMQJycP4AOa2Ej0=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF5jYIAARijNBKWZoTQAAHAABw==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYYAAAAAoAAU=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF7jBAAACgAK
</AppendedData>
</VTKFile>
Original file line number Diff line number Diff line change
@@ -0,0 +1,46 @@
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
<UnstructuredGrid>
<FieldData>
<DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="166" format="appended" RangeMin="34" RangeMax="125" offset="0" />
<DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="26" format="appended" RangeMin="45" RangeMax="121" offset="176" />
<DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="4.1569219455e-05" RangeMax="4.1569219455e-05" offset="268" />
<DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="3.6914466404e-12" RangeMax="3.6914466404e-12" offset="348" />
</FieldData>
<Piece NumberOfPoints="4" NumberOfCells="1" >
<PointData>
<DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="-6.1354140972" RangeMax="-3.0677070486" offset="432" />
<DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="-2.9321533521e-11" RangeMax="-1.466076676e-11" offset="512" />
<DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="7.7850266252e-22" RangeMax="1.3516943364e-21" offset="588" />
<DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="3" offset="696" />
<DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="3.3941125557e-06" offset="768" />
<DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="4.1569219455e-05" RangeMax="4.1569219455e-05" offset="844" />
<DataArray type="Float64" Name="ice_volume_fraction" format="appended" RangeMin="1.4869239894e-18" RangeMax="1.4869239894e-18" offset="952" />
<DataArray type="Float64" Name="pressure" format="appended" RangeMin="0" RangeMax="0" offset="1028" />
<DataArray type="Float64" Name="pressure_interpolated" format="appended" RangeMin="0" RangeMax="0" offset="1088" />
<DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="3.6914466404e-12" RangeMax="3.6914466404e-12" offset="1148" />
<DataArray type="Float64" Name="sigma_ice" NumberOfComponents="4" format="appended" RangeMin="2482606.1563" RangeMax="2482606.1563" offset="1256" />
<DataArray type="Float64" Name="temperature" format="appended" RangeMin="275.15" RangeMax="275.15" offset="1396" />
<DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="275.15" RangeMax="275.15" offset="1464" />
</PointData>
<CellData>
<DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="1532" />
<DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="0" offset="1592" />
<DataArray type="Float64" Name="fluid_density_avg" format="appended" RangeMin="1000" RangeMax="1000" offset="1652" />
<DataArray type="Float64" Name="stress_avg" NumberOfComponents="4" format="appended" RangeMin="3.6914466404e-12" RangeMax="3.6914466404e-12" offset="1716" />
<DataArray type="Float64" Name="viscosity_avg" format="appended" RangeMin="0.001" RangeMax="0.001" offset="1796" />
</CellData>
<Points>
<DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="0.14142135624" offset="1864" />
</Points>
<Cells>
<DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="1940" />
<DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="2012" />
<DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="2072" />
</Cells>
</Piece>
</UnstructuredGrid>
<AppendedData encoding="base64">
_AQAAAAAAAAAAgAAAAAAAAKYAAAAAAAAAYwAAAAAAAAA=eF6FzDEKgDAMheG7ZO4kTr2KSIgaJWCTktZBxLvb1UXH9374LhCtvDlVMcVsbSG501kgDtcrmi/sELsASokhQpEtEUqG9hxpYkdbcbaUTVlrA/o7/BCci+wtfCLj/QAlSj0gAQAAAAAAAAAAgAAAAAAAABoAAAAAAAAAIgAAAAAAAAA=eF4z0zPRM9E1NDQ21003NU9KSTFJMtVLySwqqQQAUwMHQw==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAGgAAAAAAAAA=eF5jvH2hsFjr5z5GKM0Ao6GAkcbyACxCNUk=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAHAAAAAAAAAA=eF4LC11aoLW/aW8YnFZNAdEMUBBGY3kA7kkssQ==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGwAAAAAAAAA=eF77rG39b2UXx4GvYFriwDsofRkqDgAg6xFGAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF5L4Iqaslhuw95kMH0ATidAxQHgnA4/AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAALgAAAAAAAAA=eF5jYHAI/a0405qBgcFWnKkIRJ+47NAOoh3EmZpA9Bkk/m4gHYKkfjcA4tAPOg==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF5jYIAARijNBKWZoTQAAHAABw==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAGAAAAAAAAAA=eF5jYEAFs0SXf/mqeGQfOh+XOADZJxOJAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAALwAAAAAAAAA=eF77c+tCYbHWz31/oPRvKM0ABb+gfBj9A02e4TaUj05DwU+o+p9o5sDkAdCJPgM=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF7T3Vd09kyutY0OlDaC0tpQGgAC9w9rAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAACwAAAAAAAAA=eF5jYMAPAAAgAAE=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAACwAAAAAAAAA=eF5jYMAPAAAgAAE=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAALwAAAAAAAAA=eF4LDl1aoLW/aW8wlA4JVU0B0QxQEAgVR9Co8mFQcQSNKh8AFYfRfmjyAOOJLIY=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAARwAAAAAAAAA=eF47eFQt9Os9U8eDUHoflL4rfO5Vq0iQ1R4oH0YfgNJGn96duSmTv+sIlA+jD0Ppa8uXaUyKStm1F8qH0TDzTaH6AWVFRYU=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5LSwMCo0KHNBw0AMEAC4U=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5LSwMCo0KHNBw0AMEAC4U=AQAAAAAAAAAAgAAAAAAAAAQAAAAAAAAADAAAAAAAAAA=eF5jYGBgAAAABAABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYIAAAAAIAAE=AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADgAAAAAAAAA=eF5jYAACh34HAAImARA=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGQAAAAAAAAA=eF4LC11aoLW/aW8YnFZNAdEMUAAA4ykLLQ==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF77s/LjJd+kAHsAGYAEpw==AQAAAAAAAAAAgAAAAAAAAGAAAAAAAAAAFwAAAAAAAAA=eF5jYMAOZs0EgZ32uMQJycP4AOa2Ej0=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF5jYIAARijNBKWZoTQAAHAABw==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYYAAAAAoAAU=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF7jBAAACgAK
</AppendedData>
</VTKFile>
Loading

0 comments on commit 9d6f6f4

Please sign in to comment.