Моделирование мелкого гигантского столкновения требует короткого временнóго шага в компьютерных расчетах; таким образом, нужно намного больше циклов, чтобы программа завершила работу. А теперь, со всеми этими новыми физическими соображениями, для каждого интервала требовалось вычислить не только давление, гравитацию и температуру, но и степень разрушения, уплотнение и тензор напряжений. Давление – это скаляр (то есть просто число), который соответствует величине тензора напряжений. Тензор напряжений содержит не лежащие на осях показатели, которые позволяют твердым телам сопротивляться деформации. Сопротивление сдвигу – это сопротивление движению в направлении оси y по плоскости, параллельной оси х (сопротивление скольжению). Такое сопротивление сдвигу обозначают как syx и так далее для других плоскостей. Так что в твердом теле есть девять компонентов напряжения (sxx, sxy, sxz, syx, syy, syz, szx, szy, szz). Давление определяется как все напряжение, приложенное на те же грани (то есть все, кроме сдвигающего усилия), поэтому P = (sxx + syy + szz)/3. Далее нам нужно определить реакцию вещества на напряжение и давление – это реология и уравнение состояния. Вся эта математика относится к элементарной алгебре, но это лишь аппроксимация, и поэтому в коде нарастает вероятность ошибок, а расчеты требуют все больше времени.
