I'm trying to figure out how to load in a csv file to implement some time
dependent BCs. We are trying to model contrast agent concentration within
tissue, and the mesaurements we have are concentration over time in the
blood vessels. Lets call this variable Cp. Contrast agent cannot cross into
cells, and therefore the flux at those boundaries is 0. Contrast can come
through vessels, through an empirically derived formula shown below. For
the boundary, I have written these two conditions:
[BCs]
[./cell]
type = NeumannBC
variable = diffused
value = 0
[../]
[./vessel]
type = NeumannBC
variable = diffused
value = _normals*P*(Cp_t-_u_t)
[../]
[]
Where, P is a constant and Cp is read in through the Functions block. Could
someone help me out and show me what I'd need to do to read in the Cp.csv
in the Functions block? It only has two columns: time and conc. Ideally,
I'd like to be able to interpolate between the values in this file, since
I'll be running the FEM model at a smaller resolution than the data was
taken. Also, could someone help me choose the correct BC type for the
vessel BC? I'm pretty certain this format is not correct since it changes
with time. Does it need to be FunctionNeumannBC, and have everything
defined in the Functions block, instead of the BCs block?
Thanks!
Ryan
------------------------------------------------------------------------------
What NetFlow Analyzer can do for you? Monitors network bandwidth and traffic
patterns at an interface-level. Reveals which users, apps, and protocols are
consuming the most bandwidth. Provides multi-vendor support for NetFlow,
J-Flow, sFlow and other flows. Make informed decisions using capacity planning
reports.http://sdm.link/zohodev2dev
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
