Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large deviations”

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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Optimization-Based Fast-Frequency Estimation and Control of Low-Inertia Microgrids

The lack of inertial response from non-synchronous, inverter-based generation in microgrids makes the power system vulnerable to a large rate of change of frequency (ROCOF) and frequency excursions. Energy storage systems (ESSs) can be utilized to provide fast-frequency support to prevent such large excursions in the system. However, fast-frequency support is a power-intensive application that has a significant impact on the ESS lifetime. In this paper, a framework that allows the ESS operator to provide fast-frequency support as a service is proposed. The framework maintains the desired quality-of-service (limiting the ROCOF and frequency) while taking into account the ESS lifetime and physical limits. The framework utilizes moving horizon estimation (MHE) to estimate the frequency deviation and ROCOF from noisy phase-locked loop (PLL) measurements. These estimates are employed by a model predictive control (MPC) algorithm that computes control actions by solving a finite-horizon, online optimization problem. Additionally, this approach avoids oscillatory behavior induced by delays that are common when using low-pass filters as with traditional derivative-based (virtual inertia) controllers. MATLAB/Simulink simulations on a test system from Cordova, Alaska, show the effectiveness of the MHE-MPC approach to reduce frequency deviations and ROCOF of a low-inertia microgrid.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Capturing plume behavior in complex terrain: an overview of the Nevada National Security Site Meteorological Experiment (METEX21)

METEX21 was an atmospheric tracer release experiment executed at the Department of Energy’s Nevada National Security Site (NNSS) in the southwestern U.S to study terrain-induced wind and thermodynamic conditions that influence local-scale (<5-km) plume transport under varying atmospheric forcing conditions. Meteorological observations were collected using 10-m tall meteorological towers, 2-m tall tripods with 3-d sonic anemometers, a 3-m tall eddy covariance flux tower, Doppler profiling lidars, Doppler scanning lidars, weather-balloon launched radiosondes, and a tethered balloon equipped with wind, temperature, and aerosol sensors at heights up to 800 m. A smoke tracer was released along three transects in the horizontal and vertical directions and observed with video cameras, aerosol sensors and lidars (via aerosol backscatter). The observations showed evidence of large-scale/synoptic transience as well as local-scale upslope and downslope flows, along-axis valley flows, recirculation eddies on leeward slopes, and periods of strong shear and veer aloft. The release days were classified as either synoptically-driven or locally-driven, and a single case day is presented in detail for each. Synoptically-forced days show relatively narrow smoke plumes traveling down the valley from north to south (with the predominant wind direction), with little deviation in transport direction regardless of the elevation or ground locations of the smoke releases, except near the presence of leeside recirculation eddies. Locally-forced days exhibit a wider range of plume behavior due to the combination of thermally-induced valley and slope flows, which are often flowing in different cardinal directions, and wind shear found aloft at higher altitudes and elevations. We saw evidence of smoke lofting on top of the mesas due to strong upslope flows on these days. A major finding of this experiment was the effectiveness of scanning lidars to measure 2-dimensional plume transport out to a 2–3 km distance; much farther than could be visibly observed. METEX21 was the first of three planned tracer experiments at NNSS, and future experiments will incorporate multiple tracers to improve individual plume identification so that finer resolution flow details can be attained from these measurements, as well as deploy a larger suite of meteorological instrumentation, including more temperature profiling data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Assessing boundary condition and parametric uncertainty in numerical-weather-prediction-modeled, long-term offshore wind speed through machine learning and analog ensemble

To accurately plan and manage wind power plants, not only does the time-varying wind resource at the site of interest need to be assessed but also the uncertainty connected to this estimate. Numerical weather prediction (NWP) models at the mesoscale represent a valuable way to characterize the wind resource offshore, given the challenges connected with measuring hub-height wind speed. The boundary condition and parametric uncertainty associated with modeled wind speed is often estimated by running a model ensemble. However, creating an NWP ensemble of long-term wind resource data over a large region represents a computational challenge. Here, we propose two approaches to temporally extrapolate wind speed boundary condition and parametric uncertainty using a more convenient setup in which a mesoscale ensemble is run over a short-term period (1 year), and only a single model covers the desired long-term period (20 year). We quantify hub-height wind speed boundary condition and parametric uncertainty from the short-term model ensemble as its normalized across-ensemble standard deviation. Then, we develop and apply a gradient-boosting model and an analog ensemble approach to temporally extrapolate such uncertainty to the full 20-year period, for which only a single model run is available. As a test case, we consider offshore wind resource characterization in the California Outer Continental Shelf. Both of the proposed approaches provide accurate estimates of the long-term wind speed boundary condition and parametric uncertainty across the region (R 2 >0.75), with the gradient-boosting model slightly outperforming the analog ensemble in terms of bias and centered root-mean-square error. At the three offshore wind energy lease areas in the region, we find a long-term median hourly uncertainty between 10 % and 14 % of the mean hub-height wind speed values. Finally, we assess the physical variability in the uncertainty estimates. In general, we find that the wind speed uncertainty increases closer to land. Also, neutral conditions have smaller uncertainty than the stable and unstable cases, and the modeled wind speed in winter has less boundary condition and parametric sensitivity than summer.

17 WIND ENERGY↗
NASA NTRS2003

Toward Immersed Boundary Simulation of High Reynolds Number Flows

In the immersed boundary (IB) method, the surface of an object is reconstructed with forcing terms in the underlying flow field equations. The surface may split a computational cell removing the constraint of the near wall gridlines to be aligned with the surface. This feature greatly simplifies the grid generation process which is cumbersome and expensive in particular for structured grids and complex geometries. The IB method is ideally suited for Cartesian flow solvers. The flow equations written in Cartesian coordinates appear in a very simple form and several numerical algorithms can be used for an efficient solution of the equations. In addition, the accuracy of numerical algorithms is dependent on the underlying grid and it usually deteriorates when the grid deviates from a Cartesian mesh. The challenge for the IB method lies in the representation of the wall boundaries and in providing an adequate near wall flow field resolution. The issue of enforcing no-slip boundary conditions at the immersed surface has been addressed by several authors by imposing a local reconstruction of the solution. Initial work by Verzicco et al. was based on a simple linear, one-dimensional operator and this approach proved to be accurate for boundaries largely aligned with the grid lines. Majumdar et al. used various multidimensional and high order polynomial interpolations schemes. These high order schemes, however, are keen to introduce wiggles and spurious extrema. Iaccarino & Verzicco and Kalitzin & Iaccarino proposed a tri-linear reconstruction for the velocity components and the turbulent scalars. A modified implementation that has proven to be more robust is reported in this paper. The issue of adequate near wall resolution in a Cartesian framework can initially be addressed by using a non-uniform mesh which is stretched near the surface. In this paper, we investigate an unstructured approach for local grid refinement that utilizes Cartesian mesh features. The computation of high Reynolds number wall bounded flows is particularly challenging as it requires the consideration of thin turbulent boundary layers, i.e. near wall regions with large gradients of the flow field variables. For such flows, the representation of the wall boundary has a large impact on the accuracy of the computation. It is also critical for the robustness and convergence of the flow solver.

Kalitzin, Georgi↗
DOE OSTI2023

Evidence for the charge asymmetry in $pp$ → $t\bar{t}$ production at $\sqrt{s}$ = 13 TeV with the ATLAS detector

Inclusive and differential measurements of the top–antitop ($t\bar{t}$) charge asymmetry ${A}_{C}^{t\bar{t}}$ and the leptonic asymmetry ${A}_{C}^{\ell\bar{\ell}}$ are presented in proton–proton collisions at $\sqrt{s}$ = 13 TeV recorded by the ATLAS experiment at the CERN Large Hadron Collider. The measurement uses the complete Run 2 dataset, corresponding to an integrated luminosity of 139 fb –1 , combines data in the single-lepton and dilepton channels, and employs reconstruction techniques adapted to both the resolved and boosted topologies. A Bayesian unfolding procedure is performed to correct for detector resolution and acceptance effects. The combined inclusive $t\bar{t}$ charge asymmetry is measured to be ${A}_{C}^{t\bar{t}}$ = 0.0068 ± 0.0015, which differs from zero by 4.7 standard deviations. Differential measurements are performed as a function of the invariant mass, transverse momentum and longitudinal boost of the $t\bar{t}$ system. Both the inclusive and differential measurements are found to be compatible with the Standard Model predictions, at next-to-next-to-leading order in quantum chromodynamics perturbation theory with next-to-leading-order electroweak corrections. The measurements are interpreted in the framework of the Standard Model effective field theory, placing competitive bounds on several Wilson coefficients.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A Modified Sand’s Time Incorporating Li-Ion Transport Across the SEI: Basis for Understanding Li Dendrite Formation and Li-Metal Battery Electrolyte Selection

Abstract Understanding the initiation of lithium dendrites remains elusive, largely due to the intricate role of the solid electrolyte interphase (SEI) which forms on the Li surface during electrodeposition. Many studies have utilized the classical Sand’s equation to estimate the onset time when lithium dendrites begin to form. The Sand’s equation provides the time when the cation (Li+) concentration at the electrode-electrolyte interface approaches zero under diffusion-limited conditions in galvanostatic Li electrodeposition. However, recent experimental studies have revealed that the observed lithium dendrite onset time deviates considerably from the Sand’s time. Here, we show that this deviation from classical theory is likely due to the transport of Li+ ions through the SEI - a transport limitation that is much more dominant in controlling dendrite formation. We develop a ‘modified’ Sand's equation, incorporating the SEI layer and the diffusional transport across it to predict Li dendrite onset times. To validate this approach, we conducted Li electrodeposition experiments at various current densities using two distinct organic electrolytes. Analysis of the results demonstrates that the modified Sand's equation provides a more accurate prediction of dendrite onset times, highlighting the importance of incorporating SEI into transport models of Li plating in next-generation rechargeable Li-metal batteries.

Ma, Yuanman (ORCID:0000000200444811)↗
NASA NTRS1997

Effects of static fingertip loading on carpal tunnel pressure

The purpose of this study was to explore the relationship between carpal tunnel pressure and fingertip force during a simple pressing task. Carpal tunnel pressure was measured in 15 healthy volunteers by means of a saline-filled catheter inserted percutaneously into the carpal tunnel of the nondominant hand. The subjects pressed on a load cell with the tip of the index finger and with 0, 6, 9, and 12 N of force. The task was repeated in 10 wrist postures: neutral; 10 and 20 degrees of ulnar deviation; 10 degrees of radial deviation; and 15, 30, and 45 degrees of both flexion and extension. Fingertip loading significantly increased carpal tunnel pressure for all wrist angles (p = 0.0001). Post hoc analyses identified significant increase (p < 0.05) in carpal tunnel pressure between unloaded (0 N) and all loaded conditions, as well as between the 6 and 12 N load conditions. This study demonstrates that the process whereby fingertip loading elevates carpal tunnel pressure is independent of wrist posture and that relatively small fingertip loads have a large effect on carpal tunnel pressure. It also reveals the response characteristics of carpal tunnel pressure to fingertip loading, which is one step in understanding the relationship between sustained grip and pinch activities and the aggravation or development of median neuropathy at the wrist.

NASA Center ARC↗
NASA NTRS2011

Emissions of Black Carbon, Organic, and Inorganic Aerosols From Biomass Burning in North America and Asia in 2008

Reliable assessment of the impact of aerosols emitted from boreal forest fires on the Arctic climate necessitates improved understanding of emissions and the microphysical properties of carbonaceous (black carbon (BC) and organic aerosols (OA)) and inorganic aerosols. The size distributions of BC were measured by an SP2 based on the laser-induced incandescence technique on board the DC-8 aircraft during the NASA ARCTAS campaign. Aircraft sampling was made in fresh plumes strongly impacted by wildfires in North America (Canada and California) in summer 2008 and in those transported from Asia (Siberia in Russia and Kazakhstan) in spring 2008. We extracted biomass burning plumes using particle and tracer (CO, CH3CN, and CH2Cl2) data. OA constituted the dominant fraction of aerosols mass in the submicron range. The large majority of the emitted particles did not contain BC. We related the combustion phase of the fire as represented by the modified combustion efficiency (MCE) to the emission ratios between BC and other species. In particular, we derived the average emission ratios of BC/CO = 2.3 +/- 2.2 and 8.5 +/- 5.4 ng/cu m/ppbv for BB in North America and Asia, respectively. The difference in the BC/CO emission ratios is likely due to the difference in MCE. The count median diameters and geometric standard deviations of the lognormal size distribution of BC in the BB plumes were 136-141 nm and 1.32-1.36, respectively, and depended little on MCE. These BC particles were thickly coated, with shell/core ratios of 1.3-1.6. These parameters can be used directly for improving model estimates of the impact of BB in the Arctic.

Kondo, Y.↗
NASA NTRS2019

Tropospheric Ozone Assessment Report: Tropospheric Ozone From 1877 to 2016, Observed Levels, Trends and Uncertainties

From the earliest observations of ozone in the lower atmosphere in the 19th century, both measurement methods and the portion of the globe observed have evolved and changed. These methods have different uncertainties and biases, and the data records differ with respect to coverage (space and time), information content, and representativeness. In this study, various ozone measurement methods and ozone datasets are reviewed and selected for inclusion in the historical record of background ozone levels, based on relationship of the measurement technique to the modern UV absorption standard, absence of interfering pollutants, representativeness of the well-mixed boundary layer and expert judgement of their credibility. There are significant uncertainties with the 19th and early 20th-century measurements related to interference of other gases. Spectroscopic methods applied before 1960 have likely underestimated ozone by as much as 11% at the surface and by about 24% in the free troposphere, due to the use of differing ozone absorption coefficients. There is no unambiguous evidence in the measurement record back to 1896 that typical mid-latitude background surface ozone values were below about 20 nmol mol–1, but there is robust evidence for increases in the temperate and polar regions of the northern hemisphere of 30–70%, with large uncertainty, between the period of historic observations, 1896–1975, and the modern period (1990–2014). Independent historical observations from balloons and aircraft indicate similar changes in the free troposphere. Changes in the southern hemisphere are much less. Regional representativeness of the available observations remains a potential source of large errors, which are difficult to quantify. The great majority of validation and intercomparison studies of free tropospheric ozone measurement methods use ECC ozonesondes as reference. Compared to UV-absorption measurements they show a modest (~1–5% ±5%) high bias in the troposphere, but no evidence of a change with time. Umkehr, lidar, and FTIR methods all show modest low biases relative to ECCs, and so, using ECC sondes as a transfer standard, all appear to agree to within one standard deviation with the modern UV-absorption standard. Other sonde types show an increase of 5–20% in sensitivity to tropospheric ozone from 1970–1995. Biases and standard deviations of satellite retrieval comparisons are often 2–3 times larger than those of other free tropospheric measurements. The lack of information on temporal changes of bias for satellite measurements of tropospheric ozone is an area of concern for long-term trend studies.

UV absorption standard↗
DOE OSTI2021

Evolution of the superconducting properties from binary to ternary APC-Nb 3 Sn wires

We present a study conducted on binary Tube Type and ternary powder-in-tube Nb 3 Sn wires manufactured using the artificial pinning centres-internal oxidation method. All the specimens are doped with Zr: oxide nano-particles of this element are responsible for the pinning improvement, both by refining the A-15 grain-size and their own point-pinning contribution. Low-field J c magnetometry confirms that the hadron-hadron Future Circular Collider (FCC-hh) specifications are met by one ternary doped-sample. The differences in microstructure were assessed by scanning electron microscopy/transmission electron microscopy to clarify the reasons for the pinning improvement between the two generations. The deviations from the Dew Hughes model are also discussed, underlying some non-linear addition due to competition between the two pinning mechanisms. Finally, we show how the introduction of Ta as a ternary addition influences the A-15 phase by focusing on the radial inhomogeneities, evaluating the T c distribution and Sn composition gradients. The latter are used to model the currents, enabling us to evaluate the individual weights of the pinning mechanisms and their absolute contributions at the High Luminosity-Large Hadron Collider and FCC-hh dipoles operational points.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Quantifying rupture characteristics of microearthquakes in the Parkfield Area using a high-resolution borehole network

It is well known that large earthquakes often exhibit significant rupture complexity such as well separated subevents. With improved recording and data processing techniques, small earthquakes have been found to exhibit rupture complexity as well. Studying these small earthquakes offers the opportunity to better understand the possible causes of rupture complexities. Specifically, if they are random or are related to fault properties. We examine microearthquakes (M < 3) in the Parkfield, California, area that are recorded by a high-resolution borehole network. We quantify earthquake complexity by the deviation of source time functions and source spectra from simple circular (omega-square) source models. We establish thresholds to declare complexity, and find that it can be detected in earthquakes larger than magnitude 2, with the best resolution above M2.5. Comparison between the two approaches reveals good agreement (>90 per cent), implying both methods are characterizing the same source complexity. For the two methods, 60–80 per cent (M 2.6–3) of the resolved events are complex depending on the method. The complex events we observe tend to cluster in areas of previously identified structural complexity; a larger fraction of the earthquakes exhibit complexity in the days following the M w 6 2004 Parkfield earthquake. Ignoring the complexity of these small events can introduce artefacts or add uncertainty to stress drop measurements. Focusing only on simple events however could lead to systematic bias, scaling artefacts and the lack of measurements of stress in structurally complex regions.

58 GEOSCIENCES↗
DOE OSTI2024

Analysis of DESI × DES using the Lagrangian effective theory of LSS

In this work we use Lagrangian perturbation theory to analyze the harmonic space galaxy clustering signal of the Bright Galaxy Survey (BGS) and luminous red galaxies (LRGs) targeted by the dark energy spectroscopic instrument (DESI), combined with the galaxy-galaxy lensing signal measured around these galaxies using Dark Energy Survey Year 3 source galaxies. The BGS and LRG galaxies are extremely well characterized by DESI spectroscopy and, as a result, lens galaxy redshift uncertainty and photometric systematics contribute negligibly to the error budget of our “2 × 2-point” analysis. On the modeling side, this work represents the first application of the SPINOSAURUS code, implementing an effective field theory model for galaxy intrinsic alignments, and we additionally introduce a new scheme ( MAIAR ) for marginalizing over the large uncertainties in the redshift evolution of the intrinsic alignment signal. Furthermore, this is the first application of a hybrid effective field theory model for galaxy bias based on the A emulus 𝜈 simulations. Our main result is a measurement of the amplitude of the lensing signal, 𝑆 8 = 𝜎 8 ⁢(Ω 𝑚 /0.3) 0.5 = 0.85⁢0$^{+0.042}_{−0.050}$, consistent with values of this parameter derived from the primary cosmic microwave background. This constraint is artificially improved by a factor of 51% if we assume a more standard, but restrictive parametrization for the redshift evolution and sample dependence of the intrinsic alignment signal, and 63% if we additionally assume the nonlinear alignment model. Furthermore, we show that when fixing the cosmological model to the best-fit values from Planck PR4 there is > 5⁢𝜎 evidence for a deviation of the evolution of the intrinsic alignment signal from the functional form that is usually assumed in cosmic shear and galaxy-galaxy lensing studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The role of driver head pose dynamics and instantaneous driving in safety critical events: Application of computer vision in naturalistic driving

This paper investigates the role of driver behavior especially head pose dynamics in safety–critical events (SCEs). Using a large dataset collected in a naturalistic driving study, this paper analyzes the head pose dynamics and driving behavior in moments leading up to crashes or near-crashes. The study uses advanced computer vision and mixed logit modeling techniques to identify patterns and relationships between drivers’ head pose dynamics and crash involvement. The results suggest that driver-head pose dynamics, especially poses that indicate distraction and movement volatility, are important factors that can contribute to undesirable safety outcomes. Marginal effects show that angular deviation for head pose dynamics indicated by yaw, pitch and roll increase the likelihood of crash intensity by 4.56%, 4.92% and 8.26% respectively. Furthermore, traffic flow and lane changing also contribute to increase in likelihood of crash intensity. These findings provide new insights into pre-crash factors, especially human factors and safety–critical events. The study highlights the importance of considering human factors in designing driver assistance systems and developing safer vehicles. This research contributes by examining naturalistic driving data at the microscopic level with early detection of behaviors that lead to SCEs and provides a basis for future research on automation.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Measurement of the CP properties of Higgs boson interactions with $\tau$-leptons with the ATLAS detector

A study of the charge conjugation and parity ($CP$) properties of the interaction between the Higgs boson and $\tau$-leptons is presented. The study is based on a measurement of $CP$-sensitive angular observables defined by the visible decay products of $\tau$-leptons produced in Higgs boson decays. The analysis uses 139 fb -1 of proton–proton collision data recorded at a centre-of-mass energy of $\sqrt{s}$ = 13 TeV with the ATLAS detector at the Large Hadron Collider. Contributions from $\tau$-violating interactions between the Higgs boson and $\tau$-leptons are described by a single mixing angle parameter Φ $\tau$ in the generalised Yukawa interaction. Without constraining the $H$ $\rightarrow$ $\tau$$\tau$ signal strength to its expected value under the Standard Model hypothesis, the mixing angle Φ $\tau$ is measured to be 9° ± 16°, with an expected value of 0° ± 28° at the 68% confidence level. The pure $CP$-odd hypothesis is disfavoured at a level of 3.4 standard deviations. The results are compatible with the predictions for the Higgs boson in the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2023

A Spaceborne Demonstration of P-Band Signals-of-Opportunity (SoOp) Reflectometry

Land-reflected signals from a geosynchronous communication satellite broadcasting in P-band (367.5 MHz) were captured in low Earth orbit (LEO) using a simple dipole antenna. A delay-Doppler map (DDM) was generated through autocorrelation. Estimates of the specular point delay were obtained from the lag of the second peak in the DDM with a bias of 239.4 m and a standard deviation of 44 m (12 m over a frozen lake) with respect to a predicted orbit model. Relative magnitudes of the first and second DDM peaks fell within the range of values predicted using dielectric models for the frozen ground and lake. Finally, retrievals of surface reflection coefficient were generated using a range of realistic values for the transmitter link budget G/T , and these also fell within the range of possible values for the antenna gain pattern. Given the lack of calibration and the large uncertainties in the receiver orbit and attitude, this agreement is sufficient to conclude a successful demonstration of the fundamental principle of single-antenna reflectometry in P-band. P-band reflectometry may offer a new approach to remote sensing of sub-canopy and root-zone soil moisture.

Delays↗
NASA NTRS2015

Comparison of High-Order and Low-Order Methods for Large-Eddy Simulation of a Compressible Shear Layer

The objective of this work is to compare a high-order solver with a low-order solver for performing Large-Eddy Simulations (LES) of a compressible mixing layer. The high-order method is the Wave-Resolving LES (WRLES) solver employing a Dispersion Relation Preserving (DRP) scheme. The low-order solver is the Wind-US code, which employs the second-order Roe Physical scheme. Both solvers are used to perform LES of the turbulent mixing between two supersonic streams at a convective Mach number of 0.46. The high-order and low-order methods are evaluated at two different levels of grid resolution. For a fine grid resolution, the low-order method produces a very similar solution to the highorder method. At this fine resolution the effects of numerical scheme, subgrid scale modeling, and filtering were found to be negligible. Both methods predict turbulent stresses that are in reasonable agreement with experimental data. However, when the grid resolution is coarsened, the difference between the two solvers becomes apparent. The low-order method deviates from experimental results when the resolution is no longer adequate. The high-order DRP solution shows minimal grid dependence. The effects of subgrid scale modeling and spatial filtering were found to be negligible at both resolutions. For the high-order solver on the fine mesh, a parametric study of the spanwise width was conducted to determine its effect on solution accuracy. An insufficient spanwise width was found to impose an artificial spanwise mode and limit the resolved spanwise modes. We estimate that the spanwise depth needs to be 2.5 times larger than the largest coherent structures to capture the largest spanwise mode and accurately predict turbulent mixing.

large eddy simulations↗
NASA NTRS2015

Comparison of High-Order and Low-Order Methods for Large-Eddy Simulation of a Compressible Shear Layer

The objective of this work is to compare a high-order solver with a low-order solver for performing large-eddy simulations (LES) of a compressible mixing layer. The high-order method is the Wave-Resolving LES (WRLES) solver employing a Dispersion Relation Preserving (DRP) scheme. The low-order solver is the Wind-US code, which employs the second-order Roe Physical scheme. Both solvers are used to perform LES of the turbulent mixing between two supersonic streams at a convective Mach number of 0.46. The high-order and low-order methods are evaluated at two different levels of grid resolution. For a fine grid resolution, the low-order method produces a very similar solution to the high-order method. At this fine resolution the effects of numerical scheme, subgrid scale modeling, and filtering were found to be negligible. Both methods predict turbulent stresses that are in reasonable agreement with experimental data. However, when the grid resolution is coarsened, the difference between the two solvers becomes apparent. The low-order method deviates from experimental results when the resolution is no longer adequate. The high-order DRP solution shows minimal grid dependence. The effects of subgrid scale modeling and spatial filtering were found to be negligible at both resolutions. For the high-order solver on the fine mesh, a parametric study of the spanwise width was conducted to determine its effect on solution accuracy. An insufficient spanwise width was found to impose an artificial spanwise mode and limit the resolved spanwise modes. We estimate that the spanwise depth needs to be 2.5 times larger than the largest coherent structures to capture the largest spanwise mode and accurately predict turbulent mixing.

large eddy simulations↗
DOE OSTI2024

Observation of $t\overline{t}$ production in the lepton+jets and dilepton channels in p+Pb collisions at $\sqrt{s_{\textrm{NN}}}$ = 8.16 TeV with the ATLAS detector

This paper reports the observation of top-quark pair production in proton-lead collisions in the ATLAS experiment at the Large Hadron Collider. The measurement is performed using 165 nb -1 of p+Pb data collected at $\sqrt{^sNN}$ = 8.16 TeV in 2016. Events are categorised in two analysis channels, consisting of either events with exactly one lepton (electron or muon) and at least four jets, or events with two opposite-charge leptons and at least two jets. In both channels at least one b-tagged jet is also required. Top-quark pair production is observed with a significance over five standard deviations in each channel. The top-quark pair production cross-section is measured to be σ $t\overline{t}$ =58.1 ± 2.0(stat.)$^{+4.8}_{-4.4}$(syst.) nb, with a total uncertainty of 9%. In addition, the nuclear modification factor is measured to be R pA = 1.090 ± 0.039(stat.)$^{+0.094}_{-0.087}$(syst.). The measurements are found to be in good agreement with theory predictions involving nuclear parton distribution functions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗