Skip to content

weak formulation and layer thickness #2

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

Open
jiangwen84 wants to merge 2 commits into
base: moving_interface_xfem
Choose a base branch
from
Open

weak formulation and layer thickness #2

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
97be48e
weak formulation and layer thickness
jiangwen84 Dec 17, 2019
9271a16
add moving interface mechanics problem
jiangwen84 Dec 19, 2019
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
38 changes: 38 additions & 0 deletions modules/xfem/include/postprocessors/OxideLayerThickness.h
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,38 @@
//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#ifndef OXIDELAYERTHICKNESS_H
#define OXIDELAYERTHICKNESS_H

#include "ElementPostprocessor.h"
#include "MovingLineSegmentCutSetUserObject.h"
// Forward Declarations
class OxideLayerThickness;

template <>
InputParameters validParams<OxideLayerThickness>();

class OxideLayerThickness : public ElementPostprocessor
{
public:
OxideLayerThickness(const InputParameters & parameters);

void initialize() override {}
void execute() override{};
void finalize() override{};
void threadJoin(const UserObject & user_object) override{};

virtual Real getValue() override;

protected:
const MovingLineSegmentCutSetUserObject * const _moving_line_segments;
const unsigned int _cut_data_index;
};

#endif // OXIDELAYERTHICKNESS_H
2 changes: 2 additions & 0 deletions modules/xfem/include/userobjects/MovingLineSegmentCutSetUserObject.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -37,6 +37,8 @@ class MovingLineSegmentCutSetUserObject : public LineSegmentCutSetUserObject

virtual Real cutFraction(unsigned int cut_num, Real time) const override;

virtual Real getCutDataXCoord(unsigned int i) const { return _cut_data[i * 6]; };

/// Pointer to XFEMMovingInterfaceVelocityBase object
const XFEMMovingInterfaceVelocityBase * _interface_velocity;
};
Expand Down
37 changes: 37 additions & 0 deletions modules/xfem/src/postprocessors/OxideLayerThickness.C
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,37 @@
//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#include "OxideLayerThickness.h"

registerMooseObject("XFEMApp", OxideLayerThickness);

template <>
InputParameters
validParams<OxideLayerThickness>()
{
InputParameters params = validParams<ElementPostprocessor>();
params.addRequiredParam<UserObjectName>("moving_line_segments",
"The MovingLineSegmentCutSetUserObject user object name");
params.addParam<unsigned int>("cut_data_index", 0, "The index of the cut data");
return params;
}

OxideLayerThickness::OxideLayerThickness(const InputParameters & parameters)
: ElementPostprocessor(parameters),
_moving_line_segments(
&getUserObject<MovingLineSegmentCutSetUserObject>("moving_line_segments")),
_cut_data_index(getParam<unsigned int>("cut_data_index"))
{
}

Real
OxideLayerThickness::getValue()
{
return 6e-4 - _moving_line_segments->getCutDataXCoord(_cut_data_index);
}
15 changes: 9 additions & 6 deletions modules/xfem/src/userobjects/XFEMC4OxideVelocity.C
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -49,7 +49,7 @@ XFEMC4OxideVelocity::computeMovingInterfaceVelocity(unsigned int point_id) const

Real delta = std::abs(xt - _x0);

// std::cout << "delta: " << delta << std::endl;
// std::cout << "delta: " << delta << std::endl;

// Current implementation only supports the case that the interface is moving horizontally
// return std::abs((_diffusivity_at_positive_level_set * grad_positive(0) -
Expand All @@ -72,8 +72,8 @@ XFEMC4OxideVelocity::computeMovingInterfaceVelocity(unsigned int point_id) const
const Real con_e_ox_m = 2 * con_v_ox_m;
const Real con_o_ox_m(5.24e28);
const Real con_o_m_ox(1.7078e28);
// use fixed temperature at 1200C until we add a temperature diffusion kernel

// use fixed temperature at 1200C until we add a temperature diffusion kernel
const Real temperature(1473);

const Real mobil_v = 4 * pow(migr_jp_l, 2) * migr_jp_f * 2 / (Kb * temperature) *
Expand All @@ -87,10 +87,13 @@ XFEMC4OxideVelocity::computeMovingInterfaceVelocity(unsigned int point_id) const
const Real eta = (-B - sqrt(pow(B, 2) - 4 * A * C)) / (2 * A);
const Real potential = Kb * temperature * log(eta);

const Real J_v = mobil_v * potential * (con_v_ox_w - con_v_ox_m * pow(eta, 2)) /
(1 - pow(eta, 2)) / delta;
const Real J_v =
mobil_v * potential * (con_v_ox_w - con_v_ox_m * pow(eta, 2)) / (1 - pow(eta, 2)) / delta;
// const Real J_o = zircaloy_density * Na / (zirconium_atmass * 1e-3) *
// _diffusivity_at_positive_level_set * (grad_positive(0) + grad_positive(1) +
// grad_positive(2)) / 3 * 1e-6;
const Real J_o = zircaloy_density * Na / (zirconium_atmass * 1e-3) *
_diffusivity_at_positive_level_set * (grad_positive(0) + grad_positive(1) + grad_positive(2)) / 3 * 1e-6;
_diffusivity_at_positive_level_set * (grad_positive(0)) * 1e-6;

if (delta == 0)
return sqrt(0.01126 * exp(-35890 / (1.987 * temperature)) / (2 * _t)) * (-1e-2);
Expand Down
177 changes: 177 additions & 0 deletions modules/xfem/test/tests/moving_interface/moving_oxide_C4_weak.i
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,177 @@
# test for an oxide growing on top of a ziroconium nuclear fuel cladding
# using the C4 model to compute the growth rate

[GlobalParams]
order = FIRST
family = LAGRANGE
[]

[Mesh]
type = GeneratedMesh
dim = 2
nx = 21
ny = 20
xmin = 0
xmax = 6e-4
ymin = 0
ymax = 1e-3
elem_type = QUAD4
[]

[XFEM]
qrule = volfrac
output_cut_plane = true
[]

[UserObjects]
[./velocity]
type = XFEMC4OxideVelocity
diffusivity_at_positive_level_set = 0.2978e-5 #1e-5*0.2978
diffusivity_at_negative_level_set = 0.6667e-11 #1e-11*0.6667
equilibrium_concentration_jump = 0.3689
value_at_interface_uo = value_uo
x0 = 5.9e-4
[../]
[./value_uo]
type = PointValueAtXFEMInterface
variable = 'u'
geometric_cut_userobject = 'moving_line_segments'
execute_on = 'nonlinear'
level_set_var = ls
[../]
[./moving_line_segments]
type = MovingLineSegmentCutSetUserObject
cut_data = '5.9e-4 0 5.9e-4 1e-3 0 0'
heal_always = true
interface_velocity = velocity
[../]
[]

[Variables]
[./u]
[../]
[]

[ICs]
[./ic_u]
type = FunctionIC
variable = u
function = 'if(x<5.9e-4, 0.03, 0.6667)'
[../]
[]

[AuxVariables]
[./ls]
order = FIRST
family = LAGRANGE
[../]
[]

# Need to use XFEMTwoSideDirichlet that has been removed
[Constraints]
[./u_constraint]
type = XFEMTwoSideDirichlet
geometric_cut_userobject = 'moving_line_segments'
use_displaced_mesh = false
variable = u
value_at_positive_level_set_interface = 1 #0.2978
value_at_negative_level_set_interface = 1 #0.6667
# value = 0.6667
alpha = 1e5
[../]
[]

[Kernels]
[./diff]
type = ConcentrationDiffusion
variable = u
diffusion_coefficient_name = 'diffusion_coefficient'
[../]
[./time]
type = TimeDerivative
variable = u
[../]
[]

[AuxKernels]
[./ls]
type = LineSegmentLevelSetAux
line_segment_cut_set_user_object = 'moving_line_segments'
variable = ls
[../]
[]

[Materials]
[./diffusivity_A]
type = GenericConstantMaterial
prop_names = A_diffusion_coefficient
prop_values = 1e-5
[../]
[./diffusivity_B]
type = GenericConstantMaterial
prop_names = B_diffusion_coefficient
prop_values = 1e-11
[../]
[./diff_combined]
type = LevelSetBiMaterialReal
levelset_positive_base = 'A'
levelset_negative_base = 'B'
level_set_var = ls
prop_name = diffusion_coefficient
[../]
[]

[BCs]
# Define boundary conditions
[./left_u]
type = DirichletBC
variable = u
value = 0.0450 #0.03/0.6667
boundary = left
[../]

[./right_u]
type = DirichletBC
variable = u
value = 2.2388 #0.6667/0.2978
boundary = right
[../]
[]

[Postprocessors]
[./cut_data_x]
type = OxideLayerThickness
moving_line_segments = moving_line_segments
execute_on = TIMESTEP_END
[../]
[]

[Executioner]
type = Transient
solve_type = 'PJFNK'
petsc_options_iname = '-pc_type'
petsc_options_value = 'lu'
line_search = 'none'

l_tol = 1e-3
nl_max_its = 15
nl_rel_tol = 1e-7
nl_abs_tol = 1e-7

start_time = 0.0
dt = 10
num_steps = 150
max_xfem_update = 1
[]


[Outputs]
execute_on = timestep_end
exodus = true
[./console]
type = Console
output_linear = true
[../]
csv = true
perf_graph = true
[]
Loading