Large eddy simulation of mechanical mixing in anaerobic digesters

A comprehensive study of anaerobic digestion requires an advanced turbulence model technique to accurately predict mixing flow patterns because the digestion process that involves mass transfer between anaerobes and their substrates is primarily dependent on detailed information about the fine structure of turbulence in the digesters. This study presents a large eddy simulation (LES) of mechanical agitation of non-Newtonian fluids in anaerobic digesters, in which the sliding mesh method is used to characterize the impeller rotation. The three subgrid scale (SGS) models investigated are: (i) Smagorinsky-Lilly model, (ii) wall-adapting local eddy-viscosity model, and (iii) kinetic energy transport (KET) model. The simulation results show that the three SGS models produce very similar flow fields. A comparison of the simulated and measured axial velocities indicates that the LES profile shapes are in general agreement with the experimental data but they differ markedly in velocity magnitudes. A check of impeller power and flow numbers demonstrates that all the SGS models give excellent predictions, with the KET model performing the best. Moreover, the performance of six Reynolds-averaged Navier-Stokes turbulence models are assessed and compared with the LES results.     

Optimization and scale‐up of oligonucleotide synthesis in packed bed reactors using computational fluid dynamics modeling

A computational fluid dynamics (CFD) model for the analysis of oligonucleotide synthesis in packed bed reactors was developed and used to optimize the scale up of the process. The model includes reaction kinetics data obtained under well defined conditions comparable to the situation in the packed bed. The model was validated in terms of flow conditions and reaction kinetics by comparison with experimental data. Experimental validation and the following model parameter studies by simulation were performed on the basis of a column with 0.3 g oligonucleotide capacity. The scale‐up studies based on CFD modelling were calculated on a 440 g scale (oligonucleotide capacity). © 2014 American Institute of Chemical Engineers Biotechnol. Prog., 30:1048–1056, 2014     

Cambrian cinctan echinoderms shed light on feeding in the ancestral deuterostome

Reconstructing the feeding mode of the latest common ancestor of deuterostomes is key to elucidating the early evolution of feeding in chordates and allied phyla; however, it is debated whether the ancestral deuterostome was a tentaculate feeder or a pharyngeal filter feeder. To address this, we evaluated the hydrodynamics of feeding in a group of fossil stem-group echinoderms (cinctans) using computational fluid dynamics. We simulated water flow past three-dimensional digital models of a Cambrian fossil cinctan in a range of possible life positions, adopting both passive tentacular feeding and active pharyngeal filter feeding. The results demonstrate that an orientation with the mouth facing downstream of the current was optimal for drag and lift reduction. Moreover, they show that there was almost no flow to the mouth and associated marginal groove under simulations of passive feeding, whereas considerable flow towards the animal was observed for active feeding, which would have enhanced the transport of suspended particles to the mouth. This strongly suggests that cinctans were active pharyngeal filter feeders, like modern enteropneust hemichordates and urochordates, indicating that the ancestral deuterostome employed a similar feeding strategy.     

Studies on fluid dynamics of the flow field and gas transfer in orbitally shaken tubes

Orbitally shaken cylindrical bioreactors [OrbShake bioreactors (OSRs)] without an impeller or sparger are increasingly being used for the suspension cultivation of mammalian cells. Among small volume OSRs, 50‐mL tubes with a ventilated cap (OSR50), originally derived from standard laboratory centrifuge tubes with a conical bottom, have found many applications including high‐throughput screening for the optimization of cell cultivation conditions. To better understand the fluid dynamics and gas transfer rates at the liquid surface in OSR50, we established a three‐dimensional simulation model of the unsteady liquid forms (waves) in this vessel. The studies verified that the operating conditions have a large effect on the interfacial surface. The volumetric mass transfer coefficient (k_La) was determined experimentally and from simulations under various working conditions. We also determined the liquid‐phase mass transfer coefficient (k_L) and the specific interfacial area (a) under different conditions to demonstrate that the value of a affected the gas transfer rate more than did the value of k_L. High oxygen transfer rates, sufficient for supporting the high‐density culture of mammalian cells, were found. Finally, the average axial velocity of the liquid was identified to be an important parameter for maintaining cells in suspension. Overall these studies provide valuable insights into the preferable operating conditions for the OSR50, such as those needed for cell cultures requiring high oxygen levels. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 33:192–200, 2017     

Evaluation of dentinal fluid flow behaviours: a fluid-structure interaction simulation

This study uses the fluid-structure interaction (FSI) method to investigate the fluid flow in dental pulp. First, the FSI method is used for the biomechanical simulation of dental intrapulpal responses during force loading (50, 100 and 150 N) on a tooth. The results are validated by comparison with experimental outcomes. Second, the FSI method is used to investigate an intact tooth subjected to a mechanical stimulus during loading at various loading rates. Force loading (0–100 N) is applied gradually to an intact tooth surface with loading rates of 125, 62.5, 25 and 12.5 N/s, respectively, and the fluid flow changes in the pulp are evaluated. FSI analysis is found to be suitable for examining intrapulpal biomechanics. An external force applied to a tooth with a low loading rate leads to a low fluid flow velocity in the pulp chamber, thus avoiding tooth pain.     

Early Archean planktonic mode of life: Implications from fluid dynamics of lenticular microfossils

Lenticular, and commonly flanged, microfossils in 3.0–3.4 Ga sedimentary deposits in Western Australia and South Africa are unusually large (20–80 μm across), robust, and widespread in space and time. To gain insight into the ecology of these organisms, we performed simulations of fluid dynamics of virtual cells mimicking lenticular forms of variable sizes, oblateness, flange presence, and flange thickness. Results demonstrate that (a) the flange reduces sedimentation velocity, (b) this flange function works more effectively in larger cells, and (c) modest oblateness lowers sedimentation rate. These observations support interpretations that the lenticular microbes were planktonic—a lifestyle that could have been advantageous in an early Earth harsh environment including violent volcanic activities, repeated asteroid impacts, and relatively high UV‐radiation. Although the robustness of these organisms could have provided additional protection on the early Earth, this architecture may have impeded a planktonic lifestyle by increasing cell density. However, our data suggest that this disadvantage could have been compensated by enlargement of cell volume, which could have enhanced the ability of the flange to slow sedimentation rate, especially if coupled with vacuolation. The results of this simulation study may help to explain the unique morphology and unusually large size of these Archean microfossils.     

Numerical simulation of blood flow in the left ventricle and aortic sinus using magnetic resonance imaging and computational fluid dynamics

Understanding cardiac blood flow patterns has many applications in analysing haemodynamics and for the clinical assessment of heart function. In this study, numerical simulations of blood flow in a patient-specific anatomical model of the left ventricle (LV) and the aortic sinus are presented. The realistic 3D geometry of both LV and aortic sinus is extracted from the processing of magnetic resonance imaging (MRI). Furthermore, motion of inner walls of LV and aortic sinus is obtained from cine-MR image analysis and is used as a constraint to a numerical computational fluid dynamics (CFD) model based on the moving boundary approach. Arbitrary Lagrangian–Eulerian finite element method formulation is used for the numerical solution of the transient dynamic equations of the fluid domain. Simulation results include detailed flow characteristics such as velocity, pressure and wall shear stress for the whole domain. The aortic outflow is compared with data obtained by phase-contrast MRI. Good agreement was found between simulation results and these measurements.     

Simulation of Fluid Flow in a Channel Induced by Three Types of Fin-Like Motion

One of many interesting research activities in biofluidmechanics is dedicated to investigations of locomotion in water. Some of propulsion mechanisms observed in the underwater world are used in the development process of underwater autonomic vehicles （AUV）. In order to characterise several solutions according to their manoeuvrability, influence on the surrounding fluid and energetic efficiency, a detailed analysis of fin-like movement is indispensable. In the current paper an analysis of undulatory, oscillatory and combined fin-like movements by means of numerical simulation is carried out. The conservation equation of mass and the conservation equation of momentum axe solved with the Finite Volume Method （FWM） by use of the software CFX-10.0. The undulatory and oscillatory fin movements axe modelled with an equation that is implemented within an additional subroutine and joined with the main solver. N carried out in the computational domain, in which one fin is fixed in a flow-through water duct. Simulations axe carded out in the range of the Re number up to 105. The results show significant influence of applied fin motion on the velocity distribution in the surrounding fluid.     

Slime travelers: Early evidence of animal mobility and feeding in an organic mat world

Mobility represents a key innovation in the evolution of complex animal life. The ability to move allows for the exploration of new food sources, escapes from unfavorable environmental conditions, enhanced ability to exchange genetic material, and is one of the major reasons for the diversity and success of animal life today. The oldest widely accepted trace fossils of animal mobility are found in Ediacaran‐aged rocks (635–539 Ma). The earliest definitive evidence for movement associated with exploitation of resources for feeding occurs in the White Sea assemblage of the Ediacara Biota—macroscopic, soft‐bodied fossils of Ediacaran age. Here, we evaluate potential support for mobility in dickinsoniomorphs, presenting new data regarding abundant Dickinsonia and associated trace fossils from the Ediacara Member, South Australia. Results quantitatively demonstrate that Dickinsonia was capable of mobility on relatively short, ecological timescales. This organism was bilaterally symmetrical, likely moved via muscular peristalsis, and left trace fossils due to active removal of the organic mat related to feeding. Analogous structures associated with Yorgia indicate that it was also mobile and fed in a similar manner. Morphological evidence suggests that two other modular taxa, Andiva and Spriggina, were able to move but did not feed in a manner that impacted the organic mat. Together, these data suggest that mobility was present in multiple disparate bilaterally symmetrical Ediacaran taxa.     

Simulated and experimental estimates of hydrodynamic drag from bio-logging tags

Drag force acting on swimming marine mammals is difficult to measure directly. Researchers often use simple modeling and kinematic measurements from animals, or computational fluid dynamics (CFD) simulations to estimate drag. However, studies that compare these methods are lacking. Here, computational simulation and physical experiments were used to estimate drag forces on gliding bottlenose dolphins (Tursiops truncatus). To facilitate comparison, variable drag loading (no-tag, tag, tag + 4, tag + 8) was used to increase force in both simulations and experiments. During the experiments, two dolphins were trained to perform controlled glides with variable loading. CFD simulations of dolphin/tag geometry in steady flow (1–6 m/s) were used to model drag forces. We expect both techniques will capture relative changes created by experimental conditions, but absolute forces predicted by the methods will differ. CFD estimates were within a calculated 90% confidence interval of the experimental results for all but the tag condition. Relative drag increase predicted by the simulation vs. experiment, respectively, differed by between 21% and 31%: tag, 4% vs. 33%; tag + 4, 47% vs. 68%; and tag + 8, 108% vs. 77%. The results from this work provide a direct comparison of computational and experimental estimates of drag, and provide a framework to quantify uncertainty.     

Comparative analysis of realistic CT-scan and simplified human airway models in airflow simulation

Efforts to model the human upper respiratory system have undergone many phases. Geometrical proximity to the realistic shape has been the subject of many research projects. In this study, three different geometries of the trachea and main bronchus were modelled, which were reconstructed from computed tomography (CT) scan images. The geometrical variations were named realistic, simplified and oversimplified. Realistic refers to the lifelike image taken from digital imaging and communications in medicine format CT scan images, simplified refers to the reconstructed image based on natural images without realistic details pertaining to the rough surfaces, and oversimplified describes the straight wall geometry of the airway. The characteristics of steady state flows with different flow rates were investigated, simulating three varied physical activities and passing through each model. The results agree with previous studies where simplified models are sufficient for providing comparable results for airflow in human airways. This work further suggests that, under most exercise conditions, the idealised oversimplified model is not favourable for simulating either airflow regimes or airflow with particle depositions. However, in terms of immediate analysis for the prediction of abnormalities of various dimensions of human airways, the oversimplified techniques may be used.     

Integration of mixing, heat transfer, and biochemical reaction kinetics in anaerobic methane fermentation

An extensive investigation of anaerobic methane fermentation requires identifying the relationship between the physical environment and biological process. In this study, a computational fluid dynamics (CFD) technique was used to characterize bacterial fermentation mechanisms intertwined with mixing and heat transfer in anaerobic digesters. The results demonstrate that the methane yield remains almost unchanged while the energy efficiency decreases with increasing mixing power in a complete‐mix digester, and that the energy output increases nonlinearly with the increase in heating energy in a plug‐flow digester. The CFD method can be applied to other bioreactors to gain valuable insights into their behavior as well. Integrating flow and temperature with kinetic behavior for anaerobic digestion not only solves the controversy about how mixing influences the digestive process, but also assists in optimizing the digester design and increasing the efficiency of energy conversion, and additionally, provides a reference for improving the mixing guidelines recommended by the U.S. Environmental Protection Agency. Biotechnol. Bioeng. 2012; 109: 2864–2874. © 2012 Wiley Periodicals, Inc.     

Comparative study of flow in right-sided and left-sided aortas: numerical simulations in patient-based models

A right-sided aorta is a rare malformation which may be associated with other various types of congenital heart disease. We utilised haemodynamic, echocardiographic measurements, computerised tomography and image reconstruction software packages that were integrated in a computational fluid dynamics model to determine blood flow patterns in patient-based aortas. In the left-sided aorta, a systolic clockwise rotational component was present, while helical flow was depicted in the aortic arch that was converted in the descending aorta as counter-rotating vortices with accompanying retrograde flow. The right-sided configuration has not altered the orientation of the three-dimensional vortex, but intensification of polymorphic flow patterns, alterations in wall shear stress distribution and development of a lateral pressure gradient at the area of an aneurysmal anomaly was observed. Moreover, increments of Reynolds, Womersley and Dean numbers were evident. These phenomena along with the formation of the aneurysm might influence cardiovascular risk in patients with right-sided aortas.     

Computational fluid dynamics modeling of the paddle dissolution apparatus: Agitation rate,mixing patterns,and fluid velocities

The purpose of this research was to further investigate the hydrodynamics of the United States Pharmacopeia (USP) paddle dissolution apparatus using a previously generated computational fluid dynamics (CFD) model. The influence of paddle rotational speed on the hydrodynamics in the dissolution vessel was simulated. The maximum velocity magnitude for axial and tangential velocities at different locations in the vessel was found to increase linearly with the paddle rotational speed. Path-lines of fluid mixing, which were examined from a central region at the base of the vessel, did not reveal a region of poor mixing between the upper cylin-drical and lower hemispherical volumes, as previously speculated. Considerable differences in the resulting flow patterns were observed for paddle rotational speeds between 25 and 150 rpm. The approximate time required to achieve complete mixing varied between 2 to 5 seconds at 150 rpm and 40 to 60 seconds at 25 rpm, although complete mixing was achievable for each speed examined. An analysis of CFD-generated velocities above the top surface of a cylindrical compact positioned at the base of the vessel, below the center of the rotating paddle, revealed that the fluid in this region was undergoing solid body rotation. An examination of the velocity boundary layers adjacent to the curved surface of the compact revealed large peaks in the shear rates for a region within∼3 mm from the base of the compact, consistent with a ‘grooving’ effect, which had been previously seen on the surface of compacts following dissolution, associated with a higher dissolution rate in this region.     

Flow modeling in a novel non-perfusion conical bioreactor

We have developed a bioreactor vessel design which has the advantages of simplicity and ease of assembly and disassembly, and with the appropriately determined flow rate, even allows for a scaffold to be suspended freely regardless of its weight. This article reports our experimental and numerical investigations to evaluate the performance of a newly developed non-perfusion conical bioreactor by visualizing the flow through scaffolds with 45 degrees and 90 degrees fiber lay down patterns. The experiments were conducted at the Reynolds numbers (Re) 121, 170, and 218 based on the local velocity and width of scaffolds. The flow fields were captured using short-time exposures of 60 microm particles suspended in the bioreactor and illuminated using a thin laser sheet. The effects of scaffold fiber lay down pattern and Reynolds number were obtained and correspondingly compared to results obtained from a computational fluid dynamics (CFD) software package. The objectives of this article are twofold: to investigate the hypothesis that there may be an insufficient exchange of medium within the interior of the scaffold when using our non-perfusion bioreactor, and second, to compare the flows within and around scaffolds of 45 degrees and 90 degrees fiber lay down patterns. Scaffold porosity was also found to influence flow patterns. It was therefore shown that fluidic transport could be achieved within scaffolds with our bioreactor design, being a non-perfusion vessel. Fluid velocities were generally same of the same or one order lower in magnitude as compared to the inlet flow velocity. Additionally, the 90 degrees fiber lay down pattern scaffold was found to allow for slightly higher fluid velocities within, as compared to the 45 degrees fiber lay down pattern scaffold. This was due to the architecture and pore arrangement of the 90 degrees fiber lay down pattern scaffold, which allows for fluid to flow directly through (channel-like flow).     

Numerical simulation of airflow and microparticle deposition in a synchrotron micro-CT-based pulmonary acinus model

The acinus consists of complex, branched alveolar ducts and numerous surrounding alveoli, and so in this study, we hypothesized that the particle deposition can be much influenced by the complex acinar geometry, and simulated the airflow and particle deposition (density = 1.0 g/cm³, diameter = 1 and 3 μm) numerically in a pulmonary acinar model based on synchrotron micro-CT of the mammalian lung. We assumed that the fluid–structure interaction was neglected and that alveolar flow was induced by the expansion and contraction of the acinar model with the volume changing sinusoidally with time as the moving boundary conditions. The alveolar flow was dominated by radial flows, and a weak recirculating flow was observed at the proximal side of alveoli during the entire respiratory cycle, despite the maximum Reynolds number at the inlet being 0.029. Under zero gravity, the particle deposition rate after single breathing was less than 0.01, although the particles were transported deeply into the acinus after inspiration. Under a gravitational field, the deposition rate and map were influenced strongly by gravity orientation. In the case of a particle diameter of 1 μm, the rate increased dramatically and mostly non-deposited particles remained in the model, indicating that the rate would increase further after repeated breathing. At a particle diameter of 3 μm, the rate was 1.0 and all particles were deposited during single breathing. Our results show that the particle deposition rate in realistic pulmonary acinar model is higher than in an idealized model.     

Protein aggregation and mitigation strategy in low pH viral inactivation for monoclonal antibody purification

ABSTRACT

Significant amounts of soluble product aggregates were observed during low-pH viral inactivation (VI) scale-up for an IgG4 monoclonal antibody (mAb IgG4-N1), while small-scale experiments in the same condition showed negligible aggregation. Poor mixing and product exposure to low pH were identified as the root cause. To gain a mechanistic understanding of the problem, protein aggregation properties were studied by varying critical parameters including pH, hold time and protein concentration. Comprehensive biophysical characterization of product monomers and aggregates was performed using fluorescence-size-exclusion chromatography, differential scanning fluorimetry, fluorescence spectroscopy, and dynamic light scattering. Results showed IgG4-N1 partially unfolds at about pH 3.3 where the product molecules still exist largely as monomers owing to strong inter-molecular repulsions and favorable colloidal stability. In the subsequent neutralization step, however, the conformationally changed monomers are prone to aggregation due to weaker inter-molecular repulsions following the pH transition from 3.3 to 5.5. Surface charge calculations using homology modeling suggested that intra-molecular repulsions, especially between CH2 domains, may contribute to the IgG4-N1 unfolding at ≤ pH 3.3. Computational fluid dynamics (CFD) modeling was employed to simulate the conditions of pH titration to reduce the risk of aggregate formation. The low-pH zones during acid addition were characterized using CFD modeling and correlated to the condition causing severe product aggregation. The CFD tool integrated with the mAb solution properties was used to optimize the VI operating parameters for successful scale-up demonstration. Our research revealed the governing aggregation mechanism for IgG4-N1 under acidic conditions by linking its molecular properties and various process-related parameters to macroscopic aggregation phenomena. This study also provides useful insights into the cause and mitigation of low-pH-induced IgG4 aggregation in downstream VI operation.     

Computational fluid dynamics as a tool for testing functional and ecological hypotheses in fossil taxa

Computational fluid dynamics is a method for simulating fluid flows that has been widely used in engineering for decades, and which also has applications for studying function and ecology in fossil taxa. However, despite the possible benefits of this approach, computational fluid dynamics has been used only rarely in palaeontology to date. The theoretical basis underlying the technique is outlined and the main steps involved in carrying out computer simulations of fluid flows are detailed. I also describe previous studies that have applied the method to fossils and discuss their potential for informing future research directions in palaeontology. Computational fluid dynamics can enable large‐scale comparative analyses, as well as exacting tests of hypotheses related to the function and ecology of ancient organisms. In this way, it could transform our understanding of many extinct fossil groups.