NumbatDensity

This Material calculates the density of the fluid with given concentration

The density is given by $(1)var element = document.getElementById("moose-equation-a4f7d572-4028-4b53-ac3f-5f03929939a7");katex.render("\\rho(c) = \\rho_0 + \\frac{c \\Delta \\rho}{c_0},", element, {displayMode:true,throwOnError:false});$ where $var element = document.getElementById("moose-equation-da8d8c4b-5723-43f5-abe5-fb285793ca34");katex.render("\\rho_0", element, {displayMode:false,throwOnError:false});$ is the density at zero solute concentration, $var element = document.getElementById("moose-equation-fc2d56cf-819f-4160-8eaa-e8aad094949e");katex.render("c", element, {displayMode:false,throwOnError:false});$ is the concentration of solute, $var element = document.getElementById("moose-equation-420c7a2f-7237-43d7-ab7a-d702151535c7");katex.render("\\Delta \\rho", element, {displayMode:false,throwOnError:false});$ is the density increase at full saturation, and $var element = document.getElementById("moose-equation-f06227c0-4f96-440e-b559-02a4a969eceb");katex.render("c_0", element, {displayMode:false,throwOnError:false});$ is the equilibrium solute concentration.

These parameters are specified in the input file as zero_density ( $var element = document.getElementById("moose-equation-c53af6e0-9f67-4442-a4fc-cbb88230a1c5");katex.render("\\rho_0", element, {displayMode:false,throwOnError:false});$ ), delta_density ( $var element = document.getElementById("moose-equation-22986142-e322-4021-92dd-70a31b40bc48");katex.render("\\Delta \\rho", element, {displayMode:false,throwOnError:false});$ ), and saturated_concentration ( $var element = document.getElementById("moose-equation-5bd74f32-d5c5-460c-b533-6b357488f1c3");katex.render("c_0", element, {displayMode:false,throwOnError:false});$ ).

Input Parameters

concentrationThe concentration variable
C++ Type:std::vector
Description:The concentration variable
delta_densityThe density increase due to concentration at equilibrium
C++ Type:double
Description:The density increase due to concentration at equilibrium
zero_densityThe density with zero concentration
C++ Type:double
Description:The density with zero concentration

Required Parameters

saturated_concentration1Concentration of species at full saturation. Note: should be equal to boundary condition concentration
Default:1
C++ Type:double
Description:Concentration of species at full saturation. Note: should be equal to boundary condition concentration
boundaryThe list of boundary IDs from the mesh where this boundary condition applies
C++ Type:std::vector
Description:The list of boundary IDs from the mesh where this boundary condition applies
computeTrueWhen false, MOOSE will not call compute methods on this material. The user must call computeProperties() after retrieving the Material via MaterialPropertyInterface::getMaterial(). Non-computed Materials are not sorted for dependencies.
Default:True
C++ Type:bool
Description:When false, MOOSE will not call compute methods on this material. The user must call computeProperties() after retrieving the Material via MaterialPropertyInterface::getMaterial(). Non-computed Materials are not sorted for dependencies.
blockThe list of block ids (SubdomainID) that this object will be applied
C++ Type:std::vector
Description:The list of block ids (SubdomainID) that this object will be applied

Optional Parameters

enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Description:Set the enabled status of the MooseObject.
use_displaced_meshFalseWhether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.
Default:False
C++ Type:bool
Description:Whether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.
control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector
Description:Adds user-defined labels for accessing object parameters via control logic.
seed0The seed for the master random number generator
Default:0
C++ Type:unsigned int
Description:The seed for the master random number generator
implicitTrueDetermines whether this object is calculated using an implicit or explicit form
Default:True
C++ Type:bool
Description:Determines whether this object is calculated using an implicit or explicit form
constant_onNONEWhen ELEMENT, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps.When SUBDOMAIN, MOOSE will only call computeSubdomainProperties() for the 0th quadrature point, and then copy that value to the other qps. Evaluations on element qps will be skipped
Default:NONE
C++ Type:MooseEnum
Description:When ELEMENT, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps.When SUBDOMAIN, MOOSE will only call computeSubdomainProperties() for the 0th quadrature point, and then copy that value to the other qps. Evaluations on element qps will be skipped

Advanced Parameters

output_propertiesList of material properties, from this material, to output (outputs must also be defined to an output type)
C++ Type:std::vector
Description:List of material properties, from this material, to output (outputs must also be defined to an output type)
outputsnone Vector of output names were you would like to restrict the output of variables(s) associated with this object
Default:none
C++ Type:std::vector
Description:Vector of output names were you would like to restrict the output of variables(s) associated with this object

Outputs Parameters

Input Parameters
Input Files

# Sample input file for 2D density-driven convective mixing # Note: this input file is primarily used to generate documentation # PLEASE DO NOT ALTER THIS FILE [Mesh] type = NumbatBiasedMesh dim = 2 ymax = 1.5 nx = 100 ny = 50 refined_edge = top refined_resolution = 0.001 [] [Variables] [./concentration] initial_condition = 0 scaling = 1e6 [../] [./pressure] initial_condition = 1e6 [../] [] [AuxVariables] [./noise] family = MONOMIAL order = CONSTANT [../] [] [ICs] [./noise] type = RandomIC variable = noise max = 0.003 min = -0.003 [../] [] [Materials] [./porosity] type = NumbatPorosity porosity = 0.3 noise = noise [../] [./permeability] type = NumbatPermeability permeability = '1e-11 0 0 0 1e-11 0 0 0 1e-11' [../] [./diffusivity] type = NumbatDiffusivity diffusivity = 2e-9 [../] [./density] type = NumbatDensity concentration = concentration zero_density = 995 delta_density = 10.5 saturated_concentration = 0.049306 [../] [./viscosity] type = NumbatViscosity viscosity = 6e-4 [../] [] [Kernels] [./time] type = NumbatTimeDerivative variable = concentration [../] [./diffusion] type = NumbatDiffusion variable = concentration [../] [./convection] type = NumbatConvection variable = concentration pressure = pressure [../] [./darcy] type = NumbatDarcy variable = pressure concentration = concentration [../] [] [BCs] [./conctop] type = PresetBC variable = concentration boundary = top value = 0.049306 [../] [./Periodic] [./x] variable = concentration auto_direction = x [../] [../] [] [Preconditioning] [./smp] type = SMP full = true [../] [] [Executioner] type = Transient l_max_its = 200 end_time = 3e5 solve_type = NEWTON petsc_options = -ksp_snes_ew petsc_options_iname = '-pc_type -sub_pc_type -ksp_atol' petsc_options_value = 'bjacobi ilu 1e-12' nl_abs_tol = 1e-10 nl_max_its = 25 dtmax = 2e3 [./TimeStepper] type = IterationAdaptiveDT dt = 1 [../] [] [Postprocessors] [./boundaryfluxint] type = NumbatSideFlux variable = concentration boundary = top [../] [./mass] type = NumbatTotalMass variable = concentration [../] [] [Outputs] perf_graph = true [./exodus] type = Exodus file_base = 2Dddc execute_on = 'INITIAL TIMESTEP_END' [../] [./csvoutput] type = CSV file_base = 2Dddc execute_on = 'INITIAL TIMESTEP_END' [../] []

# Sample input file for 3D density-driven convective mixing # Note: this input file is primarily used to generate documentation # PLEASE DO NOT ALTER THIS FILE [Mesh] type = NumbatBiasedMesh dim = 3 zmax = 1.5 nx = 20 ny = 20 nz = 500 refined_edge = front refined_resolution = 0.001 [] [Variables] [./concentration] initial_condition = 0 scaling = 1e6 [../] [./pressure] initial_condition = 1e6 [../] [] [AuxVariables] [./noise] family = MONOMIAL order = CONSTANT [../] [] [ICs] [./noise] type = RandomIC variable = noise max = 0.003 min = -0.003 [../] [] [Materials] [./porosity] type = NumbatPorosity porosity = 0.3 noise = noise [../] [./permeability] type = NumbatPermeability permeability = '1e-11 0 0 0 1e-11 0 0 0 1e-11' [../] [./diffusivity] type = NumbatDiffusivity diffusivity = 2e-9 [../] [./density] type = NumbatDensity concentration = concentration zero_density = 995 delta_density = 10.5 saturated_concentration = 0.049306 [../] [./viscosity] type = NumbatViscosity viscosity = 6e-4 [../] [] [Kernels] [./time] type = NumbatTimeDerivative variable = concentration [../] [./diffusion] type = NumbatDiffusion variable = concentration [../] [./convection] type = NumbatConvection variable = concentration pressure = pressure [../] [./darcy] type = NumbatDarcy variable = pressure concentration = concentration [../] [] [BCs] [./conctop] type = PresetBC variable = concentration boundary = front value = 0.049306 [../] [./Periodic] [./x] variable = concentration auto_direction = 'x y' [../] [../] [] [Preconditioning] [./smp] type = SMP full = true [../] [] [Executioner] type = Transient l_max_its = 200 end_time = 3e5 solve_type = NEWTON petsc_options = -ksp_snes_ew petsc_options_iname = '-pc_type -sub_pc_type -ksp_atol' petsc_options_value = 'bjacobi ilu 1e-12' nl_abs_tol = 1e-10 nl_max_its = 25 dtmax = 2e3 [./TimeStepper] type = IterationAdaptiveDT dt = 1 [../] [] [Postprocessors] [./boundaryfluxint] type = NumbatSideFlux variable = concentration boundary = top [../] [./mass] type = NumbatTotalMass variable = concentration [../] [] [Outputs] perf_graph = true [./exodus] type = Exodus file_base = 3Dddc execute_on = 'INITIAL TIMESTEP_END' [../] [./csvoutput] type = CSV file_base = 3Dddc execute_on = 'INITIAL TIMESTEP_END' [../] []

NumbatDensity

Input Parameters

Required Parameters

Optional Parameters

Advanced Parameters

Outputs Parameters

Input Files

tests/3D/3D_action.i

examples/3D/isotropic/3Dddc.i

tests/darcy/3Dsteady.i

tests/convection_diffusion/2DtransientDG.i

tests/actions/2D.i

tests/convection_diffusion/2Dtransient2.i

tests/2D/2D_c0_action.i

examples/inputfiles/2D.i

tests/1D/1D_error.i

tests/2D/2D_action.i

tests/convection/2Dtransient_heavy.i

tests/convection_diffusion/2Dtransient.i

tests/2D/2D_c0.i

tests/convection/3Dtransient_heavy.i

examples/2D/anisotropic/2Dddc.i

tests/darcy/2Dsteady.i

tests/jacobians/3D.i

tests/darcy/3Dsteady_gravity.i

examples/inputfiles/3D.i

tests/convection/2Dtransient.i

tests/convection/2DtransientDG.i

tests/convection/3DtransientDG.i

tests/3D/3D.i

tests/jacobians/2D.i

examples/2D/isotropic/2Dddc.i

tests/actions/3D.i

tests/jacobians/2DDG.i

tests/darcy/2Dsteady_gravity.i

tests/convection_diffusion/2DtransientDG2.i

tests/2D/2D.i

tests/convection/3Dtransient.i