Reset grid in time dependent problems

General discussions about how to formulate a script for FlexPDE.

Reset grid in time dependent problems

Postby jlsanhueza on Wed Jan 02, 2019 12:29 pm

Hello users and moderators,

I'm working in fluid mechanics with moving meshes to solve the Navier-Stokes equations and my main goal is to "reset" the grid every time step or after the script is done, so I can run another script loading the final boundary geometry but with a new mesh configuration.

In the attached script, there are 2 fluids with different viscosities and densities (that's why I added a buoyancy gravity term to a particular region). At the end of the run I want to preserve the boundaries (mainly the region 2 - yellow one) but reset the mesh configuration. The reason is because I want to study how the region 2 travels along the region 1. However, when time becomes large enough (increasing the T value in the code) some nodes collapse, that's why reseting the mesh would be useful to adress this problem.

I've tried TRANSFERMESH, RESOLVE and FRONT commands with no success at all, altough I could transfer the whole boundaries+mesh configuration at the final stage, but I couldnt re-mesh from that. Also, I haven't tried STAGES yet since I'm getting a "MEMORY PROTECTION FAULT" if I remove the time dependence.

Any suggestions?

Thanks in advance to everyone.

Notes:
1. I'm using 6.5 professional version.
2. Some units in the script are meaningless.
3. I want to use this code to solve geological problems (deformation controlled by creep), so my viscosities are up to ~10^20 [Pa s] and time ~10^13 seconds (millions of years).
Attachments
two_fluids_test.pde
(984 Bytes) Downloaded 13 times
jlsanhueza
 
Posts: 1
Joined: Thu Dec 27, 2018 1:30 pm

Re: Reset grid in time dependent problems

Postby moderator on Mon Jan 07, 2019 4:15 pm

The mesh motion equation div(grad(xm))=0 is a generic kind of motion equation that works in some problems, but in others (like yours) it allows the mesh to smash at the front of a moving object.
There are other equations that can be used to move the mesh that prevent that kind of compression.
Unfortunately, I have meetings this afternoon, and will have to address this in more detail tomorrow.
moderator
 
Posts: 749
Joined: Tue Jan 11, 2011 1:45 pm

Re: Reset grid in time dependent problems

Postby moderator on Wed Jan 09, 2019 2:03 pm

One approach in fluid problems is to use a Lagrangean formulation, and lock the mesh motion to the fluid motion.
Since you are already calculating fluid velocities, you can drag the mesh along with the fluid simply by providing the mesh motion equations
Xm: dt(Xm) = vx
Ym: dt(Ym) = vy

Unfortunately, you have conflicting boundary conditions which prevent this from working correctly.
You have provided both a uniform pressure and a uniform velocity on the left boundary, along with a fixed upper and lower boundary velocity.
A uniform pressure does not create a uniform velocity when there is drag on the walls. The velocity will tend toward a parabolic shape in the channel.
Also, forcing both a zero velocity and a non-zero velocity at the meeting of two boundaries induces an unresolvable shear at the corner point.
moderator
 
Posts: 749
Joined: Tue Jan 11, 2011 1:45 pm


Return to FlexPDE Discussion

Who is online

Users browsing this forum: No registered users and 1 guest