كيف يمكن أن نحقق الديناميكا الحرارية تحليل خلية أوتو مهندس يستخدم ماتلاب

فهم نموذج قناة "أوتو سيكل" الجوية

ويُعدّ هذا النموذج المثالي للحركة الأربع التي تعمل بالسيارات، ويحل محل المحركات المعقدة، والضغط، والحرق، وعمليات العادم بأربع خطوات قابلة للنقض تعزل آليات تحويل الطاقة الأساسية، وفي نموذج قياس الهواء، يُعالج سوائل العمل (تُعالج كغاز مثالي ذي حرارة ثابتة)، ويُضاف الاحتراق إلى الناتج الفوري المستمر.

وتتمثل العمليات الأربع التي تشكل الدورة في ما يلي:

• Isentropic compression (1 ⁇ 2): The piston compresses the charge without heat transfer, raising temperature and pressure.
• Constant —volume heat addition (2 ⁇ 3)]: An momentaneous release of energy (representing combustion) at top dead centre raises the temperature and pressure to their toptom values.
• Isentropic expansion (3 ⁇ 4)]: The high —pressure gases push the piston down, deliver net work to the crankshaft.
• Constant —volume heat rejection (4 ⁇ 1): The cylinder contents lose heat back to the initial state, closing the cycle.

ونظراً لأن النموذج يستخدم نسبة حرارة ثابتة محددة ويتجاهل تبادل الغاز، فإنه يوفر معياراً سريعاً يمكن تنقيحه بصورة منهجية، كما أن تنفيذ هذه الدورة في MATLAB يعطيك أداة تفاعلية مرنة لاستكشاف مدى تأثير نسبة الضغط، وإضافات الحرارة، وخواص السوائل على الأداء.

إقامة بيئة مهيأة

(أ) بداية من إنشاء نص جديد أو ملف عمل جديد: يتيح حيز عمل نظيف إجراء تعديلات مباشرة على التزييف والتعديلات المستقبلية، واستخدام رئيس وصفي يوثق الغرض، ومعايير المدخلات، والنواتج، والصفائف السابقة لمتغيرات الدولة عند مسح البارامترات لتجنب إعادة التوازن الدينامية للرؤوس العامة. (أ) اعتماد أسماء متغيرة ذات معنى مثل بدلاً من مجرد (البنود القياسية).

% otto_cycle_analysis.m
% Air‑standard Otto cycle analysis
% Inputs:
% r – compression ratio (dimensionless)
% T1, P1 – intake state (K, Pa)
% gamma – specific heat ratio (cp/cv)
% Q_in – specific heat addition (J/kg)
% Outputs:
% state_table – array of T, P, v at four corners
% W_net, eta – specific net work (J/kg) and thermal efficiency

(ب) استخدام صفائف MATLAB أو لتخزين الولايات الأربع يحسن من إمكانية القراءة عندما تقطع لاحقاً أو تُنتج نتائج التصدير، مثلاً، إنشاء صفيفة هيكلية مع حقول ]، ، مفهرسة من 1 إلى 4.

تحديد البارامترات الهندسية الحرجة

وتتوقف دقة المحاكاة على بيانات المدخلات، وتحدد على الأقل نسبة الضغط وضغط الاستيعاب ودرجة الحرارة المحددة ونسبة الحرارة المحددة والطاقة المضافة لكل كتلة من الوحدات، والقيم النموذجية لمحرك البنزين هي:

• Comppression ratio r]: 8-12 for modern spark —hinignition motorignition motors; higher ratios improve efficiency but increase knock risk. In research motors, ratios as high as 14 are sometimes explored with advanced fuels.
• Intake temperature T1]: 300 K (ambient) or 350 K if accounting for intake air heating from the manifold and residual gases.
• Intake pressure P1]: 101325 Pa (1 atm) for naturally aspirated; up to 200 kPa or more for turbocharged motors.
• Specific heat ratio غاما ]: 1.4 للهواء الجاف؛ 1.3-1.35 بالنسبة للمخاليط الوقودية - الهواء ومنتجات الاحتراق بسبب وجود جزيئات ثلاثية الأبعاد مثل ثاني أكسيد الكربون وه2O.
• Specific heat addition Q in]: حوالي 2.5-3 MJ/kg for gasoline based on lower heating value and stoichiometric air —fuel-fuel. For a typical gasoline, the lower heating value is about 44 MJ/kg, but mixed with air (roughly 14.7:1 airfuelfuelity ratio), the mixture dens energy

Define these as MATLAB variables at the top of the script so they can be easily changed for sensitivity studies.

r = 9.0; % compression ratio
T1 = 300; % K
P1 = 101325; % Pa
gamma = 1.35; % approximate for fuel‑air mixture
Q_in = 2500e3; % J/kg
R_air = 287; % J/(kg·K)

يمكنك أيضاً أن تُضمّن الكتلة الطرية من الهواء (28.97 غرام/ملت) إذا تحولت لاحقاً إلى حسابات مُتّسمة بالجاسوسية، والاحتفاظ بالغاز الثابت يتجنب أخطاء الوحدة.

حوسبة دول العشائر الأربع

ومع تحديد المعايير، فإن العلاقات مع مستوى الهواء تعطي الدولة نقاطا مباشرة، والتحليل هو كل كتلة وحدة، وبالتالي نعمل بحجم محدد، ويبسط حجم التطهير إلى الوحدة (من 5 إلى 1) التصعيد النسبي؛ ثم من 1 إلى 1 إلى 1 لأن نسبة الضغط من 1 إلى 5 إلى 1 إلى 1 إلى 1، أي س.

الدولة 1 (المستقبل، مركز الموتى القاتم)

v1 = r; % relative to clearance volume
% T1 and P1 are inputs

الدولة 2 (الضغط على المركز الأعلى للموتى)

بالنسبة لعملية مُتَعَدِّدة مع غاما الثابتة:

v2 = 1;
T2 = T1 * r^(gamma-1);
P2 = P1 * r^gamma;

هذه العلاقات مستمدة من القانون المثالي للغاز وعلاقة البوليتروبيك P vga = ثابت .

الدولة 3 (إضافة الحرارة الثابتة - المُعدية)

أولاً، تُحسب الحرارة المحددة في الحجم الثابت:

cv = R_air / (gamma-1);
T3 = T2 + Q_in / cv;
P3 = P2 * (T3 / T2);
v3 = v2; % no volume change

ويمكن أن تتجاوز درجة الحرارة في الولاية ٣ ٠٠٥٢ كاف بالنسبة للاحتراق النموذجي للبنزين، وكثيرا ما يصل ضغط البقاع من نوع P3 إلى ٦٠ - ٠٨ بار، وهو ما يحدد الاحتياجات الهيكلية لمبنى المحرك ورأسه.

الدولة 4 (التوسع في الحجم الأصلي)

أما نسبة التوسع من درجة حرارة التقلبات إلى نسبة التراكم البيولوجي فهي تعادل، ولكن من حيث الحجم المحدد مقابل 3 = 2 = 1 و4 = 1 = 1 = 1 = س، وبالتالي فإن نسبة الحجم هي 1/ص.

v4 = v1;
T4 = T3 * (1/r)^(gamma-1);
P4 = P3 * (1/r)^gamma;

وهذه النقاط الأربع تحدد تماما الدورة من حيث درجة الحرارة والضغط والحجم المحدد، وتخزنها في صفائف للتخطيط والحساب فيما بعد، ومن المفيد التحقق من أن درجة P4 أعلى قليلا من الرتبة ف-1؛ والفرق يدفع إلى الانهيار في العادم.

حساب صافي العمل والكفاءة الحرارية، والضغط الفعال

أما العمل الصافي المحدد فهو الفرق بين الحرارة المضافة والحرارة المرفوضة، أما بالنسبة لنموذج قياس الهواء، فإن الحرارة المرفوضة في العملية 4-1 هي:

Q_out = cv * (T4 - T1);
W_net = Q_in - Q_out;

الكفاءة الحرارية هي:

eta = W_net / Q_in; % from energy balance
eta_direct = 1 - 1/r^(gamma-1); % direct formula
% Both should match within numerical round-off

Mean effective pressure (MEP)] is a performance metric that comparisons net work per cycle to the displacement volume. For the air —standard cycle, the specific displacement volume is v1 - v2 = r - 1. Therefore:

mep = W_net / (v1 - v2); % Pa

وشركة الطاقة البيئية مستقلة عن حجم المحرك، وهي مؤشر ملائم لنوعية تصميم المحركات، والقيم النموذجية لمحركات الغازولين المطهرة طبيعياً هي 8-12 باراً (0.8-1.2 ميغابايا)، وشركة الإنتاج المتعددة الأطراف إلى جانب الكفاءة لتوفير صورة أكمل، ويمكن أيضاً أن تتحول إلى بار للقابلية للقراءة: .

بما في ذلك القفازات المحددة لتحسين الاستحقاق

ويمكن للافتراض المستمر في غاما أن يبالغ في تقدير الكفاءة لأن غازات العادم الحقيقية لديها نسب حرارة محددة عند درجات حرارة عالية، ولحساب ذلك، يستعاض عن الغاما الثابتة بترابط متصل بدرجة الحرارة، ويستخدم نموذج البوليومات في وكالة ناسا (NASA) كمعاملات متعددة المقاييس (FLT:0)(ج) _{(ب) [النموذج: 2 من النوع الثاني من المواد الكيميائية]}

ويتطلب تنفيذ المتغيرات في غاما حلا متكررا لأن المسرّع نفسه يعتمد على درجة الحرارة غير المعروفة، واستخدام ماتيلبار أو حلقة بسيطة، مثلا، أثناء الضغط على الجليد، تصبح العلاقة كما يلي:

% Use an anonymous function to find T2 such that
% integral of cp dT from T1 to T2 = R * ln(v1/v2)
% This replaces the constant‑gamma formula.
% For air, cp(T) can be approximated by:
% cp = a0 + a1*T + a2*T^2 + a3*T^3 + a4*T^4
% Using NASA 7-coefficient polynomials.
fun = @(T2) integral(@(T) polyval(cp_coeffs, T), T1, T2) - R_air * log(r);
T2 = fzero(fun, T1 * r^0.4);

ويوفر NASA البوليوميوميومية ] معاملا قابلا للتنزيل بالنسبة لأنواع الهواء والحرق، ويقلل نصك من تحميل هذه المعاملات ويتصل بها أثناء حسابات الدولة، وستلاحظون أن حرارات معينة متغيرة تقلل من الكفاءة المتوقعة بنسبة تتراوح بين 2 و5 نقاط مئوية في النسب المضغوطة النموذجية، التي تتواءم بشكل وثيق مع بيانات المحركات الحقيقية.

ولزيادة التقريب السريع دون التكامل الكامل بين التعددية، تستخدم نسبة حرارة محددة متباينة بشكل متواز مع درجة الحرارة: غاما (T) = 1.4 - 0.00005*(T - 300) وهذا الترابط البسيط يستوعب الاتجاه ويمكن تنفيذه في وقت ما إلى حين التقارب.

بناء تحليل أداء نظامي

وبدلا من النص الأحادي، يلخص منطق الدورة في وظيفة MATLAB، مما يجعله قابلا لإعادة استخدامه في عمليات مسح البارامترات ويقلل من الازدواجية في استخدام الرموز، ويرد أدناه هيكل للوظيفة:

function [states, performance] = ottoCycle(r, T1, P1, Q_in, gamma, R_air)
 % states: structure with T, P, v for states 1-4
 % performance: structure with W_net, eta, mep
 % (calculation code from previous section goes here)
end

بالنسبة للحرارات المحددة المتغيرة، تضاف علم سوائل:

function [states, performance] = ottoCycle(r, T1, P1, Q_in, fluid_type)
 % fluid_type: 'constant_gamma' or 'variable_cp'
 % (logic branches based on flag)
end

مع وجود الوظيفة، يمكن لنص السائق أن يمسح من خلال نسب الضغط:

r_values = 6:1:14;
for i = 1:length(r_values)
 [~, perf] = ottoCycle(r_values(i), 300, 101325, 2500e3, 'constant_gamma');
 eta(i) = perf.eta;
 mep(i) = perf.mep;
end
plot(r_values, eta*100, 'b-', r_values, mep/1e5, 'r--');

ويظهر هذا النهج النموذجي على الفور انخفاض العائدات من زيادة نسبة الضغط إلى ما بعد 12 عاماً، وتبرز المؤامرة المبادلات الرئيسية: فالكفاءة لا تزال ترتفع ولكن بمعدل تناقص، في حين أن برنامج الطاقة المتعدد الأطراف كثيراً ما يصل إلى ذروته نحو 10 لترات بالنسبة للمدخلات الحرارية النموذجية.

إنتاج تشخيصي للأشعة المقطعية وأشعة T - S

وتساعد التصورات على التحقق من النموذج وعلى الإبلاغ عن النتائج، وتنتج وظائف تخطيط برنامج العمل المتعدد التخصصات برامج للنشر.

مضخة الضغط (PANV)

The compression and expansion curves follow P vغاما = constant]. Generate a series of medium points for smooth curves:

n = 50;
v_comp = linspace(v2, v1, n); % reverse direction for plotting
P_comp = P2 * (v_comp ./ v2).^(-gamma);
v_exp = linspace(v3, v4, n);
P_exp = P3 * (v_exp ./ v3).^(-gamma);
figure;
hold on;
plot(v_comp, P_comp, 'b', 'LineWidth', 1.5);
plot(v_exp, P_exp, 'b', 'LineWidth', 1.5);
plot([v2 v3], [P2 P3], 'r--', 'LineWidth', 1.2);
plot([v4 v1], [P4 P1], 'r--', 'LineWidth', 1.2);
xlabel('Specific volume (relative)'); ylabel('Pressure (Pa)');
title('Air‑Standard Otto Cycle P–V Diagram');
legend('Compression', 'Expansion', 'Heat addition', 'Heat rejection');
grid on;

ويمثل المجال المرفق العمل الصافي، ويمكنكم أن تدمجوا في عدد من الحلقات باستخدام [(FLT:26] للتحقق من الشبكة الحاسوبية، وعلى سبيل المثال، دمج المنحنى الأعلى (الضغط + التوسع) وطرح المنحنى الأدنى (خط الرفض الحراري) والإحالة إلى وثائق [(MATLAB بشأن التكامل الرقمي .

مقياس درجة الحرارة - الانتروبية (TINS)

بالنسبة لغاز مثالي ذي حرارة محددة ثابتة، فإن التغير في النسخ مقارنة بالوضع 1 هو:

% s = cp * log(T / T1) - R * log(P / P1)
cp = gamma * R_air / (gamma - 1);
s1 = 0;
s2 = cp * log(T2/T1) - R_air * log(P2/P1);
s3 = cp * log(T3/T1) - R_air * log(P3/P1);
s4 = cp * log(T4/T1) - R_air * log(P4/P1);

(ب) استخدام الدمج الرقمي للترددات/الثانية في حساب الطراز، ويكشف الرسم البياني لدائرة T -S عن الطابع المائي للعمليات 1-2 و 3-4 (الخطوط الخارجية) والتدفئة/الرفض المستمرين (الممرات المكشوفة) بينما تمثل المنطقة تحت منحنى الحرارة المضافة، بينما تمثل المنطقة تحت خط التسخين الثابت.

ولتحسين مؤامرة TINS، تُتجاوز خطوط الضغط المستمر (الحدود) من العلاقة المثلى للغاز، وهذا يساعد على وضع السياق للدورة في إطار أوسع نطاقاً للدماغ الحراري.

إجراء دراسات توعية

وبعد أن يعمل نموذج خط الأساس، يُغيّر بانتظام معايير المدخلات لفهم تأثيرها، وتشمل الدراسات المشتركة ما يلي:

• Comppression ratio sweep]: Plot efficiency and MEP against r. Efficiency increases rapidly at low r and then plateaus; MEP often tops at a moderate r because of the trade —off between pressure rise and heat addition. Use subplots to show both metrics on the same figure.
• (ب) زيادة معدلات الطاقة الكهربائية/الوقود الجوي؛ وزيادة كمية الطاقة العالية/العمل على الحدود القصوى للمواد؛ وتشمل قيوداً تفرضها P3 على 150 باراً (ميغابايا) تعكس قوة المحرك النموذجية.
• Specific heat ratio effect]: أنخفض غاما الأدنى من الكفاءة، وأن تنفق على النتائج الثابتة والمتغيرة لغما لتظهر التصويب، وهذا أمر مهم بصفة خاصة عند تحليل أنواع الوقود البديلة مثل الغاز الطبيعي أو الهيدروجين.

من أجل مسح مقاسين، استخدام حلقات مُستَنَقَة، و المخزن يُنتج في مصفوفة، ثمّ يُنشئ مؤامرة أو خريطة لونية لتصوير التأثير المشترك:

r_vec = 6:1:14;
Q_vec = 2e6:0.2e6:3e6;
for i = 1:length(r_vec)
 for j = 1:length(Q_vec)
 [~, perf] = ottoCycle(r_vec(i), 300, 101325, Q_vec(j), 1.35, 287);
 eta_mat(i,j) = perf.eta;
 end
end
contourf(r_vec, Q_vec/1e6, eta_mat', 20);
xlabel('Compression ratio'); ylabel('Q_in (MJ/kg)');
colorbar; title('Thermal Efficiency (%)');

وهذه الخرائط الحساسية لا تقدر بثمن بالنسبة لتحليل المبادلات في تصميم المحركات، فهي تساعد على تحديد مناطق التشغيل التي تكون فيها الكفاءة عالية بينما يظل الضغط على ذروته مقبولا.

تقييم نتائجك

(أ) قبل استخدام نموذج القرارات الهندسية، التحقق من المعايير المعروفة، فإن صيغة الكفاءة في الهواء (غاما) = 1 - 1/rغاما (-1) توفر التحقق المباشر من الحالة المستمرة، إضافة إلى مقارنة الضغط في ذروة مدروسة من البيانات المستمدة من خرائط المحركات المنشورة أو نماذج الكتب المدرسية (مثلاً، ستون، Introduction to Internal Combustion Engines [FLT]

results = table(r_values', eta', mep', ...
 'VariableNames', {'CompressionRatio','Efficiency','MEP_bar'});
writetable(results, 'otto_sweep.csv');