GeoDict Contact Correction GeoApp

The Contact Correction GeoApp is a specialized tool developed for granular digital materials. It identifies and corrects false or oversized particle contact regions in 3D segmented structures caused by insufficient image resolution. Using two methods—Area‑Dependent Stiffness Reduction and Explicit Contact Area Reduction—the GeoApp improves the accuracy of mechanical property predictions for rocks, battery electrodes, fuel cell electrodes, ceramics, sintered materials, particle packings, and powder materials. It provides more reliable digital models for material analysis, simulation, and optimization driven by real microstructures.

Why Contact Correction?

In particle‑based microstructures, contact conditions between particles directly affect the overall mechanical response. When image resolution is insufficient:

Thin gaps between particles may not be resolved

Adjacent particles may be incorrectly connected

Contact areas are enlarged

Local structures exhibit excessive stiffness

Porosity and local volume distribution deviate

Computed elastic moduli and acoustic velocities may be overestimated

The Contact Correction GeoApp addresses these issues by selectively correcting contact regions in segmented structures, bringing digital models closer to the true contact state of the material. The method has been validated through digital rock physics experiments and is applicable to a wide range of granular materials, including battery electrodes, fuel cell electrodes, catalyst supports, ceramics, sintered materials, particle packings, and powder materials.

Particle & Contact Region Identification Based on GrainFind

Both correction methods build on particle identification results from GrainFind. GrainFind uses a watershed‑based segmentation algorithm to identify individual particles, their spatial positions, and contact regions between adjacent particles.

Before performing contact correction, users can set Reconnect Fragmented Grains – Interface Threshold to control whether contacting particle fragments are merged:

0%: All contacting particle fragments are merged

100%: No particle fragments are merged

Method 1: Area‑Dependent Stiffness Reduction

This method does not alter the original segmented structure or the geometry of particle contact regions. Instead, it adjusts the effective mechanical properties of contact regions based on the relative size of the contact area.

Users can set the Reduce Contact Properties by parameter to define the reduction ratio for the following material properties in contact regions: density, Young's modulus, and Poisson's ratio.

Key features:

Preserves the original voxel‑based segmented structure

Does not change the original porosity

Modifies local mechanical properties based on contact area

Suitable for workflows where geometric changes are not desired

After correction, each contact region receives an independent material ID and corresponding mechanical properties, determined by the relative area of that contact.

Method 2: Explicit Contact Area Reduction

This method directly erodes contact‑region voxels, geometrically reducing the contact area between particles.

Users set the Eroded Voxel Number to define how many voxel layers are removed from the contact region. This parameter determines how far the particle contact area is eroded toward the pore space.

Key features:

Directly modifies contact region geometry

Reduces the actual connection area between particles

Changes the porosity of the original structure

Brings resolution‑limited images closer to experimental conditions

Suitable for simulations requiring explicit geometric correction of particle contacts

After processing, all identified and reduced contact regions are assigned new material IDs. Their mechanical properties can remain consistent with the original particle material, while the contact area is reduced by the specified number of voxel layers.

How to Choose Between the Two Methods?

The choice between the two correction methods depends on your specific simulation goals and whether you want to preserve the original segmented geometry.

Area‑Dependent Stiffness Reduction does not change the geometry or the porosity of the structure. Instead, it reduces the local effective stiffness in contact regions based on the relative area of each contact. This method is ideal when you wish to keep the original voxel model intact and only adjust mechanical properties to account for overestimated contacts.

Explicit Contact Area Reduction directly erodes voxels in the contact regions, thereby reducing the physical contact area between particles. This changes both the geometry and the porosity of the structure, making the model closer to what would be observed at higher resolution. This method is suitable when you need a geometric correction that affects subsequent simulations of flow, transport, or mechanical response where contact area plays a critical role.

In summary, choose Area‑Dependent when preserving the original segmentation and porosity is important, and choose Explicit when you want to physically reduce false connections and allow porosity to change accordingly.

Application Examples

Digital Rock Physics

Contact Correction is used to correct false contacts in granular rock models such as sandstones and carbonates, improving prediction accuracy for: elastic moduli, rock compressibility, P‑wave and S‑wave velocities, Vp/Vs ratio, and acoustic and rock mechanical parameters.

Batteries & Fuel Cells

The tool is applied to study how contact conditions within particle electrodes affect: overall electrode stiffness, load transfer between particles, electrode compression deformation, stress and swelling during charge/discharge, and particle fracture and structural degradation.

Ceramics & Sintered Materials

It corrects oversized particle connection regions in segmented models for analyzing: sintered neck structures, local stress concentrations, and overall elastic and deformation behavior.

Powders & Particle Packings

It enables study of how particle size, shape, packing state, and contact relationships influence macroscopic mechanical properties.

Required GeoDict Licenses

According to Math2Market's official configuration, running the Contact Correction GeoApp requires:

Essential Modules: GrainFind, MatDict

Essential GeoApp: Acoustic Rock Properties

Commonly Used Supporting Modules: ImportGeo‑Vol (for importing, processing, and segmenting 3D images such as µCT), GrainGeo (for generating digital models of granular materials), ElastoDict (for computing mechanical properties, stress, strain, and deformation)

The specific module combination depends on the actual material and research objectives.

Value of Contact Correction

The Contact Correction GeoApp establishes a critical calibration step between 3D imaging and mechanical simulation:

µCT 3D Imaging → Image Processing & Segmentation → Particle & Contact Region Identification → Contact Correction → Mechanical & Acoustic Property Simulation → Comparison & Validation with Experimental Data

By correcting particle contact artifacts caused by limited image resolution, users can:

Reduce the risk of overestimating material stiffness

Improve the physical realism of digital material models

Enhance predictions of mechanical and acoustic parameters

Enable more reliable comparisons between different microstructure designs

Reduce reliance on repeated experiments and physical sample testing