Nanotechnology is rapidly developing as an enabling technology for introducing new highly performant theragnostic tools for biomedical applications. Recently, plasmonic nanostructures (PNs) have shown a great potential for providing an effective method in this field of research and applications. PNs present very peculiar optical properties related to the presence of a large optical absorption cross-section due to the collective oscillation of the conduction electrons by incident laser field known as localized surface plasmon resonance (LSPR). The PNs are also capable of enhancing the incident laser field in the form of scattered near field for even an off-resonance irradiation where the incident laser wavelength is far away from LSPR. If the intensity of the incident laser is large enough, this field enhancement capability can be employed to generate a highly localized free electron (plasma) in the vicinity of the particles. The energy deposition in generated plasma can be considered as an energy source for photothermal ablation therapy, photoacoustic tomography, and light-controlled drug release. In this presentation we will present a model that will serve as a design tool for simulating ultrashort pulsed laser interaction with PNs in water for such applications. This tool is built based on the following physical phenomena that couples: 1) optical properties of the PN obtained from simulation of electromagnetic wave interaction with PN and coupled through resistive heating due to laser pulse to; 2) two temperature model to analyze transient heat transfer in ultrafast regime and temperature increase calculation, where temperature is used as a coupling parameter to; 3) plasma dynamics and 4) structural mechanics analysis for linear thermal expansion, stress and strain calculation for the input as a boundary condition in acoustic-structure interaction for; 5) acoustic pressure wave propagation modeling.