Chemical Reactor Design Optimization

Click For Summary
The discussion revolves around challenges in determining the amount of styrene reacting in a chemical reactor design problem, particularly due to uncertainties regarding the molecular weight of polystyrene and the molar ratios of reactants. The group has 1000 kg of polystyrene exiting the reactor but lacks sufficient data to calculate the moles of styrene entering. MATLAB codes have been shared for problem-solving, but issues with file formats hinder some participants' ability to contribute effectively. The conversation highlights the importance of understanding the density and concentration of the monomer in the reactor feed, as well as the need to establish kinetic equations for the reactor's operation. Ultimately, the goal is to identify the most cost-effective conversion rate for the process.
member 392791

Homework Statement


Homework Equations


The Attempt at a Solution



For this problem, my group is having difficulty figuring out how much Styrene is reacting, since we don't know the molecular weight of polystyrene, we can't find the moles of polystyrene, which would be used to aid in finding the moles of styrene entering the reactor. We know there is 1000 kg of polystyrene coming out of the PFR, but also without knowing any molar ratios of the initiator or styrene, cannot find anything on a molar basis.

I included the MATLAB code for problem 46 that we wrote, and the 2nd MATLAB code was given to us by the instructor for our own use.
 

Attachments

Physics news on Phys.org
A web search would turn up the chemical formula for the styrene "unit cell" of polystyrene. Add up the atomic masses!
 
The problem is that the molecular weight of polystyrene depends on the concentration of the monomer
 
Woopydalan said:
The problem is that the molecular weight of polystyrene depends on the concentration of the monomer

My understanding was that polystyrene consists of linked styrene "cells", so that the monomer in this case is styrene which has a fixed molecular weight. 1000 kg of pure polystyrene should consist of a fixed number of moles of styrene.

So, I guess I don't understand your concern about "concentration of the monomer". Then again, I will admit that it's been a long while since I did any organic chemistry :blushing:
 
It's in the problem statement of problem 46
 
Woopydalan said:
It's in the problem statement of problem 46

Ah. Unfortunately I'm not able to read the docx format files in a way that gives me any more to work with. Much of the text is resolved, but the equations are not being rendered intelligibly. I blame it on the antiquated software I'm running. Sorry about that. Perhaps another Homework Helper will be able to provide some insight.
 
Woopydalan said:
For this problem, my group is having difficulty figuring out how much Styrene is reacting, since we don't know the molecular weight of polystyrene, we can't find the moles of polystyrene, which would be used to aid in finding the moles of styrene entering the reactor. We know there is 1000 kg of polystyrene coming out of the PFR, but also without knowing any molar ratios of the initiator or styrene, cannot find anything on a molar basis.

I included the MATLAB code for problem 46 that we wrote, and the 2nd MATLAB code was given to us by the instructor for our own use.
According to the problem statement, the density of the mixture is identical to the density of the polystyrene, and also to the density of the polystyrene. You know the amount of polystyrene being produced, and you know the conversion, so you know the mass flow rate of mixture exiting the reactor. If you divide by the density, you get the volumetric rate of flow, which is constant through the reactor. You know the density of monomer, so you know the mass concentration of monomer in the inlet stream. Since you know the molecular weight of monomer, you also know the molar concentration of monomer in the feed.

You have a plug flow reactor. If t represents the cumulative residence time through the reactor, then the kinetic balance equations for a plug flow reactor are, in terms of t, the same as a batch reactor. Set up the balance equations for the concentrations of the species in the reactor as a function of the cumulative residence time.
 
gneill said:
Ah. Unfortunately I'm not able to read the docx format files in a way that gives me any more to work with. Much of the text is resolved, but the equations are not being rendered intelligibly. I blame it on the antiquated software I'm running. Sorry about that. Perhaps another Homework Helper will be able to provide some insight.

I've taken the liberty of printing Woody's files as acrobats so that everyone who is interested can read them:
 

Attachments

The molecular weight expression for the polystyrene has been modified as of today, apparently something was wrong with the previous expression, and the initial concentration of the initiator has been changed. When I get home I will post the modified document, this design project is a living document.

Ultimately the goal is to find the conversion that will be the most cost effective
 
  • #10
The equation for the molecular weight of styrene has been modified, the initial concentration of the initiator has been modified, and the 95% purity constraint has been removed, and the recycle ratio is now 1.
 

Attachments

  • #11
You know the density of the reactor feed. The feed is essentially pure styrene monomer. From this, you should be able to calculate the initial concentration of monomer in moles/liter. What is that initial concentration C0? You have the kinetic equations for the rates of disappearance of initiator and monomer. So, you should be able to write equations for the time derivatives of the initiator concentration and monomer concentration, where the time parameter in these equations is the cumulative residence time through the reactor. What are those two ordinary differential equations?
 
  • #12
As of today, the feed cost is now $15 kg/hr for the feed, and here is our MATLAB code. The problem is that the cost is not making a U shape, so it appears that the cost just goes down as you increase the conversion, which doesn't make sense. There should be some intermediate value of conversion where it is a minimum. Just choose an nrange=[.2,.5];

If single pass conversion is 50%, then the overall conversion is 100%, which should not be the least expensive conversion

Should come as Fig. 16, haven't gotten it to be more efficient and not have multiple plots that are the same generated
 

Attachments

Last edited by a moderator:

Similar threads

  • · Replies 9 ·
Replies
9
Views
7K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
4
Views
3K
  • · Replies 28 ·
Replies
28
Views
8K
  • Poll Poll
  • · Replies 12 ·
Replies
12
Views
2K
  • · Replies 10 ·
Replies
10
Views
40K
Replies
2
Views
1K
  • · Replies 13 ·
Replies
13
Views
4K
Replies
7
Views
5K
  • · Replies 4 ·
Replies
4
Views
3K