Sweta Narayan Sahu, Sumit Sen, Sourav Hossain, Koeli Ghoshal
{"title":"Unsteady suspended sediment distribution in an ice-covered channel through fractional advection–diffusion equation","authors":"Sweta Narayan Sahu, Sumit Sen, Sourav Hossain, Koeli Ghoshal","doi":"10.1007/s10665-024-10380-0","DOIUrl":"https://doi.org/10.1007/s10665-024-10380-0","url":null,"abstract":"<p>Despite several applications of the fractional advection–diffusion equation (fADE) in studying sediment transport in an open channel flow, its application is limited to apprehending the non-local movement of sediment particles in an ice-covered channel with a steady, uniform flow field. An unsteady fADE is considered where the space term is non-local with a non-integer order and the mathematical model with Caputo fractional derivative is able to estimate the variation of sediment concentration along a vertical as well as with time in the ice-covered channel. An eddy viscosity expression is used, which includes the variation in roughness between the channel bed and ice cover surface. The Chebyshev collocation method and the Euler backward method are used to solve the fADE with the initial and boundary conditions and the convergence of the methods is established. The temporal variation of concentration shows that for a zero initial condition, the concentration profile first increases and then becomes stable after a certain time; for a non-zero initial concentration, the profile decreases with an increase in time and eventually a steady state is achieved. The effect of the order of the fractional derivative on the vertical variation of concentration at different times for zero and non-zero initial concentrations is studied and it is found that the order of the fractional derivative has a greater impact at smaller times. The impact of several parameters on concentration profiles is studied at different times and the validation of the model is done by comparing it with experimental studies under restricted conditions.</p>","PeriodicalId":50204,"journal":{"name":"Journal of Engineering Mathematics","volume":"32 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-06-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141502828","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}
{"title":"Perspective on the description of viscoelastic flows via continuum elastic dumbbell models","authors":"Evgeniy Boyko, Howard A. Stone","doi":"10.1007/s10665-024-10374-y","DOIUrl":"https://doi.org/10.1007/s10665-024-10374-y","url":null,"abstract":"<p>Non-Newtonian fluid mechanics and computational rheology widely exploit elastic dumbbell models such as Oldroyd-B and FENE-P for a continuum description of viscoelastic fluid flows. However, these constitutive equations fail to accurately capture some characteristics of realistic polymers, such as the steady extension in simple shear and extensional flows, thus questioning the ability of continuum-level modeling to predict the hydrodynamic behavior of viscoelastic fluids in more complex flows. Here, we present seven elastic dumbbell models, which include different microstructurally inspired terms, i.e., (i) the finite polymer extensibility, (ii) the conformation-dependent friction coefficient, and (iii) the conformation-dependent non-affine deformation. We provide the expressions for the steady dumbbell extension in shear and extensional flows and the corresponding viscosities for various elastic dumbbell models incorporating different microscopic features. We show the necessity of including these microscopic features in a constitutive equation to reproduce the experimentally observed polymer extension in shear and extensional flows, highlighting their potential significance in accurately modeling viscoelastic channel flow with mixed kinematics.</p>","PeriodicalId":50204,"journal":{"name":"Journal of Engineering Mathematics","volume":"2013 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141528794","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}
{"title":"Numerical simulations of shallow-water sloshing coupled to horizontal vessel motion in the presence of a time-dependent porous baffle","authors":"M. R. Turner","doi":"10.1007/s10665-024-10376-w","DOIUrl":"https://doi.org/10.1007/s10665-024-10376-w","url":null,"abstract":"<p>Shallow-water fluid sloshing in the Lagrangian Particle Path formulation, with the addition of an energy-extracting porous baffle, is simulated numerically using a symplectic numerical scheme which captures, in an essential way, the energy exchange. The fluid motion in a rectangular vessel is dynamically coupled to a surface-piercing porous baffle. The fluid transmission through the baffle is characterized by a nonlinear Darcy–Forchheimer model equation. The numerical scheme is symplectic, based on the implicit-midpoint rule, and thus is strategically designed to maintain the energy partition between the fluid and vessel throughout numerous time steps. Our results demonstrate the non-conservative nature of the system, with the porous baffle effectively dissipating energy from the overall system. Furthermore, we present findings that demonstrate the role of time-periodic variations in baffle porosity on energy dissipation. By manipulating the frequency and magnitude of this time-dependent variability, it is established that a greater amount of energy can be extracted from the system compared with the optimal fixed porosity baffle. These results shed new light on potential strategies for enhancing energy dissipation in such configurations.</p>","PeriodicalId":50204,"journal":{"name":"Journal of Engineering Mathematics","volume":"237 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-06-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141520754","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}
{"title":"Buckling phenomenon of vertical beam/column of variable density carrying a top mass","authors":"Mustafa Turkyilmazoglu","doi":"10.1007/s10665-024-10378-8","DOIUrl":"https://doi.org/10.1007/s10665-024-10378-8","url":null,"abstract":"<p>This study focuses on modeling ideal nonuniform standing beams and towers supporting a constant top mass. We also analyze their dynamical stability, as determining the design parameters influencing their shape and stability holds significant value for structural engineering. Initially, we employ a statical mechanics approach to balance the mechanical and gravitational forces. By solving an initial-value problem, we derive the cross-sectional areas of the columns. Our findings reveal that these areas, rather than the shapes, are the primary contributors to the engineering performance of the columns. Additionally, the top mass acts as a multiplying factor for the cross-sectional areas, and the density distribution along the column determines whether the top should be heavier or lighter. Furthermore, we demonstrate that exponential, parabolic, or linear cross-sections with significantly wider base profiles are crucial for accommodating heavier top loads. Moving on to the dynamical analysis, we consider two ideal tower configurations: FC and SC. Numerical and analytical results reveal that higher modes exhibit shorter amplitudes. FC modes necessitate higher design parameters to resist buckling phenomena, whereas SC modes show lower resistance to vibrational deflections. In terms of stability, a heavier top mass enhances the vertical beam’s stability, while towers with parabolic bases are more susceptible to instabilities.</p>","PeriodicalId":50204,"journal":{"name":"Journal of Engineering Mathematics","volume":"75 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-06-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141520755","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}
{"title":"A Nusselt number correlation for a superhydrophobic solid sphere encapsulated in a perfect plastron","authors":"Zavier Berti, M. R. Flynn","doi":"10.1007/s10665-024-10377-9","DOIUrl":"https://doi.org/10.1007/s10665-024-10377-9","url":null,"abstract":"<p>Surface-attached air bubbles are known to provide lubricating (i.e., drag reducing) benefits but their contribution to inhibiting heat transfer is not as well understood. The present theoretical study considers Stokes flow around a solid sphere and uses matched asymptotic expansions to estimate the degree of thermal insulation offered by an encapsulating air layer of uniform thickness. Key to our analysis is to derive an expression for the Nusselt number in terms of the air layer thickness and the Péclet number, <span>(text{ Pe}_w)</span>, of the surrounding liquid, here assumed to be water. This latter parameter, which characterizes advective to diffusive heat transport, is assumed to be small such that our zeroth- and first-order solutions are, respectively, proportional to <span>(text{ Pe}_w^0)</span> and <span>(text{ Pe}_w^1)</span>. Although small <span>(text{ Pe}_w)</span> favors small free stream velocities, forced convection will dominate over natural convection only if the free stream velocity (and/or the solid sphere radius) exceeds a certain threshold. This requirement constrains the solution space; on the other hand, it is straightforward to generalize our analysis so that it considers fluid pairs other than water and air.</p>","PeriodicalId":50204,"journal":{"name":"Journal of Engineering Mathematics","volume":"27 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-06-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141528911","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}
{"title":"Stabilizing an adverse density difference in the presence of phase change","authors":"Lewis Johns, Ranga Narayanan","doi":"10.1007/s10665-024-10372-0","DOIUrl":"https://doi.org/10.1007/s10665-024-10372-0","url":null,"abstract":"<p>Given two phases in equilibrium in a porous solid, the heavy phase lying above the light phase in a gravitational field, we stabilize this adverse density arrangement by heating from below and derive a formula for how steep the temperature gradient must be to do this. The input temperature gradient has two effects on the stability of our system. Its effect on the heat convection is destabilizing, its effect on the heat conduction at the surface is stabilizing. By directing our attention to the case of zero growth rate, we obtain the critical value of the input temperature gradient as it depends on the permeability of the porous solid, the density difference across the surface, the distance between the planes bounding our system, and the physical properties. Our problem makes connections to the Bénard problem where it has two, one, or no critical points, and to the Rayleigh–Taylor problem where it has no critical points.</p>","PeriodicalId":50204,"journal":{"name":"Journal of Engineering Mathematics","volume":"23 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-06-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141502831","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}
Ranusha Rajakrishnan, Seng Huat Ong, Choung Min Ng
{"title":"Alternating minimization algorithm with a probability generating function-based distance measure","authors":"Ranusha Rajakrishnan, Seng Huat Ong, Choung Min Ng","doi":"10.1007/s10665-024-10349-z","DOIUrl":"https://doi.org/10.1007/s10665-024-10349-z","url":null,"abstract":"<p>The Expectation Maximization (EM) algorithm, a popular method for maximum likelihood estimation of parameters, requires a complete data space and construction of a conditional expectation. For many statistical models, these may not be straightforward. This paper proposes a simpler Alternating Minimization (AM) algorithm using a probability generating function (pgf)-based divergence measure for estimation in univariate and bivariate distributions. The performance of the estimation method is studied for the negative binomial and Neyman Type-A distributions in the univariate setting, while for bivariate cases, the bivariate Poisson and the bivariate negative binomial distributions are considered. Comparison is made with direct optimization of pgf-based divergence measure and maximum likelihood (ML) estimates. Results produced via AM in both simulated and real-life datasets show an improvement in comparison to direct pgf optimization, especially in the bivariate setting, with the execution time showing an improvement for large sample sizes when compared to ML. Goodness-of-fit tests show that the pgf divergence measure with AM estimates mostly perform similarly to the ML estimates in terms of power of the test.</p>","PeriodicalId":50204,"journal":{"name":"Journal of Engineering Mathematics","volume":"46 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-05-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141059987","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}
{"title":"Evolution of acceleration waves in non-ideal relaxing gas subjected to the transverse magnetic field","authors":"G. Nath, V. S. Kadam","doi":"10.1007/s10665-024-10345-3","DOIUrl":"https://doi.org/10.1007/s10665-024-10345-3","url":null,"abstract":"<p>Using the characteristics of the governing quasi-linear system as the referencing coordinate system in the presence of a transverse magnetic field, the evolution of acceleration waves in a non-ideal relaxing gas has been examined along its characteristic path. It is demonstrated that a linear solution in the characteristic plane can behave non-linearly in the physical plane. We have determined the critical amplitude of the initial disturbance; if the initial amplitude of the compressive disturbance is greater than the critical value, the disturbance must culminate into a shock wave, while if it is less than this value, the disturbance will decay, and no shock formation will happen. We establish the criteria for shock generation and the transport equation that governs the development of weak shock waves. Acceleration waves having planar and cylindrical symmetry are analyzed as their steepening, or flattening is investigated as a function of the non-idealness parameter, relaxation parameter, adiabatic index, and magnetic field strength parameter. In both the planar and cylindrical symmetries, the shock formation process is slowed by increasing the relaxation parameter as well as the magnetic field parameter but accelerated by non-idealness and the adiabatic index. In the ideal gas case with adiabatic exponent <span>(gamma = 2)</span>, the magnetic field has no effect on the steepening or flattening of the wavefront in both the planar and cylindrical symmetries.</p>","PeriodicalId":50204,"journal":{"name":"Journal of Engineering Mathematics","volume":"54 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-05-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140885685","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}
Sid Becker, Stefanie Gutschmidt, Bradley Boyd, Dan Zhao
{"title":"Pulsed heating of the self-actuated cantilever: a one-dimensional exact solution investigation of non-axial temperature gradients","authors":"Sid Becker, Stefanie Gutschmidt, Bradley Boyd, Dan Zhao","doi":"10.1007/s10665-024-10359-x","DOIUrl":"https://doi.org/10.1007/s10665-024-10359-x","url":null,"abstract":"<p>Self-actuated bimorph cantilevers are implemented in a variety of micro-electro-mechanical systems. Their tip deflection relies on the unmatched coefficients of thermal expansion between layers. The thermal bimorph phenomenon is dependent on the temperature rise within the cantilever and, while previous studies have investigated variations in the thermal profile along the cantilever length, these have usually neglected variations in the thermal profile along the cantilever thickness. The current study investigates the thermal distribution across the thickness of the cantilever. The exact closed form solution to the one-dimensional problem of heat conduction in the composite (layered) domain subjected to transient volumetric heating is developed using the appropriate Green’s function. This solution is applied to a one-dimensional case study of a 3-layer cantilever with an Aluminium heater, a silicon dioxide resistive layer, and a silicon base layer. The aluminium heater experiences volumetric heating at a rate of 0.2 mW/μm<sup>3</sup> of 5 μs duration at 100 μs intervals (10 kHz with a 1/20 duty cycle). Benchmark solutions of the temperature at select times and positions are provided. It is shown that there are negligible temperature gradients across the cantilever thickness during the heating and the first ~ 5 μs afterward. These short-lived temperature differences are positively biased with the unmatched thermal expansion coefficients between the layers, though their relative influence on bending is not clear. A simple parametric analysis indicates that the relative magnitude of the temperature differences across the cantilever (compared to the overall temperature) decreases substantially with increasing duty cycle.</p>","PeriodicalId":50204,"journal":{"name":"Journal of Engineering Mathematics","volume":"108 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-04-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140831700","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}
{"title":"Droplet impact onto a porous substrate: a Wagner theory for early-stage spreading","authors":"Gavin Moreton, Richard Purvis, Mark J. Cooker","doi":"10.1007/s10665-024-10352-4","DOIUrl":"https://doi.org/10.1007/s10665-024-10352-4","url":null,"abstract":"<p>An analytical model for droplet impact onto a porous substrate is presented, based on Wagner theory. An idealised substrate boundary condition is introduced, mimicking the effect of fluid entry into a genuinely porous substrate. The asymptotic analysis yields a solution for a small porous correction with free-surfaces and pressures compared with the impermeable case. On a global scale, it is found that the impact region on the substrate grows more slowly with porosity included due to loss of mass into the substrate. The spatial distribution of liquid volume flux into the substrate is also described. Locally near the turn-over regions, the expected jetting along the surface is calculated with the same volume flux but the jet is found to be slower and thicker than for an impermeable substrate.</p>","PeriodicalId":50204,"journal":{"name":"Journal of Engineering Mathematics","volume":"15 1","pages":""},"PeriodicalIF":1.3,"publicationDate":"2024-04-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140616390","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}