Numerical analysis of geosynthetic-reinforced embankment performance under moving loads

The performance of geosynthetic-reinforced embankments under traffic moving loads is always a hotspot in the geotechnical engineering field.A three-dimensional(3D)model of a geosynthetic-reinforced embankment without drainage consolidation was established using the finite element software ABAQUS.In this model,the traffic loads were simulated by two moving loads of rectangular pattern,and their amplitude,range,and moving speed were realized by a Fortran subroutine.The embankment fill was simulated by an equivalent linear viscoelastic model,which can reflect its viscoelasticity.The geogrid was simulated by the truss element,and the geocell was simulated by the membrane element.Infinite elements were utilized to weaken the boundary effect caused by the model geometry at the boundaries.Validation of the established numerical model was conducted by comparing the predicted deformations in the cross-section of the geosynthetic-reinforced embankment with those from the existing literature.On this basis,the dynamic stress and strain distribution in the pavement structure layer of the geosynthetic-reinforced embankment under a moving load was also analyzed.Finally,a parametric study was conducted to examine the influences of the different types of reinforcement,overload,and the moving load velocity on the geosynthetic-reinforced embankment.

Numerical Analysis of the Influence of the Impinging Distance on the Scouring Efficiency of Submerged Jets

Water jet technology is widely used in submerged buried pipes as a non-traditional trenching process,often invol-ving a complex sediment response.An important adjustable and influential engineering variable in this technol-ogy is represented by the impinging distance.In this study,the FLOW-3D software was used to simulate the jet scouring of sand beds in a submerged environment.In particular,four sets of experimental conditions were con-sidered to discern the relationship between the maximum scour depth and mass and the impinging distance.As shown by the results,a critical impinging distance h0 exists by which the static scour depth can be maximized;the scour mass ratio between dynamic and static conditions decreases as the impinging distance increases.Moreover,the profile contours are similar when the erosion parameter Ec is in the range 0.35<Ec<2.Empirical equations applicable for predicting the jet trenching contour under both dynamic and static scour modes are also provided in this study.

Numerical analysis of high‑speed railway slab tracks using calibrated and validated 3D time‑domain modelling

Concrete slabs are widely used in modern railways to increase the inherent resilient quality of the tracks,provide safe and smooth rides,and reduce the maintenance frequency.In this paper,the elastic performance of a novel slab trackform for high-speed railways is investigated using three-dimensional finite element modelling in Abaqus.It is then compared to the performance of a ballasted track.First,slab and ballasted track models are developed to replicate the full-scale testing of track sections.Once the models are calibrated with the experimental results,the novel slab model is developed and compared against the calibrated slab track results.The slab and ballasted track models are then extended to create linear dynamic models,considering the track geodynamics,and simulating train passages at various speeds,for which the Ledsgard documented case was used to validate the models.Trains travelling at low and high speeds are analysed to investigate the track deflections and the wave propagation in the soil,considering the issues associated with critical speeds.Various train loading methods are discussed,and the most practical approach is retained and described.Moreover,correlations are made between the geotechnical parameters of modern high-speed rail and conventional standards.It is found that considering the same ground condition,the slab track deflections are considerably smaller than those of the ballasted track at high speeds,while they show similar behaviour at low speeds.

Numerical Simulation Analysis of Welded Joints in Arch Ribs of Large Span Steel Pipe Arch Bridges

In order to study the residual stress distribution law of welded joints of arch ribs of large-span steel pipe concrete arch bridges,numerical simulation of temperature,stress and strain fields based on ABAQUS for welded joints of arch-ribbed steel tubes using 7-,8-and 9-layer welds is carried out and its accuracy is demonstrated.The steel pipe welding temperature changes,residual stress distribution,different processes residual stress changes in the law,the prediction of post-weld residual stress distribution and deformation are studied in this paper.The results show that the temperature field values and test results are more consistent with the accuracy of numerical simulation of welding,the welding process is mainly in the form of heat transfer;Residual high stresses are predominantly distributed in the Fusion zone(FZ)and Heat-affected zone(HAZ),with residual stress levels tending to decrease from the center of the weld along the axial path,the maximum stress appears in the FZ and HAZ junction;The number of welding layers has an effect on the residual stress distribution,the number of welding layers increases,the residual stress tends to decrease,while the FZ and HAZ high stress area range shrinks;Increasing the number of plies will increase the amount of residual distortion.

Analysis of debris flow control effect and hazard assessment in Xinqiao Gully,Wenchuan M_(s)8.0 earthquake area based on numerical simulation

Xinqiao Gully is located in the area of the 2008 Wenchuan M_(s)8.0 earthquake in Sichuan province,China.Based on the investigation of the 2023"6-26"Xinqiao Gully debris flow event,this study assessed the effectiveness of the debris flow control project and evaluated the debris flow hazards.Through field investigation and numerical simulation methods,the indicators of flow intensity reduction rate and storage capacity fullness were proposed to quantify the effectiveness of the engineering measures in the debris flow event.The simulation results show that the debris flow control project reduced the flow intensity by41.05%to 64.61%.The storage capacity of the dam decreases gradually from upstream to the mouth of the gully,thus effectively intercepting and controlling the debris flow.By evaluating the debris flow of different recurrence intervals,further measures are recommended for managing debris flow events.

Numerical Modeling and Technico-Economic Analysis of a Hybrid Energy Production System for Self-Consumption: Case of Rural Area in the Comoros

This study aims to provide electricity to a remote village in the Union of Comoros that has been affected by energy problems for over 40 years. The study uses a 50 kW diesel generator, a 10 kW wind turbine, 1500 kW photovoltaic solar panels, a converter, and storage batteries as the proposed sources. The main objective of this study is to conduct a detailed analysis and optimization of a hybrid diesel and renewable energy system to meet the electricity demand of a remote area village of 800 to 1500 inhabitants located in the north of Ngazidja Island in Comoros. The study uses the Hybrid Optimization Model for Electric Renewable (HOMER) Pro to conduct simulations and optimize the analysis using meteorological data from Comoros. The results show that hybrid combination is more profitable in terms of margin on economic cost with a less expensive investment. With a diesel cost of $1/L, an average wind speed of 5.09 m/s and a solar irradiation value of 6.14 kWh/m2/day, the system works well with a proportion of renewable energy production of 99.44% with an emission quantity of 1311.407 kg/year. 99.2% of the production comes from renewable sources with an estimated energy surplus of 2,125,344 kWh/year with the cost of electricity (COE) estimated at $0.18/kWh, presenting a cost-effective alternative compared to current market rates. These results present better optimization of the used hybrid energy system, satisfying energy demand and reducing the environmental impact.

Numerical analysis of matrix swelling and its effect on microstructure of digital coal and its associated permeability during CO_(2)-ECBM process based on X-ray CT data

Matrix swelling effect will cause the change of microstructure of coal reservoir and its permeability,which is the key factor affecting the engineering effect of CO_(2)-ECBM technology.The Sihe and Yuwu collieries are taken as research objects.Firstly,visualization reconstruction of coal reservoir is realized.Secondly,the evolution of the pore/fracture structures under different swelling contents is discussed.Then,the influence of matrix phase with different swelling contents on permeability is discussed.Finally,the mechanism of swelling effect during the CO_(2)-ECBM process is further discussed.The results show that the intra-matrix pores and matrix-edge fractures are the focus of this study,and the contacting area between matrix and pore/fracture is the core area of matrix swelling.The number of matrix particles decreases with the increase of size,and the distribution of which is isolated with small size and interconnected with large size.The swelling effect of matrix particles with larger size has a great influence on the pore/fracture structures.The number of connected pores/fractures is limited and only interconnected in a certain direction.With the increase of matrix swelling content,the number,porosity,width,fractal dimension,surface area and volume of pores/fractures decrease,and their negative contribution to absolute permeability increases from 0.368% to 0.633% and 0.868%-1.404%,respectively.With the increase of swelling content,the number of intra-matrix pores gradually decreases and the pore radius becomes shorter during the CO_(2)-ECBM process.The matrix continuously expands to the connected fractures,and the width of connected fractures gradually shorten.Under the influence of matrix swelling,the bending degree of fluid flow increases gradually,so the resistance of fluid migration increases and the permeability gradually decreases.This study shows that the matrix swelling effect is the key factor affecting CBM recovery,and the application of this effect in CO_(2)-ECBM process can be discussed.

Numerical analysis of arc-droplet behavior in thin wire high current MAG welding with magnetic control

High efficiency welding has always been one of the hot research subjects at home and overseas. Aiming at the problem of unstable droplet rotating spray transfer process in 450 A thin wire metal active gas arc(MAG) welding, a three-dimensional arc-droplet coupled model is established in this paper on the base of fluid dynamics to conduct numerical analysis for droplet transfer behavior with external magnetic controlled MAG welding, and used high-speed camera system to verify the simulation results. The results indicate that without external magnetic field, the droplet is transferred by the way of rotating spray, the arc is changeable and the spatter is serious. The electromagnetic force accelerates the radial liquid beam rotation under direct current(DC) magnetic field, and increases the angle between liquid beam and welding wire axis. Meanwhile, the rotation frequency of the liquid beam reaches 400 Hz, the area of arc column and heat loss increase under DC magnetic field. Based on the alternating magnetic field, the arc is the “bell shaped”, the arc stiffness raises, its temperature achieves 26 000 K, and the maximum speed of the liquid stream beam can reach 500 m/s. And the droplet transfer is more stable, the spatter significantly reduces, and welding quality is effectively improved.

Numerical modeling and parametric sensitivity analysis of heat transfer and two-phase oil and water flow characteristics in horizontal and inclined flowlines using OpenFOAM

Estimating the oil-water temperatures in flowlines is challenging especially in deepwater and ultra-deepwater offshore applications where issues of flow assurance and dramatic heat transfer are likely to occur due to the temperature difference between the fluids and the surroundings. Heat transfer analysis is very important for the prediction and prevention of deposits in oil and water flowlines, which could impede the flow and give rise to huge financial losses. Therefore, a 3D mathematical model of oil-water Newtonian flow under non-isothermal conditions is established to explore the complex mechanisms of the two-phase oil-water transportation and heat transfer in different flowline inclinations. In this work, a non-isothermal two-phase flow model is first modified and then implemented in the InterFoam solver by introducing the energy equation using OpenFOAM® code. The Low Reynolds Number (LRN) k-ε turbulence model is utilized to resolve the turbulence phenomena within the oil and water mixtures. The flow patterns and the local heat transfer coefficients (HTC) for two-phase oil-water flow at different flowlines inclinations (0°, +4°, +7°) are validated by the experimental literature results and the relative errors are also compared. Global sensitivity analysis is then conducted to determine the effect of the different parameters on the performance of the produced two-phase hydrocarbon systems for effective subsea fluid transportation. Thereafter, HTC and flow patterns for oil-water flows at downward inclinations of 4°, and 7° can be predicted by the models. The velocity distribution, pressure gradient, liquid holdup, and temperature variation at the flowline cross-sections are simulated and analyzed in detail. Consequently, the numerical model can be generally applied to compute the global properties of the fluid and other operating parameters that are beneficial in the management of two-phase oil-water transportation.

Numerical Analysis for the Effect of Irresponsible Immigrants on HIV/AIDS Dynamics

The human immunodeficiency viruses are two species of Lentivirus that infect humans.Over time,they cause acquired immunodeficiency syndrome,a condition in which progressive immune system failure allows life-threatening opportunistic infections and cancers to thrive.Human immunodeficiency virus infection came from a type of chimpanzee in Central Africa.Studies show that immunodeficiency viruses may have jumped from chimpanzees to humans as far back as the late 1800s.Over decades,human immunodeficiency viruses slowly spread across Africa and later into other parts of the world.The Susceptible-Infected-Recovered(SIR)models are significant in studying disease dynamics.In this paper,we have studied the effect of irresponsible immigrants on HIV/AIDS dynamics by formulating and considering different methods.Euler,Runge Kutta,and a Non-standardfinite difference(NSFD)method are developed for the same problem.Numerical experiments are performed at disease-free and endemic equilibria points at different time step sizes‘ℎ’.The results reveal that,unlike Euler and Runge Kutta,which fail for large time step sizes,the proposed Non-standardfinite difference(NSFD)method gives a convergence solution for any time step size.Our proposed numerical method is bounded,dynamically con-sistent,and preserves the positivity of the continuous solution,which are essential requirements when modeling a prevalent disease.

Blast wave characteristics of multi-layer composite charge:Theoretical analysis,numerical simulation,and experimental validation

This article investigates the characteristics of shock wave overpressure generated by multi-layer composite charge under different detonation modes.Combining dimensional analysis and the explosion mechanism of the charge,a peak overpressure prediction model for the composite charge under singlepoint detonation and simultaneous detonation was established.The effects of the charge structure and initiation method on the overpressure field characteristics were investigated in AUTODYN simulation.The accuracy of the prediction model and the reliability of the numerical simulation method were subsequently verified in a series of static explosion experiments.The results reveal that the mass of the inner charge was the key factor determining the peak overpressure of the composite charge under single-point detonation.The peak overpressure in the radial direction improved apparently with an increase in the aspect ratio of the charge.The overpressure curves in the axial direction exhibited a multi-peak phenomenon,and the secondary peak overpressure even exceeded the primary peak at distances of 30D and 40D(where D is the charge diameter).The difference in peak overpressure among azimuth angles of 0-90°gradually decreased with an increase in the propagation distance of the shock wave.The coupled effect of the detonation energy of the inner and outer charge under simultaneous detonation improved the overpressure in both radial and axial directions.The difference in peak overpressure obtained from model prediction and experimental measurements was less than 16.4%.

Numerical Analysis on Temperature Distribution in a Single Cell of PEFC Operated at Higher Temperature by1D Heat Transfer Model and 3D Multi-Physics Simulation Model

This study is to understand the impact of operating conditions, especially initial operation temperature (Tini) which is set in a high temperature range, on the temperature profile of the interface between the polymer electrolyte membrane (PEM) and the catalyst layer at the cathode (i.e., the reaction surface) in a single cell of polymer electrolyte fuel cell (PEFC). A 1D multi-plate heat transfer model based on the temperature data of the separator measured using the thermograph in a power generation experiment was developed to evaluate the reaction surface temperature (Treact). In addition, to validate the proposed heat transfer model, Treact obtained from the model was compared with that from the 3D numerical simulation using CFD software COMSOL Multiphysics which solves the continuity equation, Brinkman equation, Maxwell-Stefan equation, Butler-Volmer equation as well as heat transfer equation. As a result, the temperature gap between the results obtained by 1D heat transfer model and those obtained by 3D numerical simulation is below approximately 0.5 K. The simulation results show the change in the molar concentration of O2 and H2O from the inlet to the outlet is more even with the increase in Tini due to the lower performance of O2 reduction reaction. The change in the current density from the inlet to the outlet is more even with the increase in Tini and the value of current density is smaller with the increase in Tini due to the increase in ohmic over-potential and concentration over-potential. It is revealed that the change in Treact from the inlet to the outlet is more even with the increase in Tini irrespective of heat transfer model. This is because the generated heat from the power generation is lower with the increase in Tini due to the lower performance of O2 reduction reaction.

Social-ecological perspective on the suicidal behaviour factors of early adolescents in China:a network analysis

Background In early adolescence,youth are highly prone to suicidal behaviours.Identifying modifiable risk factors during this critical phase is a priority to inform effective suicide prevention strategies.Aims To explore the risk and protective factors of suicidal behaviours(ie,suicidal ideation,plans and attempts)in early adolescence in China using a social-ecological perspective.Methods Using data from the cross-sectional project‘Healthy and Risky Behaviours Among Middle School Students in Anhui Province,China',stratified random cluster sampling was used to select 5724 middle school students who had completed self-report questionnaires in November 2020.Network analysis was employed to examine the correlates of suicidal ideation,plans and attempts at four levels,namely individual(sex,academic performance,serious physical llness/disability,history of self-harm,depression,impulsivity,sleep problems,resilience),family(family economic status,relationship with mother,relationship with father,family violence,childhood abuse,parental mental illness),school(relationship with teachers,relationship with classmates,school-bullying victimisation and perpetration)and social(social support,satisfaction with society).Results In total,37.9%,19.0%and 5.5%of the students reported suicidal ideation,plans and attempts in the past 6 months,respectively.The estimated network revealed that suicidal ideation,plans and attempts were collectively associated with a history of self-harm,sleep problems,childhood abuse,school bullying and victimisation.Centrality analysis indicated that the most influential nodes in the network were history of self-harm and childhood abuse.Notably,the network also showed unique correlates of suicidal ideation(sex,weight=0.60;impulsivity,weight=0.24;family violence,weight=0.17;relationship with teachers,weight=-0.03;school-bullying perpetration,weight=0.22),suicidal plans(social support,weight=-0.15)and suicidal attempts(relationship with mother,weight=-0.10;parental mental llness,weight=0.61).Conclusions This study identified the correlates of suicidal ideation,plans and attempts,and provided practical implications for suicide prevention for young adolescents in China.Firstly,this study highlighted the importance of joint interventions across multiple departments.Secondly,the common risk factors of suicidal ideation,plans and attempts were elucidated.Thirdly,this study proposed target interventions to address the unique influencing factors of suicidal ideation,plans and attempts.

Numerical analysis of anisotropic stiffness and strength for geomaterials

In numerical modelling,selection of the constitutive model is a critical factor in predicting the actual response of a geomaterial.The use of oversimplified or inadequate models may not be sufficient to reproduce the actual geomaterial behaviour.That selection is especially relevant in the case of aniso-tropic rocks,and particularly for shales and slates,whose behaviour may be affected,e.g.well stability in geothermal or oil and gas production operations.In this paper,an alternative anisotropic constitutive model has been implemented in the finite element method software CODE_BRIGHT,which is able to account for the anisotropy of shales and slates in terms of both deformability and strength.For this purpose,a transversely isotropic version of the generalised Hooke’s law is adopted to represent the stiffness anisotropy,while a nonuniform scaling of the stress tensor is introduced in the plastic model to represent the strength anisotropy.Furthermore,a detailed approach has been proposed to determine the model parameters based on the stressestrain results of laboratory tests.Moreover,numerical analyses are performed to model uniaxial and triaxial tests on Vaca Muerta shale,Bossier shale and slate from the northwest of Spain(NW Spain slate).The experimental data have been recovered from the literature in the case of the shale and,in the case of the slate,performed by the authors in terms of stress-strain curves and strengths.A good agreement can be generally observed between numerical and experi-mental results,hence showing the potential applicability of the approach to actual case studies.Therefore,the presented constitutive model may be a promising approach for analysing the anisotropic behaviour of rocks and its impact on well stability or other relevant geomechanical problems in aniso-tropic rocks.

A methodology for damage evaluation of underground tunnels subjected to static loading using numerical modeling

We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensities.To account for these variations,we utilized a Monte Carlo Simulation(MCS)technique coupled with the finite difference code FLAC^(3D),to conduct two thousand seven hundred numerical simulations of a horseshoe tunnel located within a rock mass with different geological strength index system(GSIs)and subjected to different states of static loading.To quantify the severity of damage within the rock mass,we selected one stress-based(brittle shear ratio(BSR))and one strain-based failure criterion(plastic damage index(PDI)).Based on these criteria,we then developed fragility curves.Additionally,we used mathematical approximation techniques to produce vulnerability functions that relate the probabilities of various damage states to loading intensities for different quality classes of blocky rock mass.The results indicated that the fragility curves we obtained could accurately depict the evolution of the inner and outer shell damage around the tunnel.Therefore,we have provided engineers with a tool that can predict levels of damages associated with different failure mechanisms based on variations in rock mass quality and in situ stress state.Our method is a numerically developed,multi-variate approach that can aid engineers in making informed decisions about the robustness of underground tunnels.

Experimental and Numerical Analysis of the Influence ofMicrochannel Size and Structure on Boiling Heat Transfer

Computational fluid dynamics was used and a numerical simulation analysis of boiling heat transfer in microchannels with three depths and three cross-sectional profiles was conducted.The heat transfer coefficient and bubble generation process of three microchannel structures with a width of 80μm and a depth of 40,60,and 80μm were compared during the boiling process,and the factors influencing bubble generation were studied.A visual test bench was built,and test substrates of different sizes were prepared using a micro-nano laser.During the test,the behavior characteristics of the bubbles on the boiling surface and the temperature change of the heated wall were collected with a high-speed camera and a temperature sensor.It was found that the microchannel with a depth of 80μm had the largest heat transfer coefficient and shortest bubble growth period,the rectangular channel had a larger peak heat transfer coefficient and a lower frequency of bubble occurrence,while the V-shaped channel had the shortest growth period,i.e.,the highest frequency of bubble occurrence,but its heat transfer coefficient was smaller than that of the rectangular channel.

Numerical manifold method for thermo-mechanical coupling simulation of fractured rock mass

As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.

Lagrangian coherent structure analysis on transport of Acetes chinensis along coast of Lianyungang,China

Spatial heterogeneity or“patchiness”of plankton distributions in the ocean has always been an attractive and challenging scientific issue to oceanographers.We focused on the accumulation and dynamic mechanism of the Acetes chinensis in the Lianyungang nearshore licensed fishing area.The Lagrangian frame approaches including the Lagrangian coherent structures theory,Lagrangian residual current,and Lagrangian particle-tracking model were applied to find the transport pathways and aggregation characteristics of Acetes chinensis.There exist some material transport pathways for Acetes chinensis passing through the licensed fishing area,and Acetes chinensis is easy to accumulate in the licensed fishing area.The main mechanism forming this distribution pattern is the local circulation induced by the nonlinear interaction of topography and tidal flow.Both the Lagrangian coherent structure analysis and the particle trajectory tracking indicate that Acetes chinensis in the licensed fishing area come from the nearshore estuary.This work contributed to the adjustment of licensed fishing area and the efficient utilization of fishery resources.

Seismic Wavelet Analysis Based on Finite Element Numerical Simulation

The practice of exploration and production has proved that explosives are excited in different surrounding rocks and the seismic wavelets collected have different characteristics. In this paper, by establishing a numerical model of the explosion in the well, using finite element analysis technology for numerical simulation, the simulation calculated the stress structure in the near-source area of the earthquake excitation, and extracted the seismic wavelet. The results show that the simulation seismic wavelet characteristics of different thin interbedded sand and mudstone structures have changed significantly. Through excitation simulation, the amplitude and spectrum information of seismic wavelets can be compared and analyzed, and the excitation parameters can be optimized. .

Numerical analysis of the effects of vesicle distribution characteristics on the engineering properties of volcanic rocks

Vesicles can be of different sizes and shapes and can be randomly distributed within vesicular volcanic rocks. This study investigates the variation of engineering properties of vesicular rocks due to the changes in vesicle distribution characteristics for different cases of bulk porosity and vesicle diameter using a systematic numerical simulation program using the finite element method-based rock failure process analysis (RFPA) software. Models with uniform-size vesicles and combinations of different proportions of different-sized vesicles were considered to resemble natural vesicular rocks more closely, and ten different random vesicle distributions were tested for each case. Increasing bulk porosity decreased the uniaxial compressive strength (UCS) and elastic modulus of the specimens, and the specimens with the lowest bulk porosity showed the greatest range of UCS values in the case of uniform-size vesicles. The effect of vesicle diameter on UCS showed an unsystematic response which was understood to be a result of different vesicle distribution patterns, some of which facilitated a shear failure. Specimens with multiple-size vesicles in different proportions revealed that the variation of UCS due to vesicle distribution characteristics is minimum when the bulk porosity is equally shared by different size vesicles. In addition, when the proportion of smaller-sized vesicles is higher, UCS showed an increase compared to that of the equal proportion of different size vesicles case at low porosities, but a decrease at higher porosities. Variation of elastic modulus showed minor, unsystematic fluctuations as a function of vesicle diameter and different proportions of different-sized vesicles, and the range for different vesicle distribution patterns was narrow in general. Overall, the findings of this study recommend cautious use of the engineering properties determined through a limited number of laboratory tests on vesicular rocks.