Numerical Modelling of Coupled Heat and Mass Transfer in Porous Materials: Application to Cinder Block Bricks

In this work, we present numerical modelling of coupled heat and mass transfer within porous materials. Our study focuses on cinder block bricks generally used in building construction. The material is assumed to be placed in air. Moisture content and temperature have been chosen as the main transfer drivers and the equations governing these transfer drivers are based on the Luikov model. These equations are solved by an implicit finite difference scheme. A Fortran code associated with the Thomas algorithm was used to solve the equations. The results show that heat and mass transfer depend on the temperature of the air in contact with the material. As this air temperature rises, the temperature within the material increases, and more rapidly at the material surface. Also, thermal conductivity plays a very important role in the thermal conduction of building materials and influences heat and mass transfer in these materials. Materials with higher thermal conductivity diffuse more heat.

Numerical modelling of structural controls on fluid flow and mineralization

This paper presents the results of a set of numerical models focussing on structural controls on hydrothermal mineralization. We first give an overview of natural phenomena of structurally-controlled ore formation and the background theory and mechanisms for such controls. We then provide the results of a group of simple 2D numerical models validated through comparison with Cu-vein structure observed near the Shilu Copper deposit （Yangchun, Guangdong Province, China） and finally a case study of 3D numerical modelling applied to the Hodgkinson Province in North Queensland （Australia）. Two modelling approaches, discrete deformation modelling and continuum coupled deformation and fluid flow modelling, are involved. The 2D model-derived patterns are remarkably consistent with the Cu-vein structure from the Shilu Copper deposit, and show that both modelling approaches can realistically simulate the mechanical behaviours of shear and dilatant fractures. The continuum coupled deformation and fluid flow model indicates that pattern of the Cu- veins near the Shilu deposit is the result of shear strain localization, development of dilation and fluid focussing into the dilatant fracture segments. The 3D case-study models （with deformation and fluid flow coupling） on the Hodgkinson Province generated a number of potential gold mineralization

3-D numerical modelling of Domino failure of hard rock pillars in Fetr6 Chromite Mine, Iran, and comparison with empirical methods

Fetr6 is an underground mine using the stope-and-pillar mining method. As there was some evidence regarding pillar failure in this mine, improving works such as roof support and replacing existing pillars with concrete pillars (CP) were carried out. During the construction of the second CP, in the space between the remaining pillars, one of the pillars failed leading to the progressive failure of other pillars until 4 000 m 2 of mine had collapsed within a few minutes. In this work, this phenomenon is described by applying both numerical and empirical methods and the respective results are compared. The results of numerical modelling are found to be closer to the actual condition than those of the empirical method. Also, a width-to-height (W/H) ratio less than 1, an inadequate support system and the absence of a detailed program for pillar recovery are shown to be the most important causes of the Domino failure in this mine.

Overhanging rock slope by design:An integrated approach using rock mass strength characterisation,large-scale numerical modelling and limit equilibrium methods

Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlights an integrated approach of designing overhanging rock slopes where the relative dimensions of the slope exceed the scale of fracturing and the rock mass failure needs to be considered rather than kinematic release of individual blocks. The key to the method is a simplified limit equilibrium(LE) tool that was used for the support design and analysis of a multi-faceted overhanging rock slope. The overhanging slopes required complex geometries with constantly changing orientations. The overhanging rock varied in height from 30 m to 66 m. Geomechanical modelling combined with discrete fracture network(DFN)representation of the rock mass was used to validate the rock mass strength assumptions and the failure mechanism assumed in the LE model. The advantage of the simplified LE method is that buttress and support design iterations(along with sensitivity analysis of design parameters) can be completed for various cross-sections along the proposed overhanging rock sections in an efficient manner, compared to the more time-intensive, sophisticated methods that were used for the initial validation. The method described presents the development of this design tool and assumptions made for a specific overhanging rock slope design. Other locations will have different geological conditions that can control the potential behaviour of rock slopes, however, the approach presented can be applied as a general guiding design principle for overhanging rock cut slope.

Numerical modelling of non-ideal detonation in ANFO explosives applying Wood-Kirkwood theory coupled with EXPLO5 thermochemical code

Ammonium nitrate and fuel oil(ANFO)based explosive is a classic example of non-ideal high explosives.Its detonation is characterized by a strong dependence of detonation parameters on explosive charge diameter,presence and characteristics of confinement,as well as incomplete consumption of explosive at the sonic point.In this work we propose a detonation model based on the Wood-Kirkwood(WK)theory coupled with the thermochemical code EXPLO5 and supplemented with reaction rate models.Our objective is to analyze the validity of the model for highly non-ideal ANFO explosives,with emphasis on effect of reaction rate models.It was found that both single-step and two-step pressure-based models can be calibrated to reproduce experimental detonation velocity-charge radius data of ANFO at radii significantly above the failure radius(i.e.for D/D_(id)>~0.6).Single-step pressure-based model,with the pressure exponent equal to 1.4,proved to be the most accurate,even in the vicinity of the failure radius.The impact of the rate models is most evident on temporal(and spatial)distribution of flow parameters in detonation driving zone,especially when it comes to the conversion and width of detonation driving zone.

Coupled Transfer of Water and Heat in Red Soil: Experiment and Numerical Modelling

Coupled transfer of soil water and heat in closed columns of homogeneous red soil was studied under laboratory conditions. A coupled model was constructed using soil physical theory, empirical equations and experimental data to predict the coupled transfer. The results show that transport of soil water was affected by temperature gradient, and the largest net water transport was found in the soil column with initial water content of 0.148 m3 m-3. At the same time, temperature changes with the transport of soil water was in a nonlinear shape as heat parameters were function of water content, and the changes of temperature were positively correlated with the net amount of water transported. Numerical modelling results show that the predicted values of temperature distribution were close to the observed values, while the predicted values of water content exhibited limited deviation at both ends of the soil column due to the slight temperature changes at both ends. It was indicated that the model proposed here was applicable.

Numerical Modellings of Properties of the Summer Quasi-Stationary Circulation Systems and Their Monthly Variations

An ocean-atmosphere and land-air coupled numerical model system is used to study the basic properties and the monthly time variations of the summer quasi-stationary circulation systems. It is found that either at the upper or at the lower levels of the atmosphere, the circulation patterns have a two-wave structure in the zonal direction at the mid and high latitudes of the Northern Hemisphere. Such a structure of circulation is totally matchable to that of the land-sea distribution there. It is proved, hence, that the land-sea distributive pattern is the fundamental cause for the summer quasi-stationary circulation pattern. The topography in the globe is the secondary factor for circulation systems. The circulation centres of the quasi-stationary systems are always located in certain areas due to the thermodynamic contrast between land and sea.From the time evolutions of the circulation systems it is seen that the change is larger at the beginning period of the time integration, it is because of using the zonally averaged mean fields as the initial values of the model. As long as the basic simulated pattern of circulations reaches the state similar to that of the real climatic fields resulting from the coeffects of the land-sea distribution and the topography, the circulation systems modelled will change slowly and tend to a quasi-stationary state. Therefore, the time integration does not need to last for a very long time, if the purpose of numerical modellings is to test sensitivities of some factors influencing the climate, 20 model days may be enough for sensitive experiments.

Cyclone and Tsunami Hazards in the Arabian Sea-A Numerical Modelling Case Study by Royal HaskoningDHV

Significant loss of life and damage to properties, ecosystems and marine facilities occur due to various natural hazards such as cyclones and tsunamis. Royal HaskoningDHV has developed regional hydrodynamic and wave models covering the Northern Arabian Sea to address these issues. Cyclone modelling was carried out on 11 major cyclones since 1945 and the tsunami modelling on an earthquake along the Makran Fault Line in 1945. Sample results from these modelling studies are presented in this paper. The methodology described in this article for modelling cyclones and tsunamis in the Arabian Sea could be applied to simulate these natural hazards at other sites around the world.

Numerical Modelling and Experimental Investigation of Microse-gregation in Al-4.45 wt pct Cu:Effect of Dendrite Joining

A numerical model Alloy/M has been applied to investigate the microsegregation in Al-4.45 wt pct Cu alloy. The calculated data were compared with the experimental measurements for samples solidified at different cooling rates. Discrepancies in solute concentration occur between the experimental observations and calculated results. Reasons for the discrepancies were discussed, and the effect of dendrite joining at a later stage of solidification on the microsegregation was investigated. Calculations that have included this effect showed a better fit with experimental results.

Numerical modelling of strength and energy release characteristics of pillar-scale coal mass

Coal burst is a manifestation of rapid energy release,which is considered as one of the most critical operational hazards in underground coal mines.This study numerically investigates the effects of discontinuities on the strength and energy release characteristics of coal mass samples under uniaxial compression.The universal distinct element code(UDEC)was used to model pillar-scale coal mass samples that were represented by an assembly of triangular deformable blocks,and pre-existing discontinuities such as bedding planes and cleats were also included in the models.It shows that cleat spacing can have a significant impact on compressive strength and energy release,with both strength and energy release(magnitude and rate)reducing as the number of cleats was increased.This work is one of the first attempts to numerically model and quantify the energy release which occurs during the failure of pillar-scale coal mass samples with varying cleat densities.The insights from the numerical modelling can help to understand the possible energy release mechanisms and associated coal burst potential in changing coal cleat conditions.

Numerical modelling of the tide - induced residual current in the East China Sea

-A two-dimensional.nonlinear numerical model is used to study the residual current generated by tides in the East China Sea (ECS)and the South Huanghai Sea (SHS). At first, the principal semidiurnal lunar tide (M2)and the tidal current are derived in these areas. The results obtained with the model are strongly supported by the observational results available. Then, the tide-induced residual flow is determined by using the currents generated by the tidal input. The main features of the residual current in ECS and SHS are presented by analyzing the calculated results. Some of the problems are discussed such as the cause of generating residual current and the contribution of the residual current to the observed current.

Numerical Modelling of Nitrogen Thermal Effects Produced by the Negative DC Corona Discharge

In this work, we simulate numerically the thermal effects in nitrogen at atmospheric pressure caused by a negative corona DC discharge of low current. A mathematical function that simulates the injection of the thermal transfer in gas is proposed. The simulated discharge is of a negative point to plane mass type, with an interelectrode distance of 12 mm and a symmetry about the axis of discharge. The spatial and temporal evolution of neutrals is analyzed based upon the equations of continuity, momentum and energy in a cylindrical geometry. For the geometry of the system, the FCT （flux corrected transport） procedure was adopted.

Experiments in numerical modelling of the Pacific sea surface temperature anomalies

By using the atmosphere-ocean coupled model (CGCM) which is composed of a 2-level global atmospheric general circulation model and a 4-layer Pacific oceanic general circulation model developed in the Institute of Atmospheric Physics of Chinese Academy of Sciences, and two model climatological fields got from the two independent models' numerical integrations respectively, the Pacific sea surface temperature anomalies (SSTA) from 1988 to 1989 are simulated in this paper with observed atmospheric general circulation data and sea surface temperature fields as initial conditions and monthly coupling scheme. In order to remove systematic biases of the model climatological fields, interaction variables between atmosphere and ocean are also corrected simultaneously. The experiments show that the simulation results can be improved effectively if these interaction variables are corrected in spite of the fact that there always exist systematic biases in independent numerical simulations of atmospheric part and oceanic part within CGCM. The basic characteristics of the observed Pacific SSTA in September and October 1988 have been simulated by using the correction scheme, such as the negative SSTA domain in the whole E-quatorial Pacific east to 150°E and the positive SSTA domain in the Western Pacific, the northern subtropical Pacific and nearly the whole Southern Pacific. Further numerical simulations show that the model can simulate not only the SSTA in the Pacific and its seasonal variations but also its interannual changes (for example, La Nino event in the Equatorial Pacific terminated after May 1989) to a certain degree. Furthermore, some problems existing in experiment processes and what we shoud do in the following stage are also discussed and analysed in this paper.

On Numerical Modelling of Industrial Powder Compaction Processes for Large Deformation of Endochronic Plasticity at Finite Strains

Compaction processes are one the most important par ts of powder forming technology. The main applications are focused on pieces for a utomotive, aeronautic, electric and electronic industries. The main goals of the compaction processes are to obtain a compact with the geometrical requirements, without cracks, and with a uniform distribution of density. Design of such proc esses consist, essentially, in determine the sequence and relative displacements of die and punches in order to achieve such goals. A.B. Khoei presented a gener al framework for the finite element simulation of powder forming processes based on the following aspects; a large displacement formulation, centred on a total and updated Lagrangian formulation; an adaptive finite element strategy based on error estimates and automatic remeshing techniques; a cap model based on a hard ening rule in modelling of the highly non-linear behaviour of material; and the use of an efficient contact algorithm in the context of an interface element fo rmulation. In these references, the non-linear behaviour of powder was adequately desc ribed by the cap plasticity model. However, it suffers from a serious deficiency when the stress-point reaches a yield surface. In the flow theory of plasticit y, the transition from an elastic state to an elasto-plastic state appears more or less abruptly. For powder material it is very difficult to define the locati on of yield surface, because there is no distinct transition from elastic to ela stic-plastic behaviour. Results of experimental test on some hard met al powder show that the plastic effects were begun immediately upon loading. In such mater ials the domain of the yield surface would collapse to a point, so making the di rection of plastic increment indeterminate, because all directions are normal to a point. Thus, the classical plasticity theory cannot deal with such materials and an advanced constitutive theory is necessary. In the present paper, the constitutive equations of powder materials will be discussed via an endochronic theory of plasticity. This theory provides a unifi ed point of view to describe the elastic-plastic behaviour of material since it places no requirement for a yield surface and a ’loading function’ to disting uish between loading an unloading. Endochronic theory of plasticity has been app lied to a number of metallic materials, concrete and sand, but to the knowledge of authors, no numerical scheme of the model has been applied to powder material . In the present paper, a new approach is developed based on an endochronic rate independent, density-dependent plasticity model for describing the isothermal deformation behavior of metal powder at low homologous temperature. Although the concept of yield surface has not been explicitly assumed in endochronic theory, it is shown that the cone-cap plasticity yield surface (Fig.1), which is the m ost commonly used plasticity models for describing the behavior of powder materi al can be easily derived as a special case of the proposed endochronic theory. Fig.1 Trace of cone-cap yield function on the meridian pl ane for different relative density As large deformation is observed in powder compaction process, a hypoelastic-pl astic formulation is developed in the context of finite deformation plasticity. Constitutive equations are stated in unrotated frame of reference that greatly s implifies endochronic constitutive relation in finite plasticity. Constitutive e quations of the endochronic theory and their numerical integration are establish ed and procedures for determining material parameters of the model are demonstra ted. Finally, the numerical schemes are examined for efficiency in the model ling of a tip shaped component, as shown in Fig.2. Fig.2 A shaped tip component. a) Geometry, boundary conditio n and finite element mesh; b) density distribution at final stage of

Numerical modelling of the sea level under the actions of the tide and strong wind in the Bohai Sea

In this paper, the numerical modelling of the tidal level and current in the Bohai Sea was carried out with ADI method, by taking the sum of four main tidal components M2,S2K2,O1 as the open boundary condition. The calculated values were consistent with the predicted ones (the observed values in the case of calm) in the Tidal Table. On the basis of the modelling of the tide, the sea level and current fields under the effects of strong wind were simulated. The calculated results were also quite satisfactory.

Application of Numerical Modelling to Optimise Design of a Coastguard Harbour

As part of the design process of the Sidab Coastguard harbour in Oman, 2-dimensional wave transformation modelling was carried out to generate inshore wave conditions （outside the harbour） using the MIKE21 Spectral Wave （SW） model of DHI （Danish Hydraulic Institute）. The MIKE21 Boussinesq Wave （BW） model was then used to simulate the wave penetration so that acceptable wave conditions were achieved inside the harbour. Tide levels and currents （speeds and directions） inside and outside the harbour were simulated using the MIKE21 Hydrodynamic （HD） model. Finally, the MIKE21 model was used to simulate the flushing mechanism to ensure good water quality inside the harbour. Both the existing and proposed harbour layouts were modelled for operational and extreme conditions. The influence of cyclones on the design was also considered. An acceptable harbour layout was determined which provided good shelter from waves and with an exchange rate of water between the harbour and the sea which was high enough to prevent eutrophication and to allow removal of man-made pollutants such as accidental oil spills. The paper compares the findings of the present study with a range of internationally used guidelines on harbour flushing. Finally, the paper demonstrates the value of using numerical modelling technique in optimising the layout of a new harbour.

Preliminary Assessment of Operational Downtime of Sea Ports and Berthing Facilities through Numerical Modelling

Downtime due to adverse wave conditions for vessels at berth is an important commercial aspect in the planning and development of a sea port or a berthing terminal. This paper describes a practical technique for preliminary assessment of operational downtime at a proposed bulk terminal. Time-series wind and wave data at an offshore location was purchased. Numerical modelling was then carried out using the MIKE21 SW （spectral wave） model developed by DHI （Deutsches Hydrographisches Institut） to transform these offshore waves into inshore in order to derive wave conditions at the berth. Both wind-waves and swell-waves were considered. Waves affecting the head and beam of a vessel were considered separately for a wide range of vessel sizes with DT （displacement tonnage） ranging from 5,000 tons to 〉 200,000 tons. The limiting wave height H5% was used. Operational downtime was also calculated using significant wave height, Hs as a criterion with limits ofHs = 1.0 m for beam seas and Hs = 1.5 m for head seas. The methodology and lessons learnt from the study can be applied for the development of any sea port worldwide.

The State Key Laboratory of Numerical Modelling for Atmospheric Science and Geophysical Fluid Dynamics

The State Key Laboratory of Numerical Modelling for Atmospheric Science and Geophysical Fluid Dynamics (LASG) was established in 1985. It was considered an Excellent Opening Laboratory in the national assessment organized by the former State Science and Technology Committee (now the Ministry of Science and Technology) in 1988. Because of this, LASG was upgraded to a State Key Laboratory in 1989. It won the honor of "Excellent

Numerical Modelling of Waves and Surge from Cyclone Mekunu (May 2018) in the Arabian Sea

Cyclone Mekunu developed in the Arabian Sea on 22 May 2018 and made landfall near the Port of Salalah(Oman)on 25 May.Wide spread damages to properties and coastal facilities and human casualties were reported in Yemen and Oman.Less information on numerical modelling of waves and surge is publicly available on this cyclone.Therefore,numerical modelling of Cyclone Mekunu was carried out in the present study to derive waves and storm surge.The MIKE21 Spectral Wave Model and the Flow Model were used in coupled mode to simulate the waves and surge from the cyclone.Model results of waves and surge are presented in this paper for illustration purposes.The methodology of the present study can be used to simulate any cyclone around the world.

Numerical modelling of the depth variation and the fluctuation of the thermocline under the effects of wind in the Bohai Sea

-In this paper, numerical modelling of the fluctuation of the thermocline in the Bohai Sea has been made using a two-dimensional nonlinear model in stratified ocean and the model for the depth of the thermocline under the effects of wind stirring. The computed results depict the variations of the fluctuation of the thermocline driven by different kinds of wind fields. The fluctuation of the thermocline in the Bohai Sea varies somewhat with different directions, paths and locations of typhoon (cyclone). Under the effects of strong wind, the thermoclines both sink due to mixing and fluctuate. Furthermore, the fluctuation of the thermocline speeds up mixing. At last, the thermoclines disappear after 12-15 h when the strong wind increases from Force 6 to Force 9.