Researcher profile

Gianluca Iaccarino

Gianluca Iaccarino contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

6 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.

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.

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.

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.

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.

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.