BSP-fields: An Exact Representation of Polygonal Objects by Differentiable Scalar Fields Based on Binary Space Partitioning

The problem considered in this work is to find a dimension independent algorithm for the generation of signed scalar fields exactly representing polygonal objects and satisfying the following requirements: the defining real function takes zero value exactly at the polygonal object boundary; no extra zero-value isosurfaces should be generated; C1 continuity of the function in the entire domain. The proposed algorithms are based on the binary space partitioning (BSP) of the object by the planes passing through the polygonal faces and are independent of the object genus, the number of disjoint components, and holes in the initial polygonal mesh. Several extensions to the basic algorithm are proposed to satisfy the selected optimization criteria. The generated BSP-fields allow for applying techniques of the function-based modeling to already existing legacy objects from CAD and computer animation areas, which is illustrated by several examples

A multiscale mechanobiological model of bone remodelling predicts site-specific bone loss in the femur during osteoporosis and mechanical disuse

We propose a multiscale mechanobiological model of bone remodelling to investigate the site-specific evolution of bone volume fraction across the midshaft of a femur. The model includes hormonal regulation and biochemical coupling of bone cell populations, the influence of the microstructure on bone turnover rate, and mechanical adaptation of the tissue. Both microscopic and tissue-scale stress/strain states of the tissue are calculated from macroscopic loads by a combination of beam theory and micromechanical homogenisation. This model is applied to simulate the spatio-temporal evolution of a human midshaft femur scan subjected to two deregulating circumstances: (i) osteoporosis and (ii) mechanical disuse. Both simulated deregulations led to endocortical bone loss, cortical wall thinning and expansion of the medullary cavity, in accordance with experimental findings. Our model suggests that these observations are attributable to a large extent to the influence of the microstructure on bone turnover rate. Mechanical adaptation is found to help preserve intracortical bone matrix near the periosteum. Moreover, it leads to non-uniform cortical wall thickness due to the asymmetry of macroscopic loads introduced by the bending moment. The effect of mechanical adaptation near the endosteum can be greatly affected by whether the mechanical stimulus includes stress concentration effects or not.Comment: 25 pages, 10 figure

Procedural function-based modelling of volumetric microstructures

We propose a new approach to modelling heterogeneous objects containing internal volumetric structures with size of details orders of magnitude smaller than the overall size of the object. The proposed function-based procedural representation provides compact, precise, and arbitrarily parameterised models of coherent microstructures, which can undergo blending, deformations, and other geometric operations, and can be directly rendered and fabricated without generating any auxiliary representations (such as polygonal meshes and voxel arrays). In particular, modelling of regular lattices and cellular microstructures as well as irregular porous media is discussed and illustrated. We also present a method to estimate parameters of the given model by fitting it to microstructure data obtained with magnetic resonance imaging and other measurements of natural and artificial objects. Examples of rendering and digital fabrication of microstructure models are presented

TauFactor: An open-source application for calculating tortuosity factors from tomographic data

TauFactor is a MatLab application for efficiently calculating the tortuosity factor, as well as volume fractions, surface areas and triple phase boundary densities, from image based microstructural data. The tortuosity factor quantifies the apparent decrease in diffusive transport resulting from convolutions of the flow paths through porous media. TauFactor was originally developed to improve the understanding of electrode microstructures for batteries and fuel cells; however, the tortuosity factor has been of interest to a wide range of disciplines for over a century, including geoscience, biology and optics. It is still common practice to use correlations, such as that developed by Bruggeman, to approximate the tortuosity factor, but in recent years the increasing availability of 3D imaging techniques has spurred interest in calculating this quantity more directly. This tool provides a fast and accurate computational platform applicable to the big datasets (>10^8 voxels) typical of modern tomography, without requiring high computational power

Darcian permeability constant as indicator for shear stresses in regular scaffold systems for tissue engineering

The shear stresses in printed scaffold systems for tissue engineering depend on the flow properties and void volume in the scaffold. In this work, computational fluid dynamics (CFD) is used to simulate flow fields within porous scaffolds used for cell growth. From these models the shear stresses acting on the scaffold fibres are calculated. The results led to the conclusion that the Darcian (k 1) permeability constant is a good predictor for the shear stresses in scaffold systems for tissue engineering. This permeability constant is easy to calculate from the distance between and thickness of the fibres used in a 3D printed scaffold. As a consequence computational effort and specialists for CFD can be circumvented by using this permeability constant to predict the shear stresses. If the permeability constant is below a critical value, cell growth within the specific scaffold design may cause a significant increase in shear stress. Such a design should therefore be avoided when the shear stress experienced by the cells should remain in the same order of magnitud

New Method for Estimating Fractal Dimension in 3D Space and Its Application to Complex Surfaces

The concept of “surface modeling” generally describes the process of representing a physical or artificial surface by a geometric model, namely a mathematical expression. Among the existing techniques applied for the characterization of a surface, terrain modeling relates to the representation of the natural surface of the Earth. Cartographic terrain or relief models as three-dimensional representations of a part of the Earth's surface convey an immediate and direct impression of a landscape and are much easier to understand than two-dimensional models. This paper addresses a major problem in complex surface modeling and evaluation consisting in the characterization of their topography and comparison among different textures, which can be relevant in different areas of research. A new algorithm is presented that allows calculating the fractal dimension of images of complex surfaces. The method is used to characterize different surfaces and compare their characteristics. The proposed new mathematical method computes the fractal dimension of the 3D space with the average space component of Hurst exponent H, while the estimated fractal dimension is used to evaluate, compare and characterize complex surfaces that are relevant in different areas of research. Various surfaces with both methods were analyzed and the results were compared. The study confirms that with known coordinates of a surface, it is possible to describe its complex structure. The estimated fractal dimension is proved to be an ideal tool for measuring the complexity of the various surfaces considered

New method for estimating fractal dimension in 3d space and its application to complex surfaces

The concept of “surface modeling” generally describes the process of representing a physical or artificial surface by a geometric model, namely a mathematical expression. Among the existing techniques applied for the characterization of a surface, terrain modeling relates to the representation of the natural surface of the Earth. Cartographic terrain or relief models as threedimensional representations of a part of the Earth's surface convey an immediate and direct impression of a landscape and are much easier to understand than two-dimensional models. This paper addresses a major problem in complex surface modeling and evaluation consisting in the characterization of their topography and comparison among different textures, which can be relevant in different areas of research. A new algorithm is presented that allows calculating the fractal dimension of images of complex surfaces. The method is used to characterize different surfaces and compare their characteristics. The proposed new mathematical method computes the fractal dimension of the 3D space with the average space component of Hurst exponent H, while the estimated fractal dimension is used to evaluate, compare and characterize complex surfaces that are relevant in different areas of research. Various surfaces with both methods were analyzed and the results were compared. The study confirms that with known coordinates of a surface, it is possible to describe its complex structure. The estimated fractal dimension is proved to be an ideal tool for measuring the complexity of the various surfaces considered

Multiscale Bone Remodelling with Spatial P Systems

Many biological phenomena are inherently multiscale, i.e. they are characterized by interactions involving different spatial and temporal scales simultaneously. Though several approaches have been proposed to provide "multilayer" models, only Complex Automata, derived from Cellular Automata, naturally embed spatial information and realize multiscaling with well-established inter-scale integration schemas. Spatial P systems, a variant of P systems in which a more geometric concept of space has been added, have several characteristics in common with Cellular Automata. We propose such a formalism as a basis to rephrase the Complex Automata multiscaling approach and, in this perspective, provide a 2-scale Spatial P system describing bone remodelling. The proposed model not only results to be highly faithful and expressive in a multiscale scenario, but also highlights the need of a deep and formal expressiveness study involving Complex Automata, Spatial P systems and other promising multiscale approaches, such as our shape-based one already resulted to be highly faithful.Comment: In Proceedings MeCBIC 2010, arXiv:1011.005