Boundary Conditions Research Articles

Impact of realistic boundary conditions on CFD simulations: A case study of vehicle ventilation

Impact of realistic boundary conditions on CFD simulations: A case study of vehicle ventilation

Fast and Efficient Mitigation of Coupled Axial-Torsional Drillstring Vibrations through State-Feedback Linear-Quadratic-Integral (LQI) Control

Summary In the exploitation of subsurface hydrocarbon and geothermal energy, more than 50% of total expenses are typically attributed to drilling and well completion. Coupled triaxial vibrations in the drilling system, excited by various mechanical, hydraulic, and geological sources, are major causes of premature downhole tool failures, excessive bit damage, compromised hole quality, reduced rate of penetration (ROP), and increased nonproductive time. Drillstring dynamic models, shock absorption tools, and controllers have been developed since the 1960s to understand and mitigate detrimental vibrations in the drilling system, aiming for smoother and more efficient operations. With advancements in measurement/logging while drilling technologies, rig equipment, mud pulse/wired pipe telemetry, and the digitalization of drilling engineering, a foundation has been established for proactive drilling vibration detection, analysis, and control to achieve efficient and safe drilling. In this study, a reduced-order axial-torsional drillstring dynamic model with soft string lateral contacts is first developed. Nonlinear drilling boundary conditions and inputs, such as bit-rock interactions (BRIs) with torque-on-bit (TOB) velocity-weakening effects, wellbore-drillstring contacts/friction, drillstring gravity, mud buoyancy, and more, are defined within the modeling framework. Equilibrium states and linearization of the model are derived for any given weight-on-bit (WOB) and rotational speed (RPM) setpoints, 3D well trajectories, and drillstring dimensions. A linear-quadratic-integral (LQI) controller is developed and tested for drillstring vibration suppression using a 6,562 ft (2000 m) drillstring model. In comparison to commercially available controllers for stick/slip mitigation, which are designed based on impedance matching and only address the decoupled drillstring torsional dynamics, the proposed state-feedback LQI controller emphasizes the coupling between the axial and torsional drillstring dynamics by simultaneously controlling WOB and RPM. Simulation results demonstrate that the LQI controller can suppress stick/slip and stabilize downhole WOB within 7 seconds, while a commercial drilling controller takes more than 35 seconds to achieve the same. Additionally, valuable insights and potential future directions for nonlinear drillstring vibration mitigation and drilling controller design are provided. With its fast drilling dysfunction stabilization time, small downhole RPM/WOB overshoots, minimal rig control inputs, robustness against drilling nonlinearities, and high efficiency, the proposed multi-input-multi-output state-feedback drilling controller holds great potential for application in proactive drilling vibration control, drilling automation, and real-time drilling optimization.

Heat Transfer and Entropy Generation in Vibrational Flow: Newtonian vs. Inelastic Non-Newtonian Fluid

A computational method is employed to solve heat transfer and entropy generation within a circular pipe. The thermal boundary condition assumes a constant wall temperature, while viscosity is taken to be dependent on temperature. A power-law type shear-thinning fluid is utilized in the analysis, with sinusoidal vibration applied horizontally perpendicular to the flow direction. Temperature distributions across the pipe are illustrated. Additionally, the entropy generation rate over the entire fluid volume under vibration was examined, comparing the results between steady flow and vibrational flow for both types of fluids. It was found that radial mixing is more pronounced in non-Newtonian fluids as vibration increases the strain rate, which is higher for low Reynolds numbers. The research provides a quantitative analysis of heat transfer and entropy generation for both Newtonian and shear-thinning fluids at different Reynolds numbers. It was observed that the effectiveness of superimposed vibrational flow is limited, especially for low Reynolds numbers and flow behavior index characteristic of shear-thinning fluids.

Corrigendum: About the Quantum‐Kinetic Derivation of Boundary Conditions for Quasiparticle Boltzmann Equations at Interfaces

Corrigendum: About the Quantum‐Kinetic Derivation of Boundary Conditions for Quasiparticle Boltzmann Equations at Interfaces

Optimization of thermal performance in hybrid nanofluids for parabolic trough solar collectors using Adams–Bashforth–Moulton method

Solar energy is an environment-friendly renewable energy source that fulfills the energy consumption requirements of the world economy. Concentrated solar power is the most advanced technology that efficiently absorbs concentrated solar energy. This article investigates the solar energy storage and entropy generation in concentrated parabolic trough solar collectors for hybrid nanofluid flow due to a spinning tube inserted at the receiver's center. The effect of copper nanoparticles and multi-walled carbon nanotube mixture is studied using water-based fluid. Also, the influence of an external magnetic field is investigated with thermophoresis and Brownian motion of nanoparticles. The governing equations for the nanofluid flow are derived using the conservation principles of mass, momentum, energy, and concentration with boundary layer assumptions and no-slip boundary conditions. The governing equations are numerically solved using Adam Bashforth and Adam Moulten's fourth-order predictor-corrector numerical scheme along with the shooting approach. The numerical outcomes for the fluid velocity, temperature, Nusselt number, and entropy formation against various influential physical parameters have been discussed using graphs. It is observed that the augmenting Eckert number, thermophoretic diffusion parameter, and Brownian motion parameter escalates the thermal profiles. Also, an escalation of the radiation parameter and Lewis number enhances the Nusselt number. Thermal energy storage in solar collectors utilizing different terminologies is crucial for improving the efficiency and performance of solar thermal collectors.

Biofilm mass transfer and thermodynamic constraints shape biofilm in trickle bed reactor syngas biomethanation

Syngas (H2, CO2, CO) produced thermochemically from lignocellulosic biomass is an underexploited source for resource recovery and valorisation through its biological conversion for the production of a wide range of chemicals and fuels. Syngas biomethanation is one such promising bioconversion pathway, displaying interesting features such as high conversion efficiency and product selectivity, as methane is the sole product of the process. The biological conversion of syngas to high purity biomethane is still typically associated with a number of challenges related to the syngas composition, CO toxicity, and gas–liquid mass transfer limitations. In this work, the syngas biomethanation process carried out in trickle bed reactors was investigated at its boundary conditions to explore the limits of the process, focusing on syngas composition and mass-transfer conditions potentially limiting process performance. The process was found to be robust when exposed to CO excess, but highly sensitive to H2 excess, which caused severe inhibition even under small amounts of excess H2. This implies that biomethane purity comparable to natural gas can be achieved by addition of renewable H2, but this requires precise control to avoid process failure. Modulating the liquid recirculation rate and gas residence time allowed for a maximum methane productivity of 9.8 ± 0.5 mmol CH4 h−1 Lreactor−1 with full conversion of H2 and CO at a gas residence time of 1 h. Nevertheless, increasing the gas–liquid mass transfer with increasing liquid recirculation rate did not lead to increased methane productivity, which suggested additional rate-limiting bottlenecks in the process. Careful investigation of other factors potentially limiting the process led to the conclusion that diffusive transport of syngas components in the biofilm was the main bottleneck of the process. This diffusive limitation leads to a scenario of severe substrate scarcity in the biofilm phase that conditions the thermodynamic feasibility of the different biochemical reactions involved in syngas biomethanation. In turn, these thermodynamic constraints were found to drive the stratification of microbial groups and sequential consumption of syngas components along the height of the reactor

Simulation and optimisation of vacuum (pressure) swing adsorption with simultaneous consideration of real vacuum pump data and bed fluidisation

Pressure swing adsorption (PSA) is a promising technology for gas separation and purification, receiving considerable attention in the past decades. Existing studies on simulating PSA involving a vacuum pump, the vacuum (pressure) swing adsorption [V(P)SA] process, typically estimate vacuum pump energy consumption using an approximate constant energy efficiency or an empirical energy efficiency correlation, leading to inaccurate representation of realistic vacuum pump performance. In this work, we propose an enhanced computational approach for simulation and optimisation of V(P)SA processes through simultaneous incorporation of realistic vacuum pump performance prediction models and adsorption bed fluidisation limits. The vacuum pump prediction models are developed using data-driven modelling techniques with realistic vacuum pump performance data. The enhanced computational approach more accurately accounts for the V(P)SA performance, without relying on an estimated vacuum pump efficiency and assumed pressure/flow rate boundary conditions at the vacuum pump end of the adsorption bed. Furthermore, this approach prevents the bed fluidisation. The computational results show that the developed prediction models accurately represent the actual performance curves of the vacuum pump. It is also demonstrated that incorporating the vacuum pump prediction models and fluidisation constraints in PSA optimisation leads to significantly different optimal solutions compared to when these factors are not considered.

Effect of Source-Sink Misalignment and Sink Cavity Aspect Ratio on the Thermal Performance of Gallium Heat Sinks

Herein, we numerically investigated methods to improve heat dissipation from the base of a phase-change material (PCM) heat sink exposed to high heat flux (15 W/cm2). The sink cavity exhibited a rectangular cross-section with two opposite vertical walls of high thermal conductivity. These walls serve as heat-dissipating fins. We used gallium, a metallic PCM with high thermal conductivity and a low melting point to improve heat dissipation by leveraging buoyancy-driven circulations within the molten gallium. We adopted a method to position the sink base on top of the heat source surface in such a way to augment temperature gradients within the liquid to induce imbalances in gallium density and consequent internal circulations. In this technique, one vertical wall is subjected to forced-convective heat exchange with the surroundings, while the other remains adiabatic. Subsequently, for the case which showed optimum performance, we relaxed the adiabatic thermal boundary condition applied at one of the vertical walls and replaced it with convective heat transfer condition, and assessed the sink performance based on that. Further, we analyzed the impact of various aspect ratios of the sink cavity on the induced internal circulations and base cooling. Sink cavities with aspect ratios of 1.33, 0.96, and 0.64 corresponding to sink heights of 20, 17, and 14 mm, respectively, were examined under a consistent heat flux applied across a 15-mm base length. We also investigated the precise positioning of the sink base relative to the heat source surface and explored its influence on temperature gradients and gallium density imbalances crucial for induced internal circulations. Results indicate that a sink cavity aspect ratio of 0.96 yields the optimal heat dissipation from the base. Further, aligning the left bottom edge of the sink base with the left edge of the heat source surface maximizes heat dissipation to the ambient surroundings, thereby prolonging latent-heat dumping into the sink before the gallium completely melts. Any misalignment between the sink base and the surface of the heat source may have a profound impact on the temporal evolution of induced internal circulations within the molten gallium, essentially affecting heat dissipation from the sink base.

Coherent flow structures and magnetic field patterns in rotating spherical shell convective dynamos: A data-driven approach

The Earth's outer core dynamics involve convective fluid motion generating an observable geomagnetic field. The velocity and magnetic fields exhibit characteristic spatiotemporal features possessing geophysical significance for which extensive datasets are available from direct observations and computational simulations. This study demonstrates the robustness of proper orthogonal decomposition (POD), a data-driven technique, in detecting prominent and relevant features in these datasets. Improvising on previous practices, the POD efficiently detects infinitesimal instabilities at the onset of convection, providing an accurate and objective methodology to determine the convective threshold, even for heterogeneous buoyancy forcing. Time evolution of paired, phase-shifted modes efficiently reconstructs the azimuthally drifting of traveling wave instabilities. Simultaneously reduced order modeling of velocity components clearly distinguish the equatorial and polar coherent flow structures. Supercritical convection-driven magnetic field data over long periods, generated using numerical simulations, produce dominant modes that are more accurately representative of time-averaged patterns than geocentric axial dipole patterns. Moreover, the quantitative significance of the dominant modes determines the extent of dimensional reduction complementing established diagnostics for dipolarity. Finally, analysis of observational geomagnetic field data reveals long-lived dominant patterns influenced by thermal core–mantle interaction consistent with numerical models employing tomographic heat flux boundary conditions in present as well as previous studies.

A study of sign-changing Poisson-type equations in two configurations

Abstract. In this paper, We are interested in specific non-coercive issues concerning electromagnetic wave propagation in the presence of metals or particular metamaterials. We focus on some non-coercive problems which cannot be studied by a classical Lax-Milgram approach. We consider the sign-changing Poisson-type equations with the homogeneous Dirichlet boundary conditions in the circular configuration and in the three square configuration. We utilize a method called T-isomorphism, which allows the conversion of non-coercive problems into coercive ones, to examine specific Poisson-type equations with sign-changing parameters. We first use the classical method to transform the equations into the corresponding variational formulations, and then apply the T-isomorphism method to show the well-posedness of these variational formulations with some parameters. By applying appropriate isomorphisms, we can conclude the well-posedness of the related variational formulations with certain parameters, which constitute the main findings of this paper. After giving appropriate isomorphisms, we can deduce the well-posedness of the corresponding variational formulations with some parameters, which are our main results in this paper.

Volumetric Drag Coefficients for Generic Urban Configurations: Insights from Canopy Flow Analysis

Alternative drag approaches for representing unresolved buildings were proposed in literature for computational fluid dynamics (CFD) simulation of macroscopic urban airflow. As a contribution, the present work derives the volumetric drag coefficient (Cd*) through canopy drag and velocity analysis and provides appropriate correlations for Cd*against urban morphological parameters. A total of 72 cases across various urban configurations are investigated, categorized by building typology, horizontal layout, height variability, and plan area density (λp, from 0.0625 to 0.57). Reynolds-Averaged Navier-Stokes (RANS) simulations with periodic boundary conditions are performed to model fully developed flows. Results for the normalized drag force and superficial velocity and their relations with λp are evaluated. Subsequent evaluation of the profiles for the sectional coefficients (Cd*(Z)) reveals four distinct types with variations in uniform-height cases and combinations in varying-height cases. A throughout correlations analysis, facilitated by data transformation, identifies the straightforward relations between Cd* and frontal area density (λf) and tortuosity (τ). The followed stepwise regression provides a recommended formula for Cd*, demonstrating a proper fit with the simulated values. These findings facilitate the understanding and appropriate estimation of Cd* and Cd*(Z), promoting the application of macroscopic turbulence models, for neighborhood-scale wind and air quality studies.

A reverse Faber-Krahn inequality for the magnetic Laplacian

We consider the first eigenvalue of the magnetic Laplacian in a bounded and simply connected planar domain, with uniform magnetic field and Neumann boundary conditions. We investigate the reverse Faber-Krahn inequality conjectured by S. Fournais and B. Helffer, stating that this eigenvalue is maximized by the disk for a given area. Using the method of level lines, we prove the conjecture for small enough values of the magnetic field (those for which the corresponding eigenfunction in the disk is radial).

The Role of Thermal Analysis in Engine Fin Design: Insights and Perspectives

Objective: The study aims to address the significance of fins in enhancing the heat transfer rate of engine cylinders and highlight their applications and advancements in various fields. Methods: Thermal analysis of fins using ANSYS Workbench is an invaluable tool for engineers and researchers in determining and optimizing the heat transfer characteristics of finned structures. ANSYS Workbench provides a comprehensive platform for evaluating the performance of fins in various applications using simulations of computational fluid dynamics (CFD) and finite element analysis (FEA). With ANSYS Workbench, it is possible to assess the thermal behaviour of fins under a variety of conditions, such as variable boundary conditions and thermal loads. Some significant recent patents on Engine fin are also discussed in this review article. Results: Fins are designed to maximize contact with the surrounding air or coolant, thereby facilitating engine heat transfer to the surrounding environment. Engineers can enhance heat dissipation capabilities by considering parameters such as fin geometry, density, material selection, and manufacturing techniques while minimizing associated drawbacks. Conclusion: The numerous applications of fins highlighted in this review article will encourage researchers to find novel strategies to improve heat transfer by employing fins for use in a variety of sectors.

Shape reconstruction of a cavity with impedance boundary condition via the reciprocity gap method

We consider an interior inverse scattering problem of reconstructing the shape of a cavity with impedance boundary condition from measured Cauchy data of the total field. The incident point sources and the measurements are distributed on two different manifolds inside the cavity. We first prove that the boundary of the cavity and the surface impedance can be uniquely determined by the scattered field data on the measurement manifold. Then we develop a reciprocity gap (RG) method to reconstruct the cavity. The theoretical analysis shows the uniquely solvability and existence of the approximate solution for the linear integral equation constructed in the RG method. We also prove that the shape of the cavity can be characterized by the blow-up property of the approximate solution of the proposed integral equation. Numerical examples are presented to verify the feasibility of the RG method.

Direct hydroacoustics analyses of pipe orifice using lattice Boltzmann method

Acoustic waves in water pipes pose a structural threat but also offer a valuable tool for non-invasive inspection. Complex pipe geometries such as bends and surface liners create complex fluid boundaries, triggering interactions between fluid flow and acoustics. The lattice Boltzmann method (LBM) excels at capturing these interactions near complex boundaries, in contrast to the finite volume method, which can result in errors. However, the application of LBM in water pipes has been limited by stability problems. This study proposes a two-step (DM-TS) collision operator based on a direct method (LBM-HA) for stable LBM simulations in water pipes. LBM-HA enables direct hydroacoustic predictions for both acoustic and dynamic flow fields in a pipe orifice. A novel acoustic-dynamic complex boundary condition was formulated from the linearized Euler's equation to account for the physical effects of the bounded domain of the LBM-HA analysis. To distinguish between dynamic and acoustic components, wavenumber and frequency decomposition techniques were applied to the pressure data obtained within the pipe. In addition, mode decomposition of the acoustic pressure was conducted to identify the dominant acoustic modes. By comparing the LBM-HA results with experimental data, we highlighted the method's potential for direct hydroacoustic analysis considering the acoustic characteristics of the water pipe.

On the uniqueness of the solution and finite-dimensional attractors for the 3D Flow with dynamic slip boundary condition

On the uniqueness of the solution and finite-dimensional attractors for the 3D Flow with dynamic slip boundary condition

Forced flow reversal in ferrofluidic Couette flow via alternating magnetic field

Time-dependent boundary conditions are very common in natural and industrial flows and by far no exception. An example of this is the movement of a magnetic fluid forced due to temporal modulations. In this study, we used numerical methods to examine the dynamics of ferrofluidic wavy vortex flows (WVF2, with dominant azimuthal wavenumber m=2) in the counter-rotating Taylor–Couette system, which was subjected to time-periodic modulation/forcing in a spatially homogeneous magnetic field. In the absence of a magnetic field, all WVF2 states move in the opposite direction to the rotation of the inner cylinder, they are retrograde. However, when strength or frequency of the alternating magnetic field increases, the motion direction of the flow pattern changes. Thus, the alternating field provides a precise and controllable key parameter for triggering the system response and controlling the flow. Aside, we also observed intermittent behavior when one solution became unstable, leading to random transitions in both, the transition time and toward the different final solutions. Our findings suggest that, in ferrofluids, flow pattern reversal can be induced by varying a magnetic field in a controlled manner, which may have applications in the development of modern fluid devices in laboratory experiments. These findings provide a framework to study other types of magnetic flows driven by time-dependent forcing.

Electro-osmotic flow and heat transfer in Jeffery fluid: A multi-membrane microchannel model

Microscale heat transfer is vital for the performance of smart thermal devices like heat sinks, thermosyphons, and microheat pipes. This study introduces a biothermal pumping flow model based on a multi-membrane pumping mechanism that leverages microscale heat transfer. The model describes rhythmic contraction and relaxation of membranes, combined with electro-osmosis in Jeffery fluid flow within a vertical microchannel of finite length. Two membranes on the microchannel walls, with varying amplitudes, diameters, and phase lags, generate pressure that moves fluid in both directions through contraction and expansion cycles. The model is based on the conservation of mass and momentum, using a low Reynolds number approximation to capture microscale transport phenomena at biomedical scales. Dimensionless conservation equations are analytically solved under no-slip boundary conditions, with results computed in MATLAB for clarity. Axial velocity results are simulated and verified using the optimal homotopy analysis method. The model explores the influence of key parameters (UHS, me, λ, Gr, β) on pressure gradient, velocity distribution, volumetric flow rates, skin friction, Nusselt number, and stream function. The findings demonstrate that pressure from membrane motion is significantly affected by thermal effects and buoyancy forces, and flow and pumping characteristics are largely determined by the fluid's rheological qualities and the geometrical features of the membrane. This study provides novel ideas for enhancing the functionality and design of smart thermal devices while also advancing microscale heat transfer technology.

Characteristics of constant pressure dehydrogenation and construction of linear diffusion model for solid hydrogen storage materials: A case study of TiMn1.5 alloy

Solid-state hydrogen storage is considered to be the most promising hydrogen storage technology because of its high volumetric capacity and good security. The process of hydrogenation and dehydrogenation of solid hydrogen storage materials can be characterized by a diffusion kinetic model, which can be used to guide the safe and efficient application of hydrogen storage materials. Currently, experimental study on the dehydrogenation characteristics under constant pressure outer boundary condition has not been carried out methodically, and a simple and reliable diffusion kinetic model under this condition should be established based on the experimental results. In this study, an experimental platform for hydrogen absorption and constant pressure dehydrogenation of solid hydrogen storage materials was independently established. The TiMn1.5 alloy and its interstitial hydrides were taken as the experimental research objects to conduct constant pressure dehydrogenation tests, and the dehydrogenation properties at different temperatures were systematically studied. The results show that the dehydrogenation rate decreases with the increase of testing time. The ultimate capacity of hydrogen release increases first and then decreases with the increasing temperature, reaching the maximum value at 30 °C. It is found that the diffusion coefficient increases linearly with time during constant pressure dehydrogenation process. On this basis, a linear diffusion kinetic model is established. The initial diffusion coefficient D0 calculated by the new model has a quadratic polynomial relationship with temperature, and the increase coefficient α of diffusion coefficient has a linear relationship with temperature. According to the mathematical expressions, the dehydrogenation process of TiMn1.5 alloy at different temperatures can be predicted. The research results can provide theoretical guidance for the safe, efficient and controllable application of hydrogen storage materials.

Dynamic modeling of air–metal plasma mixture of single-turn coil with erosion at megaGauss magnetic field

Single-turn coil (STC) is a destructive pulsed magnet aiming at 100–300 T magnetic field. Due to the high discharge current, the conductor of STC is heated rapidly and undergoes melting and vaporization, leading to the generation of supersonic air–metal vapor mixed plasma jet and the magneto-fluid effect. In this study, the mixed plasma mass-transfer and fluid dynamic characteristics are modeled at megaGauss magnetic field, high temperature, high pressure, and supersonic conductor shock deformation. The collision integral method is employed to calculate the fluid transport properties. In addition, a boundary constraint model of fluid–structure interaction (FSI) compatible with both fluid wall boundary condition and plasma jet entrance condition and a model to simultaneously solve the thermal ionization and high electric field ionization of the mixed vapor are proposed. As the result, the distributions of plasma electrical conductivity, current density, electron, heavy particles, temperature, air body load, and velocity are derived. Especially, the region of highest electrical conductivity is not the air domain near the inner surface of the conductor with the highest electron density and the highest magnetic field, but the air domain near the outer surface of the conductor with the relatively higher electron density and lower magnetic field.