A numerical study of how to simulate heavy atmospheric particle collisions with a 38-m, 1.5-MW horizontal axis wind turbine blade is discussed. Two types of particles were considered, namely hailstones and rain drops. Computations were performed by using a two-dimensional inviscid flowfield solver along with a particle position predictor code. Three blade sections were considered: at 35% span and characterized by a DU 97-W-300 airfoil, at 70% span with a DU 96-W-212 airfoil, and at 90% span using a DU 96-W-180 airfoil. The three blade sections are constituted by 8-ply carbon/epoxy panels, coated with ultra-high molecular weight polyethylene (UHMWPE). Hailstone and raindrop simulations were performed to estimate the location of the striking occurrences and the blade surface area subject to damage. Results show that the impact locations along the blade are a function of airfoil angle of attack, local relative velocity, airfoil shape, aerodynamics and mass of the particle. Hailstones were found to collide on nearly every portion of the blade section along their trajectory due to their insensitivity to the blade flowfield. The damaged surface areas were found to be small when compared to the overall impingement surface, and most of delamination damage was localized on the blade leading edge. Moreover, panel delamination occurred for outboard sections, when r/R ≥ 0.90. The damage due to raindrops was divided in an erosive and a fatigue contribution due to the impact force. It was observed that the erosive damage follows the cubic power of the blade velocity, whereas the impact force follows the square power of the blade velocity. Moreover, it was seen that the rain drops are sensitive to the blade flowfield, due to shape modifications through the Weber number. In particular, a sensitive behavior of the damage with respect to the blade angle of attack was observed. Nomenclature a = axial induction factor A = particle reference area AK = particle nondimensional mass b C = coating constant related to the fatigue curve c = airfoil chord length C = speed of sound COE = Cost of Energy C d = airfoil drag coefficient C D = particle drag coefficient C l = airfoil lift coefficient d = particle diameter D = particle drag force E = erosion rate E D = particle damage efficiency E I = particle impact efficiency F imp = particle impact force FT E = failure threshold energy FTV = failure threshold velocity g = gravitational acceleration GAEP = Gross Annual Energy Production h = airfoil projected height perpendicular to freestream k = number of coating stress wave reflections m = particle mass n i = number of droplet impacts per site during incubation period P = impact pressure r/R = blade section radial location R D = damage surface ratio R I = impingement surface ratio Re = particle Reynolds number Re ∞ = freestream Reynolds number s = impact location in airfoil arc lengths S e f f ,C = effective coating strength s tot = airfoil total arc length t = time t/c = airfoil thickness-to-chord ratio U = chordwise flowfield velocity compone...