  The azimuthal anisotropy of $\Upsilon$(1S) mesons in high-multiplicity
proton-lead collisions is studied using data collected by the CMS experiment at
a nucleon-nucleon center-of-mass energy of 8.16 TeV. The $\Upsilon$(1S) mesons
are reconstructed using their dimuon decay channel. The anisotropy is
characterized by the second Fourier harmonic coefficients, found using a
two-particle correlation technique, in which the $\Upsilon$(1S) mesons are
correlated with charged hadrons. A large pseudorapidity gap is used to suppress
short-range correlations. Nonflow contamination from the dijet background is
removed using a low-multiplicity subtraction method, and the results are
presented as a function of $\Upsilon$(1S) transverse momentum. The azimuthal
anisotropies are smaller than those found for charmonia in proton-lead
collisions at the same collision energy, but are consistent with values found
for $\Upsilon$(1S) mesons in lead-lead interactions at a nucleon-nucleon
center-of-mass energy of 5.02 TeV.
