By HeshamFS
End-to-end materials simulation workflows: plan, configure, run, verify, and validate PDE/atomistic simulations on HPC clusters with FAIR metadata, ontology mapping, and automated diagnostics for LAMMPS, VASP, QE, and MOOSE.
Perform spatial and temporal convergence analysis for solution verification — compute observed convergence orders from grid or timestep refinement studies, apply Richardson extrapolation to estimate discretization error, and calculate the Grid Convergence Index (GCI) per ASME V&V 20 standards. Use when verifying that a numerical solution converges at the expected rate, estimating the error on the finest mesh, checking whether grids are in the asymptotic range, or preparing formal verification reports, even if the user only asks "is my mesh fine enough" or "how accurate is my solution."
Select and apply numerical differentiation schemes for PDE and ODE discretization — generate finite-difference stencils at arbitrary order and accuracy, choose between central, upwind, compact (Pade), and spectral methods, handle boundary stencils, and estimate truncation error scaling. Use when discretizing spatial derivatives, picking a scheme for advection- or diffusion-dominated problems, building custom stencils for nonstandard operators, or comparing dispersion and dissipation properties of candidate schemes, even if the user just says "how do I approximate this derivative" or "my solution is too diffusive."
Select and configure linear solvers for Ax=b systems arising in numerical simulations — choose between direct (LU, Cholesky) and iterative (CG, GMRES, BiCGSTAB, MINRES) methods, analyze sparsity patterns and matrix conditioning, recommend preconditioners (AMG, ILU, IC), apply row/column scaling, and diagnose convergence stagnation from residual histories. Use when setting up a linear solve for FEM/FVM assembly, debugging slow or stalled Krylov iterations, choosing a preconditioner for SPD or nonsymmetric systems, or investigating ill-conditioning, even if the user only says "my solver is slow" or "GMRES won't converge."
Plan and evaluate mesh generation for numerical simulations — estimate grid resolution from physics scales (interface width, boundary layers, wavelengths), check aspect ratios and skewness against quality thresholds, choose between structured, unstructured, and adaptive mesh refinement strategies, and compute grid sizing for 1D/2D/3D domains. Use when setting up a new mesh, diagnosing poor solver convergence caused by mesh quality, deciding how many points to place across a phase-field interface or boundary layer, or preparing a mesh convergence study, even if the user only asks "what resolution do I need" or "why is my solver failing."
Select and configure nonlinear solvers for root-finding f(x)=0, optimization min F(x), and least-squares problems — choose among Newton, Newton-Krylov, quasi-Newton (BFGS, L-BFGS), Broyden, Anderson acceleration, and Levenberg-Marquardt methods, configure line search or trust-region globalization, diagnose convergence rate (quadratic, linear, stagnated), and assess Jacobian quality and conditioning. Use when a Newton solver converges slowly or diverges, choosing between line search and trust region, debugging nonlinear iteration failures in FEM or phase-field codes, or selecting a solver for large-scale unconstrained optimization, even if the user only says "my Newton iterations aren't converging."
Own this plugin?
Verify ownership to unlock analytics, metadata editing, and a verified badge. GitHub access is read-only (username + org membership).
Sign in to claimOwn this plugin?
Verify ownership to unlock analytics, metadata editing, and a verified badge. GitHub access is read-only (username + org membership).
Sign in to claimBased on adoption, maintenance, documentation, and repository signals. Not a security audit or endorsement.
Give your AI coding agent real expertise in numerical methods, simulation best practices, and computational materials science — so it stops guessing.
New to Agent Skills? A "skill" is a portable folder of instructions + small scripts that an AI coding agent discovers automatically and loads only when relevant. They follow the open Agent Skills standard and work across 20+ tools — Claude Code, Codex, Cursor, Antigravity, GitHub Copilot, and more. Nothing to wire up: drop them in and your agent gets smarter at the task.
Simulation engineers repeat the same guidance to AI agents constantly: "Check the CFL number before running," "Use Richardson extrapolation for grid convergence," "Exit code 2 means bad input." General-purpose agents lack the domain knowledge to run reliable numerical simulations without heavy hand-holding — so they pick wrong time steps, miss convergence checks, and misread solver failures.
This project packages that domain knowledge into skills the agent finds and runs on its own:
You: Check if dt=0.001 is stable for my advection problem with v=2.0 m/s and dx=0.01.
Agent: I'll use the numerical-stability skill to check this.
[runs cfl_checker.py --dx 0.01 --dt 0.001 --velocity 2.0 --json]
CFL number = 0.2 (limit: 1.0) -- stable.
Recommended max dt = 0.005 with safety factor 1.0.
No prompt engineering, no copy-pasting formulas. 24 skills span numerical methods, simulation workflows, HPC deployment, verification & validation, FAIR data, robustness, and materials ontologies.
Most agent-skill collections are unverified prompt text. These are measured:
sbatch spec, CMSO IRIs).How and why: evaluation methodology · harness.
The repo follows the open standard, so the ecosystem installers work out of the box:
npx skills add HeshamFS/materials-simulation-skills # any agent → ~/.agents/skills/
gh skill install HeshamFS/materials-simulation-skills --pin v1.0.0 # version-pinned
Or install a curated bundle into a specific agent with the bundled mss CLI:
mss bundles # list bundles
mss install --agent claude --bundle verification-and-validation
Claude Code users can also /plugin marketplace add HeshamFS/materials-simulation-skills.
Full per-agent guide → docs/INSTALL.md.
Once installed, just describe your task — the agent picks the right skill, runs the script, and interprets the result (as in the example above). You can also name a skill explicitly:
Use convergence-study to check if my mesh is in the asymptotic range:
h = 0.4, 0.2, 0.1 gave stress = 98.5, 99.6, 99.9 MPa.
Browse the full catalog in docs/SKILLS.md (or the machine-readable
skills_index.json).
Try it locally / develop:
git clone https://github.com/HeshamFS/materials-simulation-skills.git
cd materials-simulation-skills && pip install -e ".[dev]"
mss list # list skills
mss run numerical-stability cfl_checker -- --dx 0.01 --dt 0.001 --velocity 2.0 --json
mss eval # run the deterministic eval gate
python -m pytest tests/ # full test suite
24 skills across 8 categories (one-line summary; full tables in docs/SKILLS.md):
npx claudepluginhub heshamfs/materials-simulation-skills --plugin core-numericalJarvis-CD MCP - Pipeline Management for High-Performance Computing with comprehensive workflow operations
Skills for NVIDIAs ecosystem spans GPU acceleration, CUDA, AI agents, inference, robotics, Physical AI, Omniverse, and simulation. This plugin helps you understand the pieces, choose a path, validate your setup, and build practical NVIDIA-powered workflows.
Autonomous research loops with 10 commands. Generalizes Karpathy's autoresearch loop to any domain with mechanical evaluation, overnight persistence, and zero dependencies.
Three AI models, one synthesis — multi-model research workflow for scientific domains
Life sciences computational skills for scientific AI agents — 197 skills covering genomics, proteomics, drug discovery, biostatistics, scientific computing, and scientific writing
Commands for scenario simulation and decision modeling