			Home About us Contact

Momentum Equations (momentum + equation)

Distribution by Scientific Domains

Engineering	74%
Polymers and Materials Science	3%
5 Other Domains	23%

Selected Abstracts

Nonlinear Smoluchowski velocity for electroosmosis of Power-law fluids over a surface with arbitrary zeta potentials

ELECTROPHORESIS, Issue 5 2010
Cunlu Zhao
Abstract Electroosmotic flow of Power-law fluids over a surface with arbitrary zeta potentials is analyzed. The governing equations including the nonlinear Poisson,Boltzmann equation, the Cauchy momentum equation and the continuity equation are solved to seek exact solutions for the electroosmotic velocity, shear stress, and dynamic viscosity distributions inside the electric double layer. Specifically, an expression for the general Smoluchowski velocity is obtained for electroosmosis of Power-law fluids in a fashion similar to the classic Smoluchowski velocity for Newtonian fluids. The existing Smoluchowski slip velocities under two special cases, (i) for Newtonian fluids with arbitrary zeta potentials and (ii) for Power-law fluids with small zeta potentials, can be recovered from our derived formula. It is interesting to note that the general Smoluchowski velocity for non-Newtonian Power-law fluids is a nonlinear function of the electric field strength and surface zeta potentials; this is due to the coupling electrostatics and non-Newtonian fluid behavior, which is different from its counterpart for Newtonian fluids. This general Smoluchowski velocity is of practical significance in determining the flow rates in microfluidic devices involving non-Newtonian Power-law fluids. [source]

Mixed finite element formulation of algorithms for double-diffusive convection in a fluid-saturated porous medium

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 11 2007
J. C.-F.
Abstract The consistent splitting scheme involving a mixed finite element method for considering the influence of the Forchheimer-extended Brinkman,Darcy model in the momentum equation is applied to double-diffusive convection in a fluid-saturated porous medium. It is shown that the method is robust and can accurately predict flow, pressure distribution, temperature and concentration fields. The numerical scheme may be an alternative to some other existing methods for the solution of porous thermosolutal convection problems since its implementation is very handy. Copyright © 2006 John Wiley & Sons, Ltd. [source]

Hybrid finite element/volume method for shallow water equations

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 13 2010
Shahrouz Aliabadi
Abstract A hybrid numerical scheme based on finite element and finite volume methods is developed to solve shallow water equations. In the recent past, we introduced a series of hybrid methods to solve incompressible low and high Reynolds number flows for single and two-fluid flow problems. The present work extends the application of hybrid method to shallow water equations. In our hybrid shallow water flow solver, we write the governing equations in non-conservation form and solve the non-linear wave equation using finite element method with linear interpolation functions in space. On the other hand, the momentum equation is solved with highly accurate cell-center finite volume method. Our hybrid numerical scheme is truly a segregated method with primitive variables stored and solved at both node and element centers. To enhance the stability of the hybrid method around discontinuities, we introduce a new shock capturing which will act only around sharp interfaces without sacrificing the accuracy elsewhere. Matrix-free GMRES iterative solvers are used to solve both the wave and momentum equations in finite element and finite volume schemes. Several test problems are presented to demonstrate the robustness and applicability of the numerical method. Copyright © 2010 John Wiley & Sons, Ltd. [source]

An embedded Dirichlet formulation for 3D continua

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 5 2010
A. Gerstenberger
Abstract This paper presents a new approach for imposing Dirichlet conditions weakly on non-fitting finite element meshes. Such conditions, also called embedded Dirichlet conditions, are typically, but not exclusively, encountered when prescribing Dirichlet conditions in the context of the eXtended finite element method (XFEM). The method's key idea is the use of an additional stress field as the constraining Lagrange multiplier function. The resulting mixed/hybrid formulation is applicable to 1D, 2D and 3D problems. The method does not require stabilization for the Lagrange multiplier unknowns and allows the complete condensation of these unknowns on the element level. Furthermore, only non-zero diagonal-terms are present in the tangent stiffness, which allows the straightforward application of state-of-the-art iterative solvers, like algebraic multigrid (AMG) techniques. Within this paper, the method is applied to the linear momentum equation of an elastic continuum and to the transient, incompressible Navier,Stokes equations. Steady and unsteady benchmark computations show excellent agreement with reference values. The general formulation presented in this paper can also be applied to other continuous field problems. Copyright © 2009 John Wiley & Sons, Ltd. [source]

A two-step Taylor-characteristic-based Galerkin method for incompressible flows and its application to flow over triangular cylinder with different incidence angles

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 11 2010
Yan Bao
Abstract An alternative characteristic-based scheme, the two-step Taylor-characteristic-based Galerkin method is developed based on the introduction of multi-step temporal Taylor series expansion up to second order along the characteristic of the momentum equation. Contrary to the classical characteristic-based split (CBS) method, the current characteristic-based method does not require splitting the momentum equation, and segregate the calculation of the pressure from that of the velocity by using the momentum,pressure Poisson equation method. Some benchmark problems are used to examine the effectiveness of the proposed algorithm and to compare with the original CBS method, and the results show that the proposed method has preferable accuracy with less numerical dissipation. We further applied the method to the numerical simulation of flow around equilateral triangular cylinder with different incidence angles in free stream. In this numerical investigation, the flow simulations are carried out in the low Reynolds number range. Instantaneous streamlines around the cylinder are used as a means to visualize the wake region behind, and they clearly show the flow pattern around the cylinder in time. The influence of incidence angle on flow characteristic parameters such as Strouhal number, Drag and Lift coefficients are discussed quantitatively. Copyright © 2009 John Wiley & Sons, Ltd. [source]

Two-dimensional modeling for stability analysis of two-phase stratified flow

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 7 2010
Ghassem Heidarinejad
Abstract The effect of wavelength and relative velocity on the disturbed interface of two-phase stratified regime is modeled and discussed. To analyze the stability, a small perturbation is imposed on the interface. Growth or decline of the disturbed wave, relative velocity, and surface tension with respect to time will be discussed numerically. Newly developed scheme applied to a two-dimensional flow field and the governing Navier,Stokes equations in laminar regime are solved. Finite volume method together with non-staggered curvilinear grid is a very effective approach to capture interface shape with time. Because of the interface shape, for any time advancement, a new grid is performed separately on each stratified field, liquid, and gas regime. The results are compared with the analytical characteristics method and one-dimensional modeling. This comparison shows that solving the momentum equation including viscosity term leads to physically more realistic results. In addition, the newly developed method is capable of predicting two-phase stratified flow behavior more precisely than one-dimensional modeling. It was perceived that the surface tension has an inevitable role in dissipation of interface instability and convergence of the two-phase flow model. Copyright © 2009 John Wiley & Sons, Ltd. [source]

Friction term discretization and limitation to preserve stability and conservation in the 1D shallow-water model: Application to unsteady irrigation and river flow

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 4 2008
J. Burguete
Abstract Friction is one of the relevant forces included in the momentum equation of the 1D shallow-water model. This work shows that a pointwise discretization of the friction term unbalances this term with the rest of the terms in the equation in steady state. On the other hand, an upwind discretization of the friction term ensures the correct discrete balance. Furthermore, a conservative technique based on the limitation of the friction value is proposed in order to avoid unbounded values of the friction term in unsteady cases of advancing front over dry and rough surfaces. This limitation improves the quality of unsteady solutions in wet/dry fronts and guarantees the numerical stability in cases with dominant friction terms. The proposed discretization is validated in some test cases with analytical solution or with measured data and used in some practical cases. Copyright © 2007 John Wiley & Sons, Ltd. [source]

Numerical aspects of improvement of the unsteady pipe flow equations

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 11 2007
Romuald Szymkiewicz
Abstract The paper presents an analysis of some recently proposed improvements of the water hammer equations, which concern the friction term in the momentum equation. A comparison of the experimental data and numerical results shows that the required damping and smoothing of the pressure wave cannot be obtained by modification of the friction factor only. In order to evaluate the significance of the introduced improvements into the momentum equation, the accuracy of the numerical solution has been analysed using the modified equation approach. The analysis shows why the physical dissipation process observed in the water hammer phenomenon cannot be reproduced with the commonly used source term in Darcy,Weisbach form, representing friction force in the momentum equation. Therefore, regardless of the proposed form of the friction factor for unsteady flow, the model of water hammer improved in such a way keeps its hyperbolic character. Consequently, it cannot ensure the expected effects of damping and smoothing of the calculation head oscillations. Copyright © 2007 John Wiley & Sons, Ltd. [source]

Using vorticity to define conditions at multiple open boundaries for simulating flow in a simplified vortex settling basin

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 1 2007
A. N. Ziaei
Abstract In this paper a method is developed to define multiple open boundary (OB) conditions in a simplified vortex settling basin (VSB). In this method, the normal component of the momentum equation is solved at the OBs, and tangential components of vorticity are calculated by solving vorticity transport equations only at the OBs. Then the tangential vorticity components are used to construct Neumann boundary conditions for tangential velocity components. Pressure is set to its ambient value, and the divergence-free condition is satisfied at these boundaries by employing the divergence as the Neumann condition for the normal-direction momentum equation. The 3-D incompressible Navier,Stokes equations in a primitive-variable form are solved using the SIMPLE algorithm. Grid-function convergence tests are utilized to verify the numerical results. The complicated laminar flow structure in the VSB is investigated, and preliminary assessment of two popular turbulence models, k,, and k,,, is conducted. Copyright © 2006 John Wiley & Sons, Ltd. [source]

Pressure boundary condition for the time-dependent incompressible Navier,Stokes equations

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 6 2006
R. L. Sani
Abstract In Gresho and Sani (Int. J. Numer. Methods Fluids 1987; 7:1111,1145; Incompressible Flow and the Finite Element Method, vol. 2. Wiley: New York, 2000) was proposed an important hypothesis regarding the pressure Poisson equation (PPE) for incompressible flow: Stated there but not proven was a so-called equivalence theorem (assertion) that stated/asserted that if the Navier,Stokes momentum equation is solved simultaneously with the PPE whose boundary condition (BC) is the Neumann condition obtained by applying the normal component of the momentum equation on the boundary on which the normal component of velocity is specified as a Dirichlet BC, the solution (u, p) would be exactly the same as if the ,primitive' equations, in which the PPE plus Neumann BC is replaced by the usual divergence-free constraint (, · u = 0), were solved instead. This issue is explored in sufficient detail in this paper so as to actually prove the theorem for at least some situations. Additionally, like the original/primitive equations that require no BC for the pressure, the new results establish the same thing when the PPE approach is employed. Copyright © 2005 John Wiley & Sons, Ltd. [source]

Flux and source term discretization in two-dimensional shallow water models with porosity on unstructured grids

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 3 2006
Vincent Guinot
Abstract Two-dimensional shallow water models with porosity appear as an interesting path for the large-scale modelling of floodplains with urbanized areas. The porosity accounts for the reduction in storage and in the exchange sections due to the presence of buildings and other structures in the floodplain. The introduction of a porosity into the two-dimensional shallow water equations leads to modified expressions for the fluxes and source terms. An extra source term appears in the momentum equation. This paper presents a discretization of the modified fluxes using a modified HLL Riemann solver on unstructured grids. The source term arising from the gradients in the topography and in the porosity is treated in an upwind fashion so as to enhance the stability of the solution. The Riemann solver is tested against new analytical solutions with variable porosity. A new formulation is proposed for the macroscopic head loss in urban areas. An application example is presented, where the large scale model with porosity is compared to a refined flow model containing obstacles that represent a schematic urban area. The quality of the results illustrates the potential usefulness of porosity-based shallow water models for large scale floodplain simulations. Copyright © 2005 John Wiley & Sons, Ltd. [source]

A 2D implicit time-marching algorithm for shallow water models based on the generalized wave continuity equation

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 3 2004
Kendra M. Dresback
Abstract This paper builds upon earlier work that developed and evaluated a 1D predictor,corrector time-marching algorithm for wave equation models and extends it to 2D. Typically, the generalized wave continuity equation (GWCE) utilizes a three time-level semi-implicit scheme centred at k, and the momentum equation uses a two time-level scheme centred at k+12. It has been shown that in highly non-linear applications, the algorithm becomes unstable at even moderate Courant numbers. This work implements and analyses an implicit treatment of the non-linear terms through the use of an iterative time-marching algorithm in the two-dimensional framework. Stability results show at least an eight-fold increase in the maximum time step, depending on the domain. Studies also examined the sensitivity of the G parameter (a numerical weighting parameter in the GWCE) with results showing the greatest increase in stability occurs when 1,G/,max,10, a range that coincides with the recommended range to minimize errors. Convergence studies indicate an increase in temporal accuracy from first order to second order, while overall error is less than the original algorithm, even at higher time steps. Finally, a parallel implementation of the new algorithm shows that it scales well. Copyright © 2004 John Wiley & Sons, Ltd. [source]

An efficient finite difference scheme for free-surface flows in narrow rivers and estuaries

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 3 2003
XinJian ChenArticle first published online: 13 MAY 200
Abstract This paper presents a free-surface correction (FSC) method for solving laterally averaged, 2-D momentum and continuity equations. The FSC method is a predictor,corrector scheme, in which an intermediate free surface elevation is first calculated from the vertically integrated continuity equation after an intermediate, longitudinal velocity distribution is determined from the momentum equation. In the finite difference equation for the intermediate velocity, the vertical eddy viscosity term and the bottom- and sidewall friction terms are discretized implicitly, while the pressure gradient term, convection terms, and the horizontal eddy viscosity term are discretized explicitly. The intermediate free surface elevation is then adjusted by solving a FSC equation before the intermediate velocity field is corrected. The finite difference scheme is simple and can be easily implemented in existing laterally averaged 2-D models. It is unconditionally stable with respect to gravitational waves, shear stresses on the bottom and side walls, and the vertical eddy viscosity term. It has been tested and validated with analytical solutions and field data measured in a narrow, riverine estuary in southwest Florida. Model simulations show that this numerical scheme is very efficient and normally can be run with a Courant number larger than 10. It can be used for rivers where the upstream bed elevation is higher than the downstream water surface elevation without any problem. Copyright © 2003 John Wiley & Sons, Ltd. [source]

Computation of unsteady viscous incompressible flows in generalized non-inertial co-ordinate system using Godunov-projection method and overlapping meshes

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 11 2002
H. Pan
Abstract Time-dependent incompressible Navier,Stokes equations are formulated in generalized non-inertial co-ordinate system and numerically solved by using a modified second-order Godunov-projection method on a system of overlapped body-fitted structured grids. The projection method uses a second-order fractional step scheme in which the momentum equation is solved to obtain the intermediate velocity field which is then projected on to the space of divergence-free vector fields. The second-order Godunov method is applied for numerically approximating the non-linear convection terms in order to provide a robust discretization for simulating flows at high Reynolds number. In order to obtain the pressure field, the pressure Poisson equation is solved. Overlapping grids are used to discretize the flow domain so that the moving-boundary problem can be solved economically. Numerical results are then presented to demonstrate the performance of this projection method for a variety of unsteady two- and three-dimensional flow problems formulated in the non-inertial co-ordinate systems. Copyright © 2002 John Wiley & Sons, Ltd. [source]

A finite volume solver for 1D shallow-water equations applied to an actual river

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 1 2002
N. Gouta
Abstract This paper describes the numerical solution of the 1D shallow-water equations by a finite volume scheme based on the Roe solver. In the first part, the 1D shallow-water equations are presented. These equations model the free-surface flows in a river. This set of equations is widely used for applications: dam-break waves, reservoir emptying, flooding, etc. The main feature of these equations is the presence of a non-conservative term in the momentum equation in the case of an actual river. In order to apply schemes well adapted to conservative equations, this term is split in two terms: a conservative one which is kept on the left-hand side of the equation of momentum and the non-conservative part is introduced as a source term on the right-hand side. In the second section, we describe the scheme based on a Roe Solver for the homogeneous problem. Next, the numerical treatment of the source term which is the essential point of the numerical modelisation is described. The source term is split in two components: one is upwinded and the other is treated according to a centred discretization. By using this method for the discretization of the source term, one gets the right behaviour for steady flow. Finally, in the last part, the problem of validation is tackled. Most of the numerical tests have been defined for a working group about dam-break wave simulation. A real dam-break wave simulation will be shown. Copyright © 2002 John Wiley & Sons, Ltd. [source]

High-order filtering for control volume flow simulation

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 7 2001
G. De Stefano
Abstract A general methodology is presented in order to obtain a hierarchy of high-order filter functions, starting from the standard top-hat filter, naturally linked to control volumes flow simulations. The goal is to have a new filtered variable better represented in its high resolved wavenumber components by using a suitable deconvolution. The proposed formulation is applied to the integral momentum equation, that is the evolution equation for the top-hat filtered variable, by performing a spatial reconstruction based on the approximate inversion of the averaging operator. A theoretical analysis for the Burgers' model equation is presented, demonstrating that the local de-averaging is an effective tool to obtain a higher-order accuracy. It is also shown that the subgrid-scale term, to be modeled in the deconvolved balance equation, has a smaller absolute importance in the resolved wavenumber range for increasing deconvolution order. A numerical analysis of the procedure is presented, based on high-order upwind and central fluxes reconstruction, leading to congruent control volume schemes. Finally, the features of the present high-order conservative formulation are tested in the numerical simulation of a sample turbulent flow: the flow behind a backward-facing step. Copyright © 2001 John Wiley & Sons, Ltd. [source]

Diffusion in poro-plastic media

MATHEMATICAL METHODS IN THE APPLIED SCIENCES, Issue 18 2004
R. E. Showalter
Abstract A model is developed for the flow of a slightly compressible fluid through a saturated inelastic porous medium. The initial-boundary-value problem is a system that consists of the diffusion equation for the fluid coupled to the momentum equation for the porous solid together with a constitutive law which includes a possibly hysteretic relation of elasto-visco-plastic type. The variational form of this problem in Hilbert space is a non-linear evolution equation for which the existence and uniqueness of a global strong solution is proved by means of monotonicity methods. Various degenerate situations are permitted, such as incompressible fluid, negligible porosity, or a quasi-static momentum equation. The essential sufficient conditions for the well-posedness of the system consist of an ellipticity condition on the term for diffusion of fluid and either a viscous or a hardening assumption in the constitutive relation for the porous solid. Copyright © 2004 John Wiley & Sons, Ltd. [source]

Treatment of vector equations in deep-atmosphere, semi-Lagrangian models.

THE QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY, Issue 647 2010
I: Momentum equation
Abstract Application of the semi-Lagrangian method to the vector momentum equation in orthogonal curvilinear systems is considered, with emphasis on spherical and spheroidal coordinates for deep-atmosphere models (in which the shallow-atmosphere approximation is not made). In spherical coordinates, a certain rotation matrix allows vector components at semi-Lagrangian departure points to be expressed in the required component forms at the corresponding arrival points. This rotation matrix also allows the unit vector triad at the arrival point to be expressed in terms of the unit vector triad at the departure point. The resulting formalisation of the momentum equation applies to timesteps of arbitrary length, and in the limit of infinitesimal timestep it delivers the familiar Eulerian components. Extension of the rotation matrix technique to more general orthogonal curvilinear systems is remarkably straightforward, thus demonstrating the versatility of the semi-Lagrangian method. Explicit forms are given for systems having longitude as one coordinate and general orthogonal coordinates in the meridional plane, and in particular these forms are applicable to both confocal oblate spheroidal and similar oblate spheroidal systems. A specific factorisation of the rotation matrix in the spherical polar case involves three matrices, one of which represents rotation of the unit vector triad along a great circle arc; the non-Euclidean, shallow-atmosphere approximation of the total rotation matrix results when the great circle rotation matrix is replaced by the identity matrix (but the other two matrices involved in the factorisation are kept intact). This shallow-atmosphere rotation matrix agrees with that used by ECMWF and Météo-France. © Crown Copyright 2010. Published by John Wiley & Sons, Ltd. [source]

Balanced boundary layers used in hurricane models

THE QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY, Issue 635 2008
Roger K. Smith
Abstract We examine the formulation and accuracy of various approximations made in representing the boundary layer in simple axisymmetric hurricane models, especially those that assume strict gradient wind balance in the radial direction. Approximate solutions for a steady axisymmetric slab boundary-layer model are compared with a full model solution. It is shown that the approximate solutions are generally poor in the inner core region of the vortex, where the radial advection term in the radial momentum equation is important and cannot be neglected. These results affirm some prior work and have implications for a range of theoretical studies of hurricane dynamics, including theories of potential intensity, that employ balanced boundary-layer formulations. Copyright © 2008 Royal Meteorological Society [source]

A critique of Emanuel's hurricane model and potential intensity theory

THE QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY, Issue 632 2008
Roger K. Smith
Abstract We present a critique of Emanuel's steady-state hurricane model, which is a precursor to his theory for hurricane potential intensity (PI). We show that a major deficiency of the theory is the tacit assumption of gradient wind balance in the boundary layer, a layer that owes its existence to gradient wind imbalance in the radial momentum equation. If a more complete boundary-layer formulation is included using the gradient wind profiles obtained from Emanuel's theory, the tangential wind speed in the boundary layer becomes supergradient, invalidating the assumption of gradient wind balance. We show that the degree to which the tangential wind is supergradient depends on the assumed boundary-layer depth. The full boundary-layer solutions require a knowledge of the tangential wind profile above the boundary layer in the outer region where there is subsidence into the layer and they depend on the breadth of this profile. This effect is not considered in Emanuel's theory. We argue that a more complete theory for the steady-state hurricane would require the radial pressure gradient above the boundary layer to be prescribed or determined independently of the boundary layer. The issues raised herein highlight a fundamental problem with Emanuel's theory for PI, since that theory makes the same assumptions as in the steady-state hurricane model. Our current findings together with recent studies examining intense hurricanes suggest a way forward towards a more consistent theory for hurricane PI. Copyright © 2008 Royal Meteorological Society [source]

A self-consistent treatment of the electromotive force in magnetohydrodynamics for large diffusivities

ASTRONOMISCHE NACHRICHTEN, Issue 7 2010
A. Courvoisier
Abstract The coupled equations that describe the effect of large-scale magnetic and velocity fields on forced high-diffusivity magnetohydrodynamic flows are investigated through an extension of mean field electrodynamics. Our results generalise those of Rädler & Brandenburg (2010), who consider a similar situation but assume that the effect of the Lorentz force on the momentum equation can be neglected. New mean coupling terms are shown to appear, which can lead to large-scale growth of magnetic and velocity fields even when the usual a-effects are absent (© 2010 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim) [source]

Heat transfer for Marangoni-driven boundary layer flow

HEAT TRANSFER - ASIAN RESEARCH (FORMERLY HEAT TRANSFER-JAPANESE RESEARCH), Issue 2 2002
David M. Christopher
Abstract Marangoni convection induced by variation of the surface tension with temperature along a surface influences crystal growth melts and other processes with liquid,vapor interfaces, such as boiling in both microgravity and normal gravity in some cases. This paper presents the Nusselt number for Marangoni flow over a flat surface calculated using a similarity solution for both the momentum equations and the energy equation assuming developing boundary layer flow along a surface. Solutions are presented for the surface velocity, the total flow rate, and the Nusselt number for various temperature profiles, Marangoni numbers, and Prandtl numbers. For large bubbles, the predicted boundary layer thickness would be less than the bubble diameter, so the curvature effects could be neglected and this analysis could be used as a first estimate of the effect of Marangoni flow around a vapor bubble. © 2002 Scripta Technica, Heat Trans Asian Res, 31(2): 105,116, 2002; DOI 10.1002/htj.10019 [source]

Hybrid finite element/volume method for shallow water equations

Meshfree simulation of failure modes in thin cylinders subjected to combined loads of internal pressure and localized heat

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 8 2008
Dong Qian
Abstract This paper focuses on the non-linear responses in thin cylindrical structures subjected to combined mechanical and thermal loads. The coupling effects of mechanical deformation and temperature in the material are considered through the development of a thermo-elasto-viscoplastic constitutive model at finite strain. A meshfree Galerkin approach is used to discretize the weak forms of the energy and momentum equations. Due to the different time scales involved in thermal conduction and failure development, an explicit,implicit time integration scheme is developed to link the time scale differences between the two key mechanisms. We apply the developed approach to the analysis of the failure of cylindrical shell subjected to both heat sources and internal pressure. The numerical results show four different failure modes: dynamic fragmentation, single crack with branch, thermally induced cracks and cracks due to the combined effects of pressure and temperature. These results illustrate the important roles of thermal and mechanical loads with different time scales. Copyright © 2008 John Wiley & Sons, Ltd. [source]

Unstructured finite volume discretization of two-dimensional depth-averaged shallow water equations with porosity

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 8 2010
L. Cea
Abstract This paper deals with the numerical discretization of two-dimensional depth-averaged models with porosity. The equations solved by these models are similar to the classic shallow water equations, but include additional terms to account for the effect of small-scale impervious obstructions which are not resolved by the numerical mesh because their size is smaller or similar to the average mesh size. These small-scale obstructions diminish the available storage volume on a given region, reduce the effective cross section for the water to flow, and increase the head losses due to additional drag forces and turbulence. In shallow water models with porosity these effects are modelled introducing an effective porosity parameter in the mass and momentum conservation equations, and including an additional drag source term in the momentum equations. This paper presents and compares two different numerical discretizations for the two-dimensional shallow water equations with porosity, both of them are high-order schemes. The numerical schemes proposed are well-balanced, in the sense that they preserve naturally the exact hydrostatic solution without the need of high-order corrections in the source terms. At the same time they are able to deal accurately with regions of zero porosity, where the water cannot flow. Several numerical test cases are used in order to verify the properties of the discretization schemes proposed. Copyright © 2009 John Wiley & Sons, Ltd. [source]

An approximate-state Riemann solver for the two-dimensional shallow water equations with porosity

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 12 2010
P. Finaud-Guyot
Abstract PorAS, a new approximate-state Riemann solver, is proposed for hyperbolic systems of conservation laws with source terms and porosity. The use of porosity enables a simple representation of urban floodplains by taking into account the global reduction in the exchange sections and storage. The introduction of the porosity coefficient induces modified expressions for the fluxes and source terms in the continuity and momentum equations. The solution is considered to be made of rarefaction waves and is determined using the Riemann invariants. To allow a direct computation of the flux through the computational cells interfaces, the Riemann invariants are expressed as functions of the flux vector. The application of the PorAS solver to the shallow water equations is presented and several computational examples are given for a comparison with the HLLC solver. Copyright © 2009 John Wiley & Sons, Ltd. [source]

Numerical simulation of bubble and droplet deformation by a level set approach with surface tension in three dimensions

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 9 2010
Roberto Croce
Abstract In this paper we present a three-dimensional Navier,Stokes solver for incompressible two-phase flow problems with surface tension and apply the proposed scheme to the simulation of bubble and droplet deformation. One of the main concerns of this study is the impact of surface tension and its discretization on the overall convergence behavior and conservation properties. Our approach employs a standard finite difference/finite volume discretization on uniform Cartesian staggered grids and uses Chorin's projection approach. The free surface between the two fluid phases is tracked with a level set (LS) technique. Here, the interface conditions are implicitly incorporated into the momentum equations by the continuum surface force method. Surface tension is evaluated using a smoothed delta function and a third-order interpolation. The problem of mass conservation for the two phases is treated by a reinitialization of the LS function employing a regularized signum function and a global fixed point iteration. All convective terms are discretized by a WENO scheme of fifth order. Altogether, our approach exhibits a second-order convergence away from the free surface. The discretization of surface tension requires a smoothing scheme near the free surface, which leads to a first-order convergence in the smoothing region. We discuss the details of the proposed numerical scheme and present the results of several numerical experiments concerning mass conservation, convergence of curvature, and the application of our solver to the simulation of two rising bubble problems, one with small and one with large jumps in material parameters, and the simulation of a droplet deformation due to a shear flow in three space dimensions. Furthermore, we compare our three-dimensional results with those of quasi-two-dimensional and two-dimensional simulations. This comparison clearly shows the need for full three-dimensional simulations of droplet and bubble deformation to capture the correct physical behavior. Copyright © 2009 John Wiley & Sons, Ltd. [source]

Fluid,solid interaction problems with thermal convection using the immersed element-free Galerkin method

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 1 2010
Claudio M. Pita
Abstract In this work, the immersed element-free Galerkin method (IEFGM) is proposed for the solution of fluid,structure interaction (FSI) problems. In this technique, the FSI is represented as a volumetric force in the momentum equations. In IEFGM, a Lagrangian solid domain moves on top of an Eulerian fluid domain that spans over the entire computational region. The fluid domain is modeled using the finite element method and the solid domain is modeled using the element-free Galerkin method. The continuity between the solid and fluid domains is satisfied by means of a local approximation, in the vicinity of the solid domain, of the velocity field and the FSI force. Such an approximation is achieved using the moving least-squares technique. The method was applied to simulate the motion of a deformable disk moving in a viscous fluid due to the action of the gravitational force and the thermal convection of the fluid. An analysis of the main factors affecting the shape and trajectory of the solid body is presented. The method shows a distinct advantage for simulating FSI problems with highly deformable solids. Copyright © 2009 John Wiley & Sons, Ltd. [source]

Momentum/continuity coupling with large non-isotropic momentum source terms

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 9 2009
J. D. Franklin
Abstract Pressure-based methods such as the SIMPLE algorithm are frequently used to determine a coupled solution between the component momentum equations and the continuity equation. This paper presents a colocated variable pressure correction algorithm for control volumes of polyhedral/polygonal cell topologies. The correction method is presented independent of spatial approximation. The presence of non-isotropic momentum source terms is included in the proposed algorithm to ensure its applicability to multi-physics applications such as gas and particulate flows. Two classic validation test cases are included along with a newly proposed test case specific to multiphase flows. The classic validation test cases demonstrate the application of the proposed algorithm on truly arbitrary polygonal/polyhedral cell meshes. A comparison between the current algorithm and commercially available software is made to demonstrate that the proposed algorithm is competitively efficient. The newly proposed test case demonstrates the benefits of the current algorithm when applied to a multiphase flow situation. The numerical results from this case show that the proposed algorithm is more robust than other methods previously proposed. Copyright © 2009 John Wiley & Sons, Ltd. [source]

Incorporating spatially variable bottom stress and Coriolis force into 2D, a posteriori, unstructured mesh generation for shallow water models

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 3 2009
D. Michael Parrish
Abstract An enhanced version of our localized truncation error analysis with complex derivatives (LTEA,CD ) a posteriori approach to computing target element sizes for tidal, shallow water flow, LTEA+CD , is applied to the Western North Atlantic Tidal model domain. The LTEA + CD method utilizes localized truncation error estimates of the shallow water momentum equations and builds upon LTEA and LTEA,CD-based techniques by including: (1) velocity fields from a nonlinear simulation with complete constituent forcing; (2) spatially variable bottom stress; and (3) Coriolis force. Use of complex derivatives in this case results in a simple truncation error expression, and the ability to compute localized truncation errors using difference equations that employ only seven to eight computational points. The compact difference molecules allow the computation of truncation error estimates and target element sizes throughout the domain, including along the boundary; this fact, along with inclusion of locally variable bottom stress and Coriolis force, constitute significant advancements beyond the capabilities of LTEA. The goal of LTEA + CD is to drive the truncation error to a more uniform, domain-wide value by adjusting element sizes (we apply LTEA + CD by re-meshing the entire domain, not by moving nodes). We find that LTEA + CD can produce a mesh that is comprised of fewer nodes and elements than an initial high-resolution mesh while performing as well as the initial mesh when considering the resynthesized tidal signals (elevations). Copyright © 2008 John Wiley & Sons, Ltd. [source]