yade-users team mailing list archive

[mrhappy <question691802@xxxxxxxxxxxxxxxxxxxxx>]

New question #691802 on Yade:
https://answers.launchpad.net/yade/+question/691802

Hi,

I have a few questions on modeling hydraulic fractures using DFNFlowEngine. My model size is 40m by 40m by 40m and I have precrack (in shape of a circle) on the x-y plane located exactly at the center of the model. The center of the circle of the precrack is located at 20,20,20. Here are my questions:

1. I want to pump in 10m^3 of liquid in my precrack @ 1m^3 per second. To do this, this is what I have done:
ti = O.time
def Pumping(): 
  if O.time<(ti+10): 
    flow.imposeFlux(Vector3(20,20,20),-1) 
    print("PUMPING")
  else: 
    flow.imposeFlux(Vector3(20,20,20),0) 
    print("PUMPING STOPED")
O.engines=O.engines+[PyRunner(iterPeriod=1,command='Pumping()')] 

I am concerned about how I switch off the pumping. is it correct? Should I use clearImposedFlux?

2. I want to keep the fluid inside the precrack (or the newly developed crack). This means that there is no leakoff. To do this, I did this:
flow.permeabilityFactor=0
Is this correct? Will this keep the fluid only in the crack?

3. Going back to question 1, when I pump in fluid at 20,20,20 (the center of the precrack), I would expect that the crack would grow radially (i.e. begcome a bigger circle in the xy plane). However, the bonds (between the particles) just above and below (20,20,20) start to break first and crack mini cracks in the z direction before the precrack starts to grow radially. How can I prevent the crack from propagating in the z direction? I am aware that my imposed flux is a singular point and therefore, this is bound to happen. Should I change the way I load the flux?

4. What is flow.jointsResidualAperture ? According to the YADE manual, it is "residual aperture of joints". Is this the width of the precrack? Is so, how does this value effect the calculations?

5. What is flow.clampKValues? What is K? I read the manual, but I still don't understand what this does.

6. Should I use flow.blockHook for all the pores that are not within the precrack to increase the computations? 

Thank you so much for your help. I really appreciate everyone's time in responding to these questions

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.