add inviscid flow states, aero loads tests for flat plate

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_flow/Adjoint/.gitignore

@@ -0,0 +1 @@
+debug.getgrad

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_flow/Flow/fun3d.nml

@@ -0,0 +1,82 @@
+&project
+  project_rootname = 'miniMesh'
+/
+&raw_grid
+  grid_format = 'aflr3'
+  data_format = 'default'
+  swap_yz_axes = .false.
+/
+&governing_equations
+  eqn_type = 'compressible'
+  viscous_terms = 'inviscid'
+/
+&reference_physical_properties
+  mach_number     = 0.3
+  angle_of_attack = 10.0
+  reynolds_number = 6.91137e6
+  temperature     = 300.0
+  temperature_units = 'Kelvin'
+/
+&inviscid_flux_method
+  flux_construction       = 'roe'
+  flux_limiter            = 'none'
+  smooth_limiter_coeff    = 1.0
+  freeze_limiter_iteration = 100
+/
+&code_run_control
+  restart_write_freq = 500
+  restart_read       = 'off'
+  steps              =  500
+/
+&nonlinear_solver_parameters
+  time_accuracy = "steady"
+  time_step_nondim =   0.1
+  subiterations =          0
+  schedule_iteration(1:2) = 1,200
+  schedule_cfl(1:2) = 1.0, 100.0
+  schedule_cflturb(1:2) = 5.0,50.0
+  hanim = .false.
+/
+&hanim
+  cfl_increase_factor = 1.2
+  cfl_decrease_factor = 0.7
+/ 
+&elasticity_gmres
+  elasticity = 2
+  tol = 1.e-20
+  tol_abs = 1.e-30
+/
+&global
+  moving_grid = .true.
+  boundary_animation_freq = 1000
+  volume_animation_freq = 1000
+  recompute_turb_dist = .false.
+/
+&slice_data
+  nslices	= 1
+  slice_z(:)	= .true.
+  slice_location(:) = 2.5
+/
+&volume_output_variables
+  export_to = 'tecplot'
+  temperature = .true.
+  x = .true.
+  y = .true.
+  z = .true.
+/
+&boundary_output_variables
+  number_of_boundaries = 1
+  boundary_list = '2'
+  temperature = .true.
+  x = .true.
+  y = .true.
+  z = .true.
+/
+&massoud_output
+  funtofem_include_skin_friction = .true.
+  funtofem_internal_adjoint_test = .false.
+  funtofem_test_epsilon = 1.e-4
+  funtofem_flow_test = .true.
+  funtofem_restart_grid_forward = .false.
+  funtofem_restart_grid_adjoint = .true.
+/

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_flow/Flow/miniMesh.mapbc

@@ -0,0 +1,3 @@
+2             
+1              5000           Farfield
+2              4000           wall

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_flow/Flow/moving_body.input

@@ -0,0 +1,9 @@
+&body_definitions
+  n_moving_bodies = 1,
+  body_name(1) = 'plate',
+  parent_name(1) = '',         ! '' means motion relative to inertial ref frame
+  n_defining_bndry(1) =  1,    ! number of boundaries that define this body
+  defining_bndry(1,1) = 2,     ! index 1: boundary number index 2: body number
+  motion_driver(1) = 'funtofem', ! tells fun3d to use motion inputs from python
+  mesh_movement(1) = 'deform',  ! can use 'rigid', 'deform', 'rigid+deform' with funtofem interface
+/

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_flow/Flow/perturb.input

@@ -0,0 +1,20 @@
+   PERTURB   EPSILON GRIDPOINT
+         0    1e-30       666
+
+0 = No perturbation
+1 = Mach number
+2 = Alpha
+3 = Shape
+4 = x-rotation rate
+5 = y-rotation rate
+6 = z-rotation rate
+7 = Grid point x
+8 = Grid point y
+9 = Grid point z
+10 = Yaw
+11 = error transport (truncation error)
+12 = RCS jet plenum pressure, p0
+100+ = add an imaginary source term to equation
+         PERTURB-100 of node GRIDPOINT
+         (to verify the adjoint lambda value)
+

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_flow/clear.sh

@@ -0,0 +1,9 @@
+rm Flow/*.dat*
+rm Flow/*.flow
+rm Flow/*.forces*
+rm Flow/*.grid_info
+rm Adjoint/*getgrad
+rm Adjoint/*.tx*
+rm Adjoint/*.dat
+rm Adjoint/*.adjoint
+rm Adjoint/*.grid_info

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_quads/.gitignore

@@ -0,0 +1,3 @@
+# have to force changes, doesn't show in git status
+Flow/
+Adjoint/

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_quads/Adjoint/.gitignore

@@ -0,0 +1 @@
+*

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_quads/Flow/fun3d.nml

@@ -0,0 +1,73 @@
+&project
+  project_rootname = 'miniMesh'
+/
+&raw_grid
+  grid_format = 'aflr3'
+  data_format = 'default'
+  swap_yz_axes = .false.
+/
+&governing_equations
+  eqn_type = 'compressible'
+  viscous_terms = 'inviscid'
+/
+&reference_physical_properties
+  mach_number     = 0.3
+  angle_of_attack = 10.0
+  reynolds_number = 6.91137e6
+  temperature     = 300.0
+  temperature_units = 'Kelvin'
+/
+&inviscid_flux_method
+  flux_construction       = 'roe'
+  flux_limiter            = 'none'
+  smooth_limiter_coeff    = 1.0
+  freeze_limiter_iteration = 100
+/
+&code_run_control
+  restart_write_freq = 500
+  restart_read       = 'off'
+  steps              =  500
+/
+&nonlinear_solver_parameters
+  time_accuracy = "steady"
+  time_step_nondim =   0.1
+  subiterations =          0
+  schedule_iteration(1:2) = 1,200
+  schedule_cfl(1:2) = 1.0, 100.0
+  schedule_cflturb(1:2) = 5.0,50.0
+/
+&elasticity_gmres
+  elasticity = 2
+  tol = 1.e-20
+  tol_abs = 1.e-30
+/
+&global
+  moving_grid = .true.
+  boundary_animation_freq = 1000
+  volume_animation_freq = 1000
+  recompute_turb_dist = .false.
+/
+&slice_data
+  nslices	= 1
+  slice_z(:)	= .true.
+  slice_location(:) = 2.5
+/
+&volume_output_variables
+  export_to = 'tecplot'
+  temperature = .true.
+  x = .true.
+  y = .true.
+  z = .true.
+/
+&boundary_output_variables
+  number_of_boundaries = 1
+  boundary_list = '2'
+  temperature = .true.
+  x = .true.
+  y = .true.
+  z = .true.
+/
+&massoud_output
+  funtofem_include_skin_friction = .true.
+  funtofem_flow_test = .true.
+/

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_quads/Flow/miniMesh.mapbc

@@ -0,0 +1,3 @@
+2             
+1              5000           Farfield
+2              4000           wall

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_quads/Flow/moving_body.input

@@ -0,0 +1,9 @@
+&body_definitions
+  n_moving_bodies = 1,
+  body_name(1) = 'plate',
+  parent_name(1) = '',         ! '' means motion relative to inertial ref frame
+  n_defining_bndry(1) =  1,    ! number of boundaries that define this body
+  defining_bndry(1,1) = 2,     ! index 1: boundary number index 2: body number
+  motion_driver(1) = 'funtofem', ! tells fun3d to use motion inputs from python
+  mesh_movement(1) = 'deform',  ! can use 'rigid', 'deform', 'rigid+deform' with funtofem interface
+/

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_quads/Flow/perturb.input

@@ -0,0 +1,20 @@
+   PERTURB   EPSILON GRIDPOINT
+         0    1e-30       666
+
+0 = No perturbation
+1 = Mach number
+2 = Alpha
+3 = Shape
+4 = x-rotation rate
+5 = y-rotation rate
+6 = z-rotation rate
+7 = Grid point x
+8 = Grid point y
+9 = Grid point z
+10 = Yaw
+11 = error transport (truncation error)
+12 = RCS jet plenum pressure, p0
+100+ = add an imaginary source term to equation
+         PERTURB-100 of node GRIDPOINT
+         (to verify the adjoint lambda value)
+

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_quads/clear.sh

@@ -0,0 +1,9 @@
+rm Flow/*.dat*
+rm Flow/*.flow
+rm Flow/*.forces*
+rm Flow/*.grid_info
+rm Adjoint/*getgrad
+rm Adjoint/*.tx*
+rm Adjoint/*.dat
+rm Adjoint/*.adjoint
+rm Adjoint/*.grid_info

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_tris/.gitignore

@@ -0,0 +1,3 @@
+# have to force changes, doesn't show in git status
+Flow/
+Adjoint/

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_tris/Adjoint/.gitignore

@@ -0,0 +1 @@
+debug.getgrad

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_tris/Flow/fun3d.nml

@@ -0,0 +1,82 @@
+&project
+  project_rootname = 'miniMesh'
+/
+&raw_grid
+  grid_format = 'aflr3'
+  data_format = 'default'
+  swap_yz_axes = .false.
+/
+&governing_equations
+  eqn_type = 'compressible'
+  viscous_terms = 'inviscid'
+/
+&reference_physical_properties
+  mach_number     = 0.3
+  angle_of_attack = 10.0
+  reynolds_number = 6.91137e6
+  temperature     = 300.0
+  temperature_units = 'Kelvin'
+/
+&inviscid_flux_method
+  flux_construction       = 'roe'
+  flux_limiter            = 'none'
+  smooth_limiter_coeff    = 1.0
+  freeze_limiter_iteration = 100
+/
+&code_run_control
+  restart_write_freq = 500
+  restart_read       = 'off'
+  steps              =  500
+/
+&nonlinear_solver_parameters
+  time_accuracy = "steady"
+  time_step_nondim =   0.1
+  subiterations =          0
+  schedule_iteration(1:2) = 1,200
+  schedule_cfl(1:2) = 1.0, 100.0
+  schedule_cflturb(1:2) = 5.0,50.0
+  hanim = .false.
+/
+&hanim
+  cfl_increase_factor = 1.2
+  cfl_decrease_factor = 0.7
+/ 
+&elasticity_gmres
+  elasticity = 2
+  tol = 1.e-20
+  tol_abs = 1.e-30
+/
+&global
+  moving_grid = .true.
+  boundary_animation_freq = 1000
+  volume_animation_freq = 1000
+  recompute_turb_dist = .false.
+/
+&slice_data
+  nslices	= 1
+  slice_z(:)	= .true.
+  slice_location(:) = 2.5
+/
+&volume_output_variables
+  export_to = 'tecplot'
+  temperature = .true.
+  x = .true.
+  y = .true.
+  z = .true.
+/
+&boundary_output_variables
+  number_of_boundaries = 1
+  boundary_list = '2'
+  temperature = .true.
+  x = .true.
+  y = .true.
+  z = .true.
+/
+&massoud_output
+  funtofem_include_skin_friction = .true.
+  funtofem_internal_adjoint_test = .false.
+  funtofem_test_epsilon = 1.e-4
+  funtofem_grid_test = .false.
+  funtofem_restart_grid_forward = .false.
+  funtofem_restart_grid_adjoint = .false.
+/

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_tris/Flow/miniMesh.mapbc

@@ -0,0 +1,3 @@
+2             
+1              5000           Farfield
+2              4000           wall

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_tris/Flow/moving_body.input

@@ -0,0 +1,9 @@
+&body_definitions
+  n_moving_bodies = 1,
+  body_name(1) = 'plate',
+  parent_name(1) = '',         ! '' means motion relative to inertial ref frame
+  n_defining_bndry(1) =  1,    ! number of boundaries that define this body
+  defining_bndry(1,1) = 2,     ! index 1: boundary number index 2: body number
+  motion_driver(1) = 'funtofem', ! tells fun3d to use motion inputs from python
+  mesh_movement(1) = 'deform',  ! can use 'rigid', 'deform', 'rigid+deform' with funtofem interface
+/

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_tris/Flow/perturb.input

@@ -0,0 +1,20 @@
+   PERTURB   EPSILON GRIDPOINT
+         0    1e-30       666
+
+0 = No perturbation
+1 = Mach number
+2 = Alpha
+3 = Shape
+4 = x-rotation rate
+5 = y-rotation rate
+6 = z-rotation rate
+7 = Grid point x
+8 = Grid point y
+9 = Grid point z
+10 = Yaw
+11 = error transport (truncation error)
+12 = RCS jet plenum pressure, p0
+100+ = add an imaginary source term to equation
+         PERTURB-100 of node GRIDPOINT
+         (to verify the adjoint lambda value)
+

tests/fun3d_tests/fully_coupled_disc_deriv/inviscid/meshes/plate_tris/clear.sh

@@ -0,0 +1,9 @@
+rm Flow/*.dat*
+rm Flow/*.flow
+rm Flow/*.forces*
+rm Flow/*.grid_info
+rm Adjoint/*getgrad
+rm Adjoint/*.tx*
+rm Adjoint/*.dat
+rm Adjoint/*.adjoint
+rm Adjoint/*.grid_info