BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Lin Fu

Lin Fu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dynphysics.comp-phmath.NANumerical Analysiscond-mat.softcond-mat.stat-mechGraphicsmath-phmath.MPphysics.app-ph

Trust snapshot

Quick read

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

15 published item(s)

preprint2022arXiv

A Block-based Adaptive Particle Refinement SPH Method for Fluid-Structure Interaction Problems

The multi-resolution method, e.g., the Adaptive Particle Refinement (APR) method, has been developed to increase the local particle resolution and therefore the solution quality within a pre-defined refinement zone instead of using a globally uniform resolution for Smoothed Particle Hydrodynamics (SPH). However, sometimes, the targeted zone of interest can be varying, and the corresponding topology is very complex, thus the conventional APR method is not able to track these characteristics adaptively. In this study, a novel Block-based Adaptive Particle Refinement (BAPR) method is developed, which is able to provide the necessary local refinement flexibly for any targeted characteristic, and track it adaptively. In BAPR, the so-called activation status of the block array defines the refinement regions, where the transition and activated zones are determined accordingly. A regularization method for the generated particles in the newly activated blocks is developed to render an isotropic distribution of these new particles. The proposed method has been deployed for simulating Fluid-Structure Interaction (FSI) problems. A set of 2D FSI cases have been simulated with the proposed BAPR method, and the performance of the BAPR method is quantified and validated comprehensively. In a word, the BAPR method is viable and potential for complex multi-resolution FSI simulations by tracking any targeted characteristic of interest.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A class of new high-order finite-volume TENO schemes for hyperbolic conservation laws with unstructured meshes

The recently proposed high-order TENO scheme [Fu et al., Journal of Computational Physics, 305, pp.333-359] has shown great potential in predicting complex fluids owing to the novel weighting strategy, which ensures the high-order accuracy, the low numerical dissipation, and the sharp shock-capturing capability. However, the applications are still restricted to simple geometries with Cartesian or curvilinear meshes. In this work, a new class of high-order shock-capturing TENO schemes for unstructured meshes are proposed. Similar to the standard TENO schemes and some variants of WENO schemes, the candidate stencils include one large stencil and several small third-order stencils. Following a strong scale-separation procedure, a tailored novel ENO-like stencil selection strategy is proposed such that the high-order accuracy is restored in smooth regions by selecting the candidate reconstruction on the large stencil while the ENO property is enforced near discontinuities by adopting the candidate reconstruction from smooth small stencils. The nonsmooth stencils containing genuine discontinuities are explicitly excluded from the final reconstruction, leading to excellent numerical stability. Different from the WENO concept, such unique sharp stencil selection retains the low numerical dissipation without sacrificing the shock-capturing capability. The newly proposed framework enables arbitrarily high-order TENO reconstructions on unstructured meshes. For conceptual verification, the TENO schemes with third- to sixth-order accuracy are constructed. Without parameter tuning case by case, the performance of the proposed TENO schemes is demonstrated by examining a set of benchmark cases with broadband flow length scales.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A six-point neuron-based ENO (NENO6) scheme for compressible fluid dynamics

In this work, we introduce a deep artificial neural network (ANN) that can detect locations of discontinuity and build a six-point ENO-type scheme based on a set of smooth and discontinuous training data. While a set of candidate stencils of incremental width is constructed, the ANN instead of a classical smoothness indicator is deployed for an ENO-like sub-stencil selection. A convex combination of the candidate fluxes with the re-normalized linear weights forms the six-point neuron-based ENO (NENO6) scheme. The present methodology is inspired by the work [Fu et al., Journal of Computational Physics 305 (2016): 333-359] where contributions of candidate stencils containing discontinuities are removed from the final reconstruction stencil. The binary candidate stencil classification is performed by a well-trained ANN with high fidelity. The proposed framework shows an improved generality and robustness compared with other ANN-based schemes. The generality and performance of the proposed NENO6 scheme are demonstrated by examining one- and two-dimensional benchmark cases with different governing conservation laws and comparing to those of WENO-CU6 and TENO6-opt schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Compressible Velocity Transformations for Various Noncanonical Wall-Bounded Turbulent Flows

This work assesses several popular transformations for the velocity profile through their application to several types of non-canonical compressible wall-bounded turbulent flows. Specifically, this work explores DNS databases of high-enthalpy boundary layers with dissociation and vibrational excitation, supercritical channel and boundary-layer flows, and adiabatic boundary layers with pressure gradients. The transformations considered include the van Driest [Van Driest, J. Aeronaut. Sci., 18(1951):145-216], Zhang et al. [Zhang et al., Phys. Rev. Lett., 109(2012):054502], Trettel-Larsson [Trettel and Larsson, Phys. Fluids, 28(2016):026102], data-driven [Volpiani et al., Phys. Rev. Fluids, 5(2020):052602], and total-stress-based [Griffin et al., Proc. Natl. Acad. Sci. U.S.A., 118(2021):e2111144118] transformations. The Trettel-Larsson transformation collapses velocity profiles of high-enthalpy temporal boundary layers but not the spatial boundary layers considered. For supercritical channel flows, the Trettel-Larsson transformation also performs well over the entire inner layer. None of the transformations above works for supercritical boundary layers. For all the considered methods, the transformed velocity profiles of boundary layers with weak pressure gradients coincide well with the universal incompressible law of the wall. In summary, all these popular methods fail to deliver uniform performance for non-canonical compressible wall-bounded flows in the logarithmic region, and a more sophisticated version, which accounts for these different physics, is needed. The data-driven and total-stress-based transformations perform well in the viscous sublayer for all the considered flows.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Consistency between the attached-eddy model and the inner-outer interaction model: a study of streamwise wall-shear stress fluctuations in a turbulent channel flow

The inner-outer interaction model (Marusic, Mathis & Hutchins, Science, vol. 329, 2010, 193-196) and the attached-eddy model (Townsend, Cambridge University Press, 1976) are two fundamental models describing the multi-scale turbulence interactions and the organization of energy-containing motions in the logarithmic region of high-Reynolds number wall-bounded turbulence, respectively. In this paper, by coupling the additive description with the attached-eddy model, the generation process of streamwise wall-shear fluctuations, resulting from wall-attached eddies, is portrayed. Then, by resorting to the inner-outer interaction model, the streamwise wall-shear stress fluctuations generated by attached eddies in a turbulent channel flow are isolated. Direct comparison between the statistics from these two models demonstrates that they are consistent to and complement each other. Meanwhile, we further show that the superpositions of attached eddies follow an additive process strictly by verifying the validity of the strong and extended self similarity. Moreover, we propose a Gaussian model to characterize the instantaneous distribution of streamwise wall-shear stress, resulting from the attached-eddy superpositions. These findings are important for developing an advanced reduced-order wall model.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Prediction of aerothermal characteristics of a generic hypersonic inlet flow

The accurate prediction of aerothermal surface loading is of paramount importance for the design of high speed flight vehicles. In this work, we consider the numerical solution of hypersonic flow over a double-finned geometry, representative of the inlet of an air-breathing flight vehicle, characterized by three-dimensional intersecting shock-wave/turbulent boundary-layer interaction at Mach 8.3. High Reynolds numbers ($Re_L \approx 11.6 \times 10^6$ based on free-stream conditions) and the presence of cold walls ($T_w/T_o \approx 0.27$) leading to large near-wall temperature gradients necessitate the use of wall-modeled large-eddy simulation (WMLES) in order to make calculations computationally tractable. The comparison of the WMLES results with experimental measurements shows good agreement in the time-averaged surface heat flux and wall pressure distributions, and the WMLES predictions show reduced errors with respect to the experimental measurements than prior RANS calculations. The favorable comparisons are obtained using an LES wall model based on equilibrium boundary layer approximations despite the presence of numerous non-equilibrium conditions including three dimensionality, shock-boundary layer interactions, and flow separation. Lastly, it is also demonstrated that the use of semi-local eddy viscosity scaling (in lieu of the commonly used van Driest scaling) in the LES wall model is necessary to accurately predict the surface pressure loading and heat fluxes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Streamwise inclination angle of wall-attached eddies in turbulent channel flows

We develop a new methodology to assess the streamwise inclination angles (SIAs) of the wall-attached eddies populating the logarithmic region with a given wall-normal height. To remove the influences originating from other scales on the SIA estimated via two-point correlation, the footprints of the targeted eddies in the vicinity of the wall and the corresponding streamwise velocity fluctuations carried by them are isolated simultaneously, by coupling the spectral stochastic estimation with the attached-eddy hypothesis. Datasets produced with direct numerical simulations spanning $Re_τ \sim O(10^2)-O(10^3)$ are dissected to study the Reynolds-number effect. The present results show, for the first time, that the SIAs of attached eddies are Reynolds-number dependent in low and medium Reynolds numbers and tend to saturate at $45^{\circ}$ as the Reynolds number increases. The mean SIA reported by vast previous experimental studies are demonstrated to be the outcomes of the additive effect contributed by multi-scale attached eddies. These findings clarify the long-term debate and perfect the picture of the attached-eddy model.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

General method for determining the boundary layer thickness in nonequilibrium flows

While the computation of the boundary-layer thickness is straightforward for canonical equilibrium flows, there are no established definitions for general non-equilibrium flows. In this work, a method is developed based on a local reconstruction of the "inviscid" velocity profile $U_I[y]$ resulting from the application of the Bernoulli equation in the wall-normal direction. The boundary-layer thickness $δ_{99}$ is then defined as the location where $U/U_I = 0.99$, which is consistent with its classical definition for the zero-pressure-gradient boundary layers (ZPGBLs). The proposed local-reconstruction method is parameter free and can be deployed for both internal and external flows without resorting to an iterative procedure, numerical integration, or numerical differentiation. The superior performance of the local-reconstruction method over various existing methods is demonstrated by applying the methods to laminar and turbulent boundary layers and two flows over airfoils. Numerical experiments reveal that the local-reconstruction method is more accurate and more robust than existing methods, and it is applicable for flows over a wide range of Reynolds numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Feature-aware SPH for Isotropic Unstructured Mesh Generation

In this paper, we present a feature-aware SPH method for the concurrent and automated isotropic unstructured mesh generation. Two additional objectives are achieved with the proposed method compared to the original SPH-based mesh generator (Fu et al., 2019). First, a feature boundary correction term is introduced to address the issue of incomplete kernel support at the boundary vicinity. The mesh generation of feature curves, feature surfaces and volumes can be handled concurrently without explicitly following a dimensional sequence. Second, a two-phase model is proposed to characterize the mesh-generation procedure by a feature-size-adaptation phase and a mesh-quality-optimization phase. By proposing a new error measurement criterion and an adaptive control system with two sets of simulation parameters, the objectives of faster feature-size adaptation and local mesh-quality improvement are merged into a consistent framework. The proposed method is validated with a set of 2D and 3D numerical tests with different complexities and scales. The results demonstrate that high-quality meshes are generated with a significant speedup of convergence.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A low-dissipation shock-capturing framework with flexible nonlinear dissipation control

In this work, a framework to construct arbitrarily high-order low-dissipation shock-capturing schemes with flexible and controllable nonlinear dissipation for convection-dominated problems is proposed. While a set of candidate stencils of incremental width is constructed, each one is indicated as smooth or nonsmooth by the ENO-like stencil selection procedure proposed in the targeted essentially non-oscillatory (TENO) scheme [Fu et al., Journal of Computational Physics 305 (2016): 333-359]. Rather than being discarded directly as with TENO schemes, the nonsmooth candidates are filtered by an extra nonlinear limiter, such as a monotonicity-preserving (MP) limiter or a total variation diminishing (TVD) limiter. Consequently, high-order reconstruction is achieved by assembling candidate fluxes with optimal linear weights since they are either smooth reconstructions or filtered ones which feature good non-oscillation property. A weight renormalization procedure as with the standard TENO paradigm is not necessary. This new framework concatenates the concepts of TENO, WENO and other nonlinear limiters for shock-capturing, and provides a new insight to designing low-dissipation nonlinear schemes. Through the adaptation of nonlinear limiters, nonlinear dissipation in the newly proposed framework can be controlled separately without affecting the performance in smooth regions. Based on the proposed framework, a family of new six- and eight-point nonlinear schemes with controllable dissipation is proposed. A set of critical benchmark cases involving strong discontinuities and broadband fluctuations is simulated. Numerical results reveal that the proposed new schemes capture discontinuities sharply and resolve the high-wavenumber fluctuations with low dissipation, while maintaining the desired accuracy order in smooth regions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Shock-induced heating and transition to turbulence in a hypersonic boundary layer

The interaction between an incident shock wave and a Mach-6 undisturbed hypersonic laminar boundary layer over a cold wall is addressed using direct numerical simulations (DNS) and wall-modeled large-eddy simulations (WMLES) at different angles of incidence. At sufficiently high shock-incidence angles, the boundary layer transitions to turbulence via breakdown of near-wall streaks shortly downstream of the shock impingement, without the need of any inflow free-stream disturbances. The transition causes a localized significant increase in the Stanton number and skin-friction coefficient, with high incidence angles augmenting the peak thermomechanical loads in an approximately linear way. Statistical analyses of the boundary layer downstream of the interaction for each case are provided that quantify streamwise spatial variations of the Reynolds analogy factors and indicate a breakdown of the Morkovin's hypothesis near the wall, where velocity and temperature become correlated. A modified strong Reynolds analogy with a fixed turbulent Prandtl number is observed to perform best. Conventional transformations fail at collapsing the mean velocity profiles on the incompressible log law. The WMLES prompts transition and peak heating, delays separation, and advances reattachment, thereby shortening the separation bubble. When the shock leads to transition, WMLES provides predictions of DNS peak thermomechanical loads within $\pm 10\%$ at a computational cost lower than DNS by two orders of magnitude. Downstream of the interaction, in the turbulent boundary layer, WMLES agrees well with DNS results for the Reynolds analogy factor, the mean profiles of velocity and temperature, including the temperature peak, and the temperature/velocity correlation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Time-Accurate and highly-Stable Explicit operators for stiff differential equations

Unconditionally stable implicit time-marching methods are powerful in solving stiff differential equations efficiently. In this work, a novel framework to handle stiff physical terms implicitly is proposed. Both physical and numerical stiffness originating from convection, diffusion and source terms (typically related to reaction) can be handled by a set of predefined Time-Accurate and highly-Stable Explicit (TASE) operators in a unified framework. The proposed TASE operators act as preconditioners on the stiff terms and can be deployed to any existing explicit time-marching methods straightforwardly. The resulting time integration methods remain the original explicit time-marching schemes, yet with nearly unconditional stability. The TASE operators can be designed to be arbitrarily high-order accurate with Richardson extrapolation such that the accuracy order of original explicit time-marching method is preserved. Theoretical analyses and stability diagrams show that the $s$-stages $s$th-order explicit Runge-Kutta (RK) methods are unconditionally stable when preconditioned by the TASE operators with order $p \leq s$ and $p \leq 2$. On the other hand, the $s$th-order RK methods preconditioned by the TASE operators with order of $p \leq s$ and $p > 2$ are nearly unconditionally stable. The only free parameter in TASE operators can be determined a priori based on stability arguments. Unlike classical implicit methods, the TASE methodology allows for solving non-linear problems with arbitrary order without requiring solving a nonlinear system of equations. A set of benchmark problems with strong stiffness is simulated to assess the performance of the TASE method. Numerical results suggest that the proposed framework preserves the high-order accuracy of the explicit time-marching methods with very-large time steps for all the considered cases.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A Consistent Parallel Isotropic Unstructured Mesh Generation Method based on Multi-phase SPH

In this paper, we propose a consistent parallel unstructured mesh generator based on a multi-phase SPH method. A set of physics-motivated modeling equations are developed to achieve the targets of domain decomposition, communication volume optimization and high-quality unstructured mesh generation simultaneously. A unified density field is defined as the target function for both partitioning the geometry and distributing the mesh-vertexes. A multi-phase Smoothing Particle Hydrodynamics (SPH) method is employed to solve the governing equations. All the optimization targets are achieved implicitly and consistently by the particle relaxation procedure without constructing triangulation/tetrahedralization explicitly. The target of communication reduction is achieved by introducing a surface tension model between distinct partitioning sub-domains, which are characterized by colored SPH particles. The resulting partitioning diagram features physically localized sub-domains and optimized interface communication. The target of optimizing the mesh quality is achieved by introducing a tailored equation-of-state (EOS) and a smooth isotropic kernel function. The mesh quality near the interface of neighboring sub-domains is improved by gradually removing the surface-tension force once a steady state is achieved. The proposed method is developed basing on a new parallel environment for multi-resolution SPH to exploit both coarse- and fine-grained parallelism. A set of benchmarks are conducted to verify that all the optimization targets are achieved consistently within the current framework.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Correlation lengths in quasi-one-dimensional systems via transfer matrices

Using transfer matrices up to next-nearest-neighbour (NNN) interactions, we examine the structural correlations of quasi-one-dimensional systems of hard disks confined by two parallel lines and hard spheres confined in cylinders. Simulations have shown that the non-monotonic and non-smooth growth of the correlation length in these systems accompanies structural crossovers (Fu et al., Soft Matter, 2017, 13, 3296). Here, we identify the theoretical basis for these behaviour. In particular, we associate kinks in the growth of correlation lengths with eigenvalue crossing and splitting. Understanding the origin of such structural crossovers answers questions raised by earlier studies, and thus bridges the gap between theory and simulations for these reference models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Assembly of hard spheres in a cylinder: a computational and experimental study

Hard spheres are an important benchmark of our understanding of natural and synthetic systems. In this work, colloidal experiments and Monte Carlo simulations examine the equilibrium and out-of-equilibrium assembly of hard spheres of diameter $σ$ within cylinders of diameter $σ\leq D\leq 2.82σ$. Although in such a system phase transitions formally do not exist, marked structural crossovers are observed. In simulations, we find that the resulting pressure-diameter structural diagram echoes the densest packing sequence obtained at infinite pressure in this range of $D$. We also observe that the out-of-equilibrium self-assembly depends on the compression rate. Slow compression approximates equilibrium results, while fast compression can skip intermediate structures. Crossovers for which no continuous line-slip exists are found to be dynamically unfavorable, which is the source of this difference. Results from colloidal sedimentation experiments at high Péclet number are found to be consistent with the results of fast compressions, as long as appropriate boundary conditions are used. The similitude between compression and sedimentation results suggests that the assembly pathway does not here sensitively depend on the nature of the out-of-equilibrium dynamics.

0 citations0 reviewsNo code linkedOpen access