Calculate seepage force

FLAC3D
1

Hello guys,

quick question, would you suggest me a method to calculate the seepage force vector at gridpoint level? The model is one-way coupled fluid->mech.

I only managed to estimate seepage force module at gridpoints (using gp.flow, as float) and force vector at zones (zone.flow, as vector).

Thanks,
Gian

1 Like
2

Is there already a solution for this query? I’m trying to solve the same “problem”
Thanks!
Josef

3

Currently it cannot be directly output. But again, we can use a small simple FISH function to do the job (Thanks, FISH!). We can loop all gridpoints, calculable the gridpoint flow by averaging flows of its neighboring zones with weights of zone volumes. We can then use the user defined vectors to visualize the gridpoint flow. Below is an example

model restore "ch2a.sav"

fish define calc_gp_flow
    loop foreach local g gp.list
        local gp_flow = vector(0,0,0)
        local z_vol = 0.0
        loop foreach local z gp.zone(g)
            zvol = zone.vol(z)
            gp_flow = gp_flow + zone.flow(z)* zvol
            z_vol = z_vol + zvol
        endloop
        gp_flow = gp_flow/z_vol
        zp = data.vector.create(gp.pos(g))
        data.vector.value(zp) = gp_flow
    endloop
end
[calc_gp_flow]

Then plot the user defined data → vector:

a
a3800×1474 164 KB

2 Likes