Graph-like asymptotics for the Dirichlet Laplacian in connected tubular domains

We consider the Dirichlet Laplacian in a waveguide of uniform width and infinite length which is ideally divided into three parts: a "vertex region", compactly supported and with non zero curvature, and two "edge regions" which are semi-infinite straight strips. We make the waveguide collapse onto a graph by squeezing the edge regions to half-lines and the vertex region to a point. In a setting in which the ratio between the width of the waveguide and the longitudinal extension of the vertex region goes to zero, we prove the convergence of the operator to a selfadjoint realization of the Laplacian on a two edged graph. In the limit operator, the boundary conditions in the vertex depend on the spectral properties of an effective one dimensional Hamiltonian associated to the vertex region.Comment: Major revision. Reviewed introduction. Changes in Th. 1, Th. 2, and Th. 3. Updated references. 23 page

Stability of Attached Transonic Shocks in Steady Potential Flow past Three-Dimensional Wedges

We develop a new approach and employ it to establish the global existence and nonlinear structural stability of attached weak transonic shocks in steady potential flow past three-dimensional wedges; in particular, the restriction that the perturbation is away from the wedge edge in the previous results is removed. One of the key ingredients is to identify a "good" direction of the boundary operator of a boundary condition of the shock along the wedge edge, based on the non-obliqueness of the boundary condition for the weak shock on the edge. With the identification of this direction, an additional boundary condition on the wedge edge can be assigned to make sure that the shock is attached on the edge and linearly stable under small perturbation. Based on the linear stability, we introduce an iteration scheme and prove that there exists a unique fixed point of the iteration scheme, which leads to the global existence and nonlinear structural stability of the attached weak transonic shock. This approach is based on neither the hodograph transformation nor the spectrum analysis, and should be useful for other problems with similar difficulties.Comment: 28 Pages; 2 figure

Circular edge singularities for the Laplace equation and the elasticity system in 3-D domains

International audienceAsymptotics of solutions to the Laplace equation with Neumann or Dirichlet conditions in the vicinity of a circular singular edge in a three-dimensional domain are derived and provided in an explicit form. These asymptotic solutions are represented by a family of eigen-functions with their shadows, and the associated edge flux intensity functions (EFIFs), which are functions along the circular edge. We provide explicit formulas for a penny-shaped crack for an axisymmetric case as well as a case in which the loading is non-axisymmetric. Explicit formulas for other singular circular edges such as a circumferential crack and an external crack are also derived. The mathematical machinery developed in the framework of the Laplace operator is extended to derive the asymptotic solution (three-component displacement vector) for the elasticity system in the vicinity of a circular edge in a three-dimensional domain. As a particular case we present explicitly the series expansion for a traction free or clamped penny-shaped crack in an axisymmetric or a non-axisymmetric situation. The precise representation of the asymptotic series is required for constructing benchmark problems with analytical solutions against which numerical methods can be assessed, and to develop new extraction techniques for the edge flux/intensity functions which are of practical engineering importance in predicting crack propagation. <br

Investigation of advanced counterrotation blade configuration concepts for high speed turboprop systems, task 1: Ducted propfan analysis

The time-dependent three-dimensional Euler equations of gas dynamics were solved numerically to study the steady compressible transonic flow about ducted propfan propulsion systems. Aerodynamic calculations were based on a four-stage Runge-Kutta time-marching finite volume solution technique with added numerical dissipation. An implicit residual smoothing operator was used to aid convergence. Two calculation grids were employed in this study. The first grid utilized an H-type mesh network with a branch cut opening to represent the axisymmetric cowl. The second grid utilized a multiple-block mesh system with a C-type grid about the cowl. The individual blocks were numerically coupled in the Euler solver. Grid systems were generated by a combined algebraic/elliptic algortihm developed specifically for ducted propfans. Numerical calculations were initially performed for unducted propfans to verify the accuracy of the three-dimensional Euler formulation. The Euler analyses were then applied for the calculation of ducted propfan flows, and predicted results were compared with experimental data for two cases. The three-dimensional Euler analyses displayed exceptional accuracy, although certain parameters were observed to be very sensitive to geometric deflections. Both solution schemes were found to be very robust and demonstrated nearly equal efficiency and accuracy, although it was observed that the multi-block C-grid formulation provided somewhat better resolution of the cowl leading edge region

Extracting generalized edge flux intensity functions with the quasidual function method along circular 3-D edges

International audienceExplicit asymptotic series describing solutions to the Laplace equation in the vicinity of a circular edge in a three-dimensional domain was recently provided in Yosibash et al, Int. Jour. Fracture, 168 (2011), pp. 31-52. Utilizing it, we extend the quasidual function method (QDFM) for extracting the generalized edge flux intensity functions (GEFIFs) along circular singular edges in the cases of axisymmetric and non-axisymmetric data. This accurate and efficient method provides a functional approximation of the GEFIFs along the circular edge whose order is adaptively increased so to approximate the exact GEFIFs. It is implemented as a post-solution operation in conjunction with the p-version of the finite element method. The mathematical analysis of the QDFM is provided, followed by numerical investigations, demonstrating the efficiency, robustness and high accuracy of the proposed quasi-dual function method. The mathematical machinery developed in the framework of the Laplace operator is important to realize its possible extension for the elasticity system

Computations of supersonic viscous flow over a finite-width plate

Finite difference methods were applied to solve the parabolic Navier-Stokes equations for the flow over a finite width plate at 0 deg and 10 deg angles of attack. The methods were developed on the basis of the operator factorization concept resulting in the split of a three dimensional equation into successive two dimensional equations. Backward and centered implicit factorization schemes, were used and their results were compared. Available numerical solutions and experimental data obtained at low Reynolds number conditions were also used for comparison. The backward implicit method provides a more successful solution, which ranges from the merged layer to the strong interaction regimes. Detailed structures were revealed of the shear layer around and behind the side edge