B-Spline Solutions of General Euler-Lagrange Equations

: The Euler-Lagrange equations are useful for solving optimization problems in mechanics. In this paper, we study the B-spline solutions of the Euler-Lagrange equations associated with the general functionals. The existing conditions of B-spline solutions to general Euler-Lagrange equations are given. As part of this work, we present a general method for generating B-spline solutions of the second-and fourth-order Euler-Lagrange equations. Furthermore, we show that some existing techniques for surface design, such as Coons patches, are exactly the special cases of the generalized Partial differential equations (PDE) surfaces with appropriate choices of the constants.


Introduction
Computer-aided geometric design (CAGD) is an area where the improvement of surface generation techniques is an ongoing demand.Partial differential equations (PDEs) were introduced as a valuable tool for geometric modeling, since they offer several features from which these areas can benefit [1][2][3].The use of PDEs for shape design is conceptually different to conventional methods such as splines.The basic philosophy behind this method is that shape design is effectively treated as a mathematical boundary-value problem, i.e., shapes are produced by finding the solutions to a suitably chosen PDE that satisfies certain boundary conditions [4,5].
The importance of such a mathematical tool is that almost every physical phenomenon is modeled by a PDE.Now, PDEs have been introduced to many areas, such as computer graphics and animation and manufacturing designs, where they have been capable of solving a variety of problems efficiently [6].It is well known that both the harmonic and biharmonic operators associated with Laplacian and bi-Laplacian equations are widely used in many application areas [7].For example, the harmonic operator is widely used in physical problems, such as gravity, electromagnetism, and fluid flows.Similarly, the biharmonic operator is also associated with a variety of physical problems, such as tension in elastic membranes and the study of stress and strain in physical structures.
A Lagrangian functional defined on the space of smooth patches − → x : We take a functional to be Mathematics 2019, 7, 365; doi:10.3390/math7040365www.mdpi.com/journal/mathematics Minimizing the functional I is equivalent to requiring that the first variation of I is zero which then gives rise to the corresponding Euler-Lagrange equation [8].For instance, the general quadratic functional [9] is where a, b, c are constants and the Euler-Lagrange equation associated with the above general functional is which is a second-order partial differential equation, thereby a general form of the second-order PDEs is where a, b, c, d, e, f are constants and subscripts denote derivatives.The general Euler-Lagrange equations is an effective tool in geometrical modeling.Some of the existing techniques for surface design, such as the Coons patches [10] and harmonic surfaces [11], are the particular cases of the generalized PDE surfaces with appropriate choices of the constants in above equations.
In 1989, a method was presented for approximating surfaces which are the solutions of PDEs [2].Monterde [11] focused on harmonic and biharmonic Bézier surfaces and presented a method of surface generation from prescribed boundaries based on the elliptic partial differential operators.In 2006, Monterde [12] also presented a method for generating Bézier surfaces from the boundary information based on the general 4th-order Euler-Lagrange equations which are a generalization of harmonic and biharmonic Bézier surfaces corresponding to Laplace and the standard biharmonic equations.A sixth-order PDE was presented by Zhang and You [13], which provided enough degrees of freedom not only to accommodate tangent, but also curvature boundary conditions and offered more shape control parameters to serve as user controls for the manipulation of surface shapes.Arnal [14] presented methods for designing triangular Bézier PDE surfaces by solving the second-order and fourth-order PDEs.Wang [15] presented a general 8th-order PDE method to generate Bézier surfaces from the boundary with positions and tangent vector information and extended Monterde's results [11] to tri-harmonic Bézier surfaces.In 2015, Beltran [16] studied polynomial solutions in the Bézier form of the wave equation in dimensions one and two and determined which control points of the Bézier solution at two different times fix the solution.In various literature [17][18][19][20], they have indicated the ways in which, by a suitable choice of PDEs and boundary conditions imposed upon its solutions, surfaces satisfying a wide range of functional requirements can be created.
The work presented in this paper is geared to show the B-spline solutions of the general Euler-Lagrange equations.For this purpose, we study the existing conditions of B-spline solutions and then give a general method to compute the B-spline solutions of the general PDEs.Furthermore, we show that some of the existing techniques for surface design, such as the Coons patches and harmonic surfaces, are the special cases of the generalized PDE surfaces we present here.
The paper is organized as follows.In Section 2, we recall some notations and backgrounds about the Euler-Lagrange equations.In Section 3, we study the existing conditions of B-spline solutions to the Euler-Lagrange equations and the general method to compute the B-spline solutions are developed.As part of this work, the B-spline solutions to the 4th PDEs are discussed in Section 4. In Section 5, we indicate some of the existing methods for surface design are the special cases of the general PDE surfaces.The last Section concludes this paper.

Preliminary
For notational purposes, we assume that a parametric surface patch − → x : We also represent the usual derivatives within the parametric domain such that − → ∂v .We impose the standard harmonic and biharmonic PDE operators respectively acting on these surfaces to be, ∆ − → x = 0 and ∆ 2 − → x = 0 where We study functionals defined on the space of smooth patches − → x : Ω → R 3 , Ω ⊂ R 2 .Given a Lagrangian functional as Equation (1) displayed in the Introduction is the general first-order quadratic form of the above functional.The general second-order quadratic form is where a, b, c, d, e, f ∈ R are constants [12].The associated Euler-Lagrange equation for the general quadratic functionals given above is To simplify our formulation, let us consider the more general functionals which is divergence free [21], i.e., where a, b, c, d, e are constants, whose Euler-Lagrange equation [12] is

B-Spline Solutions of the Second-Order PDEs
In this section, we give the existing conditions of B-spline solutions to the general second-order Euler-Lagrange equations and develop a method to get solutions.The general second-order Euler-Lagrange equations are where a, b and c are constants.Given a topological rectangular set of control points denote the B-spline basis functions of degree p and degree q defined on the knots vector [22], which can be obtained by the Cox-de Boor recursion formula [23], taking N i,p (u) for example, The associated B-spline surface of degree p × q is given by The derivative of a B-spline basis function of degree p is simply a linear combination of two B-spline basis functions of degree p − 1 [24].
For the sake of simplicity, it is assumed that the knot vector is uniform spaced, namely the interval [u 0 , u m+p+1 ] is partitioned into m intervals of spacing h 1 .Similarly, the interval [v 0 , v n+q+1 ] is partitioned into n intervals of spacing h 2 .Thus, starting from Equation ( 2), substituting the B-spline surface − → x (u, v) into the above PDE, we can obtained an equation with the control points P i,j to be determined.
We can rewrite Equation ( 3) by the Cox-de Boor recursion formula as follows where Q i,j is a linear combination of the control points obtained by the recursion formula.The aim of the method is to find a B-spline solution where − → x (u, v) is the solution of PDE (2), P i,j are control points to be determined and N i,p (u) and N j,q (v) are the tensor-product B-spline basis functions.Since the degree elevation of B-spline basis functions is more complicated compared to the Bernstein basis functions, the expressions of Q i,j is complicated in the most cases.Let us look at it from another point of view.Obviously, − → x (u, v) is a piecewise polynomial with total degree p + q.Therefore, the degree of − → x (u, v) uu , − → x (u, v) uv or − → x (u, v) vv is p + q − 2 with respect to u and v. Collecting the terms of each degree, there are finite terms, specifically, we have the following results.Lemma 1.In general, there are pq + p + q − 2 terms about u and v in Equation (2).
Proof.For − → x (u, v) uu , the largest degrees of u and v are p − 2 and q respectively.Then there are (p Similarly, for − → x (u, v) uv , the largest degree of u and v are p − 1 and q − 1 respectively, thus there are (p − 1 + 1)(q − 1 + 1) terms in − → x (u, v) uv , but some of them are the same as the terms in − → x (u, v) uu .Actually, there are only q terms are new for Equation (2), to be specific, which are In the same way, − → x (u, v) vv totally has (p + 1)(q − 2 + 1) terms, in which q − 1 terms are new for (2), specifically, u p v q−2 , • • • , u p v, u p .The above analysis suggests that there are (p − 2 + 1)(q + 1) + q + q − 1 terms in Equation ( 2), and then the lemma holds.
Theorem 1.In general, while n ≥ p and m ≥ q, PDE (2) has B-spline solutions.
and only if the coefficient of each term about u and v vanishes.By Lemma 1, there are pq + p + q − 2 terms about u and v in Equation ( 2), which implies that we get pq + p + q − 2 linear equations about the control points.B-spline surface has (n + 1)(m + 1) control points to be determined, while n ≥ p and m ≥ q, (n + 1)(m + 1) > pq + p + q − 2. This suggests that the pq + p + q − 2 linear equations are solvable, that is to say, the PDE has B-spline solutions.
From the previous discussion, we also get a general method to solve the PDEs Equation ( 2) with B-spline solutions.Actually, as the coefficients of terms about u and v vanish, pq + p + q − 2 linear equations of the control points are obtained.While n ≥ p and m ≥ q, the control points to be determined are redundancy, this means we can impose some boundary conditions on the equation.For instance, the values of B-spline surface − → x (u, v) on the boundary and its normal derivative − → x (u, v) u and − → x (u, v) v at the corner points are specified.Thus, we conclude the method for computing the B-spline solutions of Equation ( 2) can be described as Algorithm 1 and explain the method in Example 1.

Input:
Euler-Lagrange Equation (2) and B-spline surfaces − → x (u, v) with the control points to be determined; Output: The B-spline solutions − → x (u, v) of Equation (2).

Step 1:
Substituting the B-spline surface − → x (u, v) into Equation (2), collect the terms about u and v.

Step 2:
Collecting the pq + p + q − 2 coefficients of each terms of equations, give pq + p + q − 2 linear equations of the control points {P i,j , 0 ≤ i ≤ n, 0 ≤ j ≤ m}.

Step 3:
Combining with the given boundary conditions or control points and solving the linear equations, the unknown control points are obtained.

Step 4:
Construction of B-spline surface − → x (u, v) by the control points.
Example 1.For our purposes, we require that the basis functions be not only B-splines, but also at least two-times differentiable.This implies that the degree of B-spline surface is at least three.Thus, as it is not desirable to add any more complication than is strictly necessary, cubic B-splines will be used as the basis functions throughout the discussion.In the interest of simplicity, assuming that the B-splines are uniformly spaced, i.e., the interval [u 0 , u m+p+1 ] is partitioned into m intervals of spacing h.Then at each point u i is centered a piecewise cubic B-spline N i,3 (u) as N i,3 (v) can be given in the same way.By the Theorem 1, in order to get the B-spline solutions of Equation (2), the control points set satisfies m ≥ 3 and n ≥ 3. Assume m = 3, n = 3, the B-spline surface to be determined can be represented by where P i,j are the unknown control points.
For simplicity, assuming a = 1, b = −2, c = 1 in Equation (2) and the space of interval h = 1.Substituting the B-spline surface − → x (u, v) into Equation (2), and collecting the terms about u and v, we get 13 terms in the equation, specifically, u 3 ) holds if and only if the coefficients of these terms vanish, which implies the following 13 equations about control points holds.
−2 P 0,0 + 6 P 0,1 − 6 P 0,2 + 2 P 0,3 + 6 P 1,0 − 18 P 1,1 + 18 P 1,2 − 6 P 1,3 − 6 P 2,0 + 18 P 2,1 − 18 P 2,2 + 6 P 2,3 + 2 P 3,0 In fact, since the coefficients of uv 3 and u 3 v are identical, the above linear system has 12 linear equations about the control points.Since m = 3, n = 3, there are 16 control points to be determined which implies there are 4 control points are free.We impose the additional values of B-spline surface or the control points on the boundary, and then we can obtain the unknown control points unique by solving the above linear equations.Supposing by solving the above equations, the remaining control points are given as

B-Spline Solutions of the 4th Order PDEs
In this section, we represent the solutions of the following fourth-order PDE with B-spline surfaces, which can be discussed using the proposed method.
To get the B-spline solutions of the above PDE, we can establish the algorithm in the similar fashion with the two-order PDEs in Section 3. Substituting the B-spline surface − → x (u, v) into the above PDE, an equation with the control points P i,j can be obtained.Obviously, − → x (u, v) is a piecewise polynomial surface with total degree p + q.Therefore, the total degree of − → x (u, v) uuuu , − → x (u, v) uuuv or − → x (u, v) uuvv is p + q − 4 with respect to u and v.We give the following results analogous to the two-order PDEs.Lemma 2. In general, there are pq + p + q − 9 terms about u and v in Equation (4).
By the proposed lemma and analysis, the following result can be given, and the proof is similar with Theorem 1, which will be omitted here.
Theorem 2. In general, while m ≥ p and n ≥ q, PDE (4) has B-spline solutions.
Thus, in similar with the method to get the B-spline solutions of the second-order PDE in Section 3, we conclude the method for computing the B-spline solutions of Equation ( 4) as Algorithm 2 and explain the method in Example 2.
Algorithm 2 Algorithm to get the B-spline solution of the Euler-Lagrange Equation ( 4)

Input:
B-spline surfaces − → x (u, v) with the control points to be determined and the general Euler-Lagrange Equation (4); Output: The B-spline solutions − → x (u, v) of Equation (4).

Step 1:
Substituting the B-spline surface − → x (u, v) into Equation ( 4), collect the terms about u and v.

Step 2:
Collecting the pq + p + q − 9 coefficients of each terms of equations, give pq + p + q − 9 linear equations of the control points {P i,j , 0 ≤ i ≤ n, 0 ≤ j ≤ m}.

Step 3:
Combining with the given boundary conditions or control points and solving the linear equations, the unknown control points are obtained.

Step 4:
Construction of B-spline surface − → x (u, v) by the control points.
Example 2. In a similar way, the fourth-order PDE implies that the B-spline solution is at least four-times differentiable.Therefore, we require the degree of B-spline surface is at least quintic.As our aim is that the solution be obtained in the form of a tensor-product B-spline surface, and then N i,p (u) is chosen to be quintic B-spline basis functions which span the space of piecewise quintic polynomials.Assuming that the knot vector of B-splines is uniformly spaced, and the space of the intervals is h.In the interest of simplicity, suppose h = 1.Then at each point u i is centered a piecewise quintic B-spline N i,5 (u) as

Applications in Geometrical Modeling
As mentioned previously, the PDE surfaces arising from the functional serve as a useful tool in geometrical modeling.In this section, we discuss some existing surfaces, such as Coons surfaces and biharmonic surface, are exactly the special cases of the general PDE surfaces.
While dealing with biharmonic surfaces [25], we often study the functionals associated with the standard biharmonic operators ∆ 2 , As the definition of the biharmonic operator itself suggests a kind of control on the shape of the resulting surface in some sense.
Other similar functionals used in CAGD are also the special cases of the previous general functional.Among them we cite the functional related to Coons patches which can be used to show that the resulting surface minimizing the surface twist [10] The modified biharmonic functional is [26] A more general form is where a, b and c are constants.This functional is applied to minimize the resulting surface area where more weights are assigned along one particular parametric coordinate direction.The examples of functionals displayed above are some of the typical functionals commonly used in CAGD.It is noteworthy that they are the special cases of the most general functionals, Please note that above surface design techniques, such as the biharmonic surface and Coons surface, are exactly the special cases of the generalized framework with appropriate choices of the constants.

Conclusions
PDEs have been introduced to many areas such as computer graphics, CAGD, and animation since they can solve a variety of problems in these areas.In this paper, we consider the B-spline solutions of the general Euler-Lagrange equations.It is well known that both the harmonic and biharmonic operators associated with Laplacian and bi-Laplacian equations are widely used in many application areas.We use a linear combination of these common operators and get a wide variety of PDEs which can be used in CAGD.The existing conditions of B-spline solutions to a general PDE are given.As part of this work, we present a general method for generating B-spline solutions to the second-and fourth-order Euler-Lagrange equations.
There are several extensions of this work.It is interesting to study how the various coefficients associated with the Euler-Lagrange equations affect the shape of the resulting surfaces.We could also carry out a study of the geometric properties of surfaces generated by the PDEs which will be useful in the developing interactive tools in future.

Figure 1 .
Figure 1.B-spline solutions of the second-order PDEs.

Figure 2 .
Figure 2. B-spline solutions of the fourth-order PDEs.