Interpebble contact radius in a comet nucleus

In recent years, the gravitational collapse of pebble clumps in the early Solar System has been regarded as a plausible scenario for the origin of comets. In this context, ``pebbles'' represent mm- to cm-sized dust aggregates composed of (sub)micron-sized dust particles, and the structure of km-sized comets is thought to be an agglomerate of pebbles. The contact radius for pebble-pebble contacts was modelled in an earlier study; however, the pressure dependence of the interpebble contact radius was not considered. Here, we revisit the interpebble contact radius in a comet nucleus. We calculated the interpebble contact radius based on JKR contact theory, and we took into consideration the effect of lithostatic pressure. We found that the interpebble contact radius varies with depth from the surface, and the earlier model underestimated it by one order of magnitude at the centre of the comet nucleus.Comment: 9 pages, 6 figures. Accepted for publication in MNRA

An efficient multi-dimensional implementation of VSIAM3 and its applications to free surface flows

We propose an efficient multidimensional implementation of VSIAM3 (volume/surface integrated average-based multi-moment method). Although VSIAM3 is a highly capable fluid solver based on a multi-moment concept and has been used for a wide variety of fluid problems, VSIAM3 could not simulate some simple benchmark problems well (for instance, lid-driven cavity flows) due to relatively high numerical viscosity. In this paper, we resolve the issue by using the efficient multidimensional approach. The proposed VSIAM3 is shown to capture lid-driven cavity flows of the Reynolds number up to Re = 7500 with a Cartesian grid of 128 × 128, which was not capable for the original VSIAM3. We also tested the proposed framework in free surface flow problems (droplet collision and separation of We = 40 and droplet splashing on a superhydrophobic substrate). The numerical results by the proposed VSIAM3 showed reasonable agreements with these experiments. The proposed VSIAM3 could capture droplet collision and separation of We = 40 with a low numerical resolution (8 meshes for the initial diameter of droplets). We also simulated free surface flows including particles toward non-Newtonian flow applications. These numerical results have showed that the proposed VSIAM3 can robustly simulate interactions among air, particles (solid), and liquid

MovieMaker: A Parallel Movie-Making Software for Large Scale Simulations

We have developed a parallel rendering software for scientific visualization of large-scale, three-dimensional, time development simulations. The goal of this software, MovieMaker, is to generate a movie, or a series of visualization images from totally one TB-scale data within one night (or less than 12 hours). The isocontouring, volume rendering, and streamlines are implemented. MovieMaker is a parallel program for the shared memory architecture with dynamic load balancing and overlapped disk I/O.Comment: 3pages, 5figures, submitted to J. Plasma Physcs (special issue for 19th ICNSP

Performance of large scaled tsunami run-up analysis using explicit ISPH method

The tsunami run-up simulation by the particle method at city level needs to huge number of particle at least 1 billion particles. The conventional particle simulation method is not easy to solve these huge problem even on the premise of using supercomputer. Then, a new particle method ’fully explicit Incompressible SPH’ is developed that takes into consideration both calculation efficiency and accuracy. Finally, we demonstrate the future plan how to use our simulation resultes for a practical ’Soft’ disaster mitigation method through the evacuation education with the Virtual Reality(VR) system

Interparticle normal force in highly porous granular matter during compression

We perform a numerical simulation of compression of a highly porous dust aggregate of monodisperse spheres. We find that the average interparticle normal force within the aggregate is inversely proportional to both the filling factor and the average coordination number, and we also derive this relation theoretically. Our findings would be applicable for granular matter of arbitrary structures, as long as the constituent particles are monodisperse spheres.Comment: 9 pages, 7 figures. Accepted for publication in PR

Control of Fault Weakening on the Structural Styles of Underthrusting-Dominated Non-Cohesive Accretionary Wedges

Underthrusting is a typical process at compressive margins responsible for nappe stacking and sediment subduction. In nature, underthrusting is often associated with weak basal faults, although static mechanical analysis (critical taper theory) suggests that weak basal faults promote accretion while strong basal faults promote underthrusting. We perform mathematical analyses and numerical simulations to determine whether permanent fault weakening promotes or inhibits underthrusting. We investigate the control of permanent fault weakening on the dynamics of a strong-based ( (1 − λ ∗ b )μ b ≈ (1 − λ ∗ )μ ) non-cohesive wedge ( μ and μ b are internal and basal friction, respectively). We control the wedge material strength by a spatially constant fluid overpressure factor ( λ ∗ ), and fault strength by a plastic strain weakening factor ( χ ). First, we use the critical taper theory to determine a mechanical mode diagram that predicts structural styles. Then, we perform numerical simulations of accretionary wedge formation to establish their dynamical structural characteristics. We determine a continuum of structural styles between three end-members which correspond to the theoretical mechanical mode transitions. Style 1 is characterized by thin tectonic slices and little to no underthrusting. Style 2 shows thick slices, nappe stacking, and shallow gravity-driven tectonics. Style 3 displays the complete underthrusting of the incoming sediments, that are exhumed when they reach the backstop. We conclude that in the condition of an initially strong wedge base, permanent fault weakening promotes underthrusting. Thus, this contribution enlightens the control of the dynamic evolution of material properties on the formation of subduction channels, slope instabilities, and antiformal nappe stacks

Implicit solution of the material transport in Stokes flow simulation: Toward thermal convection simulation surrounded by free surface

We present implicit time integration schemes suitable for modeling free surface Stokes flow dynamics with marker in cell (MIC) based spatial discretization. Our target is for example thermal convection surrounded by deformable surface boundaries to simulate the long term planetary formation process. The numerical system becomes stiff when the dynamical balancing time scale for the increasing/decreasing load by surface deformation is very short compared with the time scale associated with thermal convection. Any explicit time integration scheme will require very small time steps; otherwise, serious numerical oscillation (spurious solutions) will occur. The implicit time integration scheme possesses a wider stability region than the explicit method; therefore, it is suitable for stiff problems. To investigate an efficient solution method for the stiff Stokes flow system, we apply first (backward Euler (BE)) and second order (trapezoidal method (TR) and trapezoidal rule—backward difference formula (TR-BDF2)) accurate implicit methods for the MIC solution scheme. The introduction of implicit time integration schemes results in nonlinear systems of equations. We utilize a Jacobian free Newton Krylov (JFNK) based Newton framework to solve the resulting nonlinear equations. In this work we also investigate two efficient implicit solution strategies to reduce the computational cost when solving stiff nonlinear systems. The two methods differ in how the advective term in the material transport evolution equation is treated. We refer to the method that employs Lagrangian update as “fully implicit” (Imp), whilst the method that employs Eulerian update is referred to as “semi-implicit” (SImp). Using a finite difference (FD) method, we have performed a series of numerical experiments which clarify the accuracy of solutions and trade-off between the computational cost associated with the nonlinear solver and time step size. In comparison with the general explicit Euler method, the second order accurate Imp methods reduce total computational cost successfully through the utilization of a large time step without sacrificing accuracy and stability. Moreover, the proposed SImp method is effective in reducing the computational cost associated with evaluating the nonlinear residual while obtaining a solution similar to the Imp method.ISSN:0010-4655ISSN:1879-294