Interface Conditions Research Articles

Study of the shear bond behavior of a hybrid fiber-reinforced concrete (HyFRC) composite slab

In this study, a new composite slab consisting of a prefabricated hybrid fiber-reinforced concrete (HyFRC) formwork and a cast-in-place ordinary concrete (OC) topping was proposed. The HyFRC was comprised of steel and PVA fibers. The compressive strength and bending strength of the HyFRC were 111.4 MPa and 15.6 MPa, respectively, and the toughness coefficient was 12.1 MPa. Seven full-scale composite slabs including six simply-supported slabs and one two-span continuous slab were tested under four-point bending to study the shear bond behavior. Three different contact surface conditions were considered for the simply supported slabs. Despite the different interface conditions, results showed that the interfacial bond strength between the HyFRC formwork and the OC topping was sufficient to achieve flexural failure. Numerical simulation was used to analyze the interfacial force transfer mechanism and was validated by the experimental data. The validated numerical model was used to study the influence of concrete topping depth and amount of reinforcement on the load-carrying capacity and longitudinal shear stress at the joint surface. The numerical results identified a critical condition of positive effect between reinforcement ratio and OC topping thickness. The peak load of a properly designed composite slab was almost two times that of the corresponding slab with ordinary concrete only, and the deformation capability was also larger. A procedure for calculating the flexural capacity of the composite slab was proposed.

Developing organs-on-chips for biomedical applications.

In recent years, organs-on-chips have been arousing great interest for their bionic and stable construction of crucial human organs in vitro. Compared with traditional animal models and two-dimensional cell models, organs-on-chips could not only overcome the limitations of species difference and poor predict ability but also be capable of reappearing the complex cell-cell interaction, tissue interface, biofluid and other physiological conditions of humans. Therefore, organs-on-chips have been regarded as promising and powerful tools in diverse fields such as biology, chemistry, medicine and so on. In this perspective, we present a review of organs-on-chips for biomedical applications. After introducing the key elements and manufacturing craft of organs-on-chips, we intend to review their cut-edging applications in biomedical fields, incorporating biological analysis, drug development, robotics and so on. Finally, the emphasis is focused on the perspectives of organs-on-chips.

Identification of excessive contact pressures under hand orthosis based on finite element analysis.

Implicit magnitudes and distribution of excessive contact pressures under hand orthoses hinder clinicians from precisely adjusting them to relieve the pressures. To address this, contact pressure under a hand orthosis were analysed using finite element method. This paper proposed a method to numerically predict the relatively high magnitudes and critical distribution of contact pressures under hand orthosis through finite element analysis, to identify excessive contact pressure locations. The finite element model was established consisting of the hand, orthosis and bones. The hand and bones were assumed to be homogeneous and elastic bodies, and the orthosis was considered as an isotropic and elastic shell. Two predictions were conducted by assigning either low (fat) or high (skin) material stiffness to the hand model to attain the range of pressure magnitudes. An experiment was conducted to measure contact pressures at the predicted pressure locations. Identical pressure distributions were obtained from both predictions with relatively high pressure values disseminated at 12 anatomical locations. The highest magnitude was found at the thumb metacarpophalangeal joint with the maximum pressure range from 13 to 78 KPa. The measured values were within the predicted range of pressure magnitudes. Moreover, 6 excessive contact pressure locations were identified. The proposed method was verified by the measurement results. It renders understanding of interface conditions underneath the orthosis to inform clinicians regarding orthosis design and adjustment. It could also guide the development of 3D printed or sensorised orthosis by indicating optimal locations for perforations or pressure sensors.

Threshold Voltage Modulation and Gate Leakage Suppression of Enhancement‐Mode GaN HEMT by Metal/Insulator/p‐GaN Gate Structure

Herein, a metal/insulator/p‐GaN gate HEMT (MIP‐HEMT) with Si3N4 gate dielectric layer is fabricated. Compared to the conventional p‐GaN HEMT, the MIP‐HEMT has a higher threshold voltage (Vth) of 4.8 V and better gate leakage suppression. The mechanism of threshold voltage increase in MIP‐HEMT is elucidated through an analysis of the electric field distribution in the gate region and the proposed gate capacitance model. Furthermore, MIP‐HEMTs with gate dielectric layers of varying dielectric constants and thicknesses are designed and simulated. The obtained results lead to the derivation of an empirical formula for the threshold voltage of the MIP‐HEMT under ideal dielectric/p‐GaN interface conditions, in relation to the dielectric constant and thickness of the gate dielectric layer. Additionally, simulations are conducted on the MIP‐HEMT incorporating fixed charges at the dielectric/p‐GaN interface, and the impact of interface charges is investigated. This refinement enables a more accurate prediction of the MIP‐HEMT's threshold voltage. Conclusively, MIP‐HEMTs can effectively modulate the threshold voltage by manipulating the parameters of the gate dielectric layer, thereby extending the threshold voltage range of conventional enhancement‐mode devices.

Streaming Potential and Associated Electrokinetic Effects through a Channel Filled with Electrolyte Solution Surrounded by a Layer of Immiscible and Dielectric Liquid.

The present article deals with the streaming potential-mediated pressure-driven flow across a channel in which the electrolyte solution is surrounded by a layer of cell membrane. Such a membrane of a biological cell may be modeled as an immiscible and dielectric liquid, which may bear free lipid molecules or charged surfactants. The presence of such additional charged molecules may lead to formation of liquid-liquid interfacial charge. In addition, the dielectric gradient-mediated ion partitioning effect further plays an important role in two-phase electrokinetic motion. We aim to study the generation of streaming potential and electrokinetic conversion efficiency as well as associated electroviscous effect for the undertaken problem. The mathematical model is based on the Poisson-Boltzmann equation for electrostatic potential and the Stokes equation for fluid flow, and the problem is studied considering suitable interfacial conditions for the flow variables along the liquid-liquid interface. The explicit analytical results for velocity and streaming field, electrokinetic energy conversion efficiency, and the parameter indicating the electroviscous effect are derived under the Donnan limit and within the Debye-Hückel electrostatic framework. We further numerically calculated the aforementioned intrinsic electrokinetic parameter associated with the problem undertaken for a wide range of pertinent parameters. The results are illustrated to indicate the impact of pertinent parameters on the generation of the streaming potential and associated electrokinetic effects.

Mesoscopic simulations of a fracture process in reinforced concrete beam in bending using a 2D coupled DEM/micro-CT approach

In this study, a complex fracture process in a short rectangular concrete (RC) beam reinforced by one longitudinal bar (without vertical reinforcement) and subjected to quasi-static three-point bending was numerically explored in 2D conditions. A critical diagonal shear crack in the beam caused it to fail during the experiment. The numerical simulations were conducted with a classical particle discrete element method (DEM). A three-phase concrete description (aggregate, mortar, and interfacial transitional zones (ITZs) around aggregates) accounted for the concrete heterogeneity. In mesoscopic DEM calculations based on a 2D X-ray CT scan, the actual shape and placement of aggregate particles in concrete were taken for granted. In the study, the steel bar with ribs was replicated. ITZ was also assumed between the bar and mortar. Without imposing any bond-slip law, a geometrical bar/concrete interface condition was explicitly considered. The focus was on the force–deflection diagram, fracture process, contact forces, and stresses along the bar. A good level of agreement about the evolution of the vertical force versus the deflection and failure mechanism was attained between DEM analyses and in-house laboratory tests despite simplified 2D conditions. A strong effect of concrete mesostructure on the crack pattern was found.

An Inverse Problem for Estimating Spatially and Temporally Dependent Surface Heat Flux with Thermography Techniques

In this study, an inverse conjugate heat transfer problem is examined to estimate temporally and spatially the dependent unknown surface heat flux using thermography techniques with a thermal camera in a three-dimensional domain. Thermography techniques encompass a broad set of methods and procedures used for capturing and analyzing thermal data, while thermal cameras are specific tools used within those techniques to capture thermal images. In the present study, the interface conditions of the plate and air domains are obtained using perfect thermal contact conditions, and therefore we define the problem studied as an inverse conjugate heat transfer problem. Achieving the simultaneous solution of the continuity, Navier–Stokes, and energy equations within the air domain, alongside the heat conduction equation in the plate domain, presents a more intricate challenge compared to conventional inverse heat conduction problems. The validity of our inverse solutions was verified through numerical simulations, considering various inlet air velocities and plate thicknesses. Notably, it was found that due to the singularity of the gradient of the cost function at the final time point, the estimated results near the final time must be discarded, and exact measurements consistently produce accurate boundary heat fluxes under thin-plate conditions, with air velocity exhibiting no significant impact on the estimates. Additionally, an analysis of measurement errors and their influence on the inverse solutions was conducted. The numerical results conclusively demonstrated that the maximum error when estimating heat flux consistently remained below 3% and higher measurement noise resulted in the accuracy of the heat flux estimation decreasing. This underscores the inherent challenges associated with inverse problems and highlights the importance of obtaining accurate measurement data in the problem domain.

Interfacial thermal resistance between mechanically exfoliated nm-thick MoS2 and silicon from −60 °C to 50 °C based on ns-ET Raman technique

In the miniaturization and integration development of micro-devices, the large thermal contact resistance significantly affects the overall performance of micro-devices, which is influenced by both interfacial contact condition and material properties. However, the temperature dependence of interfacial thermal resistance of realistic materials remains unclear. In this work, the interfacial thermal resistance between nm-thick MoS2 nanosheets and Si substrate is measured by using the ns energy transport state-resolved Raman (ns ET-Raman) technique. By comparing the experimental data with finite difference numerical simulation, the real value of the interfacial thermal resistance is determined. By changing the ambient temperature from -60 °C to 50 °C, the interfacial thermal resistance against temperature curve is obtained. As the environmental temperature deviates from room temperature, a large increasement of the interfacial thermal resistance can be observed. As temperature goes back to room temperature, the unreversible interfacial thermal resistance indicates that the large thermal expansion mismatch deteriorates the contact situation. In the last, the interfacial thermal resistance evolution under varying environmental temperatures from both experiments and modeling results in literatures are compared and discussed to illustrate the effect of different interface formation methods. For epitaxial and CVD grown materials, the interfacial thermal resistance generally reduces with raising ambient temperature. While for mechanically exfoliated samples, the effect of thermal expansion mismatch and the resulting gaps dominate the interfacial thermal resistance, which leads to non-monotonous temperature dependency with a valley located at room temperature.

Fluorescence-Based Calibration Model for In-Situ Measurement of Micro-scaled Lubricant Thickness Distribution at Indentation Interface

This study proposes a model for the measurement of microscale liquid film thickness distribution using fluorescence signals. The interfacial conditions between the tool and the workpiece in mechanical machining are important for understanding these phenomena and mechanisms. In this study, indentation tests with transparent tools were used to observe interfaces; however, it was challenging to obtain the signal from a thin fluorescent liquid film on smooth and steeply inclined surfaces. Therefore, fluorescence-based measurement, such as laser-induced fluorescence, was employed. To measure the absolute thickness of the thin fluorescent film, calibration of the measurement system is necessary. Therefore, a theoretical model was proposed considering the multiple reflections of excitation light and fluorescence at the inclined surface between the indenter and workpiece. By measuring the profile of the surface topography of the indented workpiece and comparing the results with those measured by a surface profiler, the validity of the proposed calibration method and the performance of this measurement system were demonstrated. The measured surface profiles, including scratches of 2–4 µm, were in good agreement, demonstrating the validity of the proposed method.

Plastic zones for ductile layer sandwiched between two different substrates

This study is an extension of previously developed models of plasticity occurring in the vicinity of ​​a longitudinal crack in a thin elastic-plastic layer sandwiched between two identical, ideally elastic substrates: Varias–Suo–Shih [1], Rogowski-Molski [2], and Rogowski [3]. The solution presented in the current study describes the plastic zone size problem in a thin constrained layer bonded to two different substrates. A well-known two-term K-T Williams [4] stress solution for interfacial cracks is used to obtain the plastic zone size. A new equation with coefficients dependent on the layer and substrate material properties was obtained, allowing for the calculation of the plastic zone size in the thin ductile layer based on the values of the complex stress intensity factor K and T-stresses. The proposed model is tested for several arbitrarily selected interfacial crack geometries, load conditions, bonding layers, and substrate material properties. Because of the complicated oscillatory description of the local stress field, the final solution of the mathematical equation is obtained numerically, ignoring a small zone of intense stress oscillations near the crack tip. Good agreement is obtained between the results from the proposed model and the FEM results in terms of quantity and quality. A novel approach may be utilized in the future as a method for estimating the parameter of fracture criterion for bonded joints.

Distributed Physics-Informed machine learning strategies for two-phase flows

This study delves into applying Distributed Learning Machines (DLM), a subset of Physics-Informed Neural Networks (PINN), in tackling benchmark challenges within two-phase flows. Specifically, two variants of DLM, namely Distributed Physics-Informed Neural Networks (DPINN) and Transfer Physics-Informed Neural Networks (TPINN), are studied. The DLM architecture strategically divides the global domain into distinct non-overlapping sub-domains, with interconnected solutions facilitated by interface conditions embedded in the loss function. Forward and inverse benchmark problems in two-phase flows are explored: (a) bubble in a reversing vortex and (b) Bubble Rising under Buoyancy. The Volume of Fluid (VOF) method handles interface transport in both scenarios, with the inverse problem incorporating interface-position data during the training phase. The forward problem highlights the effectiveness of DPINN in capturing the interface using a simple transport equation. The distinctive contribution of this work lies in its exploration of the inverse problem, offering insights into the scalability of distributed architectures when dealing with a system of governing equations. Following the validation of an initial PINN model against Computational Fluid Dynamics (CFD) data, the study extends to DPINN and TPINN. A parametric study optimizes network hyperparameters, emphasizing the regularization of loss terms within the DPINN loss function. A self-adaptive weighting strategy based on a Gaussian probabilistic model dynamically adjusts loss weights during training to overcome challenges associated with manual parameter tuning. The evaluation of accuracy against CFD data and published results underscore the efficacy of DLMs in addressing two-phase flow problems. Additionally, the computational efficiency of distributed networks is explored compared to traditional PINNs.

Improved Smoothed Particle Hydrodynamics for the Three-Dimensional Solid Impact Problems

Total Lagrangian-Smoothed Particle Hydrodynamics (TL-SPH) method is an important method to improve the tensile instability of traditional SPH, but it still faces some defects, such as insufficient accuracy at boundary region, difficulty in obtaining damage and applying interface conditions. Therefore, the high-order TL-SPH is applied for the elastic–plastic solid. A damage model is proposed by adaptively transforming the type of kernel function. Furthermore, a Riemann contact algorithm is introduced for different solids. The advantages of the methods in improving tensile instability and simulating damage are verified by the cases of the plates and the notched beams with an impact load.

A unified-field monolithic fictitious domain-finite element method for fluid-structure-contact interactions and applications to deterministic lateral displacement problems

Based upon two overlapped, body-unfitted meshes, a type of unified-field monolithic fictitious domain-finite element method (UFMFD-FEM) is developed in this paper for moving interface problems of dynamic fluid-structure interactions (FSI), which is accompanied with high-contrast physical coefficients across the interface and with contacting collisions between the structure and fluidic channel wall in the immersed FSI case. In particular, the proposed novel numerical method consists of a monolithic, stabilized mixed finite element method within the frame of fictitious domain/immersed boundary method (FD/IBM) for generic fluid-structure-contact interaction (FSCI) problems in the Eulerian–updated Lagrangian description, while involving the no-slip type of interface conditions on the fluid-structure interface, and the repulsive contact force/stress on the structural surface when the immersed structure contacts the fluidic channel wall. The developed UFMFD-FEM for FSI/FSCI problems can deal with the structural motion with large rotational and translational displacements and/or large deformation in an accurate and efficient fashion, which is first validated by two benchmark FSI problems and one FSCI model problem, then by experimental results of a realistic FSCI scenario – the microfluidic deterministic lateral displacement (DLD) problem that is applied to isolate circulating tumor cells (CTCs) from blood cells in the blood fluid through a cascaded filter DLD microchip in practice, where a particulate fluid with the pillar obstacles effect in the fluidic channel, i.e., the effects of fluid-structure interaction and structure collision, play significant roles to sort particles (cells) of different sizes with tilted pillar arrays. The developed unified-field, monolithic fictitious domain-based mixed finite element method can be seamlessly extended to more sophisticated, high dimensional FSCI problems with contacting collisions between the moving elastic structure and fluidic channel wall.

Potential-Modulated Surface-Enhanced Raman Spectroscopy of Tolmetin at Gold Nanoparticle Film Functionalized Polarizable Liquid-Liquid Interfaces.

An aqueous colloidal suspension of gold nanoparticles (AuNPs) may be condensed into a thin fractal film at the polarizable liquid-liquid interface formed between two immiscible electrolyte solutions upon injection of millimolar concentrations of sodium chloride to the aqueous phase. By adjusting the interfacial polarization conditions (negative, intermediate, and positive open-circuit potentials), the morphology of the film is modified, resulting in unique surface plasmon properties of the film, which enable in situ surface-enhanced Raman spectroscopy (SERS). Intense SERS signals are observed at the polarizable liquid-liquid interface when micromolar concentrations of tolmetin, a nonsteroidal anti-inflammatory drug, are entrapped in the AuNP fractal film. The change in the signal intensity, averaged over multiple spectra, with respect to the concentration of tolmetin, depends on the polarization conditions and suggests the presence of chemical-induced damping effects on the surface plasmons of the gold film.

Investigation of effects of initial interface conditions on the two-dimensional single-mode compressible Rayleigh–Taylor instability: Based on the discrete Boltzmann method

The Rayleigh–Taylor (RT) instability in inertial confinement fusion implosions evolves at the unstable interface of two fluids when the light fluid is pushing the heavy one. The effects of the initial amplitude and transition layer on the compressible RT instability are investigated numerically by using the discrete Boltzmann method. On the one hand, during the RT evolution, higher initial amplitudes initially increase the global density gradient and non-equilibrium area, with a subsequent reversal. The increasing initial amplitude leads to an initial rise followed by a decline in the system’s maximum Mach number. On the other hand, the impact of the transition layer is generally opposite to the one of initial amplitude in the RT process. These findings offer significant insights into controlling and understanding RT instability in fusion implosion scenarios, emphasizing novel aspects relative to existing literature.

Viscoelastic modelling and analysis of two-dimensional woven CNT-based multiscale fibre reinforced composite material system

Carbon nanotube (CNT) has fostered research as a promising nanomaterial for a variety of applications due to its exceptional mechanical, optical, and electrical characteristics. The present article proposes a novel and comprehensive micromechanical framework to assess the viscoelastic properties of a multiscale CNT-reinforced two-dimensional (2D) woven hybrid composite. It also focuses on demonstrating the utilisation of the proposed micromechanics in the dynamic analysis of shell structure. First, the detailed constructional attributes of the proposed trans-scale composite material system are described in detail. Then, according to the nature of the constructional feature, mathematical modelling of each constituent phase or building block’s material properties is established to evaluate the homogenised viscoelastic properties of the proposed composite material system. To highlight the novelty of this study, the viscoelastic characteristics of the modified matrix are developed using the micromechanics method of Mori–Tanaka (MT) in combination with the weak viscoelastic interphase (WI) theory. In the entire micromechanical framework, the CNTs are considered to be randomly oriented. The strength of the material (SOM) approach is used to establish mathematical frameworks for the viscoelastic characteristics of yarns, whereas the unit cell method (UCM) is used to determine the viscoelastic properties of the representative unit cell (RUC). Different numerical results have been obtained by varying the CNT composition, interface conditions, agglomeration, carbon fibre volume percentage, excitation frequency, and temperature. The influences of geometrical parameters like yarn thickness, width, and the gap length to yarn width ratio on the viscoelasticity of such composite material systems are also explored. The current study also addresses the issue of resultant anisotropic viscoelastic properties due to the use of dissimilar yarn thickness. The results of this micromechanical analysis provide valuable insights into the viscoelastic properties of the proposed composite material system and suggest its potential applications in vibration damping. To demonstrate the application of developed novel micromechanics in vibration analysis, as one of the main contributions, comprehensive numerical experiments are conducted on a shell panel. The results show a significant reduction in vibration amplitudes compared to traditional composite materials in the frequency response and transient response analyses. To focus on the aspect of micromechanical behaviour on dynamic response and for the purpose of brevity, only linear strain displacement relationships are considered for dynamic analysis. These insights could inform future research and development in the field of composite materials.

A high order generalized finite difference method for solving the anisotropic elliptic interface problem in static and moving systems

In this paper, a high order meshless numerical scheme based on the Generalized Finite Difference Method (GFDM) is proposed to analyze the anisotropic elliptic interface problem. The GFDM has an advantage in dealing with some anisotropic elliptic problems with complex interfaces, interface conditions with the jump of derivatives. Four problems with eight examples are provided to verify the existence of the good performance of the GFDM for anisotropic elliptic interface problems, including that the simplicity, accuracy, and stability in static and moving systems. Especially, GFDM has the tolerance of the complexity of interface shape, the interface movement and the large jump.

An improved partitioned time stepping method based on modified characteristic FEM for the evolutionary dual-porosity-Navier–Stokes model

In this paper, we study an improved partitioned time stepping method based on the modified characteristic finite element method for the evolutionary dual-porosity-Navier–Stokes coupling model with the Beavers–Joseph interface condition. By referring to the idea of partitioned time stepping, we lead into the half-time steps to separate the original coupling problem into three subproblems, implement two-level decoupling to simplify the complexity of solving coupling systems. Before that, the matrix and microfracture equations on the half-time steps in the dual-porosity media are calculated first, with the aim of replacing the initial values of the decoupling equations. Besides, for the nonlinear convection term in Navier–Stokes equations, the modified characteristic finite element method is chosen to avoid the problem of low computational efficiency. The convergence of the decoupling algorithm is rigorously derived by using the idea of mathematical induction. Finally, the reliability of the theoretical analysis and the rationality of the scheme are verified by numerical experiments.

Microcell and macrocell corrosion kinetics of steel bars in reinforced concrete slabs subjected to alternate wet and dry environments

Electrochemical kinetics of steel bar corrosion in concrete is of complex nature, especially when it is subjected to alternate wet and dry environments with changes in conditions of steel-concrete interface and rust properties. In this paper, a quantitative electrochemical mechanism was proposed based on thermodynamics and anodic polarisation kinetics to account for variations of microcell and macrocell current densities in reinforced concrete (RC) slabs subjected to periodic wetting and drying cycles. Effects of environmental conditions on corrosion process were incorporated into anodic Tafel slope through three quantifiable parameters, namely, corrosion potential, linear polarisation resistance and ferrous ion concentration at the steel bar surface. A numerical model was established based on the proposed mechanism to quantify microcell and macrocell current densities. An experimental programme of two RC slabs with multiple steel bars arranged at different layers was conducted to verify the numerical model. It shows that the model could reasonably predict variations of both microcell and macrocell corrosion under different wet and dry environments. When concrete was moistened, current density was substantially increased compared to that in dry concrete. Macrocell contribution from uncorroded steel bars at different layers was more significant than microcell corrosion on the corroded bar itself, especially in wet environments.

Multiscale finite element analysis of layer interface effects on cracking in semi-flexible pavements at different temperatures

This paper presents comprehensive research work on the crack initiation behavior of Semi Flexible Pavements (SFP) with the help of a one-way coupled multiscale modelling approach. The study focuses on elucidating the influence of layer interface conditions on potential cracking patterns under different temperature conditions. By integrating material properties at both local (mixture) and global (pavement) scales using a local de-homogenization approach, the interlocking effect was addressed within the asphalt binder's mesoscale properties. CT-image based mesostructure was employed to reconstruct the Representative Volume Element (RVE). Cohesive Zone Model (CZM) was utilized to simulate cracking initiation within the RVEs, with simulation results validated against practical engineering inspections. Additionally, a potential macroscopic failure index was introduced for SFP. It was found that the typical macroscopic failure indexes utilized for AC pavement may not be suitable for SFP due to the higher degree of heterogeneity. Analysis showed that SFP layer's coarse aggregate and cement grout at the bottom may experience compressive stress, while asphalt binder beneath the tire load may face high tensile stresses. Regardless of the interface conditions, top-down cracking was found to be the dominating failure mode in SFP. Besides, the vulnerable area in bonded interface SFP shifts from under tire load to the adjacent area with the rise in pavement temperature, while bottom-up cracking may appear only in unbonded interface SFP under low temperature conditions. Subsequent analysis demonstrated that top-down cracking is prominent in bonded interface SFP, especially at higher temperatures, compared to unbonded interface SFP. It is expected that critical observations for SFP highlighted above will certainly help in understanding critical failure modes and, hence, designing and constructing durable SFP structures.