In comparison with engineering problems, ore-forming problems commonly have both large length-scale and large time-scale characteristics. The length-scale of an ore-forming problem is commonly measured in either tens of kilometers or even hundreds of kilometers, while the time-scale of a typical ore-forming problem is often measured in several million years or even several tens of million years. In addition, most ore-forming problems are coupled across both multiple processes and multiple scales. Particularly, complicated geochemical reactions, such as redox reactions and fluid-rock reactions, are involved in almost all ore-forming systems. Because of these significant differences between engineering problems and ore-forming ones, existing computer programs and related algorithms, which are designed for solving engineering problems, cannot be directly used to solve ore-forming problems without any modification. For this reason, extensive and systematic research has been conducted, in recent years, to develop new algorithms for solving different aspects of ore-forming problems or to modify the existing computer programs, originally designed for solving engineering problems, so as to be suitable for solving ore-forming problems. The major achievements in this emerging field include the development of: (1) the progressive asymptotic approach procedure for simulating convective pore-fluid flow in ore-forming systems; (2) the term splitting algorithm for dealing with fluid-rock interaction problems that are closely associated with mineral dissolution and precipitation in hydrothermal systems; (3) the segregated algorithm for solving chemical-dissolution front instability problems in ore-forming systems of critical and supercritical Zhao numbers; (4) the decoupling procedure for simulating fluid mixing, heat transfer and non-equilibrium redox chemical reactions in ore-forming systems; (5) the equivalent source algorithm for simulating thermal and chemical effects of intruded magma solidification problems; (6) the particle simulation method for dealing with spontaneous crack generation problems in the brittle rocks during mechanical deformation processes; and (7) the porosity-gradient replacement approach for computational simulation of chemical-dissolution front propagation in fluid-saturated porous media including pore-fluid compressibility.