Next Article in Journal
Three-Echelon Supply Chain Management with Deteriorated Products under the Effect of Inflation
Next Article in Special Issue
Remarks on the “Onsager Singularity Theorem” for Leray–Hopf Weak Solutions: The Hölder Continuous Case
Previous Article in Journal
A Review of the Key Technology in a Blockchain Building Decentralized Trust Platform
Previous Article in Special Issue
Numerical Study of a Confined Vesicle in Shear Flow at Finite Temperature
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Drag Coefficients around Objects Created Using Log-Aesthetic Curves

1
Special Interest Group on Modelling & Data Analytics, Faculty of Ocean Engineering Technology and Informatics, University Malaysia Terengganu, Kuala Nerus 21030, Malaysia
2
Graduate School of Engineering, Shizuoka University, 3-5-1 Johoku, Naka Ward, Hamamatsu 432-8011, Japan
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(1), 103; https://doi.org/10.3390/math11010103
Submission received: 7 November 2022 / Revised: 11 December 2022 / Accepted: 13 December 2022 / Published: 26 December 2022
(This article belongs to the Special Issue Mathematical Dynamic Flow Models)

Abstract

:
A fair curve with exceptional properties, called the log-aesthetic curves (LAC) has been extensively studied for aesthetic design implementations. However, its implementation in terms of functional design, particularly hydrodynamic design, remains mostly unexplored. This study examines the effect of the shape parameter α of LAC on the drag generated in an incompressible fluid flow, simulated using a semi-implicit backward difference formula coupled with P 2 P 1 Taylor–Hood finite elements. An algorithm was developed to create LAC hydrofoils that were used in this study. We analyzed the drag coefficients of 47 LAC hydrofoils of three sizes with various shapes in fluid flows with Reynolds numbers of 30, 40, and 100, respectively. We found that streamlined LAC shapes with negative α values, of which curvature with respect to turning angle are almost linear, produce the lowest drag in the incompressible flow simulations. It also found that the thickness of LAC objects can be varied to obtain similar drag coefficients for different Reynolds numbers. Via cluster analysis, it is found that the distribution of drag coefficients does not rely solely on the Reynolds number, but also on the thickness of the hydrofoil.
MSC:
65D17, 68U05, 35Q30, 76-05

1. Introduction

Log-aesthetic curve (LAC), first proposed by Yoshimoto and Harada [1], refers to a family of aesthetically pleasing curves with monotonic curvatures. Harada et al. [2] found that the manufactured and natural shapes that are deemed beautiful have a linear Logarithmic Distribution Diagram of Curvature (LDDC). Thus, Miura [3] derived the LAC fundamental equation by equating the gradient of the Logarithmic Curvature Graph (LCG), the analytical form of LDDC, to a constant α , which is the parameter determining the shape of the LAC. A designer usually interrogates curves using the curvature profiles, which involve second derivatives of curves, whereas LCG involves third derivatives, which are thus suitable for higher-order shape interrogation [4].
LAC has since then been proposed for many applications such as car design [5], modeling transitional curves [6], path planning [7], Computer-Aided Design (CAD) systems [8], and architecture design [9]. The success of LAC is due to its underlying properties; it has a sufficient degree of freedom and shape parameters to represent various spirals [10]. However, there is still much to be discovered about its fluid dynamics properties. It is unknown how fluid flow and drag changes with the change in the LAC’s shape, which is determined by its shape parameter α . To the best of our knowledge, there is not yet any literature that elucidates this relationship.
Takuma et al. [11] studied how curvature affects the energy absorption characteristic of cylindrical corrugated tubes, showing that curve properties such as curvature have more to offer than contributing solely to aesthetic appeals. Likewise, we wish to elucidate how LACs’ shape properties, e.g., LCG gradient and curvature, affect fluid dynamics of incompressible fluid flows. Hence, we constructed hydrofoil-like objects using LAC and investigated their drag coefficient in an incompressible fluid flow using numerical simulation. The information obtained from this study may aid in designing submerged structures, objects, or vehicles, especially in drag reduction.
Lift and drag evaluations were based on the volume integral formulation found in [12,13]. The volume integral formulation has been known to provide better lift and drag coefficient estimations than conventional line integrals along our streamlined-shaped object. Furthermore, volume integral formulation is known to be less sensitive to a slight change in the mesh generated around the object. Recently, [14] reported the behavior of streamlines and fluid flow around streamlined-shaped objects constructed with LAC.
The contribution of this paper is two-fold. This paper completes the work of Wo et al. [14] to report how an LAC’s shape, dictated by its shape parameter, influences the drag of an incompressible fluid flow. An effect of the Reynolds number on the drag coefficient trend is also observed. Furthermore, an algorithm was developed to create LAC that satisfies given length, height, G 1 Hermite data on one end, and G 0 Hermite data on the other. This algorithm is used to construct LAC hydrofoils used in this study.
The rest of the paper is organized as follows. Section 2 discusses the numerical method used to simulate the incompressible non-stationary fluid flow. Section 3 elucidates the creation of an LAC hydrofoil with user-specified thickness and shape parameter α . Section 4 specifies the computational and domain settings used for the simulations. The simulation results and their drag coefficient distribution are presented and discussed in Section 5 and Section 6, respectively. Finally, a conclusion is made, and future work is briefly discussed at the end of this paper. The results are expected to serve as a stepping stone in preparing LACs for hydrodynamic design.

2. Modeling Incompressible Fluid Flow and the Drag Coefficient

To solve an incompressible two-dimensional fluid flow problem using a numerical scheme, we first rewrite the Navier–Stokes Equation (NSE) in its weak formulation form [13]:
u t + u · u ν Δ u + 1 ρ p = f on Ω × ( 0 , T ) , · u = 0 on Ω × ( 0 , T ) , u ( x , 0 ) = u 0 in Ω , u ( x , t ) = g on Γ × ( 0 , T ) ,
where u is the velocity vector field, ν is the kinematic viscosity, p is the pressure, f represents any external force (which is usually zero), u 0 is the initial velocity, T is the maximum time (t), and g is the condition (velocity field) enforced on the boundary Γ of domain Ω . The additional f in the equation will not affect any of the subsequent analyses [15]. We set ρ = 1 and assume the effect of gravity is negligible for the problem.
We choose the Finite Element Method (FEM) for the space approximation of NSE. The semi-discrete form for Equation (1) is given as follows:
Ω d d t u h ^ ( t ) · v h ^ d x + Ω ν u h ^ ( t ) : v h ^ d x + Ω u h ^ ( t ) · u h ^ ( t ) · v h ^ d x Ω p h ^ ( t ) · v h ^ d x = Ω f h ^ ( t ) · v h ^ d x , v h ^ V h ^ ,
Ω q h ^ · u h ^ t d x = 0 , q h ^ M h ^ ,
u h ^ 0 = u h ^ 0 ,
where u h ^ 0 is an approximation of u 0 [13]. Our goal is to find u h ^ ( t ) V h ^ and p h ^ ( t ) M h ^ for all t ( 0 , T ) . We follow the standard notation for Sobolev spaces for the rest of this paper [16]. Using the Taylor–Hood discretization of space, the following spaces are obtained [13]:
V h ^ = X h ^ H 0 1 Ω 2 , X h ^ = v h ^ C 0 Ω ¯ 2 | v h ^ | K P 2 , K T h ^ , M h ^ = Y h ^ L 0 2 Ω , Y h ^ = q h ^ C 0 Ω ¯ | q h ^ | K P 1 , K T h ^ ,
where q h ^ is a continuous function, p h ^ , q h ^ = Ω p h ^ q h ^ d Ω , and q h ^ 0 = q h ^ , q h ^ 1 2 . The notation H 0 1 denotes the subset of the Sobolev space H 1 (of functions with at least one weak derivative), whose members are equal to zero at the boundaries of the domain. Meanwhile, C 0 denotes the space of continuous functions defined in Ω .
The Taylor–Hood elements satisfy the inf-sup condition, also known as the Ladyzhenskaya–Babuška–Brezzi condition (LBB), indicating that the system is well-posed [17,18]. The Taylor–Hood elements consist of continuous piece-wise quadratic polynomials for velocity and linear polynomials for pressure (see Figure 1). The notation T h ^ represents a regular triangulation of Ω , while P k denotes the Lagrange polynomials space of k-degree on the triangles K. The pressure, p, is taken into L 0 2 ( Ω ) for both Stokes equation and NSE. The notation L 0 2 ( Ω ) stands for the space of (generalized) functions L 2 ( Ω ), which are square-integrable and have a zero average on Ω [13,19]. This is because p is only unique up to a constant when the velocity’s Dirichlet boundary condition is imposed on all Ω . To counter this, we penalize the pressure to restrain the constant such that a non-singular algebraic system can be acquired [13]. The solution ( u ε , p ε ), where ε > 0 is the penalized parameter, produced by the penalized system, converges to the non-penalized system’s ( u , p ) [20].
The application of the second-order semi-implicit time-stepping method onto the NSE can produce the following problem [13]: given values { u j } j = 0 1 and a constant time step, Δ t , find the solution ( u n + 2 , p n + 2 ) of
1 Δ t α u n + 2 + j = i k β j u n + 2 j The discretization of time derivative ν Δ u n + 2 + j = 1 k γ j u n + 2 j · u n + 2 j Nonlinear advection extrapolated at t   =   t n + 2 + p n + 2 = f n + 2 .
The coefficients α , { β j } j = 1 2 R , where R represents real numbers, are the result of applying the backward differentiation formula on the time derivative with α + j = 1 2 β j = 0 . Meanwhile, the coefficients { γ j } j = 1 2 R generate the extrapolation formula for the nonlinear term, which satisfies j = 1 2 γ j = 0 . The terms p n + 2 and · u n + 2 used in the velocity–pressure coupling are taken implicitly to strictly enforce the discrete incompressibility condition. The diffusion term u n + 2 is taken implicitly as well, to prevent stringent stability conditions in O ( ν 1 h ^ 2 ) on the time step.
This study employs the second-order semi-implicit backward difference formula (SBDF) [13]. The SBDF does not self-start and hence requires a proper initialization for { u j } j = 1 2 , in other words, obtainin u 1 from u 0 , which can be achieved by simply employing the first-order SBDF. The time-discretized momentum equation when second-order SBDF is applied is written as [13]: given the initial solution u h 0 and a proper initialization, u h ^ 1 , we can find the solution ( u h ^ n + 2 , p h ^ n + 2 ) of
3 u h ^ n + 2 4 u h ^ n + 1 + u h ^ n 2 Δ t ν Δ u h ^ n + 2 + p h ^ n + 2 = f h ^ n + 2 2 B u h ^ n + 1 + B u h ^ n ,
for all n N , where the nonlinear advection term is denoted by B ( u h ^ ) = u h ^ · u h ^ . The notation N denotes the set of natural numbers.

3. Creating Same-Sized LAC Hydrofoils with Various Shapes

The hydrofoils are generated using the LAC Equation [21]:
P ( θ ) = 0 θ e ( Λ + i ) ψ d ψ , if α = 1 0 θ ( ( α 1 ) Λ ψ + 1 ) 1 α 1 e ψ i d ψ , otherwise
where α dictates the type and shape of the LAC generated, Λ is the rate of change of curvature with respect to the curve’s arc length, and θ is the turning angle of the curve. It is well known that LACs with α = 1 , 0, and 1 are the Euler spiral, logarithmic spiral, and Nielsen’s spiral, respectively. LACs with α < 0 are curves classified as divergent, α = 0 as neutral, and α > 0 as convergent. These classifications are based on the designers’ impression of the curve [22].
The hydrofoils used in this study were generated using Equation (8) by fixing the length (c) and height of the LAC, which makes up the upper half of the hydrofoil, through manipulating Λ and scaling the LAC while ensuring the tangent vector ( T ) at the leading edge is orthogonal to the x-axis (see Figure 2). The generated curve is then reflected along the x-axis to produce the complete hydrofoil shape. The generated shape is G 2 -continuous everywhere except at the trailing edge [23]. The thickness of the thickest section of the hydrofoil is denoted by h.
The steps for generating the hydrofoils are elucidated in Algorithm 1. Two examples of the output are shown in Figure 3 and Figure 4. The size difference in the figure refers to the difference between the size of the LAC of the input α value (orange) and the size of the LAC of α = 0 (grey).
Algorithm 1 Building LAC Hydrofoil with User-specified Size
Mathematics 11 00103 i001

4. Simulation and Domain Settings

The simulation domain comprises a rectangular boundary and the boundary of the hydrofoil-shaped object in it. The rectangular boundary has vertices ( 12.5 , 12.5 ) , ( 12.5 , 12.5 ) , ( 21 , 12.5 ) , and ( 21 , 12.5 ) . This is 33.5 times the chord length, c, of the LAC hydrofoil in length and 25 times c in height to ensure a sufficiently accurate simulation and drag coefficient computation [24,25]. The leading edge of the hydrofoil is set at ( 1.05551 , 0 ) , and the trailing edge at ( 0 , 0 ) .
The boundary conditions for the rectangular boundary’s edges ( L 1 , L 2 , L 3 , L 4 ) and the hydrofoil ( L 5 ) are shown in Figure 5. Each meshed domain has approximately 75,000 nodes. The maximum mesh spacing is 0.6 units on the rectangular boundary. The minimum mesh spacing is approximately 0.0004 units on the LAC hydrofoil. The meshed domain is illustrated in Figure 6. FreeFem++ [26] was employed to solve the incompressible Navier–Stokes equations. The space approximation was carried out using the Taylor–Hood finite element ( P 2 P 1 ) while a second-order semi-implicit backward difference formula was chosen for time integration. The combined method provides fairly accurate approximations for both velocity and pressure [13]. The drag coefficients are then recorded upon reaching a steady state or when the first derivative for velocity with respect to time (in the Navier–Stokes equation) is less than 10 5 .
Three experiments (cases) were designed in which the hydrofoils had different thicknesses at the thickest section of the hydrofoil-shaped object. The objects’ thickness h for Case 1, Case 2, and Case 3 are 0.234185, 0.3, and 0.4, respectively. In each case, the flow is simulated around hydrofoils with the same maximum thickness but different α values to find their drag coefficients C D . These hydrofoils have different shapes and leading-edge curvatures (see Figure 7, Figure 8 and Figure 9). The α values of the hydrofoils in Figure 7 (Case 1) are {−0.05, −0.025, −0, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.4, 0.5, 0.6, 0.75, 0.8, 0.9, 1}. The hydrofoil with the smallest α value is colored in blue, and the color gradually changes into yellow as α increases. The α values of the hydrofoil in Figure 8 (Case 2) are {−0.05, 0, 0.05, 0.1, 0.15, 0.25, 0.3, 0.4, 0.5, 0.6, 0.75, 0.8, 0.9, 1, 1.2} whereas for Figure 9 (Case 3), they are {−0.25, −0.05, 0, 0.05, 0.1, 0.15, 0.25, 0.3, 0.4, 0.5, 0.6, 0.75, 0.8, 0.9, 1, 1.5}. Note that the LACs with α = 0 and α = 1 are also known as the Nielsen’s spiral and Logarithmic spiral, respectively [21]. We need to fulfill the G 2 data at the leading edge, the G 0 data at the trailing edge, and the thickness of the hydrofoil to generate the hydrofoil shapes. However, due to the lack of the single-segment LACs’ degree of freedom, we can only create profiles within a specific range of α in each case, thus resulting in the difference between the range of α values used in the three cases.
Three simulations were run using different Reynolds numbers ( Re ): 30, 40, and 100. While Re = 30 was chosen arbitrarily, the other two Re values were chosen to examine how the LAC hydrofoils’ drag coefficients change when there is a small increment in Re , i.e., from 30 to 40, and a large increment from 40 to 100. It is also notable that at Re = 100 , the fluid’s viscosity is similar to that of water at 20 degrees Celsius, assuming 1 unit equals 0.1 m. The flows are assumed to be laminar, meaning the streamlines were smooth and regular. The Reynolds number is a similarity parameter that measures the ratio of inertial forces to viscous forces in a flow [27]. The fluid at Re = 100 is deemed to produce less viscous flows than Re = 30 and 40. We ran nine sets of simulations comprising 141 individual simulations in total. The results are presented and discussed in the following section.

5. Results and Discussion

The drag coefficients denoted as C D for each simulation are listed in Table 1, Table 2 and Table 3 and plotted in Figure 10, Figure 11 and Figure 12. The difference between the lowest C D and the highest is shown on the right side of Figure 10, Figure 11 and Figure 12.
It is observed that for each simulation of a different Re , the C D values for shapes of different α exhibit a similar trend, as shown in Case 1 and Case 2. This is reflected in the sudden drop in C D values for shapes of α = 0.8 . The C D values for these cases peaked around α = 0.6 to α = 0.75 . However, for Case 3, the peak of the C D graph with Re = 100 shifted left, from α = 0.8 to α = 0.6 , obviously deviating from the simulations with Re = 30 and Re = 40 . The standardized C D graphs are plotted in Figure 13. The standardization shifts and rescaled C D values of LAC hydrofoils with the same thickness but different Re have zero mean and unit sample variance for better comparison of the C D trends. The α values of the LAC shape that creates the most drag are shown in Table 4, along with the corresponding C D values. The rate of change of the drag coefficient decreases as Re increases. This statement is true for Cases 1–3.
The streamlines (flow lines), which show the local directions of the vector field of the simulations with the lowest and highest C D values at a steady state are shown in Figure 14, Figure 15 and Figure 16. The rainbow color in the background shows the scalar field, i.e., flow speed, with red indicating low speed and purple indicating higher speed. The detachment of a boundary layer from a surface is known as flow separation [28]. Flow separation is most apparent in the simulation of Case 3 with Re = 100 . More separation is seen in the shape, which creates a more substantial drag. As such, most of the drag force is made up of skin friction drag. The possible onset of wakes or flow separation in the simulation of Case 3 ( Re = 100 ) could be the reason for the change in the C D trend observed earlier.
Figure 17 presents the curvature profiles of each LAC shape and its corresponding case with different colors. The darker blue curves represent the curvature profile κ ( θ ) of the LAC shapes with lower α values, where θ is the turning angle of the curve. As α increases, the color of the κ ( θ ) curve turns yellow. The purple-, red-, and orange-colored curves are the κ ( θ ) curves for the LAC shape with the highest C D . Shapes with the lowest C D have an almost linear curvature κ ( θ ) . However, C D also decreases as the κ ( θ ) of the LAC shape bends more.
Thus, it is clear that C D values for hydrofoils constructed with neutral (Nielsen’s Spiral α = 0 ) and divergent LACs ( α < 0 ) are lower than those of convergent LACs ( α > 0 ) and decrease as α decreases. The C D for hydrofoil constructed with convergent LACs ( α > 0 ) does not increase or decrease monotonically as α increases or decreases. Instead, it gradually increases until it peaks at around α = 0.6 to 0.8 and decreases as α increases.
Figure 18, Figure 19 and Figure 20 illustrates the distribution of time taken to reach steady state for each α .

6. Cluster Analysis of Drag Distribution

In this section, we clustered drag coefficients in Table 1, Table 2 and Table 3 using agglomerative hierarchical cluster analysis [29]. This step groups a similar drag distribution among the three different thicknesses of LAC hydrofoils with different Reynolds numbers. In this analysis, we employed the complete linkage approach [30], which merges two clusters with the closest maximum distance:
d AB = max i A , j B d i j
where d i j is the distance and i and j are observations in clusters A and B . The complete linkage method is coupled with the dissimilarity distance matrix obtained from Dynamic Time Warping (DTW) [31,32]. DTW is an algorithm that measures the similarity or distance between two arrays or time series of different lengths [32]. The difference between DTW and Euclidean distance is elucidated in Figure 21. Two data vectors were connected based on their minimal distance using DTW (green) and Euclidean (gray) distance [33]. The DTW method provides more accurate results than the Euclidean distance. DTW does not require the data sets to be equal in length and is not affected by shifting, unlike Euclidean distance. The detailed algorithm of DTW can be found in [32].
Since we have computed drag coefficients with different h values from different α values, the DTW algorithm matches it in a direction that minimizes the distance of drag coefficients between various α values without losing information. Figure 22 shows a dendrogram illustrating the clusters obtained. The x-axis is labeled based on the cases, and their corresponding Reynolds number, e.g., C1_30, represents Case 1 with LAC object thickness h 0.2 and Re = 30 .
The dendrogram has two distinct clusters at the phenon line of distance 2.1 onward with a Cophenetic correlation coefficient (CCC) value equal to 0.86. The CCC measuring close to 1 indicates that the accuracy of the resulting dendrogram preserves the pair-wise distances between the drag distributions. The first cluster in orange consists of Case 3 with Re = 100 . We observed this cluster with the drag that may demonstrate the onset of unsteady flows at Re = 100 regardless of the LAC shapes. The second cluster (green) is the combination of cases with Re between 30 and 40, where the steady flow is guaranteed. Indeed, the two subclusters in green are also grouped based on the Reynolds number. The exception is C3_40, where the LAC hydrofoil has h = 0.4 and Re = 40 . The drag coefficients of C3_40 have a similar distribution to C1_30, where the LAC hydrofoil has h 0.2 and Re = 30 . Similarly, the case of C3_30 tends to move away from its own group of Re = 30 . The formation of these two subclusters indicates that the distribution of drag coefficients does not rely on the Reynolds number alone; the thickness of the LAC objects also plays a pivotal role in the distribution of drag coefficients. This shows that the evaluation of Reynolds numbers based on chord length (or characteristic length) for streamlined-shape objects may somehow need to be reformulated and is still an open problem. This is true since the chord length of Case 3 hydrofoils, c 1 , started to become much closer to its thickness of 0.4. We recall that characteristic length directly influences the Reynolds numbers besides the mean velocity and the eddy viscosity. In other words, the “more accurate” Reynolds number for C3_30 can be very close to that of C1_40. Hence, they form a subcluster. The initial effort of cluster analysis can be helpful to study the resulting flow behavior (in this case, by only accessing their drag coefficients), even if the Reynolds numbers are not properly evaluated. Further studies need to be carried out to make the DTW algorithm more practical for flow characterizations, especially in the post-processing of CFD results.

7. Conclusions

Simulations of incompressible fluid flow around streamlined shapes built using LAC, with chord length c 1 and thickness 0.234185, 0.3, and 0.4, were carried out. The results indicated that LAC shapes with negative α values, classified as divergent curves, and with almost linear curvature profiles κ ( θ ) , representing Clothoids, generate the lowest drag. As α increases, the drag coefficient increases until it reaches a maximum of around α = 0.6 to 0.8 and then decreases. Flow separation was not observed for any of the three thickness variations except for the thickest LAC shape (Case 3, h = 0.4 ) when Re = 100 . The separation may have caused the difference in C D for Case 3 when Re = 100 . It was also observed that the time used for the simulation to reach a steady state for Case 3 exhibits an entirely different trend compared to the other cases. Furthermore, the thickness of LAC objects can be varied to obtain similar drag coefficients for different Reynolds numbers. Thus, LAC with negative α values is better suited for designing submerged structures or objects that minimize drag.
The findings obtained in this study may help in decision-making in designing ship hulls or submerged bodies such as hydrofoils, underwater vehicles and structures, and marine-bio-logging tags [34]. A designer may now opt to prioritize LAC with negative values in designing such objects to minimize drag. They can also anticipate how the drag coefficient varies as the shape of the LAC changes. This shall act as a step towards implementing LACs in submerged objects or ship-hull designs. Additionally, this paper presents a new algorithm for generating the LAC shapes required to construct a LAC hydrofoil. The resulting LAC satisfies the user-specified height, length, G 1 Hermite data at one end, and G 0 Hermite data at the other.
For future research, we wish to extend our study to incompressible turbulent flow situations. We also hope to simulate incompressible flows around three-dimensional objects built with LAC or LA surfaces and examine their hydrodynamics and vector-field topology using an emerging method called Topological Data Analysis [35].

Author Contributions

Conceptualization, M.S.W., R.G., K.C.L. and K.T.M.; methodology, M.S.W. and R.G.; software, M.S.W., K.C.L. and R.G.; validation, M.S.W., R.G., K.C.L. and F.N.H.; formal analysis, M.S.W.; investigation, M.S.W.; resources, M.S.W., R.G. and K.T.M.; data curation, M.S.W.; writing—original draft preparation, M.S.W. and R.G.; writing—review and editing, R.G., K.C.L. and F.N.H.; visualization, M.S.W. and R.G.; supervision, R.G. and K.C.L.; project administration, M.S.W. and R.G.; funding acquisition, R.G. and K.T.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was initially funded by the Ministry of Education Malaysia under the Fundamental Research Grant Scheme (FRGS/1/2016/STG06/UMT/02/3). This work was supported by JST CREST Grant Number JPMJCR1911. It was also supported JSPS Grant-in-Aid for Scientific Research (B) Grant Number 19H02048, JSPS Grant-in-Aid for Challenging Exploratory Research Grant Number 26630038.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors acknowledge University Malaysia Terengganu and Japan Society for the Promotion of Science for providing the research fund, software, and facilities that were utilized for this research. Special thanks to the anonymous reviewers for the helpful feedback that improved the readability of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CADComputer-Aided Design
DTWDynamic Time Warping
FEMFinite Element Method
LACLog-Aesthetic Curve
LCGLogarithmic Curvature Graph
LDDCLogarithmic Distribution Diagram of Curvature
SBDFSemi-Implicit Backward Difference Formula

References

  1. Yoshimoto, F.; Harada, T. Analysis of the characteristics of curves in natural and factory products. In Proceedings of the 2nd IASTED International Conference on Vizualization, Imaging and Image Processing, Malaga, Spain, 9–12 September 2002; pp. 276–281. [Google Scholar]
  2. Harada, T.; Yoshimoto, F.; Moriyama, M. Aesthetic curve in the field of industrial design. In Proceedings of the IEEE Symposium on Visual Languages, Tokyo, Japan, 13–16 September 1999; pp. 38–47. [Google Scholar] [CrossRef]
  3. Miura, K.T. A general equation of aesthetic curves and its self-affinity. Comput.-Aided Des. Appl. 2006, 3, 457–464. [Google Scholar] [CrossRef]
  4. Gobithaasan, R.U.; Miura, K.T. Logarithmic curvature graph as a shape interrogation tool. Appl. Math. Sci. 2014, 8, 755–765. [Google Scholar] [CrossRef] [Green Version]
  5. Miura, K.T.; Shibuya, D.; Gobithaasan, R.U.; Usuki, S. Designing log-aesthetic splines with G2 continuity. Comput.-Aided Des. Appl. 2013, 10, 1021–1032. [Google Scholar] [CrossRef] [Green Version]
  6. Arslan, A.; Tari, E.; Ziatdinov, R.; Nabiyev, R.I. Transition curve modeling with kinematical properties: Research on Log-aesthetic curves. CAD Solut. LLC 2014, 11, 509–517. [Google Scholar] [CrossRef] [Green Version]
  7. Gobithaasan, R.U.; Yip, S.W.; Miura, K.T.; Madhavan, S. Optimal path smoothing with Log-aesthetic curves based on shortest distance, minimum bending energy or curvature variation energy. Comput.-Aided Des. Appl. 2020, 17, 639–658. [Google Scholar] [CrossRef]
  8. Kineri, Y.; Endo, S.; Maekawa, T. Surface design based on direct curvature editing. CAD Comput. Aided Des. 2014, 55, 1–12. [Google Scholar] [CrossRef]
  9. Suzuki, T. Application of Log-aesthetic curves to the eaves of a wooden house. In Proceedings of the 4th International Conference on Archi-Cultural Interactions through the Silk Road, Nishinomiya, Japan, 16–18 July 2016; pp. 67–70. [Google Scholar]
  10. Levien, R.; Séquin, C.H. Interpolating Splines: Which is the fairest of them all? Comput.-Aided Des. Appl. 2009, 6, 91–102. [Google Scholar] [CrossRef] [Green Version]
  11. Imai, T.; Shibutani, T.; Matsui, K.; Kumagai, S.; Tran, D.T.; Mu, K.; Maekawa, T. Curvature sensitive analysis of axially compressed cylindrical tubes with corrugated surface using isogeometric analysis and experiment. Comput. Aided Geom. Des. 2016, 49, 17–30. [Google Scholar] [CrossRef]
  12. John, V. Reference values for drag and lift of a two-dimensional time-dependent flow around a cylinder. Int. J. Numer. Methods Fluids 2004, 44, 777–788. [Google Scholar] [CrossRef]
  13. Loy, K.C.; Bourgault, Y. On efficient high-order semi-implicit time-stepping schemes for unsteady incompressible Navier–Stokes equations. Comput. Fluids 2017, 148, 166–184. [Google Scholar] [CrossRef]
  14. Wo, M.S.; Gobithaasan, R.U.; Miura, K.T.; Loy, K.C.; Yasmeen, S.; Harun, F.N. Log-aesthetic curves and their relation to fluid flow patterns in terms of streamlines. J. Comput. Des. Eng. 2020, 8, 55–68. [Google Scholar] [CrossRef]
  15. Quartapelle, L. Numerical Solution of the Incompressible Navier–Stokes Equations; Springer Science & Business Media: Berlin, Germany, 1993. [Google Scholar] [CrossRef]
  16. Girault, V.; Raviart, P.A. Mathematical Foundation of the Stokes Problem BT—Finite Element Methods for Navier–Stokes Equations: Theory and Algorithms; Springer: Berlin/Heidelberg, Germany, 1986; pp. 1–111. [Google Scholar] [CrossRef]
  17. Babuška, I. The finite element method with Lagrangian multipliers. Numer. Math. 1973, 20, 179–192. [Google Scholar] [CrossRef]
  18. Girault, V.; Raviart, P.A. Numerical Solution of the Stokes Problem in the Primitive Variables BT—Finite Element Methods for Navier–Stokes Equations: Theory and Algorithms; Springer: Berlin/Heidelberg, Germany, 1986; pp. 112–192. [Google Scholar] [CrossRef]
  19. Barrenechea, G.R.; Wachtel, A. The inf-sup stability of the lowest order Taylor–Hood pair on affine anisotropic meshes. IMA J. Numer. Anal. 2020, 40, 2377–2398. [Google Scholar] [CrossRef] [Green Version]
  20. Temam, R. Navier–Stokes Equations: Theory and Numerical Analysis; American Mathematical Society: Chelsea, MA, USA, 2000. [Google Scholar]
  21. Yoshida, N.; Saito, T. Interactive aesthetic curve segments. Vis. Comput. 2006, 22, 896–905. [Google Scholar] [CrossRef]
  22. Kanaya, I.; Nakano, Y.; Sato, K. Classification of Aesthetic Curves and Surfaces for Industrial Designs. Des. Discourse 2007, 2, 4. [Google Scholar]
  23. Gobithasan, R.; Ali, J. Towards G2 curve design with Timmer parametric cubic. In Proceedings of the Proceedings. International Conference on Computer Graphics, Imaging and Visualization, Penang, Malaysia, 2 July 2004; pp. 109–114. [Google Scholar] [CrossRef]
  24. Hardie, S. Drag Estimations on Experimental Aircraft Using CFD; Technical Report; Mälardalens Högskola: Västerås, Sweden, 2006. [Google Scholar]
  25. Khchine, Y.E.; Sriti, M. Boundary layer and amplified grid effects on aerodynamic performances of S809 airfoil for horizontal axis wind turbine (HAWT). J. Eng. Sci. Technol. 2017, 12, 3011–3022. [Google Scholar]
  26. Hecht, F. New development in freefem+. J. Numer. Math. 2012, 20, 251–265. [Google Scholar] [CrossRef] [Green Version]
  27. Anderson, J.D. Fundamentals of Aerodynamics, 6th ed.; McGraw-Hill Higher Education: New York, NY, USA, 2017. [Google Scholar]
  28. Anderson, J.D. Introduction to Flight, 5th ed.; McGraw-Hill Higher Education: New York, NY, USA, 2004. [Google Scholar]
  29. Everitt, B.S.; Landau, S.; Leese, M.; Stahl, D. Cluster Analysis, 5th ed.; John Wiley & Sons: Hoboken, NJ, USA, 2011. [Google Scholar]
  30. Everitt, B.S.; Dunn, G. Applied Multivariate Data Analysis, 2nd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2001. [Google Scholar] [CrossRef]
  31. Aghabozorgi, S.; Seyed Shirkhorshidi, A.; Teh, Y.W. Time-series clustering—A decade review. Inf. Syst. 2015, 53, 16–38. [Google Scholar] [CrossRef]
  32. Senin, P. Dynamic Time Warping Algorithm Review. Science 2008, 855, 40. [Google Scholar]
  33. Sobolewska, E. Dynamic Time Warping (DTW) as a Mean to Cluster Time Series. 2022. Available online: https://rstudio-pubs-static.s3.amazonaws.com/474160_0761428e8683401e941a92a4177254a4.html (accessed on 1 December 2022).
  34. Fiore, G.; Anderson, E.; Garborg, C.S.; Murray, M.; Johnson, M.; Moore, M.J.; Howle, L.; Shorter, K.A. From the track to the ocean: Using flow control to improve marine bio-logging tags for cetaceans. PLoS ONE 2017, 12, e0170962. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Bridel-Bertomeu, T.; Fovet, B.; Tierny, J.; Vivodtzev, F. Topological Analysis of High Velocity Turbulent Flow. In Proceedings of the 2019 IEEE 9th Symposium on Large Data Analysis and Visualization (LDAV), Vancouver, BC, Canada, 21 October 2019; pp. 87–88. [Google Scholar] [CrossRef]
Figure 1. Conventional representation of the Taylor–Hood P 2 P 1 element.
Figure 1. Conventional representation of the Taylor–Hood P 2 P 1 element.
Mathematics 11 00103 g001
Figure 2. An LAC that makes up half of the hydrofoil.
Figure 2. An LAC that makes up half of the hydrofoil.
Mathematics 11 00103 g002
Figure 3. Example output of Algorithm 1 (input α = 0.75 ).
Figure 3. Example output of Algorithm 1 (input α = 0.75 ).
Mathematics 11 00103 g003
Figure 4. Example output of Algorithm 1 (input α = 0.25 ).
Figure 4. Example output of Algorithm 1 (input α = 0.25 ).
Mathematics 11 00103 g004
Figure 5. FEM domain with the LAC hydrofoil ( α = 1 ) and boundary conditions.
Figure 5. FEM domain with the LAC hydrofoil ( α = 1 ) and boundary conditions.
Mathematics 11 00103 g005
Figure 6. Meshed FEM domain with the LAC hydrofoil ( α = 0.5 ).
Figure 6. Meshed FEM domain with the LAC hydrofoil ( α = 0.5 ).
Mathematics 11 00103 g006
Figure 7. Comparison of Case 1 ( h = 0.234185 ) hydrofoil shapes from α = 0.05 (blue) to α = 1 (yellow).
Figure 7. Comparison of Case 1 ( h = 0.234185 ) hydrofoil shapes from α = 0.05 (blue) to α = 1 (yellow).
Mathematics 11 00103 g007
Figure 8. Comparison of Case 2 ( h = 0.3 ) hydrofoil shapes from α = 0.05 (blue) to α = 1.2 (yellow).
Figure 8. Comparison of Case 2 ( h = 0.3 ) hydrofoil shapes from α = 0.05 (blue) to α = 1.2 (yellow).
Mathematics 11 00103 g008
Figure 9. Comparison of Case 3 ( h = 0.4 ) hydrofoil shapes from α = 0.25 (blue) to α = 1.5 (yellow).
Figure 9. Comparison of Case 3 ( h = 0.4 ) hydrofoil shapes from α = 0.25 (blue) to α = 1.5 (yellow).
Mathematics 11 00103 g009
Figure 10. The drag coefficient of objects for various α values ( Re = 30 ).
Figure 10. The drag coefficient of objects for various α values ( Re = 30 ).
Mathematics 11 00103 g010
Figure 11. Drag coefficient of objects for various α values ( Re = 40 ).
Figure 11. Drag coefficient of objects for various α values ( Re = 40 ).
Mathematics 11 00103 g011
Figure 12. C D of objects for various α values ( Re = 100 ).
Figure 12. C D of objects for various α values ( Re = 100 ).
Mathematics 11 00103 g012
Figure 13. Comparison of the C D graph shape of all simulations.
Figure 13. Comparison of the C D graph shape of all simulations.
Mathematics 11 00103 g013
Figure 14. Streamline plots ( Re = 30 ).
Figure 14. Streamline plots ( Re = 30 ).
Mathematics 11 00103 g014
Figure 15. Streamline plots ( Re = 40 ).
Figure 15. Streamline plots ( Re = 40 ).
Mathematics 11 00103 g015
Figure 16. Streamline plots ( Re = 100 ).
Figure 16. Streamline plots ( Re = 100 ).
Mathematics 11 00103 g016
Figure 17. Curvature profiles of the objects.
Figure 17. Curvature profiles of the objects.
Mathematics 11 00103 g017
Figure 18. Final time (T) for each simulation with shapes of h 0.2 (Case 1).
Figure 18. Final time (T) for each simulation with shapes of h 0.2 (Case 1).
Mathematics 11 00103 g018
Figure 19. Final time (T) for each simulation with shapes of h = 0.3 (Case 2).
Figure 19. Final time (T) for each simulation with shapes of h = 0.3 (Case 2).
Mathematics 11 00103 g019
Figure 20. Final time (T) for each simulation with shapes of h = 0.4 (Case 3).
Figure 20. Final time (T) for each simulation with shapes of h = 0.4 (Case 3).
Mathematics 11 00103 g020
Figure 21. Visual comparison of matched points based on DTW (green) and Euclidean (gray) distance.
Figure 21. Visual comparison of matched points based on DTW (green) and Euclidean (gray) distance.
Mathematics 11 00103 g021
Figure 22. Dendrogram with two distinct clusters at a phenon line of distance 2.3.
Figure 22. Dendrogram with two distinct clusters at a phenon line of distance 2.3.
Mathematics 11 00103 g022
Table 1. α and C D for simulations with Re = 30 .
Table 1. α and C D for simulations with Re = 30 .
α Case 1Case 2Case 3
−0.25--0.521258
−0.050.4466400.4753290.522634
−0.0250.446800--
00.4469690.4757340.522961
0.050.4473400.4761120.523232
0.10.4476970.4764870.523516
0.150.4480720.4768430.523750
0.20.448431--
0.250.4487610.4777210.524190
0.30.4490710.478150.524372
0.40.4495560.4784540.524683
0.50.4498600.4786410.524931
0.60.4500000.478690.525093
0.750.4499510.478660.525211
0.80.4497830.4785190.525236
0.90.4495550.4783270.525195
10.4491680.4776690.525139
1.2-0.475734-
1.5--0.524178
Table 2. α and C D for simulations with Re = 40 .
Table 2. α and C D for simulations with Re = 40 .
α Case 1Case 2Case 3
−0.25--0.448448
−0.050.379420.4058160.449585
−0.0250.379562--
00.3797120.4061680.449856
0.050.3800410.4064950.450076
0.10.3803580.4068210.450311
0.150.380690.4071270.450501
0.20.381008--
0.250.3812990.4078850.450860
0.30.3815750.4082550.451006
0.40.3820050.4085150.451255
0.50.3822760.4086740.451452
0.60.3824030.4087140.451575
0.750.3823720.4086880.451655
0.80.3822290.4085660.451671
0.90.382040.4084030.451623
10.3817120.4078440.451566
1.2-0.407665-
1.5--0.450725
Table 3. α and C D for simulations with Re = 100 .
Table 3. α and C D for simulations with Re = 100 .
α Case 1Case 2Case 3
−0.25--0.285798
−0.050.2300450.2508440.286243
−0.0250.230138--
00.2302360.2510510.286353
0.050.2304540.2512420.286429
0.10.2306600.2514310.286522
0.150.2308790.2516070.286585
0.20.231089--
0.250.2312780.2519170.286711
0.30.2314610.2520420.286755
0.40.2317480.2522500.286825
0.50.2319330.2523940.286874
0.60.2320260.2524780.286883
0.750.2320350.2524870.286845
0.80.2319500.2524690.286833
0.90.2318580.2523920.286756
10.2316780.2523010.286686
1.2-0.251984-
1.5--0.286111
Table 4. α values of the LAC shape that create the highest C D and the corresponding C D .
Table 4. α values of the LAC shape that create the highest C D and the corresponding C D .
Case 1 ( h 0.2 )Case 2 ( h = 0.3 )Case 3 ( h = 0.4 )
α C D α C D α C D
Re = 30 0.60.4500000.750.4786600.80.525236
Re = 40 0.60.3824030.750.4086880.80.451671
Re = 100 0.750.2320350.750.2524870.60.286883
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wo, M.S.; Gobithaasan, R.U.; Miura, K.T.; Loy, K.C.; Harun, F.N. Analysis of Drag Coefficients around Objects Created Using Log-Aesthetic Curves. Mathematics 2023, 11, 103. https://doi.org/10.3390/math11010103

AMA Style

Wo MS, Gobithaasan RU, Miura KT, Loy KC, Harun FN. Analysis of Drag Coefficients around Objects Created Using Log-Aesthetic Curves. Mathematics. 2023; 11(1):103. https://doi.org/10.3390/math11010103

Chicago/Turabian Style

Wo, Mei Seen, R.U. Gobithaasan, Kenjiro T. Miura, Kak Choon Loy, and Fatimah Noor Harun. 2023. "Analysis of Drag Coefficients around Objects Created Using Log-Aesthetic Curves" Mathematics 11, no. 1: 103. https://doi.org/10.3390/math11010103

APA Style

Wo, M. S., Gobithaasan, R. U., Miura, K. T., Loy, K. C., & Harun, F. N. (2023). Analysis of Drag Coefficients around Objects Created Using Log-Aesthetic Curves. Mathematics, 11(1), 103. https://doi.org/10.3390/math11010103

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop