SIMULATION METHOD, SIMULATION PROGRAM, AND SIMULATION DEVICE

摘要

A renormalization transformation process is performed for a granular system S which is a simulation target based on a renormalization factor α depending on the number of renormalizations. Position vectors and momentum vectors of grains of a renormalized granular system S′ are calculated by executing molecular dynamics calculation for the renormalized granular system S′. An interaction potential φ between grains of the granular system S is expressed as φ(r)=εf((r−r0)/σ), where ε represents an interaction coefficient having a dimension of energy, f represents a non-dimensional function, r0 and σ represent parameters characterizing grains, and r represents an inter-grain distance. When a dimensionality of a space of the granular system S is represented as d, by applying transformation laws expressed as N′=N/αd, m′=mαd, ε′=εαd, r0′=αr0, and σ′=ασ, the molecular dynamics calculation is executed based on an interaction potential of the renormalized granular system S′ expressed as φ′(r)=ε′f((r−r0′)/σ′).

主权项

1. A simulation method comprising:
a process of performing a renormalization transformation process with respect to a granular system S which is a simulation target formed of a plurality of grains based on a renormalization factor α depending on the number of renormalizations; and a process of calculating a position vector and a momentum vector of a grain of a renormalized granular system S′, by executing molecular dynamics calculation with respect to the renormalized granular system S′, wherein when an interaction potential φ between the grains of the granular system S is expressed as follows,φ(r)=ɛ·f(r-r0σ) where ε represents an interaction coefficient having a dimension of energy, f represents a non-dimensional function, r0 and σ represent parameters characterizing a grain, and r represents an inter-grain distance, and when a dimensionality of a space of the granular system S is represented as d, by applying transformation laws expressed as follows,N′=Nαdm′=m·αdɛ′=ɛ·αdr0′=α·r0σ′=α·σ the molecular dynamics calculation is executed based on an interaction potential of the renormalized granular system S′ expressed as follows:φ′(r)=ɛ′·f(r-r0′σ′)