← Back to team overview

yade-users team mailing list archive

Re: [Question #593653]: How to simulate a sinusoidal-shaped pressure condition at the top boundary of the packing

 

Question #593653 on Yade changed:
https://answers.launchpad.net/yade/+question/593653

Bruno Chareyre posted a new comment:
Hi,
1] There is no strong notion of top/bottom in the code design although most default settings are assuming the order x-,x+,y-,y+,z-,z+. In details the ordering is defined by flow.wallIds.
So, if  flow.bndCondIsPressure=[0,0,0,1,0,0] result in a change at the bottom I guess flow.wallIds[3] correspond to the id of the bottom Box.

2] It should be enough to set 
flow.waveAction=True
flow.Sinus_Amplitude = Max_Amplitude*sin(frequency*2*pi*(O.time-time_zero))

As suggested by Robert the change of boundary pressure is not
immediately reflected. However there is no need to rebuild the mesh and
the matrix of the problem every time (since changing the boundary
condition boils down to change B in "A*X=B" while keeping matrix A
unchanged). You can call flow.updateBCs(), which is _much_ faster than
updateTriangulation.

Bruno

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.