Loading [MathJax]/jax/output/SVG/jax.js
Research article

Phase Field Theory and Analysis of Pressure-Shear Induced Amorphization and Failure in Boron Carbide Ceramic

  • Received: 10 April 2014 Accepted: 11 June 2014 Published: 08 July 2014
  • A nonlinear continuum phase field theory is developed to describe amorphization of crystalline elastic solids under shear and/or pressure loading. An order parameter describes the local degree of crystallinity. Elastic coefficients can depend on the order parameter, inelastic volume change may accompany the transition from crystal to amorphous phase, and transitional regions parallel to bands of amorphous material are penalized by interfacial surface energy. Analytical and simple numerical solutions are obtained for an idealized isotropic version of the general theory, for an element of material subjected to compressive and/or shear loading. Solutions compare favorably with experimental evidence and atomic simulations of amorphization in boron carbide, demonstrating the tendency for structural collapse and strength loss with increasing shear deformation and superposed pressure.

    Citation: D. Clayton John. Phase Field Theory and Analysis of Pressure-Shear Induced Amorphization and Failure in Boron Carbide Ceramic[J]. AIMS Materials Science, 2014, 1(3): 143-158. doi: 10.3934/matersci.2014.3.143

    Related Papers:

    [1] Shouyun Cheng, Lin Wei, Xianhui Zhao, Yinbin Huang, Douglas Raynie, Changling Qiu, John Kiratu, Yong Yu . Directly catalytic upgrading bio-oil vapor produced by prairie cordgrass pyrolysis over Ni/HZSM-5 using a two stage reactor. AIMS Energy, 2015, 3(2): 227-240. doi: 10.3934/energy.2015.2.227
    [2] Wasinton Simanjuntak, Kamisah Delilawati Pandiangan, Tika Dwi Febriyanti, Aryani Putri Islami, Sutopo Hadi, Ilim Ilim . Catalytic upgrading of palm oil derived bio-crude oil for bio-hydrocarbon enrichment using protonated zeolite-Y as catalyst. AIMS Energy, 2024, 12(3): 600-616. doi: 10.3934/energy.2024028
    [3] Vincent E. Efeovbokhan, Augustine O. Ayeni, Osuvwe P. Eduvie, James A. Omoleye, Oladotun P. Bolade, Ajibola T. Ogunbiyi, Victoria N. Anyakora . Classification and characterization of bio-oil obtained from catalytic and non-catalytic pyrolysis of desludging sewage sample. AIMS Energy, 2020, 8(6): 1088-1107. doi: 10.3934/energy.2020.6.1088
    [4] Xianhui Zhao, Lin Wei, James Julson . First stage of bio-jet fuel production: non-food sunflower oil extraction using cold press method. AIMS Energy, 2014, 2(2): 193-209. doi: 10.3934/energy.2014.2.193
    [5] Ahmad Indra Siswantara, Illa Rizianiza, Tanwir Ahmad Farhan, M. Hilman Gumelar Syafei, Dyas Prawara Mahdi, Candra Damis Widiawaty, Adi Syuriadi . Analyzing temperature distribution in pyrolysis systems using an atomic model. AIMS Energy, 2023, 11(6): 1012-1030. doi: 10.3934/energy.2023048
    [6] Lihao Chen, Hu Wu, Kunio Yoshikawa . Research on upgrading of pyrolysis oil from Japanese cedar by blending with biodiesel. AIMS Energy, 2015, 3(4): 869-883. doi: 10.3934/energy.2015.4.869
    [7] Lifita N. Tande, Valerie Dupont . Autothermal reforming of palm empty fruit bunch bio-oil: thermodynamic modelling. AIMS Energy, 2016, 4(1): 68-92. doi: 10.3934/energy.2016.1.68
    [8] Sunbong Lee, Shaku Tei, Kunio Yoshikawa . Properties of chicken manure pyrolysis bio-oil blended with diesel and its combustion characteristics in RCEM, Rapid Compression and Expansion Machine. AIMS Energy, 2014, 2(3): 210-218. doi: 10.3934/energy.2014.3.210
    [9] Hasanudin Hasanudin, Wan Ryan Asri, Utari Permatahati, Widia Purwaningrum, Fitri Hadiah, Roni Maryana, Muhammad Al Muttaqii, Muhammad Hendri . Conversion of crude palm oil to biofuels via catalytic hydrocracking over NiN-supported natural bentonite. AIMS Energy, 2023, 11(2): 197-212. doi: 10.3934/energy.2023011
    [10] Boua Sidoine KADJO, Mohamed Koïta SAKO, Kouadio Alphonse DIANGO, Amélie DANLOS, Christelle PERILHON . Characterization and optimization of the heat treatment of cashew nutshells to produce a biofuel with a high-energy value. AIMS Energy, 2024, 12(2): 387-407. doi: 10.3934/energy.2024018
  • A nonlinear continuum phase field theory is developed to describe amorphization of crystalline elastic solids under shear and/or pressure loading. An order parameter describes the local degree of crystallinity. Elastic coefficients can depend on the order parameter, inelastic volume change may accompany the transition from crystal to amorphous phase, and transitional regions parallel to bands of amorphous material are penalized by interfacial surface energy. Analytical and simple numerical solutions are obtained for an idealized isotropic version of the general theory, for an element of material subjected to compressive and/or shear loading. Solutions compare favorably with experimental evidence and atomic simulations of amorphization in boron carbide, demonstrating the tendency for structural collapse and strength loss with increasing shear deformation and superposed pressure.


    With the rapid technological advancement today, the access to realistic 3D human shapes is of great importance in both computer vision and graphics, and has various applications in different industries including virtual game design, film making, bioinformatics[41], healthcare[38], and especially, those related to garment design. Some applications involve fitting predictions, virtual try-on simulations[13,16,23,34] or size recommendations[5,6,7], that help to recommend relevant clothing which would fit specific occasions or fashion trends for online customers. Such applications require a critical ingredient on digital transformation from humans bodies to digital 3D shapes, such that the shapes maintain some of the main features from human bodies.

    The traditional approaches to access reliable digital information of a human body are through laser range scanners[1], stereo reconstruction[18,22,33] or structured light methods for 3D sensing[11,20,24]. However, considering the cost of data storage, network transmission and expensive scanning equipment, it is rather unpractical to scan individuals for each application. Hence many studies have been done to generate 3D human shapes based on partial input information. These prior systems can be mainly classified into three types: marker-based systems, silhouette-based systems and measurement-based systems. Marker-based system estimates dynamic 3D human body shapes by capturing a sparse set of marker positions. These techneqiues proceed by using a single static scan or multiple scans and a marker motion capture sequence of the person[4]. For the static case, silhouette-based system estimates human body shapes based on a set of input images by fitting the silhouette in each view [3,8,9,15]. Some apporaches also combine with machine learning that build a correlation between a training dataset of 3D body shapes and a set of 2D images, and then predict a shape based on the correlation[14].

    Although marker-based systems and silhouette-based systems could yield satisfactory reconstructions on 3D human body shapes under tight dresses or naked human shapes, most of the schemes are so computationally expensive and the results are easily affected if heavy or loose clothes are worn. To overcome these difficulties, a great deal of efforts have been devoted to the investigation of simple and fast measurement-based systems [17,32,35]. Typically, one considers the landmarks or circumferences from the human structures at specific locations as characteristic values. Since such characteristic values are linear or curvilinear, they are relatively invariant to articulation changes than those silhouettes measurements. If the set of characteristic values is well selected, one can achieve meaningful estimation for both global and local body shapes. Kart et al. (2011) built a system which only requires to input some personal information, such as weight, height and age as well as a 2D photograph. The decision algorithm then determines the human shape according to the measurements and the body mass index (BMI)[12]. Seo et al. [35] presented a human body generation apporach by taking the anthropometric measurements, e.g., stature, crotch length, arm length, neck girth, chest/bust girth, underbust girth, waist girth and hip girth as input. They derived the relationship between the input characteristic values and the preprocessing database of 3D scanned data of human body models by using radial basis interpolation. At run-time, the system generates new human body shapes from the user input characteristic values by fitting the template model onto each scanned data.

    In this paper, we develop a completely novel methodology for the geometric body generation, which fulfils the following two basic requirements: (ⅰ) the geometric body generation is automatically determined by the input characteristic sets; (ⅱ) the predicted geometric shape fits for all input characteristic values and moreover it can well approximate the exact geometric body possessing the aforesaid characteristic values. The proposed method follows a machine-learning flavour that generates the inferred geometric body with the customized characteristic parameters from a training dataset. The training dataset consists of some preprocessed body shapes associated with appropriately sampled characteristic parameters. One of the critical ingredients and novelties of our method is the borrowing of inverse scattering techniques in the theory of wave propagation to the body generation. This is done by establishing a delicate one-to-one correspondence between a geometric body and the far-field pattern of a source scattering problem governed by the Helmholtz system. It in turn enables us to establish a one-to-one correspondence between the geometric body space and the function space defined by the far-field patterns. Hence, the far-field patterns can act as the shape generators. The shape generation with prescribed characteristic parameters is achieved by first manipulating the shape generators in the function space and then reconstructing the corresponding geometric body from the obtained shape generator by a stable multiple-frequency Fourier method. The proposed method is in sharp difference from the existing methodologies in the literature, which usually treat the human body as a suitable Riemannian manifold and the generation is based on non-Euclidean approximation and interpolation. In fact, in all of the literature mentioned earlier on manifold learning of body generation, one typically uses Principal Component Analysis (PCA) or Principal Geodesic Analysis (PGA). PCA and PGA are used for optimal reduction of the data and thus efficient deformation by computing statistics on Euclidean manifolds or non-Euclidean manifolds can be achieved; see [35,10] and the references therein for more relevant discussion. In our new approach, the shape generator enables us to train the learning dataset via the algebraic operations in the shape space directly without dealing with the deformation of the manifold meshes between geometric shapes.

    The rest of the paper is organized as follows. In section 2, we provide rigorous mathematical formulations of characteristic values and shape space. Section 3 introduces the notion of shape generator via the inverse source scattering associated with the Helmholtz system. In section 4, we present the mathematical setup of the geometric body generation from a machine-learning perspective. Section 5 is devoted to the development of the new method for the shape generation. In section 6, we present several two- and three-dimensional numerical examples to show the effectiveness and efficiency of our method. The paper is concluded in section 7 with some relevant discussion.

    In this section, we present some preliminary knowledge on the shape manifold theory that shall be needed in our subsequent study of body generation. Generally speaking, a geometric shape or a geometric body is a topological $ n $ -manifold, $ n\in\mathbb{N} $ , equipped with certain shape descriptors, which give the full information to describe the geometric shape. We call such shape descriptors as characteristic values. We have the following formal definition.

    Definition 2.1. Let $ D $ be a topological $ n $ -manifold with $ n\in\mathbb{N} $ . Let $ \Lambda_D: = \{\lambda^{(j)}\}_{j\in\mathscr{C}} $ be a set of parameters associated with $ D $ that are invariant with respect to isometric deformations and are independent to the parametrizations of $ D $ . Here, the cardinality $ \mathscr{C} $ might be finite or infinite. $ \Lambda_D $ is said to be a characteristic set of $ D $ if it uniquely determines $ D $ . $ D $ and its characteristic set $ \Lambda_D $ , written as the $ (D, \Lambda_D) $ is referred to as a geometric shape or a geometric body.

    Clearly, Definition 2.1 includes much general geometric objects. However, for the present study, we are mainly concerned with the case that $ D $ can be embedded into $ \mathbb{R}^d $ , $ d = 2, 3 $ , as a bounded domain. That means, we exclude some interesting cases such as $ D $ is a Riemannian surface with boundary in $ \mathbb{R}^3 $ . Nevertheless, our study is general enough to include the human body as a specific case.

    In Definition 2.1, the set of characteristic values is typically a set of measurements which gives a systematic characterization of the size, shape and composition of a geometric object for us to determine the shape of the object. For example, when considering a rectangular object, once can introduce a set of characteristic values containing its height, width and length, which provide all details to determine a unique rectangular shape. Expanding the same idea to human body shapes, one could also use characteristic sets to represent them. There are many different ways to represent a human shape. We would try to group those characteristic values into four main catagories, including Eucidean distance, geodesic distance, circumference and ratio. The Eucidean or geodisic distance is linear or curvilinear distance between two points on the human model, such as stature, crotch length, arm length, shoulder breadth etc. The circumference can be computed by the horizontal girth of the body, such as neck girth, chest/bust girth, under-bust girth, waist girth, hip girth, etc. The ratio can be information of weight, Body Mass Index, muscle and fat rate. In spite of the above characteristic values, one can also consider some pure measurements such as age or gender as characteristic values.

    The full set of characteristic values gives the complete information of a geometric shape without lossing any information. It is easy to imagine that the cardinality of a set of characteristic values depends on the complexity of a shape. Hence, the number of characteristic values required can be considered as the dimensionality of the geometric shape. For those complicated objects, like human shapes, it may require infinite set of characteristic values for accurate formulations. Due to practical reasons, one can consider a reasonable truncation of an infinite characteristic set into a finite one for a complicated geometric shape. In doing so, we can consider our study in the following product space

    $ S:=D×V, $ (2.1)

    where $ \mathscr{D} $ is composed of all the bounded domains in $ \mathbb{R}^d $ and $ \mathscr{V} $ is an $ M $ -dimensional vector space containing the characteristic values. In fact, in the present study, the characteristic values are usually real numbers and one can take $ \mathscr{V} = \mathbb{R}^{M} $ with $ M\in\mathbb{N} $ . $ \mathcal{S} $ is referred to as the geometric shape space. According to the (approximate) one-to-one correspondence between a geometric shape and its characteristic values in Definition 2.1, we readily see that all the shape information can be obtained by a single point of this $ M $ -dimensional vector space $ \mathscr{V} $ . By adjusting the characteristic values, we can obtain new geometric shapes and this is a key ingredient in our human body interpolation.

    In the previous section, we introduce the important notion of shape space for our study. We proceed to introduce another critical ingredient, shape generator, for our subsequent study of the geometric body generation. In fact, the generation of a new geometric shape shall be based on algebraic interpolation of exemplar models from the shape space. If the algebraic operations are to be conducted directly in the shape space, dealing with geometric deformations of manifolds, one would certainly encounter very complicated and tedious calculations and manipulations because of the lack of global parametrizations for the non-Eucidean shapes involved. The shape generator can overcome this challenge by bridging the geometric shape space and the function space. To that end, we next introduce the inverse scattering problem in finding an active source from its generated far-field pattern.

    Let $ f:\mathbb{R}^d\mapsto\mathbb{C} $ be a function having a compact support, $ f = \chi_D \varphi $ , where $ D\subset\mathbb{R}^d $ is a bounded domain and $ \varphi\in L^\infty(\mathbb{R}^d) $ . The set $ D $ is the external shape of $ f $ while $ \varphi $ describes the intensity of the source at various points in $ D $ . We assume that $ \varphi $ and $ D $ do not depend on the wavenumber $ k\in\mathbb{R}_+ $ . In other words we are considering monochomatic scattering. The source $ f $ produces a scattered wave $ u\in H^2_{loc}(\mathbb{R}^d) $ given by the unique solution to

    $ (Δ+k2)u=f,limrrd12(rik)u=0, $ (3.1)

    where $ r = |x| $ for $ x\in\mathbb{R}^d $ . The limit in (3.1) is known as the Sommerfeld radiation condition which characterizes the outgoing nature of the radiating wave. By the limiting absorption principle (cf. [21]), the solution to (3.1) can be computed as follows,

    $ u=(Δ+k2)1f=limε+0(Δ+(kiε)2)1f=limε+0Rdeixξˆf(ξ)|ξ|2(kiε)2 dξ, $ (3.2)

    where

    $ ˆf(ξ):=Ff(ξ)=(2π)dRnf(x)eiξx dx $ (3.3)

    signifies the Fourier transform of $ f $ . Inverting the Fourier transform in (3.2), one has the following integral representation,

    $ u=(Δ+k2)1f:=i4(k2π)d22Rd|xy|2d2H(1)d22(k|xy|)f(y) dy, $ (3.4)

    where $ H_{(d-2)/2}^{(1)} $ is the first-kind Hankel function of order $ (d-2)/2 $ . Stationary phase applied to (3.4) yields that

    $ u(x)=eik|x||x|(d1)/2Cd,kRdeikˆxyf(y) dy+O(|x|d2),|x|, $ (3.5)

    where $ \hat x: = x/|x|\in\mathbb{S}^{d-1} $ , $ x\in\mathbb{R}^d\backslash\{0\} $ , and

    $ C_{d, k} = \frac{-i}{\sqrt{8\pi}}\left(\frac{k}{2\pi} \right)^{\frac{d-2}{2}} e^{-\frac{(d-1)\pi}{4}i}. $

    The far-field pattern of $ u $ is given by

    $ u(ˆx,k;f):=Cd,kRdeikˆxyf(y) dy=(2π)dCd,kFf(kˆx)L2(Sd1). $ (3.6)

    It is obvious that $ u_\infty $ is (real) analytic in both $ \hat x $ and $ k $ . Hence, if $ u_\infty(\hat x, k) $ is known on any open portion of $ \mathbb{S}^{d-1}\times\mathbb{R}_+ $ , then it is known on the whole set by analytic continuation.

    The inverse source scattering problem is concerned with the recovery of $ f = \varphi\chi_D $ by knowledge of $ u_\infty(\hat x, k; f) $ for $ (\hat x, k)\in\Sigma $ , where $ \Sigma $ is an open subset of $ \mathbb{S}^{d-1}\times\mathbb{R}_+ $ . According to our discussion above, without loss of generality, we always assume that $ \Sigma = \mathbb{S}^{d-1}\times\mathbb{R}_+ $ in what follows. The inverse source problem arises in a variety of important applications including detection of hazardous chemicals, medical imaging, photoacoustic and thermoacoustic tomography, brain imaging, artificial intelligence in gesture computing and others. We refer to the recent articles [25,26,39,40] by two of the authors of this article for some recent developments on the inverse source problem.

    Next, let us consider a specific case by assuming a source supported in a domain $ D $ with a constant density $ 1 $ . Then clearly by (3.6), there is a one-to-one correspondence between $ D $ and $ u_\infty(\Sigma; D): = \{u_\infty(\hat x, k; 1\cdot\chi_D)\}_{(\hat x, k)\in\Sigma}\in L^2(\mathbb{S}^{d-1}\times\mathbb{R}_+) $ in the sense that for two domains $ D_1 $ and $ D_2 $ ,

    $ u(Σ;D1)=u(Σ;D2)if and only ifD1=D2. $ (3.7)

    Based on (3.7), we next introduce

    Definition 3.1. For a geometric shape $ (D, \Lambda_D)\in\mathcal{S} $ ,

    $ u_\infty(\Sigma; D): = \{u_\infty(\hat x, k; 1\cdot\chi_D)\}_{(\hat x, k)\in\Sigma}\in L^2(\mathbb{S}^{d-1}\times\mathbb{R}_+) $

    defined via the Helmholtz system (3.1) is called a shape generator for $ D $ .

    Remark 3.1. By Definition 3.1, a geometric body $ D $ can be completely determined by a shape generator $ u_\infty(\Sigma) $ . Since $ u_\infty(\Sigma) $ is from a function space, this paves the way for the new body generation through function interpolations.

    Remark 3.2. By (3.6), we know the far-field pattern is actually the Fourier transform of the source density up a dimensional constant. However, introducing the shape generator via the inverse scattering approach shall provide more physical insights in our study, and moreover it enables us to borrow ideas from the inverse scattering literature of recovering the geometric shape $ D $ from the associated far-field pattern. This also paves the way of extending the idea by using other inverse scattering models that have such one-to-one correspondence between geometric shapes and far-field patterns; see more relevant discussion in section 7.

    In this section, we introduce the mathematical formulation of the geometric body generation for our study from a machine learning perspective. For a geometric shape $ (D, \Lambda_{D}) $ with the associated shape generator $ u_\infty(\Sigma; D) $ , the pair of the high dimensional variables, written as $ \{(\Lambda_{D}, u_\infty(\Sigma; D))\} $ , is referred to as an input-output pair. Let $ \{(\Lambda_{D_i}, u_\infty(\Sigma; D_i))\}_{i\in \mathscr{N}} $ with $ \mathscr{N} = \{0, 1, \dots, N_{pair}\} $ be a set of input-output pairs associated with the characteristic sets $ \{\Lambda_{D_i}\}_{i\in \mathscr{N}} $ . Here the input characteristic sets $ \Lambda_{D_i} $ are introduced as

    $ ΛDi:={λ(j)i}jC,iN, $ (4.1)

    with $ \mathscr{C} = \{1, \dots, M\} $ . In (4.1), the notation $ \lambda^{(j)}_{i} $ represents the characteristic value of the $ i $ -th geometric shape $ D_i $ in the $ j $ -th direction of its characteristic set. Here, the cardinality $ \mathscr{N} $ is finite. The training dataset of the geometric body generation is introduced to be

    $ Z:={(ΛDi,u(Σ;Di))}iN $ (4.2)

    with

    $ Λ=(λ(1),,λ(M))RM and u(Σ)L2(Sd1×R+). $ (4.3)

    The training dataset consists of certain pre-sampled geometric shapes with statistically well selected characteristic values. The corresponding shape generator of a specific body in the training dataset can also be pre-calculated and stored. The main goal of our study is to first infer a learning model from the training dataset, $ T_\mathbf{Z}: \mathbb{R}^M\rightarrow L^2(\mathbb{S}^{d-1}\times\mathbb{R}_+) $ that fulfils the following requirements:

    1. It fits the training data well in the sense that

    $ TZ(ΛDi):=ˆu(Σ;Di)u(Σ;Di),DiZ. $ (4.4)

    2. It can be used to infer the shape generator for a given new shape with prescribed characteristic values, namely,

    $ ˆu(Σ;Dnew):=TZ(ΛDnew), $ (4.5)

    and with a statically well selected training dataset, it is justifiable to expect that

    $ ˆu(Σ;Dnew)u(Σ;Dnew), $ (4.6)

    where $ u_\infty(\Sigma; D_{\mathrm{new}}) $ is the shape generator for $ D_{\mathrm{new}} $ .

    If a learning model can be achieved that fulfils the two requirements as described above, then the body generation can be proceeded as follows. For a given new set of characteristic values, one first generates the learned shape generator as in (4.5). By a certain inverse scattering approach, one can then reconstruct the (approximate) shape $ D_{\mathrm{new}} $ from the corresponding shape generator $ \hat{u}_\infty(\Sigma; D_{\mathrm{new}}) $ . In the next section, we shall develop the two critical ingredients in the body generation procedure described above, namely, the learning model and the reconstruction method. To be more definite and specific, we first introduce the following definition from a machine learning perspective.

    Definition 4.1. (Body learning model) Given a training dataset

    $ Z:={(ΛDi,u(Σ;Di))}iN. $ (4.7)

    Let $ \mathcal{H} $ be a compact subset of $ L^2(\mathbb{S}^{d-1}\times\mathbb{R}_+) $ . $ T_\mathbf{Z}\in\mathcal{H} $ (with specified coefficients $ C $ ) is said to be the best fit learning model associated with the training dataset $ \mathbf{Z} $ if it is the minimizer of the following optimization problem,

    $ minTZH1Npair+1iNTZ(ΛDi)u(Σ;Di)2H. $ (4.8)

    According to Definition 4.1, the choice of the learning subspace $ \mathcal{H} $ plays a critical role. However, we note that the shape generator is actually (real) analytic in all of its arguments. Hence, instead of solving the computationally costly optimization problem (4.8), we can make use of the functional interpolation to produce a well-rounded shape learning model. This is one of the main advantages of introducing the shape generator through the inverse scattering model. In the next section, for a given training dataset as in (4.7), we shall derive a learning model using the cubic B-spline interpolation through the use of the high-dimensional data-points (4.4). For the reconstruction of the approximate body shape from the shape generator obtained through the learning model, we shall make use of a multiple-frequency Fourier method, and it can also produce an efficient and stable recovery. Throughout, we assume that the characteristic values in the training dataset is statistically well selected and it is not the focus of the present article.

    In this section, we develop the details of our scheme for the geometric body generation following the general discussion made in the previous section. We first derive the learning model through the functional interpolation of the high-dimensional data in the training dataset. To that end, we present some preliminary knowledge on the cubic B-spline, and we also refer to [2,19,36,37] for more relevant discussion on the cubic B-spline.

    Consider the training dataset (4.2). Let the $ M $ sets of unique grids in the directions of $ \{\lambda^{(j)}\}_{j\in \mathscr{C}} $

    $ Δj={λ(j)0,,λ(j)kj}jC,kjN, $ (5.1)

    define on the intervals $ [a_j, b_j] $ as $ M $ sets of points $ \lambda_{g_j}^{(j)}\in[a_j, b_j]\subset \mathbb{R} $ , where $ g_j\in \{0, \dots, k_j\} $ and $ a_j = \lambda_{0}^{(j)} < \lambda_{1}^{(j)} < \dots < \lambda_{k_j}^{(j)} = b_j, j = 1, \dots, M $ . Here, $ k_j $ is the greatest number of distinct characteristic values in the $ j $ -th direction of the characteristic set. We remark that if the characteristic values of the training dataset are all collected in distinct values, then $ k_j $ is actually the last index of the training dataset, $ N_{pair} $ . However, the training dataset might be collected in such a way that some body shapes may possess the same characteristic value in the $ j $ -th direction, and hence $ k_j $ is usually smaller than $ N_{pair} $ .

    With the above notation, the training dataset stored as the array in (4.2) can be represented as elements on the grid mesh corresponding to the characteristic numbers $ \lambda^{(j)} $ as described above. The interpolation data are the corresponding shape generators and are written as

    $ u(Σ;Di):=Ug1,g2,,gML2(Sd1×R+),iN, $ (5.2)

    where $ g_j = 0, \dots, k_j, j = 1, \dots, M. $ In the subsequent study, we shall stick to the same notation $ g_1, g_2, \dots, g_M $ to represent the linear indexing. The following example demonstrates a real application for the human body generation.

    Example 5.1. The training dataset consists of 20 bodies with two characteristic values as consideration, say, height and relative weight. Here, the height and relative weight are the two directions of the grids, $ \lambda^{(height)} $ and $ \lambda^{(weight)} $ . Suppose the heights of the sampled bodies are given by 1.5 m, 1.6 m, 1.7 m, 1.8 m, 1.9 m and the relative weights of the sampled bodies are given by 60%, 80%, 100%, 120%. Then the first grid $ \Delta_{height} = \{\lambda_0^{(height)}, \lambda_1^{(height)}, \lambda_2^{(height)}, \lambda_3^{(height)}, \lambda_4^{(height)}\} = \{1.5, 1.6, 1.7, 1.8, 1.9\} $ and the second grid $ \Delta_{weight} = \{\lambda_0^{(weight)}, \lambda_1^{(weight)}, \lambda_2^{(weight)}, \lambda_3^{(weight)}\} = \{0.6, 0.8, 1, 1.2\} $ . The interpolation data are actually stored as listed in Table 1; e.g., $ U_{0, 0} = u_\infty(\Sigma; D_1), U_{0, 3} = u_\infty(\Sigma; D_4), U_{1, 3} = u_\infty(\Sigma; D_{8}) $ .

    Table 1.  Training dataset with the given characteristic grids $ (\lambda^{(height)}, \lambda^{(weight)})$.
    $ 2^{nd}$ grid
    $ \lambda_0^{(weight)}$ $ \lambda_1^{(weight)}$ $ \lambda_2^{(weight)}$ $ \lambda_3^{(weight)}$
    $ \lambda_0^{(height)}$ $ (\Lambda_{D_1}, u_\infty(\Sigma, D_1))$
    $ =(1.5, 0.6, U_{0, 0})$
    $ (\Lambda_{D_2}, u_\infty(\Sigma, D_2))$
    $ =(1.5, 0.8, U_{0, 1})$
    $ (\Lambda_{D_3}, u_\infty(\Sigma, D_3))$
    $ =(1.5, 1, U_{0, 2})$
    $ (\Lambda_{D_4}, u_\infty(\Sigma, D_4))$
    $ =(1.5, 1.2, U_{0, 3})$
    $ 1^{st}$ grid$ \lambda_1^{(height)}$ $ (\Lambda_{D_5}, u_\infty(\Sigma, D_5))$
    $ =(1.6, 0.6, U_{1, 0})$
    $ (\Lambda_{D_6}, u_\infty(\Sigma, D_6))$
    $ =(1.6, 0.8, U_{1, 1})$
    $ (\Lambda_{D_7}, u_\infty(\Sigma, D_7))$
    $ =(1.6, 1, U_{1, 2})$
    $ (\Lambda_{D_8}, u_\infty(\Sigma, D_8))$
    $ =(1.6, 1.2, U_{1, 3})$
    $ \lambda_2^{(height)}$ $ (\Lambda_{D_9}, u_\infty(\Sigma, D_9))$
    $ =(1.7, 0.6, U_{2, 0})$
    $ (\Lambda_{D_{10}}, u_\infty(\Sigma, D_{10}))$
    $ =(1.7, 0.8, U_{2, 1})$
    $ (\Lambda_{D_{11}}, u_\infty(\Sigma, D_{11}))$
    $ =(1.7, 1, U_{2, 2})$
    $ (\Lambda_{D_{12}}, u_\infty(\Sigma, D_{12}))$
    $ =(1.7, 1.2, U_{2, 3})$
    $ \lambda_3^{(height)}$ $ (\Lambda_{D_{13}}, u_\infty(\Sigma, D_{13}))$
    $ =(1.8, 0.6, U_{3, 0})$
    $ (\Lambda_{D_{14}}, u_\infty(\Sigma, D_{14}))$
    $ =(1.8, 0.8, U_{3, 1})$
    $ (\Lambda_{D_{15}}, u_\infty(\Sigma, D_{15}))$
    $ =(1.8, 1, U_{3, 2})$
    $ (\Lambda_{D_{16}}, u_\infty(\Sigma, D_{16}))$
    $ =(1.8, 1.2, U_{3, 3})$
    $ \lambda_4^{(height)}$ $ (\Lambda_{D_{17}}, u_\infty(\Sigma, D_{17}))$
    $ =(1.9, 0.6, U_{4, 0})$
    $ (\Lambda_{D_{18}}, u_\infty(\Sigma, D_{18}))$
    $ =(1.9, 0.8, U_{4, 1})$
    $ (\Lambda_{D_{19}}, u_\infty(\Sigma, D_{19}))$
    $ =(1.9, 1, U_{4, 2})$
    $ (\Lambda_{D_{20}}, u_\infty(\Sigma, D_{20}))$
    $ =(1.9, 1.2, U_{4, 3})$

     | Show Table
    DownLoad: CSV

    Let $ S_3(\Delta_j), j = 1, \dots, M $ be a function subspace of $ C^2([a_j, b_j]) $ consisting of one dimensional, complex-valued functions in the direction of $ \lambda^{(j)}, j = 1, \dots, M $ on the bounded interval $ [a_j, b_j] $ . The function in $ S_3(\Delta_j), j = 1, \dots, M $ is piecewise polynomial of degree 3 on every subinterval $ [\lambda_{g_j-1}^{(l)}, \lambda_{g_j}^{(j)}] $ , where $ g_j = 1, \dots, k_j, j = 1, \dots, M $ . Then we introduce a function subspace of multidimensional and complex-valued $ C^2([a_j, b_j]) $ functions as

    $ S3(Δ1,,ΔM), $ (5.3)

    on each rectangular grid

    $ Ig1,g2,,gM:=jC[λ(j)gj,λ(j)gj+1] $ (5.4)

    for all $ 0\leq g_j\leq k_j-1, j = 1, \dots, M $ that are piecewise polynomials of degree $ 3 $ on every interval. For easy reference we provide the definition of B-splines.

    Definition 5.1. The sets of $ k_j+k, j = 1, \dots, M $ , B-spline basis functions $ \{B_{l, k}(\lambda^{(j)})\}_{l = 1}^{k_j+k} $ of degree $ k $ of the function space $ S_k(\Delta_j) $ are defined based on concurrent boundary knots vectors with Cox-deBoor recurrence[19],

    $ Bl,k(λ(j))=λ(j)λ(j)l1kλ(j)l1λ(j)l1kBl1,k1(λ(j))+λ(j)lλ(j)λ(j)lλ(j)lkBl,k1(λ(j)) $ (5.5)

    with

    $ Bl,0(λ(j))={1if λ(j)lλ(j)<λ(j)l+10else, $ (5.6)

    for $ l = 1, \dots, k_j+k $ , where $ k = 3 $ is for the cubic B-spline and $ \lambda_l^{(j)} $ are elements of the knot vectors, satisfying the relation $ \lambda_l^{(j)} < \lambda_{l+1}^{(j)} $ .

    All methods are in the following using splines with a $ k + 1 $ regular knot vector, and the interior knots are the grid points. Based on Definition 5.1, we next introduce a general learning model for the geometric shape generation through the multidimensional cubic B-spline interpolation.

    Body learning model Ⅰ. Given the training dataset $ \mathbf{Z}: = \{(\Lambda_{D_i}, u_\infty(\Sigma; D_i))\}_{i\in \mathscr{N}} $ , the learning model $ T_\mathbf{Z}\in S_3(\Delta_1, \dots, \Delta_M) $ at $ \Lambda = (\lambda^{(1)}, \dots, \lambda^{(M)}) $ for the geometric body generation associated with the sets of the grids $ \{\Delta_j\}_ {j\in \mathscr{C}} $ is defined as follows

    $ TZ(Λ)=k1+3g1=1kM+3gM=1cg1,g2,,gMjCBjgj,3(λ(j)),λ(j)[aj,bj], $ (5.7)

    which satisfies the following conditions by (4.4)

    $ TZ(ΛDi)=k1+3g1=1kM+3gM=1cg1,g2,,gMjCBjgj,3(λ(j)i)=u(Σ;Di), $ (5.8)

    where $ c_{g_1, g_2, \dots, g_{M}} $ with $ g_j = 0, \dots, k_j, j = 1, \dots, M $ are the coefficients to be determined from the training dataset $ \mathbf{Z} $ , $ B_{g_j, 3}^j(\lambda^{(j)}), j\in \mathscr{C} $ are the B-spline basis functions of degree 3 defined in (5.5), and $ k_j, j = 1, \dots, M $ is the number of different characteristic values in each direction.

    Remark 5.1. In learning model Ⅰ, (5.7) presents a general form of the learning model $ T_\mathbf{Z} $ for the geometric body generation associated with the non-uniform grids $ \{\Delta_j\}_ {j\in \mathscr{C}} $ . The learning model eventually generates a B-spline interpolation with the associated spacing for each segment. If the training dataset consists of equidistant grids, we can derive a faster and easier learning model and this shall be provided in the next subsection.

    In this subsection, we derive a learning model for a special case with the training dataset consisting of equidistant grids. By (5.1), for the $ M $ set of equidistant grids $ \Delta_j = \{\lambda_0^{(j)}, \cdots, \lambda_{k_j}^{(j)}\}_{j\in \mathscr{C}} $ with additional conditions

    $ λ(j)gj=aj+gjhj,hj=bjajkj,,gj=1,,kj, $ (5.9)

    the B-spline basis function $ \beta^k(t) $ of degree $ k $ is a symmetrical, bell-shaped function constructed from $ k+1 $ times self-convolution of the $ \beta^0(t) $ basis function of degree zero which is a centered rectangle around origin [37]

    $ β0(t)={1,12<t<1212,|x|=120,otherwise, $ (5.10)
    $ βk(t)=β0(t)β0(t)(k+1) times. $ (5.11)

    The centered symmetric B-spline of degree $ k $ has an explicit expression [36]

    $ βk(t)=1k!k+1j=0Ck+1j(1)j(t+k+12j)k+, $ (5.12)

    where the function $ x_+ $ is defined as follows

    $ x+={x,forx>0,0,otherwise. $ (5.13)

    In this paper, we are particular intereted in the cubic B-spline. By (5.12), the closed-form representation of the cubic B-spline basis function can be also expressed as

    $ β3(t)=16{(2|t|)31<|t|2,46|t|2+3|t|3,|t|1,0,elsewhere, $ (5.14)

    which is used for preforming the interpolation. Then we choose the interpolation kernels to be

    $ Lgj(λ(j))=β3(λ(j)),kj=1,,Npair+3, $ (5.15)

    as the basis of $ S_3(\Delta_j) $ in the $ \lambda^{(j)} $ -direction such that $ L = \{L_1, L_2, \dots, L_{k_j+3}\} $ is a basis of the $ (k_j+3) $ -dimensional space $ S_3(\Delta) $ and hence, the basis of the $ \prod_{j\in\mathscr{C}}(k_j+3) $ -dimensional space $ S_3(\Delta_1, \dots, \Delta_M) $ in the directions of $ \lambda^{(1)}, \dots, \lambda^{(M)} $ is given by

    $ {Lg1Lg2LgM|gj{1,,Npair+3},jC}. $ (5.16)

    Based on (5.16), (5.1) and (5.9), we next introduce the learning model for the uniform case.

    Body learning model Ⅱ. Given the training dataset $ \mathbf{Z}: = \{(\Lambda_{D_i}, u_\infty(\Sigma; D_i))\}_{i\in \mathscr{N}} $ , the learning model $ T_\mathbf{Z}\in S_3(\Delta_1, \dots, \Delta_M) $ at $ \Lambda = (\lambda^{(1)}, \dots, \lambda^{(M)}) $ for the geometric body generation associated with the sets of equidistent grids $ \{\Delta_j\}_ {j\in \mathscr{C}} $ defined in (5.9) is defined as follows

    $ TZ(Λ)=k1+3g1=1kM+3gM=1cg1,g2,,gMjCLjgj(λ(j)), $ (5.17)

    which is required to satisfy the following conditions

    $ TZ(ΛDi)=k1+3g1=1kM+3gM=1cg1,g2,,gMjCLjgj(λ(j)i)=u(Σ;Di), $ (5.18)

    where $ c_{g_1, g_2, \dots, g_{M}} $ with $ g_j = 0, \dots, k_j, j = 1, \dots, M. $ are the coefficients to be determined from the training dataset $ \mathbf{Z} $ , $ L_{g_j}^j(\lambda^{(j)}) $ are B-spline basis functions of degree 3 defined in (5.15) and $ k_j, j = 1, \dots, M $ . is the number of different characteristic values in each direction.

    In learning models Ⅰ and Ⅱ, the learning functionals for the non-uniform and uniform case are respectively considered in the $ \prod_{j\in \mathscr{C}}(k_j+3) $ -dimensional space $ S_3(\Delta_1, \dots, \Delta_M) $ . $ \prod_{j\in \mathscr{C}}(k_j+3) $ interpolation conditions are required to determine the coefficients $ c_{i_1, i_2, \dots, i_{M}} $ in the training models. However, there are only $ \prod_{j\in \mathscr{C}}(k_j+1) $ shape generators to specify $ \prod_{j\in \mathscr{C}}(k_j+1) $ conditions in (5.8) or (5.18). To obtain a unique correlation between the characteristic values and the shape generator, we need to add $ \prod_{j\in \mathscr{C}}2 = 2^{M} $ conditions, which define the second-order derivatives of the spline function at the boundary $ a_j $ and $ b_j $ to be equal to 0 and lead to a natural spline.

    With the learning models Ⅰ and Ⅱ established in the previous subsections, for an input new set of characteristic values $ \Lambda_{D_{new}} $ associated with a new geometric body $ D_{new} $ , the unknown shape generator can be generated as follows,

    $ TZ(ΛDnew)=k1+3g1=1kM+3gM=1cg1,g2,,gMjCLjgj(λ(j)new)=ˆu(Σ;Dnew)u(Σ;Dnew), $ (5.19)

    where the coefficients $ c_{g_1, g_2, \dots, g_{M}} $ in (5.19) could be determined by solving the natural spline problem.

    In this subsection, we briefly outline the Fourier method for the reconstruction of geometry shape $ D_{{new}} $ by using the shape generators $ {u}_\infty(\Sigma; D_{new}) $ .

    Define the periodic Sobolev space by

    $ Hσ(Rd)={gL2(Rd):(1+|ξ|2)σ2ˆg(ξ)L2(Rd),}, $

    where $ \sigma \geq 1 $ , $ \xi\in \mathbb{Z}^d $ and $ \widehat g(\xi) $ denote the Fourier coefficients of $ g $ . Suppose that $ f\in H^{\sigma}(\mathbb{R}^d) $ has a compact support in domain $ V_0 = (-a/2, a/2)^d, \, (a > 0) $ , then the Fourier transform of $ f $ is represented by

    $ ˆfξ=1adV0f(x)¯ϕξ(x)dx, $ (5.20)

    where the overbar stands for the complex conjugate and the Fourier basis functions are given by

    $ ϕξ(x)=exp(i2πaξx),ξZd. $

    Definition 5.2 (Admissible wavenumbers and observation directions). Let $ \mu $ be a sufficiently small positive constant such that $ 0 < \mu < 1 $ and

    $ ξ0:={(μ,0),d=2,(μ,0,0),d=3, $

    then the admissible set of wavenumbers is defined by

    $ K:={2πa|ξ|:ξZ3{0}}k0, $

    correspondingly, the admissiable set of observation directions is given by

    $ X:={ξ|ξ|:ξZ3{0}}ˆx0, $

    where $ k_0 = 2\pi|\xi_0|/a $ and $ \hat x_0 = \xi_0/|\xi_0| $ for $ \xi = 0\in \mathbb{Z}^d $ .

    Due to $ \mathrm{supp}\, S \subset\subset V_0 $ , for $ \xi\in \mathbb{Z}^{3}\backslash \{0\} $ , the far-field pattern defined in (3.6) can be written as

    $ u(ˆx,k;f)=Cd,kRdeikˆxyf(y) dy=Cd,kV0eikˆxyf(y) dy, $ (5.21)

    where $ k\in \mathbb{K} $ and $ \hat x\in \mathbb{X} $ depend on $ \xi $ . Combining (5.20) and (5.21), one has

    $ ˆfξ=1adCd,ku(ˆx,k;f),ξZ3{0}. $ (5.22)

    For $ \xi = 0 $ , using the Fourier expansion of $ f $ , we derive that

    $ u(ˆx0,k0;f)=Cd,kV0eik0ˆx0yf(y) dy,=Cd,kV0(ˆf0+ξZd{0}ˆfξϕξ)¯ϕξ0(y) dy=Cd,kV0ˆf0¯ϕξ0(y) dy+Cd,kξZd{0}ˆfξV0ϕξ(y)¯ϕξ0(y) dy=adCd,ksin(μπ)μπˆf0+Cd,kξZd{0}ˆfξV0ϕξ(y)¯ϕξ0(y) dy, $

    which implies

    $ ˆf0=μπadsinμπ{u(ˆx0,k0;f)Cd,kξZ3{0}ˆfξV0ϕξ(y)¯ϕξ0(y) dy}. $ (5.23)

    Therefore, the Fourier method is to approximate $ f $ by a truncated Fourier expansion

    $ fN=ˆf0+1|ξ|Nˆfξϕξ, $

    where $ N\in \mathbb{N}_+ $ denotes the truncation order and the Fourier coefficients are given by (5.22) and (5.23). Hence the domain $ D_{new} $ is determined since the set $ D_{new} $ is the external shape of $ f $ .

    Next, we investigate the stability of the proposed Fourier method. In practical computation, there exists some noise between the shape generators $ u_\infty(\Sigma; D_{{new}}) $ and the predictions $ \hat u_\infty(\Sigma; D_{{new}}) $ , which satisfies

    $ ˆu(Σ;Dnew)u(Σ;Dnew)L2δu(Σ;Dnew)L2, $

    where $ \delta > 0 $ denotes the noise level. Noting that $ (\hat x, k)\in \mathbb{X}\times \mathbb{K}\subset \Sigma $ , then the approximation of $ f $ from predicted shape generators is given by

    $ fδN=ˆfδ0+1|ξ|Nˆfδξϕξ, $

    where

    $ ˆfδξ=1adCd,kˆu(ˆx,k),ξZ3{0}, $ (5.24)
    $ ˆfδ0=μπadsinμπ{ˆu(ˆx0,k0)Cd,kξZ3{0}ˆfδξV0ϕξ(y)¯ϕξ0(y) dy}. $ (5.25)

    Theorem 5.1. Let $ f $ be a compactly supported function in $ H^{\sigma}(\mathbb{R}^d) $ , $ \sigma \geq 1 $ , with $ \mathrm{supp} f \subset\subset V_0 $ , then we have the following estimate

    $ fδf2L2(V0)Cδ+C(τd+τ2)δ42+d. $

    where $ d\leq \tau\in \mathbb{R}_+ $ and $ C $ is a constant which depends on $ f, a, d, \mu $ .

    Proof. Using the Plancherel theorem, we have

    $ fδf2L2(D)=fδf2L2(Rd)=1a2dRd|ˆfδξˆfξ|2 dξ=1a2d(|ξ|N|ˆfδξˆfξ|2 dξ+|ξ|>N|ˆfδξˆfξ|2 dξ), $ (5.26)

    where $ N\in \mathbb{N}_+ $ . Due to $ f\in H^{\sigma}(\mathbb{R}^d) $ , that is,

    $ (1+|ξ|2)α2ˆfξL2(Rd),|α|σ. $

    It means that both $ |\xi|\widehat f_{\xi} $ and $ |\xi|\widehat f_{\xi}^{\delta} $ are bounded in $ L^2(\mathbb{R}^d) $ , so we can find $ N > 0 $ , such that

    $ |ξ|>N|ˆfδξˆfξ|2 dξ1N2|ξ|>N|ξ|2|ˆfδξˆfξ|2 dξ<C1N2, $ (5.27)

    where $ C_1 > 0 $ is a constant. For $ 1\leq |\xi|_{\infty}\leq N $ , from (5.22) and (5.24), we have

    $ |ˆfδξˆfξ|=1adCd,k|ˆu(ˆx,k)u(ˆx,k)|δadCd,k|u(ˆx,k)|=δadCd,k|Cd,kV0f(y)eikˆxy dy|δad(V0|f(y)|2 dy)12(V0|eikˆxy|2dy)12=fL2(V0)ad/2δ, $

    which implies

    $ 1|ξ|N|ˆfδξˆfξ|2 dξC2(2N+1)dδ2, $ (5.28)

    for $ C_2 = \|f\|_{L^2(V_0)}^2/ a^d $ . Define $ \xi = (\xi_{1}, \xi_2) \in \mathbb{Z}^2 $ or $ \xi = (\xi_{1}, \xi_2, \xi_3) \in \mathbb{Z}^3 $ , by a straight forward calculation, one finds that

    $ V0ϕξ(y)¯ϕξ0(y) dy={0,|ξ||ξ1|,adcosξ1πsinλπ(ξ1λ)π,|ξ|=|ξ1|. $

    For $ \xi = 0 $ , using (5.23), (5.25), (5.27), (5.28) and the last equation, it derives that

    $ |ˆfδ0ˆf0|μπadCd,ksinμπ|ˆuδ(ˆx0,k0)uδ(ˆx0,k0)|+μπadsinμπ1|ξ|N|ˆfδξˆfξ||V0ϕξ(y)¯ϕξ0(y) dy|+μπadsinμπ|ξ|N|ˆfδξˆfξ||V0ϕξ(x)¯ϕξ0(x) dy|C3δ+C2(2N+1)dδ+C1N, $ (5.29)

    where $ C_3 = {\mu \pi\|f\|_{L^2(V_0)} }/(a^{d/2}\sin \mu \pi) $ . Hence, substituting (5.27), (5.28) and (5.29) into (5.26), it deduces that

    $ fδf2L2(V0)Cδ2+CNdδ2+CN2, $

    where $ C = \max\{2 C_1, 2^{d+1}C_2, C_3^2\}/ a^{2d} $ . Furthermore, if we take $ N = \tau\delta^{-\frac{2}{2+d}} $ with $ \tau\geq d $ in Theorem 5.1, then it holds that

    $ fδf2L2(V0)Cδ+C(τd+τ2)δ42+d. $

    Let $ N = \left[\tau\delta^{-\frac{2}{2+d}}, \, \tau\geq d\right] $ , here $ [X] $ denotes the largest integer that is smaller than $ X + 1 $ . From definition 5, the truncated wavenumbers and observation directions can be written as

    $ KN:={2πa|ξ|:1|ξ|N}k0,XN:={ξ|ξ|:1|ξ|N}ˆx0. $

    Thus, the truncated Fourier expansion of $ f $ from the predictions $ \{\hat u_\infty(\Sigma; D_{{new}_q})\}_{q\in\mathscr{Q}} $ takes the form

    $ fδN:=ˆfδ0+1|ξ|Nˆfδξϕξ(x), $ (5.30)

    where

    $ ˆfδξ=1adCd,kˆu(ˆx,k),1|ξ|N,ˆfδ0=μπadsinμπ{ˆu(ˆx0,k0)Cd,k1|ξ|NˆfδξV0ϕξ(y)¯ϕξ0(y) dy}. $ (5.31)

    Motivated by above discussion, we are ready to present our novel modeling methodology for geometric shape in $ \mathbb{R}^d, d = 2, 3 $ , see Algorithm 1.

    Algorithm 1: Inverse-scattering-based geometric body generation scheme
    1: Select the parameter $ N $ , the set of admissible wavenumbers $ \mathbb{K}_{N} $ and the set of admissible observation directions $ \mathbb{X}_{N} $ .
    2: Given a training dataset $ \mathbf{Z}:=\{(\Lambda_{D_i}, u_\infty(\hat x, k; D_i))\}_{i\in \mathscr{N}} $ for $ \hat x\in X_N, k\in K_N $ , obtain the coefficients $ c_{g_1, \dots, g_{M}} $ of the learning model $ T_\mathbf{Z} $ by solving the problem of the natural spline interpolation.
    3: Given the characteristic sets $ \Lambda_{D_{new}} $ , predict the new shape generators $ \{\hat u_\infty(\hat x, k)\} $ for $ \hat x\in X_N, k\in K_N $ with the use of the learning model $ T_\mathbf{Z} $ .
    4: Compute the Fourier coefficients $ \widehat{f}_{ 0}^{\delta} $ and $ \widehat{f}_{ \xi}^{\delta} $ defined in (5.31) for $ 1\leq| \xi|_\infty\leq N_t $ .
    5: Select a sampling mesh $ \mathcal{T}_h $ in a region $ V_0 $ . For each sampling point $ z_j\in \mathcal{T} $ , calculate the imaging function $ f_{N}^{\delta} $ defined in (5.30). $ D_{new} $ is obtained as the external shape of $ f_{N}^{\delta} $ .

     | Show Table
    DownLoad: CSV

    In this section, several numerical examples of geometric body generation are conducted to show that the proposed method is effective and efficient. Here the geometric body includes general geometric shape and human body shape.

    The proposed algorithm is implemented by using Matlab 2016. The shape generator $ \{u_\infty(\hat x, k; D_i)\}_{i\in \mathscr{N}} $ is obtained by solving the direct problem of (3.1). To avoid the inverse crime, we use the quadratic finite elements on a truncated spherical domain enclosed by a PML layer. The mesh of the forward solver is successively refined till the relative error of the successive measured scattered data is below $ 0.1\% $ . Then artificial shape generators are generated by applying the Kirchhoff integral formula to the scattered data. Thus the training dataset is given by

    $ {(ΛDi,u(ˆxj,kj;Di)):ˆxjXN,kjKN}, $

    where $ j = 1, 2, \cdots, (2N+1)^{d} $ and $ i = 1, 2, \cdots, N_t $ denotes the $ i $ -th geometry shape. In what follows, we set $ \tau = 2 $ ($ \tau = 3 $ ) for $ d = 2 $ ($ d = 3 $ ) and $ \delta = 1\% $ , then we have $ N = 20 $ ($ N = 19 $ ) for $ d = 2 $ ($ d = 3 $ ).

    Next, we present the implementation of interpolation. The characteristic value set is given by $ \{\Lambda_{D_i}\} $ , where $ \Lambda_{D_i} = \{\lambda_i^{(1)}, \lambda_i^{(2)}, \cdots, \lambda_i^{({M})}\} $ has $ {M} $ variables. For a fixed wavenumber $ k_j $ , we use cubic spline interpolation to obtain the coefficients $ c_{i_1\dots, i_{{M}}} $ from the characteristic value $ \Lambda_{D_i} $ and the shape generator $ u_{\infty}(\hat{x}_j, k_j; D_i) $ . Therefore, given a new characteristic value $ \Lambda_{D_{new}} = \{\lambda_{D_{new}}^{(1)}, \lambda_{D_{new}}^{(2)}, \cdots, \lambda_{D_{new}}^{(M)}\} $ , we obtain the predicted shape generators $ \{\hat u_{\infty}(\hat{x}_j, k_j; D_{new})\}, \, j = 1, 2, \cdots, (2N+1)^d $ .

    Finally, we specify details of recovering the geometry. As discussed above, reconstructing the geometry shape is equally to reconstructing the source function $ f $ . In the discrete formula, the domain $ V_0 $ is divided into a uniform mesh with size $ 100\times 100 $ in two dimensions and size $ 100\times 100\times 100 $ in three dimensions. Further, the approximated Fourier series $ f_{N}^{\delta} $ are computed at the mesh nodes $ \mathcal{T}_j, \, j = 1, \cdots, 100^3 $ in (5.30). Thus, the geometry shape $ D_{new} $ is approximated by the boundary of the imaging results $ f_{N}^{\delta} $ .

    In the first example, we aim to reconstruct a kite shaped domain with scale changing. The kite shaped domain is parameterized by

    $ x(t)=(β1(cost+0.65cos2t0.65), 1.5β2sint),t[0,2π], $

    where $ \beta_1 $ and $ \beta_2 $ are scale factors (characteristic values) with $ \beta_1, \beta_2 \in [0.5, \, 1.8] $ . The training dataset consists of $ 14 \times 14 $ different scale domains, i.e., $ \beta_1 $ and $ \beta_2 $ uniformly distributed on $ [0.5, \, 1.8] $ with $ M = 14^2 $ . Next, we consider four sets of different scale factors which are not covered by the training data. In this numerical experiments, the imaging results with different characteristic values are shown in Figure 1, where the black dotted lines denote the exact boundary. It is clear that the reconstructions are very closed to the exact domains.

    Figure 1.  Contour plots of reconstructed kite shape with different scale factors $ (\beta_1, \beta_2) $ : (a) $ (0.93, 1.76) $ , (b) $ (1.58, 0.87) $ , (c) $ (0.95, 0.95) $ , (d) $ (1.65, 1.65) $.

    In the second example, we aim to recover multi-domain with different scale factors. The apple shaped domain is parameterized by

    $ y(t)=β1((0.5+0.4cost+0.1sin2t)/(1+0.7cost))(cost, sint)t[0,2π], $

    and the rounded triangle shaped domain is parameterized by

    $ z(t)=β2(1+0.15cos3t)(cost, sint)t[0,2π], $

    where $ \beta_1\in [1,2] $ and $ \beta_2\in [0.5, 1.5] $ are scale factors (characteristic values) for different domains. The training dataset consists of $ 11 \times 11 $ different scale domains, that is, $ \beta_1 $ and $ \beta_2 $ uniformly distributed on $ [1,2]\times [0.5, 1] $ . Similarly, we give four sets of different scale factors which are not covered by the training data. Figure 2 shows the the reconstruction of multi-domain with different characteristic values via contour plots, where the black dotted lines denote the exact boundary. It demonstrates very good imaging performance of the approach.

    Figure 2.  Contour plots of reconstructed mixed shape with different scale factors $ (\beta_1, \beta_2) $ : (a) $ (1.13, 1.13) $ , (b) $ (1.94, 0.53) $ , (c) $ (1.88, 0.94) $ , (d) $ (1.96, 1.47) $.

    In the third example, we verify the proposed method by using a set of artificial experiments on rectangular solid. The training dataset consists of $ 125 $ rectangular solids with different height, width and length. Here the height, width and length are uniformly distributed on $ [1,2] $ with $ 5 $ amounts, i.e., $ 1, 1.25, 1.5, 1.75, 2 $ . Here, We consider four sets of different height, width and length of rectangular solids which are not covered by the training data. The imaging results with different characteristic values are shown in Figure 3, where the black dotted lines denote the shadows of the exact cube boundary. Due to discontinuities of the source, there is Gibbs phenomena on the boundary of the rectangular solids. On the whole, given the characteristic values, our proposed method is valid for determining the geometry shape.

    Figure 3.  Isosurface plots of rectangular solid with different height, width and length, where the isosurface value is $ 1 $ . The sets of height, width and length are as follows: (a) $ (1.9, 1.2, 1.4) $ , (b) $ (1.2, 1.5, 1.8) $ , (c) $ (1.6, 1.8, 1.3) $ , (d) $ (1.8, 1.8, 1.8) $.

    In last example, we consider a challenging case and verify the proposed method by using a set of synthetic experiments on 3D human body shape. The training dataset consists of $ 25 $ bodies which are generated by the MakeHuman soft [42]. MakeHuman soft is an open source software designed to create 3D virtual human models using 3D morphing technology. Providing the main parameters gender, age, muscle mass, weight, height, proportion and ethnicity, etc. of any specific character, the intermediate human shape can be reproduced with its large and own database. In this experiment, we consider two characteristic values, i.e., height and relative weight. Define the exact weight by $ EW $ and standard weight by $ SW $ , then the relative weight $ RW $ is calculated by

    $ RW=EWSW×100%. $

    Here the height of the body is given by $ 1.5\; m, 1.6\; m, 1.7\; m, 1.8\; m, 1.9\; m $ and the relative weight of the body is given by $ 60\%, 80\%, 100\%, 120\%, 140\% $ . Some human body shapes in the training dataset are presented in Figure 4. In addition, we choose two characteristic values of human body which are not covered by the training data. The first body's height is $ 1.55\; m $ and the relative weight is $ 130\% $ . The second body's height is $ 1.85\; m $ and the relative weight is $ 110\% $ . Figure 5(a) and Figure 6(a) present the exact body shape with the given characteristic value. Figure 5(b) and Figure 6(b) show the prediction of the human body shape with the given characteristic value. The results show that our method is efficient to predict the human body shape.

    Figure 4.  Isosurface plots of some training body data. Relative weight: the left column is $ 60\% $ ; the center column is $ 100\% $ , the right column is $ 140\% $ ; height: the top row is $ 1.50\; m $ , the center row is $ 1.70\; m $ , the bottom row is $ 1.90\; m $.
    Figure 5.  The relative weight is $ 130\% $ , the height is $ 1.55\; m $ . (a) Exact body, (b) reconstruction body.
    Figure 6.  The relative weight is $ 110\% $ , the height is $ 1.85\; m $ . (a) Exact body, (b) reconstruction body.

    In this paper, we develop a machine-learning method in generating a geometric body shape through prescribing a set of characteristic values of the body. The generation is mainly based on a given training dataset consisting of certain pre-selected body shapes with statistically well-sampled characteristic values. A major novelty and critical ingredient of our study is the borrowing of inverse scattering techniques in the theory of wave propagation to the geometric shape generation. We introduce the notion of shape generator which establishes a one-to-one correspondence between the geometric shape space and the function space consisting of the multiple-frequency far-field patterns associated with the time-harmonic source scattering problem. The shape generator plays an intermediate role in the geometric shape generation. First, the training dataset of geometric shapes is converted into a subset of the function space consisting of the corresponding shape generators. Then a learning model is derived through a functional interpolation of the aforementioned shape generators. For a given set of characteristic values, one then uses the learning model to obtain the shape generator of the underlying geometric body and finally reconstructs it through a multiple-frequency Fourier method.

    In the current article, the shape generator is introduced through an inverse source scattering model where we make use of the one-to-one correspondence between a geometric shape and the multiple-frequency far-field pattern associated with a compactly-supported acoustic source. One may consider to introduce the shape generator through other inverse scattering models, e.g., the inverse acoustic obstacle scattering model (cf. [28,29,31]) or the inverse electromagnetic scattering model (cf. [27,30]). In doing so, one may achieve other geometric shape generation schemes that are suitable for different applications.

    The work of H. Liu was supported by the FRG and startup grants from Hong Kong Baptist University, Hong Kong RGC General Research Funds, 12302017, 12301218 and 12302919.

    The authors declare no conflict of interest.

    [1] Gregoryanz E, Hemley R, Mao H, et al. (2000) High pressure elasticity of α-quartz: instability and ferroelastic transition.Phys Rev Lett 84: 3117-3120.
    [2] Watson G, Parker S (1995) Dynamical instabilities in α-quartz and α-berlinite: a mechanism for amorphization.Phys Rev B 52: 13306-13309.
    [3] Goel P, Mittal R, Choudhury N, et al. (2010) Lattice dynamics and Born instability in yttrium aluminum garnet, Y3Al5O12..J Phys: Condens Mat 22: 065401.
    [4] Yan X, Tang Z, Zhang L, et al. (2009) Depressurization amorphization of single-crystal boron carbide.Phys Rev Lett 102: 075505.
    [5] Born M (1940) On the stability of crystal lattices I.Math Proc Cambridge 36: 160-172.
    [6] Hill R (1975) On the elasticity and stability of perfect crystals at finite strain.Math Proc Cambridge 77: 225-240.
    [7] Morris J, Krenn C (2000) The internal stability of an elastic solid.Phil Mag A 80: 2827-2840.
    [8] Clayton J, Bliss K (2014) Analysis of intrinsic stability criteria for isotropic third-order Green elastic and compressible neo-Hookean solids.Mech Mater 68: 104-119.
    [9] Fanchini G, McCauley J, Chhowalla M (2006) Behavior of disordered boron carbide under stress.Phys Rev Lett 97: 035502.
    [10] Taylor D, Wright T, McCauley J (2011) First principles calculation of stress induced amorphization in armor ceramics. US Army Research Lab Res Rep ARL-MR-0779.Aberdeen Proving Ground, MD .
    [11] Taylor D, McCauley J, Wright T (2012) The effects of stoichiometry on the mechanical properties of icosahedral boron carbide under loading.J Phys: Condens Mat 24: 505402.
    [12] Aryal S, Rulis P, Ching W (2011) Mechanism for amorphization of boron carbide B4C under uniaxial compression.Phys Rev B 84: 184112.
    [13] Subhash G, Maiti S, Geubelle P, Ghosh D (2008) Recent advances in dynamic indentation fracture, impact damage and fragmentation of ceramics.J Am Ceram Soc 91: 2777-2791.
    [14] Chen M, McCauley J, Hemker K (2003) Shock-induced localized amorphization in boron carbide.Science 299: 1563-1566.
    [15] Grady D (2011) Adiabatic shear failure in brittle solids.Int J Impact Eng 38: 661-667.
    [16] Clayton J (2012) Towards a nonlinear elastic representation of finite compression and instability of boron carbide ceramic.Phil Mag 92: 2860-2893.
    [17] Clayton J (2013) Mesoscale modeling of dynamic compression of boron carbide polycrystals.Mech Res Commun 49: 57-64.
    [18] Koslowski M, Cuitino A, Ortiz M (2002) A phase-field theory of dislocation dynamics, strain hardening, and hysteresis in ductile single crystals.J Mech Phys Solids 50: 2597-2635.
    [19] Clayton J, Knap J (2011) A phase field model of deformation twinning: nonlinear theory and numerical simulations.Physica D 240: 841-858.
    [20] Clayton J, Knap J (2011) Phase field modeling of twinning in indentation of transparent single crystals.Model Simul Mater Sci Eng 19: 085005.
    [21] Clayton J, Knap J (2013) Phase-field analysis of fracture induced twinning in single crystals.Acta Mater 61: 5341-5353.
    [22] Kuhn C, Müller R (2012) Interpretation of parameters in phase field models for fracture.Proc Appl Math Mec 12: 161-162.
    [23] Borden M, Verhoosel C, Scott M, et al. (2012) A phase-field description of dynamic brittle fracture.Comput Method Appl Mech Eng 271: 77-95.
    [24] Voyiadjis G, Mozaffari N (2013) Nonlocal damage model using the phase field method: theory and applications.Int J Solids Structures 50: 3136-3151.
    [25] Clayton J, Knap J (2014) A geometrically nonlinear phase field theory of brittle fracture.Int J Fracture, submitted .
    [26] Clayton J (2005) Dynamic plasticity and fracture in high density polycrystals: constitutive modeling and numerical simulation.J Mech Phys Solids 53: 261-301.
    [27] Clayton J (2009) A continuum description of nonlinear elasticity, slip and twinning, with application to sapphire.Proc R Soc Lond A 465: 307-334.
    [28] Clayton J (2010) Modeling nonlinear electromechanical behavior of shocked silicon carbide.J Appl Phys 107: 013520.
    [29] Clayton J (2010) Analysis of shock compression of strong single crystals with logarithmic thermoelastic-plastic theory.Int J Eng Sci 79: 1-20.
    [30] Benallal A, Marigo J (2007) Bifurcation and stability issues in gradient theories with softening.Model Simul Mater Sci Eng 15: S283-S295.
    [31] Clayton J (2011) Nonlinear Mechanics of Crystals, Dordrecht.Springer .
    [32] Ashbee K (1971) Defects in boron carbide before and after irradiation.Acta Metall 19: 1079-1085.
    [33] Sano T, Randow C (2011) The effect of twins on the mechanical behavior of boron carbide.Metall Mater Trans A 42: 570-574.
    [34] Clayton J (2010) On anholonomic deformation, geometry, and differentiation.Math Mech Solids 17: 702-735.
    [35] Clayton J (2010) Deformation, fracture, and fragmentation in brittle geologic solids.Int J Fracture 163: 151-172.
    [36] Clayton J, McDowell D, Bammann D (2006) Modeling dislocations and disclinations with finite micropolar elastoplasticity.Int J Plasticity 22: 210-256.
    [37] Clayton J (2013) Defects in nonlinear elastic crystals: differential geometry, finite kinematics, and second-order analytical solutions.Z Angew Math Mech (ZAMM) in press .
    [38] Clayton J (2014) An alternative three-term decomposition for single crystal deformation motivated by non-linear elastic dislocation solutions.Q J Mech Appl Math 67: 127-158.
    [39] Clayton J, Bammann D, McDowell D (2005) A geometric framework for the kinematics of crystals with defects.Phil Mag 85: 3983-4010.
    [40] Clayton J (2009) A non-linear model for elastic dielectric crystals with mobile vacancies.Int J Non-Linear Mech 44: 675-688.
    [41] Clayton J (2008) A model for deformation and fragmentation in crushable brittle solids.Int J Impact Eng 35: 269-289.
    [42] Ivaneshchenko V, Shevchenko V, Turchi P (2009) First principles study of the atomic and electronic structures of crystalline and amorphous B4C.Phys Rev B 80: 235208.
    [43] Clayton J (2013) Nonlinear Eulerian thermoelasticity for anisotropic crystals.J Mech Phys Solids 61: 1983-2014.
    [44] Clayton J, McDowell D (2003) A multiscale multiplicative decomposition for elastoplasticity of polycrystals.Int J Plasticity 19: 1401-1444.
    [45] Hu S, Henager C, Chen LQ (2010) Simulations of stress-induced twinning and de-twinning: a phase field model.Acta Mater 58: 6554-6564.
    [46] Ghosh D, Subhash G, Lee C, et al. (2007) Strain-induced formation of carbon and boron clusters in boron carbide during dynamic indentation.Appl Phys Lett 91: 061910.
  • This article has been cited by:

    1. Tharifkhan Shan Ahamed, Susaimanickam Anto, Thangavel Mathimani, Kathirvel Brindhadevi, Arivalagan Pugazhendhi, Upgrading of bio-oil from thermochemical conversion of various biomass – Mechanism, challenges and opportunities, 2021, 287, 00162361, 119329, 10.1016/j.fuel.2020.119329
    2. Qun Wang, Bin Qin, Xiaohua Zhang, Xiaoling Xie, Li'e Jin, Qing Cao, Synthesis of N-doped carbon nanosheets with controllable porosity derived from bio-oil for high-performance supercapacitors, 2018, 6, 2050-7488, 19653, 10.1039/C8TA07563H
    3. Jia Wen Chong, Nishanth G. Chemmangattuvalappil, Suchithra Thangalazhy-Gopakumar, 2023, Chapter 2, 978-981-19-4846-6, 33, 10.1007/978-981-19-4847-3_2
    4. Anisha GS, Shyni Raphel M, 2023, 9780124095489, 10.1016/B978-0-323-93940-9.00057-8
    5. Zengtong Deng, Yuanjing Chen, Xun Hu, Wei Deng, Mortaza Gholizadeh, Song Hu, Sheng Su, Long Jiang, Jun Xu, Kai Xu, Yi Wang, Jun Xiang, Investigation on bio-oil pyrolysis with Ni/Al2O3 blending: Influence of the blended catalyst on coke formation, 2024, 358, 00162361, 130274, 10.1016/j.fuel.2023.130274
    6. Wasinton Simanjuntak, Kamisah Delilawati Pandiangan, Tika Dwi Febriyanti, Aryani Putri Islami, Sutopo Hadi, Ilim Ilim, Catalytic upgrading of palm oil derived bio-crude oil for bio-hydrocarbon enrichment using protonated zeolite-Y as catalyst, 2024, 12, 2333-8334, 600, 10.3934/energy.2024028
  • Reader Comments
  • © 2014 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(6036) PDF downloads(1114) Cited by(15)

Figures and Tables

Figures(3)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog