-
SNF-ROM: Projection-based nonlinear reduced order modeling with smooth neural fields
Authors:
Vedant Puri,
Aviral Prakash,
Levent Burak Kara,
Yongjie Jessica Zhang
Abstract:
Reduced order modeling lowers the computational cost of solving PDEs by learning a low-order spatial representation from data and dynamically evolving these representations using manifold projections of the governing equations. While commonly used, linear subspace reduced-order models (ROMs) are often suboptimal for problems with a slow decay of Kolmogorov $n$-width, such as advection-dominated fl…
▽ More
Reduced order modeling lowers the computational cost of solving PDEs by learning a low-order spatial representation from data and dynamically evolving these representations using manifold projections of the governing equations. While commonly used, linear subspace reduced-order models (ROMs) are often suboptimal for problems with a slow decay of Kolmogorov $n$-width, such as advection-dominated fluid flows at high Reynolds numbers. There has been a growing interest in nonlinear ROMs that use state-of-the-art representation learning techniques to accurately capture such phenomena with fewer degrees of freedom. We propose smooth neural field ROM (SNF-ROM), a nonlinear reduced modeling framework that combines grid-free reduced representations with Galerkin projection. The SNF-ROM architecture constrains the learned ROM trajectories to a smoothly varying path, which proves beneficial in the dynamics evaluation when the reduced manifold is traversed in accordance with the governing PDEs. Furthermore, we devise robust regularization schemes to ensure the learned neural fields are smooth and differentiable. This allows us to compute physics-based dynamics of the reduced system nonintrusively with automatic differentiation and evolve the reduced system with classical time-integrators. SNF-ROM leads to fast offline training as well as enhanced accuracy and stability during the online dynamics evaluation. Numerical experiments reveal that SNF-ROM is able to accelerate the full-order computation by up to $199\times$. We demonstrate the efficacy of SNF-ROM on a range of advection-dominated linear and nonlinear PDE problems where we consistently outperform state-of-the-art ROMs.
△ Less
Submitted 8 July, 2024; v1 submitted 17 May, 2024;
originally announced May 2024.
-
Projection-based reduced order modeling and data-driven artificial viscosity closures for incompressible fluid flows
Authors:
Aviral Prakash,
Yongjie Jessica Zhang
Abstract:
Projection-based reduced order models rely on offline-online model decomposition, where the data-based energetic spatial basis is used in the expensive offline stage to obtain equations of reduced states that evolve in time during the inexpensive online stage. The online stage requires a solution method for the dynamic evolution of the coupled system of pressure and velocity states for incompressi…
▽ More
Projection-based reduced order models rely on offline-online model decomposition, where the data-based energetic spatial basis is used in the expensive offline stage to obtain equations of reduced states that evolve in time during the inexpensive online stage. The online stage requires a solution method for the dynamic evolution of the coupled system of pressure and velocity states for incompressible fluid flows. The first contribution of this article is to demonstrate the applicability of the incremental pressure correction scheme for the dynamic evolution of pressure and velocity states. The evolution of a large number of these reduced states in the online stage can be expensive. In contrast, the accuracy significantly decreases if only a few reduced states are considered while not accounting for the interactions between unresolved and resolved states. The second contribution of this article is to compare three closure model forms based on global, modal and tensor artificial viscosity approximation to account for these interactions. The unknown model parameters are determined using two calibration techniques: least squares minimization of error in energy approximation and closure term approximation. This article demonstrates that an appropriate selection of solution methods and data-driven artificial viscosity closure models is essential for consistently accurate dynamics forecasting of incompressible fluid flows.
△ Less
Submitted 30 January, 2024;
originally announced January 2024.
-
A global unstructured, coupled, high-resolution hindcast of waves and storm surge
Authors:
Lorenzo Mentaschi,
Michalis Vousdoukas,
Guillermo Garcia-Sanchez,
Tomas Fernandez Montblanc,
Aron Roland,
Evangelos Voukouvalas,
Ivan Federico,
Ali Abdolali,
Yinglong J. Zhang,
Luc Feyen
Abstract:
Accurate information on waves and storm surges is essential to understand coastal hazards that are expected to increase in view of global warming and rising sea levels. Despite the recent advancement in development and application of large-scale coastal models, nearshore processes are still not sufficiently resolved due to coarse resolutions, transferring errors to coastal risk assessments and oth…
▽ More
Accurate information on waves and storm surges is essential to understand coastal hazards that are expected to increase in view of global warming and rising sea levels. Despite the recent advancement in development and application of large-scale coastal models, nearshore processes are still not sufficiently resolved due to coarse resolutions, transferring errors to coastal risk assessments and other large-scale applications. Here we developed a 50-year hindcast of waves and storm surges on an unstructured mesh of >650,000 nodes with an unprecedented resolution of 2-4 km at the global coast. Our modelling system is based on the circulation model SCHISM that is fully coupled with the WWM-V (WindWaveModel) and is forced by surface winds, pressure, and ice coverage from the ERA5 reanalysis. Results are compared with observations from satellite altimeters, tidal gauges and buoys, and show good skill for both Sea Surface Height (SSH) and Significant Wave Height (Hs), and a much-improved ability to reproduce the nearshore dynamics compared with previous, lower-resolution studies. Besides SSH, the modelling system also produces a range of other wave-related fields at each node of the mesh with a time step of 3 hours, including the spectral parameters of the first three largest energy peaks. This dataset offers the potential for more accurate global-scale applications on coastal hazard and risk
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
Search for relativistic fractionally charged particles in space
Authors:
DAMPE Collaboration,
F. Alemanno,
C. Altomare,
Q. An,
P. Azzarello,
F. C. T. Barbato,
P. Bernardini,
X. J. Bi,
M. S. Cai,
E. Casilli,
E. Catanzani,
J. Chang,
D. Y. Chen,
J. L. Chen,
Z. F. Chen,
M. Y. Cui,
T. S. Cui,
Y. X. Cui,
H. T. Dai,
A. De-Benedittis,
I. De Mitri,
F. de Palma,
M. Deliyergiyev,
A. Di Giovanni,
M. Di Santo
, et al. (126 additional authors not shown)
Abstract:
More than a century after the performance of the oil drop experiment, the possible existence of fractionally charged particles FCP still remains unsettled. The search for FCPs is crucial for some extensions of the Standard Model in particle physics. Most of the previously conducted searches for FCPs in cosmic rays were based on experiments underground or at high altitudes. However, there have been…
▽ More
More than a century after the performance of the oil drop experiment, the possible existence of fractionally charged particles FCP still remains unsettled. The search for FCPs is crucial for some extensions of the Standard Model in particle physics. Most of the previously conducted searches for FCPs in cosmic rays were based on experiments underground or at high altitudes. However, there have been few searches for FCPs in cosmic rays carried out in orbit other than AMS-01 flown by a space shuttle and BESS by a balloon at the top of the atmosphere. In this study, we conduct an FCP search in space based on on-orbit data obtained using the DArk Matter Particle Explorer (DAMPE) satellite over a period of five years. Unlike underground experiments, which require an FCP energy of the order of hundreds of GeV, our FCP search starts at only a few GeV. An upper limit of $6.2\times 10^{-10}~~\mathrm{cm^{-2}sr^{-1} s^{-1}}$ is obtained for the flux. Our results demonstrate that DAMPE exhibits higher sensitivity than experiments of similar types by three orders of magnitude that more stringently restricts the conditions for the existence of FCP in primary cosmic rays.
△ Less
Submitted 9 September, 2022;
originally announced September 2022.
-
Modeling material transport regulation and traffic jam in neurons using PDE-constrained optimization
Authors:
Angran Li,
Yongjie Jessica Zhang
Abstract:
The intracellular transport process plays an important role in delivering essential materials throughout branched geometries of neurons for their survival and function. Many neurodegenerative diseases have been associated with the disruption of transport. Therefore, it is essential to study how neurons control the transport process to localize materials to necessary locations. Here, we develop a n…
▽ More
The intracellular transport process plays an important role in delivering essential materials throughout branched geometries of neurons for their survival and function. Many neurodegenerative diseases have been associated with the disruption of transport. Therefore, it is essential to study how neurons control the transport process to localize materials to necessary locations. Here, we develop a novel optimization model to simulate the traffic regulation mechanism of material transport in complex geometries of neurons. The transport is controlled to avoid traffic jam of materials by minimizing a pre-defined objective function. The optimization subjects to a set of partial differential equation (PDE) constraints that describe the material transport process based on a macroscopic molecular-motor-assisted transport model of intracellular particles. The proposed PDE-constrained optimization model is solved in complex tree structures by using isogeometric analysis (IGA). Different simulation parameters are used to introduce traffic jams and study how neurons handle the transport issue. Specifically, we successfully model and explain the traffic jam caused by reduced number of microtubules (MTs) and MT swirls. In summary, our model effectively simulates the material transport process in healthy neurons and also explains the formation of a traffic jam in abnormal neurons. Our results demonstrate that both geometry and MT structure play important roles in achieving an optimal transport process in neuron.
△ Less
Submitted 30 August, 2021;
originally announced August 2021.
-
The divergence-conforming immersed boundary method: Application to vesicle and capsule dynamics
Authors:
Hugo Casquero,
Carles Bona-Casas,
Deepesh Toshniwal,
Thomas J. R. Hughes,
Hector Gomez,
Yongjie Jessica Zhang
Abstract:
We extend the recently introduced divergence-conforming immersed boundary (DCIB) method [1] to fluid-structure interaction (FSI) problems involving closed co-dimension one solids. We focus on capsules and vesicles, whose discretization is particularly challenging due to the higher-order derivatives that appear in their formulations. In two-dimensional settings, we employ cubic B-splines with perio…
▽ More
We extend the recently introduced divergence-conforming immersed boundary (DCIB) method [1] to fluid-structure interaction (FSI) problems involving closed co-dimension one solids. We focus on capsules and vesicles, whose discretization is particularly challenging due to the higher-order derivatives that appear in their formulations. In two-dimensional settings, we employ cubic B-splines with periodic knot vectors to obtain discretizations of closed curves with C^2 inter-element continuity. In three-dimensional settings, we use analysis-suitable bi-cubic T-splines to obtain discretizations of closed surfaces with at least C^1 inter-element continuity. Large spurious changes of the fluid volume inside closed co-dimension one solids is a well-known issue for IB methods. The DCIB method results in volume changes orders of magnitude lower than conventional IB methods. This is a byproduct of discretizing the velocity-pressure pair with divergence-conforming B-splines, which lead to negligible incompressibility errors at the Eulerian level. The higher inter-element continuity of divergence-conforming B-splines is also crucial to avoid the quadrature/interpolation errors of IB methods becoming the dominant discretization error. Benchmark and application problems of vesicle and capsule dynamics are solved, including mesh-independence studies and comparisons with other numerical methods.
△ Less
Submitted 22 January, 2020;
originally announced January 2020.
-
Material characterization and precise finite element analysis of fiber reinforced thermoplastic composites for 4D printing
Authors:
Yuxuan Yu,
Haolin Liu,
Kuanren Qian,
Humphrey Yang,
Matthew McGehee,
Jianzhe Gu,
Danli Luo,
Lining Yao,
Yongjie Jessica Zhang
Abstract:
Four-dimensional (4D) printing, a new technology emerged from additive manufacturing (3D printing), is widely known for its capability of programming post-fabrication shape-changing into artifacts. Fused deposition modeling (FDM)-based 4D printing, in particular, uses thermoplastics to produce artifacts and requires computational analysis to assist the design processes of complex geometries. Howev…
▽ More
Four-dimensional (4D) printing, a new technology emerged from additive manufacturing (3D printing), is widely known for its capability of programming post-fabrication shape-changing into artifacts. Fused deposition modeling (FDM)-based 4D printing, in particular, uses thermoplastics to produce artifacts and requires computational analysis to assist the design processes of complex geometries. However, these artifacts are weak against structural loads, and the design quality can be limited by less accurate material models and numerical simulations. To address these issues, this paper propounds a composite structure design made of two materials - polylactic acid (PLA) and carbon fiber reinforced PLA (CFPLA) - to increase the structural strength of 4D printed artifacts and a workflow composed of several physical experiments and series of dynamic mechanical analysis (DMA) to characterize materials. We apply this workflow to 3D printed samples fabricated with different printed parameters to accurately characterize the materials and implement a sequential finite element analysis (FEA) to achieve accurate simulations. The accuracy of deformation induced by the triggering process is both computationally and experimentally verified with several creative design examples, and the 95% confidence interval of the accuracy is (0.972, 0.985). We believe the presented workflow is essential to the combination of geometry, material mechanism and design, and has various potential applications.
△ Less
Submitted 6 November, 2020; v1 submitted 4 September, 2019;
originally announced September 2019.
-
Instrument Performance and Simulation Verification of the POLAR Detector
Authors:
M. Kole,
Z. H. Li,
N. Produit,
T. Tymieniecka,
J. Zhang,
A. Zwolinska,
T. W. Bao,
T. Bernasconi,
F. Cadoux,
M. Z. Feng,
N. Gauvin,
W. Hajdas,
S. W. Kong,
H. C. Li,
L. Li,
X. Liu,
R. Marcinkowski,
S. Orsi,
M. Pohl,
D. Rybka,
J. C. Sun,
L. M. Song,
J. Szabelski,
R. J. Wang,
Y. H. Wang
, et al. (10 additional authors not shown)
Abstract:
POLAR is a new satellite-born detector aiming to measure the polarization of an unprecedented number of Gamma-Ray Bursts in the 50-500 keV energy range. The instrument, launched on-board the Tiangong-2 Chinese Space lab on the 15th of September 2016, is designed to measure the polarization of the hard X-ray flux by measuring the distribution of the azimuthal scattering angles of the incoming photo…
▽ More
POLAR is a new satellite-born detector aiming to measure the polarization of an unprecedented number of Gamma-Ray Bursts in the 50-500 keV energy range. The instrument, launched on-board the Tiangong-2 Chinese Space lab on the 15th of September 2016, is designed to measure the polarization of the hard X-ray flux by measuring the distribution of the azimuthal scattering angles of the incoming photons. A detailed understanding of the polarimeter and specifically of the systematic effects induced by the instrument's non-uniformity are required for this purpose. In order to study the instrument's response to polarization, POLAR underwent a beam test at the European Synchrotron Radiation Facility in France. In this paper both the beam test and the instrument performance will be described. This is followed by an overview of the Monte Carlo simulation tools developed for the instrument. Finally a comparison of the measured and simulated instrument performance will be provided and the instrument response to polarization will be presented.
△ Less
Submitted 2 August, 2017;
originally announced August 2017.
-
The DArk Matter Particle Explorer mission
Authors:
J. Chang,
G. Ambrosi,
Q. An,
R. Asfandiyarov,
P. Azzarello,
P. Bernardini,
B. Bertucci,
M. S. Cai,
M. Caragiulo,
D. Y. Chen,
H. F. Chen,
J. L. Chen,
W. Chen,
M. Y. Cui,
T. S. Cui,
A. D'Amone,
A. De Benedittis,
I. De Mitri,
M. Di Santo,
J. N. Dong,
T. K. Dong,
Y. F. Dong,
Z. X. Dong,
G. Donvito,
D. Droz
, et al. (139 additional authors not shown)
Abstract:
The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives…
▽ More
The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.
△ Less
Submitted 14 September, 2017; v1 submitted 26 June, 2017;
originally announced June 2017.
-
Molecular Dynamics Simulations for Anisotropic Thermal Conductivity of Borophene
Authors:
Yue Jia,
Chun Li,
Jin-Wu Jiang,
Ning Wei,
Yang Chen,
Yongjie Jessica Zhang
Abstract:
The present work carries out molecular dynamics simulations to compute the thermal conductivity of the borophene nanoribbon and the borophene nanotube using the Muller-Plathe approach. We investigate the thermal conductivity of the armchair and zigzag borophenes, and show the strong anisotropic thermal conductivity property of borophene. We compare the results of the borophene nanoribbon and the b…
▽ More
The present work carries out molecular dynamics simulations to compute the thermal conductivity of the borophene nanoribbon and the borophene nanotube using the Muller-Plathe approach. We investigate the thermal conductivity of the armchair and zigzag borophenes, and show the strong anisotropic thermal conductivity property of borophene. We compare the results of the borophene nanoribbon and the borophene nanotube, and find the thermal conductivity of the borophene is structure dependent.
△ Less
Submitted 31 May, 2017;
originally announced May 2017.
-
Calibration of the Space-borne Compton Polarimeter POLAR flight model with 100% polarized X-ray beams
Authors:
H. L. Xiao,
W. Hajdas,
P. Socha,
R. Marcinkowski,
B. B. Wu,
T. W. Bao,
J. Y. Chai,
Y. W. Dong,
M. N. Kong,
L. Li,
Z. H. Li,
J. T. Liu,
H. L. Shi,
L. M. Song,
J. C. Sun,
R. J. Wang,
Y. H. Wang,
X. Wen,
S. L. Xiong,
J. Zhang,
L. Y. Zhang,
S. N. Zhang,
X. F. Zhang,
Y. J. Zhang,
F. Cadoux
, et al. (10 additional authors not shown)
Abstract:
POLAR is space-borne detector designed for a precise measurement of gamma-ray polarization of the prompt emissions of Gamma-Ray Bursts in the energy range 50 keV - 500 keV. POLAR is a compact Compton polarimeter consisting of 40$\times$ 40 plastic scintillator bars read out by 25 multi-anode PMTs. In May 2015, we performed a series of tests of the POLAR flight model with 100\% polarized x-rays bea…
▽ More
POLAR is space-borne detector designed for a precise measurement of gamma-ray polarization of the prompt emissions of Gamma-Ray Bursts in the energy range 50 keV - 500 keV. POLAR is a compact Compton polarimeter consisting of 40$\times$ 40 plastic scintillator bars read out by 25 multi-anode PMTs. In May 2015, we performed a series of tests of the POLAR flight model with 100\% polarized x-rays beams at the European Synchrotron Radiation Facility beam-line ID11 aming to study thresholds, crosstalk between channels and responses of POLAR flight model to polarized X-ray beams. In this paper we present the data analysis method and some analysis results. According the results, POLAR FM has good polarimetric capabilities.
△ Less
Submitted 24 April, 2017; v1 submitted 20 April, 2017;
originally announced April 2017.
-
Gain factor and parameter settings optimization of the new gamma-ray burst polarimeter POLAR
Authors:
X. F. Zhang,
W. Hajdas,
H. L. Xiao,
X. Wen,
B. B. Wu,
T. W. Bao,
T. Batsch,
T. Bernasconi,
F. Cadoux,
I. Cernuda,
J. Y. Chai,
Y. W. Dong,
N. Gauvin,
J. J. He,
M. Kole,
M. N. Kong,
C. Lechanoine-Leluc,
L. Li,
Z. H. Li,
J. T. Liu,
X. Liu,
R. Marcinkowski,
S. Orsi,
M. Pohl,
D. Rapin
, et al. (16 additional authors not shown)
Abstract:
As a space-borne detector POLAR is designed to conduct hard X-ray polarization measurements of gamma-ray bursts on the statistically significant sample of events and with an unprecedented accuracy. During its development phase a number of tests, calibrations runs and verification measurements were carried out in order to validate instrument functionality and optimize operational parameters. In thi…
▽ More
As a space-borne detector POLAR is designed to conduct hard X-ray polarization measurements of gamma-ray bursts on the statistically significant sample of events and with an unprecedented accuracy. During its development phase a number of tests, calibrations runs and verification measurements were carried out in order to validate instrument functionality and optimize operational parameters. In this article we present results on gain optimization togeter with verification data obtained in the course of broad laboratory and environmental tests. In particular we focus on exposures to the $^{137}$Cs radioactive source and determination of the gain dependence on the high voltage for all 1600 detection channels of the polarimeter. Performance of the instrument is described in detail with respect to the dynamic range, energy resolution and temperature dependence. Gain optimization algorithms and response non-uniformity studies are also broadly discussed. Results presented below constitute important parts for development of the POLAR calibration and operation database.
△ Less
Submitted 14 March, 2017; v1 submitted 12 March, 2017;
originally announced March 2017.
-
Calibration of Gamma-ray Burst Polarimeter POLAR
Authors:
H. L. Xiao,
W. Hajdas,
T. W. Bao,
T. Batsch,
T. Bernasconi,
I. Cernuda,
J. Y. Chai,
Y. W. Dong,
N. Gauvin,
M. Kole,
M. N. Kong,
S. W. Kong,
L. Li,
J. T. Liu,
X. Liu,
R. Marcinkowski,
S. Orsi,
M. Pohl,
N. Produit,
D. Rapin,
A. Rutczynska,
D. Rybka,
H. L. Shi,
L. M. Song,
J. C. Sun
, et al. (11 additional authors not shown)
Abstract:
Gamma Ray Bursts (GRBs) are the strongest explosions in the universe which might be associated with creation of black holes. Magnetic field structure and burst dynamics may influence polarization of the emitted gamma-rays. Precise polarization detection can be an ultimate tool to unveil the true GRB mechanism. POLAR is a space-borne Compton scattering detector for precise measurements of the GRB p…
▽ More
Gamma Ray Bursts (GRBs) are the strongest explosions in the universe which might be associated with creation of black holes. Magnetic field structure and burst dynamics may influence polarization of the emitted gamma-rays. Precise polarization detection can be an ultimate tool to unveil the true GRB mechanism. POLAR is a space-borne Compton scattering detector for precise measurements of the GRB polarization. It consists of a 40$\times$40 array of plastic scintillator bars read out by 25 multi-anode PMTs (MaPMTs). It is scheduled to be launched into space in 2016 onboard of the Chinese space laboratory TG2. We present a dedicated methodology for POLAR calibration and some calibration results based on the combined use of the laboratory radioactive sources and polarized X-ray beams from the European Synchrotron Radiation Facility. They include calibration of the energy response, computation of the energy conversion factor vs. high voltage as well as determination of the threshold values, crosstalk contributions and polarization modulation factors.
△ Less
Submitted 9 December, 2015;
originally announced December 2015.