Modeling fractional behaviors is fundamental for characterizing complex phenomena governed by power-law dynamics, long-range temporal and spatial correlations, fractal structures, and memory effects, which are commonly observed in diverse scientific and engineering applications. In this work, we provide an analytical derivation of a fractional formulation for steady-state heat-conduction modeling. The boundary-value problem is defined on an annular domain, with boundary conditions prescribed as functions of the angular coordinate. The proposed approach integrates two analytical techniques, namely separation of variables and fractional power-series expansion, within a framework based on the Caputo fractional derivative. Beyond heat conduction, the methodology presented here is also directly applicable to a broader class of boundary-value problems governed by elliptic fractional differential equations, supporting advancements in modeling heterogeneous materials, porous media, biological tissues, and systems with multiscale or memory-driven behavior. A distinguishing aspect of this study is its emphasis on analytical transparency and didactic clarity, with each mathematical step introduced and developed systematically. The resulting closed-form solution provides a useful reference for future investigations into fractional differential equations, offering both methodological guidance and conceptual insight. The results include the explicit solution in polar coordinates and its interpretation as a function of an arbitrarily selected fractional order, thereby illustrating how fractional operators modulate steady-state heat conduction under the specified boundary conditions.