BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jean-Christophe Nave

Jean-Christophe Nave appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAphysics.flu-dynNumerical Analysisphysics.comp-phmath-phmath.MPComputational Engineering, Finance, and Sciencemath.APnlin.PSphysics.soc-ph

Catalog footprint

What is connected

17works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

17 published item(s)

preprint2024arXiv

Singularity formation of vortex sheets in 2D Euler equations using the characteristic mapping method

The goal of this numerical study is to get insight into singular solutions of the two-dimensional (2D) Euler equations for non-smooth initial data, in particular for vortex sheets. To this end high resolution computations of vortex layers in 2D incompressible Euler flows are performed using the characteristic mapping method (CMM). This semi-Lagrangian method evolves the flow map using the gradient-augmented level set method (GALS). The semi-group structure of the flow map allows its decomposition into sub-maps (each over a finite time interval), and thus the precision can be controlled by choosing appropriate remapping times. Composing the flow map yields exponential resolution in linear time, a unique feature of CMM, and thus fine scale flow structures can be resolved in great detail. Here the roll-up process of vortex layers is studied varying the thickness of the layer showing its impact on the growth of palinstrophy and possible blow up of absolute vorticity. The curvature of the vortex sheet shows a singular-like behavior. The self-similar structure of the vortex core is investigated in the vanishing thickness limit. Conclusions on the non-uniqueness of weak solutions of 2D Euler for non-smooth initial data are drawn and the presence of flow singularities is revealed tracking them in the complex plane.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

High-order FDTD schemes for Maxwell's interface problems with discontinuous coefficients and complex interfaces based on the Correction Function Method

We propose high-order FDTD schemes based on the Correction Function Method (CFM) for Maxwell's interface problems with discontinuous coefficients and complex interfaces. The key idea of the CFM is to model the correction function near an interface to retain the order of a finite difference approximation. For this, we solve a system of PDEs based on the original problem by minimizing an energy functional. The CFM is applied to the standard Yee scheme and a fourth-order FDTD scheme. The proposed CFM-FDTD schemes are verified in 2-D using the transverse magnetic mode (TM$_z$). Numerical examples include scattering of magnetic and non-magnetic dielectric cylinders, and problems with manufactured solutions using various complex interfaces and discontinuous piecewise varying coefficients. Long-time simulations are also performed to provide numerical evidences of the stability of the proposed numerical approach. The proposed CFM-FDTD schemes achieve up to fourth-order convergence in $L^2$-norm and provide approximations devoid of spurious oscillations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Conservative Integrators for Many-body Problems

Conservative symmetric second-order one-step schemes are derived for dynamical systems describing various many-body systems using the Discrete Multiplier Method. This includes conservative schemes for the $n$-species Lotka-Volterra system, the $n$-body problem with radially symmetric potential and the $n$-point vortex models in the plane and on the sphere. In particular, we recover Greenspan-Labudde's conservative schemes for the $n$-body problem. Numerical experiments are shown verifying the conservative property of the schemes and second-order accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Conservative Integrators for Vortex Blob Methods

Conservative symmetric second-order one-step integrators are derived using the Discrete Multiplier Method for a family of vortex-blob models approximating the incompressible Euler's equations on the plane. Conservative properties and second order convergence are proved. A rational function approximation was used to approximate the exponential integral that appears in the Hamiltonian. Numerical experiments are shown to verify the conservative property of these integrators, their second-order accuracy, and as well as the resulting spatial and temporal accuracy of the vortex blob method. Moreover, the derived implicit conservative integrators are shown to be better at preserving conserved quantities than standard higher-order explicit integrators on comparable computation times.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Treatment of complex interfaces for Maxwell's equations with continuous coefficients using the correction function method

We propose a high-order FDTD scheme based on the correction function method (CFM) to treat interfaces with complex geometry without increasing the complexity of the numerical approach for constant coefficients. Correction functions are modeled by a system of PDEs based on Maxwell's equations with interface conditions. To be able to compute approximations of correction functions, a functional that is a square measure of the error associated with the correction functions' system of PDEs is minimized in a divergence-free discrete functional space. Afterward, approximations of correction functions are used to correct a FDTD scheme in the vicinity of an interface where it is needed. We perform a perturbation analysis on the correction functions' system of PDEs. The discrete divergence constraint and the consistency of resulting schemes are studied. Numerical experiments are performed for problems with different geometries of the interface. A second-order convergence is obtained for a second-order FDTD scheme corrected using the CFM. High-order convergence is obtained with a corrected fourth-order FDTD scheme. The discontinuities within solutions are accurately captured without spurious oscillations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A fast-marching algorithm for non-monotonically evolving fronts

The non-monotonic propagation of fronts is considered. When the speed function $F:\mathbb{R}^{n} \times [0,T]\rightarrow \mathbb{R}$ is prescribed, the non-linear advection equation $ϕ_{t}+F|\nabla ϕ|=0$ is a Hamilton-Jacobi equation known as the level-set equation. It is argued that a small enough neighbourhood of the zero-level-set $\mathcal{M}$ of the solution $ϕ: \mathbb{R}^{n} \times [0,T] \rightarrow \mathbb{R}$ is the graph of $ψ:\mathbb{R}^{n} \rightarrow \mathbb{R}$ where $ψ$ solves a Dirichlet problem of the form $H(\vec{u},ψ(\vec{u}),\nabla ψ(\vec{u}))=0$. A fast-marching algorithm is presented where each point is computed using a discretization of such a Dirichlet problem, with no restrictions on the sign of $F$. The output is a directed graph whose vertices evenly sample $\mathcal{M}$. The convergence, consistency and stability of the scheme are addressed. Bounds on the computational complexity are estimated, and experimentally shown to be on par with the Fast Marching Method. Examples are presented where the algorithm is shown to be globally first-order accurate. The complexities and accuracies observed are independent of the monotonicity of the evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Aerodynamic ground effect in fruitfly sized insect takeoff

Aerodynamic ground effect in flapping-wing insect flight is of importance to comparative morphologies and of interest to the micro-air-vehicle (MAV) community. Recent studies, however, show apparently contradictory results of either some significant extra lift or power savings, or zero ground effect. Here we present a numerical study of fruitfly sized insect takeoff with a specific focus on the significance of leg thrust and wing kinematics. Flapping-wing takeoff is studied using numerical modelling and high performance computing. The aerodynamic forces are calculated using a three-dimensional Navier--Stokes solver based on a pseudo-spectral method with volume penalization. It is coupled with a flight dynamics solver that accounts for the body weight, inertia and the leg thrust, while only having two degrees of freedom: the vertical and the longitudinal horizontal displacement. The natural voluntary takeoff of a fruitfly is considered as reference. The parameters of the model are then varied to explore possible effects of interaction between the flapping-wing model and the ground plane. These modified takeoffs include cases with decreased leg thrust parameter, and/or with periodic wing kinematics, constant body pitch angle. The results show that the ground effect during natural voluntary takeoff is negligible. In the modified takeoffs, when the rate of climb is slow, the difference in the aerodynamic forces due to the interaction with the ground is up to 6%. Surprisingly, depending on the kinematics, the difference is either positive or negative, in contrast to the intuition based on the helicopter theory, which suggests positive excess lift. This effect is attributed to unsteady wing-wake interactions. A similar effect is found during hovering.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Fourier penalty method for solving the time-dependent Maxwell's equations in domains with curved boundaries

We present a high order, Fourier penalty method for the Maxwell's equations in the vicinity of perfect electric conductor boundary conditions. The approach relies on extending the smooth non-periodic domain of the equations to a periodic domain by removing the exact boundary conditions and introducing an analytic forcing term in the extended domain. The forcing, or penalty term is chosen to systematically enforce the boundary conditions to high order in the penalty parameter, which then allows for higher order numerical methods. We present an efficient numerical method for constructing the penalty term, and discretize the resulting equations using a Fourier spectral method. We demonstrate convergence orders of up to 3.5 for the one-dimensional Maxwell's equations, and show that the numerical method does not suffer from dispersion (or pollution) errors. We also illustrate the approach in two dimensions and demonstrate convergence orders of 2.5 for transverse magnetic modes and 1.5 for the transverse electric modes. We conclude the paper with numerous test cases in dimensions two and three including waves traveling in an irregular waveguide, and scattering off of a windmill-like geometry.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A low complexity algorithm for non-monotonically evolving fronts

A new algorithm is proposed to describe the propagation of fronts advected in the normal direction with prescribed speed function F. The assumptions on F are that it does not depend on the front itself, but can depend on space and time. Moreover, it can vanish and change sign. To solve this problem the Level-Set Method [Osher, Sethian; 1988] is widely used, and the Generalized Fast Marching Method [Carlini et al.; 2008] has recently been introduced. The novelty of our method is that its overall computational complexity is predicted to be comparable to that of the Fast Marching Method [Sethian; 1996], [Vladimirsky; 2006] in most instances. This latter algorithm is O(N^n log N^n) if the computational domain comprises N^n points. Our strategy is to use it in regions where the speed is bounded away from zero -- and switch to a different formalism when F is approximately 0. To this end, a collection of so-called sideways partial differential equations is introduced. Their solutions locally describe the evolving front and depend on both space and time. The well-posedness of those equations, as well as their geometric properties are addressed. We then propose a convergent and stable discretization of those PDEs. Those alternative representations are used to augment the standard Fast Marching Method. The resulting algorithm is presented together with a thorough discussion of its features. The accuracy of the scheme is tested when F depends on both space and time. Each example yields an O(1/N) global truncation error. We conclude with a discussion of the advantages and limitations of our method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Adaptive gradient-augmented level set method with multiresolution error estimation

A space-time adaptive scheme is presented for solving advection equations in two space dimensions. The gradient-augmented level set method using a semi-Lagrangian formulation with backward time integration is coupled with a point value multiresolution analysis using Hermite interpolation. Thus locally refined dyadic spatial grids are introduced which are efficiently implemented with dynamic quadtree data structures. For adaptive time integration, an embedded Runge-Kutta method is employed. The precision of the new fully adaptive method is analysed and speed up of CPU time and memory compression with respect to the uniform grid discretization are reported.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The multiplier method to construct conservative finite difference schemes for ordinary and partial differential equations

We present the multiplier method of constructing conservative finite difference schemes for ordinary and partial differential equations. Given a system of differential equations possessing conservation laws, our approach is based on discretizing conservation law multipliers and their associated density and flux functions. We show that the proposed discretization is consistent for any order of accuracy when the discrete multiplier has a multiplicative inverse. Moreover, we show that by construction, discrete densities can be exactly conserved. In particular, the multiplier method does not require the system to possess a Hamiltonian or variational structure. Examples, including dissipative problems, are given to illustrate the method. In the case when the inverse of the discrete multiplier becomes singular, consistency of the method is also established for scalar ODEs provided the discrete multiplier and density are zero-compatible.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Sharp-Interface Active Penalty Method for the Incompressible Navier-Stokes Equations

The volume penalty method provides a simple, efficient approach for solving the incompressible Navier-Stokes equations in domains with boundaries or in the presence of moving objects. Despite the simplicity, the method is typically limited to first order spatial accuracy. We demonstrate that one may achieve high order accuracy by introducing an active penalty term. One key difference from other works is that we use a sharp, unregularized mask function. We discuss how to construct the active penalty term, and provide numerical examples, in dimensions one and two. We demonstrate second and third order convergence for the heat equation, and second order convergence for the Navier-Stokes equations. In addition, we show that modifying the penalty term does not significantly alter the time step restriction from that of the conventional penalty method.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Convecting reference frames and invariant numerical models

In the recent paper by Bernardini et al. [J. Comput. Phys. 232 (2013), 1-6] the discrepancy in the performance of finite difference and spectral models for simulations of flows with a preferential direction of propagation was studied. In a simplified investigation carried out using the viscous Burgers equation the authors attributed the poorer numerical results of finite difference models to a violation of Galilean invariance in the discretization and propose to carry out the computations in a reference frame moving with the bulk velocity of the flow. Here we further discuss this problem and relate it to known results on invariant discretization schemes. Non-invariant and invariant finite difference discretizations of Burgers equation are proposed and compared with the discretization using the remedy proposed by Bernardini et al..

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Invariant Discretization Schemes Using Evolution-Projection Techniques

Finite difference discretization schemes preserving a subgroup of the maximal Lie invariance group of the one-dimensional linear heat equation are determined. These invariant schemes are constructed using the invariantization procedure for non-invariant schemes of the heat equation in computational coordinates. We propose a new methodology for handling moving discretization grids which are generally indispensable for invariant numerical schemes. The idea is to use the invariant grid equation, which determines the locations of the grid point at the next time level only for a single integration step and then to project the obtained solution to the regular grid using invariant interpolation schemes. This guarantees that the scheme is invariant and allows one to work on the simpler stationary grids. The discretization errors of the invariant schemes are established and their convergence rates are estimated. Numerical tests are carried out to shed some light on the numerical properties of invariant discretization schemes using the proposed evolution-projection strategy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Correction Function Method for Poisson Problems with Interface Jump Conditions

In this paper we present a method to treat interface jump conditions for constant coefficients Poisson problems that allows the use of standard "black box" solvers, without compromising accuracy. The basic idea of the new approach is similar to the Ghost Fluid Method (GFM). The GFM relies on corrections applied on nodes located across the interface for discretization stencils that straddle the interface. If the corrections are solution-independent, they can be moved to the right-hand-side (RHS) of the equations, producing a problem with the same linear system as if there were no jumps, only with a different RHS. However, achieving high accuracy is very hard (if not impossible) with the "standard" approaches used to compute the GFM correction terms. In this paper we generalize the GFM correction terms to a correction function, defined on a band around the interface. This function is then shown to be characterized as the solution to a PDE, with appropriate boundary conditions. This PDE can, in principle, be solved to any desired order of accuracy. As an example, we apply this new method to devise a 4th order accurate scheme for the constant coefficients Poisson equation with discontinuities in 2D. This scheme is based on (i) the standard 9-point stencil discretization of the Poisson equation, (ii) a representation of the correction function in terms of bicubics, and (iii) a solution of the correction function PDE by a least squares minimization. Several applications of the method are presented to illustrate its robustness dealing with a variety of interface geometries, its capability to capture sharp discontinuities, and its high convergence rate.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A gradient-augmented level set method with an optimally local, coherent advection scheme

The level set approach represents surfaces implicitly, and advects them by evolving a level set function, which is numerically defined on an Eulerian grid. Here we present an approach that augments the level set function values by gradient information, and evolves both quantities in a fully coupled fashion. This maintains the coherence between function values and derivatives, while exploiting the extra information carried by the derivatives. The method is of comparable quality to WENO schemes, but with optimally local stencils (performing updates in time by using information from only a single adjacent grid cell). In addition, structures smaller than the grid size can be located and tracked, and the extra derivative information can be employed to obtain simple and accurate approximations to the curvature. We analyze the accuracy and the stability of the new scheme, and perform benchmark tests.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Self-sustained nonlinear waves in traffic flow

In analogy to gas-dynamical detonation waves, which consist of a shock with an attached exothermic reaction zone, we consider herein nonlinear traveling wave solutions, termed "jamitons," to the hyperbolic ("inviscid") continuum traffic equations. Generic existence criteria are examined in the context of the Lax entropy conditions. Our analysis naturally precludes traveling wave solutions for which the shocks travel downstream more rapidly than individual vehicles. Consistent with recent experimental observations from a periodic roadway (Sugiyama et al., New Journal of Physics, 10, 2008), our numerical calculations show that, under appropriate road conditions, jamitons are attracting solutions, with the time evolution of the system converging towards a jamiton-dominated configuration. Jamitons are characterized by a sharp increase in density over a relatively compact section of the roadway. Applications of our analysis to traffic modeling and control are examined by way of a detailed example.

0 citations0 reviewsNo code linkedOpen access