Showing papers in "arXiv: Computational Engineering, Finance, and Science in 2019"

Direct shape optimization through deep reinforcement learning

Abstract: Deep Reinforcement Learning (DRL) has recently spread into a range of domains within physics and engineering, with multiple remarkable achievements. Still, much remains to be explored before the capabilities of these methods are well understood. In this paper, we present the first application of DRL to direct shape optimization. We show that, given adequate reward, an artificial neural network trained through DRL is able to generate optimal shapes on its own, without any prior knowledge and in a constrained time. While we choose here to apply this methodology to aerodynamics, the optimization process itself is agnostic to details of the use case, and thus our work paves the way to new generic shape optimization strategies both in fluid mechanics, and more generally in any domain where a relevant reward function can be defined.

Solving high-dimensional optimal stopping problems using deep learning

Abstract: Nowadays many financial derivatives, such as American or Bermudan options, are of early exercise type. Often the pricing of early exercise options gives rise to high-dimensional optimal stopping problems, since the dimension corresponds to the number of underlying assets. High-dimensional optimal stopping problems are, however, notoriously difficult to solve due to the well-known curse of dimensionality. In this work, we propose an algorithm for solving such problems, which is based on deep learning and computes, in the context of early exercise option pricing, both approximations of an optimal exercise strategy and the price of the considered option. The proposed algorithm can also be applied to optimal stopping problems that arise in other areas where the underlying stochastic process can be efficiently simulated. We present numerical results for a large number of example problems, which include the pricing of many high-dimensional American and Bermudan options, such as Bermudan max-call options in up to 5000 dimensions. Most of the obtained results are compared to reference values computed by exploiting the specific problem design or, where available, to reference values from the literature. These numerical results suggest that the proposed algorithm is highly effective in the case of many underlyings, in terms of both accuracy and speed.

A multi-resolution SPH method for fluid-structure interactions.

Abstract: In this paper, we present a multi-resolution smoothed particle hydrodynamics (SPH) method for modeling fluid-structure interaction (FSI) problems. By introducing different smoothing lengths and time steps, the spatio-temporal discretization is applied with different resolutions for fluid and structure. To ensure momentum conservation at the fluid-structure coupling, a position-based Verlet time integration scheme is introduced. Furthermore, the time-averaged velocity and acceleration of solid particles are introduced to enhance force matching in the fluid and solid equations. A set of numerical examples including several bio-mechanical problems are considered to demonstrate the efficiency, accuracy and robustness of the present method. A open-source code for all the examples is also provided.

Model-Free Data-Driven Methods in Mechanics: Material Data Identification and Solvers

Abstract: This paper presents an integrated model-free data-driven approach to solid mechanics, allowing to perform numerical simulations on structures on the basis of measures of displacement fields on representative samples, without postulating a specific constitutive model. A material data identification procedure, allowing to infer strain-stress pairs from displacement fields and boundary conditions, is used to build a material database from a set of mutiaxial tests on a non-conventional sample. This database is in turn used by a data-driven solver, based on an algorithm minimizing the distance between manifolds of compatible and balanced mechanical states and the given database, to predict the response of structures of the same material, with arbitrary geometry and boundary conditions. Examples illustrate this modelling cycle and demonstrate how the data-driven identification method allows importance sampling of the material state space, yielding faster convergence of simulation results with increasing database size, when compared to synthetic material databases with regular sampling patterns.

Determination of Metamaterial Parameters by Means of a Homogenization Approach Based on Asymptotic Analysis

Abstract: Owing to additive manufacturing techniques, a structure at millimeter length scale (macroscale) can be produced by using a lattice substructure at micrometer length scale (microscale). Such a system is called a metamaterial at the macroscale as the mechanical characteristics deviate from the characteristics at the microscale. As a remedy, metamaterial is modeled by using additional parameters; we intend to determine them. A homogenization approach based on the asymptotic analysis establishes a connection between these different characteristics at micro- and macroscales. A linear elastic first order theory at the microscale is related to a linear elastic second order theory at the macroscale. Relation for parameters at the macroscale is derived by using the equivalence of energy at macro- and microscales within a so-called Representative Volume Element (RVE). Determination of parameters are succeeded by solving a boundary value problem with the Finite Element Method (FEM). The proposed approach guarantees that the additional parameters vanish if the material is purely homogeneous, in other words, it is fully compatible with conventional homogenization schemes based on spatial averaging techniques. Moreover, the proposed approach is reliable as it ensures that such resolved additional parameters are not sensitive to choices of RVE consisting in the repetition of smaller RVEs but depend upon the intrinsic size of the structure.

Encoding Candlesticks as Images for Patterns Classification Using Convolutional Neural Networks

Abstract: Candlestick charts display the high, low, opening, and closing prices in a specific period. Candlestick patterns emerge because human actions and reactions are patterned and continuously replicate. These patterns capture information on the candles. According to Thomas Bulkowski's Encyclopedia of Candlestick Charts, there are 103 candlestick patterns. Traders use these patterns to determine when to enter and exit. Candlestick pattern classification approaches take the hard work out of visually identifying these patterns. To highlight its capabilities, we propose a two-steps approach to recognize candlestick patterns automatically. The first step uses the Gramian Angular Field (GAF) to encode the time series as different types of images. The second step uses the Convolutional Neural Network (CNN) with the GAF images to learn eight critical kinds of candlestick patterns. In this paper, we call the approach GAF-CNN. In the experiments, our method can identify the eight types of candlestick patterns with 90.7% average accuracy automatically in real-world data, outperforming the LSTM model.

Global cracking elements: a novel tool for Galerkin-based approaches simulating quasi-brittle fracture

Abstract: Following the so-called Cracking Elements Method (CEM), recently presented in \cite{Yiming:14,Yiming:16}, we propose a novel Galerkin-based numerical approach for simulating quasi-brittle fracture, named Global Cracking Elements Method (GCEM). For this purpose the formulation of the original CEM is reorganized. The new approach is embedded in the standard framework of the Galerkin-based Finite Element Method (FEM), which uses disconnected element-wise crack openings for capturing crack initiation and propagation. The similarity between the proposed Global Cracking Elements (GCE) and the standard 9-node quadrilateral element (Q9) suggests a special procedure: the degrees of freedom of the center node of the Q9, originally defining the displacements, are "borrowed" to describe the crack openings of the GCE. The proposed approach does not need remeshing, enrichment, or a crack-tracking strategy, and it avoids a precise description of the crack tip. Several benchmark tests provide evidence that the new approach inherits from the CEM most of the advantages. The numerical stability and robustness of the GCEM are better than the ones of the CEM. However, presently only quadrilateral elements with nonlinear interpolations of the displacement field can be used.

Multi-material Topology Optimization of Lattice Structures using Geometry Projection

Abstract: This work presents a computational method for the design of architected truss lattice materials where each strut can be made of one of a set of available materials. We design the lattices to extremize effective properties. As customary in topology optimization, we design a periodic unit cell of the lattice and obtain the effective properties via numerical homogenization. Each bar is represented as a cylindrical offset surface of a medial axis parameterized by the positions of the endpoints of the medial axis. These parameters are smoothly mapped onto a continuous density field for the primal and sensitivity analysis via the geometry projection method. A size variable per material is ascribed to each bar and penalized as in density-based topology optimization to facilitate the entire removal of bars from the design. During the optimization, we allow bars to be made of a mixture of the available materials. However, to ensure each bar is either exclusively made of one material or removed altogether from the optimal design, we impose optimization constraints that ensure each size variable is 0 or 1, and that at most one material size variable is 1. The proposed material interpolation scheme readily accommodates any number of materials. To obtain lattices with desired material symmetries, we design only a reference region of the unit cell and reflect its geometry projection with respect to the appropriate planes of symmetry. Also, to ensure bars remain whole upon reflection inside the unit cell or with respect to the periodic boundaries, we impose a no-cut constraint on the bars. We demonstrate the efficacy of our method via numerical examples of bulk and shear moduli maximization and Poisson's ratio minimization for two- and three-material lattices with cubic symmetry.

FLEXI: A high order discontinuous Galerkin framework for hyperbolic-parabolic conservation laws

Abstract: High order (HO) schemes are attractive candidates for the numerical solution of multiscale problems occurring in fluid dynamics and related disciplines. Among the HO discretization variants, discontinuous Galerkin schemes offer a collection of advantageous features which have lead to a strong increase in interest in them and related formulations in the last decade. The methods have matured sufficiently to be of practical use for a range of problems, for example in direct numerical and large eddy simulation of turbulence. However, in order to take full advantage of the potential benefits of these methods, all steps in the simulation chain must be designed and executed with HO in mind. Especially in this area, many commercially available closed-source solutions fall short. In this work, we therefor present the FLEXI framework, a HO consistent, open-source simulation tool chain for solving the compressible Navier-Stokes equations in a high performance computing setting. We describe the numerical algorithms and implementation details and give an overview of the features and capabilities of all parts of the framework. Beyond these technical details, we also discuss the important, but often overlooked issues of code stability, reproducibility and user-friendliness. The benefits gained by developing an open-source framework are discussed, with a particular focus on usability for the open-source community. We close with sample applications that demonstrate the wide range of use cases and the expandability of FLEXI and an overview of current and future developments.

Communication-Efficient Jaccard Similarity for High-Performance Distributed Genome Comparisons

Abstract: The Jaccard similarity index is an important measure of the overlap of two sets, widely used in machine learning, computational genomics, information retrieval, and many other areas. We design and implement SimilarityAtScale, the first communication-efficient distributed algorithm for computing the Jaccard similarity among pairs of large datasets. Our algorithm provides an efficient encoding of this problem into a multiplication of sparse matrices. Both the encoding and sparse matrix product are performed in a way that minimizes data movement in terms of communication and synchronization costs. We apply our algorithm to obtain similarity among all pairs of a set of large samples of genomes. This task is a key part of modern metagenomics analysis and an evergrowing need due to the increasing availability of high-throughput DNA sequencing data. The resulting scheme is the first to enable accurate Jaccard distance derivations for massive datasets, using largescale distributed-memory systems. We package our routines in a tool, called GenomeAtScale, that combines the proposed algorithm with tools for processing input sequences. Our evaluation on real data illustrates that one can use GenomeAtScale to effectively employ tens of thousands of processors to reach new frontiers in large-scale genomic and metagenomic analysis. While GenomeAtScale can be used to foster DNA research, the more general underlying SimilarityAtScale algorithm may be used for high-performance distributed similarity computations in other data analytics application domains.

Data-driven discovery of free-form governing differential equations.

Abstract: We present a method of discovering governing differential equations from data without the need to specify a priori the terms to appear in the equation. The input to our method is a dataset (or ensemble of datasets) corresponding to a particular solution (or ensemble of particular solutions) of a differential equation. The output is a human-readable differential equation with parameters calibrated to the individual particular solutions provided. The key to our method is to learn differentiable models of the data that subsequently serve as inputs to a genetic programming algorithm in which graphs specify computation over arbitrary compositions of functions, parameters, and (potentially differential) operators on functions. Differential operators are composed and evaluated using recursive application of automatic differentiation, allowing our algorithm to explore arbitrary compositions of operators without the need for human intervention. We also demonstrate an active learning process to identify and remedy deficiencies in the proposed governing equations.

Deformable Models for Surgical Simulation: A Survey

Abstract: This paper presents a survey of the state-of-the-art deformable models studied in the literature concerning soft tissue deformable modeling for interactive surgical simulation. It first introduces the challenges of surgical simulation, followed by discussions and analyses on the deformable models, which are classified into three categories: the heuristic modeling methodology, continuum-mechanical methodology, and other methodologies. It also examines linear and nonlinear deformable modeling, model internal forces, and numerical time integrations, together with modeling of soft tissue anisotropy, viscoelasticity, and compressibility. Finally, various issues in the existing deformable models are discussed to outline the remaining challenges of deformable models in surgical simulation.

Three-dimensional topology optimization of auxetic metamaterial using isogeometric analysis and model order reduction

Abstract: In this work, we present an efficiently computational approach for designing material micro-structures by means of topology optimization. The central idea relies on using the isogeometric analysis integrated with the parameterized level set function for numerical homogenization, sensitivity calculation and optimization of the effective elastic properties. Design variables, which are level set values associated with control points, are updated from the optimizer and represent the geometry of the unit cell. We further improve the computational efficiency in each iteration by employing reduced order modeling when solving linear systems of the equilibrium equations. We construct a reduced basis by reusing computed solutions from previous optimization steps, and a much smaller linear system of equations is solved on the reduced basis. Two- and three-dimensional numerical results show the effectiveness of the topology optimization algorithm coupled with the reduced basis approach in designing metamaterials.

Accurate prediction of melt pool shapes in laser powder bed fusion by the non-linear temperature equation including phase changes - isotropic versus anisotropic conductivity.

Abstract: In this contribution, we validate a physical model based on a transient temperature equation (including latent heat) w.r.t. the experimental set AMB2018-02 provided within the additive manufacturing benchmark series, established at the National Institute of Standards and Technology, USA. We aim at predicting the following quantities of interest: width, depth, and length of the melt pool by numerical simulation and report also on the obtainable numerical results of the cooling rate. We first assume the laser to posses a double ellipsoidal shape and demonstrate that a well calibrated, purely thermal model based on isotropic thermal conductivity is able to predict all the quantities of interest, up to a deviation of maximum 7.3\% from the experimentally measured values. However, it is interesting to observe that if we directly introduce, whenever available, the measured laser profile in the model (instead of the double ellipsoidal shape) the investigated model returns a deviation of 19.3\% from the experimental values. This motivates a model update by introducing anisotropic conductivity, which is intended to be a simplistic model for heat material convection inside the melt pool. Such an anisotropic model enables the prediction of all quantities of interest mentioned above with a maximum deviation from the experimental values of 6.5\%. We note that, although more predictive, the anisotropic model induces only a marginal increase in computational complexity.

A mortar-type finite element approach for embedding 1D beams into 3D solid volumes

Abstract: In this work we present a novel computational method for embedding arbitrary curved one-dimensional (1D) fibers into three-dimensional (3D) solid volumes, as e.g. in fiber-reinforced materials. The fibers are explicitly modeled with highly efficient 1D geometrically exact beam finite elements, based on various types of geometrically nonlinear beam theories. The surrounding solid volume is modeled with 3D continuum (solid) elements. An embedded mortar-type approach is employed to enforce the kinematic coupling constraints between the beam elements and solid elements on non-matching meshes. This allows for very flexible mesh generation and simple material modeling procedures in the solid, since it can be discretized without having to capture for the reinforcements, while still being able to account for complex nonlinear effects due to the embedded fibers. Several numerical examples demonstrate the consistency, robustness and accuracy of the proposed method, as well as its applicability to rather complex fiber-reinforced structures of practical relevance.

A Total Lagrangian SPH Method for Modelling Damage and Failure in Solids

Abstract: An algorithm is proposed to model crack initiation and propagation within the total Lagrangian Smoothed Particle Hydrodynamics (TLSPH) framework. TLSPH avoids the two main deficiencies of conventional SPH, i.e., tensile instability and inconsistency, by making use of the Lagrangian kernel and gradient correction, respectively. In the present approach, the support domain of a particle is modified, where it only interacts with its immediately neighbouring particles. A virtual link defines the level of interaction between each particle pair. The state of the virtual link is determined by damage law or cracking criterion. The virtual link approach allows easy and natural modelling of cracking surfaces without explicit cracking treatments such as particle splitting, field enrichment or visibility criterion. The performance of the proposed approach is demonstrated via a few numerical examples of both brittle and ductile failure under impact loading.

Fleet Prognosis with Physics-informed Recurrent Neural Networks

Abstract: Services and warranties of large fleets of engineering assets is a very profitable business. The success of companies in that area is often related to predictive maintenance driven by advanced analytics. Therefore, accurate modeling, as a way to understand how the complex interactions between operating conditions and component capability define useful life, is key for services profitability. Unfortunately, building prognosis models for large fleets is a daunting task as factors such as duty cycle variation, harsh environments, inadequate maintenance, and problems with mass production can lead to large discrepancies between designed and observed useful lives. This paper introduces a novel physics-informed neural network approach to prognosis by extending recurrent neural networks to cumulative damage models. We propose a new recurrent neural network cell designed to merge physics-informed and data-driven layers. With that, engineers and scientists have the chance to use physics-informed layers to model parts that are well understood (e.g., fatigue crack growth) and use data-driven layers to model parts that are poorly characterized (e.g., internal loads). A simple numerical experiment is used to present the main features of the proposed physics-informed recurrent neural network for damage accumulation. The test problem consist of predicting fatigue crack length for a synthetic fleet of airplanes subject to different mission mixes. The model is trained using full observation inputs (far-field loads) and very limited observation of outputs (crack length at inspection for only a portion of the fleet). The results demonstrate that our proposed hybrid physics-informed recurrent neural network is able to accurately model fatigue crack growth even when the observed distribution of crack length does not match with the (unobservable) fleet distribution.

Short Overview of Early Developments of the Hardy Cross Type Methods for Computation of Flow Distribution in Pipe Networks

Abstract: Hardy Cross originally proposed a method for analysis of flow in networks of conduits or conductors in 1936. His method was the first really useful engineering method in the field of pipe network calculation. Only electrical analogs of hydraulic networks were used before the Hardy Cross method. A problem with flow resistance versus electrical resistance makes these electrical analog methods obsolete. The method by Hardy Cross is taught extensively at faculties, and it remains an important tool for the analysis of looped pipe systems. Engineers today mostly use a modified Hardy Cross method which considers the whole looped network of pipes simultaneously (use of these methods without computers is practically impossible). A method from a Russian practice published during the 1930s, which is similar to the Hardy Cross method, is described, too. Some notes from the work of Hardy Cross are also presented. Finally, an improved version of the Hardy Cross method, which significantly reduces the number of iterations, is presented and discussed. We also tested multi-point iterative methods, which can be used as a substitution for the Newton-Raphson approach used by Hardy Cross, but in this case this approach did not reduce the number of iterations. Although many new models have been developed since the time of Hardy Cross, the main purpose of this paper is to illustrate the very beginning of modelling of gas and water pipe networks and ventilation systems. As a novelty, a new multi-point iterative solver is introduced and compared with the standard Newton-Raphson iterative method.

An adaptive space-time phase field formulation for dynamic fracture of brittle shells based on LR NURBS

Abstract: We present an adaptive space-time phase field formulation for dynamic fracture of brittle shells. Their deformation is characterized by the Kirchhoff-Love thin shell theory using a curvilinear surface description. All kinematical objects are defined on the shell's mid-plane. The evolution equation for the phase field is determined by the minimization of an energy functional based on Griffith's theory of brittle fracture. Membrane and bending contributions to the fracture process are modeled separately and a thickness integration is established for the latter. The coupled system consists of two nonlinear fourth-order PDEs and all quantities are defined on an evolving two-dimensional manifold. Since the weak form requires $C^1$-continuity, isogeometric shape functions are used. The mesh is adaptively refined based on the phase field using Locally Refinable (LR) NURBS. Time is discretized based on a generalized-$\alpha$ method using adaptive time-stepping, and the discretized coupled system is solved with a monolithic Newton-Raphson scheme. The interaction between surface deformation and crack evolution is demonstrated by several numerical examples showing dynamic crack propagation and branching.

Incorporating Fine-grained Events in Stock Movement Prediction

Abstract: Considering event structure information has proven helpful in text-based stock movement prediction. However, existing works mainly adopt the coarse-grained events, which loses the specific semantic information of diverse event types. In this work, we propose to incorporate the fine-grained events in stock movement prediction. Firstly, we propose a professional finance event dictionary built by domain experts and use it to extract fine-grained events automatically from finance news. Then we design a neural model to combine finance news with fine-grained event structure and stock trade data to predict the stock movement. Besides, in order to improve the generalizability of the proposed method, we design an advanced model that uses the extracted fine-grained events as the distant supervised label to train a multi-task framework of event extraction and stock prediction. The experimental results show that our method outperforms all the baselines and has good generalizability.

Topology Optimization of Fluidic Pressure Loaded Structures and Compliant Mechanisms using the Darcy Method

Abstract: In various applications, design problems involving structures and compliant mechanisms experience fluidic pressure loads. During topology optimization of such design problems, these loads adapt their direction and location with the evolution of the design, which poses various challenges. A new density-based topology optimization approach using Darcy's law in conjunction with a drainage term is presented to provide a continuous and consistent treatment of design-dependent fluidic pressure loads. The porosity of each finite element and its drainage term are related to its density variable using a Heaviside function, yielding a smooth transition between the solid and void phases. A design-dependent pressure field is established using Darcy's law and the associated PDE is solved using the finite element method. Further, the obtained pressure field is used to determine the consistent nodal loads. The approach provides a computationally inexpensive evaluation of load sensitivities using the adjoint-variable method. To show the efficacy and robustness of the proposed method, numerical examples related to fluidic pressure loaded stiff structures and small-deformation compliant mechanisms are solved. For the structures, compliance is minimized, whereas for the mechanisms a multi-criteria objective is minimized with given resource constraints.

A state-of-knowledge review on the Endurance Time Method

Abstract: Endurance time method is a time history dynamic analysis in which structures are subjected to predesigned intensifying excitations. This method provides a tool for response prediction that correlates structural responses to the intensity of earthquakes with a considerably less computational demand as compared to conventional time history analysis. The endurance time method is being used in different areas of earthquake engineering such as performance-based assessment and design, life-cycle cost-based design, value-based design, seismic safety, seismic assessment, and multicomponent seismic analysis. Successful implementation of the endurance time method relies heavily on the quality of endurance time excitations. In this paper, a review of the endurance time method from conceptual development to its practical applications is provided. Different types of endurance time excitations are described. Features related to the existing endurance time excitations are also presented. Particular attention is given to different applications of the endurance time method in the field of earthquake engineering.

Diversification Across Mining Pools: Optimal Mining Strategies under PoW.

Abstract: Mining is a central operation of all proof-of-work (PoW) based cryptocurrencies. The vast majority of miners today participate in "mining pools" instead of "solo mining" in order to lower risk and achieve a more steady income. However, this rise of participation in mining pools negatively affects the decentralization levels of most cryptocurrencies. In this work, we look into mining pools from the point of view of a miner: We present an analytical model and implement a computational tool that allows miners to optimally distribute their computational power over multiple pools and PoW cryptocurrencies (i.e. build a mining portfolio), taking into account their risk aversion levels. Our tool allows miners to maximize their risk-adjusted earnings by diversifying across multiple mining pools which enhances PoW decentralization. Finally, we run an experiment in Bitcoin historical data and demonstrate that a miner diversifying over multiple pools, as instructed by our model/tool, receives a higher overall Sharpe ratio (i.e. average excess reward over its standard deviation/volatility).

Accelerated Steady-State Torque Computation for Induction Machines using Parallel-In-Time Algorithms

Abstract: This paper focuses on efficient steady-state computations of induction machines. In particular, the periodic Parareal algorithm with initial-value coarse problem (PP-IC) is considered for acceleration of classical time-stepping simulations via non-intrusive parallelization in time domain, i.e., existing implementations can be reused. Superiority of this parallel-in-time method is in its direct applicability to time-periodic problems, compared to, e.g, the standard Parareal method, which only solves an initial-value problem, starting from a prescribed initial value. PP-IC is exploited here to obtain the steady state of several operating points of an induction motor, developed by Robert Bosch GmbH. Numerical experiments show that acceleration up to several dozens of times can be obtained, depending on availability of parallel processing units. Comparison of PP-IC with existing time-periodic explicit error correction method highlights better robustness and efficiency of the considered time-parallel approach.

BRTSim, a general-purpose computational solver for hydrological, biogeochemical, and ecosystem dynamics.

Abstract: This paper introduces the recent release v3.1a of BRTSim (BioReactive Transport Simulator), a general-purpose multiphase and multi-species liquid, gas and heat flow solver for reaction-advection-dispersion processes in porous and non-porous media with application in hydrology and biogeochemistry. Within the philosophy of the BRTSim platform, the user can define (1) arbitrary chemical and biological species; (2) arbitrary chemical and biological reactions; (3) arbitrary equilibrium reactions; and (4) combine solvers for phases and heat flows as well as for specialized biological processes such as bioclogging and chemotaxis. These capabilities complement a suite of processes and process-feedback solvers not currently available in other general-purpose codes. Along with the flexibility to design arbitrarily complex reaction networks, and setup and synchronize solvers through one input text file, BRTSim can communicate with third-party software with ease. Here, four cases study that combine experimental observations and modeling with BRTSim are reported: (i) water table dynamics in a heterogeneous aquifer for variable hydrometeorological conditions; (ii) soil biological clogging by cells and exopolymers; (iii) biotic degradation and isotopic fractionation of nitrate; and (iv) dispersion and biodegradation of atrazine herbicide in agricultural crops.

Phase-field modeling of fracture for quasi-brittle materials

Abstract: This paper addresses the modeling of fracture in quasi-brittle materials using a phase-field approach to the description of crack topology. Within the computational mechanics community, several studies have treated the issue of modeling fracture using phase fields. Most of these studies have used an approach that implies the lack of a damage threshold. We herein explore an alternative model that includes a damage threshold and study how it compares with the most popular approach. The formulation is systematically explained within a rigorous variational framework. Subsequently, we present the corresponding three-dimensional finite element discretization that leads to a straightforward numerical implementation. Benchmark simulations in two dimensions and three dimensions are then presented. The results show that while an elastic stage and a damage threshold are ensured by the present model, good agreement with the results reported in the literature can be obtained, where such features are generally absent.

Numerical Modelling of Crack Initiation, Propagation and Branching under Dynamic Loading

Abstract: In this paper crack initiation, propagation and branching phenomena are simulated using the Pseudo-Spring Smooth Particle Hydrodynamics (SPH) in two and three-dimensional domains. The pseudo-spring analogy is used to model material damage. Here, the interaction of particles is limited to its initial immediate neighbours. The particles are connected via springs. These springs do not provide any extra stiffness in the system but only define the level of interaction between the connecting pairs. It is assumed that a crack has passed through a spring connecting a particle pair if the damage indicator of that spring becomes more than a predefined value. The crack branching of a pre-notched plate under dynamic loading and the effect of loading amplitude are studied here. The computed crack speeds, crack paths and surfaces are compared with experimental and numerical results available in the literature and are found to be in good agreement. The effect of notch location for a plate with a circular hole is studied here. The ability of the framework to model arbitrary crack paths and surfaces are demonstrated via three-dimensional simulations of chalk under torsion, Kalthoff-Winkler experiment and Taylor bullet impact.

A hybrid material-point spheropolygon-element method for solid and granular material interaction

Abstract: Capturing the interaction between objects that have an extreme difference in Young s modulus or geometrical scale is a highly challenging topic for numerical simulation. One of the fundamental questions is how to build an accurate multi-scale method with optimal computational efficiency. In this work, we develop a material-point-spheropolygon discrete element method (MPM-SDEM). Our approach fully couples the material point method (MPM) and the spheropolygon discrete element method (SDEM) through the exchange of contact force information. It combines the advantage of MPM for accurately simulating elastoplastic continuum materials and the high efficiency of DEM for calculating the Newtonian dynamics of discrete near-rigid objects. The MPM-SDEM framework is demonstrated with an explicit time integration scheme. Its accuracy and efficiency are further analysed against the analytical and experimental data. Results demonstrate this method could accurately capture the contact force and momentum exchange between materials while maintaining favourable computational stability and efficiency. Our framework exhibits great potential in the analysis of multi-scale, multi-physics phenomena

An efficient iterative method for looped pipe network hydraulics.

Abstract: Original and improved version of the Hardy Cross iterative method with related modifications are today widely used for calculation of fluid flow through conduits in loops-like distribution networks of pipes with known node fluid consumptions. Fluid in these networks is usually natural gas for distribution in the municipalities, water in waterworks or hot water in district heating system, air in the case of ventilation systems in buildings or mines, etc. Since, the resistances in these networks depend of flow, problem is not linear like in electrical circuits, and iterative procedure must be used. In both version of the Hardy Cross method, in original and in the improved one, initial results of calculation in iteration procedure is not flow, but rather the correction of flow. Unfortunately, these corrections should be added to or subtracted from a flow calculated in previous iteration according to complicate algebraic rules. After the here presented node-loop method, final results in each of the iterations is flow directly rather than flow correction. In that way complex algebraic scheme for sign of flow correction is avoided, while the final results still remain unchanged. Numbers of required iterations for the same results are comparable with the improved Hardy Cross method.

A Fractional Subgrid-scale Model for Turbulent Flows: Theoretical Formulation and a Priori Study

Abstract: Coherent structures/motions in turbulence inherently give rise to intermittent signals with sharp peaks, heavy-skirt, and skewed distributions of velocity increments, highlighting the non-Gaussian nature of turbulence. That suggests that the spatial nonlocal interactions cannot be ruled out of the turbulence physics. Furthermore, filtering the Navier-Stokes equations in the large eddy simulation of turbulent flows would further enhance the existing nonlocality, emerging in the corresponding subgrid scale fluid motions. That urges the development of new nonlocal closure models, which respect the corresponding non-Gaussian statistics of the subgrid stochastic motions. To this end and starting from the filtered Boltzmann equation, we model the corresponding equilibrium distribution function with a \textit{Levy}-stable distribution, leading to the proposed fractional-order modeling of subgrid-scale stresses. We approximate the filtered equilibrium distribution function with a power-law term, and derive the corresponding filtered Navier-Stokes equations. Subsequently in our functional modeling, the divergence of subgrid-scale stresses emerges as a single-parameter fractional Laplacian, $(-\Delta)^{\alpha}(\cdot)$, $\alpha \in (0,1]$, of the filtered velocity field. The only model parameter, i.e., the fractional exponent, appears to be strictly depending on the filter-width and the flow Reynolds number. We furthermore explore the main physical and mathematical properties of the proposed model under a set of mild conditions. Finally, the introduced model undergoes \textit{a priori} evaluations based on the direct numerical simulation database of forced and decaying homogeneous isotropic turbulent flows at relatively high and moderate Reynolds numbers, respectively. Such analysis provides a comparative study of predictability and performance of the proposed fractional model.