Using step function in kinetic reaction simulation

Topics: Microfluidics, 4.2a

Thread index  |  Previous thread  |  Next thread  |  Start a new discussion

RSS FeedRSS feed   |   Email notificationsTurn on email notifications   |   5 Replies   Last post: January 29, 2013 9:01am UTC
Pattamon Teerapanich

Pattamon Teerapanich

January 25, 2013 2:26pm UTC

Using step function in kinetic reaction simulation

Hi

I am trying to simulate the kinetic reaction of biomolecule interaction in fludic channel expecting to have association and dissociation phase (concentration of complex versus time) in the same curve. I knew that it is possible to use step function or smoothed Heaviside step function from COMSOL in order to change the concentration of analyte species at the inlet from C=x mol/m3 to C=0 mol/m3 at a certain time. I have no idea how to find this function and what kind of variables I need to defind. Anyone has any experience on this?

Regards,

Pattamon

Reply  |  Reply with Quote  |  Send private message  |  Report Abuse

Ivar Kjelberg

Ivar Kjelberg

January 26, 2013 12:38am UTC in response to Pattamon Teerapanich

Re: Using step function in kinetic reaction simulation

Hi

in v4 you find these under root-Globa-Definitions, or better: Model-Definitions-Functions, use the dynamic doc and search the HTML or PDF doc

--
Good luck
Ivar

Reply  |  Reply with Quote  |  Send private message  |  Report Abuse

Pattamon Teerapanich

Pattamon Teerapanich

January 28, 2013 2:21pm UTC in response to Ivar Kjelberg

Re: Using step function in kinetic reaction simulation

Hi,

I have tried to put the step function where I want to change the inlet concentration from c0 to 0 at time = 400s.
I defined the step function for concentration in the variables. It didn't not work for now. I guess my function is not correct or I need to add something more. I also attached my short simulation model here.

Thank you for your advice.

Pattamon

Reply  |  Reply with Quote  |  Send private message  |  Report Abuse

Ivar Kjelberg

Ivar Kjelberg

January 28, 2013 8:12pm UTC in response to Pattamon Teerapanich

Re: Using step function in kinetic reaction simulation

Hi

first of all your shape is very thin, with only 1 element across you do not resolve a lot, try to scale the vertical axis by i.e 100 in the mesh advanced tab.

Then your process is over in 4 seconds, your solver steps seem to be far too long, what about cutting your rectangle with a layer (geometry tab) at te half height of 0.3 um, then solve and look at the concentration along this central line in a line plot

--
Good luck
Ivar

Reply  |  Reply with Quote  |  Send private message  |  Report Abuse

Magnus Ringh

Magnus Ringh
COMSOL Employee
Sweden

January 29, 2013 7:57am UTC in response to Pattamon Teerapanich

Re: Using step function in kinetic reaction simulation

Hi,

The step function uses the concentration as the input argument (the switch), but that should be the time t. Also the step amplitude is the concentration c0, so I don't see why that value is multiplied with the step output in the variable intl: c0*step1((c0-0)[m^3/mol])

I think that you should use:

step1(t[1/s])[mol/m^3]

where "[1/s]" is used to make the argument unitless and the unit of the output is mol/m^3 for a concentration.

Also, perhaps the initial value should be c0 rather than 0? Now the initial value and the boundary condition for the concentration don't match at the start of the simulation.

Best regards,
Magnus Ringh, COMSOL

Reply  |  Reply with Quote  |  Send private message  |  Report Abuse

Pattamon Teerapanich

Pattamon Teerapanich

January 29, 2013 9:01am UTC in response to Magnus Ringh

Re: Using step function in kinetic reaction simulation

Hi

Thank you so much for all advices. It works now :)

Best Regards,

Pattamon

Reply  |  Reply with Quote  |  Send private message  |  Report Abuse


Rules and guidelines