Source author record

Kamran Mohseni

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

13 published item(s)

preprint2020arXiv

On the maintenance of an attached leading-edge vortex via model bird alula

Researchers have hypothesized that the post-stall lift benefit of bird's alular feathers, or alula, stems from the maintenance of an attached leading-edge vortex (LEV) over their thin-profiled, outer hand-wing. Here, we investigate the connection between the alula and LEV attachment via flow measurements in a wind tunnel. We show that a model alula, whose wetted area is 1% that of the wing, stabilizes a recirculatory aft-tilted LEV on a steadily-translating unswept wing inclined at post-stall incidences. The attached vortex is the result of the alula's ability to smoothly merge otherwise separate leading- and side-edge vortical flows. We identify two key processes that facilitate this merging: i) the steering of spanwise vorticity generated at the wing's leading edge back to the wing plane and ii) an aft-located wall-jet of high-magnitude root-to-tip spanwise flow (>80% that of the freestream velocity). The latter feature induces LEV roll-up while the former feature tilts LEV vorticity aft and evacuates this flow toward the wing tip via an outboard vorticity flux. We identify the alula's streamwise position (relative to the leading-edge of the thin wing) as important for vortex steering and the alula's cant angle as important for high-magnitude spanwise flow generation. These findings advance our understanding of the likely ways bird's leverage LEVs to augment slow flight.

preprint2016arXiv

Structure-Preserving Model Reduction of Forced Hamiltonian Systems

This paper reports a development in the proper symplectic decomposition (PSD) for model reduction of forced Hamiltonian systems. As an analogy to the proper orthogonal decomposition (POD), PSD is designed to build a symplectic subspace to fit empirical data. Our aim is two-fold. First, to achieve computational savings for large-scale Hamiltonian systems with external forces. Second, to simultaneously preserve the symplectic structure and the forced structure of the original system. We first reformulate d'Alembert's principle in the Hamiltonian form. Corresponding to the integral and local forms of d'Alembert's principle, we propose two different structure-preserving model reduction approaches to reconstruct low-dimensional systems, based on the variational principle and on the structure-preserving projection, respectively. These two approaches are proven to yield the same reduced system. Moreover, by incorporating the vector field into the data ensemble, we provided several algorithms for energy preservation. In a special case when the external force is described by the Rayleigh dissipative function, the proposed method automatically preserves the dissipativity, boundedness, and stability of the original system. The stability, accuracy, and efficiency of the proposed method are illustrated through numerical simulations of a dissipative wave equation.

preprint2015arXiv

A unified slip boundary condition for flow over a surface

Interface between two phases of matter are ubiquitous in nature and technology. Determining the correct velocity condition at an interface is essential for understanding and designing of flows over a surface. We demonstrate that both the widely used no-slip and the Navier and Maxwell slip boundary conditions do not capture the complete physics associated with complex problems, such as spreading of liquids or corner flows. Hence, we present a unified boundary condition that is applicable to a wide-range of flow problems.

preprint2015arXiv

Symplectic Model Reduction of Hamiltonian Systems

In this paper, a symplectic model reduction technique, proper symplectic decomposition (PSD) with symplectic Galerkin projection, is proposed to save the computational cost for the simplification of large-scale Hamiltonian systems while preserving the symplectic structure. As an analogy to the classical proper orthogonal decomposition (POD)-Galerkin approach, PSD is designed to build a symplectic subspace to fit empirical data, while the symplectic Galerkin projection constructs a reduced Hamiltonian system on the symplectic subspace. For practical use, we introduce three algorithms for PSD, which are based upon: the cotangent lift, complex singular value decomposition, and nonlinear programming. The proposed technique has been proven to preserve system energy and stability. Moreover, PSD can be combined with the discrete empirical interpolation method to reduce the computational cost for nonlinear Hamiltonian systems. Owing to these properties, the proposed technique is better suited than the classical POD-Galerkin approach for model reduction of Hamiltonian systems, especially when long-time integration is required. The stability, accuracy, and efficiency of the proposed technique are illustrated through numerical simulations of linear and nonlinear wave equations.

preprint2014arXiv

An Online Manifold Learning Approach for Model Reduction of Dynamical Systems

This article discusses a newly developed online manifold learning method, subspace iteration using reduced models (SIRM), for the dimensionality reduction of dynamical systems. This method may be viewed as subspace iteration combined with a model reduction procedure. Specifically, starting with a test solution, the method solves a reduced model to obtain a more precise solution, and it repeats this process until sufficient accuracy is achieved. The reduced model is obtained by projecting the full model onto a subspace that is spanned by the dominant modes of an extended data ensemble. The learning procedure is computed in the online stage, as opposed to being computed offline, which is used in many projection-based model reduction techniques that require prior calculations or experiments. After providing an error bound of the classical POD-Galerkin method in terms of the projection error and the initial condition error, we prove that the sequence of approximate solutions converge to the actual solution of the original system as long as the vector field of the full model is locally Lipschitz on an open set that contains the solution trajectory. Good accuracy of the proposed method has been demonstrated in two numerical examples, from a linear advection-diffusion equation to a nonlinear Burgers equation. In order to save computational cost, the SIRM method is extended to a local model reduction approach by partitioning the entire time domain into several subintervals and obtaining a series of local reduced models of much lower dimensionality. The accuracy and efficiency of the local SIRM are shown through the numerical simulation of the Navier--Stokes equation in a lid-driven cavity flow problem.

preprint2014arXiv

Nonlinear Model Reduction via an Adaptive Weighting of Snapshots

In this paper, we propose a new approach to model reduction of parameterized partial differential equations (PDEs) based on the concept of adaptive reduced bases. The presented approach is particularly suited for large-scale nonlinear systems characterized by parameter variations. Instead of using a global basis to construct a global reduced model, the proposed method approximates the original system by multiple lower-dimensional subspaces. Each localized reduced basis is generated by the SVD of a weighted snapshot ensemble; here, each weighting coefficient is a function of the input parameter. Compared with a global model reduction method, such as the classical POD, the adaptive model reduction method could yield a more accurate solution with a fixed subspace dimension. Moreover, we combine the adaptive reduced model with the chord iteration to solve elliptic PDEs in a computationally efficient fashion. The potential of the method for achieving large speedups, while maintaining good accuracy, is demonstrated for both elliptic and parabolic PDEs in a few numerical examples.

preprint2013arXiv

Effect of slip on circulation inside a droplet

Internal recirculation in a moving droplet plays an important role in several droplet-based microfluidic devices as it enhances mixing, chemical reaction and heat transfer. The occurrence of fluid slip at the wall, which becomes prominent at high shear rates and lower length scales, results in a significant change in droplet circulation. Using molecular dynamics (MD) simulations, the presence of circulation in droplets is demonstrated and quantified. Circulation is shown to vary inversely with slip length, which is a measure of interface wettability. A simple circulation model is established that captures the effect of slip on droplet circulation. Scaling parameters for circulation and slip length are identified from the circulation model which leads to the collapse of data for droplets with varying aspect ratio (AR) and slip length. The model is validated using continuum and MD simulations and is shown to be accurate for droplets with high AR.

preprint2013arXiv

Observations on the flow structures and transport in a warm-core ring in the Gulf of Mexico

This study presents several new observations from the study of a warm-core ring (WCR) in the Gulf of Mexico based on the ECCO2 global ocean simulation. Using Lagrangian coherent structures (LCS) techniques to investigate this flow reveals a pattern of transversely intersecting LCS in the mixed layer of the WCR which experiences consistent stretching behavior over a large region of space and time. A detailed analysis of this flow region leads to an analytical model velocity field which captures the essential elements that generate the transversely intersecting LCS. The model parameters are determined from the WCR and the resulting LCS show excellent agreement with those observed in the WCR. The three-dimensional transport behavior which creates these structures relies on the small radial outflow which is present in the mixed layer and is not seen below the pycnocline, leading to a sharp change in the character of the LCS at the bottom of the mixed layer. The flow behavior revealed by the LCS limits fluid exchange between the WCR and the surrounding ocean, contributing to the long life of WCRs. Further study of these structures and their associated transport behavior may lead to further insights into the development and persistence of such geophysical vortices as well as their transport behavior.

preprint2012arXiv

A 3D fast algorithm for computing Lagrangian coherent structures via ridge tracking

Lagrangian coherent structures (LCS) in fluid flows appear as co-dimension one ridges of the finite time Lyapunov exponent (FTLE) field. In three- dimensions this means two-dimensional ridges. A fast algorithm is presented here to locate and extract such ridge surfaces while avoiding unnecessary computations away from the LCS. This algorithm reduces the order of the computational complexity from O(1/dx^3) to about O(1/dx^2) by eliminating computations over most of the three dimensional domain and computing the FTLE only near the two-dimensional ridge surfaces. The algorithm is grid based and proofs of error bounds for ridge locations are included. The algorithm performance and error bounds are verified in several examples. The algorithm offers significant advantages in computational cost as well as later data analysis.

preprint2012arXiv

Analysis of Boundary Slip in a Flow with an Oscillating Wall

Molecular dynamic (MD) simulation is used to study slip at the fluid-solid boundary in an unsteady flow based on the Stokes second problem. An increase in slip is observed in comparison to the steady flow for shear rates below the critical shear rate of the corresponding steady flow. This increased slip is attributed to fluid inertial forces not represented in a steady flow. An unsteady mathematical model for slip is established, which estimates the increment in slip at the boundary. The model shows that slip is also dependent on acceleration in addition to the shear rate of fluid at the wall. By writing acceleration in terms of shear rate, it is shown that slip at the wall in unsteady flows is governed by the gradient of shear rate and shear rate of the fluid. Non-dimensionalizing the model gives a universal curve which can be used to find the slip boundary condition at the fluid-solid interface based on the information of shear rate and gradient of shear rate of the fluid. A governing non-dimensional number, defined as the ratio of phase speed to speed of sound, is identified to help in explaining the mechanism responsible for the transition of slip boundary condition from finite to a perfect slip and determining when this would occur. Phase lag in fluid velocity relative to wall is observed. The lag increases with decreasing time period of wall oscillation and increasing wall hydrophobicity. The phenomenon of hysteresis is seen when looking into the variation of slip velocity as a function of wall velocity and slip velocity as a function of fluid shear rate. The cause for hysteresis is attributed to the unsteady inertial forces of the fluid.

preprint2010arXiv

A regularization for the transport equations using spatial-averaging

This paper examines an averaging technique applied to the transport equations as an alternative to vanishing viscosity. Such techniques have been shown to be valid shock-regularizations of the Burgers equation and the Euler equations, but has yet to be applied to the similar transport equations. However, for this system, the classical notion of weak solutions is not always sufficient thus a more general notion of a distribution solution containing Dirac-delta functions must be introduced. Moreover, the distribution solution to the Riemann problem is known to be the weak-$*$ limit of the viscous perturbed transport equations as viscosity vanishes. In comparison to the classical method of vanishing viscosity, the Riemann problem is examined for the averaged transport equation and it is shown that the same delta-shock distribution solution is captured as filtering vanishes. Both mathematical and physical motivation are provided for the averaging techniques considered including an existence and uniqueness result of solutions within the class of smooth initial conditions.

preprint2010arXiv

An arbitrary Lagrangian-Eulerian formulation for the numerical simulation of flow patterns generated by the hydromedusa \textit{Aequorea victoria}

A new geometrically conservative arbitrary Lagrangian-Eulerian (ALE) formulation is presented for the moving boundary problems in the swirl-free cylindrical coordinates. The governing equations are multiplied with the radial distance and integrated over arbitrary moving Lagrangian-Eulerian quadrilateral elements. Therefore, the continuity and the geometric conservation equations take very simple form similar to those of the Cartesian coordinates. The continuity equation is satisfied exactly within each element and a special attention is given to satisfy the geometric conservation law (GCL) at the discrete level. The equation of motion of a deforming body is solved in addition to the Navier-Stokes equations in a fully-coupled form. The mesh deformation is achieved by solving the linear elasticity equation at each time level while avoiding remeshing in order to enhance numerical robustness. The resulting algebraic linear systems are solved using an ILU(k) preconditioned GMRES method provided by the PETSc library. The present ALE method is validated for the steady and oscillatory flow around a sphere in a cylindrical tube and applied to the investigation of the flow patterns around a free-swimming hydromedusa \textit{Aequorea victoria} (crystal jellyfish). The calculations for the hydromedusa indicate the shed of the opposite signed vortex rings very close to each other and the formation of large induced velocities along the line of interaction while the ring vortices moving away from the hydromedusa. In addition, the propulsion efficiency of the free-swimming hydromedusa is computed and its value is compared with values from the literature for several other species. The fluid dynamics video presented here shows the time variation of the instantaneous three-dimensional vorticity isosurfaces around a free-swimming hydromedusa \textit{Aequorea victoria}.

preprint2010arXiv

Observable Divergence Theorem: Evolution Equations for Inviscid Regularization of Shocks and Turbulence

The divergence theorem of Gauss plays a central role in the derivation of the governing differential equations in fluid dynamics, electrodynamics, gravitational fields, and optics. One is often interested in an evolution equation for the large scale quantities without resolving the details of the small scales. As a result, there has been a significant effort in developing time-averaged and spatially-filtered equations for large scale dynamics from the fully resolved governing differential equations. One should realize that by starting from these fully-resolved equations (e.g. the Euler or Navier-Stokes equations) to derive an averaged evolution equation one has already taken the limit of the wave-numbers approaching infinity with no regards to our observational abilities at such a limit. As a result, obtaining the evolution equations for large scale quantities (low wave-numbers) by an averaging or filtering process is done after the fact. This could explain many of the theoretical and computational difficulties with the Euler or Navier-Stokes equations. Here, a rather different approach is proposed. The averaging process in implemented before the derivation of the differential form of the transport equations. A new observable divergence concept is defined based on fluxes calculated from observable quantities at a desired averaging scale, $α$. An observable divergence theorem is then proved and applied in the derivation of the observable and regularized transport equations. We further show that the application of the observable divergence theorem to incompressible flows results in a formal derivation of the inviscid Leray turbulence model first proposed in 1934. It is argued that such a methodology in deriving fluid evolution equations removes many of the theoretical and computational difficulties in multi-scale problems such as turbulence and shocks.