Source author record

David Amsallem

David Amsallem appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2020arXiv

Gradient-based Constrained Optimization Using a Database of Linear Reduced-Order Models

A methodology grounded in model reduction is presented for accelerating the gradient-based solution of a family of linear or nonlinear constrained optimization problems where the constraints include at least one linear Partial Differential Equation (PDE). A key component of this methodology is the construction, during an offline phase, of a database of pointwise, linear, Projection-based Reduced-Order Models (PROM)s associated with a design parameter space and the linear PDE(s). A parameter sampling procedure based on an appropriate saturation assumption is proposed to maximize the efficiency of such a database of PROMs. A real-time method is also presented for interpolating at any queried but unsampled parameter vector in the design parameter space the relevant sensitivities of a PROM. The practical feasibility, computational advantages, and performance of the proposed methodology are demonstrated for several realistic, nonlinear, aerodynamic shape optimization problems governed by linear aeroelastic constraints.

preprint2016arXiv

Robust Model Reduction Of Hyperbolic Problems by $L^1$-norm Minimization and Dictionary Approximation

We propose a novel model reduction approach for the approximation of non linear hyperbolic equations in the scalar and the system cases. The approach relies on an offline computation of a dictionary of solutions together with an online $L^1$-norm minimization of the residual. It is shown why this is a natural framework for hyperbolic problems and tested on nonlinear problems such as Burgers' equation and the one-dimensional Euler equations involving shocks and discontinuities. Efficient algorithms are presented for the computation of the $L^1$-norm minimizer, both in the cases of linear and nonlinear residuals. Results indicate that the method has the potential of being accurate when involving only very few modes, generating physically acceptable, oscillation-free, solutions.

preprint2015arXiv

Linearized Aeroelastic Computations in the Frequency Domain Based on Computational Fluid Dynamics

An iterative, CFD-based approach for aeroelastic computations in the frequency domain is presented. The method relies on a linearized formulation of the aeroelastic problem and a fixed-point iteration approach and enables the computation of the eigenproperties of each of the wet aeroelastic eigenmodes. Numerical experiments on the aeroelastic analysis and design optimization of two wing configurations illustrate the capability of the method for the fast and accurate aeroelastic analysis of aircraft configurations and its advantage over classical time-domain approaches.

preprint2015arXiv

Projection-based model reduction for contact problems

To be feasible for computationally intensive applications such as parametric studies, optimization and control design, large-scale finite element analysis requires model order reduction. This is particularly true in nonlinear settings that tend to dramatically increase computational complexity. Although significant progress has been achieved in the development of computational approaches for the reduction of nonlinear computational mechanics models, addressing the issue of contact remains a major hurdle. To this effect, this paper introduces a projection-based model reduction approach for both static and dynamic contact problems. It features the application of a non-negative matrix factorization scheme to the construction of a positive reduced-order basis for the contact forces, and a greedy sampling algorithm coupled with an error indicator for achieving robustness with respect to model parameter variations. The proposed approach is successfully demonstrated for the reduction of several two-dimensional, simple, but representative contact and self contact computational models.

preprint2015arXiv

Real-time solution of computational problems using databases of parametric linear reduced-order models with arbitrary underlying meshes

A comprehensive approach for real-time computations using a database of parameterized linear reduced-order models (ROMs) is proposed. The method proceeds by sampling offline ROMs for specific values of the parameters and interpolating online the associated reduced operators. In the offline phase, a pre-processing step transforms the reduced operators into consistent sets of generalized coordinates prior to their interpolation. The present paper also introduces a consistency enforcement approach for models defined on arbitrary underlying meshes. In the online phase, the operators are interpolated on matrix manifolds. The proposed framework is illustrated on two realistic multi-physics problems: an inverse acoustic scattering problem around a submarine and flutter predictions for a wing-tank system. The second application is implemented on a mobile device, illustrating the capability of the proposed framework to operate in real-time.

preprint2015arXiv

Robust Model Reduction by $L^1$-norm Minimization and Approximation via Dictionaries: Application to Linear and Nonlinear Hyperbolic Problems

We propose a novel model reduction approach for the approximation of non linear hyperbolic equations in the scalar and the system cases. The approach relies on an offline computation of a dictionary of solutions together with an online $L^1$-norm minimization of the residual. It is shown why this is a natural framework for hyperbolic problems and tested on nonlinear problems such as Burgers' equation and the one-dimensional Euler equations involving shocks and discontinuities. Efficient algorithms are presented for the computation of the $L^1$-norm minimizer, both in the cases of linear and nonlinear residuals. Results indicate that the method has the potential of being accurate when involving only very few modes, generating physically acceptable, oscillation-free, solutions.

preprint2013arXiv

On the stability of projection-based linear reduced-order models: Descriptor vs non-descriptor forms

Two comprehensive approaches are considered for constructing projection-based reduced-order computational models for linear dynamical systems. The first one reduces the governing equations written in the descriptor form, using a Galerkin or Petrov-Galerkin projection onto a reduced-order basis or pair of them, respectively. These bases can be constructed by any preferred method such as the Proper Orthogonal Decomposition, Balanced Proper Orthogonal Decomposition, or Moment Matching method. The second approach transforms first the governing equations into their non-descriptor form before applying the same projection-based model order reduction method. For several structural and coupled fluid-structure dynamical systems, it is observed that the first reduction approach leads to reduced-order models that exhibit significantly better numerical stability and accuracy properties than those resulting from the second one. These observations are reinforced by theoretical results that anticipate and confirm the better stability properties obtained in general when reducing the descriptor rather than non-descriptor form of the equations governing a linear dynamical system.

preprint2012arXiv

High-order accurate difference schemes for the Hodgkin-Huxley equations

A novel approach for simulating potential propagation in neuronal branches with high accuracy is developed. The method relies on high-order accurate difference schemes using the Summation-By-Parts operators with weak boundary and interface conditions applied to the Hodgkin-Huxley equations. This work is the first demonstrating high accuracy for that equation. Several boundary conditions are considered including the non-standard one accounting for the soma presence, which is characterized by its own partial differential equation. Well-posedness for the continuous problem as well as stability of the discrete approximation is proved for all the boundary conditions. Gains in terms of CPU times are observed when high-order operators are used, demonstrating the advantage of the high-order schemes for simulating potential propagation in large neuronal trees.

preprint2012arXiv

The GNAT method for nonlinear model reduction: effective implementation and application to computational fluid dynamics and turbulent flows

The Gauss--Newton with approximated tensors (GNAT) method is a nonlinear model reduction method that operates on fully discretized computational models. It achieves dimension reduction by a Petrov--Galerkin projection associated with residual minimization; it delivers computational efficency by a hyper-reduction procedure based on the `gappy POD' technique. Originally presented in Ref. [1], where it was applied to implicit nonlinear structural-dynamics models, this method is further developed here and applied to the solution of a benchmark turbulent viscous flow problem. To begin, this paper develops global state-space error bounds that justify the method's design and highlight its advantages in terms of minimizing components of these error bounds. Next, the paper introduces a `sample mesh' concept that enables a distributed, computationally efficient implementation of the GNAT method in finite-volume-based computational-fluid-dynamics (CFD) codes. The suitability of GNAT for parameterized problems is highlighted with the solution of an academic problem featuring moving discontinuities. Finally, the capability of this method to reduce by orders of magnitude the core-hours required for large-scale CFD computations, while preserving accuracy, is demonstrated with the simulation of turbulent flow over the Ahmed body. For an instance of this benchmark problem with over 17 million degrees of freedom, GNAT outperforms several other nonlinear model-reduction methods, reduces the required computational resources by more than two orders of magnitude, and delivers a solution that differs by less than 1% from its high-dimensional counterpart.