# 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
Esempio n. 2
0
# 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