We introduce an unconditionally stable finite element (FE) method, the automatic variationally stable FE (AVS-FE) method for the numerical analysis of the Korteweg-de Vries (KdV) equation. The AVS-FE method is a Petrov-Galerkin method which employs the concept of optimal discontinuous test functions of the discontinuous Petrov-Galerkin (DPG) method. However, since AVS-FE method is a minimum residual method, we establish a global saddle point system instead of computing optimal test functions element-by-element. This system allows us to seek both the approximate solution of the KdV initial boundary value problem (IBVP) and a Riesz representer of the approximation error. The AVS-FE method distinguishes itself from other minimum residual methods by using globally continuous Hilbert spaces, such as H 1 , while at the same time using broken Hilbert spaces for the test. Consequently, the AVS-FE approximations are classical C 0 continuous FE solutions. The unconditional stability of this method allows us to solve the KdV equation space and time without having to satisfy a CFL condition. We present several numerical verifications for both linear and nonlinear versions of the KdV equation leading to optimal convergence behavior. Finally, we present a numerical verification of adaptive mesh refinements in both space and time for the nonlinear KdV equation.