Simulating the Hydraulic Characteristics of the Lower Yellow River By the Finite-Volume Technique

The finite-volume technique is used to solve the two-dimensional shallow-water equations on unstructured mesh consisting of quadrilateral elements. In this paper the algorithm of the finite-volume method is discussed in detail and particular attention is paid to accurately representing the complex irregular computational domain. The lower Yellow River reach from Huayuankou to Jiahetan is a typical meandering river. The generation of the computational mesh, which is used to simulate the flood, is affected by the distribution of water works in the river channel. The spatial information about the two Yellow River levee, the protecting dykes, and those roads that are obviously higher than the ground, need to be used to generate the computational mesh. As a result these dykes and roads locate the element interfaces of the computational mesh. In the model the finite-volume method is used to solve the shallow-wave equations, and the Osher scheme of the empirical function is used to calculate the flux through the interface between the neighbouring elements. The finite-volume method has the advantage of using computational domain with complex geometry, and the Osher scheme is a method based on characteristic theory and is a monotone upwind numerical scheme with high resolution. The flood event with peak discharge of 15 300 m(3)/s, occurring in the period from 30 July to 10 August 1982, is simulated. The estimated result indicates that the simulation method is good for routing the flood in a region with complex geometry. Copyright (C) 2002 John Wiley Sons, Ltd

Computational simulation of an unmanned air vehicle impacting water

Computational simulations were performed to study the splashdown of an unmanned air vehicle (UAV) falling nose-first into seawater from various heights. Solutions were generated with a time-accurate finite-volume method based on the unsteady compressible ensemble averaged Navier-Stokes equations for the air and the unsteady incompressible ensemble averaged Navier-Stokes equations for the seawater. The volume of fluid model was used to track the air-water interface and a deforming mesh algorithm was used to move the UAV through the computational domain. Computed pressure histories at four key locations on the UAV forebody were compared with experimentally measured values to validate this study. The computational simulations were shown to have accurately predicted the magnitude and character of the pressure histories, but with some discrepancies in the behavior of the pressure within the UAV inlet aperture. Results are presented for various drop heights, which simulated a range of impact velocities. Modifications were also made to the UAV geometry to examine the effect of deflecting the upper inlet lip downward. Deflection angles of 30 deg and 20 deg were analyzed for a drop height condition of 35 ft with results showing a significant decrease in impact force and pressure within the inlet

Computationally Efficient Forward Operator for Photoacoustic Tomography Based on Coordinate Transformations

IEEE Photoacoustic tomography (PAT) is an imaging modality that utilizes the photoacoustic effect. In PAT, a photoacoustic image is computed from measured data by modeling ultrasound propagation in the imaged domain and solving an inverse problem utilizing a discrete forward operator. However, in realistic measurement geometries with several ultrasound transducers and relatively large imaging volume, an explicit formation and use of the forward operator can be computationally prohibitively expensive. In this work, we propose a transformation based approach for efficient modeling of photoacoustic signals and reconstruction of photoacoustic images. In the approach, the forward operator is constructed for a reference ultrasound transducer and expanded into a general measurement geometry using transformations that map the formulated forward operator in local coordinates to the global coordinates of the measurement geometry. The inverse problem is solved using a Bayesian framework. The approach is evaluated with numerical simulations and experimental data. The results show that the proposed approach produces accurate three-dimensional photoacoustic images with a significantly reduced computational cost both in memory requirements and in time. In the studied cases, depending on the computational factors such as discretization, over 30-fold reduction in memory consumption and was achieved without a reduction in image quality compared to a conventional approach

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

Mechanistic and pathological study of the genesis, growth, and rupture of abdominal aortic aneurysms

Postprint (published version

Phase-field boundary conditions for the voxel finite cell method: surface-free stress analysis of CT-based bone structures

The voxel finite cell method employs unfitted finite element meshes and voxel quadrature rules to seamlessly transfer CT data into patient-specific bone discretizations. The method, however, still requires the explicit parametrization of boundary surfaces to impose traction and displacement boundary conditions, which constitutes a potential roadblock to automation. We explore a phase-field based formulation for imposing traction and displacement constraints in a diffuse sense. Its essential component is a diffuse geometry model generated from metastable phase-field solutions of the Allen-Cahn problem that assumes the imaging data as initial condition. Phase-field approximations of the boundary and its gradient are then employed to transfer all boundary terms in the variational formulation into volumetric terms. We show that in the context of the voxel finite cell method, diffuse boundary conditions achieve the same accuracy as boundary conditions defined over explicit sharp surfaces, if the inherent length scales, i.e., the interface width of the phase-field, the voxel spacing and the mesh size, are properly related. We demonstrate the flexibility of the new method by analyzing stresses in a human femur and a vertebral body