پایتون 88250 /

---

شرح پروژه

#کد_پروژه: 88250

موضوع: پروژه پایتون /

import numpy as np

from scipy.optimize import newton

import matplotlib.pyplot as plt

# تعریف معادله کولبروک

def colebrook(f, re, epsilon):

    return 1 / np.sqrt(f) + 2.0 * np.log10(epsilon / 3.7 + 2.51 / (re * np.sqrt(f)))

# پارامترها

reynolds_numbers = np.logspace(3.6, 8, num=100)  # محدوده رینولدز (log space)

roughness_values = [0.0, 0.00000, 0.0008, 0.02]  # مقادیر زبری

initial_guess = 0.02  # حدس اولیه برای حل نیوتن

# محاسبه ضریب اصطکاک برای هر ترکیب از Re و ε

friction_factors = {}

for epsilon in roughness_values:

    friction_factors[epsilon] = []

    for re in reynolds_numbers:

        if re < 4000:  # شرط جریان متلاطم

            friction_factors[epsilon].append(np.nan)

            continue

        try:

            # حل معادله با روش نیوتن

            f_solution = newton(colebrook, initial_guess, args=(re, epsilon))

            friction_factors[epsilon].append(f_solution)

        except RuntimeError:

            friction_factors[epsilon].append(np.nan)

# رسم نمودار

plt.figure(figsize=(12, 8))

for epsilon in roughness_values:

    plt.plot(

        reynolds_numbers,

        friction_factors[epsilon],

        label=f"ε = {epsilon}",

        linewidth=2

    )

# تنظیمات نمودار

plt.xscale("log")

plt.yscale("log")

plt.xlabel("Reynolds Number (log scale)", fontsize=14)

plt.ylabel("Friction Factor f (log scale)", fontsize=14)

plt.title("Colebrook Equation - Friction Factor vs Reynolds Number", fontsize=16)

plt.legend(title="Roughness (ε)", fontsize=12)

plt.grid(which="both", linestyle="--", linewidth=0.5)

plt.tight_layout()

# نمایش نمودار

plt.show()

خروجی ان شکلی هست 

این پروژه شامل 3 فایل مهم است، لطفا قبل از ارسال پیشنهاد حتما نسبت به بررسی این فایل اقدام فرمایید.

مهارت ها و تخصص های مورد نیاز

---

---