XiaoMao

Hi, I am trying to simulate the weathering process by toughreact and you can find the description of my question. I choose EOS9 because I want to know the rainwater infiltration process. However, actually, I have no idea how to set the CO2 partial pressure. How to keep a fixed CO2 partial pressure if I want to inject rainwater from model top? How to set a CO2 partial pressure in the soil? model.pptx

When I tried to run my model, I got the error message: Time step = 0 Time (yr) = 0.000000E+00 There is a problem with fugacity coefficient I use eos2 to simulate the CO2 injection. I set a higher partial pressure 12 for CO2 at the left bottom coner. All other cell is zero. I can not run my model and I have no idea how to solve this porblem. I can not upload the sim file because it is 69KB. the max total szie only 30KB...

Ok Thank you so much

Can I use a large density , such as 1e50 kg/m3, as a constant boundary?

My question is Why low pH caused by carbon dioxide cannot move forward

I downloaded the example from https://www.thunderheadeng.com/petrasim/petrasim-4-documentation/ Example Problem 1, Linear sorption and decay (EOS9). I want to simulate that the low pH caused by CO2 dissolution transport in the domain. Therefore, I modified the example problem 1 I set two Gas zones in the model. CO2 partial pressure is zero in the first zone and it has been set all cells. and then CO2 partial pressure is 10 in the second zone, and it has been set at middle one cell. However, I got resutls like this, low pH only occurs in the cell which is set high partial pressure and it did not move with groundwater flow. Logically, low pH should migrate to the right. But it doesn't。 Is anyone can tell me the reason ? you can find the sim file in the attachment. Thank you so much. treact_prob_1.sim

Thanks Alison

In my conceptual model, geochemical reaction will produce some clay minerals, and then the clay minerals swelling will clog or block porosity in the reservoir . so, I am not sure if Petrasim can simulate this process?

There are three ways to set a constant boundary in Petrasim, First, set the boundary cells to fixed state; Second, set the boundary cells to 1E50 or more Third, set an extra cell which has large or very large volume, and then connect the extra cell to boundary cell. Actually, I am a beginner for Petrasim, so I am not sure whether these three ways are right. Please let me know if I make mistake, THANKS! My question is that there is a large lake, and my model domain is the bottom of the lake, so what type should I use when I set the top boundary? I focus on pressure, saturation and temperature, I do not know what type boundary is more appropriately. Thank you very much!

wa, Tom Thanks a lot I will try to convert it by gwb. I am doing some model by gwb.. Thansk again. Xiao

I exported permeability vlaues from Plot Line, they should 1E-12, 1E-13 1E-15 However, All of these are Zero. How to solve this problem? Permeabi. 0 0 0 0

OH, can I calculate them by myself ? CO2, H2S and H2 consist of mixed gases. They are injected into my model. I can get them CO2(g), H2S(g), and H2(g), so can I calculate by some method to convert them to aqueous?

I have run my model which describes a process of injection of mixed gases. However, I only get the results of CO2(g), H2S(g). How can I convert them form gas to aqueous? I wish I can get the result of CO2(aq) and H2S(aq).

When I run my model, after 3.4 days, the DT is always 0.01s . Should I wait？ Or How can I adjust model parameters? Do you have any suggestions?

Recently, I would like to model the Mixed gases injection process. I have used a Extra Cell as a injection source. However, I do not know how to determine the Air Mass Fraction. Actually, I want to inject a hot mixed gases into a aquifer, and I have known the temperature, but I do not know the partial pressure. Only I known is that the concentration of gas in final water, CO2 is 17mmol/kg. So, How can I set the initial condition at the Extra Cell? When I set the Air Mass Fraction is 0.5 or 0.55, it will stop. The model can run if I set the Air Mass Fraction is 0.6, but it can not reach the final time I set.