Устройство для решения дифференциальных уравнений

д Куцен 8. Р РЕШЕНИЯ АВНЕ НИЙ числитель- пользовано ГОСУДАРСТВЕННЫЙ КОМИТЕТ СССРПО ДЕЛАМ ИЗОБРЕТЕНИЙ И ОТКРЫТИЙ К АВТОРСКОМУ СВИДЕТЕЛЬСТВУ(54) УСТРОЙСТВО ДЛЯДИФФЕРЕНЦИАЛЬНЫХ УР(57) Изобретение относится к вной технике и может быть ис ля решения дифференциальных уравнений в частных производных, Цель изобретения - повышение быстродействия устройства. Поставленная цель достигается тем, что устройство содержит блокввода-вывода, блок 2 управления, два блока 31 и 3 местного управления, % блоков 41 - 4 двунаправленной передачи данных, где У - линейный размер области определения уравнения, У решающих блоков 5 - 5 1 Ч групп по л и (У+1)-ю группу из и- блоков 6 однонаправленной передачи данных, где ив размерность уравнения,и блок 7 вычислений, 11 ил.Редактор Н. ТупицЗаказ 915/49ВНИИПИ Государ1130Производственно Смир итель Верес уска Р по делам изоушская наб.,иятие, г. Ужгор тенин4/5ул. П ткрытий ктная, 4 Состав а Техред И. Тираж 704 твенного комитета ССС 5, Москва, Ж - 35, Ра олиграфическое предпрИзобретение относится к вычислительной технике и может быть использовано для решения дифферециальных уравнений в частных производных.Цель изобретения - повышение быстродействия устройства.На фиг. 1 представлена структурная схема устройства; на фиг. 2 - схема блока управления; на фиг. 3 и 4 - схемы блока местного управления; на фиг: 5 - схема блока двунаправленной передачи данных; на фиг. 6 - схема решающего блока; на фиг. 7 - схема блока однонаправленной передачи данных; на фиг. 8 - схема блока вычислений; на фиг: 9 и 10 - алгоритмы работы устройства; на фиг. 11 - временная диаграмма работы устройства. Устройство содержит блок 1 ввода-вывода, блок 2 управления, первый З и второй 32 блоки местного управления, Л блоков 4, - 4 двунаправленной передачи данных, У решающих блоков 51 - 5 М групп блоков 6 - б однонаправленной передачи данных, блок 7 вычислений, системные выходные шины 8, системные входные шины 9, входные шины 10 блока 3 местного управления, выходные шины 11 блока 3 местного управления, выходные шины 12 решающего блока 5, входные шины 13 решающего блока 5, первые выходные шины 14 блока 7, первые входные шины 15 блока 7, вторые входные шины 16 матричного вычислителя 7, вторые выходные шины 17 матричного вычислителя 7,Блок 2 управления может быть выполнен по известной схеме, и содержит узел 18 памяти, регистр 19 адреса, буфер 20 ввода- вывода, счетчик 21 команд, регистр 22 команд, буферный регистр 23, дешифратор 24 команд, узел 25 синхронизации.Блок 3 местного управления может быть выполнен по известной схеме и содержит узел 26 памяти, регистр 27 адреса, буфер 28 ввода-вывода, счетчик 29 команд; регистр ЗО команд, буферный регистр 31, дешифратор 32 команд, узел 33 синхронизации, элемент И 34, дешифратор 35 адреса, элемент И 36, группу ключей 37, регистр 38 вывода, регистр 39 ввода, группу ключей 40, элемент И 41, элемент И 42, дешифратор 43 адреса.Блок 4 двунаправленной передачи данных может быть выполнен по известной схеме и содержит регистр 44 вывода, группу клк)чей 45, элементы И 46 и 47, группу ключей 48, регистр 49 ввода, элемент И 50, дешифратор 51 адреса, дешифратор 52 адреса, элемент И 53.Решающий блок 5 может быть выполнен по известной схеме и содержит регистр 54 результата, узел 55 синхронизации, регистр 56 второго операнда, регистр 57 первого операнда, дешифратор 58 команд, счетчик 59 команд, сумматор 60, регистр 61 команд, регистр 62 адреса, регистр 63 частичного1(ант а) а 1 ащ аФщКе дх ду ду дх д хдудФ+д Р О (б)дх ду 2+ = - фх, у, г), (в) удовлетворяющей на ЭО условиям прилипания и непроницаемости:(3) Т=То 45на границе охлаждения;2) второго родаУ,=О (а=1, 2, 3) (4)для боковых поверхностей параллелепипеда.50 Для системы (1) исходными являютсяпервое и второе уравнения, с помощью которых определяются конвективные члены (1l=- и Г=), а результатом решения зфду дхраспределение температур (Т) в области О.Введя в область О разностную сетку и заменив производные разностными отношениями, т. е. для аппроксимации оператора результата, группу ключей 64, узел 65 памяти.Блок 6 однонаправленной передачи данных может быть выполнен по известной схеме и содержит регистр 66, группу ключей 67, элементы И 68 и 69, дешифраторы 70 и 71 адреса.Блок 7 вычислений может быть выполненпо известной схеме и содержит группу решающих блоков 72 (выполненных аналогично решающему блоку 5), группу блоков 73 двунаправленной передачи данных и две группы входных и выходных шин 74, 75 и 76, 77 блоков 73.Устройство работает следующим образом.15 В исходном состоянии устройство подготовлено к работе, т. е, в узел 18 памяти блока 2 управления записаны через блок 1 ввода управляющие программы, в узел 26 памяти блоков 31 и 32 местного управления записаны блоком 2 управляющие программы, 20 массивы коэффициентов, начальные и граничные условия, а в узлах 65 памяти решающих блоков 5 и 72 записаны управляющие и прикладные программы.Рассмотрим работу устройства на приме ре решения задачи о стационарном течениивязкой жидкости, в терминах безразмерных переменных функция тока - вихрь, в ограниченной области О, которая описывается системой дифференциальных уравнений20 Алгоритм решения задачи включает следующие этапы:1. Загрузка исходных данных из блоков 3 и Зчерез блоки 4 - 4 и 6 - 6, в решающи блоки 51 - 5, и в блок 7 вычислений; 2. Решение разностного аналога уравнения движения (1 а) методом оптимальной релаксации в блоке 7 (в соответствии с формулой (9) определяется значение ь внутри области);3. Решение разностного аналога уравнения Пуассона (1 б) методом оптимальной релаксации в блоке 7 (в соответствии с формулой (10) определяется уточненное ф;4. Определение значения вихря на границе (ж,) методом оптимальной релаксации с использованием регуляризации в блоке 7 (в соответствии с формулами 11 и 12);5. Определяется невязка для каждой точки области 0 для значений . Если невязка в допустимых пределах, то переход к п. 6, в противном случае возврат к п. 2;6. Пересылаются полученные значения р внутри области и на границе, а для ь только на границе из блока 7 через блоки 6, - 6 в решающие блоки 51 - 5;7. Вычисляются в решающих блоках 51 - 5 в соответствии с формулой (6) значения / и 1 т;8. Решение разностного аналога уравнения (1 в) методом разделения переменных с использованием БПФ (в соответствии с формулами 13 и 14).Для реализации приведенного алгоритма по командам блока 2 управления блоки 3, и Зг местного управления по собственным программам параллельно загружают исходную информацию (коэффициенты, значения граничных условий) в решающие блоки 72. Для этого блок Зг местного управления выдает по шинам 11 ь 11 г, 11 д соответственно данные, управляющий сигнал Запись (ЗП) и адрес. В блоке бь к которому адресуется блок Зг местного управления, срабатывает дешифратор 70 адреса, с выхода которого активный сигнал совместно с активным сигналом ЗП поступает на входы элемента И 68, По сигналу с выхода элемента И 68 передаваемая информация с шин 111 записывается в регистр 66 и одновременно сигнал Запрос поступает в решающий блок 72 (узел 55 синхронизации), По сигналу Запрос решающий блок 72 переходит на подпрограмму чтения информации из регистра 66. Для этого решающий блок 72 выставляет на шинах 14 г адрес, который поступает на вход дешифратора 71 адреса, а по шине 14, - активный сигнал Чтение (ЧТ). По совокупности активных сигналов на входе элемента И 69 ключи 67 подключают выходы регистра 66 к шинам 151 и информация через группу ключей 64 записывается в регистр 54 результата, а затем запоминается в узле 65 памяти (фиг. 6), При чтении информации 25 30 35 40 45 50 55 из регистра 66 сигнал Запрос по шине 10 г сообщает блоку 3 местного управления, что можно передавать следующее слово. По инициативе блока 3 рассмотренный цикл передачи повторяется. Последней командой передаваемого массива слов блок Зг местного управления запускает первый решающий блок 72на выполнение прикладной программы в соответствии с формулами (9) - (12). Для этого передается в решающий блок 721код, соответствуюший передаче управления, и адрес перехода.Далее блок Зг местного управления переключается через соответствующий блок бг к решающему блоку 72 г ь загружает исходной информацией и передает ему управление. Таким образом, все п решающих блоков 72 со смещением во времени параллельно выполняют вычисления. После выполнения очередной итерации решающие блоки 72 через блоки 73 обмениваются информацией с соседними решающими блоками в соответстви и с ал гор итмом (9) и (10) .Для примера рассмотрим обмен информацией блоками 72 1 и 72 ьг через блок 73, Если обмен инициируется решающим блоком 72 и то на шинах 74, 74 г и 74 д выставляются соответственно данные, сигнал ЗП и адрес. В блоке 73 срабатывает дешифратор 52 адреса и по активному сигналу с выхода элемента И 46 передаваемая информация записывается в регистр 44. По сигналу Запрос (шина 76 г) решающий блок 72 принимает информацию через ключи 45 и шину 761 в узел памяти. Для этого выставляется по шинам 77 г,г и 77 соответственно сигнал ЧТ и адрес блока 73 ьпгл. Передача информации из решающего блока 72 гв 72 и осуществляется аналогично рассмотренному через регистр 49 и ключи 48.После получения решения с определенной точностью в блоке 7 значения ф; и ь пересылаются в решающие блоки 51 - 5, через группы блоков 6. Каждому -му столбцу решающих блоков 72 в блоке 7 соответствует -й решающий блок 5. Во время выполнения вычислений в блоке 7 блок 31 местного управления через блоки 41 - 4 загружает исходной информацией решающие блоки 5 - 5 которые после приема информации из блока 7 приступают к выполнению п. 7 и 8 описанного алгоритма. Более подробно алгоритм работы блоков 5 - 5, приведен на фиг, 10.После того, как определено распределение температур (Т) в текущей плоскости, результаты через блоки 41 - 4 пересылаются в блок 31 для последующего вывода. Во время работы блоков 51 - 5 по расчету значений искомой функции в текущей плоскости блок 7 приступает к расчету ф и ю для следующей плоскости.Формула изобретенияУстройство для решения дифференциальных уравнений, содержащее блок ввода-вы 1383391вода, блок управления, два блока местного управления, Лг блоков двунаправленной передачи данных, где Л/ - линейный размер области определения уравнения, с первого по Лг-й рещающие блоки, первый, второй и третий выходы блока управления подключены соответственно к первым, вторым и к третьим входам режима первого и второго блоков местного управления, первые выходы первого и второго блоков местного управления объединены и подключены к первому входу режима блока управления, вторые выходы первого и второго блоков местного управления объединены и подключены к второму входу режима блока управления, четвертый и пятый выходы блока управления подключены соответственно к информационному и управляющему входам блока ввода- вывода, информационный выход которого подключен к третьему входу режима блока управления, третий, четвертый и пятый выходы первого блока местного управления подключены соответственно к первым информационным, первым управляющим и к вторым информационным входам У блоков двунаправленной передачи данных, первые информационные выходы Лг блоков двунаправленной передачи данных объединены и подключены к четвертому входу режима первого блока местного управления, первые управляющие выходы М блоков двунаправленной передачи данных объединены и подключены к пятому входу режима первого блока местного управления, отличающееся тем, что, с целью повышения быстродействия устройства, оно содержит Лг групп по и блоков однонаправленной передачи данных, где и - размерность уравнения, (Л/+ 1)-ю группу из гг блоков однонаправленной передачи данных и блок вычислений, первый информационный выход К-го решающего блока (К=1, .;., Л) подключен к третьему информационному входу К-го блока двунаправленной передачи данных, второй информационный и управляющий выходы К-го решающего блока подключены, соответственно, к четвертому информационному и к второму управляющему входам каждого блока двунаправленной передачи данных, а также соответственно к первым информационным и к первым управляющим входам и блоков однонаправленной передачи данных К-й группы, первые информационные выходы и блоков однонаправленной передачи данных К-й группы и второй информа 5 10 15 20 25 30 35 40 45 ционный выход К-го блока двунаправленной передачи данных объединены и подключены к информационному входу К-го решающего блока, первые управляющие выходы г блоков однонаправленной передачи данных К-й группы и второй управляюгций выход К-го блока двунаправленной передачи данных объединены и подключены к управляющему входу К-го решающего блока, вторые управляющие выходы гг блоков однонаправленной передачи данных (Лг+1) -й группы объединены и подключены к пятому входу режима второго блока местного управления, третий, четвертый и пятый выходы второго блока местного управления подключены соответственно к вторым информационным, вторым управляющим и к третьим информационным входам и- блоков однонаправленной передачи данных (%+1) -й группы, первые информационные выходы блоков однонаправленной передачи данных с первого по и-й (Лг+1)-й группы подключены соответственно к информационным входам с первого по и- первой группы блока вычислений, первые управляющие выходы блоков однонаправленной передачи данных с первого по и-й (Л+1)-й группы подключены соответственно к управляющим входам с первого по ап первой группы блока вычислений, вторые управляющие выходы и блоков однонаправленной передачи данных К-й группы объединены и подключены к К-му управляющему входу второй группы блока вычислений, первые информационные входы блоков однонаправленной передачи данных с первого по л -й (Лг+ 1)-й группы подключены соответственно к информационным выходам с первого по гг первой группы блока вычислений, первые управляющие входы блоков однонаправленной передачи данных с первого по гг-й (%+1)-й группы подключены соответственно к управляющим выходам с первого по и-й первой группы блока вычислений, К-й информационный выход второй группы блока вычислений подключен к вторым информационным входам и блоков однонаправленной передачи данных К-й группы, К-й управляющий выход второй группы блока вычислений подключен к вторым управляющим входам и блоков однонаправленной передачи данных К-й группы, К-й информационный выход третьей группы блока вычислений подключен к третьим информационным входам п блоков однонаправленной передачи данных каждой группы.1383391 го ачало с Методрел акса граничнык услйцй- =адчд Решение разностнаго урадненця Пуассона дг ч дгчналога етоЮелаисации слодцй Определение значения дцкря на границе Метод релансаии с использоанием регуляри Я ифг = с 1.5( чс,/ч)=и =( - Г реп а М, -Решение разностного аналогау оаЬненця теплопрододности:д . дт д дт д дтГл ) + - СЛу )+ - (лг - ) +Зх дк др ду дг дг1 Т - + ч - = - а (х,у,г) с учетомдТ ддх дуграничных усладио Т =Т Р/,=0 Конец фиг.9 Решениеура днения1 дг)с учетомю =о,азностного аналога двинеиия:Югш 1 дУ ди дч ди( +дуг ду дк дх ду д и и.игМ у+г ф,о+ дг се г етоды разделеия переменных с пользоданоеи дыстого преИразодания урье