Hi Jonathan, The temperature change predicted is simply too high. This simulates sliding of a granular material under modest stress at 10 µm/s. I did find the issue though. For what it’s worth the representation should be:
eqX = TransientTerm() == ExplicitDiffusionTerm(coeff=D) + (mask * ((mu*sigma_n*v/w) / (rho*cp))) eqI = TransientTerm() == DiffusionTerm(coeff=D) + (mask * ((mu*sigma_n*v/w) / (rho*cp))) eqCN= eqI + eqX I was just worried I didn’t grasp the way I needed to interact with FiPy. An example along this vein could be a useful addition in the docs possibly? Cheers, John > On Apr 19, 2016, at 10:27 AM, Guyer, Jonathan E. Dr. (Fed) > <[email protected]> wrote: > > I don't see anything obviously wrong with the way you've posed the source. > Why do you think you have a problem? > >> On Apr 18, 2016, at 4:44 PM, John Leeman <[email protected]> wrote: >> >> Hi all, >> >> I’m working on implementing a rather simple thermal diffusion model with >> FiPy, but haven’t been able to find any examples or things in the docs to >> get me over the last challenge. >> >> The model represents a stack of several materials. One of them is a granular >> substance that is sheared and generates some heat. The shear heating >> generation is just the product of the shear stress and velocity. The >> velocity will vary with time in the model (in the end). Basically I’m >> solving equations 2,3 from >> http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.211.384&rep=rep1&type=pdf >> >> I was able to build the model and solve a simple diffusion problem, but I’m >> having trouble figuring out how to get the heat generation term into FiPy. >> You can see the simple notebook and a notebook in which I tried to add the >> heating terms on this gist: >> https://gist.github.com/jrleeman/8bcafffd512feaa77ec2946dc92460ae >> I’m probably missing something, but after looking through the source term >> docs I was unclear how to proceed. >> >> Many Thanks, >> John >> _______________________________________________ >> fipy mailing list >> [email protected] >> http://www.ctcms.nist.gov/fipy >> [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] > > > _______________________________________________ > fipy mailing list > [email protected] > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
