Соңғы нүкте әдісі - Finite point method - Wikipedia

The ақырғы нүкте әдісі (FPM) Бұл meshfree әдісі шешу үшін дербес дифференциалдық теңдеулер (PDE) нүктелердің шашыранды үлестірімдері бойынша. FPM тоқсаныншы жылдардың ортасында ұсынылды (Oñate, Idelsohn, Zienkiewicz & Taylor, 1996a),[1] (Oñate, Idelsohn, Zienkiewicz, Taylor & Sacco, 1996b)[2] және (Oñate & Idelsohn, 1998a)[3] күрделі геометрия, еркін беттер, қозғалатын шекаралар мен байланысты мәселелерді шешуді жеңілдету мақсатында адаптивті нақтылау. Содан бері FPM айтарлықтай сұйықтық пен қатты механика мәселелерін шешудің қанағаттанарлық дәлдігі мен мүмкіндіктерін көрсетіп, айтарлықтай дамыды.

Тарих

PDE-ге арналған басқа meshfree әдістеріне ұқсас, ақырлы нүктелік әдіс (FPM) шашыранды деректерді орналастыру және интерполяциялау, негізінен ең кіші квадраттар әдістер (WLSQ). Соңғысын ерекше формалары ретінде қарастыруға болады ең кіші квадраттар Ланкастер мен Салкаускас ұсынған әдіс (MLS).[4] WLSQ әдістері meshfree техникасында кеңінен қолданылды, өйткені MLS-дің көп бөлігін сақтауға мүмкіндік береді, бірақ оларды қолдану тиімді және қарапайым. Осы мақсаттарды ескере отырып, FPM дамуына себеп болған көрнекті тергеу басталды (Oñate, Idelsohn & Zienkiewicz, 1995a)[5] және (Тейлор, Зиенкевич, Oñate & Idelsohn, 1995).[6] Ұсынылған әдістеме жергілікті нүктелер бұлтындағы WLSQ жуықтауларымен және нүктелік коллокацияға негізделген теңдеулерді дискреттеу процедурасымен сипатталды (Батина шығармаларының жолында, 1989 ж.).[7] 1992[8]). FPM-дің алғашқы қосымшалары адаптивті сығылатын ағын проблемаларына бағытталған (Fischer, Onate & Idelsohn, 1995;[9] Oñate, Idelsohn & Zienkiewicz, 1995a;[5] Oñate, Idelsohn, Zienkiewicz & Fisher, 1995b[10]). Жергілікті бұлттардың жақындасуына және салмақ өлшеу функцияларына әсерлер сызықтық және квадраттық көпмүшелік негіздердің көмегімен талданды (Фишер, 1996).[11] Конвекция-диффузия және ағынның қысылмайтын мәселелері аясында қосымша зерттеулер FPM-ге неғұрлым берік негіз берді; cf. (Oñate, Idelsohn, Zienkiewicz & Taylor, 1996a)[1] және (Oñate, Idelsohn, Zienkiewicz, Taylor & Sacco, 1996b).[2] Бұл жұмыстар және (Oñate & Idelsohn, 1998)[3] қазіргі кезде қолданылатын негізгі FPM техникасын анықтады.

Санды жуықтау

FPM сандық жуықтау схемасы

FPM-дегі жуықтауды келесідей қорытындылауға болады. Әр ұпай үшін талдау доменінде (жұлдыз нүктесі), жуықталған шешім жергілікті тірек нүктелерінің ішкі жиынын қолдану арқылы салынған , проблемалық доменге жататын (жергілікті нүктелер бұлты ). Жақындау бұлттың белгісіз түйіндік мәндерінің (немесе параметрлерінің) және белгілі бір метрикалық коэффициенттердің сызықтық комбинациясы ретінде есептеледі. Бұлар түйін параметрлері мен жуықталған шешім арасындағы қашықтық LSQ мағынасында минимизацияланатын бұлт деңгейінде WLSQ есебін шешу арқылы алынады. Жақындау метрикалық коэффициенттері белгілі болғаннан кейін, PDE-ді басқаратын мәселе әрқайсысына іріктеліп алынады жұлдыз нүктесі көмегімен коллокация әдісі. Үздіксіз айнымалылар (және олардың туындылары) таңдалған теңдеулерде дискретті жуықталған формалармен ауыстырылады, ал алынған жүйенің шешімі белгісіз түйін мәндерін есептеуге мүмкіндік береді. Демек, есептің басқарушы теңдеулерін қанағаттандыратын жуықталған шешімді алуға болады. FPM-дің жергілікті сипаттамасы әдісті тиімді параллельді схемаларды жүзеге асыруға ыңғайлы ететіндігін атап өту маңызды.

Әдеттегі FPM жуықтауының құрылысы сипатталған (Oñate & Idelsohn, 1998).[3] Жақындау параметрлерін талдауға болады (Ortega, Oñate & Idelsohn, 2007)[12] және толығырақ зерттеу (Ortega, 2014) өткізіледі.[13] Басқа тәсілдер де ұсынылды, мысалы, қараңыз (Boroomand, Tabatabaei and Oñate, 2005).[14] FPM жуықтауының кеңеюі ұсынылған (Boroomand, Najjar & Oñate, 2009).[15]

Қолданбалар

Сұйықтық механикасы

Сұйық ағыны проблемаларына FPM зерттеуі мен қолданылуының алғашқы бағыттары қысқаша келтірілген (Фишер, 1996).[11] Онда конвективті-диффузиялық есептер LSQ және WLSQ полиномдық жуықтамаларын қолдану арқылы зерттелді. Зерттеу барысында нүктелер бұлтының және салмақтау функцияларының жергілікті жуықтау дәлдігіне әсері басты назарға алынды, бұл FPM негізгі мінез-құлқын түсінуге көмектесті. Нәтижелер көрсеткендей, 1D FPM жуықтауы екінші реттік дәл болатын, орталық айырмашылық жуықтамаларымен алынған дискретті туынды формаларға әкеледі. Алайда, дәлдік өлшеу функциясына байланысты симметриялы емес бұлттар үшін бірінші ретті төмендетеді. Жергілікті бұлтқа сәйкес келетін нүктелерді таңдаудың алдын-ала критерийлері минимизация проблемасының нашар күйін жақсарту мақсатында анықталды. Бұл жұмыста қолданылған ағынды шешуші екі сатылы Тейлор-Галеркин схемасына негізделген, жасанды диссипация. Сандық мысалдарда дыбыстық емес, трансоникалық және дыбыстан жоғары екі өлшемді мәселелер қамтылды, бірақ тұтқыр төмен Рейнольдс сынағының жағдайы ұсынылды. Жалпы алғанда, бұл жұмыста алынған нәтижелер қанағаттанарлық болды және LSQ минимизациясына салмақты енгізу жоғары нәтижелерге әкелетінін көрсетті (сызықтық негіз қолданылды).

Осыған ұқсас зерттеу желісінде ақырғы домендегі ағынды теңдестіру тұрғысынан алынған қалдықты тұрақтандыру әдістемесі, ақырлы өсім есебі (FIC) деп аталады (Oñate, 1996,[16] 1998[17]) енгізілді. Нәтижелер жасанды диссипациямен алынған нәтижелермен салыстырылды, бірақ FIC-тегі тұрақтандыру дәйекті түрде енгізілетіндігімен, қараңыз (Oñate, Idelsohn, Zienkiewicz, Taylor & Sacco, 1996b)[2] және (Oñate & Idelsohn, 1998a).[3]

Осы әзірлемелердің арасында нүктені қалыптастыру мәселесі бірінші кезекте қарастырылды (Löhner & Oñate, 1998).[18] Алдыңғы қатарлы техниканың негізінде авторлар торсыз есептеулерге сәйкес нүктелік дискретизацияларды әдеттегі торды құруға қажет әдеттегі сапа тексерулерінен аулақ бола отырып, тиімдірек құруға болатындығын көрсетті. Дәстүрлі мешерлермен салыстырғанда бәсекеге қабілетті ұрпақтың уақытына қол жеткізілді, бұл бірінші рет торсыз әдістер дискретизация проблемаларын жеңілдетуге болатын балама екенін көрсетті.

Қысылмайтын 2D ағындары алғаш рет зерттелген (Oñate, Sacco & Idelsohn, 2000)[19] пайдалану проекциялау әдісі FIC техникасы арқылы тұрақтандырылды. Бұл тәсілге егжей-тегжейлі талдау жасалды (Sacco, 2002).[20] Осы жұмыстағы керемет жетістіктер FPM-ге берік негіз берді; олардың ішінде жергілікті және нормаланған жуықтау негіздерін анықтау, жергілікті Delaunay триангуляциясы негізінде нүктелердің жергілікті бұлттарын құру процедурасы және нәтижеге жақындау сапасын бағалау критерийі. Ұсынылған сандық қосымшалар негізінен екі өлшемді (тұтқыр және инвисцидті) сығылмайтын ағындарға бағытталған, бірақ үш өлшемді қолдану мысалы да ұсынылған.

Лагранж шеңберінде FPM-ді алдын-ала қолдану (Idelsohn, Storti & Oñate, 2001),[21] туралы да айта кеткен жөн. Сығымдалмаған үшін алынған қызықты нәтижелерге қарамастан еркін бет ағындар, бұл зерттеу желісі FPM шеңберінде жалғаспады және кейінгі тұжырымдар тек Эйлерия ағынының сипаттамаларына негізделген.

3D қысылатын ағындарды шешуге FPM-ді алғашқы қолдану пионер жұмысында ұсынылды (Löhner, Sacco, Oñate & Idelsohn, 2002).[22] Онда нүктелердің жергілікті бұлттарын салудың сенімді және жалпы процедурасы (Delaunay техникасы негізінде) және ағын теңдеулерін шешудің қолайлы схемасы жасалды. Ұсынылған шешім схемасында ағынның дискретті туындылары конвективті тұрақтылықты қамтамасыз ететін орталық айырмашылық тәрізді өрнек және жоғары бағыттағы термиялық термин ретінде бұлт нүктелерін байланыстыратын жиектер бойынша жазылады. Ол үшін Рим мен Ван Лер ағынының векторын бөлудің шамамен Риман шешушісі қолданылды. Ұсынылған тәсіл жасанды диссипация әдістеріне қарағанда дәлірек (сонымен бірге қымбатырақ) және қосымша, жергілікті бұлттағы геометриялық өлшемдерді және проблемаларға тәуелді параметрлерді анықтауды қажет етпейді. Теңдеулердің уақыттық интеграциясы Рунге-Кутта әдісі бойынша көп сатылы айқын схема арқылы орындалды.

Бірнеше жылдан кейін, 3D FPM жуықтауларына қатысты одан әрі зерттеулер жүргізілді (Ortega, Oñate & Idelsohn, 2007).[12] Бұл жұмыс жергілікті қолдаудың сипаттамаларына қарамастан сенімді жуықтама құруға бағытталды. Осы мақсатта салмақ өлшеу функциясын және басқа жуықтау параметрлерін жергілікті автоматты түрде реттеу ұсынылды. Сығымдалатын аэродинамикамен байланысты әдісті бұдан әрі 3D қолдану адаптивті нақтылау (Ortega, Oñate & Idelsohn, 2009)[23] және қозғалатын / деформацияланатын шекара мәселелер (Ortega, Oñate & Idelsohn, 2013).[24] Бұл жұмыстарда FPM қанағаттанарлық беріктік пен дәлдікті және практикалық есептеулерді шешу мүмкіндіктерін көрсетті. Басқа жетістіктермен қатар, модельдік дискреттеудің толық қалпына келуі, тіпті үлкен модельдеу проблемаларында да қол жетімді шешім стратегиясы бола алатындығы көрсетілді. Бұл нәтиже қозғалатын / деформацияланатын домен проблемаларын торсыз талдаудың жаңа мүмкіндіктерін ұсынады. FPM адаптацияға сәтті қолданылды таяз су проблемалары (Ortega, Oñate, Idelsohn & Buachart, 2011)[25] және (Buachart, Kanok-Nukulchai, Ortega & Oñate, 2014).[26] Рейнольдстың тұтқыр ағыны проблемаларында торсыз артықшылықтарды пайдалану туралы ұсыныс келтірілген (Ortega, Oñate, Idelsohn & Flores, 2014a).[27]

Сол қосымшалар саласында FPM дәлдігі, есептеу құны және параллель өнімділігі бойынша үлкен зерттеу жүргізілді (Ortega, Oñate, Idelsohn & Flores, 2014b).[28] Онда FPM эквивалентті Finite Element негізіндегі еріткішпен салыстырылды, бұл екеуін де, торсыз еріткіштің сипаттамаларын және оның практикалық қолданбаларды шешуге жарамдылығын бағалауға стандарт берді. Бұл жұмыста FPM техникасының кейбір жеңілдетулері тиімділікті арттыру және ФЭМ-мен жұмыс алшақтықты азайту үшін ұсынылды. Содан кейін қанат корпусының конфигурациясын қолдана отырып, тордың конвергенциясы бойынша зерттеулер жүргізілді. Нәтижелер салыстырмалы дәлдік пен өнімділікті көрсетіп, FPM-ді өзінің әріптесіне қатысты бәсекеге қабілетті екенін көрсетті. Бұл өте маңызды, өйткені торсыз техникалар бастапқы іске асырудың тиімділігінің төмендігіне байланысты практикалық емес болып саналады.

FPM қолданылды аэроакустика жылы (Байко, Чермак және Джича, 2014).[29] Ұсынылған шешім схемасы сызықты Риман еріткішіне негізделген және жоғары деңгейлі FPM жуықтамаларының артықшылықтарын сәтті пайдаланады. Алынған нәтижелер FPM-дің дыбысты тарату мәселелерін шешуге арналған әлеуетін көрсетеді.

Қатты механика

Ағымдағы тергеу бағыттары

Ағымдағы күштер негізінен FPM-дің ауқымды практикалық мәселелерді шешу үшін параллельді ортада жұмыс істеу мүмкіндіктерін пайдалануға бағытталған, әсіресе торсыз процедуралар пайдалы үлес қосуы мүмкін жерлерде, мысалы, күрделі геометрия, қозғалатын / деформацияланатын домен, адаптивті нақтылау сияқты мәселелер және көпөлшемді құбылыстар.

Әдебиеттер тізімі

  1. ^ а б Онате, Э .; Иделсон, С .; Зиенкевич, О. С .; Тейлор, Р.Л (1996). «Сұйықтық механикасы мәселелерін талдаудың ақырғы нүктелік әдісі. Конвективті тасымалдауға және сұйықтық ағынына қосымшалар» Инженериядағы сандық әдістерге арналған халықаралық журнал. 39 (2): 3839–3866. Бибкод:1996IJNME..39.3839O. дои:10.1002 / (SICI) 1097-0207 (19961130) 39:22 <3839 :: AID-NME27> 3.0.CO; 2-R.
  2. ^ а б c Oñate, E .; Иделсон, С .; Зиенкевич, О. С .; Тейлор, Р.Л .; Sacco, C. (1996). «Сұйықтық механикасы мәселелерін талдауға арналған тұрақтандырылған ақырғы нүкте әдісі». Қолданбалы механика мен техникадағы компьютерлік әдістер. 139 (1): 315–346. Бибкод:1996CMAME.139..315O. дои:10.1016 / s0045-7825 (96) 01088-2.
  3. ^ а б c г. Oñate, E .; Иделсон, С. (1998). «Адвективті-диффузиялық тасымалдау және сұйықтық ағыны проблемалары үшін торсыз ақырлы нүктелік әдіс». Есептеу механикасы. 24 (4–5): 283–292. Бибкод:1998CompM..21..283O. дои:10.1007 / s004660050304.
  4. ^ Ланкастер, П .; Салкаускас, К. (1981). «Ең кіші квадраттардың қозғалуы нәтижесінде пайда болған беттер». Есептеу математикасы. 37 (155): 141–158. дои:10.2307/2007507. JSTOR  2007507.
  5. ^ а б Oñate, E .; Иделсон, С .; Zienkiewicz, O. C. (1995). «Есептеу механикасындағы ақырғы нүкте әдістері». CIMNE басылымы Nº 74: Инженериядағы сандық әдістердің халықаралық орталығы.
  6. ^ Тейлор, Р.Л .; Зиенкевич, О. С .; Oñate, E .; Иделсон, С. (1995). «Дифференциалдық теңдеулерді шешуге арналған ең кіші квадраттық жуықтамаларды жылжыту». CIMNE басылымы Nº 74 (31-бет): Инженериядағы сандық әдістердің халықаралық орталығы.
  7. ^ Батина, Дж. Т. (1989). «Компьютерлік аэроэластикалық талдауға арналған құрылымсыз динамикалық тормен тұрақсыз Эйлер алгоритмі». AIAA қағазы. 89: 1189.
  8. ^ Батина, Дж. Т. (1992). «Күрделі екі өлшемді қосымшалар үшін торсыз Эйлер / Навье-Стокс шешімінің алгоритмі». Nasa-Tm-107631.
  9. ^ Фишер, Т .; Oñate, E .; Иделсон, С. (1995). «Жоғары жылдамдықты ағындарды компьютерлік талдауға арналған торсыз әдіс». AGARD CFD әдістері мен алгоритмдеріндегі прогресс және қиындықтар туралы симпозиумда ұсынылған жұмыс, Севилья.
  10. ^ Oñate, E .; Иделсон, С .; Зиенкевич, О. С .; Фишер, Т. (1995). «Есептеу механикасындағы ақырғы нүкте әдістері». Сұйықтардағы ақырғы элементтердің әдістері бойынша конференция, Венизе, Италия, 15-21.
  11. ^ а б Фишер, Т. (1996). «Сығылатын ағынды есептердің адаптивті сандық шешіміне үлес». Докторлық диссертация, Universitat Politècnica de Catalunya.
  12. ^ а б Ортега, Е .; Oñate, E .; Иделсон, С. (2007). «Үш өлшемді потенциал ағындарының жетілдірілген ақырғы нүктелік әдісі». Есептеу механикасы. 40 (6): 949–963. Бибкод:2007CompM..40..949O. дои:10.1007 / s00466-006-0154-6.
  13. ^ Ортега, Е .; Oñate, E .; Иделсон, С. (2014). Сығылатын аэродинамикалық есептерге ақырғы нүктелік әдісті құру және қолдану (PDF). CIMNE монографиясы M143. ISBN  978-84-941686-7-3.
  14. ^ Боруманд, Б .; Табатабаей, А.А .; Oñate, E. (2005). «Ақырғы нүктелік әдісті тұрақтандыруға арналған қарапайым модификация». Инженериядағы сандық әдістерге арналған халықаралық журнал. 63 (3): 351–379. Бибкод:2005IJNME..63..351B. дои:10.1002 / nme.1278.
  15. ^ Боруманд, Б .; Наджар, М .; Oñate, E. (2009). «Жалпыланған ақырғы нүкте әдісі». Есептеу механикасы. 44 (2): 173–190. Бибкод:2009CompM..44..173B. дои:10.1007 / s00466-009-0363-x.
  16. ^ Oñate, E. (1996). «Конвективті тасымал және сұйықтық ағыны мәселелерінің сандық шешімін тұрақтандыру туралы». Nº 81 зерттеу есебі: Инженериядағы сандық әдістердің халықаралық орталығы.
  17. ^ Oñate, E. (1998). «Адвективті-диффузиялық тасымалдау және сұйықтық ағыны мәселелерін сандық шешуге арналған тұрақтандырылған теңдеулер шығару». Қолданбалы механика мен техникадағы компьютерлік әдістер. 151 (1): 233–265. Бибкод:1998CMAME.151..233O. дои:10.1016 / s0045-7825 (97) 00119-9.
  18. ^ Лёнер, Р .; Oñate, E. (1998). «Алдыңғы қатарды қалыптастыру техникасы». Техникадағы сандық әдістердегі байланыс. 14 (12): 1097–1108. дои:10.1002 / (sici) 1099-0887 (199812) 14:12 <1097 :: aid-cnm183> 3.0.co; 2-7.
  19. ^ Oñate, E .; Сакко, С .; Иделсон, С. (2000). «Ағынның қысылмайтын мәселелеріне арналған ақырғы нүктелік әдіс». Ғылымдағы есептеу және көрнекілік. 3 (1–2): 67–75. дои:10.1007 / s007910050053.
  20. ^ Sacco, C. (2002). «Desarrollo del método de puntos finitos en mecánica de fluidos». Докторлық диссертация, Universitat Politècnica de Catalunya.
  21. ^ Иделсон, С .; Сторти, М .; Oñate, E. (2001). «Сұйық ағындардың еркін бетіндегі қысылмайтын сұйықтық ағындарын шешуге арналған лагранжды формулалар». Қолданбалы механика мен техникадағы компьютерлік әдістер. 191 (6): 583–593. Бибкод:2001CMAME.191..583R. дои:10.1016 / s0045-7825 (01) 00303-6.
  22. ^ Лёнер, Р .; Сакко, С .; Oñate, E .; Иделсон, С. (2002). «Қысылатын ағынның ақырғы нүктелік әдісі». Инженериядағы сандық әдістерге арналған халықаралық журнал. 53 (8): 1765–1779. Бибкод:2002IJNME..53.1765L. дои:10.1002 / nme.334. hdl:2117/167123.
  23. ^ Ортега, Е .; Онате, Э .; Иделсон, С. (2009). «Үш өлшемді қысылатын ағынды есептеудің ақырғы нүктелік әдісі». Сұйықтықтағы сандық әдістерге арналған халықаралық журнал. 60 (9): 937–971. Бибкод:2009IJNMF..60..937O. дои:10.1002 / fld.1892. hdl:2117/24488.
  24. ^ Ортега, Е .; Онате, Э .; Иделсон, С .; Флорес, Р. (2013). «Қозғалмалы шекаралар мен адаптацияға байланысты қысылатын ағын мәселелерін үш өлшемді талдауға арналған торсыз ақырлы нүктелік әдіс». Сұйықтықтағы сандық әдістерге арналған халықаралық журнал. 73 (4): 323–343. Бибкод:2013IJNMF..73..323O. дои:10.1002 / fld.3799. hdl:2117/86276.
  25. ^ Ортега, Е .; Oñate, E .; Иделсон, С .; Buachart, C. (2011). «Таяз су теңдеулеріне арналған адаптивті ақырлы нүктелік әдіс». Инженериядағы сандық әдістерге арналған халықаралық журнал. 88 (2): 180–204. Бибкод:2011IJNME..88..180O. дои:10.1002 / nme.3171.
  26. ^ Буахарт, С .; Канок-Нукулчай, В .; Ортега, Е .; Oñate, E. (2014). «Шекті нүкте әдісімен таяз су моделі». Халықаралық есептеу әдістері журналы. 11 (1): 1350047. дои:10.1142 / S0219876213500473.
  27. ^ Ортега, Е .; Oñate, E .; Иделсон, С .; Флорес, Р. (2014). «Рейнольдс саны жоғары қысылатын ағынды есептерге ақырғы нүктелік әдісті қолдану». Сұйықтықтағы сандық әдістерге арналған халықаралық журнал. 74 (10): 732. Бибкод:2014IJNMF..74..732O. дои:10.1002 / fld.3871.
  28. ^ Ортега, Е .; Oñate, E .; Иделсон, С .; Флорес, Р. (2014). «Сығымдалатын ағынды есептердегі ақырғы нүктелік әдістің салыстырмалы дәлдігі мен өнімділігін бағалау». Компьютерлер және сұйықтықтар. 89: 53–65. дои:10.1016 / j.compfluid.2013.10.024.
  29. ^ Байко, Дж .; Кермак, Л .; Jícha, M. (2014). «Дыбысты тарату мәселелерін шешудің жоғары ретті ақырлы нүктелік әдісі». Қолданбалы механика мен техникадағы компьютерлік әдістер. 280: 157–175. Бибкод:2014CMAME.280..157B. дои:10.1016 / j.cma.2014.07.022.