BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Gianluca Iaccarino

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

ResearcherUnclaimed source record
math.NAmath.OCphysics.flu-dynComputational Engineering, Finance, and ScienceNumerical Analysisphysics.comp-phComputationDistributed, Parallel, and Cluster Computingmath-phmath.MPNeural and Evolutionary Computingphysics.app-phphysics.data-anRobotics

Catalog footprint

What is connected

16works
14topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

16 published item(s)

preprint2023arXiv

Simultaneous Identification and Denoising of Dynamical Systems

In recent years there has been a push to discover the governing equations dynamical systems directly from measurements of the state, often motivated by systems that are too complex to directly model. Although there has been substantial work put into such a discovery, doing so in the case of large noise has proved challenging. Here we develop an algorithm for Simultaneous Identification and Denoising of a Dynamical System (SIDDS). We infer the noise in the state measurements by requiring that the denoised data satisfies the dynamical system with an equality constraint. This is unlike existing work where the mismatch in the dynamical system is treated as a penalty in the objective. We assume the dynamics is represented in a pre-defined basis and develop a sequential quadratic programming approach to solve the SIDDS problem featuring a direct solution of KKT system with a specialized preconditioner. In addition, we show how we can include sparsity promoting regularization using an iteratively reweighted least squares approach. The resulting algorithm leads to estimates of the dynamical system that approximately achieve the Cramér-Rao lower bound and, with sparsity promotion, can correctly identify the sparsity structure for higher levels of noise than existing techniques. Moreover, because SIDDS decouples the data from the evolution of the dynamical system, we show how to modify the problem to accurately identify systems from low sample rate measurements. The inverse problem approach and solution framework used by SIDDS has the potential to be expanded to related problems identifying governing equations from noisy data.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An integrated heterogeneous computing framework for ensemble simulations of laser-induced ignition

An integrated computational framework is introduced to study complex engineering systems through physics-based ensemble simulations on heterogeneous supercomputers. The framework is primarily designed for the quantitative assessment of laser-induced ignition in rocket engines. We develop and combine an implicit programming system, a compressible reacting flow solver, and a data generation/management strategy on a robust and portable platform. We systematically present this framework using test problems on a hybrid CPU/GPU machine. Efficiency, scalability, and accuracy of the solver are comprehensively assessed with canonical unit problems. Ensemble data management and autoencoding are demonstrated using a canonical diffusion flame case. Sensitivity analysis of the ignition of a turbulent, gaseous fuel jet is performed using a simplified, three-dimensional model combustor. Our approach unifies computer science, physics and engineering, and data science to realize a cross-disciplinary workflow. The framework is exascale-oriented and can be considered a benchmark for future computational science studies of real-world systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Spinning-enabled Wireless Amphibious Origami Millirobot

Wireless millimeter-scale origami robots that can locomote in narrow spaces and morph their shapes have recently been explored with great potential for biomedical applications. Existing millimeter-scale origami devices usually require separate geometrical components for locomotion and functions, which increases the complexity of the robotic systems and their operation upon limited locomotion modes. Additionally, none of them can achieve both on-ground and in-water locomotion. Here we report a magnetically actuated amphibious origami millirobot that integrates capabilities of spinning-enabled multimodal locomotion, controlled delivery of liquid medicine, and cargo transportation with wireless operation. This millirobot takes full advantage of the geometrical features and folding/unfolding capability of Kresling origami, a triangulated hollow cylinder, to fulfill multifunction: its geometrical features are exploited for generating omnidirectional locomotion in various working environments, including on unstructured ground, in liquids, and at air-liquid interfaces through rolling, flipping, and spinning-induced propulsion; the folding/unfolding is utilized as a pumping mechanism for integrated multifunctionality such as controlled delivery of liquid medicine; furthermore, the spinning motion provides a sucking mechanism for targeted solid cargo transportation. This origami millirobot breaks the conventional way of utilizing origami folding only for shape reconfiguration and integrates multiple functions in one simple body. We anticipate the reported magnetic amphibious origami millirobots have the potential to serve as minimally invasive devices for biomedical diagnoses and treatments.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Task-parallel in-situ temporal compression of large-scale computational fluid dynamics data

Present day computational fluid dynamics simulations generate extremely large amounts of data, sometimes on the order of TB/s. Often, a significant fraction of this data is discarded because current storage systems are unable to keep pace. To address this, data compression algorithms can be applied to data arrays containing flow quantities of interest to reduce the overall amount of storage. Compression methods either exactly reconstruct the original dataset (lossless compression) or provide an approximate representation of the original dataset (lossy compression). The matrix column interpolative decomposition (ID) can be implemented as a type of lossy compression for data matrices that factors the original data matrix into a product of two smaller factor matrices. One of these matrices consists of a subset of the columns of the original data matrix, while the other is a coefficient matrix which approximates the columns of the original data matrix as linear combinations of the selected columns. Motivating this work is the observation that the structure of ID algorithms makes them a natural fit for the asynchronous nature of task-based parallelism; they are able to operate independently on sub-domains of the system of interest and, as a result, provide varied levels of compression. Using the task-based Legion programming model, a single-pass ID algorithm (SPID) for CFD applications is implemented. Performance studies, scalability, and the accuracy of the compression algorithms are presented for an analytical Taylor-Green vortex problem, followed by a large-scale implementation of a compressible Taylor-Green vortex using a high-order Navier-Stokes solver. In both cases, compression factors exceeding 100 are achieved with relative errors at or below 10e-3. Moreover, strong and weak scaling results demonstrate that introducing SPID to solvers leads to negligible increases in runtime.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pass-efficient methods for compression of high-dimensional turbulent flow data

The future of high-performance computing, specifically on future Exascale computers, will presumably see memory capacity and bandwidth fail to keep pace with data generated, for instance, from massively parallel partial differential equation (PDE) systems. Current strategies proposed to address this bottleneck entail the omission of large fractions of data, as well as the incorporation of $\textit{in situ}$ compression algorithms to avoid overuse of memory. To ensure that post-processing operations are successful, this must be done in a way that a sufficiently accurate representation of the solution is stored. Moreover, in situations where the input/output system becomes a bottleneck in analysis, visualization, etc., or the execution of the PDE solver is expensive, the the number of passes made over the data must be minimized. In the interest of addressing this problem, this work focuses on the utility of pass-efficient, parallelizable, low-rank, matrix decomposition methods in compressing high-dimensional simulation data from turbulent flows. A particular emphasis is placed on using coarse representation of the data -- compatible with the PDE discretization grid -- to accelerate the construction of the low-rank factorization. This includes the presentation of a novel single-pass matrix decomposition algorithm for computing the so-called interpolative decomposition. The methods are described extensively and numerical experiments on two turbulent channel flow data are performed. In the first (unladen) channel flow case, compression factors exceeding $400$ are achieved while maintaining accuracy with respect to first- and second-order flow statistics. In the particle-laden case, compression factors of 100 are achieved and the compressed data is used to recover particle velocities.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Estimated aleatoric uncertainty from initial and inlet conditions for variable density mixing

Variable density flows occur in a variety of different systems with a wide range of scales, from astrophysics to atmospheric flows to inertial confinement fusion or reacting flows. Given the inherent limitations of RANS simulations, it is important to find ways to quantify the uncertainty in the predictions. The aleatoric uncertainty of two different variable density flows to inlet and initial conditions, respectively, is studied. The two cases are the turbulent mixing of a jet in a co-flow at small Atwood number $A$, and the Rayleigh-Taylor mixing in a tilted rocket rig at medium $A$. Uncertainty estimates are made for simulations using the BHR-1 turbulence model and compared to reference data. Estimated distributions of measurement uncertainties from an experimental report are used to do Monte-Carlo sampling for some inlet parameters of the jet in a co-flow case. The results show only little sensitivity to all inlet parameter perturbations, and even with increased perturbations the estimated uncertainty band does not fully match the experimental results. Therefore, the model predictions are very robust with respect to the inlet parameter variability. For the tilted rocket rig case, a perturbation to the initial interface location is varied over an expected range of parameter values. The simulation results are significantly more sensitive to the initial conditions than the jet in co-flow simulation results are to the inlet conditions. The estimated uncertainty band of the tilted rocket rig matches quite well the experimental data at nominal perturbation amplitude. The predictions of this Rayleigh-Taylor instability case hence are less robust to the initial interface perturbations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Many physical laws are ridge functions

A ridge function is a function of several variables that is constant along certain directions in its domain. Using classical dimensional analysis, we show that many physical laws are ridge functions; this fact yields insight into the structure of physical laws and motivates further study into ridge functions and their properties. We also connect dimensional analysis to modern subspace-based techniques for dimension reduction, including active subspaces in deterministic approximation and sufficient dimension reduction in statistical regression.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A density-matching approach for optimization under uncertainty

Modern computers enable methods for design optimization that account for uncertainty in the system---so-called optimization under uncertainty. We propose a metric for OUU that measures the distance between a designer-specified probability density function of the system response the target and system response's density function at a given design. We study an OUU formulation that minimizes this distance metric over all designs. We discretize the objective function with numerical quadrature and approximate the response density function with a Gaussian kernel density estimate. We offer heuristics for addressing issues that arise in this formulation, and we apply the approach to a CFD-based airfoil shape optimization problem. We qualitatively compare the density-matching approach to a multi-objective robust design optimization to gain insight into the method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Exploiting Active Subspaces to Quantify Uncertainty in the Numerical Simulation of the HyShot II Scramjet

We present a computational analysis of the reactive flow in a hypersonic scramjet engine with focus on effects of uncertainties in the operating conditions. We employ a novel methodology based on active subspaces to characterize the effects of the input uncertainty on the scramjet performance. The active subspace identifies one-dimensional structure in the map from simulation inputs to quantity of interest that allows us to reparameterize the operating conditions; instead of seven physical parameters, we can use a single derived active variable. This dimension reduction enables otherwise infeasible uncertainty quantification, considering the simulation cost of roughly 9500 CPU-hours per run. For two values of the fuel injection rate, we use a total of 68 simulations to (i) identify the parameters that contribute the most to the variation in the output quantity of interest, (ii) estimate upper and lower bounds on the quantity of interest, (iii) classify sets of operating conditions as safe or unsafe corresponding to a threshold on the output quantity of interest, and (iv) estimate a cumulative distribution function for the quantity of interest.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Flexible Uncertainty Quantification Framework for General Multi-Physics Systems

We present a "module-based hybrid" Uncertainty Quantification (UQ) framework for general nonlinear multi-physics simulation. The proposed methodology, introduced in [\hyperlink{ref1}{1}], supports the independent development of each \emph{stochastic} linear or nonlinear physics module equipped with the most suitable probabilistic UQ method: non-intrusive, semi-intrusive or intrusive; and provides a generic framework to couple these stochastic simulation components. Moreover, the methodology is illustrated using a common "global" uncertainty representation scheme based on generalized polynomial chaos (gPC) expansions of inputs and outputs. By using thermally-driven cavity flow as the multi-physics model problem, we demonstrate the utility of our framework and report the computational gains achieved.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An Efficient Intrusive Uncertainty Propagation Method For Multi-Physics System With Random Inputs

Coupled partial differential equation (PDE) systems, which often represent multi-physics models, are naturally suited for modular numerical solution methods. However, several challenges yet remain in extending the benefits of modularization practices to the task of uncertainty propagation. Since the cost of each deterministic PDE solve can be usually expected to be quite significant, statistical sampling based methods like Monte-Carlo (MC) are inefficient because they do not take advantage of the mathematical structure of the problem, and suffer for poor convergence properties. On the other hand, even if each module contains a moderate number of uncertain parameters, implementing spectral methods on the combined high-dimensional parameter space can be prohibitively expensive due to the curse of dimensionality. In this work, we present a module-based and efficient intrusive spectral projection (ISP) method for uncertainty propagation. In our proposed method, each subproblem is separated and modularized via block Gauss-Seidel (BGS) techniques, such that each module only needs to tackle the local stochastic parameter space. Moreover, the computational costs are significantly mitigated by constructing reduced chaos approximations of the input data that enter each module. We demonstrate implementations of our proposed method and its computational gains over the standard ISP method using numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An Efficient Non-Intrusive Uncertainty Propagation Method for Stochastic Multi-Physics Models

Multi-physics models governed by coupled partial differential equation (PDE) systems, are naturally suited for partitioned, or modular numerical solution strategies. Although widely used in tackling deterministic coupled models, several challenges arise in extending the benefits of modularization to uncertainty propagation. On one hand, Monte-Carlo (MC) based methods are prohibitively expensive as the cost of each deterministic PDE solve is usually quite large, while on the other hand, even if each module contains a moderate number of uncertain parameters, implementing spectral methods on the combined high-dimensional parameter space can be prohibitively expensive. In this work, we present a reduced non-intrusive spectral projection (NISP) based uncertainty propagation method which separates and modularizes the uncertainty propagation task in each subproblem using block Gauss-Seidel (BGS) techniques. The overall computational costs in the proposed method are also mitigated by constructing reduced approximations of the input data entering each module. These reduced approximations and the corresponding quadrature rules are constructed via simple linear algebra transformations. We describe these components of the proposed algorithm assuming a generalized polynomial chaos (gPC) model of the stochastic solutions. We demonstrate our proposed method and its computational gains over the standard NISP method using numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Factorizing the Stochastic Galerkin System

Recent work has explored solver strategies for the linear system of equations arising from a spectral Galerkin approximation of the solution of PDEs with parameterized (or stochastic) inputs. We consider the related problem of a matrix equation whose matrix and right hand side depend on a set of parameters (e.g. a PDE with stochastic inputs semidiscretized in space) and examine the linear system arising from a similar Galerkin approximation of the solution. We derive a useful factorization of this system of equations, which yields bounds on the eigenvalues, clues to preconditioning, and a flexible implementation method for a wide array of problems. We complement this analysis with (i) a numerical study of preconditioners on a standard elliptic PDE test problem and (ii) a fluids application using existing CFD codes; the MATLAB codes used in the numerical studies are available online.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Optimization Under Uncertainty Using the Generalized Inverse Distribution Function

A framework for robust optimization under uncertainty based on the use of the generalized inverse distribution function (GIDF), also called quantile function, is here proposed. Compared to more classical approaches that rely on the usage of statistical moments as deterministic attributes that define the objectives of the optimization process, the inverse cumulative distribution function allows for the use of all the possible information available in the probabilistic domain. Furthermore, the use of a quantile based approach leads naturally to a multi-objective methodology which allows an a-posteriori selection of the candidate design based on risk/opportunity criteria defined by the designer. Finally, the error on the estimation of the objectives due to the resolution of the GIDF will be proven to be quantifiable

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-intrusive Low-Rank Separated Approximation of High-Dimensional Stochastic Models

This work proposes a sampling-based (non-intrusive) approach within the context of low-rank separated representations to tackle the issue of curse-of-dimensionality associated with the solution of models, e.g., PDEs/ODEs, with high-dimensional random inputs. Under some conditions discussed in details, the number of random realizations of the solution, required for a successful approximation, grows linearly with respect to the number of random inputs. The construction of the separated representation is achieved via a regularized alternating least-squares regression, together with an error indicator to estimate model parameters. The computational complexity of such a construction is quadratic in the number of random inputs. The performance of the method is investigated through its application to three numerical examples including two ODE problems with high-dimensional random inputs.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Spectral Methods for Parameterized Matrix Equations

We apply polynomial approximation methods -- known in the numerical PDEs context as spectral methods -- to approximate the vector-valued function that satisfies a linear system of equations where the matrix and the right hand side depend on a parameter. We derive both an interpolatory pseudospectral method and a residual-minimizing Galerkin method, and we show how each can be interpreted as solving a truncated infinite system of equations; the difference between the two methods lies in where the truncation occurs. Using classical theory, we derive asymptotic error estimates related to the region of analyticity of the solution, and we present a practical residual error estimate. We verify the results with two numerical examples.

0 citations0 reviewsNo code linkedOpen access