Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA Second Order Thermal and Momentum Immersed Boundary Method for Conjugate Heat Transfer in a Cartesian Finite Volume Solver
136
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
A conjugate heat transfer (CHT) immersed boundary (IB and CHTIB) method is developed for use with laminar and turbulent flows with low to moderate Reynolds numbers. The method is validated with the canonical flow of two co-annular rotating cylinders at $Re=50$ which shows second order accuracy of the $L_{2}$ and $L_{infty}$ error norms of the temperature field over a wide rage of solid to fluid thermal conductivities, $kappa_{s}/kappa_{f} = left(9-100right)$. To evaluate the CHTIBM with turbulent flow a fully developed, heated, turbulent channel $left(Re_{u_{tau}}=150text{ and } kappa_{s}/kappa_{f}=4 right)$ is used which shows near perfect correlation to previous direct numerical simulation (DNS) results. The CHTIB method is paired with a momentum IB method (IBM), both of which use a level set field to define the wetted boundaries of the fluid/solid interfaces and are applied to the flow solver implicitly with rescaling of the difference operators of the finite volume (FV) method (FVM).
rate research
Read More
We present an effective thermal open boundary condition for convective heat transfer problems on domains involving outflow/open boundaries. This boundary condition is energy-stable, and it ensures that the contribution of the open boundary will not cause an ``energy-like temperature functional to increase over time, irrespective of the state of flow on the open boundary. It is effective in coping with thermal open boundaries even in flow regimes where strong vortices or backflows are prevalent on such boundaries, and it is straightforward to implement. Extensive numerical simulations are presented to demonstrate the stability and effectiveness of our method for heat transfer problems with strong vortices and backflows occurring on the open boundaries. Simulation results are compared with previous works to demonstrate the accuracy of the presented method.
The study of viscous fluid flow coupled with rigid or deformable solids has many applications in biological and engineering problems, e.g., blood cell transport, drug delivery, and particulate flow. We developed a partitioned approach to solve this coupled Multiphysics problem. The fluid motion was solved by Palabos (Parallel Lattice Boltzmann Solver), while the solid displacement and deformation was simulated by LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator). The coupling was achieved through the immersed boundary method (IBM). The code modeled both rigid and deformable solids exposed to flow. The code was validated with the Jeffery orbits of an ellipsoid particle in shear flow, red blood cell stretching test, and effective blood viscosity flowing in tubes. It demonstrated essentially linear scaling from 512 to 8192 cores for both strong and weak scaling cases. The computing time for the coupling increased with the solid fraction. An example of the fluid-solid coupling was given for flexible filaments (drug carriers) transport in a flowing blood cell suspensions, highlighting the advantages and capabilities of the developed code.
We extend the recently introduced divergence-conforming immersed boundary (DCIB) method [1] to fluid-structure interaction (FSI) problems involving closed co-dimension one solids. We focus on capsules and vesicles, whose discretization is particularly challenging due to the higher-order derivatives that appear in their formulations. In two-dimensional settings, we employ cubic B-splines with periodic knot vectors to obtain discretizations of closed curves with C^2 inter-element continuity. In three-dimensional settings, we use analysis-suitable bi-cubic T-splines to obtain discretizations of closed surfaces with at least C^1 inter-element continuity. Large spurious changes of the fluid volume inside closed co-dimension one solids is a well-known issue for IB methods. The DCIB method results in volume changes orders of magnitude lower than conventional IB methods. This is a byproduct of discretizing the velocity-pressure pair with divergence-conforming B-splines, which lead to negligible incompressibility errors at the Eulerian level. The higher inter-element continuity of divergence-conforming B-splines is also crucial to avoid the quadrature/interpolation errors of IB methods becoming the dominant discretization error. Benchmark and application problems of vesicle and capsule dynamics are solved, including mesh-independence studies and comparisons with other numerical methods.
The dissolution of solids has created spectacular geomorphologies ranging from centimeter-scale cave scallops to the kilometer-scale stone forests of China and Madagascar. Mathematically, dissolution processes are modeled by a Stefan problem, which describes how the motion of a phase-separating interface depends on local concentration gradients, coupled to a fluid flow. Simulating these problems is challenging, requiring the evolution of a free interface whose motion depends on the normal derivatives of an external field in an ever-changing domain. Moreover, density differences created in the fluid domain induce self-generated convecting flows that further complicate the numerical study of dissolution processes. In this contribution, we present a numerical method for the simulation of the Stefan problem coupled to a fluid flow. The scheme uses the Immersed Boundary Smooth Extension method to solve the bulk advection-diffusion and fluid equations in the complex, evolving geometry, coupled to a {theta}-L scheme that provides stable evolution of the boundary. We demonstrate third-order temporal and pointwise spatial convergence of the scheme for the classical Stefan problem, and second-order temporal and pointwise spatial convergence when coupled to flow. Examples of dissolution of solids that result in high-Rayleigh number convection are numerically studied, and qualitatively reproduce the complex morphologies observed in recent experiments.
We present a novel moving immersed boundary method (IBM) and employ it in direct numerical simulations (DNS) of the closed-vessel swirling von Karman flow in laminar and turbulent regimes. The IBM extends direct-forcing approaches by leveraging a time integration scheme, that embeds the immersed boundary forcing step within a semi-implicit iterative Crank-Nicolson scheme. The overall method is robust, stable, and yields excellent results in canonical cases with static and moving boundaries. The moving IBM allows us to reproduce the geometry and parameters of the swirling von Karman flow experiments in (F. Ravelet, A. Chiffaudel, and F. Daviaud, JFM 601, 339 (2008)) on a Cartesian grid. In these DNS, the flow is driven by two-counter rotating impellers fitted with curved inertial stirrers. We analyze the transition from laminar to turbulent flow by increasing the rotation rate of the counter-rotating impellers to attain the four Reynolds numbers 90, 360, 2000, and 4000. In the laminar regime at Reynolds number 90 and 360, we observe flow features similar to those reported in the experiments and in particular, the appearance of a symmetry-breaking instability at Reynolds number 360. We observe transitional turbulence at Reynolds number 2000. Fully developed turbulence is achieved at Reynolds number 4000. Non-dimensional torque computed from simulations matches correlations from experimental data. The low Reynolds number symmetries, lost with increasing Reynolds number, are recovered in the mean flow in the fully developed turbulent regime, where we observe two tori symmetrical about the mid-height plane. We note that turbulent fluctuations in the central region of the device remain anisotropic even at the highest Reynolds number 4000, suggesting that isotropization requires significantly higher Reynolds numbers.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
