def __init__(self, problem, solver_params_fluid, solver_params_flow0d): self.pb = problem self.solver_params_fluid = solver_params_fluid self.solver_params_flow0d = solver_params_flow0d # initialize nonlinear solver class self.solnln = solver_nonlin.solver_nonlinear_constraint_monolithic( self.pb, self.pb.pbs.V_v, self.pb.pbs.V_p, self.solver_params_fluid, self.solver_params_flow0d)
def __init__(self, problem, solver_params_solid, solver_params_flow0d): self.pb = problem self.solver_params_solid = solver_params_solid self.solver_params_flow0d = solver_params_flow0d # initialize nonlinear solver class self.solnln = solver_nonlin.solver_nonlinear_constraint_monolithic( self.pb, self.pb.pbs.V_u, self.pb.pbs.V_p, self.solver_params_solid, self.solver_params_flow0d) if self.pb.pbs.prestress_initial: # add coupling work to prestress weak form self.pb.pbs.weakform_prestress_u -= self.pb.work_coupling_prestr # initialize solid mechanics solver self.solverprestr = SolidmechanicsSolver(self.pb.pbs, self.solver_params_solid)