Search results

Accurate presentation of the rock microstructure is critical to the grain-scale analysis of rock deformation and failure in numerical modelling. 3D granite microstructure modelling has only been used in limited studies with the mineral pattern often remaining poorly constructed. In this study, the authors developed a new approach for generating 2D and 3D granite microstructure models from a 2D image by combining a heterogeneous material reconstruction method (simulated annealing method) with Voronoi tessellation.

More specifically, the stochastic information in the 2D image is first extracted using the two-point correlation function (TPCF). Then an initial 2D or 3D Voronoi diagram with a random distribution of the minerals is generated and optimised using a simulated annealing method until the corresponding TPCF is consistent with that in the 2D image. The generated microstructure model accurately inherits the stochastic information (e.g. volume fraction and mineral pattern) from the 2D image. Lastly, the authors compared the topological characteristics and mechanical properties of the 2D and 3D reconstructed microstructure models with the model obtained by direct mapping from the 2D image of a real rock sample.

The good agreements between the mapped and reconstructed models indicate the accuracy of the reconstructed microstructure models on topological characteristics and mechanical properties.

The newly developed reconstruction method successfully transfers the mineral pattern from a granite sample into the 2D and 3D Voronoi-based microstructure models ready for use in grain-scale modelling.

High head gates are commonly used in hydropower plants for flow regulation and emergence closure. Hydrodynamic downpull can be a critical parameter in design of the lifting mechanism. The purpose of this paper is to show that a simplified two-dimensional (2D) computational fluid dynamics solution can be used in the prediction of the downpull force on the gate lip by comparison of computed results to experimentally measured data.

In this study, ANSYS FLUENT CFD software was used to obtain 2D numerical solution for the flow field around a generic gate model located in a power intake structure which was previously used in an experimental study. Description of the flow domain, computational grid resolution, requirements on setting appropriate boundary conditions and methodology in describing downpull coefficient are discussed. Total number of 245 simulations for variable gate lip geometry and gate openings were run. The downpull coefficient evaluated from the computed pressure field as function of gate opening and lip angle are compared with the experimental results.

The computed downpull coefficient agrees well with the previous experimental results, except one gate with small lip angle where a separation bubble forms along the lip, which is responsible from this deviation. It is observed that three-dimensional (3D) effects are confined to the large gate openings where downpull is minimum or even reversed.

In large gate openings, three dimensionality of the flow around gate slots plays an important role and departure from 2D solutions become more pronounced. In that case, one might need to perform a 3D solution instead.

This paper presents a very fast and accurate way to predict downpull force on high head gates in the absence of experimental data.

An extensive amount of simulations are run within the scope of this study. It is shown that knowing its limitations, 2D numerical models can be used to calculate downpull for a wide range of gate openings without the need of expensive experimental models.

The purpose of this paper is to investigate the applicability of the LS-DYNA software using a Lagrangian formulation in the jet formation, flight and penetration of improvised explosively formed projectiles (EFPs). Numerical results dealing with different properties of the EFPs have been validated with a significant number of field tests.

2D and 3D Lagrangian models, using different material definition, are developed to reproduce the field-measured characteristics of copper- and steel-made EFPs: projectile size and velocity. After validation, the model has been extended to analyse the penetration features. Two different plasticity models have been used to describe the steel target, Plastic-Kinematic and Johnson–Cook.

Despite the difficulty in characterizing a non-industrial artefact, the results show that both Lagrangian models (2D and 3D) are able to simulate the projectile size, velocity and penetration capability with errors less than 10 per cent when using the Johnson–Cook material model for both liner and target.

These data can be used to test the penetration ability of improvised EFP’s against different targets, i.e. light armoured vehicles.

There are no references that address the application of the Lagrangian simulation of non-industrial EFPs and its validation with field tests, including penetration assessment.

This paper aims to investigate the mechanisms lying behind the cycloidal rotor under hovering status.

Experiments were conducted to validate the numerical simulation results. The simulations were based on unsteady Reynolds-averaged Navier–Stokes (URANS) equations solver and the sliding mesh technique was used to model the blade motion. 2D and 2.5D simulations were made to investigate the 3D effects of turbulence. The effects of pressure and viscosity were compared to study the significance of the blade motion on force generation.

The 2.5D numerical simulation cannot produce more accurate results than the 2D counterpart. The pitching motion of the blade results in dynamic stall. The dynamic stall vortices induce parallel blade vortex interaction (BVI) upon downstream blades. The interactions between the blades delay the stall of the blade which is beneficial to the thrust generation. The blade pitching motion is the dominant contributor to the force generation and the turbulence is the secondary. Strong downwash in the rotor cage varied the inflow velocity as well as the effective angle of attack (AOA) of the blade.

Cycloidal rotor is a propulsion device that can provide omni-directional vectored thrust with high efficiency and low noise. To understand the mechanisms lying behind the cycloidal rotor helps the authors to design efficient cycloidal rotors for aircraft.

The authors discovered that the blade pitching motion plays primary role in force generation. The effects of the dynamic stall and BVI were studied. The reason why cycloidal rotor can be more efficient was discussed.

A program for numerical simulation of two‐dimensional semiconductor devices coupled with an external circuit is described. The circuit equations are formulated using modified nodal analysis to allow an arbitrary configuration of elements like, e.g., also semiconductor compact models. Coupling to the numerical devices is attained via their admittance matrix leading to a two‐level Newton method. To calculate this matrix two methods are compared: (a) a linearization scheme and (b) a secant method. The comparison shows a significant speed advantage of the secant method despite its lower rate of convergence. The linearization scheme, however, is the more stable and robust method and should be used in critical cases where convergence problems can occur. An efficient bypassing scheme was developed for the linearization scheme leading to a computation speed comparable to that of the secant method, but maintaining the better convergence properties. A further advantage of the two‐level Newton method used in this work is that the CPU‐time consuming solution for the numerical devices can be done in parallel on different processors. Several examples are given to demonstrate the capabilities of the developed simulator.

This paper aims to describe a proposal for an innovative method of normal shock wave–turbulent boundary layer interaction (SBLI) and shock-induced separation control.

The concept is based on the introduction of a tangentially moving wall upstream of the shock wave and in the interaction region. The SBLI control mechanism may be implemented as a closed belt floating on an air cushion, sliding over two cylinders and forming the outer skin of the suction side of the airfoil. The presented exploratory numerical study is conducted with SPARC solver (steady 2D RANS). The effect of the moving wall is presented for the NACA 0012 airfoil operating in transonic conditions.

To assess the accuracy of obtained solutions, validation of the computational model is demonstrated against the experimental data of Harris, Ladson & Hill and Mineck & Hartwich (NASA Langley). The comparison is conducted not only for the reference (impermeable) but also for the perforated (permeable) surface NACA 0012 airfoils. Subsequent numerical analysis of SBLI control by moving wall confirms that for the selected velocity ratios, the method is able to improve the shock-upstream boundary layer and counteract flow separation, significantly increasing the airfoil aerodynamic performance.

The moving wall concept as a means of normal shock wave–turbulent boundary layer interaction and shock-induced separation control has been investigated in detail for the first time. The study quantified the necessary operational requirements of such a system and practicable aerodynamic efficiency gains and simultaneously revealed the considerable potential of this promising idea, stimulating a new direction for future investigations regarding SBLI control.

The unstable dynamic propagation of multistage hydrofracturing fractures leads to uneven development of the fracture network and research on the mechanism controlling this phenomenon indicates that the stress shadow effects around the fractures are the main mechanism causing this behaviour. Further studies and simulations of the stress shadow effects are necessary to understand the controlling mechanism and evaluate the fracturing effect.

In the process of stress-dependent unstable dynamic propagation of fractures, there are both continuous stress fields and discontinuous fractures; therefore, in order to study the stress-dependent unstable dynamic propagation of multistage fracture networks, a series of continuum-discontinuum numerical methods and models are reviewed, including the well-developed extended finite element method, displacement discontinuity method, boundary element method and finite element-discrete element method.

The superposition of the surrounding stress field during fracture propagation causes different degrees of stress shadow effects between fractures and the main controlling factors of stress shadow effects are fracture initiation sequence, perforation cluster spacing and well spacing. The perforation cluster spacing varies with the initiation sequence, resulting in different stress shadow effects between fractures; for example, the smaller the perforation cluster spacing and well spacing are, the stronger the stress shadow effects are and the more seriously the fracture propagation inhibition arises. Moreover, as the spacing of perforation clusters and well spacing increases, the stress shadow effects decrease and the fracture propagation follows an almost straight pattern. In addition, the computed results of the dynamic distribution of stress-dependent unstable dynamic propagation of fractures under different stress fields are summarised.

A state-of-art review of stress shadow effects and continuum-discontinuum methods for stress-dependent unstable dynamic propagation of multiple hydraulic fractures are well summarized and analysed. This paper can provide a reference for those engaged in the research of unstable dynamic propagation of multiple hydraulic structures and have a comprehensive grasp of the research in this field.

Aims to present recent activities in numerical modeling of cold crucible melting process.

3D numerical analysis was used for electromagnetic problem and 3D large eddy simulation (LES) method was applied for fluid flow modeling.

The comparative modeling shows, that higher H/D ratio of the melt is more efficient when total power consumption is considered, but this advantage is held back by higher heat losses through the crucible walls. Also, calculations reveal that lower frequencies, which are energetically less effective, provide better mixing of the melt.

3D electromagnetic model, which allows to take into account non‐symmetrical distribution of Joule heat sources, together with transient LES fluid flow simulation gives the opportunity of accurate prediction of temperature distribution in the melt.

A nonlinear finite element (FE) beam-column model for the analysis of reinforced concrete (RC) frames with due account of shear is presented in this paper. The model is an expansion of the traditional flexural fibre beam formulations to cases where multiaxial behaviour exists, being an alternative to plane and solid FE models for the nonlinear analysis of entire frame structures. The paper aims to discuss these issues.

Shear is taken into account at different levels of the numerical model: at the material level RC is simulated through a smeared cracked approach with rotating cracks; at the fibre level, an iterative procedure guarantees equilibrium between concrete and transversal reinforcement, allowing to compute the biaxial stress-strain state of each fibre; at the section level, a uniform shear stress pattern is assumed in order to estimate the internal shear stress-strain distribution; and at the element level, the Timoshenko beam theory takes into account an average rotation due to shear.

The proposed model is validated through experimental tests available in the literature, as well as through an experimental campaign carried out by the authors. The results on the response of RC elements critical to shear include displacements, strains and crack patterns and show the capabilities of the model to efficiently deal with shear effects in beam elements.

A formulation for the nonlinear shear-bending interaction based on the fixed stress approach is implemented in a fibre beam model. Shear effects are accurately accounted during all the nonlinear path of the structure in a computationally efficient manner.

This paper aims to develop a numerical model of bead spreading architecture of a viscous polymer in fused filament fabrication (FFF) process with different nozzle geometry. This paper also focuses on the manufacturing feasibility of the nozzles and 3D printing of the molten beads using the developed nozzles.

The flow of a highly viscous polymer from a nozzle, the melt expansion in free space and the deposition of the melt on a moving platform are captured using the FLUENT volume of fluid (VOF) method based computational fluid dynamics code. The free surface motion of the material is captured in VOF, which is governed by the hydrodynamics of the two-phase flow. The phases involved in the numerical model are liquid polymer and air. A laminar, non-Newtonian and non-isothermal flow is assumed. Under such assumptions, the spreading characteristic of the polymer is simulated with different nozzle-exit geometries. The governing equations are solved on a regular stationary grid following a transient algorithm, where the boundary between the polymer and the air is tracked by piecewise linear interface construction (PLIC) to reconstruct the free surface. The prototype nozzles were also manufactured, and the deposition of the molten beads on a flatbed was performed using a commercial 3D printer. The deposited bead cross-sections were examined through optical microscopic examination, and the cross-sectional profiles were compared with those obtained in the numerical simulations.

The numerical model successfully predicted the spreading characteristics and the cross-sectional shape of the extruded bead. The cross-sectional shape of the bead varied from elliptical (with circular nozzle) to trapezoidal (with square and star nozzles) where the top and bottom surfaces are significantly flattened (which is desirable to reduce the void spaces in the cross-section). The numerical model yielded a good approximation of the bead cross-section, capturing most of the geometric features of the bead with a reasonable qualitative agreement compared to the experiment. The quantitative comparison of the cross-sectional profiles against experimental observation also indicated a favorable agreement. The significant improvement observed in the bead cross-section with the square and star nozzles is the flattening of the surfaces.

The developed numerical algorithm attempts to address the fundamental challenge of voids and bonding in the FFF process. It presents a new approach to increase the inter-bead bonding and reduce the inter-bead voids in 3D printing of polymers by modifying the bead cross-sectional shape through the modification of nozzle exit-geometry. The change in bead cross-sectional shape from elliptical (circular) to trapezoidal (square and star) cross-section is supposed to increase the contact surface area and inter-bead bonding while in contact with adjacent beads.