Alan M. Shih

Parallel unstructured mesh generation by an advancing front method

Mesh generation is a critical step in high fidelity computational simulations. High-quality and high-density meshes are required to accurately capture the complex physical phenomena. A robust approach for a parallel framework has been developed to generate large-scale meshes in a short period of time. A coarse tetrahedral mesh is generated first to provide the basis of block interfaces and then is partitioned into a number of sub-domains using METIS partitioning algorithms. A volume mesh is generated on each sub-domain in parallel using an advancing front method. Dynamic load balancing is achieved by evenly distributing work among the processors. All the sub-domains are combined to create a single volume mesh. The combined volume mesh can be smoothed to remove the artifacts in the interfaces between sub-domains. A void region is defined inside each sub-domain to reduce the data points during the smoothing operation. The scalability of the parallel mesh generation is evaluated to quantify the improvement on shared- and distributed-memory computer systems.

Multiple marching direction approach to generate high-quality hybrid meshes

This paper describes the method to generate hybrid meshes composed of triangular prisms, pyramids, hexahedra, and tetrahedra for viscous computational fluid dynamics simulations. Surface triangulation is performed using a direct advancing front method or a modified mesh-decimation method. From a surface mesh, a near-field mesh is generated using an advancing layer approach. To generate high-quality meshes and to mesh around singular points, multiple marching directions are prepared from nodes on sharp convex corners. Special placement of elements around the sharp convex corners avoids using generalized elements. Tetrahedral meshing is then performed to fill the rest of the domain using an advancing front method. The hybrid mesh generation method is applied to several models to demonstrate its capability.

Turbulent flow evaluation of the venous needle during hemodialysis.

Arteriovenous (AV) grafts and fistulas used for hemodialysis frequently develop intimal hyperplasia (IH) at the venous anastomosis of the graft, leading to flow-limiting stenosis, and ultimately to graft failure due to thrombosis. Although the high AV access blood flow has been implicated in the pathogenesis of graft stenosis, the potential role of needle turbulence during hemodialysis is relatively unexplored. High turbulent stresses from the needle jet that reach the venous anastomosis may contribute to endothelial denudation and vessel wall injury. This may trigger the molecular and cellular cascade involving platelet activation and IH, leading to eventual graft failure. In an in-vitro graft/needle model dye injection flow visualization was used for qualitative study of flow patterns, whereas laser Doppler velocimetry was used to compare the levels of turbulence at the venous anastomosis in the presence and absence of a venous needle jet. Considerably higher turbulence was observed downstream of the venous needle, in comparison to graft flow alone without the needle. While turbulent RMS remained around 0.1 m/s for the graft flow alone, turbulent RMS fluctuations downstream of the needle soared to 0.4-0.7 m/s at 2 cm from the tip of the needle and maintained values higher than 0.1 m/s up to 7-8 cm downstream. Turbulent intensities were 5-6 times greater in the presence of the needle, in comparison with graft flow alone. Since hemodialysis patients are exposed to needle turbulence for four hours three times a week, the role of post-venous needle turbulence may be important in the pathogenesis of AV graft complications. A better understanding of the role of needle turbulence in the mechanisms of AV graft failure may lead to improved design of AV grafts and venous needles associated with reduced turbulence, and to pharmacological interventions that attenuate IH and graft failure resulting from turbulence.

An Approach to Generate High Quality Unstructured Hybrid Meshes

This paper describes the method to generate hybrid meshes comprised of triangular prisms, pyramids, hexahedra and tetrahedra for viscous computational fluid dynamics (CFD) simulations. Surface triangulation is performed using a direct advancing front method or a modified mesh-decimation method. From a surface mesh, a near-field mesh is generated using an advancing layer approach. To generate high quality meshes and to mesh around singular points, multiple marching directions are prepared from nodes on sharp convex corners. Tetrahedral meshing is then performed to fill the rest of the domain using an advancing front method. The hybrid mesh generation method is applied to several models to demonstrate its capability.

Hemodynamic Analysis of a Compliant Femoral Artery Bifurcation Model using a Fluid Structure Interaction Framework

The influence of wall motion on the hemodynamic characteristics of the human femoral bifurcation and its effects on the development of peripheral artery disease has not been previously investigated. This study aimed in investigating the hemodynamics of a compliant patient-specific femoral artery bifurcation model by a fluid structure interaction (FSI) scheme. The complex physiological geometry of the femoral artery bifurcation was reproduced from sequentially obtained transverse CT scan images. Velocity waveforms derived from phase contrast MR images were extracted and mapped to define boundary conditions. Equations governing blood flow and wall motion were solved using an FSI framework that utilizes commercial codes: FLUENT for computational fluid dynamics and ANSYS for computational structural dynamics. The results showed that wall compliance decreased flow velocities at the relatively high curvature geometries including common and superficial femoral artery (SFA), and it created strong recirculation in the profunda femoris artery close to the bifurcation. In the SFA region near the apex, time averaged wall shear stress (TAWSS) differences up to 25% between compliant and rigid models were observed. The compliant model also exhibited lower TAWSS and oscillatory shear at the superior section of the common femoral artery close to the bifurcation. The presence of wall motion, however, created minor differences in the general flow-field characteristics. We conclude that wall motion does not have significant influence on the global fluid dynamic characteristics of the femoral artery bifurcation. Longer arterial segments need to be simulated to see the effect of wall motion on tortuousity which was previously cited as an important factor in the development of atherosclerosis at the femoral artery.

Finite element model development of a child pelvis with optimization-based material identification

A finite element (FE) model of a 10-years-old child pelvis was developed and validated against experimental data from lateral impacts of pediatric pelves. The pelvic bone geometry was reconstructed from a set of computed tomography images, and a hexahedral mesh was generated using a new octree-based hexahedral meshing technique. Lateral impacts to the greater trochanter and iliac wing of the seated pelvis were simulated. Sensitivity analysis was conducted to identify material parameters that substantially affected the model response. An optimization-based material identification method was developed to obtain the most favorable material property set by minimizing differences in biomechanical responses between experimental and simulation results. This study represents a pilot effort in the development and validation of age-dependent musculoskeletal FE models for children, which may ultimately serve to evaluate injury mechanisms and means of protection for the pediatric population.

Computational Fluid Dynamic Analysis of the Posterior Airway Space After Maxillomandibular Advancement for Obstructive Sleep Apnea Syndrome

PURPOSE This study evaluated the soft tissue change of the upper airway after maxillomandibular advancement (MMA) using computational fluid dynamics. MATERIALS AND METHODS Eight patients with obstructive sleep apnea syndrome who required MMA were recruited into this study. All participants underwent pre- and postoperative computed tomography and then MMA by a single oral and maxillofacial surgeon. Upper airway computed tomographic datasets for these 8 patients were created with high-fidelity 3-dimensional numerical models for computational fluid dynamics. The 3-dimensional models were simulated and analyzed to study how changes in airway anatomy affect the pressure effort required for normal breathing. Airway dimensions, skeletal changes, apnea-hypopnea index, and pressure effort of pre- and postoperative 3-dimensional models were compared and correlations were interpreted. RESULTS After MMA, laminar and turbulent air flows were significantly decreased at every level of the airway. The cross-sectional areas at the soft palate and tongue base were significantly increased. CONCLUSIONS This study showed that MMA increased airway dimensions by increasing the distance from the occipital base to the pogonion. An increase of this distance showed a significant correlation with an improvement in the apnea-hypopnea index and a decreased pressure effort of the upper airway. Decreasing the pressure effort will decrease the breathing workload. This improves the condition of obstructive sleep apnea syndrome.

Efficient Computational Fluid Dynamics Evaluation of Small-Device Locations with Automatic Local Remeshing

This paper describes a new efficient automatic remeshing method for three-dimensional hybrid meshes for viscous flow simulations to accommodate the meshes with changes of small devices quickly and easily. This remeshing method has two notable advantages. First, hybrid meshes can be generated automatically from a baseline mesh when the small devices are moved and/or deformed. This enables shape optimization techniques to find better models quickly, because tens or hundreds of meshes are usually required during the process. Second, the updated hybrid meshes are the same as the baseline mesh except for elements around the small devices. Their effect can be evaluated more accurately. Solution data from the baseline mesh can be shared with new hybrid meshes (e.g., as an initial condition to expedite the convergence of computational simulations). The remeshing method is applied to the Japan Aerospace Exploration Agencys high-lift-configuration standard model with a nacelle chine in different locations to demonstrate its capability. Computational simulations are also performed to further discuss the effectiveness of the remeshing method.

Original articles: Patient-specific geometry modeling and mesh generation for simulating Obstructive Sleep Apnea Syndrome cases by Maxillomandibular Advancement

The objective of this paper is the reconstruction of upper airway geometric models as hybrid meshes from clinically used Computed Tomography (CT) data sets in order to understand the dynamics and behaviors of the pre- and postoperative upper airway systems of Obstructive Sleep Apnea Syndrome (OSAS) patients by viscous Computational Fluid Dynamics (CFD) simulations. The selection criteria for OSAS cases studied are discussed because two reasonable pre- and postoperative upper airway models for CFD simulations may not be created for every case without a special protocol for CT scanning. The geometry extraction and manipulation methods are presented with technical barriers that must be overcome so that they can be used along with computational simulation software as a daily clinical evaluation tool. Eight cases are presented in this paper, and each case consists of pre- and postoperative configurations. The results of computational simulations of two cases are included in this paper as demonstration.

Assessment of Surgical Effects on Patients with Obstructive Sleep Apnea Syndrome Using Computational Fluid Dynamics Simulations.

Obstructive sleep apnea syndrome is one of the most common sleep disorders. To treat patients with this health problem, it is important to detect the severity of this syndrome and occlusion sites in each patient. The goal of this study is to test the hypothesis that the cure of obstructive sleep apnea syndrome by maxillomandibular advancement surgery can be predicted by analyzing the effect of anatomical airway changes on the pressure effort required for normal breathing using a high-fidelity, 3-D numerical model. The employed numerical model consists of: 1) 3-D upper airway geometry construction from patient-specific computed tomographic scans using an image segmentation technique, 2) mixed-element mesh generation of the numerically constructed airway geometry for discretizing the domain of interest, and 3) computational fluid dynamics simulations for predicting the flow field within the airway and the degree of severity of breathing obstruction. In the present study, both laminar and turbulent flow simulations were performed to predict the flow field in the upper airway of the selected patients before and after maxillomandibular advancement surgery. Patients of different body mass indices were also studied to assess their effects. The numerical results were analyzed to evaluate the pressure gradient along the upper airway. The magnitude of the pressure gradient is regarded as the pressure effort required for breathing, and the extent of reduction of the pressure effort is taken to measure the success of the surgery. The description of the employed numerical model, numerical results from simulations of various patients, and suggestion for future work are detailed in this paper.