Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear analysis”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Exploring the Effects of Population and Employment Characteristics on Truck Flows: An Analysis of NextGen NHTS Origin-Destination Data

Truck transportation remains the dominant mode of US freight transportation because of its advantages, such as the flexibility of accessing pickup and drop-off points and faster delivery. Because of the massive freight volume transported by trucks, understanding the effects of population and employment characteristics on truck flows is critical for better transportation planning and investment decisions. The US Federal Highway Administration published a truck travel origin-destination data set as part of the Next Generation National Household Travel Survey program. This data set contains the total number of truck trips in 2020 within and between 583 predefined zones encompassing metropolitan and nonmetropolitan statistical areas within each state and Washington, DC. In this study, origin-destination-level truck trip flow data was augmented to include zone-level population and employment characteristics from the US Census Bureau. Census population and County Business Patterns data were included. The final data set was used to train a machine learning algorithm-based model, Extreme Gradient Boosting (XGBoost), where the target variable is the number of total truck trips. Shapley Additive ExPlanation (SHAP) was adopted to explain the model results. Results showed that the distance between the zones was the most important variable and had a nonlinear relationship with truck flows.

Uddin, Majbah↗
DOE OSTI2020

Uncertainty amplification due to density/refractive index gradients in background-oriented schlieren experiments

Here, we theoretically analyze the effect of density/refractive index gradients on the measurement precision of background-oriented schlieren (BOS) experiments by deriving the Cramer–Rao lower bound (CRLB) for the 2D centroid estimation process. A model is derived for the diffraction limited image of a dot viewed through a medium containing density gradients that includes the effect of the experimental parameters such as the magnification and f-number. It is shown using the model that nonlinearities in the density gradient field lead to blurring of the dot image. This blurring amplifies the effect of image noise on the centroid estimation process, leading to an increase in the CRLB and a decrease in the measurement precision. The ratio of position uncertainties of a dot in the reference and gradient images is shown to be a function of the ratio of the dot diameters and dot intensities. We termed this parameter the amplification ratio (A F ), and a methodology for reporting position uncertainties in tracking-based BOS measurements is proposed. The theoretical predictions of the dot position estimation variance from the CRLB are compared to ray tracing simulations, and agreement is obtained. The uncertainty amplification is also demonstrated on experimental BOS images of flow induced by a spark discharge, where it is seen that regions of high amplification ratio correspond to regions of density gradients. This analysis elucidates the dependence of the position uncertainty on density and refractive index gradient-induced distortion parameters, provides a methodology for accounting its effect on uncertainty quantification and provides a framework for optimizing experiment design.

42 ENGINEERING↗
DOE OSTI2023

Microscopic Origins of the Nonlinear Behavior of Particle-Filled Rubber Probed with Dynamic Strain XPCS

The underlying microscopic response of filler networks in reinforced rubber to dynamic strain is not well understood due to the experimental difficulty of directly measuring filler network behavior in samples undergoing dynamic strain. This difficulty can be overcome with in situ X-ray photon correlation spectroscopy (XPCS) measurements. The contrast between the silica filler and the rubber matrix for X-ray scattering allows us to isolate the filler network behavior from the overall response of the rubber. This in situ XPCS technique probes the microscopic breakdown and reforming of the filler network structure, which are responsible for the nonlinear dependence of modulus on strain, known in the rubber science community as the Payne effect. These microscopic changes in the filler network structure have consequences for the macroscopic material performance, especially for the fuel efficiency of tire tread compounds. In this work, we elucidate the behavior with in situ dynamic strain XPCS experiments on industrially relevant, vulcanized rubbers filled (13 vol %) with novel air-milled silica of ultrahigh-surface area (UHSA) (250 m 2 /g). The addition of a silane coupling agent to rubber containing this silica causes an unexpected and counterintuitive increase in the Payne effect and decrease in energy dissipation. For this rubber, we observe a nearly two-fold enhancement of the storage modulus and virtually equivalent loss tangent compared to a rubber containing a coupling agent and conventional silica. Interpretation of our in situ XPCS results simultaneously with interpretation of traditional dynamic mechanical analysis (DMA) strain sweep experiments reveals that the debonding or yielding of bridged bound rubber layers is key to understanding the behavior of rubber formulations containing the silane coupling agent and high-surface area silica. These results demonstrate that the combination of XPCS and DMA is a powerful method for unraveling the microscale filler response to strain which dictates the dynamic mechanical properties of reinforced soft matter composites. With this combination of techniques, we have elucidated the great promise of UHSA silica when used in concert with a silane coupling agent in filled rubber. Such composites simultaneously exhibit large moduli and low hysteresis under dynamic strain.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Thermal Reservoir Networks for Modularly Expandable Thermal Microgrids

The Department of Defense (DoD) faces the substantial challenge of cost-effectively retrofitting one to two installations per month, each comprising approximately 1,000 buildings, to improve resilience, reduce energy consumption, and enhance energy supply security. Achieving these objectives requires optimal system selection and effective risk mitigation during system integration. To address this need, we introduce Platform-Based Design (PBD), a structured, hierarchical methodology adapted from other industrial sectors to the domain of energy system retrofits. We demonstrate the effectiveness of PBD through a techno-economic feasibility study comparing geothermal-coupled thermal energy networks (TENs) with conventional energy systems for heating, cooling, and powering 17 buildings at Joint Base Andrews (JBA) in Maryland. Our analysis illustrates that the PBD approach enables rigorous, data-driven, sequential decision making, resulting in a family of Pareto-optimal systems, among which the TEN emerged as the most promising solution. The selected TEN design integrates geothermal borefields, heat recovery heat pumps, photovoltaic (PV) arrays, and battery storage. Compared to the baseline system – gas heating combined with air-source chillers – the proposed TEN reduces annual imported energy by 74% and peak electricity demand by 45%, achieves a levelized cost of energy of $\$0.210$/kWh, and substantially enhances resilience. Life-cycle costs increase by approximately 6%, and initial investment costs are about 2.5 times higher than the baseline. However, if central plant infrastructure, district loops, and utility-scale PV and battery systems are privately funded and operated, the initial investment would fall below the baseline system cost. Critical to achieving these significant performance improvements were detailed nonlinear dynamic simulations coupling geothermal heat transfer, energy system operation, and realistic feedback control logic. These simulations identified essential design modifications and control strategy refinements that substantially reduced energy use, peak demand, and compressor shortcycling, thereby improving durability and reliability—issues that would have been significantly more expensive to resolve during operation. Additionally, the verification step highlighted sensitivities to key design parameters that could reduce initial investment by approximately $\$2$ million and reduce annual life-cycle costs more than $\$300,000$. We recommend adopting the PBD methodology for future feasibility studies and TEN pilot projects to gain valuable operational experience. Furthermore, we recommend that DoD invest in transferring and scaling the PBD methodology to other installations. This entails developing standardized computational frameworks and component libraries as well as training industry in conducting PBD. Such investments would enable rapid, robust, reliable, and cost-effective retrofits, supporting DoD’s ambitious energy system modernization goals.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Nonlinear photoconductivities and quantum geometry of chiral multifold fermions

Here, chiral multifold fermions are quasiparticles that appear only in chiral crystals such as transition metal silicides in the cubic B20 structure (i.e., the CoSi family), and they may show exotic physical properties. Here we study the injection and shift photoconductivities and also the related geometrical quantities for several types of chiral multifold fermions, including spin-1/2 as well as pseudospin-1 and -3/2 fermions, dubbed as Kramers Weyl, triple-point, and Rarita-Schwinger-Weyl (RSW) fermions, respectively. We utilize the minimal symmorphic model to describe the triple-point fermions (TPF). We also consider the more realistic model Hamiltonian for the CoSi family including both linear and quadratic terms. We find that injection currents due to circularly polarized light are quantized as a result of the Chern numbers carried by the multifold fermions within the linear models. Surprisingly, we discover that in the TPF model, the linear shift conductivities, responsible for the shift current generation by linearly polarized light, are proportional to the pseudo spin-orbit coupling and independent of photon frequency. In contrast, for the RSW and Kramer Weyl fermions, the linear shift conductivity is linearly proportional to photon frequency. The numerical results agree with the power-counting analysis for quadratic Hamiltonians. The frequency independence of the linear shift conductivity could be attributed to the strong resonant symplectic Christoffel symbols of the flat bands. Moreover, the calculated symplectic Christoffel symbols show significant peaks at the nodes, suggesting that the shift currents are due to the strong geometrical response near the topological nodes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

The Heavy Rainfall Mechanism Revealed by a Terrain-Resolving 4DVar Data Assimilation System—A Case Study

In this study, a newly developed terrain-resolving four-dimensional variational (4DVar)-based data assimilation system, Immersed Boundary Method_Variational Doppler Radar Analysis System (IBM_VDRAS), is applied to investigate the mechanisms leading to a heavy precipitation event that occurred in Taiwan during the Southwesterly Monsoon Experiment (SoWMEX) conducted in 2008. The multivariate analyses using IBM_VDRAS and surface observations reveal that the warm and moist southwesterly flow from the ocean decelerates after making landfall, forming a surface convergence zone along the coast, which is further strengthened during the passage of a prefrontal rainband. The flow ascends as it advances inland until reaching the mountains, producing persistent precipitation and the enhancement of evaporative cooling as well as a widespread high pressure zone. A very shallow (<0.4 km) layer of offshore flow can be identified over the southwestern plain, which helps to generate a quasi-stationary convergence zone near the coast. Furthermore, sensitivity studies are carried out to quantify the relative importance of the contributions made by topographic blockage, evaporative cooling, and their nonlinear interaction, to the evolution of this type of convective system. The influence of the topography is identified as the dominant factor in modulating the flow structure of the rainfall system. However, it is the nonlinear interaction between terrain and evaporation that determines the distribution of the temperature and pressure fields.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The saturation mechanism of thermal instability

The literature on thermal instability (TI) reveals that even for a simple homogeneous plasma, the nonlinear outcome ranges from a gentle reconfiguration of the initial state to an explosive one, depending on whether the condensations that form evolve in an isobaric or nonisobaric manner. After summarizing the recent developments on the linear and nonlinear theory of TI, here we derive several general identities from the evolution equation for entropy that reveal the mechanism by which TI saturates; whenever the boundary of the instability region (the Balbus contour) is crossed, a dynamical change is triggered that causes the comoving time derivative of the pressure to change the sign. This event implies that the gas pressure force reverses direction, slowing the continued growth of condensation. For isobaric evolution, this “pressure reversal” occurs nearly simultaneously for every fluid element in condensation and a steady state is quickly reached. For nonisobaric evolution, the condensation is no longer in mechanical equilibrium and the contracting gas rebounds with greater force during the expansion phase that accompanies the gas reaching the equilibrium curve. The cloud then pulsates because the return to mechanical equilibrium becomes wave mediated. We show that both the contraction rebound event and subsequent pulsation behavior follow analytically from an analysis of the new identities. Our analysis also leads to the identification of an isochoric TI zone and makes it clear that unless this zone intersects the equilibrium curve, isochoric modes can only become unstable if the plasma is in a state of thermal non-equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

On the attribution of the impacts of extreme weather events to anthropogenic climate change

Investigations into the role of anthropogenic climate change in extreme weather events are now starting to extend into analysis of anthropogenic impacts on non-climate (e.g. socio-economic) systems. However, care needs to be taken when making this extension, because methodological choices regarding extreme weather attribution can become crucial when considering the events’ impacts. The fraction of attributable risk (FAR) method, useful in extreme weather attribution research, has a very specific interpretation concerning a class of events, and there is potential to misinterpret results from weather event analyses as being applicable to specific events and their impact outcomes. Using two case studies of meteorological extremes and their impacts, we argue that FAR is not generally appropriate when estimating the magnitude of the anthropogenic signal behind a specific impact. Attribution assessments on impacts should always be carried out in addition to assessment of the associated meteorological event, since it cannot be assumed that the anthropogenic signal behind the weather is equivalent to the signal behind the impact because of lags and nonlinearities in the processes through which the impact system reacts to weather. Whilst there are situations where employing FAR to understand the climate change signal behind a class of impacts is useful (e.g. ‘system breaking’ events), more useful results will generally be produced if attribution questions on specific impacts are reframed to focus on changes in the impact return value and magnitude across large samples of factual and counterfactual climate model and impact simulations. We advocate for constant interdisciplinary collaboration as essential for effective and robust impact attribution assessments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Qualification and Quantification of Porosity at the Top of the Fuel Pins in Metallic Fuels Using Image Processing

Approximately 130,000 metal fuel pins were irradiated in the Experimental Breeder Reactor II (EBR-II) during its 30 years of operation to develop and characterize existing and prospective fuels. For many of the metal fuel irradiation experiments, neutron radiography imaging was performed to characterize fuel behavior, such as fuel axial expansion. While several fuel expansion results obtained from neutron radiography imaging have been published, the analysis of neutron radiography for the purpose of describing statistical properties of porous matter formed on top of the fuel pins, also referred to as fluff in previous publications, is significantly less represented in the literature with just a single paper so far. This study aims to validate and augment results reported in previous publications using automated image processing. The paper describes the statistical properties of the porous matter in terms of nine parameters derived from radiography images and correlates those parameters with such fuel properties as composition, expansion, temperature, and burnup. The reported results are based on 1097 fuel pins of eight different fuel compositions. For three major fuel types, U-10Zr, U-8Pu-10Zr, and U-19Pu-10Zr, a clear negative correlation is found between the Pu content and five parameters describing the amount of porous matter generated. The parameters describing granularity properties, however, showed either negative correlation or nonlinear dependency from fuel composition. The parameters describing the amount showed a positive correlation with fuel axial expansion, while granularity parameters showed a negative correlation with axial expansion. The dependency on cladding temperature was found to be weak. A positive correlation is demonstrated for volume parameters and fuel burnup. In general, reported results confirm and validate findings published in previous studies using a much larger number of pins and automated processing techniques, which easily lend themselves to reproducibility, thus avoiding subjective bias.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Spectral fluctuations in the Sachdev-Ye-Kitaev model

We present a detailed quantitative analysis of spectral correlations in the Sachdev-Ye-Kitaev (SYK) model. We find that the deviations from universal Random Matrix Theory (RMT) behavior are due to a small number of long-wavelength fluctuations (of the order of the number of Majorana fermions N ) from one realization of the ensemble to the next one. These modes can be parameterized effectively in terms of Q-Hermite orthogonal polynomials, the main contribution being due to scale fluctuations for which we give a simple analytical estimate. Our numerical results for N = 32 show that only the lowest eight polynomials are needed to eliminate the nonuniversal part of the spectral fluctuations. The covariance matrix of the coefficients of this expansion can be obtained analytically from low-order double-trace moments. We evaluate the covariance matrix of the first six moments and find that it agrees with the numerics. We also analyze the spectral correlations in terms of a nonlinear σ-model, which is derived through a Fierz transformation, and evaluate the one and two-point spectral correlation functions to two-loop order. The wide correlator is given by the sum of the universal RMT result and corrections whose lowest-order term corresponds to scale fluctuations. However, the loop expansion of the σ-model results in an ill-behaved expansion of the resolvent, and it gives universal RMT fluctuations not only for q= 4 or higher even q-body interactions, but also for the q= 2 SYK model albeit with a much smaller Thouless energy while the correct result in this case should have been Poisson statistics. In our numerical studies we analyze the number variance and spectral form factor for N = 32 and q = 4. We show that the quadratic deviation of the number variance for large energies appears as a peak for small times in the spectral form factor. After eliminating the long-wavelength fluctuations, we find quantitative agreement with RMT up to an exponentially large number of level spacings for the number variance or exponentially short times in the case of the spectral form factor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Deconvolving the roles of E × B shear and pedestal structure in the energy confinement quality of super H-mode experiements

Analysis of “super H-mode” experiments on DIII-D has put forward that high plasma toroidal rotation, not high pedestal, plays the essential role in achieving energy confinement quality H 98y2 >>1. Recently, super H-mode experiments with variable input torque have confirmed that high rotation shear discharges have very high levels of H 98y2 (>1.5), independent of the pedestal height, and that high pedestal discharges with low rotation shear have levels of H 98y2 only slightly above 1 (≤1.2). Although some increase in stored energy with higher pedestal occurs, the energy confinement quality mainly depends on the toroidal rotation shear, which varies according to different levels of injected neutral beam torque per particle. Here, quasi-linear gyrofluid modeling achieves a good match of the experiment when including the E×B shear; without including plasma rotation, the modeling predicts a confinement quality consistent with the empirical observation of H 98y2 ~1.2 at low rotation. Nonlinear gyrokinetic transport modeling shows that the effect of E×B turbulence stabilization is far larger than other mechanisms, such as the so-called hot-ion stabilization (T i /T e ) effect. Consistent with these experimental and modeling results are previous simulations of the ITER Baseline Scenario using a super H-mode pedestal solution, which showed the potential to exceed the Q=10 target if the pedestal density could be increased above the Greenwald limit. A close look at these simulations reveals that the predicted energy confinement quality is below 1 even at the highest pedestal pressure. The improvement in Q at higher pedestal density is due to the improved fusion power generation at the higher core density associated with higher pedestal density, not to an improved energy confinement quality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Numerical fluid dynamics for FRG flow equations: Zero-dimensional QFTs as numerical test cases. II. Entropy production and irreversibility of RG flows

Herein we demonstrate that the reformulation of renormalization group (RG) flow equations as nonlinear heat equations has severe implications on the understanding of RG flows in general. We demonstrate by explicitly constructing an entropy function for a zero-dimensional Z 2 -symmetric model that the dissipative character of generic nonlinear diffusion equations is also hard-coded in the functional RG equation. This renders RG flows manifestly irreversible, revealing the semigroup property of RG transformations on the level of the flow equation itself. Additionally, we argue that the dissipative character of RG flows, its irreversibility and the entropy production during the RG flow may be linked to the existence of a so-called C– / A-function. In total, this introduces an asymmetry in the so-called RG time—in complete analogy to the thermodynamic arrow of time—and allows for an interpretation of infrared actions as equilibrium solutions of dissipative RG flows equations. The impossibility of resolving microphysics from macrophysics is evident in this framework. Furthermore, we directly link the irreversibility and the entropy production in RG flows to an explicit numerical entropy production, which is manifest in diffusive and non-linear partial differential equations (PDEs) and a standard mathematical tool for the analysis of PDEs. Using exactly solvable zero-dimensional Z 2 -symmetric models, we explicitly compute the (numerical) entropy production related to the total variation nonincreasing property of the PDE during RG flows toward the infrared limit. Finally, we discuss generalizations of our findings and relations to the C– / A-theorem as well as how our work may help to construct truncations of RG flow equations in the future, including numerically stable schemes for solving the corresponding PDEs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Direct uranium enrichment assay in gaseous uranium hexafluoride with laser induced breakdown spectroscopy

A set of eleven U atomic emission lines, between 421.099 and 421.460 nm, were characterized for direct enrichment assay of gaseous UF 6 samples with laser induced breakdown spectroscopy (LIBS). Several iterative, multivariate nonlinear spectral-fitting algorithms were evaluated for their efficacies to extract the enrichment information from the measured LIBS spectra. Furthermore, wavelength-dependent weight factors, which take into consideration of the spectral-line position in the fitting model, the determined isotopic ratio, and the isotopic shift of the emission line, are essential for the spectral-fitting model to function adequately. The analytical accuracies and precision were typically within 0.5% in absolute [ 235 U/( 235 U + 238 U)] ratios.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Pedestal collapse by resonant magnetic perturbations

Pedestal collapse (i.e., the complete loss of the edge transport barrier (ETB)) in DIII-D H-mode plasmas occurs when resonant magnetic perturbations (RMPs) penetrate the steep gradient region at the plasma edge. Normally, RMP driven magnetic islands can occur at the top and bottom of the H-mode pedestal and these islands generate conditions consistent with edge-localized-mode (ELM) suppression and density pump-out, respectively, based on nonlinear two-fluid MHD simulations. In contrast, MHD simulations show that the steep pressure gradient region between the top and bottom of the DIII-D pedestal is generally immune to resonant field penetration due to large local E × B and diamagnetic flows. By this fortuitous circumstance, the edge-transport-barrier and H-mode confinement can be maintained while achieving ELM suppression. However, pedestal collapse can occur in DIII-D when the screening flows are inadequate to prevent field penetration in the steep gradient region of the pedestal. Non-linear two-fluid MHD simulations support the role of resonant field penetration in pedestal collapse for DIII-D H-mode plasmas with weak edge E × B and diamagnetic screening flows. ITER will likely have weaker edge screening flows than present experiments due to its much larger size, making it more susceptible to resonant field penetration in the steep gradient region of the pedestal. Analysis of model ITER equilibria demonstrates that resonant field penetration in the steep pressure gradient region is possible for RMP levels of the order required for ELM suppression. The effect of such penetration on the ITER pedestal will depend sensitively on the resulting degree of island overlap.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Artificial neural networks estimate evapotranspiration for Miscanthus × giganteus as effectively as empirical model but with fewer inputs

Estimating actual evapotranspiration (ET) is particularly crucial for addressing how vegetation affects the water balance of ecosystems. ET estimation can be complex with empirical models due to their many parameters and reliance on aridity. In contrast, artificial neural networks (ANNs) could potentially estimate ET with fewer and more common meteorological parameters. In this study, we trained two ANNs, one using a feed-forward approach (FFN) and the other a nonlinear auto-regressive network (NARX), to predict ET and compared them to the commonly used empirical model Granger and Gray (GG). We trained our models on a nine-year eddy covariance (EC) dataset for Miscanthu s × giganteus ( M . × giganteus ) from Illinois (UIEF), then tested them using out-of-sample data from both UIEF and a different location in Iowa (SABR) to compare the accuracy of FFN, NARX, and GG models in estimating daily ET. A combination of air temperature (T a ) and solar radiation (R s ) was chosen as inputs due to the highest R 2 for FFN (R 2 = 0.79, 0.81, and 0.79 for training, testing, and validation, respectively) and only T a for NARX (R 2 = 0.70 for out-of-sample validation). The predictive power of the FFN model was superior to the NARX and GG models at the UIEF site (R 2 = 0.84, 0.70, and 0.83 for out-of-sample validation, respectively). Our analysis showed that ANN approaches are as accurate as empirical approaches for estimating ET but use fewer inputs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Characterizing the effect of incrementally increasing dry bulb temperature on linear and nonlinear measures of heart rate variability in nonpregnant, mid-gestation, and late-gestation sows

Abstract Characterizing the sow physiological response to an increased heat load is essential for effective heat stress mitigation. The study objective was to characterize the effects of a 400-min heating episode on sow heart rate variability (HRV) at different reproductive stages. HRV is a commonly used noninvasive proxy measure of autonomic function. Twenty-seven sows were enrolled in the study according to their gestation stage at time of selection: 1) nonpregnant (NP; n = 7), 2) mid-gestation (MID; 57.3 ± 11.8 d gestation; n = 11), and 3) late-gestation (LATE; 98.8 ± 4.9 d gestation; n = 8). The HRV data utilized in the study were collected from each pig as the dry bulb temperature in the room increased incrementally from 19.84 ± 2.15 °C to 35.54 ± 0.43 °C (range: 17.1–37.5 °C) over a 400-min period. After data collection, one 5-min set of continuous heart rate data were identified per pig for each of nine temperature intervals (19–20.99, 21–22.99, 23–24.99, 25–26.99, 27–28.99, 29–30.99, 31–32.99, 33–34.99, and 35–36.99 °C). Mean inter-beat interval length (RR), standard deviation of r-r intervals (SDNN), root mean square of successive differences (RMSSD), high frequency spectral power (HF), sample entropy (SampEn), short-term detrended fluctuation analysis (DFAα1), and three measures (%REC, DET, LMEAN) derived from recurrence quantification analysis were calculated for each data set. All data were analyzed using the PROC GLIMMIX procedure in SAS 9.4. Overall, LATE sows exhibited lower RR than NP sows (P < 0.01). The standard deviation of r–r intervals and RMSSD differed between each group (P < 0.01), with LATE sows exhibiting the lowest SDNN and RMSSD and NP sows exhibiting the greatest SDNN and RMSSD. Late-gestation sows exhibited lower HF than both MID and NP sows (P < 0.0001), greater DFA values than NP sows (P = 0.05), and greater DET compared to MID sows (P = 0.001). Late-gestation also sows exhibited greater %REC and LMEAN compared to MID (P < 0.01) and NP sows (all P < 0.01). In conclusion, LATE sows exhibited indicators of greater autonomic stress throughout the heating period compared to MID and NP sows. However, temperature by treatment interactions were not detected as dry bulb increased. Future studies are needed to fully elucidate the effect of gestational stage and increasing dry bulb temperature on sow HRV.

Agriculture↗
DOE OSTI2024

A Comprehensive Analysis of Uncertainties in Warm-Rain Parameterizations in Climate Models Based on In Situ Measurements

Abstract Because of the coarse grid size of Earth system models (ESMs), representing warm-rain processes in ESMs is a challenging task involving multiple sources of uncertainty. Previous studies evaluated warm-rain parameterizations mainly according to their performance in emulating collision–coalescence rates for local droplet populations over a short period of a few seconds. The representativeness of these local process rates comes into question when applied in ESMs for grid sizes on the order of 100 km and time steps on the order of 20–30 min. We evaluate several widely used warm-rain parameterizations in ESM application scenarios. In the comparison of local and instantaneous autoconversion rates, the two parameterization schemes based on numerical fitting to stochastic collection equation (SCE) results perform best. However, because of Jessen’s inequality, their performance deteriorates when grid-mean, instead of locally resolved, cloud properties are used in their simulations. In contrast, the effect of Jessen’s inequality partly cancels the overestimation problem of two semianalytical schemes, leading to an improvement in the ESM-like comparison. In the assessment of uncertainty due to the large time step of ESMs, it is found that the rainwater tendency simulated by the SCE is roughly linear for time steps smaller than 10 min, but the nonlinearity effect becomes significant for larger time steps, leading to errors up to a factor of 4 for a time step of 20 min. After considering all uncertainties, the grid-mean and time-averaged rainwater tendency based on the parameterization schemes is mostly within a factor of 4 of the local benchmark results simulated by SCE.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Evolution of coupled weakly driven waves in a dissipative plasma

The nonlinear collisional dynamics of coupled driven plasma waves in the presence of background dissipation is studied analytically within kinetic theory. Sufficiently near marginal stability, phase space correlations are poorly preserved and time delays become unimportant. The system is then shown to be governed by two first-order coupled autonomous differential equations of cubic order for the wave amplitudes and two complementary first-order equations for the evolution of their phases. That system of equations can be decoupled and further simplified to a single second-order differential equation of Liénard's type for each amplitude. Numerical solutions for this equation are obtained in the general case, while analytic solutions are obtained for special cases in terms of parameters related to the spacing of the resonances of the two waves in frequency space, e.g., wave lengths and oscillation frequencies. These parameters are further analyzed to find classes of quasi-steady saturation and pulsating scenarios. To classify equilibrium points, local stability analysis is applied, and bifurcation conditions are determined. When the two waves saturate at similar amplitude levels, their combined signal is shown to invariably exhibit amplitude beating and phase jumps of nearly π. In conclusion, the obtained analytical results can be used to benchmark simulations and to interpret eigenmode amplitude measurements in fusion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗