Abstract: The Robin problem for the Helmholtz equation in normal-polar annuli is addressed by using a suitable Fourier-Hankel series technique. Attention is in particular focused on the wide class of domains whose boundaries are defined by the so-called superformula introduced by Gielis. A dedicated numerical procedure based on the computer algebra system Mathematica© is developed in order to validate the proposed methodology. In this way, highly accurate approximations of the solution, featuring properties similar to the classical ones, are obtained.

Abstract: In this work, the electroporation phenomenon induced by pulsed electric field on different nucleated biological cells is studied. A nonlinear, non-local, dispersive, and space-time multiphysics model based on Maxwell's and asymptotic Smoluchowski's equations has been developed to calculate the transmembrane voltage and pore density on both plasma and nuclear membrane perimeters. The irregular cell shape has been modeled by incorporating in the numerical algorithm the analytical functions pertaining to Gielis curves. The dielectric dispersion of the cell media has been modeled considering the multi-relaxation Debye-based relationship. Two different irregular nucleated cells have been investigated and their response has been studied applying both the dispersive and non-dispersive models. By a comparison of the obtained results, differences can be highlighted confirming the need to make use of the dispersive model to effectively investigate the cell response in terms of transmembrane voltages, pore densities, and electroporation opening angle, especially when irregular cell shapes and short electric pulses are considered. Bioelectromagnetics. 2019;40:331-342. (c) 2019 Wiley Periodicals, Inc.

Abstract: During the past decades, the poration of cell membrane induced by pulsed electric fields has been widely investigated. Since the basic mechanisms of this process have not yet been fully clarified, many research activities are focused on the development of suitable theoretical and numerical models. To this end, a nonlinear, nonlocal, dispersive, and space-time numerical algorithm has been developed and adopted to evaluate the transmembrane voltage and pore density along the perimeter of realistic irregularly shaped cells. The presented model is based on the Maxwell's equations and the asymptotic Smoluchowski's equation describing the pore dynamics. The dielectric dispersion of the media forming the cell has been modeled by using a general multirelaxation Debye-based formulation. The irregular shape of the cell is described by using the Gielis' superformula. Different test cases pertaining to red blood cells, muscular cells, cell in mitosis phase, and cancer-like cell have been investigated. For each type of cell, the influence of the relevant shape, the dielectric properties, and the external electric pulse characteristics on the electroporation process has been analyzed. The numerical results demonstrate that the proposed model is an efficient numerical tool to study the electroporation problem in arbitrary-shaped cells.

Abstract: The interior and exterior Robin problems for the Helmholtz equation in starlike planar domains are addressed by using a suitable Fourier-like technique. Attention is in particular focused on normal-polar domains whose boundaries are defined by the so-called superformula introduced by J. Gielis. A dedicated numerical procedure based on the computer algebra system Mathematica© is developed in order to validate the proposed approach. In this way, highly accurate approximations of the solution, featuring properties similar to the classical ones, are obtained. The computed results are found to be in good agreement with the theoretical findings on Fourier series expansion presented by L. Carleson.

Abstract: Plant leaves exhibit diverse shapes that enable them to utilize a light resource maximally. If there were a general parametric model that could be used to calculate leaf area for different leaf shapes, it would help to elucidate the adaptive evolutional link among plants with the same or similar leaf shapes. We propose a simplified version of the original Gielis equation (SGE), which was developed to describe a variety of object shapes ranging from a droplet to an arbitrary polygon. We used this equation to fit the leaf profiles of 53 species (among which, 48 bamboo plants, 5 woody plants, and 10 geographical populations of a woody plant), totaling 3310 leaves. A third parameter (namely, the floating ratio c in leaf length) was introduced to account for the case when the theoretical leaf length deviates from the observed leaf length. For most datasets, the estimates of c were greater than zero but less than 10%, indicating that the leaf length predicted by the SGE was usually smaller than the actual length. However, the predicted leaf areas approximated their actual values after considering the floating ratios in leaf length. For most datasets, the mean percent errors of leaf areas were lower than 6%, except for a pooled dataset with 42 bamboo species. For the elliptical, lanceolate, linear, obovate, and ovate shapes, although the SGE did not fit the leaf edge perfectly, after adjusting the parameter c, there were small deviations of the predicted leaf areas from the actual values. This illustrates that leaves with different shapes might have similar functional features for photosynthesis, since the leaf areas can be described by the same equation. The anisotropy expressed as a difference in leaf shape for some plants might be an adaptive response to enable them to adapt to different habitats.

Abstract: There is convincing evidence for a scaling relationship between leaf dry weight (DW) and leaf surface area (A) for broad-leaved plants, and most estimates of the scaling exponent of DW vs. A are greater than unity. However, the scaling relationship of leaf fresh weight (FW) vs. A has been largely neglected. In the present study, we examined whether there is a statistically strong scaling relationship between FW and A and compared the goodness of fit to that of DW vs. A. Between 250 and 520 leaves from each of 12 bamboo species within 2 genera (Phyllostachys and Pleioblastus) were investigated. The reduced major axis regression protocols were used to determine scaling relationships. The fit for the linearized scaling relationship of FW vs. A was compared with that of DW vs. A using the coefficient of determination (i.e., r2). A stronger scaling relationship between FW and A than that between DW and A was observed for each of the 12 bamboo species investigated. Among the 12 species examined, five had significantly smaller scaling exponents of FW vs. A compared to those of DW vs. A; only one species had a scaling exponent of FW vs. A greater than that of DW vs. A. No significant difference between the two scaling exponents was observed for the remaining 6 species. Researchers conducting future studies might be well advised to consider the influence of leaf fresh weight when exploring the scaling relationships of foliar biomass allocation patterns.

Abstract: According to Thompson’s principle of similarity, the area of an object should be proportional to its length squared. However, leaf area–length data of some plants have been demonstrated not to follow the principle of similarity. We explore the reasons why the leaf area–length allometry deviates from the principle of similarity and examine whether there is a general model describing the relationship among leaf area, width and length. We sampled more than 11,800 leaves from six classes of woody and herbaceous plants and tested the leaf area–length allometry. We compared six mathematical models based on root-mean-square error as the measure of goodness-of-fit. The best supported model described a proportional relationship between leaf area and the product of leaf width and length (i.e., the Montgomery model). We found that the extent to which the leaf area–length allometry deviates from the principle of similarity depends upon the extent of variation of the ratio of leaf width to length. Estimates of the parameter of the Montgomery model ranged between 1/2, which corresponds to a triangular leaf with leaf length as its height and leaf width as its base, and π/4, which corresponds to an elliptical leaf with leaf length as its major axis and leaf width as its minor axis, for the six classes of plants. The narrow range in practice of the Montgomery parameter implies an evolutionary stability for the leaf area of large-leaved plants despite the fact that leaf shapes of these plants are rather different.

Abstract: A novel class of supershaped dielectric lens antennas, whose geometry is described by the three-dimensional (3D) Gielis formula, is introduced and analysed. To this end, a hybrid modelling approach based on geometrical and physical optics is adopted in order to efficiently analyse the multiple wave reflections occurring within the lens and to evaluate the relevant impact on the radiation properties of the antenna under analysis. The developed modelling procedure has been validated by comparison with numerical results already reported in the literature and, afterwards, applied to the electromagnetic characterisation of Gielis dielectric lens antennas with shaped radiation pattern. Furthermore, a dedicated optimisation algorithm based on quantum particle swarm optimisation has been developed for the synthesis of 3D supershaped lens antennas with single feed, as well as with beamforming capabilities.

Abstract: Many cross-sectional shapes of plants have been found to approximate a superellipse rather than an ellipse. Square bamboos, belonging to the genus Chimonobambusa (Poaceae), are a group of plants with round-edged square-like culm cross sections. The initial application of superellipses to model these culm cross sections has focused on Chimonobambusa quadrangularis (Franceschi) Makino. However, there is a need for large scale empirical data to confirm this hypothesis. In this study, approximately 750 cross sections from 30 culms of C. utilis were scanned to obtain cross-sectional boundary coordinates. A superellipse exhibits a centrosymmetry, but in nature the cross sections of culms usually deviate from a standard circle, ellipse, or superellipse because of the influences of the environment and terrain, resulting in different bending and torsion forces during growth. Thus, more natural cross-sectional shapes appear to have the form of a deformed superellipse. The superellipse equation with a deformation parameter (SEDP) was used to fit boundary data. We find that the cross-sectional shapes (including outer and inner rings) of C. utilis can be well described by SEDP. The adjusted root-mean-square error of SEDP is smaller than that of the superellipse equation without a deformation parameter. A major finding is that the cross-sectional shapes can be divided into two types of superellipse curves: hyperellipses and hypoellipses, even for cross sections from the same culm. There are two proportional relationships between ring area and the product of ring length and width for both the outer and inner rings. The proportionality coefficients are significantly different, as a consequence of the two different superellipse types (i.e., hyperellipses and hypoellipses). The difference in the proportionality coefficients between hyperellipses and hypoellipses for outer rings is greater than that for inner rings. This work informs our understanding and quantifying of the longitudinal deformation of plant stems for future studies to assess the influences of the environment on stem development. This work is also informative for understanding the deviation of natural shapes from a strict rotational symmetry.

Abstract: The Lamé curve is an extension of an ellipse, the latter being a special case. Dr. Johan Gielis further extended the Lamé curve in the polar coordinate system by introducing additional parameters (n1, n2, n3; m): rφ=1Acosm4φn2+1Bsinm4φn3−1/n1, which can be applied to model natural geometries. Here, r is the polar radius corresponding to the polar angle φ; A, B, n1, n2 and n3 are parameters to be estimated; m is the positive real number that determines the number of angles of the Gielis curve. Most prior studies on the Gielis equation focused mainly on its applications. However, the Gielis equation can also generate a large number of shapes that are rotationally symmetric and axisymmetric when A = B and n2 = n3, interrelated with the parameter m, with the parameters n1 and n2 determining the shapes of the curves. In this paper, we prove the relationship between m and the rotational symmetry and axial symmetry of the Gielis curve from a theoretical point of view with the condition A = B, n2 = n3. We also set n1 and n2 to take negative real numbers rather than only taking positive real numbers, then classify the curves based on extremal properties of r(φ) at φ = 0, π/m when n1 and n2 are in different intervals, and analyze how n1, n2 precisely affect the shapes of Gielis curves.

Abstract: Electroporation is a non-thermal electromagnetic phenomenon widely used in medical diseases treatment. Different mathematical models of electroporation have been proposed in literature to study pore evolution in biological membranes. This paper presents a nonlinear dispersive multiphysic model of electroporation in irregular shaped biological cells in which the spatial and temporal evolution of the pores size is taken into account. The model solves Maxwell and asymptotic Smoluchowski equations and it describes the dielectric dispersion of cell media using a Debye-based relationship. Furthermore, the irregular cell shape has been modeled using the Gielis superformula. Taking into account the cell in mitosis phase, the electroporation process has been studied comparing the numerical results pertaining the model with variable pore radius with those in which the pore radius is supposed constant. The numerical analysis has been performed exposing the biological cell to a rectangular electric pulse having duration of 10 μs. The obtained numerical results highlight considerable differences between the two different models underling the need to include into the numerical algorithm the differential equation modeling the spatial and time evolution of the pores size.