Documents published in Scipedia

• Model-based design approach finding optimal liquid cooling flow path for electric vehicle battery

  D. Aketo, K. Ono

  WCCM2024.

  
  

  Abstract

  Today many electronic devices that generate significant heat are required to be equipped with liquid cooling systems to reduce their temperature. Since the liquid flow path in the cooling system affects cooling performance, determining flow path in the early development phase can improve the efficiency of the downstream development process and reduce the total cost. In this paper, we propose a model-based analysis system for thermo-fluid phenomena based on CFD results and demonstrate the parametric optimization of flow path.

  Abstract
  Today many electronic devices that generate significant heat are required to be equipped with liquid cooling systems to reduce their temperature. Since the liquid flow path [...]

  • 0
  •  read

• Input-output reduced order modeling for public health intervention evaluation

  A. Viguerie, E. Jacobson, C. Piazzola, H. Islam

  WCCM2024.

  
  

  Abstract

  In recent years, mathematical models have become an indispensable tool in the planning, evaluation, and implementation of public health interventions. Models must often provide detailed information for many levels of population stratification. Such detail comes at a price: in addition to the computational costs, the number of considered input parameters can be large, making effective study design difficult. To address these difficulties, we propose a novel technique to reduce the dimension of the model input space to simplify model-informed intervention planning. The method works by first applying a dimension reduction technique on the model output space. We then develop a method which allows us to map each reduced output to a corresponding vector in the input space, thereby reducing its dimension. We apply the method to the HIV Optimization and Prevention Economics (HOPE) model, to validate the approach and establish proof of concept.

  Abstract
  In recent years, mathematical models have become an indispensable tool in the planning, evaluation, and implementation of public health interventions. Models must often provide [...]

  • 0
  •  read

• Computational model for local buckling of compressively loaded omega-stringer-stiffened panels

  C. El Yaakoubi-Mesbah, C. Mittelstedt

  WCCM2024.

  
  

  Abstract

  Thin-walled composite structures are used in applications such as aircraft and spacecraft due to their low weight and corresponding high stiffness properties. To optimize the potential of these structures to the fullest extent, a complete understanding of their stability behavior is required. Thereby, uniaxial compression describes an important load case that is investigated. A closed-form analytical method based on the energy method for determining the local buckling load of omega-stringer-stiffened panels is presented. The stiffened panel under consideration consists of the skin plate with eccentrically attached stringer feet along the longitudinal sides of the panel, while the remaining part of the omega-stringer is modeled by corresponding elastically restrained edges. Due to the applied stringer feet, stiffness discontinuities occur in the stiffened panel. This is covered by the presented method, whereas in comparable studies in the literature, a homogeneous stiffness is often assumed across the entire panel. To evaluate the new analysis method, a comparison with the numerical solution of the corresponding Levy-type ´solution and the finite element analysis is being drawn.

  Abstract
  Thin-walled composite structures are used in applications such as aircraft and spacecraft due to their low weight and corresponding high stiffness properties. To optimize [...]

  • 0
  •  read

• Developments in the use of the Bonded Particle Model to study ore fracture

  T. Oladele, D. Weatherley, L. Bbosa

  WCCM2024.

  
  

  Abstract

  In mineral processing, ore fracture is an essential first step for which the objective is to increase the exposed surface area of the valuable mineral, thereby increasing the likelihood of liberation in subsequent separation stages. This process is well known to be energy-intensive, and increasing scrutiny around sustainable practices has heightened the need to examine the efficiency of current industry approaches. Factors such as mineralogical structure and inherent weakening in the form of micro cracks are known to affect ore breakage mechanisms. However, isolating and investigating individual factors under experimental conditions is challenging and typically impractical. Numerical techniques such as the Bonded Particle Model-Discrete Element Method (BPM-DEM) have been developed as a means of investigating in isolation, the effects of different factors on ore breakage behaviour under closely controlled breakage conditions. In this work, the robustness of the BPM-DEM in predicting fracture characteristics during SILC impact breakage is evaluated. Thereafter, the BPM-DEM is used to analyse the internal mechanical response of a simulated rock specimen under impact loading commensurate with that of the SILC. The method is shown to be an insightful opportunity to study intrinsic and extrinsic rock properties during dynamic loading and breakage

  Abstract
  In mineral processing, ore fracture is an essential first step for which the objective is to increase the exposed surface area of the valuable mineral, thereby increasing [...]

  • 0
  •  read

• Effect of particle-size-scaling on particle interactions in DEM-simulations of sand in the context of air pluviation

  N. Heim, S. Henke

  WCCM2024.

  
  

  Abstract

  DEM (discrete element method) is a widely used numerical simulation method, which models the behaviour of a bulk substrate based on the individual interactions of many particles. One of its possible applications is the modelling of sand behaviour in different laboratory tests, e.g., cone penetration tests [1] or direct shear tests [2]. Furthermore, DEM is specifically of interest as a modelling method for investigating air pluviation, because it models the individual inter-particle and particle-environment interactions, both friction and collisions, which determine the compaction and homogeneity of the created samples. However, one disadvantage of DEM is the relatively long computational time [3] especially with decreasing particle sizes. This makes larger particle sizes compared to reality more interesting, especially for large scale or repeating simulations. On the other hand, if the size of the chosen particles is too large, certain interactions, such as interactions with other materials and equipment, may not be simulated in a way that properly represents real behaviour. This would lead to preferring smaller sized particles, which again would lead to longer computational times. Therefore, the chosen particle size as an important aspect of DEM simulations will be discussed, as well as the effects on different simulation aspects. This includes necessary parameter calibrations, the resulting inter-particle and particle-environment interactions as well as the achieved simulation results and accuracies. Of specific interest is the largest particle size, at which accurate and realistic results concerning real-world particle interactions can be achieved. Further, the effects of graded particle sizes to better represent the sand during the pluviation process will be discussed.

  Abstract
  DEM (discrete element method) is a widely used numerical simulation method, which models the behaviour of a bulk substrate based on the individual interactions of many particles. [...]

  • 0
  •  read

• Full waveform modeling in seismic exploration based on a digital geological model using spectral element method on GPU

  A. Vershinin, A. Yury, V. Levin, K. Petrovsky

  WCCM2024.

  
  

  Abstract

  The paper considers the solution of a three-dimensional problem of modeling of all types of seismic waves propagating in real geological media. The numerical algorithm based on the spectral element method (SEM). The main advantages of SEM (high order space discretization, explicit time integration scheme) are presented in comparison with the classical approach based on the finite element method (FEM). The features of the massively parallel implementation of the algorithm on modern MultiGPU systems (based on A100 GPU) using CUDA technology are considered. The efficiency of parallelization on hybrid systems with different SEM orders and parameters of the numerical time integration scheme is analyzed. The results of solving a three-dimensional problem of modeling the propagation of seismic waves in a heterogeneous geological media with faults and sharply varying properties of layers are presented. Analysis of the numerical convergence of SEM for dispersive waves of the Rayleigh type is performed. Local and non-local non-reflective boundary conditions on the artificial boundary of the computational region are considered. The 3D computational model is constructed using a detailed digital geological model built for one of the Arctic regions. It was converted to an unstructured hexahedral mesh to perform SEM calculations using CAE FIDESYS software. The model is further generalized for typical seismic-geological conditions of Western Siberia, so that on the basis of such modeling it is possible to conduct a wide range of studies on the possibilities of seismic exploration to study the main oil and gas reservoirs in this region. The solution was sought on a hexahedral mesh consisting of 5.5 mln spectral elements of the 5th order with a total number of SEM nodes 1.2 billion. The output results of full-wave modeling are stored in the SEG-Y format, suitable for all types of industrial seismic processing. The analysis of the obtained model seismograms and wave fields is carried out. The conclusion is made about the practical significance of the conducted research.

  Abstract
  The paper considers the solution of a three-dimensional problem of modeling of all types of seismic waves propagating in real geological media. The numerical algorithm based [...]

  • 0
  •  read

• The method of boundary integral equations in the boundary value problems of dynamics of fluids and gases

  Y. Krashanytsia

  WCCM2024.

  
  

  Abstract

  Based on the created generalized apparatus of vector-tensor analysis, integral representations of the main dynamic and kinematic characteristics of the problem of viscous gas flow around force systems of arbitrary spatial shape are constructed. The boundary value problem of the interaction of such systems with a viscous gas flow is reduced to a system of linear, conditioned by physical boundary conditions, boundary integral equations regarding the kinematic and dynamic characteristics of the problem. It is proven that all the obtained characteristics depend on the newly obtained irrotational vector potential of the momentum, which significantly simplifies the integral representations of solutions and their numerical implementation. On the basis of the created generalized apparatus of vector-tensor analysis, integral representations of the main dynamic and kinematic characteristics of the problem of the flow of a viscous gas flow around supporting systems of satisfactory spatial form have been constructed. The boundary value problem of the interaction of such systems with a viscous gas flow is reduced to a system of linear, conditioned by physical boundary conditions, boundary integral equations regarding the kinematic and dynamic characteristics of the problem. It is proven that all the obtained characteristics depend on the newly obtained, vortex-free vector potential of the momentum, which significantly simplifies the integral representations of the solutions and their numerical implementation.

  Abstract
  Based on the created generalized apparatus of vector-tensor analysis, integral representations of the main dynamic and kinematic characteristics of the problem of viscous [...]

  • 0
  •  read

• Landslide run-out simulations with depth-averaged models and integration with 3D impact analysis using the Material Point Method

  M. Fois, A. Makarim, C. FALCO, A. Larese, L. Formaggia

  WCCM2024.

  
  

  Abstract

  Landslides pose a significant threat to human safety and the well-being of communities, making them one of the most challenging natural phenomena. Their potential for catastrophic consequences, both in terms of human lives and economic impact, is a major concern. Additionally, their inherent unpredictability adds to the complexity of managing the risks associated with landslides. It is crucial to continuously monitor areas susceptible to landslides. In situ detection systems like piezometers and strain gauges play a vital role in accurately monitoring internal pressures and surface movements in the targeted areas. Simultaneously, satellite surveys contribute by offering detailed topographic and elevation data for the study area. However, relying solely on empirical monitoring is insufficient for ensuring effective management of hazardous situations, especially in terms of preventive measures. This study provides advanced simulations of mudflows and fast landslides using particle depth-averaged methods, specifically employing the Material Point Method adapted for shallow water (Depth Averaged Material Point Method). The numerical method has been parallelized and validated through benchmark tests and real-world cases. Furthermore, the investigation extends to coupling the depth-averaged formulation with a three-dimensional one in order to have a detailed description of the impact phase of the sliding material on barriers and membranes. The multidimensional approach and its validation on real cases provide a robust foundation for a more profound and accurate understanding of the behavior of mudflows and fast landslides

  Abstract
  Landslides pose a significant threat to human safety and the well-being of communities, making them one of the most challenging natural phenomena. Their potential for catastrophic [...]

  • 0
  •  read

• Numerical simulation investigation on nonlinear flow characteristics of rough single fractures with different contact areas

  Y. Liang, Y. Cheng, Z. Han, C. Yan

  WCCM2024.

  
  

  Abstract

  To comprehensively understand the influence of the contact area on the flow characteristics of rough single fractures, a rough fracture surface is initially constructed using a spatial frequency domain approach. Subsequently, rough single fractures with varying contact ratios are derived by translating and displacing the fracture surface. The Navier- Stokes equation and Mass-conservation equation are solved by utilizing the laminar flow module integrated within the COMSOL software. The simulation results show that the nonlinear correlation between fluid flow velocity and pressure gradient can be described by using Forchheimer equation. Under the same flow velocity, a higher contact rate will exacerbate the nonlinear characteristics of fluid flow. In contrast to non-contact fractures, the streamlines within contact fractures exhibit increased tortuosity, accompanied by an elongation of flow pathways. Furthermore, with an expanding contact area, the complexity of the streamline pattern amplifies. The overall pressure field distinctly exhibits non-uniform characteristics, with larger pressure gradient observed within localized contact regions, consequently facilitating an increase in flow velocity.

  Abstract
  To comprehensively understand the influence of the contact area on the flow characteristics of rough single fractures, a rough fracture surface is initially constructed using [...]

  • 0
  •  read

• Numerical study of fractal-tree-generated turbulence

  Y. Yin, R. Onishi, S. Watanabe, I. Segrovets

  WCCM2024.

  
  

  Abstract

  We study the aerodynamics of fractal trees by using a simulation based on the Lattice Boltzmann Method with a cumulant collision term. We have applied L-system rules to construct self-similar fractal tree models in aerodynamic computations. We found that the drag coefficient closely matches that of previous literature at high tree-height-based Reynolds numbers (ReH ≥ 60 000). A normalization process capable of collapsing turbulence intensity for various tree models is made. This process reveals that, at the same Reynolds number, different tree models exhibit the same behaviour in the turbulence intensity of their wake region. Our assessment of global and local isotropy in the turbulence generated by fractal trees reveals that the distant wake can be considered nearly locally isotropic at a high Reynolds number (ReH ≥ 60 000). Finally, the present numerical results confirm the non-equilibrium dissipation behaviour previously observed in the case of space-filling fractal square grids[2]. In the wake region, the non-dimensional dissipation rate Cϵ = constant is not valid. Instead, it is inversely proportional to the local Taylor-microscale-based Reynolds number, Cϵ ∝ 1/Reλ.

  Abstract
  We study the aerodynamics of fractal trees by using a simulation based on the Lattice Boltzmann Method with a cumulant collision term. We have applied L-system rules to construct [...]

  • 0
  •  read