Central Difference Scheme
Theory and Physics
Overview
Teacher! Today's topic is about the central difference scheme, right? What is it like?
Second-order accurate but unbounded. Used in LES. Stabilized by Rhie-Chow interpolation.
Governing Equations
I see. So, if I can handle the central difference scheme, I'm basically okay to start, right?
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) is directly linked to the trade-off between solution accuracy and computational cost.
Matrix Solution Algorithms
What exactly are matrix solution algorithms?
Solve the simultaneous equations by direct methods (LU decomposition, Cholesky decomposition) or iterative methods (CG method, GMRES method). Preconditioned iterative methods are effective for large-scale problems.
| 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 I cut corners on the finite element method part, I'll pay for it later. I'll keep that in mind!
Implementation in Commercial Tools
So, what software can I use to do the central difference scheme?
| Tool Name | Developer/Current | Main File Formats |
|---|---|---|
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Ansys CFX | Ansys Inc. | .cfx, .def, .res, .ccl |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| OpenFOAM | Open Source (OpenCFD/ESI, OpenFOAM Foundation) | Dictionary files (blockMeshDict, etc.), .foam |
Vendor Lineage and Product Integration History
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.
Ansys CFX
Please tell me about "Ansys CFX"!
CFX developed by AEA Technology (UK). Acquired by Ansys in 2003. Features a coupled solver.
Current Affiliation: Ansys Inc.
After hearing this, I finally understand why development is important!
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
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 regarding the correspondence of element types, compatibility of material models, and differences in the representation of loads and boundary conditions. Particularly, higher-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
Are there things like "field wisdom" that aren't in textbooks?
Verifying mesh convergence, validating the appropriateness 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
なった
詳しく
報告