Euler–Lagrange simulations of particle-laden flow require hydrodynamic models of drag and lift forces for individual particles. Our goal is to develop models that can prescribe these forces for arbitrarily orientated shell objects. Here, we use computational fluid dynamics simulations of steady bottom-boundary layer flow over a series of spherical triangle shell fragments to calculate the hydrodynamic forces. The simulations explicitly resolve the wall boundary layers using grid resolution on the order of y+=1 at the shell fragment surface and use the SST k-omega turbulence closure model. These fragments cover a range of aspect ratio and flatness characteristics. The shell fragments are generated as triangular selections of a spherical shell with azimuthal and longitudinal angles proscribed based on elongation and flatness parameters (varying between 1 to 5, and 0.02 to 0.2 respectively), while characteristic length of the fragment is held constant to define the overall fragment size. Fragment orientations are considered with independently varying pitch, roll, and yaw each ranging from 0 to 180 degrees. The numerical estimates for the forces from all simulations were used to develop robust parameterizations of the drag and lift as a function of aspect ratio and flatness characteristics, as well as orientation of the shell fragments.