yade-users team mailing list archive

Thread
Date

Re: [Question #702316]: how to get the change of permeability during triaxial compression in flowengine

To: yade-users@xxxxxxxxxxxxxxxxxxx
From: Robert Caulk <question702316@xxxxxxxxxxxxxxxxxxxxx>
Date: Fri, 01 Jul 2022 11:41:08 -0000
Reply-to: question702316@xxxxxxxxxxxxxxxxxxxxx
Sender: noreply@xxxxxxxxxxxxx

Question #702316 on Yade changed:
https://answers.launchpad.net/yade/+question/702316

    Status: Open => Answered

Robert Caulk proposed the following answer:
Calling a function like this in your addData:

def getPermeability():
	
	# set new bcs for perm check
	flow.bndCondIsPressure = [0, 0, 0, 0, 1, 1]
	flow.bndCondValue = [0, 0, 0, 0, deltaP, 0]

	# force update the triangulation	
	flow.updateTriangulation = True

	flow.emulateAction()
	flow.emulateAction()

	Qin = flow.getBoundaryFlux(4)
	Qout = flow.getBoundaryFlux(5)

        perm = abs(Qin) * flow.viscosity * L / (A * deltaP)

        print("Q", Qin, 'perm', perm, "A", A, "L",L)

        # reset boundary conditions before returning to sim
	flow.bndCondIsPressure = [0, 0, 0, 1, 0, 0]
	flow.bndCondValue = [0, 0, 0, 0, 0, 0]
	flow.updateTriangulation = True

        return perm

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