Let Ω be either R n or a strongly Lipschitz domain of R n , and ω ∈ A ∞ (R n) (the class of Muckenhoupt weights). Let L be a second-order divergence form elliptic operator on L 2 (Ω) with the Dirichlet or Neumann boundary condition, and assume that the heat semigroup generated by L has the Gaussian property (G 1) with the regularity of their kernels measured by μ ∈ (0, 1]. Let Φ be a continuous, strictly increasing, subadditive, positive and concave function on (0, ∞) of critical lower type index p − Φ ∈ (0, 1]. In this paper, the authors first introduce the "geometrical" weighted local Orlicz-Hardy spaces h Φ ω, r (Ω) and h Φ ω, z (Ω) via the weighted local Orlicz-Hardy spaces h Φ ω (R n), and obtain their two equivalent characterizations in terms of the nontangential maximal function and the Lusin area function associated with the heat semigroup generated by L when p − Φ ∈ (n/(n + μ), 1]. Second, the authors furthermore establish three equivalent characterizations of h Φ ω, r (Ω) in terms of the grand maximal function, the radial maximal function and the atomic decomposition when the complement of Ω is unbounded and p − Φ ∈ (0, 1]. Third, as applications, the authors prove that the operators ∇ 2 G D are bounded from h Φ ω, r (Ω) to the weighted Orlicz space L Φ ω (Ω), and from h Φ ω, r (Ω) to itself when Ω is a bounded semiconvex domain in R n and p − Φ ∈ (n n+1 , 1], and the operators ∇ 2 G N are bounded from h Φ ω, z (Ω) to L Φ ω (Ω), and from h Φ ω, z (Ω) to h Φ ω, r (Ω) when Ω is a bounded convex domain in R n and p − Φ ∈ (n n+1 , 1], where G D and G N denote, respectively, the Dirichlet Green operator and the Neumann Green operator.