def AddVariables(fluid_model_part, structure_model_part):
import incompressible_fluid_solver
incompressible_fluid_solver.AddVariables(fluid_model_part);
structure_model_part.AddNodalSolutionStepVariable(RELAXED_DISP);
print "variables for FractionalStepCoupling added correctly"
benchmarking.Output(min_press, "minimum pressure", 0.00001)
benchmarking.Output(max_press, "maximum pressure", 0.00001)
benchmarking.Output(id_min_vel,
"Id of the node with minimum velocity norm", 0.0)
benchmarking.Output(x_min_vel, "coord x minimum velocity norm", 0.0)
benchmarking.Output(y_min_vel, "coord y minimum velocity norm", 0.0)
# defining a model part
print("before creation of the model part")
model_part = ModelPart("FluidPart")
print("after creation of the model part")
# importing the solver files and adding the variables
import incompressible_fluid_solver
incompressible_fluid_solver.AddVariables(model_part)
# adding of Variables to Model Part should be here when the "very fix container will be ready"
# reading a model
gid_mode = GiDPostMode.GiD_PostBinary
use_multifile = MultiFileFlag.MultipleFiles
deformed_print_flag = WriteDeformedMeshFlag.WriteUndeformed
write_conditions = WriteConditionsFlag.WriteConditions
gid_io = GidIO("cavity2d", gid_mode, use_multifile, deformed_print_flag,
write_conditions)
write_conditions = WriteConditionsFlag.WriteElementsOnly
# gid_io.ReadMesh(model_part.GetMesh())
gid_io.ReadModelPart(model_part)
# gid_io.WriteMesh((model_part).GetMesh(),domain_size,0.0,GiDPostMode.GiD_PostBinary);
gid_io.InitializeMesh(0.0)
def AddVariables(fluid_model_part, structure_model_part):
import incompressible_fluid_solver
incompressible_fluid_solver.AddVariables(fluid_model_part)
print("variables for FractionalStepCoupling added correctly")