Researcher profile

Richard M. McLaughlin

Richard M. McLaughlin contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 15 - UnverifiedVerification L1Unclaimed author
3works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

3 published item(s)

preprint2023arXiv

Critical density triplets for the arrestment of a sphere falling in a sharply stratified fluid

We study the motion of a rigid sphere falling in a two-layer stratified fluid under the action of gravity in the potential flow regime. Experiments at a moderate Reynolds number of approximately 20 to 450 indicate that a sphere with the precise critical density, higher than the bottom layer density, can display behaviors such as bounce or arrestment after crossing the interface. We experimentally demonstrate that such a critical sphere density increases linearly as the bottom fluid density increases with a fixed top fluid density. Additionally, the critical density approaches the bottom layer fluid density as the thickness of density transition layer increases. We propose an estimation of the critical density based on the potential energy. With assuming the zero layer thickness, the estimation constitutes an upper bound of the critical density with less than 0.043 relative difference within the experimental density regime 0.997 $g/cm^{3}$ $\sim $ 1.11 $g/cm^{3}$ under the zero layer thickness assumption. By matching the experimental layer thickness, we obtain a critical density estimation with less than 0.01 relative difference within the same parameter regime.

preprint2022arXiv

Correlation function of a random scalar field evolving with a rapidly fluctuating Gaussian process

We consider a scalar field governed by an advection-diffusion equation (or a more general evolution equation) with rapidly fluctuating, Gaussian distributed random coefficients. In the white noise limit, we derive the closed evolution equation for the ensemble average of the random scalar field by three different strategies, i.e., Feynman-Kac formula, the limit of Ornstein-Uhlenbeck process, and evaluating the cluster expansion of the propagator on an $n$-simplex. With the evolution equation of ensemble average, we study the passive scalar transport problem with two different types of flows, a random periodic flow, and a random strain flow. For periodic flows, by utilizing the homogenization method, we show that the $N$-point correlation function of the random scalar field satisfies an effective diffusion equation at long times. For the strain flow, we explicit compute the mean of the random scalar field and show that the statistics of the random scalar field have a connection to the time integral of geometric Brownian motion. Interestingly, all normalized moment (e.g., skewness, kurtosis) of this random scalar field diverges at long times, meaning that the scalar becomes more and more intermittent during its decay.

preprint2022arXiv

Numerical Algorithms for Water Waves with Background Flow over Obstacles and Topography

We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply-connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).