severel questions about the postprocess of results

[Dongxueyang]

Can anybody share your experience of postporcess of results?

I want to get the results with vtk format. Is there any way to complete the transition of formats to vtk ? At the same time I want to get the value of flow(m3/s m3/h) , can I output the flow of one section,or how to calculate the flow?

I am appreciate it if someone can share some ideas.

Best regards dong

[stoiver]

@Dongxueyang I know some people have written wrappers to create vtk output. They may be able to help. I would also suggest that you submit your questions to the anuga mailing list https://lists.sourceforge.net/lists/listinfo/anuga-user which has a larger readership than the github issues.

[Dongxueyang]

@stoiver Thanks a lot. I have post my question to the mailing list. And I have another question. How to simulate the process that water flows into the domain. Which boundary condition should I adapt? Has the example "cairns_excel" achieved the function of inlet/outlet? Best wishes dong

[stoiver]

@Dongxueyang, I tend to use an inlet operator instead of a boundary condition to control the input of water into the domain. Notebook 2 (Merewether) in the anuga clinic notebooks https://github.com/stoiver/anuga-clinic-2018 shows how to use an inlet operator.

[Dongxueyang]

@stoiver Thanks for your reply. And I have learned to complete inlet with your help successfully. Best regards dong

[Dongxueyang]

@stoiver Hi, I have another question. Now I am trying to output the section flow . And I use file "cairns_excel/flow_through_cross_sections.py" And I got some flux of the section. The value of flux is not equal the real flow. So what is the meaning of flux. And how can I calculate real flow with flux? Best regards dong

[Dongxueyang]

@stoiver Hi , what is the meaning of xmomentum and ymomentum. I see the unit of them is m^2/sec . how to get the momentum with flow velicity? dong

[stoiver]

@Dongxueyang, xmomentum is the depth integrated x velocity, and so can also be considered as the depth averaged velocity times the depth.

[Dongxueyang]

@stoiver Thanks for your explanition. I made a model of a channel with an inlet at left. The bottom of channel is 40m long and 5m wide, and the elevation of the bottom is -1m. I set the inlet boundary is Bi=anuga.Dirichlet_boundary(1.0, 1.0, 0.0). I output the flow,depth,velocity of section[(5 , 0),( 0, 0)].(the section of inlet) . I got depth=1.4m, velocity=1.93m/s, flow=1.36m^3/s. (Calculate from inlet_boundary: stage=1.0m, depth=stage-elevation=1-(-1)=2.0m. xmomentum=1.0
if flow is 1.36m^3/s, velocity(averaged depth)=1.36/(2.0*5.0)=0.136m/s. Is it right?) But how can I control the inlet flow by momentum,if I want to give a const flow Q. I do not know which value the momentum will be. Hope for your reply dong