Menu
▾
▴
Re: [Anuga-user] Actual height of wall using breaklines
|
From: Gareth D. <gar...@gm...> - 2019-01-30 07:36:26
|
Hi,
Are you passing riverWall_Par as an argument to
domain.riverwallData.create_riverwalls? It's not enough to just define
it in the script.
In the script you mention, the "riverWall_Par" line is commented out,
and nothing is passed to create_riverwalls. I guess that is because the
intention was to use the default behaviour.
If you want to change the defaults, you need to pass the argument in to
the create_riverwalls function.
For further info, see the documentation in
anuga_core/anuga/structures/riverwall.py
Cheers,
G
On 30/1/19 1:59 pm, Ananth Wuppukondur wrote:
> Hi Gareth,
>
> Thanks for promptly responding. Yes, it is modified from
> 'runMerewether.py'. I've tried both the methods now.
>
> 1. By widening the polygon to higher than mesh resolution so triangles
> can be accommodated. But the result is still the same. So i tried
> using 'riverWall'.
>
> 2. Using 'riverWall', the simulation runs and it does seem to generate
> a wall (although I'm yet to compare the results with experiments as
> i'm awaiting complete simulation). When trying to set the default
> riverWall parameters using
> "riverWall_Par={'centralWall':{'Qfactor':1.0, 's1': 0.999,
> 's2':0.9999, 'h1':100, 'h2':150}}" as in 'weir_1.py" example, the
> simulation shows default values are being used instead of the
> specified values. Could you please suggest how to set these values?
> Thanks in advance!
>
> On Tue, Jan 29, 2019 at 7:57 PM Gareth Davies
> <gar...@gm... <mailto:gar...@gm...>> wrote:
>
> Hi Ananth,
>
> From your code snippet, it looks like your "breaklines" are
> actually "polygons" with a given specified height. I suspect this
> has been modified from the runMerewether.py example? In the latter
> script, 'breaklines' are used to force edge boundaries in the mesh
> (defining buildings in that case). Then we set the elevation
> inside the buildings to a 'house-height'. It looks like your code
> is similar.
>
> I'm not 100% sure -- but I suspect this is not really what you
> want to do. If it really is what you want, then be careful that
> you are not 'accidently smoothing' the elevation (i.e. be sure to
> set the elevation quantity at centroids, not vertices). Also, be
> careful that your "wall polygons" are wide enough to include a
> continuous set of triangles.
>
> An alternative approach is to keep using "breaklines" to ensure
> that mesh-edges follow some specified line (i.e. your walls).
> Next, instead of defining the wall using a group of cells, you can
> use a "riverwall" operator to tell ANUGA what the "edge-elevation"
> should be. In this case, the "riverwall" is like a 'thin
> discontinuity in elevation' between 2 triangles. Flow over this
> "weir" can be computed using some mixture of a weir equation, and
> the Riemann flux function (for details see the documentation of
> riverWalls in the code).
>
> An example which uses the riverwall functionality like what I
> describe is here:
> https://github.com/GeoscienceAustralia/anuga_core/tree/master/validation_tests/behaviour_only/weir_1
>
> Probably that's the easiest place to look to get started. However,
> there are a bunch of other examples which use the riverwall
> functionality in the same subdirectory. It is also used in a
> generic script here:
> https://github.com/GeoscienceAustralia/anuga_core/tree/master/examples/cairns_excel.
>
>
> Cheers,
> Gareth.
>
>
> On 29/1/19 7:34 pm, Ananth Wuppukondur wrote:
>> Hi,
>>
>> I'm using the following to define elevation of a wall using
>> 'breaklines'.
>>
>> for i in range(len(breaklines)):
>> breaklines[i] = breaklines[i] + [breaklines[i][0]]
>> house_addition_poly_fun_pairs.append(
>> [ breaklines[i], house_height])
>> house_addition_poly_fun_pairs.append(['All', 0.])
>>
>> The simulation runs without any error. However, the results are
>> not as expected when comparing to experiments. I suspect this
>> difference in experiments and simulations could be because of
>> some issue in translating wall elevation from the above function.
>> I've two cases: 1. With wall defined by breaklines, 2. without
>> any wall. The second case matches well with corresponding
>> experimental results. But not the first case. Is there any
>> command to display the elevation assigned to the breaklines? Or
>> would you recommend any alternative for defining a wall? Thanks
>> in advance.
>>
>> --
>> --
>> Best wishes,
>> *Ananth Wuppukondur,*
>> PhD candidate | Coastal Research Group
>> School of Civil Engineering | The University of Queensland
>> Brisbane, Australia.
>> /Web-links:/Personal webpage
>> <https://ananthwsharma.wixsite.com/ananthwuppukondur> Google
>> scholar
>> <https://scholar.google.ca/citations?user=AROgxX4AAAAJ&hl=en>
>> ResearchGate
>> <https://www.researchgate.net/profile/Ananth_Wuppukondur>LinkedIn
>> <https://in.linkedin.com/in/ananth-wuppukondur-a37a6048>
>>
>>
>>
>> _______________________________________________
>> Anuga-user mailing list
>> Anu...@li... <mailto:Anu...@li...>
>> https://lists.sourceforge.net/lists/listinfo/anuga-user
>
>
> _______________________________________________
> Anuga-user mailing list
> Anu...@li...
> <mailto:Anu...@li...>
> https://lists.sourceforge.net/lists/listinfo/anuga-user
>
>
>
> --
> --
> Best wishes,
> *Ananth Wuppukondur,*
> PhD candidate | Coastal Research Group
> School of Civil Engineering | The University of Queensland
> Brisbane, Australia.
> /Web-links:/Personal webpage
> <https://ananthwsharma.wixsite.com/ananthwuppukondur> Google scholar
> <https://scholar.google.ca/citations?user=AROgxX4AAAAJ&hl=en>
> ResearchGate
> <https://www.researchgate.net/profile/Ananth_Wuppukondur>LinkedIn
> <https://in.linkedin.com/in/ananth-wuppukondur-a37a6048>
>
>
>
> _______________________________________________
> Anuga-user mailing list
> Anu...@li...
> https://lists.sourceforge.net/lists/listinfo/anuga-user
|
