-
Filter by Topic
Release
Popular Topics
- List all discussions
LinkedIn
Facebook
Twitter
Delicious
Digg
Stumbleupon
More services...Using step function in kinetic reaction simulation
|
Thread index | Previous thread | Next thread | Start a new discussion |
RSS feed |
Turn on email notifications |
5 Replies
Last post: January 29, 2013 9:01am UTC
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
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
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
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
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
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
