Documents published in 2021 : Scipedia

A Method for Accurate and Efficient Propagation of Satellite Orbits: A Case Study for a Molniya Orbit
R. Flores, M. Burhani

Alexandria Engineering Journal (2021). Alexandria Engineering Journal (2021) 60, 2661–2676

Abstract
Fast and precise propagation of satellite orbits is required for mission design, orbit determination and payload data analysis. We present a method to improve the computational performance of numerical propagators and simultaneously maintain the accuracy level required by any particular application. This is achieved by determining the positional accuracy needed and the corresponding acceptable error in acceleration on the basis of the mission requirements, removing those perturbation forces whose effect is negligible compared to the accuracy requirement, implementing an efficient and precise algorithm for the harmonic synthesis of the geopotential gradient (i.e., the gravitational acceleration) and adjusting the tolerance of the numerical propagator to achieve the prescribed accuracy level with minimum cost. In particular, to achieve the optimum balance between accuracy and computational performance, the number of geopotential spherical harmonics to retain is adjusted during the integration on the basis of the accuracy requirement. The contribution of high-order harmonics decays rapidly with altitude, so the minimum expansion degree meeting the target accuracy decreases with height. The optimum degree for each altitude is determined by making the truncation error of the harmonic synthesis equal to the admissible acceleration error. This paper presents a detailed description of the technique and test cases highlighting its accuracy and efficiency.
Abstract
Fast and precise propagation of satellite orbits is required for mission design, orbit determination and payload data analysis. We present a method to improve the computational [...]
- 2
- read
BEHAVIOUR OF A CRACKED GAS PIPELINE
J. Massa, A. Giudici

Revista Internacional de Desastres Naturales, Accidentes e Infraestructura Civil (2009). Rev. Int. de Desastres Naturales, Accidentes e Infraestructura Civil. Vol. 9(1-2)

Abstract
When designing systems for the transportation and distribution of gas, safety is one of the main issues to be considered. In order to minimize potential hazards, this subject is treated in international regulations. Manufacturing characteristics of gas pipes and their typical defects are described in the first part of this work. Safety concerns related to the presence of real and potential crack type defects in gas pipes are presented. An actual gas pipeline, which has been feeding natural gas to the city of Córdoba (Argentina) for more than forty years, is considered as a study case. In order to determine the properties of this gas pipeline, mechanical tests from material samples were carried out. The stress state associated with “plane” type defects is determined using a model that is based on the material toughness of the pipe, which is appropriate for studying the behavior of cracks. This model allows computing, depending on the dimensions of the pipe, the transition pressure which separates possible defects into two categories: defects which cause failure by breakage of the pipe and defects which only cause gas leaks. The pipe failure pressure for passing and non-passing defects can also be determined using this model. The system safety conditions are defined considering the size of the defects and the work pressure. The range of lengths of possible defects, the size of the defects that are critical, and the size of the defects which are tolerable for a fixed safety coefficient, are determined. Finally, the problem of determining in a quick way, the reduced pressure acceptable for the safe operation until repairs, of a system under the presence of a crack of “non-tolerable size”, is addressed.
Abstract
When designing systems for the transportation and distribution of gas, safety is one of the main issues to be considered. In order to minimize potential hazards, this subject is [...]
- 0
- read
DAMAGE BY CORROSION IN PIPELINES
J. Massa, A. Giudici

Revista Internacional de Desastres Naturales, Accidentes e Infraestructura Civil (2010). Rev. Int. de Desastres Naturales, Accidentes e Infraestructura Civil. Vol. 10(2)

Abstract
When designing systems for the transportation and distribution of gas, safety is one of the main issues to be considered. In order to minimize potential hazards, this subject is treated in international regulations. Safety concerns related to the presence of real and potential corrosion type defects in gas pipes are presented. An actual gas pipeline, which has been feeding natural gas to the city of Córdoba (Argentina) for more than forty years, is considered as a study case. In order to determine the properties of this gas pipeline, mechanical tests from material samples were carried out. The stress state associated with “volumetric” type defects is determined using a model that is based on the flux stress of the material of the pipe, which is appropriate for studying the behavior of corroded pipelines. This model allows computing, depending on the dimensions of the pipe and the length of the defect, the transition pressure which separates possible defects into two categories: those which cause failure by breakage of the pipe and those which only cause gas leaks. The pipe failure pressure for passing and non-passing defects can also be determined using this model. The system safety conditions are defined considering the size of the defects and the work pressure. The range of lengths of possible defects, the size of the defects that are critical, and the size of the defects which are tolerable for a fixed safety coefficient, are determined. Finally, the problem of determining in a quick way, the reduced pressure acceptable for the safe operation until repairs, of a system under the presence of flaws of “non-tolerable size”, is addressed.
Abstract
When designing systems for the transportation and distribution of gas, safety is one of the main issues to be considered. In order to minimize potential hazards, this [...]
- 0
- read
Effect of Covid-19 Pandemic in Construction Labor Productivity: A Quantitative and Qualitative Data Analysis
E. Quezon, A. Ibanez

American Journal of Civil Engineering and Architecture (2021). vol. 9, no. 1 (2021): 23-33. doi: 10.12691/ajcea-9-1-4.

Abstract
This research aims to identify and analyze the various factors affecting construction labor productivity covering the period from December 9, 2020, to January 31, 2021, a year after it was declared Covid-19 a pandemic. Due to the pandemic effect, the local government units in four selected study areas imposed intermittent Enhance Community Quarantine on all places to control the Coronavirus spread in compliance with the health department protocol. Fifty-five questionnaires returned out of 63 distributed are tabulated according to each group-related factor. The correlation analysis resulted in the highest coefficient value of 0.89 between the CTDEO and contractor groups. Meaning, most respondents have the same perception of the factors affecting construction's low labor productivity. Besides, results depicted that the absence of health workers on the construction site ranked 1st in the health & safety provision factor group with an RII of 0.97, followed by no safety engineers on the construction project sites with an RII of 0.93. From all the seven groups of factors, the health and safety group-related factors ranked 1st with an RII of 0.81, next, the schedule compression group with an RII of 0.78. Hypothesis testing asserted that working six days per week was one of the significant factors affecting labor productivity on the contractor's side, suggested by more than 50% of the respondents. The workforce group-related factors with an RII of 0.77 demonstrated a slight difference with the schedule compression group-related factors. Thus, the Covid-19 pandemic has a significant effect on the essential factors affecting construction's low labor productivity.
Abstract
This research aims to identify and analyze the various factors affecting construction labor productivity covering the period from December 9, 2020, to January 31, 2021, [...]
- 2
- read
Accurate Timoshenko beam elements for linear elastostatics and LPB stability
C. Felippa, E. Oñate

Archives of Computational Methods in Engineering (2021). Online First

Abstract
Several methods to derive accurate Timoshenko beam finite elements are presented and compared. Two application problems are examined: linear elastostatics and linearized prebuckling (LPB) stability analysis. Accurate elements can be derived for both problems using a well known technique that long preceeds the Finite Element Method: using homogeneous solutions of the governing equations as shape functions. An interesting question is: can accurate elements be derived with simpler assumptions? In particular, can linear-linear interpolation of displacements and rotations with one-point integration reproduce those elements? The answers are: no if standard variational tools based on classical functionals are used, but yes if modified functionals are introduced. The connection of modified functionals to newer methods, in particular templates, modified differential equations and Finite Increment Calculus (FIC) are examined. The results brings closure to a 50-year conumdrum centered on this particular finite element model. In addition, the discovery of modified functionals provides motivation for extending these methods to full geometrically nonlinear analysis while still using inexpensive numerical integration.
Abstract
Several methods to derive accurate Timoshenko beam finite elements are presented and compared. Two application problems are examined: linear elastostatics and linearized prebuckling [...]
- 1
- read
Mathematical optimization problems for particle finite element analysis applied to 2D landslide modeling
L. Wang, X. Zhang, F. Zaniboni, E. Oñate, S. Tinti

Mathematical Geosciences (2021). Vol. 53, pp. 81–103

Abstract
Notwithstanding its complexity in terms of numerical implementation and limitations in coping with problems involving extreme deformation, the finite element method (FEM) offers the advantage of solving complicated mathematical problems with diverse boundary conditions. Recently, a version of the particle finite element method (PFEM) was proposed for analyzing large-deformation problems. In this version of the PFEM, the finite element formulation, which was recast as a standard optimization problem and resolved efficiently using advanced optimization engines, was adopted for incremental analysis whilst the idea of particle approaches was employed to tackle mesh issues resulting from the large deformations. In this paper, the numerical implementation of this version of PFEM is detailed, revealing some key numerical aspects that are distinct from the conventional FEM, such as the solution strategy, imposition of displacement boundary conditions, and treatment of contacts. Additionally, the correctness and robustness of this version of PFEM in conducting failure and post-failure analyses of landslides are demonstrated via a stability analysis of a typical slope and a case study on the 2008 Tangjiashan landslide, China. Comparative studies between the results of the PFEM simulations and available data are performed qualitatively as well as quantitatively.
Abstract
Notwithstanding its complexity in terms of numerical implementation and limitations in coping with problems involving extreme deformation, the finite element method (FEM) [...]
- 2
- read
A nodal-integration based particle finite element method (N-PFEM) to model cliff recession
J. Meng, X. Zhang, S. Utili, E. Oñate

Geomorphology (2021). Vol. 381, pp. 107666

Abstract
Cliff recession poses a significant threat to the built environment, transportation infrastructure and land use. In this paper, a novel computational framework called the Nodal-integration based Particle Finite Element Method (N-PFEM) is developed for modelling the cliff recession resulting from weathering-induced landslides. The N-PFEM combines the nodal-integration technique with the PFEM in second-order cone programming and thus requires no variable mapping operation, which is essential in the classical PFEM for modelling history-dependent materials, for modelling large deformation problems such as landslides in cliff recession processes. To verify the developed N-PFEM, a series of benchmarks have been simulated including the cliff recession under both the weathering-limited and transport-limited conditions. Simulation results from the N-PFEM are validated in detail to these from the limit analysis method, well established geomorphologic models and the discrete element method. Additionally, measures for preventing cliff recession such as the construction of retaining wall structures are also investigated using the N-PFEM.
Abstract
Cliff recession poses a significant threat to the built environment, transportation infrastructure and land use. In this paper, a novel computational framework called the [...]
- 2
- read
A review of graphene reinforced Cu matrix composites for thermal management of smart electronics
S. Ali, F. Ahmad, P. Yusoff, N. Muhamad, E. Oñate, M. Raza, K. Malik

Composites Part A: Applied Science and Manufacturing (2021). Vol. 144, pp. 106357

Abstract
Heat dissipation remains a key challenge to be addressed, determining the performance and durability of smart electronic devices. Graphene reinforced metal matrix composites have been extensively studied as a thermal management material due to their high thermal conductivity and low coefficient of thermal expansion. The emphasis of this review is pivoted on the thermal conductivity enhancement of graphene reinforced Cu matrix composites developed in the recent literature. An overview of factors affecting thermal conductivity of composite namely defect processing route, density, graphene derivative, lateral size, concentration, alignment, graphene/matrix interfacial bonding and graphene modification are discussed. An extensive weightage is given to the processing route as it is the most influential factor in determining the enhancement efficiency. Furthermore, graphene based functional products such as heat spreader and heat sink developed for heat dissipation of electronic devices are also reviewed. Finally, the development and outlook for graphene based Cu composites are presented.
Abstract
Heat dissipation remains a key challenge to be addressed, determining the performance and durability of smart electronic devices. Graphene reinforced metal matrix composites [...]
- 1
- read
Methodological-technological framework for construction 4.0
F. Muñoz, J. Mora-Serrano, I. Valero, E. Oñate

Archives of Computational Methods in Engineering (2021). Vol. 28, pp. 689–711

Abstract
The construction industry has traditionally been characterised by the high diversity of its agents and processes, high resistance to change and low incorporation of technology compared to manufacturing industries. However, the construction sector is experiencing now a strong renovation process in methodology and tools due to the incorporation of the Building Information Modelling, Lean Construction and Integrated Project Delivery. Meanwhile, in production systems, “Industry 4.0” is a new paradigm that proposes automation, monitoring, sensorisation, robotisation, and digitalisation to improve production and distribution processes. In this context, some authors have proposed the concept of “Construction 4.0” as the counterpart of Industry 4.0 for the construction sector, although the methodological-technological implications are not clear. This research shows a methodological-technological framework adapted to the Architecture, Engineering, Construction, and Operations industry. This papers includes a detailed proposal for a reference frameworks and related technologies that could impact on this sector, responding to its complexities and specific challenges, such as the unique spaces for each work, which are difficult to standardise, arbitrary cost overruns and a productivity far below the average for other industries, increasing competitiveness and globalisation, as opposed to its traditionally local deployment, and an increasing demand to reduce the carbon footprint for all its activities.
Abstract
The construction industry has traditionally been characterised by the high diversity of its agents and processes, high resistance to change and low incorporation of technology [...]
- 4
- read
Evaluation of the HMA performance using combined effects of marble waste dust filler and superpave aggregate gradation
A. Mulatu, E. Quezon, O. Gudina

Australian Journal of Civil Engineering (Select a year). Adamu Mulatu , Emer Tucay Quezon & Oluma Gudina (2020): Evaluation of the HMA performance using combined effects of marble waste dust filler and superpave aggregate gradation, Australian Journal of Civil Engineering. 12(1); 1-11, Doi: 10.1080/14488353.2020.1854942

Abstract
During the fabrication of new products, a variety of wastes discharged and stored by causing harsh environmental problems. The marble processing plant is amongst the industrial sector that produces huge waste. Thus, the main purpose of this study is to evaluate the combined effects of marble waste dust filler and Superpave aggregate gradation on the performance of HMA. Based on the Marshall Mix design, seventy-two (72) specimens were produced from three below-restricted zones (BRZ) of Superpave gradation (BRZ4, BRZ5, & BRZ6) by using the 0.45 power chart. The design gradation obtained to be BRZ5 with 4.96% OBC with corresponding values 9.13 KN of stability, 2.96 mm of flow, 73.51% of VFA, 15.73% of VMA, 2.308 gm/cm3, of bulk density, and 4% of air void. The crushed stone dust filler was partially replaced by Marble Waste Dust filler (under No. 200 sieve) at a replacement rate of (1%, 1.5, 2, 2.5, 3, 3.5, 4, 4.5, and 5%) by keeping design gradation and OBC constant. The test showed satisfactory results. The Marshall mix design parameters of all mixtures are within the specification limit under international and local specifications. Marble waste dust filler when blended with a BRZ of Superpave gradation indicated a substantial effect on Marshall stability with 10.18 KN value at 3% MWD filler content.
Abstract
During the fabrication of new products, a variety of wastes discharged and stored by causing harsh environmental problems. The marble processing plant is amongst the industrial [...]
- 0
- read
Evaluation of waste engine oil rejuvenation for highly short term aged asphalt binder
B. Tadele, E. Quezon

Australian Journal of Civil Engineering (2021). Biruk Tadele & Emer Tucay Quezon (2021): Evaluation of waste engine oil rejuvenation for highly short term aged asphalt binder, Australian Journal of Civil Engineering, DOI: 10.1080/14488353.2021.1896124

Abstract
The aging property of the asphalt binder is time-dependent. A time-dependent short-term aging property of 80/100 penetration grade bitumen and the possibility of waste engine oil rejuvenation for the highly aged asphalt binder performance were investigated. Four specimens of equal weight from the penetration grade-80/100 bitumen are collected. The first specimen was checked for quality requirements. The other three specimens were aged using a rolling thin film oven for 85, 115, and 145 minutes to simulate the delay during hot mix asphalt production, hauling, and compaction. The highly aged bitumen was rejuvenated with 2%, 5%, and 10% Waste Engine Oil by weight. Results indicated that as the aging time increased, penetration and ductility decreased, softening point, flash point, fire point, and mass loss increased. A conventional test showed that highly aged bitumen from the trial period was 145 minutes and 10% waste engine oil obtained the optimum dosage. Further, multiple stress creeps recovery analyses indicated the rejuvenated binder is prone to pavement rutting above 70°C, and rejuvenation is effective for pavement temperature below 70°C. Hence, exposing the asphalt binder for temperature for a more extended period affects pavement performance.
Abstract
The aging property of the asphalt binder is time-dependent. A time-dependent short-term aging property of 80/100 penetration grade bitumen and the possibility of waste engine [...]
- 0
- read
Statistical Analysis on Shear Strength Parameter from Index Properties of Fine-grained Soils
T. Bakala, E. Quezon, M. Yasin, S. contributions

Journal of Engineering Research and Reports (2021). Tariku Tafari Bakala, Emer Tucay Quezon and Mohammed Yasin. Statistical Analysis on Shear Strength Parameter from Index Properties of Fine-grained Soils. Journal of Engineering Research and Reports. 2021; 20(4): 15-28. DOI: 10.9734/JERR/2021/v20i417291

Abstract
Shear strength is the essential engineering property of soil required to analyze and design foundations, retaining walls, bridges, embankment, and related infrastructure. The laboratory equipment and field instruments are not sufficient in developing countries to obtain soil engineering properties, especially strength properties. Thus, Geotechnical engineers usually endeavor to develop statistical models that best fit a particular area and soil type, especially for analysis and design purposes. In this research, a Statistical Analysis on the Shear Strength parameter from the Index Properties of Fine-Grained Soils was studied. For predicting the undrained shear strength parameter, single linear regression (SLR) and multiple linear regressions (MLR) analyses were developed. To develop the intended statistical models for a study, SAS JMP Pro 13, SPSS v22, and Microsoft Excel-2013 software were introduced. The results of a study indicated that undrained shear strength(Cu) was significantly correlated with liquid limit(LL), plastic limit(PL), bulk density (ρbulk), dry density(ρdry), natural moisture content(NMC), and plasticity index(PI). While it was not significantly correlated with a specific gravity (Gs) and liquidity index (LI) of study area soil. Finally, a strong Model of Cu with a coefficient of determination (R2 = 0.806), good significance level, and less Std. error was obtained from multiple linear regression (MLR) analysis. The developed model can figure undrained shear strength parameter and wide application in the construction industry to minimize the cost, effort, and time for laboratory tests of shear strength parameter of a study area.
Abstract
Shear strength is the essential engineering property of soil required to analyze and design foundations, retaining walls, bridges, embankment, and related infrastructure. [...]
- 0
- read
Model-Based Urban Road Network Performance Measurement Using Travel Time Reliability: A Case Study of Addis Ababa City, Ethiopia
S. Assen, E. Quezon

American Journal of Civil Engineering and Architecture (2019). Seada Mohammed Assen, and Emer Tucay Quezon, “Model-Based Urban Road Network Performance Measurement Using Travel Time Reliability: A Case Study of Addis Ababa City, Ethiopia.” American Journal of Civil Engineering and Architecture, vol. 7, no. 5 (2019): 202-207. doi: 10.12691/ajcea-7-5-2.

Abstract
In heterogeneous traffic conditions, the performance of the road network is described by vehicles and driver characteristics. Nowadays, traffic congestion, delay, and unreliability are terms that are most associated with present-day travel, in which transport users spend their precious time on long traffic queues. Due to this, late arrival at the workplace and appointment for social or business activities have become a perennial problem in the study area. In addition, during traffic queues, vehicle fuel emission increases in congested traffic segments affecting the environment, particularly the issue of global warming. This research study evaluated the performance of the road network in terms of travel time reliability in order to determine the main factors affecting travel time reliability. Ten road segments were selected to analyze the performance level and efficiency of the road network, considering the travel time probability distribution, and reliability of road segments in Addis Ababa City. From these road segments, nine were selected to formulate the model, and one road segment considered for validation of the result. The reliability of the road segments was analyzed using travel time reliability measures such as buffer time, buffer time index, planning time index, and the travel time index. It was used a multiple linear regression model to predict the travel time reliability of the road segments with R2
Abstract
In heterogeneous traffic conditions, the performance of the road network is described by vehicles and driver characteristics. Nowadays, traffic congestion, delay, and unreliability [...]
- 1
- read
Performance Evaluation of Grade Separation Within the Selected Major Intersection in Addis Ababa City
S. Leleisa, E. Quezon, T. Wedajo

American Journal of Traffic and Transportation Engineering (2018). Sentayehu Leleisa, Emer Tucay Quezon, Teyba Wedajo. Performance Evaluation of Grade Separation Within the Selected Major Intersection in Addis Ababa City. American Journal of Traffic and Transportation Engineering. Vol. 3, No. 3, 2018, pp. 41-49. doi: 10.11648/j.ajtte.20180303.11

Abstract
Urban City Traffic Management introduced various remedies to decongest traffic and to reduce the severity of accidents, such as providing traffic signal controls at the intersecting roads at different levels, or grade separating the traffic movement towards different directions. One type of grade separation is Road Underpass or Flyover. The construction of Road Underpass or Flyover poses some doubts whether it is effective or ineffective from the point of view of the travelers or motorists. It is for this reason that this research sought to evaluate the performance and effectiveness of the Road Underpass structure on traffic flow at the major intersection, specifically the Wellosefer intersection in Addis Ababa City. The primary data used the existing traffic volume recorded between 7:00 AM – 6:00 PM from Monday to Friday. While the secondary data obtained from the Consulting Office 'before' the improvement of Wellosefer intersection, including the previous study, which served as the control analysis for the research. From these data, the analyses provided results on vehicle travel time and delay for a case 'before' and 'after' improvement of the intersection. SIDRA Software was used to analyze the traffic conditions of the intersection. Based on the findings of the study, a case 'before' Improvement of Wellosefer Unsignalized T-Intersection indicated an average delay of vehicles of about 535.7 seconds and a travel time of vehicles of about 574.7 seconds. On the other hand, in a case 'after' Improvement of Wellosefer intersection, three alternatives considered; such as 1) Signalized T-Intersection, 2) Roundabout and 3) Road Underpass with Roundabout at the upper level. The results of the first alternative showed the Signalized T-intersection indicated an average delay and travel time of 344.3 seconds and 365.0 seconds, respectively. It means, there could be a reduction of 15.4% in average delay if signalized intersection implemented. The second alternative, a case of Roundabout using the same traffic data at current condition, the results indicated the motorists and commuters would be experiencing an average delay of about 355.8 seconds and average travel time of 390.0 seconds. The third alternative considering Road Underpass with Roundabout at the upper level, the results revealed that there was a significant reduction in an average delay of vehicles of about 42% from the Unsignalized T-Intersection. Comparing this result with the Signalized intersection using the same traffic data, showed that there was a significant decrease in average delay. Therefore, the findings of the study, a combination of Road Underpass and Roundabout at the upper level as in the case of Wellosefer intersection performed well than the other traffic control design alternatives.
Abstract
Urban City Traffic Management introduced various remedies to decongest traffic and to reduce the severity of accidents, such as providing traffic signal controls at the intersecting [...]
- 1
- read
Traffic Safety Evaluation at Unsignalized road Intersection using Proactive Approach
. Abstract, E. Quezon

Proceedings of the 2nd Annual Ethiopian Road Safety (RSC2018) (2018). Seifeddin Haji, Prof. Emer T. Quezon, Murad Mohammed, Biruk Yigezu. Traffic Safety Evaluation on Urban Road Un-signalized Intersection using Proactive Approach: A Case Study in Adama City. Proceedings of the 2nd Annual Ethiopian Road Safety (RSC2018) Addis Ababa, Ethiopia, August 23-24, 2018

Abstract
Deaths and injuries as a result of road accidents are a recognized global problem and authorities are aware of the consequences this global phenomenon might lead to. Because of the rapid development of urbanization and motorization in developing and transitional countries, like Ethiopia, actions to improve traffic safety conditions are urgently needed. Highway intersections are among the greatest safety challenges because they are some of the most common highway features. Improving the safety level of the urban intersection is an important way to increase safety. Safety evaluations can help to develop effective safety countermeasures to lower crash rates and reduce crash severity. So how to diagnose safety problems and improve urban road intersection safety effectively becomes an important issue. Now, there have been quite a number of safety analysis methods and theories in the city intersection, but these methods and theories are formed on the basis of traffic accidents. However, the data on traffic accidents is not accurate enough, these theories are not appropriate for our country. Also, it is impractical and unethical to wait for accidents to occur before being able to draw statistically sound conclusions regarding safety improvement. Hence there is a need to develop a proactive approach, a non-accident-based approach, to evaluate urban road Unsignalized intersection safety performance. A proactive approach (indirect and non-accident based approach) was attempted to that of traditional reactive approach (direct and accident based approach) which is purely based on the existing conditions that have direct relationships or impacts to traffic safety such as the geometric, traffic control & engineering, road Surface and environment characteristics. The non-accident-based approach is based on field surveys under the conditions mentioned previously, summarizes the intersection safety diagnosis in a safety index to indicate the safety performance of the intersection. The study indicates that the proposed methodology and intersection safety diagnosis technique can be used to perform an evaluation of traffic safety and improvement at an intersection even when there are few resources available. Also, the use of intersection safety diagnosis is a more resource-efficient and ethically appealing alternative for fast, reliable, and effective safety assessment.
Abstract
Deaths and injuries as a result of road accidents are a recognized global problem and authorities are aware of the consequences this global phenomenon might lead to. Because [...]
- 1
- read
Investigation of Pedestrian Safety Problems and Its Countermeasures: A Case Study in Nekemte City, Ethiopia
H. Alemgena, P. Quezon, T. Kumela, E. Quezon

International Journal of Engineering and Technical Research (2018). Alemgena, H.; Quezon, ET.; Kumela, T. Investigation of Pedestrian Safety Problems and Its Countermeasures: A Case Study in Nekemte City, Ethiopia. International Journal of Engineering and Technical Research. 2018; 8(2): 58-64

Abstract
The actual growth of urban economic activities and simple movements of peoples and goods measured through its primary transport system. The requirements of the pedestrian should be considered in the design of the urban environment and transportation facilities. Since road infrastructure improvements and designed pedestrian environments are an essential part of road safety enhancements that increase walking and decrease fatalities. This study investigated the pedestrian safety problems and its countermeasures on road segments in Nekemte city. This study examined the existing issues affecting the safety of pedestrian movement, evaluates the service quality of walking facilities, and finally identified the most significant factors affecting pedestrian safety with engineering countermeasures in the city. For the data analysis, the study used both objective measurements and subjective assessments to model PLOS using Australian methods. As per field study indicates 60% of the respondents' think the sidewalk is not user-friendly. In the city, the built road infrastructure is lacking pedestrian facility management and appropriate maintenance activities. Since the walkways in Nekemte city, in general, have almost the same problems like non-continuous, crack, and non-covered utility holes which are difficult especially for pedestrians with disabilities. The current sidewalk is encroached either by vendors, utility poles, parking carts or illegal construction materials in different places, and crosswalk facilities are not visible and not enough provided at a logical distance. The Pedestrian LOS analysis indicates that poor pedestrian conditions exist and the factors that negatively affect pedestrian LOS are wide-ranging in this city road segment. It concluded that Pedestrian facilities and all the concerned factors should be provided and maintained to more efficient for both encouraging people to walk and improve pedestrian safety.
Abstract
The actual growth of urban economic activities and simple movements of peoples and goods measured through its primary transport system. The requirements of the pedestrian [...]
- 1
- read
Analysis of StressStrain and Deflection of Flexible Pavements Using Finite Element Method Case Study on Bako-Nekemte Road
S. Wayessa, E. Quezon, T. Kumela

Journal of Civil, Construction and Environmental Engineering (2017). Shiferaw Garoma Wayessa, Emer Tucay Quezon, Tarekegn Kumela. Analysis of Stress- Strain and Deflection of Flexible Pavements Using Finite Element Method Case Study on Bako-Nekemte Road. Journal of Civil, Construction and Environmental Engineering. Vol. 2, No. 4, 2017, pp. 100-111. doi: 10.11648/j.jccee.20170204.11

Abstract
In Ethiopia, the failure of roads is before the handover period. Most roads found in Ethiopia are flexible pavements. Nowadays, the failure of the surface of flexible pavement roads is common before the expected design period. For example, the Bako-Nekemte road/ has become a critical issue in our country. The most common parameters that cause stress, strain, and deflection of the roads are loads and pressures that come from vehicles. Moreover, the modulus of elasticity, Poisson’s ratio, and thickness of each layer need to be characterized. Further, the load magnitude, contact pressure (or load radius), and location are defined for each load (wheel) considered. The finite element method (FEM) is a numerical analysis technique to obtain the stress-strain and deflection of each pavement layer. The analytical method usually uses layers thickness, loads, elastic modulus, and Poisson’s ratio of the pavement materials as design parameters. The objective of this research was to study the sensitivity of the road parameters in analyzing the major causes of failure in asphalt pavement layers fatigue cracking and rutting deformation which came due to the critical tensile strains at the bottom of the asphalt layer and the critical compressive strains on the top of subgrade using the finite element method by relating the standard specification of ERA and laboratory test result. This thesis studied the analysis of stress-strain and deflection of flexible pavements using the Everstress finite element method. The Ever stress program will take into account any stress-dependent stiffness characteristics. This thesis dealt with ways to reduce deflections by varying the design configuration, such as increasing the HMA modulus, the base modulus, sub-base modulus, the subgrade modulus, and increasing the thickness of each layer. Based on the type of materials to use the value of elastic modulus and poison’s ratio are various in each layer, in layer 1 is varied from 1500 to 3500 MPa, in layer 2 is varied from 200 to 1000MPa, in layer 3 is varied from 100 to 250 MPa and in layer 4 is varied from 20MPa to 150MPa.
Abstract
In Ethiopia, the failure of roads is before the handover period. Most roads found in Ethiopia are flexible pavements. Nowadays, the failure of the surface of flexible pavement [...]
- 1
- read
Efficiency Assessment of the Complex Process of Variable Cutting Parameters
M. Jia, Y. Su, Z. Zuo, Y. Zhan*, Y. Xu, D. Xiang, . Diamond20210

Diamond20210's personal collection (2021). 1

Abstract
Saving energy and low carbon strategy has becoming the future trend of manufacturing industry. Cutting process in mechanical manufacturing has the characteristics of heavy energy consumption and complex process. The estimation and evaluation of energy consumption and energy efficiency of cutting process are hot topics in research groups and industry. The energy consumption of cutting process is determined by the load, which is related to the cutting parameters under fixed cutting systems, while the cutting parameter might change with the cutting process and material. Based on the relationship between energy consumption and cutting characteristics and the quality requirements, the cutting specific energy u(t) = P(t)/M(t), which is varies with time is proposed to evaluate the energy consumption of various cutting parameters under complex machining processes. The influence of cutting parameters on cutting energy consumption is studied by experiments and a calculation method is proposed to calculate cutting energy consumption. Finally, a case study of complex curve machining in mold manufacturing is illustrated to prove the effectiveness of the proposed evaluation method.
Abstract
Saving energy and low carbon strategy has becoming the future trend of manufacturing industry. Cutting process in mechanical manufacturing has the characteristics of heavy [...]
- 2
- read
Offshore Oil and Gas Safety: Protection against Explosions
D. Brkić, Z. Stajić

Journal of Marine Science and Engineering (2021). 9(3), article 331.

Abstract
Offshore oil and gas operations carry a high risk of explosions, which can be efficiently prevented in many cases. The two most used approaches for prevention are: (1) the “International Electrotechnical Commission System for Certification to Standards Relating to Equipment for Use in Explosive Atmospheres” (IECEx) and (2) European “Atmosphere Explosible” (ATEX) schemes. The main shortcoming for the IECEx scheme is in the fact that it does not cover nonelectrical equipment, while for the ATEX scheme, it is due to the allowed self-certification for a certain category of equipment in areas with a low probability of explosions, as well as the fact that it explicitly excludes mobile offshore drilling units from its scope. An advantage of the IECEx scheme is that it is prescribed by the US Coast Guard for protection against explosions on foreign mobile offshore drilling units, which intend to work on the US continental shelf but have never operated there before, with an additional requirement that the certificates should be obtained through a US-based Certified Body (ExCB). Therefore, to avoid bureaucratic obstacles and to be allowed to operate with minimized additional costs both in the US and the EU/EEA’s offshore jurisdictions (and very possibly worldwide), all mobile offshore drilling units should be certified preferably as required by the US Coast Guard.
Abstract
Offshore oil and gas operations carry a high risk of explosions, which can be efficiently prevented in many cases. The two most used approaches for prevention are: (1) the [...]
- 2
- read
Simplified Indirect Estimation of Pump Flow Discharge: An Example from Serbia
M. Amon, M. Radić, Z. Stajić, D. Brkić, P. Praks

Water (2021). 13(6), article 796.

Abstract
In the absence of a flowmeter or due to its inaccuracy, the flow rate at the discharge section of the pipeline following the observed pump can be roughly estimated if the pressure can be measured instead. To use the proposed procedure two main conditions should be achieved: (1) a manometer should be installed at the discharge pipeline between the pump and the flow regulation valve, and (2) the actual curve that relates pressure and flow for the observed pump unit should be known in advance. The described example is from Serbia, but it is of interest for any water pumping station with a submersible pump (installed in wells or tanks) where a limited number of adequate places for the measuring of flow are available (if any are available at all), but where the pressure at the discharge pipe of the observed pump can be measured. This simplified method can find applicability in installations in remote rural regions where limited resources are available. The results show that the calculated values of the flow obtained by the presented method deviate greatly in relation to the measured values provided by the portable ultrasonic flowmeter, up to 60% at one of the measuring points. However, in relation to the measured values provided by the permanently installed flowmeter the discrepancy is significantly lower (0.6–6.8%).
Abstract
In the absence of a flowmeter or due to its inaccuracy, the flow rate at the discharge section of the pipeline following the observed pump can be roughly estimated if the [...]
- 1
- read
Calcined Termite Hill Clay Powder: As Partial Cement Replacement in Production of C-25 Grade Concrete
G. Danbala, E. Quezon, G. Kebede

American Journal of Civil Engineering and Architecture (2019). Guyo Konso Danbala, Emer Tucay Quezon, and Getachew Kebede, “Calcined Termite Hill Clay Powder: As Partial Cement Replacement in Production of C-25 Grade Concrete.” American Journal of Civil Engineering and Architecture, vol. 7, no. 3 (2019): 128-134. doi: 10.12691/ajcea-7-3-2.

Abstract
Cement is one of the components of concrete plays a great role in the construction industry. Nevertheless, it is the most expensive and environmentally unfriendly material. Therefore, requirements for economical and more environmentally-friendly cementing materials have extended interest in other cementing materials that can be used as partial replacement of the ordinary Portland cement. This research was, conducted to examine the suitability of calcined termite hill clay powder as a cement replacing material in the production of C-25 grade concrete as a relief for this problem. The Termite Hill Clay sample was collected from Bokuluboma vicinity. The sample was calcined at the temperature of 650°C using muffle furnace, ground to the fineness of 150 µm and its chemical composition was investigated. Normal consistency and setting time of the paste having Ordinary Portland Cement and calcined termite hill clay powder from 0% to 40% replacement in 5% increment percentage investigated. The compressive strength of eight different concrete mixes with the CTHCP replacing the ordinary Portland cement prepared for 22MPa and 34MPa (i.e., for 7th and 28th days targeted compressive strength) concrete with water cement ratio of 0.5 and 360kg/m3 cement content. The properties of the mixture assessed both at the fresh and hardened state. Results revealed that calcined termite hill clay powder found pozzolanic and can partially replace cement. It has shown that up to 11.3% replacement of the ordinary Portland cement by CTHCP, the cubes achieved a target mean compressive strength of C-25 grade concrete at 28th day 34Mpa. Although, the replacement percentage greater than 11.3% of the cement of calcined termite hill clay powder in the concrete has shown a slightly lower compressive strength. Moreover, the mix indicated workable up to 25% replacement percentage and also the setting time of the paste containing calcined termite hill clay powder found faster than that of the control mix. Hence, it is concluded that CTHCP was suitable to replace cement partially with 11.3% replacement of cement by calcined termite hill clay powder as the optimal percentage of the production of C-25 grade concrete. The fresh concrete mix was workable up to 25% replacement, and the mineral admixture served as an accelerator which is suitable for cold weather concreting works.
Abstract
Cement is one of the components of concrete plays a great role in the construction industry. Nevertheless, it is the most expensive and environmentally unfriendly material. [...]
- 3
- read
Evaluation of the Performance of Brick Dust as a Filler Material for Hot Asphalt Mix Design A Case Study in Jimma Zone
F. Wagaw, P. Quezon, A. Geremew, E. Quezon

The International Journal of Engineering and Science (2018). Fisseha Wagaw, Prof. Emer Tucay Quezon, Anteneh Geremew. Evaluation of the Performance of Brick Dust as a Filler Material for Hot Asphalt Mix Design: A Case Study in Jimma Zone. The International Journal of Engineering and Science. 2018; 7(3): 64-72. DOI:10.9790/1813-0703026472

Abstract
A bituminous paving mixture is a mixture of coarse aggregate, fine aggregate, and bitumen mixed in suitable proportion to result in a strong and durable mix in order to withstand traffic load. In this paving mix, ordinary stone dust, and cement used as filler material. One of the main problems in the construction of a bituminous paving mixture is the insufficient amount of filler from crushing stone aggregate, and cement supply is low. A study has been carried out in this research to explore the use of Brick dust as filler material for the bituminous mixture. The objective is intended to evaluate the performance of brick dust as fillers in Hot Asphalt Mix design, according to the test procedure specified by ASTM. Several specimens of hot asphalt mixture were prepared according to ASTM D 1559 using an aggregate blend with brick filler and aggregate blend without brick filler. The aggregate blending made by using Job mix formula to obtain the percentage of mixed material. For aggregate blended without brick G-1 32%, G-2 23%, and G-3 45% for Aggregate blended with brick filler G-1 30%, G-2 18%, G-3 45% and G-4 7%; where G-1 Coarse Aggregate 3/4, G-2 Coarse Aggregate 3/8, G-3 Fine Aggregate, and G-4 brick filler. It concluded that the results of the Marshall test of mix design showed satisfactory when hot asphalt mixed with these brick fillers. The Specimens blended with brick filler lead produce asphalt mixture with higher Marshall stability, lower flow, a less void filled with asphalt. Hence, brick dust can replace stone dust and cement filler in the bituminous paving mix. It is recommended to use brick dust as filler material in a bituminous paving mix may save considerable investment; as well as a reliable performance of the in-service highway can be achieved.
Abstract
A bituminous paving mixture is a mixture of coarse aggregate, fine aggregate, and bitumen mixed in suitable proportion to result in a strong and durable mix in order to withstand [...]
- 1
- read
International Journal of Scientific & Engineering Research
E. Quezon

International Journal of Scientific & Engineering Research (2016). Tarekegn Kumala, Prof. Emer T. Quezon, Bogale Shiferaw. Effect of Cycle Time and Signal Phase on Average Time Delay, Congestion and Level of Service: A Case Study at Hager Astedader Signalized Intersection in Addis Ababa. International Journal of Scientific & Engineering Research. 2016; 7(10): 524-530.

Abstract
High urbanization rate and Economic development have caused many challenges to transportation systems. Among these, long time delay and high fuel consumption of vehicles at congestion places are a few to mention. Some literatures have revealed that road traffic congestions are caused by inadequate infrastructures, long signal cycle time, and poor traffic management, such as incapable roads, inefficient public transit, and high travel demand. The research study was focused on the effect of traffic congestion on average time delay at a selected signalized road intersection in Addis Ababa. The study area was the Hager Astedadrerroad traffic signal intersection. The methodology that has been employed for the study was a quantitative descriptive research design method. The data needed were; road geometry data, signal data, traffic vehicles flow data including the pedestrian data were collected on-peak hours (with 15minutes intervals) from 7:30 to 9:30 AM and from 5:30 – 7:30 PM for the four consecutive working days. Data of traffic classes were extracted manually on a separate worksheet. The volume of each vehicle category was converted to the same vehicle category using the passenger car unit (PCU) of each vehicle class. Data analysis and processing have been performed using SIDRA (Signalized and Unsignalized Intersection Design and Research Aid) intersection software in order to know the traffic flow condition at the intersection. The result of average time delay was used to know the corresponding level of service and operational performance of the intersection. The outcomes of the research work minimize average time delay by adopting different improving strategies that range from low-cost measures such as improvements to signal timing and phase numbers, to high-cost measures such as intersection reconstruction, which be applied in the study area after knowing the amount of average time delay as well as the types of the level of service within the intersection area.
Abstract
High urbanization rate and Economic development have caused many challenges to transportation systems. Among these, long time delay and high fuel consumption of vehicles at [...]
- 0
- read
Regression Analysis of Index Properties of Soil as Strength Determinant for California Bearing Ratio (CBR)
D. Abdella, T. Abebe, P. Quezon, E. Quezon

GSJ (2017). Dino Abdella, Tadesse Abebe, Prof. Emer Tucay Quezon. Regression Analysis of Index Properties of Soil as Strength Determinant for California Bearing Ratio (CBR). GSJ. 2017; 5(6): 1-12

Abstract
Investigation of the variation in different soil types and origins is an essential task for Geotechnical engineers. To overcome the effects of this change, the Geotechnical engineers, as well as other professionals, attempted to develop empirical equations unique to a region and soil type to use the soil for its intended purpose. However, these empirical equations are more reliable for the kind of soil where there is a strong relationship of parameters is developed. Hence, it is good practice to develop empirical equations that best fit the soil available in the location that we can access. On the flexible pavement, the subgrade is considered to be an ideal layer to resist wheel load, and its CBR value is considered as the strength measuring a parameter. Conducting a CBR test is an expensive and time-consuming procedure. Also, it is challenging to mold the sample at a desired in-situ density in the laboratory. Furthermore, if the available soil is of poor quality, suitable additives are mixed with soil, and the resulting strength of the soil is assessed by the CBR value which is cumbersome. To overcome such a problem, statistical approaches such as regression-based models are used. This can be used for quick and easily determined parameters. Relative to this, the study has been conducted to develop an equation to show the relationship between the index properties and CBR values precisely located along WelkiteArekit –Hosanna Road of about 121km stretches. It was carried out using thirty samples retrieved from this road and tested in a laboratory. The test result of the regression-based statistical analysis was used to develop the predetermined relationship. The relationship development was performed in the form of an equation of CBR as a function of grain size parameters, Atterberg’s liquid limits, and compaction parameters by considering the effect of individual soil properties as well as the effect of a combination of soil properties on the CBR value. Based on the results of the study for both linear and multiple linear regression analyses, it was revealed that there was a relatively fair relationship obtained by combining plasticity index, the percentage of fine content, and maximum dry density which are strength determinants of fine-grained soils. Also, the results showed that the coefficient of determination for multiple linear regression is R2=0. 731, while for single linear regression is R2=0. 682. Therefore, it is concluded that the index properties of soils are sufficiently accurate in determining the CBR values, of which it can be utilized for preliminary characterization.
Abstract
Investigation of the variation in different soil types and origins is an essential task for Geotechnical engineers. To overcome the effects of this change, the Geotechnical [...]
- 0
- read
Effects of human behaviors and geometric factors in road traffic accidents a case study along Sululta-Fiche road, Ethiopia
E. Quezon

GJS (2017). Galeta Chala, Prof. Emer Tucay Quezon, Tarekegn Kumela. Effects of human behaviors and geometric factors in road traffic accidents: a case study along Sululta-Fiche road, Ethiopia. GSJ: 2017; 5(11). 163-172

Abstract
A road traffic accident was one of the severe problems which are affecting once’s a country and leading to death, injuries and property damage resulted in economic loss. The objective of this research was to investigate the effects of Human behaviors and Geometric factors in road traffic accidents along the road of Sululta-Fiche, which is located in the Northern Shewa of the Oromia region. Both probability sampling and non-probability sampling techniques were applied. The primary data were obtained from site observation of blackspot locations. While the secondary source of data was taken from records compiled by the North Shewa traffic office and Oromia Special Zone Surrounding Finfinne traffic Office. Descriptive and analytical design methods were used in this research, while the results were displayed by graphs, tables, and pie charts. From the study, the status of road traffic accidents showed that there were dramatically increase in accidents every year within the study period. A total of 830 crashes occurred in the past three years. From this, 237 accidents happened in 2012, 279 accidents in 2013, and 314 accidents in 2015. From the traffic police report analysis, there were 41 locations in Sululta Town, and 70 locations in Sululta Wereda, where those accidents happened. Based on the results of this research using the priority value approach, there were 12 selected locations in Sululta Town and 10 chosen locations in Sululta Wereda. Therefore, it is concluded that the significant factor in traffic accidents along Sululta –Fiche road was the driver errors such as over-speeding, night driving, and driving without attention, and it was followed by the geometric factors. These road deficiencies contributory to some errors of the drivers, which were inadequate road width, insufficient sight distance, and narrow bridge at sag curves, and improper intersection design. Based on the findings of the study, some of the countermeasures recommended reducing the road traffic accidents are the provision of speed limit signs to forewarn drivers, installation of speed barkers and rumble strips, widening the lane width of pavement, maintenance of the deteriorated road, repair of the road signs and provision of adequate sight visibility for drivers.
Abstract
A road traffic accident was one of the severe problems which are affecting once’s a country and leading to death, injuries and property damage resulted in economic loss. [...]
- 0
- read
Effectiveness of using Geosynthetic Material for Improvement of Road Construction and Performance Case Study on Adis Ababa
M. Gurara, K. Getachew, E. Quezon, E. Quezon

International Journal of Engineering Research & Technology (2017). Mikiyas. J Gurara, Kabtamu Getachew, Emer Quezon. Effectiveness of using Geosynthetic Material for Improvement of Road Construction and Performance - Case Study on Adis Ababa. International Journal of Engineering Research & Technology. 2017; (6)2: 640-642

Abstract
The research is concerned with the Effectiveness of geosynthetic to improve the performance of the pavement. It has been done on road construction projects under Addis Ababa Road Authority. Results indicated that the performance of bottom pavements structures constructed on poor subgrades with CBR values not greater than 4% was significantly enhanced using geosynthetic material. An increase of bearing capacity for subgrade and can improve the subgrade quality without excavating and filling of select material only by using geomembrane. Geosynthetic reinforcement in pavement design and construction should be widespread. Geosynthetic reinforcements are incorporated into permanent, paved roads either as base (or subbase) reinforcement - in flexible pavements to aid in the support of vehicular loads over the life of the pavement; or as subgrade restraint for construction of flexible or rigid roadways over weak subgrade conditions to aid in support of equipment loads on the unpaved base, or subbase, course during construction. Clearly, both base reinforcement and subgrade restraint with geosynthetics are proven techniques for use in pavement design and construction effectively.
Abstract
The research is concerned with the Effectiveness of geosynthetic to improve the performance of the pavement. It has been done on road construction projects under Addis Ababa [...]
- 0
- read
Cost and Benefit Analysis of Rigid and Flexible Pavement: A Case Study at Chancho –Derba-Becho Road Project
E. Quezon

International Journal of Scientific & Engineering Research (2016). Yonas Ketema, Prof. Emer T. Quezon, Getachew Kebede. Cost and Benefit Analysis of Rigid and Flexible Pavement: A Case Study at Chancho –Derba-Becho Road Project. International Journal of Scientific & Engineering Research. 2016; 7(10): 181-188

Abstract
Road construction projects have been implemented all over Ethiopia as part of the national development plan. Roads are one of the country’s basic infrastructural facilities where high amounts of the budget allocated every fiscal year planning period. Since the cost comprises of a large portion of government investment, a careful evaluation of the alternatives is of utmost importance to make the right choice for a particular project. In the history of Ethiopia's road development program, almost all of the road pavements are flexible, and it demands high foreign currency for asphalt material importing from abroad. In addition, flexible pavement needs to be maintained and rehabilitated within a few years after its initial construction. In view of the emerging cement factories and the availability of cement in Ethiopia, it is practical to consider rigid pavement as one of the alternatives. Relative to this, the research project was conducted with the main objective of identifying the cost and benefit of rigid and flexible pavements at the Chancho-Derba-Becho road project, North Showa Zone in Oromia. The research work had been focused on the specific objectives to determine and compare the life cycle costs of rigid and flexible pavements and to investigate all other qualitative merits of rigid and flexible pavement. To achieve these objectives, a review of related literature, design and specifications, observations and investigations of the actual pavement construction projects, evaluation of life cycle costs, the future value of money and present worth calculation were undertaken with an analysis period of 40 years. While the data considered was gathered through investigation at the actual rigid and flexible pavement projects, an examination of specifications, drawings and pavement design, Ethiopian Road Authority manuals, rehabilitation, and maintenance strategy. In this regard, the cost parameters investigated are initial construction cost, maintenance cost, rehabilitation cost, user’s cost, and salvage value, in addition to other qualitative and quantitative data. Based on the results of the research project, it revealed that the initial cost of rigid pavement was almost twice of the flexible pavement, but in the long run, the cost of flexible pavement per kilometer was found out to have 7.9 Million ETB more than the rigid pavement because of the incurring costs of maintenance through its design life. Therefore, it is suggested that Portland Cement Concrete Pavement (PCCP) shall be used in pavement construction to cater to local material requirements.
Abstract
Road construction projects have been implemented all over Ethiopia as part of the national development plan. Roads are one of the country’s basic infrastructural facilities [...]
- 0
- read
Comparative Study of Labor-Based Method and Machine-Based Method at Rural Road Construction and Maintenance in Jimma Zone, Oromia Region, Ethiopia
A. Dunffa, E. Quezon, E. Agon

American Journal of Civil Engineering (2017). Anwar Kassahun Dunffa, Emer Tucay Quezon, Elmer De Castro Agon. Comparative Study of Labor-Based Method and Machine-Based Method at Rural Road Construction and Maintenance in Jimma Zone, Oromia Region, Ethiopia. American Journal of Civil Engineering. Vol. 5, No. 4, 2017, pp. 212-219. doi: 10.11648/j.ajce.20170504.13

Abstract
The construction of physical infrastructure is contributory to poverty alleviation in every country. The road is one of the major components which is very significant in Ethiopia due to the lack of highly developed infrastructure, particularly in the rural road sector. Selection of the appropriate road construction and maintenance methods such as a labor-based method or machine-based method is necessary to be financially viable. Most rural road construction projects usually suffer from the problem of selecting the appropriate technology for the proposed project and implementation of the chosen technology. In Ethiopia, Universal Rural Road Access Program (URRAP) is an employment-intensive program for employment purposes as a first step, out of poverty using labor-based road construction. However, there are problems, constraints, and challenges encountered when implementing labor-based methods leading to enforce the mobilization of machines. The aim of this research study is to carry out the comparative analysis of ongoing and finished rural road projects in Jimma Zone and to establish qualitative and quantitative merit of various technologies used in rural road construction and maintenance. Results showed from financial cost comparison, the labor-based method was 47% cheaper than a machine-based method to carry out the same tasks. On the other hand, the share of unskilled average direct labor cost comprised an astounding 65% more of the total direct cost of a labor-based method. Also, the labor-based method directly employed 300 unskilled laborers per workday for a month to complete a 1.0 km section of a rural road.
Abstract
The construction of physical infrastructure is contributory to poverty alleviation in every country. The road is one of the major components which is very significant in Ethiopia [...]
- 0
- read
Combined Effects of Molasses-Lime Treatment on Poor Quality Natural Gravel Materials Used for Sub-Base and Base Course Construction
Y. Mamuye, E. Quezon, A. Geremew

GSJ (Select a year). Yibas Mamuye, Emer Tucay Quezon, Anteneh Geremew. Combined Effects of Molasses-Lime Treatment on Poor Quality Natural Gravel Materials Used for Sub-Base and Base Course Construction. GSJ: V2018; 6(7): 621-633

Abstract
At the time of industrial development, the production of large amounts of wastage required proper disposal. To reduce the disposal problem, utilization of wastage in construction works is a very important aspect. The reuse of industrial waste had gained high momentum for achieving sustainable waste management locally and globally. Keeping this in mind, research carried out to assess the utilization of molasses, lime, and molasses-lime for improving natural gravel for sub-base and base course construction from the Jimma Jiren quarry site. The natural gravel from this quarry site is known to have poor quality since it did not satisfy the ERA standard specification for sub-base and base course construction. To use Jiren natural gravel as a sub-base or base course construction material, it is necessary to improve its engineering properties. The laboratory test matrix in each test plan included variations in additive type, additive content, and curing period. It was established 8% cane molasses by weight of dry soil as the maximum for effective stabilization of natural gravel for sub-base construction, but not suitable for a base course layer. Results indicated that the natural gravel specimens from the Jiren quarry site stabilized with lime provided higher CBR values than molasses used alone. On the other hand, replacing 50% of lime content with molasses indicated better results than the performance of either individual additive. Therefore, the natural gravel blended with the molasses-lime combination is much better when it used for sub-base and base course construction to treat such an inferior quality of natural aggregates.
Abstract
At the time of industrial development, the production of large amounts of wastage required proper disposal. To reduce the disposal problem, utilization of wastage in construction [...]
- 0
- read
Geotechnical Conditions and Stability Analysis of Landslide Prone Area A Case Study in Bonga Town, South-Western Ethiopia
E. Quezon

International Journal of Scientific & Engineering Research (2017). Damtew Tsige, Prof. Emer T. Quezon, Dr. Kifle Woldearegay. Geotechnical Conditions and Stability Analysis of Landslide Prone Area: A Case Study in Bonga Town, South-Western Ethiopia. International Journal of Scientific & Engineering Research. 2017; 8(4): 239-248

Abstract
There were several slides and associated ground subsidence which brought significant impact on the cracking of walls and floor of several private and governmental buildings in Bonga Town. The principal and secondary roads were also affected by subsidence with vertical displacement up to 1m which hampered the traffic in the town. Water pipelines along the road were disturbed by the sliding which was later repaired. Cracking of the walls and floor of more than 120 private residences and more than 10 government buildings were recorded. The main highway that connects Bonga-TepiMasha via Alamo and Gatiba has been disrupted at four locations. This resulted in hampering traffic for several days. This research aimed to evaluate the cause and failure mechanism as well as the stability condition of the landslides. The study involved the investigation of the Geotechnical parameters of soil and the terrain characteristics to be used for the stability analysis of the slope, including distribution and characteristics of soils, the groundwater table, and the depth and geometry of the failures. The Slope stability analysis is supplemented by using Geo-studio 2004 software. Soil samples were collected and were tested for grain size analysis, distribution analysis (sieve & hydrometer) plastic limit, liquid limit, plasticity index, water content, unit weight of soil, specific gravity, and shear strength parameters following the ASTM procedures. Based on the findings, the landslides were triggered by heavy rainfall. Therefore, the main factors controlling the stability of the slope are soil type and characteristics, slope angle, water (surface and groundwater), and slope steepness. The design of the retaining wall is recommended to mitigate the impact of landslides in the study area.
Abstract
There were several slides and associated ground subsidence which brought significant impact on the cracking of walls and floor of several private and governmental buildings [...]
- 0
- read
Investigation of the Damping Ratio and Shear Modulus of Soil along Light Rail Transit Route in Megenagna-Hayat Road
E. Quezon

International Journal of Scientific & Engineering Research (2017). Abenezer Takele, Emer T. Quezon, Tadesse Abebe. Investigation of the Damping Ratio and Shear Modulus of Soil along Light Rail Transit Route in Megenagna-Hayat Road. International Journal of Scientific & Engineering Research. 2017; 8 (1):111-122

Abstract
Soil vibration may also cause by the movement of heavy vehicles and train, which affect adjacent structures along the route. The response of soils for incoming vibration from the train is measured by the dynamic properties of soils such as shear modulus and damping ratio. These soil parameters are very significant to study the ground motion as well as the site response of soil deposits under cyclic loading and soil-structure interaction. About this, the research study has been sought to give some information about the condition, specifically the shear modulus and damping ratio of the soil along the Addis Ababa light rail transit route from Megenagna to Hayat road. Also, this study was dealt with soil distribution along the route, which was conducted to determine soil properties. The sand replacement method was also utilized to investigate the controlling factors for the dynamic soil properties. From related literature, it was known that the dynamic soil properties are influenced by soil types and location of the soil profile. Five Test pits, designated A, B, C, D & E, were considered randomly for soil sampling to determine the properties of soil. The soil samples were taken from a depth of 2.5m below the surface of the roadside. Aside from the soil laboratory testing for common tests, another test was performed using a simple cyclic shear test for the representative samples. The values of normalized shear modulus and damping ratio are compared with already known curves from literature. The research study found out that the values of normalized shear modulus G/Gmax plotted against the shear strain, showed scattered points when it was compared with the curves of Seed and Idriss. Likewise, for the wet clay from the study area, all measured points were close to the known curve of the plot. On the other hand, the values of damping ratio were also compared to curves provided by Seed and Idriss. The Soil properties from all test pits had almost closer values as seen in the plot. However, the sample collected from test Pit D has the highest value of G/Gmax only, while samples from test Pit C has high values both for Gmax and Damping ratio. Test Pit C showed values of Dry density and liquid limit higher than the other test pits. Based on the comparison of the sand; for the strain less than 1%, all measured points of G/Gmax lying close to the extended curve given by Seed and Idriss. The strain of 0.2%, measured points of G/Gmax is closed to curve plot within the boundary, and it was below the curve given by Seed and Idriss. Likewise, for saturated clay, the measured points are close to the curve of Seed and Idriss. The values of damping ratio obtained for strain less than 1% obtained within the range of either the curves given by Seed and Idriss considers sand or clay soils. For strain, less than and equal to 1% the points lie within a range of clay. Therefore, the controlling factor that affects dynamic soil properties from a comparison of shear modulus and damping ratio curves based on consolidation pressure is the confining pressure, which shows high influence on the values of shear modulus and damping ratio based on the test results of the shear strain amplitudes. Index Terms—Damping Ratio, Dynamic Properties, Index property, Normalized Shear Modulus, Sand Replacement, Shear Modulus, Shear Strain Amplitude, Simple Cyclic Shear Test.
Abstract
Soil vibration may also cause by the movement of heavy vehicles and train, which affect adjacent structures along the route. The response of soils for incoming vibration from [...]
- 1
- read
Causes of Delays During Construction Phase of Road Projects due to the Failures of Contractor, Consultant, and Employer in Addis Ababa City Road Authority
E. Quezon

International Journal of Scientific & Engineering Research (2017). Amare, Y.; Quezon, ET.; Busier, M. Causes of Delays During Construction Phase of Road Projects due to the Failures of Contractor, Consultant, and Employer in Addis Ababa City Road Authority. International Journal of Scientific & Engineering Research. 2017; 8(3): 15-25.

Abstract
Construction delay considered as one of the most recurring problems in the implementation of construction projects. It is widely known to have an adverse impact on project success in terms of time, quality, and cost. The effect of construction delay is not only confined to the construction industry but also its influence on the overall economy of a country like Ethiopia. This research was conducted to assess the causes of excessive delays in the completion of road projects during the construction phase due to the failures of Employer, Consultant, and Contractor in Addis Ababa City Road Authority projects. Spearman rank correlation coefficient from the Relative Importance Index (RII) analysis was used to test the agreement between different groups of respondents who participated in the questionnaire survey and to rank the three construction parties according to their responsibility area and importance as perceived by the respondents which factor causing the delay of road projects. This research identified sixty-five (65) causes of delay. There were fifty-one (51) valid questionnaires received back from Contractors, Consultants, and Employer (AACRA). Based on the results, the contractors have the highest percentage of responsibility area that causes the delay of about 40%. While the second was on the part of the Employer, which comprised 26.15%, and the consultant which placed third of 23.08%. On the other hand, there were 10.77% of the respondents attested that the Shared groups (3-parties) have the responsibility area which causes project delay during the implementation phase. About this, the research study identified and ranked the top ten factors causing delays of construction projects in Addis Ababa City Road Authority. Poor ﬁnancial control of the project ranked 1st with Relative Importance Index (RII) of 0.905. Difficulties in financing projects by contractors ranked 2nd with the Relative Importance Index (RII) of 0.854. Type of project bidding and award (lowest bidder) ranked 3rd with an RII of 0.850. Poor site management and supervision of contractors ranked 4th with an RII of 0.839. Selecting inappropriate contractors ranked 5th with an RII of 0.823. Lack of high-technology mechanical equipment ranked 6th with an RII of 0.819. Inaccurate initial project scope estimate and Ineffective project scheduling ranked 7th and 8th with an RII of 0.803. Weak control of the project progress ranked 9th with an RII of 0.788. And the Contractor’s staff is not adequately trained in professional construction management techniques ranked 10th with an RII of 0.784. Therefore, it concluded that the main party Contractor did not perform properly his duties and obligations leading to the main contributory factor causing the failure of the project.
Abstract
Construction delay considered as one of the most recurring problems in the implementation of construction projects. It is widely known to have an adverse impact on project [...]
- 0
- read
A Study on the Effect of Time Duration by Vibrating or Tamping Fresh Concrete on the Compressive Strength of C-25 Concrete
A. Achalu, E. Quezon, G. Macarubbo

Effect of Time Duration by Vibrating or Tamping Fresh Concrete on the Compressive Strength (Select a year). Afework Achalu, Emer Tucay Quezon, Gerry Bangayan Macarubbo. A Study on the Effect of Time Duration by Vibrating or Tamping Fresh Concrete on the Compressive Strength of C-25 Concrete. American Journal of Construction and Building Materials. Vol. 2, No. 2, 2017, pp. 40-47. doi: 10.11648/j.ajcbm.20170202.11

Abstract
Many years ago, consolidation accomplished by laborers wielding a variety of spades, tampers, and similar tools. In Ethiopia, there is no enough knowledge by the laborers, while the site engineers ignored the correct application due to some constraints. The agreement of acceptable vibration level in fresh concrete to which tamping time duration is optimum is an issue to obtain the desired specific compressive strength. This research study aims to evaluate the maximum time length and techniques used to Vibrating or Tamping fresh concrete and its effect on the compressive strength of C-25 Concrete through internal Rodding. Therefore, to study the effects of the period by vibrating or tamping, the needle-type internal vibrator was used to test new concrete cubes based on ACI309R-96. The quality of concrete is important to make physical tests on materials used before any actual experiments to perform. The fresh concrete samples tested by using internal vibration methods and techniques of compaction with different time duration of 1sec, 2sec, 3sec, 4sec, and 5sec. These tests performed by using steel rod in fresh concrete to know the Optimum time of compaction by internal vibration on the cube test of new concrete. Laboratory experiments had been conducted to check the methods and techniques applied to construction sites, as well as site observation on the effects of the vibration time of fresh concrete on the compressive strength and flexural strength of C-25 Concrete with varying vibration time duration and steel rod tamping. Based on the results of cube tests of compressive strength of Concrete Samples have been checked on the 3rd day, 7th day, and 28th day. This research study found out that the Optimum period of compaction was 3sec to produce an acceptable right quality of the concrete mix. Likewise, the flexural strength of the Concrete Sample, when it tested in the laboratory on the 7th day and 28th day, results revealed that the flexural strength improved when vibration time increased, but the amount of water has to be reduced to obtain the desired strength. Finally, Tamping is needed to place fresh concrete in the mold correctly. Furthermore, new concrete produced should match the properties and avoid improper tamping or vibration applications. This situation must not cause differences in the compressive strength of concrete.
Abstract
Many years ago, consolidation accomplished by laborers wielding a variety of spades, tampers, and similar tools. In Ethiopia, there is no enough knowledge by the laborers, [...]
- 0
- read
Contract claim Analysis on Building Construction Project in Addis Ababa: A case study at Yeka Sub City
E. Quezon

International Journal of Scientific & Engineering Research (2016). Zenebe, E.; Quezon, ET.; Mosisa, A. Contract claim Analysis on Building Construction Project in Addis Ababa: A case study at Yeka Sub City. International Journal of Scientific & Engineering Research. 2016; 7(7): 1154-1160

Abstract
Contract claims are an inevitable fact of life in the construction industry. Contract claims are one of the problems that construction projects may face. Claims might cause many effects such as payment delay, an extension of time, and work suspension. The main aim of this research was an analysis of contract claims on building construction projects in Addis Ababa, Yeka Sub City by investigating the causes of contract claims and determining the factors to be considered in the selection of appropriate methodology for contract claims analysis. Literature and survey-based studies were made to provide an overview of the different contribution effects of contract claim analysis on building construction projects at the selected site. Interviews were also conducted with the three parties which have direct involvement in construction activities to get concrete ideas on the subject matter. A total of 75 questionnaires was distributed to respondents in Addis Ababa, Yeka Sub City building construction projects. Out of these 12 for client representatives, 35 for contractors representatives, 20 for consultant representatives, and 8 for other companies in order to identify claim causes that occurred in the study area. The quantitative data were analyzed using Excel and SPSS version 20 software, while qualitative data were analyzed by content analysis. The main factors of contact claim cause in the study area were identified. These factors further categorized into four groups; client-related, contract document-related, contractor related and consultant-related factors. While types of delays, knowing the strength and weakness of the project (project complexity) and conditions of the contract were identified as the factor to be considered in the selection of appropriate methodology for contract claim analysis. It is recommended that the three main contracting parties (Owner, Contractor, and Consultant) have responsibility for the effective accomplishment of the project. So each party would take part in their homework on time by developing technical knowledge about the effects of Contract claims on the progressive construction projects. Index Terms— causes, claim analysis, Contract claims, construction project, and risks.
Abstract
Contract claims are an inevitable fact of life in the construction industry. Contract claims are one of the problems that construction projects may face. Claims might cause [...]
- 0
- read
Evaluation of Health and Safety Practice in Building Construction: A Case Study in Addis Ababa
E. Quezon

International Journal of Scientific & Engineering Research (2016). Fekele, L.; Quezon, ET.; Macarubbo, YC. Evaluation of Health and Safety Practice in Building Construction: A Case Study in Addis Ababa. International Journal of Scientific & Engineering Research. 2016; 7(10): 122-131

Abstract
The construction industry has been seen as one of the hazardous industries. This is because the industry has a poor health and safety performance record compared to other industries all over the world. Labor law in every nation provides that it is the duty of an employer to ensure that every worker employed works under satisfactory, safe and healthy conditions. The research is focused on evaluating the health and safety issues in building construction projects in Addis Ababa, Ethiopia. The roles of the Client, Consultant and Contractor in health and safety consideration during project implementation were examined. It also determined the factors affecting health and safety performance of laborers. Whether, safety is used as criteria in selecting a contractor during the prequalification phase of public bidding was also identified. Building Construction Companies were involved in the sampling. All the information that could help in attaining the study objectives were collected, reviewed and formalized from the literature review. Data were collected and evaluated and the Relative Importance Index for each factor affecting health and safety was determined. The study determined the practical situation of building construction projects in consideration of health and safety, which was evaluated as low to medium. Relative Importance (RII) index of 12 factors affecting health and safety performance of labors were calculated. The top three factors identified were: (1) Non availability of a clear company Health and safety policy (0.85); (2) Inadequate enforcement of the existing building rules and regulations (0.82); and (3) Safety awareness of the company's top management (0.80). To improve the safety and health performance of building construction projects, the following major recommendations were proposed: each building, construction projects should have their own safety and health policy; include safety as a pay item in contract document; and allocate budget and time frames for health and safety in the contract document. Index Terms— Addis Ababa, Company Policies, Hazards, Health and Safety Practices, Rules and Regulations .
Abstract
The construction industry has been seen as one of the hazardous industries. This is because the industry has a poor health and safety performance record compared to other [...]
- 0
- read
Assessment on the effects of weather change on road construction planning at some selected projects in Addis Ababa and Oromia
E. Quezon

International Journal of Scientific & Engineering Research (2016). Befkadu, A.; Quezon, ET.; Busier, M. Assessment on the effects of weather change on road construction planning at some selected projects in Addis Ababa and Oromia. International Journal of Scientific & Engineering Research. 2016; 7(12): 134-145.

Abstract
Road construction in Ethiopia is in high demand to meet its short and long-term goals in infrastructure development programs. The road network is constructed radiating outwards North-South and East-West direction of the Capital City. Some projects are experiencing weather malady due to changing weather conditions in project locations. The study area of this research was focused on some selected projects in Addis Ababa and Oromia. The research had been addressed the effects of changes in weather in the three project planning phases, namely conception planning, design planning, and construction planning. A quantitative research design was used in this study. A literature review was undertaken and the result of which showed that the effects of weather changes are not extensively researched in the construction industry and there are limited references on the effects of weather change within the road construction industry. Sample data were obtained from sixty-three (63) construction practitioners through questionnaires, which were mailed to prospective respondents to assess their opinions. Questions were subdivided into three relevant areas covering demographics, experience with weather change factors, and the effects of weather change on road construction project planning. The study found that the effects of weather change on road construction project planning are similar to those found in other construction industries. The effects of weather change factors on road project planning vary with planning phases. These weather change factors as well as their effects were listed and ranked based on their frequency of occurrence and mean score. An independent sample t-test was used to compare respondents’ opinions on weather change effects between the two sample groups taken in Addis Ababa and part of Oromia. A partial correlation was also employed to identify the relationship between weather change factors and their effect on road construction project planning using SPSS-20. Delay of site preparation activities as a result of weather change on road project conception planning phase has a mean score of 4.30 Addis Ababa and 4.00 in for Oromia, while the influence of weather change on the selection and specification of construction materials for road project design planning phase has the highest mean score of 4.23 for Addis Ababa and 4.03 for Oromia. The effect of weather change on road workers' physiological comfort and quality of work for the road project construction planning phase has the highest mean score of 4.80 and 4.58 in Addis Ababa and Oromia, respectively. The study concludes with recommendations for potential research and with useful information and considerations on the effects of weather change on road construction project planning for construction professionals and project managers.
Abstract
Road construction in Ethiopia is in high demand to meet its short and long-term goals in infrastructure development programs. The road network is constructed radiating outwards [...]
- 0
- read
Data-sets of the model basin experiments carried out within the NICESHIP project
J. Gutiérrez Romero, S. Ruiz-Capel, J. Esteve-Perez, B. Zamora-Parra, J. Luna-Abad

NICESHIP PROJECT (2021). 1

Abstract
Photogrammetry data-sets corresponding to the model basin experiments carried out within the NICESHIP project.
- 4
- read
Immunisation for HIV-Infected People
D. Cebo

(2017). Current Topics in Virology

Abstract
Vaccines are critical components for protecting HIV-infected adults from an increasing number of preventable diseases. However, missed opportunities for vaccination among HIV-infected persons persist, likely due to concerns regarding the safety and efficacy of vaccines, as well as the changing nature of vaccine guidelines. In addition, the optimal timing of vaccination among HIV-infected adults with regard to HIV stage and receipt of antiretroviral therapy remain an important question. This article provides a review of the current recommendations regarding vaccines among HIV-infected adults and a comprehensive summary of the evidence-based literature of the benefits and risks of vaccines among this vulnerable population.
Abstract
Vaccines are critical components for protecting HIV-infected adults from an increasing number of preventable diseases. However, missed opportunities for vaccination among [...]
- 0
- read
Discordant CD4 cell count and viral load responses to antiretroviral therapy
D. Cebo

(2017). Current Trends in Immunology

Abstract
Most patients starting highly active antiretroviral therapy (HAART) show suppression of HIV replication below the level of detection (currently <50 HIV-1 RNA copies/mL in most assays), and experience a gradual rise in CD4 lymphocyte count, which may continue for several years. The CD4 count response is generally related to the degree of viral load suppression, and this typical pattern of CD4 and viral load response is associated with a marked improvement in prognosis. In some patients, however, there is discordance in the response. Either there is suppression of viral load but poor recovery of immune function, characterized by little or no CD4 cell count increase or, conversely, an improvement in CD4 cell count with incomplete or delayed viral load suppression. Little is known about the pathogenesis of discordant responses, which seems to depend on the interaction of a multitude of viral, host and treatment-related factors. Available evidence indicates that discordant responses are associated with an intermediate risk of death or clinical progression. At present, recommendations for the clinical management of patients with discordant responses to antiretroviral therapy are largely based on observational, uncontrolled data. The development of standardized and universally accepted definitions of discordant responses is necessary to allow meaningful comparisons between studies to be made, as well as to help in the design of trials of possible therapeutic interventions.
Abstract
Most patients starting highly active antiretroviral therapy (HAART) show suppression of HIV replication below the level of detection (currently <50 HIV-1 RNA copies/mL [...]
- 0
- read
PARAMETRIC MODELLING ON THE RELATIONSHIPS BETWEEN ATTERBERG LIMITS AND COMPACTION CHARACTERISTICS OF FINE-GRAINED SOILS
P. Quezon, E. Quezon

International Journal of Advanced Research in Engineering and Applied Sciences (2019). Worku Firomsa, Prof. Emer Tucay Quezon; (2019). “PARAMETRIC MODELLING ON THE RELATIONSHIPS BETWEEN ATTERBERG LIMITS AND COMPACTION CHARACTERISTICS OF FINE-GRAINED SOILS.” International Journal of Advanced Research in Engineering and Applied Sciences, 8(7), 1-20. https:// www.garph.co.uk

Abstract
Compaction is one of the essential engineering techniques, performed to assure the stability of soils based on specified strength. However, in most construction projects such as large-scale projects, obtaining the desired compaction characteristics, namely optimum moisture content (OMC) and maximum dry density (MDD), becomes time-consuming. In this case, predicting the compaction characteristics from the Atterberg's limits that involve a more straightforward and quicker method of testing becomes a vital task. This study focused on obtaining valid correlations between Atterberg's limit parameters and compaction characteristics of fine-grained soils. A series of laboratory tests for 50 samples conducted for the investigation. Statistical relationships of all the parameters were analyzed. The laboratory test results indicated that both OMC and MDD showed a strong correlation with LL and PL together using multiple linear regressions than with single parameters from single linear regression. Therefore, the study concluded that during the prediction of OMC and MDD from Atterberg's limits, the combined parameters of Atterberg's limits should be used rather than single parameters. It is recommended, the result of this research could be applied in different civil engineering practices, directly related to the parameters to be known.
Abstract
Compaction is one of the essential engineering techniques, performed to assure the stability of soils based on specified strength. However, in most construction projects such [...]
- 3
- read
Evaluation on the Performance of Lowest Responsive Bid Contract and the Quality of Materials Used on Governmental Building Projects in Jimma Town
E. Quezon

International Journal of Scientific & Engineering Research (2016). BinyamLetarge, Emer T. Quezon, Yolente C. Macarubbo. Evaluation on the Performance of Lowest Responsive Bid Contract and the Quality of Materials Used on Governmental Building Projects in Jimma Town. International Journal of Scientific & Engineering Research. 2016; 7(12): 60-73

Abstract
Construction industry participants have started recognizing that accepting the least bid price does not guarantee maximum value. Continuous problems of inferior quality of constructing facilities, high incidence of claims and litigation, frequent cost overruns, and use of poor quality of materials have become the main features of Ethiopian’s public construction work contracts. This research was undertaken to evaluate the performance of publicly owned construction projects awarded on a lowest bid awarding system, to determine the effect of advance payment on the contractor’s performance, compare lowest and average bid systems, and check the quality of local construction materials use by lowest price won contracts. A literature review was carried out to identify different practices and floated questionnaire survey and laboratory test was conducted for selected and expected materials. Two alternative bid evaluation methods were discussed and suggested the better one from a performance point of view. The questionnaire was distributed to contractors, clients, consultants, and other related professionals. Additionally, interviews were conducted with them. A total of 88 questionnaires were distributed, including laboratory test results for selected and expected materials found in Jimma town. The data were collected and 80 valid questionnaires were analyzed by using SPSS-20, Excel, and laboratory test results requirement. The study identified that the causes of the poor performance of contractors were won projects with the lowest price; therefore more of the respondents do not like the lowest bidding method. Advance payment is the main solution to support the financial problems of contractors, but instead of paying in cash, purchasing the necessary construction materials is seen as an alternative to protect contractors from using the money for other uses. Obviously known that in the current bid awarding method of Ethiopia, most of the time, in order to be the winner of the bid price should be lesser. This leads the bidder not to get adequate profit, this pushes them to use poor quality of local construction materials. Finally, this work provides valuable information to the Ethiopian government, clients, consultants and contractors, and other stakeholders who desire to improve bidding methods; usage of advance payment; ways to improve performance of contractors, and to protect the project from contractors using poor quality of construction materials.
Abstract
Construction industry participants have started recognizing that accepting the least bid price does not guarantee maximum value. Continuous problems of inferior quality of [...]
- 1
- read
Nonlinear Dynamics of Structures: Projects, Tests, Damping, Damages and Their Diagnostics
V. Kulyabko

WCCM-ECCOMAS2020.

Abstract
The paper considers the features of nonlinear modeling of problems of dynamics of structures at different stages: design, testing and operation. Recommendations are given for nonlinear dynamic calculations and dynamic diagnostics (certification, identification of structural defects, etc.).
Abstract
The paper considers the features of nonlinear modeling of problems of dynamics of structures at different stages: design, testing and operation. Recommendations are given [...]
- 28
- read
Multiphysics Modelling Of A Hybrid Rocket Engine
A. Ferrero, F. Masseni, L. Muscarà, D. Pastrone

WCCM-ECCOMAS2020.
- 16
- read
Multidisciplinary Modelling, Analysis and Optimisation for Aircraft and System Level Design and Validation
W. J. Vankan, W. Lammen, E. Baalbergen

WCCM-ECCOMAS2020.

Abstract
Against the background of the big environmental and societal challenges as formulated for example in Flightpath 2050, current developments in aircraft design are aiming at further emission reduction through integrated, unconventional propulsion, systems and airframe innovations. This requires the further integration of methods for multidisciplinary modelling, analysis and optimization for aircraft design, but also for propulsion and system level designs. Moreover, experimental validation of the methods and physical testing of critical unconventional propulsion and system designs are prerequisites for industrially relevant development processes. This paper presents some key technologies for computationally efficient collaborative MDO (multidisciplinary design and optimization) frameworks for multidisciplinary design and validation of advanced aeronautic products like aircraft and propulsion systems.
Abstract
Against the background of the big environmental and societal challenges as formulated for example in Flightpath 2050, current developments in aircraft design are aiming at [...]
- 7
- read
A Bi-Fidelity Method for Multiscale Kinetic Models with Uncertainties
L. Liu

WCCM-ECCOMAS2020.

Abstract
In this paper, we introduce a bi-fidelity stochastic collocation (SC) method for the linear transport equation with diffusive scaling and high-dimensional random inputs characterized by random variables. For the high-fidelity linear transport model, the asymptotic-preserving Discontinuous Galerkin implicit-explicit method in the micro-macro decomposition framework is used. We discuss different choices of low-fidelity models and conduct the corresponding uniform error estimates for the bi-fidelity method. Extensive numerical tests are presented to validate our bi-fidelity SC method.
Abstract
In this paper, we introduce a bi-fidelity stochastic collocation (SC) method for the linear transport equation with diffusive scaling and high-dimensional random inputs characterized [...]
- 10
- read
Numerical Analysis of a Mechanical De-Icing Process by Low Frequency Oscillation of a CFRP Layer
F. Grubert, M. Nuño, K. Schröder

WCCM-ECCOMAS2020.

Abstract
De-icing of general aviation aircraft is usually realized by chemical or thermal processes. These lead to an increase in fuel consumption or usage of electrical energy. Mechanical de-icing is a way to significantly reduce this consumption. In mechanical de-icing, a surface is deformed so that the required failure mechanisms are induced in the ice, causing it to detach. In this paper, the release behavior of ice on a CFRP layer is investigated. The CFRP layer consists of two plies, each 0.3 mm thick. A numerical calculation is performed to determine the natural frequency and the required amplitudes of the vibration. In addition, the relationship between the various failure mechanisms of ice and the ice layer thickness, as well as the control values are determined. Tests are also being conducted to demonstrate the feasibility of mechanical de-icing. For this purpose, the surface of a CFRP layer is iced with water in a climate chamber at -20°C. A modal shaker is connected to the CFRP and generates the required displacements at desired frequencies to observe the detachment of the ice.
Abstract
De-icing of general aviation aircraft is usually realized by chemical or thermal processes. These lead to an increase in fuel consumption or usage of electrical energy. Mechanical [...]
- 9
- read
A Simple Formulation for Visco-Hyperelastic Behavior for Soft Materials Suitable for Different Loading Types
O. Garcia-Bañales, L. Saucedo-Mora, J. Benítez, F. Montáns

WCCM-ECCOMAS2020.

Abstract
Soft materials, and especially soft biological tissues, have a complex highly nonlinear behavior both for quasielastic (slow) and viscous loading. In partiular, the cyclic behavior is different depending on the loading speed, number of cycles, and their magnitude. Furthermore, different soft materials and soft tissues have different particularities in their behavior. Therefore, a phenomenological proposal capable of accurately capturing all these singularities with few, easy to obtain parameters based on experimental data, is valuable. In this study, a visco-hyperelastic one-dimensional formulation to characterize different biological tissues is proposed, which has proven to be capable of capturing the response of numerous soft biological tissues (brain tissue, coronary arteries, tendons, tongue tissues, abdominal muscle, cells...) under pure and combined loading modes, including tension, compression, simple shear and the combination of the latter one with tension and compression. One of the main advantages of the proposed model is its simplicity, being that the formulation is calibrated with four simple parameters (two of them for the hyperelastic component and four for dealing with the different viscous aspects) obtained from the uniaxial loading. The formulation, based on a combination of Maxwell and Kelvin-Voigt rheological models has proven to represent to very good accuracy the behavior of a wide range of materials under different types of loadings, including effects like preconditioning and cycle stabilization. In all these cases, under different monotonic and cyclic loading, all aspects of the viscous and elastic behavior are accurately captured. Thanks to its structure, this model incorporates strain-level dependent nonequilibrium viscoelasticity and it may be easily incorporated to 3D nonlinear finite strains formulations.
Abstract
Soft materials, and especially soft biological tissues, have a complex highly nonlinear behavior both for quasielastic (slow) and viscous loading. In partiular, the cyclic [...]
- 15
- read
Higher-Order Space-Time Continuous Galerkin Methods for the Wave Equation
M. Zank

WCCM-ECCOMAS2020.

Abstract
We consider a space-time variational formulation of the second-order wave equation, where integration by parts is also applied with respect to the time variable. Conforming tensor-product finite element discretisations with piecewise polynomials of this space-time variational formulation require a CFL condition to ensure stability. To overcome this restriction in the case of piecewise multilinear, continuous ansatz and test functions, a stabilisation is well-known, which leads to an unconditionally stable space-time finite element method. In this work, we generalise this stabilisation idea from the lowestorder case to the higher-order case, i.e. to an arbitrary polynomial degree. We give numerical examples for a one-dimensional spatial domain, where the unconditional stability and optimal convergence rates in space-time norms are illustrated.
Abstract
We consider a space-time variational formulation of the second-order wave equation, where integration by parts is also applied with respect to the time variable. Conforming [...]
- 4
- read
Probabilistic Agent Based Model For Tumoral Cells, And 3d Model For Angiogenesis
J. Benítez-Baena, L. García-Mozos, J. Ramos-Lojo, M. Benítez-Aguilar, L. Saucedo-Mora, F. Montans

WCCM-ECCOMAS2020.

Abstract
The study of tumoral cells behaviour through computational models is arising. The movement of these cells is governed by physical laws; therefore, the different forces exerted between them need to be implemented. Nevertheless, biological criteria should be also considered since other factors, as oxygen level or cellular density, are decisive in the real movement. These phenomena are captured by probabilistic models such as the Agent Based Model (ABM). Following this research line, the present paper outlines a numerical model that tries to join both criteria with the aim of reproducing the behaviour of the cells that are part of a brain tumor: Glioblastoma Multiforme (GBM). The study has been carried out by the implementation of the different force equations in a Smoothed-Particle Hydrodynamic (SPH) framework. The SPH method is a meshfree lagrangian method based on the discretization of the study domain into finite particles which carry their own information, as tumoral cells do. Cohesive, viscous and pressure forces have been taking into account. Also, the possible attraction or repulsiveness between cells is considered through the implementation of a mechanical force formulation that combines the Maxwell and KelvinVoigt viscoelastic models. In addition to this forces approach, an energetic model is proposed to consider the results provided by an ABM. It evaluates the energy consumption and the associated extra-force that the cell needs to reach the ABM position, which is considered the biologically optimal one. The model has been tested under different sets of parameters, getting the logical outcome. Successful results have also been found in the evaluation of the energy consumption and, therefore, of the extra-force, finding a formulation that joins both criteria.
Abstract
The study of tumoral cells behaviour through computational models is arising. The movement of these cells is governed by physical laws; therefore, the different forces exerted [...]
- 11
- read
Numerical Prediction Of Springback In Sheet Metal Forming In Anisotropic Elastoplastic Finite Strains
R. Díaz, K. Nguyen, F. Montáns, M. Sanz

WCCM-ECCOMAS2020.

Abstract
We present an application on springback of an anisotropic elastoplastic finite strain model. The computational procedure allows the sole use of internal elastic variables as it is fully hyperelastic. Based on the Lee-type multiplicative decomposition, which permits an additive solution in a logarithmic strain approach, it resolves the 'rate issue' by the employment of a corrector rate of elastic strains and its algorithmic implementation [1]. The associative nature and Clausius-Duhem inequality recover the formulation of Simo's isotropic strain-hardening particular case. Its application to metal and soft materials has been corroborated as it is also valid for anisotropic yield functions and for any anisotropic stored energy (linear and nonlinear in logarithmic strains) [2]. The computation of the stress-point algorithm is done using a simple backward-Euler scheme. Its validation has been performed by the implementation of the procedure into a subroutine which allows the user to compute benchmark models in the commercial program ADINA [3]. This permits to explore its usefulness in application for the numerical prediction of springback in thin sheet metal forming processes. The results are compared with both experimental and other authors' results. Moreover, many simulations were performed such as the draw-bending test, the unconstrained bending problem, the square cup or the 3DS warping benchmark.
Abstract
We present an application on springback of an anisotropic elastoplastic finite strain model. The computational procedure allows the sole use of internal elastic variables [...]
- 0
- read
High Performance Optimization Computing Platform
S. Sotiropoulos, G. Kazakis, N. Kallioras, S. Xynogalas, C. Mitropoulou, S. Chatzieleftheriou, S. Damikoukas, P. Tsakalis, N. Lagaros

WCCM-ECCOMAS2020.

Abstract
Structural optimization over the past decades matured from an academic theoretical field, to an important tool in the design procedure in various engineering disciplines. Some commercial software applications provide some suites with optimization solutions, but they are focused mostly in the aeronautics, automotive and aerospace industry. High Performance Optimization Computing Platform (HP-OCP) is a software developed by the ISAAR-NTUA and provides a holistic optimization approach for civil engineering structures. More precisely, HPOCP is a computational suite that has the ability to integrate with several structural analysis and design software and provide optimization solutions. Structural optimization is mainly divided in three groups, sizing (or parametric), shape and topology optimization. All of them are integrated in HP-OCP and the appropriate algorithms are provided in each category. Considering size and shape optimization, the parametric optimization module is developed, in which the design variables of the mathematical formulation can be the dimension of the section properties, the quality of the material, the coordinates of the nodes etc. In this module plenty of derivativebased and derivative-free algorithms are provided like the Projected Quasi-Newton, Constrained Optimization by Linear Approximation, Latin Hypercube Sampling etc. [1]. Considering the topology optimization module [2], the SIMP method is applied and the mathematical algorithms that are implemented are the Optimality Criteria and Method of Moving Asymptotes. HP-OCP was developed in C# programming language, making it a powerful suite that can be integrated with any commercial software that provide Application Programming Interface, batch analysis via XML files or any other type of data exchange format. In the current work the integration of HP-OCP with the SAP2000, ETABS and SCIA Engineering software is presented. Several examples considering parametric and topology optimization problems are examined. Remarkable cost reduction is succeeded in real-world structures, validating in this way the usefulness of HP-OCP not only in the research field but also in applied civil engineering problems.
Abstract
Structural optimization over the past decades matured from an academic theoretical field, to an important tool in the design procedure in various engineering disciplines. [...]
- 17
- read
Strain-Driven Generative Design Framework Coupled With A Mimetic Metamaterial: A Process Towards Mechanical And Shape Adaptation To Observed Structures And Functionalities
H. García-Modet, L. Saucedo-Mora, G. Gómez-Carano, M. Sanz-Gómez, F. Montans

WCCM-ECCOMAS2020.

Abstract
Topology optimization has undergone tremendous development since its introduction by Bendsøe and Kikuci in 1988, especially in recent years, due to its involvement in revolutionary generative design techniques. This paper aims to lay the foundations of a generative design methodology powered by an alternative approach to the well-known density methods. Based on finite element analysis, the objective is to develop an optimization algorithm with the Young modulus of the elements as design variables. That way, while previous studies have focused on void/solid distributions, this study searches for a distribution of different E values that could be manufactured due to progress in metamaterials and additive manufacturing. A mimetic metamaterial was also developed to be coupled with the topological optimization, but will not be included in this paper. To assess the optimization algorithm, several analyses have been carried out under different load and boundary conditions. The outcome shows correlation with our initial hypothesis: elements under higher strains increase their stiffness value, while the opposite occurs for those under minor stresses. Consequently, the results present a structure with a Young modulus distribution that optimizes the strain energy, and therefore, reduces the displacements.
Abstract
Topology optimization has undergone tremendous development since its introduction by Bendsøe and Kikuci in 1988, especially in recent years, due to its involvement in revolutionary [...]
- 14
- read
THE CHARACTERIZATION OF METAMATERIALS IN VIRTUAL TESTING: MACROSCOPIC EFFECTIVE PROPERTIES FOR STRUCTURAL APPLICATIONS
I. Gonzalez, G. Rodríguez, F. Montans, L. Mora, M. Gómez

WCCM-ECCOMAS2020.

Abstract
Metamaterials research and, moreover, their obtaining by additive manufacturing is trendy because of their potential industrial applications, since lattice structures can save weight when integrated in larger parts. For this reason, it is important to draw these new material properties in order to compare them with traditional structures obtained by classic manufacturing methods. In this work we present a numerical procedure to determine the mechanical properties of a metamaterial based on periodic cells. To date, for these materials there are no standard procedures to test and validate their structural properties, and specifically to determine their mechanical bending behaviour. We apply a virtual testing procedure using finite elements to determine the homogeneous elastic modulus, Young's modulus and shear modulus, using a set of modified bending virtual test, based on the three point bending and the four point bending test [1]. Theses physical testing procedures are common in sandwich composite materials [2], and can be apply to metamaterial testing taking into account some geometrical and material aspects to determine the proper parameters for testing. The results obtained are macroscopically independent on the number of elemental cells used in the metamaterial, and permit us to consider the intrinsic defects that may appear in the 3D additive manufacturing materials.
Abstract
Metamaterials research and, moreover, their obtaining by additive manufacturing is trendy because of their potential industrial applications, since lattice structures can [...]
- 4
- read
Geometric-Topologic Description of a Complex Road Junction Considering the Requirements of Highly Automated Driving
J. Logo, V. Poto, N. Krausz, A. Barsi

WCCM-ECCOMAS2020.

Abstract
The highly automated driving research is one of the most remarkable developments in the automotive world. Still, ordinary people can also realize the milestones facing more and more built-in vehicle assistant services brought by the new car manufacturers. The navigation system is also part of the 'success story' because the maps are currently high-resolution 3D databases, and permanently better visualization supports the drivers. The creation of these map databases relies on cutting-edge field data capture techniques, like aerial surveys. The paper presents a technology based on aerial photographs' evaluation done by human operators. The captured photographs were transformed into distortion-free orthophotos; then, skilled humans evaluated the expected database elements: road and lane segments. The map creation procedure aimed to serve the content fully compatible with the standards of the industry. This evaluation procedure resulted in a road and a lane-level graph, clearly in harmony with the traffic rules. The CAD-drawings of each operator form an excellent starting point to start new artificial intelligence developments.
Abstract
The highly automated driving research is one of the most remarkable developments in the automotive world. Still, ordinary people can also realize the milestones facing more [...]
- 5
- read
SHAPE OPTIMIZATION FOR A NOISE REDUCTION PROBLEM BY NON-INTRUSIVE PARAMETRIC REDUCED MODELING
F. Bonizzoni, D. Pradovera

WCCM-ECCOMAS2020.

Abstract
We study a PDE-constrained optimization problem, where the shape and liner material of the nacelle of an aircraft engine are optimized in order to minimize the noise radiated by the engine. More precisely, the acoustic problem is modeled by the Helmholtz equation with varying wavenumber kon an exterior domain. A model reduction strategy is employed to alleviate the cost of the design optimization: the minimal rational interpolation technique is used to construct a surrogate (w.r.t. k) for the quantity of interest at fixed shape/material parameter values, and a parametric model order reduction approach is employed to combine surrogates at different shape/material designs, resulting in a nonintrusivemethodology. Numerical experiments for shape and shape/material optimization are provided, to showcase the effectiveness of the presented methodology.
Abstract
We study a PDE-constrained optimization problem, where the shape and liner material of the nacelle of an aircraft engine are optimized in order to minimize the noise radiated [...]
- 22
- read
Grid Generation and Heat Transfer in a Complex Topology Maritime Pine Tree in Forest Fire Environments
E. Conceição, J. Gomes, M. Lúcio, J. Raposo, D. Viegas, M. Viegas

WCCM-ECCOMAS2020.

Abstract
This article presents a study on the grid generation and heat transfer in a complex topology maritime pine tree in a forest fire environment. The numerical model of the maritime pine tree thermal behaviour is developed and applied. This numerical model is based on energy balance integral equations for the maritime pine tree elements. These equations are obtained for the heat exchanges by conduction within the maritime pine tree elements, by convection between the external maritime pine tree surfaces and the environment, and by radiation between the external maritime pine tree surfaces and the fire front. The geometry of the virtual maritime pine maritime pine tree was here developed applying mesh generation, which is used to evaluate the view factors. In turn, these view factors are used to calculate the heat exchanges by radiation between the virtual maritime pine tree and the fire front. The virtual maritime pine tree geometry is constituted by 8863 cylindrical elements, which represent their trunks, branches and leaves, and it is divided into four levels. The input data considered in the numerical simulation are the mean flame temperature of 500ºC, the air temperature of 20ºC, the wind speed of 10 m/s and the air relative humidity around the pine tree of 50%. The numerical simulation results are the view factors and the evolution of the surface temperature of the maritime pine tree elements obtained in a transient regime. The values of the view factors and the surface temperature of the maritime pine tree elements decrease from the first level, near the ground, to the last level, at the top of the maritime pine tree. The stabilization time of the surface temperatures of the maritime pine tree elements decreases according to the increase of the height to which these elements are located.
Abstract
This article presents a study on the grid generation and heat transfer in a complex topology maritime pine tree in a forest fire environment. The numerical model of the maritime [...]
- 9
- read
Unstructured Level-Set Method For Saturated Liquid-Vapor Phase Change
N. Balcazar, J. Rigola, A. Oliva

WCCM-ECCOMAS2020.

Abstract
A novel conservative level-set method for saturated liquid-vapor phase change on unstructured meshes is introduced. Transport equations are discretized by the finite-volume method on collocated unstructured grids. Mass transfer promoted by thermal phase change is computed using the energy jump condition at the interface, as a function of the temperature gradient. The fractional-step projection method is used for solving the pressure-velocity coupling, convective terms are discretized by unstructured flux-limiter schemes, central difference scheme is used for discretization of diffusive terms. Verification and validation cases have been undertaken to prove the accuracy and robustness of the numerical methods, including simulation of the Stefan problem, and film boiling on a cylindrical surface. Excellent agreement between numerical solutions against analytical solution and empirical correlations from the literature is reported.
Abstract
A novel conservative level-set method for saturated liquid-vapor phase change on unstructured meshes is introduced. Transport equations are discretized by the finite-volume [...]
- 36
- read
Lagrangian Magneto-Hydrodynamics Based On Curvilinear Finite Elements
J. Nikl, M. Kuchařík, S. Weber

WCCM-ECCOMAS2020.

Abstract
The magneto-hydrodynamic model is widely used for description of magnetized fluids in plasma dynamics, microfluidics, astrophysics and many other applications. In terms of modelling, the Lagrangian formulation is favourable for the rapid expansion during lasertarget interaction for example. This is the case for inertial fusion and laboratory astrophysics applications, which are our primary interest. However, the proposed numerical method remains general and can be applied elsewhere. The conservation properties and divergence-free magnetic field are crucial aspects, which are not satisfied by the traditional numerical schemes. Here, the Lagrangian hydrodynamics using curvilinear finite elements is extended to the resistive magneto-hydrodynamics. An energy-conserving numerical scheme is formulated maintaining divergence-free magnetic field. The mixed finite element formulation provides theoretically arbitrary order of the spatial convergence and application on unstructured Lagrangian grids in multiple dimensions. An example of a physically relevant numerical simulation is presented.
Abstract
The magneto-hydrodynamic model is widely used for description of magnetized fluids in plasma dynamics, microfluidics, astrophysics and many other applications. In terms of [...]
- 6
- read
Unstructured Mesh Adaptation for Moving Bodies in Immersed Boundary Methods
T. Kozubskaya, L. Kudryavtseva, V. Tsvetkova, N. Zhdanova

WCCM-ECCOMAS2020.

Abstract
The paper presents a method for simulating the flow of a viscous compressible gas around a moving obstacle, specified by the immersed boundary method on a simply connected unstructured mesh, and the features of its application for bodies of different configuration. An advantage of the developed moving-mesh method is that it keeps the topology of an initial mesh and provides an anisotropic mesh refinement to the 'immersed' surface of a streamlined body. The position and shape of the body are determined by the interpolation grid, which stores the distance and normal to the body surface. The paper also describes a special technique to improve the efficiency of adaptation. The technique is based on the prediction of the body location in some automatically chosen time. This trick allows to find the coordinates of mesh vertices during the predicted time interval by interpolation. The integral strategy combining the immersed boundary method with the developed moving-mesh adaptation is tested on several model cases in twoand three-dimensional formulations. However, the paper does not consider the numerical results of the cases, it is mostly targeted at the meshing details.
Abstract
The paper presents a method for simulating the flow of a viscous compressible gas around a moving obstacle, specified by the immersed boundary method on a simply connected [...]
- 6
- read
Numerical Investigation On Fluid Flow of Protruding Branching Duct in a Dividing Junction of a Minichannel Heat Exchanger
M. Settati, C. Oliet, A. Oliva

WCCM-ECCOMAS2020.

Abstract
This study focuses on investigating the effect of protruding branching duct on laminar fluid field in 90◦ dividing pipe junction of mini heat exchanger. The ducts under study, main and branch side, are of minichannel scale where the corresponding diameters are 3.96 and 1.6 mm respectively. First, numerical validations are given based on the experimental measurements of Cho et al. [1] for the case of non-protruding branching duct (h=0, β ≠0) and also comparing with Poiseuille flow for nonprotruding case with no flow in the branch side (h=0, β=0) where good agreements were found in each case. Second, the degree of protrusion in pipe branching duct is characterized studying and analyzing two cases of protrusions (h = D1/4 and h = D1/2). The effects of branching flow rate ratio β on local pressure changes and the velocity profile are presented and examined in both main and branching ducts. Finally, the influence of the extent of protruded branch on velocity field and the local pressure drop are documented as expected. The results shown a strong effect of protrusion on hydrodynamic flow
Abstract
This study focuses on investigating the effect of protruding branching duct on laminar fluid field in 90◦ dividing pipe junction of mini heat exchanger. The ducts under [...]
- 4
- read
Computation of Changes in Local Quantities of Interest Due to Design Modifications
D. Materna

WCCM-ECCOMAS2020.

Abstract
In this contribution a method for the prediction of the change in a quantity of interest due to structural modifications is considered. The approach is based on a goal-oriented method by using the primal and dual problems. The method can be used in order to reduce the overall computational cost in the numerical simulation. Furthermore, the proposed method is easily to implement in existing finite element programs, because no derivatives with respect to the design variables are necessary.
Abstract
In this contribution a method for the prediction of the change in a quantity of interest due to structural modifications is considered. The approach is based on a goal-oriented [...]
- 8
- read
Development of a Hybrid Rans-Les Model Based On Temporal Filtering for Natural Convection Flows
P. Bikkanahally, R. Manceau, F. Mastrippolito

WCCM-ECCOMAS2020.

Abstract
Natural convection,Hybrid RANS-LES,Turbulence modeling Transient,buoyancy-affected flows play a major role in some industries,for instance in the design of nuclear power plants or the under-hood space of automobiles. In order to avoid expensive experiments,promote rapid design,development,it is quintessential to employ efficient,robust,accurate CFD models to simulate such flows. With the aim of being cost-effective,retaining the ability to capture the necessary transient phenomena,the industry standard of using RANS models alone does not suffice,whereas LES is as of yet not affordable for industrial configurations. In this context,a Hybrid RANS-LES model based on temporal filtering,called Hybrid Temporal Large Eddy Simulation (HTLES) is developed,thus reducing computational costs by switching to RANS in the near-wall regions while being able to predict a large part of the turbulent scales in regions where it is necessary. The present work is based on the HTLES model developed for iso-thermal configurations [1],with the hybridisation applied to the standard k SST model. The extension to natural convection flows is carried out using the DNS database of a differentially heated square cavity [2]. The test case is characterised by the three non-dimensional numbers viz. Rayleigh number(Ra = 10 ),Prandtl number(Pr = 0.71),aspect ratio(A = L /L = 1),features a stratified core with a spatially developing buoyant boundary layer. The study reveals that the HTLES model is able to predict the mean profiles of temperature,velocity satisfactorily. Comparisons with URANS show that the HTLES model predicts the kinetic energy profiles more accurately,especially when considering coarser grids. However,it is noted that the current definition of the shielding function (used to enforce RANS close to the walls based on the Kolmogorov length scale) is inadvertently modifying the energy ratio far from the wall due to fluctuations in the dissipation rate,. This is further confirmed by way of a priori tests using the DNS database. It is essential to have a shielding function that is only effective in the near-wall regions. Therefore,to ward off this issue,a shielding function based on the Elliptic blending framework is developed. For the validation of the approach,flows in a differentially heated vertical channel [3],a 4:1 cavity [4] are considered,in order to investigate the ability of the model to reproduce natural convection flows with very different stratification.
Abstract
Natural convection,Hybrid RANS-LES,Turbulence modeling Transient,buoyancy-affected flows play a major role in some industries,for instance in the design of nuclear power plants [...]
- 5
- read
Contact Problems for Interface Cracks Under Harmonic Shear Loading
V. Menshykov, O. Menshykov, I. Guz

WCCM-ECCOMAS2020.

Abstract
The linear crack between two dissimilar elastic isotropic half-spaces under normal harmonic shear loading is considered. To take the crack faces interaction into account we assumed that the contact satisfies the Signorini constraints and the Coulomb friction law. The problem is solved numerically using the iterative process the solution changes until the distribution of physical values satisfying the contact constraints is found. The numerical convergence of the method with respect to the number of the Fourier coefficients and mesh size is analysed. The effects of material properties and values of the friction coefficient on the distribution of displacements and contact forces are presented and analysed. Special attention is paid to the size of the contact zone and the results are compared with the classical model solutions obtained for the static problems with and without friction.
Abstract
The linear crack between two dissimilar elastic isotropic half-spaces under normal harmonic shear loading is considered. To take the crack faces interaction into account we [...]
- 15
- read
Micromechanics-Driven Variational Method for Diffuse-to-Localised Fracture in Quasi-Brittle Solids
C. Azua-Gonzalez, I. Mihai, T. Jefferson

WCCM-ECCOMAS2020.

Abstract
A novel multiscale variational method for modelling fracture propagation is proposed. The method employs strong discontinuity kinematics enhancement, enabling macroscopic cracks to be modelled explicitly, while minimum remeshing is ensured. In addition, the response of quadrature points in the bulk is up-scaled using a Micromechanical continua, which enables the evolution of directional micro-defects (e.g. microcracks) without venturing into prohibitive computational burden. Noticeably, the method allows the Micromechanical continua to interact with macroscopic cracks. The framework is conveniently formulated as a variational setting to provide a minimum energy solution. The new computational framework has been found to diagnose realistic failure mechanisms in quasi-brittle materials.
Abstract
A novel multiscale variational method for modelling fracture propagation is proposed. The method employs strong discontinuity kinematics enhancement, enabling macroscopic [...]
- 8
- read
PDE-BASED AND SOLUTION-DEPENDENT PARAMETERIZATION FOR ISOGEOMETRIC ANALYSIS
Z. Ali, W. Ma

WCCM-ECCOMAS2020.

Abstract
This paper presents some results for PDE-based and solution-dependent parameterization of computational domains for isogeometric analysis (IGA) using non-uniform rational B-splines (NURBS). The final parameterization is produced based on the solution of a partial differential equation (PDE) that is solved using isogeometric collocation method (IGA-C) with Dirichlet boundary condition being the input boundary of the final desired computational domain for IGA, namely the IGA-C-PDE method for domain parameterization. The theory of PDE guarantees that the mapping between physical and transformed region will be one-to-one. In addition, we also apply intuitive position and ratio constraints while solving the PDE to achieve solution-dependent parameterization. While one may use any general PDE with any constraint, the PDEs and additional constraints selected in our case are such that the resulting solution can be efficiently solved through a system of linear equations with or without additional linear constraints. This approach is different from typical existing parameterization methods in IGA that are often solved through an expensive nonlinear optimization processes. The results show that the proposed method can efficiently produce satisfactory analysis-suitable parameterization.
Abstract
This paper presents some results for PDE-based and solution-dependent parameterization of computational domains for isogeometric analysis (IGA) using non-uniform rational [...]
- 22
- read
Numerical Study of the Impact on Cooling Behaviour of Vent-Holes Design of Fruit Packaging Boxes
J. Curto, A. Ilangovan, P. Gaspar, P. Silva, N. Alves

WCCM-ECCOMAS2020.

Abstract
Fruit packaging is highly impactful in temperature and humidity distribution for post-harvest fresh fruit cooling and conservation. It is also the most flexible part on fruit production and should be able to increase shelf-life by deterring adverse conditions. However, packaging thermal performance during cold storage where it is displayed to the internal airflow is influenced by the size and position of vent-holes. In this study airflow inside the fruit box and fruits thermal behaviour are evaluated in order to determine differences between three prototype package models. Models differ from each other by the holes' configuration and size. Furthermore, models' performance is compared to a commercially available box of similar dimensions. Experimental tests were carried out in three stacked boxes, equipped with temperature and humidity sensors, and placed inside a refrigeration chamber with forced air cooling. Packaging box CFD models of each prototype wall were developed to predict the airflow and heat transfer during storage in a refrigeration chamber. Models predictions were validated by experiments results. According with simulation results, all prototypes have better thermal performance than the commercial model. However, significant performance differences between prototype walls were found. Fewer, larger and strategically distributed vent-holes allow a better result in terms of temperature homogeneity inside the box and significant reduction in fruit cooling times.
Abstract
Fruit packaging is highly impactful in temperature and humidity distribution for post-harvest fresh fruit cooling and conservation. It is also the most flexible part on fruit [...]
- 10
- read
A Comparative Study of Numerical Simulation Strategies in Injection Molding
D. Anders, M. Baum, J. Alken

WCCM-ECCOMAS2020.

Abstract
Injection molding is undoubtedly one of the most widely used manufacturing process for polymers [1]. Due to its advantages over other processing techniques in terms of good surface finish, the ability to process complex parts without the need of secondary operations and low cost for mass production, injection molding has found its way into various branches of industry [2]. In order to implement process optimization a thorough understanding of the material behaviour, the filling process as well as the underlying physical phenomena is necessary. A suitable way to address this issue is the application of numerical simulation of the injection molding process. In this contribution the authors present a comparative numerical study for the mold filling of an exemplary thin-walled mold geometry. For the numerical simulation the commercial CFD software packages Cadmould 3D-F and ANSYS CFX are employed. Both software packages make use of different simulation strategies. In ANSYS CFX a threedimensional classical CFD approach is employed for the mold filling which is of course very accurate but extremely time consuming. Cadmould 3D-F pursues a generalized HeleShaw approximation (also known as so-called 2,5-D approach) which is due to some simplifying assumptions less accurate but also much less costly in terms of computing time [2, 5].
Abstract
Injection molding is undoubtedly one of the most widely used manufacturing process for polymers [1]. Due to its advantages over other processing techniques in terms of good [...]
- 8
- read
Numerical Analysis of Viscoelastic Fluid Injection Processes
A. Amani, E. Schillaci, D. Kizildag, C. Pérez-Segarra

WCCM-ECCOMAS2020.

Abstract
In this work, a numerical framework aimed at simulating high-viscous and viscoelastic liquid injection processes is presented. The study of fluid injection with viscoelastic properties, such as synthetic polymers, is of great importance in several industrial sectors, such as metallurgy, automotive, food and pharmaceutical products, ink jets and 3D printers. In this work, a Geometrical Volume-of-Fluid (VOF) method is used to represent the interface, while finite-volume discretizations of Navier-Stokes equations on collocated unstructured meshes are solved through a fractional step method. Viscoelastic constitutive models are used to resolve the non-Newtonian behaviours. The employed implementations allow integrating different types of constitutive equations and stabilization approaches. The test case proposed in the current work consist of the simulation of the discharge of a polymeric jet from an upper nozzle into an air-containing cavity with solid bottom. The behavior of the fluids under analysis is validated by observing the onset of fluid-buckling s tructures. This phenomenon consists in the appearance of toroidal oscillation (as coiling and folding patterns) after the jet hits the solid surface and begins to accumulate. First, low-Reynolds number (Re<1.0) Newtonian Jets are analyzed, with the objective of validating the appearance of buckling as a function of Re number, and comparing the results with reference works. Finally, the potentialities of the proposed numerical methods are shown by simulating buckling phenomena in viscoelastic Jets described with the Oldroyd-B constitutive model.
Abstract
In this work, a numerical framework aimed at simulating high-viscous and viscoelastic liquid injection processes is presented. The study of fluid injection with viscoelastic [...]
- 11
- read
Three-Dimensional Direct Numerical Simulation (dns) of Taylor Bubbles Rising in Non-Newtonian Environments
A. Amani, J. Castro, A. Oliva

WCCM-ECCOMAS2020.

Abstract
Three-dimensional numerical simulation of Taylor gas bubbles as primary unites of slug flow patterns rising in non-Newtonian environments is performed in the context of Direct Numerical Simulation (DNS) of the governing equations, where the whole physics of fluid motions will be taken into account. State-of-the-art numerical tools are proposed to tackle the numerical challenges in the DNS study of this problem. E.g. a coupled level-set volume-of-fluid (CLSVOF) interface capturing method is used to solve the topological changes of the interface. Physical formulations are integrated with moving-mesh (MM) technique to decrease the computational cost of 3D simulations and adaptivemesh-refinement (AMR) technique to increase the local accuracy around the interface. The governing equations are solved using High-Performance Computing (HPC) parallel approaches. To the best of the authors' knowledge, this is the first work dealing with three-dimensional direct numerical simulation of Taylor bubbles rising in non-Newtonian environments.
Abstract
Three-dimensional numerical simulation of Taylor gas bubbles as primary unites of slug flow patterns rising in non-Newtonian environments is performed in the context of Direct [...]
- 11
- read
Lattice-Particle Microstructural Model for Ion Diffusion in Graphite Electrode Batteries
J. Marin-Montin, C. Fresneda-Portillo, F. Montero-Chacón

WCCM-ECCOMAS2020.

Abstract
In this work, we propose a lattice-particle approach to study ionic diffusion across graphite electrodes. In our approach, we generate virtual representative volume elements (RVE) of the electrode material based on its composition, i.e., active particles, carbon additives, and binder. Porosity is also accounted for as an input parameter. To account for the evolution of the ionic concentration, Fickean diffusion is considered. This problem is solved within a network of one-dimensional elements, which is constructed upon the particles of the RVE, yielding a three-dimensional lattice. We use the centraldifference time-integration scheme to solve the transient problem within the framework of the finite element method for the spatial discretization. One of the main advantages of our approach is that we are able to reduce the number of degrees of freedom and thus the computational cost in comparison to the conventional continuum-based finite element simulations. For the transport simulations, we consider Li ions, although our approach can be also applied to other type of species, such as PF-
Abstract
In this work, we propose a lattice-particle approach to study ionic diffusion across graphite electrodes. In our approach, we generate virtual representative volume elements [...]
- 7
- 5/5
- read
Subgrid-Scale Model Based On the Invariants of the Gradient Model Tensor
D. Folch, F. Trias, A. Gorobets, A. Oliva

WCCM-ECCOMAS2020.

Abstract
The incompressible Navier-Stokes equations stand as the best mathematical model for turbulent flows. However, direct numerical simulations at high Reynolds numbers are not yet feasible because the convective term produces far too many relevant scales of motion, thus remaining limited to relatively low-Reynolds numbers. Dynamically less complex mathematical formulations have been developed for coarse-grain simulations, like the well known eddy-viscosity models. Most of these models are based on the combination of invariants of a symmetric tensor that depends on the gradient of the resolved velocity field, G=∇u, and should properly detect different flow configurations (laminar and 2D flows, near-wall behavior, transitional regime, etc.). Brand-new models have been constructed considering the first three invariants of the symmetric tensor GGT with excellent results;1 hence, it is formally based on the lowest-order approximation of the subgrid stress tensor, τ(u) = (∆2/12) GGT +O(∆4), i.e. the gradient model proposed by Clark.2 Furthermore, these models have been implemented on a pseudo-spectral algorithm with a fully-explicit second-order time-integration method.3 The performance of this special configuration has been successfully tested for decaying isotropic turbulence and a turbulent channel flow. It is currently being developed for a semi-infinite boundary layer with periodic conditions as a previous step to carry out wind farm simulations. Details of the implementation and numerical results will be presented.
Abstract
The incompressible Navier-Stokes equations stand as the best mathematical model for turbulent flows. However, direct numerical simulations at high Reynolds numbers are not [...]
- 8
- read
Model Order Reduced Transient Acoustic Finite Element Simulations with Impedance Boundary Conditions
M. Miller, S. Ophem, E. Deckers, W. Desmet

WCCM-ECCOMAS2020.

Abstract
The efficiency of time-domain acoustic simulations is improved immensely if the degrees of freedom imparted by the spatial discretization are reduced. Time-stable model order reduction strategies achieve this by ensuring that frequency-domain system realizations transform in a physical manner after reduction. Frequency dependent damping matrices add to the challenge and require additional consideration. Handling the associated boundary condition imposition in the time domain is one way to approach the problem. Convolution complicates this strategy but the obstacle can be surmounted using a recursive formulation. This work proposes such a method, combining Krylov subspace projection based model order reduction with an efficient time domain-impedance boundary condition implementation. The mass and stiffness matrices are frequency independent and reduced using a second-order Arnoldi algorithm. As Arnoldi iterations implicitly match the moments of the system transfer function, the complex damping matrix must still be contended with. Discussion is included on when reduced order, time-domain, simulations bear fruit. Comparison with full-order time and frequency-domain calculations demonstrate the effectiveness of the proposed algorithms with system degrees of freedom decreased from NDOF= 13125 to RDOF= 63, and simulation time reductions of 9197%.
Abstract
The efficiency of time-domain acoustic simulations is improved immensely if the degrees of freedom imparted by the spatial discretization are reduced. Time-stable model order [...]
- 24
- read
Parallel Computing In Hp-Ocp
G. Kazakis, S. Sotiropoulos, N. Kallioras, S. Xynogalas, C. Mitropoulou, S. Chatzieleftheriou, S. Damikoukas, P. Tsakalis, N. Lagaros

WCCM-ECCOMAS2020.

Abstract
The most computationally demanding part of structural design optimization is the solution of the FE equations and design of the structural model. Therefore, there is a need for the implementation of strategies that can reduce the computational cost of each iteration and thus manage to achieve the same optimized result with considerable reduction in the optimization time. High Performance Optimization Computing Platform (HP-OCP) is an optimization software developed in C# programming language by ISAAR-NTUA and OptiStructre Ltd. [1] which provides a holistic optimization approach for civil engineering structures. It combines powerful derivative-based and derivative-free optimization algorithms like the Projected Quasi-Newton (PQN), Constrained Optimization by Linear Approximation (COBYLA), Latin Hypercube (LH), Differential Evolution etc. [2] integrated with different structural analysis software's like SAP2000, ETABS & SCIA Engineer utilizing their abilities in finite element analysis and most importantly different design codes into the optimization procedure. To deal with the computational demand deriving from this coupling of optimization algorithms and commercial structural analysis software's parallel computational procedures have been implemented to HP-OCP. These procedures were tested in real world civil engineering problems and produced very good results. Parallel strategies are implemented both at the level of the optimization algorithm, by exploiting the natural parallelization features of the evolutionary algorithms, as well as at the level of the repeated structural analysis problems that are required by the optimization algorithm. The numerical tests presented demonstrate the computational advantages of the proposed parallel strategies, which become more pronounced in large-scale optimization problems. 1
Abstract
The most computationally demanding part of structural design optimization is the solution of the FE equations and design of the structural model. Therefore, there is a need [...]
- 11
- read
The Inertial Cavitation Threshold in Soft Tissue Using a Dual-Frequency Driving Signal
T. Filonets, M. Solovchuk, T. Sheu

WCCM-ECCOMAS2020.

Abstract
High Intensity Focused Ultrasound (HIFU) is a non-invasive technology that can be applied for treatment of different diseases and ablation of tumours in different parts of the body. When high intensity ultrasound propagates through the medium bubbles can be formed, a phenomenon known as acoustic cavitation. There are two different regimes of acoustic cavitation: stable cavitation when a bubble just oscillates around an equilibrium state, and inertial cavitation which is accompanied by bubble collapse. These two different regimes can be used for different biomedical applications. However, in some cases it can also make the treatment less predictable. Therefore, fundamental understanding of these effects is very important. In the current study theoretical investigation of the bubble dynamics in viscoelastic medium is performed and inertial cavitation thresholds have been calculated. To describe the bubble dynamics, Gilmore-Akulichev-Zener model has been used, which is suitable for a large bubble oscillations and high ultrasound powers. The results showed that using the dual-frequency driving signal the threshold value of inertial cavitation can be significantly reduced compared to single-frequency signal mode. Large difference between frequencies in the dual-frequency signal leads to lower threshold values. Numerical simulations also showed the dependencies of the cavitation threshold on the bubble radius.
Abstract
High Intensity Focused Ultrasound (HIFU) is a non-invasive technology that can be applied for treatment of different diseases and ablation of tumours in different parts of [...]
- 30
- read
Computation of Eigenvalues in Linear Elasticity with Least-Squares Finite Elements: Dealing with the Mixed System
L. Alzaben, F. Bertrand, D. Boffi

WCCM-ECCOMAS2020.

Abstract
In this paper we discuss some aspects related to the practical implementation of a method that has been introduced recently for the approximation of the eigenvalues of the linear elasticity problem. The scheme, based on a least-squares finite element formulation, gives rise to a non-symmetric discrete formulation that may have complex eigenvalues. Moreover the algebraic eigenvalue problem to be solved is singular, so that the theoretical estimates about the convergence of the scheme should be carefully interpreted.
Abstract
In this paper we discuss some aspects related to the practical implementation of a method that has been introduced recently for the approximation of the eigenvalues of the linear [...]
- 11
- read
On Preconditioning Variable Poisson Equation with Extreme Contrasts in the Coefficients
À. Alsalti-Baldellou, F. Trias, A. Gorobets, A. Oliva

WCCM-ECCOMAS2020.

Abstract
It is well known that the solution by means of iterative methods of very ill-conditioned systems leads to very poor convergence rates. In this context, preconditioning becomes crucial in order to modify the spectrum of the system being solved and improve the performance of the solvers. A proper balance between the reduction in the number of iterations and the overhead of the construction and application of the preconditioner needs to be sought to actually decrease the total execution time of the solvers. This is particularly important when considering variable coefficients matrices as, in general, its preconditioners will also be variable and need to be updated regularly at an affordable cost. In this work we present a family of variable preconditioners designed for the effective solution of variable Poisson equation with extreme contrasts in the coefficients, which represents a particularly challenging case as it translates into a variable and extremely ill-conditioned system arising in many situations such as with multiphase flows presenting high density ratios or in the presence of highly-stretched adaptive mesh refinements. Finally, the results of the numerical experiments performed are presented and discussed, confirming our preconditioners as extremely affordable, highly-parallelizable and easy-to-implement alternatives to the more standard (and usually unfeasible) preconditioners, still showing great improvements in the rate of convergence of the solvers without requiring the variable coefficients matrix to be explicitly rebuilt at each iteration.
Abstract
It is well known that the solution by means of iterative methods of very ill-conditioned systems leads to very poor convergence rates. In this context, preconditioning becomes [...]
- 22
- read
Full-Waveform Inversion of Sh-Wave Input Motions in a Near-Surface 2d Domain
B. Guidio, C. Jeong

WCCM-ECCOMAS2020.

Abstract
There is a need to estimate complex seismic input motions in a near-surface domain, without resorting to the hypocenter, from restricted seismic measurement data. Thus, engineers can replicate responses within structures and soils after an earthquake occurrence by using the estimated seismic inputs and evaluate the impact of an earthquake on the built environment. To date, there has been no robust numerical method that can identify complex seismic input motions in a solid, truncated by a wave-absorbing boundary condition. Existing methods are limited to large-scale seismic-source inversion approaches and deconvolution. To fill this gap, a new inversion modeling method is presented for reconstructing complex, incoherent SH wave input motions in a two-dimensional (2D) domain that is truncated by a wave-absorbing boundary condition (WABC), using a partial differential equation (PDE)constrained optimization method. In a set of numerical examples, targeted, dynamic traction at the WABC mimics seismic incident wavefield. The discretize-then-optimize (DTO) approach is used in the mathematical modeling and numerical implementation, and the finite element method (FEM) is applied to solve state and adjoint problems. The numerical results indicate that the presented inversion algorithm can reconstruct incident, inclined plane waves. Second, the accuracy of our inversion solver is not compromised by the material complexity of a background domain. Lastly, the minimizer suffers less from solution multiplicity when it identifies lower frequency traction (e.g., a realistic seismic signal).
Abstract
There is a need to estimate complex seismic input motions in a near-surface domain, without resorting to the hypocenter, from restricted seismic measurement data. Thus, engineers [...]
- 12
- read
Non-Modal Analysis of Multigrid Schemes for the High-Order Flux Reconstruction Method
A. Hurtado-de-Mendoza, J. Kou, S. Joshi, K. Puri, C. Hirsch, E. Ferrer

WCCM-ECCOMAS2020.

Abstract
The present study introduces an application of the non-modal analysis to multigrid operators with explicit Runge-Kutta smoothers in the context of Flux Reconstruction discretizations of the linear convection-diffusion equation. A dissipation curve is obtained that reflects upon the convergence properties of the multigrid operator. The number of smoothing steps, the type of cycle (V/W) and the combination of pand h-multigrid are taken into account in order to find those configurations which yield faster convergence rates. The analysis is carried out for polynomial orders up to P = 6, in 1D and 2D for varying degrees of convection (Péclet number), as well as for high aspect ratio cells. The non-modal analysis can support existing evidence on the behaviour of multigrid schemes. W-cycles, a higher number of coarse-level sweeps or the combined use of hp-multigrid are shown to increase the error dissipation, while higher degrees of convection and/or high aspect-ratio cells both decrease the error dissipation rate.
Abstract
The present study introduces an application of the non-modal analysis to multigrid operators with explicit Runge-Kutta smoothers in the context of Flux Reconstruction discretizations [...]
- 10
- read
A Wall-Modeled/Wall-Resolved Les Method For Turbulent Wall Flows
F. Vanna, M. Cogo, M. Bernardini, F. Picano, E. Benini

WCCM-ECCOMAS2020.

Abstract
This present paper reports a novel methodology to simulate wall-bounded flows in the LargeEddy Simulation framework using an automatic transition between a wall-modelled and a wall-resolved approach. The proposed technique aims at prescribing the right (modelled or resolved) wall shear stress and wall heat flux, preserving at the same time the no-slip/no-penetration conditions for the velocity and the isothermal/adiabatic conditions for the temperature fields. The approach is successfully implemented in a high-order finite-difference framework, and it is found able to adapt smoothly to the available nearwall grid spacing. Thus, the method falls into the wall-resolved case, when the near-wall dynamic is directly computed, whereas it employs the wall stress model when a full resolution of the near-wall region is not achievable. The method is tested on a nearly-incompressible turbulent channel flow and a supersonic spatially-devolving boundary layer flow. The obtained results highlight an excellent accuracy in representing the wall turbulence dynamics in terms of mean velocity profiles and fluctuations, almost independently of the near-wall spatial resolution. Thus, the proposed method results in a promising technique for analysis of high-Reynolds wall-bounded flows.
Abstract
This present paper reports a novel methodology to simulate wall-bounded flows in the LargeEddy Simulation framework using an automatic transition between a wall-modelled and [...]
- 25
- read
Investigating The Capabilities Of Cfd-Based Data-Driven Models For Indoor Environmental Design And Control
N. Morozova, F. Trias, R. Capdevila, A. Oliva

WCCM-ECCOMAS2020.

Abstract
In this work, we study the accuracy of CFD-based data-driven models, which predict comfortrelated flow parameters in a ventilated cavity with a heated floor. We compare the computational cost and accuracy of three different models, namely artificial neural network, support vector regression, and gradient boosting regression. The tested scenarios include short and long cavities with different inlet velocities. Among the studied frameworks, the artificial neural network provides the most accurate predictions for most of the tested flow configurations. However, test configurations with jet separation and a secondary vortex are more difficult to predict correctly; thus more high-fidelity data is required in order to construct a more robust and reliable model.
Abstract
In this work, we study the accuracy of CFD-based data-driven models, which predict comfortrelated flow parameters in a ventilated cavity with a heated floor. We compare the [...]
- 7
- read
Simultaneous Biomass Flame Temperature and Spectral Emissivity Estimation by Inverse Method
C. Salinas, L. Gonzales, P. Cândido, D. Marcelo-Aldana

WCCM-ECCOMAS2020.

Abstract
This work presents a numerical and experimental study in a laboratorial rig to measure local flame temperature during the biomass combustion. Most of practical situations of biomass combustion is a transient and non-stable phenomenon. A spectral analysis method in the visible spectrum range was used for measurement of local temperatures in a sugarcane bagasse flame. The two-color method considering spectral flame emissivity variation joint with an inverse numerical procedure is used to simultaneous estimation of flame temperature and spectral emissivity. Also, a sequential iterative numerical procedure is tested. It includes the application of two-color method considering grey emissivity between two near spectral intensity data in a first step. In a second step and iterative procedure considering spectral emissivity variation along whole spectral range is applied. To apply the numerical methods was used several spectral intervals in the visible range. In inverse estimation, the LevenbergMarquardt method is used and a polynomial and a cosine models for estimation of coefficients were tested. Six test of biomass combustion were made focusing the sensor to the reaction zone, and spectral data collected. The flame temperature estimated for the six spectral data were closed to 1530 to 1540 K. Also, the spectral emissivity is closed in all sixcombustion test. Instantaneous bagasse flame temperature in the range of 1400 K to 1500 K in the reactive combustion zone were found. An image processing method was used to measure the temperature in the same reactive combustion zone and the comparison show
Abstract
This work presents a numerical and experimental study in a laboratorial rig to measure local flame temperature during the biomass combustion. Most of practical situations [...]
- 6
- read
Simulation Of Fluid-Structure Interaction And Impact Force On A Reed Valve
P. Castrillo, E. Schillaci, J. Rigola

WCCM-ECCOMAS2020.

Abstract
The cyclic impact force between a reed valve and the seat plate is the main reason of the valve failure in many thermo-technical devices as compressors, engines, etc. According to experimental observations the latter is due to fatigue and usually occurs in the leading part of the valve `neck'. In this work, a complex numerical analysis is presented aimed to studying the external forces and internal stresses suffered by the valve. In particular, the impact force between the valve and the seat is studied. The numerical analysis relies on the coupled synergy of two different simulation concepts. In order to do so, two codes are used: (1) first, the in-house Computational Fluid Dynamics (CFD) code presented in [1] is employed to simulate the Fluid-Structure Interaction (FSI) between gas and valve, extracting reference data for valve displacement and external gas pressures; (2) second, the analysis of the internal structure stresses, together with the impact forces with the plate is implemented in a Computational Solid Dynamics (CSD) code developed in FreeFEM++ [2]. The impact force representation is based on the formulation presented in [3] where a conserving algorithm for frictionless dynamic contact/impact is developed. Due to the importance of obtaining an adequate impact force, an exhaustive study is carried out on its characterization in terms of numerical parameters, such as the penalty stiffness. Under this framework, the valve displacement and impact velocities are verified. Hence, impact forces are analysed in different scenarios, obtaining interesting observations about stresses distribution, with a particular focus on the points where failure is experienced.
Abstract
The cyclic impact force between a reed valve and the seat plate is the main reason of the valve failure in many thermo-technical devices as compressors, engines, etc. According [...]
- 9
- read
Stiffness and Damping Properties of Cellular Shear Walls with Viscoelastic Fillers
P. Syrimi, S. Papathanasiou, P. Tsopelas

WCCM-ECCOMAS2020.

Abstract
Aiming to enhance the stiffness and damping properties of cellular materials, the idea of inserting viscoelastic fillers into the gaps of cells is developed. A representative volume element of typical hexagonal honeycomb cells with and without inclusions is examined under dynamic shear loading, in order to define the mechanical properties of the proposed composite material. Within each cell the viscoelastic material is either through viscoelastic shell elements or through viscoelastic bars attached to the hexagonal cell apexes. The calibrated viscoelastic bar RVE model is, by extent, implemented to a shear wall, tested under dynamic shear loading, to verify the efficiency of the viscoelastic inclusions. Finally, different viscoelastic filler patterns are explored to study the proposed system's performance.
Abstract
Aiming to enhance the stiffness and damping properties of cellular materials, the idea of inserting viscoelastic fillers into the gaps of cells is developed. A representative [...]
- 7
- read
Hybrid Numerical-Analytical Scheme for Locally Inhomogeneous Elastic Waveguides
E. Glushkov, N. Glushkova, A. Evdokimov

WCCM-ECCOMAS2020.

Abstract
Numerical simulation of guided wave excitation, propagation, and diffraction in laminate structures with local inhomogeneities (obstacles) is associated with high computational cost due to the need for a mesh-based approximation of extended domains with a rigorous account for the radiation conditions at infinity. To obtain computationally efficient solutions, hybrid numerical-analytical approaches are currently being developed, based on linking a numerical solution in a local vicinity of the source and/or obstacles with an explicit analytical representation in the external semi-infinite domain. However, the developed methods are generally not widely spread because the possibility of such coupling with an external multimode wave field is generally not provided in standard finite-element (FE) software. We propose a scheme that allows the use of the FE software as a black box for the required correct matching of local numerical and global analytical solutions (FEM-An). The FEM is used to obtain a set of local numerical solutions that serve as a basis in the inner domain. These solutions satisfy the boundary conditions induced by guided wave modes so that they fit correctly with the modal expansion in the outer region. The expansion coefficients of both FE and modal decompositions are determined then from the condition of stress and displacement continuity at the interface between the inner and outer domains. This scheme was numerically validated against analytical solutions to test problems and FE solutions for long waveguide sections with perfect match layer absorbing conditions at the ends (FEM PML). Along the way, it turned out that the FEM-PML approach gives an incorrect result in the backward-wave bands and at high frequencies. The application of the FEM-An hybrid scheme is illustrated by examples of Lamb wave diffraction by elastic inclusions and delaminations.
Abstract
Numerical simulation of guided wave excitation, propagation, and diffraction in laminate structures with local inhomogeneities (obstacles) is associated with high computational [...]
- 4
- read
Large Scale Discrete Element Simulation Campaigns – Simulating Extraterrestrial Soils in Partsival
R. Lichtenheldt, S. Ono, L. Stubbig

WCCM-ECCOMAS2020.

Abstract
In planetary exploration, testing under the actual mission conditions is inherently not possible. Hence, simulation campaigns complement ground test campaigns. This specially applies to surface missions that include the complex behaviour of soils under non-terrestrial gravitation. Increasingly ambitious mission goals made large simulation campaigns with very precise particle models necessary for the simulation of soil interaction. Thus, to limit the amount of time and the computation hardware needed, DLR developed the particle simulation tool 'Sir partsival'. This tool does not only speed up simulations by usage of GPU computing, but also integrates the institute's experience in modelling of soil on Earth and beyond. Using partsival it was possible to speed up simulations by more than a factor of ten and thus conduct large simulation campaigns. Two examples are shown: a large, on-going validation campaign of DEM for wheel simulations, and the completed traction optimization for the MMX rover wheel.
Abstract
In planetary exploration, testing under the actual mission conditions is inherently not possible. Hence, simulation campaigns complement ground test campaigns. This specially [...]
- 4
- read
Predictive Modeling of Holdup in Horizontal Water-Oil Flow Using a Neural Network Approach
C. Díaz, O. González-Estrada, M. Cely

WCCM-ECCOMAS2020.

Abstract
In this work, the application of an artificial neural network (ANN) is proposed to develop a predicting model for the holdup of a two-phase flow composed of water and mineral oil in a horizontal pipe. For this, the surface velocities of each fluid and the differential pressure in the pipeline are used as input parameters of the multilayer artificial neural network with backpropagation, while the holdup of the fluids is used as the output parameter for the training. A set of 56 experimental data was obtained in the LabPetroCEPETRO-UNICAMP laboratory. The best performing results for the predictive model show a mean absolute error (AAPE) of 3.01% and a coefficient of determination 2 of 0.9964 using 15 neurons in the hidden layer of the network and the TanSig transfer function.
Abstract
In this work, the application of an artificial neural network (ANN) is proposed to develop a predicting model for the holdup of a two-phase flow composed of water and mineral [...]
- 6
- read
Free Surface Flow in Vertical Taylor-Couette System
T. Watanabe, H. Furukawa

WCCM-ECCOMAS2020.

Abstract
Vertical Taylor-Couette flow with a free surface at its top has been examined by numerical and experimental approaches. Being compared with the ideal case with infinite length cylinders, this system has effects of the finite lengths, the gravitational acceleration and the surface tension. The inner cylinder rotates and the outer cylinder and the bottom wall are stationary. The numerical approach uses the VOF method to model the free surface. In the experiment, the flow patterns are observed from the top and the meridional views, and the displacement of the free surface is measured. The flow modes are classified by the number of vortices appearing in the meridional plane and the radial flow directions at the top and the bottom. The power spectra of the displacement of the free surface and the bulk energy are numerically evaluated and the bulk energy tends to give favorable agreement with the spectra of the displacement obtained experimentally. The transition from the secondary mode flow appearing at higher Reynolds number to the primary mode flow at lower Reynolds number is examined and comparison between the numerical and experimental results are made.
Abstract
Vertical Taylor-Couette flow with a free surface at its top has been examined by numerical and experimental approaches. Being compared with the ideal case with infinite length [...]
- 7
- read
Peridynamic Modelling of Harmonic Structured Materials Under High Strain Rate Deformation
A. Ahadi, T. Sousa, D. Orlov, S. Elis

WCCM-ECCOMAS2020.

Abstract
This paper applies the peridynamic continuum mechanics theory on a new type of material known as harmonic-structured materials. Using the Peridigm software, rapid uniaxial elongation is simulated on a peridynamic model of a thin bimodal harmonic-structured metal sheet. Mechanical wave initiation, propagation, and reflection, as well as fracture initiation and propagation are successfully simulated.
Abstract
This paper applies the peridynamic continuum mechanics theory on a new type of material known as harmonic-structured materials. Using the Peridigm software, rapid uniaxial [...]
- 6
- read
Considerations On Relationship Between Each Sensitivity Term and Results of Topology Optimization in Unsteady Oscillation Problem for 3d Cantilever Beam
M. Kishida, T. Kurahashi

WCCM-ECCOMAS2020.

Abstract
In this study, we carried out the topology optimization based on the density method in the unsteady oscillation problems. This paper considers two items, the proposal of new performance function and the relationship between each sensitivity term and results of topology optimization. In other studies, the adjoint analysis must be performed for the propose of minimizing strain energy. Optimization problems are classified into the problem of maximizing work in the negative and the problem of minimizing work in the positive, and self-adjoint relationships are derived. When using this method for cantilever beam model, there is a problem that numerous high-density elements is distributed at the load point, i.e., the tip. By classifying the sensitivity, it is found that the mass term has a great influence on the problem. In our study, some results are shown by changing the penalization parameter in the SIMP method that determines the mass coefficient, damping coefficient and elastic coefficient.
Abstract
In this study, we carried out the topology optimization based on the density method in the unsteady oscillation problems. This paper considers two items, the proposal of new [...]
- 4
- read
A MODEL FOR COUPLED BELOUSOV-ZHABOTINSKY OSCILLATORS WITH DELAY
M. Budroni, D. Conte, R. D‘Ambrosio, G. Pagano, B. Paternoster, F. Rossi, A. Hassan, S. Ristori

WCCM-ECCOMAS2020.

Abstract
The purpose of this work is to introduce time delays in an Oregonator-based model [11] to describe molecular communication in a system of coupled chemical oscillators [5]. With this approach we assume that the molecules responsible for communication (signal molecules) take a non-zero time to cross the barrier (phospholipid membrane) which physically separate the oscillators. In particular, here we present numerical results and we discuss and highlight the differences between delayed and non-delayed models.
Abstract
The purpose of this work is to introduce time delays in an Oregonator-based model [11] to describe molecular communication in a system of coupled chemical oscillators [5]. [...]
- 6
- read
An Equal-Order Finite Element Framework for Incompressible Non-Newtonian Flow Problems
D. Pacheco, O. Steinbach

WCCM-ECCOMAS2020.

Abstract
Various materials of engineering and biomedical interest can be modelled as generalised Newtonian fluids, i.e., via an apparent viscosity depending locally on the flow field. In spite of the particular features of those models, they are often handled in practice by classical numerical techniques originally conceived for Newtonian fluids. Methods designed specifically for the generalised case are rather scarce in the literature, as well as their use in practical applications. As it turns out, tackling nonNewtonian problems with standard finite element formulations can have undesired consequences such as the induction of spurious pressure boundary layers and the emergence of natural boundary conditions not suitable for realistic flow scenarios. In this context, we introduce a novel framework that deals with those issues while maintaining simplicity and low computational cost. The new stabilised formulation is based on a modified system combining the continuity equation with a Poisson equation for the pressure and consistent pressure boundary conditions. A weak enforcement of the rheological law is employed to enable full consistency even for first-order finite element pairs. Simple numerical examples are provided to demonstrate the potential of our method in yielding accurate solutions for relevant problems.
Abstract
Various materials of engineering and biomedical interest can be modelled as generalised Newtonian fluids, i.e., via an apparent viscosity depending locally on the flow field. [...]
- 9
- read
Vibration Control of Offshore Wind Turbine Towers
A. Kampitsis

WCCM-ECCOMAS2020.

Abstract
In this study, the dynamic response of a Offshore Wind Turbine (OWT) tower under wind loading is investigated adopting the vibration absorption concepts. Application of vibration control systems (VCS) to OWT towers has the potential to significantly improve the dynamics of the tower and the nacelle, increasing its reliability. The first vibration mitigation approach is the classical tuned mass damped (TMD) and it is used as benchmark. In the TMD approach an additional mass is installed inside the nacelle at the top of the OWT tower. Usually, a large additional mass is required to achieve significant effective damping, which may lead to undesirable second order effects and ultimately increase the tower displacement. In the second approach, a passive vibration absorption concept is implemented, namely the Extended KDamper (EKD). The KDamper device is a special mechanical configuration consisting of an additional mass, an artificial damper and both positive and negative stiffness elements. The OWT is subjected to extreme environmental loading such as the wind and sea waves. These OWT loadings, are usually calculated based on stochastic interpolation of meteorological data or by artificial time-histories. In this study, the Blade Element Momentum (BEM) theory is taken into consideration in order to produce the wind load time histories. The improvement on the dynamic response of the OWT tower using vibration control systems is illustrated. Moreover, this research shown that the EKD presents a significant increase of the damping ratio with a minimal value of added mass at the top of the tower.
Abstract
In this study, the dynamic response of a Offshore Wind Turbine (OWT) tower under wind loading is investigated adopting the vibration absorption concepts. Application of vibration [...]
- 13
- read
Unsteady Dynamics of Free-Interface Gravitational Liquid Sheet Flows
A. Pia, M. Chiatto, L. Luca

WCCM-ECCOMAS2020.

Abstract
Numerical simulations of gravitational planar liquid sheet flows, interacting with unconfined gaseous environments located on both sides of the liquid phase, are performed through Volume-of-Fluid (VOF) technique. The global unsteady dynamics of the non-parallel flow is analyzed by perturbing the initial steady configuration by means of a Gaussian bump in the transverse velocity component of relatively very small amplitude, thereby exciting sinuous modes. Thanks to the development of a theoretical linear one-dimensional model, more physical insights are gained on the flow system. It is found that surface tension plays a stabilizing role for the gravitational sheet, and for relatively high values of density ratio rρ of gaseous-to-liquid phases it becomes unstable. An analogy is shown between the global unstable behavior exhibited by the liquid sheet as rρ increases, and the shear-induced global instability found by Tammisola et al. [“Surface tension-induced global instability of planar jets and wakes”, J. Fluid Mech. 713, 632–658 (2012)] for planar jet and wake flows of two immiscible fluids in the presence of surface tension.
Abstract
Numerical simulations of gravitational planar liquid sheet flows, interacting with unconfined gaseous environments located on both sides of the liquid phase, are performed [...]
- 9
- read
Surrogate-Assisted Optimization for Augmentation of Finite Element Techniques
S. Bagheri, U. Reinicke, D. Anders, W. Konen

WCCM-ECCOMAS2020.

Abstract
The application of finite element techniques for the analysis and optimization of complex thermo-mechanical structures typically involves highly nonlinear models for material characterization, tribological contact, large deformation, damage, etc. These nonlinearities usually call for a higher-order Spatio-temporal discretization, including a large number of elements and time-steps in order to provide good convergence and sufficiently accurate simulation results. This inevitably leads to many expensive simulations in terms of cost and time if an optimization or adaption of model parameters has to be done. In this work, a FEM simulation modeling approach is proposed, which uses radial basis function interpolations (RBF) as efficient surrogate models to save FEM simulations. Also, a surrogate-assisted optimization algorithm [3] is utilized to find the parameter setting, which would lead to maximum damage in a simple tensile testing scenario involving a notched specimen with as few FEM simulations as possible. The relatively high accuracy of the utilized surrogate models showcases promising results and indicates the potential of surrogate models in saving time-expensive simulations.
Abstract
The application of finite element techniques for the analysis and optimization of complex thermo-mechanical structures typically involves highly nonlinear models for material [...]
- 5
- read
Maneuver-Based Cross-Validation Approach for Angle-of-Attack Estimation
A. Brandl, M. Battipede

WCCM-ECCOMAS2020.

Abstract
The estimation of the Angle of Attack (AOA) and Angle of Sideslip (AOS) is crucial for flight monitoring and control. However, a gap has been identified on the data selection technique for the class of estimators based on data-driven methods, such as the synthetic sensor based on Neural Network (NN). This paper proposes a Cross Validation (CV) technique applied on a manoeuver-based partitioning method to provide evidence that a given selection of data can lead to better estimation performance, with the final aim of providing a list of manoeuvers suitable for the training phase of the estimator. Results are shown using simulated data related to the CleanSky 2 project MIDAS.
Abstract
The estimation of the Angle of Attack (AOA) and Angle of Sideslip (AOS) is crucial for flight monitoring and control. However, a gap has been identified on the data selection [...]
- 6
- read
Application of Level-Set Type Topology Optimization Analysis for Cavity Shape Estimation Problem in Structures Based On Non-Destructive Hammering Test
T. Kurahashi, Y. Murakami, S. Toyama, F. Ikeda, T. Iyama, I. Ihara

WCCM-ECCOMAS2020.

Abstract
In this study, we present application of the level-set type topology optimization analysis for the cavity shape estimation problem in structures based on the non-destructive hammering test. The cavity shape is identified so as to minimize a performance function. The performance function is defined as the square sum of the residual between computed and the observed displacements on structure surface. In this study, accuracy of identified cavity shape is investigated by changing numerical parameters in the topology optimization.
Abstract
In this study, we present application of the level-set type topology optimization analysis for the cavity shape estimation problem in structures based on the non-destructive [...]
- 4
- read
A New Boundary Element Technique For One- And Two-Temperature Models Of Biothermomechanical Behavior Of Anisotropic Biological Tissues
M. Fahmy

WCCM-ECCOMAS2020.

Abstract
The main objective of this paper is to develop a novel boundary element technique for describing the three-dimensional (3D) biothermomechanical behavior of anisotropic biological tissues. The governing equations are studied on the basis of the dual phase lag bioheat transfer and Biot's theory for oneand two-temperature models. Because of the benefits of CQBEM, such as not being restricted by the complex shape of biological tissues and not requiring discretization of the interior of the treated region, it can cope with complex bioheat models and has low use of RAM and CPU. CQBEM is therefore a flexible and efficient tool for modeling the distribution of bioheat in anisotropic biological tissues and associated deformation. The resulting linear equations arising from CQBEM are solved by the generalized modified shift-splitting (GMSS) iterative method which reduces the number of iterations and the total time of the CPU. Numerical findings show the validity, efficacy and consistency of the proposed technique.
Abstract
The main objective of this paper is to develop a novel boundary element technique for describing the three-dimensional (3D) biothermomechanical behavior of anisotropic biological [...]
- 7
- read
Comparison Study Between Using Air Data And Soil Data In Predicting Soil Temperature
S. Hamdane, L. Pires, P. Dinho, A. Moummi

WCCM-ECCOMAS2020.

Abstract
soil in shallow depths could be used as heat sink / heat source for cooling / heating devices. To predict its thermal behaviour during year, some authors use data from soil temperature records some else use data from average yearly air temperature measurements of meteorologic stations to avoid measure soil temperature in shallow depths for long time. In this research paper, we present a comparison study about the use of soil data and the use of air data for soil temperature prediction, where transient heat conduction in semi-infinite solids is utilised in the study. We measured the soil temperature of three sites for three years in Covilha -Portugal to support the current study. Later, we create reference years of soil temperature in the three sites.
Abstract
soil in shallow depths could be used as heat sink / heat source for cooling / heating devices. To predict its thermal behaviour during year, some authors use data from soil [...]
- 8
- read
Selected Problems Regarding Construction Of Lightweight Motorcycle Swing Arm
W. Nosko, M. Mrzygłód

WCCM-ECCOMAS2020.

Abstract
Manufacturers of light motorcycles have always sought the use of more and more advanced materials. The application of carbon fiber in motorcycle parts design can bring high strength, stiffness, and low weight in a manner previously unavailable. This article contains selected problems regarding the construction of a swingarm made of carbon fiber reinforced plastics (CFRP). The swingarm was made specifically for Honda NSF 100. For the analysis, series of tests such as geometry, FEM, and structural analysis was chosen. The analysis concludes that the swingarm shape that was adopted from the original steel made part is not suitable for CFRP use. The analysis showed manufacturing errors, like the inconsistent thickness of laminate and poor quality of assembling. The full potential of the CFRP material was not used.
Abstract
Manufacturers of light motorcycles have always sought the use of more and more advanced materials. The application of carbon fiber in motorcycle parts design can bring high [...]
- 3
- read
A Coupled Diffusion-Mechanical Lattice Model for the Degradation of Graphite Active Particles of Li-Ion Battery Anodes
J. Marin-Montin, F. Montero-Chacón

WCCM-ECCOMAS2020.

Abstract
The performance,durability of lithium-ion batteries (LIBs) are constrained by the degradation mechanisms that take place during charge,discharge cycles. Degradation of active particles (APs) of LIBs is a complex problem involving several physical phenomena (e.g.,diffusion,mechanical deformation,heat transfer,to cite a few). During lithium insertion,extraction cycles,volume changes in the AP result in high mechanical stresses and,consequently,mechanical damage that promotes capacity fade. In this work,we present a microscale 3D finite element model that takes into account the coupled effects between lithium diffusion,mechanical stress within the AP. Using the surface of an ellipsoid as the base for the geometrical construction,we are able to generate different shapes of APs,with both concave,convex surfaces. Porosity,other types of defects that may be present inside the AP are explicitly modeled,different volume fractions,shapes,orientations are also accounted for. In our approach,the material is discretized into a lattice of one-dimensional elements: we consider beam elements for the mechanical problem,while in the diffusive approach,the material is treated as an assembly of 'nanopipes' through which the flow of Li-ions takes place. The same lattice network is used for both simulations. We follow a classical lattice model approach to characterize the fracture behavior of a single AP of a LIB anode when subjected to charge/discharge cycles. The material of the APs analyzed in this work is graphite,which presents a brittle,disordered material structure,making it suitable for lattice modeling. The mechanical problem is solved,obtaining the crack patterns associated with specific charge,discharge strategies,potential initial defects. The simulation results correctly reproduce the experimental observations on mechanical stresses,the evolution of damage. This lattice model framework analyzing the degradation in the APs of LIBs (durability) can be used to provide more information regarding the microstructural evolution,morphological changes,mechanical degradation in APs,identify improvement strategies.
Abstract
The performance,durability of lithium-ion batteries (LIBs) are constrained by the degradation mechanisms that take place during charge,discharge cycles. Degradation of active [...]
- 14
- 5/5
- read