A Hybrid Reactive Multiphasic Mixture with a Compressible Fluid Solvent

Author(s):  
Jay J. Shim ◽  
Gerard A. Ateshian

Abstract Mixture theory is a general framework that has been used to model mixtures of solid, fluid, and solute constituents, leading to significant advances in modeling the mechanics of biological tissues and cells. Though versatile and applicable to a wide range of problems in biomechanics and biophysics, standard multiphasic mixture frameworks incorporate neither dynamics of viscous fluids nor fluid compressibility, both of which facilitate the finite element implementation of computational fluid dynamics solvers. This study formulates governing equations for reactive multiphasic mixtures where the interstitial fluid has a solvent which is viscous and compressible. This hybrid reactive multiphasic framework uses state variables that include the deformation gradient of the porous solid matrix, the volumetric strain and rate of deformation of the solvent, the solute concentrations, and the relative velocities between the various constituents. Unlike standard formulations which employ a Lagrange multiplier to model fluid pressure, this framework requires the formulation of a function of state for the pressure, which depends on solvent volumetric strain and solute concentrations. Under isothermal conditions the formulation shows that the solvent volumetric strain remains continuous across interfaces between hybrid multiphasic domains. Apart from the Lagrange multiplier-state function distinction for the fluid pressure, and the ability to accommodate viscous fluid dynamics, this hybrid multiphasic framework remains fully consistent with standard multiphasic formulations previously employed in biomechanics. With these additional features, the hybrid multiphasic mixture theory makes it possible to address a wider range of problems that are important in biomechanics and mechanobiology.

2020 ◽  
Vol 17 (170) ◽  
pp. 20200612
Author(s):  
A. Andreozzi ◽  
M. Iasiello ◽  
P. A. Netti

Macromolecules and drug delivery to solid tumours is strongly influenced by fluid flow through interstitium, and pressure-induced tissue deformations can have a role in this. Recently, it has been shown that temperature-induced tissue deformation can influence interstitial fluid velocity and pressure fields, too. In this paper, the effect of modulating-heat strategies to influence interstitial fluid transport in tissues is analysed. The whole tumour tissue is modelled as a deformable porous material, where the solid phase is made up of the extracellular matrix and cells, while the fluid phase is the interstitial fluid that moves through the solid matrix driven by the fluid pressure gradient and vascular capillaries that are modelled as a uniformly interspersed fluid point-source. Pulsating-heat generation is modelled with a time-variable cosine function starting from a direct current approach to solve the voltage equation, for different pulsations. From the steady-state solution, a step-variation of vascular pressure included in the model equation as a mass source term via the Starling equation is simulated. Dimensionless 1D radial equations are numerically solved with a finite-element scheme. Results are presented in terms of temperature, volumetric strain, pressure and velocity profiles under different conditions. It is shown that a modulating-heat procedure influences velocity fields, that might have a consequence in terms of mass transport for macromolecules or drug delivery.


2021 ◽  
pp. 108128652110545
Author(s):  
S Kiana Naghibzadeh ◽  
Noel Walkington ◽  
Kaushik Dayal

Accretion and ablation, i.e., the addition and removal of mass at the surface, are important in a wide range of physical processes, including solidification, growth of biological tissues, environmental processes, and additive manufacturing. The description of accretion requires the addition of new continuum particles to the body, and is therefore challenging for standard continuum formulations for solids that require a reference configuration. Recent work has proposed an Eulerian approach to this problem, enabling side-stepping of the issue of constructing the reference configuration. However, this raises the complementary challenge of determining the stress response of the solid, which typically requires the deformation gradient that is not immediately available in the Eulerian formulation. To resolve this, the approach introduced the elastic deformation as an additional kinematic descriptor of the added material, and its evolution has been shown to be governed by a transport equation. In this work, the method of characteristics is applied to solve concrete simplified problems motivated by biomechanics and manufacturing. Specifically, (1) for a problem with both ablation and accretion in a fixed domain and (2) for a problem with a time-varying domain, the closed-form solution is obtained in the Eulerian framework using the method of characteristics without explicit construction of the reference configuration.


Processes ◽  
2020 ◽  
Vol 9 (1) ◽  
pp. 79
Author(s):  
Minghan Luo ◽  
Wenjie Xu ◽  
Xiaorong Kang ◽  
Keqiang Ding ◽  
Taeseop Jeong

The ultraviolet photochemical degradation process is widely recognized as a low-cost, environmentally friendly, and sustainable technology for water treatment. This study integrated computational fluid dynamics (CFD) and a photoreactive kinetic model to investigate the effects of flow characteristics on the contaminant degradation performance of a rotating annular photoreactor with a vacuum-UV (VUV)/UV process performed in continuous flow mode. The results demonstrated that the introduced fluid remained in intensive rotational movement inside the reactor for a wide range of inflow rates, and the rotational movement was enhanced with increasing influent speed within the studied velocity range. The CFD modeling results were consistent with the experimental abatement of methylene blue (MB), although the model slightly overestimated MB degradation because it did not fully account for the consumption of OH radicals from byproducts generated in the MB decomposition processes. The OH radical generation and contaminant degradation efficiency of the VUV/UV process showed strong correlation with the mixing level in a photoreactor, which confirmed the promising potential of the developed rotating annular VUV reactor in water treatment.


2013 ◽  
Vol 80 (2) ◽  
Author(s):  
Minh H. Tran ◽  
Younane N. Abousleiman

The porochemoelectroelastic analytical models have been used to describe the response of chemically active and electrically charged saturated porous media such as clay soils, shales, and biological tissues. However, existing studies have ignored the anisotropic nature commonly observed on these porous media. In this work, the anisotropic porochemoelectroelastic theory is presented. Then, the solution for an inclined wellbore drilled in transversely isotropic shale formations subjected to anisotropic far-field stresses with time-dependent down-hole fluid pressure and fluid activity is derived. Numerical examples illustrating the combined effects of porochemoelectroelastic behavior and anisotropy on wellbore responses are also included. The analysis shows that ignoring either the porochemoelectroelastic effects or the formation anisotropy leads to inaccurate prediction of the near-wellbore pore pressure and effective stress distributions. Finally, wellbore responses during a leak-off test conducted soon after drilling are analyzed to demonstrate the versatility of the solution in simulating complex down-hole conditions.


Author(s):  
Gerard A. Ateshian ◽  
Brandon Zimmerman

Abstract Mixture theory models continua consisting of multiple constituents with independent motions. In constrained mixtures all constituents share the same velocity but they may have different reference configurations. The theory of constrained reactive mixtures was formulated to analyze growth and remodeling in living biological tissues. It can also reproduce and extend classical frameworks of damage mechanics and viscoelasticity under isothermal conditions, when modeling bonds that can break and reform. This study focuses on establishing the thermodynamic foundations of constrained reactive mixtures under more general conditions, for arbitrary reactive processes where temperature varies in time and space. By incorporating general expressions for reaction kinetics, it is shown that the residual dissipation statement of the Clausius-Duhem inequality must include a reactive power density, while the axiom of energy balance must include a reactive heat supply density. Both of these functions are proportional to the molar production rate of a reaction, and they depend on the chemical potentials of the mixture constituents. We present novel formulas for the classical thermodynamic concepts of energy of formation and heat of reaction, making it possible to evaluate the heat supply generated by reactive processes from the knowledge of the specific free energy of mixture constituents as well as the reaction rate. We illustrate these novel concepts with mixtures of ideal gases, and isothermal reactive damage mechanics and viscoelasticity, as well as reactive thermoelasticity. This framework facilitates the analysis of reactive tissue biomechanics and physiological and biomedical engineering processes where temperature variations cannot be neglected.


1970 ◽  
Vol 7 ◽  
pp. 60-64 ◽  
Author(s):  
Ruchi Khare ◽  
Vishnu Prasad Prasad ◽  
Sushil Kumar

The testing of physical turbine models is costly, time consuming and subject to limitations of laboratory setup to meet International Electro technical Commission (IEC) standards. Computational fluid dynamics (CFD) has emerged as a powerful tool for funding numerical solutions of wide range of flow equations whose analytical solutions are not feasible. CFD also minimizes the requirement of model testing. The present work deals with simulation of 3D flow in mixed flow (Francis) turbine passage; i.e., stay vane, guide vane, runner and draft tube using ANSYS CFX 10 software for study of flow pattern within turbine space and computation of various losses and efficiency at different operating regimes. The computed values and variation of performance parameters are found to bear close comparison with experimental results.Key words: Hydraulic turbine; Performance; Computational fluid dynamics; Efficiency; LossesDOI: 10.3126/hn.v7i0.4239Hydro Nepal Journal of Water, Energy and EnvironmentVol. 7, July, 2010Page: 60-64Uploaded date: 31 January, 2011


2018 ◽  
Vol 24 (4) ◽  
pp. 369-378 ◽  
Author(s):  
J.N.M. Batista ◽  
R.C. Brito ◽  
R. Béttega

The spouted bed presents limitations in terms of scale-up. Furthermore, its stability depends on its geometry as well as the properties of the fluid and solid phases. CFD provides an important tool to improve understanding of these aspects, enabling a wide range of information to be obtained rapidly and at low cost. In this work, CFD simulation was used to evaluate the effects of different inlet air distributors (Venturi and straight tube) and the effects of static bed height on the fluid and solid dynamics of a conical spouted bed. Simulations were performed using the two-dimensional Euler-Euler approach. In order to evaluate the fluid dynamics model, static pressure data obtained by simulation were compared with experimental data obtained with the Venturi distributor. The fluid and solid dynamics of the conical spouted bed were obtained by CFD simulation. The results showed that the pressure drop was lower for the straight tube air distributor, while the Venturi air distributor provided higher stability and a more homogenous air distribution at the bed entrance.


2021 ◽  
Vol 56 (5) ◽  
pp. 294-317
Author(s):  
A. I. Siswantara ◽  
H. Pujowidodo ◽  
M. A. Budiyanto ◽  
G. G. Ramdlan Gunadi ◽  
C. D. Widiawaty

This research aims to find the optimal standard k-e turbulence model constants (cµ, c1e, and c2e) for better predicting compressible fluid dynamics in an air jet ejector. The turbulence field in a jet flow plays an important role in influencing the performance of the momentum transfer process at a shear layer in nozzle application for momentum source and mixing process. In this research, some activities have been done before analyzing and optimizing the turbulence model constants, including preliminary turbulence modeling study for compressible flow in the air-jet ejector, verification, and validation with primary experimental data as well as by other secondary data. The preliminary studies in turbulence modeling presented that the turbulence modeling of a 3mm air jet-ejector resulted in a similar trend of the relation between entrainment ratio and motive fluid pressure. The results showed that the sensitive parameters in the standard k-emodel dissipation and diffusion terms, cµ, c1e, and c2e, strongly affected the optimum value of turbulence kinetic energy (k) and dissipation rate (e), compared to the reference model. Better k and e could be obtained by changing the c2e into positively proportional, but the cµ and c1e must be changed with opposite proportionality. It was found that the optimum standard k-e model constants in the case of air-jet ejector with 3 mm nozzle diameter for cµ, c1e, and c2e are 0.05, 1.48, and 1.88, respectively, with the error values for k being -8.88% and e being -17.44%.


Author(s):  
Andrea G. Sanvito ◽  
Giacomo Persico ◽  
M. Sergio Campobasso

Abstract This study provides a novel contribution toward the establishment of a new high-fidelity simulation-based design methodology for stall-regulated horizontal axis wind turbines. The aerodynamic design of these machines is complex, due to the difficulty of reliably predicting stall onset and poststall characteristics. Low-fidelity design methods, widely used in industry, are computationally efficient, but are often affected by significant uncertainty. Conversely, Navier–Stokes computational fluid dynamics (CFD) can reduce such uncertainty, resulting in lower development costs by reducing the need of field testing of designs not fit for purpose. Here, the compressible CFD research code COSA is used to assess the performance of two alternative designs of a 13-m stall-regulated rotor over a wide range of operating conditions. Validation of the numerical methodology is based on thorough comparisons of novel simulations and measured data of the National Renewable Energy Laboratory (NREL) phase VI turbine rotor, and one of the two industrial rotor designs. An excellent agreement is found in all cases. All simulations of the two industrial rotors are time-dependent, to capture the unsteadiness associated with stall which occurs at most wind speeds. The two designs are cross-compared, with emphasis on the different stall patterns resulting from particular design choices. The key novelty of this work is the CFD-based assessment of the correlation among turbine power, blade aerodynamics, and blade design variables (airfoil geometry, blade planform, and twist) over most operational wind speeds.


Sign in / Sign up

Export Citation Format

Share Document