Thanks for your detailed answer. PSB On Thu, Apr 30, 2009 at 6:20 PM, Daniel Wheeler <[email protected]>wrote:
> > On Wed, Apr 29, 2009 at 8:34 AM, david wende <[email protected]> > wrote: > > Hi > > Excuse my "list oriented" nature. > > > > 1 Thanks for what seems to be great software. > > Thanks for your interest. > > > 2. This is my first forray into F.E. analysis. > > Good luck! This is an FV code, but the distinction is rather pedantic. > > > 3. Need to simulate pulsed laser directed into thin silicon wafer > > and the resulting time varying heat gradient. > > 4. mayavi 2 doesn't work (known problem) > > Mayavi 2 works, it's just we don't have an interface to it right now. > However, you can extract the required data arrays and pass them to > mayavi2 like any other numpy arrays, which is the main reason we have > wrapper classes to take care of this. If you want to make a > publication quality figure or customize your figure, you would > probably have to hack our viewers or interface with Mayavi2 directly > anyway. OK - I'll try and extract array and view with mayavi2 How would I dump data into format suitable for mayavi2? I tried TSVViewer but mayavi2 does not like it. > > > > 5. mayavi 1.5 not working on my system > > What happens? Does it freeze on the screen? Are you using Windows? Mayavi 1.5 on Linux. Error messages on startup, such as: " _tkinter.Tcl_Obj has no len() " > > > > 6. So - could I do a 3D simulation and then slice a cross section and > send > > that to a 2D viewer? > > Yes, extract the data, create a 2D mesh and a new CellVariable based > on the 2D mesh and initialized with the sliced data. This sounds hard - I will try first suggestion first (data dump -> mayavi2). Do you mean something like: solution2D = solution3D(where x = 3) where solution3D is the solution variable in the eqn.solve step ? > > > > 7. The pulse of laser I would (I think) model as a time and spatially > > dependent source - could > > someone send me an example of how to do that? > > Here is a code snippet that might be helpful > > t = Variable() > x, y = mesh.getCellCenters() > sv = CellVariable(mesh=mesh, x * y) > eqn = SomeTerms + sv * t > dt = 1.0 > > for step in range(steps): > eqn.solve(var, dt=dt) > t.setValue(t + dt) > > The important thing is to remember that the time variable needs to be > updated and that the source needs to be a CellVariable object. > Unfortunately, getCellCenters does not return x and y as CellVariables > right now, which adds an additional line of code to this example. > Hopefully, this will change in future releases. I think that I understood the above about the time-dependence - thanks. > > > Cheers > > -- > Daniel Wheeler > > -- David Wende home +972-8-9353488 work +972-2-5886116 mobile +972-54-234-6479
