Let G=(V,E) be a finite weighted graph, and Ω⊆V be a domain such that Ω^∘≠∅. In this paper, we study the following initial boundary problem for the non-homogenous wave equation ∂^2_t u(t,x) − Δ_Ω u(t,x) = f(t,x), (t,x)∈[0,∞)×Ω^∘ u(0,x)=g(x), x∈Ω^∘, ∂_t u(0,x)=h(x), x∈Ω^∘, u(t,x)=0, (t,x)∈[0,∞)×∂Ω, where Δ_Ω denotes the Dirichlet Laplacian on Ω^∘. Using Rothe's method, we prove that the above wave equation has a unique solution.

The sine-Gordon (SG) equation and perturbed nonlinear Schrdinger (NLS) equations are studied numerically for modeling the propagation of two space dimensional (2D) localized pulses (the so-called <i>light bullets</i>) in nonlinear dispersive optical media. We begin with the (2+1) SG equation obtained as an asymptotic reduction in the two level dissipationless MaxwellBloch system, followed by the review on the perturbed NLS equation in 2D for SG pulse envelopes, which is globally well posed and has all the relevant higher order terms to regularize the collapse of standard critical (cubic focusing) NLS. The perturbed NLS is approximated by truncating the nonlinearity into finite higher order terms undergoing focusingdefocusing cycles. Efficient semi-implicit sine pseudospectral discretizations for SG and perturbed NLS are proposed with rigorous error estimates. Numerical comparison results between light bullet

We consider an elliptic pseudo-differential equation with a highly oscillating linear potential modeled as a stationary ergodic random field. The random field is a function composed with a centered long-range correlated Gaussian process. In the limiting of vanishing correlation length, the heterogeneous solution converges to a deterministic solution obtained by averaging the random potential. We characterize the deterministic and stochastic correctors. With proper rescaling, the mean-zero stochastic corrector converges to a Gaussian random process in probability and weakly in the spatial variables. In addition, for two prototype equations involving the Laplacian and the fractional Laplacian operators, we prove that the limit holds in distribution in some Hilbert spaces. We also determine the size of the deterministic corrector when it is larger than the stochastic corrector. Depending on the correlation structure of the random field and on the singularities of the Green’s function, we show that either the deterministic or the random part of the corrector dominates.

We introduce two exponential-type integrators for the “good” Boussinesq equation. They are of orders one and two, respectively, and they require lower spatial regularity of the solution compared to classical exponential integrators. For the first integrator, we prove first-order convergence in H^r for solutions in H^{r+1} with r>1/2. This new integrator even converges (with lower order) in H^r for solutions in H^r, i.e., without any additional smoothness assumptions. For the second integrator, we prove second-order convergence in H^r for solutions in H^{r+3} with r>1/2 and convergence in L^2 for solutions in H^3. Numerical results are reported to illustrate the established error estimates. The experiments clearly demonstrate that the new exponential-type integrators are favorable over classical exponential integrators for initial data with low regularity.

No abstract uploaded!