Professor! Today's topic is the Stefan problem, right? What exactly is it?

It's a moving boundary problem involving solidification and melting. It's a classical formulation for tracking the temporal evolution of the interface position.

How do we actually solve this equation computationally?

We use spatial discretization via the Finite Element Method (FEM). This involves assembling the element stiffness matrices and constructing the global stiffness equation.

What specifically does the matrix solution algorithm entail?

We solve the system of equations using either direct methods (LU decomposition, Cholesky decomposition) or iterative methods (Conjugate Gradient method, GMRES method). For large-scale problems, preconditioned iterative methods are particularly effective.

Stefan problem

Category: 熱解析 | Integrated 2026-04-06

ステファン問題

Theory and Physics

Overview

🧑‍🎓

Professor! Today's topic is about the Stefan problem, right? What is it?

🎓

It's a moving boundary problem for solidification and melting. A classical formulation for tracking the time evolution of the interface position.

🧑‍🎓

Wait, wait, a moving boundary for solidification and melting... so can it also be used for cases like this?

Governing Equations

$$ k_s\frac{\partial T_s}{\partial x}\bigg|_{x=s} - k_l\frac{\partial T_l}{\partial x}\bigg|_{x=s} = \rho L_f \frac{ds}{dt} $$

$$ s(t) = 2\lambda\sqrt{\alpha t} $$

🧑‍🎓

I see. So if we can describe the Stefan problem, we're basically good to start, right?

Discretization Methods

🧑‍🎓

How do we actually solve these equations on a computer?

🎓

We use spatial discretization by the Finite Element Method (FEM). We assemble the element stiffness matrices and construct the global stiffness equation.

🎓

We perform a transformation to the weak form (variational form) and use formulation by the Galerkin method using test functions and shape functions. The choice of element type (low-order elements vs. high-order elements, full integration vs. reduced integration) directly affects the trade-off between solution accuracy and computational cost.

Matrix Solution Algorithms

🧑‍🎓

What exactly are matrix solution algorithms?

🎓

We solve the simultaneous equations using direct methods (LU decomposition, Cholesky decomposition) or iterative methods (CG method, GMRES method). For large-scale problems, preconditioned iterative methods are effective.

Solver	Classification	Memory Usage	Applicable Scale
LU Decomposition	Direct Method	O(n²)	Small to Medium Scale
Cholesky Decomposition	Direct Method (Symmetric Positive Definite)	O(n²)	Small to Medium Scale
PCG Method	Iterative Method	O(n)	Large Scale
GMRES Method	Iterative Method	O(n·m)	Large Scale / Non-symmetric
AMG Preconditioner	Preprocessing	O(n)	Very Large Scale

🧑‍🎓

So if we cut corners on the finite element method part, we'll pay for it later. I'll keep that in mind!

Implementation in Commercial Tools

🧑‍🎓

So, what software can be used to tackle the Stefan problem?

Tool Name	Developer/Current	Main File Formats
Ansys Fluent	Ansys Inc.	.cas, .dat, .msh, .jou
Simcenter STAR-CCM+	Siemens Digital Industries Software	.sim, .java, .csv
COMSOL Multiphysics	COMSOL AB	.mph
Abaqus FEA (SIMULIA)	Dassault Systèmes SIMULIA	.inp, .odb, .cae, .sta, .msg

Vendor History and Product Integration Background

🧑‍🎓

Are the origins of each software quite dramatic?

Ansys Fluent

🧑‍🎓

Next is about Ansys Fluent. What's the story?

🎓

Developed by Fluent Inc. Acquired by Ansys in 2006. A general-purpose CFD solver based on unstructured grids.

Current affiliation: Ansys Inc.

Simcenter STAR-CCM+

🧑‍🎓

Next is about Simcenter STAR. What's the story?

🎓

Developed by CD-adapco. Acquired by Siemens in 2016 and integrated into the Simcenter brand. Known for polyhedral meshes.

Current affiliation: Siemens Digital Industries Software

🧑‍🎓

After hearing this, I finally understand why development is so important!

COMSOL Multiphysics

🧑‍🎓

Tell me about "COMSOL Multiphysics"!

🎓

Founded in Sweden in 1986. Started as FEMLAB with MATLAB integration, later renamed COMSOL. Strong in multiphysics.

Current affiliation: COMSOL AB

🧑‍🎓

Wow, the story of development is super interesting! Tell me more.

File Formats and Interoperability

🧑‍🎓

Are there any points to note when transferring data between different software?

Format	Extension	Type	Overview
STEP	.stp/.step	Neutral CAD	ISO 10303 compliant 3D CAD data exchange format. Supports geometry + PMI.
CGNS	.cgns	CFD Data	CFD General Notation System. Standard exchange format for CFD results.
VTK	.vtk/.vtu	Visualization	Visualization Toolkit format. Used by ParaView, etc.

🎓

When converting models between different solvers, attention is needed regarding the correspondence of element types, compatibility of material models, and differences in the representation of loads and boundary conditions. Especially high-order elements and special elements (cohesive elements, user-defined elements, etc.) often cannot be directly converted between solvers.

🧑‍🎓

I see... formats seem simple at first glance, but they're actually quite deep, aren't they?

Practical Considerations

🧑‍🎓