MeltPoolDG::Simulation::StefansProblem1WithFlowAndHeat Namespace Reference

Namespaces
namespace	AnalyticalSolution

Classes
class	InitialValuesTemperature
class	SimulationStefansProblem1WithFlowAndHeat

Functions
	MELTPOOLDG_REGISTER_CASE (MeltPoolCase, SimulationStefansProblem1WithFlowAndHeat, "stefans_problem1_with_flow_and_heat", 1, double)
	MELTPOOLDG_REGISTER_CASE (MeltPoolCase, SimulationStefansProblem1WithFlowAndHeat, "stefans_problem1_with_flow_and_heat", 2, double)
	MELTPOOLDG_REGISTER_CASE (MeltPoolCase, SimulationStefansProblem1WithFlowAndHeat, "stefans_problem1_with_flow_and_heat", 3, double)

Variables
static constexpr double	x_min = 0.0
static constexpr double	y_min = 0.0
static constexpr double	y_max = 1.e-3

Detailed Description

This example is derived from

Hardt, S., and F. Wondra. "Evaporation model for interfacial flows based on a continuum-field representation of the source terms." Journal of Computational Physics 227.11 (2008): 5871-5895.

and represents the example denoted as "Stefan's Problem 1".

The parameters listed in the paper are:

domain size [0.0, 0.001] discretized in 1d by 1000 cells

Initial interface location: 1e-6 m

boiling temperature: 373.15 K temperature at the hot wall: 383.15 K

gas (vapor) phase: – density: 1 kg/m^3 – viscosity: 0.0001 Pa/s – thermal_conductivity: 1e-2 W/(mK) – specific_heat_capacity: 1000 J/(kgK)

liquid phase: – density: 1 kg/m^3 – viscosity: 0.01 Pa/s – thermal_conductivity: 1 W/(mK) (Note: thermal diffusivity of the liquid phase was increased by order of magnitudes to obtain a constant temperature in the liquid phase) – specific_heat_capacity: 1000 J/(kgK)

Enthalpy of evaporation: 10^6 J/kg

NOTE: Due to the equal densities in the two phases, no flow velocities will be induced.

NOTE: In the publication, they did not use the evaporative mass flux calculated according to Schrage's theory. We used the model by Hardt and Wondra and calibrated the evaporation coefficient.

Function Documentation

MELTPOOLDG_REGISTER_CASE() [1/3]

MeltPoolDG::Simulation::StefansProblem1WithFlowAndHeat::MELTPOOLDG_REGISTER_CASE	(	MeltPoolCase	,
		SimulationStefansProblem1WithFlowAndHeat	,
		"stefans_problem1_with_flow_and_heat"	,
		1	,
		double
	)

MELTPOOLDG_REGISTER_CASE() [2/3]

MeltPoolDG::Simulation::StefansProblem1WithFlowAndHeat::MELTPOOLDG_REGISTER_CASE	(	MeltPoolCase	,
		SimulationStefansProblem1WithFlowAndHeat	,
		"stefans_problem1_with_flow_and_heat"	,
		2	,
		double
	)

MELTPOOLDG_REGISTER_CASE() [3/3]

MeltPoolDG::Simulation::StefansProblem1WithFlowAndHeat::MELTPOOLDG_REGISTER_CASE	(	MeltPoolCase	,
		SimulationStefansProblem1WithFlowAndHeat	,
		"stefans_problem1_with_flow_and_heat"	,
		3	,
		double
	)

Variable Documentation

x_min

constexpr double MeltPoolDG::Simulation::StefansProblem1WithFlowAndHeat::x_min = 0.0

staticconstexpr

y_max

constexpr double MeltPoolDG::Simulation::StefansProblem1WithFlowAndHeat::y_max = 1.e-3

staticconstexpr

y_min

constexpr double MeltPoolDG::Simulation::StefansProblem1WithFlowAndHeat::y_min = 0.0

staticconstexpr