This paper discusses analytical and numerical issues related to elliptic equations with random coefficients which are generally nonlinear functions of white noise. Singularity issues are avoided by using the Itô-Skorohod calculus to interpret the interactions between the coefficients and the solution. The solution is constructed by means of the Wiener Chaos (Cameron-Martin) expansions. The existence and uniqueness of the solutions are established under rather weak assumptions, the main of which...

We present a numerical algorithm to solve the micromagnetic equations based on tangential-plane minimization for the magnetization update and a homothethic-layer decomposition of outer space for the computation of the demagnetization field. As a first application, detailed results on the flower-vortex transition in the cube of Micromagnetic Standard Problem number 3 are obtained, which confirm, with a different method, those already present in the literature, and validate our method and code. We...

We present a numerical algorithm to solve the micromagnetic equations based on tangential-plane minimization for the magnetization update and a homothethic-layer decomposition of outer space for the computation of the demagnetization field. As a first application, detailed results on the flower-vortex transition in the cube of Micromagnetic Standard Problem number 3 are obtained, which confirm, with a different method, those already present in the literature, and validate our method and...

We show that while Runge-Kutta methods cannot preserve polynomial invariants in general, they can preserve polynomials that are the energy invariant of canonical Hamiltonian systems.

We prove the existence of solutions to nonlinear parabolic problems of the following type: $$\left\{\begin{array}{cc}{\displaystyle \frac{\partial b\left(u\right)}{\partial t}}+A\left(u\right)=f+\mathrm{div}\left(\Theta (x;t;u)\right)\hfill & \text{in}Q,\hfill \\ u(x;t)=0\hfill & \text{on}\partial \Omega \times [0;T],\hfill \\ b\left(u\right)(t=0)=b\left(u_0\right)\hfill & \text{on}\Omega ,\hfill \end{array}\right.$$ where $b:ℝ\to ℝ$ is a strictly increasing function of class ${𝒞}^1$, the term $$A\left(u\right)=-\mathrm{div}\left(a\right(x,t,u,\nabla u\left)\right)$$ is an operator of Leray-Lions type which satisfies the classical Leray-Lions assumptions of Musielak type, $\Theta :\Omega \times [0;T]\times ℝ\to ℝ$ is a Carathéodory, noncoercive function which satisfies the following condition: ${sup}_{\left|s\right|\le k}\left|\Theta (·,·,s)\right|\in E_{\psi }\left(Q\right)$ for all $k>0$, where $\psi$ is the Musielak complementary function of $\Theta$, and the second term $f$ belongs to $L^1\left(Q\right)$.

For contractive interval functions $\left[g\right]$ we show that $\left[g\right]\left({\left[x\right]}_{ϵ}^{k_0}\right)\subseteq \int \left({\left[x\right]}_{ϵ}^{k_0}\right)$ results from the iterative process ${\left[x\right]}^{k+1}:=\left[g\right]\left({\left[x\right]}_{ϵ}^k\right)$ after finitely many iterations if one uses the epsilon-inflated vector ${\left[x\right]}_{ϵ}^k$ as input for $\left[g\right]$ instead of the original output vector ${\left[x\right]}^k$. Applying Brouwer’s fixed point theorem, zeros of various mathematical problems can be verified in this way.