def getPfem():
global w
if w: return w
w = pfem
mshFile = runPath+os.sep+'birdImpact_deformable_panel_Mtf_Pfem_Smojver.msh'
print 'mshFile: ', mshFile
rho0 = 938.
mu = 0.001
U0 = 145.7
pbl = w.Problem()
pbl.rho0 = rho0
pbl.mu = mu
pbl.nonLinAlgorithm = 1
pbl.solScheme = 1
pbl.alpha = 1.2
pbl.beta = 0.
pbl.gravity = 0.
pbl.scalingU = U0
msh = w.MshData(pbl)
msh.load(mshFile)
print msh
scheme = w.BackwardEuler(msh, pbl)
w.Medium(msh, 13, mu, rho0, 3)
w.Medium(msh, 16, mu, rho0, 1)
# boundaries
w.Boundary(msh, 15, 3, 0.0)
w.Boundary(msh, 13, 1, 0.0)
w.Boundary(msh, 13, 2, 0.0)
#Initial velocity
bird = w.Group(msh, 16)
loadingset = w.LoadingSet(msh)
loadingset.add(1,w.InitialVelocity(msh,bird,0.,-U0,0.))
scheme.savefreq=1
scheme.nthreads=3
scheme.gamma = 0.5
scheme.omega = 0.5
scheme.addRemoveNodesOption = True
scheme.tollNLalgo = 1e-7
#Results
extManager = w.ExtractorsManager(msh)
extManager.add(1,w.PositionExtractor(msh,163))
extManager.add(6,wt.KineticEnergyExtractor(msh,pbl,16))
extManager.add(7,wt.ViscousEnergyExtractor(msh,pbl,scheme,16))
extManager.add(8,wt.PressureWorkExtractor(msh,pbl,scheme,16))
extManager.add(9,w.MassExtractor(msh,pbl,16))
gui = v.MeshViewer(msh, scheme, True)
return Module(w, msh, pbl, bird, loadingset, scheme, extManager, gui)
def getPfem():
global w
if w: return w
w = pfem
mshFile = runPath + os.sep + 'waterColoumnFallFakeSolidSolver.msh'
rho0 = 1000.
mu = 0.001
pbl = w.Problem()
pbl.rho0 = rho0
pbl.mu = mu
pbl.nonLinAlgorithm = 1
pbl.solScheme = 1
pbl.alpha = 1.2
pbl.extP = 100.
pbl.scalingU = 1.0
msh = w.MshData(pbl)
msh.load(mshFile)
print msh
scheme = w.BackwardEuler(msh, pbl)
w.Medium(msh, "Walls", mu, rho0, 3)
w.Medium(msh, "Face1", mu, rho0, 1)
# boundaries
w.Boundary(msh, "Free", 3, pbl.extP)
w.Boundary(msh, "Walls", 1, 0.0)
w.Boundary(msh, "Walls", 2, 0.0)
scheme.ttot = 1.0
scheme.savefreq = 1
scheme.nthreads = 3
scheme.gamma = 0.5
scheme.omega = 0.5
scheme.addRemoveNodesOption = True
scheme.tollNLalgo = 1e-12
#Results
extManager = w.ExtractorsManager(msh)
extManager.add(1, w.PositionExtractor(msh, 1))
extManager.add(2, w.IntForceExtractor(msh, 1))
extManager.add(3, w.ExtForceExtractor(msh, 1))
extManager.add(4, w.IneForceExtractor(msh, 1))
extManager.add(5, w.IntForceExtractor(msh, 6))
extManager.add(6, w.ExtForceExtractor(msh, 6))
extManager.add(7, w.IneForceExtractor(msh, 6))
extManager.add(8, w.PressureExtractor(msh, 1))
extManager.add(9, w.VelocityExtractor(msh, "Face1"))
extManager.add(10, w.MassExtractor(msh, pbl, "Face1"))
extManager.add(11, wt.KineticEnergyExtractor(msh, pbl, "Face1"))
extManager.add(12, wt.ViscousEnergyExtractor(msh, pbl, scheme, "Face1"))
gui = v.MeshViewer(msh, scheme, True)
return Module(w, msh, pbl, scheme, extManager, gui)
def getPfem():
global w
if w: return w
w = pfem
mshFile = runPath + os.sep + 'birdImpact_deformable_panel_Mtf_Pfem_Axisym.msh'
rho0 = 1000.
mu = 0.001
U0 = 100
pbl = w.Problem()
pbl.rho0 = rho0
pbl.mu = mu
pbl.nonLinAlgorithm = 1
pbl.solScheme = 1
pbl.alpha = 1.2
pbl.beta = 0.
pbl.gravity = 0.
pbl.scalingU = U0
pbl.isAxisymmetric = True
msh = w.MshData(pbl)
msh.load(mshFile)
print msh
scheme = w.BackwardEuler(msh, pbl)
w.Medium(msh, 13, mu, rho0, 3)
w.Medium(msh, 17, mu, rho0, 1)
# boundaries
w.Boundary(msh, 15, 7, 0.0)
w.Boundary(msh, 16, 3, 0.0)
w.Boundary(msh, 13, 1, 0.0)
w.Boundary(msh, 13, 2, 0.0)
# --- Necessary for the correct elimination of nodes along the axis of symmetry ---
n4 = msh.getNode(4)
n4.isOnAxisOfSymmetry = True
# ---
#Initial velocity
bird = w.Group(msh, 17)
loadingset = w.LoadingSet(msh)
loadingset.add(1, w.InitialVelocity(msh, bird, 0., -U0, 0.))
scheme.savefreq = 1
scheme.nthreads = 3
scheme.gamma = 0.5
scheme.omega = 0.5
scheme.addRemoveNodesOption = True
scheme.tollNLalgo = 1e-7
#Results
extManager = w.ExtractorsManager(msh)
extManager.add(1, w.PositionExtractor(msh, 4))
extManager.add(6, wt.KineticEnergyExtractor(msh, pbl, 17))
extManager.add(7, wt.ViscousEnergyExtractor(msh, pbl, scheme, 17))
extManager.add(8, wt.PressureWorkExtractor(msh, pbl, scheme, 17))
extManager.add(9, w.MassExtractor(msh, pbl, 17))
gui = v.MeshViewer(msh, scheme, True)
return Module(w, msh, pbl, bird, loadingset, scheme, extManager, gui)
def getPfem():
global w
if w: return w
w = pfem
mshFile = runPath + os.sep + 'birdImpact_deformable_panel_Mtf_Pfem_Smojver_Axisym.msh'
rho0 = 938.
mu = 0.001
U0 = 145.7
pbl = w.Problem()
pbl.rho0 = rho0
pbl.mu = mu
pbl.nonLinAlgorithm = 1
pbl.solScheme = 1
pbl.alpha = 1.2
pbl.beta = 0.
pbl.gravity = 0.
pbl.scalingU = U0
pbl.isAxisymmetric = True
msh = w.MshData(pbl)
msh.load(mshFile)
print msh
toll = 1e-6
nItMax = 10
solScheme = w.SchemeMonolithicPSPG(msh, pbl)
convCriterion = w.ForcesBalanceNormedBodyForceCriterion(msh, pbl, toll)
nonLinAlgo = w.PicardAlgorithm(solScheme, convCriterion, nItMax)
scheme = w.BackwardEuler(msh, pbl, nonLinAlgo)
w.Medium(msh, 17, mu, rho0, 3)
w.Medium(msh, 21, mu, rho0, 1)
# boundaries
w.Boundary(msh, 20, 7, 0.0)
w.Boundary(msh, 19, 3, 0.0)
w.Boundary(msh, 17, 1, 0.0)
w.Boundary(msh, 17, 2, 0.0)
# --- Necessary for the correct elimination of nodes along the axis of symmetry ---
n4 = msh.getNode(6)
n4.isOnAxisOfSymmetry = True
# ---
#Initial velocity
bird = w.Group(msh, 21)
loadingset = w.LoadingSet(msh)
loadingset.add(1, w.InitialVelocity(msh, bird, 0., -U0, 0.))
scheme.savefreq = 1
scheme.nthreads = 3
scheme.gamma = 0.5
scheme.omega = 0.5
scheme.addRemoveNodesOption = True
scheme.tollNLalgo = 1e-7
#Results
extManager = w.ExtractorsManager(msh)
extManager.add(1, w.PositionExtractor(msh, 6))
extManager.add(6, wt.KineticEnergyExtractor(msh, pbl, 21))
extManager.add(7, wt.ViscousEnergyExtractor(msh, pbl, scheme, 21))
extManager.add(8, wt.PressureWorkExtractor(msh, pbl, scheme, 21))
extManager.add(9, w.MassExtractor(msh, pbl, 21))
gui = v.MeshViewer(msh, scheme, True)
return Module(w, msh, pbl, solScheme, nonLinAlgo, convCriterion, scheme,
extManager, gui)
