Validation data for a hybrid smoothed dissipative particle dynamics (SDPD) spatial stochastic simulation algorithm (sSSA) method.

We present the validation of the hybrid sSSA-SDPD method for advection-diffusion-reaction problems coupled to discrete biochemical systems, as presented in the publication "A hybrid smoothed dissipative particle dynamics (SDPD) spatial stochastic simulation algorithm (sSSA) for advection-diffusion-reaction problems" (Drawert et al., 2019). We validate 1D diffusion, and 2D diffusion cases against their analytical solutions. We present graphs and tables of data showing the error in the simulation method.

Specifications table

Value of the data

The data provide a detailed numerical validation of the sSSA-SDPD method compared to analytic solutions for 1D and 2D diffusion problems.

This data contain normalized average concentration profiles for 1D (transient) and 2D (steady-state) diffusion, obtained using SDPD and sSSA-SDPD methods, as well as their L2 errors, which extend the reader׳s understanding of the capabilities and limitations of the methods.

Data and methods employed can be used as a benchmark for future validations and comparisons with other spatial stochastic methods, as well as other particle-based methods such as dissipative particle dynamics (DPD).

Data

Data provided in this article are based on computations performed applying the hybrid SDPD-sSSA method [1] to 1D and 2D diffusion problems. The dataset includes: schematics of the validation cases tested (Fig. 1a and Fig. 2a); the normalized concentration profiles (Figs. 1b, 2c and d); the level set curves of the concentration field (Fig. 2b) and the L2 errors of the simulations, compared with the analytical results at different particle arrangements (Table 1, Table 2). Additionally, a brief description of the partial differential equations (PDEs) and their analytical solutions is provided. Routines to reproduce the data are freely available at the Git repository: https://github.com/briandrawert/hybrid_SSA_SDPD.

Fig. 1

(a) Schematic of diffusion in a slab, with initial and boundary conditions. (b) Comparison of profiles of the normalized concentration for selected times (bottom to top: and ).

Fig. 2

(a) Schematic of diffusion in a flat plate, with initial and boundary conditions. (b) Level set curves of the normalized concentration at steady-state, for particles. (c) Horizontal and (d) vertical centerline profiles of the normalized concentration at steady-state, for the case with for particles. SSA results were averaged over realizations.

Table 1

Results: validation of 1D transient diffusion.

Δx	N	ϵ^SDPDa	ϵ^SSAb	C_0
1/16	17	6.617×10^−4	8.514×10^−4	1.6×10^5
1/32	33	4.669×10^−4	5.278×10^−4	3.2×10^5
1/64	65	3.764×10^−4	8.070×10^−4	6.4×10^5

Error computed at steady-state, .

Error computed based at the mean profile of realizations at steady-state.

Table 2

Results: diffusion in a flat plate.

Δx	N	ϵ^SDPDa	ϵ^SSAb	C_0
1/16	17^2	2.337×10^−4	2.902×10^−4	2.890×10^5
1/32	33^2	5.818×10^−5	1.294×10^−4	1.089×10^6
1/64	65^2	1.438×10^−5	2.043×10^−4	4.225×10^6

Error computed at steady-state, .

Error computed based at the mean profile of realizations at steady-state.

Experimental design, materials, and methods

Validation for one-dimensional diffusion

The first validation case to be considered consists of the one-dimensional transient, isotropic diffusion problem. Dirichlet boundary conditions are imposed on both sides of the domain, and the initial concentration of the domain is kept as zero. The physical model of the problem is depicted in Fig. 1a. All the physical parameters are expressed in SI units.

The concentration of an arbitrary species evolves along time in a domain . The mathematical model is given by:with Dirichlet boundary conditions , and initial condition . The initial-boundary value problem can be solved analytically, using separation of variables. The solution has the form:

The problem is solved using both sSSA and SDPD diffusion formulations. Boundary conditions are modeled as layers of particles with fixed imposed concentrations. The mass diffusivity is kept constant, and a total of three levels of refinement were used, with the distance between particles and . A time step was used to ensure stability.

Due the stochastic nature of the sSSA algorithm, results were rendered as averages of multiple realizations in order to be compared with the SDPD and exact solutions. For all the three cases considered, realizations were performed. Hence, the sSSA results should be interpreted as first-order statistical moments, i.e., average and standard deviation.

Temporal profiles comparing sSSA and SDPD with the exact solution are shown in Fig. 1b.

Table 1 shows the quantitative results obtained. Errors and were computed using the norm of the difference between the present method and reference solutions, normalized by the number of particles in the domain, i.e.,where is a reference concentration, e.g., the exact solution, when available, is a normalization factor and N is the number of particles considered.

Validation for two-dimensional diffusion

The next validation case consists of a two-dimensional, steady-state, isotropic diffusion problem in a square region , with a non-homogeneous Dirichlet boundary condition in the upper wall, and homogeneous conditions on the remaining boundaries, as shown in Fig. 2a. The initial condition is homogeneous, i.e., , . We seek to find the steady-state by evolving the problem in time until it reaches the steady-state. The mathematical model is given bywhere the steady-state exact solution in known to have the form

Simulations were performed using and equally-spaced particles distributed in a squared domain of length . Dirichlet boundary conditions are modeled as layers of particles with fixed imposed concentrations, avoiding the truncation of the kernel in near-boundary regions. The mass diffusivity is kept constant, and a time step was used.

A qualitative plot of the level set curves of the normalized concentration is described by Fig. 2b. Figs. 2c–d show the horizontal and vertical centerline profiles of the normalized concentration at steady-state. Profiles of sSSA-resolved concentrations, along with the standard deviation bars were obtained as averages over realizations. A list of the estimated -norm errors is summarized in Table 2.

Subject area	Computational Physics, Mechanical Engineering
More specific subject area	SDPD, Particle-based fluid dynamics. Reaction-diffusion master equation, Discrete stochastic simulation.
Type of data	Table, graph, figure
How data was acquired	SDPD simulations were performed using LAMMPS. Validation data was generated using Matlab codes.
Data format	Raw and analyzed.
Experimental factors	The data was created by simulation using the SDPD method
Experimental features	The data shows diffusion dynamics in 1D and 2D domains
Data source location	The data was generated on computers at UCSB and UNCA
Data accessibility	Data can be reproduced using the Docker container, freely available at the Git repository: https://github.com/briandrawert/hybrid_SSA_SDPD
Related research article	B. Drawert, B. Jacob, Z. Li, T.-M. Yi, L. Petzold. A hybrid smoothed dissipative particle dynamics (SDPD) spatial stochastic simulation algorithm (sSSA) for advection-diffusion-reaction problems, Journal of Computational Physics[1]