Tài liệu Chương 15: Các phương pháp số trong giải bài toán trường điện từ

Tìm hiểu lý thuyết trường điện từ trong Chapter 15 với các định luật cơ bản, phương trình sóng và bài tập vận dụng chi tiết môn vật lý.

• Chuyên ngành: Điện từ trường
• Đề tài: Phương pháp số
• Loại tài liệu: Chương sách

Tóm tắt

I. Tổng quan về Chapter 15 trường điện từ và phương pháp số

Chapter 15 trường điện từ giới thiệu các phương pháp số hiện đại trong việc phân tích trường. Đây là công cụ thiết yếu để giải quyết các bài toán kỹ thuật phức tạp. Kỹ thuật giải tích truyền thống cung cấp nghiệm dạng đóng chính xác. Tuy nhiên, giải tích chỉ áp dụng hiệu quả cho các cấu hình hình học đơn giản và lý tưởng hóa. Khi biên bài toán không khớp với các mặt tọa độ chuẩn, nghiệm giải tích trở nên bất khả thi. Sự phát triển mạnh mẽ của máy tính số thúc đẩy việc ứng dụng các phương pháp tính toán số. Các kỹ thuật này cung cấp nghiệm xấp xỉ với độ tin cậy cao và thỏa mãn hoàn toàn yêu cầu kỹ thuật thực tế. Phương pháp số không loại bỏ phương pháp giải tích mà đóng vai trò bổ trợ chặt chẽ. Mỗi phương pháp số đều thực hiện đơn giản hóa toán học để chuyển đổi phương trình vi phân hoặc tích phân thành hệ phương trình đại số tuyến tính.

1.1. Khái niệm và vai trò của phương pháp số trong điện từ học

Phương pháp số là kỹ thuật tính toán gần đúng trên máy tính nhằm giải các bài toán trường điện từ. Kỹ thuật này chuyển đổi các phương trình đạo hàm riêng và phương trình tích phân thành hệ phương trình đại số ma trận. Trong thực tế kỹ thuật, nhiều bài toán sở hữu hình học phức tạp và môi trường vật liệu không đồng nhất. Khi đó, việc tìm kiếm nghiệm giải tích là điều bất khả thi. Phương pháp số trở thành công cụ chủ đạo giúp xác định phân bố điện thế, điện trường và dòng điện. Độ chính xác của nghiệm số được kiểm soát dễ dàng thông qua việc tinh chỉnh kích thước lưới chia.

1.2. Mối quan hệ bổ trợ giữa nghiệm giải tích và nghiệm số

Nghiệm giải tích biểu diễn dưới dạng hàm toán học tường minh, giúp quan sát rõ sự biến thiên của trường theo các tham số. Tuy nhiên, phạm vi áp dụng giải tích bị giới hạn trong các bài toán đối xứng cao. Ngược lại, phương pháp số giải quyết linh hoạt mọi cấu trúc bất kỳ nhưng chỉ đưa ra giá trị số tại từng điểm cụ thể. Nghiệm giải tích đóng vai trò quan trọng làm chuẩn đối chiếu để đánh giá độ chính xác và kiểm tra tính hội tụ của các thuật toán số. Sự kết hợp giữa hai phương pháp tạo nên quy trình phân tích hoàn chỉnh.

II. Phân tích hạn chế giải tích trong Chapter 15 trường điện từ

Nhiều bài toán điện từ kinh điển gặp trở ngại lớn khi áp dụng giải tích thuần túy. Ví dụ tiêu biểu là tụ điện phẳng hữu hạn trong kỹ thuật. Phân tích giải tích thường bỏ qua hiệu ứng rìa ngoài mép bản cực. Giả định này chỉ đúng khi khoảng cách giữa hai bản cực rất nhỏ so với kích thước bề mặt. Trong thực tế, điện trường rìa làm thay đổi đáng kể điện dung và phân bố năng lượng. Phương trình Laplace và Poisson cũng đối mặt với giới hạn tương tự. Nghiệm giải tích đòi hỏi biên miền tính toán phải trùng khớp với hệ tọa độ chuẩn Đề-các, trụ hoặc cầu. Khi cấu trúc có biên cong bất kỳ, phương pháp tách biến thất bại hoàn toàn. Ngoài ra, môi trường truyền sóng không đồng nhất có hằng số điện môi biến đổi theo không gian khiến phương trình vi phân trở nên phức tạp. Do đó, việc số hóa miền tính toán là giải pháp bắt buộc.

2.1. Rào cản biên hình học phức tạp và hiệu ứng rìa vật lý

Các thiết bị thực tế trong kỹ thuật cao tần hiếm khi có biên phẳng lý tưởng. Mép kim loại sắc nhọn tạo ra mật độ điện tích tập trung cao, làm biến dạng trường xung quanh. Phương pháp giải tích buộc phải đưa ra nhiều giả thiết giản lược, dẫn đến sai số lớn khi tính toán thông số tán xạ và điện dung ký sinh. Sự sai lệch này ảnh hưởng trực tiếp đến hiệu năng của mạch tích hợp và anten. Việc mô hình hóa chính xác đường bao hình học thực tế là yêu cầu bắt buộc mà chỉ phương pháp số mới đáp ứng được.

2.2. Sự phức tạp của phương trình vi phân và tích phân điện từ

Các bài toán điện từ thường quy về hai dạng toán học chính là phương trình đạo hàm riêng và phương trình tích phân. Phương trình đạo hàm riêng mô tả trường cục bộ tại từng điểm trong không gian. Phương trình tích phân mô tả tương tác trường thông qua nguồn dòng và điện tích phân bố. Cả hai dạng toán này đều cực kỳ khó giải khi cấu trúc chứa nhiều lớp vật liệu khác nhau. Điều kiện biên phức tạp tại mặt phân cách môi trường khiến các công thức toán học giải tích tổng quát hoàn toàn không thể thiết lập.

III. Các giải pháp tính toán cốt lõi ở Chapter 15 trường điện từ

Chapter 15 trường điện từ cung cấp ba phương pháp số cốt lõi nhằm giải quyết triệt để các hạn chế trên. Phương pháp sai phân hữu hạn FDM xấp xỉ các đạo hàm riêng bằng tỉ số sai phân trên lưới điểm đều. Điện thế tại mỗi nút được xác định thông qua giá trị trung bình của các nút lân cận bằng kỹ thuật lặp. Phương pháp phần tử hữu hạn FEM vượt trội trong việc xử lý miền có biên cong phức tạp. Miền khảo sát được chia nhỏ thành các phần tử tam giác hoặc tứ diện với kích thước tùy biến. Hàm thế được xấp xỉ cục bộ trên từng phần tử, sau đó ghép nối thành hệ phương trình ma trận tổng thể. Phương pháp moment MoM chuyên dùng để giải các phương trình tích phân biên. Kỹ thuật này phân chia nguồn điện tích hoặc dòng điện trên bề mặt vật dẫn thành các đoạn nhỏ, giúp giảm số chiều bài toán và tính toán nhanh chóng bức xạ anten.

3.1. Phương pháp sai phân hữu hạn và kỹ thuật lặp ma trận

Phương pháp sai phân hữu hạn FDM thay thế phương trình đạo hàm riêng Laplace bằng sơ đồ sai phân năm điểm. Toàn bộ miền tính toán được chia thành lưới ô vuông hoặc chữ nhật. Giá trị điện thế tại từng nút nội bộ bằng trung bình cộng giá trị của bốn nút xung quanh. Kỹ thuật lặp Gauss-Seidel giúp cập nhật liên tục giá trị điện thế cho đến khi sai số đạt mức cho phép. Thuật toán FDM có cấu trúc đơn giản, dễ dàng lập trình và thực thi hiệu quả trên các phần mềm tính toán kỹ thuật.

3.2. Phương pháp phần tử hữu hạn và phương pháp moment

Phương pháp phần tử hữu hạn FEM chia nhỏ miền tính toán thành lưới tam giác linh hoạt, cho phép làm mịn lưới tại các khu vực trường biến thiên mạnh. Ma trận phần tử được xây dựng dựa trên nguyên lý biến phân và ghép nối thành hệ phương trình đại số. Trong khi đó, phương pháp moment MoM biến đổi phương trình tích phân thành phương trình ma trận bằng cách sử dụng hàm cơ sở và hàm kiểm tra. MoM đặc biệt tối ưu cho các bài toán phân tích bức xạ anten và tán xạ sóng điện từ.

IV. Ứng dụng thực tiễn của Chapter 15 trường điện từ hiện đại

Nền tảng lý thuyết trong Chapter 15 trường điện từ là cơ sở vận hành của hầu hết các phần mềm mô phỏng thương mại hiện đại. Các công cụ kỹ thuật hàng đầu như Ansys HFSS, CST Studio Suite và COMSOL Multiphysics đều phát triển dựa trên thuật toán FEM, FDM và MoM. Kỹ sư sử dụng các phần mềm này để thiết kế anten vi dải, ống dẫn sóng và vi mạch tích hợp cao tần. Phương pháp số giúp tính toán chính xác phân bố trường, xác định thông số trở kháng và tối ưu hóa hiệu suất bức xạ. Bên cạnh đó, mô phỏng số đóng vai trò then chốt trong việc phân tích độ tương thích điện từ EMC và kiểm soát can nhiễu EMI trên bo mạch tốc độ cao. Việc áp dụng mô phỏng số giúp phát hiện lỗi thiết kế sớm, giảm thiểu chi phí chế tạo mẫu thử nghiệm và đẩy nhanh tiến độ nghiên cứu sản phẩm.

4.1. Tích hợp thuật toán số trong phần mềm mô phỏng thương mại

Các phần mềm mô phỏng công nghiệp khai thác tối đa ưu thế của từng phương pháp số. Phương pháp FEM được ứng dụng chủ yếu trong phân tích vi ba và cộng hưởng khoang rỗng nhờ khả năng mô tả hình học chi tiết. Phương pháp MoM chiếm ưu thế trong mô phỏng mảng anten radar và bài toán tán xạ điện từ không gian mở. Quy trình tự động hóa chia lưới và giải ma trận lớn giúp kỹ sư thu được kết quả chính xác về đồ thị bức xạ và tham số truyền sóng trong thời gian ngắn.

4.2. Tối ưu hóa thiết kế phần cứng và tương thích điện từ

Mô phỏng số là giải pháp tiêu chuẩn để kiểm soát bức xạ và nâng cao tương thích điện từ trong các hệ thống điện tử. Trên bo mạch in mật độ cao, hiện tượng xuyên âm giữa các đường mạch dẫn gây suy hao tín hiệu nghiêm trọng. Phương pháp số hỗ trợ xác định chính xác phân bố dòng cảm ứng và vùng rò rỉ trường. Nhờ đó, đội ngũ kỹ thuật có thể tối ưu hóa vị trí linh kiện, thiết kế lớp che chắn kim loại và đường tiếp địa hiệu quả.

Trích đoạn nội dung tài liệu

Chapter 15 NUMERICAL METHODS The recipe for ignorance is: be satisfied with your opinions and content with your knowledge.1 INTRODUCTION In the preceding chapters we considered various analytic techniques for solving EM prob- lems and obtaining solutions in closed form. A closed form solution is one in the form of an explicit, algebraic equation in which values of the problem parameters can be substituted. Some of these analytic solutions were obtained assuming certain situations, thereby making the solutions applicable to those idealized situations. For example, in deriving the formula for calculating the capacitance of a parallel-plate capacitor, we assumed that the fringing effect was negligible and that the separation distance was very small compared with the width and length of the plates.

Also, our application of Laplace's equation in Chapter 6 was restricted to problems with boundaries coinciding with coordinate surfaces. Analytic solutions have an inherent advantage of being exact. They also make it easy to observe the behavior of the solution for variation in the problem parameters. However, an- alytic solutions are available only for problems with simple configurations.

When the complexities of theoretical formulas make analytic solution intractable, we resort to nonanalytic methods, which include (1) graphical methods, (2) experimental methods, (3) analog methods, and (4) numerical methods. Graphical, experimental, and analog methods are applicable to solving relatively few problems. Numerical methods have come into prominence and become more attractive with the advent of fast digital computers. The three most commonly used simple numerical techniques in EM are (1) moment method, (2) finite difference method, and (3) finite element method.

Most EM problems involve either partial differential equations or integral equations. Partial differ- ential equations are usually solved using the finite difference method or the finite element method; integral equations are solved conveniently using the moment method. Although numerical methods give approximate solutions, the solutions are sufficiently accurate for engineering purposes. We should not get the impression that analytic techniques are out- dated because of numerical methods; rather they are complementary.

As will be observed later, every numerical method involves an analytic simplification to the point where it is easy to apply the method. The Matlab codes developed for computer implementation of the concepts developed in this chapter are simplified and self-explanatory for instructional purposes.2 FIELD PLOTTING • 661 used in the programs are as close as possible to those used in the main text; some are defined wherever necessary. These programs are by no means unique; there are several ways of writing a computer program. Therefore, users may decide to modify the programs to suit their objectives.2 FIELD PLOTTING In Section 4.9, we used field lines and equipotential surfaces for visualizing an electrosta- tic field.

However, the graphical representations in Figure 4.21 for electrostatic fields and in Figures 7.16 for magnetostatic fields are very simple, trivial, and qualitative. Accurate pictures of more complicated charge distributions would be more helpful. In this section, a numerical technique that may be developed into an interactive computer program is presented. It generates data points for electric field lines and equipotential lines for arbitrary configuration of point sources.

Electric field lines and equipotential lines can be plotted for coplanar point sources with simple programs. Suppose we have N point charges located at position vectors r t , r 2 ,., rN, the electric field intensity E and potential V at position vector r are given, re- spectively, by y Qk(r-rk) (15.1) t=i Airs \r - rk\3 and (15.2) ift\ 4ire |r - rk\ If the charges are on the same plane (z = constant), eqs.2) become N E= 2 - xkf + (y - ykff2 (15.4) k=\ - xkf - ykf}} m To plot the electric field lines, follow these steps: 1. Choose a starting point on the field line. Calculate Ex and Ey at that point using eq.

Take a small step along the field line to a new point in the plane. As shown in Figure 15.1, a movement A€ along the field line corresponds to movements AJC and Ay along x- and y-directions, respectively. From the figure, it is evident that Ax Ex E [E2X 2 U2 ] 662 Numerical Methods Figure 15.1 A small displacement on a field line. field line new point point or Ax = (15.6) Move along the field line from the old point (x, y) to a new point x' = x + Ax, y' =y + Ay.

Go back to steps 2 and 3 and repeat the calculations. Continue to generate new points until a line is completed within a given range of coordinates. On completing the line, go back to step 1 and choose another starting point. Note that since there are an infi- nite number of field lines, any starting point is likely to be on a field line.

The points generated can be plotted by hand or by a plotter as illustrated in Figure 15. To plot the equipotential lines, follow these steps: 1. Choose a starting point. Calculate the electric field (Ex, Ey) at that point using eq.2 Generated points on £-field lines (shown thick) and equipotential lines (shown dotted).

Move a small step along the line perpendicular to £-field line at that point. Utilize the fact that if a line has slope m, a perpendicular line must have slope — Mm. Since an fi-field line and an equipotential line meeting at a given point are mutually or- thogonal there, Ax = (15.8) + Move along the equipotential line from the old point (x, y) to a new point (x + Ax, y + Ay). As a way of checking the new point, calculate the potential at the new and old points using eq.4); they must be equal because the points are on the same equipotential line.

Go back to steps 2 and 3 and repeat the calculations. Continue to generate new points until a line is completed within the given range of x and >>. After completing the line, go back to step 1 and choose another starting point. Join the points gener- ated by hand or by a plotter as illustrated in Figure 15.

By following the same reasoning, the magnetic field line due to various current distri- butions can be plotted using Biot-Savart law. Programs for determining the magnetic field line due to line current, a current loop, a Helmholtz pair, and a solenoid can be developed. Programs for drawing the electric and magnetic field lines inside a rectangular waveguide or the power radiation pattern produced by a linear array of vertical half-wave electric dipole antennas can also be written. Write a program to plot the electric field and equipotential lines due to: EXAMPLE 15.

Take QIAire = landA€ = 0. Consider the range — 5 < x , y < 5. Solution: Based on the steps given in Section 15.2, the program in Figure 15. Enough comments are inserted to make the program as self-explanatory as possible.

For example, to use the program to generate the plot in Figure 15.4(a), load program plotit in your Matlab directory. At the command prompt in Matlab, type plotit ([1 - 4 ] , [-1 0; 1 0], 1, 1, 0.01, 8, 2, 5) where the numbers have meanings provided in the program. Further explanation of the program is provided in the following paragraphs. Since the £"-field lines emanate from positive charges and terminate on negative charges, it seems reasonable to generate starting points (xs, ys) for the £-field lines on small circles centered at charge locations (xQ, yQ); that is, xs = xQ + r cos 0 (15.1b) 664 Numerical Methods function plotit(charges,location,ckEField,ckEq,DLE,DLV,NLE,NLV,PTS) figure; hold on; % Program for plotting the electric field lines % and equipotential lines due to coplanar point charges % the plot is to be within the range -5<x,y<5 % % This is the correct usage: % function plotit(charges, location,ckEField,ckEq,DLE,DLV,NLE,NLV,PTS) % % where, % charges = a vector containing the charges % location = a matrix where each row is a charge location % ckEField = Flag set to 1 plots the Efield lines % ckEq = Flag set to 1 plots the Equipotential lines % DLE or DLV = the increment along E & V lines % NLE = No.

of E-Field lines per charge % NLV = No. of Equipotential lines per charge % PTS => Plots every PTS point (i. if PTS = 5 then plot every 5th point) % note that constant Q/4*Pie*ErR is set equal to 1.0 % Determine the E-Field Lines % For convenience, the starting points (XS,YS) are radially distributed about charge locations Q=charges; XQ = location(:,1); YQ = location):,2); JJ=1; NQ = length(charges); if (ckEField) for K=1:NQ for I =1:NLE THETA = 2*pi*(1-1)/(NLE); XS=XQ(K) + 0.1*sin (THETA); XE=XS; YE=YS; JJ=JJ+1; if (~mod(JJ ,PTS)) plot (XE , YE) ; end while(1) % Find i ncrement and n e w point (X,Y) EX=0; EY=0; Figure 15.3 Computer program for Example 15.2 FIELD PLOTTING 665 for J=1:NQ R =sqrt((XE-XQ(J))~2 + (YE - YQ(J)) A 2 ) ; EX = EX +Q(J)*(XE-XQ(J) ) / (FT3) ; EY = EY +Q(J)*(YE-YQ(J))/(RA3); end E = sqrt(EX^2 + EY A 2); % CHECK FOR A SINGULAR POINT if (E <=.00005) break; end DX = DLE*EX/E; DY = DLE*EY/E; % FOR NEGATIVE CHARGE, NEGATE DX & DY SO THAT INCREMENT % IS AWAY FROM THE CHARGE if (Q(K) < 0) DX = -DX; DY = -DY; end XE = XE + DX; YE = YE + DY; % CHECK WHETHER NEW POINT IS WITHIN THE GIVEN RANGE OR TOO % CLOSE TO ANY OF THE POINT CHARGES - TO AVOID SINGULAR POINT if ((abs(XE) >= 5) | (abs(YE) >= 5)) break; end if (sum(abs(XE-XQ) < .05) >0) break; end JJ=JJ+1; if (~mod(JJ,PTS)) plot (XE,YE); end end % while loop end % I =1:NLE end % K = 1:NQ end % if % NEXT, DETERMINE THE EQUIPOTENTIAL LINES % FOR CONVENIENCE, THE STARTING POINTS (XS,YS) ARE % CHOSEN LIKE THOSE FOR THE E-FIELD LINES if (ckEq) JJ=1; DELTA = .2; ANGLE = 45*pi/180; Figure 15.3 (Continued) 666 Numerical Methods for K =1:NQ FACTOR = .5; for KK = 1:NLV XS = XQ(K) + FACTOR*cos(ANGLE); YS = YQ(K) + FACTOR*sin(ANGLE); if ( abs(XS) >= 5 | abs(YS) >=5) break; end DIR = 1; XV = XS; YV = YS; JJ=JJ+1; if (~mod(JJ,PTS)) plot(XV,YV); end % FIND INCREMENT AND NEW POINT (XV,YV) N=l; while (1) EX = 0; EY = 0; for J = 1:NQ R = sqrt((XV-XQ(J))"2 + (YV-YQ(J))^2); EX = EX + Q(J)*(XV-XQ(J))/(RA3); EY = EY + Q(J)*(YV-YQ(J))/(R"3); end E=sqrt(EXA2 + EY~2); if (E <= .00005) FACTOR = 2 *FACTOR; break; end DX = -DLV*EY/E; DY = DLV*EV/E; XV = XV + DIR*DX; YV = YV + DIR*DY; % CHECK IF THE EQUIPOTENTIAL LINE LOOPS BACK TO (X,YS) R0 = sqrt((XV - XS)"2 + (YV - YS)^2); if (R0 < DELTA & N < 50) FACTOR = 2 *FACTOR; break; end % CHECK WHETHER NEW POINT IS WITHIN THE GIVEN RANGE % IF FOUND OUT OF RANGE, GO BACK TO THE STARTING POINT % (S,YS)BUT INCREMENT IN THE OPPOSITE DIRECTION if (abs(XV) > 5 | abs(YV) > 5) DIR = DIR -2; XV = XS; YV = YS; Figure 15.2 FIELD PLOTTING m 667 if (abs(DlR) > 1) FACTOR = 2 *FACTOR; break; end else if (sumfabs(XV-XQ) < .005) >0) break; end end JJ=JJ+1; if (~mod(JJ,PTS)) N=N+1; plot(XV,YV); end end % WHILE loop end % KK end % K end % if Figure 15.3 (Continued) where r is the radius of the small circle (e.05), and 6 is a prescribed angle chosen for each £-field line. The starting points for the equipotential lines can be generated in different ways: along the x- and y-axes, along line y = x, and so on.

However, to make the program as general as possible, the starting points should depend on the charge loca- tions like those for the £-field lines. They could be chosen using eq.1) but with fixed 0(e.