Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter compensation”

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 109 records · Page 6

Using neural network ensembles to separate ocean biogeochemical and physical drivers of phytoplankton biogeography in Earth system models

Abstract. Earth system models (ESMs) are useful tools for predicting and understanding past and future aspects of the climate system. However, the biological and physical parameters used in ESMs can have wide variations in their estimates. Even small changes in these parameters can yield unexpected results without a clear explanation of how a particular outcome was reached. The standard method for estimating ESM sensitivity is to compare spatiotemporal distributions of variables from different runs of a single ESM. However, a potential pitfall of this method is that ESM output could match observational patterns because of compensating errors. For example, if a model predicts overly weak upwelling and low nutrient concentrations, it might compensate for this by allowing phytoplankton to have a high sensitivity to nutrients. Recently, we demonstrated that neural network ensembles (NNEs) are capable of extracting relationships between predictor and target variables within ocean biogeochemical models. Being able to view the relationships between variables, along with spatiotemporal distributions, allows for a more mechanistically based examination of ESM outputs. Here, we investigated whether we could apply NNEs to help us determine why different ESMs produce different spatiotemporal distributions of phytoplankton biomass. We tested this using three cases. The first and second case used different runs of the same ESM, except that the physical circulations differed between them in the first case, while the biological equations differed between them in the second. Our results indicated that the NNEs were capable of extracting the relationships between variables for different runs of a single ESM, allowing us to distinguish between differences due to changes in circulation (which do not change relationships) from changes in biogeochemical formulation (which do change relationships). In the third case, we applied NNEs to two different ESMs. The results of the third case highlighted the capability of NNEs to contrast the apparent relationships of different ESMs and some of the challenges it presents. Although applied specifically to the ocean components of an ESM, our study demonstrates that Earth system modelers can use NNEs to separate the contributions of different components of ESMs. Specifically, this allows modelers to compare the apparent relationships across different ESMs and observational datasets.

54 ENVIRONMENTAL SCIENCES↗

Surface chemistry on a polarizable surface: Coupling of CO with KTaO 3 (001)

Polarizable materials attract attention in catalysis because they have a free parameter for tuning chemical reactivity. Their surfaces entangle the dielectric polarization with surface polarity, excess charge, and orbital hybridization. How this affects individual adsorbed molecules is shown for the incipient ferroelectric perovskite KTaO 3 . This intrinsically polar material cleaves along (001) into KO- and TaO 2 -terminated surface domains. At TaO 2 terraces, the polarity-compensating excess electrons form a two-dimensional electron gas and can also localize by coupling to ferroelectric distortions. TaO 2 terraces host two distinct types of CO molecules, adsorbed at equivalent lattice sites but charged differently as seen in atomic force microscopy/scanning tunneling microscopy. Temperature-programmed desorption shows substantially stronger binding of the charged CO; in density functional theory calculations, the excess charge favors a bipolaronic configuration coupled to the CO. These results pinpoint how adsorption states couple to ferroelectric polarization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Global Atmosphere‐aerosol Model ICON‐A‐HAM2.3–Initial Model Evaluation and Effects of Radiation Balance Tuning on Aerosol Optical Thickness

Abstract The Hamburg Aerosol Module version 2.3 (HAM2.3) from the ECHAM6.3‐HAM2.3 global atmosphere‐aerosol model is coupled to the recently developed icosahedral nonhydrostatic ICON‐A (icon‐aes‐1.3.00) global atmosphere model to yield the new ICON‐A‐HAM2.3 atmosphere‐aerosol model. The ICON‐A and ECHAM6.3 host models use different dynamical cores, parameterizations of vertical mixing due to sub‐grid scale turbulence, and parameter settings for radiation balance tuning. Here, we study the role of the different host models for simulated aerosol optical thickness (AOT) and evaluate impacts of using HAM2.3 and the ECHAM6‐HAM2.3 two‐moment cloud microphysics scheme on several meteorological variables. Sensitivity runs show that a positive AOT bias over the subtropical oceans is remedied in ICON‐A‐HAM2.3 because of a different default setting of a parameter in the moist convection parameterization of the host models. The global mean AOT is biased low compared to MODIS satellite instrument retrievals in ICON‐A‐HAM2.3 and ECHAM6.3‐HAM2.3, but the bias is larger in ICON‐A‐HAM2.3 because negative AOT biases over the Amazon, the African rain forest, and the northern Indian Ocean are no longer compensated by high biases over the sub‐tropical oceans. ICON‐A‐HAM2.3 shows a moderate improvement with respect to AOT observations at AERONET sites. A multivariable bias score combining biases of several meteorological variables into a single number is larger in ICON‐A‐HAM2.3 compared to standard ICON‐A and standard ECHAM6.3. In the tropics, this multivariable bias is of similar magnitude in ICON‐A‐HAM2.3 and in ECHAM6.3‐HAM2.3. In the extra‐tropics, a smaller multivariable bias is found for ICON‐A‐HAM2.3 than for ECHAM6.3‐HAM2.3.

54 ENVIRONMENTAL SCIENCES↗

Two-stage nonlinear compression of high-power femtosecond laser pulses

Two-stage compression of laser pulses with a power of 250 TW is experimentally realised by broadening their spectrum during self-phase modulation in fused silica and subsequent dispersion compensation upon reflection from chirping mirrors. A five-fold decrease in the duration is demonstrated, from 75 to 15 fs, with a B-integral value of about 5 at each stage. It is possible to avoid small-scale self-focusing due to self-filtering of the laser beam during free propagation in vacuum. With optimal parameters of the dispersive mirror, the pulse can be compressed to a duration of less than 5 fs. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Detecting thermodynamic phase transition via explainable machine learning of photoemission spectroscopy

Identifying thermodynamic signatures of electronic phases, such as superconductivity, is challenging in low-dimensional materials due to strong fluctuations and low probing volume. Spectroscopic methods are often used to identify new bulk phases, but their main measurable quantity—electronic energy gaps—is no longer an effective order parameter in low-dimensional and fluctuating systems. Combining angle-resolved photoemission with a domain-adversarial neural network, we report a data-driven method to identify thermodynamic phase transitions solely based on single-particle spectra. We demonstrate 97.6% accuracy in cuprate superconductor Bi 2 Sr 2 CaCu 2 O 8+δ with strong superconducting fluctuations. This model notably compensates for the scarcity of experimental data by leveraging virtually inexhaustible simulated data. Further, its explainability reveals the crucial role of in-gap spectral weight in detecting phase fluctuations and thermodynamic transitions. Our work pinpoints the spectroscopic signatures of fluctuating orders and enables using spectroscopy for machine-learning-assisted material discovery for low-dimensional and strong coupling systems.

2D materials↗

Shear reflectivity compensation in full-waveform inversion using least-squares reverse-time migration

SUMMARY The computational cost of elastic-waveform inversion is too high for inverting PP reflections, while using acoustic full-waveform inversion (FWI) is inaccurate because it does not depend on the shear modulus/velocity/impedance that affects elastic PP wavefield amplitudes. To solve this problem, we develop a waveform inversion method that uses acoustic least-squares reverse-time migration (LSRTM) to compensate the shear reflectivity for acoustic FWI. Our method is based on the quasi-elastic-wave equation developed by Chapman et al. (2014). The quasi-elastic-wave equation uses a linearized acoustic-wave equation with shear modulus μ as a virtual source to correct the acoustic PP wavefield amplitudes toward elastic ones. Our waveform inversion method inverts for elastic parameters by minimizing the L2 norm of the difference between recorded and predicted PP reflections modelled using the quasi-elastic-wave equation. Numerical tests on synthetic and field data show that our method can properly handle the amplitudes of elastic PP reflections and provides an accurate estimate of the P- and S-wave velocities/impedances and, in some cases, the density. The method does not need the computationally expensive numerical solution to the elastic-wave equation. It also gives a better estimate of elastic parameters than a pure LSRTM method for elastic PP reflections.

Feng, Zongcai↗

Understanding water and energy fluxes in the Amazonia: Lessons from an observation‐model intercomparison

Abstract Tropical forests are an important part of global water and energy cycles, but the mechanisms that drive seasonality of their land‐atmosphere exchanges have proven challenging to capture in models. Here, we (1) report the seasonality of fluxes of latent heat (LE), sensible heat ( H ), and outgoing short and longwave radiation at four diverse tropical forest sites across Amazonia—along the equator from the Caxiuanã and Tapajós National Forests in the eastern Amazon to a forest near Manaus, and from the equatorial zone to the southern forest in Reserva Jaru; (2) investigate how vegetation and climate influence these fluxes; and (3) evaluate land surface model performance by comparing simulations to observations. We found that previously identified failure of models to capture observed dry‐season increases in evapotranspiration (ET) was associated with model overestimations of (1) magnitude and seasonality of Bowen ratios (relative to aseasonal observations in which sensible was only 20%–30% of the latent heat flux) indicating model exaggerated water limitation, (2) canopy emissivity and reflectance (albedo was only 10%–15% of incoming solar radiation, compared to 0.15%–0.22% simulated), and (3) vegetation temperatures (due to underestimation of dry‐season ET and associated cooling). These partially compensating model‐observation discrepancies (e.g., higher temperatures expected from excess Bowen ratios were partially ameliorated by brighter leaves and more interception/evaporation) significantly biased seasonal model estimates of net radiation ( R n ), the key driver of water and energy fluxes (LE ~ 0.6 R n and H ~ 0.15 R n ), though these biases varied among sites and models. A better representation of energy‐related parameters associated with dynamic phenology (e.g., leaf optical properties, canopy interception, and skin temperature) could improve simulations and benchmarking of current vegetation–atmosphere exchange and reduce uncertainty of regional and global biogeochemical models.

Restrepo‐Coupe, Natalia↗

Achieving 100 GW idler pulses from an existing petawatt optical parametric chirped pulse amplifier

Optical parametric chirped-pulse-amplification produces two broadband pulses, a signal and an idler, that can both provide peak powers >100 GW. In most cases the signal is used, but compressing the longer-wavelength idler opens up opportunities for experiments where the driving laser wavelength is a key parameter. This paper will describe several subsystems that were added to a petawatt class, Multi-Terawatt optical parametric amplifier line (MTW-OPAL) at the Laboratory for Laser Energetics to address two long-standing issues introduced by the use of the idler, angular dispersion, and spectral phase reversal. To the best of our knowledge, this is the first time that compensation of angular dispersion and phase reversal has been achieved in a single system and results in a 100 GW, 120-fs duration, pulse at 1170 nm.

47 OTHER INSTRUMENTATION↗

Boosts in leaf-level photosynthetic capacity aid Pinus ponderosa recovery from wildfire

Forests mitigate climate change by sequestering massive amounts of carbon, but recent increases in wildfire activity are threatening carbon storage. Currently, our understanding of wildfire impacts on forest resilience and the mechanisms controlling post-fire recovery remains unresolved due to a lack of empirical data on mature trees in natural settings. Here, we quantify the physiological mechanisms controlling carbon uptake immediately following wildfire in mature individuals of ponderosa pine (Pinus ponderosa), a wide-spread and canopy-dominant tree species in fire-prone forests. While photosynthetic capacity was lower in burned than unburned trees due to an overall depletion of resources, we show that within the burned trees, photosynthetic capacity increases with the severity of damage. Our data reveal that boosts in the efficiency of carbon uptake at the leaf-level may compensate for whole-tree damage, including the loss of leaf area and roots. We further show that heightened photosynthetic capacity in remaining needles on burned trees may be linked with reduced water stress and leaf nitrogen content, providing pivotal information about post-fire physiological processes. Our results have implications for Earth system modeling efforts because measurements of species-level physiological parameters are used in models to predict ecosystem and landscape-level carbon trajectories. Finally, current land management practices do not account for physiological resilience and recovery of severely burned trees. Our results suggest premature harvest may remove individuals that may otherwise survive, irrevocably altering forest carbon balance.

54 ENVIRONMENTAL SCIENCES↗

Data Assimilation using Non-invasive Monte Carlo Sensitivity Analysis of Reactor Kinetics Parameters [Slides]

Sensitivity coefficients of prompt neutron decay constant (α) and effective delayed neutron fraction (β eff ) to Pu-239 nuclear data were calculated using a newly available tool (ACEtk). Uncertainty calculations showed the following trends: (1) The prompt neutron decay constant can be used to reduce nuclear data-induced uncertainty in Pu-239(n,f) and (2) the effective delayed neutron fraction can be used to effectively reduce nuclear data-induced uncertainty in Pu-239 total fission $\tilde{v}$. Sensitivity/uncertainty analysis can be used to determine optimal experiment/detector setup to identify compensating errors in nuclear data. Upcoming work includes investigating other nuclide-reaction pairs and response sensitivities, such as Δ$\rho$, neutron leakage spectra, and reaction rate measurements to constrain nuclear data of interest.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Excitation of long electrostatic solitary waves in ion beam neutralization process

Unusually long electrostatic solitary waves (ESWs) are discovered in a particle-in-cell simulation study of the process of ion beam neutralization by electron emission from a filament. These ESWs are long because the density perturbation responses to the potential wells created by the ESWs are very small. Notably, the density perturbation is small because the trapped (positive) and untrapped (negative) electron density perturbations nearly compensate each other because of a non-Maxwellian electron velocity distribution in the beam.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

PIP-II Linac Cryogenic Distribution System Design Challenges

The PIP-II linac Cryogenic Distribution System (CDS) is characterized by extremely small heat inflows and robust mechanical design. It consists of a Distribution Valve Box (DVB), Intermediate Transfer Line, Tunnel Transfer Line, comprising 25 Bayonet Cans, and ends with a Turnaround Can. Multiple helium streams, each characterized by distinct helium parameters, flow through each of these elements. The CDS geometry allows maintaining an acceptable pressure drop for each helium stream, considering the planned flows and helium parameters in different operation modes. This is particularly crucial for the return line of helium vapors, which return from cryomodules to the cold compressors and thus have very restrictive pressure drop requirements. On both sides of the DVB there are fixed supports for process pipes. One of the DVB design challenges was to route the process pipes in such a way that their shape provided sufficient compensation for thermal shrinkage. This ensures th at the forces resulting from thermal shrinkage acting on the cryogenic valves remain at a level acceptable to the manufacturer. The required thermal budget of the CDS was achieved by thermo-mechanical optimization of its components, like process pipes fixed supports in Bayonet Cans.

43 PARTICLE ACCELERATORS↗

Unveiling the influence of selective-area-regrowth interfaces on local electronic properties of GaN p-n junctions for efficient power devices

Here, we report correlated nanoscale mapping of the structure, composition, and properties of regrown GaN p-n junctions to identify how etching and non-planar regrowth processes limit diode performance via the introduction of unintentional dopants and defect states. p-GaN was selectively regrown in n-GaN trenches with SiO 2 masks of variable mask-to-trench-width ratio. Dilute Al layers were periodically introduced during regrowth as markers of the growth interface. Correlated nanoscale mapping of doping, conductivity, and dopant complexes was achieved with atom probe tomography (APT), scanning spreading resistance microscopy (SSRM), and cathodoluminescence (CL) spectroscopy, respectively. The Al marker layers, detected by APT, enabled reconstruction of the faceted growth interface and correlation of the dopant concentration with position and time. The p-GaN growth rate is proportional to the mask-to-trench width ratio while the dopant incorporation rate is invariant. At trench edges, magnesium incorporation is suppressed, and oxygen incorporation enhanced, due to preferential incorporation on the semi-polar growth surface, leading to compensation and less abrupt p-n junctions; the SiO 2 mask is a source of oxygen. Residual etch damage below the regrowth interface induces n-type and p-type conductivity, creating leakage pathways. The non-uniform Mg incorporation is driven by crystal anisotropy and is thus inherent to non-planar regrowth, but can be mitigated by engineering the regrowth interface and process parameters. The unprecedented integration of spatially resolved mapping of dopants, impurities, conductivity, and carrier type is a powerful approach to discriminating distinct factors that limit the performance of regrown diodes, enabling the rational optimization of process and device design.

36 MATERIALS SCIENCE↗

Inertial Confinement Fusion Design Search Using Bayesian Optimization

Inertial confinement fusion (ICF) experiments rely on complex multi-physics simulation codes such as the Lawrence Livermore National Laboratory-developed HYDRA to guide design work. However, these simulations have several dozen tunable parameters and can be computationally expensive. This makes searching the parameter space challenging and time-consuming. Recently developed automated tools utilize Bayesian optimization to search these high-dimensional parameter spaces for optimal designs. The optimization tools run 2D integrated simulations in HYDRA to converge on a design that produces specified scalar or vector outputs. In this paper, we apply the Bayesian optimization tools to two common tuning scenarios. First, we tune simulation inputs to match measurements of a well-characterized experiment at the National Ignition Facility. This type of tuning is commonly performed to compensate for the use of simplified simulation settings (e.g. reduced resolution) or to account for missing physics in the simulations. Second, we search for an ICF simulation design that has a particular radiation drive profile. These optimizations replicate the kinds of tuning researchers routinely perform, but do so with significantly reduced manual effort. This approach demonstrates a powerful and efficient pathway toward autonomous, high-fidelity design optimization for future ICF experiments.

Bayesian optimization↗

Estimation of the Required Dipole Corrector Magnetic Field for PIP-II Injection Based on Beam Studies

The Fermilab Booster will accept a 600 µs beam pulse from the new superconducting Linac for PIP-II operations. The Booster is a rapid-cycling synchrotron that uses a resonant magnet circuit ramping at 15 Hz. For PIP-II, the cycle rate will increase to 20 Hz, and the injection pulse length will expand from 40 µs to 600 µs due to the lower output current from the new Linac. Because the Booste main bending field follows a sinusoidal waveform, the magnetic field is not constant during the extended injection window. The longer pulse length and higher repetition rate modify the beam orbit and can lead to increased beam losses. The Booster contains 48 dipole-corrector packages distributed across its 24 periods. Each package includes horizontal and vertical dipoles, quadrupole, sextupole, skew-quadrupole, and skew-sextupole elements. By driving the dipole correctors with an appropriately shaped sinusoidal waveform during injection, we can compensate the changing main field and create an effectively flat bending field—referred to as flat injection. Over the past several years, machine studies have been performed to characterize the required correction fields and to determine the corresponding power-supply specifications needed for PIP-II operation. In this presentation, we will summarize the study results and discuss the estimated magnetic field requirements and power supply parameters for achieving flat injection in the Booster.

Seiya, K. [Fermilab] (ORCID:0000000250576943)↗

Designing Particle Morphologies for Materials with Solid Transport Limitations: A Case Study of Lithium and Manganese Rich Cathode Oxides

A lithium and manganese rich nickel-manganese-cobalt oxide (LMR-NMC) cathode is a promising candidate for next-generation batteries due to its high specific capacity, low cost, and low cobalt content. However, the material suffers from poor rate capability due to the diffusion limitations of lithium in the cathode particles. Understanding the material performance requires careful control of the morphology of the cathode particles, taking into account the primary and agglomerated diffusion pathways and the presence of pores, some of which could be closed from electrolyte infiltration. Here, in this study, we use a microstructure-based mathematical model combined with experimental data to understand the role of the complex cathode particle morphology in the rate performance of the material. Scanning electron microscopy images of cathodes made under different synthesis conditions, which results in different agglomerate morphologies, serve as the input into the mathematical model. The model is then compared to rate data to understand the controlling parameters. The presence of intra-agglomerate closed pores results in a large agglomerate diffusion length in comparison to the ideal condition, where the primary particles are agglomerated in an open and dispersed manner such that the entire interfacial area is available for electrochemical reaction. Smaller primary and agglomerate diffusion lengths result in better electrochemical performance. This points us toward designing the morphology of the cathode particles to compensate for the diffusion limitation of LMR-NMC while maximizing the density.

Tewari, Deepti↗

Disentangling Rotational Dynamics and Ordering Transitions in a System of Self-Organizing Protein Nanorods via Rotationally Invariant Latent Representations

The dynamics of complex ordering systems with active rotational degrees of freedom exemplified by protein self-assembly is explored using a machine learning workflow that combines deep learning-based semantic segmentation and rotationally invariant variational autoencoder-based analysis of orientation and shape evolution. The latter allows for disentanglement of the particle orientation from other degrees of freedom and compensates for lateral shifts. The disentangled representations in the latent space encode the rich spectrum of local transitions that can now be visualized and explored via continuous variables. The time dependence of ensemble averages allows insight into the time dynamics of the system and, in particular, illustrates the presence of the potential ordering transition. Finally, analysis of the latent variables along the single-particle trajectory allows tracing these parameters on a single-particle level. The proposed approach is expected to be universally applicable for the description of the imaging data in optical, scanning probe, and electron microscopy seeking to understand the dynamics of complex systems where rotations are a significant part of the process.

representation learning↗

Estimating List-Mode Data Sensitivities to Nuclear Data with MCNP6

Nuclear data are a vital component of predictive simulations used in applications like experiment design, stockpile stewardship, nuclear nonproliferation/safeguards, health physics, and criticality safety. A singular simulation requires the coalescence of different areas of nuclear data such as cross sections, angular distributions, and energy distributions of emitted neutrons for different materials and energy ranges. Improving nuclear data and thus reducing the uncertainty in simulated parameters could enable smaller, better-informed safety factors and ultimately reduce operational and procedural costs. There is a constant effort to garner a better understanding of the physical quantities represented by nuclear data through experiments. Integral experiment benchmarks use simulated and measured results to validate current nuclear data values. In the past, benchmarks primarily focused on the effective multiplication factor (k eff ); however, this limited scope has caused compensating errors and areas of nuclear data that lack validation. Compensating errors are inaccuracies in nuclear data that are obfuscated by cancellation when observing integrated values such as k eff . Diverse integral benchmark experiments that look for quantities of interest other than k eff and include multiple responses minimize the possibility of compensating errors and provides validation to areas of nuclear data previously lacking experimental validation. Benchmark experiments can be optimized during the design process to be highly dependent on specific areas of nuclear data. The dependence of a response in an experiment to a specific area/type of nuclear data is defined as sensitivity. A larger sensitivity means that nuclear data uncertainties will play a larger role in the response(s) resulting in larger bias. Currently, the sensitivity capabilities of the Monte Carlo N-Particle (MCNP ®1 ) transport code are limited to responses of k eff and tallied values (e.g., flux, surface current). As a part of the EUCLID project, this work explores estimating list-mode nuclear data sensitivities that can be used to design experiments aimed to constrain and reduce compensating errors in nuclear data by focusing on responses other than k eff . Tallied values are ideal quantities that are estimated with detectors during experiments. List-mode data (a list of neutron collection times) are the direct output of detector systems in subcritical neutron noise experiments. Expanding MCNP sensitivity capabilities to include the sensitivity of responses estimated from list-mode data, such as the prompt neutron decay constant (α) and multiplicity estimates (S and D), enables more direct comparison of simulated and measured experimental quantities. Additionally, deterministic tools such as SENSMG are capable of obtaining sensitivities to a wide variety of responses; however, these tools cannot handle complex geometries due to the assumptions made in discretizing the phase-space variables of the Boltzman transport equation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗