Biomod/2013/NanoUANL/Reactor

What is a reactor?

Introduction

The CCMV capsid could be considered as a continuous stirred-tank reactor with an accumulation of the product inside the capsid. This type of reactor is a common and ideal type in chemical engineering. The open system will be treated as it operates on a steady-state assumption, where the conditions of the reactor will not change with time. The reactor proposed is a complete opposite type from tubular plug-flow and stirred batch reactors, and can be very useful when studying the behavior of a gas, liquid or solid.

The reactor is modeled by that of a Continuous Ideally Stirred-Tank Reactor (CISTR), assuming perfect mixing in the container.

Enzymatic Reaction

The general reaction scheme is described as follows:

[math]\displaystyle{ E + S \leftrightarrow ES \rightarrow E^0 + P }[/math]

With a reaction rate of:

[math]\displaystyle{ \frac{d[ES]}{dt}=k_1[E][S]-k_{-1}[ES]-k_2[ES] }[/math] . . . 1.1

This equation is affected by the constants k₁ , k_-1 and k₂.

Considerations

In order to solve the reactor, we established the following considerations in our system:

• Uniform distribution throughout the reactor
• k₁ and k₂ >> k_-1
• One enzyme per reactor/VLP
• Tortuosity approaches zero during diffusion
• Excess of specie S

Mass balance

Material balances are important, as a first step in devising a new process (or analyzing an existing one). They are almost always a prerequisite for all calculations for process engineering problems. The concept of mass balance is based on the physical principle that matter cannot be either created nor destroyed, only transformed. The law of mass transformation balances describe the mass of the inputs of the process with the output, as waste, products or emissions. This whole process is accounting for the material used in a reaction.

Applying a mass balance to our system we got:

ACCUMULATION = INPUT - OUTPUT - DISAPPEARANCE BY REACTION

where

Input= F₀

Output= F₀(1-x_S)

Disappearance = V(-r_S)

Accumulation = [math]\displaystyle{ \tfrac{d[P]}{dt} }[/math]

[math]\displaystyle{ \frac{d[P]}{dt}=F_0-F_0(1-x_S)-V(-r_S) }[/math] . . . 1.2

The inlet and outlet flow were determined by diffusion. For the simplification of the diffusion phenomenon we considered:

• Constant temperature
• Constant pressure
• Species B stays in a stationary state (it does not diffuse in A)
• The container (VLP) has a spherical shape

A mass balance, applied to a spherical envelope is described as follows:

[math]\displaystyle{ \frac{d}{dr}(r^2N_{Ar})=0 }[/math] . . . 1.3

where N_Ar represents molar flux. When N_Br=0 we obtain

[math]\displaystyle{ \frac{d}{dr}(r^2\frac{cD_{AB}}{1-x_A}\frac{dx_A}{dr})=0 }[/math] . . . 1.4

At a constant temperature the product (cD_AB) is equally constant and x_A=1-x_B, the equation can be integrated into the following expression¹:

[math]\displaystyle{ F_A=4\pi r_1^2N_{Ar}|_{r=r1}= \frac{4\pi cD_{AB}}{1/r_1-1/r_2} \ln\frac{x_{B2}}{x_{B1}} }[/math] . . . 1.5

where x are the fractions, c is the concentration and r are the respective radius.

This equation defines the nanoreactor inflow.

For the outflow we made a similar analysis, but in this case we are describing uniquely the specie P flow. We neglect the possibility of an outflow of specie S because:

• The gradient of concentration of S tends to stay inside the capsid
• The positive charge in the outside of the VLP made a repulsion of the specie S

Diffusion coefficient

For S and P being ionic silver and reduced silver respectively. The ionic silver diffusion coefficient in solution is described by Nerst's equation (1888)²:

[math]\displaystyle{ D_{AB}°= \frac{RT}{F^2} \frac{\lambda^o_+\lambda^o_-}{\lambda^o_++\lambda^o_-} \frac{|Z_-|+|Z_+|}{|Z_+Z_-|} }[/math] . . . 1.6

where

• F = Faraday's constant [A·s/g_eq]
• D_AB° = Diffusion coefficient at infinite dilution [m²/s]
• λ₊° = Cationic conductivity at infinite dilution
• λ_-° = Anionic conductivity at infinite dilution
• Z⁺ = Cation valence
• Z^- = Anionic valence
• T = Absolute temperature [K]

Boiling temperature

Via Joback's method, we estimate the normal boiling temperature:

[math]\displaystyle{ T_b=\mathbf{198} + \sum_{k} N_k(tbk) }[/math] . . . 1.7

in which N_k is the number of times that the contribution group is present in the compound.

Critical temperature

Using a similar approach, also by Joback, we estimate the critical temperature:

[math]\displaystyle{ T_c=T_b\Bigg[ \mathbf{0.584}+\mathbf{0.965} \bigg\{\sum_{k} N_k(tck) \bigg\} - \bigg\{\sum_{k} N_k(tck) \bigg\}^2 \Bigg]^{-1} }[/math] . . . 1.8

Joback's Method Contributions (Table C1. Prausnitz)

Contribution Group	tbk, K	tck, K	Contribution Group	tbk, K	tck, K	Contribution Group	tbk, K	tck, K	Contribution Group	tbk, K	tck, K
CH3	23.58	0.0141	CH2(ss)	27.15	0.01	ACOH	76.34	0.0240	NH	50.17	0.0295
CH2	22.88	0.0181	CH(ss)	21.78	0.0122	O	22.42	0.0168	NH(ss)	52.82	0.0130
CH	21.74	0.0164	C(ss)	21.32	0.0042	O(ss)	31.22	0.0098	N	11.74	0.0169
C	18.25	0.0067	=CH(ds)	26.73	0.0082	C=O	76.75	0.0380	=N-(ds)	74.60	0.0255
=CH2	18.18	0.0113	=C(ds)	31.01	0.0143	C=O(ds)	94.97	0.0284	=NH	X	X
=CH	24.96	0.0129	F	-0.03	0.0111	CH=O	72.20	0.0379	CN	125.66	0.0496
=C	24.14	0.0117	Cl	38.13	0.0105	COOH	169.09	0.0791	NO2	152.54	0.0437
=C=	26.15	0.0026	Br	66.86	0.0133	COO	81.10	0.0481	SH	63.56	0.0031
≡CH	9.20	0.0027	I	93.84	0.0068	=O	-10	0.0143	S	68.78	0.0119
≡C	27.38	0.0020	OH	92.88	0.0741	NH2	73.23	0.0243	S(ss)	52.10	0.0019

_{(ss) indicates a group in a nonaromatic ring, (ds) indicates a group in an aromatic ring, X indicates a non-available parameter}

Conductivity

The conductivity of the compound is determined by the Sastri method:

[math]\displaystyle{ \bold{\lambda_L=\lambda_ba^m} }[/math] . . . 1.9

where

• λ_L = thermic conductivity of the liquid [ W/(m·K)]
• λ_b = thermic conductivity at normal boiling point [ W/(m·K)]
• T_br = T_b/T_c = normal boiling reduced temperature
• T_r = T/T_c = reduced temperature
• T_c = critical temperature [K]

and

[math]\displaystyle{ m=1-\bigg(\frac{1-T_r}{1-T_{br}}\bigg)^n }[/math] . . . 1.10

with a = 0.16 and n = 0.2 for the compound

Sastri's Contributions (Table 10.5. Prausnitz)

Hydrocarbon Groups	Δλb	Non-Hydrocarbon Groups	Δλb	Non-Hydrocarbon Groups	Δλb
CH3	0.0545	O	0.0100	N(ring)	0.0135
CH2	0.0008	OH2	0.0830	CN	0.0645
CH	-0.0600	OH3	0.0680	NO2	0.0700
C	-0.1230	CO(ketone)	0.0175	S	0.0100
=CH2	0.0545	CHO(aldehyde)	0.0730	F4	0.0568
=CH	0.0020	COO(ester)	0.0070	F5	0.0510
=C	-0.0630	COOH(acid)	0.0650	Cl	0.0550
=C=	0.1200	NH2	0.0880	Br	0.0415
Ring1	0.1130	NH	0.0065	I	0.0245
		NH(ring)	0.0450	H6	0.0675
		N	-0.0605	3 member ring	0.1500
				Ring7(other)	0.1100

¹_{In polycyclic compounds, all rings are treated as separated rings, ²In aliphatic primary alcohols and phenols with no branch chains, ³In all alcohols except as described in², ⁴In perfluoro carbons, ⁵In all cases except as described in ⁴, ⁶This contribution is used for methane, formic acid, and formates, ⁷In polycyclic non-hydrocarbon compounds, all rings are considered as non-hydrocarbon rings}

Accumulation

Every time we propose a Matter Balance is quietly easy to assume a Steady-State, but in real life we could spec some disturbances in the out-flux stream caused by the Accumulation. The Accumulation of substance inside the reactor is highly common and when the capsid of the CCMV simulates a reactor it is not an exception. In this case an agglomeration of Silver Nano-particles of different sizes will be notorious and it is described by solving the differential equation of concentration of product in function of time present in the Mass Balance previously described.

A common method to approach the correct value is using the exponential matrix. Knowing that for an ordinary differential equation like

[math]\displaystyle{ Q(t)=\frac{dy}{dt}+P(x)y }[/math]

and

[math]\displaystyle{ \bold{y(0)=y_0} }[/math]

for a system of equations the solution is given by

[math]\displaystyle{ x=e^{At}x_0+\int_0^t e^{A(t-t)}Bu(t)\,\mathrm dt }[/math]

for this case A is a matrix so we can not solve the exponential directly so we solve it by Laplace

[math]\displaystyle{ \bold{xD} }[/math]

---

Silver density

[PENDIENTE]

---

References

1. Bird R. B, Stewart W.E. & Lightfoot E.N., Fenómenos de Trasporte, Reverté Ediciones SA de CV 2006, pages 17-9 – 17-11.
2. Anthony L. Hines & Robert N. Maddox, Mass Transfer Fundamentals and Applications, Prentince Hall PTR 1985, pages 34-35
3. Poling, Prausnitz & O’Connel, The Properties of Gases and Liquids, 5th Edition. McGrawl-Hill 2001, pages 2.26,10.45-10.46, C.2-C.4
4. Fogler H. Scott, Elements of Chemical Reaction Engineering 4th Edition, Pretince Hall 2006, pages 37-45

[AGREGAR CITAS DE FOGLER Y DE BEQUETTE O COUGHANOUR O LUYBEN Y TAL VEZ DE GEANKOPLIS]