No abstract uploaded!

It is known that the maximum diameter for the rupture-risk assessment of the abdominal aortic aneurysm is a generally good method, but not sufficient. Alternative features obtained with computational modeling may provide additional useful criteria. Though computational approaches are noninvasive, they are often time-consuming because of the high computational complexity. In this paper, we present a highly parallel algorithm for the numerical simulation of unsteady blood flows in the patient-specific abdominal aorta. We model the blood flow with the unsteady incompressible Navier-Stokes equations, and solve the discretized system with a highly scalable domain decomposition method. With this approach, the complete flow field can be obtained in less than an hour, instead of days with older methods. We show experimentally that the proposed method offers accurate solutions of the pressure, the velocity and the wall shear stress, and the parallel efficiency is higher than 70% on a parallel computer with more than 1,000 processor cores.

No abstract uploaded!

No abstract uploaded!

An analogue of the well-known $$ \frac{3}{{16}} $$ lower bound for the first eigenvalue of the Laplacian for a congruence hyperbolic surface is proven for a congruence tower associated with any non-elementary subgroup$L$of SL(2,$Z$). The proof in the case that the Hausdorff of the limit set of$L$is bigger than $$ \frac{1}{2} $$ is based on a general result which allows one to transfer such bounds from a combinatorial version to this archimedian setting. In the case that delta is less than $$ \frac{1}{2} $$ we formulate and prove a somewhat weaker version of this phenomenon in terms of poles of the corresponding dynamical zeta function. These “spectral gaps” are then applied to sieving problems on orbits of such groups.

We consider a system coupling the incompressible Navier–Stokes equations to the Vlasov–Fokker–Planck equation. Such a problem arises in the description of particulate flows. We design a numerical scheme to simulate the behavior of the system. This scheme is asymptotic-preserving, thus efficient in both the kinetic and hydrodynamic regimes. It has a numerical stability condition controlled by the non-stiff convection operator, with an implicit treatment of the stiff drag term and the Fokker–Planck operator. Yet, consistent to a standard asymptotic-preserving Fokker–Planck solver or an incompressible Navier–Stokes solver, only the conjugate–gradient method and fast Poisson and Helmholtz solvers are needed. Numerical experiments are presented to demonstrate the accuracy and asymptotic behavior of the scheme, with several interesting applications.

No abstract uploaded!

The diffusive transport of passive tracers or particles can be enhanced by incompressible, turbulent flow fields. Analyzing the effective behavior is a challenging problem with theoretical and practical importance in many areas of science and engineering, ranging from the transport of mass, heat, and pollutants in geophysical flows to sea ice dynamics and turbulent combustion. The long time, large scale behavior of such systems is equivalent to an enhanced diffusion process with an effective diffusivity tensor D*. Two different formulations of integral representations for D* were developed for the case of time-independent fluid velocity fields, involving spectral measures of bounded self-adjoint operators acting on vector fields and scalar fields, respectively. Here, we extend both of these approaches to the case of space-time periodic velocity fields, allowing for chaotic dynamics, providing rigorous integral representations for D* involving spectral measures of unbounded self-adjoint operators.We prove the different formulations are equivalent. Their correspondence follows from a one-to-one isometry between the underlying Hilbert spaces. We also develop a Fourier method for computing D*, which captures the phenomenon of residual diffusion related to Lagrangian chaos of a model flow. This is reflected in the spectral measure by a concentration of mass near the spectral origin.

No abstract uploaded!

We establish local higher integrability estimates for upper gradients of vector-valued parabolic quasi-minimizers in metric measure spaces, satisfying a doubling property and supporting a weak Poincaré inequality.