– We model numerically the partial normal contact of two elastic rough surfaces with highly correlated asperities. Facing surfaces are unmated and described as self-affine with a Hurst exponent H . The numerical algorithm is based on Fourier acceleration and allows efficient simulation of very large systems. The force, F , versus contact area, A, characteristics follows the law F ∼ A in accorda...