Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical solution”

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 397 records · Page 22

Proton deflectometry analysis in magnetized plasmas: Magnetic field reconstruction in one dimension

Proton deflectometry is used in magnetized high-energy-density plasmas to observe electromagnetic fields. Here we describe a reconstruction algorithm to recover the electromagnetic fields from proton fluence data in 1-D. The algorithm is verified against analytic solutions and applied to example data. Next, we study the role of source fluence uncertainty for 1-D reconstructions. We show that reconstruction boundary conditions can be used to constrain the source fluence profile and use this to develop a reconstruction using a specified pair of boundary conditions on the magnetic field. From these considerations, we experimentally demonstrate a hybrid mesh-fluence reconstruction technique where fields are reconstructed from fluence data in an interior region with boundary conditions supplied by direct mesh measurements at the boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Uncovering Non-Fermi-Liquid Behavior in Hund Metals: Conformal Field Theory Analysis of an SU(2)×SU(3) Spin-Orbital Kondo Model

Hund metals have attracted attention in recent years due to their unconventional superconductivity, which supposedly originates from non-Fermi-liquid (NFL) properties of the normal state. When studying Hund metals using dynamical mean-field theory, one arrives at a self-consistent “Hund impurity problem” involving a multiorbital quantum impurity with nonzero Hund coupling interacting with a metallic bath. If its spin and orbital degrees of freedom are screened at different energy scales, T sp < T orb , the intermediate energy window is governed by a novel NFL fixed point, whose nature had not yet been clarified. We resolve this problem by providing an analytical solution of a paradigmatic example of a Hund impurity problem, involving two spin and three orbital degrees of freedom. To this end, we combine a state-of-the-art implementation of the numerical renormalization group, capable of exploiting non-Abelian symmetries, with a generalization of Affleck and Ludwig’s conformal field theory (CFT) approach for multichannel Kondo models. We characterize the NFL fixed point of Hund metals in detail for a Kondo model with an impurity forming an SU ( 2 ) × SU ( 3 ) spin-orbital multiplet, tuned such that the NFL energy window is very wide. The impurity’s spin and orbital susceptibilities then exhibit striking power-law behavior, which we explain using CFT arguments. We find excellent agreement between CFT predictions and numerical renormalization group results. Our main physical conclusion is that the regime of spin-orbital separation, where orbital degrees of freedom have been screened but spin degrees of freedom have not, features anomalously strong local spin fluctuations: the impurity susceptibility increases as χ sp imp ~ ω – γ , with γ > 1 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Sachdev-Ye-Kitaev models and beyond: Window into non-Fermi liquids

In this work is a review of the Sachdev-Ye-Kitaev (SYK) model of compressible quantum many-body systems without quasiparticle excitations, and its connections to various theoretical studies of non-Fermi liquids in condensed matter physics. The review is placed in the context of numerous experimental observations on correlated electron materials. Strong correlations in metals are often associated with their proximity to a Mott transition to an insulator created by the local Coulomb repulsion between the electrons. The phase diagrams of a number of models of such a local electronic correlation are explored, employing a dynamical mean-field theory in the presence of random spin exchange interactions. Numerical analyses and analytical solutions, using renormalization group methods and expansions in large spin degeneracy, lead to critical regions that display SYK physics. The models studied include the single-band Hubbard model, the t–J model, and the two-band Kondo-Heisenberg model in the presence of random spin exchange interactions. Also examined are non-Fermi liquids obtained by considering each SYK model with random four-fermion interactions to be a multiorbital atom, with the SYK atoms arranged in an infinite lattice. Connections are made to theories of sharp Fermi surfaces without any low-energy quasiparticles in the absence of spatial disorder, obtained by coupling a Fermi liquid to a gapless boson; a systematic large-N theory of such a critical Fermi surface, with SYK characteristics, is obtained by averaging over an ensemble of theories with random boson-fermion couplings. Finally, an overview of the links between the SYK model and quantum gravity is presented, and the review ends with an outlook on open questions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Robust qubit interactions mediated by photonic topological edge states

Here, we investigate the coupling of two spatially separated qubits via topologically protected edge states in a two-dimensional Hofstadter lattice. In this hybrid platform, the qubits are coupled to distinct edge sites of the lattice, enabling long-range interactions mediated by topological edge modes. We solve the full system Hamiltonian and analyze the resulting eigenstate structure to uncover the conditions under which coherent qubit interactions emerge. Our analysis reveals that the effective coupling is highly sensitive to the qubit placement, energy detuning, and the topological character of the edge spectrum. We obtain an analytical solution that goes beyond the perturbative regime, capturing the full interplay between the qubits and edge modes. These results provide a foundation for exploring information transport and many-body effects in engineered quantum systems where interactions are mediated by topological edge modes.

Chern insulators↗

Magneto-opto-phononic inverse Faraday effect

Nonlinear frequency conversion processes, such as optical rectification, difference-frequency generation, and sum-frequency generation, are fundamental for producing electromagnetic radiation at diverse frequencies. In this work, we demonstrate that coherently excited infrared-active phonons can act as transducers for generating nonlinear magnetizations through phonon-magnon interactions, analogous to nonlinear optical frequency conversion. We derive analytical solutions for the time-dependent magnetizations arising from the second-order response to the electric field component of an ultrashort laser pulse. These solutions enable us to define second-order nonlinear magneto-electric susceptibilities, which describe rectification, impulsive excitation, and sum-frequency excitation of coherent magnons. Our theoretical framework naturally incorporates the conventional magneto-optic and phonon inverse Faraday effects and predicts a hybrid magneto-opto-phononic inverse Faraday effect involving photon-phonon-magnon scattering. This work highlights nonlinear phononics as a pathway for controlling magnetization in solids.

Landau-Lifschitz-Gilbert equation↗

Additive Manufactured Compact Microwave Absorbers

A high-performance, compact microwave absorber was created using Fused Deposition Modeling (FDM) 3D printing. Both a narrowband and a broadband absorber were created. The narrowband absorber was designed at 4.9 GHz, mid-band in WR-187 waveguide. The broadband absorber tried to achieve the best attenuation across the entire 3.95 to 5.85 GHz band. Two types of carbon loaded polylactic acid (PLA) plastic and one type of unloaded PLA were 3D printed with variable percentages of air to achieve different values of effective dielectric constant and loss tangent. The absorber comprised five or six rectangular pieces of these plastic materials. The thickness and fill factor values for each piece were optimized to minimize reflection through fast analytic modeling in MATLAB®. The results were then verified by HFSS® simulation as well. The stack progressed from the lowest loss and lowest dielectric constant to the highest at the shorting end. The final narrowband load had simulated return loss of 87 dB at 4.9 GHz with an analytic solution in MATLAB. The measured return loss of the 3D printed attenuator was 73 dB at 4.929 GHz. The total length of the absorber was 2.44 inches. A commercial absorber for WR-187 with return loss of 40 dB has length of 13 inches. The experiment proves that an effective and compact microwave absorber can be created using 3D printing.

36 MATERIALS SCIENCE↗

Bayesian Physics Informed Spatio-Temporal Network for Streamflow Data Imputation

Reliable reconstruction of incomplete streamflow records is critical for improving hydrological forecasting, flood preparedness, and water resource management. However, large observational gaps and uncertainties in governing physical parameters limit the accuracy of traditional statistical and machinelearning imputation frameworks. To address these challenges, we develop a Bayesian Physics-Informed Spatio-Temporal Network (BPI-STNet) that jointly captures spatial and temporal dependencies while enforcing hydrologic consistency through embedded physical constraints. The framework integrates a GraphSAGE-LSTM architecture to model spatial connectivity across gauges and temporal flow dynamics, coupled with a Bayesian update mechanism to estimate uncertain parameters in a simplified water-balance framework. Unlike conventional physics-informed networks that rely on sampling-based posterior estimation, BPI-STNet derives an analytic solution to the inverse problem, allowing closed-form Bayesian updates of uncertain parameters Λ={α,β,k} using Gaussian priors and likelihoods. Applied to daily observations from the Susquehanna River Basin (1980-2022), BPI-STNet achieves substantial improvements over a purely data-driven RGNN baseline, which reduced RMSE by 23 % and MAE by 9 %, and achieving an average NSE values up to 0.96. The results demonstrate that coupling Bayesian inference with physics-informed learning yields physically consistent, uncertainty-aware reconstructions that preserve the temporal persistence and statistical distribution of observed flows. The proposed framework establishes a generalizable paradigm for data-sparse hydrologic systems where both data fidelity and physical interpretability are essential.

Krishnan Kutty Ambika, Anukesh [ORNL] (ORCID:00000↗

Taylor-Expansion-Based Robust Power Flow in Unbalanced Distribution Systems: A Hybrid Data-Aided Method

Traditional power flow methods often adopt certain assumptions designed for passive balanced distribution systems, thus lacking practicality for unbalanced operation. moreover, their computation accuracy and efficiency are heavily subject to unknown errors and bad data in measurements or prediction data of distributed energy resources (ders). to address these issues, this paper proposes a hybrid data-aided robust power flow algorithm in unbalanced distribution systems, which combines taylor series expansion knowledge with a data-driven regression technique. the proposed method initiates a linearization power flow model to derive an explicitly analytical solution by modified taylor expansion. to mitigate the approximation loss that surges due to the der integration and bad data, we further develop a data-aided robust support vector regression approach to estimate the errors efficiently. comparative analysis in the 13-bus and 123-bus ieee unbalanced feeders shows that the proposed hybrid algorithm achieves superior computational efficiency, with guaranteed accuracy and robustness against outliers.

data-driven↗

Small-sized specimen design with the provision for high-frequency bending-fatigue testing

The testing and study of emerging materials—such as additively manufactured materials—demands for specimen designs that are cost effective and time saving. Here the design of a small-sized bending-fatigue test specimen for an ultrasonic fatigue testing system is reported in this paper. The design is optimized based on the finite element analysis and analytical-solution results to achieve the proper vibration shape and stress distribution. The proposed design is evaluated in the high- and very-high-cycle fatigue regimes under 20-kHz frequency. Both simulation and testing results confirm that the desirable vibration mode occurs and the specimen fails at the designated test (gauge) section, where the maximum stress exists. The stress–life (S–N) curve is obtained for Inconel alloy 718 and indicates an expected trend.

36 MATERIALS SCIENCE↗

Code Verification of Multiple Physics-Fidelity Models in Hypersonic Aerodynamics

Hypersonic aerodynamics models exist across a range of physics fidelities with associated computational expenses. These models may be run independently or in a multifidelity framework that leverages their complementary strengths of speed for lower-fidelity and accuracy for higher-fidelity models. This work presents applied code verification of two lower-fidelity models contained within the Sandia hypersonic aerodynamics code. Each model has a different form that requires individualized verification approaches, including comparison to analytical solutions as well as manufactured solutions with order-of-accuracy testing. In conclusion, results of this effort include the identification and resolution of code errors and shortcomings, as well as the demonstration of code correctness and consistency for both models.

Aerodynamics↗

Implementation and Verification of Cable Bending Stiffness in MoorDyn

The relatively large motions experienced by floating wind turbines and wave energy converters pose a challenge for power cables, whose internal components provide significant bending resistance and are sensitive to deformation. The behavior and associated design considerations of power cables in these highly dynamic applications make coupled analysis relevant for design. Bending stiffness capabilities have recently been added to the lumped-mass mooring dynamics model MoorDyn to enable simulation of dynamic power cables. MoorDyn is a common modeling choice for floating wind energy simulation (often coupled with OpenFAST) and floating wave energy converter simulation (often coupled with WEC-Sim) but the model’s previous line elasticity formulation only considered axial stiffness. To properly capture the dynamics of power cables, a bending stiffness model has been added that approximates cable curvature based on the difference in tangent vectors of adjacent elements. The resulting bending moment is realized by applying forces on adjacent nodes, enabling cable modeling while leaving the underlying lumped-mass formulation unchanged. In this paper, the new bending stiffness implementation is verified in static conditions against analytical solutions and then in a dynamic power cable scenario in comparison with the commercial simulator OrcaFlex. The dynamic scenario uses prescribed motions and includes wave loadings on the cable. Results indicate correct implementation of bending stiffness and show close agreement with OrcaFlex.

50 EE - Wind and Water Power Program - Water (EE-4↗

Practical considerations in using a binary pseudorandom array for instrument transfer function calibration

Binary pseudo-random array (BPRA) artifacts are useful devices for calibrating the instrument transfer function (ITF) of interferometric microscopes and other optical and non-optical surface and wavefront measurement instruments. The intrinsic white noise character of the power spectral density function of the artifact simplifies the deconvolution of the ITF from the measured power spectral density (PSD). However, resampling of the BPRA intrinsic artifact features with the measurement tool's specific sampling pattern modifies the white noise character of the intrinsic spectrum and needs to be accounted for in the ITF-based data deconvolution process. We have developed an analytic solution to the spectrum of a resampled one- and two- dimensional BPRA. The resultant nominal PSD function is a simple twoparameter cosine function with a period equal to the resampled pixel width. A transfer function model for interferometric microscopes that incorporates this function, along with an ITF that includes aliasing effects and variable numerical aperture (NA), wavelength, and obscuration factor, is used to fit to the BPRA PSDs measured by an interference microscope for a range of objective and zoom lens magnification combinations.

Takacs, PZ↗

The weakly nonlinear response and nonaffine interpretation of the Johnson–Segalman/Gordon–Schowalter model

We derive new analytical solutions for the non-affine Johnson-Segalman/Gordon-Schowalter (JS/GS) constitutive equation with a general relaxation kernel in medium-amplitude oscillatory shear (MAOS) deformation. The results show time-strain separable (TSS) nonlinearity, therefore providing new physically-meaningful interpretation to the heuristic TSS nonlinear parameter in MAOS. The upper-convected, lower-convected, and corotational Maxwell models are all subsets of the results presented here. The model assumes that the microscale elements causing stress in the material slip compared to the continuum deformation. We introduce a visualization of the non-affine deformation field that acts on stress-generating elements to reinforce the hysical interpretation of the JS/GS class of models. Lastly, a case study is presented where previously published results, from fitting TSS models to MAOS data, can be re-interpreted based on the concept of non-affine motion of the JS/GS framework.

36 MATERIALS SCIENCE↗

Dating submarine landslides using the transient response of gas hydrate stability

Abstract Submarine landslides are prevalent on the modern-day seafloor, yet an elusive problem is constraining the timing of past slope failure. We present a novel age-dating technique based on perturbations to underlying gas hydrate stability caused by slide-impacted seafloor changes. Using three-dimensional (3-D) seismic data, we mapped an irregular bottom simulating reflection (BSR) underneath a submarine landslide in the Orca Basin, Gulf of Mexico. The irregular BSR mimics the pre-slide seafloor geometry rather than the modern bathymetry. Therefore, we suggest that the gas hydrate stability zone (GHSZ) is still adjusting to the post-slide sediment temperature. We applied transient conductive heat-flow modeling to constrain the response of the GHSZ to the slope failure, which yielded a most likely age of ca. 8 ka, demonstrating that gas hydrate can respond to landslides even on multimillennial time scales. We further provide a generalized analytical solution that can be used to remotely date submarine slides in the absence of traditional dating techniques.

Geology↗

Quantifying Uncertainty in HPC Job Queue Time Predictions

High Performance Computing (HPC) has developed at an unprecedented pace in recent decades. This growth has demanded corresponding development in the area of HPC Operational Data Analytics (ODA), which encompasses a wide range of data analysis techniques, ML/AI efforts, tools, and visualizations. Published studies in ODA offer a variety of practical ways to inform HPC users, administrators, procurement managers, and other stakeholders. Uncertainty analysis, however, is rare in the related published literature. For instance, we identify only 1 out of 14 existing studies focused on job queue time prediction that investigates the uncertainty aspect of their proposed predictions. We recognize the utmost importance uncertainty quantification can have in such predictive analytics solutions, with consequences in how users interpret information they receive, and attempt to bridge this gap. With the goal of improving access to such insights, we develop a process for determining upper and lower bounds of the predicted queue times of a regression model at a specified confidence level. Our current research is focused on the uncertainty in predicting job queue times, yet our approach may be employed in predicting other metrics.

HPC↗

Phase Field Modeling of Crystallographic Corrosion Pits

The modeling of corrosion to understand and predict corrosion behavior is a topical issue. Here, a 3D phase field model is developed to simulate the pit morphology, primarily focusing on crystallographic pits. A crystallographic function is employed to incorporate different corrosion rates for different crystallographic planes. The model is benchmarked and validated against an analytical solution for a simple case. 3D crystallographic pits in a single crystal are simulated and the effect of substrate orientation on the pit morphology is studied. The crystallographic pit morphology changes significantly with the substrate orientation and these morphologies have a symmetry consistent with the substrate orientation. This first 3D phase field model of crystallographic pits will help in predicting the intricate shapes of pits, thereby, pushing the frontiers of pitting corrosion modeling.

36 MATERIALS SCIENCE↗

Parametric dependencies of resonant layer responses across linear, two-fluid, drift-MHD regimes

Non-axisymmetric magnetic fields arising in a tokamak either by external or internal perturbations can induce complex non-ideal MHD responses in their resonant surfaces while remaining ideally evolved elsewhere. This layer response can be characterized in a linear regime by a single parameter called the inner-layer Delta, which enables outer-layer matching and the prediction of torque balance to non-linear island regimes. Here, we follow strictly one of the most comprehensive analytic treatments including two-fluid and drift MHD effects and keep the fidelity of the formulation by incorporating the numerical method based on the Riccati transformation when quantifying the inner-layer Delta. The proposed scheme reproduces not only the predicted responses in essentially all asymptotic regimes but also with continuous transitions as well as improved accuracies. In particular, the Delta variations across the inertial regimes with viscous or semi-collisional effects have been further resolved, in comparison with additional analytic solutions. The results imply greater shielding of the electromagnetic torque at the layer than what would be expected by earlier work when the viscous or semi-collisional effects can compete against the inertial effects, and also due to the intermediate regulation by kinetic Alfven wave resonances as rotation slows down. These are important features that can alter the nonaxisymmetric plasma responses including the field penetration by external fields or island seeding process in rotating tokamak plasmas.

MHD instability↗

RAPID

Parallel computer code for the simulator for dynamics of power systems which has the capability to initiate the system and create different faults for the dynamic analysis. The code is based on time-parallel method (Parareal) with Adaptive Method Reduction (AMR). The coarse solvers for the Parareal algorithm include several Semi Analytical Solution methods. Also, Integrated simulation of coupled transmission and distribution systems can be studied.

Simunovic, Srdjan [Oak Ridge National Lab. (ORNL),↗