Miklós HOFFMANN, Imre JUHáSZ, Ede TROLL

1Institute of Mathematics and Computer Science, Eszterházy Károly University, Eger 3300, Hungary

2Department of Computer Graphics and Image Processing, University of Debrecen, Debrecen 4010, Hungary

3Department of Descriptive Geometry, University of Miskolc, Miskolc 3515, Hungary

Abstract: While considering a mirror and light rays coming either from a point source or from infinity,the reflected light rays may have an envelope, called a caustic curve. In this paper, we study developable surfaces as mirrors.These caustic surfaces, described in a closed form, are also developable surfaces of the same type as the original mirror surface. We provide efficient, algorithmic computation to find the caustic surface of each of the three types of developable surfaces (cone, cylinder, and tangent surface of a spatial curve). We also provide a potential application of the results in contemporary free-form architecture design.

Key words: Caustics; Developable surface; Reflected light rays; Curve of regression

1 Introduction

When light rays are reflected from a curved mirror, the following optical phenomenon may be observed: the reflected light rays may possess an envelope,called a caustic curve or surface. These caustics can not only appear in our everyday experience,such as on the surface of coffee in our coffee cup, but also play an important role in the sciences, from physics to computer graphics(Arnold et al., 1985;Lock and Andrews, 1992). Scientists, from the ancient Greeks through Huygens to contemporary opticians, engineers, and geometers, have studied reflected and refracted light rays from the theoretical point of view as well as through various applications.

In recent years,these surfaces have been usually described as parametric free-form surfaces,typically as Bézier or B-spline surfaces in various applications,to provide more freedom for users in the interactive design process (Tang et al., 2016). Optical studies frequently apply these free-form surfaces in lens design (Liu P et al., 2012; Ponce-Hernández et al.,2020)or light-emitting diode(LED)illumination research(Wu et al.,2013).

In engineering and architecture,this problem is especially relevant in terms of developable surfaces,i.e., curved surfaces that can be unfolded to (and therefore created from)a planar shape. There are basically three different types of developable surfaces.Two of them are the well-known cone and cylinder,while the third one is a more general type, namely,the tangent surface of spatial curves. From the computational point of view, this latter type is the most challenging one in engineering applications, but at the same time,this type provides much more freedom in engineering design than the classical cones and cylinders (Seguin et al., 2021). These surfaces are used, among other applications, for creating special developable mechanisms,e.g.,for cylindrical(Greenwood et al.,2019)and for conical(Hyatt et al.,2020).Non-developable surfaces are also frequently approximated by developable ones, for instance, to ease fabrication from sheet metal (Liu XH et al., 2016).For an excellent overview of developable Bézier surfaces, readers can refer to Zhang and Wang (2006).In this paper, we also follow this construction.

One of the most spectacularly evolving application of developable surfaces can be found in architecture, precisely, in the so-called free-form architecture. Due to evident mechanical and material restrictions,special attention has been paid to developable surfaces in this field (Pottmann et al., 2015;Martín-Pastor,2019).

The caustics of classical planar curves are wellknown and widely studied (Yates, 1947; Lockwood,1967). In the case of surfaces,theoretical results are also known. For a given surface,a somewhat similar notion is the focal surface,i.e.,the surface formed by the centers of the curvature spheres. The relationship between the focal and caustic surfaces has been established(Pottmann and Wallner,2000)and it has been proven that the focal surface of a developable surface will be developable of the same type as the original one(Pottmann and Wallner,2000). This result,in theory,yields the same consequence in terms of caustic surfaces, but these theoretical outcomes do not provide exact, constructive, and algorithmic solutions to compute and display these surfaces in practical applications. To calculate the caustics of a given surface, numerical solutions have been provided(Schwartzburg et al., 2014).

Instead of numerical calculation, in this paper,we provide the exact computation and closed formulae for the caustics of developable surfaces. These caustic surfaces are of utmost importance in contemporary architecture (Pottmann et al., 2015), where caustics may appear, e.g., as an outcome of the reflected sunshine beams.

2 Developable surfaces as mirrors

As is well-known,considering the curve

and directiong(t),t∈[a,b], the surface

is a ruled one. The given curve is the directrix of the surface, while for any fixedt0∈[a,b], the linesr(t0)+ug(t0)are called the generators(or rulings).The surface is developable if the normals of the tangent planes along the generators are of constant direction;i.e., considering the partial derivatives

the normal vector

does not depend onufor anyt. In other words,the tangent planes of the surface along its generators coincide.

Now let us consider a developable surface as a mirror. Here, we assume that the light source is point-like, either being at infinity or not, and that none of the generators of the developable surface go through the light source. For any generatorr(t0)+ug(t0) of the surface, incoming light rays meeting the mirror surface along this generator are coplanar,and due to the fixed tangent plane along the generator, the reflected light rays are also coplanar. The plane of reflected light rays is to be computed first.This plane is fully determined by the selected generator and one single reflected light ray. In the following subsections, the computations are presented separately for the case when the light source is at infinity(yielding parallel light rays)and for the case when the source of light is a point of the affine space.

2.1 Light source at infinity

Without loss of generality of the forthcoming computation, we can assume that the direction of parallel light rays isd=[1 0 0]T. We further assume that none of the generators are parallel to the given direction. Let the selected generator of the surface ber(t0)+ug(t0).

The direction(t0) of the reflected rays can be computed by reflectingdwith respect to the tangent plane along the selected generator with normal vectorn(t0). Since the endpoint ofdis on the unit sphere,

the endpoint of the reflected ray(t0) will also be on this sphere. Thus, we have to compute the intersection point of the unit sphere and the line:

whereλ∈R. Solving the equation

one can find

Therefore, the direction of the reflected rays along this generator is

2.2 Light source at an affine point

Without loss of generality of the foregoing computation, we can assume that the source of light is at the origin of the coordinate system. We further assume that none of the generators pass through this point. Again,let the selected generator of the surface ber(t0)+ug(t0).

As we have observed,light rays intersecting this generator form a plane, and the reflected light rays will also form a plane passing through this generator.To determine this plane, it is enough to reflect one single light ray,e.g.,the one intersecting the directrix curver(t) at this generator. In other words, we have to reflect the vectorr(t0) (this is the direction of a light ray coming from the origin, i.e.,r(t0) =d(t0))with respect to the normal vectorn(t0)of the tangent plane.

The endpoint ofr(t0)=d(t0)is on the sphere:and therefore the endpoint of the reflected vector(t0)also has to be on this sphere.

Thus,we have to compute the intersection point of the line passing throughr(t0) with the direction vectorn(t0), i.e., of the line

and sphere(11). Solving the equation

one obtains

The reflected vector is

which, together with the selected generator, determines the plane of the reflected rays.

2.3 The family of planes of reflected rays and their envelope surface

In the preceding subsections,we have computed the vector(t0), and this computation is necessary for determining the plane of reflected rays at each generatorr(t0)+ug(t0)of the developable surface.Forthcoming computations are independent of the actual position of the light source.

The reflected rays form a one-parameter family of planes:

wheret∈[a,b] is the family parameter andu,v∈R. It is easy to see that the envelope surface of this one-parameter family of reflected light planes is also developable (Do Carmo, 2016). This envelope surface is the caustic surface of the original mirror surface.

In what follows, we will determine this caustic surface in the standard form of developable surfaces analogous to Eq. (2), i.e., in the following form:

whereq(t) is the curve andf(t) is the direction of the generators passing through the points of this curve. We apply the computations as

and therefore we have

from which we can obtain

The normal vector of the planeP(u,v,t)is

the derivative of which is

To simplify this expression, we apply the formula of

Accordingly, the two terms on the right-hand side of Eq. (25)can be written as

and

Therefore,the direction of the generators of the envelope surface can simply be written in the following form:

Now,if one can fixuandvin a way that determinant (22) vanishes?t∈[a,b], then by substituting these values ofuandvinto Eq. (16), we obtain a curve (or a constant vector)q(t), based on which the caustic envelope surface can be written in the classic form of developable surfaces in Eq. (17).

3 Caustics of developable surfaces

As is well-known, there are three types of developable surfaces: (generalized)cone, (generalized)cylinder,and the tangent surface of spatial curves. In this section, we specify the general computations in Section 2 to the three different types of developable surfaces.

3.1 Cones as mirrors

Considering a cone as a mirror,defined by curver(t), apexp, and light rays with directiond(t), the caustic surface, i.e., the envelope surface of the reflected light rays with direction(t), is also a cone(Pottmann and Wallner, 2000). Herein, we present the exact computation of the caustic cone surface and provide the solution in a closed form.

In the case of a conic mirror surface, the curver(t)of Eq. (2)is an arbitrary planar or spatial curve,while the generators can be written as

wherep, the apex of the cone, is an arbitrary point out of curver. The family of planes in Eq. (16)will be of the specific form:

Finally, they entered the clothing section and headed for a rack full of Ralph Lauren Polo shirts. Yes! I could already picture myself standing12 in the fairway, contemplating13 my approach shot, while the others in my group commented on my impeccable taste.

wheret∈[a,b] andu,v∈R.

In this case,determinant(22)will be

which vanishes?t∈[a,b] ifu= 1 andv= 0. Substituting these parameter values into Eq. (31), we obtain

That is, the caustic envelope surface is also a cone with apexp. The direction of the generators can be written as

Based on the calculations above,the caustic envelope surface of the cone can be written as

wherepis the apex andf(t) denotes the direction of generators.

An example is shown in Fig.1.

Fig.1 The family of planes of reflected rays(in yellow)and their envelope caustic surface (in orange) in the case of a conic mirror (in blue) defined by a cubic Bézier curve (in red)

3.2 Cylinders as mirrors

Considering a cylinder as a mirror, defined by curver(t), directiona,and light rays with directiond(t),it is known that the caustic surface,i.e.,the envelope surface of the reflected light rays with direction(t), is also a cylinder (Pottmann and Wallner,2000). In this subsection, we present the computation of this caustic surface and provide a closed form of it.

If the mirror surface in Eq.(2)is a cylinder,then curver(t)is an arbitrary planar or spatial curve,and the direction of the generatorsg(t) is constant,i.e.,

whereais a vector (not parallel to the plane ofrwhenr(t)is planar). The family of planes in Eq. (16)of the reflected rays can be written as

The direction of generators of the caustic envelope surface will be

That is,f(t)is parallel toa,?t∈[a,b]. This means that the caustic surface of a cylindrical mirror is also a cylinder,and that the generators of these two cylinders are parallel.

Determinant(22)becomes

from which one can expressvas

Substituting the calculatedvandu= 0 into Eq. (22), it vanishes?t∈[a,b]. Substituting the samevandu= 0 into the family of planes in Eq. (37),the following curve is obtained:

Based on the calculations presented above, the caustic surface of the cylinder can be written as

A cylindrical mirror and its caustic surface are shown in Fig. 2.

Fig.2 The family of planes of reflected rays(in yellow)and their envelope caustic surface (in orange) in the case of a cylindric mirror (in blue)

Fig. 4 The caustic surface (in orange) along with its regression curve in the case of a tangent surface of a spatial curve as a mirror (in blue) defined by a cubic Bézier curve (in red)

3.3 Tangent surface of a spatial curve as a mirror

Considering a general tangent surface of a spatial curver(t) as a mirror and light rays with directiond(t), the caustic surface, i.e., the envelope surface of the reflected light rays with directiond(t),has been proven to be a tangent surface (Pottmann and Wallner,2000),but the exact construction of this surface is not known. Now,we present the computation of the caustic surface and provide the solution in a closed form.

If the developable mirror surface in Eq. (2) is a general tangent surface of a spatial curve, thenr(t)is an arbitrary spatial curve, while the direction of the generators fulfills the relationship as

In this case,determinant(22)is of the form as

and it vanishes?t∈[a,b] whenu=v=0. This immediately yields the somewhat surprising fact that the caustic envelope surface contains the original curver(t).

The normal vector of planeP(u,v,t)is

while its derivative is

In this case, the direction of the generators of the caustic envelope surface can be written as

and the caustic envelope surface is of the form as

It is clear from this expression that the caustic surface is a ruled one. However, our aim is to prove that it is a developable surface,specifically,that it is a tangent developable surface of a spatial curve.

To reach this aim, we have to find the curve of regression,i.e., a curvec(t),t∈[a,b], the tangents of which are parallel to the directionsf(t) of the generators of the surface,?t∈[a,b].

Since each curve on the ruled surface in Eq.(49)can be considered a functional translation of curver(t) along the rulings, we search for the curvec(t)in the following form:

where the functionλ(t), t∈[a,b], is to be determined. Moreover, the derivative of this curve with respect totis

which must be parallel tof(t);i.e., ˙c(t)must be orthogonal to vectors (46) and (47). These conditions yield the following system of equations:

From Eq.(52),we obtain the following equality:

and from Eq. (54),we obtain

which further yields

Therefore,

whenever

holds. Here,we prove that this scalar product cannot be equal to zero.

The first factor of the scalar product in inequality (58) is not vanishing, i.e.,(t)/= 0, since the generators of the surface are not parallel (it is not a cylinder). Thus,(t)can be written in the following form:

Since expression (59) is the linear combination of(t) and(t), this term is parallel to the plane spanned by(t)and(t),which is actually the planeP(u,v,t).

Applying identity (26) to expression (60), one can find

However,this cross product is evidently orthogonal to one of its factors,i.e.,

Therefore, expression (60) is also parallel to the plane(orthogonal to the normal vector of the plane)spanned by(t) and(t), i.e.,P(u,v,t). We obtain that both terms of(t) (expressions (59) and(60)) and(t) are parallel to the planeP(u,v,t)and orthogonal to its normal vector, i.e.,(t)⊥

For the second factor of the scalar product in inequality(58),the equality

holds;therefore,this term describes the change of the normal vector(t)(t)of the planeP(u,v,t).

As a consequence, the scalar product, i.e., the left side of inequality (58), could be equal to zero only if either the second factor, i.e., the derivative of the normal vector(t)(t), is the null vector or the normal vector is parallel to its derivative.Both would yield that the planesP(u,v,t)coincide or that they are parallel?t; i.e., curver(t) is a planar curve. However,this contradicts our assumption that we are studying the tangent surfaces of spatial curves.

Back to the system of Eqs. (52) and (53), we have expressedλ(t) from Eq. (53), but this must fulfill Eq. (52)as well.

In other words,

must hold. However, this equation indeed holds,since we have just seen thati.e., the scalar product above equals zero,?t.

Based on the computation presented above,the caustic surface of a tangent surface can be computed as a tangent surface of the spatial curvec(t)=r(t)+λ(t)f(t), where

and

Examples of this type of mirror and its caustic surface are shown in Figs. 3 and 4,respectively.

Fig.3 The family of planes of reflected rays(in yellow)and their envelope caustic surface (in orange) in the case of a tangent surface of a spatial curve as a mirror(in blue)

4 An application: developable surfaces in architecture

Developable surfaces,especially the tangent surfaces of a curve,are commonly used in modern architecture(Glaeser and Gruber,2007;Pottmann et al.,2015). If the surface of such a construction is made from a reflective material (e.g., metal sheets at the Guggenheim Museum in Bilbao), then the surface may behave like a mirror. When a free-form building or a sculpture reflects light rays, some places around the structure, where light rays are concentrated, can be of high temperature. This especially holds around the cusp of the intersection curve of the caustic surface of the building and the ground plane.Our aim is to find those “hot” points of the ground plane. Without loss of generality, we can assume that this plane is the [x,y] plane.

The intersection of the caustic surfacee(t,λ)and the [x,y] plane is a curve,expressed as follows:

This curve can be described as

whereλ(t)is the solution to the equatione(t,λ)=0 for parametert. This curve usually has a cuspe(t0,λ0),which can be very hot due to physical reasons. Applying standard methods,in most cases,this cusp cannot be directly computed exactly. Based on our results, however, this cusp can simply be computed as the intersection of the curve of regression of the caustic surface and the[x,y]plane(Fig.5),as follows:

Fig. 5 The caustic surface (in orange) and its regression curve in the case of a tangent surface of a spatial curve as a mirror (in blue) defined by a cubic Bézier curve (in red)

5 Conclusions

The caustics of developable surfaces have been studied in this paper. In the theory of developable surfaces, it is known that the caustics of each of the three different types of developable surfaces are developable surfaces of the same type. These caustic surfaces have been expressed in an exact, closed form. In the case of the tangent surface of a spatial curve, the curve of regression of the caustic surface has also been computed. This curve has been applied in a practical computation of finding the cusp of the intersection curve of a caustic surface and a plane,which may appear in applications such as free-form architecture.

Contributors

Miklós HOFFMANN, Imre JUHáSZ, and Ede TROLL contributed equally, on a shared basis to the research, jointly developing the mathematical computations and creating the experimental figures.

Compliance with ethics guidelines

Miklós HOFFMANN, Imre JUHáSZ, and Ede TROLL declare that they have no conflict of interest.

Open access

This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing,adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.