Первое крупное вычисление
Первая крупномасштабная симуляция Ричарда Хэмминга: Лос-Аламос, 1945. Цель — спроектировать рабочую атомную бомбу.
Проблема, которая сделала симуляцию необходимой: не существует маломасштабного эксперимента. Критическая масса бинарна. Либо расщепляющийся материал превышает критическую массу и происходит цепная реакция, либо нет. Вы не можете провести уменьшенную версию.
Конструкция сферического взрыва
Один из проектов использовал сферическую симметрию — имплозию. Инженеры разделили материал на множество концентрических оболочек. Для каждой оболочки они написали уравнения для сил на обеих сторонах, плюс уравнение состояния, связывающее давление с плотностью.
Время было дискретизировано на интервалы 10⁻⁸ секунд — называемые 'встряски' (от 'a shake of a lamb's tail'). На каждой встряске вычисление продвигалось вперед: куда движется каждая оболочка? Какие силы на нее действуют?
Три условия, которые требуют симуляции
Хэмминг выявил ситуации, где симуляция заменяет физический эксперимент:
1. Невозможные эксперименты — критическая масса не может быть протестирована в меньшем масштабе
2. Опасные эксперименты — вы не можете взорвать бомбу для данных калибровки
3. Слишком дорогие или слишком медленные — блокировка атмосферы, предсказание погоды, траектории ракет
Цель: получить эквивалентные результаты, а не точно копировать физический процесс. Симуляция не должна соответствовать реальности атом за атомом. Она должна производить те же наблюдаемые результаты в пределах точности, которая требуется для проекта.
Эквивалентные результаты
Ключное открытие Хэмминга в Лос-Аламосе: данные уравнения состояния были неточными. Соотношения давления и плотности поступали из лабораторий высокого давления, оценок землетрясений, моделей звездных ядер — все с значительной неопределенностью.
Инженеры читали эти кривые до трех десятичных мест, а затем табулировали их до пяти цифр. На первый взгляд, мусор на входе.
Но проект бомбы сработал.
Почему? Потому что вычисление использовало вторые разности значений на соседних оболочках. Любая локальная ошибка в уравнении состояния усреднялась на протяжении истории оболочки при прохождении кривой. Имела значение: кривизна уравнения состояния, и только в среднем.
Обратная связь в вычислении компенсировала неточные входные данные.
Повторяющееся ядро
Хэмминг выявил универсальную структурную особенность больших симуляций: весьма повторяющийся внутренний цикл.
В Лос-Аламосе: одни и те же уравнения силы работали для каждой оболочки на каждом временном шаге. Код для одной оболочки выполнялся тысячи раз. Без этой повторяющейся структуры стоимость программирования была бы запредельной.
Этот принцип обобщается: предсказание погоды разбивает атмосферу на блоки; одни и те же физические уравнения обновляют каждый блок. Симуляции ракет пошагово проходят одни и те же уравнения траектории на каждом временном приращении. Проектирование транзисторов вычисляет одни и те же полевые уравнения в каждой точке пространственной сетки.
Совет Хэмминга: рано обратите внимание на повторяющиеся части любой предложенной симуляции. Сформулируйте симуляцию в форме, которая использует повторение. Симуляция без плотного внутреннего цикла, вероятно, плохо структурирована.
Экспертное знание как обязательное условие
Хэмминг неоднократно возвращался к правилу, которое считал неоспоримым: только эксперт в области может знать, что имеет значение.
Эксперт по симуляции может структурировать код, выбрать численные методы, отладить повторяющийся цикл. Но только человек, свободно владеющий физикой, химией или инженерией в области, может определить:
- Какие эффекты должны появиться в модели
- Какие можно безопасно опустить
- Является ли необычный результат сигналом физической истины или ошибки моделирования
В Лос-Аламосе Хэмминг был экспертом по вычислениям. Физики были экспертами в области. Ни один не мог заменить другого.
Жаргон как барьер и инструмент
Одно из самых сильных операционных правил Хэмминга: выучите жаргон области, которую вы моделируете.
Его история: задача перехвата ВМФ с 28 одновременными дифференциальными уравнениями. Он настаивал, чтобы предложивший — его друг-физик — прошелся вместе с ним по каждой строке двоичного машинного кода перед запуском вычисления.
В середине разбора физик сказал: 'Дик, это ограничение плавника, а не ограничение напряжения.' Один и тот же математический символ, одно и то же формальное уравнение — но два разных физических толкования, дающие существенно разные результаты.
Урок: обе стороны понимали математику. Ни у одного не было отказа в коммуникации в обычном смысле. Но физический смысл операции ограничения был недостаточно определен одними уравнениями.
Без разбора симуляция работала бы с неправильным толкованием. Никакой ошибки времени выполнения. Никакого очевидно плохого результата. Просто неправильные ответы на важную проблему.
Устойчивые и неустойчивые задачи
Хэмминг провел четкую границу между задачами, которые симуляция хорошо решает, и задачами, где симуляция почти невозможна.
Атомная бомба: небольшие различия в траекториях оболочек по пути не сильно повлияли на финальный выход. Симуляция была устойчива — ошибки в промежуточных шагах не усиливались.
Предсказание погоды: противоположное. Небольшое возмущение — 'машет ли бабочка крыльями в Японии' — может, в принципе, определить, попадет ли шторм на континент. Чувствительность к начальным условиям делает день-за-днем прогнозирование погоды ненадежным за пределами коротких горизонтов.
Обе задачи используют одинаковую математическую структуру: дискретизировать пространство на ячейки, дискретизировать время на шаги, маршировать вперед. Разница заключается в поле направлений — растут ли (неустойчивы) или сжимаются (устойчивы) небольшие отклонения от траектории со временем.
Погода показывает оба: краткосрочная неустойчивость (день-за-днем хаос), долгосрочная устойчивость (сезоны идут своим чередом) и очень долгосрочная неустойчивость (ледниковые периоды).
Правило Хэмминга: перед началом любой симуляции проверьте, является ли задача в фундаментальном смысле устойчивой или неустойчивой. Если неустойчивая, определите, является ли неустойчивость фундаментальной для нужного вам ответа или артефактом масштаба или граничных условий. Не обнаруживайте это после инвестирования месяцев работы.
Сначала просто, потом полностью
Предпочтительный метод Хэмминга для подхода к новой симуляции:
1. Начните просто — включите только основные эффекты. Правильно получите доминирующее поведение.
2. Получите понимание рано — простая симуляция раскрывает структуру задачи до того, как вы вложитесь в полную детализацию.
3. Развивайте в сторону полноты — добавляйте вторичные эффекты постепенно, проверяя каждое дополнение против более простого базового варианта.
Он проиллюстрировал это проектом ракеты Nike. Его ранние симуляции использовали простую экспоненциальную модель атмосферы. Позже симуляции добавили профили плотности, зависящие от высоты, члены перекрестного ветра и нелинейное сопротивление. Но ранние понимания — вертикальный запуск снижает сопротивление в плотной нижней атмосфере; более крупные крылья стоят больше скорости, чем они получают в поздней маневренности — пришли из простой модели.
Предупреждение: в конце заморозьте проект, используя полную симуляцию. Простая симуляция получает понимание; полная симуляция получает обязательство.