  This paper integrates nonlinear-manifold reduced order models (NM-ROMs) with
domain decomposition (DD). NM-ROMs approximate the FOM state in a
nonlinear-manifold by training a shallow, sparse autoencoder using FOM snapshot
data. These NM-ROMs can be advantageous over linear-subspace ROMs (LS-ROMs) for
problems with slowly decaying Kolmogorov n-width. However, the number of NM-ROM
parameters that need to be trained scales with the size of the FOM. Moreover,
for "extreme-scale" problems, the storage of high-dimensional FOM snapshots
alone can make ROM training expensive. To alleviate the training cost, this
paper applies DD to the FOM, computes NM-ROMs on each subdomain, and couples
them to obtain a global NM-ROM. This approach has several advantages: Subdomain
NM-ROMs can be trained in parallel, involve fewer parameters to be trained than
global NM-ROMs, require smaller subdomain FOM dimensional training data, and
can be tailored to subdomain-specific features of the FOM. The shallow, sparse
architecture of the autoencoder used in each subdomain NM-ROM allows
application of hyper-reduction (HR), reducing the complexity caused by
nonlinearity and yielding computational speedup of the NM-ROM. This paper
provides the first application of NM-ROM (with HR) to a DD problem. In
particular, it details an algebraic DD reformulation of the FOM, trains a
NM-ROM with HR for each subdomain, and develops a sequential quadratic
programming (SQP) solver to evaluate the coupled global NM-ROM. Theoretical
convergence results for the SQP method and a priori and a posteriori error
estimates for the DD NM-ROM with HR are provided. The proposed DD NM-ROM with
HR approach is numerically compared to a DD LS-ROM with HR on the 2D
steady-state Burgers' equation, showing an order of magnitude improvement in
accuracy of the proposed DD NM-ROM over the DD LS-ROM.
