On Jan 22, 2008 1:41 PM, asreeve <[EMAIL PROTECTED]> wrote: > > Hello,
Hi, hope all is well. > I've having some difficulty figuring out how to use the convection terms > in a 3-D problem. I haven't solved a lot of 3D convection problems so I am not sure how well tested everything is. Should work though. > I've calculated a velocity term for a 3-D simulations: > vel=dot(head.getFaceGrad(),kond.getHarmonicFaceValue()) I'm confused about which axis is getting dotted. What are the shapes of "kond" and "vel"? > and would like to pass these velocities to a convection term. I assume that "vel" is rank 1. Is it? > I've simplified my problem to the point that there are non-zero velocities > only in the z direction (but I need to solve a problem where vel varies in > space and direction). Shouldn't be a problem. > > I've tried using : > > PowerLawConvectionTerm(coeff=vel()) > > and > > PowerLawConvectionTerm(coeff=(vel()[0],vel()[1],vel()[2])) > > in the equation term, but both of these yield an incorrect result. > > However, when I assign the following (values taken from vel variable), > > PowerLawConvectionTerm(coeff=(0,0,2.5e-4)) > > the output is correct. > > What's the proper syntax to assign the calculated 'vel' to a convection > term? Have you tried? PowerLawConvectionTerm(coeff=vel) Cheers -- Daniel Wheeler
