Assessment of Artificial Dissipation Models for Three-Dimensional Incompressible Flow Solutions

1997 ◽  
Vol 119 (2) ◽  
pp. 331-340 ◽  
Author(s):  
F. B. Lin ◽  
F. Sotiropoulos

Various approaches for constructing artificial dissipation terms for three-dimensional artificial compressibility algorithms are presented and evaluated. Two, second-order accurate, central-differencing schemes, with explicitly added scalar and matrix-valued fourth-difference artificial dissipation, respectively, and a third-order accurate flux-difference splitting upwind scheme are implemented in a multigrid time-stepping procedure and applied to calculate laminar flow through a strongly curved duct. Extensive grid-refinement studies are carried out to investigate the grid sensitivity of each discretization approach. The calculations indicate that even the finest mesh employed, consisting of over 700,000 grid nodes, is not sufficient to establish grid independent solutions. However, all three schemes appear to converge toward the same solution as the grid spacing approaches zero. The matrix-valued dissipation scheme introduces the least amount of artificial dissipation and should be expected to yield the most accurate solutions on a given mesh. The flux-difference splitting upwind scheme, on the other hand, is more dissipative and, thus, particularly sensitive to grid resolution, but exhibits the best overall convergence characteristics on grids with large aspect ratios.

Author(s):  
Mohammad K Hassanzadeh-Aghdam ◽  
Mohammad J Mahmoodi ◽  
Reza Ansari

A three-dimensional micromechanics-based analytical model is developed to investigate the influence of interphase on the thermo-mechanical properties of three-phase composites. The representative volume element (RVE) of composites is extended to c × r × h cells in three dimensions and the RVE consists of three phases including filler, matrix and interphase. The arrangement state of filler within the matrix materials is assumed to be random with uniform distribution. Fillers are surrounded by the interphase in the whole composite. The effects of interphase such as its thickness and stiffness on the thermo-mechanical properties of composite with various aspect ratios of filler are studied. The results illustrate that while the effects of interphase is significant for composites with randomly distributed spherical particles, it turns to be less effective as the aspect ratio of filler of composite increases. Moreover, the results demonstrate that the effect of interphase on the thermo-mechanical properties of fibrous composites in the transverse direction is more significant than that of fiber composites in the longitudinal direction.


1994 ◽  
Vol 263 ◽  
pp. 207-226 ◽  
Author(s):  
Y. Zeng ◽  
S. Weinbaum

This paper develops a three-dimensional infinite series solution for the Stokes flow through a parallel walled channel which is obstructed by a thin planar barrier with periodically spaced rectangular orifices of arbitrary aspect ratio B’/d’ and spacing D’. Here B’ is the half-height of the channel and d’ is the half-width of the orifice. The problem is motivated by recent electron microscopic studies of the intercellular channel between vascular endothelial cells which show a thin junction strand barrier with discontinuities or breaks whose spacing and width vary with the tissue. The solution for this flow is constructed as a superposition of Hasimoto's (1958) general solution for the two-dimensional flow through a periodic slit array in an infinite plane wall and a new three-dimensional solution which corrects for the top and bottom boundaries. In contrast to the well-known solutions of Sampson (1891) and Hasimoto (1958) for the flow through zero-thickness orifices of circular or elliptic cross-section or periodic slits in an infinite plane wall, which exhibit characteristic viscous velocity profiles, the present bounded solutions undergo a fascinating change in behaviour as the aspect ratio B’/d’ of the orifice opening is increased. For B’/d’ [Lt ] 1 and (D’ –- d’)/B’ of O(1) or greater, which represents a narrow channel, the velocity has a minimum at the orifice centreline, rises sharply near the orifice edges and then experiences a boundary-layer-like correction over a thickness of O(B’) to satisfy no-slip conditions. For B’/d’ of O(1) the profiles are similar to those in a rectangular duct with a maximum on the centreline, whereas for B’/d’ [Gt ] 1, which describes widely separated channel walls, the solution approaches Hasimoto's solution for the periodic infinite-slit array. In the limit (D’ –- d’)/B’ [Lt ] 1, where the width of the intervening barriers is small compared with the channel height, the solutions exhibit the same behaviour as Lee & Fung's (1969) solution for the flow past a single cylinder. The drag on the zero-thickness barriers in this case is nearly the same as for the cylinder for all aspect ratios.


2017 ◽  
Vol 834 ◽  
pp. 385-408 ◽  
Author(s):  
T. S. van den Bremer ◽  
B. R. Sutherland

We examine the wave-induced flow of small-amplitude, quasi-monochromatic, three-dimensional, Boussinesq internal gravity wavepackets in a uniformly stratified ambient. It has been known since Bretherton (J. Fluid Mech., vol. 36 (4), 1969, pp. 785–803) that one-, two- and three-dimensional wavepackets induce qualitatively different flows. Whereas the wave-induced mean flow for compact three-dimensional wavepackets consists of a purely horizontal localized circulation that translates with and around the wavepacket, known as the Bretherton flow, such a flow is prohibited for a two-dimensional wavepacket of infinite spanwise extent, which instead induces a non-local internal wave response that is long compared with the streamwise extent of the wavepacket. One-dimensional (horizontally periodic) wavepackets induce a horizontal, non-divergent unidirectional flow. Through perturbation theory for quasi-monochromatic wavepackets of arbitrary aspect ratio, we predict for which aspect ratios which type of induced mean flow dominates. We compose a regime diagram that delineates whether the induced flow is comparable to that of one-, two- or compact three-dimensional wavepackets. The predictions agree well with the results of fully nonlinear three-dimensional numerical simulations.


Author(s):  
LI CAI ◽  
JIAN-HU FENG ◽  
YU-FENG NIE ◽  
WEN-XIAN XIE

In this paper, we present a third-order central weighted essentially nonoscillatory (CWENO) reconstruction for computations of hyperbolic conservation laws in three space dimensions. Simultaneously, as a Godunov-type central scheme, the CWENO-type central-upwind scheme, i.e., the semi-discrete central-upwind scheme based on our third-order CWENO reconstruction, is developed straightforwardly to solve 3D systems by the so-called componentwise and dimensional-by-dimensional technologies. The high resolution, the efficiency and the nonoscillatory property of the scheme can be verified by solving several numerical experiments.


Author(s):  
David B. Benson ◽  
Tom I.-P. Shih ◽  
David O. Davis

CFD simulations were performed to investigate boundary-layer control through bleed patches in an axisymmetric mixed-compression inlet in which the bleed patches are modeled by two global bleed boundary conditions (BCs). In one bleed BC, the locations of the bleed holes are discerned. In the other bleed BC, each row of bleed holes is modeled as a porous surface, where the number of bleed holes in each row is accounted for by adjusting the discharge coefficient to give the correct bleed rate. Results are presented for the predicted bleed rates, pressure on the cowl and centerbody surfaces, and the flow field. Comparisons were made with available experimental data. Also presented is a method based on one-dimensional isentropic and normal shock solutions to get the flow “started” in CFD simulations of critical flow in mixed-compression inlets. This computational study is based on the ensemble-averaged conservation equations of mass (continuity), momentum (compressible Navier-Stokes), and total energy closed by shear-stress transport (SST) turbulence model, where integration is to the wall. Solutions were generated by a cell-centered finite-volume method that uses third-order accurate flux-difference splitting of Roe with limiters, multigrid acceleration of a diagonalized ADI scheme with local time stepping, and patched/overlapped structured grids.


2019 ◽  
Vol 34 (01n03) ◽  
pp. 2040022
Author(s):  
Qiang Xu ◽  
Yiyi She ◽  
Li Li

Porous electrodes in microfluidic fuel cell (MFC) operate with nonuniform reaction rate. It is intriguing to improve the utilization degree of the porous electrodes. In this work, a three-dimensional computational model is developed for MFC with flow-through porous electrodes. Characteristics of the reaction rate distributions under different electrode geometries are examined. The results show that reaction rate varies noticeably along the electrode width direction, but minimally along the electrode length direction. High reaction rate region locates in the vicinity of the interface between the porous electrode and the middle channel. A relatively high aspect ratio, defined as the ratio of the electrode length to width, is beneficial to improve the utilization degree of the porous electrodes. Yet, concentration losses increase due to the decreased fluid velocity. Considering the cell performance, optimal electrode aspect ratios are derived for the anode and cathode, respectively.


1968 ◽  
Vol 72 (691) ◽  
pp. 618-623 ◽  
Author(s):  
J. H. Argyris ◽  
I. Fried ◽  
D. W. Scharpf

Following on from the LUMINA and HERMES elements discussed in TN’s 11 and 12, this note analyses two further three-dimensional elements outlined in ref. 1. These new elements possess the desirable features of completeness and consequent invariance of the polynomials for the displacement fields. It has been shown in ref. 2 and stated in ref. 1 that third order complete polynomials yielding a second order or parabolic strain distribution fit into a tetrahedron element with 20 nodal points and 60 degrees of freedom, denoted within ASKA as TET 20 (Fig. 1).


Author(s):  
Pengchong Zhang ◽  
Chengzhi Qi ◽  
Hongyuan Fang ◽  
Xu Sun

A three dimensional elasticity analysis for the transverse free vibration characteristics of functionally graded magneto-electro-elastic plates based on the scaled boundary finite element method (SBFEM) incorporated with the precise integration algorithm (PIA) is presented. The material properties of magneto-electro-elastic plates are changing along the in-plane direction with arbitrary mathematical functions. In the proposed methodology, the strategies of only discretizing the in-plane surface and utilizing the two-dimensional spectral elements to construct diagonal coefficient matrices are adopted, which contributes to decreasing the calculation effort. The derivation process begins with the three dimensional governing equations of magneto-electro-elastic materials. Neither the plate mechanical kinematics nor invoking assumptions on the spatial distributions of electric and magnetic quantities are adopted. Built upon the introduced scaled boundary coordinates, the principle of virtual work and the technique of dual vectors, a first order ordinary differential SBFEM matrix equation for the in-plane functionally graded magneto-electro-elastic plates is obtained. Its general solution is analytically denoted as the matrix exponent. To improve the computation accuracy of the matrix exponent, the PIA is utilized to form the stiffness matrix. By virtue of the kinetic energy technique, it is convenient to construct the mass matrix of the in-plane functionally graded magneto-electro-elastic plates based on the SBFEM for the first time. Finally, comparisons of flexural frequency parameters with those from exact solutions and other numerical methods are provided. The accuracy, effectiveness, and versatility of the employed technique are validated. Moreover, additional numerical exercises are conducted to exhibit the influences of boundary conditions, material gradation functions, and aspect ratios on the free vibration behaviors of in-plane functionally graded magneto-electro-elastic rectangular, circular, and perforated plates.


2003 ◽  
Vol 70 ◽  
pp. 201-212 ◽  
Author(s):  
Hideaki Nagase ◽  
Keith Brew

The tissue inhibitors of metalloproteinases (TIMPs) are endogenous inhibitors of the matrix metalloproteinases (MMPs), enzymes that play central roles in the degradation of extracellular matrix components. The balance between MMPs and TIMPs is important in the maintenance of tissues, and its disruption affects tissue homoeostasis. Four related TIMPs (TIMP-1 to TIMP-4) can each form a complex with MMPs in a 1:1 stoichiometry with high affinity, but their inhibitory activities towards different MMPs are not particularly selective. The three-dimensional structures of TIMP-MMP complexes reveal that TIMPs have an extended ridge structure that slots into the active site of MMPs. Mutation of three separate residues in the ridge, at positions 2, 4 and 68 in the amino acid sequence of the N-terminal inhibitory domain of TIMP-1 (N-TIMP-1), separately and in combination has produced N-TIMP-1 variants with higher binding affinity and specificity for individual MMPs. TIMP-3 is unique in that it inhibits not only MMPs, but also several ADAM (a disintegrin and metalloproteinase) and ADAMTS (ADAM with thrombospondin motifs) metalloproteinases. Inhibition of the latter groups of metalloproteinases, as exemplified with ADAMTS-4 (aggrecanase 1), requires additional structural elements in TIMP-3 that have not yet been identified. Knowledge of the structural basis of the inhibitory action of TIMPs will facilitate the design of selective TIMP variants for investigating the biological roles of specific MMPs and for developing therapeutic interventions for MMP-associated diseases.


2012 ◽  
Vol 9 (1) ◽  
pp. 142-146
Author(s):  
O.A. Solnyshkina

In this work the 3D dynamics of two immiscible liquids in unbounded domain at low Reynolds numbers is considered. The numerical method is based on the boundary element method, which is very efficient for simulation of the three-dimensional problems in infinite domains. To accelerate calculations and increase the problem size, a heterogeneous approach to parallelization of the computations on the central (CPU) and graphics (GPU) processors is applied. To accelerate the iterative solver (GMRES) and overcome the limitations associated with the size of the memory of the computation system, the software component of the matrix-vector product


Sign in / Sign up

Export Citation Format

Share Document