Velocity Vector Field Reconstruction on Unstructured Grids
Find Innovative SolutionsGenerate Solutions
Solution Overview
Problem
Existing methods for simulating fluid flow in porous media struggle to accurately convert fluxes into a consistent velocity vector field, especially on unstructured grids, which is crucial for validating simulation codes and evaluating grid quality, due to the complexity of comparing fluxes across different grids and the need for geometric information.
Innovation Solution
A method is developed to reconstruct a velocity vector field from given fluxes on unstructured grids by subdividing cells into triangular prisms, using a divergence-free correction method, and solving a steady-state pressure equation with Mixed Finite Element Method (MFE) to minimize complementary energy, ensuring mass conservation and consistency with fluxes.
Engineering Contradictions & Design Principles
Engineering Contradiction Analysis
1Measurement precision
If fluxes are converted to velocity vector field on unstructured grids using conventional methods, then the conversion can be performed, but the accuracy and consistency of the velocity field with respect to flux constraints is poor
Solution Approach 1:
The method segments each unstructured grid cell into multiple subcells (tetrahedra or triangular prisms) that share common edges. This segmentation allows the velocity field to be constructed by solving flux balance equations on shared edges, ensuring local mass conservation and consistency with the given fluxes while maintaining computational tractability through the structured subcell decomposition.
Solution Approach 2:
The invention changes the mathematical formulation from directly computing velocity components to solving for edge fluxes first, then deriving velocity from those fluxes. This parameter transformation ensures that the velocity field automatically satisfies the flux constraints by construction, improving accuracy without proportionally increasing computational complexity.
2Measurement precision
If the grid resolution is increased to improve simulation accuracy, then the precision of flux and velocity computation improves, but the computational cost increases significantly
Solution Approach 1:
The method enables the velocity field computation to serve multiple purposes simultaneously: it satisfies flux constraints, ensures mass conservation, provides streamline information, and validates grid quality. This self-service capability reduces the need for separate validation and post-processing computations, lowering overall computational cost while maintaining high resolution.
Solution Approach 2:
The invention computes the velocity field locally on each cell using a simplified subcell decomposition approach rather than solving the full global system. This partial action approach provides sufficient accuracy for most applications without the excessive computational cost of global high-resolution solutions, enabling resolution increases without proportional cost increases.
3Reliability
If a velocity vector field is computed to ensure mass conservation and flux consistency, then the physical realism improves, but the computational complexity and time increase
Solution Approach 1:
The method performs preliminary decomposition of cells into subcells and pre-computation of geometric parameters (edge lengths, areas, normals) before solving the flux balance equations. This preliminary action ensures that the main computation step efficiently enforces mass conservation and flux consistency without iterative corrections, reducing computation time while maintaining reliability.
Data Source
AI summary
A method for constructing a velocity vector field from a grid and a set of fluxes for each face of the grid cells. The cells are first subdivided and internal fluxes are calculated for each cell subject to the constraints of the flux for each cell and to achieve the minimum energy state for the each cell. The minimum energy state is computed efficiently using a divergence-free correction method without introducing a pressure variable. Then, the velocity vector field is constructed from the subcell fluxes using mixed finite element interpolation.