Fluids, Volume 6, Issue 8 (August 2021) – 38 articles

Wavelet transform has become a common tool for processing non-stationary signals in many different fields. The present paper reports a review of some applications of wavelet in aeroacoustics with a special emphasis on the analysis of experimental data taken in compressible jets. The focus is on three classes of wavelet-based signal processing procedures: (i) conditional statistics; (ii) acoustic and hydrodynamic pressure separation; (iii) stochastic modeling. The three approaches are applied to an experimental database consisting of pressure time series measured in the near field of a turbulent jet. Future developments and possible generalization to other applications, e.g., airframe or propeller noise, are also discussed. Full article

In this paper, we consider the poroelasticity problem in fractured and heterogeneous media. The mathematical model contains a coupled system of equations for fluid pressures and displacements in heterogeneous media. Due to scale disparity, many approaches have been developed for solving detailed fine-grid problems on a coarse grid. However, some approaches can lack good accuracy on a coarse grid and some corrections for coarse-grid solutions are needed. In this paper, we present a coarse-grid approximation based on the generalized multiscale finite element method (GMsFEM). We present the construction of the offline and online multiscale basis functions. The offline multiscale basis functions are precomputed for the given heterogeneity and fracture network geometry, where for the construction, we solve a local spectral problem and use the dominant eigenvectors (appropriately defined) to construct multiscale basis functions. To construct the online basis functions, we use current information about the local residual and solve coupled poroelasticity problems in local domains. The online basis functions are used to enrich the offline multiscale space and rapidly reduce the error using residual information. Only with appropriate offline coarse-grid spaces can one guarantee a fast convergence of online methods. We present numerical results for poroelasticity problems in fractured and heterogeneous media. We investigate the influence of the number of offline and online basis functions on the relative errors between the multiscale solution and the reference (fine-scale) solution. Full article

The paper describes a distributed temperature sensing system that was used to monitor the artificial freezing of soils during the construction of a potash mine shaft. The technique of reconstructing the temperature field by solving the inverse problem in the entire volume of frozen soils using the measured temperatures in four thermal monitoring (TM) wells is described. Two local anomalies in temperature distributions in TM wells are described and analyzed theoretically using thermo-hydraulic modeling. The first anomaly concerns the asymmetric temperature distribution in one of the soil layers and is associated with the influence of natural groundwater flow in the horizontal direction. The second anomaly consists of a sharp decrease in water temperature in the section of the TM well located inside the freezing contour. Calculations showed that it is most likely associated with the entry of cold groundwater from the overlying layers of soils through a well filter at a depth of 160 m and the subsequent movement of the water up the well. Full article

Geometrically parametrized partial differential equations are currently widely used in many different fields, such as shape optimization processes or patient-specific surgery studies. The focus of this work is some advances on this topic, capable of increasing the accuracy with respect to previous approaches while relying on a high cost–benefit ratio performance. The main scope of this paper is the introduction of a new technique combining a classical Galerkin-projection approach together with a data-driven method to obtain a versatile and accurate algorithm for the resolution of geometrically parametrized incompressible turbulent Navier–Stokes problems. The effectiveness of this procedure is demonstrated on two different test cases: a classical academic back step problem and a shape deformation Ahmed body application. The results provide insight into details about the properties of the architecture we developed while exposing possible future perspectives for this work. Full article

Ferrofluid sloshing vibration energy harvesters use ferrofluid sloshing movement as a moving magnet between a fixed coil to induce current and, in turn, harvest energy from external excitations. A symmetric ferrofluid sloshing vibration energy harvester configuration is introduced in this study which utilizes four external, symmetrically placed, permanent magnets to magnetize a ferrofluid inside a tank. An external sinusoidal excitation of amplitude 1 m/s² is imparted, and the whole system is studied numerically using a level-set method to track the sharp interface between ferrofluid and air. The system is studied for two significant length scales of 0.1 m and 0.05 m while varying the four external magnets’ polarity arrangements. All of the system configuration dimensions are parametrized with the length scale to keep the system configuration invariant with the length scale. Finally, a frequency sweep is performed, encompassing the structure’s first modal frequency and impedance matching to obtain the system’s energy harvesting characteristics. Full article

In this study, a fully coupled electro-thermo-mechanical model of radiofrequency (RF)-assisted cardiac ablation has been developed, incorporating fluid–structure interaction, thermal relaxation time effects and porous media approach. A non-Fourier based bio-heat transfer model has been used for predicting the temperature distribution and ablation zone during the cardiac ablation. The blood has been modeled as a Newtonian fluid and the velocity fields are obtained utilizing the Navier–Stokes equations. The thermal stresses induced due to the heating of the cardiac tissue have also been accounted. Parametric studies have been conducted to investigate the effect of cardiac tissue porosity, thermal relaxation time effects, electrode insertion depths and orientations on the treatment outcomes of the cardiac ablation. The results are presented in terms of predicted temperature distributions and ablation volumes for different cases of interest utilizing a finite element based COMSOL Multiphysics software. It has been found that electrode insertion depth and orientation has a significant effect on the treatment outcomes of cardiac ablation. Further, porosity of cardiac tissue also plays an important role in the prediction of temperature distribution and ablation volume during RF-assisted cardiac ablation. Moreover, thermal relaxation times only affect the treatment outcomes for shorter treatment times of less than 30 s. Full article

This work extends the characteristic-based volume penalization method, originally developed and demonstrated for compressible subsonic viscous flows in (J. Comput. Phys. 262, 2014), to a hyperbolic system of partial differential equations involving complex domains with moving boundaries. The proposed methodology is shown to be Galilean-invariant and can be used to impose either homogeneous or inhomogeneous Dirichlet, Neumann, and Robin type boundary conditions on immersed boundaries. Both integrated and non-integrated variables can be treated in a systematic manner that parallels the prescription of exact boundary conditions with the approximation error rigorously controlled through an a priori penalization parameter. The proposed approach is well suited for use with adaptive mesh refinement, which allows adequate resolution of the geometry without over-resolving flow structures and minimizing the number of grid points inside the solid obstacle. The extended Galilean-invariant characteristic-based volume penalization method, while being generally applicable to both compressible Navier–Stokes and Euler equations across all speed regimes, is demonstrated for a number of supersonic benchmark flows around both stationary and moving obstacles of arbitrary shape. Full article

Analytical and numerical investigations were performed to study the influence of the Soret and Dufour effects on double-diffusive convection in a vertical porous layer filled with a binary mixture and subject to horizontal thermal and solute gradients. In particular, the study was focused on the effect of Soret and Dufour diffusion on bifurcation types from the rest state toward steady convective state, and then toward oscillatory convective state. The Brinkman-extended Darcy model and the Boussinesq approximation were employed to model the convective flow within the porous layer. Following past laboratory experiments, the investigations dealt with the particular situation where the solutal and thermal buoyancy forces were equal but acting in opposite direction to favor the possible occurrence of the rest state condition. For this situation, the onset of convection could be either supercritical or subcritical and occurred at given thresholds and following various bifurcation routes. The analytical investigation was based on the parallel flow approximation, which was valid only for a tall porous layer. A numerical linear stability analysis of the diffusive and convective states was performed on the basis of the finite element method. The thresholds of supercritical, $R_{TC}^{sup}$, and overstable, $R_{TC}^{over}$, convection were computed. In addition, the stability of the established convective flow, predicted by the parallel flow approximation, was studied numerically to predict the onset of Hopf’s bifurcation, $R_{TC}^{Hopf}$, which marked the transition point from steady toward unsteady convective flows; a route towards the chaos. To support the analytical analyses of the convective flows and the numerical stability methodology and results, nonlinear numerical solutions of the full governing equations were obtained using a second-order finite difference method. Overall, the Soret and Dufour effects were seen to affect significantly the thresholds of stationary, overstable and oscillatory convection. The Hopf bifurcation was marked by secondary convective flows consisting of superposed vertical layers of opposite traveling waves. A good agreement was found between the predictions of the parallel flow approximation, the numerical solution and the linear stability results. Full article

The Alber equation is a phase-averaged second-moment model used to study the statistics of a sea state, which has recently been attracting renewed attention. We extend it in two ways: firstly, we derive a generalized Alber system starting from a system of nonlinear Schrödinger equations, which contains the classical Alber equation as a special case but can also describe crossing seas, i.e., two wavesystems with different wavenumbers crossing. (These can be two completely independent wavenumbers, i.e., in general different directions and different moduli.) We also derive the associated two-dimensional scalar instability condition. This is the first time that a modulation instability condition applicable to crossing seas has been systematically derived for general spectra. Secondly, we use the classical Alber equation and its associated instability condition to quantify how close a given nonparametric spectrum is to being modulationally unstable. We apply this to a dataset of 100 nonparametric spectra provided by the Norwegian Meteorological Institute and find that the vast majority of realistic spectra turn out to be stable, but three extreme sea states are found to be unstable (out of 20 sea states chosen for their severity). Moreover, we introduce a novel “proximity to instability” (PTI) metric, inspired by the stability analysis. This is seen to correlate strongly with the steepness and Benjamin–Feir Index (BFI) for the sea states in our dataset (>85% Spearman rank correlation). Furthermore, upon comparing with phase-resolved broadband Monte Carlo simulations, the kurtosis and probability of rogue waves for each sea state are also seen to correlate well with the PTI (>85% Spearman rank correlation). Full article

A new visualization technology is presented, which was used in applied research when observing and modeling the dynamics of the flow of gaseous environments. In the process of developing and improving the technology, a set of experimental results was compiled to study the phenomenon of combustion and detonation of a hydrogen-oxygen mixture, as well as the phenomena of propagation, action, and interaction of shock waves and gas-dynamic structures. On the example of analyzing data on the dynamics of the formation of a vortex ring, the possibilities of verifying the computational model of the implemented physical process are shown. The presented results reflect the level of information content when using technology to carry out tests. Full article

This paper presents a numerical study on the impact of partial leeward porous coatings on the drag of circular cylinders in cross-flow. Porous coatings are receiving increasing attention for their potential in passive flow control. An unsteady Reynolds-averaged Navier–Stokes model was developed that agreed well with the numerical and experimental literature. Using the two-equation shear stress transport $k-\omega$ turbulence model, 2D flow around a circular cylinder was simulated at $Re$ = $4.2\times {10}^4$ with five different angles of partial leeward porous coatings and a full porous coating. For coating angles below ${130}^{\circ }$, the coating resulted in an increase in pressure on the leeward side of the cylinder. There was a significant reduction in the fluctuation of the pressure and aerodynamic forces and a damping effect on vortex shedding. Flow separation occurred earlier; the wake was widened; and there was a decrease in turbulence intensity at the outlet. A reduction of drag between 5 and 16% was measured, with the maximum at a 70${}^{\circ }$ coating angle. The results differed greatly for a full porous coating and a 160${}^{\circ }$ coating, which were found to cause an increase in drag of 42% and 43%, respectively. The results showed that leeward porous coatings have a clear drag-reducing potential, with possibilities for further research into the optimum configuration. Full article

The availability of reliable and efficient turbulent flow simulation methods is highly beneficial for wind energy and aerospace developments. However, existing simulation methods suffer from significant shortcomings. In particular, the most promising methods (hybrid RANS-LES methods) face divergent developments over decades, there is a significant waste of resources and opportunities. It is very likely that this development will continue as long as there is little awareness of conceptional differences of hybrid methods and their implications. The main purpose of this paper is to contribute to such clarification by identifying a basic requirement for the proper functioning of hybrid RANS-LES methods: a physically correct communication of RANS and LES modes. The state of the art of continuous eddy simulations (CES) methods (which include the required mode communication) is described and requirements for further developments are presented. Full article

Calcification and bicuspid valve formation are important aortic valve disorders that disturb the hemodynamics and the valve function. The detailed analysis of aortic valve hemodynamics would lead to a better understanding of the disease’s etiology. We computationally modeled the aortic valve using simplified three-dimensional geometry and inlet velocity conditions obtained via echocardiography. We examined various calcification severities and bicuspid valve formation. Fluid-structure interaction (FSI) analyses were adapted using ANSYS Workbench to incorporate both flow dynamics and leaflet deformation accurately. Simulation results were validated by comparing leaflet movements in B-mode echo recordings. Results indicate that the biomechanical environment is significantly changed for calcified and bicuspid valves. High flow jet velocities are observed in the calcified valves which results in high transvalvular pressure difference (TPG). Wall shear stresses (WSS) increased with the calcification on both fibrosa (aorta side) and ventricularis (left ventricle side) surfaces of the leaflet. The WSS distribution is regular on the ventricularis, as the WSS values proportionally increase from the base to the tip of the leaflet. However, WSS patterns are spatially complex on the fibrosa side. Low WSS levels and spatially complex WSS patterns on the fibrosa side are considered as promoting factors for further calcification and valvular diseases. Full article

The preferential organisation of coherent vortices in a turbulent boundary layer in relation to local large-scale streamwise velocity features was investigated. Coherent vortices were identified in the wake region using the Triple Decomposition Method (originally proposed by Kolář) from 2D particle image velocimetry (PIV) data of a canonical turbulent boundary layer. Two different approaches, based on conditional averaging and quantitative statistical analysis, were used to analyze the data. The large-scale streamwise velocity field was first conditionally averaged on the height of the detected coherent vortices and a change in the sign of the average large scale streamwise fluctuating velocity was seen depending on the height of the vortex core. A correlation coefficient was then defined to quantify this relationship between the height of coherent vortices and local large-scale streamwise fluctuating velocity. Both of these results indicated a strong negative correlation in the wake region of the boundary layer between vortex height and large-scale velocity. The relationship between vortex height and full large-scale velocity isocontours was also studied and a conceptual model based on the findings of the study was proposed. The results served to relate the hairpin vortex model of Adrian et al. to the scale interaction results reported by Mathis et al., and Chung and McKeon. Full article

Two-dimensional computational fluid dynamics (CFD) simulations are employed to investigate the effects of adding a circular-arc-shaped windshield on the performance of a Banki type vertical axis wind turbine (VAWT), particularly to the application where the VAWT is harnessing wind energy in highway caused by the passing vehicles. Unsteady Reynolds-Averaged Navier-Stokes (URANS) is the computational approach used to calculate the turbulent flow within the domain. Two sets of simulation cases based on two different vehicles (i.e., car and a bus) are performed with varying arc-length of the windshield. The results show that the windshield provides an increase in the energy captured by the VAWT by up to 16.14% compared to no windshield case when the car model is used. In contrast, windshield in all the simulation cases using a bus model gives a negative effect to VAWT performance where the worst case yields −64.77%. Full article

Cardiovascular diseases are a major cause of death and disability worldwide and they are commonly associated with the occurrence of atherosclerotic plaque deposition in the vessel walls, a process denoted as atherosclerosis. This is a chronic and progressive inflammatory disease of large-/medium-sized blood vessels that affects blood flow profiles, with the abdominal aorta and its branches being one of the locations prone to the development of this pathology, due to their curvatures and bifurcations. In this regard, the effect of flow patterns was studied and compared for both a simplified three-dimensional model of aorta bifurcation on the iliac arteries and a realistic model of iliac bifurcation, which was constructed from a computational tomography medical image. The flow patterns were analyzed in terms of velocity and wall shear stress distribution, but a special focus was given to the size and location of the recirculation zone. The simulations were performed using the Computational Fluid Dynamics software, FLUENT, taking into account the cardiac cycle profile at the infrarenal aorta. The shear stress and the velocity distribution observed for both models indicated that higher shear stress occurred along the flow divider wall (inner wall) and low shear stress occurred along the outer walls. In addition, the results demonstrated that the wall shear stress profiles were deeply affected by the transient profile of the cardiac cycle, with the deceleration phase being the most critical phase to the occurrence of backflow. Full article

We remark on the use of regularized Stokeslets in the slender body theory (SBT) approximation of Stokes flow about a thin fiber of radius $ϵ>0$. Denoting the regularization parameter by $\delta$, we consider regularized SBT based on the most common regularized Stokeslet plus a regularized doublet correction. Given sufficiently smooth force data along the filament, we derive $L^{\infty }$ bounds for the difference between regularized SBT and its classical counterpart in terms of $\delta$, $ϵ$, and the force data. We show that the regularized and classical expressions for the velocity of the filament itself differ by a term proportional to $log(\delta /ϵ)$; in particular, $\delta =ϵ$ is necessary to avoid an $O\left(1\right)$ discrepancy between the theories. However, the flow at the surface of the fiber differs by an expression proportional to $log(1+{\delta }^2/{ϵ}^2)$, and any choice of $\delta \propto ϵ$ will result in an $O\left(1\right)$ discrepancy as $ϵ\to 0$. Consequently, the flow around a slender fiber due to regularized SBT does not converge to the solution of the well-posed slender body PDE which classical SBT is known to approximate. Numerics verify this $O\left(1\right)$ discrepancy but also indicate that the difference may have little impact in practice. Full article

Nonlinear dynamics of patterns near the threshold of long-wave monotonic Marangoni instability of conductive state in a heated thin layer of liquid covered by insoluble surfactant is considered. Pattern selection between roll and square planforms is analyzed. The dependence of pattern stability on the heat transfer from the free surface of the liquid characterized by Biot number and the gravity described by Galileo number at different surfactant concentrations is studied. Using weakly nonlinear analysis, we derive a set of amplitude equations governing the large-scale roll distortions in the presence of the surface deformation and the surfactant redistribution. These equations are used for the linear analysis of modulational instability of stationary rolls. Full article

We provide a new framework for analyzing the flow of an axisymmetric liquid film flowing down a vertical fiber, applicable to fiber coating flows and those in similar geometries in heat exchangers, water treatment, and desalination processes. The problem considered is that of a viscous liquid film falling under the influence of gravity and surface tension on a solid cylindrical fiber. Our approach is different from existing ones in that we derive our mathematical model by using a control-volume approach to express the conservation of mass and axial momentum in simple and intuitively appealing forms, resulting in a pair of equations that are reminiscent of the Saint-Venant shallow-water equations. Two versions of the model are obtained, one assuming a plug-flow velocity profile with a linear drag force expression, and the other using the fully-developed laminar velocity profile for a locally uniform film to approximate the drag. These can, respectively, model high- and low-Reynolds number regimes of flow. Linear stability analyses and fully nonlinear numerical simulations are presented that show the emergence of traveling wave solutions representing chains of identical droplets falling down the fiber. Physical experiments with safflower oil on a fishing line are also undertaken and match the theoretical predictions from the laminar flow model well when machine learning methods are used to estimate the parameters. Full article

We propose a new hyper-reduction method for a recently introduced nonlinear model reduction framework based on dynamically transformed basis functions and especially well-suited for transport-dominated systems. Furthermore, we discuss applying this new method to a wildland fire model whose dynamics feature traveling combustion waves and local ignition and is thus challenging for classical model reduction schemes based on linear subspaces. The new hyper-reduction framework allows us to construct parameter-dependent reduced-order models (ROMs) with efficient offline/online decomposition. The numerical experiments demonstrate that the ROMs obtained by the novel method outperform those obtained by a classical approach using the proper orthogonal decomposition and the discrete empirical interpolation method in terms of run time and accuracy. Full article

Wall jets are applied to reduce flow separation and recirculation of the airflow entering the inner space of a laboratory fume hood through its front opening. The flow separation in the hood was further reduced by introducing a self-induced oscillatory motion using fluidic oscillators. The design of the oscillators integrated in the predetermined contour are based on numerical simulations. The effect of the steady and unsteady wall jet was investigated experimentally using flow visualization, particle image velocimetry (PIV), and containment measurements. The oscillatory wall-jet led to reduction of flow separation and recirculation even at lower injection volume flows. In consequence, the usage of fluidic oscillators for a laboratory fume hood increases the energy efficiency of the system without reducing the safety of the laboratory fume hood. Full article

A numerical simulation of the gas-dynamic processes in the thrust vectorable nozzle of the solid rocket motor is considered. Construction of a geometric model and a generation of computational mesh, and reconstruction of model and mesh at each time step are discussed. Calculations of the flowfield of combustion products in the pre-nozzle chamber and nozzle block are carried out for various angles of nozzle rotation. The distributions of the gas dynamic quantities in the pre-nozzle volume corresponding to the outflow of the combustion products from the cylindrical channel and star-shaped channel are compared, as well as the solutions of the problem obtained with quasi-stationary and unsteady formulations. The effects of the channel shape on the distribution of flow quantities and formation of a vortical flow structure in the nozzle block are discussed. Full article

Although recent advances of four-dimensional (4D) flow magnetic resonance imaging (MRI) has introduced a new way to measure Reynolds stress tensor (RST) in turbulent flows, its measurement accuracy and possible bias have remained to be revealed. The purpose of this study was to compare the turbulent flow measurement of 4D flow MRI and particle image velocimetry (PIV) in terms of velocity and turbulence quantification. Two difference flow rates of 10 and 20 L/min through a 50% stenosis were measured with both PIV and 4D flow MRI. Not only velocity through the stenosis but also the turbulence parameters such as turbulence kinetic energy and turbulence production were quantitatively compared. Results shows that 4D flow MRI velocity measurement well agreed with the that of PIV, showing the linear regression slopes of two methods are 0.94 and 0.89, respectively. Although turbulence mapping of 4D flow MRI was qualitatively agreed with that of PIV, the quantitative comparison shows that the 4D flow MRI overestimates RST showing the linear regression slopes of 1.44 and 1.66, respectively. In this study, we demonstrate that the 4D flow MRI visualize and quantify not only flow velocity and also turbulence tensor. However, further optimization of 4D flow MRI for better accuracy might be remained. Full article

An unstable discharge rate occurs during dry fine powder discharge from a hopper because of the significant two-phase solid/gas interactions that occur in powder flows. In addition, the air bubble phenomenon may occur in a silo during fine powder discharge. In this study, we conducted experiments using a semi-conical dual-structure hopper, and examined the effects on the hopper internal flow structure, cavity fluid pressure, pressure inside the airtight cavity section, and the powder discharge rate when changes are made in the position of the supplied air injection port and the solenoid valve open/close timing. From the experimental results, it was confirmed that an appropriate pressure supply port position exists, and the change in expansion/contraction of the flexible container due to air vibration is determined by the balance between the amount of air inserted and the amount of air discharged, and does not affect the presence or absence of powder so much. Furthermore, as the pressure value in the airtight void is directly related to the change in the expansion and contraction of the flexible container, the maximum amplitude value of the pressure in the airtight void can be kept high and constant at the time of opening and closing the solenoid valve. Full article

One essential utilization of phase change materials as energy storage materials is energy saving and temperature control in air conditioning and indirect solar air drying systems. This study presents an experimental investigation evaluating the characteristics and energy savings of multiple phase change materials subjected to internal flow in an air heating system during charging and discharging cycles. The experimental tests were conducted using a test rig consisting of two main parts, an air supply duct and a room model equipped with phase change materials (PCMs) placed in rectangular aluminum panels. Analysis of the results was based on three test cases: PCM1 (Paraffin wax) placed in the air duct was used alone in the first case; PCM2 (RT–42) placed in the room model was used alone in the second case; and in the third case, the two PCMs (PCM1 and PCM2) were used at the same time. The results revealed a significant improvement in the energy savings and room model temperature control for the air heating system incorporated with multiple PCMs compared with that of a single PCM. Complete melting during the charging cycle occurred at temperatures in the range of 57–60 °C for PCM1 and 38–43 °C for PCM2, respectively, thereby validating the reported PCMs’ melting–solidification results. Multiple PCMs maintained the room air temperature at the desired range of 35–45.2 °C in the air heating applications by minimizing the air temperature fluctuations. The augmentation in discharging time and improvement in the room model temperature using multiple PCMs were about 28.4% higher than those without the use of PCMs. The total energy saving using two PCMs was higher by about 29.5% and 46.7% compared with the use of PCM1 and PCM2, respectively. It can be concluded that multiple PCMs have revealed higher energy savings and thermal stability for the air heating system considered in the current study. Full article

This paper presents the Eulerian–Lagrangian approach for numerical modeling of high-speed gas-droplet flows and aeroacoustics. The proposed hybrid approach is implemented using the OpenFOAM library and two different methods. The first method is based on a hybrid convective terms approximation method employing a Kurganov–Tadmor and PIMPLE scheme. The second method employs the regularized or quasi-gas dynamic equations. The Lagrangian part of the flow description uses the OpenFOAM cloud model. Within this model, the injected droplets are simulated as packages (parcels) of particles with constant mass and diameter within each parcel. According to this model, parcels moving in the gas flow could undergo deceleration, heating, evaporation, and breakup due to hydrodynamic instabilities. The far-field acoustic noise is predicted using Ffowcs Williams and Hawking’s analogy. The Lagrangian model is verified using the cases with droplet evaporation and motion. Numerical investigation of water microjet injection into the hot ideally expanded jet allowed studying acoustic properties and flow structures, which emerged due to the interaction of gas and liquid. Simulation results showed that water injection with a mass flow rate equal to 13% of the gas jet mass flow rate reduced the noise by approximately 2 dB. This result was in good coincidence with the experimental observations, where maximum noise reduction was about 1.6 dB. Full article

With the development of mature Computational Fluid Dynamics (CFD) tools for fluids (air and liquid) and Finite Element Methods (FEM) for solids and structures, many approaches have been proposed to tackle the so-called Fluid–Structure Interaction or Fluid–Solid Interaction (FSI) problems. Traditional partitioned iterations are often used to link available FEM codes with CFD codes in the study of FSI systems. Although these procedures are convenient, fluid mesh adjustments according to the motion and finite deformation of immersed solids or structures can be challenging or even prohibitive. Moreover, complex dynamic behaviors of coupled FSI systems are often lost in these iterative processes. In this paper, the author would like to review the so-called monolithic approaches for the solution of coupled FSI systems as a whole in the context of the immersed boundary method. In particular, the focus is on the implicit monolithic algorithm for compressible solids immersed inside a compressible liquid. Notice here the main focus of this paper is on liquid or more precisely liquid phase of water as working fluid. Using the word liquid, the author would like to emphasize the consideration of the compressibility of the fluid and the assumption of constant density and temperature. It is a common practice to assume that the pressure variations are not strong enough to alter the liquid density in any significant fashion for acoustic fluid–solid interactions problems. Although the algorithm presented in this paper is not directly applicable to aerodynamics in which the density change is significant along with its relationship with the pressure and the temperature, the author did revisit his earlier work on merging immersed boundary method concepts with a fully-fledged compressible aerodynamic code based on high-order compact scheme and energy conservative form of governing equations. In the proposed algorithm, on top of a uniform background (ghost) mesh, a fully implicit immersed method is implemented with mixed finite element methods for compressible liquid as well as immersed compressible solids with a matrix-free Newton–Krylov iterative solution scheme. In this monolithic approach, with the simple modulo function, the immersed solid or structure points can be easily located and thus the displacement projections and force distributions stipulated in the immersed boundary method can be effectively and efficiently implemented. This feature coupled with the key concept of the immersed boundary method helps to avoid topologically challenging mesh adjustments and to incorporate parallel processing commands such as Message Passing Interface (MPI) and further vectorization of the numerical operation. Once these high-performance procedures are implemented coupled with the monolithic implicit matrix-free Newton–Krylov iterative scheme with immersed methods, effective and efficient reduced order modeling techniques can then be employed to explore phase and parametric spaces. The in-house developed programs are at the moment two-dimensional. Furthermore, based on the same approach implemented in one-dimensional test example with one continuum immersed in another continuum, such monolithic implicit matrix-free Newton–Krylov iterative approach can be extended for the study of composites with deformable aggregates and matrix. Full article

This work discusses in vivo experiments that were performed to evaluate whether local or whole-body heating to 40 °C reduced interstitial fluid pressures (IFPs) and enhanced nanoparticle delivery to subcutaneous PC3 human prostate cancer xenograft tumors in mice. After heating, 0.2 mL of a previously developed nanofluid containing gold nanoparticles (10 mg Au/mL) was injected via the tail vein. The induced whole-body hyperthermia led to increases in tumor and mouse body blood perfusion rates of more than 50% and 25%, respectively, while the increases were much smaller in the local heating group. In the whole-body hyperthermia groups, the IFP reduction from the baseline at the tumor center immediately after heating was found to be statistically significant when compared to the control group. The 1 h of local heating group showed IFP reductions at the tumor center, while the IFPs increased in the periphery of the tumor. The intratumoral gold nanoparticle accumulation was quantified using inductively coupled plasma mass spectrometry (ICP-MS). Compared to the control group, 1 h or 4 h of experiencing whole-body hyperthermia resulted in an average increase of 51% or 67% in the gold deposition in tumors, respectively. In the 1 h of local heating group, the increase in the gold deposition was 34%. Our results suggest that 1 h of mild whole-body hyperthermia may be a cost-effective and readily implementable strategy for facilitating nanoparticle delivery to PC3 tumors in mice. Full article

This paper presents an analysis of linear viscous stress Favre averaged turbulence models for computational fluid dynamics (CFD) of fully turbulent round jets with a long straight tube geometry in the near field. Although similar work has been performed in the past with very relevant solutions, considerations were not given for the issues and limitations involved with coupling between an Eulerian and Lagrangian phase, such as in fully two-way coupled CFD-DEM. These issues include limitations on solution domain, mesh cell size, wall modelling, and momentum coupling between the two phases in relation to the particles size. Therefore, within these considerations, solutions are provided to the Navier–Stokes equations with various turbulence models using a three-dimensional wedge long straight tube geometry for fully developed turbulence flow. Simulations are performed with a Reynolds number of 13,000 and 51,000 using two different tube diameters. It is found that a modified k-$\epsilon$ turbulence model achieved the most agreeable results for both the velocity and turbulent flow fields between these two flow regimes, while a modified k-$\omega$ SST/BSL also provided suitable results. Full article

This study presents an evaluation of semi-empirical single-point wall pressure spectrum models by comparing model predictions with wind tunnel and flight test data. The mean squared error was used to compare the power spectral density of the wall pressure fluctuations predicted by semi-empirical models with a large amount of experimental data. Results show that the models proposed by Goody and Smol’yakov have the lowest mean squared error when predicting the power spectral density for wind tunnel experiments and the Rackl and Weston model has the lowest mean squared error when predicting the power spectral density for flight test data. In addition, although current studies of the power spectra obtained in the wind tunnel are similar, they are not generally an accurate representation of flight test experiments. Full article