formatted the code

Author: krishnatejas12

examples/00-systemcoupling/cht_pipe.py

@@ -43,7 +43,7 @@
 
 
 The flow is smooth inside the pipe and the outer wall of the pipe is adiabatic. Fluid enters
-at an initial temperature of 300K while the outside pipe of the wall is at 400K. 
+at an initial temperature of 300K while the outside pipe of the wall is at 400K.
 The setup is simulated for a few iterations to allow the examination
 of the temperature field in the pipe.
 
@@ -58,8 +58,9 @@
 # ~~~~~~~~~~~~~~~~~~~~~~~~
 # Import ``ansys-systemcoupling-core``, ``ansys-fluent-core``
 
-import ansys.mapdl.core as pymapdl
 import ansys.fluent.core as pyfluent
+import ansys.mapdl.core as pymapdl
+
 import ansys.systemcoupling.core as pysyc
 from ansys.systemcoupling.core import examples
 
@@ -69,101 +70,96 @@
     "fluid_domain.msh", "pysystem-coupling/cht_pipe"
 )
 
-mapdl=pymapdl.launch_mapdl()
+mapdl = pymapdl.launch_mapdl()
 mapdl.clear()
 mapdl.prep7()
 
-#%%
-#Define material properties
+# %%
+# Define material properties
 mapdl.mp("EX", 1, 69e9)
 mapdl.mp("NUXY", 1, 0.33)
 mapdl.mp("DENS", 1, 2700)
 mapdl.mp("ALPX", 1, 23.6e-6)
 mapdl.mp("KXX", 1, 237)
-mapdl.mp("C", 1 , 900)
+mapdl.mp("C", 1, 900)
 
-#%%
-#Set element type to SOLID279
+# %%
+# Set element type to SOLID279
 mapdl.et(1, 279)
 mapdl.keyopt(1, 2, 1)
 print(mapdl)
 
-#%%
-#Parameter of the pipe
-r_in=0.025
-r_out=0.035
-l=0.2
+# %%
+# Parameter of the pipe
+r_in = 0.025
+r_out = 0.035
+l = 0.2
 
-#%%
-#Create hollow pipe
-mapdl.cyl4(0, 0, rad1=r_in, rad2=r_out, depth=l)  
+# %%
+# Create hollow pipe
+mapdl.cyl4(0, 0, rad1=r_in, rad2=r_out, depth=l)
 mapdl.esize(0.002)
 mapdl.vsweep(1)
 print(mapdl.geometry.anum)
 
+
 # %%
 # Biot number prediction
-def biot_number(
-        rho=1000,
-        mu=1e-3,
-        cp=4180,
-        k_f=0.6,
-        k_s=237,
-        L_c=r_out-r_in,
-        U=0.1
-):
-    Re= (rho*U*L_c)/mu
-    Pr= (mu*cp)/ k_f
-    Nu= 0.023 * Re**0.8 * Pr**0.4
-    h= Nu * k_f / L_c
-    Bi= h *L_c / k_s
+def biot_number(rho=1000, mu=1e-3, cp=4180, k_f=0.6, k_s=237, L_c=r_out - r_in, U=0.1):
+    Re = (rho * U * L_c) / mu
+    Pr = (mu * cp) / k_f
+    Nu = 0.023 * Re**0.8 * Pr**0.4
+    h = Nu * k_f / L_c
+    Bi = h * L_c / k_s
     return Bi
-Bi=biot_number()
+
+
+Bi = biot_number()
 print("The Biot number is ", Bi)
 
 # %%
 # Creating the regions from the geometry for named selections
-#Inner wall NS
-mapdl.asel('S', 'AREA', '', 5, 6)
-mapdl.nsla('S', 1)
+# Inner wall NS
+mapdl.asel("S", "AREA", "", 5, 6)
+mapdl.nsla("S", 1)
 mapdl.cm("FSIN_1", "NODE")
 mapdl.allsel()
 
-#Outer wall NS
-mapdl.asel('S', 'AREA', '', 3, 4)
+# Outer wall NS
+mapdl.asel("S", "AREA", "", 3, 4)
 mapdl.cm("Outer_wall", "AREA")
 mapdl.allsel()
 
-#Outlet NS
-mapdl.asel('S', 'AREA', '', 2)
+# Outlet NS
+mapdl.asel("S", "AREA", "", 2)
 mapdl.cm("Outlet", "AREA")
 mapdl.allsel()
 
-#Inlet NS
-mapdl.asel('S', 'AREA', '', 1)
+# Inlet NS
+mapdl.asel("S", "AREA", "", 1)
 mapdl.cm("Inlet", "AREA")
 mapdl.allsel()
 
-#%%
-#Boundary conditions
-mapdl.cmsel('S', 'Outer_wall')
-mapdl.d('Outer_wall', 'TEMP', 77)
+# %%
+# Boundary conditions
+mapdl.cmsel("S", "Outer_wall")
+mapdl.d("Outer_wall", "TEMP", 77)
 mapdl.allsel()
 
-mapdl.cmsel('S', 'Inlet')
-mapdl.sf('ALL', 'HFLUX', 0)
+mapdl.cmsel("S", "Inlet")
+mapdl.sf("ALL", "HFLUX", 0)
 mapdl.allsel()
 
-mapdl.cmsel('S', 'Outlet')
-mapdl.sf('ALL', 'HFLUX', 0)
+mapdl.cmsel("S", "Outlet")
+mapdl.sf("ALL", "HFLUX", 0)
 mapdl.allsel()
 
-mapdl.cmsel('S', 'FSIN_1')
-mapdl.sf('FSIN_1', 'FSIN', 1)
+mapdl.cmsel("S", "FSIN_1")
+mapdl.sf("FSIN_1", "FSIN", 1)
 mapdl.allsel()
 
-#%%
-#Setup the rest of the analysis
+# %%
+# Setup the rest of the analysis
 mapdl.run("/SOLU")
 mapdl.antype(0)
 
@@ -174,24 +170,26 @@ def biot_number(
 # %%
 # Read the pre-created mesh file
 
-fluent=pyfluent.launch_fluent(start_transcript=False)
+fluent = pyfluent.launch_fluent(start_transcript=False)
 fluent.file.read(file_type="mesh", file_name=fluent_msh_file)
 
-#%%
-#Define the fluid solver settings
-fluent.setup.models.energy.enabled =True
+# %%
+# Define the fluid solver settings
+fluent.setup.models.energy.enabled = True
 
-#%%
-#Add the material
+# %%
+# Add the material
 fluent.setup.materials.database.copy_by_name(type="fluid", name="water-liquid")
 
 fluent.setup.cell_zone_conditions.fluid["fff_fluiddomain"].material = "water-liquid"
 
-#%%
-#Define boundary conditions
+# %%
+# Define boundary conditions
 fluent.setup.boundary_conditions.velocity_inlet["inlet"].momentum.velocity = 0.1
 fluent.setup.boundary_conditions.velocity_inlet["inlet"].thermal.temperature = 300
-fluent.setup.boundary_conditions.wall["inner_wall"].thermal.thermal_bc = "via System Coupling"
+fluent.setup.boundary_conditions.wall["inner_wall"].thermal.thermal_bc = (
+    "via System Coupling"
+)
 
 fluent.solution.run_calculation.iter_count = 20
 
@@ -201,33 +199,35 @@ def biot_number(
 # System Coupling setup involves adding the structural and fluid
 # participants, adding coupled interfaces and data transfers,
 # and setting other coupled analysis properties.
-syc=pysyc.launch(start_output= True)
+syc = pysyc.launch(start_output=True)
 
-#%%
+# %%
 # Add participants by passing session handles to System Coupling.
-fluid_name=syc.setup.add_participant(participant_session= fluent)
-solid_name=syc.setup.add_participant(participant_session= mapdl)
+fluid_name = syc.setup.add_participant(participant_session=fluent)
+solid_name = syc.setup.add_participant(participant_session=mapdl)
 
-syc.setup.coupling_participant[fluid_name].display_name= "Fluid"
-syc.setup.coupling_participant[solid_name].display_name= "Solid"
+syc.setup.coupling_participant[fluid_name].display_name = "Fluid"
+syc.setup.coupling_participant[solid_name].display_name = "Solid"
 
-#%%
-#Add coupling face and data transfers
-interface_name=syc.setup.add_interface(
-    side_one_participant= fluid_name, side_one_regions= ["inner_wall"],
-    side_two_participant= solid_name, side_two_regions= ["FSIN_1"]
+# %%
+# Add coupling face and data transfers
+interface_name = syc.setup.add_interface(
+    side_one_participant=fluid_name,
+    side_one_regions=["inner_wall"],
+    side_two_participant=solid_name,
+    side_two_regions=["FSIN_1"],
 )
 
-#%%
-#Set up 2-way thermal FSI coupling
-syc.setup.add_thermal_data_transfers(interface = interface_name)
+# %%
+# Set up 2-way thermal FSI coupling
+syc.setup.add_thermal_data_transfers(interface=interface_name)
 
 # %%
 # Time step size, end time, output controls
 syc.setup.solution_control.time_step_size = "0.1 [s]"  # time step is 0.1 [s]
 syc.setup.solution_control.end_time = 10  # end time is 10.0 [s]
 
-syc.setup.solution_control.maximum_iterations=100
+syc.setup.solution_control.maximum_iterations = 100
 
 # %%
 # Solution
@@ -239,4 +239,3 @@ def biot_number(
 # ----
 syc.end_output()
 syc.exit()
-