Background: Understanding the transport of nanoparticles within blood vessels and their distribution in tumor tissues is crucial for the successful implementation of nanotechnological strategies in clinical practice. Although numerous studies have examined nanoparticle transport in blood flow, none have comprehensively investigated all the sequential steps a nanoparticle must undergo prior to internalization by target cells. Methods: A computational framework was developed in COMSOL Multiphysics to simulate nanoparticle (NP) transport from systemic administration through to tumor cell internalization. The model integrates three coupled stages: (1) NP movement within a non-Newtonian blood flow; (2) trans-endothelial transport; and (3) NP motion within the tumor stroma, incorporating affinity forces to capture ligand–receptor interactions. The tumor geometry was reconstructed, including cancer cells and fibroblasts, to reproduce physiological porosity. Multiple case studies were conducted to evaluate the impact of particle density, injection velocity, and size on NP biodistribution. Results: The computational model effectively simulates nanoparticle transport across all stages. Notably, it is the first model in the literature to incorporate the affinity of functionalized nanoparticles, which facilitates ligand–receptor interactions for targeted delivery. Simulation outcomes indicate that a low Stokes number is critical for ensuring a higher percentage of particles reach the end of the capillary network. Furthermore, surface modification of nanoparticles with ligands promotes more specific distribution within the stroma, reducing the percentage of nanoparticles that fail to reach target cells by approximately 50% Conclusions: A novel and comprehensive computational model has been developed to include the entire process of nanoparticle distribution following systemic administration, including specific recognition by cellular receptors.

Magnetic hyperthermia is a promising cancer treatment method that involves complex multiphysics phenomena, including interstitial tissue fluid flow, magnetic nanoparticle (MNP) transport, and temperature evolution. However, these intricate processes have rarely been studied simultaneously, primarily due to the lack of a comprehensive simulation tool. To address this issue, we develop a comprehensive numerical framework in this study. Using this framework, we simulate a circular-shaped tumor embedded in healthy tissue. The treatment process is examined under two scenarios: one considering gravity and the other neglecting it. Without gravity, the interstitial tissue flow remains stationary, and hence MNP transport and temperature evolution are determined solely by diffusion. The optimal treatment time, when the tumor cells are completely ablated, decreases with both the Lewis number and the heat source number, following a power law. When gravity is considered, treatment efficacy deteriorates due to buoyancy-induced MNP movement, significantly extending the time required to completely ablate the tumor cells. This required time increases with both the buoyancy ratio and the Darcy ratio, also following a power law. The results from this study could provide valuable guidelines for practical magnetic hyperthermia treatment.

We report recent advances in the development of a computationally efficient and accurate timedependent RANS-based (T-RANS) method for predicting heavy gas turbulent dispersion in complex urban areas. We analyze two cases: (i) the simplified Kit Fox setup at the field scale, and (ii) a complex urban area (the city of Beijing). The present approach, which models the two-way coupled transport of turbulent airflow and heavy gas concentration, combined with a passive-element method for representing obstacles and buildings, has proven to be both numerically efficient and accurate in predicting the temporal evolution of local concentrations.

We present recent results from large eddy simulations (LES) of bounded, turbulent, double-diffusive convection in seawater, focusing on small-aspect-ratio domains (2:2:1) across an extensive parameter range. We focus on presenting the main instantaneous and time-averaged features of flow, heat, and concentration fields. In addition to analyzing the power spectral density (PSD) and probability density function (PDF) at characteristic locations within the simulated domain, special attention is devoted to a detailed analysis of the budgets of the governing second-order correlations for different strengths of the imposed stable thermal stratification gradient.

We report on numerical studies of bounded double-diffusive turbulent convection, which involves the combined effects of concentration/solutal and thermal buoyancy forces. Our study focuses on an intermediate range of the characteristic non-dimensional numbers, specifically 107≤Rac≤109, and 0≤Raθ≤106. We use fixed values for the concentration and temperature Prandtl numbers (i.e. Prc=700, Prθ = 7), which approximately correspond to seawater properties. We apply wall-resolved Large Eddy Simulations (LES) and compare the obtained results with available Direct Numerical Simulations (DNS) in the literature. Our findings show an overall good agreement in predicting the global wall mass and heat transfer coefficients, achieved with significantly reduced computational costs. Furthermore, the local mass and heat transfer distributions reveal a high sensitivity to the strength of the vertically imposed stable thermal stratification. Finally, we present the vertical profiles of the long-term time-averaged first and second moments.

<Note to JM: Errors in article heads, particularly titles, author names, affiliations must not be repeated. In such cases, use the correct title, authors, affiliation(s), etc. If the corrigendum is the result of authorship changes, include the correct author details in the article head.>

The authors regret that the final version of the Graphical Abstract was not included with the original manuscript submission. The Graphical Abstract is now available.

We apologize for any inconvenience this may have caused.

A novel approach to generate left ventricular (LV) pressure–volume (PV) loops from combined 4D Flow MRI and computational fluid dynamics (CFD) is presented. Pressure was calculated from person-specific three-dimensional (3D) CFD models created from LV segmentations and peak-systolic pressure from the one-dimensional 111-artery CFD model, with aortic flow from 4D Flow MRI as input. Ten healthy volunteers underwent scan–rescan MRI. Additionally, one patient without cardiovascular abnormalities underwent MRI and invasive catheter measurement for single-case comparison. Scan–rescan reproducibility was very good overall, with no significant differences in any parameters and ICCs of all parameters but minimum pressure were significant and high (0.55–0.99). Aortic flow of 3D CFD model correlated well with 4D Flow (ICC = 0.74) and stroke volume of LV segmentation (ICC = 0.90). Segmentation volume variability resulted in 12% difference in stroke work and mean external power, while aortic flow variability resulted in 10–11% difference in most parameters. Single-case comparison is promising, with only 1.8 mmHg and 0.005 mmHg/mL difference in minimum pressure and EDPVR, and <10% differences for other parameters. Noninvasive pressure–volume loops can therefore reproducibly be generated from only aortic flow, cine short axis MRI, and brachial pressure measurement. Single-case comparison shows promise, but larger validation studies are needed.

In this work, we developed deep neural networks for the fast and comprehensive estimation of the most salient features of aortic blood flow. These features include velocity magnitude and direction, 3D pressure, and wall shear stress. Starting from 40 subject-specific aortic geometries obtained from 4D Flow MRI, we applied statistical shape modeling to generate 1,000 synthetic aorta geometries. Complete computational fluid dynamics (CFD) simulations of these geometries were performed to obtain ground-truth values. We then trained deep neural networks for each characteristic flow feature using 900 randomly selected aorta geometries. Testing on remaining 100 geometries resulted in average errors of 3.11% for velocity and 4.48% for pressure. For wall shear stress predictions, we applied two approaches: (i) directly derived from the neural network-predicted velocity, and, (ii) predicted from a separate neural network. Both approaches yielded similar accuracy, with average error of 4.8 and 4.7% compared to complete 3D CFD results, respectively. We recommend the second approach for potential clinical use due to its significantly simplified workflow. In conclusion, this proof-of-concept analysis demonstrates the numerical robustness, rapid calculation speed (less than seconds), and good accuracy of the CFD-based machine learning approach in predicting velocity, pressure, and wall shear stress distributions in subject-specific aortic flows.

Targeted drug delivery to the deep lung improves therapeutic outcomes, but respiratory system variability complicates drug spray design. Numerical simulations offer insights for individualized treatments but are computationally intensive, highlighting the need for surrogate models for real-time deposition prediction. This study comprehensively explores the multi-task predictive capability of regression models, including Linear regression (LR), Bayesian regression (BR), Back Propagation Neural Network (BPNN), Support Vector Machine (SVM), Random Forest (RF), XGBoost, and CatBoost, for predicting total and regional deposition rates of inhaled particles in airway. A training dataset is obtained from well-validated CFD simulations with realistic human airway model using Euler-Lagrangian method. The results indicate that LR, BR, and SVM yield unsatisfactory predictive accuracy, with average R² values in range of 0.21 to 0.73. Comparatively, BPNN and decisiontree-based models show great potential in predicting total deposition rate in the upper and central airway. However, for regional deposition rate prediction, BPNN did not consistently yield high accuracy, particularly for oral deposition (R² = 0.538). Comparatively, XGBoost emerges as optimal model, achieving an R² approximately close to 1 on both the training and testing datasets, with predictive errors within the range of ±0.5. The overall results demonstrate that decision-tree-based models, particularly XGBoost, have superior performance in accurately predicting both total and regional deposition rates of inhaled particles within airway. Despite limitations like geometry complexity and data quantity, the workflow developed in this study is expected to pave the way for future research integrating ML models into drug delivery device design and evaluation.

Several studies have recommended the use of hydrogels for localized targeted delivery of chemotherapeutic drugs following tumor removal surgery. This approach aims to both fill the cavity and prevent cancer recurrence. The use of Multiphysics-based simulation emerges as a valuable strategy for minimizing experimental work, providing detailed insights into how drug release occurs in the tissue, and enabling the optimization of the design. In this study, we introduced a mathematical model, utilizing experimental data, to investigate the transport of liposomes carrying MZ1 from a thermosensitive hydrogel and their impact on the viability of breast cancer cells. The proposed comprehensive model considers not just the transport within the interstitial tissue, represented as a porous medium, but also the uptake by cells and its influence on cell viability, along with the potential lymphatic drainage. The six real patient-specific tumor shapes extracted from MRI scans were used to investigate how the size and form of the tumor can modify the transport pattern. The computational results revealed that the concentration of liposomes in the tissue is significantly influenced by their release from the hydrogel, which proved to be the limiting step. Liposome concentrations of approximately 0.1 % weight were found to be sufficient in ensuring minimal cell survival in the vicinity of the tumor.

Numerous studies in the literature have proposed the use of thermo-responsive hydrogels for filling cavities after tumor resection. However, optimizing the injection process is challenging due to the complex interplay of various multi-physics phenomena, such as the coupling of flow and heat transfer, multi-phase interactions, and phase-change dynamics. Therefore, gaining a fundamental understanding of these processes is crucial. In this study, we introduce a thermo-sensitive hydrogel formulated with poloxamer 407 and Gellan gum as a promising filling agent, offering an ideal phase-transition temperature along with suitable elastic and viscous modulus properties.

We performed multi-physics simulations to predict the flow and temperature distributions during hydrogel injection. The results suggested that the hydrogel should be kept at 4 °C and injected within 90 s to avoid reaching the transition temperature. Cavity filling simulations indicated a symmetric distribution of the hydrogel, with minimal influence from the syringe's position.

The temperature gradient at the cavity edge delays gelation during injection, which is essential to guarantee its administration as a liquid. The hydrogel's viscosity follows a sigmoidal function relative to temperature, taking five minutes to reach its maximum value. In summary, the multi-physics simulations carried out in this study confirm the potential of thermo-responsive hydrogels for use in post-tumor surgery treatment and define the conditions for a proper administration. Furthermore, the proposed model can be widely applied to other thermo-responsive hydrogels or under different conditions.

BACKGROUND: Kidney disease is the most important predictor of death in patients with a Fontan circulation, yet its clinical and hemodynamic correlates have not been well established. METHODS AND RESULTS: A total of 53 ambulatory patients with a Fontan circulation (median age, 16.2 years, 52.8% male patients) underwent advanced cardiovascular magnetic resonance assessment, including 4-dimensional flow imaging and computational fluid dynamics. Estimated glomerular filtration rate (eGFR) <90 mL/min per 1.73 m2 was observed in 20.8% and albumin-to-creatinine ratio >3 mg/mmol in 39.6%. The average eGFR decline rate was -1.83 mL/min per 1.73 m2 per year (95% CI, -2.67 to -0.99; P<0.001). Lower eGFR was associated with older age, larger body surface area at examination, longer time since Fontan procedure, and lower systemic ventricular ejection fraction. Higher albumin-to-creatinine ratio was associated with absence of fenestration at the Fontan operation, and older age and lower systemic ventricular ejection fraction at the assessment. Lower cross-sectional area of the Fontan conduit indexed to flow (r=0.32, P=0.038), higher inferior vena cava-conduit velocity mismatch factor (r=-0.35, P=0.022), higher kinetic energy indexed to flow in the total cavopulmonary connection (r=-0.59, P=0.005), and higher total cavopulmonary connection resistance (r=-0.42, P=0.005 at rest; r=-0.43, P=0.004 during exercise) were all associated with lower eGFR but not with albuminuria. CONCLUSIONS: Kidney dysfunction and albuminuria are common among clinically well adolescents and young adults with a Fontan circulation. Advanced cardiovascular magnetic resonance-derived metrics indicative of declining Fontan hemodynamics are associated with eGFR and might serve as targets to improve kidney health. Albuminuria might be driven by other factors that need further investigation.

Background
Properly understanding the origin and progression of the thoracic aortic aneurysm (TAA) can help prevent its growth and rupture. For a better understanding of this pathogenesis, the aortic blood flow has to be studied and interpreted in great detail. We can obtain detailed aortic blood flow information using magnetic resonance imaging (MRI) based computational fluid dynamics (CFD) with a prescribed motion of the aortic wall.

Methods
We performed two different types of simulations—static (rigid wall) and dynamic (moving wall) for healthy control and a patient with a TAA. For the latter, we have developed a novel morphing approach based on the radial basis function (RBF) interpolation of the segmented 4D-flow MRI geometries at different time instants. Additionally, we have applied reconstructed 4D-flow MRI velocity profiles at the inlet with an automatic registration protocol.

Results
The simulated RBF-based movement of the aorta matched well with the original 4D-flow MRI geometries. The wall movement was most dominant in the ascending aorta, accompanied by the highest variation of the blood flow patterns. The resulting data indicated significant differences between the dynamic and static simulations, with a relative difference for the patient of 7.47±14.18% in time-averaged wall shear stress and 15.97±43.32% in the oscillatory shear index (for the whole domain).

Conclusions
In conclusion, the RBF-based morphing approach proved to be numerically accurate and computationally efficient in capturing complex kinematics of the aorta, as validated by 4D-flow MRI. We recommend this approach for future use in MRI-based CFD simulations in broad population studies. Performing these would bring a better understanding of the onset and growth of TAA.

The treatment for asthma and chronic obstructive pulmonary disease relies on forced inhalation of drug particles. Their distribution is essential for maximizing the outcomes. Patient-specific computational fluid dynamics (CFD) simulations can be used to optimize these therapies. In this regard, this study focuses on creating a parametric model of the human respiratory tract from which synthetic anatomies for particle deposition analysis through CFD simulation could be derived. A baseline geometry up to the fourth generation of bronchioles was extracted from a CT dataset. Radial basis function (RBF) mesh morphing acting on a dedicated tree structure was used to modify this baseline mesh, extracting 1000 synthetic anatomies. A total of 26 geometrical parameters affecting branch lengths, angles, and diameters were controlled. Morphed models underwent CFD simulations to analyze airflow and particle dynamics. Mesh morphing was crucial in generating high-quality computational grids, with 96% of the synthetic database being immediately suitable for accurate CFD simulations. Variations in wall shear stress, particle accretion rate, and turbulent kinetic energy across different anatomies highlighted the impact of the anatomical shape on drug delivery and deposition. The study successfully demonstrates the potential of tree-structure-based RBF mesh morphing in generating parametric airways for drug delivery studies.

In the present study, we have applied a combined wall-resolving dynamic Large-Eddy Simulation (LES) (for the velocity field) and Direct Numerical Simulation (DNS) (for the temperature field) approach for mixing of parallel triple-jets with different temperatures of liquid sodium in a turbulent forced convection regime. Because of the high thermal conductivity of sodium (a low-Prandtl fluid), we adopted the dynamic Smagorinsky subgrid closure for the unresolved velocity scales, while the thermal scales are fully resolved. Furthermore, the Time-dependent Reynolds-Averaged Navier-Stokes (T-RANS) approach with the high-Reynolds number variant (i.e. with the wall functions as boundary conditions along solid boundaries) of the four-equation eddy viscosity model (k−ε−k_θ−ε_θ) was applied. The fine-mesh LES/DNS provided a close agreement with the experimental data for both velocity and temperature fields (for both first- and second-moments). In contrast, the coarse-mesh LES/DNS overestimated the turbulent kinetic energy profiles at different distances from the inlet plane. The T-RANS results confirmed a good agreement with the mean streamwise velocity and turbulent kinetic energy, as well as the mean temperature profiles. Finally, the analysis of power spectral density distributions of the temperature signal revealed that all simulation techniques captured a dominant flow frequency originating from the induced Kelvin-Helmholtz instabilities between the side and central jets. The presented combined dynamic LES/DNS approach is recommended for future simulations of the turbulent forced convection flows of low Prandtl fluids, especially if thermal fatigue effects need to be predicted correctly.

The present study addresses numerical simulations of the double-diffusive convection in a turbulent regime. The characteristic concentration and temperature Prandtl numbers (PrC = 700 and PrT = 7) correspond to the typical seawater properties. Here, we applied the Large-Eddy Simulation (LES) approach, with relatively simple subgrid turbulence closures of momentum, concentration, and temperature for the unresolved scales. The wall-resolved LES approach proved to work well on significantly coarser numerical mesh than used in recent Direct Numerical Simulation (DNS) studies of Yang et al. (2016) in the intermediate range of working parameters, 107 ≤ RaC ≤ 109, 0 ≤ RaT ≤ 106, which covers the quasi-Rayleigh-Bénard, fingering, and damping flow regimes. A good agreement between concentration and temperature Nusselt numbers is obtained. The instantaneous Nusselt numbers distribution revealed a significant impact of the imposed strong thermal stratification (RaT = 106) in comparison to the neutral case (RaT = 0).

The paper reports on the numerical turbulence model in predicting mass, momentum and heat transfer in a partially divided cavity heated from the side using buoyancy-extended eddy-viscosity and elliptic relaxation approach with the algebraic expressions for the Reynold stress tensor and turbulent heat flux vector. The CDS (central differencing scheme) and LUDS (linear upwind differencing scheme) were used as the discretization method and the governing equations were solved using the finite volume method and Navier-Stokes solver. Validation of the model has been carried out by experimental data of convective flow in the cavity as well as by numerical data DNS (direct numerical simulation). The model agrees very well with the experiment and DNS and it is also able to demonstrate the performance which is comparable to that of the previous advanced second-moment closure model (SMC) in the literature. The results show that the model is suitable for use in simulations of the turbulent convective flow in a cavity with partition and it has the potential to be applied to more complex cavities and a wide range of turbulence levels.

The turbulent penetrative convection into a stable convective boundary layer represents an important phenomenon in environmental engineering and atmospheric science. In the present study, we present a series of numerical simulations performed by two modeling approaches: the high-fidelity Large-Eddy Simulations (LES), and the less computationally demanding transient Reynolds-Averaged Approach (TRANS), but with an advanced sub-scale turbulent heat flux model. By simulating different localized heat sources over the ground, and by performing a direct comparative assessment of results obtained by LES and TRANS, we confirmed an overall good agreement in predicting the time evolution of the horizontally averaged temperature profiles. Similarly, the morphology of instantaneous thermal plumes and large convective structures predicted by TRANS were in reasonable agreement with the referent LES predictions.

Thoracic aortic aneurysm (TAA) is a local dilation of the thoracic aorta. Although universally used, aneurysm diameter alone is a poor predictor of major complications such as rupture. There is a need for better biomarkers for risk assessment that also reflect the aberrant flow patterns found in TAAs. Furthermore, hypertension is often present in TAA patients and may play a role in progression of aneurysm. The exact relation between TAAs and hypertension is poorly understood. This study aims to create a numerical model of hypertension in the aorta by using computational fluid dynamics. First, a normotensive state was simulated in which flow and resistance were kept unaltered. Second, a hypertensive state was modeled in which blood inflow was increased by 30%. Third, a hypertensive state was modeled in which the proximal and peripheral resistances and capacitance parameters from the three-element Windkessel boundary condition were adjusted to mimic an increase in resistance of the rest of the cardiovascular system. One patient with degenerative TAA and one healthy control were successfully simulated at hypertensive states and were extensively analyzed. Furthermore, three additional TAA patients and controls were simulated to validate our method. Hemodynamic variables such as wall shear stress, oscillatory shear index, endothelial cell activation potential (ECAP), vorticity and helicity were studied to gain more insight on the effects of hypertension on flow patterns in TAAs. By comparing a TAA patient and a control at normotensive state at peak-systole, helicity and vorticity were found to be lower in the TAA patient throughout the entire domain. No major changes in flow and flow derived quantities were observed for the TAA patient and control when resistance was increased. When flow rate was increased, regions with high ECAP values were found to reduce in TAA patients in the aneurysm region which could reduce the risk of thrombogenesis. Thus, it may be important to assess cardiac output in patients with TAA.

This paper presents a study of developing a numerical turbulent model in a cavity heated from below using eddy viscosity combined with elliptic relaxation approach. The model uses a set of differential equations that consist of kinetic energy, its dissipation, variance of temperature, velocity scale and elliptic relaxation parameter, which are solved using a finite-volume and Navier-Stokes solver. The unresolved stress tensors and heat flux vectors are modelled with an algebraic formula. The discretization method is carried out by CDS, or second-order differencing scheme, and LUDS, or second-order linear upwind scheme. The model is applied to the natural convection heated from below, known as Rayleigh-Benard convection, in a two-dimensional cavity with a height-to-length aspect ratio of 1:1.5, 1:4, and 1:8. The model has been validated using numerical data from DNS (direct numerical simulation) and experiments. The model produced similar results with both DNS and experiments. It was also shown that the model can visualize the main feature of turbulent convective flow in the enclosure for various aspect ratio.