Re: [Libmesh-users] Defining two mesh objects (mesh_pos+ mesh_negative)

[Xiao Ma <xi...@il...>]

Hi Roy,


Thank you for your reply.
I think I will go with two separate equation systems one on each mesh.
The problem I am trying to solve is a two half plane and in between it is
governing by a friction model, and the friction contact is enforce by a
method "TSN" (Traction Separation nodes).
Say I have two half plane (mesh_pos,mesh_neg),  so the mesh_pos 's bottom
layer is gonna to be applied some traction , and the mesh_neg' top layer is
gonna be applied some traction to enforce the jump condition and also the
friction law.

So what I need is to have these two separate EquationsSystems(System for
mesh_pos, System for mesh_neg), two parts are solved separately but at the
same time step, the only interaction is through the traction at the
separate nodes.


Best,
Xiao

On Mon, Oct 30, 2017 at 9:53 AM, Roy Stogner <roy...@ic...>
wrote:

>
> On Sun, 29 Oct 2017, Xiao Ma wrote:
>
> > I am wondering if there is a way to define two mesh objects , which each
> > represents half of the domain, and solve the two part at the same time ,
> > mesh_pos and mesh_neg , two meshes are separated, they are not connected
> .
>
> Depends what you mean by "at the same time".  If you want to solve a
> fully coupled problem, then you can use
> ReplicatedMesh::stitch_meshes() to combine them into a single mesh,
> and solve the problem on that.  If you want to solve a weakly coupled
> problem, then you can create a System on each mesh... either as part
> of the same EquationSystems but using subdomain-specific variables, or
> as part of two separate EquationSystems, one on each mesh.
> ---
> Roy
>



-- 
-------------------------------------------------------------------------------------
*Xiao Ma*
Graduate Research Assistant
PhD student in Structural Engineering
University of Illinois at Urbana-Champaign
E-mail: *xi...@il... <xi...@il...>*
--------------------------------------------------------------------------------------