'

“ Моделирование кодом FiDAP свободной поверхности плоской турбулентной струи, формируемой соплом SHIMA”.

Понравилась презентация – покажи это...





Слайд 0

“ Моделирование кодом FiDAP свободной поверхности плоской турбулентной струи, формируемой соплом SHIMA”. ______________________________________________


Слайд 1

______________________________________________ ЗАЧЕМ ЭТО НАДО ? Цель расчетного исследования заключалась: в моделировании контурной гидродинамики литиевой петли в исследовании устойчивости свободной поверхности изотермической плоской струи при истечении ее из сопла SHIMA в вакуум [4?15] в начальном этапе оптимизации геометрии кромки сопла и выяснении причин образования аэрозолей


Слайд 2

______________________________________________ ЧТО СДЕЛАНО ? С помощью лицензионного кода FiDAP v.8.7.2 [1?3] выполнено 2D расчетное моделирование части петли литиевого стенда. Исследован стационарный изотермический, турбулентный режим в условиях полного смачивания и отсутствия шероховатости конструкций. С помощью лицензионного пре-процессора GAMBIT cмоделированы расчетные сетки и исследованы 3 варианта угловых конфигураций кромок сопла SHIMA на четырехугольной криволинейной сетке: 22?, 62?30? (проектное значение) и 67?30?. Количество ячеек расчетной сетки варьировалась в пределах 108284 (для проектного варианта) до 92342(для вариантов: 22? и 67?30?). При этом количество ячеек для области струи ( по глубине) для проектного варианта составило 12 (размер ячейки 0.83 мм) и 4 (размер ячейки 2.5 мм) для остальных вариантов.


Слайд 3

______________________________________________ ЧТО ПРЕДСТАВЛЯЕТ СОБОЙ ФРАГМЕНТ ПЕТЛИ ? В состав фрагмента петли вошли: камера -успокоитель (на входе) без решеток модельное сопло SHIMA с габаритами: 70x97 мм сечение на входе в первую ступень 70x10 мм сечение на выходе второй ступени 177 мм –общая длина . прямолинейный 90 мм участок, наклоненный под углом 67?30' к горизонту. криволинейный участок (R=250 мм, с угловым растром в 45?) прямолинейный 198 мм участок на сливе приемный успокоительный бак на сливе.


Слайд 4

______________________________________________ РАСЧЕТНАЯ ОБЛАСТЬ


Слайд 5

______________________________________________ ЧТО ОЗНАЧАЕТ ТЕРМИН “ FIDAP ” ?


Слайд 6

______________________________________________ МАТЕМАТИЧЕСКАЯ МОДЕЛЬ V.O.F. для FiDAP уравнение для турбулентной энергии уравнение для диссипации турбулентной энергии: уравнение метода (Volume Of Fluid) кода FIDAP ([1] стр.13-2): Система уравнений расчетной модели решена с помощью штатного алгоритма кода : ‘SEGREGATED FREE-SURFACE ALGORITM’.


Слайд 7

______________________________________________ СТРУКТУРА РЕШАТЕЛЯ FiDAP


Слайд 8

______________________________________________ ВЗАИМОДЕЙСТВИЕ FiDAP С ПРЕ-ПРОЦЕССОРОМ GAMBIT


Слайд 9

______________________________________________ ЧТО ОЗНАЧАЕТ ТЕРМИН ”GAMBIT” ?


Слайд 10

______________________________________________ ГРАНИЧНЫЕ УСЛОВИЯ


Слайд 11

______________________________________________ РАСЧЕТНЫЕ РЕЗУЛЬТАТЫ ДЛЯ варианта кромки сопла 62?30? (проектное значение).


Слайд 12

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). 2D расчетная урезанная расчетная сетка в области второй ступени сопла, кромки и прямолинейного 90 мм участка. УВЕЛИЧЕНИЕ


Слайд 13

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). 2D расчетная урезанная расчетная сетка в области первой и второй ступеней сопла, кромки и прямолинейного 90 мм участка. УВЕЛИЧЕНИЕ


Слайд 14

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). 2D расчетная урезанная расчетная сетка в области действия пучка дейтронов. УВЕЛИЧЕНИЕ


Слайд 15

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). 2D векторное поле скорости в области второй ступени сопла, кромки и прямолинейного 90 мм участка. УВЕЛИЧЕНИЕ


Слайд 16

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). 2D векторное поле скорости в области первой и второй ступеней сопла, кромки и прямолинейного 90 мм участка. УВЕЛИЧЕНИЕ


Слайд 17

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). 2D поле вертикальной составляющей скорости в области второй ступени сопла, кромки и прямолинейного 90 мм участка. УВЕЛИЧЕНИЕ


Слайд 18

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). 2D поле вертикальной составляющей скорости внутри сопла, кромки и прямолинейного 90 мм участка


Слайд 19

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). 2D векторное поле в области кромки (с наложением на расчетную сетку). УВЕЛИЧЕНИЕ. (построено с помощью пост-процессора FIELDVIEW v.8.1.f [3]).


Слайд 20

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). 2D векторное поле в области кромки (с наложением на расчетную сетку). УВЕЛИЧЕНИЕ. (построено с помощью пост-процессора FIELDVIEW v.8.1.f [3]).


Слайд 21

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). 2D поле турбулентного числа Рейнольдса в области второй ступени, кромки и прямолинейного участка. УВЕЛИЧЕНИЕ.


Слайд 22

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). 2D поле сдвига в области второй ступени, кромки и прямолинейного участка. УВЕЛИЧЕНИЕ.


Слайд 23

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). 2D поле турбулентной вязкости в области второй ступени, кромки и прямолинейного участка. УВЕЛИЧЕНИЕ.


Слайд 24

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). 2D поле завихренности в области второй ступени, кромки и прямолинейного участка. УВЕЛИЧЕНИЕ.


Слайд 25

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). 2D поле турбулентной кинетической энергии в области второй ступени, кромки и прямолинейного участка. УВЕЛИЧЕНИЕ.


Слайд 26

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). 1D распределение давления (Pa) в области сопряжения прямолинейного наклонного участка и и криволинейной части ( угол 22?30? от оси симметрии пучка). Отсчет от свободной поверхности. Примечание : вследствие гидродинамического прыжка с повышением уровня толщина струи в сечении составляет 11.0 мм


Слайд 27

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). 1D распределение давления (Pa) в области оси симметрии пучка ( угол 0?). Отсчет от свободной поверхности. Примечание : вследствие гидродинамического прыжка с повышением уровня толщина струи в сечении составляет 11.5 мм


Слайд 28

______________________________________________ Вариант кромки сопла 62?30? (проектное значение)


Слайд 29

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Векторное поле и функция FILLING для сечения струи лития


Слайд 30

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Векторное поле и функция FILLING для сечения струи лития


Слайд 31

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Векторное поле и функция FILLING для сечения струи лития


Слайд 32

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Векторное поле и функция FILLING для сечения струи лития


Слайд 33

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Векторное поле и функция FILLING для сечения струи лития


Слайд 34

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Векторное поле и функция FILLING для сечения струи лития


Слайд 35

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Векторное поле на выходе из второй ступени сопла SHIMA


Слайд 36

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Векторное поле на выходе из второй ступени сопла SHIMA


Слайд 37

______________________________________________ Вариант кромки сопла 62?30? (проектное значение)


Слайд 38

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Вертикальная компонента поля скорости.


Слайд 39

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Вертикальная компонента поля скорости.


Слайд 40

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Горизонтальная компонента поля скорости.


Слайд 41

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Поле скорости.


Слайд 42

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Поле турбулентной кинетической энергии.


Слайд 43

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Распределение диссипации турбулентной энергии.


Слайд 44

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Распределение турбулентной вязкости.


Слайд 45

______________________________________________ Вариант кромки сопла 62?30? (проектное значение).Распределение поля сдвига.


Слайд 46

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Распределение числа Рейнольдса


Слайд 47

______________________________________________ РАСЧЕТНЫЕ РЕЗУЛЬТАТЫ ДЛЯ варианта кромки сопла 22? ( НЕ проектное значение).


Слайд 48

______________________________________________ Вариант кромки сопла 22? ( НЕ проектное значение). 2D расчетная сетка в области первой –второй ступеней, кромки и прямолинейного 90 мм участка. УВЕЛИЧЕНИЕ


Слайд 49

______________________________________________ Вариант кромки сопла 22? ( НЕ проектное значение). 2D векторное поле скорости в области первой –второй ступеней, кромки и прямолинейного 90 мм участка. УВЕЛИЧЕНИЕ


Слайд 50

______________________________________________ Вариант кромки сопла 22? ( НЕ проектное значение). 2D поле скорости в области первой –второй ступеней, кромки и прямолинейного 90 мм участка. УВЕЛИЧЕНИЕ


Слайд 51

______________________________________________ Вариант кромки сопла 22? ( НЕ проектное значение). 2D поле давления в области первой –второй ступеней, кромки и прямолинейного 90 мм участка. УВЕЛИЧЕНИЕ


Слайд 52

______________________________________________ Вариант кромки сопла 22? ( НЕ проектное значение). 2D поле вертикальной составляющей скорости в области первой –второй ступеней, кромки и прямолинейного 90 мм участка. УВЕЛИЧЕНИЕ


Слайд 53

______________________________________________ Вариант кромки сопла 22? ( НЕ проектное значение). 2D поле горизонтальной составляющей скорости в области первой –второй ступеней, кромки и прямолинейного 90 мм участка. УВЕЛИЧЕНИЕ


Слайд 54

______________________________________________ Вариант кромки сопла 22? ( НЕ проектное значение). 2D поле завихренности в области кромки и прямолинейного 90 мм участка. УВЕЛИЧЕНИЕ


Слайд 55

______________________________________________ Вариант кромки сопла 22? ( НЕ проектное значение). 2D поле турбулентной вязкости в области кромки и прямолинейного 90 мм участка. УВЕЛИЧЕНИЕ


Слайд 56

______________________________________________ Вариант кромки сопла 22? ( НЕ проектное значение). 2D поле в области кромки и прямолинейного 90 мм участка. УВЕЛИЧЕНИЕ


Слайд 57

______________________________________________ Вариант кромки сопла 22? ( НЕ проектное значение). 2D поле турбулентной кинетической энергии в области кромки и прямолинейного 90 мм участка. УВЕЛИЧЕНИЕ


Слайд 58

______________________________________________ Вариант кромки сопла 22? ( НЕ проектное значение). 2D поле диссипации турбулентной энергии в области кромки и прямолинейного 90 мм участка. УВЕЛИЧЕНИЕ


Слайд 59

______________________________________________ МУЛЬТИПЛИКАЦИЯ . Вариант кромки сопла 62?30? (проектное значение).


Слайд 60

______________________________________________ МУЛЬТИПЛИКАЦИЯ . Вариант кромки сопла 62?30? (проектное значение). Грубая сетка.


Слайд 61

______________________________________________ МУЛЬТИПЛИКАЦИЯ . Вариант кромки сопла 62?30? (проектное значение). Грубая сетка.


Слайд 62

______________________________________________ МУЛЬТИПЛИКАЦИЯ . Вариант кромки сопла 62?30? (проектное значение). Грубая сетка.


Слайд 63

______________________________________________ МУЛЬТИПЛИКАЦИЯ . Вариант кромки сопла 62?30? (проектное значение). Урезанная (подробная) сетка сетка.


Слайд 64

______________________________________________ РЕЗУЛЬТАТЫ Выявлена причина образования аэрозолей, которые будут присутствовать в рабочей области пучка всегда и , возможно, будут оседать на стенках конструкций и уноситься в тракт укорителя. 1.1. Возможно лишь снизить интенсивность их возникновения (если это позволяет технология и требования к эрозионному износу кромок). 1.2. В области краевого угла второй ступени образуется зона максимальной турбулентности, которая образует турбулентный след линейным размером порядка 10 мм, в котором происходит зарождений аэрозолей. 1.3. Прослеживается закономерность размера краевого угла ? максимума турбулентной кинетической энергии? его расположения относительно угла ? интенсивности возникновения аэрозолей: 1.3.1. C уменьшением краевого угла со значения 62?30? до 22? максимальное значение кинетической турбулентной энергии (?) увеличивается со значения 54.17м2/c3 до 56.04 м2/c3. Это сопровождается смещением максимума турбулентности от краевого угла (приблизительно с 2 мм до 8 мм) вниз по потоку и снижением интенсивности возникновения аэрозолей! Однако полного устранения аэрозолей, по всей видимости, достичь невозможно из-за технологических требований к надежности конструкции сопла. 2. Получена исчерпывающая информация гидродинамики сопла SHIMA для штатного расположения 67?30? к горизонту в виде двумерных распределений полей (и компонентов) скорости, давления и турбулентности, как для внутренней части (1 и 2 ступень), так для выхода струи в вакуум на прямолинейный 90 мм участок. И ЧТО ДАЛЬШЕ ? (смотри далее…..)


Слайд 65

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Пробный вариант сетки . 40 ячеек по глубине струи. Размер ячейки по высоте 0.25 мм. На область аэрозолей отведено дополнительно 14 ячеек с таким же размером по высоте.


Слайд 66

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Пробный вариант сетки . 40 ячеек по глубине струи. Размер ячейки по высоте 0.25 мм. На область аэрозолей отведено дополнительно 14 ячеек с таким же размером по высоте. ОБЛАСТЬ СОПРЯЖЕНИЯ.


Слайд 67

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Пробный вариант сетки . 40 ячеек по глубине струи. Размер ячейки по высоте 0.25 мм. На область аэрозолей отведено дополнительно 14 ячеек с таким же размером по высоте. ОБЛАСТЬ ПУЧКА. ПОВЕРНУТО.


Слайд 68

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Пробный вариант сетки . 40 ячеек по глубине струи. Размер ячейки по высоте 0.25 мм. На область аэрозолей отведено дополнительно 14 ячеек с таким же размером по высоте. ОБЛАСТЬ СЛИВА. ПОВЕРНУТО.


Слайд 69

______________________________________________ International Association for Hydraulic Engineering and Research (IAHR) / Section on Energy Exchange and Fluid Phenomena Group Meeting July 5 - 10 , 2004, Obninsk, Russia Вариант кромки сопла 62?30? (проектное значение). Пробный вариант сетки . ОБЛАСТЬ ВТОРОЙ СТУПЕНИ. ПОВЕРНУТО.


Слайд 70

______________________________________________ International Association for Hydraulic Engineering and Research (IAHR) / Section on Energy Exchange and Fluid Phenomena Group Meeting July 5 - 10 , 2004, Obninsk, Russia Вариант кромки сопла 62?30? (проектное значение). Пробный вариант сетки . ОБЛАСТЬ ПЕРВОЙ И ВТОРОЙ СТУПЕНЕЙ. ПОВЕРНУТО.


Слайд 71

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Пробный вариант сетки . ОБЛАСТЬ ПЕРВОЙ , ВТОРОЙ СТУПЕНЕЙ И ПОРИСТОЙ ОБЛАСТИ. ПОВЕРНУТО.


Слайд 72

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Вариант кромки сопла 62?30? (проектное значение). Пробный вариант сетки . ОБЛАСТЬ ВТОРОЙ СТУПЕНИ И ПОРИСТОЙ ОБЛАСТИ. ПОВЕРНУТО.


Слайд 73

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Пробный вариант сетки . ОБЛАСТЬ ПЕРВОЙ , ВТОРОЙ СТУПЕНЕЙ И ПОРИСТОЙ ОБЛАСТИ. ПОВЕРНУТО.


Слайд 74

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Пробный вариант сетки .


Слайд 75

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). Граничные условия в среде кода FIDAP.


Слайд 76

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). ROBAST- вариант сетки . 40 ячеек по глубине струи. Размер ячейки по высоте 0.25 мм. На область аэрозолей отведено дополнительно 14 ячеек с таким же размером по высоте. ОБЛАСТЬ СОПРЯЖЕНИЯ.


Слайд 77

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). ROBAST- вариант сетки . 40 ячеек по глубине струи. Размер ячейки по высоте 0.25 мм. На область аэрозолей отведено дополнительно 14 ячеек с таким же размером по высоте. ВТОРАЯ СТУПЕНЬ СОПЛА.


Слайд 78

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). ROBAST- вариант сетки . 40 ячеек по глубине струи. Размер ячейки по высоте 0.25 мм. На область аэрозолей отведено дополнительно 14 ячеек с таким же размером по высоте. 1и 2 СТУПЕНИ СОПЛА.


Слайд 79

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). ROBAST- вариант сетки . 40 ячеек по глубине струи. Размер ячейки по высоте 0.25 мм. На область аэрозолей отведено дополнительно 14 ячеек с таким же размером по высоте. 1 СТУПЕНЬ И ПОРИСТАЯ СТРУКТУРА


Слайд 80

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). ROBAST- вариант сетки . 40 ячеек по глубине струи. Размер ячейки по высоте 0.25 мм. На область аэрозолей отведено дополнительно 14 ячеек с таким же размером по высоте.


Слайд 81

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). ROBAST- вариант сетки . 40 ячеек по глубине струи. Размер ячейки по высоте 0.25 мм. На область аэрозолей отведено дополнительно 14 ячеек с таким же размером по высоте. ОБЛАСТЬ ПУЧКА


Слайд 82

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). ROBAST- вариант сетки . 40 ячеек по глубине струи. Размер ячейки по высоте 0.25 мм. На область аэрозолей отведено дополнительно 14 ячеек с таким же размером по высоте. ОБЛАСТЬ СЛИВА


Слайд 83

______________________________________________ Вариант кромки сопла 62?30? (проектное значение). ROBAST- вариант сетки . 40 ячеек по глубине струи. Размер ячейки по высоте 0.25 мм. На область аэрозолей отведено дополнительно 14 ячеек с таким же размером по высоте. ОБЛАСТЬ СЛИВА В БАК.


Слайд 84

______________________________________________ 1.FIDAP v. 8.7.2. Theory Manual. ( 1998. Fluent Inc., USA, Lebanon, NH, 03766) pp. 7-24, 10-14, 13-2?13-3,10.13?10.15 2.Gambit v.2.1 Theory Manual. ( 2003. Fluent Inc., USA, Lebanon, NH,03766) 3. FIELDVIEW v.8.1.f . User Gude. ( 2002. Intelligent Light Inc., USA, Lyndhurst, NJ,07071) 4. A.Shima, Theory of Direct and Inverse Methods to obtain Nozzle Shape, Mem.Inst.High Sp. Mech., Japan. Vol.17,№164,(1961/1962) 61-86. 5. M.Ida et al. Thermal-Hydraulic Characteristics of IFMIF Liquid Lithium Target. ISNT-6. April 8-12,2002,San Diego,USA. ISTC Project 2036. “Lithium Circulation Loop and Lithium Neutron Target” http://www.tech-db.ru/istc/db/projects.nsf/wu/2036 7. Hiroo Nakamura et al. Status of Lithium Target System for International Fusion Materials Irradiation Facility (IFMIF) Fusion Engineereing and Design 58-59 (2001) 919-923. 8. Hideo Nakamura et al. Water Experiment of High-speed, Free-surface, Plane Jet along Cocave wall , Proceeding of the Eighth International Topical Meeting on Nuclear Reactor Thermal-Hydraulics (NURETH-8),Kyoto,Japan,September 30-October 4,(1997), 1268-1275. 9. Hiroshi Nakamura et al. Experimental Study on the Effect of Nozzle surface Finishing in High-speed water Jet,JAERI Report,JAERI-Research 2000-068, Feb.2001 (in Japanese). 10. M.Ida et al. Water Jet Flow Simulation Experiment and Plan of Lithium free Surface Flow Experiment for IFMIF Target, Proceeding of the 10 th International Conference on Fusion reactor Materials (ICFRM 10), Baden- baden,Germany,October,14-19, 2001. 11. M.Ida et al.,Fluid Stability Analysis for IFMIF Target Proceeding of the International Topical Meeting on Nuclear Applications of Accelerator technolology (AccApp’98), Catlinburg, Tennessee, USA, September 22-23 (1998),546-555. 12. IFMIF Team (edited by A.Moslang ), IFMIF International Fusion Materials Irradiation Facility Conceptual Design Evaluation Report, A Supplement to the CDA by the IFMIF Team ,FZK Report, FZKA6199, January 1999. 13. M.Ida et al., “Thermal and Fluid Analysis for the IFMIF Lithium Target Jet Flow,”Proc. Of NURETH-8,Sept.30- Oct.4,1997,p.1276.Atomic Energy Society of Japan, (1997) 14. M.Ida et al., JAERI-Research98-0022, (1997) 15. S.Celovani, “Thermal Hydraulics of Liquid Metal Jets”,Proc.NURETH-8,Sept.30- Oct.4,1997,p.1276,Atomic Energy Socity of Japan,(1997) Литература


×

HTML:





Ссылка: