On Thu, 31 Jul 2008, luyi wrote: > Some one had asked you about the multi-physical problems using libmesh, > and you tell us you are working the conjugate heat transfer problem. > > As the question "Associating more than one mesh to EquationSystems", I > want to perform fully implicit coupling(without operator splitting) and > use the Newton-Krylov method to solve conjugate heat transfer problem, > the difference is that I am not using stagger mesh, only one mesh > departed to domain A and domain B, the variables in domain A are > rho,u,v,w,E,k,w and in domain B only T, and I want to assemble the > Jacobian matrix and the rhs together, can I using only one > EquationSystems and one mesh to implement this?
At the moment I think the answer is "not efficiently". Although our meshes can be flagged with subdomain ids, our DofMap still assumes that every subdomain in the mesh has the same variables defined on it. What we need is an API for the user to say "this variable only exists on subdomains 1, 5, and 13", as well as the DofMap code to respect those requests. Unless you want to dig into the library source code and add that feature (for which a patch would be much appreciated, naturally), I think the best you can do for a fully coupled solve is define all variables everywhere, then manually restrict T to be 0 on domain A and all the other variables to be 0 on domain B. --- Roy ------------------------------------------------------------------------- This SF.Net email is sponsored by the Moblin Your Move Developer's challenge Build the coolest Linux based applications with Moblin SDK & win great prizes Grand prize is a trip for two to an Open Source event anywhere in the world http://moblin-contest.org/redirect.php?banner_id=100&url=/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
