Skip to content

CO2_sequestration_v1.0 #84

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
wants to merge 14 commits into
base: develop
Choose a base branch
from

Conversation

askoldvilkha
Copy link

@askoldvilkha askoldvilkha commented Nov 19, 2024

This pull request adds the basic implementation of the CO2 sequestration problem (in 2d) to IGA-ADS, based on the existing heat2d simulation.

Current functionality:

  • CLI requires the user to provide the mesh dimensions (x (horizontal) and y (vertical)), the number of iteration steps, and the timestep size
  • The user also can specify the simulation parameters (density, viscosity, etc.)
  • The simulation parameters have the hard-coded default values but can be changed via the CLI (use ./co2_sequestration --help for more information)
  • Verbose added to update the user with the execution progress
  • Porosity maps are implemented as well. In order to use them, one should provide the path to the file after --porosity-map. The format should correspond to the ones in the example files (in the co2_sequestration directory)

Auxiliary codes:

  • co2_movie.py - makes an animation using the outputs from co2_sequestration executable
  • poros_and_perm_reader.py - reads the source data for the porosity and permeability maps and wraps it in ASCII format suitable for the co2_sequestration executable.

Other updates

  • Updated the dependencies script for the newer GCC compiler versions (GCC>=11)
  • Added a script to load the necessary modules for compiling the code in CYFRONET
  • Added a test .sub script for checking if the job submission system works

…This version is working but will need many corrections in the next updates
…ode; enforced proper 2D mesh (user-defined rectangle instead of a square)
@askoldvilkha askoldvilkha force-pushed the develop branch 5 times, most recently from c622025 to d0fac9e Compare February 9, 2025 15:08
askoldvilkha and others added 10 commits February 9, 2025 16:09
    1) added scripts for easier setup of the IGA-ADS software on Cyfronet cluster along with a sample test `.sub` script
    2) updated `CMakeLists.txt` to resolve the conflict with main repo
….py` script for making the viewable animation from the `co2_sequestration` executable output
…YFRONET ARES

Changes:
1) the gas injection location and rate are now user defined but not hardcoded. If no input is provided - location will be the mesh center, rate - 1e-6.
2) This version is still operating within uniform porosity/permeability space - porosity and permeability maps integration is still to be implemented
3) `co2_movie.py` now allows the user to define fps rate.
4) `ares_test.sub` - specified test script for submission on ARES. The following changes will be local unless they are major. Should be in the same directory with the executable `co2_sequestration` and the script `co2_movie.py`
…th the code to allow simple intrepolation and ASCII wrapping
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant