HOMOGENIZATION OF A HYPERBOLIC EQUATION WITH HIGHLY CONTRASTING DIFFUSIVITY COEFFICIENTS

We study a hyperbolic problem in the framework of periodic homogenization assuming a high contrast between the diffusivity coefficients of the two components M-epsilon and B-epsilon of the heterogeneous medium. There are three regimes depending on the ratio between the size of the period and the amplitude a, of the diffusivity in B-epsilon. For the critical regime alpha(epsilon) similar or equal to epsilon, the limit problem is a strongly coupled system involving both the macroscopic and the microscopic variables. We also include the results in the non critical case.

HOMOGENIZATION OF A HYPERBOLIC EQUATION WITH HIGHLY CONTRASTING DIFFUSIVITY COEFFICIENTS

We study a hyperbolic problem in the framework of periodic homogenization assuming a high contrast between the diffusivity coefficients of the two components M-epsilon and B-epsilon of the heterogeneous medium. There are three regimes depending on the ratio between the size of the period and the amplitude a, of the diffusivity in B-epsilon. For the critical regime alpha(epsilon) similar or equal to epsilon, the limit problem is a strongly coupled system involving both the macroscopic and the microscopic variables. We also include the results in the non critical case.

High-Resolution Finite Compact Difference Schemes for Hyperbolic Conservation Laws

A finite compact (FC) difference scheme requiring only bi-diagonal matrix inversion is proposed by using the known high-resolution flux. Introducing TVD or ENO limiters in the numerical flux, several high-resolution FC-schemes of hyperbolic conservation law are developed, including the FC-TVD, third-order FC-ENO and fifth-order FC-ENO schemes. Boundary conditions formulated need only one unknown variable for third-order FC-ENO scheme and two unknown variables for fifth-order FC-ENO scheme. Numerical test results of the proposed FC-scheme were compared with traditional TVD, ENO and WENO schemes to demonstrate its high-order accuracy and high-resolution.

Coherent Structures in Transition of a Flat-Plate Boundary Layer at Ma=0.7

Direct numerical simulation (DNS) of a spatially evolving flat-plate boundary layer transition process at free stream Mach number 0.7 is performed. Tollmien-Schlichting (T-S) waves are added on the inlet boundary as the disturbances before transition. Typical coherent structures in the transition process are investigated based on the second invariant of velocity gradient tensor. The instantaneous shear stress and the mean velocity profile in the transition region are studied. In our view, the fact that the peak value of shear stress in the stress concentration area increases and exceeds a threshold value during the later stage of the transition process plays an important role in the laminar breakdown process.

On a Degenerate Parabolic Problem in Holder Spaces

In this paper, we study some degenerate parabolic equation with Cauchy-Dirichlet boundary conditions. This problem is considered in little Holder spaces. The optimal regularity of the solution v is obtained and is specified in terms of those of the second member when some conditions upon the Holder exponent with respect to the degeneracy are satisfied. The proofs mainly use the sum theory of linear operators with or without density of domains and the results of smoothness obtained in the study of some abstract linear differential equations of elliptic type.

Flow-Induced Vibrations Of A Cylinder With Two Degrees Of Freedom Near Rigid Plane Boundary

Based on similarity analyses, the flow-induced vibrations of a near-wall cylinder with 2 degrees of freedom are investigated experimentally by employing a hydroelastic apparatus in conjunction with a flume. The cylinder's vibration amplitude, vibration frequency and vortex shedding frequency were measured and analyzed. The effects of gap-to-diameter ratio (e,ID) upon the vibration responses are further investigated. The experimental results indicate that, when the reduced velocity (Vr) is small (e.g. Vr = 1.2 similar to 2.6), only streamwise vibration occurs, and its frequency is quite close to its natural frequency in still water. When increasing Vr (e.g. Vr > 3.4), both streamwise and transverse vibrations of the near-wall cylinder may occur. In the examined range of gap-to-diameter ratio (0.42 < e(0)/D < 2.68), 2 vibration stages (in terms of Vr) of streamwise vibrations usually exist: First Streamwise Vibration (FSV) and Second Streamwise Vibration (SSV). In the SSV stage, the vortex shedding frequency may either undergo a jump to that of the streamwise vibration, or stay consistent with that of the transverse vibration. The amplitudes of transverse vibration are usually much larger than those of streamwise vibration for the same value of e(0)/D. The maximum amplitudes of both streamwise and transverse vibration get larger with the increase of e(0)/D (0.42 < e(0)/D < 2.68).

Direct Numerical Simulation Of Hypersonic Boundary-Layer Transition Over A Blunt Cone

Direct numerical simulation of transition How over a blunt cone with a freestream Mach number of 6, Reynolds number of 10,000 based on the nose radius, and a 1-deg angle of attack is performed by using a seventh-order weighted essentially nonoscillatory scheme for the convection terms of the Navier-Stokes equations, together with an eighth-order central finite difference scheme for the viscous terms. The wall blow-and-suction perturbations, including random perturbation and multifrequency perturbation, are used to trigger the transition. The maximum amplitude of the wall-normal velocity disturbance is set to 1% of the freestream velocity. The obtained transition locations on the cone surface agree well with each other far both cases. Transition onset is located at about 500 times the nose radius in the leeward section and 750 times the nose radius in the windward section. The frequency spectrum of velocity and pressure fluctuations at different streamwise locations are analyzed and compared with the linear stability theory. The second-mode disturbance wave is deemed to be the dominating disturbance because the growth rate of the second mode is much higher than the first mode. The reason why transition in the leeward section occurs earlier than that in the windward section is analyzed. It is not because of higher local growth rate of disturbance waves in the leeward section, but because the growth start location of the dominating second-mode wave in the leeward section is much earlier than that in the windward section.

Receptivity To Free-Stream Disturbance Waves For Hypersonic Flow Over A Blunt Cone

A high-order shock-fitting finite difference scheme is studied and used to do direction numerical simulation (DNS) of hypersonic unsteady flow over a blunt cone with fast acoustic waves in the free stream, and the receptivity problem in the blunt cone hypersonic boundary layers is studied. The results show that the acoustic waves are the strongest disturbance in the blunt cone hypersonic boundary layers. The wave modes of disturbance in the blunt cone boundary layers are first, second, and third modes which are generated and propagated downstream along the wall. The results also show that as the frequency decreases, the amplitudes of wave modes of disturbance increase, but there is a critical value. When frequency is over the critial value, the amplitudes decrease. Because of the discontinuity of curvature along the blunt cone body, the maximum amplitudes as a function of frequencies are not monotone.

Computer Simulation Of Random Sphere Packing In An Arbitrarily Shaped Container

Most simulations of random sphere packing concern a cubic or cylindric container with periodic boundary, containers of other shapes are rarely studied. In this paper, a new relaxation algorithm with pre-expanding procedure for random sphere packing in an arbitrarily shaped container is presented. Boundaries of the container are simulated by overlapping spheres which covers the boundary surface of the container. We find 0.4 similar to 0.6 of the overlap rate is a proper value for boundary spheres. The algorithm begins with a random distribution of small internal spheres. Then the expansion and relaxation procedures are performed alternately to increase the packing density. The pre-expanding procedure stops when the packing density of internal spheres reaches a preset value. Following the pre-expanding procedure, the relaxation and shrinking iterations are carried out alternately to reduce the overlaps of internal spheres. The pre-expanding procedure avoids the overflow problem and gives a uniform distribution of initial spheres. Efficiency of the algorithm is increased with the cubic cell background system and double link data structure. Examples show the packing results agree well with both computational and experimental results. Packing density about 0.63 is obtained by the algorithm for random sphere packing in containers of various shapes.

A contact crack problem in an infinite plate

The problem of an infinite plate with crack of length 2a loaded by the remote tensile stress P and a pair of concentrated forces Q is discussed. The value of the force Q for the initial contact of crack face is investigated and the contact length elevated, while the Q force increases. The problem is solved assuming that the stress intensity factor vanishes at the end point of the contact portion. By the Fredholm integral equation for the multiple cracks, the reduction of stress intensity factor due to Q is found. (C) 1999 Elsevier Science Ltd. All rights reserved.

A BIEM optimization method for fracture dynamics inverse problem

In the present paper, based on the theory of dynamic boundary integral equation, an optimization method for crack identification is set up in the Laplace frequency space, where the direct problem is solved by the author's new type boundary integral equations and a method for choosing the high sensitive frequency region is proposed. The results show that the method proposed is successful in using the information of boundary elastic wave and overcoming the ill-posed difficulties on solution, and helpful to improve the identification precision.

Molecular dynamics studies on the dislocation gliding near a tilt boundary

The gliding behavior of edge dislocation near a grain boundary(QB) in copper under pure shear stresses is simulated by using molecular dynamics(MD) method. Many-body potential incorporating the embedded atom method (EAM) is used. The critical shear stresses for a single disocation to pass across GB surface are obtained at values of sigma(c)=23MPa similar to 68 MPa and 137 MPa similar to 274 MPa for Sigma=165 small angle tilt GB at 300 K and 20 K, respectively. The first result agrees with the experimental yield stress sigma(y)(=42 MPa) quite well. It suggests that there might be one of the reasons of initial plastic yielding caused by single dislocation gliding across GB. In addition, there might be possibility to obtain yield strength from microscopic analysis. Moreover, the experimental value of sigma(y) at low temperature is generally higher than that at room temperature. So, these results are in conformity qualitatively with experimental fact. On the other hand, the Sigma=25 GB is too strong an obstacle to the dislocation. In this case, a dislocation is able to pass across GB under relatively low stress only when it is driven by other dislocations. This is taken to mean that dislocation pile-up must be built up in front of this kind of GB, if this GB may take effect on the process of plastic deformation.

Mechanical modeling of grain boundary sliding of polycrystals

Using a variational method, a general three-dimensional solution to the problem of a sliding spherical inclusion embedded in an infinite anisotropic medium is presented in this paper. The inclusion itself is also a general anisotropic elastic medium. The interface is treated as a thin interface layer with interphase anisotropic properties. The displacements in the matrix and the inclusion are expressed as polynomial series of the cartesian coordinate components. Using the virtual work principle, a set of linear algebraic equations about unknown coefficients are obtained. Then the general sliding spherical inclusion problem is accurately solved. Based on this solution, a self-consistent method for sliding polycrystals is proposed. Combining this with a two-dimensional model of an aggregate polycrystal, a systematic analysis of the mechanical behaviour of sliding polycrystals is given in detail. Numerical results are given to show the significant effect of grain boundary sliding on the overall mechanical properties of aggregate polycrystals.

Receptivity to free-stream disturbance waves for blunt cone axial symmetry hypersonic boundary layer

Based on high-order compact upwind scheme, a high-order shock-fitting finite difference scheme is studied to simulate the generation of boundary layer disturbance waves due to free-stream waves. Both steady and unsteady flow solutions of the receptivity problem are obtained by resolving the full Navier-Stokes equations. The interactions of bow-shock and free-stream disturbance are researched. Direct numerical simulation (DNS) of receptivity to free-stream disturbances for blunt cone hypersonic boundary layers is performed.

Consumption and portfolio rules whit stochastic hyperbolic discounting

We extend the classic Merton (1969, 1971) problem that investigates the joint consumption-savings and portfolio-selection problem under capital risk by assuming sophisticated but time-inconsistent agents. We introduce stochastic hyperbolic preferences as in Harris and Laibson (2013) and find closed-form solutions for Merton's optimal consumption and portfolio selection problem in continuous time. We find that the portfolio rule remains identical to the time-consistent solution with power utility and no borrowing constraints. However,the marginal propensity to consume out of wealth is unambiguously greater than the time-consistent, exponential case and,importantly, it is also more responsive to changes in risk. These results suggest that hyperbolic discounting with sophisticated agents offers promise for contributing to explaining important aspects of asset market data.