Hi Robin, thank you. yes, it works fine without restart the *.map files. On Tue, Mar 10, 2020, 2:47 AM Basso, Robin G <r.basso17@imperial.ac.uk> wrote:
Hi Wang,
I would suggest to try the restart first in 2D, from a 2D simulation. Once you are confident with that you can try with your homogeneous direction. Does your session file running fine without restart from file?
Best, Robin.
On 9 Mar 2020, at 20:37, Wang, Zhicheng <zhicheng_wang@brown.edu> wrote:
Hi Robin, Thank you very much for your kind reply. I have tried the approach you mentioned, but it complained that something was wrong for the 3DH1D simulation, please see the errors below,
************************************begin error***************************************************************** Field z not found. Field z not found. Field z not found. Field z not found. Fatal : Level 0 assertion violation 3DH1D does not support mapping in the z-direction. Fatal : Level 0 assertion violation 3DH1D does not support mapping in the z-direction. ****************************************end error ****************************************************************
*Also, I copy all the condition parameters here, The case I'm running is 2D VIV of a rigid cylinder. Maybe I have misused some parameters ?*
<CONDITIONS> <SOLVERINFO> <I PROPERTY="SolverType" VALUE="VCSMapping" /> <I PROPERTY="EQTYPE" VALUE="UnsteadyNavierStokes" /> <I PROPERTY="AdvectionForm" VALUE="Convective" /> <I PROPERTY="Projection" VALUE="Galerkin" /> <I PROPERTY="HOMOGENEOUS" VALUE="1D"/> <I PROPERTY="USEFFT" VALUE="FFTW" /> <I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder2" /> <I PROPERTY="EntropyViscosityMethod" VALUE="True" /> <I PROPERTY="GlobalSysSoln" VALUE="IterativeStaticCond" /> <I PROPERTY="Preconditioner" VALUE="Block"/> <I PROPERTY="DynamicEVM" VALUE="True" /> <I PROPERTY="VibrationType" VALUE="Free"/> <I PROPERTY="SupportType" VALUE="Free-Free"/> <I PROPERTY="RigidVibration" VALUE="True"/> <I PROPERTY="FictitiousMassMethod" VALUE="True"/> </SOLVERINFO>
<PARAMETERS> <P> TimeStep = 0.0025 </P> <P> FinalTime = 1000 </P> <P> NumSteps = FinalTime/TimeStep </P> <P> IO_CheckSteps = 500 </P> <P> IO_InfoSteps = 25 </P> <P> IO_CFLSteps = 25 </P> <P> Re = 500 </P> <P> Kinvis = 1/Re </P> <P> HomModesZ = 2 </P> <P> LZ = 0.01 </P> <P> evm_alpha = 0.1 </P> <P> evm_beta = 0.1 </P> <P> evm_start_step = 2 </P> <P> evm_reduced_order = 1 </P> <P> StructFreqY = 0.17 </P> <P> StructFreqX = 0.34 </P> <P> StructStiff = 0.0 </P> <P> StructRho = 1.5708 </P> <P> CableTension = 0.0 </P> <P> BendingStiff = 0.0 </P> <P> FictDamp = 0.0 </P> <P> FictMass = 1.5 </P> </PARAMETERS>
<VARIABLES> <V ID="0"> u </V> <V ID="1"> v </V> <V ID="2"> w </V> <V ID="3"> p </V> </VARIABLES>
<BOUNDARYREGIONS> <B ID="0"> C[5] </B> <!-- Wall --> <B ID="1"> C[3] </B> <!-- y_min --> <B ID="2"> C[4] </B> <!-- y_max --> <B ID="3"> C[2] </B> <!-- Outlet --> <B ID="4"> C[1] </B> <!-- Inflow --> </BOUNDARYREGIONS>
<BOUNDARYCONDITIONS> <REGION REF="0"> <D VAR="u" USERDEFINEDTYPE="MovingBody" VALUE="0" /> <D VAR="v" USERDEFINEDTYPE="MovingBody" VALUE="0" /> <D VAR="w" VALUE="0" /> <N VAR="p" USERDEFINEDTYPE="H" VALUE="0" /> </REGION> <REGION REF="1"> <P VAR="u" VALUE="[2]" /> <P VAR="v" VALUE="[2]" /> <P VAR="w" VALUE="[2]" /> <P VAR="p" VALUE="[2]" /> </REGION> <REGION REF="2"> <P VAR="u" VALUE="[1]" /> <P VAR="v" VALUE="[1]" /> <P VAR="w" VALUE="[1]" /> <P VAR="p" VALUE="[1]" /> </REGION> <REGION REF="3"> <N VAR="u" VALUE="0" /> <N VAR="v" VALUE="0" /> <N VAR="w" VALUE="0" /> <D VAR="p" VALUE="0" /> </REGION> <REGION REF="4"> <D VAR="u" VALUE="1" /> <D VAR="v" VALUE="0" /> <D VAR="w" VALUE="0" /> <N VAR="p" USERDEFINEDTYPE="H" VALUE="0" /> </REGION> </BOUNDARYCONDITIONS>
<FUNCTION NAME="InitialConditions"> <F FILE="cyl.rst" /> </FUNCTION>
<FUNCTION NAME="ExactSolution"> <E VAR="u" VALUE="1" /> <E VAR="v" VALUE="0" /> <E VAR="w" VALUE="0" /> <E VAR="p" VALUE="0" /> </FUNCTION>
<FUNCTION NAME="Mapping"> <F FILE="cyl.map.rst" /> </FUNCTION>
<FUNCTION NAME="MappingVel"> <F FILE="cyl.map.rst" /> </FUNCTION>
<FUNCTION NAME="displacements"> <E VAR="x" VALUE="0.0" /> <E VAR="y" VALUE="A*cos(omega*t)*cos(2*PI*z/LZ)" /> </FUNCTION> <FUNCTION NAME="velocities"> <E VAR="x" VALUE="0.0" /> <E VAR="y" VALUE="-1.0*omega*A*sin(omega*t)*cos(2*PI*z/LZ)" /> </FUNCTION> <FUNCTION NAME="accelerations"> <E VAR="x" VALUE="0.0" /> <E VAR="y" VALUE="-1.0*omega*omega*A*cos(omega*t)*cos(2*PI*z/LZ)" /> </FUNCTION>
<FUNCTION NAME="MappingFile"> <F FILE="cyl.map.rst" /> </FUNCTION>
</CONDITIONS>
<MAPPING TYPE="XYofZ"> <COORDS>Mapping</COORDS> <VEL>MappingVel</VEL> <TIMEDEPENDENT>True</TIMEDEPENDENT> </MAPPING>
<FORCING> <FORCE TYPE="MovingBody"> <DISPLACEMENTS> displacements </DISPLACEMENTS> <VELOCITIES> velocities </VELOCITIES> <ACCELERATIONS> accelerations </ACCELERATIONS> <PARAM NAME="OutputFile">TensionedCable</PARAM> <PARAM NAME="OutputFrequency">25</PARAM> <PARAM NAME="Boundary"> B[0] </PARAM> </FORCE> </FORCING>
<FILTER TYPE="HistoryPoints"> <PARAM NAME="OutputFile">TimeValues</PARAM> <PARAM NAME="OutputFrequency">10</PARAM> <PARAM NAME="Points"> 3.0 0.0 0.0 5.0 0.0 0.0 7.0 0.0 0.0 </PARAM> </FILTER> </FILTERS>
Best, Zhicheng
On Mon, Mar 9, 2020 at 4:18 PM Basso, Robin G <r.basso17@imperial.ac.uk> wrote:
Hi,
I attached an example of a simulation of two cylinders restarted from file.
You need to create two functions, namely mapping and mappingVel for instance, where you call your .map file that contains the mesh deformation and velocity. Then you need to specify in your mapping tag that you restart your simulation from this file, calling these two function (defining coords and velocity).
I hope this helps, tell me if you can’t make it work.
Best, Robin.
<?xml version="1.0" encoding="utf-8" ?> <NEKTAR> <EXPANSIONS> <E COMPOSITE="C[119]" NUMMODES="5" TYPE="MODIFIED" FIELDS="u,v,p" /> </EXPANSIONS> <CONDITIONS> <SOLVERINFO> <I PROPERTY="SolverType" VALUE="VCSMapping" /> <I PROPERTY="EQTYPE" VALUE= "UnsteadyNavierStokes" /> <I PROPERTY="EvolutionOperator" VALUE="Nonlinear" /> <I PROPERTY="Projection" VALUE="Galerkin" /> <I PROPERTY="GlobalSysSoln" VALUE= "XxtMultiLevelStaticCond"/> <I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder1" /> <I PROPERTY="MappingImplicitPressure" VALUE="True" /> <I PROPERTY="MappingImplicitViscous" VALUE="True" /> </SOLVERINFO>
<PARAMETERS> <P> TimeStep = 0.004 </P> <P> T = 40 </P> <P> NumSteps = T/TimeStep </P> <P> IO_CheckSteps = NumSteps/200 </P> <P> IO_InfoSteps = 10 </P> <P> IO_CFLSteps = 10 </P> <P> Re = 150 </P> <P> Kinvis = 1/Re </P>
<P> forcedAmplitudeY = 0.5 </P> <P> forcedAmplitudeX = 0.2*forcedAmplitudeY </P> <P> forcedFreqY = 1.5*0.2 </P> <P> forcedFreqX = 2*forcedFreqY </P>
<P> MappingPressureTolerance = 1e-6 </P> <P> MappingViscousTolerance = 1e-8 </P> <P> MappingPressureRelaxation = 0.4 </P> <P> MappingViscousRelaxation = 0.4 </P>
<!-- Structural parameters --> <P> Mr = 4.0 </P> <P> Cr = 0.007 </P> <P> Vr = 10.0 </P>
<P> fn = 1.0/Vr </P> <P> M = Mr * PI / 4.0 </P> <P> K = M* ( 2*PI *fn ) * ( 2*PI *fn ) </P> <P> C = Cr * 2 * sqrt(K*M) </P> </PARAMETERS>
<VARIABLES> <V ID="0"> u </V> <V ID="1"> v </V> <V ID="2"> p </V> </VARIABLES>
<BOUNDARYREGIONS> <B ID="0"> C[125-126] </B> <!-- Downstream cylinder --> <B ID="1"> C[120-121] </B> <!-- Far Field --> <B ID="2"> C[122] </B> <!-- Outlet --> <B ID="3"> C[123-124] </B> <!-- Upstream cylinder --> </BOUNDARYREGIONS>
<BOUNDARYCONDITIONS> <REGION REF="0"> <D VAR="u" USERDEFINEDTYPE="MovingBody" VALUE="0" /> <D VAR="v" USERDEFINEDTYPE="MovingBody" VALUE="0" /> <N VAR="p" USERDEFINEDTYPE="H" VALUE="0" /> <D VAR="x" VALUE="0" /> <D VAR="y" VALUE="0" /> </REGION> <REGION REF="1"> <D VAR="u" VALUE="1" /> <D VAR="v" VALUE="0" /> <N VAR="p" USERDEFINEDTYPE="H" VALUE="0" /> <D VAR="x" VALUE="0" /> <D VAR="y" VALUE="0" /> </REGION> <REGION REF="2"> <N VAR="u" VALUE="0" /> <N VAR="v" VALUE="0" /> <D VAR="p" VALUE="0" /> <D VAR="x" VALUE="0" /> <D VAR="y" VALUE="0" /> </REGION> <REGION REF="3"> <D VAR="u" USERDEFINEDTYPE="MovingBody" VALUE="0" /> <D VAR="v" USERDEFINEDTYPE="MovingBody" VALUE="0" /> <N VAR="p" USERDEFINEDTYPE="H" VALUE="0" /> <D VAR="x" VALUE="0" /> <D VAR="y" VALUE="0" /> </REGION> </BOUNDARYCONDITIONS>
<FUNCTION NAME="Displacement"> <E VAR="x" VALUE="forcedAmplitudeX*sin(2*PI*forcedFreqX*t)" /> <E VAR="y" VALUE="forcedAmplitudeY*sin(2*PI*forcedFreqY*t)" /> </FUNCTION>
<FUNCTION NAME="InitialConditions"> <F FILE="2cyl.fld" /> </FUNCTION>
<FUNCTION NAME="Mapping"> <F FILE="2cyl.map" /> </FUNCTION>
<FUNCTION NAME="MappingVel"> <F FILE="2cyl.map" /> </FUNCTION>
</CONDITIONS>
<MAPPING TYPE="XYofXY"> <COORDS> Mapping </COORDS> <VEL> MappingVel </VEL> </MAPPING>
<FSI TYPE="General" > <BODY TYPE="Rigid" > <PARAM NAME="Boundary"> B[0] </PARAM> <PARAM NAME="OutputFile"> downCyl </PARAM> <PARAM NAME="OutputFrequency"> 10 </PARAM> <PARAM NAME="TranslationDOFs"> 1 </PARAM> <PARAM NAME="Direction1"> 0.0 1.0 0.0</PARAM> <PARAM NAME="StartTime"> 0.0 </PARAM> <PARAM NAME="M"> M </PARAM> <PARAM NAME="C"> C </PARAM> <PARAM NAME="K"> K </PARAM> </BODY> <BODY TYPE="Forced" > <PARAM NAME="Boundary"> B[3] </PARAM> <PARAM NAME="DisplacementFCN"> Displacement </PARAM> </BODY> </FSI>
<FILTERS> <FILTER TYPE="AeroForces"> <PARAM NAME="OutputFile"> forcesDownCyl </PARAM> <PARAM NAME="OutputFrequency"> 10 </PARAM> <PARAM NAME="Boundary"> B[0] </PARAM> </FILTER> <FILTER TYPE="AeroForces"> <PARAM NAME="OutputFile"> forcesUpCyl </PARAM> <PARAM NAME="OutputFrequency"> 10 </PARAM> <PARAM NAME="Boundary"> B[3] </PARAM> </FILTER> </FILTERS>
<GEOMETRY DIM="2" SPACE="2"> <VERTEX COMPRESSED="B64Z-LittleEndian" BITSIZE="64"> eJxtm3dgT9cbxr8SRCLIED+1d+0Zm9yrlBqNaM2qUZuiiqhSI2IlSI1WSKVUiYpYIcS8MWpTu2YJsXcUIcavbd73vfc+99x/XsnHOd/zfe5z3vPec09yurKutg2DPqmUVChFy3Xu29S1vjr92pXNZbtSNp0OH2Pl2SnuzxwYGXvO18HdKUYWeuI5/qpHCv0oHC4HzwH9I3eD/vHzPeD7Yftcto8v4/j+nsCxvReMv+oknwvb2xYVnhvGh9wbvh/yPDD+zUl9zlt5Xhgf8nw
On 9 Mar 2020, at 19:41, Sherwin, Spencer J <s.sherwin@imperial.ac.uk> wrote:
Hi Robin,
Could you help Zhicheng with this request.
Thanks, Spencer.
On 5 Mar 2020, at 22:59, Wang, Zhicheng <zhicheng_wang@brown.edu> wrote:
This email from zhicheng_wang@brown.edu originates from outside Imperial. Do not click on links and attachments unless you recognise the sender. If you trust the sender, add them to your safe senders list <https://spam.ic.ac.uk/SpamConsole/Senders.aspx> to disable email stamping for this address.
Hi nektar++ users, Is there anyone can give me a hint how to restart a simulation using the VCSMapping solver ? I see folders named *.chk and *.map, but how to use the *.map file to restart a new simulation ?
Best, Zhicheng _______________________________________________ Nektar-users mailing list Nektar-users@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/nektar-users
Spencer Sherwin FREng, FRAeS Head of Aerodynamics Section, Director of Research Computing Service, Professor of Computational Fluid Mechanics, Department of Aeronautics, s.sherwin@imperial.ac.uk South Kensington Campus, Phone: +44 (0)20 7594 5052 Imperial College London, Fax: +44 (0)20 7594 1974 London, SW7 2AZ, UK http://www.imperial.ac.uk/people/s.sherwin/