Least-squares scheme

For simplicity we assume that = (d /dx ) + a(d/dx), where o is a constant. The basic steps in the least squares scheme are ... 

Working equations of the least-squares scheme in Cartesian coordinate systems... 

Retaining all of the terms in the w eight function a least-squares scheme corresponding to a second-order Petrov-Galerkin formulation will be obtained. 

Regularization methods, to invert Equation 11.31 numerically by a modified least-squares scheme... 

The standard least-squares approach provides an alternative to the Galerkin method in the development of finite element solution schemes for differential equations. However, it can also be shown to belong to the class of weighted residual techniques (Zienkiewicz and Morgan, 1983). In the least-squares finite element method the sum of the squares of the residuals, generated via the substitution of the unknown functions by finite element approximations, is formed and subsequently minimized to obtain the working equations of the scheme. The procedure can be illustrated by the following example, consider... 

The basic procedure for the derivation of a least squares finite element scheme is described in Chapter 2, Section 2.4. Using this procedure the working equations of the least-squares finite element scheme for an incompressible flow are derived as follows ... 

Extension of the streamline Petrov -Galerkin method to transient heat transport problems by a space-time least-squares procedure is reported by Nguen and Reynen (1984). The close relationship between SUPG and the least-squares finite element discretizations is discussed in Chapter 4. An analogous transient upwinding scheme, based on the previously described 0 time-stepping technique, can also be developed (Zienkiewicz and Taylor, 1994). 

Nguen, N. and Reynen, J., 1984. A space-time least-squares finite element scheme for advection-diffusion equations. Cornput. Methods Appl Mech. Eng. 42, 331- 342. 

Least-square.s and streamline upwind Petrov-Galerkin (SUPG) schemes... 

In steady-state problems 6/S.l = 1 and the time-dependent term in the residual is eliminated. The steady-state scheme will hence be equivalent to the combination of Galerkin and least-squares methods. 

Sample data from the literature18 are shown in Fig. 2-10. The curve shows the least-squares fit. A further development of this scheme is presented in Section 3.5. 

Toward these ends, the kinetics of a wider set of reaction schemes is presented in the text, to make the solutions available for convenient reference. The steady-state approach is covered more extensively, and the mathematics of other approximations ( improved steady-state and prior-equilibrium) is given and compared. Coverage of data analysis and curve fitting has been greatly expanded, with an emphasis on nonlinear least-squares regression. 

Firstly, it has been found that the estimation of all of the amplitudes of the LI spectrum cannot be made with a standard least-squares based fitting scheme for this ill-conditioned problem. One of the solutions to this problem is a numerical procedure called regularization [55]. In this method, the optimization criterion includes the misfit plus an extra term. Specifically in our implementation, the quantity to be minimized can be expressed as follows [53] ... 

Different baseline correction methods vary with respect to the both the properties of the baseline component d and the means of determining the constant k. One of the simpler options, baseline ojfset correction, nses a flat-line baseline component (d = vector of Is), where k can be simply assigned to a single intensity of the spectrum x at a specific variable, or the mean of several intensities in the spectrum. More elaborate baseline correction schemes allow for more complex baseline components, such as linear, quadratic or user-defined functions. These schemes can also utilize different methods for determining k, such as least-squares regression. 

---