diff --git a/examples/laser_welding_problems/physics/multi-scale-grain-growth/grain_growth_2D_bottom.i b/examples/laser_welding_problems/physics/multi-scale-grain-growth/grain_growth_2D_bottom.i index 82dce52..4f84fcb 100644 --- a/examples/laser_welding_problems/physics/multi-scale-grain-growth/grain_growth_2D_bottom.i +++ b/examples/laser_welding_problems/physics/multi-scale-grain-growth/grain_growth_2D_bottom.i @@ -9,11 +9,11 @@ dim = 2 # Problem dimension nx = 11 # Number of elements in the x-direction ny = 11 # Number of elements in the y-direction - xmin = 0 # minimum x-coordinate of the mesh + xmin = 0 # minimum x-coordinate of the mesh xmax = 1e5 # maximum x-coordinate of the mesh - ymin = 0 # minimum y-coordinate of the mesh + ymin = 0 # minimum y-coordinate of the mesh ymax = 1e5 # maximum y-coordinate of the mesh - elem_type = QUAD4 # Type of elements used in the mesh + elem_type = QUAD4 # Type of elements used in the mesh uniform_refine = 3 # Initial uniform refinement of the mesh parallel_type = replicated # Periodic BCs @@ -27,123 +27,123 @@ [Variables] # Variable block, where all variables in the simulation are declared - [./PolycrystalVariables] - [../] + [PolycrystalVariables] + [] [] [UserObjects] - [./voronoi] + [voronoi] type = PolycrystalVoronoi grain_num = 100 # Number of grains rand_seed = 20 - [../] - [./grain_tracker] + [] + [grain_tracker] type = GrainTracker threshold = 0.2 connecting_threshold = 0.08 compute_halo_maps = true # Only necessary for displaying HALOS - [../] + [] [] [ICs] - [./PolycrystalICs] - [./PolycrystalColoringIC] + [PolycrystalICs] + [PolycrystalColoringIC] polycrystal_ic_uo = voronoi - [../] - [../] + [] + [] [] [AuxVariables] # Dependent variables - [./bnds] + [bnds] # Variable used to visualize the grain boundaries in the simulation - [../] - [./unique_grains] + [] + [unique_grains] order = CONSTANT family = MONOMIAL - [../] - [./var_indices] + [] + [var_indices] order = CONSTANT family = MONOMIAL - [../] - [./ghost_regions] + [] + [ghost_regions] order = CONSTANT family = MONOMIAL - [../] - [./halos] + [] + [halos] order = CONSTANT family = MONOMIAL - [../] - [./T_1000K] + [] + [T_1000K] initial_condition = 0.298 - [../] - [./T] + [] + [T] initial_condition = 298 - [../] + [] [] [Kernels] # Kernel block, where the kernels defining the residual equations are set up. - [./PolycrystalKernel] + [PolycrystalKernel] # Custom action creating all necessary kernels for grain growth. All input parameters are up in GlobalParams - [../] + [] [] [AuxKernels] # AuxKernel block, defining the equations used to calculate the auxvars - [./bnds_aux] + [bnds_aux] # AuxKernel that calculates the GB term type = BndsCalcAux variable = bnds execute_on = 'initial timestep_end' - [../] - [./unique_grains] + [] + [unique_grains] type = FeatureFloodCountAux variable = unique_grains flood_counter = grain_tracker field_display = UNIQUE_REGION execute_on = 'initial timestep_end' - [../] - [./var_indices] + [] + [var_indices] type = FeatureFloodCountAux variable = var_indices flood_counter = grain_tracker field_display = VARIABLE_COLORING execute_on = 'initial timestep_end' - [../] - [./ghosted_entities] + [] + [ghosted_entities] type = FeatureFloodCountAux variable = ghost_regions flood_counter = grain_tracker field_display = GHOSTED_ENTITIES execute_on = 'initial timestep_end' - [../] - [./halos] + [] + [halos] type = FeatureFloodCountAux variable = halos flood_counter = grain_tracker field_display = HALOS execute_on = 'initial timestep_end' - [../] - [./temp_conversion] + [] + [temp_conversion] type = ParsedAux coupled_variables = 'T_1000K' expression = '1000 * T_1000K' variable = T - [../] + [] [] [BCs] # Boundary Condition block - [./Periodic] - [./top_bottom] + [Periodic] + [top_bottom] auto_direction = 'x y' - [../] - [../] + [] + [] [] [Materials] - [./CuGrGr] + [CuGrGr] # Material properties type = GBEvolution T = T # Specified temperature profile @@ -151,18 +151,18 @@ GBmob0 = 3.19e-4 #m^4(Js) for copper from Schoenfelder1997 Q = 0.23 #eV for copper from Schoenfelder1997 GBenergy = 0.708 #J/m^2 from Schoenfelder1997 - [../] + [] [] [Postprocessors] # Scalar postprocessors - [./dt] + [dt] # Outputs the current time step type = TimestepSize - [../] - [./num_dofs] + [] + [num_dofs] type = NumDOFs - [../] + [] [] [Preconditioning] @@ -191,26 +191,26 @@ start_time = 0.0 end_time = 4800 - [./TimeStepper] + [TimeStepper] type = IterationAdaptiveDT dt = 25 # Initial time step. In this simulation it changes. optimal_iterations = 8 # Time step will adapt to maintain this number of nonlinear iterations - [../] + [] - [./Adaptivity] + [Adaptivity] # Block that turns on mesh adaptivity. Note that mesh will never coarsen beyond initial mesh (before uniform refinement) initial_adaptivity = 2 # Number of times mesh is adapted to initial condition refine_fraction = 0.7 # Fraction of high error that will be refined coarsen_fraction = 0.1 # Fraction of low error that will coarsened max_h_level = 4 # Max number of refinements used, starting from initial mesh (before uniform refinement) - [../] + [] [] [Outputs] exodus = true # Exodus file will be outputted csv = true - [./console] + [console] type = Console max_rows = 20 # Will print the 20 most recent postprocessor values to the screen - [../] + [] [] diff --git a/examples/laser_welding_problems/physics/multi-scale-grain-growth/grain_growth_2D_middle.i b/examples/laser_welding_problems/physics/multi-scale-grain-growth/grain_growth_2D_middle.i index 38276db..ec5adc7 100644 --- a/examples/laser_welding_problems/physics/multi-scale-grain-growth/grain_growth_2D_middle.i +++ b/examples/laser_welding_problems/physics/multi-scale-grain-growth/grain_growth_2D_middle.i @@ -9,11 +9,11 @@ dim = 2 # Problem dimension nx = 11 # Number of elements in the x-direction ny = 11 # Number of elements in the y-direction - xmin = 0 # minimum x-coordinate of the mesh + xmin = 0 # minimum x-coordinate of the mesh xmax = 1e5 # maximum x-coordinate of the mesh - ymin = 0 # minimum y-coordinate of the mesh + ymin = 0 # minimum y-coordinate of the mesh ymax = 1e5 # maximum y-coordinate of the mesh - elem_type = QUAD4 # Type of elements used in the mesh + elem_type = QUAD4 # Type of elements used in the mesh uniform_refine = 3 # Initial uniform refinement of the mesh parallel_type = replicated # Periodic BCs @@ -27,123 +27,123 @@ [Variables] # Variable block, where all variables in the simulation are declared - [./PolycrystalVariables] - [../] + [PolycrystalVariables] + [] [] [UserObjects] - [./voronoi] + [voronoi] type = PolycrystalVoronoi grain_num = 100 # Number of grains rand_seed = 10 - [../] - [./grain_tracker] + [] + [grain_tracker] type = GrainTracker threshold = 0.2 connecting_threshold = 0.08 compute_halo_maps = true # Only necessary for displaying HALOS - [../] + [] [] [ICs] - [./PolycrystalICs] - [./PolycrystalColoringIC] + [PolycrystalICs] + [PolycrystalColoringIC] polycrystal_ic_uo = voronoi - [../] - [../] + [] + [] [] [AuxVariables] # Dependent variables - [./bnds] + [bnds] # Variable used to visualize the grain boundaries in the simulation - [../] - [./unique_grains] + [] + [unique_grains] order = CONSTANT family = MONOMIAL - [../] - [./var_indices] + [] + [var_indices] order = CONSTANT family = MONOMIAL - [../] - [./ghost_regions] + [] + [ghost_regions] order = CONSTANT family = MONOMIAL - [../] - [./halos] + [] + [halos] order = CONSTANT family = MONOMIAL - [../] - [./T_1000K] + [] + [T_1000K] initial_condition = 0.298 - [../] - [./T] + [] + [T] initial_condition = 298 - [../] + [] [] [Kernels] # Kernel block, where the kernels defining the residual equations are set up. - [./PolycrystalKernel] + [PolycrystalKernel] # Custom action creating all necessary kernels for grain growth. All input parameters are up in GlobalParams - [../] + [] [] [AuxKernels] # AuxKernel block, defining the equations used to calculate the auxvars - [./bnds_aux] + [bnds_aux] # AuxKernel that calculates the GB term type = BndsCalcAux variable = bnds execute_on = 'initial timestep_end' - [../] - [./unique_grains] + [] + [unique_grains] type = FeatureFloodCountAux variable = unique_grains flood_counter = grain_tracker field_display = UNIQUE_REGION execute_on = 'initial timestep_end' - [../] - [./var_indices] + [] + [var_indices] type = FeatureFloodCountAux variable = var_indices flood_counter = grain_tracker field_display = VARIABLE_COLORING execute_on = 'initial timestep_end' - [../] - [./ghosted_entities] + [] + [ghosted_entities] type = FeatureFloodCountAux variable = ghost_regions flood_counter = grain_tracker field_display = GHOSTED_ENTITIES execute_on = 'initial timestep_end' - [../] - [./halos] + [] + [halos] type = FeatureFloodCountAux variable = halos flood_counter = grain_tracker field_display = HALOS execute_on = 'initial timestep_end' - [../] - [./temp_conversion] + [] + [temp_conversion] type = ParsedAux coupled_variables = 'T_1000K' expression = '1000 * T_1000K' variable = T - [../] + [] [] [BCs] # Boundary Condition block - [./Periodic] - [./top_bottom] + [Periodic] + [top_bottom] auto_direction = 'x y' - [../] - [../] + [] + [] [] [Materials] - [./CuGrGr] + [CuGrGr] # Material properties type = GBEvolution T = T # Specified temperature profile @@ -151,18 +151,18 @@ GBmob0 = 3.19e-4 #m^4(Js) for copper from Schoenfelder1997 Q = 0.23 #eV for copper from Schoenfelder1997 GBenergy = 0.708 #J/m^2 from Schoenfelder1997 - [../] + [] [] [Postprocessors] # Scalar postprocessors - [./dt] + [dt] # Outputs the current time step type = TimestepSize - [../] - [./num_dofs] + [] + [num_dofs] type = NumDOFs - [../] + [] [] [Preconditioning] @@ -191,26 +191,26 @@ start_time = 0.0 end_time = 4800 - [./TimeStepper] + [TimeStepper] type = IterationAdaptiveDT dt = 25 # Initial time step. In this simulation it changes. optimal_iterations = 8 # Time step will adapt to maintain this number of nonlinear iterations - [../] + [] - [./Adaptivity] + [Adaptivity] # Block that turns on mesh adaptivity. Note that mesh will never coarsen beyond initial mesh (before uniform refinement) initial_adaptivity = 2 # Number of times mesh is adapted to initial condition refine_fraction = 0.7 # Fraction of high error that will be refined coarsen_fraction = 0.1 # Fraction of low error that will coarsened max_h_level = 4 # Max number of refinements used, starting from initial mesh (before uniform refinement) - [../] + [] [] [Outputs] exodus = true # Exodus file will be outputted csv = true - [./console] + [console] type = Console max_rows = 20 # Will print the 20 most recent postprocessor values to the screen - [../] + [] [] diff --git a/examples/laser_welding_problems/physics/multi-scale-grain-growth/larry-2d-one-spot-scaled-1st-order.i b/examples/laser_welding_problems/physics/multi-scale-grain-growth/larry-2d-one-spot-scaled-1st-order.i index a6e47ed..92d172f 100644 --- a/examples/laser_welding_problems/physics/multi-scale-grain-growth/larry-2d-one-spot-scaled-1st-order.i +++ b/examples/laser_welding_problems/physics/multi-scale-grain-growth/larry-2d-one-spot-scaled-1st-order.i @@ -1,15 +1,15 @@ -length_unit_exponent=-4 -temperature_unit_exponent=3 -mass_unit_exponent=-6 -time_unit_exponent=-5 - -endtime=${fparse .048 / 10^time_unit_exponent} -timestep=${fparse endtime / 1000} -surfacetemp=${fparse 300 / 10^temperature_unit_exponent} -bottomtemp=${fparse 300 / 10^temperature_unit_exponent} -width=${fparse 8e-4 / 10^length_unit_exponent} -half_width=${fparse width / 2} -pooldepth=${width} +length_unit_exponent = -4 +temperature_unit_exponent = 3 +mass_unit_exponent = -6 +time_unit_exponent = -5 + +endtime = '${fparse .048 / 10^time_unit_exponent}' +timestep = '${fparse endtime / 1000}' +surfacetemp = '${fparse 300 / 10^temperature_unit_exponent}' +bottomtemp = '${fparse 300 / 10^temperature_unit_exponent}' +width = '${fparse 8e-4 / 10^length_unit_exponent}' +half_width = '${fparse width / 2}' +pooldepth = ${width} [GlobalParams] gravity = '0 0 0' @@ -26,9 +26,9 @@ pooldepth=${width} [Mesh] type = GeneratedMesh dim = 2 - xmin = ${fparse -half_width} + xmin = '${fparse -half_width}' xmax = ${half_width} - ymin = ${fparse -pooldepth} + ymin = '${fparse -pooldepth}' ymax = 0 nx = 4 ny = 4 @@ -38,110 +38,110 @@ pooldepth=${width} [] [MeshModifiers] - [./corner_node] + [corner_node] type = AddExtraNodeset new_boundary = 'pinned_node' nodes = '0' - [../] + [] [] [Variables] - [./vel_x] - [./InitialCondition] + [vel_x] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] + [] scaling = 1e6 - [../] + [] - [./vel_y] - [./InitialCondition] + [vel_y] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] + [] scaling = 1e6 - [../] + [] - [./T] + [T] scaling = 1e1 - [../] + [] - [./p] + [p] order = FIRST scaling = 1e4 - [../] - [./disp_x] + [] + [disp_x] scaling = 1e-3 - [../] - [./disp_y] + [] + [disp_y] scaling = 1e-3 - [../] + [] [] [ICs] - [./T] + [T] type = FunctionIC variable = T function = '(${surfacetemp} - ${bottomtemp}) / ${pooldepth} * z + ${surfacetemp}' - [../] + [] [] [Kernels] [] [ADKernels] - [./disp_x] + [disp_x] type = ADStressDivergence variable = disp_x component = 0 - [../] - [./disp_y] + [] + [disp_y] type = ADStressDivergence variable = disp_y component = 1 - [../] + [] -[./mesh_x] + [mesh_x] type = INSConvectedMesh variable = vel_x disp_x = disp_x disp_y = disp_y use_displaced_mesh = true - [../] - [./mesh_y] + [] + [mesh_y] type = INSConvectedMesh variable = vel_y disp_x = disp_x disp_y = disp_y use_displaced_mesh = true - [../] - [./mesh_T] + [] + [mesh_T] type = INSTemperatureConvectedMesh variable = T disp_x = disp_x disp_y = disp_y use_displaced_mesh = true - [../] + [] # mass - [./mass] + [mass] type = INSADMass variable = p u = vel_x v = vel_y p = p use_displaced_mesh = true - [../] + [] # x-momentum, time - [./x_momentum_time] + [x_momentum_time] type = INSADMomentumTimeDerivative variable = vel_x use_displaced_mesh = true - [../] + [] # x-momentum, space - [./x_momentum_space] + [x_momentum_space] type = INSADMomentumBase variable = vel_x u = vel_x @@ -149,17 +149,17 @@ pooldepth=${width} p = p component = 0 use_displaced_mesh = true - [../] + [] # y-momentum, time - [./y_momentum_time] + [y_momentum_time] type = INSADMomentumTimeDerivative variable = vel_y use_displaced_mesh = true - [../] + [] # y-momentum, space - [./y_momentum_space] + [y_momentum_space] type = INSADMomentumBase variable = vel_y u = vel_x @@ -167,70 +167,70 @@ pooldepth=${width} p = p component = 1 use_displaced_mesh = true - [../] - - # temperature - [./temperature_time] - type = INSADTemperatureTimeDerivative - variable = T - use_displaced_mesh = true - [../] - - [./temperature_space] - type = INSADTemperature - variable = T - u = vel_x - v = vel_y - p = p - use_displaced_mesh = true - [../] + [] + + # temperature + [temperature_time] + type = INSADTemperatureTimeDerivative + variable = T + use_displaced_mesh = true + [] + + [temperature_space] + type = INSADTemperature + variable = T + u = vel_x + v = vel_y + p = p + use_displaced_mesh = true + [] [] [BCs] - [./x_no_disp] + [x_no_disp] type = DirichletBC variable = disp_x boundary = 'left right bottom' value = 0 - [../] - [./y_no_disp] + [] + [y_no_disp] type = DirichletBC variable = disp_y boundary = 'left right bottom' value = 0 - [../] + [] - [./x_no_slip] + [x_no_slip] type = DirichletBC variable = vel_x boundary = 'bottom left right' value = 0.0 - [../] + [] - [./y_no_slip] + [y_no_slip] type = DirichletBC variable = vel_y boundary = 'bottom left right' value = 0.0 - [../] + [] - [./T_cold] + [T_cold] type = DirichletBC variable = T boundary = 'bottom' value = ${bottomtemp} - [../] + [] - [./pressure_pin] + [pressure_pin] type = DirichletBC variable = p boundary = 'pinned_node' value = 0 - [../] + [] [] [ADBCs] - [./radiation_flux] + [radiation_flux] type = RadiationEnergyFluxBC variable = T boundary = 'top' @@ -238,36 +238,36 @@ pooldepth=${width} sb_constant = 'sb_constant' absorptivity = 'abs' use_displaced_mesh = true - [../] + [] - [./weld_flux] + [weld_flux] type = GaussianWeldEnergyFluxBC variable = T boundary = 'top' reff = 0.6 - F0 = ${fparse 5e8 / 10^mass_unit_exponent * 10^(3 * time_unit_exponent)} - R = ${fparse 1e-4 / 10^length_unit_exponent} + F0 = '${fparse 5e8 / 10^mass_unit_exponent * 10^(3 * time_unit_exponent)}' + R = '${fparse 1e-4 / 10^length_unit_exponent}' x_beam_coord = '${fparse -3e-4 / 10^length_unit_exponent} + ${fparse 1 / 10^length_unit_exponent * 10^time_unit_exponent / 60} * t' y_beam_coord = 0 z_beam_coord = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_x] + [vapor_recoil_x] type = VaporRecoilPressureMomentumFluxBC variable = vel_x boundary = 'top' component = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_y] + [vapor_recoil_y] type = VaporRecoilPressureMomentumFluxBC variable = vel_y boundary = 'top' component = 1 use_displaced_mesh = true - [../] + [] # [./surface_x] # type = SurfaceTensionBC @@ -285,7 +285,7 @@ pooldepth=${width} # use_displaced_mesh = true # [../] - [./displace_x_top] + [displace_x_top] type = PenaltyDisplaceBoundaryBC boundary = 'top' variable = 'disp_x' @@ -293,8 +293,8 @@ pooldepth=${width} disp_x = 'disp_x' vel_y = 'vel_y' disp_y = 'disp_y' - [../] - [./displace_y_top] + [] + [displace_y_top] type = PenaltyDisplaceBoundaryBC boundary = 'top' variable = 'disp_y' @@ -302,11 +302,11 @@ pooldepth=${width} disp_x = 'disp_x' vel_y = 'vel_y' disp_y = 'disp_y' - [../] + [] [] [ADMaterials] - [./kc_fits] + [kc_fits] type = DemonstrationPlantFits temperature = T beta = 1e7 @@ -314,8 +314,8 @@ pooldepth=${width} temperature_unit_exponent = ${temperature_unit_exponent} mass_unit_exponent = ${mass_unit_exponent} time_unit_exponent = ${time_unit_exponent} - [../] - [./boundary] + [] + [boundary] type = DemonstrationPlantFitsBoundary boundary = 'top' temperature = T @@ -324,28 +324,28 @@ pooldepth=${width} temperature_unit_exponent = ${temperature_unit_exponent} mass_unit_exponent = ${mass_unit_exponent} time_unit_exponent = ${time_unit_exponent} - [../] - [./stress] + [] + [stress] type = PseudoSolidStress disp_x = disp_x disp_y = disp_y - [../] + [] [] [Materials] - [./const] + [const] type = GenericConstantMaterial prop_names = 'abs sb_constant' prop_values = '1 ${fparse 5.67e-8 / 10^mass_unit_exponent * 10^(4 * temperature_unit_exponent) * 10^(3 * time_unit_exponent)}' - [../] + [] [] [Preconditioning] - [./SMP] + [SMP] type = SMP full = true solve_type = 'NEWTON' - [../] + [] [] [Executioner] @@ -360,18 +360,18 @@ pooldepth=${width} # nl_rel_tol = 1e-6 line_search = 'none' nl_max_its = 12 - l_tol = 1e-2 - [./TimeStepper] + l_tol = 1e-2 + [TimeStepper] type = IterationAdaptiveDT optimal_iterations = 8 dt = ${timestep} linear_iteration_ratio = 1e6 growth_factor = 1.1 - [../] + [] [] [MultiApps] - [./grain_growth] + [grain_growth] type = TransientMultiApp app_type = MalamuteApp positions = '-1 -2 0 @@ -382,28 +382,28 @@ pooldepth=${width} max_catch_up_steps = 4 # sub_cycling = true # library_path = '../../../../../phase_field/lib' - [../] + [] [] [Transfers] - [./to_sub] + [to_sub] source_variable = T direction = to_multiapp variable = T_1000K type = MultiAppVariableValueSampleTransfer multi_app = grain_growth - [../] + [] [] [Outputs] print_linear_residuals = false csv = true perf_graph = true - [./exodus] + [exodus] type = Exodus output_material_properties = true show_material_properties = 'mu surface_term_curvature surface_term_gradient1 surface_term_gradient2' - [../] + [] # [./dofmap] # type = DOFMap # execute_on = 'timestep_end' @@ -415,101 +415,100 @@ pooldepth=${width} show_var_residual_norms = true [] - [Adaptivity] marker = combo max_h_level = 5 initial_steps = 0 - [./Indicators] - [./error_x] + [Indicators] + [error_x] type = GradientJumpIndicator variable = vel_x - [../] - [./error_y] + [] + [error_y] type = GradientJumpIndicator variable = vel_y - [../] + [] # [./error_p] # type = GradientJumpIndicator # variable = p # [../] - [./error_T] + [error_T] type = GradientJumpIndicator variable = T - [../] - [./error_dispx] + [] + [error_dispx] type = GradientJumpIndicator variable = disp_x - [../] - [./error_dispy] + [] + [error_dispy] type = GradientJumpIndicator variable = disp_y - [../] - [../] + [] + [] - [./Markers] - [./errorfrac_x] + [Markers] + [errorfrac_x] type = ErrorFractionMarker refine = 0.8 coarsen = 0.05 indicator = error_x - [../] - [./errorfrac_y] + [] + [errorfrac_y] type = ErrorFractionMarker refine = 0.8 coarsen = 0.05 indicator = error_y - [../] + [] # [./errorfrac_p] # type = ErrorFractionMarker # refine = 0.8 # coarsen = 0.05 # indicator = error_p # [../] - [./errorfrac_T] + [errorfrac_T] type = ErrorFractionMarker refine = 0.8 coarsen = 0.05 indicator = error_T - [../] - [./errorfrac_dispx] + [] + [errorfrac_dispx] type = ErrorFractionMarker refine = 0.8 coarsen = 0.05 indicator = error_dispx - [../] - [./errorfrac_dispy] + [] + [errorfrac_dispy] type = ErrorFractionMarker refine = 0.8 coarsen = 0.05 indicator = error_dispy - [../] - [./combo] + [] + [combo] type = ComboMarker markers = 'errorfrac_x errorfrac_y errorfrac_T errorfrac_dispx errorfrac_dispy' # markers = 'errorfrac_x' - [../] - [../] + [] + [] [] [Postprocessors] - [./num_dofs] + [num_dofs] type = NumDOFs system = 'NL' - [../] - [./nl] + [] + [nl] type = NumNonlinearIterations - [../] - [./tot_nl] + [] + [tot_nl] type = CumulativeValuePostprocessor postprocessor = 'nl' - [../] - [./linear] + [] + [linear] type = NumLinearIterations - [../] - [./tot_linear] + [] + [tot_linear] type = CumulativeValuePostprocessor postprocessor = 'linear' - [../] + [] [] diff --git a/examples/laser_welding_problems/physics/no-surface-tension/cold-no-surface-tension.i b/examples/laser_welding_problems/physics/no-surface-tension/cold-no-surface-tension.i index 20a4a2b..fc92d99 100644 --- a/examples/laser_welding_problems/physics/no-surface-tension/cold-no-surface-tension.i +++ b/examples/laser_welding_problems/physics/no-surface-tension/cold-no-surface-tension.i @@ -1,7 +1,7 @@ -period=1.25e-3 -endtime=2.5e-3 -timestep=1.25e-5 -surfacetemp=300 +period = 1.25e-3 +endtime = 2.5e-3 +timestep = 1.25e-5 +surfacetemp = 300 [GlobalParams] gravity = '0 0 0' @@ -31,91 +31,91 @@ surfacetemp=300 [] [Variables] - [./vel_x] - [./InitialCondition] + [vel_x] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] + [] + [] - [./vel_y] - [./InitialCondition] + [vel_y] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] + [] + [] - [./vel_z] - [./InitialCondition] + [vel_z] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] - - [./T] - [../] - - [./p] - [../] - [./disp_x] - [../] - [./disp_y] - [../] - [./disp_z] - [../] + [] + [] + + [T] + [] + + [p] + [] + [disp_x] + [] + [disp_y] + [] + [disp_z] + [] [] [ICs] - [./T] + [T] type = FunctionIC variable = T function = '(${surfacetemp} - 300) / .7e-3 * z + ${surfacetemp}' - [../] + [] [] [Kernels] - [./disp_x] + [disp_x] type = Diffusion variable = disp_x - [../] - [./disp_y] + [] + [disp_y] type = Diffusion variable = disp_y - [../] - [./disp_z] + [] + [disp_z] type = Diffusion variable = disp_z - [../] + [] [] [ADKernels] - [./mesh_x] + [mesh_x] type = INSConvectedMesh variable = vel_x disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] - [./mesh_y] + [] + [mesh_y] type = INSConvectedMesh variable = vel_y disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] - [./mesh_z] + [] + [mesh_z] type = INSConvectedMesh variable = vel_z disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] + [] # mass - [./mass] + [mass] type = INSADMass variable = p u = vel_x @@ -123,17 +123,17 @@ surfacetemp=300 w = vel_z p = p use_displaced_mesh = true - [../] + [] # x-momentum, time - [./x_momentum_time] + [x_momentum_time] type = INSADMomentumTimeDerivative variable = vel_x use_displaced_mesh = true - [../] + [] # x-momentum, space - [./x_momentum_space] + [x_momentum_space] type = INSADMomentumBase variable = vel_x u = vel_x @@ -142,17 +142,17 @@ surfacetemp=300 p = p component = 0 use_displaced_mesh = true - [../] + [] # y-momentum, time - [./y_momentum_time] + [y_momentum_time] type = INSADMomentumTimeDerivative variable = vel_y use_displaced_mesh = true - [../] + [] # y-momentum, space - [./y_momentum_space] + [y_momentum_space] type = INSADMomentumBase variable = vel_y u = vel_x @@ -161,17 +161,17 @@ surfacetemp=300 p = p component = 1 use_displaced_mesh = true - [../] + [] # z-momentum, time - [./z_momentum_time] + [z_momentum_time] type = INSADMomentumTimeDerivative variable = vel_z use_displaced_mesh = true - [../] + [] # z-momentum, space - [./z_momentum_space] + [z_momentum_space] type = INSADMomentumBase variable = vel_z u = vel_x @@ -180,85 +180,85 @@ surfacetemp=300 p = p component = 2 use_displaced_mesh = true - [../] - - # temperature - [./temperature_time] - type = INSADTemperatureTimeDerivative - variable = T - use_displaced_mesh = true - [../] - - [./temperature_space] - type = INSADTemperature - variable = T - u = vel_x - v = vel_y - w = vel_z - p = p - use_displaced_mesh = true - [../] - [mesh_T] - type = INSTemperatureConvectedMesh - variable = T - disp_x = disp_x - disp_y = disp_y - disp_z = disp_z - use_displaced_mesh = true - [] + [] + + # temperature + [temperature_time] + type = INSADTemperatureTimeDerivative + variable = T + use_displaced_mesh = true + [] + + [temperature_space] + type = INSADTemperature + variable = T + u = vel_x + v = vel_y + w = vel_z + p = p + use_displaced_mesh = true + [] + [mesh_T] + type = INSTemperatureConvectedMesh + variable = T + disp_x = disp_x + disp_y = disp_y + disp_z = disp_z + use_displaced_mesh = true + [] [] [BCs] - [./x_no_disp] + [x_no_disp] type = DirichletBC variable = disp_x boundary = 'back' value = 0 - [../] - [./y_no_disp] + [] + [y_no_disp] type = DirichletBC variable = disp_y boundary = 'back' value = 0 - [../] - [./z_no_disp] + [] + [z_no_disp] type = DirichletBC variable = disp_z boundary = 'back' value = 0 - [../] + [] - [./x_no_slip] + [x_no_slip] type = DirichletBC variable = vel_x boundary = 'bottom right left top back' value = 0.0 - [../] + [] - [./y_no_slip] + [y_no_slip] type = DirichletBC variable = vel_y boundary = 'bottom right left top back' value = 0.0 - [../] + [] - [./z_no_slip] + [z_no_slip] type = DirichletBC variable = vel_z boundary = 'bottom right left top back' value = 0.0 - [../] + [] - [./T_cold] + [T_cold] type = DirichletBC variable = T boundary = 'back' value = 300 - [../] + [] [] [ADBCs] - [./radiation_flux] + [radiation_flux] type = RadiationEnergyFluxBC variable = T boundary = 'front' @@ -266,9 +266,9 @@ surfacetemp=300 sb_constant = 'sb_constant' absorptivity = 'abs' use_displaced_mesh = true - [../] + [] - [./weld_flux] + [weld_flux] type = GaussianWeldEnergyFluxBC variable = T boundary = 'front' @@ -279,79 +279,79 @@ surfacetemp=300 y_beam_coord = '2e-4 * sin(t * 2 * pi / ${period})' z_beam_coord = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_x] + [vapor_recoil_x] type = VaporRecoilPressureMomentumFluxBC variable = vel_x boundary = 'front' component = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_y] + [vapor_recoil_y] type = VaporRecoilPressureMomentumFluxBC variable = vel_y boundary = 'front' component = 1 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_z] + [vapor_recoil_z] type = VaporRecoilPressureMomentumFluxBC variable = vel_z boundary = 'front' component = 2 use_displaced_mesh = true - [../] + [] -[./displace_x_top] + [displace_x_top] type = DisplaceBoundaryBC boundary = 'front' variable = 'disp_x' velocity = 'vel_x' - [../] - [./displace_y_top] + [] + [displace_y_top] type = DisplaceBoundaryBC boundary = 'front' variable = 'disp_y' velocity = 'vel_y' - [../] - [./displace_z_top] + [] + [displace_z_top] type = DisplaceBoundaryBC boundary = 'front' variable = 'disp_z' velocity = 'vel_z' - [../] + [] [] [ADMaterials] - [./kc_fits] + [kc_fits] type = DemonstrationPlantFits temperature = T beta = 1e7 - [../] - [./boundary] + [] + [boundary] type = DemonstrationPlantFitsBoundary boundary = 'front' temperature = T - [../] + [] [] [Materials] - [./const] + [const] type = GenericConstantMaterial prop_names = 'abs sb_constant' prop_values = '1 5.67e-8' - [../] + [] [] [Preconditioning] - [./SMP] + [SMP] type = SMP full = true solve_type = 'NEWTON' - [../] + [] [] [Executioner] @@ -366,26 +366,26 @@ surfacetemp=300 line_search = 'none' nl_max_its = 12 l_max_its = 100 - [./TimeStepper] + [TimeStepper] type = IterationAdaptiveDT optimal_iterations = 7 dt = ${timestep} linear_iteration_ratio = 1e6 growth_factor = 1.5 - [../] + [] [] [Outputs] print_linear_residuals = false - [./exodus] + [exodus] type = Exodus output_material_properties = true show_material_properties = 'mu' - [../] - [./dofmap] + [] + [dofmap] type = DOFMap execute_on = 'initial' - [../] + [] checkpoint = true [] @@ -393,12 +393,11 @@ surfacetemp=300 show_var_residual_norms = true [] - [Adaptivity] marker = combo max_h_level = 4 - [./Indicators] + [Indicators] # [./error_x] # type = GradientJumpIndicator # variable = vel_x @@ -415,10 +414,10 @@ surfacetemp=300 # type = GradientJumpIndicator # variable = p # [../] - [./error_T] + [error_T] type = GradientJumpIndicator variable = T - [../] + [] # [./error_dispx] # type = GradientJumpIndicator # variable = disp_x @@ -427,13 +426,13 @@ surfacetemp=300 # type = GradientJumpIndicator # variable = disp_y # [../] - [./error_dispz] + [error_dispz] type = GradientJumpIndicator variable = disp_z - [../] - [../] + [] + [] - [./Markers] + [Markers] # [./errorfrac_x] # type = ErrorFractionMarker # refine = 0.7 @@ -458,12 +457,12 @@ surfacetemp=300 # coarsen = 0.3 # indicator = error_p # [../] - [./errorfrac_T] + [errorfrac_T] type = ErrorFractionMarker refine = 0.4 coarsen = 0.2 indicator = error_T - [../] + [] # [./errorfrac_dispx] # type = ErrorFractionMarker # refine = 0.7 @@ -476,36 +475,36 @@ surfacetemp=300 # coarsen = 0.3 # indicator = error_dispy # [../] - [./errorfrac_dispz] + [errorfrac_dispz] type = ErrorFractionMarker refine = 0.4 coarsen = 0.2 indicator = error_dispz - [../] - [./combo] + [] + [combo] type = ComboMarker markers = 'errorfrac_T errorfrac_dispz' - [../] - [../] + [] + [] [] [Postprocessors] - [./num_dofs] + [num_dofs] type = NumDOFs system = 'NL' - [../] - [./nl] + [] + [nl] type = NumNonlinearIterations - [../] - [./tot_nl] + [] + [tot_nl] type = CumulativeValuePostprocessor postprocessor = 'nl' - [../] - [./linear] + [] + [linear] type = NumLinearIterations - [../] - [./tot_linear] + [] + [tot_linear] type = CumulativeValuePostprocessor postprocessor = 'linear' - [../] + [] [] diff --git a/examples/laser_welding_problems/physics/no-surface-tension/hot-no-surface-tension.i b/examples/laser_welding_problems/physics/no-surface-tension/hot-no-surface-tension.i index 2463c53..65f9292 100644 --- a/examples/laser_welding_problems/physics/no-surface-tension/hot-no-surface-tension.i +++ b/examples/laser_welding_problems/physics/no-surface-tension/hot-no-surface-tension.i @@ -1,5 +1,5 @@ -period=4e-5 -timestep=4e-7 +period = 4e-5 +timestep = 4e-7 [GlobalParams] gravity = '0 0 0' @@ -29,91 +29,91 @@ timestep=4e-7 [] [Variables] - [./vel_x] - [./InitialCondition] + [vel_x] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] + [] + [] - [./vel_y] - [./InitialCondition] + [vel_y] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] + [] + [] - [./vel_z] - [./InitialCondition] + [vel_z] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] - - [./T] - [../] - - [./p] - [../] - [./disp_x] - [../] - [./disp_y] - [../] - [./disp_z] - [../] + [] + [] + + [T] + [] + + [p] + [] + [disp_x] + [] + [disp_y] + [] + [disp_z] + [] [] [ICs] - [./T] + [T] type = FunctionIC variable = T function = '(3000 - 300) / .7e-3 * z + 3000' - [../] + [] [] [Kernels] - [./disp_x] + [disp_x] type = Diffusion variable = disp_x - [../] - [./disp_y] + [] + [disp_y] type = Diffusion variable = disp_y - [../] - [./disp_z] + [] + [disp_z] type = Diffusion variable = disp_z - [../] + [] [] [ADKernels] - [./mesh_x] + [mesh_x] type = INSConvectedMesh variable = vel_x disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] - [./mesh_y] + [] + [mesh_y] type = INSConvectedMesh variable = vel_y disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] - [./mesh_z] + [] + [mesh_z] type = INSConvectedMesh variable = vel_z disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] + [] -# mass - [./mass] + # mass + [mass] type = INSADMass variable = p u = vel_x @@ -121,17 +121,17 @@ timestep=4e-7 w = vel_z p = p use_displaced_mesh = true - [../] + [] # x-momentum, time - [./x_momentum_time] + [x_momentum_time] type = INSADMomentumTimeDerivative variable = vel_x use_displaced_mesh = true - [../] + [] # x-momentum, space - [./x_momentum_space] + [x_momentum_space] type = INSADMomentumBase variable = vel_x u = vel_x @@ -140,17 +140,17 @@ timestep=4e-7 p = p component = 0 use_displaced_mesh = true - [../] + [] # y-momentum, time - [./y_momentum_time] + [y_momentum_time] type = INSADMomentumTimeDerivative variable = vel_y use_displaced_mesh = true - [../] + [] # y-momentum, space - [./y_momentum_space] + [y_momentum_space] type = INSADMomentumBase variable = vel_y u = vel_x @@ -159,17 +159,17 @@ timestep=4e-7 p = p component = 1 use_displaced_mesh = true - [../] + [] # z-momentum, time - [./z_momentum_time] + [z_momentum_time] type = INSADMomentumTimeDerivative variable = vel_z use_displaced_mesh = true - [../] + [] # z-momentum, space - [./z_momentum_space] + [z_momentum_space] type = INSADMomentumBase variable = vel_z u = vel_x @@ -178,85 +178,85 @@ timestep=4e-7 p = p component = 2 use_displaced_mesh = true - [../] - - # temperature - [./temperature_time] - type = INSADTemperatureTimeDerivative - variable = T - use_displaced_mesh = true - [../] - - [./temperature_space] - type = INSADTemperature - variable = T - u = vel_x - v = vel_y - w = vel_z - p = p - use_displaced_mesh = true - [../] - [mesh_T] - type = INSTemperatureConvectedMesh - variable = T - disp_x = disp_x - disp_y = disp_y - disp_z = disp_z - use_displaced_mesh = true - [] + [] + + # temperature + [temperature_time] + type = INSADTemperatureTimeDerivative + variable = T + use_displaced_mesh = true + [] + + [temperature_space] + type = INSADTemperature + variable = T + u = vel_x + v = vel_y + w = vel_z + p = p + use_displaced_mesh = true + [] + [mesh_T] + type = INSTemperatureConvectedMesh + variable = T + disp_x = disp_x + disp_y = disp_y + disp_z = disp_z + use_displaced_mesh = true + [] [] [BCs] - [./x_no_disp] + [x_no_disp] type = DirichletBC variable = disp_x boundary = 'back' value = 0 - [../] - [./y_no_disp] + [] + [y_no_disp] type = DirichletBC variable = disp_y boundary = 'back' value = 0 - [../] - [./z_no_disp] + [] + [z_no_disp] type = DirichletBC variable = disp_z boundary = 'back' value = 0 - [../] + [] - [./x_no_slip] + [x_no_slip] type = DirichletBC variable = vel_x boundary = 'bottom right left top back' value = 0.0 - [../] + [] - [./y_no_slip] + [y_no_slip] type = DirichletBC variable = vel_y boundary = 'bottom right left top back' value = 0.0 - [../] + [] - [./z_no_slip] + [z_no_slip] type = DirichletBC variable = vel_z boundary = 'bottom right left top back' value = 0.0 - [../] + [] - [./T_cold] + [T_cold] type = DirichletBC variable = T boundary = 'back' value = 300 - [../] + [] [] [ADBCs] - [./radiation_flux] + [radiation_flux] type = RadiationEnergyFluxBC variable = T boundary = 'front' @@ -264,9 +264,9 @@ timestep=4e-7 sb_constant = 'sb_constant' absorptivity = 'abs' use_displaced_mesh = true - [../] + [] - [./weld_flux] + [weld_flux] type = GaussianWeldEnergyFluxBC variable = T boundary = 'front' @@ -277,79 +277,79 @@ timestep=4e-7 y_beam_coord = '2e-4 * sin(t * 2 * pi / ${period})' z_beam_coord = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_x] + [vapor_recoil_x] type = VaporRecoilPressureMomentumFluxBC variable = vel_x boundary = 'front' component = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_y] + [vapor_recoil_y] type = VaporRecoilPressureMomentumFluxBC variable = vel_y boundary = 'front' component = 1 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_z] + [vapor_recoil_z] type = VaporRecoilPressureMomentumFluxBC variable = vel_z boundary = 'front' component = 2 use_displaced_mesh = true - [../] + [] - [./displace_x_top] + [displace_x_top] type = DisplaceBoundaryBC boundary = 'front' variable = 'disp_x' velocity = 'vel_x' - [../] - [./displace_y_top] + [] + [displace_y_top] type = DisplaceBoundaryBC boundary = 'front' variable = 'disp_y' velocity = 'vel_y' - [../] - [./displace_z_top] + [] + [displace_z_top] type = DisplaceBoundaryBC boundary = 'front' variable = 'disp_z' velocity = 'vel_z' - [../] + [] [] [ADMaterials] - [./kc_fits] + [kc_fits] type = DemonstrationPlantFits temperature = T beta = 1e7 - [../] - [./boundary] + [] + [boundary] type = DemonstrationPlantFitsBoundary boundary = 'front' temperature = T - [../] + [] [] [Materials] - [./const] + [const] type = GenericConstantMaterial prop_names = 'abs sb_constant' prop_values = '1 5.67e-8' - [../] + [] [] [Preconditioning] - [./SMP] + [SMP] type = SMP full = true solve_type = 'NEWTON' - [../] + [] [] [Executioner] @@ -364,26 +364,26 @@ timestep=4e-7 line_search = 'none' nl_max_its = 12 l_max_its = 100 - [./TimeStepper] + [TimeStepper] type = IterationAdaptiveDT optimal_iterations = 7 dt = ${timestep} linear_iteration_ratio = 1e6 growth_factor = 1.5 - [../] + [] [] [Outputs] print_linear_residuals = false - [./exodus] + [exodus] type = Exodus output_material_properties = true show_material_properties = 'mu' - [../] - [./dofmap] + [] + [dofmap] type = DOFMap execute_on = 'initial' - [../] + [] checkpoint = true [] @@ -391,12 +391,11 @@ timestep=4e-7 show_var_residual_norms = true [] - [Adaptivity] marker = combo max_h_level = 3 - [./Indicators] + [Indicators] # [./error_x] # type = GradientJumpIndicator # variable = vel_x @@ -413,10 +412,10 @@ timestep=4e-7 # type = GradientJumpIndicator # variable = p # [../] - [./error_T] + [error_T] type = GradientJumpIndicator variable = T - [../] + [] # [./error_dispx] # type = GradientJumpIndicator # variable = disp_x @@ -425,13 +424,13 @@ timestep=4e-7 # type = GradientJumpIndicator # variable = disp_y # [../] - [./error_dispz] + [error_dispz] type = GradientJumpIndicator variable = disp_z - [../] - [../] + [] + [] - [./Markers] + [Markers] # [./errorfrac_x] # type = ErrorFractionMarker # refine = 0.7 @@ -456,12 +455,12 @@ timestep=4e-7 # coarsen = 0.3 # indicator = error_p # [../] - [./errorfrac_T] + [errorfrac_T] type = ErrorFractionMarker refine = 0.5 coarsen = 0.2 indicator = error_T - [../] + [] # [./errorfrac_dispx] # type = ErrorFractionMarker # refine = 0.7 @@ -474,36 +473,36 @@ timestep=4e-7 # coarsen = 0.3 # indicator = error_dispy # [../] - [./errorfrac_dispz] + [errorfrac_dispz] type = ErrorFractionMarker refine = 0.5 coarsen = 0.2 indicator = error_dispz - [../] - [./combo] + [] + [combo] type = ComboMarker markers = 'errorfrac_T errorfrac_dispz' - [../] - [../] + [] + [] [] [Postprocessors] - [./num_dofs] + [num_dofs] type = NumDOFs system = 'NL' - [../] - [./nl] + [] + [nl] type = NumNonlinearIterations - [../] - [./tot_nl] + [] + [tot_nl] type = CumulativeValuePostprocessor postprocessor = 'nl' - [../] - [./linear] + [] + [linear] type = NumLinearIterations - [../] - [./tot_linear] + [] + [tot_linear] type = CumulativeValuePostprocessor postprocessor = 'linear' - [../] + [] [] diff --git a/examples/laser_welding_problems/physics/surface-tension/2d-one-spot-scaled-2nd-order.i b/examples/laser_welding_problems/physics/surface-tension/2d-one-spot-scaled-2nd-order.i index c509028..1e8ecdd 100644 --- a/examples/laser_welding_problems/physics/surface-tension/2d-one-spot-scaled-2nd-order.i +++ b/examples/laser_welding_problems/physics/surface-tension/2d-one-spot-scaled-2nd-order.i @@ -1,14 +1,14 @@ -length_unit_exponent=-4 -temperature_unit_exponent=3 -mass_unit_exponent=-6 -time_unit_exponent=-5 - -endtime=${fparse 500e-5 / 10^time_unit_exponent} -timestep=${fparse 5e-7 / 10^time_unit_exponent} -surfacetemp=${fparse 300 / 10^temperature_unit_exponent} -bottomtemp=${fparse 300 / 10^temperature_unit_exponent} -pooldepth=${fparse 2e-4 / 10^length_unit_exponent} -half_width=${fparse 4e-4 / 10^length_unit_exponent} +length_unit_exponent = -4 +temperature_unit_exponent = 3 +mass_unit_exponent = -6 +time_unit_exponent = -5 + +endtime = '${fparse 500e-5 / 10^time_unit_exponent}' +timestep = '${fparse 5e-7 / 10^time_unit_exponent}' +surfacetemp = '${fparse 300 / 10^temperature_unit_exponent}' +bottomtemp = '${fparse 300 / 10^temperature_unit_exponent}' +pooldepth = '${fparse 2e-4 / 10^length_unit_exponent}' +half_width = '${fparse 4e-4 / 10^length_unit_exponent}' [GlobalParams] gravity = '0 0 0' @@ -25,9 +25,9 @@ half_width=${fparse 4e-4 / 10^length_unit_exponent} [Mesh] type = GeneratedMesh dim = 2 - xmin = ${fparse -half_width} + xmin = '${fparse -half_width}' xmax = ${half_width} - ymin = ${fparse -pooldepth} + ymin = '${fparse -pooldepth}' ymax = 0 nx = 4 ny = 1 @@ -41,100 +41,100 @@ half_width=${fparse 4e-4 / 10^length_unit_exponent} [] [Variables] - [./vel_x] - [./InitialCondition] + [vel_x] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] + [] scaling = 1e6 - [../] + [] - [./vel_y] - [./InitialCondition] + [vel_y] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] + [] scaling = 1e6 - [../] + [] - [./T] - [../] + [T] + [] - [./p] + [p] order = FIRST - [../] - [./disp_x] + [] + [disp_x] scaling = 1e-3 - [../] - [./disp_y] + [] + [disp_y] scaling = 1e-3 - [../] + [] [] [ICs] - [./T] + [T] type = FunctionIC variable = T function = '(${surfacetemp} - ${bottomtemp}) / ${pooldepth} * z + ${surfacetemp}' - [../] + [] [] [Kernels] [] [ADKernels] - [./disp_x] + [disp_x] type = ADStressDivergence variable = disp_x component = 0 - [../] - [./disp_y] + [] + [disp_y] type = ADStressDivergence variable = disp_y component = 1 - [../] + [] -[./mesh_x] + [mesh_x] type = INSConvectedMesh variable = vel_x disp_x = disp_x disp_y = disp_y use_displaced_mesh = true - [../] - [./mesh_y] + [] + [mesh_y] type = INSConvectedMesh variable = vel_y disp_x = disp_x disp_y = disp_y use_displaced_mesh = true - [../] - [./mesh_T] + [] + [mesh_T] type = INSTemperatureConvectedMesh variable = T disp_x = disp_x disp_y = disp_y use_displaced_mesh = true - [../] + [] # mass - [./mass] + [mass] type = INSADMass variable = p u = vel_x v = vel_y p = p use_displaced_mesh = true - [../] + [] # x-momentum, time - [./x_momentum_time] + [x_momentum_time] type = INSADMomentumTimeDerivative variable = vel_x use_displaced_mesh = true - [../] + [] # x-momentum, space - [./x_momentum_space] + [x_momentum_space] type = INSADMomentumBase variable = vel_x u = vel_x @@ -142,17 +142,17 @@ half_width=${fparse 4e-4 / 10^length_unit_exponent} p = p component = 0 use_displaced_mesh = true - [../] + [] # y-momentum, time - [./y_momentum_time] + [y_momentum_time] type = INSADMomentumTimeDerivative variable = vel_y use_displaced_mesh = true - [../] + [] # y-momentum, space - [./y_momentum_space] + [y_momentum_space] type = INSADMomentumBase variable = vel_y u = vel_x @@ -160,63 +160,63 @@ half_width=${fparse 4e-4 / 10^length_unit_exponent} p = p component = 1 use_displaced_mesh = true - [../] - - # temperature - [./temperature_time] - type = INSADTemperatureTimeDerivative - variable = T - use_displaced_mesh = true - [../] - - [./temperature_space] - type = INSADTemperature - variable = T - u = vel_x - v = vel_y - p = p - use_displaced_mesh = true - [../] + [] + + # temperature + [temperature_time] + type = INSADTemperatureTimeDerivative + variable = T + use_displaced_mesh = true + [] + + [temperature_space] + type = INSADTemperature + variable = T + u = vel_x + v = vel_y + p = p + use_displaced_mesh = true + [] [] [BCs] - [./x_no_disp] + [x_no_disp] type = DirichletBC variable = disp_x boundary = 'bottom left right' value = 0 - [../] - [./y_no_disp] + [] + [y_no_disp] type = DirichletBC variable = disp_y boundary = 'bottom left right' value = 0 - [../] + [] - [./x_no_slip] + [x_no_slip] type = DirichletBC variable = vel_x boundary = 'bottom left right' value = 0.0 - [../] + [] - [./y_no_slip] + [y_no_slip] type = DirichletBC variable = vel_y boundary = 'bottom left right' value = 0.0 - [../] + [] - [./T_cold] + [T_cold] type = DirichletBC variable = T boundary = 'bottom' value = ${bottomtemp} - [../] + [] [] [ADBCs] - [./radiation_flux] + [radiation_flux] type = RadiationEnergyFluxBC variable = T boundary = 'top' @@ -224,54 +224,54 @@ half_width=${fparse 4e-4 / 10^length_unit_exponent} sb_constant = 'sb_constant' absorptivity = 'abs' use_displaced_mesh = true - [../] + [] - [./weld_flux] + [weld_flux] type = GaussianWeldEnergyFluxBC variable = T boundary = 'top' reff = 0.6 - F0 = ${fparse 2.546e9 / 10^mass_unit_exponent * 10^(3 * time_unit_exponent)} - R = ${fparse 1e-4 / 10^length_unit_exponent} + F0 = '${fparse 2.546e9 / 10^mass_unit_exponent * 10^(3 * time_unit_exponent)}' + R = '${fparse 1e-4 / 10^length_unit_exponent}' x_beam_coord = 0 y_beam_coord = 0 z_beam_coord = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_x] + [vapor_recoil_x] type = VaporRecoilPressureMomentumFluxBC variable = vel_x boundary = 'top' component = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_y] + [vapor_recoil_y] type = VaporRecoilPressureMomentumFluxBC variable = vel_y boundary = 'top' component = 1 use_displaced_mesh = true - [../] + [] - [./surface_x] + [surface_x] type = SurfaceTensionBC variable = vel_x boundary = 'top' component = 0 use_displaced_mesh = true - [../] + [] - [./surface_y] + [surface_y] type = SurfaceTensionBC variable = vel_y boundary = 'top' component = 1 use_displaced_mesh = true - [../] + [] - [./displace_x_top] + [displace_x_top] type = PenaltyDisplaceBoundaryBC boundary = 'top' variable = 'disp_x' @@ -279,8 +279,8 @@ half_width=${fparse 4e-4 / 10^length_unit_exponent} disp_x = 'disp_x' vel_y = 'vel_y' disp_y = 'disp_y' - [../] - [./displace_y_top] + [] + [displace_y_top] type = PenaltyDisplaceBoundaryBC boundary = 'top' variable = 'disp_y' @@ -288,11 +288,11 @@ half_width=${fparse 4e-4 / 10^length_unit_exponent} disp_x = 'disp_x' vel_y = 'vel_y' disp_y = 'disp_y' - [../] + [] [] [ADMaterials] - [./kc_fits] + [kc_fits] type = DemonstrationPlantFits temperature = T beta = 1e7 @@ -302,8 +302,8 @@ half_width=${fparse 4e-4 / 10^length_unit_exponent} time_unit_exponent = ${time_unit_exponent} Tl = 500 T90 = 400 - [../] - [./boundary] + [] + [boundary] type = DemonstrationPlantFitsBoundary boundary = 'top' temperature = T @@ -312,34 +312,34 @@ half_width=${fparse 4e-4 / 10^length_unit_exponent} temperature_unit_exponent = ${temperature_unit_exponent} mass_unit_exponent = ${mass_unit_exponent} time_unit_exponent = ${time_unit_exponent} - [../] - [./stress] + [] + [stress] type = PseudoSolidStress disp_x = disp_x disp_y = disp_y - [../] + [] [] [Materials] - [./const] + [const] type = GenericConstantMaterial prop_names = 'abs sb_constant' prop_values = '1 ${fparse 5.67e-8 / 10^mass_unit_exponent * 10^(4 * temperature_unit_exponent) * 10^(3 * time_unit_exponent)}' - [../] + [] [] [Preconditioning] - [./SMP] + [SMP] type = SMP full = true solve_type = 'NEWTON' - [../] + [] [] [Executioner] type = Transient end_time = ${endtime} - dtmin = ${fparse 1e-8 / 10^time_unit_exponent} + dtmin = '${fparse 1e-8 / 10^time_unit_exponent}' num_steps = 300 petsc_options = '-snes_converged_reason -ksp_converged_reason -options_left -ksp_monitor_singular_value' petsc_options_iname = '-ksp_max_it -ksp_gmres_restart -pc_type -snes_max_funcs -sub_pc_factor_levels' @@ -348,23 +348,23 @@ half_width=${fparse 4e-4 / 10^length_unit_exponent} line_search = 'none' nl_max_its = 12 l_tol = 1e-3 - [./TimeStepper] + [TimeStepper] type = IterationAdaptiveDT optimal_iterations = 8 dt = ${timestep} linear_iteration_ratio = 1e6 growth_factor = 1.5 - [../] + [] [] [Outputs] print_linear_residuals = false csv = true - [./exodus] + [exodus] type = Exodus output_material_properties = true show_material_properties = 'mu surface_term_curvature surface_term_gradient1 surface_term_gradient2' - [../] + [] # [./dofmap] # type = DOFMap # execute_on = 'timestep_end' @@ -376,7 +376,6 @@ half_width=${fparse 4e-4 / 10^length_unit_exponent} show_var_residual_norms = true [] - # [Adaptivity] # marker = combo # max_h_level = 6 @@ -455,22 +454,22 @@ half_width=${fparse 4e-4 / 10^length_unit_exponent} # [] [Postprocessors] - [./num_dofs] + [num_dofs] type = NumDOFs system = 'NL' - [../] - [./nl] + [] + [nl] type = NumNonlinearIterations - [../] - [./tot_nl] + [] + [tot_nl] type = CumulativeValuePostprocessor postprocessor = 'nl' - [../] - [./linear] + [] + [linear] type = NumLinearIterations - [../] - [./tot_linear] + [] + [tot_linear] type = CumulativeValuePostprocessor postprocessor = 'linear' - [../] + [] [] diff --git a/examples/laser_welding_problems/physics/surface-tension/2d-one-spot-unscaled-2nd-order.i b/examples/laser_welding_problems/physics/surface-tension/2d-one-spot-unscaled-2nd-order.i index 1e79a8b..3655450 100644 --- a/examples/laser_welding_problems/physics/surface-tension/2d-one-spot-unscaled-2nd-order.i +++ b/examples/laser_welding_problems/physics/surface-tension/2d-one-spot-unscaled-2nd-order.i @@ -1,8 +1,8 @@ -endtime=500e-5 -timestep=.5e-5 -surfacetemp=300 -bottomtemp=300 -pooldepth=2e-4 +endtime = 500e-5 +timestep = .5e-5 +surfacetemp = 300 +bottomtemp = 300 +pooldepth = 2e-4 [GlobalParams] gravity = '0 0 0' @@ -35,93 +35,93 @@ pooldepth=2e-4 [] [Variables] - [./vel_x] - [./InitialCondition] + [vel_x] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] + [] + [] - [./vel_y] - [./InitialCondition] + [vel_y] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] + [] + [] - [./T] - [../] + [T] + [] - [./p] + [p] order = FIRST - [../] - [./disp_x] - [../] - [./disp_y] - [../] + [] + [disp_x] + [] + [disp_y] + [] [] [ICs] - [./T] + [T] type = FunctionIC variable = T function = '(${surfacetemp} - ${bottomtemp}) / ${pooldepth} * z + ${surfacetemp}' - [../] + [] [] [Kernels] - [./disp_x] + [disp_x] type = Diffusion variable = disp_x - [../] - [./disp_y] + [] + [disp_y] type = Diffusion variable = disp_y - [../] + [] [] [ADKernels] - [./mesh_x] + [mesh_x] type = INSConvectedMesh variable = vel_x disp_x = disp_x disp_y = disp_y use_displaced_mesh = true - [../] - [./mesh_y] + [] + [mesh_y] type = INSConvectedMesh variable = vel_y disp_x = disp_x disp_y = disp_y use_displaced_mesh = true - [../] - [./mesh_T] + [] + [mesh_T] type = INSTemperatureConvectedMesh variable = T disp_x = disp_x disp_y = disp_y use_displaced_mesh = true - [../] + [] # mass - [./mass] + [mass] type = INSADMass variable = p u = vel_x v = vel_y p = p use_displaced_mesh = true - [../] + [] # x-momentum, time - [./x_momentum_time] + [x_momentum_time] type = INSADMomentumTimeDerivative variable = vel_x use_displaced_mesh = true - [../] + [] # x-momentum, space - [./x_momentum_space] + [x_momentum_space] type = INSADMomentumBase variable = vel_x u = vel_x @@ -129,17 +129,17 @@ pooldepth=2e-4 p = p component = 0 use_displaced_mesh = true - [../] + [] # y-momentum, time - [./y_momentum_time] + [y_momentum_time] type = INSADMomentumTimeDerivative variable = vel_y use_displaced_mesh = true - [../] + [] # y-momentum, space - [./y_momentum_space] + [y_momentum_space] type = INSADMomentumBase variable = vel_y u = vel_x @@ -147,63 +147,63 @@ pooldepth=2e-4 p = p component = 1 use_displaced_mesh = true - [../] - - # temperature - [./temperature_time] - type = INSADTemperatureTimeDerivative - variable = T - use_displaced_mesh = true - [../] - - [./temperature_space] - type = INSADTemperature - variable = T - u = vel_x - v = vel_y - p = p - use_displaced_mesh = true - [../] + [] + + # temperature + [temperature_time] + type = INSADTemperatureTimeDerivative + variable = T + use_displaced_mesh = true + [] + + [temperature_space] + type = INSADTemperature + variable = T + u = vel_x + v = vel_y + p = p + use_displaced_mesh = true + [] [] [BCs] - [./x_no_disp] + [x_no_disp] type = DirichletBC variable = disp_x boundary = 'bottom' value = 0 - [../] - [./y_no_disp] + [] + [y_no_disp] type = DirichletBC variable = disp_y boundary = 'bottom' value = 0 - [../] + [] - [./x_no_slip] + [x_no_slip] type = DirichletBC variable = vel_x boundary = 'bottom right left' value = 0.0 - [../] + [] - [./y_no_slip] + [y_no_slip] type = DirichletBC variable = vel_y boundary = 'bottom right left' value = 0.0 - [../] + [] - [./T_cold] + [T_cold] type = DirichletBC variable = T boundary = 'bottom' value = ${bottomtemp} - [../] + [] [] [ADBCs] - [./radiation_flux] + [radiation_flux] type = RadiationEnergyFluxBC variable = T boundary = 'top' @@ -211,9 +211,9 @@ pooldepth=2e-4 sb_constant = 'sb_constant' absorptivity = 'abs' use_displaced_mesh = true - [../] + [] - [./weld_flux] + [weld_flux] type = GaussianWeldEnergyFluxBC variable = T boundary = 'top' @@ -224,82 +224,82 @@ pooldepth=2e-4 y_beam_coord = 0 z_beam_coord = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_x] + [vapor_recoil_x] type = VaporRecoilPressureMomentumFluxBC variable = vel_x boundary = 'top' component = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_y] + [vapor_recoil_y] type = VaporRecoilPressureMomentumFluxBC variable = vel_y boundary = 'top' component = 1 use_displaced_mesh = true - [../] + [] - [./surface_x] + [surface_x] type = SurfaceTensionBC variable = vel_x boundary = 'top' component = 0 use_displaced_mesh = true - [../] + [] - [./surface_y] + [surface_y] type = SurfaceTensionBC variable = vel_y boundary = 'top' component = 1 use_displaced_mesh = true - [../] + [] - [./displace_x_top] + [displace_x_top] type = DisplaceBoundaryBC boundary = 'top' variable = 'disp_x' velocity = 'vel_x' - [../] - [./displace_y_top] + [] + [displace_y_top] type = DisplaceBoundaryBC boundary = 'top' variable = 'disp_y' velocity = 'vel_y' - [../] + [] [] [ADMaterials] - [./kc_fits] + [kc_fits] type = DemonstrationPlantFits temperature = T beta = 1e7 - [../] - [./boundary] + [] + [boundary] type = DemonstrationPlantFitsBoundary boundary = 'top' temperature = T use_displaced_mesh = true - [../] + [] [] [Materials] - [./const] + [const] type = GenericConstantMaterial prop_names = 'abs sb_constant' prop_values = '1 5.67e-8' - [../] + [] [] [Preconditioning] - [./SMP] + [SMP] type = SMP full = true solve_type = 'NEWTON' - [../] + [] [] [Executioner] @@ -314,26 +314,26 @@ pooldepth=2e-4 line_search = 'none' nl_max_its = 12 l_tol = 1e-3 - [./TimeStepper] + [TimeStepper] type = IterationAdaptiveDT optimal_iterations = 7 dt = ${timestep} linear_iteration_ratio = 1e6 growth_factor = 1.5 - [../] + [] [] [Outputs] print_linear_residuals = false - [./exodus] + [exodus] type = Exodus output_material_properties = true show_material_properties = 'mu surface_term_curvature surface_term_gradient1 surface_term_gradient2' - [../] - [./dofmap] + [] + [dofmap] type = DOFMap execute_on = 'initial' - [../] + [] checkpoint = true [] @@ -341,100 +341,99 @@ pooldepth=2e-4 show_var_residual_norms = true [] - [Adaptivity] marker = combo max_h_level = 1 initial_steps = 0 - [./Indicators] - [./error_x] + [Indicators] + [error_x] type = GradientJumpIndicator variable = vel_x - [../] - [./error_y] + [] + [error_y] type = GradientJumpIndicator variable = vel_y - [../] + [] # [./error_p] # type = GradientJumpIndicator # variable = p # [../] - [./error_T] + [error_T] type = GradientJumpIndicator variable = T - [../] - [./error_dispx] + [] + [error_dispx] type = GradientJumpIndicator variable = disp_x - [../] - [./error_dispy] + [] + [error_dispy] type = GradientJumpIndicator variable = disp_y - [../] - [../] + [] + [] - [./Markers] - [./errorfrac_x] + [Markers] + [errorfrac_x] type = ErrorFractionMarker refine = 0.4 coarsen = 0.2 indicator = error_x - [../] - [./errorfrac_y] + [] + [errorfrac_y] type = ErrorFractionMarker refine = 0.4 coarsen = 0.2 indicator = error_y - [../] + [] # [./errorfrac_p] # type = ErrorFractionMarker # refine = 0.7 # coarsen = 0.3 # indicator = error_p # [../] - [./errorfrac_T] + [errorfrac_T] type = ErrorFractionMarker refine = 0.4 coarsen = 0.2 indicator = error_T - [../] - [./errorfrac_dispx] + [] + [errorfrac_dispx] type = ErrorFractionMarker refine = 0.4 coarsen = 0.2 indicator = error_dispx - [../] - [./errorfrac_dispy] + [] + [errorfrac_dispy] type = ErrorFractionMarker refine = 0.4 coarsen = 0.2 indicator = error_dispy - [../] - [./combo] + [] + [combo] type = ComboMarker markers = 'errorfrac_x errorfrac_y errorfrac_T errorfrac_dispx errorfrac_dispy' - [../] - [../] + [] + [] [] [Postprocessors] - [./num_dofs] + [num_dofs] type = NumDOFs system = 'NL' - [../] - [./nl] + [] + [nl] type = NumNonlinearIterations - [../] - [./tot_nl] + [] + [tot_nl] type = CumulativeValuePostprocessor postprocessor = 'nl' - [../] - [./linear] + [] + [linear] type = NumLinearIterations - [../] - [./tot_linear] + [] + [tot_linear] type = CumulativeValuePostprocessor postprocessor = 'linear' - [../] + [] [] diff --git a/examples/laser_welding_problems/physics/surface-tension/3d-one-spot-unscaled-1st-order.i b/examples/laser_welding_problems/physics/surface-tension/3d-one-spot-unscaled-1st-order.i index 5884bda..50f5f25 100644 --- a/examples/laser_welding_problems/physics/surface-tension/3d-one-spot-unscaled-1st-order.i +++ b/examples/laser_welding_problems/physics/surface-tension/3d-one-spot-unscaled-1st-order.i @@ -1,8 +1,8 @@ # period=5e-4 -endtime=5e-3 -timestep=5e-6 -surfacetemp=300 -pooldepth=.2e-3 +endtime = 5e-3 +timestep = 5e-6 +surfacetemp = 300 +pooldepth = .2e-3 [GlobalParams] gravity = '0 0 0' @@ -36,98 +36,98 @@ pooldepth=.2e-3 [] [Variables] - [./vel_x] + [vel_x] scaling = 1e6 - [./InitialCondition] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] + [] + [] - [./vel_y] + [vel_y] scaling = 1e6 - [./InitialCondition] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] + [] + [] - [./vel_z] + [vel_z] scaling = 1e6 - [./InitialCondition] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] + [] + [] - [./T] - [../] + [T] + [] - [./p] + [p] scaling = 1e12 - [../] - [./disp_x] + [] + [disp_x] scaling = 1e13 - [../] - [./disp_y] + [] + [disp_y] scaling = 1e13 - [../] - [./disp_z] + [] + [disp_z] scaling = 1e14 - [../] + [] [] [ICs] - [./T] + [T] type = FunctionIC variable = T function = '(${surfacetemp} - 300) / ${pooldepth} * z + ${surfacetemp}' - [../] + [] [] [Kernels] - [./disp_x] + [disp_x] type = Diffusion variable = disp_x - [../] - [./disp_y] + [] + [disp_y] type = Diffusion variable = disp_y - [../] - [./disp_z] + [] + [disp_z] type = Diffusion variable = disp_z - [../] + [] [] [ADKernels] - [./mesh_x] + [mesh_x] type = INSConvectedMesh variable = vel_x disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] - [./mesh_y] + [] + [mesh_y] type = INSConvectedMesh variable = vel_y disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] - [./mesh_z] + [] + [mesh_z] type = INSConvectedMesh variable = vel_z disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] + [] # mass - [./mass] + [mass] type = INSADMass variable = p u = vel_x @@ -135,17 +135,17 @@ pooldepth=.2e-3 w = vel_z p = p use_displaced_mesh = true - [../] + [] # x-momentum, time - [./x_momentum_time] + [x_momentum_time] type = INSADMomentumTimeDerivative variable = vel_x use_displaced_mesh = true - [../] + [] # x-momentum, space - [./x_momentum_space] + [x_momentum_space] type = INSADMomentumBase variable = vel_x u = vel_x @@ -154,17 +154,17 @@ pooldepth=.2e-3 p = p component = 0 use_displaced_mesh = true - [../] + [] # y-momentum, time - [./y_momentum_time] + [y_momentum_time] type = INSADMomentumTimeDerivative variable = vel_y use_displaced_mesh = true - [../] + [] # y-momentum, space - [./y_momentum_space] + [y_momentum_space] type = INSADMomentumBase variable = vel_y u = vel_x @@ -173,17 +173,17 @@ pooldepth=.2e-3 p = p component = 1 use_displaced_mesh = true - [../] + [] # z-momentum, time - [./z_momentum_time] + [z_momentum_time] type = INSADMomentumTimeDerivative variable = vel_z use_displaced_mesh = true - [../] + [] # z-momentum, space - [./z_momentum_space] + [z_momentum_space] type = INSADMomentumBase variable = vel_z u = vel_x @@ -192,86 +192,86 @@ pooldepth=.2e-3 p = p component = 2 use_displaced_mesh = true - [../] - - # temperature - [./temperature_time] - type = INSADTemperatureTimeDerivative - variable = T - use_displaced_mesh = true - [../] - - [./temperature_space] - type = INSADTemperature - variable = T - u = vel_x - v = vel_y - w = vel_z - p = p - use_displaced_mesh = true - [../] - -[mesh_T] - type = INSTemperatureConvectedMesh - variable = T - disp_x = disp_x - disp_y = disp_y - disp_z = disp_z - use_displaced_mesh = true - [] + [] + + # temperature + [temperature_time] + type = INSADTemperatureTimeDerivative + variable = T + use_displaced_mesh = true + [] + + [temperature_space] + type = INSADTemperature + variable = T + u = vel_x + v = vel_y + w = vel_z + p = p + use_displaced_mesh = true + [] + + [mesh_T] + type = INSTemperatureConvectedMesh + variable = T + disp_x = disp_x + disp_y = disp_y + disp_z = disp_z + use_displaced_mesh = true + [] [] [BCs] - [./x_no_disp] + [x_no_disp] type = DirichletBC variable = disp_x boundary = 'back' value = 0 - [../] - [./y_no_disp] + [] + [y_no_disp] type = DirichletBC variable = disp_y boundary = 'back' value = 0 - [../] - [./z_no_disp] + [] + [z_no_disp] type = DirichletBC variable = disp_z boundary = 'back' value = 0 - [../] + [] - [./x_no_slip] + [x_no_slip] type = DirichletBC variable = vel_x boundary = 'bottom right left top back' value = 0.0 - [../] + [] - [./y_no_slip] + [y_no_slip] type = DirichletBC variable = vel_y boundary = 'bottom right left top back' value = 0.0 - [../] + [] - [./z_no_slip] + [z_no_slip] type = DirichletBC variable = vel_z boundary = 'bottom right left top back' value = 0.0 - [../] + [] - [./T_cold] + [T_cold] type = DirichletBC variable = T boundary = 'back' value = 300 - [../] + [] [] [ADBCs] - [./radiation_flux] + [radiation_flux] type = RadiationEnergyFluxBC variable = T boundary = 'front' @@ -279,9 +279,9 @@ pooldepth=.2e-3 sb_constant = 'sb_constant' absorptivity = 'abs' use_displaced_mesh = true - [../] + [] - [./weld_flux] + [weld_flux] type = GaussianWeldEnergyFluxBC variable = T boundary = 'front' @@ -292,104 +292,104 @@ pooldepth=.2e-3 y_beam_coord = 0 z_beam_coord = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_x] + [vapor_recoil_x] type = VaporRecoilPressureMomentumFluxBC variable = vel_x boundary = 'front' component = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_y] + [vapor_recoil_y] type = VaporRecoilPressureMomentumFluxBC variable = vel_y boundary = 'front' component = 1 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_z] + [vapor_recoil_z] type = VaporRecoilPressureMomentumFluxBC variable = vel_z boundary = 'front' component = 2 use_displaced_mesh = true - [../] + [] - [./surface_x] + [surface_x] type = SurfaceTensionBC variable = vel_x boundary = 'front' component = 0 use_displaced_mesh = true - [../] + [] - [./surface_y] + [surface_y] type = SurfaceTensionBC variable = vel_y boundary = 'front' component = 1 use_displaced_mesh = true - [../] + [] - [./surface_z] + [surface_z] type = SurfaceTensionBC variable = vel_z boundary = 'front' component = 2 use_displaced_mesh = true - [../] + [] -[./displace_x_top] + [displace_x_top] type = DisplaceBoundaryBC boundary = 'front' variable = 'disp_x' velocity = 'vel_x' - [../] - [./displace_y_top] + [] + [displace_y_top] type = DisplaceBoundaryBC boundary = 'front' variable = 'disp_y' velocity = 'vel_y' - [../] - [./displace_z_top] + [] + [displace_z_top] type = DisplaceBoundaryBC boundary = 'front' variable = 'disp_z' velocity = 'vel_z' - [../] + [] [] [ADMaterials] - [./kc_fits] + [kc_fits] type = DemonstrationPlantFits temperature = T beta = 1e7 - [../] - [./boundary] + [] + [boundary] type = DemonstrationPlantFitsBoundary use_displaced_mesh = true boundary = 'front' temperature = T - [../] + [] [] [Materials] - [./const] + [const] type = GenericConstantMaterial prop_names = 'abs sb_constant' prop_values = '1 5.67e-8' - [../] + [] [] [Preconditioning] - [./SMP] + [SMP] type = SMP full = true solve_type = 'NEWTON' - [../] + [] [] [Executioner] @@ -402,26 +402,26 @@ pooldepth=.2e-3 line_search = 'none' nl_max_its = 12 - [./TimeStepper] + [TimeStepper] type = IterationAdaptiveDT optimal_iterations = 7 dt = ${timestep} linear_iteration_ratio = 1e6 growth_factor = 1.5 - [../] + [] [] [Outputs] print_linear_residuals = false - [./exodus] + [exodus] type = Exodus output_material_properties = true show_material_properties = 'mu surface_tension grad_surface_tension' - [../] - [./dofmap] + [] + [dofmap] type = DOFMap execute_on = 'initial' - [../] + [] checkpoint = true [] @@ -429,119 +429,118 @@ pooldepth=.2e-3 show_var_residual_norms = true [] - [Adaptivity] marker = combo max_h_level = 6 - [./Indicators] - [./error_x] + [Indicators] + [error_x] type = GradientJumpIndicator variable = vel_x - [../] - [./error_y] + [] + [error_y] type = GradientJumpIndicator variable = vel_y - [../] - [./error_z] + [] + [error_z] type = GradientJumpIndicator variable = vel_z - [../] + [] # [./error_p] # type = GradientJumpIndicator # variable = p # [../] - [./error_T] + [error_T] type = GradientJumpIndicator variable = T - [../] - [./error_dispx] + [] + [error_dispx] type = GradientJumpIndicator variable = disp_x - [../] - [./error_dispy] + [] + [error_dispy] type = GradientJumpIndicator variable = disp_y - [../] - [./error_dispz] + [] + [error_dispz] type = GradientJumpIndicator variable = disp_z - [../] - [../] + [] + [] - [./Markers] - [./errorfrac_x] + [Markers] + [errorfrac_x] type = ErrorFractionMarker refine = 0.4 coarsen = 0.2 indicator = error_x - [../] - [./errorfrac_y] + [] + [errorfrac_y] type = ErrorFractionMarker refine = 0.4 coarsen = 0.2 indicator = error_y - [../] - [./errorfrac_z] + [] + [errorfrac_z] type = ErrorFractionMarker refine = 0.4 coarsen = 0.2 indicator = error_z - [../] + [] # [./errorfrac_p] # type = ErrorFractionMarker # refine = 0.7 # coarsen = 0.3 # indicator = error_p # [../] - [./errorfrac_T] + [errorfrac_T] type = ErrorFractionMarker refine = 0.4 coarsen = 0.2 indicator = error_T - [../] - [./errorfrac_dispx] + [] + [errorfrac_dispx] type = ErrorFractionMarker refine = 0.4 coarsen = 0.2 indicator = error_dispx - [../] - [./errorfrac_dispy] + [] + [errorfrac_dispy] type = ErrorFractionMarker refine = 0.4 coarsen = 0.2 indicator = error_dispy - [../] - [./errorfrac_dispz] + [] + [errorfrac_dispz] type = ErrorFractionMarker refine = 0.4 coarsen = 0.2 indicator = error_dispz - [../] - [./combo] + [] + [combo] type = ComboMarker markers = 'errorfrac_x errorfrac_y errorfrac_z errorfrac_T errorfrac_dispx errorfrac_dispy errorfrac_dispz' - [../] - [../] + [] + [] [] [Postprocessors] - [./num_dofs] + [num_dofs] type = NumDOFs system = 'NL' - [../] - [./nl] + [] + [nl] type = NumNonlinearIterations - [../] - [./tot_nl] + [] + [tot_nl] type = CumulativeValuePostprocessor postprocessor = 'nl' - [../] - [./linear] + [] + [linear] type = NumLinearIterations - [../] - [./tot_linear] + [] + [tot_linear] type = CumulativeValuePostprocessor postprocessor = 'linear' - [../] + [] [] diff --git a/examples/laser_welding_problems/physics/surface-tension/3d-scaled-1st-order.i b/examples/laser_welding_problems/physics/surface-tension/3d-scaled-1st-order.i index a966b8d..f4ef05e 100644 --- a/examples/laser_welding_problems/physics/surface-tension/3d-scaled-1st-order.i +++ b/examples/laser_welding_problems/physics/surface-tension/3d-scaled-1st-order.i @@ -1,16 +1,16 @@ -length_unit_exponent=-4 -temperature_unit_exponent=3 -mass_unit_exponent=-6 -time_unit_exponent=-5 - -endtime=${fparse 500e-5 / 10^time_unit_exponent} -period=${fparse 5e-5 / 10^time_unit_exponent} -timestep=${fparse period / 20} -surfacetemp=${fparse 300 / 10^temperature_unit_exponent} -bottomtemp=${fparse 300 / 10^temperature_unit_exponent} -width=${fparse 8e-4 / 10^length_unit_exponent} -half_width=${fparse width / 2} -pooldepth=${width} +length_unit_exponent = -4 +temperature_unit_exponent = 3 +mass_unit_exponent = -6 +time_unit_exponent = -5 + +endtime = '${fparse 500e-5 / 10^time_unit_exponent}' +period = '${fparse 5e-5 / 10^time_unit_exponent}' +timestep = '${fparse period / 20}' +surfacetemp = '${fparse 300 / 10^temperature_unit_exponent}' +bottomtemp = '${fparse 300 / 10^temperature_unit_exponent}' +width = '${fparse 8e-4 / 10^length_unit_exponent}' +half_width = '${fparse width / 2}' +pooldepth = ${width} [GlobalParams] gravity = '0 0 0' @@ -27,11 +27,11 @@ pooldepth=${width} [Mesh] type = GeneratedMesh dim = 3 - xmin = ${fparse -half_width} + xmin = '${fparse -half_width}' xmax = ${half_width} - ymin = ${fparse -half_width} + ymin = '${fparse -half_width}' ymax = ${half_width} - zmin = ${fparse -pooldepth} + zmin = '${fparse -pooldepth}' zmax = 0 nx = 2 ny = 2 @@ -41,120 +41,120 @@ pooldepth=${width} [] [MeshModifiers] - [./corner_node] + [corner_node] type = AddExtraNodeset new_boundary = 'pinned_node' nodes = '0' - [../] + [] [] [Variables] - [./vel_x] - [./InitialCondition] + [vel_x] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] + [] scaling = 1e6 - [../] + [] - [./vel_y] - [./InitialCondition] + [vel_y] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] + [] scaling = 1e6 - [../] + [] - [./vel_z] - [./InitialCondition] + [vel_z] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] + [] scaling = 1e6 - [../] + [] -[./T] + [T] scaling = 1e1 - [../] + [] - [./p] + [p] order = FIRST scaling = 1e4 - [../] - [./disp_x] + [] + [disp_x] scaling = 1e-3 - [../] - [./disp_y] + [] + [disp_y] scaling = 1e-3 - [../] - [./disp_z] + [] + [disp_z] scaling = 1e-3 - [../] + [] [] [ICs] - [./T] + [T] type = FunctionIC variable = T function = '(${surfacetemp} - ${bottomtemp}) / ${pooldepth} * z + ${surfacetemp}' - [../] + [] [] [Kernels] [] [ADKernels] - [./disp_x] + [disp_x] type = ADStressDivergence variable = disp_x component = 0 - [../] - [./disp_y] + [] + [disp_y] type = ADStressDivergence variable = disp_y component = 1 - [../] - [./disp_z] + [] + [disp_z] type = ADStressDivergence variable = disp_z component = 2 - [../] + [] -[./mesh_x] + [mesh_x] type = INSConvectedMesh variable = vel_x disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] - [./mesh_y] + [] + [mesh_y] type = INSConvectedMesh variable = vel_y disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] - [./mesh_z] + [] + [mesh_z] type = INSConvectedMesh variable = vel_z disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] - [./mesh_T] + [] + [mesh_T] type = INSTemperatureConvectedMesh variable = T disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] + [] # mass - [./mass] + [mass] type = INSADMass variable = p u = vel_x @@ -162,17 +162,17 @@ pooldepth=${width} w = vel_z p = p use_displaced_mesh = true - [../] + [] # x-momentum, time - [./x_momentum_time] + [x_momentum_time] type = INSADMomentumTimeDerivative variable = vel_x use_displaced_mesh = true - [../] + [] # x-momentum, space - [./x_momentum_space] + [x_momentum_space] type = INSADMomentumBase variable = vel_x u = vel_x @@ -181,17 +181,17 @@ pooldepth=${width} p = p component = 0 use_displaced_mesh = true - [../] + [] # y-momentum, time - [./y_momentum_time] + [y_momentum_time] type = INSADMomentumTimeDerivative variable = vel_y use_displaced_mesh = true - [../] + [] # y-momentum, space - [./y_momentum_space] + [y_momentum_space] type = INSADMomentumBase variable = vel_y u = vel_x @@ -200,17 +200,17 @@ pooldepth=${width} p = p component = 1 use_displaced_mesh = true - [../] + [] # z-momentum, time - [./z_momentum_time] + [z_momentum_time] type = INSADMomentumTimeDerivative variable = vel_z use_displaced_mesh = true - [../] + [] # z-momentum, space - [./z_momentum_space] + [z_momentum_space] type = INSADMomentumBase variable = vel_z u = vel_x @@ -219,84 +219,84 @@ pooldepth=${width} p = p component = 2 use_displaced_mesh = true - [../] - - # temperature - [./temperature_time] - type = INSADTemperatureTimeDerivative - variable = T - use_displaced_mesh = true - [../] - - [./temperature_space] - type = INSADTemperature - variable = T - u = vel_x - v = vel_y - w = vel_z - p = p - use_displaced_mesh = true - [../] + [] + + # temperature + [temperature_time] + type = INSADTemperatureTimeDerivative + variable = T + use_displaced_mesh = true + [] + + [temperature_space] + type = INSADTemperature + variable = T + u = vel_x + v = vel_y + w = vel_z + p = p + use_displaced_mesh = true + [] [] [BCs] - [./x_no_disp] + [x_no_disp] type = DirichletBC variable = disp_x boundary = 'left right' value = 0 - [../] - [./y_no_disp] + [] + [y_no_disp] type = DirichletBC variable = disp_y boundary = 'bottom top' value = 0 - [../] - [./z_no_disp] + [] + [z_no_disp] type = DirichletBC variable = disp_z boundary = 'back' value = 0 - [../] + [] - [./x_no_slip] + [x_no_slip] type = DirichletBC variable = vel_x boundary = 'bottom left right top back' value = 0.0 - [../] + [] - [./y_no_slip] + [y_no_slip] type = DirichletBC variable = vel_y boundary = 'bottom left right top back' value = 0.0 - [../] + [] - [./z_no_slip] + [z_no_slip] type = DirichletBC variable = vel_z boundary = 'bottom left right top back' value = 0.0 - [../] + [] - [./T_cold] + [T_cold] type = DirichletBC variable = T boundary = 'back' value = ${bottomtemp} - [../] + [] - [./pressure_pin] + [pressure_pin] type = DirichletBC variable = p boundary = 'pinned_node' value = 0 - [../] + [] [] [ADBCs] - [./radiation_flux] + [radiation_flux] type = RadiationEnergyFluxBC variable = T boundary = 'front' @@ -304,70 +304,70 @@ pooldepth=${width} sb_constant = 'sb_constant' absorptivity = 'abs' use_displaced_mesh = true - [../] + [] - [./weld_flux] + [weld_flux] type = GaussianWeldEnergyFluxBC variable = T boundary = 'front' reff = 0.6 - F0 = ${fparse 2.546e9 / 10^mass_unit_exponent * 10^(3 * time_unit_exponent)} - R = ${fparse 1e-4 / 10^length_unit_exponent} + F0 = '${fparse 2.546e9 / 10^mass_unit_exponent * 10^(3 * time_unit_exponent)}' + R = '${fparse 1e-4 / 10^length_unit_exponent}' x_beam_coord = '${fparse 2e-4 / 10^length_unit_exponent} * cos(2 * pi / (50 * ${timestep}) * t)' y_beam_coord = '${fparse 2e-4 / 10^length_unit_exponent} * sin(2 * pi / (50 * ${timestep}) * t)' z_beam_coord = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_x] + [vapor_recoil_x] type = VaporRecoilPressureMomentumFluxBC variable = vel_x boundary = 'front' component = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_y] + [vapor_recoil_y] type = VaporRecoilPressureMomentumFluxBC variable = vel_y boundary = 'front' component = 1 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_z] + [vapor_recoil_z] type = VaporRecoilPressureMomentumFluxBC variable = vel_z boundary = 'front' component = 2 use_displaced_mesh = true - [../] + [] - [./surface_x] + [surface_x] type = SurfaceTensionBC variable = vel_x boundary = 'front' component = 0 use_displaced_mesh = true - [../] + [] - [./surface_y] + [surface_y] type = SurfaceTensionBC variable = vel_y boundary = 'front' component = 1 use_displaced_mesh = true - [../] + [] - [./surface_z] + [surface_z] type = SurfaceTensionBC variable = vel_z boundary = 'front' component = 2 use_displaced_mesh = true - [../] + [] - [./displace_x_top] + [displace_x_top] type = PenaltyDisplaceBoundaryBC boundary = 'front' variable = 'disp_x' @@ -377,8 +377,8 @@ pooldepth=${width} disp_y = 'disp_y' vel_z = 'vel_z' disp_z = 'disp_z' - [../] - [./displace_y_top] + [] + [displace_y_top] type = PenaltyDisplaceBoundaryBC boundary = 'front' variable = 'disp_y' @@ -388,8 +388,8 @@ pooldepth=${width} disp_y = 'disp_y' vel_z = 'vel_z' disp_z = 'disp_z' - [../] - [./displace_z_top] + [] + [displace_z_top] type = PenaltyDisplaceBoundaryBC boundary = 'front' variable = 'disp_z' @@ -399,11 +399,11 @@ pooldepth=${width} disp_y = 'disp_y' vel_z = 'vel_z' disp_z = 'disp_z' - [../] + [] [] [ADMaterials] - [./kc_fits] + [kc_fits] type = DemonstrationPlantFits temperature = T beta = 1e7 @@ -411,8 +411,8 @@ pooldepth=${width} temperature_unit_exponent = ${temperature_unit_exponent} mass_unit_exponent = ${mass_unit_exponent} time_unit_exponent = ${time_unit_exponent} - [../] - [./boundary] + [] + [boundary] type = DemonstrationPlantFitsBoundary boundary = 'front' temperature = T @@ -421,59 +421,59 @@ pooldepth=${width} temperature_unit_exponent = ${temperature_unit_exponent} mass_unit_exponent = ${mass_unit_exponent} time_unit_exponent = ${time_unit_exponent} - [../] - [./stress] + [] + [stress] type = PseudoSolidStress disp_x = disp_x disp_y = disp_y disp_z = disp_z - [../] + [] [] [Materials] - [./const] + [const] type = GenericConstantMaterial prop_names = 'abs sb_constant' prop_values = '1 ${fparse 5.67e-8 / 10^mass_unit_exponent * 10^(4 * temperature_unit_exponent) * 10^(3 * time_unit_exponent)}' - [../] + [] [] [Preconditioning] - [./SMP] + [SMP] type = SMP full = true solve_type = 'NEWTON' - [../] + [] [] [Executioner] type = Transient end_time = ${endtime} - dtmin = ${fparse 1e-8 / 10^time_unit_exponent} + dtmin = '${fparse 1e-8 / 10^time_unit_exponent}' petsc_options = '-snes_converged_reason -ksp_converged_reason -options_left -ksp_monitor_singular_value' petsc_options_iname = '-ksp_max_it -ksp_gmres_restart -pc_type -snes_max_funcs -sub_pc_factor_levels' petsc_options_value = '100 100 asm 1000000 1' line_search = 'none' nl_max_its = 12 - l_tol = 1e-2 - [./TimeStepper] + l_tol = 1e-2 + [TimeStepper] type = IterationAdaptiveDT optimal_iterations = 8 dt = ${timestep} linear_iteration_ratio = 1e6 growth_factor = 1.1 - [../] + [] [] [Outputs] print_linear_residuals = false csv = true - [./exodus] + [exodus] type = Exodus output_material_properties = true show_material_properties = 'mu surface_term_curvature surface_term_gradient1 surface_term_gradient2' - [../] + [] # [./dofmap] # type = DOFMap # execute_on = 'timestep_end' @@ -485,120 +485,119 @@ pooldepth=${width} show_var_residual_norms = true [] - [Adaptivity] marker = combo max_h_level = 4 initial_steps = 0 - [./Indicators] - [./error_x] + [Indicators] + [error_x] type = GradientJumpIndicator variable = vel_x - [../] - [./error_y] + [] + [error_y] type = GradientJumpIndicator variable = vel_y - [../] - [./error_z] + [] + [error_z] type = GradientJumpIndicator variable = vel_z - [../] + [] # [./error_p] # type = GradientJumpIndicator # variable = p # [../] - [./error_T] + [error_T] type = GradientJumpIndicator variable = T - [../] - [./error_dispx] + [] + [error_dispx] type = GradientJumpIndicator variable = disp_x - [../] - [./error_dispy] + [] + [error_dispy] type = GradientJumpIndicator variable = disp_y - [../] - [./error_dispz] + [] + [error_dispz] type = GradientJumpIndicator variable = disp_z - [../] - [../] + [] + [] - [./Markers] - [./errorfrac_x] + [Markers] + [errorfrac_x] type = ErrorFractionMarker refine = 0.8 coarsen = 0.05 indicator = error_x - [../] - [./errorfrac_y] + [] + [errorfrac_y] type = ErrorFractionMarker refine = 0.8 coarsen = 0.05 indicator = error_y - [../] - [./errorfrac_z] + [] + [errorfrac_z] type = ErrorFractionMarker refine = 0.8 coarsen = 0.05 indicator = error_z - [../] + [] # [./errorfrac_p] # type = ErrorFractionMarker # refine = 0.8 # coarsen = 0.05 # indicator = error_p # [../] - [./errorfrac_T] + [errorfrac_T] type = ErrorFractionMarker refine = 0.8 coarsen = 0.05 indicator = error_T - [../] - [./errorfrac_dispx] + [] + [errorfrac_dispx] type = ErrorFractionMarker refine = 0.8 coarsen = 0.05 indicator = error_dispx - [../] - [./errorfrac_dispy] + [] + [errorfrac_dispy] type = ErrorFractionMarker refine = 0.8 coarsen = 0.05 indicator = error_dispy - [../] - [./errorfrac_dispz] + [] + [errorfrac_dispz] type = ErrorFractionMarker refine = 0.8 coarsen = 0.05 indicator = error_dispz - [../] - [./combo] + [] + [combo] type = ComboMarker markers = 'errorfrac_x errorfrac_y errorfrac_z errorfrac_T errorfrac_dispx errorfrac_dispy errorfrac_dispz' - [../] - [../] + [] + [] [] [Postprocessors] - [./num_dofs] + [num_dofs] type = NumDOFs system = 'NL' - [../] - [./nl] + [] + [nl] type = NumNonlinearIterations - [../] - [./tot_nl] + [] + [tot_nl] type = CumulativeValuePostprocessor postprocessor = 'nl' - [../] - [./linear] + [] + [linear] type = NumLinearIterations - [../] - [./tot_linear] + [] + [tot_linear] type = CumulativeValuePostprocessor postprocessor = 'linear' - [../] + [] [] diff --git a/examples/laser_welding_problems/physics/surface-tension/cold-start-twice-around-circle.i b/examples/laser_welding_problems/physics/surface-tension/cold-start-twice-around-circle.i index d5a1861..458e364 100644 --- a/examples/laser_welding_problems/physics/surface-tension/cold-start-twice-around-circle.i +++ b/examples/laser_welding_problems/physics/surface-tension/cold-start-twice-around-circle.i @@ -1,7 +1,7 @@ -period=1.25e-3 -endtime=2.5e-3 -timestep=1.25e-5 -surfacetemp=300 +period = 1.25e-3 +endtime = 2.5e-3 +timestep = 1.25e-5 +surfacetemp = 300 [GlobalParams] gravity = '0 0 0' @@ -31,91 +31,91 @@ surfacetemp=300 [] [Variables] - [./vel_x] - [./InitialCondition] + [vel_x] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] + [] + [] - [./vel_y] - [./InitialCondition] + [vel_y] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] + [] + [] - [./vel_z] - [./InitialCondition] + [vel_z] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] - - [./T] - [../] - - [./p] - [../] - [./disp_x] - [../] - [./disp_y] - [../] - [./disp_z] - [../] + [] + [] + + [T] + [] + + [p] + [] + [disp_x] + [] + [disp_y] + [] + [disp_z] + [] [] [ICs] - [./T] + [T] type = FunctionIC variable = T function = '(${surfacetemp} - 300) / .7e-3 * z + ${surfacetemp}' - [../] + [] [] [Kernels] - [./disp_x] + [disp_x] type = Diffusion variable = disp_x - [../] - [./disp_y] + [] + [disp_y] type = Diffusion variable = disp_y - [../] - [./disp_z] + [] + [disp_z] type = Diffusion variable = disp_z - [../] + [] [] [ADKernels] - [./mesh_x] + [mesh_x] type = INSConvectedMesh variable = vel_x disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] - [./mesh_y] + [] + [mesh_y] type = INSConvectedMesh variable = vel_y disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] - [./mesh_z] + [] + [mesh_z] type = INSConvectedMesh variable = vel_z disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] + [] # mass - [./mass] + [mass] type = INSADMass variable = p u = vel_x @@ -123,17 +123,17 @@ surfacetemp=300 w = vel_z p = p use_displaced_mesh = true - [../] + [] # x-momentum, time - [./x_momentum_time] + [x_momentum_time] type = INSADMomentumTimeDerivative variable = vel_x use_displaced_mesh = true - [../] + [] # x-momentum, space - [./x_momentum_space] + [x_momentum_space] type = INSADMomentumBase variable = vel_x u = vel_x @@ -142,17 +142,17 @@ surfacetemp=300 p = p component = 0 use_displaced_mesh = true - [../] + [] # y-momentum, time - [./y_momentum_time] + [y_momentum_time] type = INSADMomentumTimeDerivative variable = vel_y use_displaced_mesh = true - [../] + [] # y-momentum, space - [./y_momentum_space] + [y_momentum_space] type = INSADMomentumBase variable = vel_y u = vel_x @@ -161,17 +161,17 @@ surfacetemp=300 p = p component = 1 use_displaced_mesh = true - [../] + [] # z-momentum, time - [./z_momentum_time] + [z_momentum_time] type = INSADMomentumTimeDerivative variable = vel_z use_displaced_mesh = true - [../] + [] # z-momentum, space - [./z_momentum_space] + [z_momentum_space] type = INSADMomentumBase variable = vel_z u = vel_x @@ -180,86 +180,86 @@ surfacetemp=300 p = p component = 2 use_displaced_mesh = true - [../] - - # temperature - [./temperature_time] - type = INSADTemperatureTimeDerivative - variable = T - use_displaced_mesh = true - [../] - - [./temperature_space] - type = INSADTemperature - variable = T - u = vel_x - v = vel_y - w = vel_z - p = p - use_displaced_mesh = true - [../] - - [mesh_T] - type = INSTemperatureConvectedMesh - variable = T - disp_x = disp_x - disp_y = disp_y - disp_z = disp_z - use_displaced_mesh = true - [] + [] + + # temperature + [temperature_time] + type = INSADTemperatureTimeDerivative + variable = T + use_displaced_mesh = true + [] + + [temperature_space] + type = INSADTemperature + variable = T + u = vel_x + v = vel_y + w = vel_z + p = p + use_displaced_mesh = true + [] + + [mesh_T] + type = INSTemperatureConvectedMesh + variable = T + disp_x = disp_x + disp_y = disp_y + disp_z = disp_z + use_displaced_mesh = true + [] [] [BCs] - [./x_no_disp] + [x_no_disp] type = DirichletBC variable = disp_x boundary = 'back' value = 0 - [../] - [./y_no_disp] + [] + [y_no_disp] type = DirichletBC variable = disp_y boundary = 'back' value = 0 - [../] - [./z_no_disp] + [] + [z_no_disp] type = DirichletBC variable = disp_z boundary = 'back' value = 0 - [../] + [] - [./x_no_slip] + [x_no_slip] type = DirichletBC variable = vel_x boundary = 'bottom right left top back' value = 0.0 - [../] + [] - [./y_no_slip] + [y_no_slip] type = DirichletBC variable = vel_y boundary = 'bottom right left top back' value = 0.0 - [../] + [] - [./z_no_slip] + [z_no_slip] type = DirichletBC variable = vel_z boundary = 'bottom right left top back' value = 0.0 - [../] + [] - [./T_cold] + [T_cold] type = DirichletBC variable = T boundary = 'back' value = 300 - [../] + [] [] [ADBCs] - [./radiation_flux] + [radiation_flux] type = RadiationEnergyFluxBC variable = T boundary = 'front' @@ -267,9 +267,9 @@ surfacetemp=300 sb_constant = 'sb_constant' absorptivity = 'abs' use_displaced_mesh = true - [../] + [] - [./weld_flux] + [weld_flux] type = GaussianWeldEnergyFluxBC variable = T boundary = 'front' @@ -280,104 +280,104 @@ surfacetemp=300 y_beam_coord = '2e-4 * sin(t * 2 * pi / ${period})' z_beam_coord = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_x] + [vapor_recoil_x] type = VaporRecoilPressureMomentumFluxBC variable = vel_x boundary = 'front' component = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_y] + [vapor_recoil_y] type = VaporRecoilPressureMomentumFluxBC variable = vel_y boundary = 'front' component = 1 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_z] + [vapor_recoil_z] type = VaporRecoilPressureMomentumFluxBC variable = vel_z boundary = 'front' component = 2 use_displaced_mesh = true - [../] + [] - [./surface_x] + [surface_x] type = SurfaceTensionBC variable = vel_x boundary = 'front' component = 0 use_displaced_mesh = true - [../] + [] - [./surface_y] + [surface_y] type = SurfaceTensionBC variable = vel_y boundary = 'front' component = 1 use_displaced_mesh = true - [../] + [] - [./surface_z] + [surface_z] type = SurfaceTensionBC variable = vel_z boundary = 'front' component = 2 use_displaced_mesh = true - [../] + [] -[./displace_x_top] + [displace_x_top] type = DisplaceBoundaryBC boundary = 'front' variable = 'disp_x' velocity = 'vel_x' - [../] - [./displace_y_top] + [] + [displace_y_top] type = DisplaceBoundaryBC boundary = 'front' variable = 'disp_y' velocity = 'vel_y' - [../] - [./displace_z_top] + [] + [displace_z_top] type = DisplaceBoundaryBC boundary = 'front' variable = 'disp_z' velocity = 'vel_z' - [../] + [] [] [ADMaterials] - [./kc_fits] + [kc_fits] type = DemonstrationPlantFits temperature = T beta = 1e7 - [../] - [./boundary] + [] + [boundary] type = DemonstrationPlantFitsBoundary boundary = 'front' temperature = T use_displaced_mesh = true - [../] + [] [] [Materials] - [./const] + [const] type = GenericConstantMaterial prop_names = 'abs sb_constant' prop_values = '1 5.67e-8' - [../] + [] [] [Preconditioning] - [./SMP] + [SMP] type = SMP full = true solve_type = 'NEWTON' - [../] + [] [] [Executioner] @@ -392,26 +392,26 @@ surfacetemp=300 line_search = 'none' nl_max_its = 12 l_max_its = 100 - [./TimeStepper] + [TimeStepper] type = IterationAdaptiveDT optimal_iterations = 7 dt = ${timestep} linear_iteration_ratio = 1e6 growth_factor = 1.5 - [../] + [] [] [Outputs] print_linear_residuals = false - [./exodus] + [exodus] type = Exodus output_material_properties = true show_material_properties = 'mu' - [../] - [./dofmap] + [] + [dofmap] type = DOFMap execute_on = 'initial' - [../] + [] checkpoint = true [] @@ -419,12 +419,11 @@ surfacetemp=300 show_var_residual_norms = true [] - [Adaptivity] marker = combo max_h_level = 4 - [./Indicators] + [Indicators] # [./error_x] # type = GradientJumpIndicator # variable = vel_x @@ -441,10 +440,10 @@ surfacetemp=300 # type = GradientJumpIndicator # variable = p # [../] - [./error_T] + [error_T] type = GradientJumpIndicator variable = T - [../] + [] # [./error_dispx] # type = GradientJumpIndicator # variable = disp_x @@ -453,13 +452,13 @@ surfacetemp=300 # type = GradientJumpIndicator # variable = disp_y # [../] - [./error_dispz] + [error_dispz] type = GradientJumpIndicator variable = disp_z - [../] - [../] + [] + [] - [./Markers] + [Markers] # [./errorfrac_x] # type = ErrorFractionMarker # refine = 0.7 @@ -484,12 +483,12 @@ surfacetemp=300 # coarsen = 0.3 # indicator = error_p # [../] - [./errorfrac_T] + [errorfrac_T] type = ErrorFractionMarker refine = 0.4 coarsen = 0.2 indicator = error_T - [../] + [] # [./errorfrac_dispx] # type = ErrorFractionMarker # refine = 0.7 @@ -502,36 +501,36 @@ surfacetemp=300 # coarsen = 0.3 # indicator = error_dispy # [../] - [./errorfrac_dispz] + [errorfrac_dispz] type = ErrorFractionMarker refine = 0.4 coarsen = 0.2 indicator = error_dispz - [../] - [./combo] + [] + [combo] type = ComboMarker markers = 'errorfrac_T errorfrac_dispz' - [../] - [../] + [] + [] [] [Postprocessors] - [./num_dofs] + [num_dofs] type = NumDOFs system = 'NL' - [../] - [./nl] + [] + [nl] type = NumNonlinearIterations - [../] - [./tot_nl] + [] + [tot_nl] type = CumulativeValuePostprocessor postprocessor = 'nl' - [../] - [./linear] + [] + [linear] type = NumLinearIterations - [../] - [./tot_linear] + [] + [tot_linear] type = CumulativeValuePostprocessor postprocessor = 'linear' - [../] + [] [] diff --git a/examples/laser_welding_problems/physics/surface-tension/gold-2d-one-spot-scaled-1st-order.i b/examples/laser_welding_problems/physics/surface-tension/gold-2d-one-spot-scaled-1st-order.i index ebc5273..b3e67f4 100644 --- a/examples/laser_welding_problems/physics/surface-tension/gold-2d-one-spot-scaled-1st-order.i +++ b/examples/laser_welding_problems/physics/surface-tension/gold-2d-one-spot-scaled-1st-order.i @@ -1,16 +1,16 @@ -length_unit_exponent=-4 -temperature_unit_exponent=3 -mass_unit_exponent=-6 -time_unit_exponent=-5 - -endtime=${fparse 500e-5 / 10^time_unit_exponent} -period=${fparse 5e-5 / 10^time_unit_exponent} -timestep=${fparse period / 20} -surfacetemp=${fparse 300 / 10^temperature_unit_exponent} -bottomtemp=${fparse 300 / 10^temperature_unit_exponent} -width=${fparse 8e-4 / 10^length_unit_exponent} -half_width=${fparse width / 2} -pooldepth=${width} +length_unit_exponent = -4 +temperature_unit_exponent = 3 +mass_unit_exponent = -6 +time_unit_exponent = -5 + +endtime = '${fparse 500e-5 / 10^time_unit_exponent}' +period = '${fparse 5e-5 / 10^time_unit_exponent}' +timestep = '${fparse period / 20}' +surfacetemp = '${fparse 300 / 10^temperature_unit_exponent}' +bottomtemp = '${fparse 300 / 10^temperature_unit_exponent}' +width = '${fparse 8e-4 / 10^length_unit_exponent}' +half_width = '${fparse width / 2}' +pooldepth = ${width} [GlobalParams] gravity = '0 0 0' @@ -27,9 +27,9 @@ pooldepth=${width} [Mesh] type = GeneratedMesh dim = 2 - xmin = ${fparse -half_width} + xmin = '${fparse -half_width}' xmax = ${half_width} - ymin = ${fparse -pooldepth} + ymin = '${fparse -pooldepth}' ymax = 0 nx = 4 ny = 4 @@ -43,102 +43,102 @@ pooldepth=${width} [] [Variables] - [./vel_x] - [./InitialCondition] + [vel_x] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] + [] scaling = 1e6 - [../] + [] - [./vel_y] - [./InitialCondition] + [vel_y] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] + [] scaling = 1e6 - [../] + [] - [./T] + [T] scaling = 1e1 - [../] + [] - [./p] + [p] order = FIRST scaling = 1e4 - [../] - [./disp_x] + [] + [disp_x] scaling = 1e-3 - [../] - [./disp_y] + [] + [disp_y] scaling = 1e-3 - [../] + [] [] [ICs] - [./T] + [T] type = FunctionIC variable = T function = '(${surfacetemp} - ${bottomtemp}) / ${pooldepth} * z + ${surfacetemp}' - [../] + [] [] [Kernels] [] [ADKernels] - [./disp_x] + [disp_x] type = ADStressDivergence variable = disp_x component = 0 - [../] - [./disp_y] + [] + [disp_y] type = ADStressDivergence variable = disp_y component = 1 - [../] + [] -[./mesh_x] + [mesh_x] type = INSConvectedMesh variable = vel_x disp_x = disp_x disp_y = disp_y use_displaced_mesh = true - [../] - [./mesh_y] + [] + [mesh_y] type = INSConvectedMesh variable = vel_y disp_x = disp_x disp_y = disp_y use_displaced_mesh = true - [../] - [./mesh_T] + [] + [mesh_T] type = INSTemperatureConvectedMesh variable = T disp_x = disp_x disp_y = disp_y use_displaced_mesh = true - [../] + [] # mass - [./mass] + [mass] type = INSADMass variable = p u = vel_x v = vel_y p = p use_displaced_mesh = true - [../] + [] # x-momentum, time - [./x_momentum_time] + [x_momentum_time] type = INSADMomentumTimeDerivative variable = vel_x use_displaced_mesh = true - [../] + [] # x-momentum, space - [./x_momentum_space] + [x_momentum_space] type = INSADMomentumBase variable = vel_x u = vel_x @@ -146,17 +146,17 @@ pooldepth=${width} p = p component = 0 use_displaced_mesh = true - [../] + [] # y-momentum, time - [./y_momentum_time] + [y_momentum_time] type = INSADMomentumTimeDerivative variable = vel_y use_displaced_mesh = true - [../] + [] # y-momentum, space - [./y_momentum_space] + [y_momentum_space] type = INSADMomentumBase variable = vel_y u = vel_x @@ -164,63 +164,63 @@ pooldepth=${width} p = p component = 1 use_displaced_mesh = true - [../] - - # temperature - [./temperature_time] - type = INSADTemperatureTimeDerivative - variable = T - use_displaced_mesh = true - [../] - - [./temperature_space] - type = INSADTemperature - variable = T - u = vel_x - v = vel_y - p = p - use_displaced_mesh = true - [../] + [] + + # temperature + [temperature_time] + type = INSADTemperatureTimeDerivative + variable = T + use_displaced_mesh = true + [] + + [temperature_space] + type = INSADTemperature + variable = T + u = vel_x + v = vel_y + p = p + use_displaced_mesh = true + [] [] [BCs] - [./x_no_disp] + [x_no_disp] type = DirichletBC variable = disp_x boundary = 'left right' value = 0 - [../] - [./y_no_disp] + [] + [y_no_disp] type = DirichletBC variable = disp_y boundary = 'bottom' value = 0 - [../] + [] - [./x_no_slip] + [x_no_slip] type = DirichletBC variable = vel_x boundary = 'bottom left right' value = 0.0 - [../] + [] - [./y_no_slip] + [y_no_slip] type = DirichletBC variable = vel_y boundary = 'bottom left right' value = 0.0 - [../] + [] - [./T_cold] + [T_cold] type = DirichletBC variable = T boundary = 'bottom' value = ${bottomtemp} - [../] + [] [] [ADBCs] - [./radiation_flux] + [radiation_flux] type = RadiationEnergyFluxBC variable = T boundary = 'top' @@ -228,54 +228,54 @@ pooldepth=${width} sb_constant = 'sb_constant' absorptivity = 'abs' use_displaced_mesh = true - [../] + [] - [./weld_flux] + [weld_flux] type = GaussianWeldEnergyFluxBC variable = T boundary = 'top' reff = 0.6 - F0 = ${fparse 2.546e9 / 10^mass_unit_exponent * 10^(3 * time_unit_exponent)} - R = ${fparse 1e-4 / 10^length_unit_exponent} + F0 = '${fparse 2.546e9 / 10^mass_unit_exponent * 10^(3 * time_unit_exponent)}' + R = '${fparse 1e-4 / 10^length_unit_exponent}' x_beam_coord = '${fparse 2e-4 / 10^length_unit_exponent} * sin(2 * pi / (50 * ${timestep}) * t)' y_beam_coord = 0 z_beam_coord = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_x] + [vapor_recoil_x] type = VaporRecoilPressureMomentumFluxBC variable = vel_x boundary = 'top' component = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_y] + [vapor_recoil_y] type = VaporRecoilPressureMomentumFluxBC variable = vel_y boundary = 'top' component = 1 use_displaced_mesh = true - [../] + [] - [./surface_x] + [surface_x] type = SurfaceTensionBC variable = vel_x boundary = 'top' component = 0 use_displaced_mesh = true - [../] + [] - [./surface_y] + [surface_y] type = SurfaceTensionBC variable = vel_y boundary = 'top' component = 1 use_displaced_mesh = true - [../] + [] - [./displace_x_top] + [displace_x_top] type = PenaltyDisplaceBoundaryBC boundary = 'top' variable = 'disp_x' @@ -283,8 +283,8 @@ pooldepth=${width} disp_x = 'disp_x' vel_y = 'vel_y' disp_y = 'disp_y' - [../] - [./displace_y_top] + [] + [displace_y_top] type = PenaltyDisplaceBoundaryBC boundary = 'top' variable = 'disp_y' @@ -292,11 +292,11 @@ pooldepth=${width} disp_x = 'disp_x' vel_y = 'vel_y' disp_y = 'disp_y' - [../] + [] [] [ADMaterials] - [./kc_fits] + [kc_fits] type = DemonstrationPlantFits temperature = T beta = 1e7 @@ -304,8 +304,8 @@ pooldepth=${width} temperature_unit_exponent = ${temperature_unit_exponent} mass_unit_exponent = ${mass_unit_exponent} time_unit_exponent = ${time_unit_exponent} - [../] - [./boundary] + [] + [boundary] type = DemonstrationPlantFitsBoundary boundary = 'top' temperature = T @@ -314,58 +314,58 @@ pooldepth=${width} temperature_unit_exponent = ${temperature_unit_exponent} mass_unit_exponent = ${mass_unit_exponent} time_unit_exponent = ${time_unit_exponent} - [../] - [./stress] + [] + [stress] type = PseudoSolidStress disp_x = disp_x disp_y = disp_y - [../] + [] [] [Materials] - [./const] + [const] type = GenericConstantMaterial prop_names = 'abs sb_constant' prop_values = '1 ${fparse 5.67e-8 / 10^mass_unit_exponent * 10^(4 * temperature_unit_exponent) * 10^(3 * time_unit_exponent)}' - [../] + [] [] [Preconditioning] - [./SMP] + [SMP] type = SMP full = true solve_type = 'NEWTON' - [../] + [] [] [Executioner] type = Transient end_time = ${endtime} - dtmin = ${fparse 1e-8 / 10^time_unit_exponent} + dtmin = '${fparse 1e-8 / 10^time_unit_exponent}' petsc_options = '-snes_converged_reason -ksp_converged_reason -options_left -ksp_monitor_singular_value' petsc_options_iname = '-ksp_max_it -ksp_gmres_restart -pc_type -snes_max_funcs -sub_pc_factor_levels' petsc_options_value = '100 100 asm 1000000 1' line_search = 'none' nl_max_its = 12 - l_tol = 1e-2 - [./TimeStepper] + l_tol = 1e-2 + [TimeStepper] type = IterationAdaptiveDT optimal_iterations = 8 dt = ${timestep} linear_iteration_ratio = 1e6 growth_factor = 1.1 - [../] + [] [] [Outputs] print_linear_residuals = false csv = true - [./exodus] + [exodus] type = Exodus output_material_properties = true show_material_properties = 'mu surface_term_curvature surface_term_gradient1 surface_term_gradient2' - [../] + [] # [./dofmap] # type = DOFMap # execute_on = 'timestep_end' @@ -377,101 +377,100 @@ pooldepth=${width} show_var_residual_norms = true [] - [Adaptivity] marker = combo max_h_level = 5 initial_steps = 0 - [./Indicators] - [./error_x] + [Indicators] + [error_x] type = GradientJumpIndicator variable = vel_x - [../] - [./error_y] + [] + [error_y] type = GradientJumpIndicator variable = vel_y - [../] + [] # [./error_p] # type = GradientJumpIndicator # variable = p # [../] - [./error_T] + [error_T] type = GradientJumpIndicator variable = T - [../] - [./error_dispx] + [] + [error_dispx] type = GradientJumpIndicator variable = disp_x - [../] - [./error_dispy] + [] + [error_dispy] type = GradientJumpIndicator variable = disp_y - [../] - [../] + [] + [] - [./Markers] - [./errorfrac_x] + [Markers] + [errorfrac_x] type = ErrorFractionMarker refine = 0.8 coarsen = 0.05 indicator = error_x - [../] - [./errorfrac_y] + [] + [errorfrac_y] type = ErrorFractionMarker refine = 0.8 coarsen = 0.05 indicator = error_y - [../] + [] # [./errorfrac_p] # type = ErrorFractionMarker # refine = 0.8 # coarsen = 0.05 # indicator = error_p # [../] - [./errorfrac_T] + [errorfrac_T] type = ErrorFractionMarker refine = 0.8 coarsen = 0.05 indicator = error_T - [../] - [./errorfrac_dispx] + [] + [errorfrac_dispx] type = ErrorFractionMarker refine = 0.8 coarsen = 0.05 indicator = error_dispx - [../] - [./errorfrac_dispy] + [] + [errorfrac_dispy] type = ErrorFractionMarker refine = 0.8 coarsen = 0.05 indicator = error_dispy - [../] - [./combo] + [] + [combo] type = ComboMarker markers = 'errorfrac_x errorfrac_y errorfrac_T errorfrac_dispx errorfrac_dispy' # markers = 'errorfrac_x' - [../] - [../] + [] + [] [] [Postprocessors] - [./num_dofs] + [num_dofs] type = NumDOFs system = 'NL' - [../] - [./nl] + [] + [nl] type = NumNonlinearIterations - [../] - [./tot_nl] + [] + [tot_nl] type = CumulativeValuePostprocessor postprocessor = 'nl' - [../] - [./linear] + [] + [linear] type = NumLinearIterations - [../] - [./tot_linear] + [] + [tot_linear] type = CumulativeValuePostprocessor postprocessor = 'linear' - [../] + [] [] diff --git a/examples/laser_welding_problems/physics/surface-tension/hot-once-around-circle.i b/examples/laser_welding_problems/physics/surface-tension/hot-once-around-circle.i index dda8ad6..b7ab28f 100644 --- a/examples/laser_welding_problems/physics/surface-tension/hot-once-around-circle.i +++ b/examples/laser_welding_problems/physics/surface-tension/hot-once-around-circle.i @@ -1,5 +1,5 @@ -period=4e-5 -timestep=4e-7 +period = 4e-5 +timestep = 4e-7 [GlobalParams] gravity = '0 0 0' @@ -29,91 +29,91 @@ timestep=4e-7 [] [Variables] - [./vel_x] - [./InitialCondition] + [vel_x] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] + [] + [] - [./vel_y] - [./InitialCondition] + [vel_y] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] + [] + [] - [./vel_z] - [./InitialCondition] + [vel_z] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] - - [./T] - [../] - - [./p] - [../] - [./disp_x] - [../] - [./disp_y] - [../] - [./disp_z] - [../] + [] + [] + + [T] + [] + + [p] + [] + [disp_x] + [] + [disp_y] + [] + [disp_z] + [] [] [ICs] - [./T] + [T] type = FunctionIC variable = T function = '(3000 - 300) / .7e-3 * z + 3000' - [../] + [] [] [Kernels] - [./disp_x] + [disp_x] type = Diffusion variable = disp_x - [../] - [./disp_y] + [] + [disp_y] type = Diffusion variable = disp_y - [../] - [./disp_z] + [] + [disp_z] type = Diffusion variable = disp_z - [../] + [] [] [ADKernels] - [./mesh_x] + [mesh_x] type = INSConvectedMesh variable = vel_x disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] - [./mesh_y] + [] + [mesh_y] type = INSConvectedMesh variable = vel_y disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] - [./mesh_z] + [] + [mesh_z] type = INSConvectedMesh variable = vel_z disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] + [] -# mass - [./mass] + # mass + [mass] type = INSADMass variable = p u = vel_x @@ -121,17 +121,17 @@ timestep=4e-7 w = vel_z p = p use_displaced_mesh = true - [../] + [] # x-momentum, time - [./x_momentum_time] + [x_momentum_time] type = INSADMomentumTimeDerivative variable = vel_x use_displaced_mesh = true - [../] + [] # x-momentum, space - [./x_momentum_space] + [x_momentum_space] type = INSADMomentumBase variable = vel_x u = vel_x @@ -140,17 +140,17 @@ timestep=4e-7 p = p component = 0 use_displaced_mesh = true - [../] + [] # y-momentum, time - [./y_momentum_time] + [y_momentum_time] type = INSADMomentumTimeDerivative variable = vel_y use_displaced_mesh = true - [../] + [] # y-momentum, space - [./y_momentum_space] + [y_momentum_space] type = INSADMomentumBase variable = vel_y u = vel_x @@ -159,17 +159,17 @@ timestep=4e-7 p = p component = 1 use_displaced_mesh = true - [../] + [] # z-momentum, time - [./z_momentum_time] + [z_momentum_time] type = INSADMomentumTimeDerivative variable = vel_z use_displaced_mesh = true - [../] + [] # z-momentum, space - [./z_momentum_space] + [z_momentum_space] type = INSADMomentumBase variable = vel_z u = vel_x @@ -178,86 +178,86 @@ timestep=4e-7 p = p component = 2 use_displaced_mesh = true - [../] - - # temperature - [./temperature_time] - type = INSADTemperatureTimeDerivative - variable = T - use_displaced_mesh = true - [../] - - [./temperature_space] - type = INSADTemperature - variable = T - u = vel_x - v = vel_y - w = vel_z - p = p - use_displaced_mesh = true - [../] - - [mesh_T] - type = INSTemperatureConvectedMesh - variable = T - disp_x = disp_x - disp_y = disp_y - disp_z = disp_z - use_displaced_mesh = true - [] + [] + + # temperature + [temperature_time] + type = INSADTemperatureTimeDerivative + variable = T + use_displaced_mesh = true + [] + + [temperature_space] + type = INSADTemperature + variable = T + u = vel_x + v = vel_y + w = vel_z + p = p + use_displaced_mesh = true + [] + + [mesh_T] + type = INSTemperatureConvectedMesh + variable = T + disp_x = disp_x + disp_y = disp_y + disp_z = disp_z + use_displaced_mesh = true + [] [] [BCs] - [./x_no_disp] + [x_no_disp] type = DirichletBC variable = disp_x boundary = 'back' value = 0 - [../] - [./y_no_disp] + [] + [y_no_disp] type = DirichletBC variable = disp_y boundary = 'back' value = 0 - [../] - [./z_no_disp] + [] + [z_no_disp] type = DirichletBC variable = disp_z boundary = 'back' value = 0 - [../] + [] - [./x_no_slip] + [x_no_slip] type = DirichletBC variable = vel_x boundary = 'bottom right left top back' value = 0.0 - [../] + [] - [./y_no_slip] + [y_no_slip] type = DirichletBC variable = vel_y boundary = 'bottom right left top back' value = 0.0 - [../] + [] - [./z_no_slip] + [z_no_slip] type = DirichletBC variable = vel_z boundary = 'bottom right left top back' value = 0.0 - [../] + [] - [./T_cold] + [T_cold] type = DirichletBC variable = T boundary = 'back' value = 300 - [../] + [] [] [ADBCs] - [./radiation_flux] + [radiation_flux] type = RadiationEnergyFluxBC variable = T boundary = 'front' @@ -265,9 +265,9 @@ timestep=4e-7 sb_constant = 'sb_constant' absorptivity = 'abs' use_displaced_mesh = true - [../] + [] - [./weld_flux] + [weld_flux] type = GaussianWeldEnergyFluxBC variable = T boundary = 'front' @@ -278,104 +278,104 @@ timestep=4e-7 y_beam_coord = '2e-4 * sin(t * 2 * pi / ${period})' z_beam_coord = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_x] + [vapor_recoil_x] type = VaporRecoilPressureMomentumFluxBC variable = vel_x boundary = 'front' component = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_y] + [vapor_recoil_y] type = VaporRecoilPressureMomentumFluxBC variable = vel_y boundary = 'front' component = 1 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_z] + [vapor_recoil_z] type = VaporRecoilPressureMomentumFluxBC variable = vel_z boundary = 'front' component = 2 use_displaced_mesh = true - [../] + [] - [./surface_x] + [surface_x] type = SurfaceTensionBC variable = vel_x boundary = 'front' component = 0 use_displaced_mesh = true - [../] + [] - [./surface_y] + [surface_y] type = SurfaceTensionBC variable = vel_y boundary = 'front' component = 1 use_displaced_mesh = true - [../] + [] - [./surface_z] + [surface_z] type = SurfaceTensionBC variable = vel_z boundary = 'front' component = 2 use_displaced_mesh = true - [../] + [] - [./displace_x_top] + [displace_x_top] type = DisplaceBoundaryBC boundary = 'front' variable = 'disp_x' velocity = 'vel_x' - [../] - [./displace_y_top] + [] + [displace_y_top] type = DisplaceBoundaryBC boundary = 'front' variable = 'disp_y' velocity = 'vel_y' - [../] - [./displace_z_top] + [] + [displace_z_top] type = DisplaceBoundaryBC boundary = 'front' variable = 'disp_z' velocity = 'vel_z' - [../] + [] [] [ADMaterials] - [./kc_fits] + [kc_fits] type = DemonstrationPlantFits temperature = T beta = 1e7 - [../] - [./boundary] + [] + [boundary] type = DemonstrationPlantFitsBoundary use_displaced_mesh = true boundary = 'front' temperature = T - [../] + [] [] [Materials] - [./const] + [const] type = GenericConstantMaterial prop_names = 'abs sb_constant' prop_values = '1 5.67e-8' - [../] + [] [] [Preconditioning] - [./SMP] + [SMP] type = SMP full = true solve_type = 'NEWTON' - [../] + [] [] [Executioner] @@ -390,26 +390,26 @@ timestep=4e-7 line_search = 'none' nl_max_its = 12 l_max_its = 100 - [./TimeStepper] + [TimeStepper] type = IterationAdaptiveDT optimal_iterations = 7 dt = ${timestep} linear_iteration_ratio = 1e6 growth_factor = 1.5 - [../] + [] [] [Outputs] print_linear_residuals = false - [./exodus] + [exodus] type = Exodus output_material_properties = true show_material_properties = 'mu' - [../] - [./dofmap] + [] + [dofmap] type = DOFMap execute_on = 'initial' - [../] + [] checkpoint = true [] @@ -417,12 +417,11 @@ timestep=4e-7 show_var_residual_norms = true [] - [Adaptivity] marker = combo max_h_level = 3 - [./Indicators] + [Indicators] # [./error_x] # type = GradientJumpIndicator # variable = vel_x @@ -439,10 +438,10 @@ timestep=4e-7 # type = GradientJumpIndicator # variable = p # [../] - [./error_T] + [error_T] type = GradientJumpIndicator variable = T - [../] + [] # [./error_dispx] # type = GradientJumpIndicator # variable = disp_x @@ -451,13 +450,13 @@ timestep=4e-7 # type = GradientJumpIndicator # variable = disp_y # [../] - [./error_dispz] + [error_dispz] type = GradientJumpIndicator variable = disp_z - [../] - [../] + [] + [] - [./Markers] + [Markers] # [./errorfrac_x] # type = ErrorFractionMarker # refine = 0.7 @@ -482,12 +481,12 @@ timestep=4e-7 # coarsen = 0.3 # indicator = error_p # [../] - [./errorfrac_T] + [errorfrac_T] type = ErrorFractionMarker refine = 0.5 coarsen = 0.2 indicator = error_T - [../] + [] # [./errorfrac_dispx] # type = ErrorFractionMarker # refine = 0.7 @@ -500,36 +499,36 @@ timestep=4e-7 # coarsen = 0.3 # indicator = error_dispy # [../] - [./errorfrac_dispz] + [errorfrac_dispz] type = ErrorFractionMarker refine = 0.5 coarsen = 0.2 indicator = error_dispz - [../] - [./combo] + [] + [combo] type = ComboMarker markers = 'errorfrac_T errorfrac_dispz' - [../] - [../] + [] + [] [] [Postprocessors] - [./num_dofs] + [num_dofs] type = NumDOFs system = 'NL' - [../] - [./nl] + [] + [nl] type = NumNonlinearIterations - [../] - [./tot_nl] + [] + [tot_nl] type = CumulativeValuePostprocessor postprocessor = 'nl' - [../] - [./linear] + [] + [linear] type = NumLinearIterations - [../] - [./tot_linear] + [] + [tot_linear] type = CumulativeValuePostprocessor postprocessor = 'linear' - [../] + [] [] diff --git a/examples/laser_welding_problems/profiling/kc-profiling.i b/examples/laser_welding_problems/profiling/kc-profiling.i index 74b7821..854dfdb 100644 --- a/examples/laser_welding_problems/profiling/kc-profiling.i +++ b/examples/laser_welding_problems/profiling/kc-profiling.i @@ -24,62 +24,62 @@ [] [Variables] - [./vel_x] - [./InitialCondition] + [vel_x] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] + [] + [] - [./vel_y] - [./InitialCondition] + [vel_y] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] + [] + [] - [./vel_z] - [./InitialCondition] + [vel_z] + [InitialCondition] type = ConstantIC value = 1e-15 - [../] - [../] + [] + [] - [./T] - [./InitialCondition] + [T] + [InitialCondition] type = ConstantIC value = 300 - [../] - [../] - - [./p] - [../] - [./disp_x] - [../] - [./disp_y] - [../] - [./disp_z] - [../] + [] + [] + + [p] + [] + [disp_x] + [] + [disp_y] + [] + [disp_z] + [] [] [Kernels] - [./disp_x] + [disp_x] type = Diffusion variable = disp_x - [../] - [./disp_y] + [] + [disp_y] type = Diffusion variable = disp_y - [../] - [./disp_z] + [] + [disp_z] type = Diffusion variable = disp_z - [../] + [] [] [ADKernels] # mass - [./mass] + [mass] type = INSADMass variable = p u = vel_x @@ -87,17 +87,17 @@ w = vel_z p = p use_displaced_mesh = true - [../] + [] # x-momentum, time - [./x_momentum_time] + [x_momentum_time] type = INSADMomentumTimeDerivative variable = vel_x use_displaced_mesh = true - [../] + [] # x-momentum, space - [./x_momentum_space] + [x_momentum_space] type = INSADMomentumBase variable = vel_x u = vel_x @@ -106,17 +106,17 @@ p = p component = 0 use_displaced_mesh = true - [../] + [] # y-momentum, time - [./y_momentum_time] + [y_momentum_time] type = INSADMomentumTimeDerivative variable = vel_y use_displaced_mesh = true - [../] + [] # y-momentum, space - [./y_momentum_space] + [y_momentum_space] type = INSADMomentumBase variable = vel_y u = vel_x @@ -125,17 +125,17 @@ p = p component = 1 use_displaced_mesh = true - [../] + [] # z-momentum, time - [./z_momentum_time] + [z_momentum_time] type = INSADMomentumTimeDerivative variable = vel_z use_displaced_mesh = true - [../] + [] # z-momentum, space - [./z_momentum_space] + [z_momentum_space] type = INSADMomentumBase variable = vel_z u = vel_x @@ -144,102 +144,102 @@ p = p component = 2 use_displaced_mesh = true - [../] - - # temperature - [./temperature_time] - type = INSADTemperatureTimeDerivative - variable = T - use_displaced_mesh = true - [../] - - [./temperature_space] - type = INSADTemperature - variable = T - u = vel_x - v = vel_y - w = vel_z - p = p - use_displaced_mesh = true - [../] - - [./mesh_x] + [] + + # temperature + [temperature_time] + type = INSADTemperatureTimeDerivative + variable = T + use_displaced_mesh = true + [] + + [temperature_space] + type = INSADTemperature + variable = T + u = vel_x + v = vel_y + w = vel_z + p = p + use_displaced_mesh = true + [] + + [mesh_x] type = INSConvectedMesh variable = vel_x disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] - [./mesh_y] + [] + [mesh_y] type = INSConvectedMesh variable = vel_x disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] - [./mesh_z] + [] + [mesh_z] type = INSConvectedMesh variable = vel_x disp_x = disp_x disp_y = disp_y disp_z = disp_z use_displaced_mesh = true - [../] + [] [] [BCs] - [./x_no_disp] + [x_no_disp] type = DirichletBC variable = disp_x boundary = 'back' value = 0 - [../] - [./y_no_disp] + [] + [y_no_disp] type = DirichletBC variable = disp_y boundary = 'back' value = 0 - [../] - [./z_no_disp] + [] + [z_no_disp] type = DirichletBC variable = disp_z boundary = 'back' value = 0 - [../] + [] - [./x_no_slip] + [x_no_slip] type = DirichletBC variable = vel_x boundary = 'bottom right left top back' value = 0.0 - [../] + [] - [./y_no_slip] + [y_no_slip] type = DirichletBC variable = vel_y boundary = 'bottom right left top back' value = 0.0 - [../] + [] - [./z_no_slip] + [z_no_slip] type = DirichletBC variable = vel_z boundary = 'bottom right left top back' value = 0.0 - [../] + [] - [./T_cold] + [T_cold] type = DirichletBC variable = T boundary = 'back' value = 300 - [../] + [] [] [ADBCs] - [./radiation_flux] + [radiation_flux] type = RadiationEnergyFluxBC variable = T boundary = 'front' @@ -247,9 +247,9 @@ sb_constant = 'sb_constant' absorptivity = 'abs' use_displaced_mesh = true - [../] + [] - [./weld_flux] + [weld_flux] type = GaussianWeldEnergyFluxBC variable = T boundary = 'front' @@ -260,79 +260,79 @@ y_beam_coord = 0 z_beam_coord = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_x] + [vapor_recoil_x] type = VaporRecoilPressureMomentumFluxBC variable = vel_x boundary = 'front' component = 0 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_y] + [vapor_recoil_y] type = VaporRecoilPressureMomentumFluxBC variable = vel_y boundary = 'front' component = 1 use_displaced_mesh = true - [../] + [] - [./vapor_recoil_z] + [vapor_recoil_z] type = VaporRecoilPressureMomentumFluxBC variable = vel_z boundary = 'front' component = 2 use_displaced_mesh = true - [../] + [] - [./displace_x_top] + [displace_x_top] type = DisplaceBoundaryBC boundary = 'front' variable = 'disp_x' velocity = 'vel_x' - [../] - [./displace_y_top] + [] + [displace_y_top] type = DisplaceBoundaryBC boundary = 'front' variable = 'disp_y' velocity = 'vel_y' - [../] - [./displace_z_top] + [] + [displace_z_top] type = DisplaceBoundaryBC boundary = 'front' variable = 'disp_z' velocity = 'vel_z' - [../] + [] [] [ADMaterials] - [./kc_fits] + [kc_fits] type = DemonstrationPlantFits temperature = T - [../] - [./boundary] + [] + [boundary] type = DemonstrationPlantFitsBoundary use_displaced_mesh = true boundary = 'front' temperature = T - [../] + [] [] [Materials] - [./const] + [const] type = GenericConstantMaterial prop_names = 'abs sb_constant' prop_values = '1 5.67e-8' - [../] + [] [] [Preconditioning] - [./SMP] + [SMP] type = SMP full = true solve_type = 'NEWTON' - [../] + [] [] [Executioner] @@ -347,24 +347,24 @@ line_search = 'none' nl_max_its = 12 l_max_its = 100 - [./TimeStepper] + [TimeStepper] type = IterationAdaptiveDT optimal_iterations = 6 dt = 1e-5 linear_iteration_ratio = 1e6 - [../] + [] [] [Outputs] - [./exodus] + [exodus] type = Exodus output_material_properties = true show_material_properties = 'mu' - [../] - [./dofmap] + [] + [dofmap] type = DOFMap execute_on = 'initial' - [../] + [] checkpoint = true [] @@ -373,15 +373,15 @@ [] [Postprocessors] - [./num_dofs] + [num_dofs] type = NumDOFs system = 'NL' - [../] - [./nl] + [] + [nl] type = NumNonlinearIterations - [../] - [./tot_nl] + [] + [tot_nl] type = CumulativeValuePostprocessor postprocessor = 'nl' - [../] + [] []