Algorithms
rep-no: UCRL-LR--121751
Viscosity
Matter
United States. Department of Energy.
Lawrence Livermore National Laboratory
Mathematical Models
Nonequilibrium flows with smooth particle applied mechanics
Other Information: DN: Thesis submitted by Oyeon Kum to University of California, Davis; TH: Thesis (Ph.D.)
Classical Mechanics
1995-07-01
Kum, O.
Instability
ark: ark:/67531/metadc623301
other: DE95017266
grantno: W-7405-ENG-48
99 Mathematics, Computers, Information Science, Management, Law, Miscellaneous
Turbulence
Shock Waves
Smooth particle methods are relatively new methods for simulating solid and fluid flows through they have a 20-year history of solving complex hydrodynamic problems in astrophysics, such as colliding planets and stars, for which correct answers are unknown. The results presented in this thesis evaluate the adaptability or fitness of the method for typical hydrocode production problems. For finite hydrodynamic systems, boundary conditions are important. A reflective boundary condition with image particles is a good way to prevent a density anomaly at the boundary and to keep the fluxes continuous there. Boundary values of temperature and velocity can be separately controlled. The gradient algorithm, based on differentiating the smooth particle expression for (u{rho}) and (T{rho}), does not show numerical instabilities for the stress tensor and heat flux vector quantities which require second derivatives in space when Fourier`s heat-flow law and Newton`s viscous force law are used. Smooth particle methods show an interesting parallel linking to them to molecular dynamics. For the inviscid Euler equation, with an isentropic ideal gas equation of state, the smooth particle algorithm generates trajectories isomorphic to those generated by molecular dynamics. The shear moduli were evaluated based on molecular dynamics calculations for the three weighting functions, B spline, Lucy, and Cusp functions. The accuracy and applicability of the methods were estimated by comparing a set of smooth particle Rayleigh-Benard problems, all in the laminar regime, to corresponding highly-accurate grid-based numerical solutions of continuum equations. Both transient and stationary smooth particle solutions reproduce the grid-based data with velocity errors on the order of 5%. The smooth particle method still provides robust solutions at high Rayleigh number where grid-based methods fails.
Hydrodynamics
osti: 105054
199 p.
66 Physics
Fluid Mechanics