def __init__(self,
boundary,
btype,
value,
porosity=None,
description=None):
"""Boundary condition initialisation with :
- one or several boundary(ies) (object Region). The objet is a region in order to be able to treat several materials
- a boundary condition type.
It can be:
a displacement
a normal force ( or a pressure)
In the case of a pressure, the surface of the boundary is retrieved to derive the normal force
- A boundary condition value. Value depends of boundary condition type.
"""
print("b", boundary, boundary.__class__.__name__)
print("k", btype)
print("v", value)
boundaryConditionDico = makeDico(
Dirichlet=[Displacement, NormalForce, Pressure])
CommonBoundaryCondition.__init__(self, boundary, btype, value,
boundaryConditionDico, porosity,
description)
return None
def __init__(self,
boundary,
btype,
value,
porosity=None,
description=None):
"""Boundary condition initialisation with :
- one or several boundaries (object Body or StructuredMesh
or list of Supports)
- a boundary condition type.
It can be Dirichlet,Neumann and in the future Flux
- a boundary condition value. Value depends of boundary condition type.
If Dirichlet, value can be Head, Pressure or SaturationLevel.
If Neumann, value can be HeadGradient or PressureGradient.
If Flux, value can be HydraulicFlux.
"""
#print("b",boundary)
#print("k",btype)
#print("v",value)
boundaryConditionDico = makeDico(Dirichlet=[Head, Pressure, SaturationLevel],\
Neumann=[HeadGradient, PressureGradient],\
Flux=[HydraulicFlux])
CommonBoundaryCondition.__init__(self, boundary, btype, value,
boundaryConditionDico, porosity,
description)
return None
def __init__(self, boundary, type,value, description = None):
"""Boundary condition initialisation with :
- one or several boundaries
- a boundary condition type.
It can be Dirichlet,Neumann,Mixed or Flux
- a boundary condition value. Value depend of boundary condition type.
If Dirichlet, value can be Pressure.
If Neumann, value can be PressureGradient.
If Flux, value can be HydraulicFlux.
If Mixed, value can be PressureMixed.
"""
from datamodel import Pressure, PressureGradient, PressureMixed
from datamodel import HydraulicFlux
bcdict=makeDico(Dirichlet=[Pressure],Neumann=[PressureGradient],\
Mixed=[PressureMixed],Flux=[HydraulicFlux])
CommonBoundaryCondition.__init__(self, boundary, type, value, bcdict, description)
def __init__(self, boundary, type, value, description = None):
"""Boundary condition initialisation with :
- one or several boundaries (object Support or StructuredMesh
or list of Supports)
- a boundary condition type.
It can be Dirichlet,Neumann or Flux
- a boundary condition value. Value depend of boundary condition type.
If Dirichlet, value can be Head, Pressure or SaturationDegree.
If Neumann, value can be HeadGradient or PressureGradient.
If Flux, value can be HydraulicFlux.
"""
bcdict = makeDico (Dirichlet=[Head,Pressure,SaturationDegree],\
Neumann=[HeadGradient,PressureGradient],\
Flux=[HydraulicFlux])
CommonBoundaryCondition.__init__(self,boundary, type, value, bcdict, description)
return
def __init__(self, boundary, type, value, description=None):
"""Boundary condition initialisation with :
- one or several boundaries (object Support or StructuredMesh
or list of Supports)
- a boundary condition type.
It can be Dirichlet,Neumann or Flux
- a boundary condition value. Value depend of boundary condition type.
If Dirichlet, value can be Head, Pressure or SaturationDegree.
If Neumann, value can be HeadGradient or PressureGradient.
If Flux, value can be HydraulicFlux.
"""
bcdict = makeDico (Dirichlet=[Head,Pressure,SaturationDegree],\
Neumann=[HeadGradient,PressureGradient],\
Flux=[HydraulicFlux])
CommonBoundaryCondition.__init__(self, boundary, type, value, bcdict,
description)
return
def __init__(self, boundary, type, value, description=None):
"""Boundary condition initialisation with :
- one or several boundaries
- a boundary condition type.
It can be Dirichlet, Neumann, Mixed or Flux.
- a boundary condition value. Value depend of boundary condition type.
If Dirichlet, value can be Pressure.
If Neumann, value can be PressureGradient.
If Flux, value can be HydraulicFlux.
If Mixed, value can be PressureMixed.
"""
from datamodel import Pressure, PressureGradient, PressureMixed
from datamodel import HydraulicFlux
bcdict = makeDico(Dirichlet=[Head, Pressure, Temperature], Neumann=[PressureGradient],\
Mixed=[PressureMixed], Flux=[HydraulicFlux])
CommonBoundaryCondition.__init__(self, boundary, type, value, bcdict,
description)
def __init__(self, boundary, btype, value, description=None):
"""
Boundary condition initialisation with :
- one or several boundaries (a body or a structured mesh element)
- a boundary condition type.
It can be Dirichlet,Neumann or Flux
- a boundary condition value. Value depends of boundary condition type.
If Dirichlet, value can be Head or Pressure.
If Neumann, value can be HeadGradient or PressureGradient.
If Flux, value can be HydraulicFlux.
"""
boundaryConditionDico = makeDict(
Dirichlet=[Head, Pressure, Temperature],
Neumann=[HeadGradient, PressureGradient],
Flux=[HydraulicFlux])
CommonBoundaryCondition.__init__(self,boundary, btype, value,\
boundaryConditionDico, porosity = None, description = None)
def __init__(self, boundary, btype, value, porosity = None, description = None):
"""Boundary condition initialisation with :
- one or several boundaries (object Body or StructuredMesh
or list of Supports)
- a boundary condition type.
It can be Dirichlet,Neumann and in the future Flux
- a boundary condition value. Value depends of boundary condition type.
If Dirichlet, value can be Head, Pressure or SaturationLevel.
If Neumann, value can be HeadGradient or PressureGradient.
If Flux, value can be HydraulicFlux.
"""
print("b",boundary)
print("k",btype)
print("v",value)
boundaryConditionDico = makeDico(Dirichlet=[Head,Pressure,SaturationLevel],\
Neumann=[HeadGradient,PressureGradient],\
Flux=[HydraulicFlux])
CommonBoundaryCondition.__init__(self,boundary, btype, value, boundaryConditionDico, porosity, description)
return None
def __init__(self, boundary, btype, value, porosity = None, description = None):
"""Boundary condition initialisation with :
- one or several boundary(ies) (object (Body) Region). The objet is a region in order to be able to treat several materials
- a boundary condition type.
It can be:
a displacement
a normal force ( or a pressure)
In the case of a pressure, the surface of the boundary is retrieved to derive the normal force
- A boundary condition value. Value depends of boundary condition type.
"""
print("b",boundary,boundary.__class__.__name__)
print("k",btype)
print("v",value)
boundaryConditionDico = makeDico(Dirichlet=[Displacement, NormalForce, Pressure, ChemicalState])
#print boundaryConditionDico.keys()
#raw_input()
CommonBoundaryCondition.__init__(self, boundary, btype, value, boundaryConditionDico, description)
BoundaryCondition.porosity = porosity
ok = 0
for val in value:
if val.__class__.__name__ == "ChemicalState":
ok = 1
self.chemicalState = val
break
pass
if ok == 0:
raise Warning(" the association of a chemical state to a boundary is mandatory ")
ok = 0
for val in value:
if val.__class__.__name__ in ["Displacement","NormalForce"]:
ok = 1
self.value = val
break
pass
if ok == 0:
raise Warning(" check the definition of boundary conditions")
return None
def __init__(self, boundary, btype, value = None, massTCoef = None, velocity = None, flowRate = None, porosity = None, timeVariation = None,
description = None):
"""
Constructor with :
- boundary : a mesh part element of type Cartesian or Unstructured ( made of bodies)
- btype : is a string and should be "Dirichlet", "Flux", "Mixed", "Neumann"
For a "symmetry", a Neumann boundary condition with g = 0 must be specified
- OPTIONAL :
--> value : a PhysicalQuantity or a list of tuples (PhysicalQuantity,species)
or a ChemicalState
--> massTCoef : float : mass transfer coefficient or set to zero
--> velocity : object Velocity
--> porosity : a scalar.
--> flowRate : a Flowrate, see PhysicalQuantities
--> timeVariation : a list of tuples [(time,chemical state)] , [(time,(list of species and eventually temperature))];
the temperature can also be introduced through a file.
-- description a string which will be eventually set as a support for the model comprehension
"""
bcDico = makeDico(Dirichlet = [ChemicalState, Head, Displacement, NormalForce],\
Flux = [ChemicalState, HeadGradient],\
Neumann = [ChemicalState, HeadGradient])
CommonBoundaryCondition.__init__(self,boundary, btype, value, bcDico, description)
# print "dbg commonmodel CommonBoundaryCondition1"
if type(boundary) is types.ListType:
# print "type of boundary is list type "
#raw_input("type of boundary is list type ")
verifyClassList(boundary,[ CartesianMesh, Body])
pass
else:
memberShip(boundary,[ CartesianMesh, Body])
#raw_input("membership ")
pass
#raw_input("dbg commonmodel CommonBoundaryCondition2")
self.boundary = boundary
if type(btype) != types.StringType:
raise TypeError, " problem on the definition of the boundary type "
if btype.lower() not in ["dirichlet","symmetry","flux","mixed","neumann","noflux"]: raise Exception, " check the boundary condition kind"
self.btype = btype
self.chemicalStateValue = None
self.headValue = None
self.massTCoef = 0.
self.value_species = None
self.value_property = None
self.value = None
#
# the next ones are linked to a well sim.
#
self.enthalpyBoundaryCondition = None
self.wellMassFlowBoundaryCondition = None
self.wellPressureBoundaryCondition = None
#
# We treat B.C.
# by default, a chemical state is introduced
# and in the case of a transient flow, eventually a list
# made of a chemical state, a displacement, a head.
#
if type(value) is types.ListType:
#
# useful for debugging
#
#for i in value:
# print "dbg commonmodel",type(i)
# pass
verifyClassList(value, [ Head, ChemicalState, Displacement, NormalForce, TupleType])
for bc in value:
if isinstance(bc, Head):
self.headValue = bc # it should be the charge
pass
elif isinstance(bc, NormalForce):
self.normalForceValue = bc # it should be NormalForce
pass
elif isinstance(bc, Displacement):
self.displacementValue = bc # it should be Displacement
pass
elif isinstance(bc, ChemicalState):
self.value = bc
self.chemicalStateValue = bc # it should be ChemicalState
pass
elif bc[0].lower() == "enthalpy": # it can also be an enthalpy in the
# case of a well
#
if type(bc[1]) == types.StringType:
self.enthalpyBoundaryCondition = refindall(recompile(r'([xyzXYZ0-9.*/+-])'),bc[1])
pass
elif type(bc[1]) in [types.FloatType,types.IntType]:
self.enthalpyBoundaryCondition = bc[1]
pass
elif bc[0].lower() == "wellpressure": # it can also be the pressure in the
# case of a well
#
if type(bc[1]) == types.StringType:
self.wellPressureBoundaryCondition = refindall(recompile(r'([xyzXYZ0-9.*/+-])'),bc[1])
pass
elif type(bc[1]) in [types.FloatType,types.IntType]:
self.wellPressureBoundaryCondition = bc[1]
#print("commonmodel well pressure debug yes\n")
#raw_input()
pass
pass
elif bc[0].lower() == "wellmassflow": # it can also be the mass flow in the
# case of a well
#
if type(bc[1]) == types.StringType:
self.wellMassFlowBoundaryCondition = refindall(recompile(r'([$mfunction()ifelse{} ><_=xyzXYZ0-9.*/+-])'),bc[1])
pass
elif type(bc[1]) in [types.FloatType,types.IntType]:
self.wellMassFlowBoundaryCondition = bc[1]
pass
elif type(bc[1]) is tuple:
self.wellMassFlowBoundaryCondition = bc[1]
pass
pass
else:
#self.value = bc # it should be chemistry
pass
pass
pass
else:
memberShip(value,[PhysicalQuantity, ChemicalState, Displacement, NormalForce])
if (isinstance(value, PhysicalQuantity) or
type(value) is types.ListType):
self.value_species, self.value_property = createList(value, PhysicalQuantity)
pass
else:
self.value = value
self.chemicalStateValue = value
pass
pass
print "massTCoef",massTCoef,type(massTCoef)
if massTCoef:
memberShip(massTCoef,[types.FloatType])
if (type(massTCoef) is types.FloatType):
self.massTCoef = massTCoef
pass
else:
self.massTCoef = 0.0
pass
print " common model mass transfer coefficient ",self.massTCoef
pass
if porosity:
self.porosity = porosity
pass
if velocity:
memberShip(velocity,Velocity)
pass
self.velocity = velocity
if flowRate:
if flowRate.__class__.__name__=="FlowRate":
pass
else:
flowrate = FlowRate(flowrate,"m**3/s") # the flow rate is supposed to be in m**3/s
pass
self.flowRate = flowRate
if timeVariation:
if type(timeVariation) != types.ListType:
raise typeError, " Time variation should be a list"
for item in timeVariation:
if type(item[0]) not in [types.FloatType,types.IntType]:
raise typeError, "item[@] should be a list"
memberShip(item[1],[ChemicalState])
pass
pass
self.timeVariation = timeVariation
return None
def __init__(self, boundary, btype, value = None, massTCoef = None, velocity = None, flowRate = None, porosity = None, timeVariation = None,
description = None):
"""
Constructor with :
- boundary : a mesh part element of type Cartesian or Unstructured ( made of bodies)
- btype : is a string and should be "Dirichlet", "Flux", "Mixed", "Neumann"
For a "symmetry", a Neumann boundary condition with g = 0 must be specified
- OPTIONAL :
--> value : a PhysicalQuantity or a list of tuples (PhysicalQuantity,species)
or a ChemicalState
--> massTCoef : float : mass transfer coefficient or set to zero
--> velocity : object Velocity
--> porosity : a scalar.
--> flowRate : a Flowrate, see PhysicalQuantities
--> timeVariation : a list of tuples [(time,chemical state)] , [(time,(list of species and eventually temperature))];
the temperature can also be introduced through a file.
-- description a string which will be eventually set as a support for the model comprehension
"""
bcDico = makeDico(Dirichlet = [ChemicalState, Head, Displacement, NormalForce],\
Flux = [ChemicalState, HeadGradient],\
Neumann = [ChemicalState, HeadGradient])
CommonBoundaryCondition.__init__(self,boundary, btype, value, bcDico, description)
# print "dbg commonmodel CommonBoundaryCondition1"
if type(boundary) is types.ListType:
# print "type of boundary is list type "
#raw_input("type of boundary is list type ")
verifyClassList(boundary,[ CartesianMesh, Body])
pass
else:
memberShip(boundary,[ CartesianMesh, Body])
#raw_input("membership ")
pass
#raw_input("dbg commonmodel CommonBoundaryCondition2")
self.boundary = boundary
if type(btype) != types.StringType:
raise TypeError, " problem on the definition of the boundary type "
if btype.lower() not in ["dirichlet","symmetry","flux","mixed","neumann","noflux"]: raise Exception, " check the boundary condition kind"
self.btype = btype
self.chemicalStateValue = None
self.headValue = None
self.massTCoef = 0.
self.value_species = None
self.value_property = None
self.value = None
#
# the next ones are linked to a well sim.
#
self.enthalpyBoundaryCondition = None
self.wellMassFlowBoundaryCondition = None
self.wellPressureBoundaryCondition = None
#
# We treat B.C.
# by default, a chemical state is introduced
# and in the case of a transient flow, eventually a list
# made of a chemical state, a displacement, a head.
#
if type(value) is types.ListType:
#
# useful for debugging
#
#for i in value:
# print "dbg commonmodel",type(i)
# pass
verifyClassList(value, [ Head, ChemicalState, Displacement, NormalForce, TupleType])
for bc in value:
if isinstance(bc, Head):
self.headValue = bc # it should be the charge
pass
elif isinstance(bc, NormalForce):
self.normalForceValue = bc # it should be NormalForce
pass
elif isinstance(bc, Displacement):
self.displacementValue = bc # it should be Displacement
pass
elif isinstance(bc, ChemicalState):
self.value = bc
self.chemicalStateValue = bc # it should be ChemicalState
pass
elif bc[0].lower() =="enthalpy": # it can also be an enthalpy in the
# case of a well
#
if type(bc[1]) == types.StringType:
self.enthalpyBoundaryCondition = refindall(recompile(r'([xyzXYZ0-9.*/+-])'),bc[1])
pass
elif type(bc[1]) in [types.FloatType,types.IntType]:
self.enthalpyBoundaryCondition = bc[1]
pass
elif bc[0].lower() =="wellpressure": # it can also be the pressure in the
# case of a well
#
if type(bc[1]) == types.StringType:
self.wellPressureBoundaryCondition = refindall(recompile(r'([xyzXYZ0-9.*/+-])'),bc[1])
pass
elif type(bc[1]) in [types.FloatType,types.IntType]:
self.wellPressureBoundaryCondition = bc[1]
#print("commonmodel well pressure debug yes\n")
#raw_input()
pass
pass
elif bc[0].lower() =="wellmassflow": # it can also be the mass flow in the
# case of a well
#
if type(bc[1]) == types.StringType:
self.wellMassFlowBoundaryCondition = refindall(recompile(r'([xyzXYZ0-9.*/+-])'),bc[1])
pass
elif type(bc[1]) in [types.FloatType,types.IntType]:
self.wellMassFlowBoundaryCondition = bc[1]
pass
pass
else:
#self.value = bc # it should be chemistry
pass
pass
pass
else:
memberShip(value,[PhysicalQuantity, ChemicalState, Displacement, NormalForce])
if (isinstance(value, PhysicalQuantity) or
type(value) is types.ListType):
self.value_species, self.value_property = createList(value, PhysicalQuantity)
pass
else:
self.value = value
self.chemicalStateValue = value
pass
pass
print "massTCoef",massTCoef,type(massTCoef)
if massTCoef:
memberShip(massTCoef,[types.FloatType])
if (type(massTCoef) is types.FloatType):
self.massTCoef = massTCoef
pass
else:
self.massTCoef = 0.0
pass
print " common model mass transfer coefficient ",self.massTCoef
pass
if porosity:
self.porosity = porosity
pass
if velocity:
memberShip(velocity,Velocity)
pass
self.velocity = velocity
if flowRate:
if flowRate.__class__.__name__=="FlowRate":
pass
else:
flowrate = FlowRate(flowrate,"m**3/s") # the flow rate is supposed to be in m**3/s
pass
self.flowRate = flowRate
if timeVariation:
if type(timeVariation) != types.ListType:
raise typeError, " Time variation should be a list"
for item in timeVariation:
if type(item[0]) not in [types.FloatType,types.IntType]:
raise typeError, "item[@] should be a list"
memberShip(item[1],[ChemicalState])
pass
pass
self.timeVariation = timeVariation
return None