Re: [Anuga-user] Anuga infiltration fuction

 Re: [Anuga-user] Anuga infiltration fuction From: Stephen Roberts - 2014-02-06 21:25:14 Attachments: Message as HTML ```Hi Mike, We should be able to put together an operator (or a forcing term) to do what you want. We have a rate operator which takes a quantity or a function of x,y,t as a rate. See run_rate_spatial_operator.py in the anuga/operators Should be able to use that as a base for your operator. By the way, in your formula, cumulative infiltration Z = k * T^a + f0 * T + c, what is T, time? For the rate operator, we will need an instantaneous rate, which I assume would be the function R = ka T^{a-1} + f0 Will have to think about how to deal with c, but maybe we need a quantity to store this "crack fill volume" and just remove water until this quantity reaches c. We could build this into your operator. Cheers Steve On 5/02/2014 4:52 PM, Mike.Morris@... wrote: > Hello all > > Quite some time ago I posted a query to the forum on the suitability > of anuga for modelling surface irrigation events, with the aim of > using anuga to develop better irrigation bay designs. We've since > been successful in attracting funding for this work, and have made > good progress with implementing and testing anuga for the purpose. > > In an irrigation context, infiltration is a major flux (in fact its > the flux we are most interested in). To date we've been using the > Rainfall forcing term with a negative rate (thanks for the tip, > Steve), which effectively provides us with a constant infiltration > rate, but we are now at the stage where an approach is required that > can represent time varying infiltration rates that start when each > cell in the mesh becomes wet. The thing we are unsure of is how best > to implement a function that will have a different start time in each > mesh cell. > > 1D surface irrigation models in common use implement infiltration with > an empirical Kostiakov equation - cumulative infiltration Z = k * T^a > + f0 * T + c, where k and a are empirical , f0 is a final infiltration > rate and c is an initial crack-fill/rapid adsorption term. I imagine > that the same approach could work in anuga, and that we could use our > 1D model calibration routine to provide estimates for the Kostiakov > coefficients for input to anuga. > > Has anyone attempted developing an infiltration routine for anuga in > the 10 months since I last posted? If not, do you have advice on its > development? > > Thanks in advance! > > Mike > > Please note that my working days are Monday, Tuesday, Wednesday. > > *Mike Morris | Senior Research Scientist, Irrigation Hydrology | > Agriculture Research* > Department of Environment and Primary Industries | 255 Ferguson Rd, > Tatura, VIC 3616 > T: 03 5833 5283 | M: 0419 317 454 | E: mike.morris@... > | DX: 219130 > Notice: > > This email and any attachments may contain information that is personal, > confidential, legally privileged and/or copyright. No part of it should be > reproduced, adapted or communicated without the prior written consent of the > copyright owner. > > It is the responsibility of the recipient to check for and remove viruses. > > If you have received this email in error, please notify the sender by return > email, delete it from your system and destroy any copies. You are not authorised > to use, communicate or rely on the information contained in this email. > > Please consider the environment before printing this email. > > > ------------------------------------------------------------------------------ > Managing the Performance of Cloud-Based Applications > Take advantage of what the Cloud has to offer - Avoid Common Pitfalls. > Read the Whitepaper. > http://pubads.g.doubleclick.net/gampad/clk?id=121051231&iu=/4140/ostg.clktrk > > > _______________________________________________ > Anuga-user mailing list > Anuga-user@... > https://lists.sourceforge.net/lists/listinfo/anuga-user -- +--------------------------------------------------------------------+ | Steve Roberts | My Office: (61)(2) 6125 4445 | | Department of Mathematics | Math Office: (61)(2) 6125 2908 | | Mathematical Sciences Institute| Fax: (61)(2) 6125 4984 | | John Dedman Building #27 | mailto:stephen.roberts@... | | Australian National University | http://www.maths.anu.edu.au/~steve| | Canberra ACT 0200 AUSTRALIA | ANU CRICOS # 00120C | +--------------------------------------------------------------------+ ```

 [Anuga-user] Anuga infiltration fuction From: Mike.M - 2014-02-05 06:06:57 Attachments: Message as HTML ```Hello all Quite some time ago I posted a query to the forum on the suitability of anuga for modelling surface irrigation events, with the aim of using anuga to develop better irrigation bay designs. We've since been successful in attracting funding for this work, and have made good progress with implementing and testing anuga for the purpose. In an irrigation context, infiltration is a major flux (in fact its the flux we are most interested in). To date we've been using the Rainfall forcing term with a negative rate (thanks for the tip, Steve), which effectively provides us with a constant infiltration rate, but we are now at the stage where an approach is required that can represent time varying infiltration rates that start when each cell in the mesh becomes wet. The thing we are unsure of is how best to implement a function that will have a different start time in each mesh cell. 1D surface irrigation models in common use implement infiltration with an empirical Kostiakov equation - cumulative infiltration Z = k * T^a + f0 * T + c, where k and a are empirical , f0 is a final infiltration rate and c is an initial crack-fill/rapid adsorption term. I imagine that the same approach could work in anuga, and that we could use our 1D model calibration routine to provide estimates for the Kostiakov coefficients for input to anuga. Has anyone attempted developing an infiltration routine for anuga in the 10 months since I last posted? If not, do you have advice on its development? Thanks in advance! Mike Please note that my working days are Monday, Tuesday, Wednesday. Mike Morris | Senior Research Scientist, Irrigation Hydrology | Agriculture Research Department of Environment and Primary Industries | 255 Ferguson Rd, Tatura, VIC 3616 T: 03 5833 5283 | M: 0419 317 454 | E: mike.morris@... | DX: 219130 Notice: This email and any attachments may contain information that is personal, confidential, legally privileged and/or copyright. No part of it should be reproduced, adapted or communicated without the prior written consent of the copyright owner. It is the responsibility of the recipient to check for and remove viruses. If you have received this email in error, please notify the sender by return email, delete it from your system and destroy any copies. You are not authorised to use, communicate or rely on the information contained in this email. Please consider the environment before printing this email.```
 Re: [Anuga-user] Anuga infiltration fuction From: Stephen Roberts - 2014-02-06 21:25:14 Attachments: Message as HTML ```Hi Mike, We should be able to put together an operator (or a forcing term) to do what you want. We have a rate operator which takes a quantity or a function of x,y,t as a rate. See run_rate_spatial_operator.py in the anuga/operators Should be able to use that as a base for your operator. By the way, in your formula, cumulative infiltration Z = k * T^a + f0 * T + c, what is T, time? For the rate operator, we will need an instantaneous rate, which I assume would be the function R = ka T^{a-1} + f0 Will have to think about how to deal with c, but maybe we need a quantity to store this "crack fill volume" and just remove water until this quantity reaches c. We could build this into your operator. Cheers Steve On 5/02/2014 4:52 PM, Mike.Morris@... wrote: > Hello all > > Quite some time ago I posted a query to the forum on the suitability > of anuga for modelling surface irrigation events, with the aim of > using anuga to develop better irrigation bay designs. We've since > been successful in attracting funding for this work, and have made > good progress with implementing and testing anuga for the purpose. > > In an irrigation context, infiltration is a major flux (in fact its > the flux we are most interested in). To date we've been using the > Rainfall forcing term with a negative rate (thanks for the tip, > Steve), which effectively provides us with a constant infiltration > rate, but we are now at the stage where an approach is required that > can represent time varying infiltration rates that start when each > cell in the mesh becomes wet. The thing we are unsure of is how best > to implement a function that will have a different start time in each > mesh cell. > > 1D surface irrigation models in common use implement infiltration with > an empirical Kostiakov equation - cumulative infiltration Z = k * T^a > + f0 * T + c, where k and a are empirical , f0 is a final infiltration > rate and c is an initial crack-fill/rapid adsorption term. I imagine > that the same approach could work in anuga, and that we could use our > 1D model calibration routine to provide estimates for the Kostiakov > coefficients for input to anuga. > > Has anyone attempted developing an infiltration routine for anuga in > the 10 months since I last posted? If not, do you have advice on its > development? > > Thanks in advance! > > Mike > > Please note that my working days are Monday, Tuesday, Wednesday. > > *Mike Morris | Senior Research Scientist, Irrigation Hydrology | > Agriculture Research* > Department of Environment and Primary Industries | 255 Ferguson Rd, > Tatura, VIC 3616 > T: 03 5833 5283 | M: 0419 317 454 | E: mike.morris@... > | DX: 219130 > Notice: > > This email and any attachments may contain information that is personal, > confidential, legally privileged and/or copyright. No part of it should be > reproduced, adapted or communicated without the prior written consent of the > copyright owner. > > It is the responsibility of the recipient to check for and remove viruses. > > If you have received this email in error, please notify the sender by return > email, delete it from your system and destroy any copies. You are not authorised > to use, communicate or rely on the information contained in this email. > > Please consider the environment before printing this email. > > > ------------------------------------------------------------------------------ > Managing the Performance of Cloud-Based Applications > Take advantage of what the Cloud has to offer - Avoid Common Pitfalls. > Read the Whitepaper. > http://pubads.g.doubleclick.net/gampad/clk?id=121051231&iu=/4140/ostg.clktrk > > > _______________________________________________ > Anuga-user mailing list > Anuga-user@... > https://lists.sourceforge.net/lists/listinfo/anuga-user -- +--------------------------------------------------------------------+ | Steve Roberts | My Office: (61)(2) 6125 4445 | | Department of Mathematics | Math Office: (61)(2) 6125 2908 | | Mathematical Sciences Institute| Fax: (61)(2) 6125 4984 | | John Dedman Building #27 | mailto:stephen.roberts@... | | Australian National University | http://www.maths.anu.edu.au/~steve| | Canberra ACT 0200 AUSTRALIA | ANU CRICOS # 00120C | +--------------------------------------------------------------------+ ```