US20140025353A1 - Algorithm and a Method for Characterizing Fractal Volumes - Google Patents
Algorithm and a Method for Characterizing Fractal Volumes Download PDFInfo
- Publication number
- US20140025353A1 US20140025353A1 US13/915,684 US201313915684A US2014025353A1 US 20140025353 A1 US20140025353 A1 US 20140025353A1 US 201313915684 A US201313915684 A US 201313915684A US 2014025353 A1 US2014025353 A1 US 2014025353A1
- Authority
- US
- United States
- Prior art keywords
- dimensional
- physical domain
- software product
- cluster
- computer software
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
Images
Classifications
-
- G06F17/5018—
-
- C—CHEMISTRY; METALLURGY
- C09—DYES; PAINTS; POLISHES; NATURAL RESINS; ADHESIVES; COMPOSITIONS NOT OTHERWISE PROVIDED FOR; APPLICATIONS OF MATERIALS NOT OTHERWISE PROVIDED FOR
- C09D—COATING COMPOSITIONS, e.g. PAINTS, VARNISHES OR LACQUERS; FILLING PASTES; CHEMICAL PAINT OR INK REMOVERS; INKS; CORRECTING FLUIDS; WOODSTAINS; PASTES OR SOLIDS FOR COLOURING OR PRINTING; USE OF MATERIALS THEREFOR
- C09D163/00—Coating compositions based on epoxy resins; Coating compositions based on derivatives of epoxy resins
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01L—SEMICONDUCTOR DEVICES NOT COVERED BY CLASS H10
- H01L23/00—Details of semiconductor or other solid state devices
- H01L23/28—Encapsulations, e.g. encapsulating layers, coatings, e.g. for protection
- H01L23/29—Encapsulations, e.g. encapsulating layers, coatings, e.g. for protection characterised by the material, e.g. carbon
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01L—SEMICONDUCTOR DEVICES NOT COVERED BY CLASS H10
- H01L2924/00—Indexing scheme for arrangements or methods for connecting or disconnecting semiconductor or solid-state bodies as covered by H01L24/00
- H01L2924/0001—Technical content checked by a classifier
- H01L2924/0002—Not covered by any one of groups H01L24/00, H01L24/00 and H01L2224/00
Definitions
- the present invention is directed to a method for obtaining an analytical representation of an internal structure and spatial properties distribution of a selected physical domain, and in particular to applying an inverse algorithm to calculate the W-M fractal model and thereby determine parameters defining an analytical and continuous Weierstrass-Mandelbrot (W-M) representation.
- a lower dimension approach that is useful for determining the mechanical, thermal, electrical and fluid properties of the static and the sliding contact between two different conductors of heat and electric current involves the full identification of the parameters of the 2D Weirstrass-Mandelbrot (W-M) function, described in (Michopoulos, J. G., and Iliopoulos, A. P., 2011 . “Complete High Dimensional Inverse Characterization of Fractal Surfaces”, 2011 ASME International Design Engineering Technical Conferences & Computers and Information in Engineering Conference, Vol. CD-ROM, ASME) (“Michopoulos et al.”) and in U.S. patent application Ser. No. 13/507,968 filed Aug. 8, 2012 and incorporated herein by reference. The present invention represents a generalization of this methodology for higher dimensional fractal volumes.
- a method for obtaining an analytical representation of an internal structure and spatial properties distribution of a selected physical domain includes identifying d-dimensional correspondences of measured spatial properties or field distributions; and applying an inverse algorithm to the d-dimensional spatial properties or field distributions to calculate the W-M fractal model to thereby determine parameters defining an analytical and continuous fractal Weierstrass-Mandelbrot (W-M) representation.
- Exemplary physical domains for application of the method of the invention are generally any dense or coarse material where the method is applied to its internal structure as enclosed in its defined volume and surface such as a metal alloy, a carbon epoxy composite, a biological tissue, and any physical or virtual entity with a varying distribution of a scalar field.
- the physical domain may also be a man-made structure, like a microprocessor, an aircraft, a ship, an automobile, a building or similar technological objects.
- the physical domain alternatively may be a live organism or part of a live organism, or a digital representation like a 1 dimensional, 2 dimensional a 3 dimensional or a d-dimensional image of a real world or of an artificial object—the image can be of any shape.
- the physical domain may be a landform, like a mountain, or a lake or a part of a landform, a planet, or a part of a planet, a solar system, or a part of a solar system, a cluster of solar systems, or a part of a cluster of solar systems, a galaxy, or a part of a galaxy, a cluster of galaxies, or a part of a cluster of galaxies.
- the invention provides an efficient method for determining an analytical representation of the internal structure and spatial properties distribution of various materials such as composite, granular, and porous materials that have been loaded under various conditions as well as any physical field such as electric and/or magnetic fields, temperature, etc.
- various materials such as composite, granular, and porous materials that have been loaded under various conditions as well as any physical field such as electric and/or magnetic fields, temperature, etc.
- these materials are carbon epoxy composites for which is identified the parameters of the d-dimensional W-M fractal function from measured data.
- the invention provides a volume parameterization as a multivariate analytic generalization of its uni-variate version.
- the invention provides exact and full characterization of d-dimensional fractal volumes. Applications that require or are facilitated by having such an analytical representation of internal microstructure or spatial distribution of inhomogeneous material properties are enabled by this approach.
- the invention can be extended to the use of a homogeneous mesh for performing finite element analysis where the properties are represented by a W-M function, thereby bypassing the need for geometric tailoring of the mesh such as it explicitly represents the volumetric structure and replaces it with a generalized W-M fractal function.
- FIG. 1 is a graph showing the uniform sampling of a sphere (dots) and four typical sampling vectors;
- FIG. 2 is a a block diagram showing the forward problem perspective and the inverse problem perspective according to the invention
- FIG. 3 is a pseudo-colored isosurface combined with associated density plots of synthetic and inversely identified volumes of a W-M fractal according to the invention, demonstrating the success of the algorithm in this patent;
- FIG. 4 is a graph showing the mean error of inverse identification through singular value decomposition vs parameter ⁇ according to the invention.
- FIG. 5 compares noisy data with respective identified fractal volumes according to the invention in order to demonstrate the robustness against noise
- FIG. 6 is a graph showing the result of noise in the error of the fractal inverse identification from synthetic data according to the invention.
- FIG. 7 shows 8 of the 243 micro tomography slices according to the invention
- FIG. 8 compares original and inversely identified micro-tomographic volumes according to the invention in order to demonstrate the success of the method on real data
- FIG. 9 is a graph showing the mean absolute error of the inversion as a function of ⁇ according to the invention.
- fractal surface represents a very rich analytical representation that can model surface topographies of materials at small scales. This seems to be true especially because of the properties of continuity, non-differentiability and self affinity of specific types of fractals, that are also desired properties of surface topographies.
- the surface power spectra obeys a power-law relationship over a wide range of frequencies, because the surface topographies resemble a random process.
- Such a surface can be represented in polar coordinates (r, ⁇ ) by a complex function W as
- L is the linear size of the domain under consideration
- G is the fractal roughness
- D is the fractal dimension
- ⁇ mn is a table of (usually) random phases (with ⁇ mn ⁇ [0,2 ⁇ ])
- ⁇ is a parameter that controls the density of the frequencies
- M is the number of superimposed ridges
- N is a number that is chosen based on the desired highest sampling frequency. The following equation can be used to identify an appropriate value for N,
- L c being a cut-off wavelength, typically defined either by the highest sampling frequency, or by a physical barrier.
- the parameter ⁇ must be greater than 1 and it usually takes values in the vicinity of 1.5 because of surface flatness and frequency distribution density considerations. The latter though has been recently debated and only the requirement ⁇ >1 was considered as a valid assumption.
- n m are unit vectors spaced uniformly over a unit d-dimensional hypersphere.
- k 0 can be substituted by 2 ⁇ /L
- parameter A can be substituted by 2 ⁇ (2 ⁇ /G) 2-D and n can be bound between 0 and N.
- (3) can be written as:
- x 1 r ⁇ ⁇ cos ⁇ ⁇ ⁇ 1
- x 2 r ⁇ ⁇ sin ⁇ ⁇ ⁇ 1 ⁇ cos ⁇ ⁇ ⁇ 2
- x 3 r ⁇ ⁇ sin ⁇ ⁇ ⁇ 1 ⁇ sin ⁇ ⁇ ⁇ 2 ⁇ cos ⁇ ⁇ ⁇ 3
- x d r ⁇ ⁇ sin ⁇ ⁇ ⁇ 1 ⁇ sin ⁇ ⁇ ⁇ 2 ⁇ ⁇ ... ⁇ ⁇ sin ⁇ ⁇ ⁇ d - 2 ⁇ sin ⁇ ⁇ ⁇ d - 1 ⁇ , ( 5 )
- ⁇ i are the angular coordinates.
- a uniform sampling of the hypersphere can be achieved by setting,
- FIG. 2 The systemic view of both the forward and the inverse problems are shown in FIG. 2 .
- the inverse problem it is shown that its solution produces the parameters of the fractal function provided that field distribution of the volume z(r) has been experimentally determined and therefore is known.
- the system of Eq. 17 is an overdetermined system of M(N+1) equations. Since the right hand side vector z contains experimental measurements, it also contains noise; the system cannot, in general, have an exact solution. Nevertheless, one can seek a p, such as ⁇ p-z ⁇ is minimized, where ⁇ p ⁇ . is the vector norm. Such a p is known as the least squares solution to the over-determined system. It should be noted that the left hand side expression of Eq. 16 yields results in the complex domain, but as long as a minimal solution is achieved for real numbers on the right hand side, the imaginary parts will be close to 0. A solution can be given by the following equation,
- V is calculated by the Singular Value Decomposition (SVD) of C as,
- d i is the i th entry of the diagonal of D.
- z i e is the value of the scalar field of the experimentally measured (or reference) points
- P is the number of those points.
- the first experiment was based on synthetic data and is aimed at inversely identifying only the phases of a volume constructed by the fractal itself.
- the original volume is shown in FIG. 3 and was constructed using Eq. (4) with random phases and the parameters as shown in Table 1.
- the inversely identified volume using the phases resulting from Eq. (19) is presented in FIG. 3 . It should be noted that this difference is very small compared to the magnitude of the volume and any discrepancies should be considered as the numerical error of the SVD algorithm.
- the second synthetic experiment involved the identification of both the phases and the ⁇ parameter.
- An exhaustive search approach was adopted in this case, as the sensitivity of the SVD inversion relative to the value of ⁇ is also of interest.
- the inversion of the phases was executed and the value of the error function (Eq. 22) was calculated.
- the error for various values of ⁇ is presented on FIG. 4 .
- AS4/3501-6 is a carbon/epoxy composite used in general purpose structural applications and the associated prepreg is manufactured by Hexcel Corporation.
- the data set used consisted of 243 tomographic slices of the specimen stored as digital images of resolution 1256 ⁇ 588 pixels each.
- the volume characterization was done on a subset of the volumetric data shown in FIG. 7 .
- the subset size was 16 ⁇ 16 ⁇ 16 voxels and was taken from the region depicted by the red square in FIG. 8 .
- the inversion of the volume gives very satisfactory results as shown in FIG. 8 .
- the mean absolute error as calculated by Eq. (22) was 4.74%.
- Non-volatile computer readable media capable of storing and/or transferring computer programs or computer-readable instructions for execution by one or more computers.
- Volatile computer readable media that can be used can include a compact disk, hard disk, floppy disk, tape, magneto-optical disk, PROM (EPROM, EEPROM, flash EPROM), DRAM, SRAM, SDRAM, or any other magnetic medium; punch card, paper tape, or any other physical medium.
- Non-volatile media can include a memory such as a dynamic memory in a computer.
- computer readable media that can be used to store and/or transmit instructions for carrying out methods described herein can include non-physical media such as an electromagnetic carrier wave, acoustic wave, or light wave such as those generated during radio wave and infrared data communications.
Abstract
A computer implemented method for obtaining an analytical representation of an internal structure and spatial properties distribution of a selected physical domain includes identifying d-dimensional correspondences of measured spatial properties or field distributions; and applying an inverse algorithm to the d-dimensional spatial properties or field distributions to calculate the Weierstrass-Mandelbrot (W-M) fractal model to thereby determine parameters defining an analytical and continuous Weierstrass-Mandelbrot (W-M) representation.
Description
- This application claims the benefit of U.S. Provisional Application 61/671,378 filed on Jul. 13, 2012 and incorporated herein by reference.
- The present invention is directed to a method for obtaining an analytical representation of an internal structure and spatial properties distribution of a selected physical domain, and in particular to applying an inverse algorithm to calculate the W-M fractal model and thereby determine parameters defining an analytical and continuous Weierstrass-Mandelbrot (W-M) representation.
- The need to adequately represent and model the structure and the behavior of physical or virtual volumes containing physical property distributions or physical scalar fields as they are relevant to various applications has motivated the utilization of appropriate fractal functions that can be used to represent such volumes.
- A lower dimension approach that is useful for determining the mechanical, thermal, electrical and fluid properties of the static and the sliding contact between two different conductors of heat and electric current involves the full identification of the parameters of the 2D Weirstrass-Mandelbrot (W-M) function, described in (Michopoulos, J. G., and Iliopoulos, A. P., 2011. “Complete High Dimensional Inverse Characterization of Fractal Surfaces”, 2011 ASME International Design Engineering Technical Conferences & Computers and Information in Engineering Conference, Vol. CD-ROM, ASME) (“Michopoulos et al.”) and in U.S. patent application Ser. No. 13/507,968 filed Aug. 8, 2012 and incorporated herein by reference. The present invention represents a generalization of this methodology for higher dimensional fractal volumes.
- According to the invention, a method for obtaining an analytical representation of an internal structure and spatial properties distribution of a selected physical domain includes identifying d-dimensional correspondences of measured spatial properties or field distributions; and applying an inverse algorithm to the d-dimensional spatial properties or field distributions to calculate the W-M fractal model to thereby determine parameters defining an analytical and continuous fractal Weierstrass-Mandelbrot (W-M) representation.
- Exemplary physical domains for application of the method of the invention are generally any dense or coarse material where the method is applied to its internal structure as enclosed in its defined volume and surface such as a metal alloy, a carbon epoxy composite, a biological tissue, and any physical or virtual entity with a varying distribution of a scalar field. The physical domain may also be a man-made structure, like a microprocessor, an aircraft, a ship, an automobile, a building or similar technological objects. The physical domain alternatively may be a live organism or part of a live organism, or a digital representation like a 1 dimensional, 2 dimensional a 3 dimensional or a d-dimensional image of a real world or of an artificial object—the image can be of any shape. The physical domain may be a landform, like a mountain, or a lake or a part of a landform, a planet, or a part of a planet, a solar system, or a part of a solar system, a cluster of solar systems, or a part of a cluster of solar systems, a galaxy, or a part of a galaxy, a cluster of galaxies, or a part of a cluster of galaxies.
- The invention provides an efficient method for determining an analytical representation of the internal structure and spatial properties distribution of various materials such as composite, granular, and porous materials that have been loaded under various conditions as well as any physical field such as electric and/or magnetic fields, temperature, etc. In particular, an example of these materials are carbon epoxy composites for which is identified the parameters of the d-dimensional W-M fractal function from measured data. The invention provides a volume parameterization as a multivariate analytic generalization of its uni-variate version.
- The invention provides exact and full characterization of d-dimensional fractal volumes. Applications that require or are facilitated by having such an analytical representation of internal microstructure or spatial distribution of inhomogeneous material properties are enabled by this approach. The invention can be extended to the use of a homogeneous mesh for performing finite element analysis where the properties are represented by a W-M function, thereby bypassing the need for geometric tailoring of the mesh such as it explicitly represents the volumetric structure and replaces it with a generalized W-M fractal function.
-
FIG. 1 is a graph showing the uniform sampling of a sphere (dots) and four typical sampling vectors; -
FIG. 2 is a a block diagram showing the forward problem perspective and the inverse problem perspective according to the invention; -
FIG. 3 is a pseudo-colored isosurface combined with associated density plots of synthetic and inversely identified volumes of a W-M fractal according to the invention, demonstrating the success of the algorithm in this patent; -
FIG. 4 is a graph showing the mean error of inverse identification through singular value decomposition vs parameter γ according to the invention; -
FIG. 5 compares noisy data with respective identified fractal volumes according to the invention in order to demonstrate the robustness against noise; -
FIG. 6 is a graph showing the result of noise in the error of the fractal inverse identification from synthetic data according to the invention; -
FIG. 7 shows 8 of the 243 micro tomography slices according to the invention; -
FIG. 8 compares original and inversely identified micro-tomographic volumes according to the invention in order to demonstrate the success of the method on real data; and -
FIG. 9 is a graph showing the mean absolute error of the inversion as a function of γ according to the invention. - d-Dimensional Modeling
- It has been established that a fractal surface represents a very rich analytical representation that can model surface topographies of materials at small scales. This seems to be true especially because of the properties of continuity, non-differentiability and self affinity of specific types of fractals, that are also desired properties of surface topographies. The surface power spectra obeys a power-law relationship over a wide range of frequencies, because the surface topographies resemble a random process. Such a surface can be represented in polar coordinates (r, θ) by a complex function W as
-
- where L is the linear size of the domain under consideration, G is the fractal roughness, D is the fractal dimension, φmn is a table of (usually) random phases (with φmnε[0,2π]), γ is a parameter that controls the density of the frequencies, M is the number of superimposed ridges and N is a number that is chosen based on the desired highest sampling frequency. The following equation can be used to identify an appropriate value for N,
-
N=n max =int[log(L/Lc)/log γ], (2) - with Lc being a cut-off wavelength, typically defined either by the highest sampling frequency, or by a physical barrier.
- It should be noted that the parameter γ must be greater than 1 and it usually takes values in the vicinity of 1.5 because of surface flatness and frequency distribution density considerations. The latter though has been recently debated and only the requirement γ>1 was considered as a valid assumption.
- Although the above mentioned representation is very useful in various physical contexts, its generalization to higher dimensions can provide a convenient representation of fields in those dimensions. Eq. (1) can be generalized to d-dimensions by the formula
-
- where nm are unit vectors spaced uniformly over a unit d-dimensional hypersphere. k0 can be substituted by 2 π/L, parameter A can be substituted by 2π(2 π/G)2-D and n can be bound between 0 and N. Then (3) can be written as:
-
- The construction of the unit vectors nm can be achieved with the help of the hyperspherical coordinates
-
- where θi are the angular coordinates. For unit vectors, one can use r=1. It should be noted that in order to avoid traversing the hypersphere twice we can choose to traverse θ1 . . . θd-2 over the range [0,π]. A uniform sampling of the hypersphere can be achieved by setting,
-
- where pεN−{0} defines the density of the winding around the surface of the hypersphere. For p=1 the hypershpere will be sampled with only a single traversal, while for higher numbers of p, denser windings can be achieved. Using (5) and (6) the unit vector sampling can be written as:
-
- An example sampling of a 3-sphere is shown in
FIG. 1 for p=8 and M=64. - The Inverse Problem
- The systemic view of both the forward and the inverse problems are shown in
FIG. 2 . In the case of the inverse problem, it is shown that its solution produces the parameters of the fractal function provided that field distribution of the volume z(r) has been experimentally determined and therefore is known. - More specifically, given an array of measurements zijk e, i=1 . . . I, j=1 . . . J, k=1 . . . K over a region of size Ld one may then identify the parameters γ and φmn of a d-dimensional volume z(r) that best fits those measurements. In previous studies the identification of those parameters was done only for the 2 Dimensional surface. In order to achieve the inverse identification of the volume parameters the approach is similar to that described in Michopoulos et al.
- The inner sum in Eq. 4 can be written as,
-
- where cm(r) and φm are vectors in P(N+1) given by,
-
- By setting,
-
- Eq. 4 can be written as:
-
- By coalescing the vectors cm and φm into larger vectors in PM(N+1) such as,
-
c(r)={c 1 T ,c 2 T , . . . ,c M T}T, (13) -
and -
φ={φ1 T,φ2 T, . . . ,φM T}T, (14) - Eq. 12 can be written as:
-
W(r)=*Qc(r)Tφ. (15) - For the needs of the inverse problem it is assumed that a number of measurements at points rt e exist for a volume represented as ze(rt e), t=1 . . . T, T≧M(N+1). We seek to identify a volume that is described by Eq. 15 and approximates the experimental points ze(rt e).
- To solve this problem we first form the following linear system,
-
- By expanding the vectors in Eq. 16, the system can be written as,
-
- The system of Eq. 17 is an overdetermined system of M(N+1) equations. Since the right hand side vector z contains experimental measurements, it also contains noise; the system cannot, in general, have an exact solution. Nevertheless, one can seek a p, such as Πp-zΠ is minimized, where ΠpΓ. is the vector norm. Such a p is known as the least squares solution to the over-determined system. It should be noted that the left hand side expression of Eq. 16 yields results in the complex domain, but as long as a minimal solution is achieved for real numbers on the right hand side, the imaginary parts will be close to 0. A solution can be given by the following equation,
-
p=Vy, (19) - where V is calculated by the Singular Value Decomposition (SVD) of C as,
-
UDV T =C, (20) - where y is a vector defined as yt=bt′/dt, b=bt is a vector given by,
-
b=U T p), (21) - and di is the i th entry of the diagonal of D.
- The solution of the inverse problem as described by the overdetermined system of Eq. 16 gives the phases φmn given known values of the other parameters. In a general volume the only other parameter that is unknown is γ. As was pointed out in Michopoulos et al parameters G and D don't need to be considered as unknowns to be determined in this optimization. This is because for any combination of the phases it is always possible to find new values for φmn, that result in generating the same surface.
- Numerical Experiments
- In order to assess the quality of the high dimensional fractal characterization results of the numerical examples that follow, we define the following error function,
-
- where zi e is the value of the scalar field of the experimentally measured (or reference) points, P is the number of those points. In the following examples P is set as P=I×J×K=40×40×40=64000 points. zi d is the value of the scalar field of the inversely identified fractal and is equal to the real part of the truncated W-M function of Eq. (4) (zi d=Re{W(r)}).
-
TABLE 1 PARAMETERS OF THE REFERENCE MODEL. Parameter Value L 1 G 1 × 106 γ 1.5 D 2.3 M 10 N 16 d 3 - To study the feasibility of the invention, a few numerical experiments were designed in the 3D space. The first experiment was based on synthetic data and is aimed at inversely identifying only the phases of a volume constructed by the fractal itself. The original volume is shown in
FIG. 3 and was constructed using Eq. (4) with random phases and the parameters as shown in Table 1. The inversely identified volume using the phases resulting from Eq. (19) is presented inFIG. 3 . It should be noted that this difference is very small compared to the magnitude of the volume and any discrepancies should be considered as the numerical error of the SVD algorithm. - The second synthetic experiment involved the identification of both the phases and the γ parameter. An exhaustive search approach was adopted in this case, as the sensitivity of the SVD inversion relative to the value of γ is also of interest. For a range of the possible values for parameter γ the inversion of the phases was executed and the value of the error function (Eq. 22) was calculated. The error for various values of γ is presented on
FIG. 4 . The smallest value for the error is at γ=1.5, which is the one used originally for the generation of the volume (Table 1). - In order to investigate the behavior of the invention under noisy data, we also applied the inversion in the original data of
FIG. 3 , but for various levels of injected noise. The noise levels varied relative to the range of the scalar field values from 0% to 300%. The results of the inverse identification fornoise levels 20%, 40% and 300% are presented inFIG. 5 and show impressively robust response. Even for noise levels of 300% (FIG. 5 ), where the volume characteristics are visually non-existent, the fractal inversion methodology resembles (FIG. 5 ) the actual field very well. InFIG. 6 the mean absolute error of the identification for various levels of noise is presented, where it can be observed that as the noise levels increase, the mean absolute error increases linearly. - Although the previous analysis demonstrates the consistency and performance of the invention, it is much more useful when applied to actual data. For this reason, a numerical test is performed based on micro-tomography scans of a failed AS4/3501-6 composite specimen used for 3 degrees of freedom in-plane mechanical loading. AS4/3501-6 is a carbon/epoxy composite used in general purpose structural applications and the associated prepreg is manufactured by Hexcel Corporation. The data set used consisted of 243 tomographic slices of the specimen stored as digital images of resolution 1256×588 pixels each.
- Because the data set is very large, the volume characterization was done on a subset of the volumetric data shown in
FIG. 7 . The subset size was 16×16×16 voxels and was taken from the region depicted by the red square inFIG. 8 . The fractal parameters were chosen to be: γ=1.5,D=2.3 μL=1,G=1*106 μM=80 and N=50. The inversion of the volume gives very satisfactory results as shown inFIG. 8 . For this test the mean absolute error as calculated by Eq. (22) was 4.74%. - To investigate the effect the choice of γ has on the fractal inversion, a parametric study was performed for various values of γ ranging from 1.02 to 2. For each of those values the fractal phases where identified and the values of the mean absolute error from Eq. (22) were calculated. The results of this study are plotted in
FIG. 9 . Although the original choice of γ in the previous example resulted in a very satisfactory inverse identification of the failed composite specimen volume, it can be seen fromFIG. 9 , that even better results can be achieved for a value of γ equal to 1.24. In that particular case the mean absolute error of the identified volume relative to the original was about 0.84% and the resulting volume is depicted inFIG. 7 . This result underlines that for optimum characterization purposes on has to consider that the set of parameters to be identified should include both the phases and γ. - It should be noted that the present invention can be accomplished by executing one or more sequences of one or more computer-readable instructions read into a memory of one or more computers from volatile or non-volatile computer-readable media capable of storing and/or transferring computer programs or computer-readable instructions for execution by one or more computers. Volatile computer readable media that can be used can include a compact disk, hard disk, floppy disk, tape, magneto-optical disk, PROM (EPROM, EEPROM, flash EPROM), DRAM, SRAM, SDRAM, or any other magnetic medium; punch card, paper tape, or any other physical medium. Non-volatile media can include a memory such as a dynamic memory in a computer. In addition, computer readable media that can be used to store and/or transmit instructions for carrying out methods described herein can include non-physical media such as an electromagnetic carrier wave, acoustic wave, or light wave such as those generated during radio wave and infrared data communications.
- While specific embodiments of the present invention have been shown and described, it should be understood that other modifications, substitutions and alternatives are apparent to one of ordinary skill in the art. Such modifications, substitutions and alternatives can be made without departing from the spirit and scope of the invention, which should be determined from the appended claims.
Claims (18)
1. A computer implemented method for obtaining an analytical representation of an internal structure and spatial properties distribution of a selected physical domain, comprising:
identifying d-dimensional correspondences of measured spatial properties or field distributions; and applying an inverse algorithm to the d-dimensional spatial properties or field distributions to calculate the Weierstrass-Mandelbrot (W-M) fractal model to thereby determine parameters defining an analytical and continuous Weierstrass-Mandelbrot (W-M) representation.
2. The method of claim 1 , wherein the physical domain is selected from a metal alloy, a carbon epoxy composite, or a biological tissue.
3. The method of claim 1 , wherein the physical domain is a man-made structure selected from a microprocessor, an aircraft, a ship, an automobile, a darn, a road layer, a bridge, or a building.
4. The method of claim 1 , wherein the physical domain is a live organism or part of a live organism or a in vitro version of them (e.g wood).
5. The method of claim 1 , wherein the physical domain is a digital representation of a real or artificial object.
6. The method of claim 5 , wherein the digital representation is a spatial distribution of a scalar field within a volume such as magnetic resonance slice images, x-ray tomography data, ultrasonic data trough a volume, or any simulated entity distributed in a volume.
7. The method of claim 1 , wherein the physical domain is a landform selected from a mountain, a lake, part of a landform, a planet, part of a planet, a solar system, part of a solar system, a cluster of solar systems, part of a cluster of solar systems, a galaxy, part of a galaxy, a cluster of galaxies, or part of a cluster of galaxies.
8. The method of claim 1 , wherein the d-dimensional modeling is obtained by applying a formula
where nm are d-dimensional unit.
9. The method of claim 8 , wherein based on an array of measurements over a selected region of the domain, parameters of the d-dimensional W-M fractal model are identified that best fit the array of measurements by decomposing the problem so as to apply singular value decomposition for determining unknown phases of the W-M fractal model.
10. A computer software product comprising a physical computer-readable medium including stored instructions that, when executed by a computer, cause the computer to:
identify d-dimensional correspondences of measured spatial properties or field distributions; and
apply an inverse algorithm to the d-dimensional spatial properties or field distributions to calculate the Weierstrass-Mandelbrot (W-M) fractal model to thereby determine parameters defining an analytical and continuous Weierstrass-Mandelbrot (W-M) representation.
11. The computer software product of claim 10 , wherein the physical domain is selected from a metal alloy, a carbon epoxy composite, or a biological tissue.
12. The computer software product of claim 10 , wherein the physical domain is a man-made structure selected from a microprocessor, an aircraft, a ship, an automobile, a dam, a road layer, a bridge, or a building.
13. The computer software product of claim 10 , wherein the physical domain is a live organism or part of a live organism or a in vitro version of them (e.g wood).
14. The computer software product of claim 10 , wherein the physical domain is a digital representation of a real or artificial object.
15. The computer software product of claim 14 , wherein the digital representation is a spatial distribution of a scalar field within a volume such as magnetic resonance slice images, x-ray tomography data, ultrasonic data trough a volume, or any simulated entity distributed in a volume.
16. The computer software product of claim 10 , wherein the physical domain is a landform selected from a mountain, a lake, part of a landform, a planet, part of a planet, a solar system, part of a solar system, a cluster of solar systems, part of a cluster of solar systems, a galaxy, part of a galaxy, a cluster of galaxies, or part of a cluster of galaxies.
17. The computer software product of claim 10 , wherein the d-dimensional modeling is obtained by applying a formula
where nm are d-dimensional unit.
18. The computer software product of claim 17 , wherein based on an array of measurements over a selected region of the domain, parameters of the d-dimensional W-M fractal model are identified that best fit the array of measurements by decomposing the problem so as to apply singular value decomposition for determining unknown phases of the W-M fractal model.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/915,684 US20140025353A1 (en) | 2012-07-13 | 2013-06-12 | Algorithm and a Method for Characterizing Fractal Volumes |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201261671378P | 2012-07-13 | 2012-07-13 | |
US13/915,684 US20140025353A1 (en) | 2012-07-13 | 2013-06-12 | Algorithm and a Method for Characterizing Fractal Volumes |
Publications (1)
Publication Number | Publication Date |
---|---|
US20140025353A1 true US20140025353A1 (en) | 2014-01-23 |
Family
ID=49947277
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/723,336 Active US8847415B1 (en) | 2012-07-13 | 2012-12-21 | Liquid compression molding encapsulants |
US13/915,684 Abandoned US20140025353A1 (en) | 2012-07-13 | 2013-06-12 | Algorithm and a Method for Characterizing Fractal Volumes |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/723,336 Active US8847415B1 (en) | 2012-07-13 | 2012-12-21 | Liquid compression molding encapsulants |
Country Status (1)
Country | Link |
---|---|
US (2) | US8847415B1 (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104091033A (en) * | 2014-07-25 | 2014-10-08 | 哈尔滨工业大学 | Bridge static force finite element model correcting method based on super-element combination virtual deformation method |
CN106354975A (en) * | 2016-09-23 | 2017-01-25 | 清华大学 | Finite element method for acquiring misalignment quantity of planetary gear |
CN107273651A (en) * | 2017-08-23 | 2017-10-20 | 中铁二院工程集团有限责任公司 | A kind of super-high pier construction design method |
CN107292100A (en) * | 2017-06-15 | 2017-10-24 | 内蒙古科技大学 | A kind of computational methods of rock structural face fractal dimension |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120029888A1 (en) * | 2010-07-30 | 2012-02-02 | Jeng-Haur Horng | Topography shaping apparatus for forming surfaces of low friction coefficient |
US8208698B2 (en) * | 2007-12-14 | 2012-06-26 | Mela Sciences, Inc. | Characterizing a texture of an image |
US8884954B2 (en) * | 2011-08-25 | 2014-11-11 | The United States Of America, As Represented By The Secretary Of The Navy | Algorithm and a method for characterizing surfaces with fractal nature |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1305702A (en) | 1969-02-20 | 1973-02-07 | ||
US4528366A (en) | 1982-09-28 | 1985-07-09 | The Dow Chemical Company | Production of polytriazines from aromatic polycyanates with cobalt salt of a carboxylic acid as catalyst |
US4477629A (en) | 1983-07-27 | 1984-10-16 | The Dow Chemical Company | Cyanate-containing polymers |
US4751323A (en) | 1983-11-16 | 1988-06-14 | The Dow Chemical Company | Novel polyaromatic cyanates |
JP5232645B2 (en) * | 2006-06-26 | 2013-07-10 | パナソニック株式会社 | Thermosetting resin composition and circuit board mounting method and repair process using the same |
EP2185621A1 (en) * | 2007-09-27 | 2010-05-19 | Nippon Shokubai Co., Ltd. | Curable resin composition for molded bodies, molded body, and production method thereof |
JP5475223B2 (en) * | 2007-10-09 | 2014-04-16 | 株式会社Adeka | One-component cyanate-epoxy composite resin composition, cured product thereof and method for producing the same, and sealing material and adhesive using the same |
US20120123053A1 (en) * | 2010-11-16 | 2012-05-17 | Makoto Kibayashi | Carbon fiber |
JP5838646B2 (en) * | 2011-08-10 | 2016-01-06 | セイコーエプソン株式会社 | Projection type image device and method of manufacturing optical element of projection type image device |
-
2012
- 2012-12-21 US US13/723,336 patent/US8847415B1/en active Active
-
2013
- 2013-06-12 US US13/915,684 patent/US20140025353A1/en not_active Abandoned
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8208698B2 (en) * | 2007-12-14 | 2012-06-26 | Mela Sciences, Inc. | Characterizing a texture of an image |
US20120029888A1 (en) * | 2010-07-30 | 2012-02-02 | Jeng-Haur Horng | Topography shaping apparatus for forming surfaces of low friction coefficient |
US8884954B2 (en) * | 2011-08-25 | 2014-11-11 | The United States Of America, As Represented By The Secretary Of The Navy | Algorithm and a method for characterizing surfaces with fractal nature |
Non-Patent Citations (5)
Title |
---|
GESchwarz "Re:6061 vs 20214" HomebuiltAirplanes.com (2008) <http://www.homebuiltairplanes.com/forums/aircraft-design-aerodynamics-new-technology/5003-6061-vs-2024-a.html>. * |
Michopoulos, .T. G., & Iliopoulos. A., "Four Parameter Inverse Characterization of Fractal Surfaces" In Proceedings of the ASME Int'l Design Engineering Technical Conf. and Computers and Information in Engineering Conf., IDETC/CIE 2010, (2010). * |
Michopoulos, John G. & Iliopoulos, Athanasios "Complete High Dimensional Inverse Characterization of Fractal Surfaces" Proceedings of ASME 2011 Int'l Design Engineering Technical Conf. & Computers and Information in Engineering Conf., pp. 447-456 (2011). * |
Wang, Shao, et al. "Determination of the Fractal Scaling Parameter From Simulated Fractal-Regular Surface Profiles Based on the Weierstrass-Mandelbrot Function" J. Tribology, vol. 129, no. 4, pp. 952-956 (2007). * |
Yan, W. & Komvopoulos, K. "Contact Analysis of Elastic-Plastic Fractal Surfaces" J. Appl. Phys., vol. 84, no. 7 (1998). * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104091033A (en) * | 2014-07-25 | 2014-10-08 | 哈尔滨工业大学 | Bridge static force finite element model correcting method based on super-element combination virtual deformation method |
CN106354975A (en) * | 2016-09-23 | 2017-01-25 | 清华大学 | Finite element method for acquiring misalignment quantity of planetary gear |
CN107292100A (en) * | 2017-06-15 | 2017-10-24 | 内蒙古科技大学 | A kind of computational methods of rock structural face fractal dimension |
CN107273651A (en) * | 2017-08-23 | 2017-10-20 | 中铁二院工程集团有限责任公司 | A kind of super-high pier construction design method |
Also Published As
Publication number | Publication date |
---|---|
US8847415B1 (en) | 2014-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Donnelly et al. | Tomographic reconstruction of a three-dimensional magnetization vector field | |
Strobelt et al. | Rolled‐out wordles: A heuristic method for overlap removal of 2d data representatives | |
Lang et al. | On the analysis of spatial binary images | |
US8884954B2 (en) | Algorithm and a method for characterizing surfaces with fractal nature | |
Bentel et al. | Different radial basis functions and their applicability for regional gravity field representation on the sphere | |
Javaran et al. | A dual reciprocity BEM approach using new Fourier radial basis functions applied to 2D elastodynamic transient analysis | |
US20140025353A1 (en) | Algorithm and a Method for Characterizing Fractal Volumes | |
Cakoni et al. | Qualitative methods in inverse electromagnetic scattering theory: Inverse scattering for anisotropic media | |
Griesmaier et al. | A multifrequency MUSIC algorithm for locating small inhomogeneities in inverse scattering | |
Li et al. | A three-dimensional model-based inversion algorithm using radial basis functions for microwave data | |
Bondarenko et al. | The factorization method for inverse acoustic scattering in a layered medium | |
He et al. | Microwave imaging of 3D dielectric-magnetic penetrable objects based on integral equation method | |
Saunders et al. | Constrained resistivity inversion using seismic data | |
Barman et al. | A three‐dimensional statistical model for imaged microstructures of porous polymer films | |
Vecherin et al. | Incorporating source directionality into outdoor sound propagation calculations | |
Kirisits et al. | Fourier reconstruction for diffraction tomography of an object rotated into arbitrary orientations | |
Hu et al. | A synergy of the velocity gradients technique and the probability density functions for identifying gravitational collapse in self-absorbing media | |
Kunyansky | Inversion of the spherical means transform in corner-like domains by reduction to the classical Radon transform | |
Jia et al. | 3-D model-based inversion using supervised descent method for aspect-limited microwave data of metallic targets | |
Chen et al. | A direct imaging method for inverse electromagnetic scattering problem in rectangular waveguide | |
Yesylevskyy et al. | Determination of the shape and curvature of nonplanar lipid bilayers that are bent in a single plane in molecular dynamics simulations | |
JP6784984B2 (en) | Measuring device and measuring method | |
Liu et al. | Overlapping Yee FDTD method on nonorthogonal grids | |
Michopoulos et al. | Complete high dimensional inverse characterization of fractal surfaces and volumes | |
Chase | Michopoulos et al.(43) Pub. Date: Jan. 23, 2014 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: THE GOVERNMENT OF THE UNITED STATES OF AMERICA, AS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:MICHOPOULOS, JOHN G;REEL/FRAME:031009/0553 Effective date: 20130206 |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |