Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dimensional 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 541 records · Page 30

Nonlinear mode coupling and energetics of driven magnetized shear-flow turbulence

To comprehensively understand the saturation of two-dimensional (2D) magnetized Kelvin–Helmholtz-instability-driven turbulence, energy transfer analysis is extended from the traditional interaction between scales to include eigenmode interactions, by using the nonlinear couplings of linear eigenmodes of the ideal instability. While both kinetic and magnetic energies cascade to small scales, a significant fraction of turbulent energy deposited by unstable modes in the fluctuation spectrum is shown to be re-routed to the conjugate-stable modes at the instability scale. They remove energy from the forward cascade at its inception. The remaining cascading energy flux is shown to attenuate exponentially at a small scale, dictated by the large-scale stable modes. Guided by a widely used instability-saturation assumption, a general quasi-linear model of instability is tested by retaining all nonlinear interactions except those that couple to the large-scale stable modes. These complex interactions are analytically removed from the magnetohydrodynamic equations using a novel technique. Observations are an explosive large-scale vortex separation instead of the well-known merger of 2D, a dramatic enhancement in turbulence level and spectral energy fluxes, and a reduced small-scale dissipation length scale. These show the critical role of the stable modes in instability saturation. Possible reduced-order turbulence models are proposed for fusion and astrophysical plasmas, based on eigenmode-expanded energy transfer analyses.

Physics↗

DAmodel: hierarchical Bayesian modelling of DA white dwarfs for spectrophotometric calibration

We use hierarchical Bayesian modelling to calibrate a network of 32 all-sky faint DA white dwarf (DA WD) spectrophotometric standards (⁠16.5 < V , 19.5⁠) alongside three CALSPEC standards, from 912 Å to 32 μm. The framework is the first of its kind to jointly infer photometric zero points and WD parameters (surface gravity log g⁠, effective temperature T eff ⁠, extinction A V ⁠, dust relation parameter R V ) by simultaneously modelling both photometric and spectroscopic data. We model panchromatic Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) UVIS and IR photometry, HST/STIS UV spectroscopy, and ground-based optical spectroscopy to sub-per cent precision. Photometric residuals for the sample are the lowest yet yielding < 0.004 mag RMS on average from the UV to the NIR, achieved by jointly inferring time-dependent changes in system sensitivity and WFC3/IR count-rate nonlinearity. Our GPU-accelerated implementation enables efficient sampling via Hamiltonian Monte Carlo, critical for exploring the high-dimensional posterior space. The hierarchical nature of the model enables population analysis of intrinsic WD and dust parameters. Inferred spectral energy distributions from this model will be essential for calibrating the James Webb Space Telescope as well as next-generation surveys, including Vera Rubin Observatory’s Legacy Survey of Space and Time and the Nancy Grace Roman Space Telescope.

methods: statistical↗

Effective field theory of Stückelberg vector bosons

We explore the effective field theory of a vector field $X^μ$ that has a Stückelberg mass. The absence of a gauge symmetry for $X^μ$ implies Lorentz-invariant operators are constructed directly from Xμ. Beyond the kinetic and mass terms, allowed interactions at the renormalizable level include $X_μX^μH^†H, (X_μX^μ)^2$, and $X_μj^μ$, where $j^μ$ is a global current of the SM or of a hidden sector. We show that all of these interactions lead to scattering amplitudes that grow with powers of $\sqrt{s}/m_X$, except for the case of $X_μj^μ$ where $j^μ$ is a nonanomalous global current. The latter is well-known when $\textit{X}$ is identified as a dark photon coupled to the electromagnetic current, often written equivalently as kinetic mixing between $\textit{X}$ and the photon. The power counting for the energy growth of the scattering amplitudes is facilitated by isolating the longitudinal enhancement. We examine in detail the interaction with an anomalous global vector current $X_μj^μ_{\text{anom}}$, carefully isolating the finite contribution to the fermion triangle diagram. We calculate the longitudinally-enhanced observables $Z → X_γ$ (when $m_X < m_Z), f\bar{f} → X_γ$, and $Z_γ → Z_γ$ when $\textit{X}$ couples to the baryon number current. Introducing a "fake" gauge-invariance by writing $X^μ = A^μ – ∂^μπ/m_X$, the would-be gauge anomaly associated with $A_μj^μ_{\text{anom}}$ is canceled by $j^μ_{\text{anom}}∂_μπ/m_X$; this is the four-dimensional Green-Schwarz anomaly-cancellation mechanism at work. Our analysis demonstrates there is a much larger set of possible interactions that an EFT with a Stückelberg vector field can have, revealing scattering amplitudes that grow with energy. The growth of these amplitudes can be tamed by a dark Higgs sector, but this requires dark Higgs boson interactions (and reintroduces fine-tuning in the dark Higgs sector) that can be separated from X interactions only in the limit $g \ll 1$ .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Tensor network simulations of quasi-GPDs in the massive Schwinger model

Generalized parton distribution functions (GPDs) are off-diagonal light-cone matrix elements that encode the internal structure of hadrons in terms of quark and gluon degrees of freedom. In this work, we present the first nonperturbative study of quasi-GPDs in the massive Schwinger model, quantum electrodynamics in 1+1 dimensions (QED 2 ), within the Hamiltonian formulation of lattice field theory. Quasidistributions are spatial correlation functions of boosted states, which approach the relevant light-cone distributions in the luminal limit. Using tensor networks, we prepare the first excited state in the strongly coupled regime and boost it to close to the light-cone on lattices of up to 400 lattice sites. We compute both quasiparton distribution functions and, for the first time, quasi-GPDs, and study their convergence for increasingly boosted states. In addition, we perform analytic calculations of GPDs in the two-particle Fock-space approximation and in the Reggeized limit, providing qualitative benchmarks for the tensor network results. Our analysis establishes computational benchmarks for accessing partonic observables in low-dimensional gauge theories, offering a starting point for future extensions to higher dimensions, non-Abelian theories, and quantum simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Identification of Dynamic Force Coefficients for an Additively Manufactured Hermetic Squeeze Film Bearing Support Damper Utilizing a Pass-Through Channel

Abstract The following paper presents breakthrough experimental results for a new hermetic squeeze film damper (HSFD) concept that is integrally designed within an externally pressurized tilting-pad radial gas bearing support. The flexibly damped gas bearing module was designed for a 7.2″ (183 mm) diameter shaft and fabricated using direct metal laser melting (DMLM); also known as additive manufacturing. The bearing and HSFD were sized based on ongoing studies for oil-free supercritical carbon dioxide (sCO2) power turbines in the 8.5 MW–10 MW power range. The development of the new damper concept was motivated by past dynamic testing on HSFD, which generated frequency-dependent stiffness and damping force coefficients. In efforts to eliminate the frequency dependency, a new HSFD architecture was conceived that adds accumulator volumes and a pass-through channel to previously conceived HSFD flow network designs. The other motivation for the work is the need to develop a cost-effective and reliable oil-free bearing technology that is scalable to large power turbomachinery applications. There were several objectives for the following work. The first objective was to successfully design and fabricate a single piece bearing-damper using additive manufacturing, while dimensionally controlling critical design features. The paper discusses the manufacturing steps and shows cut-ups that reveal adequate clearance control capability with internal damper clearances. The second objective was to perform experimental testing with the new HSFD design in efforts to extract stiffness and damping coefficients for excitation frequencies within 20–160 Hz and peak vibration amplitudes between 0.25 mils (6.35 microns) to 1 mil (25.4 microns). The test results for a single HSFD bearing module indicated that the design modifications to the HSFD architecture were successful in eliminating nearly all the frequency dependencies for the stiffness force coefficient. The dynamic tests yielded a stiffness coefficient that varied between 112 klb/in. (19.6 MN/m) and 96 klb/in. (16.8 MN/m). The damping force coefficient however, exhibited relatively more variation with frequency with values residing between 175 lb-s/in. (31 kN-s/m) to 214 lb-s/in. (37 kN-s/m). Finally, the paper advances a three-dimensional fluid-structure interaction (FSI) model using transient finite element analysis (FEA) coupled to a computational fluid dynamics (CFD) model. The FSI analysis performed between 20 Hz and 80 Hz was used to predict the stiffness and damping of the HSFD using a quarter-section model of the damper. The FSI analysis was able to support test results by showing only a 6–7.4% change in the magnitude of force coefficients. Stiffness predictions agree reasonably well with experiments whereas damping is underpredicted.

Engineering↗

Multi-Doppler radar analysis from CSAPR, CHIVO, COW, and RMA-1 radars during the CACTI/RELAMPAGO experiments in Argentina in 2018

This data set contains multi-Doppler radar analysis from CSAPR-2, CSU-CHIVO, COW, and RMA-1 radars. These radars were collecting dual-polarization data during the CACTI/RELAMPAGO experiments in Argentina in 2018. Doppler analysis is systematically conducted for 31 days with convection. Three-dimensional wind fields are retrieved using the PyDDA algorithm. Dual-polarization information is also included in the data set.

3D cartesian gridded corrected mean Doppler veloci↗

National direct-drive program on OMEGA and the National Ignition Facility

A major advantage of the laser direct-drive (DD) approach to ignition is the increased fraction of laser drive energy coupled to the hot spot and relaxed hot-spot requirements for the peak pressure and convergence ratios relative to the indirect-drive approach at equivalent laser energy. With the goal of a successful ignition demonstration using DD, the recently established national strategy has several elements and involves multiple national and international institutions. These elements include the experimental demonstration on OMEGA cryogenic implosions of hot-spot conditions relevant for ignition at MJ-scale energies available at the National Ignition Facility (NIF) and developing an understanding of laser-plasma interactions and laser coupling using DD experiments on the NIF. DD designs require reaching central stagnation pressures in excess of 100 Gbar. The current experiments on OMEGA have achieved inferred peak pressures of 56 Gbar (Regan et al 2016 Phys. Rev. Lett. 117 025001). Extensive analysis of the cryogenic target experiments and two- and three-dimensional simulations suggest that power balance, target offset, and target quality are the main limiting factors in target performance. In addition, cross-beam energy transfer (CBET) has been identified as the main mechanism reducing laser coupling. Reaching the goal of demonstrating hydrodynamic equivalence on OMEGA includes improving laser power balance, target position, and target quality at shot time. CBET must also be significantly reduced and several strategies have been identified to address this issue. (2016-01-01)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Imaging Anisotropic Waveguide Exciton Polaritons in Tin Sulfide

In recent years, novel materials supporting in-plane anisotropic polaritons have attracted a great deal of research interest due to their capability of shaping nanoscale field distributions and controlling nanophotonic energy flows. Here we report a nano-optical imaging study of waveguide exciton polaritons (EPs) in tin sulfide (SnS) in the near-infrared (near-IR) region using scattering-type scanning near-field optical microscopy (s-SNOM). With s-SNOM, we mapped in real space the propagative EPs in SnS, which show sensitive dependence on the excitation energy and sample thickness. Moreover, we found that both the polariton wavelength and propagation length are anisotropic in the sample plane. In particular, in a narrow spectral range from 1.32 to 1.44 eV, the EPs demonstrate quasi-one-dimensional propagation, which is rarely seen in natural polaritonic materials. Here, a further analysis indicates that the observed polariton anisotropy originates from the different optical band gaps and exciton binding energies along the two principal crystal axes of SnS.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Three-flavor, full momentum space neutrino spin oscillations in neutron star mergers

In the presence of anisotropic neutrino and antineutrino fluxes, the quantum kinetic equations drive coherent oscillations in neutrino helicity, frequently referred to as spin oscillations. These oscillations depend directly on the absolute mass scale and Majorana phase, but are usually too transient to produce important effects. In this paper we present a full momentum-space analysis of Majorana neutrino spin oscillations in a snapshot of a three-dimensional neutron star merger simulation. We find an interesting angular dependence that allows for that resonant and adiabatic oscillations to occur along specific directions in a large volume of the merger remnant. The solid angle spanned by these directions is extremely narrow in general. We then analyze spin transformation in the presence of flavor transformation by characterizing how the effect's resonance and timescale change during a fast flavor instability. For this analysis, we derive a generalized resonance condition that poses a restrictive requirement for resonance to exist in any flavor channel. Finally, we determine that spin oscillations at all locations in the merger snapshot have a length scale that is too large for significant oscillations to be expected even where there exist resonant and adiabatic directions.

79 ASTRONOMY AND ASTROPHYSICS↗

Clustering Analysis of Commercial Vehicles Using Automatically Extracted Features from Time Series Data

Standard of practice approaches to time series cluster analysis involve careful feature engineering, often utilizing expert input to tune and select features by hand. In many cases, expert input may not be readily available, or there may not yet exist a community consensus on the ideal features for a given application. This paper compares the results of several cluster analysis methods, using both hand selected features and those extracted automatically, when applied to large geospatial time series telematics data from commercial trucking fleets. The impacts of feature selection, dimensionality reduction, and choice of clustering algorithm on the quality of clustering results are explored. Results from this analysis confirm prior results that domain agnostic features are competitive with the hand engineered features with respect to clustering quality metrics. These results also provide new insight into the most successful strategies for identifying structure in large unstructured vehicle telematics data, and suggest that time series clustering using automatic feature extraction can be an effective approach to extract structure from large scale geospatial time series data in cases when hand selected features are not available.

33 ADVANCED PROPULSION SYSTEMS↗

An Analysis of System Balance and Architectural Trends Based on Top500 Supercomputers

Supercomputer design is a complex, multi-dimensional optimization process, wherein several subsystems need to be reconciled to meet a desired figure of merit performance for a portfolio of applications and a budget constraint. However, overall, the HPC community has been gravitating towards ever more FLOPS, at the expense of many other subsystems. To draw attention to overall system balance, in this paper, we analyze balance ratios and architectural trends in the world's most powerful supercomputers. Specifically, we have collected performance characteristics of systems between1993 and 2018 based on the Top500 lists, and then analyzed their architectures from diverse system design perspectives. Notably, our analysis studies the performance balance of the machines, across a variety of subsystems such as compute, memory, I/O, interconnect, intra-node connectivity and power. Our analysis reveals that balance ratios of the various subsystems need to be considered carefully alongside the application workload portfolio to provision the subsystem capacity and bandwidth specifications, which can help achieve optimal performance.

97 MATHEMATICS AND COMPUTING↗

An Analysis of System Balance and Architectural Trends Based on Top500 Supercomputers

Supercomputer design is a complex, multi-dimensional optimization process, wherein several subsystems need to be reconciled to meet a desired figure of merit performance for a portfolio of applications and a budget constraint. However, overall, the HPC community has been gravitating towards ever more Flops, at the expense of many other subsystems. To draw attention to overall system balance, in this paper, we analyze balance ratios and architectural trends in the world’s most powerful supercomputers. Specifically, we have collected the performance characteristics of systems between 1993 and 2019 based on the Top500 lists and then analyzed their architectures from diverse system design perspectives. Notably, our analysis studies the performance balance of the machines, across a variety of subsystems such as compute, memory, I/O, interconnect, intra-node connectivity and power. Our analysis reveals that balance ratios of the various subsystems need to be considered carefully alongside the application workload portfolio to provision the subsystem capacity and bandwidth specifications, which can help achieve optimal performance.

Khan, Awais↗

Scalability Analysis of Quantum Models for Stress and Emotion Detection

Stress and emotion detection from high-dimensional physiological signals is a challenging task, particularly when aiming for accurate classification across diverse behavioral states. Quantum machine learning (QML) is promising for modeling such high-dimensional data, but scalability is limited by qubit resources and the exponential cost of classical statevector simulation. This work studies the scalability of quantum support vector machines (QSVMs) for binary stress detection and three-class emotion recognition (Negative/Neutral/Positive) under varying qubit counts and angle-encoding strategies. We also present a comparison study with one-feature-per-qubit (1:1) and two-features-per-qubit (2:1) mappings. Experiments are executed on HPC infrastructure using NVIDIA CUDA-Q to evaluate performance, variance, and class-dependent separability at higher-qubit setups. Results show that larger Hilbert spaces can improve peak accuracy but may increase instability. At the same time, dense 2:1 encoding yields more consistent stress detection performance. For emotion recognition, scaling improves discrimination for classes like Negative and Positive more than Neutral. We find that effective QML scaling is task-dependent and benefits more from encoding design than simply increasing qubit count.

Onim, Md. Saif Hassan [University of Tennessee, Kn↗

The virtual element method for linear elastodynamics models: Design, analysis, and implementation

We design the conforming virtual element method for the numerical simulation of two dimensional time-dependent elastodynamics problems. We investigate the performance of the method both theoretically and numerically. We prove the stability and the convergence of the semi-discrete approximation in the energy norm and derive optimal error estimates. We also show the convergence in the L 2 norm. The performance of the virtual element method is assessed on a set of different computational meshes, including non-convex cells up to order four in the h-refinement setting. Exponential convergence is also experimentally seen in the p-refinement setting.

97 MATHEMATICS AND COMPUTING↗

Involute Working Group – FSI Analysis of Fuel Plates Using Finite Volume and Finite Element Methods

The three involute plate research reactors RHF, HFIR, and FRM II have expressed an interest in using computational software to carry their steady-state safety analysis. Since these tools represent a significant departure from the methods used currently (one-dimensional), the acceptability of the new approach by regulators requires thorough verification and validation of these tools. Therefore, Argonne National Laboratory and the three involute-plate reactors formed an informal group called the Involute Working Group aiming at qualifying computational tools to perform steady-state safety analysis. The present report focuses on a comparison of finite volume and finite element methods to model solids in fluid-structure interaction problems with the goal to estimate the coolant flow-induced fuel plate deflections obtained with the two methods. The finite volume method will be obsoleted in STARCCM+ by the end of 2021, nevertheless, this evaluation is important because the method was used by ANL researchers to model the response of the fuel plates, despite its drawbacks, which are discussed in the report. It was essential to check how those estimates compare to the results obtained with the finite element method that is considered superior for structural analysis. Various geometries, i.e., flat, cylindrical and circle-involute fuel plates, as well as coolant flow speed, were considered. The comparison shows that, independently of the plate geometry, the finite volume method significantly underestimates the deflection as compared to finite element method for coarser meshes. When the discretization is developed as a result of a mesh sensitivity study using finite element method, the result obtained using finite volume method can be a few times smaller than the corresponding finite element method solution. A code-to-code comparison , between STAR-CCM+ and LS-DYNA was included in the analysis. Within the LS-DYNA models, two types of finite element formulations were used: solid and shell finite elements. Mesh sensitivity study showed that both approaches converge to a similar value that was obtained with STAR-CCM+ finite element solver. The evaluation of the computational solvers was extended by adding two benchmark cases from the STAR-CCM+ Verification Suite and presented in the Appendix A. The selected cases are: (1) bending of a cantilever beam under external load, and (2) cylindrical shell deformation analysis, known in the literature as ‘Scordelis-Lo roof’. The problems were solved with finite volume, and finite element methods, and the results confirmed the previously discussed findings. The analysis shows that the finite element solver is superior to the finite volume solver in terms of representation of model geometry and estimating the structural behavior of fuel plates. Depending on the ratio of the load to the flexibility of the plate, the finite volume solver can greatly under- or overestimate the structural response if a very carefully selected mesh is not used.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A general framework for quantifying uncertainty at scale

Abstract In many fields of science, comprehensive and realistic computational models are available nowadays. Often, the respective numerical calculations call for the use of powerful supercomputers, and therefore only a limited number of cases can be investigated explicitly. This prevents straightforward approaches to important tasks like uncertainty quantification and sensitivity analysis. This challenge can be overcome via our recently developed sensitivity-driven dimension-adaptive sparse grid interpolation strategy. The method exploits, via adaptivity, the structure of the underlying model (such as lower intrinsic dimensionality and anisotropic coupling of the uncertain inputs) to enable efficient and accurate uncertainty quantification and sensitivity analysis at scale. Here, we demonstrate the efficiency of this adaptive approach in the context of fusion research, in a realistic, computationally expensive scenario of turbulent transport in a magnetic confinement tokamak device with eight uncertain parameters, reducing the effort by at least two orders of magnitude. In addition, we show that this refinement method intrinsically provides an accurate surrogate model that is nine orders of magnitude cheaper than the high-fidelity model.

Farcaş, Ionuţ-Gabriel (ORCID:0000000238841580)↗

A two-stage optical fusion framework for wildfire severity mapping across the conterminous United States

Accurate wildfire severity mapping (WSM) is essential for post-fire recovery planning, erosion risk assessment, ecosystem monitoring, and disaster risk reduction. Although Landsat and Sentinel optical imagery have been widely used for burn severity assessment, the added value of fusing multiple optical sensors has not been sufficiently quantified across diverse fire events, particularly since the launch of Landsat-9. This study evaluates whether multisensor optical fusion improves wildfire severity mapping relative to single-sensor baselines using Sentinel-2, Landsat-8, and Landsat-9 imagery across 40 wildfire events in the conterminous United States. We tested a two-stage fusion framework that combines feature-level fusion with pixel-level dimensionality reduction. First, feature-level fused datasets were created through early fusion by combining standardized post-fire bands from each sensor into a single predictor stack. Both raw reflectance bands and pairwise spectral transforms were retained to capture within- and cross-sensor spectral interactions. Second, Linear Discriminant Analysis was applied to both single-sensor and fused datasets to produce comparable low-dimensional feature spaces. Six machine-learning classifiers were then used to benchmark model performance with repeated spatially buffered train–test splits. Results show that Landsat-9 was the strongest single-sensor baseline. Among the fusion strategies, Sentinel-2 + Landsat-9 produced the most consistent improvement and reduced performance variability. Landscape-condition analysis further showed that this fusion was most beneficial in shrubland-dominated and high-terrain fires, where it achieved the highest overall mean accuracy and the fewest failures. In contrast, its benefits were less reliable in evergreen forests, mixed vegetation, and low- to moderate-elevation terrain. In operational settings, the Sentinel-2 + Landsat-9 configuration offers a practical solution for post-fire recovery planning, erosion-risk assessment, watershed management, and ecological monitoring when field observations are available and timely satellite-based information is needed.

Landsat↗