  We derive the Enskog equation utilizing orthonormal vielbein fields, enabling
the utilization of arbitrary coordinate systems to characterize spatial
geometry. Additionally, we employ an adapted coordinate system in the momentum
space, connected to the physical space through vielbeins. Within this
framework, the momentum component perpendicular to a curved boundary can be
treated as an independent one, facilitating the application of half-range
Gauss-Hermite quadratures. We develop an appropriate finite-difference Lattice
Boltzmann model and validate it against a DSMC-like particle-based method for
solving the Enskog equation in curvilinear geometries. Our test scenarios
include cylindrical Couette flow, cylindrical Fourier flow between coaxial
cylinders, and spherical Fourier flow between concentric spheres. Excellent
agreement between the two approaches is observed throughout the parameter range
and curvature-specific effects are well captured.
