My ChemE Journey



Tutorial questions and suggested solutionsNumber of bacteria living with body of the average healthy adult human is estimated to outnumber human cells 10:1. Estimate mass of all the bacteria in your body (or body of your favourite pop-star). Use internet for additional information. Explain your quantitate reasoning. You can use internet for additional information that you might require (but not to get the exact number).Assumptions:Average adult human weighs in the order of 102 kg.Average diameter of a human cell is about 10-5 m.Average diameter of bacteria is 10-6 m. Cell density of human cell and bacteria is similar to that of water 1kg/m3Volume of one cell = 10-15 m3/cellMass of one cell = 10-15 kg/cellNumber of cells = 102 kg 10-15 kg/cell = 1013 cellSince number of bacteria outnumber human cells 10:1 1014 bacteriaVolume of one bacteria = 10-18 m3/bacteriaMass of one bacteria = 10-18 kg/bacteriaTotal mass of bacteria = 10-18 kg/bacteria 1014 bacteria = 10-4 kg = 0.1gCalculate the energy yield per liter of oxygen consumed for each food type, and note that it is roughly constant. Thus we can determine a person’s metabolic rate simply by measuring her rate of oxygen consumption. In contrast, the CO2/O2 ratios are different for the different food groups; this circumstance allows us to estimate what is actually being used as the energy source, by comparing oxygen intake to carbon dioxide output.Foodkcal/gliters O2/gliters CO2/gkcal/liters O2CO2/O2Carbohydrate4.10.810.815.061Fat9.31.961.394.740.71Protein4.00.940.754.260.80Alcohol7.11.460.974.860.66Based on the data given, the average energy yield per liter of O2 is 4.73 0.34 kcal.It is noted that for the CO2/O2 ratios, the average value is 0.79 0.15 which seem to be more constant as compared to the average energy yield per liter of O2.b. An average adult at rest uses about 16 liters of O2 per hour. The corresponding heat release is called the “basal metabolic rate” (BMR). Find it, in kcal/hour and in kcal/day.16 liters/hr 4.73kcal/liter O2 = 75.7kcal/hr75.7kcal/hr 24 hr/day = 1820kcal/dayc. What power output does this correspond to in watts?75.7kcal/hr 4184J/kcal 1hr/3600s = 88Wd. Typically, the CO2 output rate might be 13.4 liters/hour. What, if anything, can you say about the type of food materials being consumed?13.4/16 = 0.8375. Although this value lies somewhere between carbohydrates and proteins, we cannot just simply assume that it is a combination of both. In actual fact for a typical meal, we often consume a mixture of carbohydrates, fat and protein (alcohol sometimes).e. During exercise, the metabolic rate increases. Someone performing hard labor for 10 hours a day might need about 3500 kcal of food per day. Suppose the person does mechanical work at a steady rate of 50W over 10 hours. We can define the body’s efficiency as the ratio of mechanical work done to excess energy intake (beyond the BMR calculated in (b)). Find this efficiency.The BMR calculated was found to be about 2000kcal/day. Thus the excess energy intake would be about 1500kcal. Body efficiency:(50J/s 10hr 3600s/hr) (1500kcal 4184J/kcal) = 0.287 0.3Binomial distribution is defined as:Pk=N!k!N-k!pk1-pN-kProve that it is normalised, i.e.k=0NP(k)=1Using the binomial theorem:(a+b)N=k=0NNka(N-k)bk=k=0NN!k!N-k!a(N-k)bkThus:k=0NP(k)=k=0NN!k!N-k!pk1-pN-k=k=0NNkpk1-pN-kLetting q=1-p,k=0NNkpkqN-k=(p+q)N=1The binomial distribution is normalized. The continuous Gaussian distribution is given by:ρx=1σ2πe-x-μ22σ2where x∈[-∞,∞]. Gaussian integral formula is given by:Ib=-∞∞e-by2dx=πbPlease show that:ddbIb=-∞∞y2e-by2dxUsing Leibniz integral rule: ddxabfx,tdt=ab??xfx,tdtddbIb=-∞∞??be-by2dy=-∞∞-y2e-by2dy=--∞∞y2e-by2dxProve that:-∞∞ρxdx=1 and x=μLet:z=x-μσ2 and dzdx=1σ2Then:-∞∞ρxdx=-∞∞1σ2πe-x-μ22σ2dx=-∞∞1πe-z2dz=1π-∞∞e-z2dzUsing:-∞∞e-x2dx2=-∞∞e-x2dx?-∞∞e-y2dy=-∞∞-∞∞e-x2+y2dxdyConverting into polar coordinates where x2+y2=r2 and dxdy=rdrdθ:-∞∞-∞∞e-x2+y2dxdy=02π0∞e-r2rdrdθ=2π?0∞r?e-r2drLet:u=-r2 and dudr=-2r2π?0∞r?e-r2dr=2π?-∞012eudu=π?e0-e-∞=πThus:-∞∞e-z2dz2=π and-∞∞e-z2dx=π∴-∞∞ρxdx=1π-∞∞e-z2dz=1π?π=1The mean of a continuous distribution function is given by:fx=abf(x)ρxdxwhere ρxis the probability distribution density function of x. Substituting appropriate values: x=-∞∞x?1σ2πe-x-μ22σ2dxLet z=x-μ:x=-∞∞z+μ?1σ2πe-z22σ2dz=1σ2π-∞∞z+μ?e-z22σ2dz=1σ2π?-∞∞z?e-z22σ2dz+μσ2π?-∞∞e-z22σ2dz=0 (Odd Function)+μσ2π?-∞∞e-z22σ2dzSubstituting back:x=μ-∞∞1σ2π?e-x-μ22σ2dx=μ-∞∞ρxdx=μPlease calculate x2, x-x2 and x2-x2x2=-∞∞x2?ρxdx=-∞∞x2?1σ2πe-x-μ22σ2dxLet z=x-μσ and dz=dxσ:x2=-∞∞μ+σz2?12πe-z22dz=-∞∞μ2+2μσz+σ2z2?12πe-z22dz=μ2-∞∞12πe-z22dz+2μσ-∞∞z?12πe-z22dz+σ2-∞∞z2?12πe-z22dzUsing:ddbIb=-∞∞y2e-by2dx=ddbπb=-12bπbσ212π-∞∞z2?e-z22dz=σ2 where b=12∴x2=μ2+σ2.For x-x2 and x2-x2:x-x2=x2-2xx+x2=x2-x2Using above results:x2-x2=μ2+σ2-μ2=σ2We have a drunken sailor, who moves in 1D with the step size of L. At each time period, she could randomly take a step forward with probability P+ = 0.5, take a step backward with probability P- = 0.25 or stay in position with the probability P0 = 0.25.Find the sailor’s mean position after Nth step.Find the variance of the position after Nth step.What is her diffusion constant D?Let Pk be the probability of taking steps of kL where k=0 or±1, and L is the step length. Pk is given by the following:Pk=0.50.250.25 k=1 k=0 k=-1Let the displacement of the jth step be ?xj=kjL and the initial position be x0=0. Let u be the mean value of kj: u=kj=kkPk=0.5×1+0.25×0+0.25×-1=0.25 and thus the mean position of the walker after the Nth step would be:xN=xN-1+LkN=xN-1+uL=NuL=0.25NLThe variance is given by:σ2xN=xN-xN2=xN-1+kNL-NuL2=xN-1-uN-1L+kNL-uL2=xN-1-uN-1L2+2xN-1-uN-1LkNL-uL+L2kN-u2=xN-1-xN-12+L2kN-kN2=VarxN-1+L2×Vark=NL2+Var(k)For Vark:σ2k=k-k2=k2-k2=kk2Pk-kkPk=0.5×12+0.25×02+0.25×-12-0.25=0.5The diffusion constant is given byD=L22?t×Vark=L22?t×0.5The chapter asserted that tiny objects stop moving at once when we stop pushing them.Consider a bacterium, idealized as a sphere of radius 1m, propelling itself at 1m/s. At time zero, the bacterium suddenly stops swimming and coasts to a stop, following Newton’s Law of motion with the Stokes drag force. How far does it travel before it stops? Comment.Our discussion of Brownian motion assumed that each random step was independent of the previous one; thus, for example, we neglected the possibility of a residual drift speed left over from the previous step. In the light of (a), would you say that this assumption is justified for a bacterium?Assume viscosity of water, = 0.001Pas, density of bacteria is same as water, = 1000kg/m3 and low Reynolds number where Stokes law assumption is reasonable. For a spherical bacterium of radius 1m moving at 1m/s:Mass of bacterium, m=ρV=1000×43π(1×10-6)3=4×10-15kgUsing Newton’s Second Law:F=mdvdtThe only force present is the viscous drag force which follows Stokes law: F=-6πμrvDoing a force balance:-6πμrv=mdvdt-t=0t6πμrmdt=v=0vdvv-6πμrmt=lnvv0v=v0?e-6πμrmtDistance travelled until stationary:v=dxdtv0?e-6πμrmt=dxdtt=0t=∞v0?e-6πμrmtdt=0xdxm6πμrv0=xx=4×10-156π(0.001)(10-6)10-6=2×10-13mThe bacterium travels about 20pm before stopping almost instantaneously as the travelling distance is almost 7 times in magnitude smaller than the bacterium.During the deceleration, we ignored any Brownian motion in the calculations as assumed that the flow due to propulsion is dominant. D=kBT6πμr=1.38×10-23×3006π(0.001)(10-6)=2×10-13m/sFrom the calculation of the diffusion coefficient and assuming that diffusion occurs, it is noted that in one second the distance moved by the bacterium would be same as the coasting distance. Thus we cannot ignore Brownian motion in our calculations. a) The initial position is given by f/k while the new equilibrium position is given by f1/k. Thus the displacement is:?x=f-f1kThe work done against the external force would then be:W=f1?x=f1f-f1kb) To find the maximal work, we need to optimize W, the work done:?W?f1=f-f1k-f1k=0The solution to the above equation is f=2f1 and the corresponding work done is:W=f1f-f1k=f2f-f2k=f24kThe useful work output is only half of what was stored in the spring. c) To make the process more efficient, the difference between the external force f1 and the original force fshould be infinitesimally small for every value of x. This would ensure that all the stored energy in the spring is retrieved without losing any to friction.For the polymer bead tied to the membrane, its microstate can be described by its position x and velocity v. Modelling it as a “spring”, the total energy would then be:ETx,v=12(mv2+kx2)Using the Boltzmann distribution, the probability is given as follow where C is a constant:Px,v=C×e-12(mv2+kx2)/kBTUsing the fact that the sum of probabilities is 1:PTx,v=C-∞∞e-12(mv2+kx2)/kBTdvdx=C-∞∞e-12mv2/kBT dv?-∞∞e-12kx2/kBT dx=1Using the following relation:-∞∞e-by2dy=πbPTx,v=C?2πkBTm2πkBTk=C?2πkBTmk=1C=mk2πkBT∴Px,v=mk2πkBTe-12(mv2+kx2)/kBTFor the average energy:ETx,v=-∞∞ETx,v?P(x,v)dvdx=mk2πkBT-∞∞12mv2?e-12(mv2+kx2)/kBTdvdx+-∞∞12kx2?e-12(mv2+kx2)/kBTdvdx=mk2πkBT-∞∞12mv2?e-12mv2/kBT dv?-∞∞e-12kx2/kBT dx+-∞∞e-12mv2/kBT dv?-∞∞12kx2?e-12kx2/kBT dx=mk2πkBT2πkBTm?m2?12π2kBT3m3+2πkBTm?k2?12π2kBT3m3 =kBT2+kBT2 =kBTTherefore at room temperature, the mean energy of the bead is:ET=kBT=1.38×10-23×298=4.14×10-21 JThe mean energy is related to the potential energy in a spring by the following equation:ET=12kx2Equating the two equations and solving:12kx2=4.14×10-2112k35×10-92=4.14×10-21∴k=6.8×10-6N/mThe hypothesis of the trap door suggests a two state system where the stiffness of the hair bundle can vary between two states. This is shown in the graph where the “active” hair bundle would behave like an ordinary spring with constant stiffness and thus we see a constant gradient for displacements of -80<x<-20 and 20<x<80. For displacements of -20<x<20, similarly a constant stiffness value was observed but this time it is negative. Thus we can attribute the observations to the trap door hypothesis where the hair bundle can vary between two states. The bump in the curve is rounded and not sharp likely due to cooperative effects of each hair in the bundle. Each hair would move individually until it becomes unstable, and then the entire bundle would jump to a thermodynamically stable position, resulting in a smooth curve and not a sharp jump.It is likely that displacement of 0nm is an unstable equilibrium point and thus any random displacement from the surrounding fluid would cause it to move towards the other two stable equilibrium points of -20 and 20nm. Evaluate:-∞∞e-p22mkBT?d3p=-∞∞e-px22mkBTdpx-∞∞e-py22mkBTdpy-∞∞e-pz22mkBTdzpz=-∞∞e-px2+py2+pz22mkBTdpxdpydpzThe osmotic pressure due to the solutes inside the bag is balanced by the osmotic pressure due to the 1mM salt solution outside. If 2mM salt solution is used, there will be a larger osmotic pressure inwards and this will result in the bags collapsing as water would seep out of the bag until the volume is reduced by half.Assuming that the salt completely ionized and is monovalent. AB(s)→A+(aq)+B-(aq)For a mole of salt, there would be two moles of ions. For glucose which does not ionize in water, two moles of glucose would be needed for every mole of salt to balance the osmotic pressure. The following relation can be used to calculate the average molar mass: p=cRT28 mmHg ×101325 Pa760 mmHg=3733 Pac=pRT=37338.314×303=1.48 mol/m3 M=60gL×1000L1m3×1m31.48 mol=40540gmol=40.5kDaFlux across capillary walls - since blood plasma proteins have depleted by 10%, the osmotic pressure should reduce by 10%:J=Lp?p=7×10-6×373.3Pa×1atm101325Pa=2.58×10-8cm/s Total flowrate: Q=JA=2.58×10-8cm/s ×1m100cm×250m2=6.45×10-8m3/s 6.45×10-8m3s×1000Lm3×3600shr×24hrday=5.57L/day Swollen bellies due to 5.57L of water going into the interstitial space per day and being retained in the body. 5.57L is almost equivalent to drinking 11 bottles of water. The depletion zone thickness is approximately equal to the size of small objects. If the radius of the protein is about 10nm, then the radius of the small object would probably be 0.1nm to 1nm. The total free energy reduction is given by the following relation:?F=ckBT×2R×A=VckBTWhen surfaces stick, the overlapping regions have a volume of about:V=10nm×10μm2=0.1μm3The free energy reduction is then given by:?F=VckBT=0.1μm3×c×kBTFor spherical proteins with volume fraction of 0.3:0.3=c×4π310-83c≈7×1022m-3∴?F=0.1μm3×7×1022m-3×kBT=7000kBTThis is significant!Using the Taylor Series for ex:ex≈1+x+x22!+x33!+…Taking the first two terms of the series expansion:eV≈1+V and e-V≈1-V Substituting: d2Vdx2=-12λD-2e-V-eV≈-12λD-21-V-1+V≈-12λD-2-2V≈λD-2VSolving the resulting linear second order differential equation: d2Vdx2-λD-2V=0Solution is given as:V=C1eλD-2x+C2e-λD-2x=C1ex/λD +C2e-x/λD Using boundary condition V()=0:V=C2e-x/λD If x?λD, potential will be very large but if x?λD, there would be almost no potential. Thus charges on the surface in a salt solution would have an influence over a distance λD, but beyond this border the surface and the cloud are essentially neutral.3.5 H-bonds per molecule:3.5×9kBTr=31.5×1.38×10-23×298=1.295×10-19J/moleculeEnergy for one mole of water:1.295×10-19J/molecule ×6.023×1023molecule/mol =78kJ/molComparing with the actual value of 40.7kJ/mol our estimated value is fairly close.As temperature rises, there would be an increase in molecular vibrations and this causes water to ionise more, forming more hydrogen ions which reduces the pH. Another way to explain this phenomenon would be that when a water molecule ionises to form two ions, the entropy increases. Since the chemical potential is defined as: μ=-TdSdN where in this case dS/dN is positive, a higher temperature will cause the chemical potential to be more negative, favouring dissociation. At equilibrium:Keq=H+[OH-][H2O]=e-?G0/kBTAs pH is defined as -log10H+:-log10H+∝-1/TThus an increase in T would lead to a decrease in pH, which results in an increase in the [H+]Concentration of acetic solution: 1mol10L=0.1MLet x be the concentration of acetic acid that ionizes:CH3COOH + H2O H3O+ + CH3COO-[CH3COOH][H3O+][CH3COO-]Initial0.100Change-x+x+xEquilibrium0.1-xxxKa=10-4.76=x20.1-xSolving: x=0.00131, thus the pH of the solution is -log10[0.00131]=2.88a) Consider a spherical protein with radius 3nm and an electrical charge of 10e. Force balance:Electrostatic force on particle = Fluid drag on particleqE=6πμrvv=qE6πμr=10×1.6×10-19C×200V/m6π×0.001kg/m?s×3×10-9m=5.66×10-6m/s b) ?v=q2-q1E6πμr=1.6×10-19C×470V/m6π×0.001kg/m?s×3×10-9m=1.33×10-6m/s =4.79×10-3m/hIf they ran the gel for 20hours, the difference in distance is about 4.79×10-3m/h×20h=0.1m = 10cm. This would allow for clear separation of the two bands. However if we just consider the distance travelled by the normal hemoglobin in 20h, it would be about 1m and thus the gel which they have made should be rather large. Running the gel for 20 hours consecutively would also generate a lot of heat but cooling system should have been in place for the experiment to work. For the vector function: rs=fs,g(s)=RcossR,RsinsRf's=-sinsR,g's=cossRThe arc length is given by:L=0sf'(s)2+g'(s)2ds=0s-sinsR2+cossR2ds=sThe unit tangent vector:ts=r'(s)r'(s)=-sinsR,cossR1=-sinsR,cossRThe derivative of the unit tangent vector:t's=-1RcossR,-1RsinsRThe bend is defined as:βs=dtds=-1RcossR,-1RsinsRThus the elastic energy cost for a 90 bend is:E=12kBT0142πRdsAβ2=12kBTA0142πRdtds2ds=12kBTA12πR1R2=πA4RkBTThe observed optical rotation of is a linear function of the fraction of amino acids in an -helix. Thus in an all-or-none model where the chain exists as either a random coil or a helix, we would be able to observe the optical rotation for the alpha-helix. By varying the number of subunits N and temperature where only the alpha-helix is observed, we would be able to observe the change in optical rotation. Letting the weight of the random coil be 1, the weight of the alpha helix conformation would then be:e-(N?G)/kBT =e-β(N?G)Thus the probability of being in the alpha helix state:e-β(N?G)1+e-β(N?G)=e-N?Ebond-T?Sbond-T?Sconf/kBT1+e-N?Ebond-T?Sbond-T?Sconf/kBT=e-N/kB?EbondT-?Sbond-?Sconf1+e-N/kB?EbondT-?Sbond-?SconfThe sharpness of the transition would increase with N. Since this is an all-or-none model, only two states of the chain would be present at equilibrium where either all residues are in random coil or in the alpha helix state. In this state of maximal cooperativity, the higher the number of residues, the easier the transition would be. Thus if there were an infinite number of residues, we would observe a step function. The hydrogen bond is being formed between the oxygen atom in the carbonyl group of monomer k and the hydrogen atom I the amide group on monomer k+4 when in the helical conformation. As such the last 2 residues are unable to form hydrogen bonding.Tm should generally decrease with increasing N.The diffusion coefficient is:D=kBT6πμr=1.38×10-23×3006π×0.001×0.37×10-6=5.936×10-13m2/sFor 2D diffusion:<x2>=4Dt=4×5.936×10-13×30=7.12×10-11m2 <x>=4Dt=7.12×10-11=8.44×10-6m≈8.4μmSince he calculated the root-mean-square displacement to be 7.84m, the measured value is within error of the predicted value.The probability that the displacement of the particles has a magnitude between r and r+dr is:Pr,tdr=14πDtexp-r24Dt2πrdr=rdr2Dtexp-r24DtThe number of particles within the first circle should therefore be:N007.84×10-64Pr,tdr=508?12Dt07.84×10-64r?exp-r24Dtdr≈26.39Counting the number of particles in the first circle, there are about 30 and thus the prediction is in line with the actual number. Similar calculations can be done for the larger rings by replacing the integration limits with d/4 and 2d/4 etc… ................
................

In order to avoid copyright disputes, this page is only a partial summary.

Google Online Preview   Download