Ask an Expert Help-24/7

ASK QUESTION & GET ANSWERED ONTIME!

High Quality, Fast Delivery, Plagiarism Free - Just in 3 Steps

Submit Your Question Here

Email Address
Word Limit
No. of Pages
(yy-mm-dd)
Multiple Attach File
Professional Assignment Help Services

MATLAB Assignment Based on Numerical Methods in Biomedical Engineering

Book Your Assignments Today!
&
Securing Higher Grades at Least Amount Of Effort?

Get Assignment Answers from Top-notch Tutors – specialize in your subject areas

Hire Your Assignment Writing Expert for  Collage/University

Question 1: [7 points]

Stryker produces three different knee replacement implants (“Triathlon”, “Restoris”, and “Tritanium”), each with different numbers of components inside. https://www.stryker.com/content/stryker/us/en/portfolios/orthopaedics/joint-replacement/knee.html

Stryker produces1

The table below shows the number of each component (pins, screws, and hinges) required to build each model.

Component
Model name: Pins Screws Hinges
Triathlon 4 4 2
Restoris 6 2 1
Tritanium 5 3 3

 

In one week, 1110 pins, 690 screws, and 480 hinges are shipped to the company.

(a) Write down a system of equations that model the production run in both conventional and matrix form [4 points]

(b) Write a Matlab script to calculate how many of each prosthetic knee model the company can produce that week. Compare the solutions obtained for two direct solution methods (use built-in Matlab functions). [2 points]

(c) How long does your code take to find a solution when each of these Matlab direct solution methods are used? [1 points]

Question 2 [12 points]:

12 points-2

The process to extract purified penicillin from a fermentation broth uses a series of reaction chambers with the raw broth solution flowing in the opposite direction to a solvent extraction stream (amyl acetate).  The goal is to separate the penicillin from the broth by extracting the drug in a pure solvent.  In each reaction chamber, the broth and the solvent do not mix, but some of the penicillin diffuses into the solvent. The exiting solvent is called the extract.

In the extraction system shown below, the input flow rate of solvent (amyl acetate) is S cm3/hr and of penicillin broth is F cm3/hr. The concentration of penicillin entering stage “n” in the feed stream is xn+1 g/cm3 and the concentration of penicillin entering stage “n” in the extraction stream is yn-1 g/cm3.

12 points-2-1

(a)  [2 points] Assume that the mass flow rates (S, F) don’t change over the length of the reaction system. Use steady state mass balance (mass in per hr = mass out per hr) to write an equation for each stage when we have n = 6 stages.

(b)  [3 points] Use the distribution coefficient K = yi / xi in each stage to eliminate the yi term in each of your equations in part (a). Express these equations in matrix form, [A][x] = [b].

What physical part(s) of the system is represented by the [b] vector?

Hint:  Gather terms in “x” on the left side of your equations in (a) in order x1, x2, x3, … and move constant terms to the right side. Matrix [A] will contain combinations of S, F, and K. Column vector [x] will contain the unknown x1, x2, x3,…

(c) [2 points] Use one of Matlab’s direct algorithms to determine the concentration of penicillin entering each stage in the feed stream when S = 125 cm3/hr, F = 265 cm3/hr, xf = 0.028 g/cm3, ys = 0 g/cm3, and K = 2. Plot these

concentrations as a function of stage #.

(d) [2 points] How many iterations does Matlab’s iterative solution algorithm lsqr need to reach a solution with a relative residual (error) below 1 x 10-10?

(d) [3 points]  As the number of stages in the extraction system increases, more of the penicillin is transferred over to the extract stream. How many stages are needed to achieve an extract concentration above 0.051 g/cm3 under the same conditions as (c)?  Plot the concentration of penicillin in the extract stream as a function of stage # for this case.

Question 2: [10 points]

The Krogh tissue cylinder is a simplified model for evaluating transport of metabolites from the capillary vessels to the surrounding tissue.   The Krogh model geometry is shown below.   The capillary is modeled as a cylinder of constant radius (rc), with wall thickness (tm), oriented in length along the positive z-axis.  As a solute moves through the capillary lumen, its concentration decreases due to radial diffusion of molecules through the capillary wall and into the surrounding tissue.

10 points-1

At some radial distance into the tissue (rcrit), the concentration of any metabolite will have reduced to zero. Solution of the Krogh tissue cylinder model leads to a nonlinear expression describing rcrit as a function of distance z along the capillary:

 

10 points-2
(i) Write a Matlab script to calculate the critical tissue radius (rcrit) where glucose is no longer supplied to cells, at a distance z = 0.025 cm under the conditions listed in the table.  Use the Newton-Raphson method with ea = 0.1%. Generate plots showing your root estimates and approximate relative error at each iteration.

Hint:  First find R2, then use this to calculate rcrit.

[6 points]

(ii) Modify your Matlab script from (i) to generate a plot of rcrit versus z as z ranges from 10 mm to 0.15 cm.

[4 points]

 

======

é

R 2 ln(R 2 )- R 2  + 1 – ê

4DT C0

ù        D

ú + 4     T

[R 2  – 1]z + 2DT

[R 2  – 1] = 0

 

2
2

ê R (rc  tm )  ú

Vrc

rc K 0

 

 

 

 

R =    rcrit

ë   0                                û

 

where

.

rc  tm

This can be written in more compact form as:

 

 

R2  ln ( R2 ) = A E R2  -1)

 

 

 

A =      4DT C0

B = 1 – 2DT

D = 4DT

E = ( B Dz )

 

c

with

R0 (rc

t m )

rc K 0

Vr 2

 

 

2

(i) Write a Matlab script to calculate the critical tissue radius (rcrit) where glucose is no longer supplied to cells, at a distance z = 0.025 cm under the conditions listed in the table.  Use the Newton-Raphson method with ea = 0.1%. Generate plots showing your root estimates and approximate relative error at each iteration.

 

Hint:  First find R2, then use this to calculate rcrit.

 

[6 points]

 

(ii) Modify your Matlab script from (i) to generate a plot of rcrit versus z as z ranges from 10 mm to 0.15 cm.

 

[4 points]

 

 

2017 Numerical Modeling – Lecture 7

2017 Homework 2 2017

Numerical Modeling – Lecture 6

Get This Answer with 24/7 Subject-Matter Experts

Get help instantly with No1 Assignment Help and Earn better grades
with 24/7 homework help

Chat Now

Top