Rigid Wall (rigid + wall)

Distribution by Scientific Domains


Selected Abstracts


Optimal transportation meshfree approximation schemes for fluid and plastic flows

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 12 2010
B. Li
Abstract We develop an optimal transportation meshfree (OTM) method for simulating general solid and fluid flows, including fluid,structure interaction. The method combines concepts from optimal transportation theory with material-point sampling and max-ent meshfree interpolation. The proposed OTM method generalizes the Benamou,Brenier differential formulation of optimal mass transportation problems to problems including arbitrary geometries and constitutive behavior. The OTM method enforces mass transport and essential boundary conditions exactly and is free from tension instabilities. The OTM method exactly conserves linear and angular momentum and its convergence characteristics are verified in standard benchmark problems. We illustrate the range and scope of the method by means of two examples of application: the bouncing of a gas-filled balloon off a rigid wall; and the classical Taylor-anvil benchmark test extended to the hypervelocity range. Copyright © 2010 John Wiley & Sons, Ltd. [source]


Numerical simulation of a single bubble by compressible two-phase fluids

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 6 2010
Siegfried Müller
Abstract The present work deals with the numerical investigation of a collapsing bubble in a liquid,gas fluid, which is modeled as a single compressible medium. The medium is characterized by the stiffened gas law using different material parameters for the two phases. For the discretization of the stiffened gas model, the approach of Saurel and Abgrall is employed where the flow equations, here the Euler equations, for the conserved quantities are approximated by a finite volume scheme, and an upwind discretization is used for the non-conservative transport equations of the pressure law coefficients. The original first-order discretization is extended to higher order applying second-order ENO reconstruction to the primitive variables. The derivation of the non-conservative upwind discretization for the phase indicator, here the gas fraction, is presented for arbitrary unstructured grids. The efficiency of the numerical scheme is significantly improved by employing local grid adaptation. For this purpose, multiscale-based grid adaptation is used in combination with a multilevel time stepping strategy to avoid small time steps for coarse cells. The resulting numerical scheme is then applied to the numerical investigation of the 2-D axisymmetric collapse of a gas bubble in a free flow field and near to a rigid wall. The numerical investigation predicts physical features such as bubble collapse, bubble splitting and the formation of a liquid jet that can be observed in experiments with laser-induced cavitation bubbles. Opposite to the experiments, the computations reveal insight to the state inside the bubble clearly indicating that these features are caused by the acceleration of the gas due to shock wave focusing and reflection as well as wave interaction processes. While incompressible models have been used to provide useful predictions on the change of the bubble shape of a collapsing bubble near a solid boundary, we wish to study the effects of shock wave emissions into the ambient liquid on the bubble collapse, a phenomenon that may not be captured using an incompressible fluid model. Copyright © 2009 John Wiley & Sons, Ltd. [source]


Existence of a weak solution to the Navier,Stokes equation in a general time-varying domain by the Rothe method

MATHEMATICAL METHODS IN THE APPLIED SCIENCES, Issue 6 2009
í Neustupa
Abstract We assume that ,t is a domain in ,3, arbitrarily (but continuously) varying for 0,t,T. We impose no conditions on smoothness or shape of ,t. We prove the global in time existence of a weak solution of the Navier,Stokes equation with Dirichlet's homogeneous or inhomogeneous boundary condition in Q[0,,T) := {(x,,t);0,t,T, x,,t}. The solution satisfies the energy-type inequality and is weakly continuous in dependence of time in a certain sense. As particular examples, we consider flows around rotating bodies and around a body striking a rigid wall. Copyright © 2008 John Wiley & Sons, Ltd. [source]


Comparison of two wave element methods for the Helmholtz problem

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 1 2009
T. Huttunen
Abstract In comparison with low-order finite element methods (FEMs), the use of oscillatory basis functions has been shown to reduce the computational complexity associated with the numerical approximation of Helmholtz problems at high wave numbers. We compare two different wave element methods for the 2D Helmholtz problems. The methods chosen for this study are the partition of unity FEM (PUFEM) and the ultra-weak variational formulation (UWVF). In both methods, the local approximation of wave field is computed using a set of plane waves for constructing the basis functions. However, the methods are based on different variational formulations; the PUFEM basis also includes a polynomial component, whereas the UWVF basis consists purely of plane waves. As model problems we investigate propagating and evanescent wave modes in a duct with rigid walls and singular eigenmodes in an L-shaped domain. Results show a good performance of both methods for the modes in the duct, but only a satisfactory accuracy was obtained in the case of the singular field. On the other hand, both the methods can suffer from the ill-conditioning of the resulting matrix system. Copyright © 2008 John Wiley & Sons, Ltd. [source]