Reliable Validation Based on Optical Flow Visualization for CFD Simulations

A reliable validation based on the optical flow visualization for numerical simulations of complex flowfields is addressed in this paper. Several test cases, including two-dimensional, axisymmetric and three-dimensional flowfields, were presented to demonstrate the effectiveness of the validation and gain credibility of numerical solutions of complex flowfields. In the validation, images of these flowfields were constructed from numerical results based on the principle of the optical flow visualization, and compared directly with experimental interferograms. Because both experimental and numerical results are of identical physical representation, the agreement between them can be evaluated effectively by examining flow structures as well as checking discrepancies in density. The study shows that the reliable validation can be achieved by using the direct comparison between numerical and experiment results without any loss of accuracy in either of them.

Flows with Inertia in a Three-Dimensional Random Fiber Network

A fiber web is modeled as a three-dimensional random cylindrical fiber network. Nonlinear behavior of fluid flowing through the fiber network is numerically simulated by using the lattice Boltzmann (LB) method. A nonlinear relationship between the friction factor and the modified Reynolds number is clearly observed and analyzed by using the Fochheimer equation, which includes the quadratic term of velocity. We obtain a transition from linear to nonlinear region when the Reynolds numbers are sufficiently high, reflecting the inertial effect of the flows. The simulated permeability of such fiber network has relatively good agreement with the experimental results and finite element simulations.

Curvature and wrinkling of premixed flame Kernels comparison between planar laser induced fluorescence (PLIF) and 3D direct numerical simulations (DNS)

A direct comparison between time resolved PLIF measurements of OH and two dimensional slices from a full three dimensional DNS data set of turbulent premixed flame kernels in lean methane/air mixture was presented. The local flame structure and the degree of flame wrinkling were examined in response to differing turbulence intensities and turbulent Reynolds numbers. Simulations were performed using the SEGA DNS code, which is based on the solution of the compressible Navier Stokes, species, and energy equations for a lean hydrocarbon mixture. For the OH PLIF measurements, a cluster of four Nd:YAG laser was fired sequentially at high repetition rates and used to pump a dye laser. The frequency doubled laser beam was formed into a sheet of 40 mm height using a cylindrical telescope. The combination of PLIF and DNS has been demonstrated as a powerful tool for flame analysis. This research will form the basis for the development of sub-grid-scale (SGS) models for LES of lean-premixed combustion systems such as gas turbines. This is an abstract of a paper presented at the 30th International Symposium on Combustion (Chicago, IL 7/25-30/2004).

Low-Dimensional Dynamical Systems for a Wake-Type Shear Flow with Vortex Dislocations

Low-dimensional systems are constructed to investigate dynamics of vortex dislocations in a wake-type shear flow. High-resolution direct numerical simulations are employed to obtain flow snapshots from which the most energetic modes are extracted using proper orthogonal decomposition (POD). The first 10 modes are classified into two groups. One represents the general characteristics of two-dimensional wake-type shear flow, and the other is related to the three-dimensional properties or non-uniform characteristics along the span. Vortex dislocations are generated by these two kinds of coherent structures. The results from the first 20 three-dimensional POD modes show that the low- dimensional systems have captured the basic properties of the wake-type shear flow with vortex dislocation, such as two incommensurable frequencies and their beat frequency.

Simulations of Mechanical Behavior of Polycrystalline Copper with Nano-Twins

Mechanical behavior and microstructure evolution of polycrystalline copper with nano-twins were investigated in the present work by finite element simulations. The fracture of grain boundaries are described by a cohesive interface constitutive model based on the strain gradient plasticity theory. A systematic study of the strength and ductility for different grain sizes and twin lamellae distributions is performed. The results show that the material strength and ductility strongly depend on the grain size and the distribution of twin lamellae microstructures in the polycrystalline copper.

Determination Ofinterfacial Mechanical Parameters For An Al/Epoxy/Al2O3 System Byusing Peel Test Simulations

Peel test measurements and simulations of the interfacial mechanical parameters for the Al/Epoxy/Al2O3 system are performed in the present investigation. A series of Al film thicknesses between 20 and 250 microns and three peel angles of 90, 135 and 180 degrees are considered. Two types of epoxy adhesives are adopted to obtain both strong and weak interface adhesions. A finite element model with cohesive zone elements is used to identify the interfacial parameters and simulate the peel test process. By simulating and recording normal stress near the crack tip, the separation strength is obtained. Furthermore, the cohesive energy is identified by comparing the simulated steady-state peel force and the experimental result. It is found from the research that both the cohesive energy and the separation strength can be taken as the intrinsic interfacial parameters which are dependent on the thickness of the adhesive layer and independent of the film thickness and peel angle.

Comparative Study Of The Influence Of Natural Convection On Directional Solidification Of Al-3.5 Wt% Ni And Al-7 Wt% Si Alloys

We present numerical simulations of thermosolutal convection for directional solidification of Al-3.5 wt% Ni and Al-7 wt% Si. Numerical results predict that fragmentation of dendrite arms resulting from dissolution could be favored in Al-7 wt% Si, but not in Al-3.5 wt% Ni. Corresponding experiments are in qualitative agreement with the numerical predictions. Distinguishing the two fragmentation mechanisms, namely dissolution and remelting, is critical during experiments on earth, when fluid flow is dominant. (C) 2007 COSPAR. Published by Elsevier Ltd. All rights reserved.

Reconsideration On The Role Of The Specific Heat Ratio In Arrhenius Law Applications

Arrhenius law implicates that only those molecules which possess the internal energy greater than the activation energy E-a can react. However, the internal energy will not be proportional to the gas temperature if the specific heat ratio gamma and the gas constant R vary during chemical reaction processes. The varying gamma may affect significantly the chemical reaction rate calculated with the Arrhenius law under the constant gamma assumption, which has been widely accepted in detonation and combustion simulations for many years. In this paper, the roles of variable gamma and R in Arrhenius law applications are reconsidered, and their effects on the chemical reaction rate are demonstrated by simulating one-dimensional C-J and two-dimensional cellular detonations. A new overall one-step detonation model with variable gamma and R is proposed to improve the Arrhenius law. Numerical experiments demonstrate that this improved Arrhenius law works well in predicting detonation phenomena with the numerical results being in good agreement with experimental data.

Two-Dimensional Kinetics Of Beta(2)-Integrin And Icam-1 Bindings Between Neutrophils And Melanoma Cells In A Shear Flow

Cell adhesion, mediated by specific receptor-ligand interactions, plays an important role in biological processes such as tumor metastasis and inflammatory cascade. For example, interactions between beta(2)-integrin ( lymphocyte function-associated antigen-1 and/or Mac-1) on polymorphonuclear neutrophils (PMNs) and ICAM-1 on melanoma cells initiate the bindings of melanoma cells to PMNs within the tumor microenvironment in blood flow, which in turn activate PMN-melanoma cell aggregation in a near-wall region of the vascular endothelium, therefore enhancing subsequent extravasation of melanoma cells in the microcirculations. Kinetics of integrin-ligand bindings in a shear flow is the determinant of such a process, which has not been well understood. In the present study, interactions of PMNs with WM9 melanoma cells were investigated to quantify the kinetics of beta(2)-integrin and ICAM-1 bindings using a cone-plate viscometer that generates a linear shear flow combined with a two-color flow cytometry technique. Aggregation fractions exhibited a transition phase where it first increased before 60 s and then decreased with shear durations. Melanoma-PMN aggregation was also found to be inversely correlated with the shear rate. A previously developed probabilistic model was modified to predict the time dependence of aggregation fractions at different shear rates and medium viscosities. Kinetic parameters of beta(2)-integrin and ICAM-1 bindings were obtained by individual or global fittings, which were comparable to respectively published values. These findings provide new quantitative understanding of the biophysical basis of leukocyte-tumor cell interactions mediated by specific receptor-ligand interactions under shear flow conditions.

Numerical exploration of new features on three-dimensional transition of a cylinder wake

The three-dimensional transition of the wake flow behind a circular cylinder is studied in detail by direct numerical simulations using 3D incompressible N-S equations for Reynolds number ranging from 200 to 300. New features and vortex dynamics of the 3D transition of the wake are found and investigated. At Re = 200, the flow pattern is characterized by mode A instability. However, the spanwise characteristic length of the cylinder determines the transition features. Particularly for the specific spanwise characteristic length linear stable mode may dominate the wake in place of mode A and determine the spanwise phase difference of the primary vortices shedding. At Re = 250 and 300 it is found that the streamwise vortices evolve into a new type of mode - "dual vortex pair mode" downstream. The streamwise vortex structures switch among mode A, mode B and dual vortex pair mode from near wake to downstream wake. At Re = 250, an independent low frequency f(m) in addition to the vortex shedding frequency f(s) is identified. Frequency coupling between f(m) and f(s) occurs. These result in the irregularity of the temporal signals and become a key feature in the transition of the wake. Based on the formation analysis of the streamwise vorticity in the vicinity of cylinder, it is suggested that mode A is caused by the emergence of the spanwise velocity due to three dimensionality of the incoming flow past the cylinder. Energy distribution on various wave numbers and the frequency variation in the wake are also described.

Molecular dynamics studies on the local disordering of σ3 and σ11 grain boundaries in aluminium bicrystal

Molecular dynamics (MD) simulations using Morse interaction potential are performed in studies of [110] symmetrical tilt grain boundary (GB) structures with mis-orientation angles 50.5 degrees(Sigma 11), 129.5 degrees(Sigma 11), 70.5 degrees(Sigma 3) and 109.5 degrees(Sigma 3) at various tempratures. The GB structures are found to start local disordering at about 0.5T(m)(T-m is the melting point of aluminium) for 50.5 degrees(Sigma 11), 0.32T(m) for 129.5 degrees(Sigma 11) and 0.38T(m) for 70.5 degrees(Sigma 3), respectively. These results agree with conclusions deduced from the anelastic measurements. But, for twin-boundary structure 109.5 degrees(Sigma 3), this disordering has not been found even when temperature increases up to 0.9T(m).

Character of simplified Navier-Stokes (SNS) equations near separation point for two-dimensional laminar flow over a flat plate

It is proved that the simplified Navier-Stokes (SNS) equations presented by Gao Zhi[1], Davis and Golowachof-Kuzbmin-Popof (GKP)[3] are respectively regular and singular near a separation point for a two-dimensional laminar flow over a flat plate. The order of the algebraic singularity of Davis and GKP equation[2,3] near the separation point is indicated. A comparison among the classical boundary layer (CBL) equations, Davis and GKP equations, Gao Zhi equations and the complete Navier-Stokes (NS) equations near the separation point is given.

Nonequilibrium statistical mechanics of one-dimensional turbulence

The method of statistical mechanics is applied to the study of the one-dimensional model of turbulence proposed in an earlier paper. The closure problem is solved by the variational approach which has been developed for the three-dimensional case, yielding two integral equations for two unknown functions. By solving the two integral equations, the Kolmogorov k−5/3 law is derived and the (one-dimensional) Kolmogorov constant Ko is evaluated, obtaining Ko=0.55, which is in good agreement with the result of numerical experiments on one-dimensional turbulence.