Solve_FSI = True, # enter "True" if you want to solve for fluid-structure interactionĬompressible_solid = True, # enter "True" if compressible: Also remember to specify compressibility (Ld) # enter "False" if incompressible Solve_temperature = True, # enter "True" if you want to solve for temperature Restart = False # Restart parameter # Physics of the problem #. solve_solid_temperature_lagrange_multiplier (a7, LmTs_, b7 ) assemble_solid_temperature_lagrange_multiplier (Ts_, uf_, dt ) start ( ) # print(BLUE % "7: Solid temperature based lagrange multiplier step", flush = True) if problem_physics and problem_physics = True :Ī7, b7 = solid_temp. assign (interpolate_nonmatching_mesh_delta (fsi_interpolation, T_, FS, interpolation_fx, "S" ) ) stop ( ) # - if problem_physics and problem_physics = True : solve_lagrange_multiplier (a6, Lm_, b6 ) assemble_lagrange_multiplier (Lm_, us_, uf_, dt ) start ( ) # print(BLUE % "6: Lagrange multiplier (fictitious force) step", flush = True) if problem_physics = True :Ī6, b6 = lagrange. mesh, problem_physics, a5, Dp_, mix, ps_, p_, bcs )ĭp_. mesh, problem_physics, a5, Dp_, mix, ps_, p_, bcs ) except : assemble_solid_problem (problem_physics, Dp_, mix, uf_, Lm_, dt ) try : step", flush = True) if problem_physics = True :Ī5 = solid. start ( ) # print(BLUE % "5: Solid momentum eq. assign (interpolate_nonmatching_mesh_delta (fsi_interpolation, uv, FS, interpolation_fx, "S" ) ) assemble_temperature (T_, uv, LmTf_, dt )įlow_temp. start ( ) # print(BLUE % "4: Energy conservation step", flush = True) if problem_physics = True :Ī4, b4 = flow_temp. assign (interpolate_nonmatching_mesh_delta (fsi_interpolation, LmTs_, FS, interpolation_fx, "F" ) ) assign (uv, for ui in range (u_components ) ] ) # - if problem_physics and problem_physics = True : solve_velocity_correction (u_, b3, bcs )Īssigner_uv. assemble_velocity_correction (u_, p_, dt )įlow. start ( ) # print(BLUE % "3: Velocity correction step", flush = True)ī3 = flow. solve_pressure_correction (p_, b2, bcs ) assemble_pressure_correction (u_, p_, Lm_f, dt )įlow. start ( ) # print(BLUE % "2: Pressure correction step", flush = True)ī2 = flow. solve_tentative_velocity (A1, u_, b1, bcs ) assemble_tentative_velocity (u_, p_, Lm_f, dt )įlow. start ( ) # print(BLUE % "1: Predict tentative velocity step", flush = True)Ī1, b1 = flow. assign (interpolate_nonmatching_mesh_delta (fsi_interpolation, Lm_, FS, interpolation_fx, "F" ) ) # Update boundary conditions : only if time-dependent # parabolic_profile.t = t tim.t = t num_cycle.cycle = int(t / t_period) # for ui, value in ems(): # inflow.v = evaluate_boundary_val(param_LSPV) inflow.v = evaluate_boundary_val(param_LIPV) # inflow.v = evaluate_boundary_val(param_RSPV) inflow.v = evaluate_boundary_val(param_RIPV) if problem_physics = True : Update_counter (counters ) # Update current time # Time loop try : while T > tsp and t < T :
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |