# <Problem_Description>
#
# Test: mesh_deformation_file
#
# Runs a mesh deformation case where deformations are specified in an
# Exodus file with a coarse mesh, interpolated to the fluid tet mesh, then
# the Lighthill tensor divergence is output to yet another medium mesh.
#
# </Problem_Description>


BEGIN SIERRA FUEGO


  TITLE Flow in a channel

$=====================================================================
$ Define the fluid mesh. This is a tet mesh here
$=====================================================================

  BEGIN FINITE ELEMENT MODEL box
    Database Name = fluid_mesh.g
    Database Type = EXODUSII
    Decomposition Method = RCB
    BEGIN PARAMETERS FOR BLOCK block_1
      MATERIAL air
    END   PARAMETERS FOR BLOCK block_1
  END   FINITE ELEMENT MODEL box

$=====================================================================
$ Define what file the displacements are saved in. This needs to either
$ be a single file or be decomposed on the same number of processors
$ as the current run.
$=====================================================================

  BEGIN FINITE ELEMENT MODEL displacements_mesh
     Database Name = displacements.e
     Decomposition Method = RCB
     Database Type = EXODUSII
  END   FINITE ELEMENT MODEL displacements_mesh

$=====================================================================

  BEGIN PROPERTY SPECIFICATION FOR FUEGO MATERIAL air
    DENSITY   = 1.1
    VISCOSITY = 1.8e-5
  END   PROPERTY SPECIFICATION FOR FUEGO MATERIAL air

$=====================================================================
  BEGIN TPETRA EQUATION SOLVER  CONTINUITY
    BEGIN GMRES SOLVER
      BEGIN MUELU PRECONDITIONER
      END
      MAXIMUM ITERATIONS = 200
      RESTART ITERATIONS = 200
      RESIDUAL SCALING = R0
      CONVERGENCE TOLERANCE = 1.000000e-10
    END
  END TPETRA EQUATION SOLVER


  BEGIN TPETRA EQUATION SOLVER  SCALAR
    BEGIN GMRES SOLVER
      BEGIN SGS PRECONDITIONER
        NUMBER OF SWEEPS = 1
      END
      MAXIMUM ITERATIONS = 50
      RESTART ITERATIONS = 50
      RESIDUAL SCALING = R0
      CONVERGENCE TOLERANCE = 1.000000e-10
    END
  END TPETRA EQUATION SOLVER


$=====================================================================
 BEGIN FUEGO PROCEDURE fuego_procedure
$=====================================================================
   Begin Solution Control Description
     Use System Main

     Begin System Main
       Begin Transient Transient1
         Advance displ_io_region
         Transfer displ_to_fluid
         Advance Fuego_region
       End   Transient Transient1
     End   System Main

     Begin Parameters For Transient Transient1
       Start Time  = 0.0
       Initial Deltat = 0.001e-4
       Termination Time Is 2.5e-3

       Begin Parameters for Fuego Region Fuego_region
          TRANSIENT STEP TYPE IS automatic
          cfl limit = 0.5
          time step change factor = 1.25
       End   Parameters for Fuego Region Fuego_region
     End   Parameters For Transient Transient1
   End   Solution Control Description

   begin transfer displ_to_fluid
     interpolate volume nodes from displ_io_region to fuego_region
     send block block_1 to block_1
     send field D state none to mesh_displacement state none
   end transfer displ_to_fluid

$=====================================================================
    BEGIN FUEGO REGION fuego_region
$=====================================================================

      BEGIN SOLUTION OPTIONS
         ACTIVATE EQUATION Continuity
         ACTIVATE EQUATION X_Momentum
         ACTIVATE EQUATION Y_Momentum
         ACTIVATE EQUATION Z_Momentum

         ACTIVATE LIGHTHILL TENSOR POSTPROCESSING
$=====================================================================
$ Define what mesh the Lighthill tensor divergence will be saved on
$=====================================================================
         BEGIN Acoustic Transfer Output acoustic_output
           INPUT MESH NAME = acoustic_mesh.g
           OUTPUT MESH NAME = acoustic.e
           FORCE SEARCH IN MODEL COORDINATES
           SEND BLOCK block_1 to block_1
         END

         PROJECTION METHOD = fourth_order SMOOTHING WITH timestep SCALING

         MINIMUM NUMBER OF NONLINEAR ITERATIONS = 1
         MAXIMUM NUMBER OF NONLINEAR ITERATIONS = 3

         USE EQUATION SOLVER continuity FOR EQUATION Continuity
         USE EQUATION SOLVER scalar FOR EQUATION X_Momentum
         USE EQUATION SOLVER scalar FOR EQUATION Y_Momentum
         USE EQUATION SOLVER scalar FOR EQUATION Z_Momentum

         FIRST ORDER UPWIND FACTOR = 0.0
         HYBRID UPWIND FACTOR = 0.0
         UPWIND METHOD IS UPW

         BEGIN MESH MOTION SPECIFICATION flex_mesh
           MESH MOTION MODEL ON all_blocks = TRANSFERRED
         END MESH MOTION SPECIFICATION flex_mesh
      END SOLUTION OPTIONS

$=====================================================================

      USE FINITE ELEMENT MODEL box

$=====================================================================

      Begin Results Output Label output
        DATABASE Name = flow.e
        At Step 0, Increment = 1
        TITLE Flow in a channel
        NODAL Variables = pressure AS P
        NODAL Variables = x_velocity AS Ux
        NODAL Variables = y_velocity AS Uy
        NODAL Variables = z_velocity AS Uz
        NODAL Variables = mesh_displacement AS D
        NODAL Variables = mesh_velocity AS UM
        NODAL Variables = control_volume
        NODAL Variables = div_mesh_velocity
        NODAL Variables = div_lighthill_tensor AS divT
      End   Results Output Label output

$=====================================================================

      Begin Initial Condition Block blah
        volume is block_1

        pressure = 0.0
        x_velocity = 0.0
        y_velocity = 0.0
        z_velocity = 0.0
      End Initial Condition Block blah

$=====================================================================

      begin wall boundary condition on surface surface_1
        x_velocity = 0.0
        y_velocity = 0.0
        z_velocity = 0.0
        VELOCITY IS RELATIVE TO MESH    # makes this a moving no-slip wall
      end   wall boundary condition on surface surface_1

      begin open boundary condition on surface surface_2
        Momentum Entrainment = Computed
        pressure = 0.0
      end   open boundary condition on surface surface_2

      begin wall boundary condition on surface surface_3
        x_velocity = 0.0
        y_velocity = 0.0
        z_velocity = 0.0
        VELOCITY IS RELATIVE TO MESH    # makes this a moving no-slip wall
      end   wall boundary condition on surface surface_3

      begin wall boundary condition on surface surface_4
        x_velocity = 0.0
        y_velocity = 0.0
        z_velocity = 0.0
        VELOCITY IS RELATIVE TO MESH    # makes this a moving no-slip wall
      end   wall boundary condition on surface surface_4


    END   FUEGO REGION fuego_region

$=====================================================================
$ IO Region to read displacements from
$  This handles interpolation in time when the time steps between the
$  saved displacements and fluid time steps are different
$=====================================================================

    BEGIN INPUT_OUTPUT REGION displ_io_region
      USE FINITE ELEMENT MODEL displacements_mesh
    END   INPUT_OUTPUT REGION displ_io_region

$=====================================================================

  END   FUEGO PROCEDURE fuego_procedure

END   SIERRA FUEGO