Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter sampling”

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 1,333 records · Page 74

DESI peculiar velocity survey – Fundamental Plane

The Dark Energy Spectroscopic Instrument (DESI) peculiar velocity survey aims to measure the peculiar velocities of early- and late-type galaxies within the DESI footprint using both the Fundamental Plane and optical Tully–Fisher relations. Direct measurements of peculiar velocities can significantly improve constraints on the growth rate of structure, reducing uncertainty by a factor of approximately 2.5 at redshift 0.1 compared to the DESI Bright Galaxy Survey’s redshift space distortion measurements alone. We assess the quality of stellar velocity dispersion measurements from DESI spectroscopic data. These measurements, along with photometric data from the Legacy Survey, establish the Fundamental Plane relation and determine distances and peculiar velocities of early-type galaxies. During survey validation, we obtain spectra for 6698 unique early-type galaxies, up to a photometric redshift of 0.15. 64 per cent of observed galaxies (4267) have relative velocity dispersion errors below 10 per cent. This percentage increases to 75 per cent if we restrict our sample to galaxies with spectroscopic redshifts below 0.1. We use the measured central velocity dispersion, along with photometry from the DESI Legacy Imaging Surveys, to fit the Fundamental Plane parameters using a 3D Gaussian maximum likelihood algorithm that accounts for measurement uncertainties and selection cuts. In addition, we conduct zero-point calibration using the absolute distance measurements to the Coma cluster, leading to a value of the Hubble constant, H 0 = 76.05 ± 0.35 (statistical) ±0.49 (systematic Fundamental Plane) ±4.86 (statistical due to calibration) km s –1 Mpc –1 ⁠. This H 0 value is within 2σ of Planck cosmic microwave background results and within 1σ of other low-redshift distance indicator-based measurements.

cosmological parameters↗

Electron-electron interactions in low-dimensional Si:P delta layers

Scientists have long studied the physics of highly disordered conducting systems, seeking to understand the multitude of quantum phenomena that govern how electrons move through material systems. Recently, research into silicon-based quantum computing has made disordered conducting systems, such as Si:P monolayers embedded in isotopically pure Si, technically relevant. Consequently, applying and advancing the theoretical frameworks developed to describe electron behavior in disordered systems is a necessary objective in this field of research. This study investigates key components of dopant-based Si quantum computing devices: embedded regions of highly doped delta layers (δ layers). We examine the transport behavior and the electron-electron interaction (EEI) physics in embedded Si:P δ layers by means of self-consistent magnetotransport measurements. Parameters associated with the electronic transport offer a meaningful quantitative characterization of δ-layer quality and dopant diffusion. In addition, by examining EEI behaviors in a set of samples with embedded Si:P δ layers produced with different PH3 exposure procedures prior to Si encapsulation, we show how details of material synthesis affect the dimensionality of charge carrier interactions in embedded Si:P δ layers. Furthermore, the relationship between δ-layer confinement and EEI screening lengths is established here. This understanding will help validate important models used for device simulation and design and lead to improvements in the control of electrostatic gating of and tunneling transport through Si:P single atom transistors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Infrared Optical Constants from Pressed Pellets of Powders: I. Improved n and k Values of (NH4)(2)SO(4)from Single-Angle Reflectance

In combination with other parameters, the real, n(?), and imaginary, k(?), components of the complex refractive index, = n + ik, can be used to simulate the optical properties of a material in different forms, e.g. its infrared spectra. But obtaining reliable values of the optical constants n/k for solid materials is challenging due to the lack of optical quality specimens, usually crystals, large enough to measure. An alternative to crystals is to press the powder into a uniform disk. We have produced pellets from ammonium sulfate, (NH4)2SO4, powder and derived the pellets’ n and k values via single-angle reflectance using a specular reflectance device in combination with a Fourier transform infrared spectrometer. The single-angle technique measures amplitude of light reflected from the material as a function of wavelength over a wide spectral domain; the optical constants are determined from the reflectance data using the Kramers-Kronig relationship. We investigate several parameters associated with the pellets and pellet formation and their effects upon delivering the most reliable n/k values. Parameters studied include pellet diameter, mass and density (void space), drying, grinding, sieving and particle size in the pellet formation, as well as pressing pressure and duration. Of these parameters, using size-selected mixtures of dried, small (< 100 µm) particles and pressing at = 10 tons for at least 30 minutes were found key to forming highly reflective samples. Comparison to two sets of previous literature n(?) and k(?) values obtained from crystalline (NH4)2SO4 both as crystal reflectance as well as extinction spectra of aerosols measured in a flow tube show reasonable agreement, but suggests the present values, as confirmed from two independent techniques, represent a substantial improvement for n/k values for (NH4)2SO4, also demonstrating promise to measure the optical constants of other materials.

Optical constants, index of refraction, single-ang↗

Stress-Drop and Source Scaling of the 2019 Ridgecrest, California, Earthquake Sequence

Stress drop, while difficult to measure reliably and at scale, is a key source parameter for understanding the earthquake rupture process and its relationship to strong ground motion. Here, we use a P-wave spectral decomposition approach, designed for large and densely sampled datasets, to measure earthquake stress drop in the region surrounding the 2019 Ridgecrest, California, earthquake sequence. With more than 11,000 measurements of earthquake stress drop in the 20-yr time period from 2000 through 2019, this dataset provides an opportunity to understand how coseismic stress changes and how other geophysical factors relate to the distribution of stress drop and its evolution in space and time. We observe a mild but persistent deviation from self-similar scaling, with larger events having systematically higher stress drops, though this trend depends on the assumption of an omega-square source spectral model. Earthquake stress drop increases with hypocentral depth in this study region, and the Ridgecrest aftershocks tend to have higher stress drops than the pre-event seismicity. This is in part due to their deeper hypocenters. Coherent spatial patterns of stress drop in the aftershock sequence correlate with the slip distribution of the M 7.1 mainshock, whose northwest rupture tip terminated in a long-lived zone of enervated stress drop. Although physical interpretation of these results is complicated by the trade-offs between the timing, depth, and location of these earthquakes, the observations provide new insight into the physics of the earthquake source in an area of renewed seismic activity in southern California.

58 GEOSCIENCES↗

Weighted greedy-optimal design of computer experiments for kernel-based and Gaussian process model emulation and calibration

This article is concerned with the approximation of high-dimensional functions by kernel-based methods. Motivated by uncertainty quantification, which often necessitates the construction of approximations that are accurate with respect to a probability density function of random variables, we aim at minimizing the approximation error with respect to a weighted $L^p$-norm. We present a greedy procedure for designing computer experiments based upon a weighted modification of the pivoted Cholesky factorization. The method successively generates nested samples with the goal of minimizing error in regions of high probability. Numerical experiments validate that this new importance sampling strategy is superior to other sampling approaches, especially when used with non-product probability density functions. We also show how to use the proposed algorithm to efficiently generate surrogates for inferring unknown model parameters from data.

97 MATHEMATICS AND COMPUTING↗

DEVELOPMENT OF A QUANTUM ELECTRON BEAM DIAGNOSTIC APPARATUS

Spatial properties of electron beams are the essential parameters needed in particle physics and accelerator research. Beam characterization using optical detection methods through the interaction between photons and electrons such as Compton scattering, or electro-optic sampling of the electric field of electrons have been extensively explored and successfully implemented as non-invasive diagnostics at various accelerator facilities. However, such methods often suffer from inherently low sensitivity. Here we present the study of a new type of electron beam diagnostic device for direct optical imaging of electron beams at various energy levels. The concept is based on the extremely high sensitivity of atoms, prepared in a specific "dark superposition" quantum state, to the external electric or magnetic perturbations induced by the passing charged particle, which allows atoms to change polarization of the probe light or absorb probe light and to fluoresce, enabling direct 3D imaging of the charged particles with high resolution. We report our recent experiment results and the design effort on a compact apparatus intended to be tested with the relativistic electron beams at Jefferson Laboratory.

Zhang, S.↗

Manufacturing Dissolvable Alloy Components using ShAPE: A Study on Solid Phase Processing of Mg-W Composites and WJ11 Alloys

Dissolvable alloys and composites are an emerging class of materials that demonstrate tailored sorptivity in corrosive environments. This work explored the viability of solid phase processing (SPP) techniques, namely the shear assisted processing and extrusion (ShAPETM) method, to manufacture dissolvable magnesium composites and alloys components such as rods and tubes. The study also determined the effects of material composition and manufacturing process conditions on material microstructures developed during SPP and the resulting performance of the components. Explicit properties of interest included ultimate tensile strength, ultimate compressive strength, and yield strength. The application fields of interest for dissolvable alloys were identified as oil-and-gas operations and biomedical implants. Accordingly, PNNL collaborated with industrial and academia partners, Fortek Industries, Houston and University of Pittsburgh, during the course of this project for material and ShAPE process development. PNNL determined the viability of the using ShAPE to process the custom-designed materials provided by the project partners. This was done by identifying process parameters such as tool rotation rate, tool plunge rate, and process cooling optimal for manufacturing components with consolidated microstructures. Subsequently, PNNL also performed microstructural characterization on unprocessed and processed samples and mechanical property testing. One of the major findings from this work is that ShAPE was able to manufacture dissolvable magnesium composite and alloy rods and tubes with consolidated homogeneous microstructures and minimal macroscale defects on component surfaces. It was observed that samples processed via ShAPE demonstrated an average grain size < 2 - 5 µm, which was markedly lower than that of the precursors with 8 – 13 µm grains. Results showed that the ShAPE synthesized dissolvable magnesium composite and alloy strength was on par with or better than those components that were manufactured using traditional manufacturing processes such as forging and extrusion. It is noteworthy that the ShAPE samples were made with considerably lower number of process steps.

36 MATERIALS SCIENCE↗

Scalable deep learning for watershed model calibration

Watershed models such as the Soil and Water Assessment Tool (SWAT) consist of high-dimensional physical and empirical parameters. These parameters often need to be estimated/calibrated through inverse modeling to produce reliable predictions on hydrological fluxes and states. Existing parameter estimation methods can be time consuming, inefficient, and computationally expensive for high-dimensional problems. In this paper, we present an accurate and robust method to calibrate the SWAT model (i.e., 20 parameters) using scalable deep learning (DL). We developed inverse models based on convolutional neural networks (CNN) to assimilate observed streamflow data and estimate the SWAT model parameters. Scalable hyperparameter tuning is performed using high-performance computing resources to identify the top 50 optimal neural network architectures. We used ensemble SWAT simulations to train, validate, and test the CNN models. We estimated the parameters of the SWAT model using observed streamflow data and assessed the impact of measurement errors on SWAT model calibration. We tested and validated the proposed scalable DL methodology on the American River Watershed, located in the Pacific Northwest-based Yakima River basin. Our results show that the CNN-based calibration is better than two popular parameter estimation methods (i.e., the generalized likelihood uncertainty estimation [GLUE] and the dynamically dimensioned search [DDS], which is a global optimization algorithm). For the set of parameters that are sensitive to the observations, our proposed method yields narrower ranges than the GLUE method but broader ranges than values produced using the DDS method within the sampling range even under high relative observational errors. The SWAT model calibration performance using the CNNs, GLUE, and DDS methods are compared using R 2 and a set of efficiency metrics, including Nash-Sutcliffe, logarithmic Nash-Sutcliffe, Kling-Gupta, modified Kling-Gupta, and non-parametric Kling-Gupta scores, computed on the observed and simulated watershed responses. The best CNN-based calibrated set has scores of 0.71, 0.75, 0.85, 0.85, 0.86, and 0.91. The best DDS-based calibrated set has scores of 0.62, 0.69, 0.8, 0.77, 0.79, and 0.82. The best GLUE-based calibrated set has scores of 0.56, 0.58, 0.71, 0.7, 0.71, and 0.8. The scores above show that the CNN-based calibration leads to more accurate low and high streamflow predictions than the GLUE and DDS sets. Our research demonstrates that the proposed method has high potential to improve our current practice in calibrating large-scale integrated hydrologic models.

54 ENVIRONMENTAL SCIENCES↗

Using Direct Numerical Simulation of Pore-Level Events to Improve Pore-Network Models for Prediction of Residual Trapping of CO2

Direct numerical simulation and pore-network modeling are common approaches to study the physics of two-phase flow through natural rocks. For assessment of the long-term performance of geological sequestration of CO 2 , it is important to model the full drainage-imbibition cycle to provide an accurate estimate of the trapped CO 2 . While direct numerical simulation using pore geometry from micro-CT rock images accurately models two-phase flow physics, it is computationally prohibitive for large rock volumes. On the other hand, pore-network modeling on networks extracted from micro-CT rock images is computationally efficient but utilizes simplified physics in idealized geometric pore elements. This study uses the lattice-Boltzmann method for direct numerical simulation of CO 2 -brine flow in idealized pore elements to develop a new set of pore-level flow models for the pore-body filling and snap-off events in pore-network modeling of imbibition. Lattice-Boltzmann simulations are conducted on typical idealized pore-network configurations, and the interface evolution and local capillary pressure are evaluated to develop modified equations of local threshold capillary pressure of pore elements as a function of shape factor and other geometrical parameters. The modified equations are then incorporated into a quasi-static pore-network flow solver. The modified model is applied on extracted pore-network of sandstone samples, and saturation of residual trapped CO 2 is computed for a drainage-imbibition cycle. The modified model yields different statistics of pore-level events compared with the original model; in particular, the occurrence of snap-off in pore-throats is reduced resulting in a more frontal displacement pattern along the main injection direction. Compared to the original model, the modified model is in closer agreement with the residual trapped CO 2 obtained from core flow experiments and direct numerical simulation.

Kohanpur, Amir H.↗

Back Trajectories during SAIL

This data set contains air mass back trajectories generated during the Surface Atmosphere Integrated Field Laboratory (SAIL) campaign using the Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT, Stein et al. 2015, Rolph et al. 2017). For each hour from January 1, 2021 to June 30, 2023, a 96-hour back trajectory was initiated at the SAIL sampling site, starting 100 m above ground level. The calculations used the GDAS meteorological data set with model vertical velocity. The files within this dataset show location and meteorological parameters of the back trajectory analysis. There are 21 columns within each data file. The first 2 columns pertain to back trajectory parameters: the trajectory number (which is always 1) and the meteorological grid number. The next 9 columns pertain to spatiotemporal information: the year, month, day, hour (0-23), and minute of the trajectory, the previous time of the trajectory (in negative hours), and the airmass location (latitude, longitude, and height in m AGL). The rest of the parameters within each file pertains to surface and meteorological characteristics of the airmass: pressure (in hPa), potential temperature (in K), temperature (in K), precipitation rate (in mm/hr), mixing depth (in m), relative humidity (in % relative to liquid water), specific humidity (in g/kg), H2O mixing ratio (in g/kg), underlying terrain height (in meters above sea level), and incoming radiation (in W/m^2).

back trajectory↗

Optimizing Grain Boundary Structures with LAMMPS Using Evolutionary Algorithms

Grain boundary structure optimization is an important part of materials modeling. Current methods for grain boundary structure optimization involve inefficient, time-consuming processes that do not fully explore the interface parameter space. Evolutionary algorithms have recently been demonstrated to be effective at determining both stable and metastable grain boundary interface structures. In this work, we demonstrate the use of GBOpt, a grain boundary structure optimization software designed to use the Large-scale Atomic/Molecular Massively Parallel Simulation (LAMMPS) software to efficiently determine grain boundary structures. We demonstrate that a only a few manipulations, namely atom insertion, atom removal, and relative grain displacement, are sufficient to explore much of the grain boundary structure parameter space. The efficacy of this approach is demonstrated on an FCC Ni system, and a BCC Fe system. The computational cost is compared against the gamma-surface sampling approach to demonstrate performance improvement.

Evolutionary algorithms↗

ZTF-SEDm Type Ia supernova sample for Twins Embedding spectrophotometric standardization

Aims. This paper has two aims: the first aim is to build a large homogeneous spectrophotometric sample of Type Ia supernovae (SNe Ia) from the second data release of the Zwicky Transient Facility (ZTF DR2). We used the spectrum sample from the low-resolution ( R ∼ 100) SEDmachine (SEDm) Integral Field Spectrograph (IFS) that gathered 3069 spectra. This is one of the largest samples of such collections that can attempt to reproduce the Twins Embedding (TE) spectrophotometric standardization method. This is our second objective. The method was developed based on high-quality spectra from 200 SNe Ia of the Nearby Supernova factory (SNfactory) and led to an exceptionally low value of 0.073 mag for the intrinsic scatter. Methods. As the SEDm is not designed as a spectrophotometric instrument, we first improved the flux-calibration accuracy of the SN Ia spectrum sample using the ZTF photometric data, which were calibrated at the percent level. We corrected the spectra for second-order polynomials, fitted by comparing the synthetic photometry in the ZTF g , r , i filters with the light-curve (LC) data. We then applied the three steps of the TE parameterization to a subset of 783 ZTF SN spectra near maximum light while comparing results from SNfactory and ZTF. We finally analyzed the standardization methods based on the TE parameters. Results. The precision of the phase-correction model, which is the first step of the TE, is estimated at 0.01 mag in g band based on ZTF data. Despite the challenge posed by the spectrum-extraction pipeline associated with the SEDm (flux calibration, leftover host signal, low signal-to-noise ratio, and low resolution), we applied a first standardization in color based on the second step of the TE, called read between the lines (RBTL), to the ZTF sample. We reached a Hubble residual scatter of 0.153 mag, all in normalized median absolute deviation, which is to be compared to the ∼0.11 mag obtained with the SNfactory data. The SALT color and stretch standardization reaches a scatter of 0.164 mag for the same ZTF SN Ia sample, and its host steps are ∼0.1 mag and zero for RBTL. When considering the scatter due to the redshift error and flux calibration error, we estimated a RBTL scatter of ∼0.129 mag for this ZTF sample as an upper limit because we identified an additional contribution from a systematic error in color. We tested the standardization based on the nonlinear TE parameters, and, as expected from the low spectrum quality, it did not improve the overall dispersion. Conclusions. We release 1897 flux calibrated spectra of 1607 SNe Ia with an estimated photometric accuracy of 0.07 mag. We further demonstrate that some amount of spectrophotometric SN Ia standardization can be done with limited-quality spectra. The RBTL standardization is more efficient than that of SALT with one parameter less, and the resulting host steps are consistent with zero. This makes it less prone to astrophysical bias. For future spectroscopic surveys, targeting the extraction pipeline for a thorough flux calibration and good signal-to-noise ratio would enable us to compute the full TE standardization, which would further reduce the scatter in the distance estimate.

Ganot, C↗

Parsec-scale Properties of Steep- and Flat-spectrum Extragalactic Radio Sources from a VLBA Survey of a Complete North Polar Cap Sample

We observed with the Very Long Baseline Array (VLBA) at 2.3 and 8.6 GHz a complete flux-density-limited sample of 482 radio sources with decl. > +75° brighter than 200 mJy at 1.4 GHz drawn from the NVSS catalog. A total of 34% of the sources show parsec-scale emission above the flux density detection limit of 30 mJy; their accurate positions and parsec-scale structure parameters are determined. Among all the sources detected at least at the shortest VLBA baselines, the majority, or 72%, have a steep single-dish spectrum. The fraction of the sources with a detectable parsec-scale structure is above 95% among the flat-spectrum objects and close to 25% among the steep-spectrum objects. We identified 82 compact steep-spectrum source candidates, which make up 17% of the sample; most of them are reported for the first time. The compactness and the brightness temperature of the sources in our sample show a positive correlation with single-dish and VLBA spectral indices. All the sources with a significant 8 GHz variability were detected by the VLBA snapshot observations, which independently confirmed their compactness. We demonstrated that 54% of the sources detected by the VLBA at 2.3 GHz in our sample have a steep VLBA spectrum. The compact radio emission of these sources is likely dominated by optically thin jets or mini-lobes, not by an opaque jet core. These results show that future VLBI surveys aimed at searching for new sources with parsec-scale structure should include not only flat-spectrum sources but also steep-spectrum ones in order to reach an acceptable level of completeness.

47 OTHER INSTRUMENTATION↗

Analysis of Spectral Lines in Large Databases of Synthetic Spectra for Massive Stars

In this paper, we describe a program that identifies in the optical spectrum the main parameters of a spectral line, namely the initial and final wavelengths, and the line depth. Moreover, using numerical calculations, it identifies and removes adjacent lines. Next, the program calculates the equivalent width and the FWHM. The software was tested in a sample of 300 lines in two databases of synthetic spectra generated by the CMFGEN and PoWR codes, and 300 lines in observed spectra from the IACOB database, showing a Gaussian distribution of relative errors, from which it is inferred that 80% of the measured lines have errors less than 17% and only 5% of the lines have errors greater than 26%. The program was also run on the entire database of 45,000 CMFGEN and 202 POWR synthetic spectra, generating a library of H i, He i, and He ii lines necessary to feed the FITspec code for the derivation of stellar parameters: effective temperature, surface gravity, and luminosity.

47 OTHER INSTRUMENTATION↗

The self-annealing of irradiation induced defects in magnetite Fe 3 O 4 : Revealing reversible irradiation-induced disorder transformation through in situ TEM

This work reports heavy ion-irradiation effects in polycrystalline Fe 3 O 4 . For this matter, Fe/Fe 3 O 4 heterostructures were irradiated in situ in a transmission electron microscopy with 1 MeV Kr ions at 50 K. Evidence of cubic to monoclinic transformation (a.k.a Verwey transition) was recorded in some magnetite grains upon cooling the sample (around 90 K); however, most of the oxide grains retain their cubic spinel structure. Irradiation effects were analyzed in the cubic phase up to a maximum dose of 38 dpa without the sign of amorphization. The extinction of first-order reflections was recorded at doses below 1 dpa, indicating the formation of a new (metastable) phase with half of the lattice parameters compared to the unirradiated Fe 3 O 4 crystal. The formation of the new crystalline phase, which also presents a high resistance to amorphization, is related to the disordering of the cation lattice and the high mobility of the cation interstitials. The metastable phase readily recovers around 225 K during the natural warming of the sample from 50 K to room temperature.

Atom irradiation effects↗

A novel methodology for assessing the hygroscopicity of aerosol filter samples

Abstract. Due to US regulations, concentrations of hygroscopic inorganic sulfate and nitrate have declined in recent years, leading to an increased importance of the hygroscopic nature of organic matter (OM). The hygroscopicity of OM is poorly characterized because only a fraction of the multitude of organic compounds in the atmosphere is readily measured, and there is limited information on their hygroscopic behaviors. Hygroscopicity of aerosol is traditionally measured using a humidified tandem differential mobility analyzer (HTDMA) or electrodynamic balance (EDB). EDB measures water uptake by a single particle. For ambient and chamber studies, HTDMA measurements provide water uptake and particle size information but not chemical composition. To fill this information gap, we developed a novel methodology to assess the water uptake by particles collected on Teflon filters. This method uses the same filter sample for both hygroscopicity measurements and chemical characterization, thereby providing an opportunity to link the measured hygroscopicity with ambient particle composition. To test the method, hygroscopic measurements were conducted in the laboratory for ammonium sulfate, sodium chloride, glucose, and malonic acid, which were collected on 25 mm Teflon filters using an aerosol generator and sampler. Constant-humidity solutions (CHSs), including potassium chloride, barium chloride dihydrate, and potassium sulfate, were employed in a saturated form to maintain the relative humidity (RH) at approximately 84 %, 90 %, and 97 % in small chambers. Our preliminary experiments revealed that, without the pouch, water uptake measurements were not feasible due to rapid water loss during weighing. Additionally, we observed some absorption by the aluminum pouch itself. To account for this, concurrent measurements were conducted for both the loaded and the blank filters at each RH level. Thus, the dry loaded and blank Teflon filters were placed in aluminum pouches with one side open and in RH-controlled chambers for more than 24 h. The wet loaded samples and wet blanks were then weighed using an ultramicrobalance to determine the water uptake by the respective compound and the blank Teflon filter. The net amount of water absorbed by each compound was calculated by subtracting the water uptake of the blank filter from that of the wet loaded filter. Hygroscopic parameters, including the water-to-solute (W / S) ratio, molality, mass fraction solute (mfs), and growth factors (GFs), were calculated from the measurements. The results obtained are consistent with those reported by the Extended Aerosol Inorganics Model (E-AIM) and previous studies utilizing HTDMA and EDB for these compounds, highlighting the accuracy of this new methodology. This new approach enables the hygroscopicity and chemical composition of individual filter samples to be assessed so that in complex mixtures, such as chamber and ambient samples, the total water uptake can be parsed between the inorganic and organic components of the aerosol.

54 ENVIRONMENTAL SCIENCES↗

Effective data sampling strategies and boundary condition constraints of physics-informed neural networks for identifying material properties in solid mechanics

Material identification is critical for understanding the relationship between mechanical properties and the associated mechanical functions. However, material identification is a challenging task, especially when the characteristic of the material is highly nonlinear in nature, as is common in biological tissue. In this work, we identify unknown material properties in continuum solid mechanics via physics-informed neural networks (PINNs). To improve the accuracy and efficiency of PINNs, we develop efficient strategies to nonuniformly sample observational data. We also investigate different approaches to enforce Dirichlet-type boundary conditions (BCs) as soft or hard constraints. Finally, we apply the proposed methods to a diverse set of time-dependent and time-independent solid mechanic examples that span linear elastic and hyperelastic material space. The estimated material parameters achieve relative errors of less than 1%. As such, this work is relevant to diverse applications, including optimizing structural integrity and developing novel materials.

97 MATHEMATICS AND COMPUTING↗

Modeling Dust Production, Growth, and Destruction in Reionization-era Galaxies with the CROC Simulations: Methods and Parameter Exploration

We introduce a model for the explicit evolution of interstellar dust in a cosmological galaxy formation simulation. We post-process a simulation from the Cosmic Reionization on Computers project, integrating an ordinary differential equation for the evolution of the dust-to-gas ratio along pathlines in the simulation sampled with a tracer particle technique. This model incorporates the effects of dust grain production in asymptotic giant branch star winds and supernovae, grain growth due to the accretion of heavy elements from the gas phase of the interstellar medium, and grain destruction due to thermal sputtering in the high-temperature gas of supernova remnants due to thermal sputtering. A main conclusion of our analysis is the importance of a carefully chosen dust destruction model, for which different reasonable parameterizations can predict very different values at the ∼100 pc resolution of the interstellar medium (ISM) in our simulations. We run this dust model on the single most massive galaxy in a 10 h$^{−1}$ comoving megaparsec box, which attains a stellar mass of ∼2 × 10$^{9}$ M $_{⊙}$ by z = 5. We find that the model is capable of reproducing dust masses and dust-sensitive observable quantities broadly consistent with existing data from high-redshift galaxies. The total dust mass in the simulated galaxy is somewhat sensitive to parameter choices for the dust model, especially the timescale for grain growth due to accretion in the ISM. Consequently, observable quantity observations that can constrain galaxy dust masses at these epochs are potentially useful for placing constraints on dust physics.

79 ASTRONOMY AND ASTROPHYSICS↗