حل معادله غير خطي ريچاردز با تركيب الگوريتم¬ فوق ابتكاري الكترومغناطيس و الگوريتم نيوتن

پيش بيني حركت سيال در خاك¬هاي غير اشباع، يكي از مسايل مهم در بسياري از شاخه¬هاي علوم مهندسي است. بطوريكه تقريبا در همه مطالعات انجام شده فرض مي¬شود، حركت سيال در ناحيه¬ي غير اشباع از معادلات ريچاردز پيروي مي¬كند. اين معادله براي پيش¬بيني و شبيه¬سازي حركت سيال در محيط غيراشباع، استفاده مي¬شود. يكي از روش¬هاي حل معادله ريچاردز، استفاده از غالب حجم محدود براي گسسته¬سازي، سپس حل معادلات حاصل، توسط روش¬هاي تكراري نيوتن يا پيكارد است. اگرچه روش نيوتن روشي كارآمد، آسان و سريع براي حل معادلات غير خطي است، اما همگرايي آن وابسته به حدس اوليه مناسب است. از اينرو در اين پژوهش سعي گرديد، با معرفي يك الگوريتم ابتكاري تركيبي، حدس اوليه مناسب براي پتانسيل فشار خاك توليد نموده، سپس با استفاده از اين حدس اوليه همگرايي روش نيوتن تضمين گردد. براي صحت¬سنجي جواب¬هاي حاصل از اين الگوريتم، معادله يك بعدي شده ريچاردز با حدس¬هاي اوليه گوناگوني بوسيله الگوريتم نيوتن حل شد و مشاهده گرديد، حدس¬هاي نامناسب در الگوريتم نيوتن با بيش از 2000تكرار بازهم به همگرايي نرسيد درحاليكه با هر حدس اوليه دلخواه جواب هاي حاصل از الگوريتم تركيبي ارائه شده براي مسايل مختلف همگرا شدند و با حل تحليلي مقايسه شدند كه مطابقت بسيار عالي داشتند.