Groundwater Flow
Theory and Physics
Overview
Teacher! Today's topic is about groundwater flow, right? What is it like?
Groundwater flow simulation. Contaminant transport. Well pumping.
Now I understand what my senior meant when he said, "Make sure you do groundwater flow simulation properly."
Governing Equations
Discretization Methods
How do you actually solve this equation on a computer?
We use spatial discretization by the Finite Element Method (FEM). We assemble the element stiffness matrix and construct the global stiffness equation.
We perform 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. higher-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 you cut corners on the finite element method part, you'll pay for it later. I'll keep that in mind!
Implementation in Commercial Tools
So, what software can be used for groundwater flow?
| 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 |
| OpenFOAM | Open Source (OpenCFD/ESI, OpenFOAM Foundation) | Dictionary files (blockMeshDict, etc.), .foam |
Vendor History and Product Integration Background
Do the origins of each software have dramatic stories?
Ansys Fluent
Next is about Ansys Fluent. What's the content?
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 content?
Developed by CD-adapco. Acquired by Siemens in 2016 and integrated into the Simcenter brand. Features polyhedral meshes.
Current affiliation: Siemens Digital Industries Software
After hearing this, I finally understand why development is so important!
COMSOL Multiphysics
Please 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! Please tell me more.
File Formats and Interoperability
When converting models between different solvers, attention is needed to the correspondence of element types, compatibility of material models, and differences in the representation of loads and boundary conditions. Especially, higher-order elements and special elements (cohesive elements, user-defined elements, etc.) often cannot be directly converted between solvers.
I see... Formats may seem simple at first glance, but they are actually very deep.
Practical Considerations
Are there things like "field wisdom" that aren't in textbooks?
Verifying mesh convergence, validating the reasonableness of boundary conditions, and performing sensitivity analysis of material parameters are extremely important.
- Mesh dependency verification: Confirm convergence with at least 3 levels of mesh density.
- Boundary condition validity: Setting physically meaningful constraint conditions.
- Result verification: Comparison with theoretical solutions, experimental data, and known benchmark problems.
Related Topics
なった
詳しく
報告