Impingement Location Research Articles

Inertial collapse of a gas bubble in a shear flow near a rigid wall

Despite the extensive research on bubble collapse near rigid walls, the bubble collapse dynamics in the presence of shear flow near a rigid wall is poorly understood. We conduct direct simulations of the Navier–Stokes equations to explore the bubble dynamics and pressures during bubble collapse near a rigid, flat wall under linear shear flow conditions. We examine the dependence of the bubble collapse morphology and wall pressures on the initial bubble location and shear rate. We find that shear distorts the bubble, generating two re-entrant jets – one developing from the side opposite to the mean flow and the other from the far end toward the wall. Upon impact of the jet on the opposite side of the bubble, water-hammer shocks are produced, which propagate outward and interact with the convoluted bubble shape. The shock stretches the bubble towards the wall, resulting in a closer impact location for the jet originating from the far end compared with the case with no shear flow. The water-hammer pressure location can be approximated as the theoretical distance travelled by a particle initialised at the bubble centre with the corresponding constant shear flow velocity. The maximum wall pressures can thus be predicted by considering the distance between the far jet impingement location and the wall along the wall-normal direction. As the shear rate is increased, the maximum wall pressure increases, although only marginally. We determine the critical initial stand-off distance from the wall at which the bubble morphology is shear dominated, i.e. characterised by converging re-entrant jets.

Nonlinear aeroelastic behavior of a two-dimensional heated panel by irregular shock reflection considering viscoelastic damping

This paper investigates the aeroelastic stability and nonlinear aeroelastic behavior of a two-dimensional heated panel in irregular shock reflection and extends prior work to include the effects of viscoelasticity. The aeroelastic model is formulated using the von Kármán large deflection plate theory and the Kelvin–Voigt damping model, accompanied by the quasi-steady thermal stress theory. The unsteady aerodynamic pressure is evaluated through the piston theory and the compressibility-corrected potential theory. The Galerkin approach is used to discretize the governing equation. The Lyapunov indirect method is applied to conduct theoretical analysis, obtaining the aeroelastic stability boundary. Also, the nonlinear aeroelastic response is numerically simulated via the fourth-order Runge–Kutta method. The proper orthogonal decomposition is applied to the panel deflection to manifest the influence of various system parameters. It is demonstrated that the shock wave aggravates the aerodynamic heating, lowering the critical buckling temperature. The viscoelastic damping restricts the impact of shock impingement location and shock strength on the stability boundary and also transforms the chaotic motions into periodic LCOs.

Dynamic analysis of breaking wave impact on a floating offshore wind turbine via smoothed particle hydrodynamics

This work leverages Lagrangian smoothed particle hydrodynamics (SPH) to explore the structural and hydrodynamic response of floating offshore wind turbines (FOWT) subject to impulsive breaking waves. The SPH formulation was first validated against breaking wave impact on a model tension leg platform (TLP) which demonstrated good consistency with experimental results. Following validation, wave focusing was utilized to generate both breaking and nonbreaking extreme waves impacting a moored semi-submersible FOWT at full scale. Impulsive forces and accelerations resulting from the plunging breaker were observed to exceed that of nonbreaking waves by up to 70 % and 230 %, respectively, and were highly sensitive to the wave impingement location relative to the FOWT. However, wave breaking did not appear to significantly influence rigid body motions and yielded lower mooring tensions than its nonbreaking counterpart due to the short duration of impact. This work ultimately demonstrates the applicability of SPH for the simulation of breaking wave interactions with floating bodies and provides further impetus towards the study of FOWTs under such conditions.

Numerical Study of the Effect of the Port Angle of the Superior Vena Cava Supplying Cannula on Hemodynamics in the Right Atrium in VV-ECMO.

Objective: To elucidate the pattern of the influence of the port angle of the superior vena cava supplying cannula (SVCS) on hemodynamics within the right atrium in VV-ECMO. Methods: A three-dimensional model of the right atrium was established based on CT images of a real patient. The 3D models of the SVCS and inferior vena cava draining cannula (IVCD) were established based on the Edwards 18Fr and Medos 22Fr real intubation models, respectively. Based on these models, three-dimensional models of the SVCS ports with bending angles of -90°, -60°, -30°, 0°, 30°, 60°, and 90° in the plane formed by the centerline of the SVCS and the center point of the tricuspid valve (TV) were established. Transient-state computational fluid dynamics (CFD) was performed to clarify the right atrium blood flow pattern and hemodynamic states at different SVCS port orientation angles. The velocity clouds, wall pressure, wall shear stress (WSS), relative residence time (RRT), and recirculation fraction (RF) were calculated to assess hemodynamic changes in the right atrium at different angles of the port of the SVCS. Results: As the angle of the port of the superior chamber cannula changed, the location of the high-velocity blood impingement from the SVCS changed, and the pattern of blood flow within the right atrium was dramatically altered. The results for the maximum right atrial wall pressure were 13,472 pa, 13,424 pa, 10,915 pa, 7680.2 pa, 5890.3 pa, 5597.6 pa, and 7883.5 pa (-90° vs. -60° vs. -30° vs. 0° vs. 30° vs. 60° vs. 90°), and the results for the mean right atrial wall pressure were 6788.9 pa, 8615.1 pa, 8684.9 pa, 6717.2 pa, 5429.2 pa, 5455.6 pa, and 7117.8 pa ( -90° vs. -60° vs. -30° vs. 0° vs. 30° vs. 60° vs. 90°). The results of the maximum right atrial wall WSS in the seven cases were 63.572 pa, 55.839 pa, 31.705 pa, 39.531 pa, 40.11 pa, 28.474 pa, and 35.424 (-90° vs. -60° vs. -30° vs. 0° vs. 30° vs. 60° vs. 90°), respectively, and the results of the mean right atrial wall WSS results were 3.8589 pa, 3.6706 pa, 3.3013 pa, 3.2487 pa, 2.3995 pa, 1.3304 pa, and 2.0747 pa (-90° vs. -60° vs. -30° vs. 0° vs. 30° vs. 60° vs. 90°), respectively. The results for the area percentage of high RRT in the seven cases were 3.44%, 2.23%, 4.24%, 1.83%, 3.69%, 7.73%, and 3.68% (-90° vs. -60° vs. -30° vs. 0° vs. 30° vs. 60° vs. 90°), and the results for the RF were 21.57%, 23.24%, 19.78%, 12.57%, 10.24%, 5.07%, and 8.05% (-90° vs. -60° vs. -30° vs. 0° vs. 30° vs. 60° vs. 90°). Conclusions: The more the port of the SVCS is oriented toward the TV, the more favorable it is for reducing RF and the impingement of blood flow in the right atrial wall, but there may be an increased risk of RRT. The opposite orientation of the SVCS port to the TV is not conducive to reducing flow impingement on the right atrial wall and RF.

Interpreting Heat Flux Measurements in a Vitiated Backward-Facing Step Flow

The investigation of the relative impact of convective and radiative heat flux from vitiated gases to a combustor wall was performed in a backward-facing step combustor. Backward-facing step flows have many flow features of gas turbine combustors (recirculation, impingement, and boundary-layer development), but in a two-dimensional, nonproprietary geometry. Reynolds number, gas temperature, and plate temperature were varied as heat flux was measured in different regions of the flow. Measurements were made using a heat flux sensor and radiometer, which measure total and radiative heat flux, respectively. The highest total heat flux was measured at the impingement location and the lowest in the recirculation zone. Separating the total heat flux into the radiative and convective components is not straightforward and requires simulations of the experiment. Spectrally resolved Monte Carlo ray-tracing simulations are used to scale the radiative heat flux measurements for the separation of convective and radiative heat flux from the total heat flux measurements. These calculations show that only about 30% of incident radiation is captured by the radiometer due to spectral differences between the two sensors. Once the radiation measurement is corrected, experimental results show higher radiative than convective heat transfer at low Reynolds numbers (Reh=2000), highlighting the impact of radiative heat transfer in hot gas environments.

Multivariate linear-mixed analysis of changes in anterior inferior iliac spine impingement incidence with posterior pelvic tilt: a computer simulation study.

It is well known that increased posterior tilt of the pelvis is an effective strategy for avoiding impingement of the femur with the pelvis during movement. Daily repetitive collisions become mechanical loads, and the more frequently they occur, the more tissue damage and pain they cause. Therefore, reducing the rate of occurrence of impingement is important to avoid aggravation of symptoms. This study aimed to evaluate the effects of changes in posterior pelvic tilt on the risk of impingement between the femur and the anterior inferior iliac spine (AIIS)/subspine in various functional postures. Patients with femoroacetabular impingement syndrome (FAIS) who were candidates for hip arthroscopic osteochondroplasty between October 2013 and June 2020 were included. A three-dimensional reconstructed model was used to simulate the incidence of impingement at 12 hip positions required for activities of daily living. We predicted value of the spatial incidence of impingement assumed that hip motion should exceed 130/30 degrees without impingement. Impingement was measured at three pelvic positions: an anterior tilt of 10°, in the functional pelvic plane and a posterior tilt of 10°. Multivariate linear-mixed models were used to assess the effect of covariate-adjusted posterior pelvic tilt on the impingement incidence in the AIIS region. AIIS type, center-edge angle, acetabular version and femoral version were used as covariates. The impingement rates and locations of the three pelvic tilt postures were assessed. Seventy-eight patients (60 males and 18 females; average age, 46 ± 15.1 years) with FAIS were analyzed. A multivariate linear-mixed model revealed a coefficient of -0.8% (95% confidence interval -0.9 to -0.7%; P < 0.001) for posterior pelvic tilt. Thus, posterior pelvic tilt affects AIIS impingement incidence. After adjusting for anatomical effects, the posterior pelvic tilt should be addressed to avoid impingement.

On the evolution of turbulent boundary layers during flame–wall interaction investigated by highly resolved laser diagnostics

The turbulent boundary layer behavior in the presence of flame–wall interactions (FWI) has an important role on the mass and energy transfer at the gas/solid interface. Detailed experiments resolving the turbulent boundary layer evolution in the presence of FWI are lacking, which impedes knowledge. This work presents a combination of particle image velocimetry (flow field), dual-pump coherent anti-Stokes Raman spectroscopy (gas temperature), and OH laser induced fluorescence (flame topology) measurements to study the evolution of the boundary layer structure in the presence of FWI. Experiments are conducted in a side-wall quenching (SWQ) burner. Findings reveal that the reacting boundary layer flow adheres to the linear scaling law u+=y+ in the viscous sublayer until y+ = 5. Beyond y+ = 5, the flame modifies the velocity and temperature field such that the uz+ streamwise velocity deviates from the viscous sublayer and the law-of-the-wall scaling in the log-layer with uz+ being smaller than that of the non-reacting flow (the subscript z refers to the streamwise coordinate and is used throughout this manuscript). As the fluid approaches the flame impingement location at the wall, the gas temperature increases significantly, causing a threefold increase in kinematic viscosity, ν. Although the near-wall streamwise velocity gradient d〈Uz〉/dy|y=0mm decreases, the larger increase in ν reduces uz+ and leads to the deviation from the law-of-the-wall. Downstream the flame impingement location, ν is relatively constant and uz+ values begin to approach those of the law-of-the-wall. Trends are presented for SWQ and head-on quenching flame topologies, and are intended to help development of more accurate wall models.

The Structure and Stability of Premixed CH4, H2, and NH3/H2 Flames in an Axially Staged Can Combustor

Abstract An experimental investigation of flame structure, stability, and emissions performance was conducted in a two-stage lab-scale generic combustor design operated with CH4, H2, and NH3/H2 fuel blends. The main flame zone features a premixed bluff body stabilized flame, with a secondary flame zone initiated downstream by injecting premixed air and fuel using two opposing radial jets. The total power and air flowrate are kept constant between the different fueling cases, while the air split between stages and equivalence ratios are varied to explore conditions relevant to gas turbine operation. Given the relative novelty of the configuration, special emphasis is given to analyzing the structure of the opposing jet flames in the secondary stage. In contrast to previous literature on reacting jets in cross flow, these interact significantly due to their proximity, leading to a merged flame zone at the impingement location in the center of the combustion chamber, and some flame propagation upstream of the jet location. As the jet-to-crossflow momentum ratio increases, the merged flame zone changes shape, reaching close to the walls for the methane cases but remaining very compact when operating with almost pure hydrogen. For the hydrogen flames, diverting more air to the second stage allows higher total thermal power conditions to be reached, while avoiding flashback, and eliminates thermoacoustic instabilities. For ammonia-hydrogen flames, air is diverted to the second stage, while a constant fuel flow is sent to the primary stage, resulting in some locally rich conditions in the primary flame. A local minima in terms of NOX occurs when the primary flame is operated at an equivalence ratio of 1.15. Analysis of the flame structure suggests that this state corresponds to almost complete combustion or pyrolysis of NH3 in the main flame, with the remaining hydrogen burned in an inverse diffusion flame in the secondary zone.

Hip Impingement of severe SCFE patients after in situ pinning causes decreased flexion and forced external rotation in flexion on 3D-CT.

In situ pinning is an accepted treatment for stable slipped capital femoral epiphysis. However, residual deformity of severe slipped capital femoral epiphysis can cause femoroacetabular impingement and forced external rotation. The aim of this study was to evaluate the (1) hip external rotation and internal rotation in flexion, (2) hip impingement location, and (3) impingement frequency in early flexion in severe slipped capital femoral epiphysis patients after in situ pinning using three-dimensional computed tomography. A retrospective Institutional Review Board-approved study evaluating 22 patients (26 hips) with severe slipped capital femoral epiphysis (slip angle > 60°) using postoperative three-dimensional computed tomography after in situ pinning was performed. Mean age at slipped capital femoral epiphysis diagnosis was 13 ± 2 years (58% male, four patients bilateral, 23% unstable, 85% chronic). Patients were compared to contralateral asymptomatic hips (15 hips) with unilateral slipped capital femoral epiphysis (control group). Pelvic three-dimensional computed tomography after in situ pinning was used to generate three-dimensional models. Specific software was used to determine range of motion and impingement location (equidistant method). And 22 hips (85%) underwent subsequent surgery. (1) Severe slipped capital femoral epiphysis patients had significantly (p < 0.001) decreased hip flexion (43 ± 40°) and internal rotation in 90° of flexion (-16 ± 21°, IRF-90°) compared to control group (122 ± 9° and 36 ± 11°). (2) Femoral impingement in maximal flexion was located anterior to anterior-superior (27% on 3 o'clock and 27% on 1 o'clock) of severe slipped capital femoral epiphysis patients and located anterior to anterior-inferior (38% on 3 o'clock and 35% on 4 o'clock) in IRF-90°. (3) However, 21 hips (81%) had flexion < 90° and 22 hips (85%) had < 10° of IRF-90° due to hip impingement and 21 hips (81%) had forced external rotation in 90° of flexion (< 0° of IRF-90°). After in situ pinning, patient-specific three-dimensional models showed restricted flexion and IRF-90° and forced external rotation in 90° of flexion due to early hip impingement and residual deformity in most of the severe slipped capital femoral epiphysis patients. This could help to plan subsequent hip preservation surgery, such as hip arthroscopy or femoral (derotation) osteotomy.

Biorthogonal decomposition of the disturbance flow field generated by particle impingement on a hypersonic boundary layer

The disturbance flow field in a hypersonic boundary layer excited by particle impingement was investigated with a focus on the first stage of the laminar-to-turbulent transition process, namely the receptivity process. A previously validated direct numerical simulation approach adopting disturbance flow tracking is used to simulate the particle-induced transition process. Particle impingement generates a highly complex disturbance flow field that can be characterised by a wide range of frequencies and wavenumbers. After providing some insight about the spectral characteristics of the disturbance flow field in the frequency and wavenumber domains, biorthogonal decomposition is employed to reveal the composition of the disturbance flow field consisting of different continuous and discrete eigenmodes that are triggered through particle impingement. The disturbance flow characteristics for different frequency and wavenumber pairs are discussed where large contributions in the disturbance flow spectrum are observed in the vicinity of the impingement location. A significant amount of the disturbance energy is diverted into the free stream leading to large coefficients of projection for the slow and fast acoustic branches while contributions to the entropy and vorticity branches are negligible. In addition to the continuous acoustic spectra, the first-, second- and other higher-order Mack modes are activated and provide large contributions to the disturbance flow field inside the boundary layer. Finally, it is demonstrated that the disturbance flow field in the vicinity of the impingement location can be reconstructed with a maximum relative error of $2.3\,\%$ by employing a theoretical biorthogonal eigenfunction system expansion and by considering contributions from fast and slow acoustic waves and at most four discrete modes only.

Poster 214: The Effect of Femoral Retroversion on Hip Range of Motion – A 3D Computational Analysis

Objectives: A subset of patients undergoing hip arthroscopy suffers from persistent pain and dysfunction which can be secondary to a number of factors. Femoral rotational deformities are often overlooked and/or providers are unsure when to address femoral deformities. Currently, there is limited to no data guiding surgeons when to consider a femoral osteotomy at the time of hip arthroscopy for FAI. Through previously validated computer-assisted 3D modeling software, this study sought to answer two questions: 1) How much femoral retroversion is too much to render cam femoroplasty ineffective in improving range of motion (ROM)? 2) In the setting of an osteotomy, what degree of femoral version correction would be sufficient to improve impingement free ROM? Methods: Ten patients were selected with symptomatic focal cam impingement lesions who underwent hip arthroscopy and had pre- and post-operative computed tomography (CT) scans. These patients were screened for normal acetabular morphology and demonstrated no radiographic signs of osteoarthritis. Pre- and post-operative 3D modeling of the involved hip was performed and subsequently dynamized with the femoral version set at (15, 5, -5 and -15 degrees) to define the impingement free ROM and location of impingement with maximum hip flexion and maximum hip internal rotation (IR) at 90 degrees of hip flexion. Results: Mean alpha angle pre and post cam resection were 66°±7.6 ° and 44°±4.9° respectively and were significantly different (p=0.003). Location of impingement shifted during hip flexion from the medial to anterior femoral neck with decreasing femoral version while the IR impingement location remained on the anterior femoral neck and was unaffected by femoral version. Kruskal-Wallis testing demonstrated decreasing femoral version significantly decreased hip flexion and IR ROM in both the pre and post cam resection groups (p&lt;0.0003). Impingement free ROM, as defined by hip flexion and IR greater than 90° and 10° respectively, was only identified in the 5 and 15 degree femoral anteversion models irrespective of whether cam femoroplasty had been performed. In the post cam resection group, hip IR ROM was significantly less with any level of decreasing femoral version (p&lt;0.002) while hip flexion ROM was significantly less with any level of decreased femoral version when compared to 15 degrees of anteversion (p&lt;0.04). When assessing the degree of change in ROM between the pre and post cam resection groups, femoral version had contrasting effects on hip ROM with decreasing femoral version showing less improvement with hip IR ROM while the femoral retroversion models had significantly greater improvements in hip flexion ROM when compared to the 15 degree anteversion model (p=0.02). Conclusions: Decreasing femoral version demonstrates significantly decreased hip flexion and IR both in the presence and absence of a CAM lesion. After cam femoroplasty, any level of decreased femoral version has a significant effect on hip ROM when compared to normal femoral anteversion. Femoral version seems to be the predominant factor in determining the amount of hip IR ROM in flexion and to a certain extent hip flexion ROM irrespective of the cam FAI. Impingement free ROM was only identified within the femoral anteversion models even after cam femoroplasty suggesting that recurrent mechanical collision and intra-articular damage seen with FAI may still be present with any level of true femoral retroversion. This study provides insight into the changes in hip flexion and IR that can be expected following both cam resection and simulated femoral osteotomy. [Table: see text]

Posterior Hip Impingement at Maximal Hip Extension in Female Patients With Increased Femoral Version or Increased McKibbin Index and Its Effect on Sports Performance.

The location of posterior hip impingement at maximal extension in patients with posterior femoroacetabular impingement (FAI) is unclear. To investigate the frequency and area of impingement at maximal hip extension and at 10° and 20° of extension in female patients with increased femoral version (FV) and posterior hip pain. Cross-sectional study; Level of evidence, 3. Osseous patient-specific 3-dimensional (3D) models were generated of 50 hips (37 female patients, 3D computed tomography) with a positive posterior impingement test and increased FV (defined as >35°). The McKibbin index (combined version) was calculated as the sum of FV and acetabular version (AV). Subgroups of patients with an increased McKibbin index >70° (24 hips) and FV >50° (20 hips) were analyzed. A control group of female participants (10 hips) had normal FV, normal AV, and no valgus deformity (neck-shaft angle, <139°). Validated 3D collision detection software was used for simulation of osseous impingement-free hip extension (no rotation). The mean impingement-free maximal hip extension was significantly lower in patients with FV >35° compared with the control group (15° ± 15° vs 55° ± 19°; P < .001). At maximal hip extension, 78% of patients with FV >35° had osseous posterior extra-articular ischiofemoral hip impingement. At 20° of extension, the frequency of posterior extra-articular ischiofemoral impingement was significantly higher for patients with a McKibbin index >70° (83%) and for patients with FV >35° (76%) than for controls (0%) (P < .001 for both). There was a significant correlation between maximal extension (no rotation) and FV (r = 0.46; P < .001) as well as between impingement area at 20° of extension (external rotation [ER], 0°) and McKibbin index (0.61; P < .001). Impingement area at 20° of extension (ER, 0°) was significantly larger for patients with McKibbin index >70° versus <70° (251 vs 44 mm2; P = .001). The limited hip extension found in our study could theoretically affect the performance of sports activities such as running, ballet dancing, or lunges. Therefore, although not examined directly in this study, these activities are not advisable for these patients. Preoperative evaluation of FV and the McKibbin index is important in female patients with posterior hip pain before hip preservation surgery (eg, hip arthroscopy).

Numerical study on the nonlinear characteristics of shock induced two-dimensional panel flutter in inviscid flow

For an aeroelastic system of a two-dimensional elastic panel subjected to an impinging inviscid oblique shockwave, the nonlinear flutter characteristics are affected by many factors such as shock impingement location, cavity pressure and initial perturbation. The effects of the above factors on the variation of system bifurcation type and dynamic behaviors are investigated numerically. A low-fidelity computational method coupled with local piston theory and van Karman plate model, and a high-fidelity computational method coupled with Euler equations and finite element model are used for fluid-structure interaction simulations. Two sets of new findings are unveiled. First, either the variation of shock impingement location or cavity pressure can induce the aeroelastic system to transition between a subcritical bifurcation and a supercritical bifurcation. For some cases, the system bifurcation characteristics exhibit strong sensitivity to these two factors. Second, it is found that in addition to the limit cycle oscillation in the form of a combination of the second and third structural modes, multiple stable limit cycle oscillations due to the coupling of higher-order modes can be triggered by suitable initial perturbations. These limit cycle oscillations have high frequencies and amplitudes, which suggest a higher risk of structural fatigue failure.

MRI 3D simulation of hip motion in female patients with and without ischiofemoral impingement.

To utilize hip MRI 3D models for demonstration of location and frequency of impingement during simulated range-of-motion in ischiofemoral impingement (IFI) compared to non-IFI hips. Sixteen hips (N = 7 IFI, 9 non-IFI) from 8 females were examined with high-resolution MRI. We performed image segmentation and generated 3D bone models and simulated hip range-of-motion and impingement. We examined the frequency and location of bone contact in early external rotation and early extension (0-20°), isolated maximum external rotation, and isolated maximum extension. Frequency and location of impingement at varied combinations of external rotation and extension and areas of simulated bone impingement at early external rotation and extension were compared between IFI and non-IFI. Higher frequency of bony impingement occurred more often in IFI hips at each simulated range-of-motion combination (P < 0.05). Impingement involved the lesser trochanter more often in IFI hips (P < 0.001) and occurred at early degrees of external rotation and extension. In isolated maximum external rotation, only the greater trochanter, intertrochanteric area, or both combined were involved, in 14%, 57%, and 29% in IFI hips. In isolated maximum extension, the lesser trochanter, intertrochanteric area, or both combined were involved in 71%, 14%, and 14% in IFI hips. The simulated area of bone impingement was significantly higher in IFI hips (P = 0.02). Hip MRI 3D models are feasible for simulated range-of-motion and show a higher frequency of extra-articular impingement at early stages of external rotation and extension in IFI compared to non-IFI hips.

Characterization of nonequilibrium shock interaction in CO2-N2 flows over double-wedges with respect to Mach number and geometry

The characterization of the shock interaction mechanism originating from the high-Mach nonequilibrium flow over double-wedge geometries is key to the design of hypersonic vehicles. The impact of changes in the freestream Mach number and double-wedge geometry on the patterns of shock interaction is investigated by means of numerical simulation in the case of CO2-N2 flows. The extended laminar Navier–Stokes equations with a two-temperature model to account for translational-to-vibrational internal energy transfer are considered the physical model of this type of flow. Simulations show that reducing the freestream Mach number leads to an increase in the separation region, both in the compression corner and in the locations of shock impingement. The impact of the size of the separation region on the patterns of interaction is such that it causes variations in the type of shock interaction. From the point of view of the flow physics near the wedges, decreasing the freestream Mach number has an equivalent effect to increasing the angle of the second wedge and an opposite effect to increasing the freestream temperature on the pattern of interaction. Results show that decreasing the freestream Mach number leads to an overall reduction in pressure and heating loads along the surface of the wedges and smaller regions of thermal equilibrium behind the bow shock.

Particle-based approach for modeling phase change and drop/wall impact at thermal spray conditions

The ability to model the impact of small, high-velocity molten droplets on a cold substrate is a key to designing the thermal spray process. Such drop-wall interactions involve complex phenomena such as droplet splash, phase change, solidification, substrate melting, and re-solidification. In this study, a particle-based approach, Smoothed Particle Hydrodynamics (SPH), is used to simulate the impact of molten droplets on a cold substrate at the thermal spray conditions, i.e., small droplets with high temperature and velocities. A solidification model is developed to simulate the phase change of the droplet and substrate. The present SPH method is validated against the experimental results using tin droplets on stainless-steel and aluminum substrate, and copper droplet on copper substrate. The histories of the spread factor, the substrate temperature, and the splat height at the impingement point are validated. A parametric study on the impact of the high melting point molybdenum droplets on different substrate materials (including tin, stainless steel, zinc, yttrium stabilized zirconia, and aluminum) is performed. The temporal evolution of the solidification interface and the height/thickness of the solidified splat are reported. Temperature distributions across the splat, substrate, and the corresponding melting/re-solidification are investigated. It was found that the cooling rates for the same impingement velocity were nearly the same for the same substrate regardless of the initial substrate temperatures but were higher for higher impingement velocities. A modified Biot number and thermal diffusivity were used to describe the heat transfer characteristics of the splat substrate combination. Results show that a higher modified Biot number for the splat-substrate combination indicates a faster cooling of the splat and a higher maximum substrate temperature. Furthermore, a substrate with high thermal diffusivity also causes heat at the impingement location to be dissipated quickly, resulting in a fast-cooling rate and accelerated solidification.

Large Hip Impingement Area and Subspine Hip Impingement in Patients With Absolute Femoral Retroversion or Decreased Combined Version.

It remains unclear if femoral retroversion is a contraindication for hip arthroscopy in patients with femoroacetabular impingement (FAI). To compare the area and location of hip impingement at maximal flexion and during the FADIR test (flexion, adduction, internal rotation) in FAI hips with femoral retroversion, hips with decreased combined version, and asymptomatic controls. Cross-sectional study; Level of evidence, 3. Twenty-four symptomatic patients (37 hips) with anterior FAI were evaluated. All patients had femoral version (FV) <5° according to the Murphy method. Two subgroups were analyzed: 13 hips with absolute femoral retroversion (FV <0°) and 29 hips with decreased combined version (McKibbin index <20°). All patients were symptomatic and had anterior groin pain and a positive anterior impingement test ; all had undergone pelvic computed tomography (CT) scans to measure FV. The asymptomatic control group consisted of 26 hips. Dynamic impingement simulation of maximal flexion and FADIR test at 90° of flexion was performed with patient-specific CT-based 3-dimensional models. Extra- or intra-articular hip impingement area and location were compared between the subgroups and with control hips using nonparametric tests. Impingement area was significantly larger for hips with decreased combined version (<20°) versus combined version (≥20°) (mean ± SD; 171 ± 140 vs 78 ± 55 mm2; P = .012) and was significantly larger for hips with FV <0° (absolute femoral retroversion) vs FV >0° (P = .025). Hips with absolute femoral retroversion had a significantly higher frequency of extra-articular subspine impingement versus controls (92% vs 0%; P < .001), compared to 84% of patients with decreased combined version. Intra-articular femoral impingement location was most often (95%) anterosuperior and anterior (2-3 o'clock). Anteroinferior femoral impingement location was significantly different at maximal flexion (anteroinferior [4-5 o'clock]) versus the FADIR test (anterosuperior and anterior [2-3 o'clock]) (P < .001). Patients with absolute femoral retroversion (FV <0°) had a larger hip impingement area, and most exhibited extra-articular subspine impingement. Preoperative FV assessment with advanced imaging (CT/magnetic resonance imaging) could help to identify these patients (without 3-dimensional modeling). Femoral impingement was located anteroinferiorly at maximal flexion and anterosuperiorly and anteriorly during the FADIR test.

Can a Computational Model Predict the Effect of Lesion Location on Cam-type Hip Impingement?

The Warwick consensus defined femoroacetabular impingement syndrome as a motion-related clinical disorder of the hip with a triad of symptoms, clinical signs, and imaging findings representing symptomatic premature contact between the proximal femur and acetabulum. Several factors appear to cause labral and cartilage damage, including joint shape and orientation and patient activities. There is a lack of tools to predict impingement patterns in a patient across activities. Current computational modeling tools either measure pure ROM of the joint or include complexity that reduces reliability and increases time to achieve a solution. The purpose of this study was to examine the efficacy of a low computational cost approach to combining cam-type hip shape and multiple hip motions for predicting impingement. Specifically, we sought to determine (1) the potential to distinguish impingement in individual hip shapes by analyzing the difference between a cam lesion at the anterior femoral neck and one located at the superior femoral neck; (2) sensitivity to three aspects of hip alignment, namely femoral neck-shaft angle, femoral version angle, and pelvic tilt; and (3) the difference in impingement measures between the individual activities in our hip motion dataset. A model of the shape and alignment of a cam-type impinging hip was created and used to describe two locations of a cam lesion on the femoral head-neck junction (superior and anterior) based on joint shape information available in prior studies. Sensitivity to hip alignment was assessed by varying three aspects from a baseline (typical alignment described in prior studies), namely, femoral neck-shaft angle, femoral version, and pelvic tilt. Hip movements were selected from an existing database of 18 volunteers performing 13 activities (10 male, eight female; mean age 44 ± 19 years). A subset was selected to maximize variation in the range of joint angles and maintain a consistent number of people performing each activity, which resulted in nine people per activity, including at least three of each sex. Activities included pivoting during walking, squatting, and golf swing. All selected hip motion cases were applied to each hip shape model. For the first part of the study, the number of motion cases in which impingement was predicted was recorded. Quantitative analyses of the depth of penetration of the cam lesion into the acetabular socket and qualitative observations of impingement location were made for each lesion location (anterior and superior). In the second part of the study, in which we aimed to test the sensitivity of the findings to hip joint orientation, full analysis of both cam lesion locations was repeated for three modified joint orientations. Finally, the results from the first part of the analysis were divided by activity to understand how the composition of the activity dataset affected the results. The two locations of cam lesion generated impingement in a different percentage of motion cases (anterior cam: 56% of motion cases; superior cam: 13% of motion cases) and different areas of impingement in the acetabulum, but there were qualitatively similar penetration depths (anterior cam: 6.8° ± 5.4°; superior cam: 7.9° ± 5.8°). The most substantial effects of changing the joint orientation were a lower femoral version angle for the anterior cam, which increased the percentage of motion cases generating impingement to 67%, and lower neck-shaft angle for the superior cam, which increased the percentage of motion cases generating impingement to 37%. Flexion-dominated activities (for example, squatting) only generated impingement with the anterior cam. The superior cam generated impingement during activities with high internal-external rotation of the joint (for example, the golf swing). This work demonstrated the capability of a simple, rapid computational tool to assess impingement of a specific cam-type hip shape (under 5 minutes for more than 100 motion cases). To our knowledge, this study is the first to do so for a large set of motion cases representing a range of activities affecting the hip, and could be used in planning surgical bone removal. The results of this study imply that patients with femoroacetabular impingement syndrome with cam lesions on the superior femoral head-neck junction may experience impinging during motions that are not strongly represented by current physical diagnostic tests. The use of this tool for surgical planning will require streamlined patient-specific hip shape extraction from imaging, model sensitivity testing, evaluation of the hip activity database, and validation of impingement predictions at an individual patient level.

Supersonic turbulent flows over sinusoidal rough walls

Direct numerical simulations were performed to characterize fully developed supersonic turbulent channel flows over isothermal rough walls. The effect of roughness was incorporated using a level-set/volume-of-fluid immersed boundary method. Turbulence statistics of five channel flows are compared, including one reference case with both walls smooth and four cases with smooth top walls and bottom walls with two-dimensional (2-D) and three-dimensional (3-D) sinusoidal roughnesses. Results reveal a strong dependence of the turbulence on the roughness topography and the associated shock patterns. Specifically, the 2-D geometries generate strong oblique shock waves that propagate across the channel and are reflected back to the rough-wall side. These strong shocks are absent in the smooth-wall channel and are significantly weaker in cases with 3-D roughness geometries, replaced by weak shocklets. At the impingement locations of the shocks on the top wall in the 2-D roughness cases, localized augmentations of turbulence shear production are observed. Such regions of augmented production also exist for the 3-D cases, at a much weaker level. The oblique shock waves are thought to be responsible for a more significant entropy generation for cases with 2-D surfaces than those with 3-D ones, leading to a higher irreversible heat generation and consequently higher temperature values in 2-D roughness cases. In the present supersonic channels, the effects of roughness extend beyond the near-wall layer due to the shocks. This suggests that outer layer similarity may not fully apply to a rough-wall supersonic turbulent flow.

Femoral impingement in maximal hip flexion is anterior-inferior distal to the cam deformity in femoroacetabular impingement patients with femoral retroversion

Femoroacetabular impingement (FAI) patients report exacerbation of hip pain in deep flexion. However, the exact impingement location in deep flexion is unknown. The aim was to investigate impingement-free maximal flexion, impingement location, and if cam deformity causes hip impingement in flexion in FAI patients. A retrospective study involving 24 patients (37 hips) with FAI and femoral retroversion (femoral version (FV) < 5° per Murphy method) was performed. All patients were symptomatic (mean age 28 years (SD 9)) and had anterior hip/groin pain and a positive anterior impingement test. Cam- and pincer-type subgroups were analyzed. Patients were compared to an asymptomatic control group (26 hips). All patients underwent pelvic CT scans to generate personalized CT-based 3D models and validated software for patient-specific impingement simulation (equidistant method). Mean impingement-free flexion of patients with mixed-type FAI (110° (SD 8°)) and patients with pincer-type FAI (112° (SD 8°)) was significantly (p < 0.001) lower compared to the control group (125° (SD 13°)). The frequency of extra-articular subspine impingement was significantly (p < 0.001) increased in patients with pincer-type FAI (57%) compared to cam-type FAI (22%) in 125° flexion. Bony impingement in maximal flexion was located anterior-inferior at femoral four and five o'clock position in patients with cam-type FAI (63% (10 of 16 hips) and 37% (6 of 10 hips)), and did not involve the cam deformity. The cam deformity did not cause impingement in maximal flexion. Femoral impingement in maximal flexion was located anterior-inferior distal to the cam deformity. This differs to previous studies, a finding which could be important for FAI patients in order to avoid exacerbation of hip pain in deep flexion (e.g. during squats) and for hip arthroscopy (hip-preservation surgery) for planning of bone resection. Hip impingement in flexion has implications for daily activities (e.g. putting on shoes), sports, and sex.Cite this article: Bone Joint Res2023;12(1):22-32.