Let's break down and solve **Question 1**, step by step, as requested: --- # 1. **Biochemical Oxygen Demand (BOD) Dynamics in a Stratified Reservoir** ## **System Description & Notation** - Two layers: - **Epilimnion** (upper, well-mixed): volume \( V_{epi} \), concentration \( C_{epi} \) - **Hypolimnion** (lower, isolated): volume \( V_{hypo} \), concentration \( C_{hypo} \) - Inflow: \( Q_{in}, C_{in} \) - Outflow: \( Q_{out}, C_{epi} \) - **Interfacial diffusion** rate: \( E_{diff} \) - **Decay constants**: \( k_{d1} \) (epilimnion), \( k_{d2} \) (hypolimnion) - **Settling/decay**: \( k_s \) (hypolimnion) --- ## **Step 1: Epilimnion Mass Balance ODE** ### **a. Write the Mass Balance:** The general equation for the epilimnion is: \[ \frac{d}{dt}(V_{epi} C_{epi}) = \text{In} - \text{Out} - \text{Decay} - \text{Diffusive Loss to Hypolimnion} \] ### **b. Term-by-Term Expansion:** - **Inflow:** \( Q_{in} C_{in} \) - **Outflow:** \( Q_{out} C_{epi} \) - **First-order Decay:** \( k_{d1} V_{epi} C_{epi} \) - **Diffusion to Hypolimnion:** \( E_{diff}(C_{epi} - C_{hypo}) \) So, \[ \frac{d}{dt}(V_{epi} C_{epi}) = Q_{in} C_{in} - Q_{out} C_{epi} - k_{d1} V_{epi} C_{epi} - E_{diff}(C_{epi} - C_{hypo}) \] Assuming \( V_{epi} \) is constant, \[ V_{epi} \frac{dC_{epi}}{dt} = Q_{in} C_{in} - Q_{out} C_{epi} - k_{d1} V_{epi} C_{epi} - E_{diff}(C_{epi} - C_{hypo}) \] \[ \boxed{ \frac{dC_{epi}}{dt} = \frac{Q_{in}}{V_{epi}} C_{in} - \frac{Q_{out}}{V_{epi}} C_{epi} - k_{d1} C_{epi} - \frac{E_{diff}}{V_{epi}} (C_{epi} - C_{hypo}) } \] --- ## **Step 2: Hypolimnion Mass Balance ODE** ### **a. Main Terms:** - **Diffusive Gain from Epilimnion:** \( E_{diff}(C_{epi} - C_{hypo}) \) - **First-order Decay:** \( k_{d2} V_{hypo} C_{hypo} \) - **Settling Loss:** \( k_s V_{hypo} C_{hypo} \) ### **b. Mass Balance Equation:** \[ \frac{d}{dt}(V_{hypo} C_{hypo}) = E_{diff}(C_{epi} - C_{hypo}) - k_{d2} V_{hypo} C_{hypo} - k_s V_{hypo} C_{hypo} \] Assuming \( V_{hypo} \) is constant, \[ V_{hypo} \frac{dC_{hypo}}{dt} = E_{diff}(C_{epi} - C_{hypo}) - (k_{d2} + k_s) V_{hypo} C_{hypo} \] \[ \boxed{ \frac{dC_{hypo}}{dt} = \frac{E_{diff}}{V_{hypo}} (C_{epi} - C_{hypo}) - (k_{d2} + k_s) C_{hypo} } \] --- ## **Step 3: Steady-State Analytical Solution** Set temporal derivatives to zero for steady state: \[ \frac{dC_{epi}}{dt} = ,\quad \frac{dC_{hypo}}{dt} = \] ### **a. Epilimnion Equation:** \[ = \frac{Q_{in}}{V_{epi}} C_{in} - \frac{Q_{out}}{V_{epi}} C_{epi} - k_{d1} C_{epi} - \frac{E_{diff}}{V_{epi}} (C_{epi} - C_{hypo}) \] \[ = Q_{in} C_{in} - Q_{out} C_{epi} - k_{d1} V_{epi} C_{epi} - E_{diff}(C_{epi} - C_{hypo}) \] ### **b. Hypolimnion Equation:** \[ = \frac{E_{diff}}{V_{hypo}} (C_{epi} - C_{hypo}) - (k_{d2} + k_s) C_{hypo} \] \[ = E_{diff}(C_{epi} - C_{hypo}) - (k_{d2} + k_s) V_{hypo} C_{hypo} \] --- ### **c. Solve the Coupled Algebraic System** #### **From Hypolimnion Equation:** \[ E_{diff}(C_{epi} - C_{hypo}) = (k_{d2} + k_s) V_{hypo} C_{hypo} \] \[ E_{diff} C_{epi} = [E_{diff} + (k_{d2} + k_s) V_{hypo}] C_{hypo} \] \[ C_{hypo} = \frac{E_{diff}}{E_{diff} + (k_{d2} + k_s) V_{hypo}} C_{epi} \] --- #### **Substitute \( C_{hypo} \) into Epilimnion Equation:** \[ Q_{in} C_{in} - Q_{out} C_{epi} - k_{d1} V_{epi} C_{epi} - E_{diff}(C_{epi} - C_{hypo}) = \] \[ Q_{in} C_{in} - Q_{out} C_{epi} - k_{d1} V_{epi} C_{epi} - E_{diff}[C_{epi} - \frac{E_{diff}}{E_{diff} + (k_{d2} + k_s) V_{hypo}} C_{epi}] = \] \[ Q_{in} C_{in} - Q_{out} C_{epi} - k_{d1} V_{epi} C_{epi} - E_{diff} C_{epi} + \frac{E_{diff}^2}{E_{diff} + (k_{d2} + k_s) V_{hypo}} C_{epi} = \] \[ Q_{in} C_{in} = C_{epi} \left[Q_{out} + k_{d1} V_{epi} + E_{diff} - \frac{E_{diff}^2}{E_{diff} + (k_{d2} + k_s) V_{hypo}} \right] \] \[ C_{epi,ss} = \frac{Q_{in} C_{in}}{Q_{out} + k_{d1} V_{epi} + E_{diff} - \frac{E_{diff}^2}{E_{diff} + (k_{d2} + k_s) V_{hypo}}} \] --- \[ C_{hypo,ss} = \frac{E_{diff}}{E_{diff} + (k_{d2} + k_s) V_{hypo}} C_{epi,ss} \] --- ## **Summary of Final Answers** ### **Steady-State Epilimnion:** \[ \boxed{ C_{epi,ss} = \frac{Q_{in} C_{in}}{Q_{out} + k_{d1} V_{epi} + E_{diff} - \frac{E_{diff}^2}{E_{diff} + (k_{d2} + k_s) V_{hypo}}} } \] ### **Steady-State Hypolimnion:** \[ \boxed{ C_{hypo,ss} = \frac{E_{diff}}{E_{diff} + (k_{d2} + k_s) V_{hypo}} C_{epi,ss} } \] --- ## **Explanation** 1. **Write mass balances** for each layer, including all relevant inflow, outflow, decay, and transfer terms. 2. **Set time derivatives to zero** for steady-state. 3. **Express \( C_{hypo} \) in terms of \( C_{epi} \)** using the hypolimnion equation. 4. **Substitute back** into the epilimnion equation and solve algebraically for \( C_{epi,ss} \). 5. **Back-calculate \( C_{hypo,ss} \)** from the value of \( C_{epi,ss} \). If you need the solution steps for **Question 2** or have additional questions, let me know!