Abstract: | This paper presents a p-version least-squares finite element formulation for unsteady fluid dynamics problems where the effects of space and time are coupled. The dimensionless form of the differential equations describing the problem are first cast into a set of first-order differential equations by introducing auxiliary variables. This permits the use of C° element approximation. The element properties are derived by utilizing p-version approximation functions in both space and time and then minimizing the error functional given by the space–time integral of the sum of squares of the errors resulting from the set of first-order differential equations. This results in a true space–time coupled least-squares minimization procedure. A time marching procedure is developed in which the solution for the current time step provides the initial conditions for the next time step. The space–time coupled p-version approximation functions provide the ability to control truncation error which, in turn, permits very large time steps. What literally requires hundreds of time steps in uncoupled conventional time marching procedures can be accomplished in a single time step using the present space–time coupled approach. For non-linear problems the non-linear algebraic equations resulting from the least-squares process are solved using Newton's method with a line search. This procedure results in a symmetric Hessian matrix. Equilibrium iterations are carried out for each time step until the error functional and each component of the gradient of the error functional with respect to nodal degrees of freedom are below a certain prespecified tolerance. The generality, success and superiority of the present formulation procedure is demonstrated by presenting specific formulations and examples for the advection–diffusion and Burgers equations. The results are compared with the analytical solutions and those reported in the literature. The formulation presented here is ideally suited for space–time adaptive procedures. The element error functional values provide a mechanism for adaptive h, p or hp refinements. The work presented in this paper provides the basis for the extension of the space–time coupled least-squares minimization concept to two- and three-dimensional unsteady fluid flow. |