Computational Hemodynamics-Based Growth Prediction for Small Abdominal Aortic Aneurysms: Laminar Simulations Versus Large Eddy Simulations.

Prior studies have shown that computational fluid dynamics (CFD) simulations help assess patient-specific hemodynamics in abdominal aortic aneurysms (AAAs); patient-specific hemodynamic stressors are frequently used to predict an AAA's growth. Previous studies have utilized both laminar and turbulent simulation models to simulate hemodynamics. However, the impact of different CFD simulation models on the predictive modeling of AAA growth remains unknown and is thus the knowledge gap that motivates this study. Specifically, CFD simulations were performed for 70 AAA models derived from 70 patients' computed tomography angiography (CTA) data with known growth status (i.e., fast-growing [> 5 mm/yr] or slowly growing [< 5 mm/yr]). We used laminar and large eddy simulation (LES) models to obtain hemodynamic parameters to predict AAAs' growth status. Predicting the growth status of AAAs was based on morphological, hemodynamic, and patient health parameters in conjunction with three classical machine learning (ML) classifiers, namely, support vector machine (SVM), K-nearest neighbor (KNN), and generalized linear model (GLM). Our preliminary results estimated aneurysmal flow stability and wall shear stress (WSS) were comparable in both laminar and LES flow simulations. Moreover, computed WSS and velocity-related hemodynamic variables obtained from the laminar and LES simulations showed comparable abilities in differentiating the growth status of AAAs. More importantly, the predictive modeling performance of the three ML classifiers mentioned above was similar, with less than a 2% difference observed (p-value > 0.05). In closing, our findings suggest that two different flow simulations investigated did not significantly affect outcomes of computational hemodynamics and predictive modeling of AAAs' growth status, given the data investigated.

Large-eddy simulation of aerosol concentrations in a realistic urban environment: Model validation and transport mechanism.

Air pollution in urban environments exhibits large spatial and temporal variations due to high heterogeneous air flow and emissions. To address the complexity of local air pollutant dynamics, a comprehensive large-eddy simulation using the PALM model system v6.0 was conducted. The distribution of flow and vehicle emitted aerosol particles in a realistic urban environment in Malmö, Sweden, was studied and evaluated against on-site measurements made using portable instrumentation on a spring morning in 2021. The canyon transport mechanisms were investigated, and the convective and turbulent mass-transport rates compared to clarify their role in aerosol transport. The horizontal distribution of aerosols showed acceptable evaluation metrics for both mass and number. Flow and pollutant concentrations were more complex than those in idealized street canyon networks. Vertical turbulent mass-transport rate was found to dominate the mass transport process compared with the convective transport rate, contributing more than 70% of the pollutant transport process. Our findings highlight the necessity of examining various aerosol metric due their distinct dispersion behaviour. This study introduces a comprehensive high-resolution modelling framework that accounts for dynamic meteorological and aerosol background boundary conditions, real-time traffic emission, and detailed building features, offering a robust toll for local urban air quality assessment.

Assessment of longitudinal dispersion and flow resistance near floodplains through riparian woodlands.

Riparian woodlands prevent bank erosions, recycle minerals, sustain biodiversity, act as flow resistance on floodplains, and filter pollutants. The emergent trees characterize woodlands with different spacing arrangements that dictate flow resistance and longitudinal dispersion of the pollutants in compound channel flow. The single- and multistage compound channels exist in urban and natural watercourses with riparian and transplanted trees on different stages of the floodplain. This study numerically validates the planting of vegetation in lines on single- and multistage floodplains using a wall-modeled large-eddy simulation model. Post-validation, the focus of the study was to assess the hydrodynamic behavior and mixing around the floodplain and main channel section of different tested configurations. The approximation of flow structures for the various configurations of tree plantations shows stronger vortices with significant characteristic length scales for floodplains closer to the main channel. The intensity of the secondary current is higher for denser planted trees at junctions of floodplains. For higher flow events, drag force contributions for staged floodplains with trees on both stages are 45-41%, and trees on the top stage contribute 27-22% to the total frictional force budget. The subsequent investigation shows that the in-line trees geometrical configuration and spacing arrangement on the floodplain dictates flow resistance and longitudinal dispersion of the pollutants and contamination in channel flow. The results show that the overall reduction in discharge for floodplains with tree planting is 19.8-36.2% for single-stage and 10.4-23.6% for multistage compound channels. The longitudinal dispersion coefficients for each multi-zone model predict a 61% and 41% dispersion reduction, respectively, in single- and multistage floodplains with planted trees. Floodplains with denser tree spacing have a maximum zonal discharge reduction of 45% for a single-stage and 27.2% and 28.0% for multistage channels. These findings strongly suggest that the planting parameters of spacing-to-diameter ratio and floodplain geometry play a pivotal role in floodplain management from the perspective of contaminant dispersion and flood risk reduction during high-flow events.

Turbulence-induced droplet grouping and augmented rain formation in cumulus clouds.

This paper provides the first observational analysis of how droplet separation is impacted by the flinging action of microscale vortices in turbulent clouds over a select radii range and how they vary over cloud cores and along the peripheral edges. It is premised that this mechanism initiates droplet separation within a cloud volume soon after condensational growth, largely in the cloud core, and operates until the cloud droplet radii exceed 20-30 µm when this effect fades rapidly. New observations are presented showing how microscale vortices also impact the settling rates of droplets over a critical size range (6-18 µm) causing them to sediment faster than in still air affecting swept volumes and thereby impacting the rain initiation and formation. Large-scale atmospheric models ignore these microscale effects linked to rapid droplet growth during the early stages of cloud conversion. Previous studies on droplet spatial organization along the cloud edges and inside the deep core have shown that homogeneous Poisson statistics, indicative of the presence of a vigorous in-cloud mixing process at small scales obtained, in contrast to an inhomogeneous distribution along the edges. In this paper, it is established that this marked core region, homogeneity can be linked to microscale vortical activity which flings cloud droplets in the range of 6-18 µm outward. The typical radius of the droplet trajectories or the droplet flung radii around the vortices correlates with the interparticle distance strongly. The correlation starts to diminish as one proceeds from the central core to the cloud fringes because of the added entrainment of cloud-free air. These first results imply that droplet growth in the core is first augmented with this small-scale interaction prior to other more large-scale processes involving entrainment mixing. This first study, combining these amplified velocities are included in a Weather Research and Forecasting- LES case study. Not only are significant differences observed in the cloud morphology when compared to a baseline case, but the 'enhanced' case also shows early commencement of rainfall along with intense precipitation activity compared to the 'standard' baseline case. It is also shown that the modelled equilibrium raindrop spectrum agrees better with observations when the enhanced droplet sedimentation rates mediated by microscale vortices are included in the calculations compared to the case where only still-air terminal velocities are used.

Velocity field and turbulence structure of the meandering flow produced by alternating deflectors.

Meandering flow can be formed during the advance of natural rivers by the scouring of river banks. However, this phenomenon is not common in artificial cement channels. This study used experimental scouring terrain data for a numerical simulation to study the meandering flow pattern formed between double alternating deflectors in a straight channel. The numerical results showed that the path of the accelerated flow generated by the upstream deflector was changed by installing a downstream deflector while the flow rate remained unchanged. Thus, a meandering flow formed, and a stable, narrow, high-speed zone formed in the downstream area. The accelerated flow between the two deflectors hit the channel bank soon after its direction changed. Then, a strong downward flow formed in this area, which may have produced an elliptical scour hole. A large-scale vortex structure was formed in the elliptical scour hole, which was influenced by the horseshoe vortex system before the downstream deflector.

A three-dimensional model of terrain-induced updrafts for movement ecology studies.

BACKGROUND: Spatially explicit simulation models of animal movements through the atmosphere necessarily require a representation of the spatial and temporal variation of atmospheric conditions. In particular, for movements of soaring birds that rely extensively on vertical updrafts to avoid flapping flight, accurate and reliable estimation of the vertical component of wind is critical. The interaction between wind and complex terrain shapes both the horizontal and vertical wind fields, highlighting the need to model the coupling between local terrain features and atmospheric conditions at scales relevant to animal movement. METHODS: In this work, we propose a new empirical model for estimating the orographic updraft field. The model is developed using computational fluid dynamics simulations of canonical atmospheric conditions over moderately complex terrain. To isolate buoyancy and thermal effects, and focus on terrain-induced effects, we use only simulations of a neutrally stratified atmosphere to develop the model. The model, which we name Engineering Vertical Velocity Estimator (EVVE), is simple to implement and is a function of the underlying terrain elevation map, the desired height above ground level (AGL), and wind conditions at a reference height (80 m). We validate the model with data from the Alaiz mountain (Spain) field campaign. RESULTS: Compared to observations, the proposed improved model estimates the updrafts at 120 m AGL with a mean error of 0.11 m/s ( σ = 0.28 m/s), compared to 0.85 m/s ( σ = 0.58 m/s) for its baseline. For typical land-based wind turbine hub heights of 80 m AGL, the proposed model has a mean error of 0.04 m/s ( σ = 0.25 m/s), compared to baseline 0.54 m/s ( σ = 0.45 m/s) estimations. We illustrate an application of the model in movement ecology by comparing simulated tracks and presence maps of golden eagles (Aquila chrysaetos) moving across two distinct landscapes. The tracks and presence maps are obtained using a simple heuristic-based movement model, with the updraft field given by the proposed model and a wind vector-based estimation approach that is currently in wide use in movement ecology studies of raptors and other soaring birds. CONCLUSIONS: We highlight that movement model results can be sensitive to the underlying orographic updraft model, especially in studies of fine-scale movements in regions of complex topography. We suggest adopting the proposed model rather than the wind vector estimation method for studies of soaring bird movements.

Computational Rhinology: Unraveling Discrepancies between In Silico and In Vivo Nasal Airflow Assessments for Enhanced Clinical Decision Support.

Computational rhinology is a specialized branch of biomechanics leveraging engineering techniques for mathematical modelling and simulation to complement the medical field of rhinology. Computational rhinology has already contributed significantly to advancing our understanding of the nasal function, including airflow patterns, mucosal cooling, particle deposition, and drug delivery, and is foreseen as a crucial element in, e.g., the development of virtual surgery as a clinical, patient-specific decision support tool. The current paper delves into the field of computational rhinology from a nasal airflow perspective, highlighting the use of computational fluid dynamics to enhance diagnostics and treatment of breathing disorders. This paper consists of three distinct parts-an introduction to and review of the field of computational rhinology, a review of the published literature on in vitro and in silico studies of nasal airflow, and the presentation and analysis of previously unpublished high-fidelity CFD simulation data of in silico rhinomanometry. While the two first parts of this paper summarize the current status and challenges in the application of computational tools in rhinology, the last part addresses the gross disagreement commonly observed when comparing in silico and in vivo rhinomanometry results. It is concluded that this discrepancy cannot readily be explained by CFD model deficiencies caused by poor choice of turbulence model, insufficient spatial or temporal resolution, or neglecting transient effects. Hence, alternative explanations such as nasal cavity compliance or drag effects due to nasal hair should be investigated.

Validation of numerical simulations and experiments on impulse characteristics induced by self-excited oscillation.

The high-frequency pulse flow, equivalent to the natural frequency of rocks, is generated by a self-excited oscillating cavity to achieve resonance rock-breaking. The flow field and oscillating mechanism of the self-excited oscillating cavity were simulated using the large eddy simulation method of Computational Fluid Dynamics (CFD). A field-scale testing apparatus was developed to investigate the impulse characteristics and verify the simulation results. The results show that the fluid at the outlet at the tool is deflected due to the pulse oscillation of the fluid. The size and shape of low-pressure vortices constantly change, leading to periodic changes in fluid impedance within the oscillating cavity. The impulse frequency reaches its highest point when the length-diameter ratio is 0.67. As the length-diameter ratio increases, the tool pressure loss also increases. Regarding the cavity thickness, the impulse frequency of the oscillating cavity initially decreases, then increases, and finally decreases again. Moreover, both the impulse frequency and pressure loss increase with an increase in displacement. The numerical simulation findings align with the experimental results, thus confirming the validity of the theoretical model. This research provides theoretical guidance for the practical application of resonance rock-breaking technology.

Aortic valve neocuspidization and bioprosthetic valves: Evaluating turbulence haemodynamics.

Aortic valve disease is often treated with bioprosthetic valves. An alternative treatment is aortic valve neocuspidization which is a relatively new reparative procedure whereby the three aortic cusps are replaced with patient pericardium or bovine tissues. Recent research indicates that aortic blood flow is disturbed, and turbulence effects have yet to be evaluated in either bioprosthetic or aortic valve neocuspidization valve types in patient-specific settings. The aim of this study is to better understand turbulence production in the aorta and evaluate its effects on laminar and turbulent wall shear stress. Four patients with aortic valve disease were treated with either bioprosthetic valves (n=2) or aortic valve neocuspidization valvular repair (n=2). Aortic geometries were segmented from magnetic resonance images (MRI), and 4D flow MRI was used to derive physiological inlet and outlet boundary conditions. Pulsatile large-eddy simulations were performed to capture the full range of laminar, transitional and turbulence characteristics in the aorta. Turbulence was produced in all aortas with highest levels occurring during systolic deceleration. In the ascending aorta, turbulence production is attributed to a combination of valvular skew, valvular eccentricity, and ascending aortic dilation. In the proximal descending thoracic aorta, turbulence production is dependent on the type of arch-descending aorta connection (e.g., a narrowing or sharp bend) which induces flow separation. Laminar and turbulent wall shear stresses are of similar magnitude throughout late systolic deceleration and diastole, although turbulent wall shear stress magnitudes exceed laminar wall shear stresses between 27.3% and 61.1% of the cardiac cycle. This emphasises the significance of including turbulent wall shear stress to improve our comprehension of progressive arterial wall diseases. The findings of this study recommend that aortic valve treatments should prioritise minimising valvular eccentricity and skew in order to mitigate turbulence generation.

Control mechanisms of different bionic structures for hydrofoil cavitation.

Cavitation limits the efficient and stable operation of rotating machinery. The exploration of control methods for hydrofoil cavitation is important for improving the performance of hydraulic machinery. The leading-edge protuberances of the humpback flipper and the spine structure of the tail fin of sailfish are two common bionic structures for cavitation control; however, the control effects of both have limitations. Accordingly, in this study, a passive control method for hydrofoil cavitation was developed by combining the two bionic structures. With the large eddy simulation method, the cavitation processes of wavy leading-edge hydrofoil, bionic fin spine structure hydrofoil, and novel bionic combined structure hydrofoil were studied under a cavitation number of σ = 0.8. The control mechanisms of the three bionic structures for the hydrofoil cavitation were investigated. The results indicated that the novel bionic combined hydrofoil realised the superposition and complementation of the control effects of the two single bionic structures and achieved a better cavitation inhibition effect, reducing the total volume of cavitation by 43 %. In addition, it enhanced the stability of the flow field and reduced the standard deviation of the pressure coefficient on the suction surface by up to 46.55 %. This research provides theoretical support for the optimisation and modification of the blades of hydraulic machinery, such as propellers and pump turbines.

Study on the flow mechanism and frequency characteristics of rales in lower respiratory tract.

The frequency characteristics of lung sounds have great significance for noninvasive diagnosis of respiratory diseases. The rales in the lower respiratory tract region that can provide rich information about symptoms of respiratory diseases are not clear. In this paper, a three-dimensional idealized bifurcated lower respiratory tract geometric model, which contains 3rd to 13th generation (G3-G13) bronchi is constructed, where Re ∼ 10 1 - 10 3 , and then the large eddy simulation and volume of fluid are used to study the fluid flow characteristics. Ffowcs Williams and Hawkings model are subsequently used to study the frequency characteristics of rale of different generations of bronchi. The results showed that bronchial blockage and sputum movement will enhance the turbulence intensity and vortex shedding intensity of flow. The dominant frequency and highest value of sound pressure level (SPL) of rhonchi/moist crackles decrease with the increase of bronchial generation. The change rates of dominant frequency of rhonchi / moist crackles in adjacent generations were 5.0 ± 0.1 ~ 9.1 ± 0.2% and 3.1 ± 0.1 ~ 11.9 ± 0.3%, respectively, which is concentrated in 290 ~ 420 Hz and 200 ~ 300 Hz, respectively. The change rates of SPL of rhonchi/moist crackles were 8.8 ± 0.1 ~ 15.7 ± 0.1% and 7.1 ± 0.1 ~ 19.5 ± 0.2%, respectively, which is concentrated in 28 ~ 50 dB and 16 ~ 32 dB, respectively. In the same generation of bronchus (e.g., G8, G9) with the same degree of initial blockage, the dominant frequency and SPL of moist crackles can be 3.7 ± 0.2% and 4.5 ± 0.3% slightly higher than that of rhonchi, respectively. This research is conducive to the establishment of a rapid and accurate noninvasive diagnosis system for respiratory diseases.

Numerical study on particle behavior in a Y-junction mixer for supercritical water hydrolysis.

In the continuous-type supercritical water hydrolysis process, rapid mixing of supercritical and subcritical streams is important to maximize yield and minimize degradation from over-reaction. This work investigated the particle behavior in a Y-junction mixer using large eddy simulation coupled with a discrete phase model, aiming to optimize the supercritical hydrolysis process for biomass conversion. A series of numerical simulations analyzed the influence of the mixer's orientation, flow directions, and flow rates on effective mixing and backflow prevention. The results demonstrated that the most effective mixing occurred in a vertically oriented Y-junction mixer with an upward-directed supercritical water inlet, aligning the momentum direction of natural and forced convection effectively. Consequently, over 80% of particles reached the temperatures close to the mixing temperature of supercritical and subcritical water within the Y-junction mixing zone, indicating enhanced mixing effectiveness and potential for efficient hydrolysis. This configuration also minimized backflow.

Assessing airflow unsteadiness in the human respiratory tract under different expiration conditions.

To enhance the understanding of airflow characteristics in the human respiratory system, the expiratory airflow in a human respiratory tract model was simulated using large eddy simulation and dynamic mesh under different expiration conditions aligned with clinically measured data. The airflow unsteadiness was quantitatively assessed using power spectral density (PSD) and spectral entropy (SE). The following findings were obtained: (1) The airflow is highly turbulent in the mouth-pharynx region during expiration, with its dynamic characteristics being influenced by both the transient expiration flow pattern at mouth piece and the glottis motion. (2) PSD analysis reveals that the expiratory airflow is very unsteady, exhibiting a broad-band attenuation spectrum in the pharynx-trachea region. When only transient expiration or glottis motion is considered, the PSD spectrum changes slightly. When both are ignored, however, the change is significant, with the peak frequency reduced to 10% of the real expiration condition. (3) SE analysis indicates that the airflow transitions into turbulence in the trachea, and there may be multiple transitions in the region of soft palate. The transient expiration or glottis motion alone increases turbulence intensity by 2%-15%, while ignoring both reduces turbulence intensity by 10%-20%. This study implies that turbulence characteristics can be significantly different under different expiratory conditions, and therefore it is necessary to determine the expiratory flow characteristics using clinically measured expiratory data.

Investigating the Hydrodynamics of Intravenous Drug Infusions.

One of the most common reported adverse events for intravenous (IV) infusions are infusion site reactions, ranging from redness and pain at the site of infusion to thrombophlebitis.  The connection between drug infusion and what drives these adverse events is not well understood. To aid in understanding these phenomena, it is crucial to accurately characterize the evolving hemodynamic environment of the infusion site when developing new intravenous formulations, as too rapid dilution may cause precipitation in the vein, while too little dilution might contribute to phlebitis. In this study, a Large-Eddy Simulation (LES) turbulence modeling inside a Computational Fluid Dynamics (CFD) framework has been used to simulate the flow and mixing characteristics of an infusion entering the bloodstream. This work represents the first such study reporting transient flow fields for intravenous infusions using LES CFD simulations with a realistic non-Newtonian blood model. The output of the CFD model closely resembled the flow and mixing patterns generated in benchtop tests for infusions into a blood analogue and water as the venous fluid across a wide range of flow rates.  These models were then investigated further to compare how changes to the fluid rheology model, needle orientation and needle position within the vein resulted in altered mixing regimes at different flow rates.

Trailing-edge fringes enable robust aerodynamic force production and noise suppression in an owl wing model.

As one of the unique owl-wing morphologies, trailing-edge (TE) fringes are believed to play a critical role in the silent flight of owls and have been widely investigated using idealized single/tandem airfoils. However, the effect of TE fringes and associated mechanisms on the aeroacoustics of owl wings, which feature curved leading edges, wavy TEs, and several feather slots at the wingtips, have not yet been addressed. In this study, we constructed two 3D owl wing models, one with and one without TE fringes, based on the geometric characteristics of a real owl wing. Large-eddy simulations and the Ffowcs Williams‒Hawkings analogy were combined to resolve the aeroacoustic characteristics of the wing models. Comparisons of the computed aerodynamic forces and far-field acoustic pressure levels demonstrate that the fringes on owl wings can robustly suppress aerodynamic noise while sustaining aerodynamic performance comparable to that of a clean wing. By visualizing the near-field flow dynamics in terms of flow and vortex structures as well as flow fluctuations, the mechanisms of TE fringes in owl wing models are revealed. First, the TE fringes on owl wings are reconfirmed to robustly suppress flow fluctuations near the TE by breaking up large TE vortices. Second, the fringes are observed to effectively suppress the shedding of wingtip vortices by mitigating the flow interaction between feathers (feather-slot interaction). These complementary mechanisms synergize to enhance the robustness and effectiveness of the TE fringe effects in owl wing models, in terms of aerodynamic force production and noise suppression. This study thus deepens our understanding of the role of TE fringes in real owl flight gliding and points to the validity and feasibility of employing owl-inspired TE fringes in practical applications of low-noise fluid machinery.

Assessing turbulent effects in ascending aorta in presence of bicuspid aortic valve.

Aortic valves with bicuspids have two rather than three leaflets, which is a congenital heart condition. About 0.5-2% of people have a bicuspid aortic valve. Blood flow through the aorta is commonly believed to be laminar, although aortic valve disorders can cause turbulent transitions. Understanding the impact of turbulence is crucial for foreseeing how the disease will progress. The study's objective was use large eddy simulation to provide a thorough analysis of the turbulence in bicuspid aortic valve dysfunction. Using a large eddy simulation, the blood flow patterns of the bicuspid and tricuspid aortic valves were compared, and significant discrepancies were found. The velocity field in flow in bicuspid configurations was asymmetrically distributed toward the ascending aorta. In tricuspid aortic valve (TAV) the flow, on the other hand, was symmetrical within the same aortic segment. Moreover, we looked into standard deviation, Q-criterion, viscosity ratio and wall shear stresses for each cases to understand transition to turbulence. Our findings indicate that in the bicuspid aortic valve (BAV) case, the fluid-dynamic abnormalities increase. The global turbulent kinetic energy and time-averaged wall shear stress for the TAV and BAV scenarios were also examined. We discovered that the global turbulent kinetic energy was higher in the BAV case compared to TAV, in addition to the increased wall shear stress induced by the BAV in the ascending aorta.

Examining the effects of soil entrainment during nuclear cloud rise on fallout predictions using a multiscale atmospheric modeling framework.

Current operational models for nuclear cloud rise over land were developed and validated using observations from shallow-buried or surface detonations, where lofted soil quickly mixed with fission products from the detonation. These models poorly predict fallout from elevated detonations near the fallout-free height of burst (FFHOB), where interactions with the ground are limited and the mixing of fission products and lofted soil is incomplete. Fallout-free is a misnomer at this HOB, as fallout was observed in these cases, but was below the levels of concern, especially off-grounds of the nuclear test site. To correctly characterize and model fallout from detonations near the FFHOB, models must be developed which can capture the stratified nature of the particle and activity-size distributions within the cloud. Previously, it was shown that the Weather Research and Forecasting (WRF) model can accurately simulate nuclear cloud rise for airbursts with little to no ground interactions (Arthur et al., 2021). That work is expanded here by (1) using a radiation-hydrodynamics code to improve the fireball initialization in WRF, (2) further developing an aerosol package from WRF-Chem to simulate lofted soil, and (3) combining the WRF cloud rise simulations with the operational models used at the National Atmospheric Release Advisory Center (NARAC) for fallout modeling. Using this combination of codes, the Upshot-Knothole Grable detonation, which was just below the FFHOB, is simulated from seconds after detonation through cloud rise and fallout, and results are compared to historical test data. The results show improved prediction of dose rate and highlight the need to correctly characterize the entrainment of material into the cloud and the subsequent mixing of fission products with entrained material.

[Experimental and numerical investigation of fluid-particle-interactions in water]. / Experimentelle und numerische Untersuchung von Fluid-Partikel-Interaktionen in Wasser.

For the development of improved sediment transport models, the basic understanding of the interaction between the solid particle and the moving fluid (water) is important. In this article, current developments in the field of fluid-particle interaction are presented based on two research articles by Gold et al. (2023) and Worf et al. (2022). One presented in this article uses state of the art measurement methods to investigate the flow around spheres of different densities that oscillate in initially resting body of water. For the spherical pendulum a similar vortex shedding characteristic was observed for all investigated fluid density ratios (m*=ρS/ρF=1.14,14.95, density ratio between solid and fluid). A new object tracking method (DOT) is also presented, which enables temporally and spatially resolved analysis of flow structures in the fluid field. The experimental results of Gold et al. (2023) show, that vortex shedding occurs during the first period. This vortex propagates downward and eventually dissipates. Furthermore, a damping optimum of the spherical pendulum in the range of m*=2.50 was observed. Additionally, an experiment with a cylindrical pendulum with m∗=4.98 was investigated numerically utilizing an immersed boundary method. The process of creation and separation up to the dissipation of a vortex ring was described. Furthermore, this investigation by Worf et al. (2022) described the creation of tip vortices. These were connected with the development of the three-dimensional flow and added mass coefficient.

Lagrangian large eddy simulations via physics-informed machine learning.

High-Reynolds number homogeneous isotropic turbulence (HIT) is fully described within the Navier-Stokes (NS) equations, which are notoriously difficult to solve numerically. Engineers, interested primarily in describing turbulence at a reduced range of resolved scales, have designed heuristics, known as large eddy simulation (LES). LES is described in terms of the temporally evolving Eulerian velocity field defined over a spatial grid with the mean-spacing correspondent to the resolved scale. This classic Eulerian LES depends on assumptions about effects of subgrid scales on the resolved scales. Here, we take an alternative approach and design LES heuristics stated in terms of Lagrangian particles moving with the flow. Our Lagrangian LES, thus L-LES, is described by equations generalizing the weakly compressible smoothed particle hydrodynamics formulation with extended parametric and functional freedom, which is then resolved via Machine Learning training on Lagrangian data from direct numerical simulations of the NS equations. The L-LES model includes physics-informed parameterization and functional form, by combining physics-based parameters and physics-inspired Neural Networks to describe the evolution of turbulence within the resolved range of scales. The subgrid-scale contributions are modeled separately with physical constraints to account for the effects from unresolved scales. We build the resulting model under the differentiable programming framework to facilitate efficient training. We experiment with loss functions of different types, including physics-informed ones accounting for statistics of Lagrangian particles. We show that our L-LES model is capable of reproducing Eulerian and unique Lagrangian turbulence structures and statistics over a range of turbulent Mach numbers.

Computational Fluid Dynamics Turbulence Model and Experimental Study for a Fontan Cavopulmonary Assist Device.

Head-flow HQ curves for a Fontan cavopulmonary assist device (CPAD) were measured using a blood surrogate in a mock circulatory loop and simulated with various computational fluid dynamics (CFD) models. The tests benchmarked the CFD tools for further enhancement of the CPAD design. Recommended Reynolds-Averaged Navier-Stokes (RANS) CFD approaches for the development of conventional ventricular assist devices (VAD) were found to have shortcomings when applied to the Fontan CPAD, which is designed to neutralize off-condition obstruction risks that could contribute to a major adverse event. The no-obstruction condition is achieved with a von Karman pump, utilizing large clearances and small blade heights, which challenge conventional VAD RANS-based CFD hemodynamic simulations. High-fidelity large eddy simulation (LES) is always recommended; however, this may be cost-inhibitive for optimization studies in commercial settings, thus the reliance on RANS models. This study compares head and power predictions of various RANS turbulence models, employing experimental measurements and LES results as a basis for comparison. The models include standard k-ϵ, re-normalization group k-ϵ, realizable k-ϵ, shear stress transport (SST) k-ω, SST with transitional turbulence, and Generalized k-ω. For the pressure head predictions, it was observed that the standard k-ϵ model provided far better agreement with experiment. For the rotor torque, k-ϵ predictions were 30% lower than LES, while the SST and LES torque values were near identical. For the Fontan CPAD, the findings support using LES for the final design simulations, k-ϵ model for head and general flow simulation, and SST for power, shear stress, hemolysis, and thrombogenicity predictions.