A conservative coupling algorithm between a compressible flow and a rigid body using an Embedded Boundary method

This paper deals with a new solid-fluid coupling algorithm between a rigid body and an unsteady compressible fluid flow, using an Embedded Boundary method. The coupling with a rigid body is a first step towards the coupling with a Discrete Element method. The flow is computed using a Finite Volume approach on a Cartesian grid. The expression of numerical fluxes does not affect the general coupling algorithm and we use a one-step high-order scheme proposed by Daru and Tenaud [Daru V,Tenaud C., J. Comput. Phys. 2004]. The Embedded Boundary method is used to integrate the presence of a solid boundary in the fluid. The coupling algorithm is totally explicit and ensures exact mass conservation and a balance of momentum and energy between the fluid and the solid. It is shown that the scheme preserves uniform movement of both fluid and solid and introduces no numerical boundary roughness. The effciency of the method is demonstrated on challenging one- and two-dimensional benchmarks

A Moving Boundary Flux Stabilization Method for Cartesian Cut-Cell Grids using Directional Operator Splitting

An explicit moving boundary method for the numerical solution of time-dependent hyperbolic conservation laws on grids produced by the intersection of complex geometries with a regular Cartesian grid is presented. As it employs directional operator splitting, implementation of the scheme is rather straightforward. Extending the method for static walls from Klein et al., Phil. Trans. Roy. Soc., A367, no. 1907, 4559-4575 (2009), the scheme calculates fluxes needed for a conservative update of the near-wall cut-cells as linear combinations of standard fluxes from a one-dimensional extended stencil. Here the standard fluxes are those obtained without regard to the small sub-cell problem, and the linear combination weights involve detailed information regarding the cut-cell geometry. This linear combination of standard fluxes stabilizes the updates such that the time-step yielding marginal stability for arbitrarily small cut-cells is of the same order as that for regular cells. Moreover, it renders the approach compatible with a wide range of existing numerical flux-approximation methods. The scheme is extended here to time dependent rigid boundaries by reformulating the linear combination weights of the stabilizing flux stencil to account for the time dependence of cut-cell volume and interface area fractions. The two-dimensional tests discussed include advection in a channel oriented at an oblique angle to the Cartesian computational mesh, cylinders with circular and triangular cross-section passing through a stationary shock wave, a piston moving through an open-ended shock tube, and the flow around an oscillating NACA 0012 aerofoil profile.Comment: 30 pages, 27 figures, 3 table

New ideas for teaching natural resource management: From the long-term realities of national forest management

Research and study of 90 years of managing multiple uses on national forests has revealed three new ideas or understandings about the nature of forest management (Fedkiw 1997a). The first idea is a new definition that describes the task of forest management and the role of forest managers. The second emphasizes the critical, continuous role of the learning experience that accompanies resource management and its relationship to both the adaptive and holistic ecological approaches to resource management. The third establishes that forest management has been on a pathway toward a holistic ecological approach from the beginning of American forestry. It also describes how forest management advanced, and continues to advance, incrementally and adaptively on that pathway in response to intensifying and diversifying uses and services; improving experience, technology, and science; changing markets and social preferences, and Nature\u27s unexpected responses to use and management and her own random vagaries. These ideas have a large potential fro improving the knowledge, teaching, communication, and progress of forest management in the classroom, in the field, and with the general public and its interest groups. To be effective, however, these ideas must be communicated, discussed, debated, researched, tested, refined, and written about, not only among resource professionals but also with students, interest groups, stakeholders, landowners, policymakers, and the public-at-large. New ideas tend to roll off like water off a duck\u27s back unless they are communicated, discussed, and debated; highlighted in their newness; packaged in a familiar context, and presented in a user/audience friendly way with graphic images (Perry 1993)

A survey of partial differential equations in geometric design

YesComputer aided geometric design is an area where the improvement of surface generation techniques is an everlasting demand since faster and more accurate geometric models are required. Traditional methods for generating surfaces were initially mainly based upon interpolation algorithms. Recently, partial differential equations (PDE) were introduced as a valuable tool for geometric modelling since they offer a number of features from which these areas can benefit. This work summarises the uses given to PDE surfaces as a surface generation technique togethe

Characterization of reaction kinetics in a porous electrode

A continuum-model approach, analogous to porous electrode theory, was applied to a thin-layer cell of rectangular and cylindrical geometry. A reversible redox couple is assumed, and the local reaction current density is related to the potential through the formula of Hubbard and Anson for a uniformily accessible thin-layer cell. The placement of the reference electrode is also accounted for in the analysis. Primary emphasis is placed on the effect of the solution-phase ohmic potential drop on the voltammogram characteristics. Correlation equations for the peak-potential displacement from E(sup 0 prime) and the peak current are presented in terms of two dimensionless parameters

MLPG_R method for modelling 2D flows of two immiscible fluids

This is a first attempt to develop the Meshless Local Petrov-Galerkin method with Rankine source solution (MLPG_R method) to simulate multiphase flows. In this paper, we do not only further develop the MLPG_R method to model two-phase flows but also propose two new techniques to tackle the associated challenges. The first technique is to form an equation for pressure on the explicitly identified interface between different phases by considering the continuity of the pressure and the discontinuity of the pressure gradient (i.e. the ratio of pressure gradient to fluid density), the latter reflecting the fact that the normal velocity is continuous across the interface. The second technique is about solving the algebraic equation for pressure, which gives reasonable solution not only for the cases with low density ratio but also for the cases with very high density ratio, such as more than 1000. The numerical tests show that the results of the newly developed two-phase MLPG_R method agree well with analytical solutions and experimental data in the cases studied. The numerical results also demonstrate that the newly developed method has a second-order convergent rate in the cases for sloshing motion with small amplitudes