A geometry model for tortuosity of flow path in porous media

A simple geometry model for tortuosity of flow path in porous media is proposed based on the assumption that some particles in a porous medium are unrestrictedly overlapped and the others are not. The proposed model is expressed as a function of porosity and there is no empirical constant in this model. The model predictions are compared with those from available correlations obtained numerically and experimentally, both of which are in agreement with each other. The present model can also give the tortuosity with a good approximation near the percolation threshold. The validity of the present tortuosity model is thus verified.

Numerical experiments on one-dimensional model of turbulence

The initial-value problem of a forced Burgers equation is numerically solved by the Fourier expansion method. It is found that its solutions finally reach a steady state of 'laminar flow' which has no randomness and is stable to disturbances. Hence, strictly speaking, the so-called Burgers turbulence is not a turbulence. A new one-dimensional model is proposed to simulate the Navier-Stokes turbulence. A series of numerical experiments on this one-dimensional turbulence is made and is successful in obtaining Kolmogorov's (1941) k exp(-5/3) inertial-range spectrum. The (one-dimensional) Kolmogorov constant ranges from 0.5 to 0.65.

Improvement of different source function expressions in SWAN model for the Bohai Sea

The simulating wave nearshore (SWAN) wave model has been widely used in coastal areas, lakes and estuaries. However, we found a poor agreement between modeling results and measurements in analyzing the chosen four typical cases when we used the default parameters of the source function formulas of the SWAN to make wave simulation for the Bohai Sea. Also, it was found that at the same wind process the simulated results of two wind generation expressions (Komen, Janssen) demonstrated a large difference. Further study showed that the proportionality coefficient alpha in linear growth term of wave growth source function plays an unperceived role in the process of wave development. Based on experiments and analysis, we thought that the coefficient alpha should change rather than be a constant. Therefore, the coefficient alpha changing with the variation of friction velocity U (*) was introduced into the linear growth term of wave growth source function. Four weather processes were adopted to validate the improvement in the linear growth term. The results from the improved coefficient alpha agree much better with the measurements than those from the default constant coefficient alpha. Furthermore, the large differences of results between Komen wind generation expression and Janssen wind generation expression were eliminated. We also experimented with the four weather processes to test the new white-capping mechanisms based on the cumulative steepness method. It was found that the parameters of the new white-capping mechanisms are not suitable for the Bohai Sea, but Alkyon's white-capping mechanisms can be applicable to the Bohai Sea after amendments, demonstrating that this improvement of parameter alpha can improve the simulated results of the Bohai Sea.

Investigation of kerosene combustion characteristics with pilot hydrogen in model supersonic combustors

Experimental investigations on the ignition and combustion stabilization of kerosene with pilot hydrogen in Mach 2.5 airflows were conducted using two test combustors, with cross sections of 30.5 x 30 and 51 x 70 mm, respectively. Various integrated modules, including the combinations of different pilot injection schemes and recessed cavity flameholders with different geometries, were designed and tested. The stagnation pressure of vitiated air varied within the range of 1.1-1.8 NiPa, while the stagnation temperature varied from 1500 to 1900 K. Specifically, effects of the pilot hydrogen injection scheme, cavity geometry, and combustor scaling on the minimally required pilot hydrogen equivalence ratio were systematically examined. Results indicated that the cavity depth and length had significant effects on the ignition and flameholding, whereas the slanted angle of the aft wall was relatively less important. Two cavities in tandem were shown to be a more effective flameholding mechanism than that with a single cavity. The minimally required pilot hydrogen equivalence ratio for kerosene ignition and stable combustion was found to be as low as 0.02. Furthermore, combustion efficiency of 80% was demonstrated to be achievable for kerosene with the simultaneous use of pilot hydrogen and a recessed cavity to promote the ignition and global burning.

The development of an inhomogeneous damage field in an elastic-plastic model

By making use of the evolution equation of the damage field as derived from the statistical mesoscopic damage theory, we have preliminarily examined the inhomogeneous damage field in an elastic-plastic model under constant-velocity tension. Three types of deformation and damage field evolution are presented. The influence of the plastic matrix is examined. It seems that matrix plasticity may defer the failure due to damage evolution. A criterion for damage localization is consistent with the numerical results.

A new model of thermocapillary convection in floating half zone

A new model of thermocapillary convection in floating half zone was suggested in the present paper. The liquid bridge floats between two co-axis rods, the lower rod consists of metal with constant temperature and the upper rod consists of thermal insulating materials, where the normal gradient of temperature is nearly zero. In this case, the new model is relatively closer to simulate a half part of floating full zone in microgravity environment in comparison with the usual model of floating half zone. Basic features of the new model were studied by both numerical simulation and experiments, and the comparisons with the usual model were also discussed.

A probabilistic model for fatigue crack propagation analysis

A simple probabilistic model for predicting crack growth behavior under random loading is presented. In the model, the parameters c and m in the Paris-Erdogan Equation are taken as random variables, and their stochastic characteristic values are obtained through fatigue crack propagation tests on an offshore structural steel under constant amplitude loading. Furthermore, by using the Monte Carlo simulation technique, the fatigue crack propagation life to reach a given crack length is predicted. The tests are conducted to verify the applicability of the theoretical prediction of the fatigue crack propagation.

Equivalent Model Of Expansion Of Cement Mortar Under Sulphate Erosion

The expansion property of cement mortar under the attack of sulfate ions is studied by experimental and theoretical methods. First, cement mortars are fabricated with the ratio of water to cement of 0.4, 0.6, and 0.8. Secondly, the expansion of specimen immerged in sulphate solution is measured at different times. Thirdly, a theoretical model of expansion of cement mortar under sulphate erosion is suggested by virtue of represent volume element method. In this model, the damage evolution due to the interaction between delayed ettringite and cement mortar is taken into account. Finally, the numerical calculation is performed. The numerical and experimental results indicate that the model perfectly describes the expansion of the cement mortar.

A Multimoment Finite-Volume Shallow-Water Model On The Yin-Yang Overset Spherical Grid

A numerical model for shallow-water equations has been built and tested on the Yin-Yang overset spherical grid. A high-order multimoment finite-volume method is used for the spatial discretization in which two kinds of so-called moments of the physical field [i.e., the volume integrated average ( VIA) and the point value (PV)] are treated as the model variables and updated separately in time. In the present model, the PV is computed by the semi-implicit semi-Lagrangian formulation, whereas the VIA is predicted in time via a flux-based finite-volume method and is numerically conserved on each component grid. The concept of including an extra moment (i.e., the volume-integrated value) to enforce the numerical conservativeness provides a general methodology and applies to the existing semi-implicit semi-Lagrangian formulations. Based on both VIA and PV, the high-order interpolation reconstruction can only be done over a single grid cell, which then minimizes the overlapping zone between the Yin and Yang components and effectively reduces the numerical errors introduced in the interpolation required to communicate the data between the two components. The present model completely gets around the singularity and grid convergence in the polar regions of the conventional longitude-latitude grid. Being an issue demanding further investigation, the high-order interpolation across the overlapping region of the Yin-Yang grid in the current model does not rigorously guarantee the numerical conservativeness. Nevertheless, these numerical tests show that the global conservation error in the present model is negligibly small. The model has competitive accuracy and efficiency.

Experimental investigation on the cavity-based scramjet model

The present work focused on improving the engine performance with different fuel equivalence ratios and fuel injections. A scramjet model with strut/cavity integrated configurations was tested under Mach 5.8 flows. The results showed that the strut may sreve as an effective tool in a kerosene-fueled scramjet. The integration of strut/cavities also had great effect on stablizing the combustion in a wide range of fuel equivalence ratio. The one-sdimensional analysis method was used to analyze the main characteristics of the model. The two-stage fuel injection should have better performance in increasing the chemical reaction rate in the first cavity region.

Nonlinear fatigue damage model based on the residual strength degradation law

In this paper, a logarithmic expression to describe the residual strength degradation process is developed in order to fatigue test results for normalized carbon steel. The definition and expression of fatigue damage due to symmetrical stress with a constant amplitude are also given. The expression of fatigue damage can also explain the nonlinear properties of fatigue damage. Furthermore, the fatigue damage of structures under random stress is analyzed, and an iterative formula to describe the fatigue damage process is deduced. Finally, an approximate method for evaluating the fatigue life of structures under repeated random stress blocking is presented through various calculation examples.

A generalized model for the effective thermal conductivity of porous media based on self-similarity

A generalized model for the effective thermal conductivity of porous media is derived based on the fact that statistical self-similarity exists in porous media. The proposed model assumes that porous media consist of two portions: randomly distributed non-touching particles and self-similarly distributed particles contacting each other with resistance. The latter are simulated by Sierpinski carpets with side length L = 13 and cutout size C = 3, 5, 7 and 9, respectively, depending upon the porosity concerned. Recursive formulae are presented and expressed as a function of porosity, ratio of areas, ratio of component thermal conductivities and contact resistance, and there is no empirical constant and every parameter has a clear physical meaning. The model predictions are compared with the existing experimental data, and good agreement is found in a wide range of porosity of 0.14-0.80, and this verifies the validity of the proposed model.

A non-equilibrium sediment transport model for rill erosion

Sediment transport in rill flows exhibits the characteristics of non-equilibrium transport, and the sediment transport rate of rill flow gradually recovers along the flow direction by erosion. By employing the concept of partial equilibrium sediment transport from open channel hydraulics, a dynamic model of rill erosion on hillslopes was developed. In the model, a parameter, called the restoration coefficient of sediment transport capacity, was used to express the recovery process of sediment transport rate, which was analysed by dimensional analysis and determined from laboratory experimental data. The values of soil loss simulated by the model were in agreement with observed values. The model results showed that the length and gradient of the hillslope and rainfall intensity had different influences on rill erosion. Copyright (c) 2006 John Wiley & Sons, Ltd.

Numerical investigation on a simulation model of floating zone convection

A simulation model with adiabatic condition at the upper rod and constant temperature at the lower rod is studied numerically in this paper. The temperature distribution in a simulation model is closer to the one in the half part of a floating full zone in comparison with the one in a usual floating half zone model with constant temperature at both rods, because the temperature distribution of a floating full zone is symmetric for the middle plane in a microgravity environment. The results of the simulation model show that the temperature profiles and the how patterns are different from those of the usual floating half zone model. Another type of half zone model, with a special non-uniform temperature distribution at the upper rod and constant temperature at the lower rod, has been suggested by recent experiments. The temperature boundary condition of the upper rod has a maximum value in the center and a lower value near the free surface. This modified simulation model is also simulated numerically in the present paper. Copyright (C)1996 Elsevier Science Ltd.

Lattice instability at a fast moving crack tip

A molecular dynamics method is used to analyze the dynamic propagation of an atomistic crack tip. The simulation shows that the crack propagates at a relatively constant global velocity which is well below the Rayleigh wave velocity. However the local propagation velocity oscillates violently, and it is limited by the longitudinal wave velocity. The crack velocity oscillation is caused by a repeated process of crack tip blunting and sharpening. When the crack tip opening displacement exceeds a certain critical value, a lattice instability takes place and results in dislocation emissions from the crack tip. Based on this concept, a criterion for dislocation emission from a moving crack tip is proposed. The simulation also identifies the emitted dislocation as a source for microcrack nucleation. A simple method is used to examine this nucleation process. (C) 1996 American Institute of Physics.