def METHOD_NAME(e, nu, inn_radius, out_radius, press, aesize):
    """Create 2D cross section modeling a pipe."""
    # reset mapdl
    mapdl.clear()
    mapdl.prep7()
    # Define element attributes
    # Quad 4 node 182 with keyoption 3 = 2 (plain strain formulation)
    mapdl.et(1, "PLANE182", kop3=2)
    # Create geometry
    # create a quadrant of the pressure vessel
    # We perform plane strain analysis on one quadrant (0deg - 90deg) of the
    # pressure vessel
    mapdl.pcirc(inn_radius, out_radius, theta1=0, theta2=90)
    mapdl.cm("PIPE_PROFILE", "AREA")
    # Define material properties
    mapdl.mp("EX", 1, e)  # Youngs modulus
    mapdl.mp("PRXY", 1, nu)  # Poissons ratio
    # Define mesh controls
    mapdl.aesize("ALL", aesize)
    mapdl.mshape(0, "2D")  # mesh the area with 2D Quad elements
    mapdl.mshkey(1)  # free mesh
    mapdl.cmsel("S", "PIPE_PROFILE")  # Select the area component to be meshed
    mapdl.amesh("ALL")
    # Create components for defining loads and constraints
    mapdl.nsel("S", "LOC", "X", 0)  # Select nodes on top left edge
    mapdl.cm("X_FIXED", "NODES")  # Create nodal component
    mapdl.nsel("S", "LOC", "Y", 0)  # Select nodes on bottom right edge
    mapdl.cm("Y_FIXED", "NODES")  # Create nodal component
    mapdl.allsel()
    mapdl.lsel("S", "RADIUS", vmin=rad1)  # Select the line along inner radius
    mapdl.cm("PRESSURE_EDGE", "LINE")  # Create a line component
    mapdl.allsel()
    # Define solution controls
    mapdl.slashsolu()  # Enter solution
    mapdl.antype("STATIC", "NEW")  # Specify a new static analysis (Optional)
    mapdl.d("X_FIXED", "UX", 0)  # Fix the selected nodes in X direction
    mapdl.d("Y_FIXED", "UY", 0)  # Fix the selected nodes in Y direction
    # Change the active Cartesian Coordinate system to Cylindrical Coordinate system
    mapdl.csys(1)
    # Apply uniform pressure load to the selected edge
    mapdl.sfl("PRESSURE_EDGE", "PRES", press)
    # Solve the model
    mapdl.allsel()
    mapdl.solve()
    mapdl.finish()
    # Enter post-processor
    mapdl.post1()
    mapdl.set(1, 1)  # Select the first load step
    max_eqv_stress = np.max(mapdl.post_processing.nodal_eqv_stress())
    all_dof = mapdl.mesh.nnum_all
    num_dof = all_dof.size
    return num_dof, max_eqv_stress