Разработан двумерный в координатном пространстве и трехмерный в пространстве скоростей численный код, основанный на методе частиц в ячейке, для моделирования мультипакторного разряда на диэлектрике, помещенном в плоскопараллельный волновод, с учетом вторичной электронной эмиссии с поверхности диэлектрика и стенок волновода, конечной температуры вторичных электронов, пространственного заряда электронов, а также упругого и неупругого отражения электронов от диэлектрика и стенок волновода. Код позволяет моделировать различные стадии мультипакторного разряда, начиная с возникновения электронной лавины и заканчивая стадией насыщения. Показано, что из-за утечки электронов на стенки волновода порог развития одностороннего мультипактора на диэлектрике, помещенном в зауженный волновод с абсорбирующими стенками, увеличивается по сравнению с порогом в случае неограниченного диэлектрика. Обнаружено, что, в зависимости от напряженности СВЧ-поля и вторично-эмиссионных характеристик стенок волновода, могут реализоваться два режима мультипакторного разряда. В первом режиме, который реализуется при относительно низких напряженностях СВЧ-поля, мультипакторный разряд развивается только на диэлектрике, поверхность которого заряжается положительно относительно стенок волновода. Во втором режиме, который реализуется при достаточно высоких напряженностях СВЧ-поля, происходит одновременное возбуждение одностороннего мультипактора на поверхности диэлектрика и двустороннего мультипактора между стенками волновода. В этом случае поверхность диэлектрика и пространство между стенками волновода приобретают отрицательный потенциал. Показано, что учет отражения электронов от поверхности диэлектрика и стенок волновода приводит к появлению высокоэнергетических хвостов функции распределения электронов.