******************************************************************************** Calling TAU-Python constructors ******************************************************************************** ******************************************************************************** Executing TAU preprocessor ******************************************************************************** Implemented boundary treatments: -------------------------------- axisymmetry wedge plane symmetry plane axisymmetry axis farfield supersonic inflow supersonic outflow reservoir-pressure inflow exit-pressure outflow dirichlet euler wall sharp edge viscous wall engine exhaust engine inflow heat exchanger inflow heat exchanger outflow actuator exhaust actuator inflow chimera actuation mirror plane periodic plane Implemented grid metrics: Cell_Vertex Cell_Centered Implemented types of agglomeration: private MGridGen Directional_MGridGen Reading parameter file 'crm_test.para' Preprocessor parameters: Files/IO ------------------------------------------------------: - Boundary mapping filename: crm_test.para Primary grid filename: tau.grid Grid prefix: edgedata Old grid prefix: (none) : - Read partitioning filename: (none) Write graph filename: (none) Parameter -----------------------------------------------------: - Number of multigrid levels: 1 Output level: 5 Grid metric: Cell_Vertex Runtime optimisation ------------------------------------------: - Cache-coloring (0/max_faces in color): 20000 Bandwidth optimisation (0/1): 1 2D offset vector (0 / x=1,y=2,z=3): 2 Compute lusgs mapping (0/1): 1 Agglomeration -------------------------------------------------: - Type of agglomeration: private : - Agglomeration version: 1 Enhanced agglomeration (0/1): 1 Point fusing reward: 1.2 Structured grid coarsening: 0.9 Surface agglomeration ratio: 0 Wallnormal agglomeration ratio: 0 : - Type of coarsening for MGridGen: 4 Type of refinement for MGridGen: 6 Minsize for the coarse graph for MGridGen: 1 Maxsize for the coarse graph for MGridGen: 4 Extras --------------------------------------------------------: - Sharp edge angle (degrees): 0 Control volume edge weight (0/1): 1 Consider wall roughness: global Global wall roughness: 0 Roughness filename: (none) Compute exact whirlflux (0/1): 0 Axisymmetry (0/1): 0 Project wall distance: 0 Line search (0/1): 0 Write lines into a file (0/1): 0 Threshold parameter for aspect ratio: 4 Partitioning --------------------------------------------------: - Number of domains: 64 Type of partitioning (name): private Theta not for Tau ---------------------------------------------: - Preprocessing for incompressible solver (0/1): 0 Output shifted points grid (0/1): 0 Correct metric for boundary control volumes (0/1): 1 Project boundary control volumes coordinates (0/1): 1 Translation factor for shifted boundary points: 1 Compute parent faces (0/1): 0 Logfile control -----------------------------------------------: - Enable logfile output on all domains (0/1): 1 ----------------------------------------------------- Reading boundary mapping: file 'crm_test.para' Preprocessor boundary parameters ---------------------------------------------------------------: - Map: 0 Type: symmetry plane Subtype: Markers: 1 Partname: Name: SYMMETRY Multigrid priority: 0 block end ---------------------------------------------------------------: - Map: 1 Type: farfield Subtype: Markers: 2 Partname: Name: farfield Multigrid priority: 0 block end ---------------------------------------------------------------: - Map: 2 Type: viscous wall Subtype: turbulent Constant wall roughness: 0 Markers: 3 Partname: Name: FUSELAGE Multigrid priority: 1 block end ---------------------------------------------------------------: - Map: 3 Type: viscous wall Subtype: turbulent Constant wall roughness: 0 Markers: 4 Partname: Name: FAIRING Multigrid priority: 1 block end ---------------------------------------------------------------: - Map: 4 Type: viscous wall Subtype: turbulent Constant wall roughness: 0 Markers: 5 Partname: Name: WING Multigrid priority: 1 block end ----------------------------------------------------- Check primary grid input 33133 surface quads orientated marker info list contains 2 invalid entries 2 marker info entries removed ----------------------------------------------------- Boundary mapping for this grid: treatment-name: map triangles/quads symmetry plane: 0 No matching element! farfield: 1 29923 viscous wall: 2 1430 viscous wall: 3 No matching element! viscous wall: 4 1791 ----------------------------------------------------- Computing dual grid for the cell vertex grid metric No Tau-like 2D grid: odd number of grid points NOTE: No Tau-like 2D grid for offset direction 2 Only one dualgrid level requested. Skipping agglomeration step. Wall distances for grid level 1 Check dual grids ----------------------------------------------------- Multigrid level 1 points ................................: 879799 faces .................................: 2264269 min point neighbors ...................: 3 max point neighbors ...................: 7 volume ................................: 7.42623e+09 max. inner surface integral ...........: 6.03785e-12 at point / position ...................: 320243 (47.9513 29.3523 7.61243) max. boundary surface integral ........: 5.47056e-12 at point / position ...................: 320195 (47.9572 29.3638 7.61408) min, max surface ratio ................: 1 1 faces in farfield ..: 29571 faces in viscous wall ..: 1299 faces in viscous wall ..: 1692 ----------------------------------------------------- Domain decomposition ----------------------------------------------------- Coloring and sorting the inner faces free preprocessing data ******************************************************************************** Initializing TAU solver with mode ******************************************************************************** Signal handling Signal 15: Write data & exit after time step Signal 12: Write data after time step ----------------------------------------------------- Implemented solver boundary treatments: 0: axisymmetry wedge plane 1: symmetry plane 2: axisymmetry axis 3: farfield 4: supersonic outflow 5: supersonic inflow 6: reservoir-pressure inflow 7: exit-pressure outflow 8: dirichlet 9: actuation 10: engine exhaust 11: engine inflow 12: euler wall 13: sharp edge 14: viscous wall 15: chimera 16: actuator inflow 17: actuator exhaust 18: heat exchanger inflow 19: heat exchanger outflow 20: mirror plane 21: periodic plane ----------------------------------------------------- Implemented 7-Equation Turbulence Model(s): RSM: Reynolds stress model Implemented initialization types: walldist direct Implemented re-distribution models: Wilcox_RSM Wilcox_RSM.2010 SSG/LRR-w SSG/LRR-w.2010 LowRe_JH_linear LowRe_JH_nonlinear SSG/LRR-g.2010 Implemented dissipation models: isotropic anisotropic_JHh Implemented diffusion models: SGDH GGDH GGDH.2010 HL Implemented length scale equations: Wilcox_omega Kok_TNT_omega Menter_BSL_omega Hellsten_omega JHh_epsilon_hom JHhV1_epsilon_hom JHhV2_epsilon_hom Menter_BSL_g Implemented transition treatments: Point_transition Trip_isotropic Implemented vortical correction functions: none ----------------------------------------------------- Implemented relaxation solvers: Runge_Kutta Backward_Euler Implemented linear relaxation solvers: Runge_Kutta Runge_Kutta_adjoint Backward_Euler Implemented inviscid flux discretization types: Upwind Central Implemented upwind fluxes: Van_Leer AUSMDV AUSMP Roe AUSM_Van_Leer EFM MAPS+ Implemented central dissipation fluxes: Scalar_dissipation Matrix_dissipation Implemented central convective meanflow fluxes: Average_of_flux Flux_of_average Skew_symmetric_kok Implemented central convective turbulence fluxes: AUSMDV Roe Roe2nd Implemented versions of cell stretching coefficient: TAU HyperFlex HyperFlex_II Implemented smoothers: Point_explicit Upwind_implicit Implemented grid movement types: static: Grid does not move rigid: Rigid body motion deformed: Deforming grid deformed+rigid: Deforming grid combined with rigid body motion Implemented unsteady time stepping schemes: (none): Unsteady time stepping is off dual: Dual time stepping scheme global: Global time stepping scheme Implemented functions for computation of flow statistics: mean: compute mean values of rho, u, v, w, p, mue_turb, mue_lam variance: compute , , , , , , meanvelgrad: compute mean velocity gradients and mean pressure gradient meanturb: compute mean values of turbulence variables (depends on turbulence model) des: compute mean des state, i.e. fd, hybrid scale, l_rans/l_hybrid ----------------------------------------------------- Implemented "Field output values": gid for output of the global point id cp for output of the pressure coefficient Ptot for output of the total pressure apsimP for output of the apsim pressure temp for output of the temperature mach for output of the Mach number mean for output of the mean flow variables meanvelgrad for output of the mean velocity gradient meanturb for output of the mean turbulence values mean-des for output of the mean des values mean-cp for output of the mean and rms pressure coefficient residual-stress for output of the residual turbulent stresses variance for output of the variances of flow variables imach for output of the isentropic Mach number fluxblend for output of the flux para blending factor muetmue for output of the ratio muet mue rs-boussinesq for output of the specific Reynolds stresses (Boussinesq) des for output of the DES mode, length and length ratio des-func for output of the various DES functions (see user guide) alg-sens for output of the algebraic DES bl-sensors turbfilter for output of the LES/DES filter widths rotcorr for output of the rotational correction indicators sas for output of the Von Karman Length divergence for output of the divergence of vel testgrad for output of the test gradients dscale for output of the dissipation scaling tu for output of the turbulence intensity macht for output of the turbulent Mach number Rrho for output of the density residual Rrhov for output of the momentum residual RrhoE for output of the residual of rho times energy Rrhor11 for output of the residual of rho r11 (RSM) Rrhor22 for output of the residual of rho r22 (RSM) Rrhor33 for output of the residual of rho r33 (RSM) Rrhor12 for output of the residual of rho r12 (RSM) Rrhor13 for output of the residual of rho r13 (RSM) Rrhor23 for output of the residual of rho r23 (RSM) Rrhowrsm for output of the residual of rho omega (RSM) Rrhoepsrsm for output of the residual of rho epsilon (RSM) Rrhogrsm for output of the residual of rho g (RSM) krsm for output of the spec. kin. turb. energy for RSM grsm-limiter for output of the grsm limiter strain for output of the strainrate vort for output of the vorticity Nk for output of the kinematic vorticity heli for output of the helicity Hn for output of the normalized helicity l2 for output of the lambda2-criterion Q for output of the second invariant Q rotxyz for output of the rotation components wdist for output of the wall-distance wrough for output of the wall-roughness wnormal for output of the wall-normal domain for output of the point domain no xyzgeod for output of the coordinates in geodesic system vxyzgeod for output of the velocity in geodesic system xyzbody for output of the coordinates in body-fixed system vxyzbody for output of the velocity in body-fixed system xyzaero for output of the coordinates in aerodynamic system vxyzaero for output of the velocity in aerodynamic system vxyzmoving for output of the velocity in moving system vxyzmvgeod for output of the velocity in moving geod. system xyz for output of the grid coordinates xyzold for output of the old grid coordinates xyzoold for output of the oold grid coordinates xyzooold for output of the ooold grid coordinates volume for output of the cv volume oldvolume for output of the old cv volume ooldvolume for output of the oold cv volume oooldvolume for output of the ooold cv volume wxyz for output of the grid velocity wxyzgeod for output of the grid velocity in geod. system entropy for output of the normalized entropy difference gradrho for output of the density gradient gradu for output of the x_velocity gradient gradv for output of the y_velocity gradient gradw for output of the z_velocity gradient gradp for output of the pressure gradient visc for output of the laminar viscosity kappa for output of the thermal conductivity momentum for output of the momentum H for output of the total enthalpy rhoE for output of the density times total energy Kn-rho for output of the Knudsen number (gradrho) Kn-tau-q for output of the Knudsen number (stress, heat flux) mfp for output of the mean free path carbuncle for output of the carbuncle switch indicator blank for output of the chimera blanked points ipzone for output of the chimera interpolated points notfound for output of the chimera points not found hole for output of the chimera blanking state holeold for output of the old chimera blanking state holeoold for output of the oold chimera blanking state holeooold for output of the ooold chimera blanking state ----------------------------------------------------- Implemented "Surface output values": xyz for output of the grid coordinates restart for output of the variables used to restart rho for output of the density v for output of the velocity p for output of the pressure temp for output of the temperature mean for output of the mean flow variables mean-cp for output of the mean and rms pressure coefficient mean-cf for output of the mean skin friction mean-delta-plus for output of the mean max. normalized grid spacing variance for output of the variances of flow variables mach for output of the Mach number cp for output of the pressure coeff Ptot for output of the total pressure clydl for output of the local force coeff. vector cmxyzl for output of the local moment coeff. vector entropy for output of the normalized entropy difference apsimP for output of the apsim pressure acd-raysc for output of the acd ray scaling factor acd-chlen for output of the acd chord length acd-twist for output of the acd blade twist acd-vrot for output of the acd rotational velocity acd-alphai for output of the acd ind angle of attack acd-alpha for output of the acd eff angle of attack acd-cl for output of the acd coefficient of lift acd-cd for output of the acd coefficient of drag acd-secload for output of the acd tau sectional load acd-dsk-xyz for output of the acd disk point coords acd-dsk-secload for output of the acd disk sectional load acd-dsk-srfload for output of the acd disk surface load acd-dsk-force for output of the acd disk discrete force acd-bld-secload for output of the acd blade sectional load acd-bld-vel for output of the acd blade velocity htc for output of the heat transfer coefficient imach for output of the isentropic Mach number ch for output of the Stanton number heatfl for output of the wall heat flux visc for output of the laminar viscosity fxyz for output of the force vector fxyzgeod for output of the force vector in geodesic system cf for output of the skin friction cfxyz for output of the skin friction vector cfxyzgeod for output of the skin friction vector in geod. system delta-plus for output of the max. normalized grid spacing yplus for output of the y+ values krplus for output of the k_r+ values wrough for output of the wall-roughness cpharm for output of the harmonics of pressure coeff. cfharm for output of the harmonics of skin friction eddy for output of the eddy viscosity r11 for output of the Reynolds stress component R11 r22 for output of the Reynolds stress component R22 r33 for output of the Reynolds stress component R33 r12 for output of the Reynolds stress component R12 r13 for output of the Reynolds stress component R13 r23 for output of the Reynolds stress component R23 krsm for output of the Spec. kin. turb. energy (RSM) wrsm for output of the turb. length scale omega (RSM) epsrsm for output of the turb. length scale epsilon (RSM) grsm for output of the turb. length scale g (RSM) xyzgeod for output of the coordinates in geodesic system vxyzgeod for output of the velocity in geodesic system xyzbody for output of the coordinates in body-fixed system vxyzbody for output of the velocity in body-fixed system xyzaero for output of the coordinates in aerodynamic system vxyzaero for output of the velocity in aerodynamic system vxyzmoving for output of the velocity in moving system vxyzmvgeod for output of the velocity in moving geod. system volume for output of the cv volume wxyz for output of the grid velocity wxyzgeod for output of the grid velocity in geod. system momentum for output of the momentum H for output of the total enthalpy rhoE for output of the density times total energy massflux for output of the mass flux blank for output of the chimera blanked points ipzone for output of the chimera interpolated points vxyzedge for output of the boundary layer edge velocities bldatasw for output of the boundary layer data streamwise bldatacf for output of the boundary layer data cross flow ltflag for output of the laminar turbulent flag blinfo for output of the boundary layer information bldelta for output of the boundary layer thicknesses fnormal for output of the boundary normal vector ----------------------------------------------------- Implemented "Plane output values": xyz for output of the grid coordinates temp for output of the temperature v for output of the velocity components p for output of the pressure rho for output of the density cp for output of the pressure coeff Ptot for output of the total pressure pl for output of the pressure loss entropy for output of the normalized entropy difference cf for output of the skin friction cfxyz for output of the skin friction vector cfabs for output of the abs. skin friction yplus for output of the y+ values psi1 for output of the adjoint variable 1 (density) psi2 for output of the adjoint variable 2 (u) psi3 for output of the adjoint variable 3 (v) psi4 for output of the adjoint variable 4 (w) psi5 for output of the adjoint variable 5 (pressure) ----------------------------------------------------- Implemented monitoring values (to be separated by '_'). Monitoring values that can be used for Cauchy convergence control (CC) are indicated by an "X" in the second column: KEY CC COMMENT Rrho root mean square of density residual Rrhovx root mean square of x-component of momentum residual Rrhovy root mean square of y-component of momentum residual Rrhovz root mean square of z-component of momentum residual RrhoE root mean square of density times energy residual Rrhor11 root mean square of rho R11 residual Rrhor22 root mean square of rho R22 residual Rrhor33 root mean square of rho R33 residual Rrhor12 root mean square of rho R12 residual Rrhor13 root mean square of rho R13 residual Rrhor23 root mean square of rho R23 residual Rrhowrsm root mean square of rho omega residual Rrhoepsrsm root mean square of rho epsilon residual Rrhogrsm root mean square of rho g residual Max-Rrho maximum norm of rho residual X-max-Rrho grid x-coordinate of Max-Rrho Y-max-Rrho grid y-coordinate of Max-Rrho Z-max-Rrho grid z-coordinate of Max-Rrho Max-Rrho-global global maximum norm of density residual over all domains Max-Rrhor11 maximum norm of rho R11 residual X-max-Rrhor11 grid x-coordinate of Max-Rrhor11 Y-max-Rrhor11 grid y-coordinate of Max-Rrhor11 Z-max-Rrhor11 grid z-coordinate of Max-Rrhor11 Max-Rrhor22 maximum norm of rho R22 residual X-max-Rrhor22 grid x-coordinate of Max-Rrhor22 Y-max-Rrhor22 grid y-coordinate of Max-Rrhor22 Z-max-Rrhor22 grid z-coordinate of Max-Rrhor22 Max-Rrhor33 maximum norm of rho R33 residual X-max-Rrhor33 grid x-coordinate of Max-Rrhor33 Y-max-Rrhor33 grid y-coordinate of Max-Rrhor33 Z-max-Rrhor33 grid z-coordinate of Max-Rrhor33 Max-Rrhor12 maximum norm of rho R12 residual X-max-Rrhor12 grid x-coordinate of Max-Rrhor12 Y-max-Rrhor12 grid y-coordinate of Max-Rrhor12 Z-max-Rrhor12 grid z-coordinate of Max-Rrhor12 Max-Rrhor13 maximum norm of rho R13 residual X-max-Rrhor13 grid x-coordinate of Max-Rrhor13 Y-max-Rrhor13 grid y-coordinate of Max-Rrhor13 Z-max-Rrhor13 grid z-coordinate of Max-Rrhor13 Max-Rrhor23 maximum norm of rho R23 residual X-max-Rrhor23 grid x-coordinate of Max-Rrhor23 Y-max-Rrhor23 grid y-coordinate of Max-Rrhor23 Z-max-Rrhor23 grid z-coordinate of Max-Rrhor23 Max-Rrhowrsm maximum norm of rho omega residual X-max-Rrhowrsm grid x-coordinate of Max-Rrhowrsm Y-max-Rrhowrsm grid y-coordinate of Max-Rrhowrsm Z-max-Rrhowrsm grid z-coordinate of Max-Rrhowrsm Max-Rrhoepsrsm maximum norm of rho epsilon residual X-max-Rrhoepsrsm grid x-coordinate of Max-Rrhoepsrsm Y-max-Rrhoepsrsm grid y-coordinate of Max-Rrhoepsrsm Z-max-Rrhoepsrsm grid z-coordinate of Max-Rrhoepsrsm Max-Rrhogrsm maximum norm of rho g residual X-max-Rrhogrsm grid x-coordinate of Max-Rrhogrsm Y-max-Rrhogrsm grid y-coordinate of Max-Rrhogrsm Z-max-Rrhogrsm grid z-coordinate of Max-Rrhogrsm Res-GMRES total linear residual from GMRES solver C-lift X lift-coefficient C-lift-p X pressure part of C-lift C-lift-v X viscous part of C-lift C-drag X drag-coefficient C-drag-p X pressure part of C-drag C-drag-v X viscous part of C-drag C-sidef X coeff of sideforce C-sidef-p X pressure part of C-sidef C-sidef-v X viscous part of C-sidef C-mx X rolling moment coefficient C-mx-p X pressure part of rolling moment coefficient C-mx-v X viscous part of rolling moment coefficient C-my X pitching moment coefficient C-my-p X pressure part of pitching moment coefficient C-my-v X viscous part of pitching moment coefficient C-mz X yawing moment coefficient C-mz-p X pressure part of yawing moment coefficient C-mz-v X viscous part of yawing moment coefficient C-fx X force-coefficient (x-component) C-fx-p X pressure part of force-coefficient (x-component) C-fx-v X viscous part of force-coefficient (x-component) C-fy X force-coefficient (y-component) C-fy-p X pressure part of force-coefficient (y-component) C-fy-v X viscous part of force-coefficient (y-component) C-fz X force-coefficient (z-component) C-fz-p X pressure part of force-coefficient (z-component) C-fz-v X viscous part of force-coefficient (z-component) Fx X x-component of force [N] Fx-p X pressure part of x-component of force [N] Fx-v X viscous part of x-component of force [N] Fy X y-component of force [N] Fy-p X pressure part of y-component of force [N] Fy-v X viscous part of y-component of force [N] Fz X z-component of force [N] Fz-p X pressure part of z-component of force [N] Fz-v X viscous part of z-component of force [N] Mx X x-component of moment [Nm] Mx-p X pressure part of x-component of moment [Nm] Mx-v X viscous part of x-component of moment [Nm] My X y-component of moment [Nm] My-p X pressure part of y-component of moment [Nm] My-v X viscous part of y-component of moment [Nm] Mz X z-component of moment [Nm] Mz-p X pressure part of z-component of moment [Nm] Mz-v X viscous part of z-component of moment [Nm] Min-mach X minimum Machnumber in field X-min-mach grid x-coordinate of Min-mach Y-min-mach grid y-coordinate of Min-mach Z-min-mach grid z-coordinate of Min-mach Max-mach X maximum Machnumber in field X-max-mach grid x-coordinate of Max-mach Y-max-mach grid y-coordinate of Max-mach Z-max-mach grid z-coordinate of Max-mach Heatflow X total heatflow [W] Res-heat heatflux-increment Min-y+ X minimum value of y+ X-min-y+ grid x-coordinate of Min-y+ Y-min-y+ grid y-coordinate of Min-y+ Z-min-y+ grid z-coordinate of Min-y+ Max-y+ X maximum value of y+ X-max-y+ grid x-coordinate of Max-y+ Y-max-y+ grid y-coordinate of Max-y+ Z-max-y+ grid z-coordinate of Max-y+ Min-kr minimum value of kr X-min-kr grid x-coordinate of Min-kr Y-min-kr grid y-coordinate of Min-kr Z-min-kr grid z-coordinate of Min-kr Max-kr maximum value of kr X-max-kr grid x-coordinate of Max-kr Y-max-kr grid y-coordinate of Max-kr Z-max-kr grid z-coordinate of Max-kr Min-kr+ minimum value of kr+ X-min-kr+ grid x-coordinate of Min-kr+ Y-min-kr+ grid y-coordinate of Min-kr+ Z-min-kr+ grid z-coordinate of Min-kr+ Max-kr+ maximum value of kr+ X-max-kr+ grid x-coordinate of Max-kr+ Y-max-kr+ grid y-coordinate of Max-kr+ Z-max-kr+ grid z-coordinate of Max-kr+ act-area Area actuation bdry act-targ-mf Actuation target massflow act-pres Actuation bdry pressure Min-eddyv X min dim-less eddy viscosity X-min-eddyv grid x-coordinate of Min-eddyv Y-min-eddyv grid y-coordinate of Min-eddyv Z-min-eddyv grid z-coordinate of Min-eddyv Max-eddyv X max dim-less eddy viscosity X-max-eddyv grid x-coordinate of Max-eddyv Y-max-eddyv grid y-coordinate of Max-eddyv Z-max-eddyv grid z-coordinate of Max-eddyv Min-k X min kinetic energy X-min-k grid x-coordinate of Min-k Y-min-k grid y-coordinate of Min-k Z-min-k grid z-coordinate of Min-k Max-k X max kinetic energy X-max-k grid x-coordinate of Max-k Y-max-k grid y-coordinate of Max-k Z-max-k grid z-coordinate of Max-k Min-k2 X min 2nd turb. quantity X-min-k2 grid x-coordinate of Min-k2 Y-min-k2 grid y-coordinate of Min-k2 Z-min-k2 grid z-coordinate of Min-k2 Max-k2 X max 2nd turb. quantity X-max-k2 grid x-coordinate of Max-k2 Y-max-k2 grid y-coordinate of Max-k2 Z-max-k2 grid z-coordinate of Max-k2 Max-mtm X max eddy-viscosity/laminar-viscosity X-max-mtm grid x-coordinate of Max-mtm Y-max-mtm grid y-coordinate of Max-mtm Z-max-mtm grid z-coordinate of Max-mtm Angle-a angle of attack alpha Angle-b yaw angle beta Sideslip sideslip angle beta Farf-vx farfield vx Farf-vz farfield vz Mtm-2t number of pseudo time steps (unsteady) Alpha(t) pitch amplitude H(t) heave amplitude [grid unit] dalpha-dt pitch angular velocity dh-dt heave velocity [grid unit/s] S-k-e structural kinetic energy T/tperiod time normalized with period length Phi rotation angle phi Psi rotation angle psi Xi rotation angle xi Trans-x translation in x direction [grid unit] Trans-y translation in y direction [grid unit] Trans-z translation in z direction [grid unit] Farf-mach effective farfield Mach number at origin Farf-alpha effective farfield angle of attack alpha at origin [deg] Farf-beta effective farfield angle of attack beta at origin [deg] les-dt X timescale resolution of des les-percent X % of des points hybcentpercent X % of fluxes computed with pure central vort-center-x center of vorticity (x) vort-center-y center of vorticity (y) vort-center-z center of vorticity (z) vort-total X total vorticity dudx X dudx dvdy X dvdy dwdz X dwdz total-ke X total ke mean-dissip1 X molecular dissipation rmsu X rms u rmsv X rms v rmsw X rms w rmsp X rms pressure skew X velocity derivative skewness kurt X normalized velocity derivative kurtosis Taylor-f X Longitudinal Taylor length scale Taylor-g X Transverse Taylor length scale Taylor-Re X Re based on Taylor scale L2n-vc X L2 norm of vortical correction Max-vc X Max value of vortical correction Min-vc X Min value of vortical correction Res-lift cl-increment Res-drag cd-increment R-time Total CPU time (real) U-time Total CPU time (user) Rhs-time CPU time for one compute_residual() call. Norm-time Normalized CPU time (real) dJ/da-1 Adjoint & Primal: Partial dJ/da dJ/da-2 Adjoint: psi * partial dR/da, Primal: Partial dJ/dW * dW/da dJ/da Adjoint & Primal: derivative of cost fn wrt design var. Sens X Functional sensitivity Freqdom-Cl-Amp Freqdom: Amplitude of CL Freqdom-Cl-Phase Freqdom: Phase shift of CL Freqdom-Cd-Amp Freqdom: Amplitude of CD Freqdom-Cd-Phase Freqdom: Phase shift of CD Freqdom-Cmx-Amp Freqdom: Amplitude shift of CM_X Freqdom-Cmx-Phase Freqdom: Phase shift of CM_X Freqdom-Cmy-Amp Freqdom: Amplitude shift of CM_Y Freqdom-Cmy-Phase Freqdom: Phase shift of CM_Y Freqdom-Cmz-Amp Freqdom: Amplitude shift of CM_Z Freqdom-Cmz-Phase Freqdom: Phase shift of CM_Z ----------------------------------------------------- Solver parameters: Files/IO ------------------------------------------------------: - Restart-data prefix: (none) Boundary mapping filename: crm_test.para Output files prefix: crm_wb Profile output description file: (none) Plane output description file: (none) Surface output description file: crm_test.para Field output description file: crm_test.para : - Automatic parameter update (0/1): 1 Automatic parameter update mode (0/1): 0 Accumulate queue time (0/1): 1 Output level: 5 Timestepping Start/Stop ---------------------------------------: - Output period: 5000 Maximal time step number: 0 Minimum residual: 1e-10 Memory management ---------------------------------------------: - Increase memory (0/1): 1 Geometry ------------------------------------------------------: - Grid scale: 1 Reference relation area: 191.84 Reference length (pitching momentum): 7.0005 Reference length (rolling/yawing momentum): 7.0005 Origin coordinate x: 33.68 Origin coordinate y: 11.91 Origin coordinate z: 4.52 Perfect gas thermodynamic -------------------------------------: - Gas constant R: 287 Gas constant gamma: 1.4 Transport coefficients ----------------------------------------: - Reynolds number: 5000000 Prandtl number: 0.72 Sutherland constant: 110.4 Sutherland reference viscosity: 1.7874e-05 Sutherland reference temperature: 273 References ----------------------------------------------------: - Reference density: 0.0469987660369619 Reference temperature: 310.928 Reference pressure: 4193.99767765972 Reynolds length: 7.0005 Reference Mach number: 0.85 Reference velocity: 300.437271096647 Reference outer pressure: 4193.99767765972 Variables -----------------------------------------------------: - Fix negative values (0/1): 0 Minimal density: 1e-12 Minimal pressure: 1e-12 Minimal energy: 1e-12 Infinite swept wing with 2D offset vector (0/1): 0 Flowfield initialization --------------------------------------: - Init total conditions (0/1): 0 Solver IO parameters ------------------------------------------: - Write pointdata dimensionless (0/1): 0 Reference system of forces and moments (tau/ln9300): tau Output scaling factor of translation (xyz): 1 1 1 Universal -----------------------------------------------------: - Solver type: Flow Monitoring ----------------------------------------------------: - Monitor history (0/1): 1 Monitoring values: Rrho_Rrhogrsm_Rrhor11_Rrhor12_Rrhor13_Rrhor22_Rrhor33_Rrhor23_C-lift_C-drag_C-my_Angle-a_Res-lift_X-max-Rrho_Y-max-Rrho_Z-max-Rrho Residuals to normalize: Density Global monitoring (0/1/2): 2 Residual evaluation for monitoring (0/1): 1 Recompute freestream residual (0/1): 0 Flux-consistent force integrals (0/1): 0 Monitoring significant figures: 8_8_8_8_8_4_8_8_4_4_4_4_4_4 Extended coefficient monitoring (0/1): 0 Logfile control -----------------------------------------------: - Enable logfile output on all domains (0/1): 1 Flux main -----------------------------------------------------: - Inviscid flux discretization type: Central Viscous flux type TSL/Full (0/1): 1 Flux parameter blending (0/1): 0 Central flux --------------------------------------------------: - Central dissipation scheme: Matrix_dissipation Central convective meanflow flux: Average_of_flux Scaling of central scheme dissipation: Standard Central convective turbulence flux: Roe2nd 2nd order dissipation coefficient: 0.5 Inverse 4th order dissipation coefficient: 64 Version of cell stretching coefficient: HyperFlex Use modified dissipation for 2D (0/1): 0 Use temperature sensor (0/1): 0 Temperature switch weighting factor: 2 Matrix Dissipation --------------------------------------------: - Pressure switch weighting factor: 1 Minimum artificial dissipation for acoustic waves: 0.2 Minimum artificial dissipation for velocity: 0.2 Couple mean flow and turbulence equations (0/1): 0 Hartens entropy fix (0/1): 0 FILE 'parameters.c', LINE 2046: WARNING: Old parameter 'Matrix dissipation terms coefficient' found in parameter file! Parameter has been removed. See README for replacement. Limiter -------------------------------------------------------: - Limiter freezing convergence: 0 Venkatakrishnan limiter constant: -1 SRR limiter active (0/1): 0 SRR limiter radius relaxation constant: 0.03 Mach number limit for limiter: 0 Flux carbuncle switch version: 0 Gradients -----------------------------------------------------: - Reconstruction of gradients: Green_Gauss Relaxation ----------------------------------------------------: - Relaxation solver: Backward_Euler Hold static velocity field (0/1): 0 Backward Euler ------------------------------------------------: - Linear solver: Lusgs Implicit overrelaxation omega: 1 Implicit overrelaxation beta: 1 Turbulence model equations eigenvalue correction: 0.15 LUSGS ---------------------------------------------------------: - Sgs stages maximum: 3 Lusgs treat rotating frame of reference source terms implicitly (0/1): 0 Multigrid -----------------------------------------------------: - MG description filename: sg Multigrid indicator (0/1): 0 SG start up steps (fine grid): 1000 SG start up steps (fine grid) in initial unsteady phase: 0 Perform SG start up steps (fine grid) every unsteady timestep (0/1): 0 Turbulence equations use multigrid (0/1): 0 Use source terms in coarse grid turbulence equations (0/1): 0 Use face gradient correction on coarse grid (0/1): 0 Set boundary fluxes for forcing function (0/1): 1 Skip dissipation update for forcing function (0/1): 0 Full multigrid ------------------------------------------------: - Multigrid start level: 1 Maximal time step number (coarse grids): 50 Minimum residual (coarse grids): 0.0001 Full multigrid central scheme first-order (0/1): 1 Timestepsize --------------------------------------------------: - CFL number: 5 CFL number (coarse grids): 4 CFL number (large grad p): 5 Limit for residual-based time step limitation: 0.1 Residual-based time step limitation: 0 Residual-based time step limitation factor: 5 Time step smoothing factor: 0 Smoother ------------------------------------------------------: - Residual smoother: Point_explicit Correction smoother: Point_explicit Correction smooth epsilon: 0.2 Residual smooth epsilon: 0.2002 Correction smoothing steps: 2 Residual smoothing steps: 2 Implicit Smoother ---------------------------------------------: - Initial CFL number: 5 Final CFL number: 25 Start up steps with initial CFL number: 10 Increase factor for CFL number: 2 Relaxation coefficient for implicit systems: 1 Evaluate preconditioner first stage only (0/1): 0 Threshold parameter for convective waves: 1 Linear solution method for implicit smoother: Block_symmetric_gauss_seidel Number of sweeps for iterative solution method: 5 Solve turbulence equations weakly coupled (0/1): 0 Ratio solving steps turbulence/mean: 5 Initial CFL number turbulence: 5 Final CFL number turbulence: 25 Start up steps with initial CFL number turbulence: 10 Increase factor for CFL number turbulence: 2 Relaxation coefficient for implicit systems turbulence: 1 Evaluate preconditioner first stage only turbulence (0/1): 0 Linear solution method for implicit smoother turbulence: Symmetric_gauss_seidel Number of sweeps for iterative solution method turbulence: 5 Preconditioning -----------------------------------------------: - Preconditioning: (none) Cut-off value: 4 Spatial cut-off value: 4 Turbulence ----------------------------------------------------: - Turbulence mode: RANS Ratio Prandtl lam/turb: 0.8 General ratio mue-t/mue-l: 0.1 Maximum limit mue-t/mue-l: -1 General turbulent intensity: 0.001 Reference bl-thickness: 1e+22 Turbulence model version: RSM Positivity scheme: 0 EARSM expansion order: 1 Vortical flow correction (0/1): 0 Turbulence initialization type: direct Vortex generator modeling (0/1): 0 RSM -----------------------------------------------------------: - Rsm re-distribution model: SSG/LRR-g.2010 Rsm implementation version: FLOWer Rsm dissipation model: isotropic Rsm diffusion model: GGDH.2010 Rsm length scale equation: Menter_BSL_g Rsm normal stress limiting factor: 0 Rsm omega limiting factor: 1e-05 Rsm epsilon limiting factor: 1e-05 Rsm g limiting factor: 1e-08 Use RSM g minimum limiter (0/1): 1 Rsm transition treatment: Point_transition Rsm transition trip update period: 100 Rsm wall boundary condition type: smooth Cauchy convergence control ------------------------------------: - Use Cauchy convergence control: off Cauchy convergence control variables: (none) Error for Cauchy convergence control: 0 Number of samples for Cauchy convergence: 20 Factor for dynamic Cauchy convergence: 0.05 Chimera -------------------------------------------------------: - 2D offset vector (0 / x=1,y=2,z=3): 2 Primary grid filename: tau.grid Suppress error on orphaned points (0/1): 0 Apply Chimera nearest neighbor search (0/1): 0 Extend valid Chimera nearest neighbors (0/1): 0 Smooth interpolation data for Chimera nearest neighbors (0/1): 0 Apply Chimera wall projection (0/1): 0 Mismatch of overlapping walls: 0 Remove overlap for loads integration (0/1): 0 Surface element tolerance: 0 Surface kink detection angle: 30 Surface kink detection angle tolerance: 15 Chimera new implementation (0/1): 1 Apply Chimera automatic holecutting (0/1): 0 Apply Chimera automatic shrink overlap (0/1): 0 Chimera automatic shrink maximum value: 15 Chimera automatic shrink overlap width: 5 Apply Chimera automatic forced surface blanking (0/1): 0 Dual time -----------------------------------------------------: - Unsteady time stepping: dual Unsteady allow external control over progress in time (0/1): 0 Unsteady show pseudo time steps (0/1): 1 Unsteady physical time step size: 0.0005 Unsteady physical time offset: 0 Unsteady computational time step size: -1 Unsteady normalized time step size: -1 Unsteady physical time steps: 5 Unsteady inner iterations per time step: 100 Minimum number of inner iterations per time step: 0 Unsteady implicit scheme order: 2 Unsteady extrapolation order: 0 Compute harmonics on surface (0/1): 0 Moving grid ---------------------------------------------------: - Type of grid movement: static Motion description filename: (none) Motion hierarchy filename: (none) Extended motion monitoring (0/1): 0 Geometric conservation law (0/1): 1 Initialize deformation (0/1): 0 Flow time averaging -------------------------------------------: - Compute flow statistics: (none) Reset flow statistics: (none) Moving average sampling window size: 0 Extra field output values -------------------------------------: - Field output values: cp_mach_rho_vort_mean_meanvelgrad_vxyzgeod_l2_Q_Hn Surface output parameters -------------------------------------: - Surface output values: xyz_rho_cp_p_cfxyz_cf_yplus_mean-cp_mean-cf_mean Surface output period: 5000 Memory management ---------------------------------------------: - Increase memory (0/1): 1 APSIM parameters ----------------------------------------------: - APSIM input file: (none) APSIM output prefix: (none) Nearest wall data parameters ----------------------------------: - Compute nearest wall data (0/1): 0 Project wall distance: 0 Primary grid filename: tau.grid Transition module parameters ----------------------------------: - Transition module description file: crm_test.para Set transition at solver start (0/1): 0 Transition prescription (0/1): 0 Transition prediction (0/1): 0 Solver boundary parameter: ---------------------------------------------------------------: - Type: symmetry plane Markers: 1 Write surface data (0/1): 1 Cutting plane allowed (0/1): 0 ---------------------------------------------------------------: - Type: farfield Markers: 2 Name: farfield Unified flux formulation (0/1): 1 : - Constant alpha/clift (0/1): 0 Targeted clift: 0 Adjustment factor: 0.01 Lift iteration period: 20 Lift iteration start: 200 : - Vortex correction (0/1): 0 Chord length: 0 : - Monitor farfield state (0/1): 1 Monitor mass flow (0/1): 0 Monitor momentum (0/1): 0 Set gradients (0/1): 0 Write surface data (0/1): 0 Cutting plane allowed (0/1): 0 Modify farfield file: (none) : - Density: 0.0469987660369619 Mach number: 0.85 Reference Mach number for turbulence: 0.85 Velocity: 300.437271096647 Angle beta (degree): 0 Angle alpha (degree): 6 Sideslip angle (degree): 0 Temperature: 310.928 Turbulent intensity: 0.001 Ratio mue-t/mue-l: 0.1 ---------------------------------------------------------------: - Markers: 3 Type: viscous wall Subtype: turbulent Name: FUSELAGE : - Use wall function (0/1): 0 Temperature filename: (none) : - Moving wall (0/1): 0 Moving wall omega: 0 0 0 Moving wall trans: 0 0 0 Moving wall origin: 0 0 0 : - Monitor forces (0/1): 1 Monitor farfield state (0/1): 0 Monitor mass flow (0/1): 0 Monitor heat flow (0/1): 0 Set gradients (0/1): 0 Write surface data (0/1): 1 Cutting plane allowed (0/1): 1 Monitor drag forces (0/1): 0 Monitor lift forces (0/1): 0 Calculate fluxes (0/1): 0 : - Apply mass flux (0/1): 0 Effusion mass flux: 0 Coolant reservoir temperature: -310.928 : - Heat flux: adiabatic Heat flux evaluation: gradient Heat flux correction (0/1): 0 Temperature: 310.928 Prescribed heat flux: 0 Emissivity: 0.8 Reservoir temperature: 310.928 Conduction coefficient: 0 ---------------------------------------------------------------: - Markers: 4 Type: viscous wall Subtype: turbulent Name: FAIRING : - Use wall function (0/1): 0 Temperature filename: (none) : - Moving wall (0/1): 0 Moving wall omega: 0 0 0 Moving wall trans: 0 0 0 Moving wall origin: 0 0 0 : - Monitor forces (0/1): 1 Monitor farfield state (0/1): 0 Monitor mass flow (0/1): 0 Monitor heat flow (0/1): 0 Set gradients (0/1): 0 Write surface data (0/1): 1 Cutting plane allowed (0/1): 1 Monitor drag forces (0/1): 0 Monitor lift forces (0/1): 0 Calculate fluxes (0/1): 0 : - Apply mass flux (0/1): 0 Effusion mass flux: 0 Coolant reservoir temperature: -310.928 : - Heat flux: adiabatic Heat flux evaluation: gradient Heat flux correction (0/1): 0 Temperature: 310.928 Prescribed heat flux: 0 Emissivity: 0.8 Reservoir temperature: 310.928 Conduction coefficient: 0 ---------------------------------------------------------------: - Markers: 5 Type: viscous wall Subtype: turbulent Name: WING : - Use wall function (0/1): 0 Temperature filename: (none) : - Moving wall (0/1): 0 Moving wall omega: 0 0 0 Moving wall trans: 0 0 0 Moving wall origin: 0 0 0 : - Monitor forces (0/1): 1 Monitor farfield state (0/1): 0 Monitor mass flow (0/1): 0 Monitor heat flow (0/1): 0 Set gradients (0/1): 0 Write surface data (0/1): 1 Cutting plane allowed (0/1): 1 Monitor drag forces (0/1): 0 Monitor lift forces (0/1): 0 Calculate fluxes (0/1): 0 : - Apply mass flux (0/1): 0 Effusion mass flux: 0 Coolant reservoir temperature: -310.928 : - Heat flux: adiabatic Heat flux evaluation: gradient Heat flux correction (0/1): 0 Temperature: 310.928 Prescribed heat flux: 0 Emissivity: 0.8 Reservoir temperature: 310.928 Conduction coefficient: 0 ----------------------------------------------------- Active turbulence model: Reynolds stress model Active re-distribution model: SSG/LRR-g.2010 Coefficient C_1 (SSG): 3.400 Coefficient C_1^(*) (SSG): 1.800 Coefficient C_2 (SSG): 4.200 Coefficient C_3 (SSG): 0.800 Coefficient C_3^(*) (SSG): 1.300 Coefficient C_4 (SSG): 1.250 Coefficient C_5 (SSG): 0.400 Coefficient C_1 (LRR): 3.600 Coefficient C_2 (LRR): 0.520 Active dissipation model: isotropic Active diffusion model: GGDH.2010 GGDH 2-layer Reynolds-stress diffusion: inner layer: C_s,i = 0.068 <=> sigma_s,i = 0.750 outer layer: C_s,o = 0.220 <=> sigma_s,o = 2.444 Active length scale equation: Menter_BSL_g Omega production, inner layer, alpha_w,i: 0.556 Omega production, outer layer, alpha_w,o: 0.440 g destruction, inner layer, beta_w,i: 0.075 g destruction, outer layer, beta_w,o: 0.083 g diffusion, inner layer, sigma_w,i: 0.500 g diffusion, outer layer, sigma_w,o: 0.856 g cross-diffusion, inner layer, sigmad_w,i: 0.000 g cross-diffusion, outer layer, sigmad_w,o: 1.712