Simulation of Effective Thermal Conductivity of Metal Hydride Packed Beds

Several mathematical models are available for estimation of effective thermal conductivity of nonreactive packed beds. Keeping in view the salient differences between metal hydride beds in which chemisorption of hydrogen takes place and conventional nonreactive packed beds, modified models are proposed here to predict the effective thermal conductivity. Variation in properties such as solid thermal conductivity and porosity during hydrogen absorption and desorption processes are incorporated. These extended models have been applied to simulate the effective thermal conductivity of the MmNi(4.5)Al(0.5) hydride bed and are compared with the experimental results. Applicability of the extended models for estimation of the effective thermal conductivity at different operating conditions such as pressure, temperature, and hydrogen concentration is discussed.

Investigations of aerosol black carbon from a semi-urban site in the Indo-Gangetic Plain region

Long-term (2009-2012) data from ground-based measurements of aerosol black carbon (BC) from a semi-urban site, Pantnagar (29.0 degrees N, 79.5 degrees E, 231 m amsl), in the Indo-Gangetic Plain (IGP) near the Himalayan foothills are analyzed to study the regional characterization. Large variations are seen in BC at both diurnal and seasonal scales, associated with the mesoscale and synoptic meteorological processes, and local/regional anthropogenic activities. BC diurnal variations show two peaks (morning and evening) arising from the combined effects of the atmospheric boundary layer (ABL) dynamics and local emissions. The diurnal amplitudes as well as the rates of diurnal evolution are the highest in winter season, followed by autumn, and the lowest in summer-monsoon. BC exhibits nearly an inverse relation with mixing layer depth in all seasons; being strongest in winter (R-2 = 0.89) and weakest (R-2 = 0.33) in monsoon (July-August). Unlike BC, co-located aerosol optical depths (AOD) and aerosol absorption are highest in spring over IGP, probably due to the presence of higher abundances of aerosols (including dust) above the ABL (in the free troposphere). AOD (500 nm) showed annual peak (>0.6) in May-June, dominated by coarse mode, while fine mode aerosols dominated in late autumn and early winter. Aerosols profiles from CALIPSO show highest values close to the surface in winter/autumn, similar to the feature seen in surface BC, whereas at altitudes > 2 km, the extinction is maximum in spring/summer. WRF-Chem model is used to simulate BC temporal variations and then compared with observed BC. The model captures most of the important features of the diurnal and seasonal variations but significantly underestimated the observed BC levels, suggesting improvements in diurnal and seasonal varying BC emissions apart from the boundary layer processes. (C) 2015 Elsevier Ltd. All rights reserved.

Energy conversion in shape memory alloy heat engine - Part II: Simulation

In this paper, a theoretical model proposed in Part I (Zhu et al., 2001a) is used to simulate the behavior of a twin crank NiTi SMA spring based heat engine, which has been experimentally studied by Iwanaga et al. (1988). The simulation results are compared favorably with the measurements. It is found that (1) output torque and heat efficiency decrease as rotation speed increase; (2) both output torque and output power increase with the increase of hot water temperature; (3) at high rotation speed, higher water temperature improves the heat efficiency. On the contrary, at low rotation speed, lower water temperature is more efficient; (4) the effects of initial spring length may not be monotonic as reported. According to the simulation, output torque, output power and heat efficiency increase with the decrease of spring length only in the low rotation speed case. At high rotation speed, the result might be on the contrary.

Numerical Investigation of Richtmyer-Meshkov Instability Driven by Cylindrical Shocks

In this paper, a numerical method with high order accuracy and high resolution was developed to simulate the Richtmyer-Meshkov(RM) instability driven by cylindrical shock waves. Compressible Euler equations in cylindrical coordinate were adopted for the cylindrical geometry and a third order accurate group control scheme was adopted to discretize the equations. Moreover, an adaptive grid technique was developed to refine the grid near the moving interface to improve the resolution of numerical solutions. The results of simulation exhibited the evolution process of RM instability, and the effect of Atwood number was studied. The larger the absolute value of Atwood number, the larger the perturbation amplitude. The nonlinear effect manifests more evidently in cylindrical geometry. The shock reflected from the pole center accelerates the interface for the second time, considerably complicating the interface evolution process, and such phenomena of reshock and secondary shock were studied.

Molecular dynamics simulation of the role of dislocations in microcrack healing

The molecular dynamics method is used to simulate microcrack healing during heating or/and under compressive stress. A centre microcrack in Cu crystal would be sealed under compressive stress or by heating. The role of compressive stress and heating in crack healing was additive. During microcrack healing, dislocation generation and motion occurred. When there were pre-existing dislocations around the microcrack, the critical temperature or compressive stress necessary for microcrack healing would decrease, and, the higher the number of dislocations, the lower the critical temperature or compressive stress. The critical temperature necessary for microcrack healing depended upon the orientation of the crack plane. For example, the critical temperature for the crack along the (001) plane was the lowest, i.e. 770K.

Monte-Carlo simulation of surface crack growth rate for offshore structural steel E36-Z35

The Monte- Carlo method is used to simulate the surface fatigue crack growth rate for offshore structural steel E36-Z35, and to determine the distributions and relevance of the parameters in the Paris equation. By this method, the time and cost of fatigue crack propagation testing can be reduced. The application of the method is demonstrated by use of four sets of fatigue crack propagation data for offshore structural steel E36-Z35. A comparison of the test data with the theoretical prediction for surface crack growth rate shows the application of the simulation method to the fatigue crack propagation tests is successful.

Intensity Distribution Design for Laser-Induced Thermal Loading Based on Numerical Simulation

To accomplish laser-induced thermal loading simulation tests for pistons,the Gaussian beam was modulated into multi-circular beam with specific intensity distribution.A reverse method was proposed to design the intensity distribution for the laser-induced thermal loading based on finite element(FE) analysis.Firstly,the FE model is improved by alternating parameters of boundary conditions and thermal-physical properties of piston material in a reasonable range,therefore it can simulate the experimental resul...

A Nanomechanical Device Based on Light-Driven Proton Pumps

In this paper, a hybrid device based on a microcantilever interfaced with bacteriorhodopsin (bR) is constructed. The microcantilever, on which the highly oriented bR film is self-assembled, undergoes controllable and reversible bending when the light-driven proton pump protein, bR, on the microcantilever surface is activated by visible light. Several control experiments are carried out to preclude the influence of heat and photothermal effects. It is shown that the nanomechanical motion is induced by the resulting gradient of protons, which are transported from the KCl solution on the cytoplasmic side of the bR film towards the extracellular side of the bR film. Along with a simple physical interpretation, the microfabricated cantilever interfaced with the organized molecular film of bR can simulate the natural machinery in converting solar energy to mechanical energy.

An Analytical Solution for Three-Dimensional Diffraction of Plane P-Waves by a Hemispherical Alluvial Valley with Saturated Soil Deposits

An analytical solution for the three-dimensional scattering and diffraction of plane P-waves by a hemispherical alluvial valley with saturated soil deposits is developed by employing Fourier-Bessel series expansion technique. Unlike previous studies, in which the saturated soil deposits were simulated with the single-phase elastic theory, in this paper, they are simulated with Biot's dynamic theory for saturated porous media, and the half space is assumed as a single-phase elastic medium. The effects of the dimensionless frequency, the incidence angle of P-wave and the porosity of soil deposits on the surface displacement magnifications of the hemispherical alluvial valley are investigated. Numerical results show that the existence of a saturated hemispherical alluvial valley has much influence on the surface displacement magnifications. It is more reasonable to simulate soil deposits with Biot's dynamic theory when evaluating the displacement responses of a hemispherical alluvial valley with an incidence of P-waves.

Rossby Vortex Simulation on a Paraboloidal Coordinate System Using the Lattice Boltzmann Method

In this paper, we apply our compressible lattice Boltzmann model to a rotating parabolic coordinate system to simulate Rossby vortices emerging in a layer of shallow water flowing zonally in a rotating paraboloidal vessel. By introducing a scaling factor, nonuniform curvilinear mesh can be mapped to a flat uniform mesh and then normal lattice Boltzmann method works. Since the mass per unit area on the two-dimensional (2D) surface varies with the thickness of the water layer, the 2D flow seems to be "compressible" and our compressible model is applied. Simulation solutions meet with the experimental observations qualitatively. Based on this research, quantitative solutions and many natural phenomena simulations in planetary atmospheres, oceans, and magnetized plasma, such as the famous Jovian Giant Red Spot, the Galactic Spiral-vortex, the Gulf Stream, and the Kuroshio Current, etc,, can be expected.

Acoustic Emission Experiments of Rock Failure under Load Simulating the Hypocenter Condition

A series of acoustic emission (AE) experiments of rock failure have been conducted under cyclic load in tri-axial stress tests. To simulate the hypocenter condition the specimens are loaded by the combined action of a constant stress, intended to simulate

Effect of Rock Mass Structure and Block Size on the Slope Stability-Physical Modeling and Discrete Element Simulation

This paper studies the stability of jointed rock slopes by using our improved three-dimensional discrete element methods (DEM) and physical modeling. Results show that the DEM can simulate all failure modes of rock slopes with different joint configurations. The stress in each rock block is not homogeneous and blocks rotate in failure development. Failure modes depend on the configuration of joints. Toppling failure is observed for the slope with straight joints and sliding failure is observed for the slope with staged joints. The DEM results are also compared with those of limit equilibrium method (LEM). Without considering the joints in rock masses, the LEM predicts much higher factor of safety than physical modeling and DEM. The failure mode and factor of safety predicted by the DEM are in good agreement with laboratory tests for any jointed rock slope.

Numerical Simulation of Nox Formation in Coal Combustion with Inlet Natural Gas Burning

A full two-fluid model of reacting gas-particle flows and coal combustion is used to simulate coal combustion with and without inlet natural gas added in the inlet. The simulation results for the case without natural gas burning is in fair agreement with the experimental results reported in references. The simulation results of different natural gas adding positions indicate that the natural gas burning can form lean oxygen combustion enviroment at the combustor inlet region and the NOz concentration is reduced. The same result can be obtained from chemical equilibrium analysis.

Lattice Boltzmann method for compressible flows with high Mach numbers

In this paper we present a lattice Boltzmann model to simulate compressible flows by introducing an attractive force. This scheme has two main advantages: one is to soften sound speed effectively, which greatly raises the Mach number (up to 5); another is its relative simple procedure. Simulations of the March cone and the comparison between theoretical expectations and simulations demonstrate that the scheme is effective in the simulation of compressible flows with high Mach numbers, which would create many new applications.

Three-Dimensional Network Model And Its Parameter Calibration

A material model, whose framework is parallel spring-bundles oriented in 3-D space, is proposed. Based on a discussion of the discrete schemes and optimum discretization of the solid angles, a 3-D network cell consisted of one-dimensional components is developed with its geometrical and physical parameters calibrated. It is proved that the 3-D network model is able to exactly simulate materials with arbitrary Poisson ratio from 0 to 1/2, breaking through the limit that the previous models in the literature are only suitable for materials with Poisson ratio from 0 to 1/3. A simplified model is also proposed to realize high computation accuracy within low computation cost. Examples demonstrate that the 3-D network model has particular superiority in the simulation of short-fiber reinforced composites.