In this contribution, a finite element scheme to impose mixed boundary conditions without introducing Lagrange multipliers is presented for hyperbolic systems described as port-Hamiltonian systems. The strategy relies on finite element exterior calculus and domain decomposition to interconnect two systems with dual input-output behavior. The spatial domain is split into two parts by introducing an arbitrary interface. Each subdomain is discretized with a mixed finite element formulation that introduces a uniform boundary condition in a natural way as the input. In each subdomain the finite element spaces are selected from a finite element subcomplex to obtain a stable discretization. The two systems are then interconnected together by making use of a feedback interconnection. This is achieved by discretizing the boundary inputs using appropriate spaces that couple the two formulations. The final systems include all boundary conditions explicitly and do not contain any Lagrange multiplier. Time integration is performed using the implicit midpoint or Störmer-Verlet scheme. The method can also be applied to semilinear systems containing algebraic nonlinearities. The proposed strategy is tested on different examples: geometrically exact intrinsic beam model, the wave equation, membrane elastodynamics and the Mindlin plate. Numerical tests assess the conservation properties of the scheme, the effectiveness of the methodology and its robustness against shear locking phenomena.

High computational costs are encountered in topology optimization problems of geometrically nonlinear structures since intensive use has to be made of incremental-iterative finite element simulations. To alleviate this computational intensity, reduced-order models (ROMs) are explored in this paper. The proposed method targets ROM bases consisting of a relatively small set of base vectors while accuracy is still guaranteed. For this, several fully automated update and maintenance techniques for the ROM basis are investigated and combined. In order to remain effective for flexible structures, path derivatives are added to the ROM basis. The corresponding sensitivity analysis (SA) strategies are presented and the accuracy and efficiency are examined. Various geometrically nonlinear examples involving both solid as well as shell elements are studied to test the proposed ROM techniques. Test cases demonstrates that the set of degrees of freedom appearing in the nonlinear equilibrium equation typically reduces to several tenth. Test cases show a reduction of up to 6 times fewer full system updates.

We introduce a mimetic dual-field discretization which conserves mass, kinetic energy and helicity for three-dimensional incompressible Navier-Stokes equations. The discretization makes use of a conservative dual-field mixed weak formulation where two evolution equations of velocity are employed and dual representations of the solution are sought for each variable. A temporal discretization, which staggers the evolution equations and handles the nonlinearity such that the resulting discrete algebraic systems are linear and decoupled, is constructed. The spatial discretization is mimetic in the sense that the finite dimensional function spaces form a discrete de Rham complex. Conservation of mass, kinetic energy and helicity in the absence of dissipative terms is proven at the discrete level. Proper dissipation rates of kinetic energy and helicity in the viscous case are also proven. Numerical tests supporting the method are provided.