{"id":625,"date":"2021-08-15T17:50:02","date_gmt":"2021-08-15T17:50:02","guid":{"rendered":"http:\/\/summergeometry.org\/sgi2021\/?p=625"},"modified":"2021-08-15T17:50:04","modified_gmt":"2021-08-15T17:50:04","slug":"geometric-modelling-iga-and-cool-simulations","status":"publish","type":"post","link":"https:\/\/summergeometry.org\/sgi2021\/geometric-modelling-iga-and-cool-simulations\/","title":{"rendered":"Geometric Modelling, IGA and cool simulations"},"content":{"rendered":"\n

Isogeometric analysis (IGA) is an analysis technique that combines computer aided design (CAD) with traditional finite element analysis (FEA). It uses the same spline basis functions to construct the geometry and the solution space, which is beneficial as traditional FEA requires geometric approximation that can lead to inaccuracies [1].<\/p>\n\n\n\n

For our week 3 project, my group and I used IGA simulation software that has been developed for modelling material transport in neurons. The software solves Navier-Stokes equations to obtain the velocity field and models the transport process by reaction-diffusion-transport equations. For our purposes, we used the solver to simulate material transport in complex neurons and heat transfer processes in various geometries including a simple block model and a rod model [2].\u00a0<\/p>\n\n\n\n

Our project was split into two main parts: geometric modeling and analysis. For geometric modelling, we used two open source software packages (HexGen and Hex2Spline [3]) for the construction of geometries, and we then used the aforementioned IGA software for simulation purposes. I was in charge of using the IGA software to run the simulations and visualizing the results with Paraview (a visualization application).\u00a0\u00a0<\/p>\n\n\n\n

In this post, I\u2019ll demonstrate our results using the simple model the team created and also share some of the results I got using the IGA software on some additional models.<\/p>\n\n\n\n

Here is the example model we used:<\/p>\n\n\n\n

\"\"<\/figure><\/div>\n\n\n\n

After the geometric modeling stage, I received 3 input files that contained the control mesh, an initial velocity field and the simulation parameters (particle concentrations, diffusion coefficient, velocities of material, etc.). These files were used as input for the IGA solver [2], which consists of four stages: spline construction, mesh partition for parallel computing, solving Navier-Stokes equations, and final transport simulation.\u00a0<\/p>\n\n\n\n

Due to the high computational needs of the last two stages, we set up the simulation environment at the Pittsburgh Supercomputer Center<\/a> (PSC). We ran the simulation by connecting to the remote supercomputer\u00a0via an ssh client. Finally, after getting the results, I was able to visualize them using Paraview.\u00a0<\/p>\n\n\n\n

\"\"<\/figure><\/div>\n\n\n\n

In the above example, the color bar shows the transport velocity magnitude. I also added streamlines that show the direction of flow.\u00a0<\/p>\n\n\n\n

Here is an animation of our results that depicts the movement of the particles:<\/p>\n\n\n\n

I also ran the solver on a rod model (that Angran Li<\/a>, our mentor, provided):<\/p>\n\n\n\n

\"\"<\/figure><\/div>\n\n\n\n

The model shows the flow of heat through the rod as it is input from the right end of the model. Again, the colorbar represents the magnitude of the velocity of particles at that point. <\/p>\n\n\n\n

Finally, I also ran the IGA solver on a neuron model (found on NeuroMorpho.org<\/em><\/a> and edited with the help of Angran). <\/p>\n\n\n\n

\"\"
The geometry of the neuron (visualised using Paraview)<\/figcaption><\/figure><\/div>\n\n\n\n