Researcher profile

Yulin Pan

Yulin Pan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
3topics
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

6 published item(s)

preprint2022arXiv

Do surface gravity waves have a frozen turbulence state?

We study the energy transfer by exact resonances for surface gravity waves in a finite periodic spatial domain. Based on a kinematic model simulating the generation of active wave modes in a finite discrete wavenumber space $\mathcal{S}_R$, we examine the possibility of direct and inverse energy cascades. More specifically, we set an initially excited region which iteratively spreads energy to wave modes in $\mathcal{S}_R$ through exact resonances. At each iteration, we first activate new modes from scale resonances (which generate modes with new lengths), then consider two bounding situations for angle resonances (which transfer energy at the same length scale): the lower bound where no angle resonance is included and the upper bound where all modes with the same length as any active mode are excited. Such a strategy is essential to enable the computation for a large domain $\mathcal{S}_R$ with the maximum wavenumber $R\sim 10^3$. We show that for both direct and inverse cascades, the energy cascade to the boundaries of $\mathcal{S}_R$ can be established when the initially excited region is sufficiently large, otherwise a frozen turbulence state occurs, with a sharp transition between the two regimes especially for the direct cascade. Through a study on the structure of resonant quartets, the mechanism associated with the sharp transition and the role of angular energy transfer in the cascades are elucidated.

preprint2022arXiv

Effects of varying initial conditions of ship encountering wave groups in computing extreme motion statistics

In computing ship motion statistics (e.g., exceeding probability) in an irregular wave field, it is a common practice to represent the irregular waves by a large number of wave groups and compute the motion statistics from the distribution of these wave groups. While this procedure significantly reduces the computational cost, the uncertainties introduced in this reduced-order computation have not been quantified. In general, the representation of a continuous wave field by separated groups loses information about wave phases, frequency modulation and initial conditions of ship when encountering the wave groups, among which the last one is arguably the most influential factor for the ship motion statistics. In this paper, we test three methods to incorporate the ship initial conditions into the computation of roll motion exceeding probability, namely the methods of natural initial condition, pre-computed initial condition and constant initial condition. For different input wave spectra and ship motion dynamics (in terms of the parametric excitation level modeled by a nonlinear roll equation), the performances of the three methods are carefully benchmarked and the superiority of the natural initial condition method is demonstrated. We finally show that the computation using the natural initial condition method can be greatly accelerated through a sequential sampling algorithm making use of the variational heteroscedastic Gaussian process regression.

preprint2022arXiv

Efficient computation of temporal exceeding probability of ship responses in a random wave field

In this work, we develop a computational framework to efficiently quantify the temporal exceeding probability of ship responses in a random wave field, i.e., the fraction of time that ship responses exceed a given threshold. In particular, we consider large thresholds so that the response exceedance needs to be treated as rare events. Our computational framework builds on the parameterization of wave field into wave groups and efficient sampling in the group parameter space through Bayesian experimental design (BED). Previous works following this framework exclusively studied extreme statistics of group-maximum response, i.e., the maximum response within a wave group, which is however not straightforward to interpret in practice (e.g., its value depends on the subjective definition of the wave group). In order to adapt the framework to the more robust measure by the temporal exceeding probability, novel developments need to be introduced to the two components of BED (on surrogate model and acquisition function), which are detailed in the paper. We validate our framework for a case of ship roll motion calculated by a nonlinear roll equation, in terms of the agreement of our results to the true solution and the independence of our results to the criterion of defining groups. Finally, we demonstrate the coupling of the framework with CFD models to handle more realistic and general ship response problems.

preprint2022arXiv

Multi-fidelity Bayesian experimental design to quantify extreme-event statistics

In this work, we develop a multi-fidelity Bayesian experimental design framework to efficiently quantify the extreme-event statistics of an input-to-response (ItR) system with given input probability and expensive function evaluations. The key idea here is to leverage low-fidelity samples whose responses can be computed with a cost of a certain fraction of that for high-fidelity samples, in an optimized configuration to reduce the total computational cost. To accomplish this goal, we employ a multi-fidelity Gaussian process as the surrogate model of the ItR function, and develop a new acquisition based on which the optimized next sample can be selected in terms of its location in the sample space and the fidelity level. In addition, we develop an inexpensive analytical evaluation of the acquisition and its derivative, avoiding numerical integrations that are prohibitive for high-dimensional problems. The new method is tested in a bi-fidelity context for a series of synthetic problems with varying dimensions, low-fidelity model accuracy and computational costs. Comparing with the single-fidelity method and the bi-fidelity method with a pre-defined fidelity hierarchy, our method consistently shows the best (or among the best) performance for all the test cases. Finally, we demonstrate the superiority of our method in solving an engineering problem of estimating the extreme ship motion statistics in irregular waves, using computational fluid dynamics (CFD) with two different grid resolutions as the high and low fidelity models.

preprint2022arXiv

Sequential Bayesian experimental design for estimation of extreme-event probability in stochastic dynamical systems

We consider an input-to-response (ItR) system characterized by (1) parameterized input with a known probability distribution and (2) stochastic ItR function with heteroscedastic randomness. Our purpose is to efficiently quantify the extreme response probability when the ItR function is expensive to evaluate. The problem setup arises often in physics and engineering problems, with randomness in ItR coming from either intrinsic uncertainties (say, as a solution to a stochastic equation) or additional (critical) uncertainties that are not incorporated in a low-dimensional input parameter space (as a result of dimension reduction applied to the original high-dimensional input space). To reduce the required sampling numbers, we develop a sequential Bayesian experimental design method leveraging the variational heteroscedastic Gaussian process regression (VHGPR) to account for the stochastic ItR, along with a new criterion to select the next-best samples sequentially. The validity of our new method is first tested in two synthetic problems with the stochastic ItR functions defined artificially. Finally, we demonstrate the application of our method to an engineering problem of estimating the extreme ship motion probability in irregular waves, where the uncertainty in ItR naturally originates from standard wave group parameterization, which reduces the original high-dimensional wave field into a two-dimensional parameter space.

preprint2021arXiv

On the Properties of Energy Flux in Wave Turbulence

We study the properties of energy flux in wave turbulence via the Majda-McLaughlin-Tabak (MMT) equation with a quadratic dispersion relation. One of our purposes is to resolve the inter-scale energy flux $P$ in the stationary state to elucidate its distribution and scaling with spectral level. More importantly, we perform a quartet-level decomposition of $P=\sum_ΩP_Ω$, with each component $P_Ω$ representing the contribution from quartet interactions with frequency mismatch $Ω$, in order to explain the properties of $P$ as well as study the wave-turbulence closure model. Our results show that time series of $P$ closely follows a Gaussian distribution, with its standard deviation several times its mean value $\overline{P}$. This large standard deviation is shown to mainly result from the fluctuation (in time) of the quasi-resonances, i.e., $P_{Ω\neq 0}$. The scaling of spectral level with $\overline{P}$ exhibits $\overline{P}^{1/3}$ and $\overline{P}^{1/2}$ at high and low nonlinearity, consistent with the kinetic and dynamic scalings respectively. The different scaling laws in the two regimes are explained through the dominance of quasi-resonances ($P_{Ω\neq 0}$) and exact resonances ($P_{Ω= 0}$) in the former and latter regimes. Finally, we investigate the wave-turbulence closure model, which connects fourth-order correlators to products of pair correlators through a broadening function $f(Ω)$, sometimes argued to be a $sinc$ function in the theory. Our numerical data show that consistent behavior of $f(Ω)$ can only be observed upon averaging over a large number of quartets, but with $f(Ω)$ showing $f\sim 1/Ω^β$ dependence with $β$ taking values between $1.3$ and $1.6$.