Source author record

Spencer J. Sherwin

Spencer J. Sherwin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.flu-dyn math.NA Numerical Analysis physics.comp-ph Computational Geometry math.DS math.OC Mathematical Software

Catalog footprint

What is connected

8works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

preprint2022arXiv

Open-source Framework for Transonic Boundary Layer Natural Transition Analysis over Complex Geometries in Nektar++

We introduce an open-source and unified framework for transition analysis for laminar boundary layer natural transition at transonic conditions and over complex geometries, where surface irregularities may be present. Different computational tools are integrated in the framework, and therefore overcomes the difficulties of two separate and usually quite disparate processes when using $e^N$ method for transition analysis. To generate a baseflow with desired pressure distribution, appropriate pressure compatible inflow boundary condition needs to be developed and enforced. We first derive the system for 1D numerical stability analysis for boundary conditions, and construct three types of pressure compatible inflow. We demonstrate that the entropy-pressure compatible inflow is stable unlike other choices. Compared with the steady baseflow computation, the unsteady simulation for the disturbance field is more challenging for compressible flows because of complex wave reflections, which can easily contaminate the results. We therefore introduce the two main sources of wave decontamination and corresponding methods to obtain clean signal. The workflow within the framework is then verified by computing the disturbance development in 2D flat plate boundary layer flows at Mach $0.8$. The $N$-factors over a clean flat plate and a flat plate with a forward-facing step are generated, and agree well with the results from the reference. Following the verified workflow, We then analyze the disturbance growth on a wing section of the CRM-NLF model. The N-factor on a 2D simulation is generated and studied.

preprint2020arXiv

Nektar++: Design and implementation of an implicit, spectral/$hp$ element, compressible flow solver using a Jacobian-free Newton Krylov approach

At high Reynolds numbers, the use of explicit in time compressible flow simulations with spectral/$hp$ element discretization can become significantly limited by time step. To alleviate this limitation we extend the capability of the spectral/$hp$ element open-source software framework, Nektar++, to include an implicit discontinuous Galerkin compressible flow solver. The integration in time is carried out by a singly diagonally implicit Runge-Kutta method. The non-linear system arising from the implicit time integration is iteratively solved by the Jacobian-free Newton Krylov (JFNK) method. A favorable feature of the JFNK approach is its extensive use of the explicit operators available from the previous explicit in time implementation. The functionalities of different building blocks of the implicit solver are analyzed from the point of view of software design and placed in appropriate hierarchical levels in the C++ libraries. In the detailed implementation, the contributions of different parts of the solver to computational cost, memory consumption, and programming complexity are also analyzed. A combination of analytical and numerical methods is adopted to simplify the programming complexity in forming the preconditioning matrix. The solver is verified and tested using cases such as manufactured compressible Poiseuille flow, Taylor-Green vortex, turbulent flow over a circular cylinder at $\text{Re}=3900$ and shock wave boundary-layer interaction. The results show that the implicit solver can speed-up the simulations while maintaining good simulation accuracy.

preprint2019arXiv

Naturally curved quadrilateral mesh generation using an adaptive spectral element solver

We describe an adaptive version of a method for generating valid naturally curved quadrilateral meshes. The method uses a guiding field, derived from the concept of a cross field, to create block decompositions of multiply connected two dimensional domains. The a priori curved quadrilateral blocks can be further split into a finer high-order mesh as needed. The guiding field is computed by a Laplace equation solver using a continuous Galerkin or discontinuous Galerkin spectral element formulation. This operation is aided by using $p$-adaptation to achieve faster convergence of the solution with respect to the computational cost. From the guiding field, irregular nodes and separatrices can be accurately located. A first version of the code is implemented in the open source spectral element framework Nektar++ and its dedicated high order mesh generation platform NekMesh.

preprint2019arXiv

Nektar++: enhancing the capability and application of high-fidelity spectral/$hp$ element methods

Nektar++ is an open-source framework that provides a flexible, high-performance and scalable platform for the development of solvers for partial differential equations using the high-order spectral/$hp$ element method. In particular, Nektar++ aims to overcome the complex implementation challenges that are often associated with high-order methods, thereby allowing them to be more readily used in a wide range of application areas. In this paper, we present the algorithmic, implementation and application developments associated with our Nektar++ version 5.0 release. We describe some of the key software and performance developments, including our strategies on parallel I/O, on in situ processing, the use of collective operations for exploiting current and emerging hardware, and interfaces to enable multi-solver coupling. Furthermore, we provide details on a newly developed Python interface that enables a more rapid introduction for new users unfamiliar with spectral/$hp$ element methods, C++ and/or Nektar++. This release also incorporates a number of numerical method developments - in particular: the method of moving frames, which provides an additional approach for the simulation of equations on embedded curvilinear manifolds and domains; a means of handling spatially variable polynomial order; and a novel technique for quasi-3D simulations to permit spatially-varying perturbations to the geometry in the homogeneous direction. Finally, we demonstrate the new application-level features provided in this release, namely: a facility for generating high-order curvilinear meshes called NekMesh; a novel new AcousticSolver for aeroacoustic problems; our development of a 'thick' strip model for the modelling of fluid-structure interaction problems in the context of vortex-induced vibrations. We conclude by commenting some directions for future code development and expansion.

preprint2016arXiv

Influence of localised smooth steps on the instability of a boundary layer

We consider a smooth forward facing step defined by the Gauss error function of height 4-30\% and four times the width of the local boundary layer thickness $δ_{99}$. The boundary layer flow over a smooth forward-facing stepped plate is studied with particular emphasis on stabilisation and destabilisation of the Tollmien-Schlichting (TS) waves and subsequently on transition. The interaction between TS waves at a range of frequencies and a base flow over a single/two forward facing smooth steps is conducted by linear analysis. The results indicate that for a high frequency TS wave, the amplitude of the TS wave is attenuated in the unstable regime of the neutral stability curve corresponding to a flat plate boundary layer. Furthermore, it is observed that two smooth forward facing steps lead to a more acute reduction of the amplitude of the TS wave. When the height of a step is increased to more than 20\% of the local boundary layer thickness for a fixed width parameter, the TS wave is amplified and thereby a destabilisation effect is introduced. Therefore, stabilisation or destabilisation effect of a smooth step is typically dependent on its shape parameters. To validate the results of the linear stability analysis, where a high-frequency TS wave is damped by the forward facing smooth steps direct numerical simulation (DNS) is performed. The results of the DNS correlate favorably with the linear analysis and show that for the investigated high frequency TS wave, the K-type transition process is altered whereas the onset of the H-type transition is postponed. The results of the DNS suggest that for a high-frequency perturbation $\mathcal{F}=150$ and in the absence of other external perturbations, two forward facing steps of height 5\% and 12\% of the boundary layer thickness delayed H-type transition scenario and completely suppresses it for the K-type transition.

preprint2015arXiv

Implicit Large Eddy Simulation of a wingtip vortex at $Re_c = 1.2\cdot 10^6$

In this article we present recent developments in numerical methods for performing a Large Eddy Simulation (LES) of the formation and evolution of a wingtip vortex. The development of these vortices in the near wake, in combination with the large Reynolds numbers present in these cases, make these types of test cases particularly challenging to investigate numerically. We first give an overview of the Spectral Vanishing Viscosity--implicit LES (SVV-iLES) solver that is used to perform the simulations, and highlight techniques that have been adopted to solve various numerical issues that arise when studying such cases. To demonstrate the method's viability, we present results from numerical simulations of flow over a NACA 0012 profile wingtip at $Re_c = 1.2\cdot 10^6$ and compare them against experimental data, which is to date the highest Reynolds number achieved for a LES that has been correlated with experiments for this test case. Our model correlates favorably with experiment, both for the characteristic jetting in the primary vortex and pressure distribution on the wing surface. The proposed method is of general interest for the modeling of transitioning vortex dominated flows over complex geometries.

preprint2014arXiv

An adaptive selective frequency damping method

The selective frequency damping (SFD) method is an alternative to classical Newton's method to obtain unstable steady-state solutions of dynamical systems. However this method has two main limitations: it does not converge for arbitrary control parameters; and when it does converge, the time necessary to reach the steady-state solution may be very long. In this paper we present an adaptive algorithm to address these two issues. We show that by evaluating the dominant eigenvalue of a "partially converged" steady flow, we can select a control coefficient and a filter width that ensure an optimum convergence of the SFD method. We apply this adaptive method to several classical test cases of computational fluid dynamics and we show that a steady-state solution can be obtained without any a priori knowledge of the flow stability properties.

preprint2013arXiv

Encapsulated formulation of the Selective Frequency Damping method

We present an alternative "encapsulated" formulation of the Selective Frequency Damping method for finding unstable equilibria of dynamical systems, which is particularly useful when analysing the stability of fluid flows. The formulation makes use of splitting methods, which means that it can be wrapped around an existing time-stepping code as a "black box". The method is first applied to a scalar problem in order to analyse its stability and highlight the roles of the control coefficient $χ$ and the filter width $Δ$ in the convergence (or not) towards the steady-state. Then the steady-state of the incompressible flow past a two-dimensional cylinder at $Re=100$, obtained with a code which implements the spectral/hp element method, is presented.