辅导BMET9960、辅导MATLAB程序语言
AMME2000/BMET2960/BMET9960 - Assignment 2, 2023
Due: 11:59 pm Friday 19th May (Week 12) 2023
Assignment Informa on
Assignment 2 focuses on your understanding of the analy cal and numerical solu on to the Wave Equa on and
Laplace Equa on, and also the analy cal solu on to the Heat Equa on without boundary condi ons.
• Your report must be your own work. Please see the FAQ on Academic Integrity for informa on on
what is appropriate and what is not appropriate for Assignments.
Any students who use third-party services such as Chegg to post or view assignment ques ons will
receive 0 for the assignment and be reported. These services provide us with detailed data on user
viewing/pos ng because they are hos ng our copyrighted material. For more informa on, see the
first lecture where academic integrity, including penal es, was discussed at length.
• Present your assignment as a concise report in PDF format;
• The report must be typed, including all mathema cal working; handwriten assignments will
receive a mark of 0 and marks will be deducted for screenshots of equa ons and/or figures;
• Your MATLAB code must be included as an appendix to the assignment report. 10% of the
assignment marks are allocated for the readability and quality of your MATLAB code;
• The report should not exceed 10 pages (this limit does NOT include the code in the appendix);
• Structure the report using the appropriate Sec on and Ques on numbers (e.g. 1.1, 2.2). Ensure all
relevant working and final answers are in that sec on;
• All figures and tables in your report must be numbered and must be referred to and discussed;
• Submit your report to Turni n by the due date. Late submissions incur a penalty of 5%/day.
We have prepared an exemplary assignment to guide you in wri ng your own.
QUESTION 1: WAVE EQUATION [40 marks]
Relevant parts of the course for this ques on:
• Week 6 Lecture 1
• Week 6 Lecture 2
• Week 7 Pre-work video: ‘The starting problem’
• Week 7 tutorial and Week 8 tutorial
A torsional wave can propagate along a stiff rod if the rod experiences a torque. The resulting wave can be
modelled using the wave equation:
??2??
????2 = ??2 ??2??
????2 (1.1)
where ?? = ??(??,??) is the torsion angle of the wave in radians at location ?? (along the length of the rod) and time ??,
and ?? is given by:
?? = ??
?? (1.2)
where ?? is the shear modulus of the rod in units of Pa, and ?? is the density of the rod in units of kg.m-3. If the
boundary conditions are ??(0,??) = ??(??,??) = 0, the general solution to Eqn (1.1) is given by:
??(??,??) = [???? cos(??????) + ???? sin(??????)]
∞
??=0
sin
??????
?? (1.3)
where ???? = ??????
?? are the eigenvalues, ?? is the length of the rod, and ???? and ???? are coefficients which depend on the
initial conditions. Consider a thin metallic filament of length L being used as a sensor in an engineering application
(Figure 1.1):
The filament is prone to torsional waves, which in turn impact the sensor function. To fully understand the impact,
you want to characterise the torsion ?? = ??(??,??) in the sensor filament. You may assume the boundary conditions
are ??(0,??) = ??(??,??) = 0 so that the solution is given by Eqn (1.3).
1.1. If the initial torsion angle along the length of the rod is 0 everywhere and the initial torsional velocity along
the rod is given by ??(??), shown in Figure 1.2, write down the displacement initial condition and the velocity
initial condition mathematically. [4 marks]
1.2. Starting from Eqn (1.3), show why for this particular problem the solution reduces to:
??(??,??) = ???? sin(??????)
∞
??=1
sin
??????
?? (1.4)
[2 marks]
1.3. Given the initial velocity profile (Figure 1.2), solve for the ???? coefficients in Eqn (1.4). [4 marks]
Note: You may use the fact that the Fourier sine series for a triangular function of the form:
is given by:
??(??) = 2????2??????
????
??
(?? − 1)??2??2
∞
??=1
??????
??????
?? (1.5)
Therefore, you do NOT need to perform any integrations yourself for this question.
1.4. Implement your analytical solution for the torsional wave in MATLAB. For this implementation use the
following parameter values:
• ?? = 0.1 m
• ?? = 30 GPa
• ?? = 9 × 103 kg. m−3
• ??0 = 50 rad. s−1
• ???? = 101 (number of points in spatial domain)
In your report, provide a plot of your solution at 4 different time points, all overlaid on the same plot. Make
sure you indicate the time corresponding to each curve and provide a brief discussion of your results. Include
your MATLAB code in the appendix of your report. [10 marks]
1.5. Write down the form of the CTCS scheme to solve this problem numerically. (Note: You do not need to
derive the CTCS scheme.) Briefly justify why the CTCS scheme is a suitable choice to solve this problem.
(Hint: Consider the properties of the scheme and the expected physics of the problem. You do not need to
include any derivations for this part.) [2 marks]
1.6. Briefly describe what is meant by the ‘starting problem’ in the CTCS scheme and explain how you will
handle it for this problem. You should outline the specific steps and how you plan to implement those steps
in MATLAB. [4 marks]
1.7. Now implement the CTCS numerical scheme in MATLAB to solve the problem numerically. For the
numerical implementation you should use the following parameters:
• ∆?? = the maximum stable time step
• ???? = 101 (number of points in spatial domain)
In your report you should include the following:
• Explain how you computed the maximum stable time step. (Note that you are NOT required to
perform a von Neumann stability analysis.) [2 marks]
• Provide a plot showing the analytical and numerical solutions overlaid at the same 4 time points used
in 1.4. Briefly discuss your results. [6 marks]
Note: Make sure your implementation compares the analytical and numerical solutions at the same (not different)
time at each step.
1.8. Plot the maximum absolute percentage error between the analytical and numerical solution as a function of
time over one complete cycle of the wave up and down the rod. On the same plot overlay the mean absolute
percentage error as a function of time. Use 51 evenly spaced time points for your horizontal axis. (Note: In
your computation of the mean and maximum absolute percentage error, make sure you avoid any time points
or spatial locations that will introduce a divide-by-0 error.) Briefly explain how you generated the plot and
discuss the results. [6 marks]
QUESTION 2: LAPLACE EQUATION [40 marks]
Relevant parts of the course for this ques on:
• Week 7 Lecture 1
• Week 7 Lecture 2
• Week 9 tutorial
An engineering control experiment involves exposing a small rectangular wafer to several lasers. The laser light
results in a ‘photon pressure’ distribution on the wafer which we can model using the Laplace equation:
??2??
????2 +
??2??
????2 = 0 (2.1)
where ?? = ??(??, ??) is the 2-dimensional pressure distribution on the wafer in units of Pa. Figure 2.1 shows the
wafer dimensions and the pressure distribution at the boundary corresponding to ?? = ??. The other boundaries have
no pressure applied.
2.1 For the scenario shown in Figure 2.1, write down the 4 boundary conditions for the pressure S
mathematically. [4 marks]
2.2 Use separation of variables to show that the steady-state pressure distribution for the scenario in Figure 2.1 is
given by:
??(??, ??) = ????
∗ sinh
??????
??
∞
??=1
sin
??????
?? (2.2)
where ????
∗ are coefficients yet to be determined. [4 marks]
2.3 Apply the non-zero boundary condition to solve for the ????
∗ coefficients and hence write down the complete
solution for ??(??, ??). [4 marks]
Note: For this question you may use Eqn (1.5) to help you solve for ????
∗ – i.e. you should NOT perform any
integrations for this part.
2.4 Suppose the final arrangement of the lasers results in the boundary conditions shown in Figure 2.2:
Briefly explain how you would obtain the steady-state solution for ??(??, ??) for this new scenario, and write
down this solution. Note: You do NOT need to perform any lengthy derivations, you should simply explain
what mathematical principles allow you to adapt your solution from 2.3 to this new solution. [4 marks]
2.5 Implement your analytical solution from 2.4 in MATLAB. For your implementation use the following
parameters:
• ?? = ?? = 0.1 m
• ???? = ???? = 101 (number of points in each dimension of the spatial domain)
• No. of Fourier terms = 100
• ??1 = ??2 = 2 × 10−3 Pa
In your report, provide a clearly labelled contour plot of the steady-state solution. Note that you can use the
pcolor function in MATLAB to generate this plot. Be sure to provide a brief discussion of your plot. Include
your MATLAB code in the appendix of your report. [8 marks]
2.6 Now implement the Gauss-Jacobi numerical scheme from Week 7 Lecture 2 to solve the same problem
numerically. In your report:
• Briefly outline the steps you used to implement the Gauss-Jacobi method, including your stopping
criteria. [4 marks]
• Present 2 plots side-by-side – the left plot should be your numerical solution to the steady-state
pressure distribution; the right plot should be 代写AMME2000the signed difference image between the numerical and
analytical solutions. Discuss your results. [8 marks]
2.7 The scientists want to know how much they can shorten the wafer before the pressure at the dead-centre
becomes too great. Using your analytical solution, determine at what wafer width, ????ℎ, the pressure at the
centre first exceeds 75% of ??1. You should determine ????ℎ using a tolerance of 1 µm. Repeat your estimate
using your numerical solution. In your report:
• Briefly outline the method you used to compute ????ℎ. [2 marks]
• Present your results using the analytical and numerical solutions and briefly discuss any differences.
[2 marks]
QUESTION 3: Fourier Integral [10 marks]
Relevant parts of the course for this ques on:
• Week 8 Lecture 2
• Week 9 Lecture 1
• Week 10 tutorial
A breach in a factory pipeline has caused a large amount of the insecticide permethrin to be released into a
long narrow waterway (Figure 3.1). The concentration C(x, t) of permethrin at location x and time t along
the waterway can be modelled using Fick’s second law, as the heat equation:
????
???? = ??
??2??
????2 (3.1)
where D is the mass diffusivity. Because the waterway is very long, there are no boundary conditions.
3.1 State the boundary conditions mathematically and briefly explain why you cannot solve this
problem analytically using Fourier series. [2 marks]
3.2 C(x, t) can be obtained using the Fourier Integral form of the heat equation solution, as shown in Eqn
(3.2):
??(??,??) =
1
?? ??(??) cos(???? − ????)??−????2??
∞
−∞
???? ???? (3.2)
∞
0
where p is a continuous variable and f(x) = f(v) is the initial condition. Using the known integral on
the Data Sheet:
??−??2
cos(2????) .
∞
0
???? = √??
2 ??−??2
(3.3)
show how Eqn (3.2) can be simplified to:
??(??,??) = 1
2√?????? ??(??)??
−
(??−??)2
4????
????
∞
−∞
(3.4)
[2 marks]
3.3 Assume the initial condition is given by:
??(??, 0) = ??(??) = ??0 if − 10 < ?? < 10
0 otherwise
as shown in Figure 3.2. By making a suitable variable substitution, solve the remaining integral in
Eqn (3.4) to show that the final solution for C(x,t) can be expressed as the following sum of error
functions:
??(??,??) = ??0
2 erf
(10 + ??)
2√???? + erf
(10 − ??)
2√???? (3.5)
[4 marks]
3.4 There is an endangered frog species 300 m from the initial site of deposition of the permethrin. The
frogs will not survive if they encounter a concentration of permethrin greater than 1% of C0. Assuming
that D=0.38 m2
/s and using the table of error function values in the Data Sheet to help you, determine
if the frogs will be able to survive for 24 hours before the habitat relocation team can rescue them.
Show all working and explain all reasoning. [2 marks]
WX：codehelp