ניתוח מרחבי-זמני של שדות תרמו-פוטואלסטיים מקושרים בסיליקון אנאיזוטרופי מחוזק בסיבים באמצעות שיטת ערך עצמי

סיכום

מחקר זה מנתח שדות תרמו-פוטואלסטיים מקושרים בסיליקון מחוזק בסיבים אניזוטרופיים באמצעות מוד נורמלי ושיטת ערך עצמי. התוצאות מראות דעיכה מרחבית והתפתחות שדה תלויה בזמן, עם רגישות חזקה לאניזוטרופיה. מפות חום ממחישות את התפלגות השדות ולוקליזציה.

תקציר

מחקר זה בוחן מערכת תרמו-פוטואלסטית מקושרת במדיום סיליקון למחצה מחוזק בסיבים אניזוטרופית, במטרה ללכוד את האינטראקציה בין שדות תרמיים, נשאים ומכניים. חומרים אניזוטרופיים מחוזקים בסיבים ממלאים תפקיד מרכזי ביישומים הנדסיים מודרניים, כולל מכשירים מיקרואלקטרוניים ואופטואלקטרוניים, טכנולוגיות מבוססות לייזר, חיישנים ומבנים מרוכבים מתקדמים, שבהם נדרשות תכונות כיווניות וביצועים מכניים משופרים. הם חשובים במיוחד בתכנון רכיבי מוליכים למחצה העומסים תרמית ואופטית, כאשר חיזוי מדויק של התנהגות השדה המקושר חיוני לאמינות ולאופטימיזציה של ביצועים. המשוואות השולטות מנוסחות על בסיס המודל הפיזיקלי המקושר ומוסבות לאחר מכן לצורה חסרת ממד כדי לפשט את הניתוח ולהדגיש את ההשפעה היחסית של הפרמטרים המעורבים. הבעיה נפתרת באמצעות טכניקת מצב נורמלי ומופחתת למערכת דיפרנציאלית של מטריצת וקטורים מסדר ראשון, ואחריה גישה של ערכי עצמי לקבלת פתרונות אנליטיים המקיימים את תנאי הגבול המוטלים בתוך תחום חצי-אינסופי. ניתוח נומרי מתבצע כדי לבחון את השפעת השונות בזמן על כל השדות הפיזיקליים, וחושף החלשה מרחבית חזקה והתנהגות מקושרת הנשלטת על ידי אניזוטרופיה וחיזוק סיבים. ייצוגי מפת חום מרחבית-זמנית משמשים להמחשת ההתפתחות והלוקליזציה של השדות, המספקים תובנות פיזיקליות על האינטראקציות הרב-פיזיות ומדגימים את יעילות הגישה האנליטית.

מבוא

תרמו-פוטואלסטיות הפכה לתחום רב-תחומי חשוב המתאר את האינטראקציה בין השפעות תרמיות, מכניות ואופטיות בחומרים מוליכים למחצה תחת עירור פוטותרמי. הקישור בין שדות אלו הופך לבולט במיוחד ביישומים מודרניים הכוללים חימום לייזר ועירוור אופטי. לדוגמה, סעיד¹ חקר אינטראקציות תרמו-פוטואלסטיות במוליכים למחצה באמצעות מודלים היפרבוליים של שתי טמפרטורות, והדגים את חשיבות אפקטי ההרפיה התרמית בחיזוי מדויק של התנהגות המערכת. התגובה התרמומכנית של קומפוזיטים מחוזקים בסיבים זכתה גם היא לתשומת לב רבה בשל חוזקם המכני המוגבר ותכונותיהם האנזוטרופיות. לי ולמברוס² ניתחו את ההתנהגות התרמומכנית הדינמית של קומפוזיטים כאלה, והדגישו את התאמתם ליישומים הנדסיים מתקדמים. באופן דומה, Kalkal ואח^' חקרו עיוותים דו-ממדיים (2D) במדיה סיבית מדורגת פונקציונלית תחת שדות מגנטיים, והראו את ההשפעה החזקה של אניזוטרופיה והשפעות חיצוניות על תגובת המערכת. יתרה מזאת, Pitarresi ^{ואחרים.4} בחנו את תפקיד ההטרוגניות המאקרוסקופית בהתנהגות תרמואלסטית, והדגישו את הצורך במידול מדויק של חומרים מרוכבים. פיתוחים אחרונים הרחיבו את המחקרים הללו למדיה מוליכים למחצה עם אפקטים מולטיפיזיקליים משולבים. מונדל ואח^' חקרו התפשטות גלים במוליכים למחצה מחוזקים תוך התחשבות בתגובות זיכרון ושדות מגנטיים, וחשפו אינטראקציות מקושרות מורכבות. מחקרים ניסיוניים, כמו אלו של Akai ואחרים^, העריכו נזק עייפות בחומרים מחוזקים בסיבים באמצעות שינויים בטמפרטורה תרמואלסטית, מה שמאשר את חשיבותם המעשית. בנוסף, נעשה שימוש בגישות מיקרומכניות להערכת תכונות תרמואלסטיות אפקטיביות, כפי שמודגם על ידי לו ואח^{' 7}. השפעת העומס החיצוני ותנאי הסביבה נחקרה גם היא באופן נרחב. ברק ודנקר⁸ ניתחו עומס מוטה במדיומים מחוזקים בסיבים מדורגים פונקציונלית, בעוד קונדו וקאלקל⁹ בחנו אינטראקציות פוטותרמיות תחת עומסים תרמיים בכוח הכבידה ובתנועה. צ'אודרי ואח^{' חקרו} תכונות תלויות טמפרטורה באמצעות מודלים של השהיה בפאזה כפולה, ופנדיט ואח^' יישמו מודלים של עיוות מסדר חלקי כדי ללכוד התנהגות עיוותים לא מקומית. מחקרים אלו מדגישים את החשיבות של התחשבות בתנאי עומס וחומרים ריאליסטיים בניתוח תרמואלסטי.

תופעות דינמיות וקשורות לרעידות נחקרו גם במערכות מוליכים למחצה. סונג ואח^{' חקרו} רטטים פוטותרמיים במבני מוליכים למחצה, בעוד מונדל וסור¹³ ניתחו התפשטות גלים במדיה אורתוטרופית עם השפעות זיכרון. השפעות ויסקואלסטיות ומיקרוסטרוקטוריות נבדקו על ידי אבואלרגל ואח' ¹⁴, והשפעות חימום מסוג רמפה נחקרו על ידי הוביני ואח' ¹⁵. יתרה מזאת, מודלים תרמואלסטיים מוכללים כמו תורת השהיית שלוש פאזות פותחו על ידי זנקור¹⁶ לשיפור דיוק התחזית. מודלים מולטיפיזיים מתקדמים המשלבים אפקטים אלקטרומגנטיים ומיקרופולריים העמיקו עוד יותר את ההבנה של התנהגות תרמו-פוטואלסטית. אל-חזאם ואח' ¹⁷ חקרו עירור פוטו-אלקטרו-מגנטו-תרמואלסטי במדיות מוליכים למחצה מסתובבות, בעוד שנאזיר וקומאר¹⁸ ניתחו אינטראקציות תרמואלסטיות מיקרופולריות. סונג ואח' ^בחנו גם אינטראקציות תרמואלסטיות לא מפזרות, ונסר ואבואלרגל²⁰ חקרו תהליכי ספיגת אור במוליכים למחצה עם חללים. מודלים לא-מקומיים ומודלים בסדר חלקי מילאו תפקיד מרכזי בהתפתחויות האחרונות. גופטה ואח^' חקרו עירור פוטותרמי במדיה לא מקומית נקבובית, בעוד הוביני ועבאס²² ניתחו התפשטות גלים בסדר שברירי במוליכים למחצה. חאפד וזנקור²³ חקרו השפעות עומס מוטה, ואוליניק ואח' ²⁴ בחנו השפעות תרמו-פוטואלסטיות לא סטציונריות. התנהגות מוליכים למחצה מדורגת פונקציונלית תחת עירור לייזר נחקרה גם היא על ידי Awwad ואח' ²⁵, בעוד גופטה ואחרים ^חקרו קישור תרמו-פייזו-פוטואלקטרי עם מודלים תלויי זיכרון. גישות ניסיוניות קלאסיות המשלבות טכניקות תרמיות ואופטיות הוקמו על ידי גרין ופטרסון²⁷ וברון ופטרסון²⁸, ומספקות שיטות אמינות לניתוח מאמצים. יתרה מזאת, קאור וסינג²⁹ פיתחו מודלים תלויי זיכרון לא מקומיים לתהודות מוליכים למחצה, ועבאס ואח^{' ניתחו} אינטראקציות פוטותרמיות עם מוליכות תרמית משתנה. מחקרים ניסיוניים ומספריים של ממברנות מרוכבות מחוזקות נערכו על ידי Lu ואח' ³¹, בעוד שפורקייט וקנורי³² חקרו תגובות זיכרון במדיות מחוזקות סיבים מסתובבות. אבו-דהאב ואח^' בחנו בהמשך את השתקפות הגלים במדיה תרמואלסטית מחוזקת בסיבים תחת תנאי עומס. לאחרונה הוצעו מודלים מתקדמים של תרמואלסטיים חלקיים ולא מקומיים לתיאור חומרים מורכבים. אבואלרגל ואח^' חקרו תגובות תרמיות ברקמות ביולוגיות באמצעות מודלים חלקיים, בעוד שסלוואמני ^{ואח' 35,36} חקרו התפשטות גלים לא-מקומית והתנהגות רטט בננו-קרניים. בנוסף, מודלים של השהיה פאזה כפולה ווויסקואלסטי יושמו על מיקרו-מבנים על ידי אבואלרגל ואח' ³⁷, ופותחו פורמולציות תרמואלסטיות חלקיות עם גרעיני זיכרון^38,39 כדי ללכוד תופעות מקושרות מורכבות.

למרות ההתפתחויות הנרחבות הללו, רוב המחקרים הקיימים מתמקדים בעיקר בפתרונות אנליטיים או נומריים מבלי לספק הדמיה מרחבית-זמנית מפורטת של השדות הפיזיקליים. ביישומים מעשיים רבים, במיוחד במדיות למחצה מחוזקות בסיבים אניזוטרופיים, תגובת המערכת תלויה מאוד הן בשינויים מרחביים והן בזמן, מה שהופך ויזואליזציה לחיונית לפרשנות מדויקת. יתרה מזאת, אף על פי שגישות נומריות וניסיוניות מספקות תובנות חשובות, שיטות אנליטיות המבוססות על טכניקות מצב נורמלי וערכי עצמי מציעות יתרונות משמעותיים בבעיות הכוללות תחומים חצי-אינסופיים ומערכות מולטיפיזיקליות מקושרות. גישות אלו מאפשרות פתרונות בצורה סגורה, המספקות הבנה פיזיקלית עמוקה יותר של מנגנוני התפשטות גלים, החלשה וקישור, ומשמשות כמדדים אמינים לאימות תוצאות מספריות וניסיוניות.

בעבודה הנוכחית מוצגת חקירה מקיפה של ההתנהגות התרמו-פוטואלסטית המקושרת במדיום מוליך למחצה מחוזק על ידי סיבים אניזוטרופיים. החידוש במחקר זה טמון בשילוב גישה אנליטית מבוססת ערכים עצמיים עם ויזואליזציה של מפת חום מרחבית-זמנית, כדי לספק תובנות עמוקות יותר להתפתחות וללוקליזציה של שדות פיזיקליים. מטרת עבודה זו היא לנתח את השפעת האניזוטרופיה וחיזוק הסיבים על האינטראקציה בין שדות תרמיים, מכניים ושדות נשאי, ולהדגים את היישום של השיטה המוצעת בפרשנות תופעות מורכבות הרלוונטיות ליישומים הנדסיים מודרניים כגון מכשירי מוליכים למחצה, חיישנים וטכנולוגיות מבוססות לייזר.

פרוטוקול

מחקר זה מבוסס כולו על מודלים תיאורטיים וסימולציות נומריות ואינו כולל משתתפים אנושיים, נבדקים בעלי חיים או דגימות ביולוגיות. לכן, לא נדרשו אישור אתי והסכמה מדעת.

ניסוח מתמטי של פוטו-תרמואלסטיות במדיה אנאיזוטרופית מחוזקת על ידי סיבים
המחקר הנוכחי בחן חצי-מרחב מוליך אניזוטרופי מחוזק בסיבים דו-ממדיים שנחשף לעוררות אופטית על פני השטח. המדיום תפס את האזור x ≥ 0, כאשר הגבול ב-x = 0 מייצג את המשטח החשוף. מערכת הקואורדינטות הוגדרה כך שציר ה-x מתפשט לתוך התווך, בעוד שציר ה-y נמצא לאורך פני השטח ותיאר את ההתנהגות במישור. החומר נחשב הומוגני אך אניזוטרופי בשל נוכחות סיבי חיזוק מיושרים, שיצרו תלות כיוונית בתכונות האלסטיות והקישור. ספיגת אופטית על פני השטח יצרה נשאי חום מקומיים ונשיאי מטען עודפים, מה שהוביל לאינטראקציה משולבת מלאה בין שדות תרמיים, מכניים ונשאים. בהתאם לכך, מצב המערכת תואר על ידי הטמפרטורה θ(x, y, t) (K), צפיפות הנשאים N (x, y, t) (^m-3), ורכיבי ההעתקה u (x, y, t) ו-v (x, y, t)(m), בהנחה של עיוותים קטנים. סכימה של התחום הפיזיקלי, מערכת הקואורדינטות, כיוון הסיבים והעירור האופטי המיושם מוצגת באיור 1. כל החישובים הסימבוליים והמספריים בוצעו באמצעות Wolfram Mathematica (גרסה 12.0).

איור 1. ייצוג סכמטי של מדיום מוליך למחצה מחוזק בסיבים חצי-אינסופי הנתון לעוררות אופטית בגבול x = 0. מערכת הקואורדינטות (x, y) מוצגת, כאשר כיוון הסיבים מיושר לאורך כיוון ה-x (a = (1, 0)), וממחיש את התצורה הגיאומטרית ואת האניזוטרופיה התלויה בכיוון של התווך. אנא לחצו כאן כדי לצפות בגרסה מוגדלת של הדמות הזו.

הקשר המרכיב לטנזור המתח בתווך למחצה תרמו-אלסטי אניזוטרופי מחוזק בסיבים הובע בצורה הכללית באמצעות משוואה 1, ^1,5. בניסוח זה, θ מציין את עליית הטמפרטורה ביחס לטמפרטורת הייחוס T₀, בעוד ש-T מייצג את הטמפרטורה המוחלטת כאשר הדבר רלוונטי.

.   (1)

כאן, _{C ijkl} הם מקדמי הקשיחות האלסטית, e_kl הוא טנזור העיוות, _{ו-β ij} ו-η_ij מייצגים את טנזורי התרמואלסטי והקישור הנשא, בהתאמה. בנוכחות חיזוק סיבים, תגובת החומר הפכה תלויה בכיוון ונשלטה על ידי וקטור כיוון הסיב a = (a_i), שהכניס תרומות אניזוטרופיות הן למונחי האלסטיות והן למונחי הקישור. בהתאם לכך, היחס המרכיב הורחב כדי לשלב במפורש את השפעת חיזוק הסיבים^כ-2,3:

. (2)

כאן, λ ו-μ_τ הם קבועי למה, _{ו-μ L} הוא מודולוס הגזירה האורכי לאורך כיוון הסיב. הפרמטר α מייצג השפעות חיזוק סיבים והוא שונה_{מ-α ij}, שמסמן מקדמי התפשטות תרמית. וקטור היחידה הגדיר את כיוון הסיבים והכניס תלות כיוונית בתגובת המתח-עיוות. בניסוח הדו-ממדי הנוכחי, הניחו שהסיבים מיושרים לאורך ציר ה-x; לכן, וקטור הכיוון נלקח במפורש כ-a = (1, 0). מפרט זה סיפק פרמטריזציה ברורה של כיוון הסיב והבטיח שהתרומות האניזוטרופיות ישולבו בעקביות במשוואות המנהלות, תוך התייחסות ישירה להתנהגות הכיוונית הנגרמת על ידי חיזוק הסיבים. בתצורה דו-ממדית הנוכחית, רכיבי המאמץ השולטים מצטמצמים ל:

, (3)

, (4)

. (5)

משוואות אלו ממחישות את ההשפעה המשולבת של אניזוטרופיה, חיזוק סיבים ואפקטי קישור מולטיפיזיקלי. המקדמים β_ij ו-η_ij הוגדרו לפי פרמטרי החומר כדלקמן:

,

,

,

.

כאן, המקדמים A_ij מייצגים את הקבועים האלסטיים האפקטיביים של התווך האנסוטרופי המחוזק בסיבים והוגדרו כך:

. (6)

כאן, λ, μ_L ו-μ_T הם הקבועים האלסטיים של התווך האניזוטרופי המחוזק על ידי סיבים, בעוד α_ij ו-ξ_ij מייצגים את מקדמי ההתפשטות התרמית והתפשטות הנשא, בהתאמה. התפשטות גלים אלסטיים במדיה למחצה תרמו-פוטואלסטית נשלטה על ידי עקרון שימור התנע הליניארי, שהיווה את הבסיס לניתוח תרמואלסטי דינמי. בהיעדר כוחות גוף, משוואת התנועה הכללית עבור רצף ניתן לעיוות מבוטאת כך, בהתבסס על ^1,15:

. (7)

כאן, ρ היא צפיפות המסה _{ו-σ ij} הוא טנזור המאמץ. במחקר הנוכחי, הניסוח הוגבל לתצורה דו-ממדית במישור x - y, ושדה ההזזה יוצג על ידי u(x, y, t) ו-v(x, y, t). בהתאם לניסוחים סטנדרטיים במדיה תרמו-פוטואלסטית, משוואות התנועה השולטות בשני ממדים נכתבו כך:

, (8)

. (9)

בהחלפת יחסי הקונסטיטוטיב האניזוטרופיים המחוזקים על ידי סיבים במשוואות הנ"ל, התקבלה מערכת משוואות דיפרנציאליות חלקיות (DEs) המקושרת שהתקבלה כך:

, (10)

. (11)

כאן, תת-כתובות מציינות התבגרות חלקית ביחס למשתנים מרחביים וזמניים. משוואות אלו מדגישות את ההשפעה המשולבת של אניזוטרופיה, חיזוק סיבים, גרדיאנטים של טמפרטורה ודיפוזיית נשאים על התגובה הדינמית של התווך. בנוכחות עירור אופטי, השדה התרמי בתוך המוליך למחצה הושפע מאוד מהאינטראקציה עם צפיפות הנשא ועיוות מכני, מה שהוביל לתהליך העברת אנרגיה מקושר במלואו. בניגוד להולכת חום קלאסית, התפתחות הטמפרטורה במדיה כזו נשלטה על ידי מונחי מקור נוספים שנבעו מרקומבינציה של נשאים ואפקטים תרמואלסטיים, אשר שינו משמעותית את מאפייני הפצת החום. משוואת הולכת החום במסגרת התרמואלסטיות המוכללת הובאה כך: ^16,20:

. (12)

כאן, C_E הוא החום הסגולי בעיוות קבוע, המייצג את הקיבולת התרמית של החומר, ו-T₀ מציין את הטמפרטורה המוחלטת של התווך במצב שיווי המשקל שלו. בתצורה דו-ממדית הנוכחית, משוואה זו הצטמצמה^ל-16,20:

. (13)

משוואה זו מראה ששדה הטמפרטורה הושפע לא רק ממוליכות תרמית כיוונית אלא גם מרקומבינציה של נשאים דרך האיבר , וכן מעיוות תלוי זמן דרך מונחי הקישור התרמואלסטי. ניסוח זה תפס את האינטראקציות הרב-פיזיקליות החיוניות השולטות בהעברת חום במוליכים למחצה מחוזקים בסיבים אניזוטרופיים והדגיש את תפקיד הדינמיקה של הנשאים והתגובה המכאנית בשינוי ההתנהגות התרמית של המערכת. כאשר מדיום מוליך למחצה היה חשוף לגירוי אופטי, נוצרו מספר משמעותי של נשאי מטען בעקבות ספיגת הקרינה הפוגעת. נשאים אלו עברו תהליכי הובלה שכללו דיפוזיה מרחבית, רקומבינציה ויצירה מונעת תרמית, שכולם היו קשורים באופן מהותי לשדה הטמפרטורה בתוך החומר. כתוצאה מכך, צפיפות הנשאים הפכה לאחד המשתנים המרכזיים השולטים בתגובה התרמו-פוטואלסטית המקושרת.

בנוסחה הנוכחית, התפתחות ריכוז הנשא תוארה באמצעות איזון בין מנגנוני דיפוזיה, השפעות דעיכה ותהליכי הפעלה תרמית, מה שהוביל ליחס השלטון הבא ^1,5"

. (14)

כאן, D_E מייצג את מקדם הדיפוזיה של הנשא והוא אופרטור לפלסיאן דו-ממדי במישור x - y. המונח מתייחס להשפעות רקומבינציה עם זמן הרפיה τ, בעוד ש-k הוא מקדם הקישור התרמו-נשא המוגדר כ-, שמאפיין את רגישות ריכוז נושא שיווי המשקל N0 לשינויים בטמפרטורה. קשר זה מדגיש את תפקיד הטמפרטורה כמנגנון מניע ליצירת נשאים ומבסס חיבור ישיר בין השדות התרמיים והאלקטרוניים במדיום המוליך למחצה המחוזק בסיבים אניזוטרופיים.

הוקמו המשוואות השולטות והניסוח המתמטי של מערכת הנשאים הפוטתרמואלסטית המקושרת. הפרמטרים הפיזיקליים והחומריים המתאימים למדיום הסיליקון (Si) מסוכמים בטבלה 1, יחד עם הערכים המספריים, היחידות וההפניות המתאימות שלהם. פרמטרים אלו משמשים לאחר מכן בחישובים נומריים ובתהליך הלא-ממדיזציה.

סמל	ערך	יחידה	מקור
λ	3.64 × 10¹⁰	N/m²	12
μT	5.46 × 10¹⁰	N/m²	12
μL	3.20 × 10¹⁰	N/m²	12
ρ	2330	ק"ג/מ³	13
CE	695	J/(kg· K)	30
K11	0.0921 × 10³	עם (m·K)	30
K22	0.0963 × 10³	עם (m·K)	30
DE	2.5 × 10⁻³	מ²/ש	22
τ	5 × 10⁻⁵	s	15
T₀	300	K	15
לדוגמה	1.11 × 10⁻¹⁹	J	12
α11	3.1 × 10⁻⁶	ק⁻¹	30
α22	3.5 × 10⁻⁶	ק⁻¹	30
ξ11	−7 × 10⁻³¹	מ³	21
ξ22	−9 × 10⁻³¹	מ³	21
κ	2.16 × 10²¹	מ⁻³·ס⁻¹· ק⁻¹	21
α	−1.28 × 10¹⁰	N/m²	28
β	220.90 × 10¹⁰	N/m²	28
ω	2.95 + 1i	ס⁻¹	12
a	1	— (ללא מימד)	13
y	0.6	m	13
θ₀	1	— (ללא מימד)	15
N₀	1	— (ללא מימד)	15

טבלה 1. תכונות החומר והפרמטרים המשמשים בניתוח נומרי של מדיום מוליך למחצה מחוזק בסיבים אניזוטרופיים. כל הכמויות מבוטאות ביחידות SI אלא אם צוין אחרת. פרמטרים חסרי ממד מסומנים בהתאם. הערכים המפורטים מתאימים לתכונות חומר מבוססות סיליקון ולפרמטרי המודל המשמשים בחישובים הנוכחיים, כפי שהתקבלו מהמקורות המצוטטים. מקדם הקישור התרמו-נושא κ מוגדר כ-κ = (∂N₀/∂T)(1/τ), בהתאם לניסוחים סטנדרטיים במודלים של מוליכים למחצה תרמו-פוטואלסטיים.

ניסוח לא-ממדי של מודל הפוטו-תרמואלסטי האנאיזוטרופי המקושר
כדי לפשט את המשוואות השולטות ולקבל ייצוג לא-ממדי עקבי של מערכת התרמו-פוטואלסטית המקושרת, הוכנסו סולמות מאפיינים מתאימים לקואורדינטות המרחביות x,y, זמן t, רכיבי הזזה u, v, טמפרטורה T, צפיפות נשאים N ומתח σ. פרמטרי קנה מידה אלו נבחרו בעקביות בהתבסס על התכונות הפיזיקליות הפנימיות של התווך ומנגנוני הקישור בין שדות תרמיים, מכניים ונשאים, בהתאם לניסוחים מבוססים שדווחו בספרות^16,21. בהתאם לכך, המשתנים חסרי הממד הוגדרו כך:

, , , , , , .

טרנספורמציה זו הפחיתה את מספר פרמטרי החומר הבלתי תלויים וסיפקה ייצוג מנורמל של המערכת המקושרת. על ידי החלפת המשתנים חסרי הממד לעיל במשוואות השלטון שנגזרו קודם, המערכת נכתבה מחדש בצורה לא-ממדית. לשם הפשטות, הסימון הראשוני הקשור למשתנים חסרי הממד הושמט לאחר מכן. הליך זה הניב קבוצה קומפקטית של DEs חלקיים חסרי ממד, שניתן לכתוב בצורתה הבאה:

, (15)

, (16)

, (17)

. (18)

לאחר יישום הטרנספורמציה הלא-ממדית, רכיבי המתח של המערכת נכתבו בצורה המנורמלת הבאה:

, (19)

, (20)

. (21)

הפרמטרים חסרי הממד a_i הוצגו כדי לייצג שילובים קומפקטיים של תכונות פיזיקליות וחומריות השולטות בהתנהגות הפוטו-תרמואלסטית האנאיזוטרופית המקושרת. כל מקדם שיקף מנגנון אינטראקציה ספציפי בתוך המערכת וסיפק תובנות לגבי ההשפעה היחסית של התהליכים הפיזיקליים הבסיסיים. מייצג את היחס בין קשיחות הקישור הנורמלית לבין הקשיחות האלסטית העיקרית, ומשקף את מידת האינטראקציה האניזוטרופית בין שני רכיבי ההזזה. מאפיין את התרומה היחסית של עיוות רוחבי לרכיב המאמץ הנורמלי. מודד את השונות הכיוונית של קישור תרמואלסטי, ומצביע על אניזוטרופיה בהשפעות התפשטות תרמית. מתאר את ההשפעה האניזוטרופית של צפיפות הנשאים על העיוות האלסטי המושרה. מייצג את קשיחות הגזירה המנורמלת ומכמת את תרומת עיוות הגזירה ביחס לעיוות נורמלי. מסביר את הקישור המשולב בין עיוות נורמלי לעיוות גזירה במשוואות ההזזה השולטות. מבטא את היחס בין קשיחות רוחבית לקשיחות גזירה, ומדגיש התנהגות עיוותים אניזוטרופית. מייצג את הפרמטר האינרציאלי המנורמל, המקשר בין השפעות התפשטות גלים לקשיחות גזירה. מאפיין את הקישור בין גרדיאנטים של הזזה בכיוונים מרחביים שונים. מכמת את התרומה היחסית של השפעות תרמיות לשדה ההעתקה בכיוון הרוחבי. מודד את השפעת העיוות המושרה על ידי הנשאים ביחס לקשיחות גזירה. : מייצג את האניזוטרופיה במוליכות תרמית בכיוונים מרחביים שונים. מאפיין את השפעת רקומבינציה של נשאים על יצירת חום בתוך המדיום. מייצג את הקישור בין אפקטים תרמיים לעיוות אלסטי תלוי זמן. מסביר את ההשפעה המשולבת של התפשטות תרמית אניזוטרופית בשני הכיוונים המרחביים. מייצג את פרמטר הדיפוזיה המנורמל השולט בקצב העברת הנשאים. מאפיין את העוצמה היחסית של השפעות הרקומבינציה של נשאי. מתאר את הקישור בין שינויים תרמיים לתהליכי יצירת נשאי.

פתרון אנליטי באמצעות טכניקת מצב נורמלי
כדי לקבל פתרונות אנליטיים למערכת תרמו-פוטואלסטית אניזוטרופית מקושרת, שימשה טכניקת מצב נורמלי בשל יעילותה בצמצום ה-DEs החלקיים השולטים למערכת נוחה יותר של DE רגילים. גישה זו משמשת באופן נרחב בניתוח תופעות התפשטות גלים, כולל פיזור והחלשה. בהתאם לכך, הונחו שינויים הרמוניים במשתני השדה הן בזמן והן בכיוון המרחב הרוחבי ^1,12,23. לכן, רכיבי ההזזה, הטמפרטורה, צפיפות הנשאים והמאמץ הובאו לידי ביטוי בצורה מעריכית כך:

. (22)

כאן, ω מציין את התדר המרוכב השולט בהתנהגות הזמן של השדות, בעוד ש-a מייצג את מספר הגל הקשור לשונות המרחבית לאורך כיוון y. פרמטרים אלו נבחרו כדי לעמוד בדרישות היציבות ולהבטיח פתרונות חסומים פיזיקליים בקביל בתוך התחום החצי-אינסופי. על ידי החלפת הצורות המניחות לעיל למשוואות השלטון הלא-ממדיות שנגזרו קודם ופישוט הביטויים המתקבלים, המערכת המקורית המקושרת של DEs חלקיים צומצמה למערכת של DEs רגילים ביחס לקואורדינטות המרחביות , שניתן לכתוב כך:

, (23)

, (24)

, (25)

. (26)

יתרה מזאת, רכיבי המאמץ המתאימים בתחום הטרנספורמציה נכתבו כך:

, (27)

, (28)

. (29)

כאן, D מציין את האופרטור הדיפרנציאלי . משוואות אלו מייצגות את הצורה המופחתת של מערכת הממשל בתחום המצב הנורמלי ומספקות את הבסיס לגזירת משוואת האופייניות ולבניית הפתרון האנליטי הכללי בשלבים הבאים. המקדמים הוגדרו כך: , , , , , .

ניסוח DE מטריצתי וניתוח ערכים עצמיים
בעקבות יישום טרנספורמציית המצב הנורמלי, מערכת הממשל שניתנה במשוואות 23–26 צומצמה לקבוצת DEs רגילים מסדר שני ביחס לקואורדינטות המרחביות. כדי להקל על פתרון שיטתי, מערכת זו הוסבה למערכת מסדר ראשון שקולה על ידי הכנסת משתנים עזר התואמים לנגזרות הראשונות של כמויות השדה. באופן ספציפי, הוגדרו המשתנים הבאים:

, . (30)

באמצעות הגדרות אלו, משוואות 23–26 נכתבו מחדש כמערכת הבאה של שמונה DEs מסדר ראשון:

, (31)

, (32)

, (33)

, (34)

. (35)

המערכת לעיל הובאה בצורת מטריצה קומפקטית A כך:

. (36)

וקטור המצב ניתן כך:

. (37)

ומטריצת המערכת קיבלה את הצורה המפורשת:

. (38)

ניסוח זה הפך את המערכת המקורית לבעיה עם ערכים עצמיים ^1,15. משוואת האופייניות התקבלה מ-

. (39)

מה שמניב פולינום מסדר שמינית השולט בערכים העצמיים. בצורה מצומצמת, ניתן לכתוב את הפולינום האופייני כך

. (40)

כאשר Z_i המקדמים הם פונקציות של פרמטרי המערכת ומוגדרים במפורש להלן. הערכים העצמיים המתקבלים קובעים את ההתנהגות המרחבית של הפתרון, כולל מאפייני החלשה והפצה. רק ערכי עצמי שמקיימים Re(m) > 0 נשמרים כדי להבטיח פתרונות פיזיקליים שמתדרדרים באופן אקספוננציאלי עם x → ∞.

. (41)

שורשי הפולינום האופייני מגדירים את הערכים העצמיים m, השולטים בהתנהגות המרחבית של הפתרון. ערכים עצמיים אלו חושבו נומרית באמצעות Mathematica על ידי בניית הפולינום המאפיין באמצעות פונקציית CharacteristicPolynomial ופתרון המשוואה האלגברית המתקבלת באמצעות NSolve. מכיוון שהבעיה מנוסחת בתחום חצי-אינסופי (x ≥ 0), נבחנים רק פתרונות פיזיקליים שנשארים חסומים כ-x → ∞. בהתאם לכך, נשמרו רק ערכי עצמי שמקיימים Re(m) > 0, מה שמבטיח דעיכה מעריכית של פתרונות בצורה exp(−mx) כ-x → ∞. השורשים הנותרים הושלכו כיוון שהם מתאימים לפתרונות שאינם מתפרקים או לא חסומים שאינם תואמים לדרישות הפיזיקליות של המודל.

עבור כל ערך עצמי שנשמר m, הושג וקטור עצמי מתאים מהמערכת האלגברית המתאימה

, (42)

והוא הובע בצורה הבאה:

. (43)

בהרחבת משוואת המטריצות לעיל, התקבלה מערכת המשוואות הליניאריות הבאה:

, (44)

, (45)

, (46)

, (47)

. (48)

עקב ההומוגניות של בעיית הערכים העצמיים, הווקטורים העצמיים הוגדרו עד לקבוע כפל שרירותי. כדי לקבל ייצוג ייחודי ועקבי, הוטל תנאי נרמול על ידי קיבוע רכיב אחד של וקטור העצמי. בעבודה הנוכחית, הרכיב הראשון נבחר כך ש-q₁ = 1, והרכיבים הנותרים נקבעו ברצף מתוך מערכת המשוואות הנ"ל. מבחינה חישובית, נורמליזציה זו יושמה על ידי הקצאת ערך יחידה לרכיב אחד ופתרון מערכת המשוואות הליניאריות שנוצרה כדי להעריך את שאר הרכיבים. הליך זה סיפק דרך שיטתית וניתנת לשחזור לחישוב הווקטורים העצמיים הקשורים לכל ערך עצמי קביל.

. (49)

והרכיבים הנותרים נובעים בהתאם מיחסי המערכת. וקטורים עצמיים אלו מתארים את התרומה היחסית של טמפרטורה, צפיפות נשאים ושדות ההזזה בתוך כל מצב. כתוצאה מכך, הפתרון הכללי של הבעיה נבנה כשילוב ליניארי של המצבים העצמיים המותרים, שכל אחד מהם קשור לערך עצמי ולוקטור העצמי המתאים לו, ובכך סיפק תיאור אנליטי מלא של ההתנהגות הפוטו-תרמואלסטית האנאיזוטרופית המקושרת במדיום חצי-מרחב. לכן, הפתרון הכללי של המערכת נכתב כך:

. (50)

כאן, C_i הם קבועים שנקבעו מתנאי הגבול. על ידי הרחבת הביטוי הווקטורי הנ"ל, משתני השדה התקבלו כך:

, (51)

, (52)

, (53)

. (54)

ייצוג זה מראה שהפתרון מורכב מסופרפוזיציה של מצבים אקספוננציאליים, כאשר כל זוג ערך עצמי-וקטור עצמי תורם באופן עצמאי לתגובה הפיזיקלית הכוללת. הערכים העצמיים המותרים נבחרים כך שהחלקים הממשיים שלהם חיוביים, מה שמבטיח פתרונות חסומים ובעלי משמעות פיזיקלית כ-x → ∞.

תנאי גבול ואילוצים פיזיקליים
על ידי החלפת הפתרון הכללי בתנאי הגבול שנקבעו ב-x = 0, התקבלה מערכת של משוואות אלגבריות ליניאריות במונחי הקבועים C_i. באופן ספציפי, כל תנאי שפה (טמפרטורה, צפיפות נשאים ואילוצי הזזה) הוצג במונחים של התרחבות המצב העצמי, מה שהוביל לקבוצת משוואות הקשורה למקדמים C_i. הליך זה הוביל למערכת ליניארית שניתן לכתוב אותה בצורת מטריצה כ-BC = D, כאשר B היא מטריצת המקדמים שנבנתה מרכיבי הווקטורים העצמיים המוערכים בגבול, C = (C₁, C₂, C_{3,C 4})^T הוא וקטור של קבועים לא ידועים, ונקבע מערכי הגבול המוטלים כגון θ₀, N₀, ומגבלות ההעתעה. המערכת הליניארית שהתקבלה נפתרה חישובית באמצעות Mathematica, כאשר מטריצת המקדמים והוקטור הימני הורכבו במפורש, והקבועים הלא ידועים הושגו באמצעות שגרת LinearSolve. קבועים אלו הוחלפו חזרה לפתרון הכללי כדי לבנות את הביטויים השלמים לשדות הפיזיקליים, אשר שימשו לאחר מכן בהערכה נומרית ובייצוג גרפי של התוצאות.

תנאי הגבול שהוטלו היו כדלקמן:

הגבלת טמפרטורה:

. (55)

תנאי זה מייצג טמפרטורת פני שטח משתנה הרמונית הנגרמת על ידי חימום אופטי תקופתי. הוא פועל כגירוי תרמי ראשוני שמניע את תהליכי ההובלה התרמואלסטית והנשיאה המקושרים בתוך התווך. המשרעת θ₀ מאפיינת את עוצמת העומס התרמי המופעל.

מגבלת צפיפות הנשאים:

. (56)

תנאי גבול זה מתאר את צפיפות הנשאים שנוצרת על ידי הפוט הנובעת מהאורה האופטית. הוא משקף את ההתרגשות האלקטרונית הנובעת מספיעת פוטונים ואת המודולציה ההרמונית שלה התואמת את השדה האופטי הנכנס.

מגבלת הזזה:

. (57)

תנאי זה מציין שהגבול מוגבל מכנית בכיוון הרוחבי. לכן, אין הזזה לאורך כיוון ה-v על פני השטח.

מגבלת מאמץ גזירה:

. (58)

תנאי זה מתאים לגבול חופשי מאחיזה ביחס למאמץ גזירה. הוא מבטיח שאין כוחות משיקים לפעול על פני השטח, וזה תואם לגבול מכני חופשי בכיוון המשיק. בנוסף לתנאי הגבול ב-x = 0, הדרישה הפיזיקלית באינסוף הוטלה כך: הבטחת פתרונות פיזיקליים חסומים בתוך התחום החצי-אינסופי. לפני הצגת התוצאות המספריות, ההליך החישובי הכולל שאומץ במחקר זה מסוכם באיור 2. הערכים המספריים של פרמטרי העירור θ₀, N₀, התדר המרוכב ω ומספר הגל a המשמשים בחישובים מפורטים בטבלה 1. הפרמטרים המפורטים בטבלה 1 כוללים גם קבועי חומר ממדיים וגם פרמטרים לא-ממדיים המשמשים בניסוח המנורמל. להערכה נומרית, התחום המרחבי הוגדר כך , הקואורדינטה הרוחבית נקבעה ב-y = 0.6, והתחום הזמני נלקח בחשבון בתוך . טווחים אלו שימשו לכל החישובים המספריים והייצוגים הגרפיים.

איור 2. זרימת עבודה חישובית של השיטה המוצעת. האיור ממחיש את רצף השלבים מהניסוח לתוצאות נומריות: משוואות שולטות, אי-ממדיזציה, יישום טכניקת המצב הנורמלי, המרה למערכת מסדר ראשון, ניסוח מטריצות, ניתוח ערך עצמי וקטורים עצמיים, יישום תנאי גבול, קביעת קבועים ויצירת גרפים נומריים. אנא לחצו כאן כדי לצפות בגרסה מוגדלת של הדמות הזו.

תוצאות

תוצאות מספריות
בחלק זה בוצעו חישובים נומריים לניתוח ההתנהגות של מערכת הנושא התרמו-פוטואלסטי המקושרת בתווך מוליך למחצה מחוזק בסיבים אניזוטרופיים. החומר שנבחן היה סיליקון (Si), והפרמטרים הפיזיים והחומריים שלו מפורטים בטבלה 1. קבועי החומר הללו הוחלפו ישירות במשוואות השולטות ויושמו בחישובים מספריים להערכת משתני השדה. כל הפרמטרים שולבו בעקביות באמצעות קנה מידה לא-ממדי מוגדר. החישובים המספריים בוצעו באמצעות הגדרות הדיוק וההתכנסות המוגדרות כברירת מחדל ב-...

דיון

התוצאות שהושגו מספקות תובנה פיזיקלית ברורה לגבי ההתנהגות התרמו-פוטואלסטית המקושרת במדיה למחצה מחוזקת בסיבים אניזוטרופיים. המחקר הנוכחי מציע מסגרת אנליטית מבוססת ערכים עצמיים לחקירת האינטראקציה בין עומס תרמי, יצירת נשאים ועיוות אלסטי במדיה כזו. התגובה הנצפית נשלטת באופן יסודי על ידי הקישור החזק בין התהליכים הפיזיקליים הללו. ספיגת האנרגיה האופטית בגבול מובילה לחימום מקומי ולעירור נשאי, אשר בתורם גורמים למאמצים והזזות תרמואלסטיות בתוך התווך. מנגנון הקישור הזה מסביר את ההתפתחות הסימולטנית של...

חומרים

רשימת החומרים שנעשה בהם שימוש במאמר זה

שם	חברה	מספר קטלוג	הערות
תוכנה חישובית (ניתוח סימבולי ונומרי)	וולפרם ריסרץ'	נעשה שימוש ב-Wolfram Mathematica (גרסה 12.0)	משמש לחישוב ערכים עצמיים, יישום פתרונות אנליטיים והערכה נומרית
כלי ויזואליזציה של נתונים (יצירת קונטור ומפת חום)	וולפרם ריסרץ'	Wolfram Mathematica (גרסה 12.0) שימשה ליצירת גרפים דו-ממדיים ומפות חום מרחב-זמניות	משמש ליצירת גרפים דו-ממדיים ומפות חום מרחב-זמניות
מערך נתונים של פרמטרי חומר (תכונות סיליקון למחצה)	מקורות ספרותיים שונים	לא זמין	קבועים פיזיקליים (אלסטיים, תרמיים, קשורים לנשאים) המשמשים בחישובים (טבלה 1)
מחשב אישי/תחנת עבודה	HP	לא זמין	החישובים בוצעו על מחשב אישי סטנדרטי שרץ Windows OS עם זיכרון מספק לסימולציות מספריות
עורך משוואות	מיקרוסופט  Word ו-MathType	לא זמין	משמש לעיצוב והצגת ביטויים מתמטיים בכתב היד
תוכנת ניהול הפניות	אלסוויר	לא זמין	משמש לניהול הפניות ועיצוב ציטוטים (בסגנון ונקובר)