CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   CFX (http://www.cfd-online.com/Forums/cfx/)
-   -   FSI problems (http://www.cfd-online.com/Forums/cfx/94939-fsi-problems.html)

lingdeer November 30, 2011 12:33

FSI problems
 
Hi, I am simulating blood flow through blood vessel and have been facing some problems that I can't reference to in the CFX manual documents.

1) The first time I used coupling convergence of 1e-2, source term of 300 (This is from trial and error until my simulation won't crash). However when I monitor two supposed-to-be symmetric points, it gives me unequal mesh displacement. Also, the solution fail to converge after running for 30 cycles (for each cycle I imposed sinusoidal velocity inflow profile. Outlet I imposed static pressure = 0Pa). Fluid convergence = RMS 1e-5
Therefore, I realize I need to tighten my coupling convergence to lower value to improve convergence.

Now, I changed to coupling convergence to 1e-3, the two monitor points now have equal values. And the mesh displacement waveform of those two points make more sense. Also, since they are at a junction, it will have the biggest displacement in the entire domain. These points are downstream of the flow, so biggest mesh displacement happen later in the cycle. However, when I run with this new coupling convergence, it crashed at the 60th timestep (each cycle time = 1s, time step = 0.01s), which is the point it reaches the biggest mesh displacement and start to retract to its original position.

Since it crash at the middle of the simulation, it's hard to check if it the source term/parameters are good at the beginning of the simulation.

I wonder for anyone who simulate blood flow in blood vessel with CFX can give me some hints on the parameters that I should use to ensure it run smoothly. e.g. what source term do you use? coefficient loop?

2) FSI simulation occupies diskspace. Until it reaches transient convergence (low cycle-to-cycle variation), I don't need results of earlier timestep, but only the last few cycles. So for fluid I know I can delete the previous .trn files to save diskspace. For solid, (.ansys folder), anyone knows how to save space? What files are necessary to keep are what are not (fluid and ansys)?

Thanks Thanks Thanks!

lingdeer November 30, 2011 12:45

P.S. the reason of the crash is because of thickness change too big, extreme distortion in the solid solver, which shows at fatal overflow error in the CFX fluid part.

stumpy November 30, 2011 16:42

For 1) posting the monitor plots showing force/displacement within each timestep close to the crash would be useful.
For 2) you need to keep all the files while it's running. After the run finishes you can delete the .rst file if you don't want to view the results. You can use MFOU and OUTRES to reduce the frequency and amount of data that is written to the rst file, but you'd need to stop and restart to increase the frequency or data output once you reached the point of interest. Make sure you write at least:
OUTRES,NSOL,LAST
otherwise you won't be able to restart at all.

lingdeer December 4, 2011 18:29

Thanks Stumpy

1) After it crashed I don't know why I cannot view the monitor points anymore.. Maybe because the .res is not created? (I am running in background so I don't have it opened till it crashed)

2) I have serious trouble with restarting FSI cases.
The ideal case for me is to backup fluid and solid at every timestep instead of coupling step, so I can start at the timestep before something crashed.

I know:
MFOU,10
MFRC,50,1

means writing ANSYS results every 10 time steps
and MFRC will clear the results every 50 time steps and save the new ones again.

The thing is I don't write fluid results every coupling step but every time step instead to save space.

Would the command you suggested write at every time step or every coupling step?

I read from the command reference that 'LAST' means every last substep of each load step. Here is each load step means coupling step or time step?

What if I want to save ANSYS results every time step instead of every coupling step? MFOU and MFRC maybe useful, but what if I am not certain who many coupling steps it would actually run in each time step?

Thanks!

stumpy December 5, 2011 12:16

A time step and a coupling step are the same thing. Perhaps you're mixing up coupling iteration with coupling step?
MFOU write data to the rst file, as you said.
MFRC creates backup points and does not remove or add any data from the rst file. A backup point in Mechanical writes data to r001, r002, etc files. If you only keep 1 backup point then you'll only get the r001 file. When restarting data from the rst file is also needed, so when MFRC writes data you have to make sure MFOU is also writing data at the same timestep. So if you want to backup points every timestep you would need:
MFOU,1
MFRC,1,1
then in CFX ask for backup files every timestep. I never use the coupling step frequency options in CFX for reasons I won't get into now. Every timestep is probably a bit excessive, and you might already to too close to failing to have "good' results. Backups every 10 timesteps would be:
MFOU,10 (or 5 or 2 or 1 if you want more frequent data for post-processing)
MFRC,10,1
Then in CFX ask for every 10th timestep.
For FSI a load step = coupling step = time step. So LAST means once per time step, but this gets over-ridden by MFOU.

lingdeer December 5, 2011 14:30

Worked! Thank you so much! :)


All times are GMT -4. The time now is 08:22.