शांत रूप में बहुत सारी जानकारी प्राप्त की जा सकती है, मुझे अनुभव हो रहा है कि अगर कोई पहले भी ऐसा कुछ देख चुका है तो मैं सोच रहा था। प्लॉट एड्विसिएशन-डिफ्यूजन समीकरण के लिए प्रारंभिक स्थिति (हरा) दिखाता है, फिर पुनरावृत्ति 200 (नीला) पर समाधान और फिर पुनरावृत्ति 400 (लाल) पर।
उत्तोलन-विसरण समीकरण का हल कुछ पुनरावृत्तियों के बाद फूटता है। Péclet संख्या , और सीएफएल हालत संतुष्ट हो जाता है, सी ≈ 0.0015 , तो समीकरण चाहिए स्थिर हो। मुझे अनुमान है कि मेरे पास संख्यात्मक कोड में एक बग है।
पृष्ठभूमि। विवेकीकरण दोनों प्रसार और प्रसार शर्तों के लिए केंद्रीय अंतर है। मेरा मानना है कि यह संधि का पहला क्रम है और प्रसार का दूसरा क्रम है। मैंने इसे एक परिमित मात्रा दृष्टिकोण (पहली बार) का उपयोग करके लागू किया है जिसमें कोशिका के चेहरे पर गुणांक (वेग और प्रसार गुणांक) मान सेल औसत से रैखिक प्रक्षेप द्वारा पाए जाते हैं। मैं रॉबिन सीमा स्थिति को बाईं और दाईं सतहों पर लागू करता हूं और फ्लक्स को सीमाओं पर शून्य पर सेट करता हूं।
आप अपने संख्यात्मक कोड को कैसे डिबग करते हैं? क्या किसी ने पहले ऐसा कुछ देखा है, जहां एक अच्छी जगह की तलाश शुरू होगी?
अपडेट करें
- यहाँ मेरी व्यक्तिगत "लैब बुक" शैली नोट है, जो परिवेदना-प्रसार समीकरण के लिए एक परिमित मात्रा विधि को लागू करने के लिए है, http://danieljfarrell.github.io/FVM/
- पायथन स्रोत कोड यहाँ उपलब्ध है, http://github.com/danieljfarrell/FVM.git
अपडेट करें
समाधान अधिक सरल नहीं हो सकता है! मैंने केवल प्रसार शब्द पर एक संकेत त्रुटि की है। यह अजीब है, मुझे यकीन है कि मैंने इसे पोस्ट नहीं किया था मुझे यह त्रुटि नहीं मिली होगी! अगर कोई इस बारे में सुझाव देना चाहता है कि वे अपने संख्यात्मक कोड को कैसे डिबग करें, तो मैं अभी भी दिलचस्पी रखता हूं। मेरे पास कोई विधि नहीं है, यह थोड़ा हिट और मिस है, मैं सुराग पाने के लिए सामान रखने की कोशिश करता रहता हूं, लेकिन इस प्रक्रिया में कई हफ्ते लग सकते हैं (कभी-कभी)।
