shocktube/shockTube1DSod_AirWater.yaml

---
test:
  # a unique test name for this integration tests
  name: shockTube1DSod_AirWater
  # create a default assert that compares the log file
  assert: "inputs/shocktube/shockTube1DSod_AirWater.txt"

environment:
  title: SOD_120pts_AirWater
  tagDirectory: true
arguments: { }
timestepper:
  name: theMainTimeStepper
  arguments:
    ts_type: rk
    ts_max_time: 0.00024
    ts_dt: 1e-10
    ts_max_steps: 10
  domain: !ablate::domain::BoxMesh
    name: simpleBoxField
    faces: [ 120 ]
    lower: [ 0 ]
    upper: [ 1 ]
    boundary: [ NONE ]
    modifiers:
      - !ablate::domain::modifiers::DistributeWithGhostCells
      - !ablate::domain::modifiers::GhostBoundaryCells
    fields:
      - !ablate::finiteVolume::CompressibleFlowFields
        eos: !ablate::eos::PerfectGas
          parameters:
            gamma: 0
            Rgas: 0
      - name: densityvolumeFraction
        type: FVM
      - name: volumeFraction
        type: FVM
      - name: pressure
        location: AUX
        type: FVM
  initialization:
    - &eulerField
      fieldName: "euler" # air/water changliou2007
      field: >-
        x < 0.5 ? 11307.222318647702 : 967.7979066316877,
        x < 0.5 ? 2500000000.0 : 2414070815.450644,
        0
    - &rhoAlpha
      fieldName: densityvolumeFraction
      field: "x < 0.5 ? 11307.222318647702 : 0.0"
    - &alpha
      fieldName: volumeFraction
      field: " x < 0.5 ? 1.0 : 0.0 "
solver: !ablate::finiteVolume::FiniteVolumeSolver
  id: SOD Problem
  processes:
    - !ablate::finiteVolume::processes::TwoPhaseEulerAdvection
      eos: !ablate::eos::TwoPhase
        eos1: !ablate::eos::PerfectGas &eosAir
          parameters:
            gamma: 1.4
            Rgas: 287.0
        eos2: !ablate::eos::StiffenedGas &eosWater
          parameters:
            gamma: 1.932
            Cp: 8095.08
            p0: 1164500000.0
      fluxCalculatorGasGas: !ablate::finiteVolume::fluxCalculator::RiemannStiff
        eosL: *eosAir
        eosR: *eosAir
      fluxCalculatorGasLiquid: !ablate::finiteVolume::fluxCalculator::RiemannStiff
        eosL: *eosAir
        eosR: *eosWater
      fluxCalculatorLiquidGas: !ablate::finiteVolume::fluxCalculator::RiemannStiff
        eosL: *eosWater
        eosR: *eosAir
      fluxCalculatorLiquidLiquid: !ablate::finiteVolume::fluxCalculator::RiemannStiff
        eosL: *eosWater
        eosR: *eosWater

  boundaryConditions:
    - !ablate::finiteVolume::boundaryConditions::EssentialGhost
      boundaryName: "air walls euler"
      labelIds: [ 1, 2 ]
      boundaryValue: *eulerField
    - !ablate::finiteVolume::boundaryConditions::EssentialGhost
      boundaryName: "vf walls"
      labelIds: [ 1, 2 ]
      boundaryValue: *rhoAlpha
    - !ablate::finiteVolume::boundaryConditions::EssentialGhost
      boundaryName: "alpha walls"
      labelIds: [ 1, 2 ]
      boundaryValue: *alpha
  monitors:
    - !ablate::monitors::TimeStepMonitor
    - !ablate::monitors::CurveMonitor
      interval: 0