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)