A dynamic coupling model for hybrid atomistic-continuum computations

A dynamic coupling model is developed for a hybrid atomistic-continuum computation in micro- and nano-fluidics. In the hybrid atomistic-continuum computation, a molecular dynamics (MD) simulation is utilized in one region where the continuum assumption breaks down and the Navier-Stokes (NS) equations are used in another region where the continuum assumption holds. In the overlapping part of these two regions, a constrained particle dynamics is needed to couple the MD simulation and the NS equations. The currently existing coupling models for the constrained particle dynamics have a coupling parameter, which has to be empirically determined. In the present work, a novel dynamic coupling model is introduced where the coupling parameter can be calculated as the computation progresses rather than inputing a priori. The dynamic coupling model is based on the momentum constraint and exhibits a correct relaxation rate. The results from the hybrid simulation on the Couette flow and the Stokes flow are in good agreement with the data from the full MD simulation and the solutions of the NS equations, respectively. (c) 2007 Elsevier Ltd. All rights reserved.

Multiscale Method for Turbulent-Flow Computations

A main method of predicting turbulent flows is to solve LES equations, which was called traditional LES method. The traditional LES method solves the motions of large eddies of size larger than filtering scale An while modeling unresolved scales less than Delta_n. Hughes et al argued that many shortcomings of the traditional LES approaches were associated with their inabilities to successfully differentiate between large and small scales. One may guess that a priori scale-separation would be better, because it can predict scale-interaction well compared with posteriori scale-separation. To this end, a multi-scale method was suggested to perform scale-separation computation. The primary contents of the multiscale method are l) A space average is used to differentiate scale. 2) The basic equations include the large scale equations and fluctuation equations. 3) The large-scale equations and fluctuation equations are coupled through turbulent stress terms. We use the multiscale equations of n=2, i.e., the large and small scale (LSS) equations, to simulate 3-D evolutions of a channel flow and a planar mixing layer flow Some interesting results are given.

Combination Of Cfd And Csd Packages For Fluid-Structure Interaction

In this article the UDF script file in the Fluent software was rewritten as the "connecting file" for the Fluent and the ANSYS/ABAQUS in order that the joined file can be used to do aero-elastic computations. In this way the fluid field is computed by solving the Navier-Stokes equations and the structure movement is integrated by the dynamics directly. An analysis of the computed results shows that this coupled method designed for simulating aero-elastic systems is workable and can be used for the other fluid-structure interaction problems.

Fast Computing For Lurr Of Earthquake Prediction

The LURR theory is a new approach for earthquake prediction, which achieves good results in earthquake prediction within the China mainland and regions in America, Japan and Australia. However, the expansion of the prediction region leads to the refinement of its longitude and latitude, and the increase of the time period. This requires increasingly more computations, and the volume of data reaches the order of GB, which will be very difficult for a single CPU. In this paper, a new method was introduced to solve this problem. Adopting the technology of domain decomposition and parallelizing using MPI, we developed a new parallel tempo-spatial scanning program.

A Note On The Numerical Simulations Of Flow Past A Wavy Square-Section Cylinder

The flow past a square-section cylinder with a geometric disturbance is investigated by numerical simulations. The extra terms, due to the introduction of mapping transformation simulating the effect of disturbance into the transformed Navier-Stokes equations, are correctly derived, and the incorrect ones in the previous literature are pointed out and analyzed. Furthermore, the relationship between the vorticity, especially on the cylinder surface, and the disturbance is derived and explained theoretically. The computations are performed at two Reynolds numbers of 100 and 180 and three amplitudes of waviness of 0.006, 0.025 and 0.167 with another aim to explore the effects of different Reynolds numbers and disturbance on the vortex dynamics in the wake and forces on the body. Numerical results have shown that, at the mild waviness of 0.025, the Karman vortex shedding is suppressed completely for Re = 100, while the forced vortex dislocation is appeared in the near wake at the Reynolds number of 180. The drag reduction is up to 21.6% at Re = 100 and 25.7% at Re = 180 for the high waviness of 0.167 compared with the non-wavy cylinder. The lift and the Strouhal number varied with different Reynolds numbers and the wave steepness are also obtained.

An investigation of electromagnetic wave propagation in plasma by shock tube

This paper presents the electromagnetic wave propagation characteristics in plasma and the attenuation coefficients of the microwave in terms of the parameters n(e), v, w, L, w(b). The phi800 mm high temperature shock tube has been used to produce a uniform plasma. In order to get the attenuation of the electromagnetic wave through the plasma behind a shock wave, the microwave transmission has been used to measure the relative change of the wave power. The working frequency is f = (2 similar to 35) GHz (w = 2pif, wave length lambda = 15 cm similar to 8 mm). The electron density in the plasma is n(e) = (3 x 10(10) similar to 1 x 10(14)) cm(-3). The collision frequency v = (1 x 10(8) similar to 6 x 10(10)) Hz. The thickness of the plasma layer L = (2 similar to 80) cm. The electron circular frequency w(b) = eB(0)/m(e), magnetic flux density B-0 = (0 similar to 0.84) T. The experimental results show that when the plasma layer is thick (such as L/lambda greater than or equal to 10), the correlation between the attenuation coefficients of the electromagnetic waves and the parameters n(e), v, w, L determined from the measurements are in good agreement with the theoretical predictions of electromagnetic wave propagations in the uniform infinite plasma. When the plasma layer is thin (such as when both L and lambda are of the same order), the theoretical results are only in a qualitative agreement with the experimental observations in the present parameter range, but the formula of the electromagnetic wave propagation theory in an uniform infinite plasma can not be used for quantitative computations of the correlation between the attenuation coefficients and the parameters n(e), v, w, L. In fact, if w < w(p), v(2) much less than w(2), the power attenuations K of the electromagnetic waves obtained from the measurements in the thin-layer plasma are much smaller than those of the theoretical predictions. On the other hand, if w > w(p), v(2) much less than w(2) (just v approximate to f), the measurements are much larger than the theoretical results. Also, we have measured the electromagnetic wave power attenuation value under the magnetic field and without a magnetic field. The result indicates that the value measured under the magnetic field shows a distinct improvement.

An improvement on the theoretical model for the large-scale coherent structures in the outer region of a turbulent boundary layer

According to the experimental results, there exist large-scale coherent structures in the outer region of a turbulent boundary layer, which have been studied by many authors.As experimental results, Antonia (1990) showed the phase- aver aged streamlines and isovorticity lines of the large-scale coherent structures in a turbulent boundary layer for different Reynolds numbers. Based on the hydrodynamic stability theory, the 2-D theoretical model for the large-scale structures was proposed by Luo and Zhou, in which the eddy viscosity was defined as a complex function of the position in the normal direction. The theoretical results showed in ref. were in agreement with those in ref. However, there were two problems in the results. One is that in the experimental results, there were divergent focuses between two saddle points in the streamlines, but in the theoretical results, there were centers. The other is that the stretched parts of the isovorticity lines appear at the location of centers in the theoretical results, while in the experimental results they located somewhere between the focuses and saddle points. The reason is that the computations were based on a 2-D model.

Plastic stress-strain fields and blunting effects during large deformations

Plastic stress-strain fields of two types of steel specimens loaded to large deformations are studied. Computational results demonstrate that, owing to the fact that the hardening exponent of the material varies as strain enlarges and the blunting of the crack tip, the well known HRR stress field in the plane strain model can only hold for the stage of a small plastic strain. Plastic dilatancy is shown to have substantial effects on strain distributions and blunting. To justify the constitutive equations used for analysis and to check the precision of computations, the load-deflection of a three-point bend beam and the load-elongation of an axisymmetric bar notched by a V-shaped cut were tested and recorded. The computed curves are in good accordance with experimental data.

Crack tip behavior and crack-propagation in ductile materials

Dilatational plastic equations, which can include the effects of ductile damage, are derived based on the equivalency in expressions for dissipated plastic work. Void damage developed internally at the large-strain stage is represented by an effective continuum being strain-softened and plastically dilated. Accumulation of this local damage leads to progressive failure in materials. With regard to this microstructural background, the constitutive parameters included for characterizing material behaviour have the sense of internal variables. They are not able to be determined explicitly by macroscopic testing but rather through computer simulation of experimental curves and data. Application of this constitutive model to mode-I cracking examples demonstrates that a huge strain concentration accompanied by a substantial drop of stress does occur near the crack tip. Eventually, crack propagation is simulated by using finite elements in computations. Two numerical examples show good accordance with experimental data. The whole procedure of study serves as a justification of the constitutive formulation proposed in the text.

A Multiscale-Domain Decomposition Method for High Reynolds Number Flows

The high Reynolds number flow contains a wide range of length and time scales, and the flow domain can be divided into several sub-domains with different characteristic scales. In some sub-domains, the viscosity dissipation scale can only be considered in a certain direction; in some sub-domains, the viscosity dissipation scales need to be considered in all directions; in some sub-domains, the viscosity dissipation scales are unnecessary to be considered at all. For laminar boundary layer region, the characteristic length scales in the streamwise and normal directions are L and L Re-1/ 2 , respectively. The characteristic length scale and the velocity scale in the outer region of the boundary layer are L and U, respectively. In the neighborhood region of the separated point, the length scale l<computations for high Reynolds number flows, an idea of solving the conservation equations for discrete cells was proposed and named the discrete fluid dynamics (DFD) algorithm. Analysis shows that the basic conservative equations for discrete cells are the Euler equations, NS- and diffusion parabolized (DP) NS equations. In this paper, a new multiscale-domain decomposition method is developed for the high Reynolds number flow. First, the whole domain is decomposed to different sub-domains with the different characteristic scales. Then the different dominant equation of all sub-domains is defined according to the diffusion parabolized (DP) theory of viscous flow. Finally these different equations are solved simultaneously in whole computational region. For numerical tests of high Reynolds numerical flows, two-dimensional supersonic flows over rearward and frontward steps as well as an interaction flow between shock wave and boundary layer were solved numerically. The pressure distributions and local coefficients of skin friction on the wall are given. The numerical results obtained by the multiscale-domain decomposition algorithm are well agreement with those by NS equations. Comparing with the usual method of solving the Navier-Stokes equations in the whole flow, under the same numerical accuracy, the present multiscale domain decomposition method decreases CPU consuming about 20% and reflects the physical mechanism of practical flow more accurately.

Numerical Analyses of Transonic Fluid/Structure interaction in Aircraft Engineering

Through the coupling between aerodynamic and structural governing equations, a fully implicit multiblock aeroelastic solver was developed for transonic fluid/stricture interaction. The Navier-Stokes fluid equations are solved based on LU-SGS (lower-upper symmetric Gauss-Seidel) Time-marching subiteration scheme and HLLEW (Harten-Lax-van Leer-Einfeldt-Wada) spacing discretization scheme and the same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Transfinite interpolation (TFI) is used for the grid deformation of blocks neighboring the flexible surfaces. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between fluid and structure. The developed code was fort validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. In the subsonic and transonic range, the calculated flutter speeds and frequencies agree well with experimental data, however, in the supersonic range, the present calculation overpredicts the experimental flutter points similar to other computations. Then the flutter character of a complete aircraft configuration is analyzed through the calculation of the change of structural stiffness. Finally, the phenomenon of aileron buzz is simulated for the weakened model of a supersonic transport wing/body model at Mach numbers of 0.98 and l.05. The calculated unsteady flow shows, on the upper surface, the shock wave becomes stronger as the aileron deflects downward, and the flow behaves just contrary on the lower surface of the wing. Corresponding to general theoretical analysis, the flow instability referred to as aileron buzz is induced by a stronger shock alternately moving on the upper and lower surfaces of wing. For the rigid structural model, the flow is stable at all calculated Mach numbers as observed in experiment

Influence of Contact Angle and Tube Size on Capillary-Driven Flow Under Microgravity

The influence of contact angle and tube radius on the capillary-driven flow for circular cylindrical tubes is studied systematically by microgravity experiments using the drop tower. Experimental results show that the velocity of the capillary flow decreases monotonically with an increase in the contact angle. However, the time-evolution of the velocity of the capillary flow is different for different sized tubes. At the beginning of the microgravity period, the capillary flow in a thinner tube moves faster than that in a thicker tube, and then the latter overtakes the former. Therefore, there is an intersection between the curves of meniscus velocity vs microgravity time for two differently sized tubes. In addition, for two given sized tubes this intersection is delayed when the contact angle increases. The experimental results are analyzed theoretically and also supported by numerical computations.