MatForge — Open Knowledge Hub for Scientific Simulation and Modeling

Reading Time: 10 minutesUnit testing is non-negotiable for trustworthy scientific software. Unlike commercial applications, research code often lacks formal testing, leading to irreproducible results and wasted effort. This guide covers pytest strategies specifically for scientific Python projects: handling numerical precision with pytest.approx, isolating external dependencies with mocking, using fixtures and parametrization efficiently, and integrating tests into continuous integration […]

Reading Time: 11 minutesFiPy simulations can achieve 10x to 100x speedups by moving compute-intensive operations to the GPU using CuPy (drop-in NumPy replacement) or Numba (JIT compilation). CuPy excels at array operations and requires minimal code changes, while Numba accelerates Python loops and custom functions. However, GPU acceleration isn’t always beneficial—small problems, memory-bound operations, and complex data structures […]

Reading Time: 10 minutesFiPy does not have built-in adaptive mesh refinement (AMR). Current approaches involve external mesh generation with Gmsh (inefficient for dynamic problems), integration with dedicated AMR libraries like libmesh (architecturally challenging), or switching to alternative phase-field codes that support AMR natively (MOOSE, PRISMS-PF). AMR provides significant speedups (often 2–10×) for phase-field problems with localized interfaces, but […]

Reading Time: 7 minutesTL;DR Maxwell’s equations describe how electric and magnetic fields evolve and interact. While FiPy was designed for diffusion-type problems, you can simulate electromagnetic waves by treating Maxwell’s equations as a coupled system of transient hyperbolic PDEs. The key is coupling Faraday’s and Ampère’s laws using FiPy’s TransientTerm and custom curl implementations. However, FiPy has limitations […]

Reading Time: 9 minutesTL;DR Multi-physics simulations often require coupling multiple specialized solvers. preCICE is a mature, open-source coupling library that enables partitioned multi-physics simulations in Python. This tutorial shows how to couple FiPy with another solver using preCICE, covering installation, adapter implementation, configuration, and common pitfalls. You’ll learn when to use partitioned coupling, how to set up data […]

Reading Time: 8 minutesTL;DR Battery electrochemistry modeling using partial differential equations (PDEs) provides high-fidelity simulation of lithium-ion cell behavior. The Single Particle Model (SPM) offers computational efficiency for real-time applications, while the Doyle-Fuller-Newman (DFN) model captures full electrolyte dynamics for high-power scenarios. This guide covers the governing equations, implementation with FiPy, and validation techniques for both approaches. Introduction […]

Reading Time: 10 minutesHDF5 is the de facto standard for storing large-scale scientific simulation data. Its hierarchical structure, parallel I/O capabilities via MPI, and built-in compression make it ideal for high-performance computing environments. However, improper use—especially poor chunking choices and incorrect parallel access patterns—can lead to severe performance degradation or data corruption. This guide covers HDF5 architecture, parallel […]

Reading Time: 7 minutesTL;DR ParaView is the go-to open-source tool for visualizing materials science simulation data. Whether you’re working with phase-field models, molecular dynamics (LAMMPS), or electronic structure calculations (VASP), ParaView transforms raw numerical output into publication-ready 3D visuals. This guide covers data import, essential filters, rendering techniques, and export workflows specifically for materials science applications. Why ParaView […]

Reading Time: 9 minutesValidation and verification (V&V) are essential quality assurance processes for PDE simulations. Verification ensures your code solves the equations correctly (solving the equations right). Validation confirms your model accurately represents real-world physics (solving the right equations). A robust V&V framework combines code verification via methods like the Method of Manufactured Solutions, solution verification with mesh […]

Reading Time: 7 minutesEditor’s note: This restored technical overview has been reconstructed from archival project traces, package descriptions, and related scientific computing references to preserve the historical and workflow context of the Illuminator library. Illuminator was not built to make simulation output merely look better. It addressed a more practical problem: how to inspect, store, and move field […]