We present a set of simulations of relativistic jets from accretion disk initial setup with numerical solutions of a system of general-relativistic magnetohydrodynamics (GRMHD) partial differential equations in a fixed black hole (BH) spacetime which is able to show substructures formations inside the jet as well as lobe formation on the jet head. For this, we used a central scheme of finite volume method without dimensional split and with no Riemann solvers namely the Nessyahu-Tadmor method. Thus, we were able to obtain stable numerical solutions with spurious oscillations under control and with no excessive numerical dissipation. Therefore, we developed some setups for initial conditions capable of simulating the formation of relativistic jets from the accretion disk falling onto central black hole until its ejection, both immersed in a magnetosphere. In our simulations, we were able to observe some substructure of a jet created from an accretion initial disk, namely, jet head, knots, cocoon, and lobe. Also, we present an explanation for cocoon formation and lobe formation. Each initial scenario was determined by ratio between disk density and magnetosphere density, showing that this relation is very important for the shape of the jet and its substructures.