From a543aa7895c096f607afd0181bb57321a0329eae Mon Sep 17 00:00:00 2001 From: Stephen DeWitt Date: Wed, 23 Aug 2017 16:41:53 -0400 Subject: [PATCH] improved the non-uniform dirichlet BC description in ICs_and_BCs.h --- applications/CHAC_anisotropy/ICs_and_BCs.h | 16 ++++++++--- .../CHAC_anisotropyRegularized/ICs_and_BCs.h | 15 +++++++--- applications/CHiMaD_benchmark1a/ICs_and_BCs.h | 13 +++++++-- applications/CHiMaD_benchmark2a/ICs_and_BCs.h | 14 ++++++++-- applications/CHiMaD_benchmark6a/ICs_and_BCs.h | 28 ++++++++++++------- applications/CHiMaD_benchmark6b/ICs_and_BCs.h | 28 ++++++++++++------- applications/allenCahn/ICs_and_BCs.h | 13 +++++++-- applications/allenCahn_pfield/ICs_and_BCs.h | 13 +++++++-- applications/anisotropyFacet/ICs_and_BCs.h | 17 +++++++---- applications/cahnHilliard/ICs_and_BCs.h | 12 ++++++-- .../cahnHilliardWithAdaptivity/ICs_and_BCs.h | 13 +++++++-- .../ICs_and_BCs.h | 16 ++++++++--- .../dendriticSolidification/ICs_and_BCs.h | 16 ++++++++--- applications/eshelbyInclusion/ICs_and_BCs.h | 12 ++++++-- applications/fickianDiffusion/ICs_and_BCs.h | 12 ++++++-- applications/grainGrowth/ICs_and_BCs.h | 12 ++++++-- applications/mechanics/ICs_and_BCs.h | 13 +++++++-- applications/nucleationModel/ICs_and_BCs.h | 17 +++++++---- .../precipitateEvolution/ICs_and_BCs.h | 12 ++++++-- .../ICs_and_BCs.h | 12 ++++++-- .../ICs_and_BCs.h | 16 ++++++++--- .../singlePrecipitateKKS/ICs_and_BCs.h | 12 ++++++-- 22 files changed, 250 insertions(+), 82 deletions(-) diff --git a/applications/CHAC_anisotropy/ICs_and_BCs.h b/applications/CHAC_anisotropy/ICs_and_BCs.h index 457f1b18f..2542085f7 100644 --- a/applications/CHAC_anisotropy/ICs_and_BCs.h +++ b/applications/CHAC_anisotropy/ICs_and_BCs.h @@ -63,10 +63,14 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // -------------------------------------------------------------------------- // ENTER THE NON-UNIFORM DIRICHLET BOUNDARY CONDITIONS HERE FOR SCALAR FIELDS // -------------------------------------------------------------------------- - // Enter the function describing conditions for the fields at point "p". + // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- @@ -80,10 +84,14 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // -------------------------------------------------------------------------- // ENTER THE NON-UNIFORM DIRICHLET BOUNDARY CONDITIONS HERE FOR VECTOR FIELDS // -------------------------------------------------------------------------- - // Enter the function describing conditions for the fields at point "p". + // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/CHAC_anisotropyRegularized/ICs_and_BCs.h b/applications/CHAC_anisotropyRegularized/ICs_and_BCs.h index c16b901e4..314c45e72 100644 --- a/applications/CHAC_anisotropyRegularized/ICs_and_BCs.h +++ b/applications/CHAC_anisotropyRegularized/ICs_and_BCs.h @@ -67,11 +67,14 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // -------------------------------------------------------------------------- // ENTER THE NON-UNIFORM DIRICHLET BOUNDARY CONDITIONS HERE FOR SCALAR FIELDS // -------------------------------------------------------------------------- - // Enter the function describing conditions for the fields at point "p". + // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. - + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- return scalar_BC; @@ -87,7 +90,11 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/CHiMaD_benchmark1a/ICs_and_BCs.h b/applications/CHiMaD_benchmark1a/ICs_and_BCs.h index 8f119a13c..6dd82abec 100644 --- a/applications/CHiMaD_benchmark1a/ICs_and_BCs.h +++ b/applications/CHiMaD_benchmark1a/ICs_and_BCs.h @@ -58,7 +58,11 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- @@ -75,8 +79,11 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. - + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/CHiMaD_benchmark2a/ICs_and_BCs.h b/applications/CHiMaD_benchmark2a/ICs_and_BCs.h index bd6d47564..295c5773f 100644 --- a/applications/CHiMaD_benchmark2a/ICs_and_BCs.h +++ b/applications/CHiMaD_benchmark2a/ICs_and_BCs.h @@ -90,10 +90,14 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // -------------------------------------------------------------------------- // ENTER THE NON-UNIFORM DIRICHLET BOUNDARY CONDITIONS HERE FOR SCALAR FIELDS // -------------------------------------------------------------------------- - // Enter the function describing conditions for the fields at point "p". + // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- @@ -110,7 +114,11 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/CHiMaD_benchmark6a/ICs_and_BCs.h b/applications/CHiMaD_benchmark6a/ICs_and_BCs.h index f0192c74b..dba266584 100644 --- a/applications/CHiMaD_benchmark6a/ICs_and_BCs.h +++ b/applications/CHiMaD_benchmark6a/ICs_and_BCs.h @@ -18,11 +18,11 @@ double InitialCondition::value (const dealii::Point &p, const unsigned double y=p[1]; double c0=0.5; double c1=0.04; - + double t1=std::cos(0.2*x)*std::cos(0.11*y); double t2=std::cos(0.13*x)*std::cos(0.087*y)*std::cos(0.13*x)*std::cos(0.087*y); double t3=std::cos(0.025*x-0.15*y)*std::cos(0.07*x-0.02*y); - + scalar_IC = c0 + c1*(t1+t2+t3); } else { @@ -62,8 +62,12 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. - + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). + if (index == 2){ if (direction == 1){ double x=p[0]; @@ -78,16 +82,20 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns template void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, dealii::Vector &vector_BC) const { - + // -------------------------------------------------------------------------- // ENTER THE NON-UNIFORM DIRICHLET BOUNDARY CONDITIONS HERE FOR VECTOR FIELDS // -------------------------------------------------------------------------- // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. - - + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). + + // ------------------------------------------------------------------------- - -} \ No newline at end of file + +} diff --git a/applications/CHiMaD_benchmark6b/ICs_and_BCs.h b/applications/CHiMaD_benchmark6b/ICs_and_BCs.h index f0192c74b..dba266584 100644 --- a/applications/CHiMaD_benchmark6b/ICs_and_BCs.h +++ b/applications/CHiMaD_benchmark6b/ICs_and_BCs.h @@ -18,11 +18,11 @@ double InitialCondition::value (const dealii::Point &p, const unsigned double y=p[1]; double c0=0.5; double c1=0.04; - + double t1=std::cos(0.2*x)*std::cos(0.11*y); double t2=std::cos(0.13*x)*std::cos(0.087*y)*std::cos(0.13*x)*std::cos(0.087*y); double t3=std::cos(0.025*x-0.15*y)*std::cos(0.07*x-0.02*y); - + scalar_IC = c0 + c1*(t1+t2+t3); } else { @@ -62,8 +62,12 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. - + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). + if (index == 2){ if (direction == 1){ double x=p[0]; @@ -78,16 +82,20 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns template void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, dealii::Vector &vector_BC) const { - + // -------------------------------------------------------------------------- // ENTER THE NON-UNIFORM DIRICHLET BOUNDARY CONDITIONS HERE FOR VECTOR FIELDS // -------------------------------------------------------------------------- // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. - - + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). + + // ------------------------------------------------------------------------- - -} \ No newline at end of file + +} diff --git a/applications/allenCahn/ICs_and_BCs.h b/applications/allenCahn/ICs_and_BCs.h index 0e5f139d5..24c7ba6a4 100644 --- a/applications/allenCahn/ICs_and_BCs.h +++ b/applications/allenCahn/ICs_and_BCs.h @@ -64,7 +64,11 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- @@ -81,8 +85,11 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. - + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/allenCahn_pfield/ICs_and_BCs.h b/applications/allenCahn_pfield/ICs_and_BCs.h index 0e5f139d5..d8a6c589d 100644 --- a/applications/allenCahn_pfield/ICs_and_BCs.h +++ b/applications/allenCahn_pfield/ICs_and_BCs.h @@ -64,8 +64,11 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. - + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- return scalar_BC; @@ -81,7 +84,11 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/anisotropyFacet/ICs_and_BCs.h b/applications/anisotropyFacet/ICs_and_BCs.h index c16b901e4..2c23d57c1 100644 --- a/applications/anisotropyFacet/ICs_and_BCs.h +++ b/applications/anisotropyFacet/ICs_and_BCs.h @@ -67,10 +67,14 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // -------------------------------------------------------------------------- // ENTER THE NON-UNIFORM DIRICHLET BOUNDARY CONDITIONS HERE FOR SCALAR FIELDS // -------------------------------------------------------------------------- - // Enter the function describing conditions for the fields at point "p". + // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- @@ -84,11 +88,14 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // -------------------------------------------------------------------------- // ENTER THE NON-UNIFORM DIRICHLET BOUNDARY CONDITIONS HERE FOR VECTOR FIELDS // -------------------------------------------------------------------------- - // Enter the function describing conditions for the fields at point "p". + // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. - + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/cahnHilliard/ICs_and_BCs.h b/applications/cahnHilliard/ICs_and_BCs.h index a5168a60d..dd9d7997d 100644 --- a/applications/cahnHilliard/ICs_and_BCs.h +++ b/applications/cahnHilliard/ICs_and_BCs.h @@ -71,7 +71,11 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- @@ -88,7 +92,11 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/cahnHilliardWithAdaptivity/ICs_and_BCs.h b/applications/cahnHilliardWithAdaptivity/ICs_and_BCs.h index 04b3c704a..c8cf36af0 100644 --- a/applications/cahnHilliardWithAdaptivity/ICs_and_BCs.h +++ b/applications/cahnHilliardWithAdaptivity/ICs_and_BCs.h @@ -70,8 +70,11 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. - + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- return scalar_BC; @@ -87,7 +90,11 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/coupledCahnHilliardAllenCahn/ICs_and_BCs.h b/applications/coupledCahnHilliardAllenCahn/ICs_and_BCs.h index 2f2818986..df066d3cc 100644 --- a/applications/coupledCahnHilliardAllenCahn/ICs_and_BCs.h +++ b/applications/coupledCahnHilliardAllenCahn/ICs_and_BCs.h @@ -71,10 +71,14 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // -------------------------------------------------------------------------- // ENTER THE NON-UNIFORM DIRICHLET BOUNDARY CONDITIONS HERE FOR SCALAR FIELDS // -------------------------------------------------------------------------- - // Enter the function describing conditions for the fields at point "p". + // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- @@ -88,10 +92,14 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // -------------------------------------------------------------------------- // ENTER THE NON-UNIFORM DIRICHLET BOUNDARY CONDITIONS HERE FOR VECTOR FIELDS // -------------------------------------------------------------------------- - // Enter the function describing conditions for the fields at point "p". + // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/dendriticSolidification/ICs_and_BCs.h b/applications/dendriticSolidification/ICs_and_BCs.h index 6ac127832..7d4badf43 100644 --- a/applications/dendriticSolidification/ICs_and_BCs.h +++ b/applications/dendriticSolidification/ICs_and_BCs.h @@ -70,10 +70,14 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // -------------------------------------------------------------------------- // ENTER THE NON-UNIFORM DIRICHLET BOUNDARY CONDITIONS HERE FOR SCALAR FIELDS // -------------------------------------------------------------------------- - // Enter the function describing conditions for the fields at point "p". + // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- @@ -87,10 +91,14 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // -------------------------------------------------------------------------- // ENTER THE NON-UNIFORM DIRICHLET BOUNDARY CONDITIONS HERE FOR VECTOR FIELDS // -------------------------------------------------------------------------- - // Enter the function describing conditions for the fields at point "p". + // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/eshelbyInclusion/ICs_and_BCs.h b/applications/eshelbyInclusion/ICs_and_BCs.h index ae6255357..2cefdf8f5 100644 --- a/applications/eshelbyInclusion/ICs_and_BCs.h +++ b/applications/eshelbyInclusion/ICs_and_BCs.h @@ -47,7 +47,11 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- @@ -64,7 +68,11 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/fickianDiffusion/ICs_and_BCs.h b/applications/fickianDiffusion/ICs_and_BCs.h index a054e971d..b115d3637 100644 --- a/applications/fickianDiffusion/ICs_and_BCs.h +++ b/applications/fickianDiffusion/ICs_and_BCs.h @@ -47,7 +47,11 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- @@ -64,7 +68,11 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/grainGrowth/ICs_and_BCs.h b/applications/grainGrowth/ICs_and_BCs.h index 17998a5ad..e8866459e 100644 --- a/applications/grainGrowth/ICs_and_BCs.h +++ b/applications/grainGrowth/ICs_and_BCs.h @@ -61,7 +61,11 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- @@ -78,7 +82,11 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/mechanics/ICs_and_BCs.h b/applications/mechanics/ICs_and_BCs.h index 01addd42e..9312ad115 100644 --- a/applications/mechanics/ICs_and_BCs.h +++ b/applications/mechanics/ICs_and_BCs.h @@ -49,7 +49,11 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- @@ -66,8 +70,11 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. - + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/nucleationModel/ICs_and_BCs.h b/applications/nucleationModel/ICs_and_BCs.h index ea69511e3..c4a06f876 100644 --- a/applications/nucleationModel/ICs_and_BCs.h +++ b/applications/nucleationModel/ICs_and_BCs.h @@ -67,11 +67,14 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // -------------------------------------------------------------------------- // ENTER THE NON-UNIFORM DIRICHLET BOUNDARY CONDITIONS HERE FOR SCALAR FIELDS // -------------------------------------------------------------------------- - // Enter the function describing conditions for the fields at point "p". + // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. - + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- return scalar_BC; @@ -84,10 +87,14 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // -------------------------------------------------------------------------- // ENTER THE NON-UNIFORM DIRICHLET BOUNDARY CONDITIONS HERE FOR VECTOR FIELDS // -------------------------------------------------------------------------- - // Enter the function describing conditions for the fields at point "p". + // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/precipitateEvolution/ICs_and_BCs.h b/applications/precipitateEvolution/ICs_and_BCs.h index 97adc3f3c..1c07f8188 100644 --- a/applications/precipitateEvolution/ICs_and_BCs.h +++ b/applications/precipitateEvolution/ICs_and_BCs.h @@ -75,7 +75,11 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- @@ -92,7 +96,11 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/precipitateEvolution_pfunction/ICs_and_BCs.h b/applications/precipitateEvolution_pfunction/ICs_and_BCs.h index 97adc3f3c..1c07f8188 100644 --- a/applications/precipitateEvolution_pfunction/ICs_and_BCs.h +++ b/applications/precipitateEvolution_pfunction/ICs_and_BCs.h @@ -75,7 +75,11 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- @@ -92,7 +96,11 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/preferential_nucleationModel/ICs_and_BCs.h b/applications/preferential_nucleationModel/ICs_and_BCs.h index 6c4eea8a0..04406c9fa 100644 --- a/applications/preferential_nucleationModel/ICs_and_BCs.h +++ b/applications/preferential_nucleationModel/ICs_and_BCs.h @@ -66,10 +66,14 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // -------------------------------------------------------------------------- // ENTER THE NON-UNIFORM DIRICHLET BOUNDARY CONDITIONS HERE FOR SCALAR FIELDS // -------------------------------------------------------------------------- - // Enter the function describing conditions for the fields at point "p". + // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- @@ -83,10 +87,14 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // -------------------------------------------------------------------------- // ENTER THE NON-UNIFORM DIRICHLET BOUNDARY CONDITIONS HERE FOR VECTOR FIELDS // -------------------------------------------------------------------------- - // Enter the function describing conditions for the fields at point "p". + // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- diff --git a/applications/singlePrecipitateKKS/ICs_and_BCs.h b/applications/singlePrecipitateKKS/ICs_and_BCs.h index 53cc31ea8..35298b1fe 100644 --- a/applications/singlePrecipitateKKS/ICs_and_BCs.h +++ b/applications/singlePrecipitateKKS/ICs_and_BCs.h @@ -81,7 +81,11 @@ double NonUniformDirichletBC::value (const dealii::Point &p, const uns // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // ------------------------------------------------------------------------- @@ -98,7 +102,11 @@ void NonUniformDirichletBCVec::vector_value (const dealii::Point &p, d // Enter the function describing conditions for the fields at point "p". // Use "if" statements to set the boundary condition for each variable // according to its variable index. This function can be left blank if there - // are no non-uniform Dirichlet boundary conditions. + // are no non-uniform Dirichlet boundary conditions. For BCs that change in + // time, you can access the current time through the variable "time". The + // boundary index can be accessed via the variable "direction", which starts + // at zero and uses the same order as the BC specification in parameters.in + // (i.e. left = 0, right = 1, bottom = 2, top = 3, front = 4, back = 5). // -------------------------------------------------------------------------