Plenary Lecture

Prof. Joannes J. Westerink (University of Notre Dame, USA)

Towards Heterogeneous Process, Scale, and Model Coupling in Simulating the Hydrodynamics of the Coastal Ocean

Computational models for wind waves and storm driven surge and currents in the coastal ocean and adjacent floodplain must provide a high level of grid resolution, fully couple the wind wave and long wave processes, and perform quickly for risk assessment, flood mitigation on system design, and forecasting purposes. We have developed a high performance unstructured grid computational framework that couples circulation and non-phase resolving wave models while scaling efficiently up to 32,000 cores. Current development is focused on incorporating a wider range of physics affecting coastal and inland water levels as well as forces on infrastructure including large scale baroclinically driven processes, rainfall runoff in upland areas and on the coastal floodplain, and wave run-up.  This is accomplished with an interleafing framework in which heterogeneous models focused on a select range of processes are coupled over the same domain and/or specific targeted equations that are dynamically assigned to changing portions of the domain as appropriate to the prevailing flow conditions.  This is all done in a dynamically load balanced framework.  Algorithmic development is focused on DG solvers, ideally suited for the associated strongly advective flows, allow super-parametric elements for p=1 and p=2 and iso-parametric elements for p=3 in order to achieve improved convergence rates and overall run me efficiency, and allow for the selection of localized physics on the elemental level.

 

Prof. Manolis Papadrakakis (National Technical University of Athens, Greece)

Seismic Assessment of Reinforced Concrete Structures based on State-of-the-art 3D Detailed Nonlinear Finite Element Simulations

The nonlinear dynamic numerical simulation of reinforced concrete structures is characterized by instabilities, which are mainly caused by the cracking of concrete and the rapture of steel reinforcement. When dealing with this numerically unstable and computationally demanding problem, the numerical solution procedure becomes extremely cumbersome, thus leading to convergence issues and the inability to capture the ultimate bearing capacity of the structure. Additionally, the lack of objectivity, when using 1D and 2D models, does not allow the study of the nonlinear dynamic response of R/C structures without introducing significant simplified assumptions in-terms of material behavior and the exact discretization of the structural geometry.

In light of these well-known modeling limitations, the main objective of this research work is to alleviate the above-mentioned numerical constraints, by developing a state-of-the-art 3D detailed modeling approach that will provide the computational tools to perform dynamic nonlinear analysis and design on large-scale reinforced concrete structures, by accounting for soil-foundation-structure interaction phenomena as well. In order to achieve this objective, the numerical handling of the solution instabilities is addressed herein, while the use of the HYMOD surrogate modeling approach is discussed as a potential more practical solution to the overall modeling problem. Furthermore, the results from a developed parallel solver for both the generation of embedded rebar elements within the concrete mesh and the handling of the resulting discretized numerical model in parallel and distributed computing environments will also be presented.

 

Prof. Seiichi Koshizuka (The University of Tokyo, Japan)

Numerical Simulation for Nuclear Plant Safety in Tsunami using Particle Method

Particle methods have been developed for analyzing multiphase flows with phase change as well as violent free surface flows. Meshless discretization is a significant advantage for such complex phenomena. Spreading of the molten nuclear reactor core is analyzed by the particle method in the case of severe accidents. Solidification is modeled by fixing the relative motion of the moving particles. Large-scale tsunami run-up simulation is carried out on the nuclear plant site using the tsunami wave of the Great East Japan Earthquake in 2011. The calculated flooded area agrees well with the observation. Inundation in a turbine building is also analyzed because the blackout of emergency power is caused by the internal flooding. Floating objects are considered as fluid-rigid body coupling problems. The trajectory of the floating object is extremely sensitive to the initial position, the coefficient of restitution, etc. The effect of the floating objects should be assessed by statistical approach, which is expected to study more in future.

 

Semi-Plenary Lectures

Prof. Olivier Allix (ENS Cachan, France)

Computational Damage and Failure Analysis of Laminates across the Scales: Progress and Challenges

The precise sizing of composite laminates requires taking into account various deterioration scenarios that happen at the micro or meso scale.  The presentation will discuss those different scenarios and the associated modeling and computational strategies to take them into account. Two objectives will be mainly considered: the virtual testing of composite material and the virtual testing of composite structure.

 

Prof. Yuri Bazilevs (Brown University, USA)

Recent Developments in Immersed IGA-Meshfree Methods for Extreme-Event Simulation

This presentation is focused on Isogeometric Analysis (IGA) and RKPM Meshfree method with applications to extreme-events simulation. A novel framework for air-blast-structure interaction (ABSI) based on an immersed approach coupling IGA and RKPM is presented and verified on a set of challenging examples. Several numerical challenges exist for carrying out the aforementioned simulations, and these are addressed in the present work. The challenges include shock capturing in both the fluid and solid parts of the problem, and addressing near incompressibility, which is important in the presence of plastic deformations. Extension of the proposed ABSI framework to handle energetic materials is also presented.

 

Prof. Kwok Fai Cheung (University of Hawaii at Manoa, USA)

Non-hydrostatic Modeling of Tsunami Waves

Tsunami waves are weakly dispersive and are often modeled using the nonlinear shallow-water equations. Inclusion of non-hydrostatic terms in the governing equations and a shock-capturing scheme in the numerical formulation greatly enhances the modeling capability for academic research and real-world application. The resulting non-hydrostatic model can describe tsunami generation from time-varying seafloor deformation making the calculation compatible with finite-fault modeling of earthquake rupture. Dispersion is a non-hydrostatic property responsible for the spatiotemporal resolution of trans-oceanic propagation that influences the subsequent coastal inundation processes. Shock-capturing is instrumental for modeling tsunami bores near the coast and conserving flow volume in flood hazard mapping. My talk will give an overview of the non-hydrostatic model NEOWAVE developed at the University of Hawaii and highlight its implementation in a series of cross-disciplinary studies from megathrust rupture mechanisms to tsunami loads on buildings.

 

Prof. John E. Dolbow (Duke University, USA)

Gradient-Based Damage Methods for Cohesive Models of Dynamic Fracture and Fragmentation

Recently, much attention has focused on gradient-based damage and phase-field models for fracture problems. In these methods, sharp fracture surfaces are approximated with a scalar damage field that varies continuously throughout the domain. The evolution of the damage field is determined by a secondary equation that incorporates a length scale for regularization. These models have enabled simulations of complex fracture problems in three dimensions and demonstrated robustness for simulating challenging phenomena such as crack bifurcations and coalescence. Many of these approaches have been based on a variational formulation for Griffithtype fracture models. While these approaches have seen considerable success, they have also suffered from a number of shortcomings when applied to dynamic fracture and fragmentation. These include, for example, the computational cost of a global reaction-diffusion auxiliary equation and challenges associated with introducing critical thresholds that trigger the onset of damage. In this talk, we describe an alternative approach that is based on recent work establishing links between gradient damage methods and cohesive-type models of fracture. The approach naturally introduces a threshold for the onset of damage and allows for the fracture properties to remain fixed as the regularization length scale vanishes. We will discuss strategies for enforcing irreversibility in these approaches, modifications for anisotropic failure problems, and methods to transition from continuous to discontinuous representations of the fracture surfaces. Finally, applications of these models to a range of problems in quasi-static and dynamic fracture in quasi-brittle systems will be presented.

 

Prof. P. Benson Shing (University of California, San Diego, USA)

Advancing the Seismic Design of Reinforced Concrete and Masonry Structures using Computational Models and Large-scale Experiments

To design new civil infrastructure systems or evaluate the safety of existing structures for extreme seismic events, engineers have been increasingly relying on nonlinear computational tools. Detailed nonlinear computational models have been used to develop design and detailing requirements in code provisions and to assess the collapse resilience of structural systems. However, the nonlinear behavior of a structural system involves many complicated mechanisms developed in the material, structural component, and system. For reinforced concrete structures, critical mechanisms include the cracking of concrete, the plastic and fracture behavior of reinforcing bars, the buckling of reinforcing bars, and the interaction between the bars and the surrounding concrete. How the structural components behave eventually depends on the reinforcing details and their interaction with other components in the structural system. This presentation will cover some recent work of the presenter’s research group on the development of practical physics-based computational models with application to the seismic design of reinforced concrete bridge structures and reinforced masonry buildings. The validation of these models with large-scale experiments conducted on structural components and systems will also be shown.

 

Prof. Kenichi Soga (University of California, Berkeley, USA)

Large Deformation Modeling and Simulations of Landslides

Traditional geotechnical analyses for landslides involve failure prediction (i.e. onset of failure) and the design of structures that can safely withstand the applied loads. But the analyses provide limited information on the post-failure behavior such as failure geometry and the rate of movement. Modern numerical methods for large deformation simulations are now emerging and some of them are started to be adopted by geotechnical engineers to simulate large mass movements. There is also a broader impact because of their potential ability to evaluate the risks of catastrophic damage if a landslide occurs. In this talk, various large-deformation analysis methods are introduced and their applicability for solving landslide problems is discussed. In particular, a technique called the material point method (MPM) is attractive because it allows numerical implementation of history-dependent soil constitutive models and boundary conditions commonly used in geotechnical analysis in a relatively straightforward manner. The recent theoretical development on the multi-phase soil-fluid coupled MPM framework is also providing an opportunity to simulate catastrophic landslides involve seepage forces. On the other hand, further development is required to build confidence in the engineering community to use large deformation simulation methods in engineering practice. This includes better appreciation of failure development processes such as softening induced shear band formation and tensile cracking, identification of energy dissipation mechanisms that are responsible for runout distance and rate, and the role of thermo-hydro-mechanical interaction on these processes and mechanisms.

 

Prof. Zhuo Zhuang (Tsinghua University, China)

Data-driving-based Theoretical and Numerical Fracking Models to Optimize Recovery Efficiency in Shale

Hydraulic fracture (fracking) technology in gas shale field engineering is highly developed last decades in North America and also recent years in China, but the knowledge of actual fracturing process is mostly empirical and makes the mechanician wonder. In this work, the data-driving-based theoretical and numerical fracking model is proposed to predict and optimize recovery efficiency in shale. Shale is a typical layered and anisotropic material whose properties are characterized primarily by locally oriented anisotropic clay minerals and naturally formed bedding planes. The debonding of bedding planes will greatly influence the shale fracking to form a large-scale highly permeable fracture network, which is the stimulate reform volume (SRV). Both theoretical and numerical models are developed to quantitatively predict the growth of debonding zone in layered shale under fracking. Some parameters are proposed to characterize the corresponding conditions of tensile and shear debonding of bedding planes. It is found that debonding is mainly caused by the shear failure of bedding planes in the actual reservoir. Then the theoretical model is applied to design the perforation cluster spacing to optimize SRV. The SRV and optimal perforation cluster spacing range can be quantitatively calculated to guide the fracking design. These results are comparable with the data from real-time signal evolution by micro earthquake monitor in field engineering.