Numerical simulation by a domain-specific language for finite element method, FreeFem++


FreeFem++ is a software package developed by Professor Frédéric Hecht at Laboratoire Jacques-Louis Lions, Universite Pierre et Marie Curie and by others.
Since the package contains capability of mesh generation of finite element, interface to linear solvers for discretized equations, and visualization, user of FreeFem++ can concentrate on building of mathematical models, time discretization, and nonlinear solution.
The most advantage feature compared to general finite element analysis software is that FreeFem++ has a domain-specific language to describe discretization process of weak formulation of partial differential equations with appropriate data structure.
FreeFem++ was originally developed for educational purpose, but nowadays it has been improved to treat three dimensional finite elements and parallel computation with MPI.
FreeFem++ is the most advanced software to carry out prototype programming with finite element in the real world problems.
This seminar consists of a lecture by Prof. F. Hecht, presentations of example and some useful information by researchers who use FreeFem++ in Japan, and introduction of usage on supercomputers VCC and Octopus in Cybermedia Center.

Recommended to whom

  • wants to performe numerical simulation of fluid dynamics, structural analysis, and electromagnetic analysis
  • is intersted in migrating from general finite element analysis software
  • has expereice to use C/C++ language



13:30-14:30 Brain imaging with FreeFem++
Frédéric Hecht LJLL, Sorbonne Universités
    The abstract is here.
14:45-15:10 Numerical simulation of crack growth by FreeFem++
Takeshi Takaishi Dept. of Mathematical Engineering, Musashino University
     Numerical results of the crack growth model with phase field approach that is derived by Kimura and the author show the advantage of the adaptive mesh method.
15:10-15:35 Shape design of unsteady heat-convection fields
Eiji Takamine Dept. of Mechanical Engineering, National Institute of Technology, Gifu College
     The authors describe the shape design for controlling the temperature distribution history in the unsteady heat convection field to the target temperature distribution history. The square error integral between the actual temperature distributions and the target temperature distributions is used as the objective functional.Reshaping is carried out by the traction method proposed as an approach to solving shape optimization problems.The validity of proposed method is confirmed by results of 2D numerical analyses.
15:45-16:10 Shape optimization problem for suppressing time periodic flows
Takashi Nakazawa Center for Mathematical Modeling and Data Science, Osaka University
     In this presentation, the new shape optimization method is developed, where the eigenvalue of Snapshot Proper Orthogonal Decomposition(POD) is defined as the cost function and the nonstationary Navier-Stokes problem and the eigenvalue problem of snapshot POD are used as the constraint functions. Numerical results reveal that time periodic flow is suppressed by the shape optimization method.
16:10-16:35 Shape optimization analyses by FreeFem++
Hideyuki Azegami Graduate School of Informatics, Nagoya University
     Problems finding shapes of domains in boundary value problems of partial differential equations are solved by FreeFem++. The outline of formulation and solution of the shape optimization problem is introduced and applied to some engineering problems.
16:45-17:00 Usage of FreeFem++ on VCC/Octopus systems in Cybermedia Center
Atsushi Suzuki Cybermedia Center, Osaka University
     FreeFem++ is installed on VCC and Octopus systems linked with Intel Math Kernel Library to utilize latest multi-core CPU. FreeFem++ script to use direct linear solvers, Intel Pardiso, MUMPS, and Dissection for multi-core shared memory computing node is shown.


Date : 12 June, 2018, 1:30 p.m. - 5:00 p.m.
Instructor: Osaka University Cybermedia Center
Venue: Cybermedia Commons, Cybermedia Center Main Building 1F, Osaka University Suita Campus
     how to access
Type : classroom lecture
Quota: 30 people
Application deadline: 11 June, 2018, 5:00 p.m.