Abstract Axisymmetric groundwater models are used for simulating radially symmetric conditions. Groundwater simulators built specifically to model axisymmetric conditions are most commonly used for simulating aquifer tests. Although some numerical models capable of simulating flow and solute transport that are developed in the cartesian coordinate system framework offer flexibility to simulate axisymmetric conditions, most of the numerical groundwater models, such as the MODFLOW family of codes, are based on structured grids in which axisymmetric flows cannot be directly simulated. Researchers in the past have provided methods to manipulate aquifer properties to mimic axisymmetric conditions. This study presents a methodology that takes advantage of the unstructured grids of MODFLOW‐USG to simulate axisymmetric models within the MODFLOW framework. To develop axisymmetric models, the intercell interface area arrays of MODFLOW‐USG were calculated to accurately represent coaxial cylindrical model cells. Three examples are presented to demonstrate the application of MODFLOW‐USG for axisymmetric modeling: a pumping well with delayed yield effects, a vadose zone flow model simulating an infiltration basin, and a density‐dependent saltwater intrusion problem for a circular island. Results were verified against analytical solutions and published numerical codes.