search

FLO-2DSoftware / FLO-2DMapCrafter

The QGIS FLO-2D MapCrafter plugin repository
0 stars 0 forks source link

BOR Hazard Maps for two phase using wrong data #66

Open FLO-2DKaren opened 7 months ago

FLO-2DKaren commented 7 months ago

Hi Robson,

You need to change the data you read for MapCrafter Hazard when looking at Two-Phase. I think you should use the VxD or the Combined Max V * the Combined Max D.

I think they are only using the water map.

FLO-2DKaren commented 7 months ago

https://github.com/FLO-2DSoftware/FLO-2DMapCrafter/assets/20424460/c786e352-06a5-4a82-a89b-53c3656b2af6

rpachaly commented 1 month ago

Yes, we need to address these issues but the problem is not MapCrafter.

USBR and SWISS methods:

We have the DEPTHMAX_2PHASE_COMBINED.OUT but we don't have a maximum velocity of both mudflow and water. We only have max velocity of water and max velocity of mudflow. That's why I did this considering only water.

FLO-2DNoemi commented 1 month ago

Combining the velocities is not straightforward. We can allow users to select the velocities they want to use. A conservative approach would be to use the mudflow velocities for this calculation, while a less conservative approach would involve using the water velocities. We should ask Jim for his input.

On Wed, Jul 31, 2024 at 1:08 PM Robson Pachaly @.***> wrote:

Yes, we need to address these issues but the problem is not MapCrafter.

USBR and SWISS methods:

We have the DEPTHMAX_2PHASE_COMBINED.OUT but we don't have a maximum velocity of both mudflow and water. We only have max velocity of water and max velocity of mudflow. That's why I did this considering only water.

— Reply to this email directly, view it on GitHub https://github.com/FLO-2DSoftware/FLO-2DMapCrafter/issues/66#issuecomment-2260980725, or unsubscribe https://github.com/notifications/unsubscribe-auth/AE32O3QSKHMJGMF6XJQYW3LZPEKYBAVCNFSM6AAAAABLY4T33KVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDENRQHE4DANZSGU . You are receiving this because you are subscribed to this thread.Message ID: @.***>

FLO-2DNoemi commented 1 month ago

Maybe a combination of the velocities at each cell. I think we could see grids with water flow depths but mudflow velocity equal to zero, so if the user selects the mudflow velocities for this calculation, the velocity will be zero. I am not sure if I explain myself here.

On Wed, Jul 31, 2024 at 1:26 PM Noemi Gonzalez @.***> wrote:

Combining the velocities is not straightforward. We can allow users to select the velocities they want to use. A conservative approach would be to use the mudflow velocities for this calculation, while a less conservative approach would involve using the water velocities. We should ask Jim for his input.

On Wed, Jul 31, 2024 at 1:08 PM Robson Pachaly @.***> wrote:

Yes, we need to address these issues but the problem is not MapCrafter.

USBR and SWISS methods:

We have the DEPTHMAX_2PHASE_COMBINED.OUT but we don't have a maximum velocity of both mudflow and water. We only have max velocity of water and max velocity of mudflow. That's why I did this considering only water.

— Reply to this email directly, view it on GitHub https://github.com/FLO-2DSoftware/FLO-2DMapCrafter/issues/66#issuecomment-2260980725, or unsubscribe https://github.com/notifications/unsubscribe-auth/AE32O3QSKHMJGMF6XJQYW3LZPEKYBAVCNFSM6AAAAABLY4T33KVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDENRQHE4DANZSGU . You are receiving this because you are subscribed to this thread.Message ID: @.***>

rpachaly commented 1 month ago

In my opinion, it does not make sense to add the velocities or to make an average of it. I'm not sure if taking the max of the velocity or the mudflow for that specific grid element is also a good approach.

MAYBE a weighted average makes sense, considering the depth and velocity of mudflow and water.

FLO-2DNoemi commented 1 month ago

I agree, it is not straightforward to combine these velocities, using the water velocities for the product of H x V is the worst scenario.

On Wed, Jul 31, 2024 at 1:34 PM Robson Pachaly @.***> wrote:

In my opinion, it does not make sense to add the velocities or to make an average of it. I'm not sure if taking the max of the velocity or the mudflow for that specific grid element is also a good approach.

MAYBE a weighted average makes sense, considering the depth and velocity of mudflow and water.

— Reply to this email directly, view it on GitHub https://github.com/FLO-2DSoftware/FLO-2DMapCrafter/issues/66#issuecomment-2261022531, or unsubscribe https://github.com/notifications/unsubscribe-auth/AE32O3WBHNIPLJ65IYSPZYTZPEN3RAVCNFSM6AAAAABLY4T33KVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDENRRGAZDENJTGE . You are receiving this because you commented.Message ID: @.***>

FLO-2DKaren commented 1 month ago

What does the Vel_X_Depth.out show use? Maybe we can extrapolate the velocity from that file?

FLO-2DNoemi commented 1 month ago

It reports the product of depth X Vel according to the following:

2 PHASE FLOW MAXIMUM VEL X DEPTH FOR FLUID AND MUDFLOW

VEL_X_DEPTH=0. !.....MAXIMUM VEL_X_DEPTH...COMPARE WITH FLOODPLAIN VALUE

    IF(MUD.LE.1)THEN

IF(VEL_X_DEPTH(J).LT.ABS(FCFD(I)VELOC(I)))VEL_X_DEPTH(J)=ABS(FCFD(I)VELOC(I)) ENDIF

    IF(MUD.EQ.2)THEN

IF(VEL_X_DEPTH(J).LT.(FCFD(I)VELOC(I)))VEL_X_DEPTH(J)=(FCFD(I)VELOC(I))

IF(VEL_X_DEPTH_MUD(J).LT.ABS(FCFD_MUD(I)VELOC_MUD(I)))VEL_X_DEPTH_MUD(J)=ABS(FCFD_MUD(I)VELOC_MUD(I))

IF(VEL_X_DEPTH_MUD(J).GT.VEL_X_DEPTH(J))VEL_X_DEPTH(J)=VEL_X_DEPTH_MUD(J) ENDIF

!......MAXIMUM VELOCITY X DEPTH

    No Mudflow
      IF(VEL(J,JDIR).LT.0.)THEN        Only when water velocity is less

than zero IF(VEL_X_DEPTH(J).LT.ABS(FPD(J)VEL(J,JDIR)))THEN VEL_X_DEPTH(J)=ABS(FPD(J)VEL(J,JDIR)) ENDIF ENDIF

  2 PHASE FLOW MAXIMUM VEL X DEPTH FOR FLUID AND MUDFLOW
     IF(VEL(J,JDIR).LT.0.)THEN       Only when water velocity is less

than zero

IF(VEL_X_DEPTH(J).LT.ABS(FPD(J)VEL(J,JDIR)))VEL_X_DEPTH(J)=ABS(FPD(J)VEL(J,JDIR)) ENDIF IF(VEL_MUD(J,JDIR).LT.0.)THEN Only when mud velocity is less than zero

IF(VEL_X_DEPTH_MUD(J).LT.ABS(FPD_MUD(J)VEL_MUD(J,JDIR)))VEL_X_DEPTH_MUD(J)=ABS(FPD_MUD(J)VEL_MUD(J,JDIR))

IF(VEL_X_DEPTH_MUD(J).GT.VEL_X_DEPTH(J))VEL_X_DEPTH(J)=VEL_X_DEPTH_MUD(J)

     ENDIF

On Thu, Aug 1, 2024 at 9:39 AM Karen @.***> wrote:

What does the Vel_X_Depth.out show use? Maybe we can extrapolate the velocity from that file?

— Reply to this email directly, view it on GitHub https://github.com/FLO-2DSoftware/FLO-2DMapCrafter/issues/66#issuecomment-2263077804, or unsubscribe https://github.com/notifications/unsubscribe-auth/AE32O3QELMELUL364G3B7DLZPI3BFAVCNFSM6AAAAABLY4T33KVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDENRTGA3TOOBQGQ . You are receiving this because you commented.Message ID: @.***>