International Journal for Numerical Methods in Fluids

Wong, K. L.; Baker, A. J.

doi: 10.1002/fld.204pmid: N/A

The velocity–vorticity formulation is selected to develop a time‐accurate CFD finite element algorithm for the incompressible Navier–Stokes equations in three dimensions.The finite element implementation uses equal order trilinear finite elements on a non‐staggered hexahedral mesh. A second order vorticity kinematic boundary condition is derived for the no slip wall boundary condition which also enforces the incompressibility constraint. A biconjugate gradient stabilized (BiCGSTAB) sparse iterative solver is utilized to solve the fully coupled system of equations as a Newton algorithm. The solver yields an efficient parallel solution algorithm on distributed‐memory machines, such as the IBM SP2. Three dimensional laminar flow solutions for a square channel, a lid‐driven cavity, and a thermal cavity are established and compared with available benchmark solutions. Copyright © 2002 John Wiley & Sons, Ltd.

Kim, Kyoungyoun; Baek, Seung‐Jin; Sung, Hyung Jin

doi: 10.1002/fld.205pmid: N/A

An efficient numerical method to solve the unsteady incompressible Navier–Stokes equations is developed. A fully implicit time advancement is employed to avoid the Courant–Friedrichs–Lewy restriction, where the Crank–Nicolson discretization is used for both the diffusion and convection terms. Based on a block LU decomposition, velocity–pressure decoupling is achieved in conjunction with the approximate factorization. The main emphasis is placed on the additional decoupling of the intermediate velocity components with only nth time step velocity. The temporal second‐order accuracy is preserved with the approximate factorization without any modification of boundary conditions. Since the decoupled momentum equations are solved without iteration, the computational time is reduced significantly. The present decoupling method is validated by solving several test cases, in particular, the turbulent minimal channel flow unit. Copyright © 2002 John Wiley & Sons, Ltd.

Sanders, Brett F.; Bradford, Scott F.

doi: 10.1002/fld.206pmid: N/A

A monotone, second‐order accurate numerical scheme is presented for solving the differential form of the adjoint shallow‐water equations in generalized two‐dimensional coordinates. Fluctuation‐splitting is utilized to achieve a high‐resolution solution of the equations in primitive form. One‐step and two‐step schemes are presented and shown to achieve solutions of similarly high accuracy in one dimension. However, the two‐step method is shown to yield more accurate solutions to problems in which unsteady wave speeds are present. In two dimensions, the two‐step scheme is tested in the context of two parameter identification problems, and it is shown to accurately transmit the information needed to identify unknown forcing parameters based on measurements of the system response. The first problem involves the identification of an upstream flood hydrograph based on downstream depth measurements. The second problem involves the identification of a long wave state in the far‐field based on near‐field depth measurements. Copyright © 2002 John Wiley & Sons, Ltd.

Heggelund, Yngve; Berntsen, Jarle

doi: 10.1002/fld.215pmid: N/A

A method for analysing different nesting techniques for the linearized shallow water equations is presented. The problem is formulated as an eigenvector–eigenvalue problem. A necessary condition for stability is that the spectral radius of the propagation matrix is less than or equal to one. Two test cases are presented. The first test case is analysed, and effects of enforcing volume conservation and nudging in time are studied. A nesting technique is found that causes no growth of any eigenvectors for reasonable time steps. This nesting technique is then used on both test cases, and results are compared to an everywhere refined model and a coarse grid model. Copyright © 2002 John Wiley & Sons, Ltd.

Li, L.; Sherwin, S. J.; Bearman, P. W.

doi: 10.1002/fld.216pmid: N/A

A mathematical and numerical formulation is derived for fluid/structure interaction problems involving arbitrary geometries relevant to the simulation of bridge deck instabilities due to cross winds. A translating and rotating moving frame of reference is attached to the body to utilize an efficient fixed mesh spectral/hp element solver. The formulation is validated against experiments with flow simulations of circular cylinders at Reynolds numbers of 100–400 undergoing free and forced motion in the transverse and in‐line directions. The well‐documented phenomena of vortex lock‐in is captured. The formulation is then applied a rectangular body at Re=250 under forced and free motion the latter of which demonstrates torsional galloping. Copyright © 2002 John Wiley & Sons, Ltd.