חשב הבדלי לחץ על פני שטח נוזלים מעוקלים באמצעות משוואת יאנג-לאפלאס. הזן מתח שטח ורדיוסים עיקריים של עיקול כדי לנתח טיפות, בועות ותופעות קפילריות.
ΔP = γ(1/R₁ + 1/R₂)
ΔP = 0.072 × (1/0.001 + 1/0.001)
ΔP = 0.072 × (1000.00 + 1000.00)
ΔP = 0.072 × 2000.00
ΔP = 0.00 Pa
ההדמיה הזו מציגה את הממשק המעוקל עם רדיוסי העיקול העיקריים R₁ ו-R₂. החיצים מצביעים על הפרש הלחצים לאורך הממשק.
משוואת יאנג-לאפלאס היא נוסחה יסודית במכניקת נוזלים שמתארת את ההפרש בלחץ על פני ממשק מעוקל בין שני נוזלים, כגון ממשק נוזל-גז או נוזל-נוזל. הפרש הלחץ הזה נובע ממתח פני השטח ומהעיקול של הממשק. ה- פותר משוואת יאנג-לאפלאס שלנו מספק דרך פשוטה ומדויקת לחשב את הפרש הלחץ הזה על ידי הזנת מתח פני השטח ורדיוסי העיקול הראשיים. בין אם אתם לומדים טיפות, בועות, פעולה קפילרית או תופעות פני שטח אחרות, כלי זה מציע פתרונות מהירים לבעיות מורכבות של מתח פני השטח.
הנוסחה, שנקראת על שם תומס יאנג ופייר-סימון לפלאס שפיתחו אותה בתחילת המאה ה-19, חיונית במספר יישומים מדעיים והנדסיים, ממיקרופלואידיקה ומדעי החומרים ועד מערכות ביולוגיות ותהליכים תעשייתיים. על ידי הבנת הקשר בין מתח פני השטח, העיקול והפרש הלחץ, חוקרים ומהנדסים יכולים לעצב ולנתח טוב יותר מערכות המעורבות בממשקי נוזלים.
משוואת יאנג-לאפלאס מקשרת בין ההפרש בלחץ על פני ממשק נוזלי לבין מתח פני השטח ורדיוסי העיקול הראשיים:
כאשר:
לממשק סגלגל (כגון טיפה או בועה), כאשר , הנוסחה מתפשטת ל:
מתח פני השטח ():
רדיוסי העיקול הראשיים ( ו-):
הפרש הלחץ ():
סימן ההגדרה עבור משוואת יאנג-לאפלאס הוא חשוב:
משטח שטוח: כאשר אחד הרדיוסים מתקרב לאינסוף, התרוממותו להפרש הלחץ מתקרבת לאפס. עבור משטח שטוח לחלוטין (), .
משטח צילינדרי: עבור משטח צילינדרי (כגון נוזל בצינור קפילרי), רדיוס אחד סופי () בעוד השני אינסופי (), מה שמוביל ל-.
רדיוסים קטנים מאוד: בקנה מידה מיקרוסקופי (למשל, טיפות ננומטריות), השפעות נוספות כמו מתח קו עשויות להיות משמעותיות, והמשוואה הקלאסית של יאנג-לאפלאס עשויה לדרוש תיקון.
השפעות טמפרטורה: מתח פני השטח בדרך כלל יורד עם עליית טמפרטורה, מה שמשפיע על הפרש הלחץ. קרוב לנקודת הקריטית, מתח פני השטח מתקרב לאפס.
סורפקטנטים: נוכחותם של סורפקטנטים מפחיתה את מתח פני השטח ולכן את הפרש הלחץ על פני הממשק.
המחשב שלנו מספק דרך פשוטה לקבוע את ההפרש בלחץ על פני ממשקים מעוקלים של נוזלים. עקבו אחרי הצעדים הבאים כדי לקבל תוצאות מדויקות:
הזינו את מתח פני השטח ():
הזינו את רדיוס העיקול הראשי הראשון ():
הזינו את רדיוס העיקול הראשי השני ():
צפו בתוצאה:
העתיקו או שתפו את התוצאות:
משוואת יאנג-לאפלאס יש לה יישומים רבים בתחומים מדעיים והנדסיים שונים:
המשוואה היא יסודית להבנת ההתנהגות של טיפות ובועות. היא מסבירה מדוע טיפות קטנות יותר יש להן לחץ פנימי גבוה יותר, מה שמניע תהליכים כמו:
משוואת יאנג-לאפלאס עוזרת להסביר ולכמת עלייה קפילרית או ירידה:
במדע הרפואה והביולוגיה, המשוואה משמשת ל:
יישומים בפיתוח חומרים כוללים:
רבים מהיישומים התעשייתיים מסתמכים על הבנת הפרשי הלחץ בממשקים:
נשקול טיפה סגלגלה של מים עם רדיוס של 1 מ"מ ב-20°C:
זה אומר שלחץ הפנימי בטיפה גבוה ב-144 Pa מהלחץ באוויר הסובב.
בעוד שמשוואת יאנג-לאפלאס היא יסודית, ישנן גישות וארכות חלופיות למצבים ספציפיים:
משוואת קלווין: מקשרת בין לחץ האדים על פני משטח נוזלי מעוקל לזה על פני שטח שטוח, שימושית לחקר התעבות ואידוי.
אפקט גיבס-תומסון: מתאר כיצד גודל חלקיקים משפיע על מסיסות, נקודת התכה ותכונות תרמודינמיות אחרות.
מודל הלפריך: מאריך את הניתוח לממברנות אלסטיות כמו ממברנות ביולוגיות, כולל קפיצים.
סימולציות נומריות: עבור גיאומטריות מורכבות, שיטות חישוביות כמו שיטת נפח הנוזל (VOF) או שיטות קו רמה עשויות להיות מתאימות יותר מאשר פתרונות אנליטיים.
דינמיקה מולקולרית: בקנה מידה קטן מאוד (ננומטרים), ההנחות הקונטינואיטיביות מתפרקות, וסימולציות דינמיקה מולקולרית מספקות תוצאות מדויקות יותר.
פיתוח משוואת יאנג-לאפלאס מייצג אבן דרך משמעותית בהבנת תופעות פני שטח וקפילריות.
לימוד פעולה קפילרית מתוארך לעבר, אך חקירה מדעית שיטתית החלה בתקופת הרנסנס:
המשוואה כפי שאנו מכירים אותה היום צמחה מעבודתם של שני מדענים שפעלו באופן עצמאי:
תומס יאנג (1805): פרסם "מאמר על קוהזיה של נוזלים" בטרנסאקציות הפילוסופיות של החברה המלכותית, והציג את המושג של מתח פני השטח ואת הקשר שלו להפרשי לחצים על פני ממשקים מעוקלים.
פייר-סימון לפלאס (1806): בספרו המונומנטלי "מכניקה שמימית", לפלאס פיתח מסגרת מתמטית לפעולה קפילרית, והפיק את המשוואה המקשרת בין הפרש הלחץ לעיקול פני השטח.
השילוב של תובנות פיזיקליות של יאנג ורצינות מתמטית של לפלאס הוביל למה שאנו מכנים היום משוואת יאנג-לאפלאס.
בעקבות המאה ה-19, המשוואה שופרה והורחבה:
היום, משוואת יאנג-לאפלאס נשארת אבן יסוד במדעי הממשקים, וממשיכה למצוא יישומים חדשים ככל שהטכנולוגיה מתקדמת למיקרו וננו סקאלות.
הנה יישומים של משוואת יאנג-לאפלאס בשפות תכנות שונות:
1' נוסחת Excel עבור משוואת יאנג-לאפלאס (ממשק סגלגל)
2=2*B2/C2
3
4' כאשר:
5' B2 מכיל את מתח פני השטח ב-N/m
6' C2 מכיל את הרדיוס ב-m
7' התוצאה היא ב-Pa
8
9' עבור המקרה הכללי עם שני רדיוסים ראשיים:
10=B2*(1/C2+1/D2)
11
12' כאשר:
13' B2 מכיל את מתח פני השטח ב-N/m
14' C2 מכיל את הרדיוס הראשון ב-m
15' D2 מכיל את הרדיוס השני ב-m
16
1def young_laplace_pressure(surface_tension, radius1, radius2):
2 """
3 Calculate pressure difference using the Young-Laplace equation.
4
5 Parameters:
6 surface_tension (float): Surface tension in N/m
7 radius1 (float): First principal radius of curvature in m
8 radius2 (float): Second principal radius of curvature in m
9
10 Returns:
11 float: Pressure difference in Pa
12 """
13 if radius1 == 0 or radius2 == 0:
14 raise ValueError("Radii must be non-zero")
15
16 return surface_tension * (1/radius1 + 1/radius2)
17
18# Example for a spherical water droplet
19surface_tension_water = 0.072 # N/m at 20°C
20droplet_radius = 0.001 # 1 mm in meters
21
22# For a sphere, both radii are equal
23pressure_diff = young_laplace_pressure(surface_tension_water, droplet_radius, droplet_radius)
24print(f"Pressure difference: {pressure_diff:.2f} Pa")
25
1/**
2 * Calculate pressure difference using the Young-Laplace equation
3 * @param {number} surfaceTension - Surface tension in N/m
4 * @param {number} radius1 - First principal radius of curvature in m
5 * @param {number} radius2 - Second principal radius of curvature in m
6 * @returns {number} Pressure difference in Pa
7 */
8function youngLaplacePressure(surfaceTension, radius1, radius2) {
9 if (radius1 === 0 || radius2 === 0) {
10 throw new Error("Radii must be non-zero");
11 }
12
13 return surfaceTension * (1/radius1 + 1/radius2);
14}
15
16// Example for a water-air interface in a capillary tube
17const surfaceTensionWater = 0.072; // N/m at 20°C
18const tubeRadius = 0.0005; // 0.5 mm in meters
19// For a cylindrical surface, one radius is the tube radius, the other is infinite
20const infiniteRadius = Number.MAX_VALUE;
21
22const pressureDiff = youngLaplacePressure(surfaceTensionWater, tubeRadius, infiniteRadius);
23console.log(`Pressure difference: ${pressureDiff.toFixed(2)} Pa`);
24
1public class YoungLaplaceCalculator {
2 /**
3 * Calculate pressure difference using the Young-Laplace equation
4 *
5 * @param surfaceTension Surface tension in N/m
6 * @param radius1 First principal radius of curvature in m
7 * @param radius2 Second principal radius of curvature in m
8 * @return Pressure difference in Pa
9 */
10 public static double calculatePressureDifference(double surfaceTension, double radius1, double radius2) {
11 if (radius1 == 0 || radius2 == 0) {
12 throw new IllegalArgumentException("Radii must be non-zero");
13 }
14
15 return surfaceTension * (1/radius1 + 1/radius2);
16 }
17
18 public static void main(String[] args) {
19 // Example for a soap bubble
20 double surfaceTensionSoap = 0.025; // N/m
21 double bubbleRadius = 0.01; // 1 cm in meters
22
23 // For a spherical bubble, both radii are equal
24 // Note: For a soap bubble, there are two interfaces (inner and outer),
25 // so we multiply by 2
26 double pressureDiff = 2 * calculatePressureDifference(surfaceTensionSoap, bubbleRadius, bubbleRadius);
27
28 System.out.printf("Pressure difference across soap bubble: %.2f Pa%n", pressureDiff);
29 }
30}
31
1function deltaP = youngLaplacePressure(surfaceTension, radius1, radius2)
2 % Calculate pressure difference using the Young-Laplace equation
3 %
4 % Inputs:
5 % surfaceTension - Surface tension in N/m
6 % radius1 - First principal radius of curvature in m
7 % radius2 - Second principal radius of curvature in m
8 %
9 % Output:
10 % deltaP - Pressure difference in Pa
11
12 if radius1 == 0 || radius2 == 0
13 error('Radii must be non-zero');
14 end
15
16 deltaP = surfaceTension * (1/radius1 + 1/radius2);
17end
18
19% Example script to calculate and plot pressure vs. radius for water droplets
20surfaceTension = 0.072; % N/m for water at 20°C
21radii = logspace(-6, -2, 100); % Radii from 1 µm to 1 cm
22pressures = zeros(size(radii));
23
24for i = 1:length(radii)
25 % For spherical droplets, both principal radii are equal
26 pressures(i) = youngLaplacePressure(surfaceTension, radii(i), radii(i));
27end
28
29% Create log-log plot
30loglog(radii, pressures, 'LineWidth', 2);
31grid on;
32xlabel('Droplet Radius (m)');
33ylabel('Pressure Difference (Pa)');
34title('Young-Laplace Pressure vs. Droplet Size for Water');
35
1#include <iostream>
2#include <stdexcept>
3#include <cmath>
4#include <iomanip>
5
6/**
7 * Calculate pressure difference using the Young-Laplace equation
8 *
9 * @param surfaceTension Surface tension in N/m
10 * @param radius1 First principal radius of curvature in m
11 * @param radius2 Second principal radius of curvature in m
12 * @return Pressure difference in Pa
13 */
14double youngLaplacePressure(double surfaceTension, double radius1, double radius2) {
15 if (radius1 == 0.0 || radius2 == 0.0) {
16 throw std::invalid_argument("Radii must be non-zero");
17 }
18
19 return surfaceTension * (1.0/radius1 + 1.0/radius2);
20}
21
22int main() {
23 try {
24 // Example for a mercury droplet
25 double surfaceTensionMercury = 0.485; // N/m at 20°C
26 double dropletRadius = 0.002; // 2 mm in meters
27
28 // For a spherical droplet, both radii are equal
29 double pressureDiff = youngLaplacePressure(surfaceTensionMercury, dropletRadius, dropletRadius);
30
31 std::cout << "Pressure difference inside mercury droplet: "
32 << std::fixed << std::setprecision(2) << pressureDiff
33 << " Pa" << std::endl;
34
35 // Example for a cylindrical interface (like in a capillary tube)
36 double tubeRadius = 0.0001; // 0.1 mm
37 double infiniteRadius = std::numeric_limits<double>::max();
38
39 double capillaryPressure = youngLaplacePressure(surfaceTensionMercury, tubeRadius, infiniteRadius);
40
41 std::cout << "Pressure difference in mercury capillary: "
42 << std::fixed << std::setprecision(2) << capillaryPressure
43 << " Pa" << std::endl;
44 }
45 catch (const std::exception& e) {
46 std::cerr << "Error: " << e.what() << std::endl;
47 return 1;
48 }
49
50 return 0;
51}
52
1#' Calculate pressure difference using the Young-Laplace equation
2#'
3#' @param surface_tension Surface tension in N/m
4#' @param radius1 First principal radius of curvature in m
5#' @param radius2 Second principal radius of curvature in m
6#' @return Pressure difference in Pa
7#' @examples
8#' young_laplace_pressure(0.072, 0.001, 0.001)
9young_laplace_pressure <- function(surface_tension, radius1, radius2) {
10 if (radius1 == 0 || radius2 == 0) {
11 stop("Radii must be non-zero")
12 }
13
14 return(surface_tension * (1/radius1 + 1/radius2))
15}
16
17# Example: Compare pressure differences for different liquids with the same geometry
18liquids <- data.frame(
19 name = c("Water", "Ethanol", "Mercury", "Benzene", "Blood plasma"),
20 surface_tension = c(0.072, 0.022, 0.485, 0.029, 0.058)
21)
22
23# Calculate pressure for a 1 mm radius spherical droplet
24droplet_radius <- 0.001 # m
25liquids$pressure <- sapply(liquids$surface_tension, function(st) {
26 young_laplace_pressure(st, droplet_radius, droplet_radius)
27})
28
29# Create a bar plot
30barplot(liquids$pressure, names.arg = liquids$name,
31 ylab = "Pressure Difference (Pa)",
32 main = "Laplace Pressure for 1 mm Droplets of Different Liquids",
33 col = "lightblue")
34
35# Print the results
36print(liquids[, c("name", "surface_tension", "pressure")])
37
משוואת יאנג-לאפלאס משמשת לחישוב ההפרש בלחץ על פני ממשק נוזלי מעוקל עקב מתח פני השטח. היא חיונית להבנת תופעות כמו פעולה קפילרית, יצירת טיפות, יציבות בועות ומגוון יישומים מיקרופלואידיים. המשוואה עוזרת למהנדסים ולחוקרים לעצב מערכות המעורבות בממשקי נוזלים ולחזות כיצד יתנהגו בתנאים שונים.
לטיפות קטנות יש לחץ פנימי גבוה יותר בגלל העיקול שלהן. לפי משוואת יאנג-לאפלאס, ההפרש בלחץ הפוך פרופורציונלי לרדיוס העיקול. ככל שהרדיוס קטן יותר, העיקול (1/R) גדל, מה שמוביל להפרש לחץ גבוה יותר. זה מסביר מדוע טיפות מים קטנות מתאדות מהר יותר מטיפות גדולות יותר ולמה בועות קטנות יותר בקצף נוטות להתכווץ בזמן שבועות גדולות יותר גדלות.
הטמפרטורה משפיעה בעיקר על משוואת יאנג-לאפלאס דרך השפעתה על מתח פני השטח. עבור רוב הנוזלים, מתח פני השטח יורד בקירוב ליניארי עם עליית הטמפרטורה. זה אומר שההפרש בלחץ על פני ממשק מעוקל גם יפחת ככל שהטמפרטורה עולה, בהנחה שהגיאומטריה נשארת קבועה. קרוב לנקודת הקריטית של נוזל, מתח פני השטח מתקרב לאפס, והאפקט של יאנג-לאפלאס הופך להיות זניח.
כן, הצורה הכללית של משוואת יאנג-לאפלאס חלה על כל ממשק מעוקל, לא רק על סגלגלים. המשוואה משתמשת בשני רדיוסים ראשיים של עיקול, שעשויים להיות שונים עבור משטחים לא סגלגלים. עבור גיאומטריות מורכבות, רדיוסים אלו עשויים להשתנות מנקודה לנקודה לאורך פני השטח, ודורשים טיפול מתמטי מתקדם יותר או שיטות נומריות כדי לפתור את כל צורת הממשק.
משוואת יאנג-לאפלאס מסבירה ישירות את העלייה הקפילרית. בצינור צר, המניסקוס המעוקל יוצר הפרש לחץ לפי המשוואה. הפרש הלחץ הזה מניע את הנוזל כלפי מעלה נגד כוח הכבידה עד שהשוויון מושג. גובה העלייה הקפילרית נגזר על ידי קביעת ההפרש בלחץ ממשוואת יאנג-לאפלאס שווה ללחץ ההידרוסטטי של עמוד הנוזל המוגבה (ρgh), מה שמוביל לנוסחה הידועה h = 2γcosθ/(ρgr).
משוואת יאנג-לאפלאס בדרך כלל מדויקת עד לקני מידה מיקרוסקופיים (מיקרומטרים), אך בקני מידה ננומטריים, השפעות נוספות נעשות משמעותיות. אלו כוללות מתח קו (בקו המגע של שלוש הפאזות), לחץ דיזוייני (בסרטים דקים) ואינטראקציות מולקולריות. בקני מידה אלו, ההנחה הקונטינואיטיבית מתחילה להתפרק, והמשוואה הקלאסית של יאנג-לאפלאס עשויה לדרוש תיקון או החלפה בגישות דינמיקה מולקולרית.
בעוד ששתי המשוואות קשורות, הן מתארות היבטים שונים של ממשקי נוזלים. משוואת יאנג-לאפלאס מקשרת בין הפרש לחץ לעיקול ומתח פני השטח. משוואת יאנג (לעיתים נקראת גם יחס יאנג) מתארת את זווית המגע הנוצרת כאשר ממשק נוזל-אוויר נפגש עם משטח מוצק, מקשרת אותה למתחים הפנימיים בין שלוש הפאזות (מוצק-אוויר, מוצק-נוזל ונוזל-אוויר). שתי המשוואות פותחו על ידי תומס יאנג והן יסודיות בהבנת תופעות פני שטח.
סורפקטנטים מפחיתים את מתח פני השטח על ידי פרסום על פני הממשק הנוזלי. לפי משוואת יאנג-לאפלאס, זה מפחית ישירות את ההפרש בלחץ על פני הממשק. בנוסף, סורפקטנטים יכולים ליצור גרדיאנטים של מתח פני השטח (אפקטי מרנגוני) כאשר הם לא מפוזרים באופן אחיד, מה שמוביל לזרימות מורכבות והתנהגויות דינמיות שלא נתפסות על ידי משוואת יאנג-לאפלאס הסטטית. זו הסיבה שסורפקטנטים מייצבים קצפים ואמולסיות - הם מפחיתים את ההפרש בלחץ המניע את ההתאגדות.
כן, משוואת יאנג-לאפלאס, בשילוב עם השפעות כבידה, יכולה לחזות את הצורה של טיפה תלויה. עבור מקרים כאלה, המשוואה נכתבת בדרך כלל במונחים של העיקול הממוצע ונפתרת נומרית כבעיה של ערך גבול. גישה זו היא הבסיס לשיטת טיפות תלויה למדידת מתח פני השטח, שבה צורת הטיפה הנצפית מותאמת לפרופילים תיאורטיים המחושבים ממשוואת יאנג-לאפלאס.
למטרות תוצאות עקביות, השתמשו ביחידות SI עם משוואת יאנג-לאפלאס:
אם אתם משתמשים במערכות יחידות אחרות, הקפידו על עקביות. לדוגמה, ביחידות CGS, השתמשו בדין/ס"מ עבור מתח פני השטח, ס"מ עבור רדיוסים ודין/ס"מ² עבור לחץ.
דה גנס, פ.ג., ברושארד-וויירט, פ., & קוּרֵה, ד. (2004). תופעות קפילריות ומתח פני השטח: טיפות, בועות, פנינים, גלים. ספרינגר.
אדמסון, א.וו., & גסט, א.פ. (1997). כימיה פיזיקלית של פני השטח (מהדורה 6). ויילי-אינטרסיינס.
ישראלצ'ווילי, ג'.נ. (2011). כוחות בין מולקולריים ופני שטח (מהדורה 3). אקדמית פרס.
רולינסון, ג'.ס., & ווידום, ב. (2002). תיאוריית מולקולות של קפילריות. דובר פרס.
יאנג, ת. (1805). "מאמר על קוהזיה של נוזלים". טרנסאקציות פילוסופיות של החברה המלכותית של לונדון, 95, 65-87.
לפלאס, פ.ס. (1806). טראיט דה מכניק צלסט (תוספת לספר 10).
דפאי, ר., & פריג'וגין, א. (1966). מתח פני השטח ופרסום. לונגמן.
פין, ר. (1986). פני קפילר שווי משקל. ספרינגר-ורלאג.
דרג'גין, ב.ו., צ'ורייב, נ.ו., & מולר, ו.מ. (1987). כוחות פני השטח. קונסולטנטס ביורו.
לאוטרופ, ב. (2011). פיזיקה של חומר רציף: תופעות אקזוטיות ויומיומיות בעולם המאקרוסקופי (מהדורה 2). CRC Press.
מוכנים לחשב הפרשי לחצים על פני ממשקים מעוקלים? נסו את פותר משוואת יאנג-לאפלאס שלנו עכשיו וקבלו תובנות על תופעות מתח פני השטח. למידע נוסף על כלים ומחשבים במכניקת נוזלים, חקרו את המשאבים האחרים שלנו.
גלה עוד כלים שעשויים להיות שימושיים עבור זרימת העבודה שלך