- Two dimensional Melting
- Milk-crown
- Collision/Destruction
- Bursting Droplets
- Condensation
- Bubble Nucleation
- Ostwald Ripening
- Self-avoiding Random Walk
- Reaction-Diffusion Systems
- Structural Analysis

Red particles are disordered and blue ones are oredred.

Fluid | Critical | Sodid |

- Fluid (monochrome) [mpg, 4.0MB] / Fluid (color) [mpg, 3.5MB]
- Critical (monochrome) [mpg, 3.5MB] / Critical (color) [mpg, 3.4MB]
- Solid (monochrome) [mpg, 3.2MB] / Solid (color) [mpg, 3.4MB]

Simulating of Milk-crown with Lennard-Jones particles.

Milk-crown [mpg, 2.9MB]Collision of two Lennard-Jones solids. Two-dimensional Lennard-Jones system. MPI parallelized (4 PU). Different colors denote different processes..

- Offset Collision [mpg, 2.4MB] 4PU Parallelized.
- Onset Collision [mpg, 2.4MB] 8PU Parallelized.

Collision of two Lennard-Jones crystals (FCC). Three-dimensional Lennard-Jones system.

- Onset Collision [mpg, 1.2MB]
- Onset Collision [mpg, 1.2MB]
- Offset Collision [mpg, 0.9MB]
- Offset Collision [mpg, 2.4MB]

- Low Speed Collision [mpg, 1.2MB]
- Middle Speed Collision [mpg, 1.2MB]
- Ditto (different angle) [mpg, 1.2MB]
- High Speed Collision [mpg, 0.9MB]

Collision of two Lennard-Jones liquid droplets. 107294 particles, 250000 steps with the time step 0.002.

Collision of two Lennard-Jones droplets. 33454 particles, 16000 steps with the time step 0.005. Images are made with POV-Ray.

Bursting of droplets. This work is inspired by Edward Horsford's artworks.

187634 particles, 8000 steps with the time step 0.005. Truncation length is 2.5. Initial density and temperature of droplet are 0.7 and 0.9, respecrively. The radius of the droplet is 40 and the system size is 400 in LJ units.

268307 particles, 8000 steps with the time step 0.005. Truncation length is 2.5. Initial density and temperature of droplet are 1.0 and 1.5, respecrively. The radius of the droplet is 40 and the system size is 400 in LJ units.

The result of molecular dynamics simulations with 300,000 Lennard-Jones particles. At first, all particles are heated enough to be gas. After thermalization, temperature is decreased quickly below the boiling point, and then, hot steam is condensed to liquid.

The result of molecular dynamics simulations with 1,468,006 Lennard-Jones particles. At first, the pure-liquid phase is achieved. After thermalization, the pressure is reduced below the boiling point, and then, a bubble starts to grow after some waiting time.

While time evolutions of the two movies are identical, different values of density threshold are used to indentify the gas state.

The result of molecular dynamics simulations with 9,663,597 Lennard-Jones particles. At first, the pure-liquid phase is achieved. After thermalization, the pressure is reduced below the boiling point. The system is in the spinodal region after the depressurization, and therefore, multi-bubble nucleation is observed.

After multi-bubble nucleation, the ostwald ripening happens because of the interaction between bubbles, and finally, the single bubble remains.

Result 1

ă€€Simulation Details: Time = 1840 in LJ unit. L=128, Initial density = 0.7, Initial Temperature = 0.9. 1,468,006 particles. Cutoff 3.0 sigma. Visualized with ParaView.

Result 2

ă€€Simulation Details: System Size = 320*320*320, Initial Density 0.7 (N = 22,937,600), Temperature 0.9, TimeStep 0.005, Cutoff 3.0: After 10000-steps thermalization, system is expanded by 1.025. Put snapshots every 1000 steps, total 619 frames (t = 3095 in LJ Units). Flat mpi, 1024 processes at ISSP System B (SGI Altix Ace 8400EX). Visualized with ParaView.

Result 3

The same conditions as Result 2. Put snapshots every 1000 steps, total 500 frames. Visualized with Pov-Ray.

Result 4

A movie created by my coworker Dr. Inaoka.

Not for research, but just for fun!

Self-avoiding random walk on 200 times 200 square lattice with fixing starting and ending point.

Self-avoiding random walk with an idea of fractal. The size of lattice starts from 4x4 and ends at 128x128.

Not for research, but just for fun!

Creating maze pattern.

Maze Pattern [gif, 1.2MB]MPEG format

- Small [mpg, 1.2MB]
- Large [mpg, 1.2MB] (With different initial condition)

Creating labyrinth, self-replicating, etc, depending on parameters.

Labyrinth [gif, 0.8MB] | Self-replicating [gif, 2.5MB] |

Visualized by ParaView.

Parameters: L=80*80; F=0.04; k = 0.06075; Du = 0.05; Dv = 0.1; dt = 0.2; 1frame/100steps; 200 frames; Total 20000 steps.

Not for research, but just for fun!

A structural analysis of a Chinese character "hito" which means "people". Dynamic and static loads are visualized for a system containing particles and springs under uniform gravity. Red parts denote high internal stresses. 772 particles, 2788 bonds (springs), 20000 steps, and 200 frames.

Source codes are available at GitHub.

A structural analysis of a Japanese character (Hiragana) "wo". Dynamic and static loads are visualized for a system containing particles and springs under uniform gravity. Red parts denote high internal stresses. 11971 particles, 46275 bonds (springs), 20000 steps, and 200 frames.

Back to the Top page