In this contribution, we present the results of a model for the simulation of three dimensional unsteady nonlinear water waves. The model, designed to evaluate the wave drag of arbitrarily shaped ship hulls, is based on potential flow theory. The governing Laplace equation is complemented with fully nonlinear boundary conditions on the water free surface. The spatial discretization of the resulting time dependent boundary value problem is carried out by means of a collocated and iso-parametric Boundary Element Method (BEM) implemented making use of the open source library deal.II. The time integration of the nonlinear differential-algebraic system resulting from the spatial discretization is performed by means of an implicit Backward Differentiation Formula method implemented in the open source library Sundials . The classes of the OpenCASCADE library have been used to interface the model with CAD data structures. At the start of each simulation, the CAD surface describing the hull is imported and used to automatically generate the unstructured and non conformal quadrilateral grid needed by the BEM solver . The model also accounts for hulls with a transom stern. In such case, a specific treatment of the free surface nodes on the stern edge has been implemented. At low speeds, when the transom stern is partially immersed, a pressure patch is applied on the water surface detaching from the transom stern, to recover the gravity effect of the recirculating water on the underlying irrotational flow domain. The pressure patch is computed making use of experimental correlations obtained by Doctors et al. The test cases considered are that of the the US Navy Combatant DTMB-5415 and the National Physical Laboratory (NPL) hull. Numerical results are compared with experiments available in the literature.