On the application of slip boundary condition on curved boundaries

Hydrodynamic simulations of sloshing phenomena often involve the application of slip boundary condition at the wetted surfaces. If these surfaces are curved, the ambiguous nature of the normal vector in the discretized problem can interfere with the application of such a boundary condition. Even the use of consistent normal vectors, preferred from the point of view of conservation, does not assure good approximation of the continuum slip condition in the discrete problem, and non‐physical recirculating flow fields may be observed. As a remedy, we consider the Navier slip condition, and more successfully, the so‐called BC‐free boundary condition. Copyright © 2004 John Wiley & Sons, Ltd.     

A time‐accurate pseudo‐compressibility approach based on an unstructured hybrid finite volume technique applied to unsteady turbulent premixed flame propagation

A preconditioning approach based on the artificial compressibility formulation is extended to solve the governing equations for unsteady turbulent reactive flows with heat release, at low Mach numbers, on an unstructured hybrid grid context. Premixed reactants are considered and a flamelet approach for combustion modelling is adopted using a continuous quenched mean reaction rate. An overlapped cell‐vertex finite volume method is adopted as a discretisation scheme. Artificial dissipation terms for hybrid grids are explicitly added to ensure a stable, discretised set of equations. A second‐order, explicit, hybrid Runge–Kutta scheme is applied for the time marching in pseudo‐time. A time derivative of the dependent variable is added to recover the time accuracy of the preconditioned set of equations. This derivative is discretised by an implicit, second‐order scheme. The resulting scheme is applied to the calculation of an infinite planar (one‐dimensional) turbulent premixed flame propagating freely in reactants whose turbulence is supposed to be frozen, homogeneous and isotropic. The accuracy of the results obtained with the proposed method proves to be excellent when compared to the data available in the literature. Copyright © 2004 John Wiley & Sons, Ltd.     

Reference values for drag and lift of a two‐dimensional time‐dependent flow around a cylinder

This paper presents a numerical study of a two‐dimensional time‐dependent flow around a cylinder. Its main objective is to provide accurate reference values for the maximal drag and lift coefficient at the cylinder and for the pressure difference between the front and the back of the cylinder at the final time. In addition, the accuracy of these values obtained with different time stepping schemes and different finite element methods is studied. Copyright © 2004 John Wiley & Sons, Ltd.     

Assessment of two‐equation turbulence modelling for high Reynolds number hydrofoil flows

This paper presents an evaluation of the capability of turbulence models available in the commercial CFD code FLUENT 6.0 for their application to hydrofoil turbulent boundary layer separation flow at high Reynolds numbers. Four widely applied two‐equation RANS turbulence models were assessed through comparison with experimental data at Reynolds numbers of 8.284×106 and 1.657×107. They were the standard k–εmodel, the realizable k–εmodel, the standard k–ωmodel and the shear‐stress‐transport (SST) k–ωmodel. It has found that the realizable k–εturbulence model used with enhanced wall functions and near‐wall modelling techniques, consistently provides superior performance in predicting the flow characteristics around the hydrofoil. Copyright © 2004 John Wiley & Sons, Ltd.     

A virtual boundary method with improved computational efficiency using a multi‐grid method

The flow around spherical, solid objects is considered. The boundary conditions on the solid boundaries have been applied by replacing the boundary with a surface force distribution on the surface, such that the required boundary conditions are satisfied. The velocity on the boundary is determined by extrapolation from the flow field. The source terms are determined iteratively, as part of the solution. They are then averaged and are smoothed out to nearby computational grid points. A multi‐grid scheme has been used to enhance the computational efficiency of the solution of the force equations. The method has been evaluated for flow around both moving and stationary spherical objects at very low and intermediate Reynolds numbers. The results shows a second order accuracy of the method both at creeping flow and at Re=100. The multi‐grid scheme is shown to enhance the convergence rate up to a factor 10 as compared to single grid approach. Copyright © 2004 John Wiley & Sons, Ltd.     

Least‐squares meshfree method for incompressible Navier–Stokes problems

A least‐squares meshfree method based on the first‐order velocity–pressure–vorticity formulation for two‐dimensional incompressible Navier–Stokes problem is presented. The convective term is linearized by successive substitution or Newton's method. The discretization of all governing equations is implemented by the least‐squares method. Equal‐order moving least‐squares approximation is employed with Gauss quadrature in the background cells. The boundary conditions are enforced by the penalty method. The matrix‐free element‐by‐element Jacobi preconditioned conjugate method is applied to solve the discretized linear systems. Cavity flow for steady Navier–Stokes problem and the flow over a square obstacle for time‐dependent Navier–Stokes problem are investigated for the presented least‐squares meshfree method. The effects of inaccurate integration on the accuracy of the solution are investigated. Copyright © 2004 John Wiley & Sons, Ltd.     

A subdomain boundary element method for high‐Reynolds laminar flow using stream function‐vorticity formulation

The paper presents a new formulation of the integral boundary element method (BEM) using subdomain technique. A continuous approximation of the function and the function derivative in the direction normal to the boundary element (further ‘normal flux’) is introduced for solving the general form of a parabolic diffusion‐convective equation. Double nodes for normal flux approximation are used. The gradient continuity is required at the interior subdomain corners where compatibility and equilibrium interface conditions are prescribed. The obtained system matrix with more equations than unknowns is solved using the fast iterative linear least squares based solver. The robustness and stability of the developed formulation is shown on the cases of a backward‐facing step flow and a square‐driven cavity flow up to the Reynolds number value 50 000. Copyright © 2004 John Wiley & Sons, Ltd.     

Molecular dynamics and information on possible sites of interaction of intramyocellular metabolites in vivo from resolved dipolar couplings in localized H NMR spectra

Proton NMR resonances of the endogenous metabolites creatine and phosphocreatine ((P)Cr), taurine (Tau), and carnosine (Cs, β-alanyl-l-histidine) were studied with regard to residual dipolar couplings and molecular mobility. We present an analysis of the direct ¹H–¹H interaction that provides information on motional reorientation of subgroups in these molecules in vivo. For this purpose, localized ¹H NMR experiments were performed on m. gastrocnemius of healthy volunteers using a 1.5-T clinical whole-body MR scanner. We evaluated the observable dipolar coupling strength SD₀ (S = order parameter) of the (P)Cr-methyl triplet and the Tau-methylene doublet by means of the apparent line splitting. These were compared to the dipolar coupling strength of the (P)Cr-methylene doublet. In contrast to the aliphatic protons of (P)Cr and Tau, the aromatic H2 (δ = 8 ppm) and H4 (δ = 7 ppm) protons of the imidazole ring of Cs exhibit second-order spectra at 1.5 T. This effect is the consequence of incomplete transition from Zeeman to Paschen-Back regime and allows a determination of SD₀ from H2 and H4 of Cs as an alternative to evaluating the multiplet splitting which can be measured directly in high-resolution ¹H NMR spectra. Experimental data showed striking differences in the mobility of the metabolites when the dipolar coupling constant D₀ (calculated with the internuclear distance known from molecular geometry in the case of complete absence of molecular dynamics and motion) is used for comparison. The aliphatic signals involve very small order parameters S ≈ (1.4 − 3) × 10⁻⁴ indicating rapid reorientation of the corresponding subgroups in these metabolites. In contrast, analysis of the Cs resonances yielded S ≈ (113 − 137) × 10⁻⁴. Thus, the immobilization of the Cs imidazole ring owing to an anisotropic cellular substructure in human m. gastrocnemius is much more effective than for (P)Cr and Tau subgroups. Furthermore, ¹H NMR experiments on aqueous model solutions of histidine and N-acetyl-l-aspartate (NAA) enabled the assignment of an additional signal component at δ = 8 ppm of Cs in vivo to the amide group at the peptide bond. The visibility of this proton could result from hydrogen bonding which would agree with the anticipated stronger motional restriction of Cs. Referring to the observation that all dipolar-coupled multiplets resolved in localized in vivo ¹H NMR spectra of human m. gastrocnemius collapse simultaneously when the fibre structure is tilted towards the magic angle (θ ≈ 55°), a common model for molecular confinement in muscle tissue is proposed on the basis of an interaction of the studied metabolites with myocellular membrane phospholipids.     

Synthesis and association of poly(N,N‐dimethylacrylamide) copolymers with perfluorocarbon pendent groups connected through polyethylene‐glycol tethers

The synthesis is reported of copolymers of N,N‐dimethylacrylamide (DMA) and methacrylates containing 2,2′‐dihydroperfluorodecanoyl (RF) groups separated from the methacrylate by long polyethylene glycol (PEG) tether groups (between 1000 and 14,000 Da). At concentrations of between 1 and 8 wt % the copolymers with macromonomer contents of 1 mol % or less give gels in organic solvents such as dioxane, THF, or methanol, as well as in water. Given the low molecular weights, this indicates very efficient association of very low numbers of RF groups. Association and gel formation is enormously enhanced in the presence of longer PEG tethers. This is consistent with smaller poly(N,N,‐dimethylacrylamide) (PDMA) intermolecular excluded volume effects that are mediated by the longer PEG tethers and possibly by the incompatibility of PEG and PDMA that may lead to the formation of PEG microdomains. This increases the local concentrations of the RF groups in the PEO domains that are not diluted by the PDMA chains, as would be the case in the absence of PEG tethers. © 2003 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 360–373, 2004     

UV‐ignited frontal polymerization of an epoxy resin

By combining frontal polymerization and radical‐induced cationic polymerization, it was possible to cure thick samples of an epoxy monomer bleached by UV light. The effect of the relative amounts of cationic photoinitiator and radical initiator was thoroughly investigated and was related to the front's velocity and its maximum temperature. The materials obtained were characterized by quantitative conversion also in the deeper layers, not reached by UV light. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 2066–2072, 2004