Questions tagged [stochastic]
For questions regarding the numerical treatment of processes whose behaviors are determined by both deterministic (predictable) and non-deterministic (random) actions.
61
questions
0
votes
0
answers
55
views
Why is Magnetic Susceptibility not showing the expected transition in 2D&3D Ising Model?
I'm trying to code the Ising Model with the metropolis algorithm to study the ferromagnetic-paramagnetic transitions. The code seems to work ; the equilibration happens. While equilibrating, the ...
1
vote
2
answers
263
views
Simulate Jump-Diffusion $dX_t = \mu(X_t)dt + \sigma(X_t)dW(t) + \int_{\{|c| <1 \}}g(X_t,c)\tilde{N}(dt,dc) + \int_{\{|c| \ge 1 \}}h(X_t,c)N(dt,dc)$
I would like to be able to model an SDE having the form
$$dX_t = \mu(X_t)dt + \sigma(X_t)dW(t) + \int_{\{|c| <1 \}}g(X_t,c)\tilde{N}(dt,dc) + \int_{\{|c| \ge 1 \}}h(X_t,c)N(dt,dc).$$
where $W$ is a ...
4
votes
1
answer
112
views
Stochastic cellular automata - algorithm limited by 1 cell per timestep
Context
Let's say I am trying to model the spread of mold in a petri dish, using a stochastic cellular automata approach. The petri dish can be thought of as a grid of 1mm x 1mm squares, each called ...
1
vote
0
answers
115
views
Stochastic differential equation system (SDE) : overflow encountered in double-scalars
I'm trying to integrate the following SDE system from Dekker et al. [1]
$$\begin{cases}
\frac{dx}{dt}=a_1x^3+a_2x+\phi+\zeta_x\\
\frac{dy}{dt}=b_1z+b_2(\kappa(x)-(y^2+z^2))y+\zeta_y\\
\frac{dz}{dt}=...
1
vote
1
answer
83
views
Rate of convergence - Stochastic Euler Method
The absolute error criterion of the pathwise approximation of an Ito process $X$ by an Euler approximation $Y$ is:
$$
\epsilon=E\left(\left|X_{T}-Y(T)\right|\right)
$$
We shall say that a time-...
1
vote
1
answer
98
views
Numerical stability of taking the `mean` of outputs from the simulation of a discrete stochastic dynamical system
I am writing a simulation for a discrete stochastic dynamical system. Since the simulation is stochastic, I need to run the simulation multiple times and then average the values of each timestep. I ...
4
votes
1
answer
51
views
Guidelines for publishing data from a stochastic simulation
So, my question is if one should ideally keep a record of all seeds that are used when publishing numerical work that involves one or more random number generators (e.g. a stochastic simulation), and ...
2
votes
1
answer
140
views
Adaptive Runge-Kutta for Stochastic (Projected) Gross-Pitaevskii Equation
I am using the XMDS library for solving the stochastic (projected) Gross-Pitaevskii equation
$$i \hbar \partial \Phi\left(\mathbf{r},t\right)_t=\hat{\mathcal{P}}\left\{(1-i \gamma)\left(\hat{H}_{\...
0
votes
1
answer
186
views
Numerically solving nonlinear parabolic stochastic PDEs
For a project I'm doing, I have to numerically solve a nonlinear parabolic stochastic partial differential equation, of the form
$$
u_t = u_{xx} + f(u)(u_x)^2 + a(u) + b(u)W(t, x),
$$
where primes ...
1
vote
0
answers
30
views
Effective fitness formula in the Moran process on a game
I was recently reading some literature about evolutionary game theory and I got particularly interested in Moran process linked to prisoner's dilemma as a model of evolution of two sub-populations.
...
5
votes
0
answers
121
views
Probability approximation: monte carlo VS sde
I have a probability measure $\mu$ (say, in $\mathbb{R}^{d}$, with density) and I want to approximate it numerically. Today I noticed that my measure is ergotic for a certain Stochastic Differential ...
2
votes
0
answers
146
views
Numerical integration of SDE: choice of $dt$ and algorithm
I am working on the following Stochastic Differential Equation (SDE) in the Quantum Mechanics context:
$$dX_{t} = a X_{t} dt + b X_{t} dW$$
where $X_{t}$ is my stochastic varible, $dt$ is my ...
1
vote
1
answer
66
views
What is the right way to set up two random tensor fields which have an identical average diffusivity
I want to compare some properties of traveling waves through two randomly diffusive media. The traveling waves follow the fisher equation:
$$\frac{du}{dt} = \nabla(\mathbf{D}_{\gamma} \nabla u) + u(1-...
3
votes
2
answers
226
views
How to optimize sampling for parameter estimation
I have a computer model with a number of parameters that need to be calibrated based on experimental results. It's also important to understand the sensitivity of the results to each parameter ...
2
votes
1
answer
117
views
How to get started with numerically solving a Stochastic Navier Stokes equation
I originally posted the question on the math stackexchange, and was told I should try here.
I’m researching Stochastic PDE, in particular the Navier Stokes Equation, and would like to estimate the ...