Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data processing and standardization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48

Physics-Informed Generative Adversarial Networks for Stochastic Differential Equations

We developed a new class of physics-informed generative adversarial networks (PI-GANs) to solve forward, inverse, and mixed stochastic problems in a unified manner based on a limited number of scattered measurements. Unlike standard GANs relying solely on data for training, here we encoded into the architecture of GANs the governing physical laws in the form of stochastic differential equations (SDEs) using automatic differentiation. In particular, we applied Wasserstein GANs with gradient penalty (WGAN-GP) for its enhanced stability compared to vanilla GANs. We first tested WGAN-GP in approximating Gaussian processes of different correlation lengths based on data realizations collected from simultaneous reads at sparsely placed sensors. We obtained good approximation of the generated stochastic processes to the target ones even if there is a mismatch between the input noise dimensionality and the effective dimensionality of the target stochastic processes. We also studied the overfitting issue for both the discriminator and the generator, and we found that overfitting occurs also in the generator in addition to the discriminator as previously reported. Subsequently, we considered the solution of elliptic SDEs requiring approximations of three stochastic processes, namely the solution, the forcing, and the diffusion coefficient. Here again, we assumed data collected from simultaneous reads at a limited number of sensors for the multiple stochastic processes. Furthermore, three generators were used for the PI-GANs: two of them were feed forward deep neural networks (DNNs), while the other one was the neural network induced by the SDE. For the case where we have one group of data, we employed one feed forward DNN as the discriminator, while for the case of multiple groups of data we employed multiple discriminators in PI-GANs. We solved forward, inverse, and mixed problems without changing the framework of PI-GANs, obtaining both the means and the standard deviations of the stochastic solution and the diffusion coefficient in good agreement with benchmarks. In this work, we have demonstrated the effectiveness of PI-GANs in solving SDEs for about 120 dimensions. In principle, PI-GANs could tackle very high dimensional problems given more sensor data with low-polynomial growth in computational cost.

97 MATHEMATICS AND COMPUTING↗

Three-dimensional anisotropic adaptive filtering of projection data for noise reduction in cone beam CT

Purpose: The combination of quickly rotating C-arm gantry with digital flat panel has enabled the acquisition of three-dimensional data (3D) in the interventional suite. However, image quality is still somewhat limited since the hardware has not been optimized for CT imaging. Adaptive anisotropic filtering has the ability to improve image quality by reducing the noise level and therewith the radiation dose without introducing noticeable blurring. By applying the filtering prior to 3D reconstruction, noise-induced streak artifacts are reduced as compared to processing in the image domain. Methods: 3D anisotropic adaptive filtering was used to process an ensemble of 2D x-ray views acquired along a circular trajectory around an object. After arranging the input data into a 3D space (2D projections + angle), the orientation of structures was estimated using a set of differently oriented filters. The resulting tensor representation of local orientation was utilized to control the anisotropic filtering. Low-pass filtering is applied only along structures to maintain high spatial frequency components perpendicular to these. The evaluation of the proposed algorithm includes numerical simulations, phantom experiments, and in-vivo data which were acquired using an AXIOM Artis dTA C-arm system (Siemens AG, Healthcare Sector, Forchheim, Germany). Spatial resolution and noise levels were compared with and without adaptive filtering. A human observer study was carried out to evaluate low-contrast detectability. Results: The adaptive anisotropic filtering algorithm was found to significantly improve low-contrast detectability by reducing the noise level by half (reduction of the standard deviation in certain areas from 74 to 30 HU). Virtually no degradation of high contrast spatial resolution was observed in the modulation transfer function (MTF) analysis. Although the algorithm is computationally intensive, hardware acceleration using Nvidia's CUDA Interface provided an 8.9-fold speed-up of the processing (from 1336 to 150 s). Conclusions: Adaptive anisotropic filtering has the potential to substantially improve image quality and/or reduce the radiation dose required for obtaining 3D image data using cone beam CT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Novel Gamma Attenuation Methods Utilized for Waste Characterization of UO{sub 2}F{sub 2}-Bearing Materials at the Paducah GDP - 17327

The Paducah Gaseous Diffusion Plant (PGDP) is currently under deactivation with demolition in the decades to follow. UO{sub 2}F{sub 2}-bearing waste is generated at PGDP by UF{sub 6} hydrolysing to materials such as pipes, valves, PPE, etc. when UF{sub 6} is exposed to humid air. In addition, Al{sub 2}O{sub 3} and NaF trap materials contain significant amounts of Uranium. The UO{sub 2}F{sub 2} deposits within the generated waste can be of sub-gram levels to several hundreds of localized grams and are much denser than the surrounding interfering matrix materials and may be of very complex geometric and chemical forms such as with the Al{sub 2}O{sub 3} and NaF trap materials. The main objective is to characterize the fissile (U-235) content of the UO{sub 2}F{sub 2}-bearing waste so as to meet the Data Quality Objectives (DQO) for proper disposal. If a waste container cannot meet the DQO then it must be subjected to costly remediation, repackaging and re-characterization processes. We have developed, at PGDP, several novel model-based methods and techniques that vastly improve the standard gamma attenuation correction methods and do not require Working Reference Materials (WRM) for calibration. These new methods involve rotating and non-rotating 18.9 liter (5 U.S. gallon) waste containers with and without a transmission modality. A new model-based variant of differential attenuation has been developed that provides a high-quality self-absorption correction that does not require knowledge of the Uranium enrichment or isotopic ratios. For more challenging waste containers, a method using multi-view (non-rotating) 18.9 liter waste containers has been considered using a recently developed Advanced ISOCS{sup C} Uncertainty Estimator (AIUE) package. The new methods suffice to significantly improve accuracy and precision and to lower the Total Measurement Uncertainty (TMU) leading to satisfactory DQO's that were not previously achievable using the standard gamma attenuation correction methods. In this paper we present real 18.9 liter waste container and Working Reference Material (WRM) data subjected to the new methods and compared to standard methods depicting the significant improvement in the DQO that is realized. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

3d Wave Propagation and Site Effects in the Humboldt Bay Area Using Strong Ground Motion Records from the M6.5 2010 Ferndale Earthquake

In this study we simulate and analyze strong ground motion waveform data recorded in the Humboldt Bay and Eureka areas during the M6.5 Ferndale earthquake of January 2010. The scope of the presented work was two-fold. First, we investigated the main aspects of seismic wave generation and propagation, including kinematic rupture process and 3D wave propagation scattering. Our goal is to analyze their potential effects on seismic motion recorded at free field stations across Humboldt Bay and Eureka, California, and test the performance of a standard broadband strong ground motion simulation technique. Second, we evaluated the implication of the non-linear response of shallow sediments on ground motion recorded at the Humboldt Bay geotechnical array. Using synthetic and recorded data we compared the performance of numerical and empirical techniques in simulating local non-linear site response. Our study provides insight into the composition of the wave field during the earthquake and an improved understanding of how the strong ground motion is affected by the 3D structure and non-linear response of shallow sedimentary layers in the Humboldt Bay.

58 GEOSCIENCES↗

Standard Ruler at Cosmic Dawn

The matter in our Universe comes in two flavors: dark and baryonic. Of these, only the latter couples to photons, giving rise to the well-known baryon acoustic oscillations and, in the process, generating supersonic relative velocities between dark matter and baryons. These velocities—imprinted with the acoustic scale in their genesis—impede the formation of the first stars during cosmic dawn (z ~ 20), modulating the expected 21-cm signal from this era. In a companion paper we showed, combining numerical simulations and analytic models, that this modulation takes the form of robust velocity-induced acoustic oscillations (VAOs), with a well-understood shape that is frozen at recombination, and unaffected by the unknown astrophysics of star formation. Here we propose using these VAOs as a standard ruler at cosmic dawn. We find that three years of 21-cm power-spectrum data from the upcoming HERA interferometer should be able to measure the Hubble expansion rate H(z) at z = 15–20 to percent-level precision, ranging from 0.3% to 11% depending on the strength of astrophysical feedback processes and foregrounds. Furthermore, this would provide a new handle on the expansion rate of our Universe during an otherwise unprobed epoch, opening a window to the mysterious cosmic-dawn era.

79 ASTRONOMY AND ASTROPHYSICS↗

Unpacking model inadequacy: The quantification of silver release from TRISO fuel by considering empirical and mechanistic approaches

Increasing adoption of the proposed tristructural isotropic (TRISO) particle fuel for both advanced and existing reactors makes it critical to assess and address any uncertainties and inadequacies of TRISO fission product release models. Model inadequacy stems from simplifications made to the computational model when compared to the experiments. The modeling and simulation efforts conducted using the BISON fuel performance code, along with the experimental campaigns carried out under the Advanced Gas Reactor Fuel Development and Qualification Program, afford a unique opportunity to conduct a rigorous modeling inadequacy assessment within the Bayesian uncertainty quantification (UQ) framework. Here, this study compares the standard Bayesian framework against the Kennedy-O'Hagan (KOH) framework, which explicitly represents modeling inadequacy, in regard to UQ for TRISO silver release models. For this purpose, both the traditional Arrhenius equation fitted to experimental data and the more advanced lower-length-scale (LLS)-informed model, which considers microstructure information, are independently considered. Applying the inverse UQ process on the AGR-2 and -3/4 datasets revealed modeling inadequacy to be the most dominant source of uncertainty. Experimental noise uncertainty is also significant; however, model parameter uncertainty can be considered negligible. Interestingly, both the Arrhenius equation and the LLS-informed model demonstrated similar levels of modeling inadequacy. For the forward predictive UQ, the KOH framework improved both the accuracy and quality of quantified uncertainties in comparison to the standard Bayesian framework. This is true for both the Arrhenius equation and the LLS-informed model. In comparing these modeling approaches, both demonstrated similar performance at the engineering scale, while the LLS-informed model expectedly outperformed the Arrhenius equation at the mesoscale. These conclusions highlight the importance of explicitly accounting for modeling inadequacy in the UQ process, and reinforce the need for continuous refinement of physics-based models in order to address the modeling inadequacy.

Advanced reactors↗

A deep neural network to search for new long-lived particles decaying to jets

A tagging algorithm to identify jets that are significantly displaced from the proton-proton (pp) collision region in the CMS detector at the LHC is presented. Displaced jets can arise from the decays of long-lived particles (LLPs), which are predicted by several theoretical extensions of the standard model. The tagger is a multiclass classifier based on a deep neural network, which is parameterised according to the proper decay length $\mathrm{c}\tau_0$ of the LLP. A novel scheme is defined to reliably label jets from LLP decays for supervised learning. Samples of pp collision data, recorded by the CMS detector at a centre-of-mass energy of 13 TeV, and simulated events are used to train the neural network. Domain adaptation by backward propagation is performed to improve the simulation modelling of the jet class probability distributions observed in pp collision data. The potential performance of the tagger is demonstrated with a search for long-lived gluinos, a manifestation of split supersymmetric models. The tagger provides a rejection factor of 10 000 for jets from standard model processes, while maintaining an LLP jet tagging efficiency of 30-80% for gluinos with 1 mm $\leq$ $c\tau_0$ $\leq$ 10 m. The expected coverage of the parameter space for split supersymmetry is presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Search for microscopic black holes in a like-sign dimuon final state using large track multiplicity with the ATLAS detector

A search is presented for microscopic black holes in a like-sign dimuon final state in proton-proton collisions at $ \sqrt{s}$=8 TeV. The data were collected with the ATLAS detector at the Large Hadron Collider in 2012 and correspond to an integrated luminosity of 20.3 fb -1 . Using a high track multiplicity requirement, 0.6±0.2 background events from Standard Model processes are predicted and none observed. This result is interpreted in the context of low-scale gravity models and 95% C.L. lower limits on microscopic black hole masses are set for different model assumptions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measurement of the direct C P -violating parameter A C P in the decay D + → K - π + π +

We measure the direct C P -violating parameter A C P for the decay of the charged charm meson, D + → K - π + π + (and charge conjugate), using the full 10.4 fb - 1 sample of p p ¯ collisions at s = 1.96 TeV collected by the D0 detector at the Fermilab Tevatron collider. We extract the raw reconstructed charge asymmetry by fitting the invariant mass distributions for the sum and difference of charge-specific samples. This quantity is then corrected for detector-related asymmetries using data-driven methods and for possible physics asymmetries (from B → D processes) using input from Monte Carlo simulation. We measure A C P = [ - 0.16 ± 0.15 ( stat ) ± 0.09 ( syst ) ] % , which is consistent with zero, as expected from the standard model prediction of C P conservation, and is the most precise measurement of this quantity to date.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effective Tooling for Linked Data Publishing in Scientific Research

Challenges that make it difficult to find, share, and combine published data, such as data heterogeneity and resource discovery, have led to increased adoption of semantic data standards and data publishing technologies. To make data more accessible, interconnected and discoverable, some domains are being encouraged to publish their data as Linked Data. Consequently, this trend greatly increases the amount of data that semantic web tools are required to process, store, and interconnect. In attempting to process and manipulate large data sets, tools–ranging from simple text editors to modern triplestores– eventually breakdown upon reaching undefined thresholds. This paper offers a systematic approach that data publishers can use to categorize suitable tools to meet their data publishing needs. We present a real-world use case, the Resource Discovery for Extreme Scale Collaboration (RDESC), which features a scientific dataset(maximum size of 1.4 billion triples) used to evaluate a toolbox for data publishing in climate research. This paper also introduces a semantic data publishing software suite developed for the RDESC project.

Purohit, Sumit↗

FAIRification, Quality Assessment, and Missingness Pattern Discovery for Spatiotemporal Photovoltaic Data

The growth of the photovoltaic market has pushed the demand for power forecasting and performance evaluation for a huge population of PV power plants. Many of these power plants have spatiotemporal coherence that can be utilized for improving model accuracy. We have demonstrated in this paper the FAIRification of spatiotemporal PV time series data. Through the creation of a solar power plant ontology, we propose standards for the naming and structure of metadata used to describe the data from these power plants. Using the structure from this ontology, we have developed both R and Python packages for the automation of the FAIRification process. Going further, we have also developed an R package that automates the analysis of the quality of a data set through the designation of letter grades. To solve the issue of data missingness, we propose the use of St-GNN autoencoders to detect and impute missing values from a data set by utilizing data from power plants nearby.

14 SOLAR ENERGY↗

A Recommended Pasquill-Gifford Stability Classification Method for Safety Basis Atmospheric Dispersion Modeling at SRS

Several of the most common methods for estimating Pasquill-Gifford (PG) stability (turbulence) class were evaluated for use in modeling the radiological consequences of SRS accidental releases using the MELCOR Accident Consequence Code System, Ver. 2 (MACCS2). Evaluation criteria included: (1) the ability of the method to represent diffusion characteristics above a predominantly forested landscape at SRS, (2) suitability of the method to provide data consistent with the formulation of the MACCS2 model, and (3) the availability of onsite meteorological data to support implementation of the method The evaluation resulted in a recommendation that PG stability classification for regulatory applications at SRS should be based on measurements of the standard deviation of the vertical component of wind direction fluctuations, σ e , collected from the 61-m level of the SRS meteorological towers, and processed in full accordance with EPA-454/R-99-005 (EPA, 2000). This approach provides a direct measurement that is fundamental to diffusion and captures explicitly the turbulence generated by both mechanical and buoyant forces over the characteristic surface (forested) of SRS. Furthermore, due to the potentially significant enhancement of horizontal fluctuations in wind direction from the occurrence of meander at night, the use of σ e will ensure a reasonably conservative estimate of PG stability class for use in dispersion models that base diffusion calculations on a single value of PG stability class. Furthermore, meteorological data bases used as input for MACCS2 calculations should contain hourly data for five consecutive annual periods from the most recent 10 years.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Analysis of Data Fusion Between Waveform Events and Radionuclide Detections Reported in 2021 by the International Data Centre

During this project, we assessed the percentage of events reported by the International Data Centre (IDC) that are candidates for evaluation using data fusion methods. We calculated the percentage of waveform (seismic, hydroacoustic, infrasound, SHI) events that may be associated with radionuclide (RN) detections and vice versa. For this task, data fusion means identifying RN detections that may be associated with SHI events. The task of producing a fused SHI and RN bulletin is included in the IDC operational manual. As of 2022, the production of a fused event bulletin remains incomplete, largely because methods for producing the bulletin are under development. Nonetheless, simple methods for data fusion have been developed and preliminary software tools are available at the IDC to associate SHI events and RN detections. This report is a first look at the percentage of SHI events that may be candidates for fusion with RN detections or bulletins. One year of IDC bulletin data (2021) is considered in this report. It is important to note that there were no known nuclear explosions during this time period. Hence, we expect that events during this time period would generate SHI and RN signals. If many SHI events are candidates for fusion with RN data, greater priority should be placed on establishing screening procedures before and/or during the fusion process. By contrast, if the number of candidate events is diminishingly small, then data fusion is not a driver for further development of event or detection screening at the IDC. LLNL reviewed the events reported in the Standard Event Bulletin (SEB) and determined which events are candidates for fusion with the Reviewed Radionuclide Reports (RRR). Much of the work involves establishing a metric or metrics for determining whether a Reviewed Event Bulletin (REB) event is a candidate for data fusion.

58 GEOSCIENCES↗

KS4A-Omics1.0_FspDS682

Soil fungi facilitate the translocation of inorganic nutrients from soil minerals to other microorganisms and plants. This ability is particularly advantageous in impoverished soils, because fungal mycelial networks can bridge otherwise spatially disconnected and inaccessible nutrient hotspots. However, the molecular mechanisms underlying fungal mineral weathering and transport through soil remains poorly understood. Here, we addressed this knowledge gap by directly visualizing nutrient acquisition and transport through fungal hyphae in a mineral doped soil micromodel using a multimodal imaging approach. Here, we observed how a representative of common saprotrophic soil fungi, Fusarium sp. DS 682, exhibited a mechanosensory response (thigmotropism) around obstacles and through pore spaces (~12 μm) in the presence of minerals.This study establishes the significance of fungal biology and nutrient translocation mechanisms in maintaining fungal growth under water and nutrient limitations in a soil-like microenvironment, using a high-throughput multi-omic analysis approach. Data package KS4A-Omics.1.0_FspDS682 (Publication: Fungal Mineral Weathering Mechanisms Revealed Through Direct Molecular Visualization) contents reported here are the first version (1.0) and contain pre- and post-processed data using high throughput data capture technologies for multi-omic analysis and integration, this data package contains raw and post-processed experimental data for X-Ray Absorption Near Edge Structure Spectroscopy (XANES/XRF), Optical Microscopy, Proteomics, Scanning Electron Microscope (SEM), Time-of-Flight Secondary Ion Mass Spectroscopy (ToF-SIMS), X-Ray Diffraction Spectroscopy (XRD) files, and X- Ray Photoelectron Spectroscopy (XPS) using EMSL capabilities. This data package DOI contains a comprehensive collection of high-throughput multi-omics data and process metadata catalog. Support files include additional data download contents “Read Me” with dataset descriptor information and data source method application ontologies (see data dictionary section). Reported data download content is structured for compliance with reported guidelines provided by community standard initiatives and publisher stakeholder policies supporting FAIR data principles.

47 OTHER INSTRUMENTATION↗

The R -process Alliance: A Nearly Complete R -process Abundance Template Derived from Ultraviolet Spectroscopy of the R -process-enhanced Metal-poor Star HD 222925 *

We present a nearly complete rapid neutron-capture process (r-process) chemical inventory of the metal-poor ([Fe/H] = -1.46 ± 0.10) r-process-enhanced ([Eu/Fe] = +1.32 ± 0.08) halo star HD 222925. This abundance set is the most complete for any object beyond the solar system, with a total of 63 metals detected and seven with upper limits. It comprises 42 elements from 31 ≤ Z ≤ 90, including elements rarely detected in r-process-enhanced stars, such as Ga, Ge, As, Se, Cd, In, Sn, Sb, Te, W, Re, Os, Ir, Pt, and Au. We derive these abundances from an analysis of 404 absorption lines in ultraviolet spectra collected using the Space Telescope Imaging Spectrograph on the Hubble Space Telescope and previously analyzed optical spectra. A series of appendices discusses the atomic data and quality of fits for these lines. The r-process elements from Ba to Pb, including all elements at the third r-process peak, exhibit remarkable agreement with the solar r-process residuals, with a standard deviation of the differences of only 0.08 dex (17%). In contrast, deviations among the lighter elements from Ga to Te span nearly 1.4 dex, and they show distinct trends from Ga to Se, Nb through Cd, and In through Te. The r-process contribution to Ga, Ge, and As is small, and Se is the lightest element whose production is dominated by the r-process. The lanthanide fraction, log X La = -1.39 ± 0.09, is typical for r-process-enhanced stars and higher than that of the kilonova from the GW170817 neutron-star merger event. We advocate adopting this pattern as an alternative to the solar r-process-element residuals when confronting future theoretical models of heavy-element nucleosynthesis with observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Muon Campus at Fermilab

The Muon Campus at Fermilab provides world class accelerator infrastructure supporting the next generation intensity frontier experiments. The anti-proton source from the Tevatron era was converted to the present day Muon Campus at the end of the collider program in 2011. Currently, the Muon Campus delivers highly polarized muon beams to the Muon g-2 experiment and in a few years, it will provide muon beams to the Mu2e experiment. The Muon g-2 experiment is a high precision test of the standard model; it is precise enough to measure QED, Weak, and QCD Standard Model contributions.The Muon g-2 experiment looks for a significant deviation from the Standard Model whereas the Mu2e experiment looks for rare processes using high intensity muon beams. The Muon g-2 experiment has measured the anomalous magnetic moment of muons to an unprecedented 460 ppb precision with Run 1 data. The experiment has completed 5 runs, accumulating 19$\times$ Brookhaven E821 data set. The Mu2e experiment has completed the construction phase and first beam to the Diagnostic Absorber was run on April 14, 2022. Installation and commissioning of the Mu2e experiment will begin in 2023 and beam commissioning and physics run is currently planned for 2026.

43 PARTICLE ACCELERATORS↗

CHANDRA PHASE-RESOLVED X-RAY SPECTROSCOPY OF THE CRAB PULSAR

We present a new study of the X-ray spectral properties of the Crab Pulsar. The superb angular resolution of the Chandra X-Ray Observatory enables distinguishing the pulsar from the surrounding nebulosity. Analysis of the spectrum as a function of pulse phase allows the least-biased measure of interstellar X-ray extinction due primarily to photoelectric absorption and secondarily to scattering by dust grains in the direction of the Crab Nebula. We modify previous findings that the line of sight to the Crab is underabundant in oxygen and provide measurements with improved accuracy and less bias. Using the abundances and cross sections from Wilms et al. we find [O/H] = (5.28 {+-} 0.28) Multiplication-Sign 10{sup -4} (4.9 Multiplication-Sign 10{sup -4} is solar abundance). We also measure for the first time the impact of scattering of flux out of the image by interstellar grains. We find {tau}{sub scat} = 0.147 {+-} 0.043. Analysis of the spectrum as a function of pulse phase also measures the X-ray spectral index even at pulse minimum-albeit with increasing statistical uncertainty. The spectral variations are, by and large, consistent with a sinusoidal variation. The only significant variation from the sinusoid occurs over the same phase range as some rather abrupt behavior in the optical polarization magnitude and position angle. We also compare these spectral variations to those observed in gamma-rays and conclude that our measurements are both a challenge and a guide to future modeling and will thus eventually help us understand pair cascade processes in pulsar magnetospheres. The data are also used to set new, and less biased, upper limits to the surface temperature of the neutron star for different models of the neutron star atmosphere. We discuss how such data are best connected to theoretical models of neutron star cooling and neutron star interiors. The data restrict the neutrino emission rate in the pulsar core and the amount of light elements in the heat-blanketing envelope. The observations allow the pulsar, irrespective of the composition of its envelope, to have a neutrino emission rate higher than {approx}0.2 of the standard rate of a non-superfluid star cooling via the modified Urca process. The observations also allow the rate to be lower but now with a limited amount of accreted matter in the envelope.

79 ASTRONOMY AND ASTROPHYSICS↗

A framework for data-driven solution and parameter estimation of PDEs using conditional generative adversarial networks

We employ and adapt the image-to-image translation concept based on conditional generative adversarial networks (cGAN) for learning a forward and an inverse solution operator of partial differential equations (PDEs). We focus on steady-state solutions of coupled hydromechanical processes in heterogeneous porous media and present the parameterization of the spatially heterogeneous coefficients, which is exceedingly difficult using standard reduced-order modeling techniques. We show that our framework provides a speed-up of at least 2,000 times compared to a finite-element solver and achieves a relative root-mean-square error (r.m.s.e.) of less than 2% for forward modeling. For inverse modeling, the framework estimates the heterogeneous coefficients, given an input of pressure and/or displacement fields, with a relative r.m.s.e. of less than 7%, even for cases where the input data are incomplete and contaminated by noise. The framework also provides a speed-up of 120,000 times compared to a Gaussian prior-based inverse modeling approach while also delivering more accurate results.

97 MATHEMATICS AND COMPUTING↗