� � � � Problems for Quasi-Linear PDEs 18.303 Linear Partial Differential Equations Matthew J. Hancock 1 Problem 1 : Traffic flow Solve the traffic flow problem ∂u ∂u + (1 − 2u) = 0, u (x, 0) = f (x)∂t ∂x for an initial traffic group 14, |x| ≥ 2 f (x) = 1 −38|x| , |x| < 2 (a) At what time ts and position xs does a shock first form? (b) Identify the important values of s and find the corresponding characteristics. Sketch the characteristics and indicate the region in the xt-plane in which the solution is well-defined (i.e. does not break down). (c) Construct tables as in class for the x and u values at the important values of s for times t = 0, 2/3, 4/3. Use these tables to sketch the density profile u = u (x, t) vs. x for t = 0, 2/3, 4/3. 2 Problem 2 : Water waves The surface displacement for shallow water waves is governed by (in scaled coordi-nates), 3 ∂h ∂h 1 + h + = 0 2 ∂x ∂t Here, h = 0 is the mean free surface of the water. Consider the initial water wave profile a sin x, 0 ≤ x ≤ π h (x, 0) = f (x) = (1) 0, x < 0, x > π 1 Fall 2006where a > 0 is constant. (a) Find the parametric solution and characteristic curves. (b) Show that two characteristics starting at s = s1 and s = s2 where s1, s2 ∈ (π/2, π) intersect at time 2 s1 − s2 tint = − 3a sin s1 − sin s2 Show that 2 , for all s1, s2 ∈ (π/2, π)tint ≥ 3a and 2 , as s1, s2 → πtint → 3a Thus the solution breaks down along the characteristics starting at s = π , when t = tc = 2/ (3a). At what position x does the solution break down? (c) Calculate ∂h/∂x using implicitly differentiation (the solution cannot be found explicitly) and hence show that along the characteristic starting at s = π, ∂h lim = −∞ − ∂x t→tc Thus the wave slope becomes vertical. (d) Plot the wave profile h (x, tc) with a = 1. Since the initial waveform is not piecewise linear, we must choose more than a few points in order to plot the solution. Choose the following 11 values of s equally spaced from 0 to π: n sn = π, n = 0, 1, 2, ..., 10 (2) 10 and construct a table at t = tc, as we did in the notes, for s, h and x. Then plot h vs. x and draw a smooth curve through the points to obtain the wave profile h (x, tc). Obviously if you use Matlab you can use 100 points or more to make the plot very smooth. Be sure to label where the wave is vertical and where the maximum displacement occurs. Plot the initial profile h (x, 0) on the same plot for comparison. 3 Problem 3 Consider the quasi-linear PDE and initial condition ut + 2u ux + 3u = 0, t > 0, −∞ < x < ∞ u (x, 0) = b sin x, −∞ < x < ∞ where b > 0 is constant. 2� � (a) Find the parametric solution and characteristic curves. (b) Show that the solution u can be written in the following implicit form 2 � � 1 3tsin x − u e 3t − 1 = ue . (3) 3 b (c) For b = 2, show that the solution first breaks down at t = tc = (1/3) ln 4 = (2/3) ln 2. Along the characteristic through (x, t) = (π, 0), find an expression for ux and show that the left limit of ux as t → t− c is lim ux = −∞. − t→tc Hint: differentiate (3) implicitly with respect to x and set (x, t) = (π, 0). (d) For b = 2, sketch the characteristics and the solution profile u (x, tc) at time tc. Note: since the initial condition is 2π-periodic, so too will be the solution, so you only have to plot one period of s values. Use the interval s ∈ [0, 2π]. When sketching the solution choose the following values for s, n sn = π, n = 0, 1, 2, ..., 10 5 and construct a table at t = tc, as we did in the notes, for s, u and x. Then plot u vs. x and draw a smooth curve through the points to obtain the wave profile u (x, tc). Obviously if you use Matlab you can use 100 points or more to make the plot very smooth. Be sure to label where the wave is vertical and where the maximum displacement occurs. Plot the initial profile u (x, 0) on the same plot for comparison. (e) Show that the solution exists for all time if 0 < b ≤ 3/2.
View Full Document
Unlocking...