  We study the evolution and oscillations of fixed massive neutrinos
interacting with stochastic gravitational waves (GWs). The energy spectrum of
these GWs is Gaussian, with the correlator of the amplitudes being arbitrary.
We derive the equation for the density matrix for flavor neutrinos in this
case. In the two flavors approximation, this equation can be solved
analytically. We find the numerical solution for the density matrix in the
general case of three neutrino flavors. We consider merging binary black holes
as sources of stochastic GWs with realistic spectra. Both normal and inverted
mass orderings are analyzed. We discuss the relaxation of the neutrino fluxes
in stochastic GWs emitted mainly by supermassive black holes. In this
situation, we obtain the range of energies and the propagation lengths for
which the relaxation process is the most efficient. We discuss the application
of our results for the observation of fluxes of astrophysical neutrinos.
