fracture.density

Syntax

fdenc = fracture.density(<vlower><,vupper><,dfn>)

Get the fracture mass density. The fracture mass density is defined as the {P21 value in 2D; P32 value in 3D} as measured in the specified region. {P21 in 2D; P32 in 3D} is the cumulated fracture length (surface) in the specified region per unit area (volume). If dfn is not given the calculation is undertaken over all fractures of all DFNs.

Returns:

fdenc - fracture mass density value

Arguments:

vlower - optional lower bound of the calculation region

vupper - optional upper bound of the calculation region

dfn - optional DFN pointer

Usage Example

The following example illustrates how fracture.density can be used get the mass density of various DFNs, or of all fractures.

model new
model random 10000
model domain extent -5 5
;Generate fractures using the default template
fracture generate dfn 'p1' mass-density 0.1
fracture generate dfn 'p2' mass-density 0.2
fracture generate dfn 'p3' mass-density 0.3
fracture generate dfn 'p4' mass-density 0.4

; fish access to densities
fish define access_density
   local dfn1 = dfn.find('p1')
   local dfn2 = dfn.find('p2')
   local dfn3 = dfn.find('p3')
   local dfn4 = dfn.find('p4')
    lower = vector(-5,-5,-5)
    upper = vector(5,5,5)
   io.out('mass density of dfn p1    = ' + ... 
           string(fracture.density(lower,upper,dfn1)))
   io.out('mass density of dfn p2    = ' + ... 
           string(fracture.density(lower,upper,dfn2)))
   io.out('mass density of dfn p3    = ' + ... 
           string(fracture.density(lower,upper,dfn3)))
   io.out('mass density of dfn p4    = ' + ... 
           string(fracture.density(lower,upper,dfn4)))
   io.out('mass density of all dfns  = ' + ... 
           string(fracture.density(lower,upper)))
end
[access_density]

program return