Large Eddy Simulation of Pseudo-Shock Waves Using Wall Model

A supersonic turbulent flowfield involving the pseudo-shock waves in an isolator of a supersonic combustion ramjet is computed using two different LES codes which are a high-order upwind finite volume scheme, and a sixth order compact differencing scheme utilizing the localized artificial diffusivity method for stabilizing shock waves and employing a wall model to enable the use of coarse mesh. In the validation study where a supersonic turbulent boundary layer flow over a flat plate is examined, both LES codes are well validated using velocity profile in the boundary layer given by the hot-wire anemometry and normal stress given by the laser Doppler anemometry. In particular, the sixth order compact differencing scheme gives closer agreements with these experimental data. Then, the validated LES codes are applied to solve the Mach number 2.5 supersonic turbulent flowfield involving the pseudo-shock waves. It is shown that typical features of unsteady flowfield of the pseudo-shock waves are well obtained by both schemes. Again, it is indicated that the sixth order compact differencing scheme gives closer agreements with the existing velocity data obtained by particle image velocimetry and pressure fluctuation data on the wall surface. Besides, the computational cost of the compact differencing scheme is found to be 1/7 of that for the upwind finite volume scheme, even though a wall model is solved at each grid point on the wall surface. Therefore, the obtained results in the present study allow recommending the sixth order compact differencing scheme with a wall model for simulating supersonic turbulent flowfield in an isolator involving the pseudo-shock waves.

SPH Particle Collisions for the Reduction of Particle Clustering, Interface Stabilisation and Wall Modelling

The pair-wise forces in the SPH momentum equation guarantee the conservation of momentum, but they cannot prevent particle clustering and wall penetration. Particle clustering may occur for several reasons. A fundamental issue is the tensile instability, which is caused by negative numerical pressures. Clustering may also occur due to certain properties of the kernel gradient. Discontinuities in the pressure and its gradient, due to surface tension and gravity, may cause particle instabilities near the interface between two fluids. Wall penetration is also a form of particle clustering. In this paper the particle collision concept is introduced to suppress particle clustering. Here, the use of kinematic conditions (motion) rather than dynamic conditions (forces) is explored. These kinematic conditions are obtained from kinetic collision theory. Conservation of momentum is maintained, and under elastic conditions conservation of energy as well. The particle collision model only becomes active when needed. It may be seen as a particle shifting method, in the sense that the velocities are changed, and as a consequence of that the particle positions change. It is demonstrated in several case studies that the particle collision model allows for realistic (low) viscosities. It was also found to stabilise the interface between two fluids up to high, realistic density ratios (1000:1) in typical liquid-gas applications. As such it can be used as a multi-fluid model. The concept allows for real wave speed ratios (and far beyond), which, as well as real viscosities, are essential in the modelling of heat transfer applications. The collisions with walls allow for no-slip conditions at real viscosities while wall penetration is suppressed. In summary, the particle collision model makes SPH more robust for engineering.

Scalar glueball in a soft-wall model of AdS/QCD

The scalar glueball is investigated in a soft-wall model of AdS/QCD. Constraints of the mass of the scalar glueball are given through an analysis of a relation between the bulk mass and the anomalous dimension. The mass of the ground scalar glueball is located at 0.96＋-0.07 ＋0.04 GeV 〈 MG 〈 1.36＋0.05 -0.10 GeV. In terms of a background dilaton field θ（z）= cz2, the two-point correlation function for the scalar gluon operator is obtained. The two-point correlation function at A=4 gives a different behavior compared with the one in QCD.

Wall-modeling for large-eddy simulation of flows around an axisymmetric body using the diffuse-interface immersed boundary method

A novel method is proposed to combine the wall-modeled large-eddy simulation(LES) with the diffuse-interface direct-forcing immersed boundary(IB) method.The new developments in this method include:(i) the momentum equation is integrated along the wall-normal direction to link the tangential component of the effective body force for the IB method to the wall shear stress predicted by the wall model;(ii) a set of Lagrangian points near the wall are introduced to compute the normal component of the effective body force for the IB method by reconstructing the normal component of the velocity. This novel method will be a classical direct-forcing IB method if the grid is fine enough to resolve the flow near the wall. The method is used to simulate the flows around the DARPA SUBOFF model. The results obtained are well comparable to the measured experimental data and wall-resolved LES results.

A NEW UNSTEADY THREE DIMENSIONAL MODEL FOR MACROMOLECULAR TRANSPORT AND WATER FILTRATION ACROSS THE ARTERIAL WALL

A new unsteady three-dimensional convective-diffusive mathematical model for the transportation of macromolecules and water across the arterial wall was proposed . After the formation of leaky junctions due to the mitosis of endothelial cell of the arterial wall, the macromolecular transport happens surrounding the leaky cells. The arterial wall was divided into four layers: the endothelial layer, the subendothelial intima, the internal elastic lamina and the media for the convenience of research. The time-dependent concentration growth, the effect of the shape of endothelial cell and the effect of physiological parameters were analyzed. The analytical solution of velocity field and pressure field of water flow across the arterial wall were obtained; and concentration distribution of three macromolecules; LDL, HRP and Albumin, were calculated with numerical simulation method. The new theory predicts, the maximum and distribution areas of time dependent concentration with round-shape endothelial cell are both larger than that with ellipse-shape endothelial cell. The model also predicts the concentration growth is much alike that of a two-dimensional model and it shows that the concentration reaches its peak at the leaky junction where atherosclerotic formation frequently occurs and falls down rapidly in a limited area beginning from its earlier-time growth to the state when macromolecular transfer approaches steadily. These predictions of the new model are in agreement with the experimental observation for the growth and concentration distribution of LDL and Albumin.

Calibration of a hysteretic model for glass fiber reinforced gypsum wall panels

Glass fi ber reinforced gypsum(GFRG) wall panels are prefabricated panels with hollow cores, originally developed in Australia and subsequently adopted by India and China for use in buildings. This paper discusses identifi cation and calibration of a suitable hysteretic model for GFRG wall panels fi lled with reinforced concrete. As considerable pinching was observed in the experimental results, a suitable hysteretic model with pinched hysteretic rule is used to conduct a series of quasi-static as inelastic hysteretic response analyses of GFRG panels with two different widths. The calibration of the pinching model parameters was carried out to approximately match the simulated and experimental responses up to 80% of the peak load in the post peak region. Interestingly, the same values of various parameters(energy dissipation and pinching related parameters) were obtained for all fi ve test specimens.

New perspective in statistical modeling of wall-bounded turbulence

Despite dedicated effort for many decades,statistical description of highly technologically important wall turbulence remains a great challenge.Current models are unfortunately incomplete,or empirical,or qualitative.After a review of the existing theories of wall turbulence,we present a new framework,called the structure ensemble dynamics （SED）,which aims at integrating the turbulence dynamics into a quantitative description of the mean flow.The SED theory naturally evolves from a statistical physics understanding of non-equilibrium open systems,such as fluid turbulence, for which mean quantities are intimately coupled with the fluctuation dynamics.Starting from the ensemble-averaged Navier-Stokes（EANS） equations,the theory postulates the existence of a finite number of statistical states yielding a multi-layer picture for wall turbulence.Then,it uses order functions（ratios of terms in the mean momentum as well as energy equations） to characterize the states and transitions between states.Application of the SED analysis to an incompressible channel flow and a compressible turbulent boundary layer shows that the order functions successfully reveal the multi-layer structure for wall-bounded turbulence, which arises as a quantitative extension of the traditional view in terms of sub-layer,buffer layer,log layer and wake. Furthermore,an idea of using a set of hyperbolic functions for modeling transitions between layers is proposed for a quantitative model of order functions across the entire flow domain.We conclude that the SED provides a theoretical framework for expressing the yet-unknown effects of fluctuation structures on the mean quantities,and offers new methods to analyze experimental and simulation data.Combined with asymptotic analysis,it also offers a way to evaluate convergence of simulations.The SED approach successfully describes the dynamics at both momentum and energy levels, in contrast with all prevalent approaches describing the mean velocity profile only.Moreover,the SED theoretical framework is general,independent of the flow system to study, while the actual functional form of the order functions may vary from flow to flow.We assert that as the knowledge of order functions is accumulated and as more flows are analyzed, new principles（such as hierarchy,symmetry,group invariance,etc.） governing the role of turbulent structures in the mean flow properties will be clarified and a viable theory of turbulence might emerge.

Effects of Magnetic Shear on the Resistive Wall Mode Instability in a Cylindrical Model

The linear evolution of a resistive wall mode(RWM)with magnetic shears was analytically studied in a cylindrical geometry.The incompressible magnetohydrodynamic(MHD)equations were solved by the Fourier analysis method,and various equilibrium magnetic fields were considered.The shear in the magnetic field had an unstable effect on the linear evolution of the RWM.The linear growth rate increased obviously with increase of the magnetic shear rate for higher magnetic shears.Slow plasma rotation could stabilize the RWM with low magnetic shears,but the plasma rotation could not stabilize the RWM with high magnetic shears.The stabilizing effects of the wall conductivity on the RWM are more efficient for lager magnetic shear.

Residual stress modeling of thin wall by laser cladding forming

The residual stress generated in the laser cladding could lead to undesirable distortions or even crack formation. In order to better understand the evolution/yielding process of stress field,a 3 D finite-element thermo-mechanical model was established for the laser cladding formation of thin wall with the 17-4 PH powder on the FV520( B) steel. The temperature field was firstly analyzed,based on which the stress field and strain field of the laser cladding forming process were analyzed.In order to validate the prediction,the final residual stress field in the obtained thin wall was tested by X-ray diffraction in comparison with the predicted results.

Cavitation Bubble Collapse near a Curved Wall by the Multiple-Relaxation-Time Shan-Chen Lattice Boltzmann Model

The cavitation bubble collapse near a cell can cause damage to the cell wall. This effect has received increasing attention in biomedical supersonics. Based on the lattice Boltzmann method, a multiple-relaxation-time Shan–Chen model is built to study the cavitation bubble collapse. Using this model, the cavitation phenomena induced by density perturbation are simulated to obtain the coexistence densities at certain temperature and to demonstrate the Young–Laplace equation. Then, the cavitation bubble collapse near a curved rigid wall and the consequent high-speed jet towards the wall are simulated. Moreover, the influences of initial pressure difference and bubble-wall distance on the cavitation bubble collapse are investigated.

SIMULATED MATHEMATICAL MODEL FOR HIGH PRESSURE CYLINDER WALL TEMPERATURE OF HOMEMADE TURBINE

In order to balance the contradiction between the demand of high precision and that of short time interval of model computing for the power plant simulator, a set of simulated mathematical models are constructed. The model describes the cylinder wall temperature located at four key positions of the high pressure cylinder. The simulated model is confirmed to be not only simple but also precise via comparison between the simulated results and the autoptic data of a power plant.

Different Ablation Models for the Wall-plasma Interaction Process in Pulsed Plasma Thruster

Ablation excited by current pulses is a very critical physical process in pulse plasma thrusters(PPT).Its effects on wall-plasma interaction directly determine the PPT performances.In order to reveal the process of the ablated wall interaction with the discharge plasma in PPT,ablation models formulated by three different boundary conditions at the wall-plasma interface are studied.These are the two widely used high-speed evaporation models(Model-L and Model-M),and the recently developed Keida-Zaghloul model(Model-K)of the Knudsen layer that takes into account the internal degrees of freedom on the energy flux conservation.First,fundamental mechanisms of the three ablation models are clarified by comparative analysis in order to gain a comprehensive understanding of the wall-plasma interaction.Then,the applicability of different ablation models with the numerical solutions of LES-6 PPT is investigated in detail using magnetohydrodynamic(MHD)modeling.Results show that Model-L and Model-M are actually special cases of Model-K when a simplified jump conditions limited by high velocity at the vapor/plasma interface is used;A ratio of ablation rate in Model-L to that in Model-M is about 0.8at the same wall surface temperature,while it rises to 1 at different surface temperature determined by Model-L and Model-M in PPT.Even though Model-K solution requires significant computational time,it shows more accurate ablation feature for the wall-plasma interaction and possesses better computing precision of impulse bit during post-pulse which is useful for future studies of the late time ablation.

Finite element modeling of power spinning of thin-walled shell with hoop inner rib

A 3D elasto-plastic finite element(FE)model of power spinning of thin-walled aluminum alloy shell with hoop inner rib was established under software ABAQUS.Key technologies were dealt with reasonably.The reliability of the FE model was verified theoretically and experimentally.The forming process was simulated and studied.The distribution of the thickness and stress,and the variations of spinning force were obtained.The workpiece springback was analyzed with ABAQUS/Standard.The results show that the FE model considering elastic deformation can not only be used to analyze the workpiece springback in the complex spinning process,but also serve as a significant guide to study the local deformation mechanism and choose the reasonable parameters.

Modeling for Collapsing Cavitation Bubble near Rough Solid Wall by Mulit-Relaxation-Time Pseudopotential Lattice Boltzmann Model

Cavitation bubble collapse near rough solid wall is modeled by the multi-relaxation-time (MRT) pseudopotential lattice Boltzmann (LB) model. The modified forcing scheme, which can achieve LB model’s thermodynamic consistency by tuning a parameter related with the particle interaction range, is adopted to achieve desired stability and density ratio. The bubble collapse near rough solid wall was simulated by the improved MRT pseudopotential LB model. The mechanism of bubble collapse is studied by investigating the bubble profiles, pressure field and velocity field evolution. The eroding effects of collapsing bubble are analyzed in details. It is found that the process and the effect of the interaction between bubble collapse and rough solid wall are affected seriously by the geometry of solid boundary. At the same time, it demonstrates that the MRT pseudopotential LB model is a potential tool for the investigation of the interaction mechanism between the collapsing bubble and complex geometry boundary.

A new attenuation model of near-fault ground motions with consideration of the hanging wall effect in the Wenchuan earthquake

The hanging wall effect is an important factor that impacts the characteristics of strong ground motions in near-fault areas. Based on a residual analysis of ground motion parameters characterizing the hanging wall effect and in recognition of the nature of the effect, many models have been developed. In this study, after a comprehensive analysis of two existing models, a new model is proposed and used to model the hanging wall effect in horizontal peak ground acceleration （PGAH） and spectral acceleration （SAH） at a period of 0.1s in the Wenchuan earthquake. Finally, comparisons between the modeling results of the hanging wall effect in the Wenchuan earthquake and the results predicted by using Abrahamson and Silva＇s NGA model （AS NGA） indicate that the AS NGA model predicts a much higher hanging wall effect than the model developed in this paper. Furthermore, the AS NGA model predicts a large hanging wall effect even at great distances, while the proposed model more accurately captures the trend of the effect.

Drinfeld twist and the domain wall partition function of the eight-vertex model

With the help of the F-basis provided by the Drinfeld twist or factorising F-matrix for the spatial optical soliton model associated with the eight-vertex model, we calculate the partition function for the eight-vertex model on an N ×N square lattice with domain wall boundary condition.

A NEW UNSTEADY THREE DIMENSIONAL MODEL FOR MACROMOLECULAR TRANSPORT AND WATER FILTRATION ACROSS THE ARTERIAL WALL

A new unsteady three-dimensional convective-diffusive mathematical model for the transportation of macromolecules and water across the arterial wall was proposed . After the formation of leaky junctions due to the mitosis of endothelial cell of the arterial wall, the macromolecular transport happens surrounding the leaky cells. The arterial wall was divided into four layers: the endothelial layer, the subendothelial intima, the internal elastic lamina and the media for the convenience of research. The time-dependent concentration growth,the effect of the shape of endothelial cell and the effect of physiological parameters were analyzed. The analytical solution of velocity field and pressure field of water flow across the arterial wall were obtained; and concentration distribution of three macromolecules ; LDL,HRP and Albumin, were calculated with numerical simulation method. The new theory predicts, the maximum and distribution areas of time dependent concentration with round shape endothelial cell are both larger than that with ellipse-shape endothelial cell. The model also predicts the concentration growth is much alike that of a two-dimensional model and it shows that the concentration reaches its peak at the leaky junction where atherosclerotic formation frequently occurs and falls down rapidly in a limited area beginning from its earlier time growth to the state when macromolecular transfer approaches steadily. These predictions of the new model are in agreement with the experimental observation for the growth and concentration distribution of LDL and Albumin.

Theoretical Model and Random Simulation for Nonlinear Interaction of Irregular Waves with Vertical Wall

According to the theoretical solutions for the nonlinear three-dimensional gravity surface waves and their interactions with vertical wall previously proposed by the lead author, in this paper an exact second-order random model of the unified wave motion process for nonlinear irregular waves and their interactions with vertical wall in uniform current is formulated, the corresponding theoretical nonlinear spectrum is derived, and the digital simulation model suitable to the use of the FFT (Fast Fourier Transform) algorithm is also given. Simulations of wave surface, wave pressure, total wave pressure and its moment are performed. The probability properties and statistical characteristics of these realizations are tested, which include the verifications of normality for linear process and of non-normality for nonlinear process; the consistencies of the theoretical spectra with simulated ones; the probability properties of apparent characteristics, such as amplitudes, periods, and extremes (maximum and minimum, positive and negative extremes). The statistical analysis and comparisons demonstrate that the proposed theoretical and computing models are realistic and effective, and estimated spectra are in good agreement with the theoretical ones, and the probability properties of the simulated waves are similar to those of the sea waves. At the same time, the simulating computation can be completed rapidly and easily.

Fault detection and identification based on combining logic and model in a wall-climbing robot

A combined logic- and model-based approach to fault detection and identification （FDI） in a suction foot control system of a wall-climbing robot is presented in this paper. For the control system, some fault models are derived by kinematics analysis. Moreover, the logic relations of the system states are known in advance. First, a fault tree is used to analyze the system by evaluating the basic events （elementary causes）, which can lead to a root event （a particular fault）. Then, a multiple-model adaptive estimation algorithm is used to detect and identify the model-known faults. Finally, based on the system states of the robot and the results of the estimation, the model-unknown faults are also identified using logical reasoning. Experiments show that the proposed approach based on the combination of logical reasoning and model estimating is efficient in the FDI of the robot.

Field measurements for calibration of simplified models of the stiffening effect of infill masonry walls in high-rise RC framed and shear-wall buildings

As a type of nonstructural component, infill walls play a significant role in the seismic behavior of high-rise buildings. However, the stiffness of the infill wall is generally either ignored or considered by simplified empirical criteria that lead to a period shortening. The difference can be greatly decreased by using a structural identification methodology. In this study, an ambient vibration test was performed on four on-site reinforced concrete high-rise buildings, and the design results were compared with the PKPM models using corresponding finite element(FE) models. A diagonal strut model was used to simulate the behavior of the infill wall, and the identified modal parameters measured from the on-site test were employed to calibrate the parameters of the diagonal strut in the FE models. The SAP2000 models with calibrated elastic modulus were used to evaluate the seismic response in the elastic state. Based on the load-displacement relationship of the infill wall, nonlinear dynamic analysis models were built in PERFORM-3 D and calibrated using the measured modal periods. The analysis results revealed that the structural performance under small/large earthquake records were both strengthened by infill walls, and the contribution of infill walls should be considered for better accuracy in the design process.