skip to main content

Explore Research Products in the NSF-PAR

Title: An Isogeometric Analysis Computational Platform for Material Transport Simulations in Complex Neurite Networks
Neurons exhibit remarkably complex geometry in their neurite networks. So far, how materials are transported in the complex geometry for survival and function of neurons remains an unanswered question. Answering this question is fundamental to understanding the physiology and disease of neurons. Here, we have developed an isogeometric analysis (IGA) based platform for material transport simulation in neurite networks. We modeled the transport process by reaction-diffusion-transport equations and represented geometry of the networks using truncated hierarchical tricubic B-splines (THB-spline3D). We solved the Navier-Stokes equations to obtain the velocity field of material transport in the networks. We then solved the transport equations using the streamline upwind/Petrov-Galerkin (SU/PG) method. Using our IGA solver, we simulated material transport in three basic models of the network geometry: a single neurite, a neurite bifurcation, and a neurite tree with three bifurcations. In addition, the robustness of our solver is illustrated by simulating material transport in three representative and complex neurite networks. From the simulation we discovered several spatial patterns of the transport process. Together, our simulation provides key insights into how material transport in neurite networks is mediated by their complex geometry.  more » « less
Award ID(s):
1804929
NSF-PAR ID:
10100412
Author(s) / Creator(s):
Date Published:
Journal Name:
Molecular & cellular biomechanics
Volume:
16
Issue:
2
ISSN:
1556-5297
Page Range / eLocation ID:
123-140
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ( , Scientific Reports)
    null (Ed.)
    Abstract Neurons exhibit complex geometry in their branched networks of neurites which is essential to the function of individual neuron but also brings challenges to transport a wide variety of essential materials throughout their neurite networks for their survival and function. While numerical methods like isogeometric analysis (IGA) have been used for modeling the material transport process via solving partial differential equations (PDEs), they require long computation time and huge computation resources to ensure accurate geometry representation and solution, thus limit their biomedical application. Here we present a graph neural network (GNN)-based deep learning model to learn the IGA-based material transport simulation and provide fast material concentration prediction within neurite networks of any topology. Given input boundary conditions and geometry configurations, the well-trained model can predict the dynamical concentration change during the transport process with an average error less than 10% and $$120 \sim 330$$ 120 ∼ 330 times faster compared to IGA simulations. The effectiveness of the proposed model is demonstrated within several complex neurite networks. 
    more » « less
  2. ; ; ( , Scientific reports)
    Neurons exhibit complex geometry in their branched networks of neurites which is essential to the function of individual neuron but also brings challenges to transport a wide variety of essential materials throughout their neurite networks for their survival and function. While numerical methods like isogeometric analysis (IGA) have been used for modeling the material transport process via solving partial differential equations (PDEs), they require long computation time and huge computation resources to ensure accurate geometry representation and solution, thus limit their biomedical application. Here we present a graph neural network (GNN)-based deep learning model to learn the IGA-based material transport simulation and provide fast material concentration prediction within neurite networks of any topology. Given input boundary conditions and geometry configurations, the well-trained model can predict the dynamical concentration change during the transport process with an average error less than 10% and 120∼330 times faster compared to IGA simulations. The effectiveness of the proposed model is demonstrated within several complex neurite networks. 
    more » « less
  3. ; ( , Scientific Reports)
    Abstract The intracellular transport process plays an important role in delivering essential materials throughout branched geometries of neurons for their survival and function. Many neurodegenerative diseases have been associated with the disruption of transport. Therefore, it is essential to study how neurons control the transport process to localize materials to necessary locations. Here, we develop a novel optimization model to simulate the traffic regulation mechanism of material transport in complex geometries of neurons. The transport is controlled to avoid traffic jam of materials by minimizing a pre-defined objective function. The optimization subjects to a set of partial differential equation (PDE) constraints that describe the material transport process based on a macroscopic molecular-motor-assisted transport model of intracellular particles. The proposed PDE-constrained optimization model is solved in complex tree structures by using isogeometric analysis (IGA). Different simulation parameters are used to introduce traffic jams and study how neurons handle the transport issue. Specifically, we successfully model and explain the traffic jam caused by reduced number of microtubules (MTs) and MT swirls. In summary, our model effectively simulates the material transport process in healthy neurons and also explains the formation of a traffic jam in abnormal neurons. Our results demonstrate that both geometry and MT structure play important roles in achieving an optimal transport process in neuron. 
    more » « less
  4. ; ( , Journal of Mechanics)
    Abstract

    The intracellular transport process plays an important role in delivering essential materials throughout branched geometries of neurons for their survival and function. Many neurodegenerative diseases have been associated with the disruption of transport. Therefore, it is essential to study how neurons control the transport process to localize materials to necessary locations. Here, we develop a novel optimization model to simulate the traffic regulation mechanism of material transport in three-dimensional complex geometries of neurons. The transport is controlled to avoid traffic jams of materials by minimizing a predefined objective function. The optimization subjects to a set of partial differential equation (PDE) constraints that describe the material transport process based on a macroscopic molecular-motor-assisted transport model of intracellular particles. The proposed PDE-constrained optimization model is solved in complex tree structures by using the isogeometric analysis. Different simulation parameters are used to introduce traffic jams and study how neurons handle the transport issue. Specifically, we successfully model and explain the traffic jam caused by the reduced number of microtubules (MTs) and MT swirls. In summary, our model effectively simulates the material transport process in healthy neurons and also explains the formation of a traffic jam in abnormal neurons. Our results demonstrate that both geometry and MT structure play important roles in achieving an optimal transport process in neurons.

     
    more » « less
  5. ; ( , ASME 2021 Fluids Engineering Division Summer Meeting)
    Abstract

    Insects rely on their olfactory system to forage, prey, and mate. They can sense odorant plumes emitted from sources of their interests with their bilateral odorant antennae, and track down odor sources using their highly efficient flapping-wing mechanism. The odor-tracking process typically consists of two distinct behaviors: surging upwind and zigzagging crosswind. Despite the extensive numerical and experimental studies on the flying trajectories and wing flapping kinematics during odor tracking flight, we have limited understanding of how the flying trajectories and flapping wings modulate odor plume structures. In this study, a fully coupled three-way numerical solver is developed, which solves the 3D Navier-Stokes equations coupled with equations of motion for the passive flapping wings, and the odorant convection-diffusion equation. This numerical solver is applied to investigate the unsteady flow field and the odorant transport phenomena of a fruit fly model in both surging upwind and zigzagging crosswind cases. The unsteady flow generated by flapping wings perturbs the odor plume structure and significantly impacts the odor intensity at the olfactory receptors (i.e., antennae). During zigzagging crosswind flight, the differences in odor perception time and peak odor intensity at the receptors potentially help create stereo odorant mapping to track odor source. Our simulation results will provide new insights into the mechanism of how fruit flies perceive odor landscape and inspire the future design of odor-guided micro aerial vehicles (MAVs) for surveillance and detection missions.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email