Project Description
Fierro (LANL code number O4659) is a modern C++ code that offers novel capabilities to simulate multiscale multiphysics engineering and material science applications and for autonomous multiphysics design optimization. Fierro contains a suite of novel solvers for quasi-static solid mechanics problems and transient, compressible material dynamic problems with Lagrangian methods, which have meshes with constant mass elements that move with the material. Multiscale material models are coupled to world-leading numerical methods inside Fierro to simulate the performance of parts accounting for the underlying microstructure in the part. Fierro supports hybrid parallelism, using MPI for coarse-grained parallelism plus a range of hardware-specific languages for fine-grained parallelism to leverage diverse GPUs for efficient simulation runtimes. The combined capabilities in the Fierro code are key to realizing the full-potential of diverse manufacturing processes by leveraging modern high-performance computing machines to: (1) simulate large-deformation, high-shear manufacturing processes; (2) assess the performance of parts across physical regimes accounting for the manufacturing processes, and (3) design optimal parts autonomously that are additively manufactured. The parallel performance of Fierro also enables assessments of applications that require meshes with millions to billions of elements, as well as, quantification of uncertainty and margins, or sensitivity studies, that may require 100's to 1000's of simulations.
Fierro is advantageous for:
- Multiphysics multiscale assessments
- Multiphysics multi-constraint topology optimization
- Multiscale modeling
- Material model research
- Numerical methods research
- Computer science research
Fierro is built on the ELEMENTS library that supports a diverse suite of element types, including high-order elements, and quadrature rules. The mesh class within the ELEMENTS library is designed for efficient calculations on unstructured meshes and to minimize memory usage. Fierro is designed to readily accommodate a range of numerical methods including continuous finite element, finite volume, and discontinuous Galerkin methods. Fierro is designed to support explicit and implicit time integration methods as well as optimization methods.
The linear Lagrangian finite element method for material dynamics in Fierro supports user developed material models and multiscale models for microstructure-aware simulations. These multiscale models can also be run in a stand-alone manner to establish structure-property relationships. This linear Lagrangian finite element method is coupled to optimization solvers for 3D transient multi-physics topology optimzation.
Fierro uniquely offers several compact-stencil, arbitrary-order Lagrangian finite element method for more efficient simulations. These schemes use meshes that edges that can bend to accurately track large deformations in material dynamics.
A physics solver is supported in Fierro for simulating static or quasistatic thermal and mechanical applications. Higher-order optimization methods work in concert with the thermal-mechanical solvers for 3D multiphysics static topology optimization that satisfies multiple constraints.
Applications
The Fierro collection of solvers includes separate tools to help users setup diverse simulations. Some solvers and work flows have Graphical User Interfaces for ease of use.
Micro-mechanical solvers
Our small strain and large strain Elasto-Viscoplastic Fast Fourier Transform (EVPFFT) solvers predict the micro-mechanical response and microstructure evolution of polycrystalline materials. This tool can be used to simulate the bulk response of a polycrystalline material under desired loading conditions. The EVPFFT solvers can also be used to calculate the bulk-scale properties of 3D printed lattices.
Static and quasi-static physics solver
Our static and quasi-static finite element physics solver can solve boundary value problems over a meshed geometry with provided boundary conditions (Neumann or Dirichlet).
Material dynamics physics solvers
Our finite element material dynamics physics solver can simulate elastic structural dynamics to shock-driven flows given a range of boundary conditions, external forces, and initial conditions.