The Robustness Assessment for Event Drıven Molecular Dynamıcs By Calculatıng Speed of Sound

Abstract

Event-driven molecular dynamics (EDMD) is a special application of Molecular Dynamics (MD) derived from kinetic theory of gases. While classical solution of Navier-Stokes equations fails at high Knudsen (Kn) number flows, EDMD is valid on entire regime. Interaction potentials are considered discrete and exist only at the moment of impact. Hence, molecule trajectories are linear. Unlike the classical MD, this helps to simulate bigger systems. Molecular interactions, interaction times and partners can be predicted deterministically. Diatomic and polyatomic molecules are handled by an implemented energy relaxation scheme. Calculation of possible event times and determination of the earliest are the most timeconsuming steps of the simulation. In order to improve computational speed, a cell partitioning methodology and a priority queue structure are implemented in this study. The effect of the implementations on the performance is investigated and optimum simulation parameters are determined; when using PQ, number of cells must be greater than the number of molecules for maximum computational speed. Robustness assessments for the implementations are conducted with a realworld problem. Extreme density difference in confined geometries has vast usage in engineering and is also a good example of stress test because of its complex nature. This paper addresses the calculation of sound speed in a shock tube filled with a diatomic gas by using EDMD simulations. The robustness is validated since the results agrees perfectly with the theoretical values.

Olay güdümlü moleküler dinamik (OGMD) gazların kinetik teorisinden türetilen moleküler dinamiğin (MD) özel bir uygulamasıdır. Navier-Stokes denklemlerinin klasik çözümü yüksek Knudsen (Kn) akışlarında başarısız olurken, OGMD tüm akış rejimi için geçerlidir. Etkileşim potansiyellerinin süreksiz ve sadece temas anında mevcut olduğu kabul edilir. Bu sayede molekül yörüngeleri doğrusaldır. Bu klasik MD’nin aksine büyük sistemlerin simülasyonunu mümkün kılar. Molekül etkileşimleri, etkileşim zamanları ve çiftleri deterministik olarak öngörülebilirdir. Diatomik ve poliatomik moleküller uyarlanan bir enerji gevşeme düzeni ile modellenir. Muhtemel olay zamanlarının hesabı ve en erken olanın seçimi simülasyonun en çok zaman alan adımlarıdır. Hesaplama hızını geliştirmek amacıyla bu çalışmada bir hücre bölümleme metodolojisi ve öncelik kuyruğu yapısı simülasyona uyarlanmıştır. Uyarlamaların performans üzerine etkileri incelenmiş ve optimum süreç parametreleri belirlenmiştir; öncelik kuyruğunun kullanılması durumunda maksimum hesaplama performansı için hücre sayısının molekül sayısından fazla olması gerekmektedir. Uyarlamaların sağlamlık değerlendirmesi gerçek bir problem ile yapılmıştır. Kapalı geometrilerdeki aşırı yoğunluk farkı mühendislikte geniş bir kullanım alanı bulmakla birlikte karmaşık doğasından dolayı stres testi için uygun bir örnektir. Bu çalışma OGMD simülasyonları kullanarak diatomik bir gaz ile doldurulmuş şok tüpünde ses hızının hesaplanmasını ele almaktadır. Sonuçların teorik değerlerle mükemmel olarak uyum göstermiş olması sağlamlığı doğrulamaktadır.