Pub Date : 2025-06-12DOI: 10.1140/epje/s10189-025-00496-1
Priyam Chakraborty, Rahul Roy, Shubhadeep Mandal
Artificial microswimmers, such as active colloids, have the potential to revolutionize targeted drug delivery, but controlling their motion under imposed flow conditions remains challenging. In this work, we implement reinforcement learning (RL) to control the navigation of a microswimmer in a plane Poiseuille flow, with applications in targeted drug delivery. With RL, the swimmer learns to efficiently reach its target by continuously adjusting its swinging or tumbling behavior depending upon its self-propulsion strength, chirality and the imposed flow strength. This RL-based approach enables precise control of the particle’s path, achieving reliable targeting even in stringent scenarios such as upstream motion in high bulk flow, thus advancing the design of intelligent in vivo medical microrobots.
{"title":"Smart navigation of microswimmers in Poiseuille flow via reinforcement learning","authors":"Priyam Chakraborty, Rahul Roy, Shubhadeep Mandal","doi":"10.1140/epje/s10189-025-00496-1","DOIUrl":"10.1140/epje/s10189-025-00496-1","url":null,"abstract":"<p>Artificial microswimmers, such as active colloids, have the potential to revolutionize targeted drug delivery, but controlling their motion under imposed flow conditions remains challenging. In this work, we implement reinforcement learning (RL) to control the navigation of a microswimmer in a plane Poiseuille flow, with applications in targeted drug delivery. With RL, the swimmer learns to efficiently reach its target by continuously adjusting its swinging or tumbling behavior depending upon its self-propulsion strength, chirality and the imposed flow strength. This RL-based approach enables precise control of the particle’s path, achieving reliable targeting even in stringent scenarios such as upstream motion in high bulk flow, thus advancing the design of intelligent <i>in vivo</i> medical microrobots.</p>","PeriodicalId":790,"journal":{"name":"The European Physical Journal E","volume":"48 6-7","pages":""},"PeriodicalIF":2.2,"publicationDate":"2025-06-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144273923","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-05-27DOI: 10.1140/epje/s10189-025-00494-3
Supti Sadhukhan, Tapati Dutta
This study investigates the impact of clay content and temperature variation on the electrical conductivity of three-dimensional fluid-filled porous rocks. The role of varying pore throat radii has been included in the course of clay fraction variation in the conducting channels of the rock samples. The research identifies a critical ratio of clay conductance to fluid conductance that dictates the regime of electrical conductance behaviour. A nonlinear increase in electrical conductance is observed when the clay-to-fluid conductance ratio exceeds the critical ratio, whereas a linear relationship is maintained below this critical ratio. A modified form of Archie’s law relating effective conductivity and porosity has been proposed for the clay coated channels. The intricate relationship between Peclet number, pore throat size, and temperature on the electrical conductivity of fluid-filled straight channels in three dimensions has also been investigated. Results revealed a quadratic increase in conductance with porosity under steady-state conditions across all Peclet number ranges examined. While the conductivity remained constant with porosity for each Peclet number, the rate of increase in conductivity diminished with it. Nonlinear increase in conductivity was observed with temperature in the transient flow regime with a threshold temperature marking the onset of conductivity. Conductivity was augmented with increase in observation time in the transient state for the entire temperature range considered. Close to the attainment of saturation in electrical conductivity, the conductivity changed linearly with temperature until a steady value was reached.
{"title":"A simulation study of electrical conductivity of porous rocks: effect of clay, porosity, temperature and Peclet number","authors":"Supti Sadhukhan, Tapati Dutta","doi":"10.1140/epje/s10189-025-00494-3","DOIUrl":"10.1140/epje/s10189-025-00494-3","url":null,"abstract":"<p>This study investigates the impact of clay content and temperature variation on the electrical conductivity of three-dimensional fluid-filled porous rocks. The role of varying pore throat radii has been included in the course of clay fraction variation in the conducting channels of the rock samples. The research identifies a critical ratio of clay conductance to fluid conductance that dictates the regime of electrical conductance behaviour. A nonlinear increase in electrical conductance is observed when the clay-to-fluid conductance ratio exceeds the critical ratio, whereas a linear relationship is maintained below this critical ratio. A modified form of Archie’s law relating effective conductivity and porosity has been proposed for the clay coated channels. The intricate relationship between Peclet number, pore throat size, and temperature on the electrical conductivity of fluid-filled straight channels in three dimensions has also been investigated. Results revealed a quadratic increase in conductance with porosity under steady-state conditions across all Peclet number ranges examined. While the conductivity remained constant with porosity for each Peclet number, the rate of increase in conductivity diminished with it. Nonlinear increase in conductivity was observed with temperature in the transient flow regime with a threshold temperature marking the onset of conductivity. Conductivity was augmented with increase in observation time in the transient state for the entire temperature range considered. Close to the attainment of saturation in electrical conductivity, the conductivity changed linearly with temperature until a steady value was reached.</p>","PeriodicalId":790,"journal":{"name":"The European Physical Journal E","volume":"48 4-5","pages":""},"PeriodicalIF":2.2,"publicationDate":"2025-05-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144155488","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-05-24DOI: 10.1140/epje/s10189-025-00489-0
Nadine du Toit, Kristian K. Müller-Nedebock
A novel field theoretical approach towards modelling dynamic networking in complex systems is presented. An equilibrium networking formalism which utilises Gaussian fields is adapted to model the dynamics of particles that can bind and unbind from one another. Here, networking refers to the introduction of instantaneous co-localisation constraints and does not necessitate the formation of a well-defined transient or persistent network. By combining this formalism with Martin–Siggia–Rose generating functionals, a weighted generating functional for the networked system is obtained. The networking formalism introduces spatial and temporal constraints into the Langevin dynamics, via statistical weights, thereby accounting for all possible configurations in which particles can be networked to one another. A simple example of Brownian particles which can bind and unbind from one another demonstrates the tool and that this leads to results for physical quantities in a collective description. Applying the networking formalism to model the dynamics of cross-linking polymers in a mixture, we can calculate the average number of networking instances. As expected, the dynamic structure factors for each type of polymer show that the system collapses once networking is introduced, but that the addition of a repulsive time-dependent potential above a minimum strength prevents this. The examples presented in this paper indicate that this novel approach towards modelling dynamic networking could be applied to a range of synthetic and biological systems to obtain theoretical predictions for experimentally verifiable quantities.
{"title":"Dynamical networking using Gaussian fields","authors":"Nadine du Toit, Kristian K. Müller-Nedebock","doi":"10.1140/epje/s10189-025-00489-0","DOIUrl":"10.1140/epje/s10189-025-00489-0","url":null,"abstract":"<p>A novel field theoretical approach towards modelling dynamic networking in complex systems is presented. An equilibrium networking formalism which utilises Gaussian fields is adapted to model the dynamics of particles that can bind and unbind from one another. Here, <i>networking</i> refers to the introduction of instantaneous co-localisation constraints and does not necessitate the formation of a well-defined transient or persistent network. By combining this formalism with Martin–Siggia–Rose generating functionals, a weighted generating functional for the networked system is obtained. The networking formalism introduces spatial and temporal constraints into the Langevin dynamics, via statistical weights, thereby accounting for all possible configurations in which particles can be networked to one another. A simple example of Brownian particles which can bind and unbind from one another demonstrates the tool and that this leads to results for physical quantities in a collective description. Applying the networking formalism to model the dynamics of cross-linking polymers in a mixture, we can calculate the average number of networking instances. As expected, the dynamic structure factors for each type of polymer show that the system collapses once networking is introduced, but that the addition of a repulsive time-dependent potential above a minimum strength prevents this. The examples presented in this paper indicate that this novel approach towards modelling dynamic networking could be applied to a range of synthetic and biological systems to obtain theoretical predictions for experimentally verifiable quantities.</p>","PeriodicalId":790,"journal":{"name":"The European Physical Journal E","volume":"48 4-5","pages":""},"PeriodicalIF":1.8,"publicationDate":"2025-05-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1140/epje/s10189-025-00489-0.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144125868","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-05-24DOI: 10.1140/epje/s10189-025-00492-5
Touvia Miloh, Eldad J. Avital
A unified ‘weak-field’ formulation is provided for calculating the combined nonlinear effect of dielectrophoresis and the induced-charge electrophoresis (dipolophoresis) of polarized rigid hydrophobic spherical colloids freely suspended in an electrolyte-saturated Brinkman-hydrogel (porous) medium under a general (direct or alternating currents) non-uniform electric forcing. Explicit expressions for the modified total dipolophoretic mobility of a conducting (metallic) spherical colloid are given in terms of the Brinkman (Darcy), Navier slip, and Debye (screening) length scales. Also presented is a rigorous derivation of the Helmholtz–Smoluchowski slip velocity in terms of these three length scales, including the induced electroosmotic flow field around a hydrophobic rigid colloid embedded in a Brinkman medium that is forced by an arbitrary (non-uniform) ambient electric field. The available solutions for a free (non-porous) electrolyte solution under a uniform forcing and no-slip surface are obtained as limiting cases. For the purpose of illustration, we present and analyse some newly explicit solutions for the mobility and the associated induced-charge electroosmotic velocity field of a slipping colloid set in an effective (hydrogel) porous medium, which is exposed to an ambient ‘sinusoidal’ travelling-wave excitation depending on frequency and wave number.
{"title":"Travelling-wave gel dipolophoresis of hydrophobic conducting colloids","authors":"Touvia Miloh, Eldad J. Avital","doi":"10.1140/epje/s10189-025-00492-5","DOIUrl":"10.1140/epje/s10189-025-00492-5","url":null,"abstract":"<div><p>A unified ‘weak-field’ formulation is provided for calculating the combined nonlinear effect of dielectrophoresis and the induced-charge electrophoresis (dipolophoresis) of polarized rigid hydrophobic spherical colloids freely suspended in an electrolyte-saturated Brinkman-hydrogel (porous) medium under a general (direct or alternating currents) non-uniform electric forcing. Explicit expressions for the modified total dipolophoretic mobility of a conducting (metallic) spherical colloid are given in terms of the Brinkman (Darcy), Navier slip, and Debye (screening) length scales. Also presented is a rigorous derivation of the Helmholtz–Smoluchowski slip velocity in terms of these three length scales, including the induced electroosmotic flow field around a hydrophobic rigid colloid embedded in a Brinkman medium that is forced by an arbitrary (non-uniform) ambient electric field. The available solutions for a free (non-porous) electrolyte solution under a uniform forcing and no-slip surface are obtained as limiting cases. For the purpose of illustration, we present and analyse some newly explicit solutions for the mobility and the associated induced-charge electroosmotic velocity field of a slipping colloid set in an effective (hydrogel) porous medium, which is exposed to an ambient ‘sinusoidal’ travelling-wave excitation depending on frequency and wave number.</p><h3>Graphical abstract</h3><div><figure><div><div><picture><source><img></source></picture></div></div></figure></div></div>","PeriodicalId":790,"journal":{"name":"The European Physical Journal E","volume":"48 4-5","pages":""},"PeriodicalIF":1.8,"publicationDate":"2025-05-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1140/epje/s10189-025-00492-5.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144131660","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-05-22DOI: 10.1140/epje/s10189-025-00488-1
Laura Lavacchi, Benjamin A. Dalton, Roland R. Netz
Barrier-crossing processes in nature are often non-Markovian and typically occur over an asymmetric double-well free-energy landscape. However, most theories and numerical studies on barrier-crossing rates assume symmetric free-energy profiles. Here, we use a one-dimensional generalized Langevin equation (GLE) to investigate non-Markovian reaction kinetics in asymmetric double-well potentials. We derive a general formula, confirmed by extensive simulations, that accurately predicts mean first-passage times from well to barrier top in an asymmetric double-well potential with arbitrary memory time and reaction coordinate mass. We extend our formalism to non-equilibrium non-Markovian systems, confirming its broad applicability to equilibrium and non-equilibrium systems in biology, chemistry, and physics.
{"title":"Non-Markovian equilibrium and non-equilibrium barrier-crossing kinetics in asymmetric double-well potentials","authors":"Laura Lavacchi, Benjamin A. Dalton, Roland R. Netz","doi":"10.1140/epje/s10189-025-00488-1","DOIUrl":"10.1140/epje/s10189-025-00488-1","url":null,"abstract":"<p>Barrier-crossing processes in nature are often non-Markovian and typically occur over an asymmetric double-well free-energy landscape. However, most theories and numerical studies on barrier-crossing rates assume symmetric free-energy profiles. Here, we use a one-dimensional generalized Langevin equation (GLE) to investigate non-Markovian reaction kinetics in asymmetric double-well potentials. We derive a general formula, confirmed by extensive simulations, that accurately predicts mean first-passage times from well to barrier top in an asymmetric double-well potential with arbitrary memory time and reaction coordinate mass. We extend our formalism to non-equilibrium non-Markovian systems, confirming its broad applicability to equilibrium and non-equilibrium systems in biology, chemistry, and physics.</p>","PeriodicalId":790,"journal":{"name":"The European Physical Journal E","volume":"48 4-5","pages":""},"PeriodicalIF":1.8,"publicationDate":"2025-05-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1140/epje/s10189-025-00488-1.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144117701","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-05-09DOI: 10.1140/epje/s10189-025-00490-7
Zhiping Li, Long He, Tianyu Pan, Yao Yin, Shaobin Li, Wei Yuan, Bo Meng
Previous studies revealed the skin-friction drag reduction properties induced by transverse grooves. However, the effects of unsteady characteristics of vortices within the grooves on the drag reduction properties have not been investigated. A hypothesis that the unsteady motion of vortices may reduce the friction drag-reduction rate induced by transverse V-grooves is proposed in this paper. To verify this hypothesis, we use the LES (large eddy simulation) method to investigate the flow field in the range of Reynolds number 0.5E5 to 7.5E5 over the different profiles of symmetric V-grooves, whose depths are 0.2 mm and AR’s are 0.5, 1, 2, 5, and 8. The results show that the AR (aspect ratio of a transverse groove) affects the stability of boundary vortices, thus driving the variation of total viscous drag and pressure drag. With the increase of AR, the boundary vortices tend to be stable at first and then gradually become unstable. When AR is 2, the boundary vortices are stable within the grooves, corresponding to optimal drag reduction. In this case, the slip velocities induced by boundary vortices are the largest, and the Reynolds shear stress is the least, suggesting that the grooves have the strongest abilities to reduce the total viscous drag. When the stability of the boundary vortices is broken, a larger area containing high pressure and low pressure is formed in the groove, and the difference also becomes greater between the high pressure and low pressure. The results provide improved understandings of the drag reduction mechanism of transverse grooves.