40 Separable Differential Equations
Learning Objectives
By the end of this section, you should be able to:
Provide initial conditions for well-mixed separable transient single-unit processes.
Solve separable transient balances to find a property of a system at any given time.
Separable Differential Equations
Separable differential equations are differential equations where the variables can be isolated to one side of the equation. Take the following differential equations:
1 – dxdy=(x3+x)∗(y−y2)
This equation is separable because you can completely isolate the x and y variables as follows:
dxx3+x=dy∗(y−y2)
2 – dxdy=xx+y
This equation is non-separable because you cannot completely isolate the x and y variables:
dy∗(x+y)=dx∗x
These non-separable equations we will discuss later.
Example: Chemical Reactor
Consider a “continuous stirred-tank reactor” (CSTR). CSTRs are reactors with continuous feed and exit streams and some kind of mixer.
Say we know the following information about this CSTR:
- The feed enters at a constant volumetric flowrate of ˙V0 in L/s
- The volume of the tank is V in L.
- Initially (t=0), the tank is filled to Vi in L
- The exit stream flows at a constant rate of ˙V in L/s
- We can assume that the density of all streams in the system is constant at ρ in g/L
We want to write a balance for the total (overall) mass in the system under transient conditions.
We start off by writing out the overall balance:
IN−OUT+GEN−CON=ACC
Mass is not being consumed or generated, just changed from one substance to another. This means the GEN and CON terms are negligible. We get:
IN=˙V0∗ρ
OUT=˙V∗ρ
ACC=dMdt=d(V∗ρ)dt=ρ∗dVdt
The units for the IN, OUT, and ACC terms are kg/s. Simplifying the balance, we get:
˙V0∗ρ−˙V∗ρ=ρ∗dVdt
Since the densities are all constant, we can cancel them out:
˙V0−˙V=dVdt
Using separation of variables (from calculus), we can integrate both sides. We want to find a given value about our system at a specific final time (tf), starting from an initial time (say t0=0). We know the initial volume Vi and want to find the final volume Vf.
(˙V0−˙V)dt=dV
∫tft0(˙V0−˙V)dt=∫VfVidV
(˙V0−˙V)∗(tf−t0)=Vf−Vi
Vf=Vi+(˙V0−˙V)∗(tf−t0)
Let’s try substituting in some numbers to this equation. Say the rate of flow in is 5 L/min, and the flow out is 6 L/min with 300 L initially in the tank. How much water remains in the tank after 1 hour? Using the formula we derived, we can find this:
Vf=Vi+(˙V0−˙V)∗(tf−t0)
[latex]V_{f} = 300 L+(5 \frac{L}{min} - 6 \frac{L}{min})*(60 min-0 min)$$
Vf=300L−60L
Vf=240L
Note, you may say this is obvious and there is no need to derive our equation, and in this case our equation was relatively simple. However, we want to get in the practice of formulating these equations for when things start getting more complicated. In this class, we'll focus on formulating these equations rather than finding the solution.
Exercise: Transient Mass Balance

Solution
Step 1: Determine what terms in the general balance are zero or negligible.
IN−OUT+GEN−CON=ACC
Since there are no feed streams and no reactions (no mass generated or consumed), the IN, GEN, and CON terms can be omitted from the balance.
−OUT=ACC
Step 2: Write out the mass balance for species A in the system.
(−˙mout,A)=dMAdt
Step 3: Expand the mass terms.
Recall that Mass=Volume∗density, so we can use that to express both terms in the mass balance:
(−˙Vout,A∗ρ)=ρ∗dVAdt
Since the density is constant, the density terms on both sides cancel each other out
(−˙Vout,A)=dVAdt
The volume terms can be further expanded. Recall that the volume of a cylinder Area∗Height or A∗h. In differential terms, dVdt=Adhdt+hdAdt. Since the area is constant in this case, this simplifies to: dVdt=Adhdt.
(−˙Vout,A)=Atank∗dhdt
Atank=π4∗Dtank2=π4∗2m2=3.14m2
Step 4: Solve the integral from t = 0 to t = 10 mins
(−˙Vout,A)dt=Atank∗dh
∫t=10min0min(−˙Vout,A)dt=∫hfinalhintialAtank∗dh
∫t=10min0min(−0.10m3min)dt=∫hfinalhintial3.14m2∗dh
(−0.1m3min)∗t|t=10min0min=(3.14m2)∗h|hfinal5m
(−0.1m3min)∗(10min−0min)=(3.14m2)∗(hfinal−5m)
−1m3∗13.14m2=hfinal−5m
−0.32m=hfinal−5m
(−.32+5)m=hfinal
hfinal=4.68m
Exercise: Transient Mass Balance
Suppose we have a tank with an outlet at the bottom as shown below. Water is the only thing in the tank (we'll call this species A). The water flows out of the tank at a rate of √19.6×hm/s, where h is the height of water in m at any specific time (note that expression like this for outlet flow come from conservation of energy in a tank with gravity driving flow out of the tank). The tank has a cross-sectional area of 1m2 and the outlet pipe has a cross-sectional area 10cm2. If the initial height of the water in the tank is 1 m, using a transient mass balance, calculate the height of water in the tank after 5 mins.
Solution
Step 1: Determine what terms in the general balance are zero or negligible.
IN−OUT+GEN−CON=ACC
Since there is no inlet stream and no chemical reactions in the tank, there is no IN, GEN or CON terms:
−OUT=ACC
Step 2: Write out the mass balance for species A in the system.
(−˙mout)=dMdt
Step 3: Express each term using the given quantities:
Replace Mass by Volume×density for both terms:
(−˙Vout×ρ)=dV×ρdt
Because ρ are multiplied by both sides and we assume it is constant throughout the system, we can cancel ρ in this step:
(−˙Vout)=dVdt
The volume of fluid flowing out of a pipe at any instant can be calculated using ˙V=u×A, where u is the instantaneous velocity of the fluid and A is the cross-sectional area of the pipe. We can use this to express the outlet flow. For the accumulation term, it is not beneficial to do this because we don't know the expression for u.
For the accumulation term, we can espress dV as A1×dh because we know that A1 is constant, and this leaves only one variable (h) for us to solve for in the mass balance:
−u×A2=A1dhdt
We can replace u using u=√19.6×h given in the question. T
−√19.6×h×A2=A1dhdt
Step 4: Solve the integral to find hf:
∫tft0−√19.6×h×A2dt=∫hfh0A1dh
Separate the equation by moving all terms involving h to one side:
∫tft0−A2A1dt=∫hfh01√19.6×hdh
Replace the variables by the given information: A2=10cm2=10×10−4m2, A1=1m2, t0=0s, tf=5min=300s.
∫t=3000−10×10−41dt=∫hfh0=11√19.6×h−12dh−10×10−4×t|t=3000=2√19.6×h12|hfh0=1−10×10−4×(300−0)=2√19.6×h12f−2√19.6×(1)12(−10×10−4×(300−0)+2√19.6×(1)122√19.6)2=hfhf=0.11m
Feedback/Errata