Skip to content
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.

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

Fix .prm file in step-60 #14940

Merged
merged 1 commit into from Mar 25, 2023
Merged

Conversation

fdrmrc
Copy link
Contributor

@fdrmrc fdrmrc commented Mar 22, 2023

With 4 refinement cycles near the interface the constraint on the ratio between the diameter of the two grids is not fulfilled in the case of a flower-like interface.

@luca-heltai I could push here the routine we are using to adjust that automatically. What do you think?

@marcfehling
Copy link
Member

Thanks for the patch! Can you also update the corresponding tutorial documentation?

<h3> Test case 2 and 3: </h3>
If we use the following parameter file:
@code
subsection Distributed Lagrange<1,2>
set Coupling quadrature order = 3
set Embedded configuration finite element degree = 1
set Embedded space finite element degree = 1
set Embedding space finite element degree = 1
set Homogeneous Dirichlet boundary ids = 0,1,2,3
set Initial embedded space refinement = 8
set Initial embedding space refinement = 4
set Local refinements steps near embedded domain = 4
set Use displacement in embedded interface = false
set Verbosity level = 10
subsection Embedded configuration
set Function constants = R=.3, Cx=.5, Cy=.5, r=.1, w=12
set Function expression = (R+r*cos(w*pi*x))*cos(2*pi*x)+Cx; (R+r*cos(w*pi*x))*sin(2*pi*x)+Cy
set Variable names = x,y,t
end
subsection Embedded value
set Function constants =
set Function expression = x-.5
set Variable names = x,y,t
end
subsection Embedding Dirichlet boundary conditions
set Function constants =
set Function expression = 0
set Variable names = x,y,t
end
subsection Embedding rhs function
set Function constants =
set Function expression = 0.0
set Variable names = x,y,t
end
subsection Schur solver control
set Log frequency = 1
set Log history = false
set Log result = true
set Max steps = 1000
set Reduction = 1.e-12
set Tolerance = 1.e-12
end
end
@endcode

@fdrmrc
Copy link
Contributor Author

fdrmrc commented Mar 22, 2023

I updated the relevant lines in the the .dox file in this patch.

@marcfehling
Copy link
Member

I updated the relevant lines in the the .dox file in this patch.

Oh, my bad! For some reason I thought you updated a stand-alone *.prm file...

@luca-heltai
Copy link
Member

Marco, yes, I think we should update the code with the routine we had developed to fix the issue.

@drwells drwells merged commit 31826f3 into dealii:master Mar 25, 2023
11 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

5 participants