Skip to content

Bug: obstacle does not work together with Pressure solver #209

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Draft
wants to merge 11 commits into
base: master
Choose a base branch
from
Draft

Bug: obstacle does not work together with Pressure solver #209

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
11 commits
Select commit Hold shift + click to select a range
b674681
start
LinhW Mar 24, 2022
63758d5
working xml
LinhW Mar 24, 2022
6856d4e
revert changes
LinhW Mar 28, 2022
aad4756
added test for even and odd size of obstacle
LinhW Mar 28, 2022
8cfdbed
added testcases with 64x64x64 grid
LinhW Mar 28, 2022
9eec106
sladvect gtest
LinhWuerzburger Mar 28, 2022
fa153a7
Merge branch 'WIP' into bug_obstacle
LinhWuerzburger Mar 28, 2022
6c7751d
sladvect gtests
LinhWuerzburger Mar 28, 2022
d0044da
fix of obstacle test
LinhWuerzburger Mar 28, 2022
6f750c9
fix sl advect unit test
LinhWuerzburger Mar 29, 2022
1fe34dd
Merge branch 'master' into bug_obstacle
LinhWuerzburger Jul 28, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
85 changes: 85 additions & 0 deletions debug/bug_obstacle/Test_Flow_not_working.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,85 @@
<?xml version="1.0" encoding="UTF-8" ?>
<ARTSS>
<physical_parameters>
<t_end> 1. </t_end> <!-- simulation end time -->
<dt> 0.01 </dt> <!-- time stepping, caution: CFL-condition dt < 0.5*dx^2/nu -->
<nu> 0.00001 </nu> <!-- kinematic viscosity -->
</physical_parameters>

<solver description="NSTurbSolver" >
<advection type="SemiLagrangian" field="u,v,w">
</advection>
<diffusion type="Explicit" field="u,v,w">
<max_iter> 100 </max_iter> <!-- max number of iterations -->
<tol_res> 1e-07 </tol_res> <!-- tolerance for residuum/ convergence -->
<w> 1 </w> <!-- relaxation parameter -->
</diffusion>
<turbulence type="ConstSmagorinsky">
<Cs> 0.2 </Cs>
</turbulence>
<source type="ExplicitEuler" force_fct="Zero" dir="xyz"> <!--Direction of force (x,y,z or combinations xy,xz,yz, xyz) -->
</source>
<pressure type="VCycleMG" field="p">
<n_level> 1 </n_level> <!-- number of restriction levels -->
<n_cycle> 1 </n_cycle> <!-- number of cycles -->
<max_cycle> 4 </max_cycle> <!-- maximal number of cycles in first time step -->
<tol_res> 1e-07 </tol_res> <!-- tolerance for residuum/ convergence -->
<n_relax> 4 </n_relax> <!-- number of iterations -->
<diffusion type="Jacobi" field="p">
<max_iter> 100 </max_iter> <!-- maximal number of iterations in solving at lowest level -->
<tol_res> 1e-07 </tol_res> <!-- tolerance for residuum/ convergence -->
<w> 0.6666666667 </w> <!-- relaxation parameter -->
</diffusion>
</pressure>
<solution available="No">
</solution>
</solver>

<domain_parameters enable_computational_domain="No">
<X1> -3. </X1> <!-- physical domain -->
<X2> 7. </X2>
<Y1> 0. </Y1>
<Y2> 2. </Y2>
<Z1> -3.5 </Z1>
<Z2> 3.5 </Z2>
<nx> 256 </nx> <!-- grid resolution (number of cells excl. ghost cells) -->
<ny> 128 </ny>
<nz> 256 </nz>
</domain_parameters>

<adaption dynamic="No" data_extraction="No"> </adaption>

<data_assimilation enabled="No"/>

<boundaries>
<boundary field="u" patch="left,right" type="dirichlet" value="0.4" />
<boundary field="v,w" patch="left,right" type="dirichlet" value="0.0" />
<boundary field="u,v,w" patch="front,back,bottom,top" type="neumann" value="0.0" />
<boundary field="p" patch="front,back,bottom,top,left,right" type="neumann" value="0.0" />
</boundaries>

<obstacles enabled="Yes">
<obstacle name="cube">
<geometry ox1="0.0273" ox2="0.964" oy1="0.0078" oy2="0.992" oz1="-0.492" oz2="0.4785"/>
<!-- <geometry ox1="0.0273" ox2="0.965" oy1="0.0078" oy2="0.994" oz1="-0.492" oz2="0.4795"/> -->
<boundary field="u,v,w" patch="front,back,left,right,bottom,top" type="dirichlet" value="0.0" />
<boundary field="p" patch="front,back,left,right,bottom,top" type="neumann" value="0.0" />
</obstacle>
</obstacles>

<surfaces enabled="No"/>

<initial_conditions usr_fct="Drift" random="No">
<u_lin> 0.4 </u_lin> <!-- background velocity -->
<v_lin> 0.0 </v_lin> <!-- background velocity -->
<w_lin> 0.0 </w_lin> <!-- background velocity -->
<pa> 0. </pa> <!-- ambient pressure -->
</initial_conditions>

<visualisation save_vtk="Yes" save_csv="No">
<vtk_nth_plot> 1 </vtk_nth_plot>
</visualisation>

<logging file="output_bug_flow.log" level="info">
</logging>
</ARTSS>
84 changes: 84 additions & 0 deletions debug/bug_obstacle/Test_Flow_working.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,84 @@
<?xml version="1.0" encoding="UTF-8" ?>
<ARTSS>
<physical_parameters>
<t_end> 1. </t_end> <!-- simulation end time -->
<dt> 0.01 </dt> <!-- time stepping, caution: CFL-condition dt < 0.5*dx^2/nu -->
<nu> 0.00001 </nu> <!-- kinematic viscosity -->
</physical_parameters>

<solver description="NSTurbSolver" >
<advection type="SemiLagrangian" field="u,v,w">
</advection>
<diffusion type="Explicit" field="u,v,w">
<max_iter> 100 </max_iter> <!-- max number of iterations -->
<tol_res> 1e-07 </tol_res> <!-- tolerance for residuum/ convergence -->
<w> 1 </w> <!-- relaxation parameter -->
</diffusion>
<turbulence type="ConstSmagorinsky">
<Cs> 0.2 </Cs>
</turbulence>
<source type="ExplicitEuler" force_fct="Zero" dir="xyz"> <!--Direction of force (x,y,z or combinations xy,xz,yz, xyz) -->
</source>
<pressure type="VCycleMG" field="p">
<n_level> 1 </n_level> <!-- number of restriction levels -->
<n_cycle> 1 </n_cycle> <!-- number of cycles -->
<max_cycle> 4 </max_cycle> <!-- maximal number of cycles in first time step -->
<tol_res> 1e-07 </tol_res> <!-- tolerance for residuum/ convergence -->
<n_relax> 4 </n_relax> <!-- number of iterations -->
<diffusion type="Jacobi" field="p">
<max_iter> 100 </max_iter> <!-- maximal number of iterations in solving at lowest level -->
<tol_res> 1e-07 </tol_res> <!-- tolerance for residuum/ convergence -->
<w> 0.6666666667 </w> <!-- relaxation parameter -->
</diffusion>
</pressure>
<solution available="No">
</solution>
</solver>

<domain_parameters enable_computational_domain="No">
<X1> -3. </X1> <!-- physical domain -->
<X2> 7. </X2>
<Y1> 0. </Y1>
<Y2> 2. </Y2>
<Z1> -3.5 </Z1>
<Z2> 3.5 </Z2>
<nx> 256 </nx> <!-- grid resolution (number of cells excl. ghost cells) -->
<ny> 128 </ny>
<nz> 256 </nz>
</domain_parameters>

<adaption dynamic="No" data_extraction="No"> </adaption>

<data_assimilation enabled="No"/>

<boundaries>
<boundary field="u" patch="left,right" type="dirichlet" value="0.4" />
<boundary field="v,w" patch="left,right" type="dirichlet" value="0.0" />
<boundary field="u,v,w" patch="front,back,bottom,top" type="neumann" value="0.0" />
<boundary field="p" patch="front,back,bottom,top,left,right" type="neumann" value="0.0" />
</boundaries>

<obstacles enabled="Yes">
<obstacle name="cube">
<geometry ox1="0.0273" ox2="0.965" oy1="0.0078" oy2="0.994" oz1="-0.492" oz2="0.4795"/>
<boundary field="u,v,w" patch="front,back,left,right,bottom,top" type="dirichlet" value="0.0" />
<boundary field="p" patch="front,back,left,right,bottom,top" type="neumann" value="0.0" />
</obstacle>
</obstacles>

<surfaces enabled="No"/>

<initial_conditions usr_fct="Drift" random="No">
<u_lin> 0.4 </u_lin> <!-- background velocity -->
<v_lin> 0.0 </v_lin> <!-- background velocity -->
<w_lin> 0.0 </w_lin> <!-- background velocity -->
<pa> 0. </pa> <!-- ambient pressure -->
</initial_conditions>

<visualisation save_vtk="Yes" save_csv="No">
<vtk_nth_plot> 1 </vtk_nth_plot>
</visualisation>

<logging file="output_working_flow.log" level="info">
</logging>
</ARTSS>
84 changes: 84 additions & 0 deletions debug/bug_obstacle/floating/flow_simple.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,84 @@
<?xml version="1.0" encoding="UTF-8" ?>
<ARTSS>
<physical_parameters>
<t_end> 10. </t_end> <!-- simulation end time -->
<dt> 0.01 </dt> <!-- time stepping, caution: CFL-condition dt < 0.5*dx^2/nu -->
<nu> 0.00001 </nu> <!-- kinematic viscosity -->
</physical_parameters>

<solver description="NSTurbSolver" >
<advection type="SemiLagrangian" field="u,v,w">
</advection>
<diffusion type="Explicit" field="u,v,w">
<max_iter> 100 </max_iter> <!-- max number of iterations -->
<tol_res> 1e-07 </tol_res> <!-- tolerance for residuum/ convergence -->
<w> 1 </w> <!-- relaxation parameter -->
</diffusion>
<turbulence type="ConstSmagorinsky">
<Cs> 0.2 </Cs>
</turbulence>
<source type="ExplicitEuler" force_fct="Zero" dir="xyz"> <!--Direction of force (x,y,z or combinations xy,xz,yz, xyz) -->
</source>
<pressure type="VCycleMG" field="p">
<n_level> 2 </n_level> <!-- number of restriction levels -->
<n_cycle> 1 </n_cycle> <!-- number of cycles -->
<max_cycle> 4 </max_cycle> <!-- maximal number of cycles in first time step -->
<tol_res> 1e-07 </tol_res> <!-- tolerance for residuum/ convergence -->
<n_relax> 4 </n_relax> <!-- number of iterations -->
<diffusion type="Jacobi" field="p">
<max_iter> 100 </max_iter> <!-- maximal number of iterations in solving at lowest level -->
<tol_res> 1e-07 </tol_res> <!-- tolerance for residuum/ convergence -->
<w> 0.6666666667 </w> <!-- relaxation parameter -->
</diffusion>
</pressure>
<solution available="No">
</solution>
</solver>

<domain_parameters enable_computational_domain="No">
<X1> -4 </X1> <!-- physical domain -->
<X2> 4 </X2>
<Y1> -4 </Y1>
<Y2> 4 </Y2>
<Z1> -4 </Z1>
<Z2> 4 </Z2>
<nx> 64 </nx> <!-- grid resolution (number of cells excl. ghost cells) -->
<ny> 64 </ny>
<nz> 64 </nz>
</domain_parameters>

<adaption dynamic="No" data_extraction="No"> </adaption>

<data_assimilation enabled="No"/>

<boundaries>
<boundary field="u" patch="left,right" type="dirichlet" value="0.4" />
<boundary field="v,w" patch="left,right" type="dirichlet" value="0.0" />
<boundary field="u,v,w" patch="front,back,bottom,top" type="neumann" value="0.0" />
<boundary field="p" patch="front,back,bottom,top,left,right" type="neumann" value="0.0" />
</boundaries>

<obstacles enabled="Yes">
<obstacle name="cube">
<geometry ox1="-1" ox2="1" oy1="-1" oy2="1" oz1="-1" oz2="1"/>
<boundary field="u,v,w" patch="front,back,left,right,bottom,top" type="dirichlet" value="0.0" />
<boundary field="p" patch="front,back,left,right,bottom,top" type="neumann" value="0.0" />
</obstacle>
</obstacles>

<surfaces enabled="No"/>

<initial_conditions usr_fct="Drift" random="No">
<u_lin> 0.4 </u_lin> <!-- background velocity -->
<v_lin> 0.0 </v_lin> <!-- background velocity -->
<w_lin> 0.0 </w_lin> <!-- background velocity -->
<pa> 0. </pa> <!-- ambient pressure -->
</initial_conditions>

<visualisation save_vtk="Yes" save_csv="No">
<vtk_nth_plot> 1 </vtk_nth_plot>
</visualisation>

<logging file="output_flow_simple.log" level="info">
</logging>
</ARTSS>
84 changes: 84 additions & 0 deletions debug/bug_obstacle/floating/flow_simple2.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,84 @@
<?xml version="1.0" encoding="UTF-8" ?>
<ARTSS>
<physical_parameters>
<t_end> 10. </t_end> <!-- simulation end time -->
<dt> 0.01 </dt> <!-- time stepping, caution: CFL-condition dt < 0.5*dx^2/nu -->
<nu> 0.00001 </nu> <!-- kinematic viscosity -->
</physical_parameters>

<solver description="NSTurbSolver" >
<advection type="SemiLagrangian" field="u,v,w">
</advection>
<diffusion type="Explicit" field="u,v,w">
<max_iter> 100 </max_iter> <!-- max number of iterations -->
<tol_res> 1e-07 </tol_res> <!-- tolerance for residuum/ convergence -->
<w> 1 </w> <!-- relaxation parameter -->
</diffusion>
<turbulence type="ConstSmagorinsky">
<Cs> 0.2 </Cs>
</turbulence>
<source type="ExplicitEuler" force_fct="Zero" dir="xyz"> <!--Direction of force (x,y,z or combinations xy,xz,yz, xyz) -->
</source>
<pressure type="VCycleMG" field="p">
<n_level> 2 </n_level> <!-- number of restriction levels -->
<n_cycle> 1 </n_cycle> <!-- number of cycles -->
<max_cycle> 4 </max_cycle> <!-- maximal number of cycles in first time step -->
<tol_res> 1e-07 </tol_res> <!-- tolerance for residuum/ convergence -->
<n_relax> 4 </n_relax> <!-- number of iterations -->
<diffusion type="Jacobi" field="p">
<max_iter> 100 </max_iter> <!-- maximal number of iterations in solving at lowest level -->
<tol_res> 1e-07 </tol_res> <!-- tolerance for residuum/ convergence -->
<w> 0.6666666667 </w> <!-- relaxation parameter -->
</diffusion>
</pressure>
<solution available="No">
</solution>
</solver>

<domain_parameters enable_computational_domain="No">
<X1> -4 </X1> <!-- physical domain -->
<X2> 4 </X2>
<Y1> -4 </Y1>
<Y2> 4 </Y2>
<Z1> -4 </Z1>
<Z2> 4 </Z2>
<nx> 64 </nx> <!-- grid resolution (number of cells excl. ghost cells) -->
<ny> 64 </ny>
<nz> 64 </nz>
</domain_parameters>

<adaption dynamic="No" data_extraction="No"> </adaption>

<data_assimilation enabled="No"/>

<boundaries>
<boundary field="u" patch="left,right" type="dirichlet" value="0.4" />
<boundary field="v,w" patch="left,right" type="dirichlet" value="0.0" />
<boundary field="u,v,w" patch="front,back,bottom,top" type="neumann" value="0.0" />
<boundary field="p" patch="front,back,bottom,top,left,right" type="neumann" value="0.0" />
</boundaries>

<obstacles enabled="Yes">
<obstacle name="cube">
<geometry ox1="-1" ox2="1.125" oy1="-1" oy2="1" oz1="-1" oz2="1"/>
<boundary field="u,v,w" patch="front,back,left,right,bottom,top" type="dirichlet" value="0.0" />
<boundary field="p" patch="front,back,left,right,bottom,top" type="neumann" value="0.0" />
</obstacle>
</obstacles>

<surfaces enabled="No"/>

<initial_conditions usr_fct="Drift" random="No">
<u_lin> 0.4 </u_lin> <!-- background velocity -->
<v_lin> 0.0 </v_lin> <!-- background velocity -->
<w_lin> 0.0 </w_lin> <!-- background velocity -->
<pa> 0. </pa> <!-- ambient pressure -->
</initial_conditions>

<visualisation save_vtk="Yes" save_csv="No">
<vtk_nth_plot> 1 </vtk_nth_plot>
</visualisation>

<logging file="output_flow_simple2.log" level="info">
</logging>
</ARTSS>
Loading