sandbox/lbattershill/README

    Introduction

    This folder is for examples in Basilisk compiled by Lily Battershill (Lily.Battershill@niwa.co.nz).

    The tsunami generation potential of pyroclastic density currents (PDCs) entering the sea is poorly understood, due to limited data and observations. Thus far, tsunami generation by PDCs has been modeled in a similar manner to tsunami generation associated with landslides or debris flows, using two-layer depth-averaged approaches. Using the adaptive partial differential equation solver Basilisk and benchmarking with published laboratory experiments, this work explores some of the important parameters not yet accounted for in numerical models of PDC-generated tsunamis. We use assumptions derived from experimental literature to approximate the granular, basal flow component of a PDC as a dense Newtonian fluid flowing down an inclined plane. This modeling provides insight into how the boundary condition of the slope and the viscosity of the dense granular-fluid influence the characteristics of the waves generated.

    The setup is based on the experimental geometry of Bougouin et al., 2020, where a laboratory fluidised granular flow runs down a ramp and into water. We assume the fluidized grains from the experiments of Bougouin et al., 2020 to behave as a continuum. This takes the form of a dense, viscous and incompressible Newtonian fluid. The dense fluid and the water are assumed to be miscible with one another, but immiscible with air: they are separated from the air by a sharp interface. Surface tension is assumed to have negligible effect on interaction dynamics and wave propagation, due to the contrast of scales. Please see fluidised_flow.c for details of implementation.

    We update conserving.h (see myconserving.h) to implement momentum-conserving VOF advection of the velocity components for the two-phase Navier–Stokes solver, with the addition of a density tracer on one side of the interface. This provides an alternative to a more classic “three phase” approach, where the phases are seperated by interfaces.

    Multi-Layer model: pyroclastic density current generated tsunami

    The setup is dimensionless and based on the experimental geometry of Bougouin et al., 2020, where a laboratory fluidised granular flow runs down a ramp and into water. We use assumptions from experimental literature to approximate the fluidized grains as a continuum, which takes the form of a dense, incompressible fluid. This assumption is confirmed to be adequate in the context of wave generation, by the numerical work of Battershill et al., 2021. This non-hydrostatic multi-layer numerical model assumes the dense fluid and the water to be miscible with one another. We consider only small variations in density between the granular-flow and the water, and we therefore justify the Boussinesq buoyancy approximation.

    • Layered density version We simulate a dense layer at base of flow, overlayed by a lighter layer. Note: the reference to density layering is different from the vertical disretisation of the multilayer scheme.
    • Single density version We simulate a dense flow (of a single density) propogating into water.

    References

    [battershill2021]

    Lily Battershill, Colin Whittaker, Emily Lane, Stephane Popinet, James White, William Power, and P Nomikou. Numerical simulations of a fluidized granular flow entry into water: insights into modeling tsunami generation by pyroclastic density currents. 2021.

    [bougouin2020]

    Alexis Bougouin, Raphael Paris, and Olivier Roche. Impact of fluidized granular flows into water: implications for tsunamis generated by pyroclastic flows. Journal of Geophysical Research: Solid Earth, 125(5):e2019JB018954, 2020. [ DOI ]