Jump to content
RockWare Support Forum

T.C. Onstott

  • Content count

  • Joined

  • Last visited

Community Reputation

0 Neutral

About T.C. Onstott

  • Rank
  1. N2O

    We would like to use our measured N2O concentrations with GWB to calculate microbial nitrogen cycling and was wondering if anyone has produced a set of logK values for an N2O(g) reaction???
  2. Dear Tom, I would like to grow biomass using the microbe reaction and monitor depletion of critical carbon substrate, N (i.e. NH3(a) and P (HPO4--) to determine when an environment becomes substrate limited as opposed to energy limited. I am not sure how to do this. I created a mineral phase bioacetate(s) with the formula C36O40H102N9P1 and comprised of 4 species, Acetic_acid(aq), NH3(aq), HPO4-- and H+. this work fine. I then have a microbe reaction that uses Acetic_acid(aq) and calculates biomass. What I then tried to do is to create a custom rate law for C36O40H102N9P1 where the rate is equal to (biomass-biomass0)/Deltat. I thought that the custom rate laws could import these values like they do in the ratelaw_custom4.bas example, but no luck. I just get a can't evaluate the rate error message. Before I beat my head against a wall, I thought I would check in to see if you can offer some advice. If my code reads as follows: kinetic microbe-24 rxn = "Acetate + H2O -> Methane(aq) + HCO3-" rate_law = ratelaw_example5.bas rate_con = 9e-10 kinetic microbe-24 biomass = 2e-14 kinetic microbe-24 growth_yield = 1 kinetic microbe-24 ATP_energy = -6 kinetic microbe-24 ATP_number = 1 kinetic microbe-24 Decay_con = 1e-11 kinetic microbe-24 power(Acetate) = 1 kinetic microbe-24 powerD(Acetate) = 1 kinetic microbe-24 order1 = 1 kinetic microbe-24 order2 = 1 kinetic microbe-24 KD = 8.6e-4 react 0.1 grams of "bioacetate(s)" kinetic "bioacetate(s)" rate_law = ratelaw_bioacetate2.bas rate_con = 3.18e-8 surface = 1 do I need to explicitly state something like kinetic "bioacetate(s)" biomass0 = ??? kinetic "bioacetate(s)" biomass = ??? in order for this ratelaw to be evaluated? Ideally I want to use the same biomass given by microbe-24. Thanks, tco
  3. Custom Rate Laws

    sure thing Tom, attached is the React code. I have highlighted the call command down about 3/4 of the way through the program.
  4. Custom Rate Laws

    Dear Tom, I have tried unsuccessfully to implement a very simple rate law command by embedding it in my react code using brackets, quotation marks and single quotes and by trying to call it as a basic program. Nothing seems to work. The problem I suspect with trying to call the basic subroutine is that it cannot recognize my text file. I've tried using the PRINT command to print the "rate" value being returned and it doesn't seem to recognize the "rate =" even though that is how the manual says I should format the command. Could you give me some guidance on how to implement my custom rate laws? thanks, tco I've attached the short rate law program and also the longer REACT program that calls it at react H2(s). This species is one that I have defined in my customize thermo database and works just fine with REACT. I'm using Vers. 8 Standard running on Windows on a MAC (as a separate sector on my hard drive). no problems otherwise. ratelaw_H2.bas.txt
  5. missing H2O

    Greetings, I've been running a combined flush with reaction program and have come across a perplexing problem. I noticed my porosity plummeting and crashing the simulation. At first I thought it must be due to mineral precipitation, but as I examined the output more closely I noticed that my solvent mass was decreasing, even though I was flushing, which I've assumed meant I was constant replenishing solvent mass. I couldn't identify any mineral dissolution or formation reaction that could account for the solvent mass loss. So I'm flummoxed. I'm using React 7.04, PC platform. see attached code and output. React_output.txt newGullfaksSRB100.txt
  6. flush, reactant times and execution time

    Thanks Tom, I was afraid you were going to say that. If I could just get someone to pony up the $$ for it, I'd purchase it in a New York second. cheers, tco
  7. flush, reactant times and execution time

    I have been using the flush command to model flow through an aquifer and in order to simulate flow rates of .01 m/yr I have had to increase my reactant times value 100000 or higher. the execution time increases because the step size decreases as I increase reactant times so that the simulation requires 24 hours on my laptop. my question: is there any way that I can increase the step size in order to reduce the execution time. the delxi command has no affect on the step size and change the volume of either fluid also has no affect. thanks, tco I'm using version 7 of React on a PC. newhighGullfaksIRB70.txt