Parallel-in-time disturbance region update method for dynamic characteristics of flight vehicles

A parallel-in-time disturbance region update method for dynamic characteristics of flight vehicles, comprising:

• step 1:
• initializing the computation;
• step 1-1:
• reading-in data:
• reading-in coordinates of each vertex and boundary conditions of a grid, final-state time, number of physical time periods and physical time layers;
• step 1-2:
• obtaining components of Chebyshev transformation matrix:
• computing Chebyshev transformation matrix based on the number of physical time layers, thereby obtaining the components of Chebyshev transformation matrix;
• step 1-3:
• initializing flow field:
• computing grid cell parameters and assigning initial values to the conservative variables of all grid cells on the physical time layers;
• Initializing the flow field from freestream condition or from specified flow field;
• step 1-4:
• establishing a dynamic computational domain based on the initialization of the flow field;
• establishing an advective dynamic computational domain and a viscous dynamic computational domain;
• step 2:
• Solving the flow-governing equations of all time layers simultaneously in the dynamic computational domains;
• step 2-1:
• allocating memory;
• step 2-2:
• treating boundary conditions;
• step 2-3:
• residual estimation;
• Estimating the residual term of a flow-governing equations in the advective dynamic domain of each physical time layer;
• step 2-4:
• time integration;
• In the advective dynamic domain of each physical time layer, computing the correction of conservation variables of a current time layer in the current iteration step:
• step 3:
• updating the dynamic computational domain;
• step 3-1:
• extending the advective dynamic domain;
• step 3-2:
• contracting the advective dynamic domain;
• step 3-3:
• extending the viscous dynamic domain;
• step 3-4:
• contracting the viscous dynamic domain;
• step 4:
• judging whether the parallel computation in the current time period converges;
• Whether the iteration of the current time period converges can be judged in terms of the maximum 2-norm of the corrections of the conservation variables over all grid cells in the dynamic computational domain of all time layers;
• if the maximum 2-norm is lower than or equal to a preassigned threshold, which means that the iteration of the current time period converges, entering step 5, and if not, entering step 2 and performing the next iteration of the current time period;
• step 5:
• judging whether the computation completes;
• If maximum time of the current time period reaches the final-state time, outputting the computation results, and if not, entering step 1-3 and performing computation of the next time period.