# NODE HISTORY RESULTS BEGIN scalar_vars = [] vector_vars = [DISPLACEMENT] for node in spheres_model_part.Nodes: node_to_follow_id = node.Id node.SetSolutionStepValue(VELOCITY_Z, terminal_velocity) results_creator = swim_proc.ResultsFileCreator(spheres_model_part, node_to_follow_id, scalar_vars, vector_vars) # NODE HISTORY RESULTS END # CHANDELLIER END # GRADIENT RECOVERY TEST BEGIN import gradient_recovery_analyser as gra p_field = gra.sinus_field() v_field = gra.sinus_vector_field() analyser = gra.GradientRecoveryAnalyser(p_field, v_field) # GRADIENT RECOVERY TEST END while (time <= final_time): time = time + Dt step += 1 fluid_model_part.CloneTimeStep(time) print("\n", "TIME = ", time) sys.stdout.flush() if DEM_parameters.coupling_scheme_type == "UpdatedDEM": time_final_DEM_substepping = time + Dt
# NODE HISTORY RESULTS BEGIN scalar_vars = [] vector_vars = [DISPLACEMENT] for node in spheres_model_part.Nodes: node_to_follow_id = node.Id node.SetSolutionStepValue(VELOCITY_Z, terminal_velocity) results_creator = swim_proc.ResultsFileCreator(spheres_model_part, node_to_follow_id, scalar_vars, vector_vars) # NODE HISTORY RESULTS END # CHANDELLIER END # GRADIENT RECOVERY TEST BEGIN import taylor_green_cell as tgc import gradient_recovery_analyser as gra p_field = gra.sinus_field() v_field = tgc.taylor_green_vector_field() analyser = gra.GradientRecoveryAnalyser(p_field, v_field) # GRADIENT RECOVERY TEST END while (time <= final_time): time = time + Dt step += 1 fluid_model_part.CloneTimeStep(time) print("\n", "TIME = ", time) sys.stdout.flush() if DEM_parameters.coupling_scheme_type == "UpdatedDEM": time_final_DEM_substepping = time + Dt