Search

Numerical simulations of pulsatile flow in a feline aorta for hypertrophic cardiomyopathy (HCM) heart conditions are performed to predict flow details and to evaluate possible thrombus trajectory patterns. The study demonstrates that average flow rate boundary conditions (QBC) at artery outlets act as a resistance-type boundary condition for pulsatile flow. For simulations when the exact artery outflows are not known, specification of estimated values from physiological conditions is a plausible approach. This boundary condition is further improved using an iterative method (I-QBC) to accurately satisfy outflow conditions when expected outflow is known. The approach is validated against experimental data for the prediction of iliac artery flow and wall stresses in a human abdominal aorta. The feline aorta simulations including Lagrangian particle transport are performed on grids with up to 11M cells for a generalized feline aorta. It is found that LES on larger grids performs significantly better than URANS for the prediction of vortical structures. Simulations for both healthy and HCM conditions show similar flow patterns in the upper abdominal aorta. However, the HCM condition shows the presence of large recirculation regions in the thoracic aorta resulting in 50% lower flow through the iliac arteries and increased entrapment of fluid-borne particles near the trifurcation region compared to the healthy condition.

A verification and validation study was performed using the open source computational fluid dynamics solver OpenFOAM version 2.0.0 for incompressible bluff body fluid flows. This includes flow over a backward facing step, a sphere in the subcritical regime, and delta wing with sharp leading edge. The study investigates solver scalability, and accuracy of numerical methods and turbulence models available in the solver. Grid verification study shows mostly monotonic convergence with averaged grid uncertainty

Robust automated vortex detection algorithms are needed to facilitate the exploration of large-scale turbulent fluid flow simulations. Unfortunately, robust non-local vortex detection algorithms are computationally intractable for large data sets and local algorithms, while computationally tractable, lack robustness. We argue that the deficiencies inherent to the local definitions occur because of two fundamental issues: the lack of a rigorous definition of a vortex and the fact that a vortex is an intrinsically non-local phenomenon. As a first step towards addressing this problem, we demonstrate the use of machine learning techniques to enhance the robustness of local vortex detection algorithms. We motivate the presence of an expert-in-the-loop using empirical results based on machine learning techniques. We employ adaptive boosting to combine a suite of widely used, local vortex detection algorithms, which we term weak classifiers, into a robust compound classifier. Fundamentally, the training phase of the algorithm, in which an expert manually labels small, spatially contiguous regions of the data, incorporates non-local information into the resulting compound classifier. We demonstrate the efficacy of our approach by applying the compound classifier to two data sets obtained from computational fluid dynamical simulations. Our results demonstrate that the compound classifier has a reduced misclassification rate relative to the component classifiers.

The performance of a pseudo-spectral solver, which uses FFT in the streamwise and spanwise directions and Chebyshev polynomials in the wall-normal direction, is improved by implementing MPI domain decomposition for the Chebyshev collocation method using the influence matrix method. The parallel code is validated for LES of plane channel flow at Reτ=590 and 2003, and applied for temporally developing flat plate transition simulations. The domain decomposition using the influence matrix method is implemented to achieve C1 continuity of the velocity across the domain interface. The method provides super linear scalability of solver computational performance. The predictions show kinks in the turbulent statistics on coarser grids, but agree very well with DNS data for finer grids. It is found that accurate turbulent predictions require a minimum of 48 collocation points per subdomain for Re ~ 105, with higher numbers required for Re ~ 106. The flat plate transitional flow simulation is performed for a...

Predictive capabilities of hybrid RANS/LES models are compared for single-phase flow over a surface combatant at Re = 5.3 × 106 and an appended DARPA Suboff model at Re = 1.2 × 107. The turbulence models used in the study are: k–ω shear stress transport (SST)-URANS; Spalart Allmaras based detached eddy simulation (SA-DDES); k–ω based improved delayed detached eddy simulation (KW-IDDES); and a dynamic hybrid RANS/LES (DHRL) model coupling SST and implicit LES. For the surface combatant case, both SA-DDES and KW-IDDES predicted 100% and 28%, respectively. Overall, the DHRL model performed best among the turbulence models tested, and KW-IDDES performed worst. The study indicates that the DHRL approach has the potential to provide accurate mean flow predictions while resolving small-scale turbulent structures. Results also highlight the importance of the wall function formulation for accurately resolving mean skin friction coefficient, especially over smooth regions of the hull.

This paper investigates a recently proposed dynamic hybrid RANS-LES framework using a general-purpose finite-volume flow solver. The new method is highly generalized, allowing coupling of any selected RANS model with any selected LES model and containing no explicit grid dependence in its formulation. Selected results are presented for three test cases: two-dimensional channel flow, backward facing step, and a nozzle flow relevant to biomedical applications. Comparison with experimental and DNS data, and with other hybrid RANS-LES approaches, highlights the advantages of the new method and suggests that further investigation is warranted.

In large-eddy simulation, the subgrid-stress model is included to represent the energy-transfer mechanisms between the resolved and unresolved scales of motion. Although most models represent forward scatter (dissipation) of energy, they fail to account for backscatter properly. The algebraic model, obtained directly from the closure of the subgrid stresses, has an additional term over the mixed model that accounts for the backscatter of energy explicitly. This model has been previously documented in the open literature and shows promising results in accounting for the subgrid-scale energy-transfer mechanisms. Detailed documentation is provided of the backscatter capability of the algebraic model, for a series of test cases including isotropic turbulence, plane channel flow, and two free shear flows. The model performace in predicting turbulent flow quantities is also compared with commonly used subgrid-stress models and with the available direct-numerical-simulation data.

One of the key challenges for computational fluid dynamics (CFD) simulations of human lung airflow is the sheer size and complexity of the complete, multiscale geometry of the bronchopulmonary tree. Since 3-D CFD simulations of the full airway tree are currently intractable, researchers have proposed reduced geometry models in which multiple airway paths are truncated downstream of the first few generations. This paper investigates a recently proposed method for closing the CFD model by application of physiologically correct boundary conditions at truncated outlets. A realistic, reduced geometry model of the lung airway based on CT data has been constructed up to generation 18, including extrathoracic, bronchi, and bronchiole regions. Results indicate that the new method yields reasonable results for pressure drop through the airway, at a small fraction of the cost of fully resolved simulations.

Subjects with active stretch reflexes responded to an imposed sinusoidal movement of the ankle joint with a reflex force whose amplitude and timing varied widely with changes in the frequency of movement. At some frequency between 6 and 8 Hz, the reflex force tended to offset the non-reflex component of resistance, and thus to reduce the total resistance to movement. At this frequency the reflex response was particularly vigorous, with a deep modulation of electromyogram (e.m.g.) activity and a displacement of the joint stiffness vectors far from their high frequency values. The total resistance to movement might then be small, or it might be zero, or the reflex might actually assist the movement. As the frequency of movement was decreased through this critical range, the timing of the reflex response to movement changed rapidly with an abrupt advancement of the triceps surae e.m.g. signal, and a wide separation of the joint stiffness vectors as they passed close to the origin. This result was attributed to a changing distribution of the movement between the muscle fibres and an elastic Achilles tendon. It was assumed that at most frequencies the muscle fibres resisted extension, so that a major part of the imposed movement went into stretching the tendon; when, however, at 6-8 Hz, the reflex response was so timed as to reduce or abolish the resistance of the muscle fibres, more of the movement would take place in them. The muscle spindles would 'see' this larger movement of the muscle fibres, and generate correspondingly more reflex activity. A simplified model of the muscle-tendon combination behaves in a way that supports this view, and the available information about the human Achilles tendon indicates that it is sufficiently compliant for such an explanation. Therefore, movements imposed on the ankle joint would not necessarily be 'seen' by the muscle spindles, since they would be modified by transmission through a compliant tendon. By assuming a value for the tendon stiffness, it was possible to calculate the course of movements that actually occurred in the muscle fibres and spindles. Records of these spindle movements indicated how some non-linearities might arise.(ABSTRACT TRUNCATED AT 400 WORDS)

Ankle joints were subjected to sinusoidal movements at a range of amplitudes and frequencies. Records were made of electromyograms (e.m.g.s) in calf muscles, and of the forces at the joints. When the leg is relaxed, the ankle joint resists an imposed sinusoidal movement with a small approximately sinusoidal force. It is stiffer in its resistance to small movements than to large ones, and this resistance is greater when the joint is dorsiflexed than when it is plantarflexed. If the subject exerts a steady mean flexing force, the imposed sinusoidal movement generates reflex activity which may be recorded as a modulation of the gastrocnemius and soleus e.m.g.s. The e.m.g. response to the sinusoidal movement occurs later in cycles of movement at high than at low frequencies, as one could expect of a reflex pathway that involves a delay. The results suggest that this delay is between 50 and 60 ms, and we conclude that under these circumstances spinal stretch reflexes are playing the important part. The relation of the resisting force to the movement has been displayed as a vector. As the frequency changes, this vector describes the circular path that is characteristic of a system which includes delays or lags; this path enables one to draw conclusions about the amplitude and timing of the reflex resistance to the movement. When a subject exerts a moderate flexing force against the sinusoidal movement for some minutes, the reflex response becomes progressively potentiated. A subject whose reflex responses are normally slight may then exhibit a vigorous reflex response to the movement of that ankle. This enhancement of spinal reflex activity was accompanied by an increase in the myotatic reflex response at the ankle. Reflex responses to sinusoidal movement were most clearly seen when the subject exerted a mean flexing force that amounted to about one-fifth of his maximum. Very small movements (+/- 0.5 degrees) generated little or no reflex response. With large amplitudes of movement there was more reflex activity, but at some amplitude (which varied from subject to subject and from time to time) the reflex mechanism appeared to 'saturate', and further increases in amplitude were not accompanied by comparable increases in the reflex response. With movements at 10-15 Hz the e.m.g. response often became large in alternate cycles, with less activity in the intervening cycles.(ABSTRACT TRUNCATED AT 400 WORDS)