The idea of analytical modeling of complex geometric
objects by implicit functions (as an inverse problem of analytical geometry)
has existed for a long time [1-5] and has a number of significant advantages
over surface models described by a parametric form for obtaining the
coordinates of boundary points.
First of all, in implicit modeling a geometric object
is defined by a zero boundary range of values, which adds to the applicability
of such approach in engineering calculations. Secondly, it can be noted that there
is no limit for the specified dimensionality of space.
However, further application possibilities in computer
technologies require an analytical function given by a complicated expression to
be additionally simplified by discretization on a given domain. First of all,
computer information is discrete and requires a transition from functional
continuity to discrete continuity. Discrete continuity should provide properties
and identity retention of functional continuity, i.e. it should also be
represented by a local function of an implicit form at a specified point.
Research papers [6-10] reflect the basic principles of Functional
Voxel modeling method (FV-method), which is based on constructing a local
function with an angular metric at each point
(
)
on a given domain:
|
|
(1)
|
The computer storage of angular metrics, referred to in
the FV-method as “local geometric characteristics”, is carried out through converting
the colour palette into numerical values, which enables to obtain special images
of local geometric characteristics (M-images) on a given domain, as well as to reassemble
them by inverse transformation into cosines to obtain a local function at the considered
point of the M-image.
For example, M-images for a circle:
|
or
|
(2)
|
could be represented by a local function
.
M-images describing the region
are shown in Figure 1.
Fig 1. M-images of a circle in RGB colour format, which displays 16777215 shades of colour
Figure 2 shows semi-tone images
of local geometric characteristics for convenient visual assessment, but with a
limited accuracy of representation up to 255 shades of gray.
Fig 2. M-images of a circle. Monochromatic colour palette displaying 256 shades of gray colour
The main advantage of a local
function is that it can be quite simply represented as a sum of arguments
multiplied by local geometric characteristics, which makes it easy to express
the desired variable through the remaining arguments of a function, providing
simple calculations for a point selected on a given domain. For example:
|
|
(3)
|
Figure 3 illustrates
in monochrome the region representing values calculated for variable
. Figure 3. a) shows the normalized values
of variable
,
and Figure 3.b) highlights in blue color the negative values of
outside the object, and positive values
inside the object in white.
|
|
|
|
a)
|
b)
|
Fig. 3. Values of
distributed over a given region for the function (2): a)
values are
normalized from white to black b) region separation
Many computer applications
based on the Functional-Voxel modelling method (FV-method) have been
implemented to solve problems in various areas of mathematical modeling [3-7],
resulting in a visual solution that is unlimited in dimension.
At the same time, the proposed
FV-method is based on a discrete representation of data, referred to as
M-images, obtained by a linear approximation of three nearest neighbour vertices
in a regular mesh, which obviously leads to limitations in the accuracy of data
representation, as well as to a limited area for calculations.
In order to
overcome these limitations, we apply the differential principles of
constructing a local function.
Let us consider again an example of a local
function defined by an implicit equation and describing a two-dimensional
domain of the function on the
plane, and
representing its value:
|
|
(4)
|
Let's reduce equation (4) to
the differential form:
|
|
(5)
|
In essence, formula (5) leads
to the following expression:
|
|
(6)
|
which implies that
,
,
– analytically defined surfaces
of partial derivative functions, where
equals one.
Let's look at an example with
an analytical description of a circle, describing how a discrete local function
defined for each point of a given domain can be expanded to a general
analytical form.
We consider equation (2) as the
initial equation for describing a circle. Taking first order partial
derivatives with respect to three arguments, we obtain:
|
|
(7)
|
Based on the equality (6), the
fourth derivative can be expressed:
|
|
(8)
|
We obtain local geometric characteristics
by multiplying all four differential components (derivatives) by the third
component
|
|
(9)
|
and then, normalizing by the length of the
homogeneous normal vector, we get the cosine values for the components:
|
|
(10)
|
|
|
Figure 4 represents M-images
displaying local geometric characteristics distributed over the same region for
the equation of a circle as those shown in Figure 1, but obtained without
linear approximation [6-10].
If we
compare both images, we see some similarities but with the obvious difference
in color patterns. This is due to the increased accuracy of the obtained
representation compared to the mesh approximation.
Fig. 4. M-images of
a circle in RGB colour format, which displays 16777215 shades of colour
Note that Figure 5 is quite
comparable to Figure 2, since a black-and-white palette and only 256 tone
gradations were used to construct the image, i.e. a lower colour resolution.
Fig. 5. M-images of a circle equation in monochrome
format, displaying 256 shades of colour
Given the fact that most
analytical functions describing objects of analytical geometry can be
differentiable at all points of the domain and along each axis, we can say with
confidence that the proposed approach is fully consistent.
This representation
significantly expands the scope of application of the local function with all
its remarkable properties.
First, we can unequivocally
assert an increase in the accuracy of local geometric characteristics. Second,
here we have the advantage of the most compact description of local geometric
characteristics by analytical expressions, which eliminates the problem of
limitations arising from a predefined domain of arguments and the number of
gradations of the color palette, which expands its applicability in
multi-object modeling.
Of course, it can be argued
that, for example, the R-function is difficult to differentiate and the
transition to representing a geometry by a local function may, on the contrary,
further complicate calculations. Let's try to solve this problem using the
developed tools of local geometry, where R-function is adapted to the
calculation of local characteristics [10].
To begin with, using the
experience gained in modeling a local function for a circle, let's
R-functionally model a three-dimensional cylinder figure with a height of
as a local function in a three-dimensional
region:
|
|
(11)
|
To do this, we use the
intersection of the regions of functions describing two spatial figures
and
:
an infinite horizontal line with a height equal to
and an infinite vertical cylinder with the radius
.
|
and
|
(12)
|
To obtain an analytically
described implicit function for a cylinder, an R-intersection is usually used:
|
|
(13)
|
Using partial derivatives, the
figure of an infinite cylinder is described by analogy with the equation of a
circle, where there is no influence of the
-coordinate:
|
|
(14)
|
Figures 6 and 7 demonstrate M-images
for the region
.
Fig.6. M-images of an infinite cylinder in RGB
colour format, which displays 16777215 shades of colour
Fig.7. M-images of an infinite cylinder in monochrome
format, displaying 256 shades of colour
Restrictions along the
-axis
with the height of
:
|
|
(15)
|
Figures 8 and 9 demonstrate M-images
for this example, the region:
.
Fig.8. M-images of a
horizontal line function in RGB colour format, which displays 16777215 shades of colour
Fig.9. M-images of a horizontal line function
in monochrome format, displaying 256 shades of colour
Let us apply the formula given
in [8] for the R-intersection of differentials. For convenience, we will apply
the previously defined notations, shortening the differential equation, and
transmit it into local geometric characteristics:
M-images representing a plane section of a cylinder in the
-plane
in RGB color format for computer
storage are shown in Figure 10.
Fig.10. M-images of a cylindric section in RGB
colour format, which displays 16777215 shades of colour
For convenient visual
evaluation by a person, Figure 11 shows these M-images in monochrome format.
Fig.11. M-images of a cylindric section in monochrome
format, displaying 256 shades of colour
The presented principle of
discretization of continuous space by continuous local functions allows us to
formulate further problems related to the analytical description of the
geometry of a complex technogenic environment by implicit functions. The idea
is to identify a group of basic differentials that at the local level enables description
of the geometry of the function space of any dimension and complexity of the
formulated object, which allows us to think about expanding the tools of the
modern graphics kernel, leading them to a multidimensional representation for
the implementation of design and control tasks.
The research was carried out within the framework of the
scientific program of the National Center for Physics and Mathematics,
direction No. 9 "Artificial intelligence and big data in technical,
industrial, natural and social systems".
1. Kaledin V.O., Reshetnikova E.V., Gridchina V.B. Theory of R-functions: study guide for Higher Education Institutions in the field of Applied Mathematics and Computer Science. Kemerovo State University, Novokuznetsk Institute. 2017. - 119 pp.
2. Kolyada R.Ђ., Maximenko-Sheiko K.V., Sheiko ’.I. R-Functions Method in the Mathematical Modeling of Convective Heat Exchange in an Octahedral Fuel Assembly with 37 Fuel Elements // Journal of Mathematical Sciences 238, 2019. Р P.154 -164.
3. Uspensky B., Avramov K., Derevianko I., Maksymenko-Sheiko K. Vibrations of cylindrical sandwich shell with fused deposition processed honeycomb core and carbon nanotubes reinforced composite faces sheets // Journal of Vibration Engineering & Technologies 12 (2), 2024. Р P. 2003-2023.
4. Maximenko-Sheiko K.V., Sheiko T.I. Mathematical modeling of geometric fractals using R-functions // Cybernetics and Systems Analysis 48, 2012. Р P. 614-620.
5. Litvinova J.S., Maximenko-Sheiko K.V., Sheiko T.I. Mathematical and computer modeling of construction structures based on R-functions [In Russian]. Problems of mechanical engineering, Ь17 (3), 2014. Р P. 12-18.
6. Tolok, A.V. Functional-voxel method in computer modeling [In Russian]. Moscow: Fizmatlit. 2016. 112 p.
7. Tolok, A.V. and Tolok, N.B. Local computer geometry [In Russian]. Moscow: IPR Media, 2025. 228 p.
8. Tolok, A.V. and Tolok, N.B. Functional-Voxel Modeling of The Cauchy Problem // Scientific Visualization, 2024, volume 16, number 1, pages 105 - 111, DOI: 10.26583/sv.16.1.09.
9. Tolok, A.V. and Tolok, N.B., Modeling the Solution of an Ordinary Differential Equation by the Functional Voxel Method // Scientific Visualization, 2024, volume 16, number 3, pages 37 - 47, DOI: 10.26583/sv.16.3.04.
10. Tolok, A.V. and Tolok, N.B., Functional Voxel Modeling of a Path Planning Algorithm to a Target Based on R-Functions. Control Sciences 5, 41Р47 (2024).