Interface Boundary Conditions Research Articles

Investigating surface and internal waves and viscosity effects in two-layer fluids with a streamlined moving body using fully nonlinear simulation and experiment

Waves generated by a moving submerged moving body were analyzed in a two-layer density fluid. The experiments were performed in a two-dimensional towing tank and compared to the numerical results from a two-dimensional, fully nonlinear numerical towing tank (NTT) based on potential flow. The experimental conditions were determined by the relationship between the depth-dependent Froude number and the critical Froude number in baroclinic mode. The study was performed with various body velocities, changing the upper and lower fluid depth while maintaining the total water depth. An artificial damping coefficient was applied to the entire free surface and interface boundary condition to simulate the viscosity effect of the fluid interface. When the baroclinic mode was dominant, internal precursor soliton waves and depression waves were similar to the experimental results with strong nonlinearity. In contrast, the trailing waves were significantly different from the experimental results without applying an artificial damping scheme. A comparison with the experimental data showed that the effects of the fluid viscosity can be simulated with an appropriate damping coefficient. As the body velocity increases, the depressed or rising waves associated with the baroclinic mode are mixed with the trailing waves of the barotropic mode. As the depth of the upper fluid decreases, the characteristics of the baroclinic mode decrease, and the barotropic mode becomes dominant under the same velocity conditions.

Convergence analysis and applicability of a domain decomposition method with nonlocal interface boundary conditions

In the past, the domain decomposition method was developed successfully for solving large-scale linear systems. However, the problems with significant nonlocal effect remain a major challenger for applying the method efficiently. In order to sort out the problem, a non-overlapping domain decomposition method with nonlocal interface boundary conditions was recently proposed and studied both theoretically and numerically. This paper is the report on the further development of the method, aiming to provide a comprehensive convergence analysis of the method, with supplementary numerical tests to support the theoretical result. The nonlocal effect of the problem is found to be reflected in both the governing equation and boundary conditions, and the effect of the latter was never taken into account, although playing a significant role in affecting the convergence. In addition, the paper extends the applicability of the analysis result drawn from the Poisson’s equation to more complicated problems by examining the symbols of the Steklov–Poincaré operators. The extended application includes a model equation arising from fluid dynamics and the high performance of the domain decomposition method in solving this equation is better elaborated.

Acoustic Radiation Characteristics of a Forced Vibrating Elastic Panel Under Thermal Environments

This paper presents the thermo-acoustic frequency response of an un-baffled rectangular panel subjected to an external excitation load. A boundary element method BEM has been employed taking into account the Kirchhoff-Helmholtz K-H integral equation for the acoustic pressure and with the fluid-plate interface boundary condition the acoustic pressure jump over the panel is calculated. The thermal effects are considered regarding in the form of a uniform increment of temperature of the panel and are analysed in order to prevent the buckling phenomena. The excitation force considered is in the form of a concentrated load at some point of the panel and the deformation modes correspond to the vacuum case. Applying a collocation method for the panel equation, a frequency transfer function is obtained that relates the deflexion of the panel with the applied load. The effect of several geometric parameters, different thermal loads and location of the load applied on the acoustic radiated power and the acoustic efficiency spectrum are evaluated. Furthermore, the influence of the excitation frequency on the sound directivity is evaluated. The verification of the method is proven with other works.

Heat transfer in shear flow over a still nanofluid: A direct integration approach

Heat transfer in a two-layer fluid system, including shear-driven viscous fluid in the upper layer and still nanofluid in the lower layer, is investigated in this work. Following single phase model for nanofluid, three different types of nanofluids are considered. Different from existing literatures on two-layer fluid flow problem, this work introduces the heat transfer of shear-driven viscous fluid over a still nanofluid which makes this investigation new. The main equations which are partial in nature are reduced to a set of coupled nonlinear ordinary differential equation (ODE) by employing suitable similarity transformations. These ODEs together with the coupled interfacial boundary conditions are supplied to MATLAB software to get the solution numerically. Comparing the results obtained in this study with the previous authoritative research results, the accuracy and correctness of the present research methodology are ensured. The results show that the temperature of the upper fluid and the lower fluid are controlled by the thermophysical properties of different types of nanofluids. Also, opposite behavior of temperature distribution for the upper and the lower fluid are observed.

Magnetohydrodynamic mixed convection of nanofluid within an annular space filled partially with a porous layer between a lid-driven wall and a discrete heater

Heat and mass transfer inside a partially porous medium is a complex coupling process in which the interface boundary conditions between free and porous space are critical for a wide range of engineering applications, including natural and engineered fractures in oil and gas extraction. In this study, we looked into the numerical simulation of magnetohydrodynamic mixed convection of nanofluid in an annular partially porous space between two coaxial cylinders with a permeable interface, in which we investigated the contribution of the interface in order to better understand the heat transfer processes. The inner cylinder is subjected to a discrete uniform heat flux, while the external cylinder is kept at a uniform cold temperature and moving at a constant speed. The base walls are designed to be impermeable and insulated. A finite difference-based vorticity-stream function method is used to solve the nonlinear coupled conservation equations using the Successive Over Relaxation approach (SOR). The acquired numerical results in terms of streamlines, isotherms, and Nusselt numbers are shown to demonstrate the effect of various control factors such as the Rayleigh number, Darcy number, Hartmann number, Reynold number, nanoparticle concentration, and the direction of the outer wall. The results of this numerical simulation show that the upward or downward direction of the wall under these control factors plays an important role in improving the rate of heat transfer or replacing the magnetic effect by limiting and controlling heat transfer

Revealing molecular insights into surface charge and local viscosity in electroosmotic flows

The limitations of the continuum theory in predicting osmotic response at the nanoscale stem from its lack of molecular-level insight into local fluid properties and the interfacial structure of fluid and electrolyte solutions. To overcome this challenge, our study integrates molecular dynamics (MD) simulation with the continuum framework to explore how surface charge and various hydrodynamic properties impact electroosmotic flow (EOF). The failure of continuum theories to account for molecular interactions and geometric boundaries leads to significant disparities between MD simulations and continuum predictions, influenced by local fluid properties and the electric field. Emphasizing the importance of incorporating appropriate local hydrodynamic properties and atomic interface boundary conditions, our findings bridge the gap between MD simulations and continuum EOF predictions. Our computational results and theoretical model, considering surface charge, atomic interface boundaries, and dynamic structure-based hydrodynamic properties, provide crucial insights and guidance for EOF investigations.

A network model to predict ionic transport in porous materials

Understanding the dynamics of electric-double-layer (EDL) charging in porous media is essential for advancements in next-generation energy storage devices. Due to the high computational demands of direct numerical simulations and a lack of interfacial boundary conditions for reduced-order models, the current understanding of EDL charging is limited to simple geometries. Here, we present a network model to predict EDL charging in arbitrary networks of long pores in the Debye-Hückel limit without restrictions on EDL thickness and pore radii. We demonstrate that electrolyte transport is described by Kirchhoff's laws in terms of the electrochemical potential of charge (the valence-weighted average of the ion electrochemical potentials) instead of the electric potential. By employing the equivalent circuit representation suggested by these modified Kirchhoff's laws, our methodology accurately captures the spatial and temporal dependencies of charge density and electric potential, matching results obtained from computationally intensive direct numerical simulations. Our network model provides results up to six orders of magnitude faster, enabling the efficient simulation of a triangular lattice of five thousand pores in 6 min. We employ the framework to study the impact of pore connectivity and polydispersity on electrode charging dynamics for pore networks and discuss how these factors affect the time scale, energy density, and power density of capacitive charging. The scalability and versatility of our methodology make it a rational tool for designing 3D-printed electrodes and for interpreting geometric effects on electrode impedance spectroscopy measurements.

Numerical Reconstruction of Time-Dependent Boundary Conditions to 2D Heat Equation on Disjoint Rectangles at Integral Observations

In this paper, two-dimensional (2D) heat equations on disjoint rectangles are considered. The solutions are connected by interface Robin’s-type internal conditions. The problem has external Dirichlet boundary conditions that, in the forward (direct) formulation, are given functions. In the inverse problem formulation, the Dirichlet conditions are unknown functions, and the aim is to be reconstructed upon integral observations. Well-posedness both for direct and inverse problems is established. Using the given 2D integrals of the unknown solution on each of the domains and the specific interface boundary conditions, we reduce the 2D inverse problem to a forward heat 1D one. The resulting 1D problem is solved using the explicit Saul’yev finite difference method. Numerical test examples are discussed to illustrate the efficiency of the approach.

A New Solution of Drag for Newtonian Fluid Droplets in a Power-Law Fluid

Understanding flow behaviors of multiple droplets in complex non-Newtonian fluids is crucial in many science and engineering applications. In this study, a new and improved analytical solution is developed based on the free surface cell model for the flow drag of swamp of Newtonian fluid drops through a power-law fluid. The developed solution is accurate and compares well to the numerical solutions. The improvement involves a new quantification of shear stress boundary condition at the interface and a more consistent approximation in linearizing the power-law fluid flow governing equation. The Newtonian fluid solutions can be reasonably used to linearize the flow governing equation. The approximation of the boundary conditions at the interface, however, has a major impact on the model prediction. The main improvement in the new solution is observed under the condition of comparable viscosities of the Newtonian drops and the outside power-law fluid when the results are sensitive to the interface boundary condition. Under the two extreme conditions of high viscosity ratio (approaching particles) and low ratio (approaching bubbles), the present and existing solutions converge.

THE ANALYSIS OF THE FLOW OF BLOOD IN A STENOSED ARTERY THROUGH SIMULATION: A COMPARISON AMONG VARIOUS NON-NEWTONIAN MODELS

This paper focuses on the dual quality of blood, Newtonian and non-Newtonian, in particular by exploring the energy curves. Careful investigation of the dual property of blood has been made by considering two different geometries to represent a stenosed arterial segment. We present a cautious assessment of non-Newtonian blood rheology impacts in arterial stream simulations by coupling the Newtonian and non-Newtonian models. The flow of energy through the two flow dimensions is meticulously investigated using velocity (kinetic energy), pressure, and wall shear stress (pressure energy). Besides, the proper implementation of an interface boundary condition (IBC) was emphasized to ensure consistency with the flow conditions downstream of a backward-facing step. The integration of the Newtonian and non-Newtonian models adjoins the novelty of the current research. The energy curves are obtained by implementing five different non-Newtonian models to designate a suitable non-Newtonian model for blood flow investigations. The combination of the non-Newtonian models enforced in this research is novel and particular attention is paid to the energy curves obtained. The conclusion was to elect the Carreau model as a suitable non-Newtonian rheological model for the blood flow study. This study was able to finalize the fact that the coupling of Newtonian and non-Newtonian models is necessary to obtain accurate results. For the sinusoidal waveform considered for the velocity, Carreau and the Power law models yield better results, eliminating the other non-Newtonian models from the list. With a better inlet condition imposed in the form of the Fourier series for pressure and velocity, the Carreau model yields the best results.

The effect of shear rate and overconsolidation ratios on undrained clay–suction caisson interface behaviors

The interface resistance during installation is crucial for the stability and safety of suction caisson in offshore geotechnical engineering, which is strongly affected by the penetration rate and soil–structure interface mechanical properties. This research conducts a series of clay–structure interface shear tests using modified direct simple shear device to fully study the mechanical behavior of clay–suction caisson interface. The effect of shear rate, over consolidation ratios (OCRs), interface boundary conditions, stress levels, and interface roughness were considered. Results show that as the OCR increases, the strength of both the clay and interface increase but show distinct patterns under constant volume (CV) and constant normal load (CNL) boundary condition. It was found that the interface strength is positively related to interface roughness and shear rate impact both the clay and corresponding interface strength. Under CNL conditions, the strength of normally consolidated (NC) clay decreases with rising shear rate, while the over consolidated (OC) clay demonstrate a opposite trend. In contrast, the effect of shear rate on interface behavior gets complicated owing to the combination of roughness, stress levels, and OCRs. Under CV conditions, the shear strength of clay and interface exhibits a logarithmic growth relationship with shear rates. The result of this work can provide a basis for interface resistance evaluation for suction caisson installation in clay.

Air-water interface boundary condition for the numerical evaporation rate prediction of a horizontal water span under different convection regimes

This paper conducts a comprehensive analysis of the existing numerical models addressing the dynamics of water vapor flux across the air-water interface. Additionally, it introduces a novel model based on the empirical friction velocity of the air over a water surface. This new model is used to predict evaporation rates, especially in the context of wind tunnel experiments where water tanks are subjected to controlled drying conditions. Then, these predictions are compared to both the empirically measured evaporation rates and those generated by other numerical models. A wide range of configurations are simulated to cover the different convection regimes. Overall, these simulations show good agreement between the model's predictions and the observed experimental measurements of evaporation rates, underscoring its robustness and reliability. One notable advantage of this model is its remarkable independence from the convection regime, setting it apart from traditional empirical correlations that typically exhibit such dependencies.

Photosynthetically-powered phototactic active nematic liquid crystal fluids and gels

One of the most ancient forms of life dating to ~3.5 billion years ago, cyanobacteria are highly abundant organisms that convert light into energy and motion, often within conjoined filaments and larger colonies that attract a great deal of interest but their active nematic behavior remained unexplored. Here we demonstrate how light causes a spontaneous self-assembly of two- and three-dimensional active nematic states of cyanobacterial filaments, with a plethora of topological defects. We quantify light-controlled evolutions of orientational and velocity order parameters during the transition between disordered and orientationally ordered states of photosynthetic active matter, as well as the subsequent active nematic’s fluid-gel transformation. Patterned illumination and foreign inclusions with different shapes interact with cyanobacterial active nematics in nontrivial ways while inducing interfacial boundary conditions and fractional boojum defects. Our phototactic model system promises opportunities to systematically explore fundamental properties and technological utility of the liquid crystalline active matter.

Multi-species kinetic-fluid coupling for high-energy density simulations

Many physics problems are subject to a mix of continuum and non-equilibrium flows or conditions where one transition into the other, as a function of time and/or space. Numerically, such flows could be described with a purely kinetic method which, in the limit of small particle mean free paths, reproduces a continuum regime. However, from a computational perspective, such approaches are usually very expensive or simply not feasible, even on modern computational architectures. A solution is the construction of hybrid methods which connect the kinetic and hydrodynamic system of evolution equations. Such coupling usually requires specific boundary conditions between non-equilibrium and continuum regimes which come with their own numerical challenges and can introduce unphysical effects. Here, we present a hybrid model which uses a buffer region to transition from a kinetic into a fluid description following the original work of Degond et al. [1]. In the buffer region, both, the kinetic and hydrodynamic equations are solved simultaneously while being coupled via a so-called transition function. The latter also ensures a smooth conversion from the coupled model to either the kinetic or continuum approach at the interfaces of the buffer region. With that, the method avoids the need to find direct interface boundary conditions and allows one to localize the use of a high dimensional kinetic model only where it is needed. In our work, we extend the original method of Degond to flows with multiple particle species in 3D. We derive the coupled equations for the multispecies Vlasov-Bhatnagar-Gross-Krook (VBGK) model and its limiting Euler or Navier-Stokes hydrodynamic equations. To validate our model numerically, we simulate a Sod shock problem and study the effect of kinetic multi-species mixing in the preheat phase of a high energy-density plasma physics experiment.

Influence of hydrophobic particles on the film drainage during bubble–solid interaction

In this article, we investigate the complex drainage behavior of liquid (“dimple”) films entrapped between hydrophilic glass substrates and air bubbles in water and aqueous suspensions of polystyrene (PS) particles. The film drainage was monitored by capturing the evolution of spatial-temporal thicknesses using interferometry. Faster drainage of the entrapped film is observed in the PS colloid suspensions compared to the water. The film drainage strongly depends on the interface boundary conditions. Our experiments reveal that the air–liquid interface of the entrapped film becomes partially mobile in the PS suspensions, which results in faster drainage. The hydrophobic PS particles tend to migrate toward the air–water interface and form an ordered layer. We argue that the adsorbed layer of PS particles makes the bubble surface elastic, delaying the formation of the immobile interface.

Application of rigorous interface boundary conditions in mesoscale plasticity simulations

The interactions between dislocations and interface/grain boundaries, including dislocation absorption, transmission, and reflection, have garnered significant attention from the research community for their impact on the mechanical properties of materials. However, the traditional approaches used to simulate grain boundaries lack physical fidelity and are often incompatible across different simulation methods. We review a new mesoscale interface boundary condition based on Burgers vector conservation and kinetic dislocation reaction processes. The main focus of the paper is to demonstrate how to unify this boundary condition with different plasticity simulation approaches such as the crystal plasticity finite element (CPFEM), continuum dislocation dynamics (CDD), and discrete dislocation dynamics (DDD) methods. In DDD and CDD, plasticity is simulated based on dislocation activity; in the former, dislocations are described as discrete lines while in the latter in terms of dislocation density. CPFEM simulates plasticity in terms of slip on each slip system, without explicit treatment of dislocations; it is suitable for larger scale simulations. To validate our interface boundary condition, we implemented simulations using both the CPFEM method and a two-dimensional CDD model. Our results show that our compact and physically realistic interface boundary condition can be easily integrated into multiscale simulation methods and yield novel results consistent with experimental observations.

Small Darcy number limit of the Navier–Stokes–Darcy system

We study the small Darcy number behavior of the Navier–Stokes–Darcy system with the conservation of mass, Beavers–Joseph–Saffman–Jones condition, and the Lions balance of the normal-force interface boundary conditions imposed on the interface separating the Navier–Stokes flow and Darcy flow. We show that the asymptotic behavior of the coupled system, at small Darcy number, can be captured by two semi-decoupled Darcy number independent sequences: a sequence of (linearized) Navier–Stokes equations, and a sequence of Darcy equations with appropriate initial and boundary data. Approximate solutions to any order of the small parameter (Darcy number) can be constructed via the two sequences. The local in time validity of the asymptotic expansion up to second order is presented. And the global in time convergence is derived under the assumption that the Reynolds number is below a threshold value.

An approach to the global well-posedness of a coupled 3-dimensional Navier-Stokes-Darcy model with Beavers-Joseph-Saffman-Jones interface boundary condition

<abstract><p>This study focused on investigating the global well-posedness of a coupled Navier-Stokes-Darcy model with the Beavers-Joseph-Saffman-Jones interface boundary condition in the three-dimensional Euclidean space. By utilizing this approach, we successfully obtained the global strong solution of the system in the three-dimensional space. Furthermore, we demonstrated the exponential stability of this strong solution. The significance of such coupled systems lies in their pivotal role in the analysis of subsurface flow problems, particularly in the context of karst aquifers.</p></abstract>

Global well-posedness of 2D incompressible Navier–Stokes–Darcy flow in a type of generalized time-dependent porosity media

<p>This study investigates the global well-posedness of a coupled Navier–Stokes–Darcy model incorporating the Beavers–Joseph–Saffman–Jones interface boundary condition in two-dimensional Euclidean space. We establish the existence of global strong solutions for the system in both linear and nonlinear cases where porosity depends on pressure. When dealing with the time-dependent porous media, the primary challenge in obtaining closed prior estimates arises from the presence of complex, sharp interfaces. To address this issue, we employ the classical Trace Theorem. Such space-time variable coupled systems are crucial for understanding underground fluid flow.</p>

Terminal-coupling induced critical eigenspectrum transition in closed non-Hermitian loops

A hallmark feature of non-Hermitian (NH) systems is the non-Hermitian skin effect (NHSE), in which the eigenenergy spectra of the system under open boundary conditions (OBC) and periodic boundary conditions (PBC) differ markedly from each other. In particular, the critical NHSE occurs in systems consisting of multiple non-Hermitian chains coupled in parallel where even an infinitesimally small inter-chain coupling can cause the thermodynamic-limit eigenenergy spectrum of the system to deviate significantly from the OBC spectra of the individual component chains. We overturn the conventional wisdom that multiple chains are required for such critical transitions by showing that such a critical effect can also be induced in a single finite-length non-Hermitian chain where its two ends are connected together by a weak terminal coupling to form a closed loop. An infinitesimally small terminal coupling can induce the thermodynamic-limit energy spectrum of the closed loop to switch from the OBC to the PBC spectrum of the chain. Similar to the critical NHSE, this switch occurs abruptly when the chain length exceeds a critical size limit. We explain analytically the underlying origin of the effect in a Hatano–Nelson chain system, and demonstrate its generality in more complex one-dimensional non-Hermitian chains. Our findings illustrate the generality of critical size-dependent effects in finite NH systems that arise from the interplay between the interfacial boundary conditions and the influence of edge localization.