Below you will find the discussion and relevant materials for the ES 105 homework assignment due on January 30, 2009.
Discussion

MaterialsSource Code


PlotsPlease consult the main Fortran source code for details. Basically, the amplitude was calculated as before in HW#2 using hte Elliptic equation. The Parabolic equation was solved using the CrankNicholson scheme with the following parameters: r = 0.75; q = 0.999; RMSE = 0.001; h =k = L/35; Dt = 137.76. RMSE was calculated between the steady state concentration values and the transient concentration values. The spectral radius was also calculated between each time step as it approached 1. Other parameters used in this calculation were: L(Length)=150*E+3; Total # of Nodes = 3780; Half Bandwidth (NHB) = 105; N+1=36 (# of nodes in Y direction); 3*N=105 (# of nodes in X direction); w = 1.41*E4; t = 20.0*E5/Height; Amplitude, Status quo (H1=20,H2=3,H3=20)
Amplitude, Status quo (H1=20,H2=3,H3=20)
Concentration,Status quo (H1=20,H2=3,H3=20) Concentration,Status quo (H1=20,H2=3,H3=20)
Sedimentation Rate: Hot Spot and Cold Spot (H1=20,H2=3,H3=20) (HOT SPOT= +1.0727*E5; COLD SPOT= 1.1243*E4)
Transient to Steady State (H1=20,H2=3,H3=20) Transient to Steady State (H1=20,H2=3,H3=20)
Comparison: Transient to Steady State VERSUS Starting @ Steady State Overlay of Transient Results vs. Steady State Results (H1=20,H2=3,H3=20)
Hot Spot Halflife~0.2 (used the following data: hot & cold) Cold Spot Halflife~0.3
Thanks to Amir for help and suggestions.
