Re: [Libmesh-users] Coupling variables with different dimensions

2012-02-08 Thread Roy Stogner
On Wed, 8 Feb 2012, David Knezevic wrote: > On 02/08/2012 06:09 PM, Mauro Werder wrote: >> I've never done it but here my thoughts: It should be possible to >> assemble the 1D system when looping over the boundary elements of the >> 2D mesh (as is usually done for the BCs) as all the machinery is

Re: [Libmesh-users] Coupling variables with different dimensions

2012-02-08 Thread David Knezevic
On 02/08/2012 06:09 PM, Mauro Werder wrote: > I've never done it but here my thoughts: It should be possible to > assemble the 1D system when looping over the boundary elements of the > 2D mesh (as is usually done for the BCs) as all the machinery is in > place to do this. All the machinery is the

Re: [Libmesh-users] Coupling variables with different dimensions

2012-02-08 Thread Mauro Werder
I've never done it but here my thoughts: It should be possible to assemble the 1D system when looping over the boundary elements of the 2D mesh (as is usually done for the BCs) as all the machinery is in place to do this. Mauro At Wed, 8 Feb 2012 22:07:16 +, Sylvain Vallaghe wrote: > > Hi, >

Re: [Libmesh-users] Q about libMesh features

2012-02-08 Thread Ali Roustaei
Thanks Jed! very nice explanations Best Regards, Ali -- Keep Your Developer Skills Current with LearnDevNow! The most comprehensive online learning library for Microsoft developers is just $99.99! Visual Studio, SharePoi

Re: [Libmesh-users] Coupling variables with different dimensions

2012-02-08 Thread David Knezevic
On 02/08/2012 05:07 PM, Sylvain Vallaghe wrote: > I'm trying to implement a 2d problem in libmesh where I have two scalar > variables : one is 2d in the whole domain, the other is 1d and mapped on one > of the boundaries of my domain. These two variables are coupled due to the > global variati

[Libmesh-users] Coupling variables with different dimensions

2012-02-08 Thread Sylvain Vallaghe
Hi, I'm trying to implement a 2d problem in libmesh where I have two scalar variables : one is 2d in the whole domain, the other is 1d and mapped on one of the boundaries of my domain. These two variables are coupled due to the global variational formulation. I don't know how to define my 1d va

Re: [Libmesh-users] Q about libMesh features

2012-02-08 Thread Jed Brown
On Wed, Feb 8, 2012 at 20:38, Roy Stogner wrote: > > 3- Support for an Uzawa solver (stokes problem) > > This is doable in libMesh application code, but it wouldn't be as easy > as many other solvers, and it isn't built in to the library. > For simple Stokes problems, PETSc can construct this au

Re: [Libmesh-users] problem with EquationSystems in a loop

2012-02-08 Thread Roy Stogner
On Wed, 8 Feb 2012, water force wrote: > Yes this algorithm seems uncommon because the material properties > associated with these meshes are different. Is there any other > elegant way to do this? It should have been obvious from your code, but no - if your material properties affect the Jacobi

Re: [Libmesh-users] Q about libMesh features

2012-02-08 Thread Ali Roustaei
Thanks for info Best Wishes, Ali -- Keep Your Developer Skills Current with LearnDevNow! The most comprehensive online learning library for Microsoft developers is just $99.99! Visual Studio, SharePoint, SQL - plus HTML5

Re: [Libmesh-users] problem with EquationSystems in a loop

2012-02-08 Thread water force
Thanks for your comments! A bug in my assemble function did result in the segfault. Actually I used attach_assemble_object in my code but made a mistake when passing the string sys_name to the constructor in the assembly class. Yes this algorithm seems uncommon because the material properties assoc

Re: [Libmesh-users] Q about libMesh features

2012-02-08 Thread Kirk, Benjamin (JSC-EG311)
Nothing particularly tricky about the FVM interface so long as you are (i) cell centered and (ii) interested only in nearest face neighbors. In that sense it is a natural subset of the DG support. If you want a node centered scheme with dual-mesh control volumes there will be many tricks - I wo

Re: [Libmesh-users] Q about libMesh features

2012-02-08 Thread Ali Roustaei
Thanks you very much Roy, For periodic BC that's enough for me. For FVM computations do we have a clean interface? or it is done with tricks? I guess this would be the more time consuming part Regards, Ali -- Keep Your

Re: [Libmesh-users] Q about libMesh features

2012-02-08 Thread Roy Stogner
On Wed, 8 Feb 2012, Derek Gaston wrote: On Feb 8, 2012, at 10:38 AM, Roy Stogner wrote: libMesh can do FVM with small flux stencils by using discontinuous shape functions for those variables; for larger stencils you'd be out of luck, since libMesh wouldn't be able to extend the parallel ghost

Re: [Libmesh-users] Q about libMesh features

2012-02-08 Thread Derek Gaston
On Feb 8, 2012, at 10:38 AM, Roy Stogner wrote: > libMesh can do FVM with small flux stencils by using discontinuous > shape functions for those variables; for larger stencils you'd be out > of luck, since libMesh wouldn't be able to extend the parallel > ghosting or sparsity pattern as far as you

Re: [Libmesh-users] Q about libMesh features

2012-02-08 Thread Roy Stogner
On Wed, 8 Feb 2012, Ali Roustaei wrote: > I'm looking for a flexible computational environment, currently my needs > right now are: > > 1- Support for P1isoP2 elements, P0,P1,P2 and Q0,Q1,Q2 and Q0isoQ2 > 2- Support for discontinuous elements P1disc, Q1disc All of this is easy enough in libMes

[Libmesh-users] Q about libMesh features

2012-02-08 Thread Ali Roustaei
Hi, I'm looking for a flexible computational environment, currently my needs right now are: 1- Support for P1isoP2 elements, P0,P1,P2 and Q0,Q1,Q2 and Q0isoQ2 2- Support for discontinuous elements P1disc, Q1disc 3- Support for an Uzawa solver (stokes problem) 4- Periodic Boundary Conditions 5-