13.07.2015 Views

Modeling bone regeneration around endosseous implants

Modeling bone regeneration around endosseous implants

Modeling bone regeneration around endosseous implants

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Summaryviiis not changed within this time period. Hence during the first four days thesources of osteogenic cells and growth factors are modeled to be located atthe old <strong>bone</strong> surface and at the implant surface. Therefore, the evolutionarydifferentiation approach is an essential part of the present model. Theclassical immediate differentiation formalism would lead to an immediatemovement of the ossification front and to not being able to define the initialsources of cells and growth factors for the current moving boundary problem,since the old <strong>bone</strong> surface and the implant surface would not be theboundary of the physical (i.e. soft tissue) domain.The numerical method for the present mathematical model is describedin Chapter 5. The current problem has a number of characteristics whichmake the algorithm elaborate and complex. The constructed numericalmethod provides stable and non-negative solutions of the nonlinearly coupledsystem of the time dependent taxis-diffusion-reaction equations. Theadditional challenges, which are faced with at the stage of the construction ofthe algorithm, consist in the need to discretize the model equations withinthe irregular and time-evolving physical domain, and in the sensitivity ofthe mathematical model with respect to negative solutions. The presentlyproposed numerical approach is based on such methods as: the method oflines, finite volume method, level set method and the embedded boundarymethod. For coarse meshes, patterns are observed to develop. These patternsshould not be there and result from numerical errors. Since we observethat convergence is only reached at high mesh resolutions, we think that furtherresearch should point into the direction of adaptive mesh refinement,improved time integration (operator splitting methods) and/or higher-ordermethods, such as spectral or DG-methods (however the requirement of thesolution positivity is still essential for higher-order methods). The series ofnumerical simulations demonstrate the ability of the present osseointegrationmodel to predict various paths of new <strong>bone</strong> formation depending on thechosen parameter values.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!