Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “predictive analytics”

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 217 records · Page 12

A vortex sheet based analytical model of the curled wake behind yawed wind turbines

Motivated by the need for compact descriptions of the evolution of non-classical wakes behind yawed wind turbines, we develop an analytical model to predict the shape of curled wakes. Interest in such modelling arises due to the potential of wake steering as a strategy for mitigating power reduction and unsteady loading of downstream turbines in wind farms. We first estimate the distribution of the shed vorticity at the wake edge due to both yaw offset and rotating blades. By considering the wake edge as an ideally thin vortex sheet, we describe its evolution in time moving with the flow. Vortex sheet equations are solved using a power series expansion method, and an approximate solution for the wake shape is obtained. The vortex sheet time evolution is then mapped into a spatial evolution by using a convection velocity. Apart from the wake shape, the lateral deflection of the wake including ground effects is modelled. Our results show that there exists a universal solution for the shape of curled wakes if suitable dimensionless variables are employed. For the case of turbulent boundary layer inflow, the decay of vortex sheet circulation due to turbulent diffusion is included. Finally, we modify the Gaussian wake model by incorporating the predicted shape and deflection of the curled wake, so that we can calculate the wake profiles behind yawed turbines. Model predictions are validated against large-eddy simulations and laboratory experiments for turbines with various operating conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Unified Analytical Method to Quantify Three Types of Fast Frequency Response from Inverter-Based Resources

With more inverter-based resources (IBRs), our power systems have lower frequency nadirs following N-1 contingencies, and undesired under-frequency load shedding (UFLS) can occur. To address this challenge, IBRs can be programmed to provide at least three types of fast frequency response (FFR), e.g., step response, proportional response (P/f droop response), and derivative response (synthetic inertia). However, these heterogeneous FFR challenge the study of power system frequency dynamics. Thus, this paper develops an analytical frequency nadir prediction method that allows for the consideration of all three potential forms of FFR provided by IBRs. The proposed method provides fast and accurate frequency nadir estimation after N-1 generation tripping contingencies. Our method is grounded on the closed-form solution for the frequency nadir, which is solved from the second-order system frequency response model considering the governor dynamics and three types of FFR. The simulation results in the IEEE 39-bus system with different types of FFR demonstrate that the proposed method provides an accurate and fast prediction of the frequency nadir under various disturbances.

fast frequency response↗

Black hole-neutron star binaries near neutron star disruption limit in the mass regime of event GW230529

In May 2023, the LIGO-Virgo-KAGRA (LVK) Collaboration reported the likely black hole-neutron star (BHNS) merger GW230529_181500. The signal was observed with high significance in only one detector, limiting constraints on the black hole spin and motivating our study of disruption in this mass regime. That event is expected to be the merger of a 2.5–4.5 M ⊙ primary with a secondary compact object of mass between 1.2–2.0 M ⊙ . This makes it the first BHNS merger with a significant potential for the production of electromagnetic (EM) counterparts, and provides further evidence for compact objects existing within the suspected lower mass gap. To produce post-merger EM transients, the component of the black hole spin aligned with the orbital angular momentum must be sufficiently high, allowing the neutron star to be tidally disrupted. The disrupting BHNS binary may then eject a few percent of a solar mass of matter, leading to an observable kilonova driven by radioactive decays in ejecta, and/or a compact-binary gamma-ray burst (cbGRB) resulting from the formation of an accretion disk and relativistic jet. Determining which mergers lead to disruption of the neutron star is necessary to predict the prevalence of EM signals from BHNS mergers, yet most BHNS simulations so far have been performed far from the minimum spin required for tidal disruption. Here, we use the Spectral Einstein Code to explore the behavior of BHNS mergers in a mass range consistent with GW230529_181500 close to that critical spin, and compare our results against the mass remnant model currently used by the LVK Collaboration to predict the probability of tidal disruption. Furthermore, our numerical results reveal the emergence of non-zero accretion disks even below the predicted NS disruption limit, of low mass but capable of powering cbGRBs. Our results also demonstrate that the remnant mass model underpredicts the disk mass for the DD2 equation of state, while they are within expected modeling errors for SFHo. The disruption limit itself, however, is not found to significantly differ from the predictions of the analytical model, unless remnant masses M rem ≲ 0.001 M ⊙ prove interesting observationally. In all of our simulations, any kilonova signal would be dim and most likely dominated by post-merger disk outflows.

GW230529_181500↗

Solidification dynamics in metal additive manufacturing: analysis of model assumptions

Solidification characteristics in the meltpool drive the process-microstructure relationship which helps dictate the material properties of as-built parts in additive manufacturing; therefore, being able to accurately and quickly simulate the size, shape, and solidification characteristics in the melt pool is of great interest to the field. In this study, we investigate various important physical phenomena (dynamic material properties, fluid-flow, radiation and vaporization) which can either be included or neglected in a continuum finite volume model (FVM) and their effect on the solidification conditions. Additionally, since the simplest form of such a model (conduction only) has an analytic solution which is much faster, its viability is also considered. Since the inclusion of some of these physical phenomena will inherently change the net energy input as well as the amount of energy needed to achieve melting of a control volume, each set of included phenomena had an effective absorption efficiency which was calibrated to closely match the dimensions of the melt pool to that of the ground truth data. The ground truth data for this study was defined to be the output of the FVM which included all the physical phenomena (OF). This study then goes on to compare the effects on solidification conditions each of these calibrated models has. It was found that most of the change in solidification conditions comes from the inclusion of latent heat. A posterior correlation factor (PCF) is then introduced to enable an analytic model to predict similar solidification conditions to OF model.

36 MATERIALS SCIENCE↗

The effect of curvature induced broken potential vorticity conservation on drift wave turbulences

Abstract The correlation theory of turbulence suppression (Zhang and Mahajan 1993 Phys. Fluids B5 2000) by velocity shear was constructed by invoking the ansatz of potential vorticity conservation (PVC) that holds for relatively simple (slab) models of drift wave turbulence. It is, therefore, surprising that a detailed modern simulation of the H mode pedestal, using the gyrokinetic code GENE (Hatch et al 2018 Plasma Phys. Control. Fusion 60 084003), found ‘striking agreement’ with the predictions of the analytic model. To understand the reasons for this remarkable agreement, an extended theory that contains finite (magnetic) curvature, and which does not conserve potential vorticity, is developed and ‘solved’ by calculating an inhomogeneous Green function reflecting the fact that the new system has a potential vorticity source. It is, then, demonstrated that the effect of the broken PVC is insignificant for the normal operation parameters in tokamaks; the correction due to curvature is at the order of 2 L 0 / R , where L 0 is the scale length of the local gradient and R is the major radius. The excellent agreement between simulation and slab model 1993 theory is therefore not accidental; the latter can be applied with confidence to the tokamak pedestal.

Physics↗

Photon phase-space dynamics in a plasma wakefield accelerator

Frequency up-shifting of laser light in a beam-driven plasma wakefield has the potential to provide high-intensity sources of short wavelength radiation. Simulations have demonstrated that a laser pulse can undergo large frequency shifts, limited only by the drive beam energy, when the plasma density is tailored to match the accelerating phase of the wake to the group velocity of the pulse. Here, we study the dynamical evolution of photons in the 1D1P phase-space vicinity of the plasma wake-phase matching condition. Numerical calculations using a photon kinetic model are validated by direct comparison with 1D and quasi-3D particle-in-cell simulations. These calculations form the basis of a linear theory of the photon dynamics which reveals several important results, including scalings for the properties of the witness pulse and a self-similar solution for the photon phase-space dynamics. One prediction of the analytic theory is that the pulse can be compressed indefinitely, though the amount of compression would ultimately be limited by practical constraints. These results suggest that photon acceleration can provide a novel source of sub-femtosecond, short wavelength radiation.

XUV generation↗

Comparison between MAST-U conventional and Super-X configurations through SOLPS-ITER modelling

Abstract MAST-U has recently started operating with a Super-X divertor, designed to increase total flux expansion and neutral trapping, both predicted through simple analytic models and SOLPS calculations to reduce the plasma and impurity density detachment thresholds. In this study, utilising the SOLPS-ITER code, we are quantifying the possible gain allowed by the MAST-U Super-X and neutral baffling geometry, in terms of access to detachment. We show that a significant reduction of the upstream density detachment threshold (up to a factor 1.6) could be achieved in MAST-U, for the Super-X, as opposed to conventional divertor geometry, mainly through an increased total flux expansion, neutral trapping being found very similar between the different configurations. We also show that variations of the strike-point angle are complex to interpret in such a tightly baffled geometry, and that a case in which the target normal points more towards the separatrix does not necessarily imply a lower detachment threshold. As in previous calculations for TCV, we quantify the role of neutral effects through developing and applying a quantitative definition of neutral trapping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Infernal instabilities in negative-triangularity plasmas with negative central shear

Abstract A systematic numerical investigation is carried out to understand magnetohydrodynamic stability of the ideal infernal-kink instability in tokamak plasmas with both negative triangularity (neg-D) shaping and negative central shear for the equilibrium safety factor profile. The latter is motivated by the desire to form the internal transport barrier in the neg-D configuration, which is known to have difficulty in forming the edge transport barrier. The infernal-kink mode is generally found to be more unstable in neg-D plasmas as compared to their positive D-shaped (pos-D) counterpart. This is mainly due to less favorable (or even unfavorable) average magnetic curvature near the radial location of the minimum safety factor ( q min ) as compared to the pos-D configuration. The larger Shafranov shift associated with the neg-D shape helps the mode stabilization but is not sufficient to overcome the destabilizing effect due to bad curvature. Strong poloidal mode coupling due to plasma shaping (toroidicity, elongation, triangularity, etc.) helps explain the slight shift with respect to that predicted by the analytic theory of the peak location of the computed mode growth versus q min .

Zhang, X. M.↗

Heavy flavor jet substructure at lepton colliders

We provide a detailed analysis of event-shape observables, namely the energy correlation function and jet angularity, for heavy-flavor jets produced in electron-positron collisions, focusing on quantum chromodynamics (QCD) interactions. Using modern jet substructure techniques, we investigate the dead-cone effect, where QCD radiation is suppressed around a heavy quark within an angle proportional to its mass. Our analysis achieves next-to-leading logarithmic accuracy, combined with partial fixed-order contributions, to improve the description of the transition near the dead-cone threshold. To ensure a comprehensive perspective, we compare our analytical results with predictions from the pythia, herwig, and sherpa Monte Carlo simulations at past and future lepton colliders.

Perturbative QCD↗

Buckling-induced interaction between circular inclusions in an infinite thin plate

Design of slender artificial materials and morphogenesis of thin biological tissues typically involve stimulation of isolated regions (inclusions) in the growing body. These inclusions apply internal stresses on their surrounding areas that are ultimately relaxed by out-of-plane deformation (buckling). Here, we utilize the Föppl-von Kármán model to analyze the interaction between two circular inclusions in an infinite plate that their centers are separated a distance of 2ℓ. In particular, we investigate a region in phase space where buckling occurs at a narrow transition layer of length ℓ D around the radius of the inclusion, $R (ℓ_D \ll R)$. We show that the latter length scale defines two regions within the system, the close separation region, $ℓ – R \sim ℓ_D$, where the transition layers of the two inclusions approximately coalesce, and the far separation region, $ℓ – R \gg ℓ_D$. While the interaction energy decays exponentially in the latter region, $E_{\text{int}} ∝ e^{– (ℓ – R)/ℓ_D}$, it presents nonmonotonic behavior in the former region. While this exponential decay is predicted by our analytical analysis and agrees with the numerical observations, the close separation region is treated only numerically. In particular, we utilize the numerical investigation to explore two different scenarios within the final configuration: The first where the two inclusions buckle in the same direction (up-up solution) and the second where the two inclusions buckle in opposite directions (up-down solution). We show that the up-down solution is always energetically favorable over the up-up solution. In addition, we point to a curious symmetry breaking within the up-down scenario; we show that this solution becomes asymmetric in the close separation region.

42 ENGINEERING↗

Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma

Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Learning Global Proliferation Expertise Evolution Using AI-Driven Analytics and Public Information

Detecting and anticipating global proliferation expertise and capability evolution from unstructured, noisy, and incomplete public data streams is a highly desired, but extremely challenging task. Here, in this article, we present our pioneering data-driven approach to support the non-proliferation mission to detect and explain the evolution of proliferation expertise and capability development globally from terabytes of publicly available information (PAI), focusing on our knowledge extraction pipeline and descriptive analytics. We first discuss how we fuse nine open-source data streams, including multilingual data, to convert 4 TB of unstructured data to structured knowledge and encode dynamically evolving proliferation expertise representations—content and context graphs. For this, we rely on natural language processing (NLP) and deep learning (DL) models to perform information extraction, topic modeling, and distributed text representation (aka embedding) learning. We then present interactive, usable, and explainable descriptive analytics to refine domain knowledge and present it in a human-understandable form. Finally, we introduce future work avenues that will leverage our dynamic knowledge representations and descriptive analytics to enable predictive and prescriptive inferences to achieve real-time domain understanding and contextual reasoning about global proliferation expertise and capability evolution.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Extension of the bright high-harmonic photon energy range via nonadiabatic critical phase matching

The concept of critical ionization fraction has been essential for high-harmonic generation, because it dictates the maximum driving laser intensity while preserving the phase matching of harmonics. In this work, we reveal a second, nonadiabatic critical ionization fraction, which substantially extends the phase-matched harmonic energy, arising because of the strong reshaping of the intense laser field in a gas plasma. We validate this understanding through a systematic comparison between experiment and theory for a wide range of laser conditions. In particular, the properties of the high-harmonic spectrum versus the laser intensity undergoes three distinctive scenarios: (i) coincidence with the single-atom cutoff, (ii) strong spectral extension, and (iii) spectral energy saturation. We present an analytical model that predicts the spectral extension and reveals the increasing importance of the nonadiabatic effects for mid-infrared lasers. These findings are important for the development of high-brightness soft x-ray sources for applications in spectroscopy and imaging.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A New Analytical Expression for Estimating the Adiabatic Temperature Rise in Lithium-Ion Batteries During High-Power Pulses

A coupled, thermal-electrochemical model is used to explain why the Joule heating assumption (i.e., I 2 R) does not provide a good representation of the temperature rise during high power pulses in lithium-ion batteries, even in cases where the reversible heat generation can be neglected. The poor agreement occurs because the internal resistance changes during the pulse due to the opposing effects of mass transport in the electrolyte, which raises the resistance, and heat generation (temperature rise), which lowers the resistance. These insights are used to propose a new analytical expression for predicting the temperature rise during adiabatic pulses with limited experimental and physical data. The expression accounts for thermal effects using a Taylor series expansion of an Arrhenius-type equation. It accounts for transport effects using an approximate solution to the one-dimensional diffusion equation. The new expression is shown to accurately estimate the simulated, adiabatic temperature rise across a range of loadings (1 to 4 mAh cm -2 ) and C-rates (1C to 10C) for cells containing LiNi 0.5 Mn 0.3 Co 0.2 O 2 positive electrodes and graphite negative electrodes. In conclusion, it is also shown to accurately estimate the experimental, adiabatic temperature rise measured for cells with LiFePO 4 positive electrodes and mesocarbon negative electrodes with minimal changes to the fitted parameters.

25 ENERGY STORAGE↗

Modeling the Vertical Response of an Aeroseismometer Balloon System

Aeroseismometery is a novel, cutting edge capability that involves balloon based systems for detecting and geolocating sources of infrasound. The incident infrasound from a range of sources such as volcanos, earthquakes, explosions, supersonic aircraft impinges upon the balloon system causing it to respond dynamically. The dynamic response is post-processed to locate the infrasound source. This report documents the derivation of an analytical model that predicts the balloon dynamics. Governing equations for the system are derived as well as a transfer function relating the infrasound signal to the net force on the balloon components. Experimental measurements of the infrasound signals are convolved with the transfer function and the governing equations numerically time integrated to obtain predictions of the displacement, velocity and acceleration of the balloon system. The predictions are compared to the experimental measurements with good agreement observed. The derivation focuses only on the vertical dynamics of the balloon system. Future work will develop governing equations for the swinging response of the balloon to the incident infrasound.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Improving Additive Manufactured Component Performance through Multi-Scale Microstructure Simulation and Process Optimization

The purpose of this project was to utilize computational tools to understand the relationships between processing, microstructure, and properties for additively manufactured (AM) aluminum alloys for automotive applications, and to provide an engineering solution for helping to optimize process conditions. The project leverages ORNL developments in computational modeling, including AM process modeling, phase-field based microstructure evolution predictions, and data analytics techniques for mapping process conditions to material outcomes. The project utilized an Al-Cu-Mn-Zr alloy as a model material for studying formation of defects and microstructural features in response to variations in process conditions. Based on both pre-existing experimental data and simulation results, statistical process maps were constructed to identify regions of process space with minimal defect formation and advantageous microstructures and properties. The software tools used for this purpose were successful disseminated to GM, who were able to successful compile the relevant HPC codes within their own computing ecosystem and perform initial calculations to reproduce ORNL results.

36 MATERIALS SCIENCE↗

Melt-Pool Dynamics and Microstructure of Mg Alloy WE43 under Laser Powder Bed Fusion Additive Manufacturing Conditions

Magnesium-based alloy WE43 is a state-of-the-art bioresorbable metallic implant material. There is a need for implants with both complex geometries to match the mechanical properties of bone and refined microstructure for controlled resorption. Additive manufacturing (AM) using laser powder bed fusion (LPBF) presents a viable fabrication method for implant applications, as it offers near-net-shape geometrical control, allows for geometry customization based on an individual patient, and fast cooling rates to achieve a refined microstructure. In this study, the laser–alloy interaction is investigated over a range of LPBF-relevant processing conditions to reveal melt-pool dynamics, pore formation, and the microstructure of laser-melted WE43. In situ X-ray imaging reveals distinct laser-induced vapor depression morphology regimes, with minimal pore formation at laser-scan speeds greater than 500 mm/s. Optical and electron microscopy of cross-sectioned laser tracks reveal three distinct microstructural regimes that can be controlled by adjusting laser-scan parameters: columnar, dendritic, and banded microstructures. These regimes are consistent with those predicted by the analytic solidification theory for conduction-mode welding, but not for keyhole-mode tracks. The results provide insight into the fundamental laser–material interactions of the WE43 alloy under AM-processing conditions and are critical for the successful implementation of LPBF-produced WE43 parts in biomedical applications.

36 MATERIALS SCIENCE↗

Resonant Metasurfaces with Van Der Waals Hyperbolic Nanoantennas and Extreme Light Confinement

This work reports on a metasurface based on optical nanoantennas made of van der Waals material hexagonal boron nitride. The optical nanoantenna made of hyperbolic material was shown to support strong localized resonant modes stemming from the propagating high-k waves in the hyperbolic material. An analytical approach was used to determine the mode profile and type of cuboid nanoantenna resonances. An electric quadrupolar mode was demonstrated to be associated with a resonant magnetic response of the nanoantenna, which resembles the induction of resonant magnetic modes in high-refractive-index nanoantennas. The analytical model accurately predicts the modes of cuboid nanoantennas due to the strong boundary reflections of the high-k waves, a capability that does not extend to plasmonic or high-refractive-index nanoantennas, where the imperfect reflection and leakage of the mode from the cavity complicate the analysis. In the reported metasurface, excitations of the multipolar resonant modes are accompanied by directional scattering and a decrease in the metasurface reflectance to zero, which is manifested as the resonant Kerker effect. Van der Waals nanoantennas are envisioned to support localized resonances and can become an important functional element of metasurfaces and transdimensional photonic components. By designing efficient subwavelength scatterers with high-quality-factor resonances, this work demonstrates that this type of nanoantenna made of naturally occurring hyperbolic material is a viable substitute for plasmonic and all-dielectric nanoantennas in developing ultra-compact photonic components.

Chemistry↗