Seminar in Computational Engineering

Thursday May 1, 2014 │ 4:00 PM │ 4-237

Hurricane Storm Surge Models Using Integrated Ocean Basin to Shelf to Inland Floodplain Unstructured Grids
Joannes J. Westerink
Joseph and Nona Ahearn Professor in Computational Science and Engineering, Henry J. Massman Chairman, Department of Civil & Environmental Engineering & Earth Sciences
University of Notre Dame

Hurricane wind wave, storm surge, and current environments in the coastal ocean and adjacent coastal floodplain are characterized by their high energy and by their spatial variability. These processes impact offshore energy assets, navigation, ports and harbors, deltas, wetlands, and coastal communities. The potential for an enormous catastrophic impact in terms of loss of life and economic losses is substantial.

Computational models for wind waves and storm driven currents and surge must provide a high level of grid resolution, fully couple the energetic processes, and perform quickly for risk assessment, flood mitigation system design, and forecasting purposes. In order to accomplish this, high performance scalable codes are essential. To this end, we have developed an MPI based domain decomposed unstructured grid framework that minimizes global communications, efficiently handles localized sub-domain to sub-domain communication, applies a local inter-model paradigm with all model to model communications being kept on identical cores for sub-domains, and carefully manages output by assigning specialized cores for this purpose. Continuous Galerkin (CG) and Discontinuous Galerkin (DG) implementations are examined. Performance of explicit and implicit implementations of the wave-current coupled system on up to 32,000 cores for various platforms is evaluated.

The system has been extensively validated with an ever increasing amount of wave, water level and current data that has being collected for recent storms including Hurricanes Katrina (2005), Rita (2005), Gustav (2008), Ike (2008), and Sandy (2012). The modeling system helps understand the physics of hurricane storm surges including processes such as geostrophically driven forerunner, shelf waves that propagate far away from the storm, wind wave – surge interaction, surge capture and propagation by protruding deltaic river systems, the influence of storm size and forward speed, and frictionally controlled inland penetration.

These models are being applied by the US Army Corps of Engineers (USACE) in the development of the recently completed hurricane risk reduction system in Southern Louisiana as well as for the development of FEMA Digital Flood Insurance Rate Maps (DFIRMS) for Texas, Louisiana, Mississippi, and other Gulf and Atlantic coast states. NOAA applies the models in extra-tropical and tropical storm surge forecasting.

Current algorithmic development is focused on DG solvers, ideally suited for the associated strongly advective flows. Due to the larger numbers of degrees of freedom for a specific grid, DG solutions have traditionally been more costly than CG solutions. It is demonstrated that high order implementations of DG leads to several orders of magnitude improvement in cost per accuracy performance as compared to lower order methods. In addition, loop level optimization further improves the efficiency of DG solutions by a factor of 4 to 5. It is noted that curved boundaries must be treated using super-parametric elements for p=1 and p=2 and iso-parametric elements for p=3 in order to achieve anticipated convergence rates.

SIAM & CCE Student Seminar Series

Thursday, April 10th at 4 PM in 4-237

Topological Microfluidics: Exploring anisotropic fluids in microfluidic environment
Anupam Sengupta, Roman Stocker group, CEE

Liquid crystals (LCs) are complex anisotropic fluids, well-known for display applications. Their properties are in contrast to the isotropic fluids which we typically encounter in state-of- the-art microfluidic science and technology. This allows us to explore LCs as a suitable functional material for developing a microfluidic technology which harnesses the coupling interactions between the flow, the molecular orientation, and the ordering or topology of the system. The anisotropic coupling, along with the spontaneous but mandatory evolution of the topological constraints, makes liquid crystal microfluidics distinct from its conventional counterpart. A feature which is particularly promising for future applications is that topological defects in liquid crystals can serve as effective traps for colloidal particles, aqueous droplets or other microfluidic cargo. By tailoring the boundary conditions and the overall geometry, topological defects in the liquid crystal bulk phase were created at will, and were guided towards a specific target. The structure of defects emerging in the system, which are otherwise considered nuisance in applications, can be used as a self-assembled system of rails, along which droplets containing the materials of interest could be transported. Harnessing the fluid anisotropy in micro-scale flows provides a novel pathway to create a range of potential optofluidic and metamaterial applications.

Reference: A. Sengupta, ‘Topological Microfluidics’, Springer International Publishing, Switzerland, 2013. ISBN: 978-3-319-00857-8.