2 A. CASTRO AND V. PADR ´ or (0.8) u (1) = 0 no flux. For k = −1 taking v = uk+1 we see that the steady configurations of the thermal structure are given by solutions to the semilinear equation (0.9) v (r) + N 1 r v (r) + λ(vp(r) vq(r)) = 0, u (0) = 0, r 0. with λ λ(k + 1), p = m k+1 , q = n k+1 if k + 1 0 and λ −λ(k + 1), p = n k+1 , q = m k+1 if k + 1 0. If k = −1, setting v = ln(u), equation (0.3) is equivalent to (0.10) v (r) + N 1 r v (r) + λ(emv(r) env(r)) = 0, u (0) = 0, r 0. This case requires independent analysis and is presented in Chapter 6. We classify (see Theorems 1, 2, 3 and 4) the positive solutions to the equation (0.9) for N 2, q p and λ 0 subject to either (0.7) or (0.8). We emphasize that p or q may exceed the critical exponent 2∗ 1 = (N + 2)/(N 2). In the conservative case, N = 1, the results are similar and can be obtained directly from the phase plane analysis. This is presented in Appendix 1. We also provide information on the stability-unstability of the radial steady states (see Theorems 6 and 8). The reader is referred to [10, 11, 12, 13, 14] for additional information on the stability of the steady state solution given by the thermal equilibrium temperature. Remark 1. Our results answer questions such as: What is the radius of a spherical star that supports a radial tempera- ture distribution with a given number of hot (radiating, u 1) and cold (absorbing, u 1) zones? What values of parameters k, m, n allow stable steady states with arbi- trarily large temperature in the center? What values of k, m, n allow steady states with arbitrarily large temper- ature and absorbing zones? When is a given temperature distribution stable? Remark 2. Assuming the thermal structure to have constant pressure, in dimensionless form (0.1) becomes (0.11) ∂u ∂t = u ∇(uk∇u) + λ(um un) in B(0, 1) × R+. and its radial steady states are also solutions to (0.9). Thus our classification of radial steady states extends to the constant pressure case. To study the problem (0.9) subject to (0.7) or to (0.8) we analyze the positive solutions to the initial value problem v (t) + N−1 t v (t) + λ(vp(t) vq(t)) = 0, t 0 v(0) = d, v (0) = 0. (0.12) We will denote by v(·, λ, d) the solution to the initial value problem (0.12). This function satisfies the re-scaling property (0.13) v(λt, 1, d) = v(t, λ2, d).
Previous Page Next Page