Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TRAJECTORIES”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Computing the solubility of argon and xenon in molten sodium chloride and potassium chloride salts

Molten salt reactors (MSRs) offer significant advancements in nuclear reactor safety and efficiency by operating at higher temperatures and lower pressures compared to traditional reactors. A critical aspect of MSR operation involves understanding the solubility of fission byproducts, particularly noble gases, in the molten salts used. This study employs molecular dynamics (MD) simulations to compute Henry’s law constants and enthalpies of solvation for argon and xenon in molten sodium chloride (NaCl) and potassium chloride (KCl). We developed a new pairwise potential for the noble gas and salt interactions based on first principles calculations. We then used this potential to calculate Henry’s law constants of the two gases in the molten salts, which were modeled using both a rigid ion model (RIM) and a polarizable ion model (PIM). The solubility calculations, performed using the Widom insertion method, show qualitative agreement with limited experimental data, highlighting the temperature dependence and greater solubility of both gases in KCl compared to NaCl. Additionally, free volume analysis elucidated the role of available space within the molten salts in governing solubility trends. Our findings suggest that PIM trajectories provide more reliable predictions for noble gas solubility than RIM due to their accurate density representation. Furthermore, these results enhance understanding of gas solubility in MSR environments, and the methods can be readily extended to other systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Range-wide population assessments for subalpine fir indicate widespread disturbance-driven decline

Subalpine forests in western North America are threatened by rapid climate change, increased activity by endemic and exotic insects and diseases, and changing wildfire regimes. The interactive effects of these stressors have resulted in pronounced population declines in many subalpine tree species; however, a systematic assessment of the status and trends of subalpine forests is lacking. Subalpine fir (Abies lasiocarpa) is a widespread species across the western United States, with documented population declines in many parts of its distribution. Here we use subalpine fir as an initial leverage point to build a more complete understanding of subalpine forest baseline conditions and responses to environmental change. Specifically, we leverage the USDA Forest Service Forest Inventory and Analysis (FIA) database to (1) ask how subalpine fir populations are changing across the species’ distribution in the western US, (2) assess the drivers of recent subalpine fir population trends, and (3) explore whether those changes imply generalized species-wide and/or system-wide decline. We found that subalpine fir abundance and basal area are declining concurrently across ~ 62% of the species’ distribution, and increasing across ~ 19%. Range-wide, we estimated 25.02 ± 2.74 % subalpine fir mortality between 2000 and 2009 and 2010–2019 FIA inventory periods, with higher mortality concentrated in the eastern Oregon Cascades, central Idaho, and parts of southern Colorado. High regeneration density did not predict positive population trajectories, which were instead associated with higher rates of adult recruitment. While the importance of different mortality agents varied substantially between ecoregions, 83.4% of total range-wide mortality was related to fire or biological disturbance. Declining subalpine fir basal area coincided with declines in the basal area of other co-occurring tree species in 39% of subalpine forest area, and with increases in conspecific basal area in 22% of forest area. Fire disturbance was the single largest cause of subalpine fir mortality; however, even where subalpine fir fire mortality was high, mortality among other species was primarily caused by insects. In conclusion, our results suggest that subalpine fir declines across large portions of the western United States are driven by forest disturbance, and that declines in subalpine fir populations may be indicative of negative change in subalpine forest systems broadly.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Assessing the fidelity of shallow-water carbonates as records of the Ni isotope composition of surface seawater

Nickel is a bioessential metal that is used in enzymes important to the C, N, and O cycles, and changes in its marine abundance and bioavailability may have affected the evolutionary trajectory of early life. Changes over time in the Ni isotope composition (δ 60 Ni) of surface seawater, which reflects biological demand for Ni, could allow for the reconstruction of the dynamics of Ni demand over Earth’s history, but this approach would require geologic records of surface seawater. Here, we investigate the fidelity of shallow-water carbonates as a record of the Ni isotope composition of surface seawater by determining how Ni is first partitioned into natural carbonates and then how post-depositional processes influence the Ni signal. Our samples come from the Great Bahama Bank, which is a well-studied, modern carbonate platform often used to study ancient platforms. We found that Ni is fractionated from seawater upon incorporation into carbonates capturing shallow (<18 cm), recent deposition (0.1 ‰–0.4 ‰ lighter than seawater). Variation among these [Ni] and δ 60 Ni values may be controlled by variation in mineral proportions. Meteoric diagenesis shifts δ 60 Ni to lower values, which we attribute to isotopically light meteoric fluids. In contrast, carbonates that experienced sediment-buffered marine diagenesis with respect to Ca isotopes and Sr/Ca ratios do not appear to differ in δ 60 Ni values from sediments generally representative of their initial deposition. The sensitivity of δ 60 Ni to diagenetic reset in these samples appears comparable to the sensitivities of Ca isotopes and Sr/Ca ratios, to first order. Thus, in general, carbonates that experienced sediment-buffered marine diagenesis with respect to these elements may hold the most promise as a record of the δ 60 Ni of coeval surface seawater. Additionally, we use our results to infer that the fraction of Ni removed from seawater into carbonates is less than 10 % of the total Ni output from the global oceans and incorporation of this Ni sink into global biogeochemical models will only have a minor impact on the modeled modern Ni budget.

Great Bahama Bank carbonates↗
DOE OSTI2024

Interpreting the spatial distribution of soil properties with a physically-based distributed hydrological model

Digital soil maps are commonly data-driven as the development of physically-based models for soil mapping is difficult due to the complexity of soils. However, physically-based hydrologic models have been successful in simulating water dynamics. Since water movement is a major driver of pedogenesis, the physical rules that govern water movement might help explain and predict the spatial variation of soil properties. Here, we demonstrate the novel use of a physically-based, distributed hydrologic model to inform the spatial distribution of soil properties. The Distributed Hydrology Soil Vegetation Model (DHSVM) was utilized to simulate soil moisture content (SM) and water table depth (WTD) in two hillslope catchments under pasture and forest management wherein hydrologic model outputs were then compared with soil properties measured in situ. SM sensors and wells were installed in both catchments to validate simulations of soil water movement via Nash-Sutcliffe Efficiency (E). In-situ observations were made at 87 sites within both catchments to study the connection between simulated water movement (SM and WTD) and observed soil properties, namely the depth and thickness of the argillic (Bt), fragic (Btx), and C horizons, and the depth of redoximorphic features. The simulated time series of SM and WTD were also clustered per season using Dynamic Time Warping (DTW), which identified similarity among time series at varying timescales. Model validation suggested that simulations of surficial SM (0–20 cm) were reasonable (E = 0.45), however, simulated subsurface SM (45–60 cm) and WTD were not sufficiently accurate. The thickness of Btx horizons were spatially grouped into different populations by SM clusters from every season except spring. For the other properties, only SM dynamics of specific seasons grouped into significantly different populations, suggesting that the explanatory power of simulated water movement varies seasonally and was greater during winter. Here, we show clusters of simulated SM separated soil properties into statistically different populations, showing that hydrologic models could inform areas that followed different water dynamics related to pedogenic trajectories and related biogeochemical processes not necessarily simulated by the model. As such, physically-based modeling of water dynamics can, therefore, inform and advance digital soil mapping by linking water movement patterns stemming from hydrologic model outputs to spatial patterns of soil properties and pedogenesis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

3D xRAGE simulation of inertial confinement fusion implosion with imposed mode 2 laser drive asymmetry

Low-mode asymmetries represent an important obstacle to achieving high-gain inertial confinement fusion implosions. As a step in learning how to control such effects, an OMEGA experiment with imposed mode 2 laser drive asymmetries was done to study the expected signatures of this type of asymmetry [M. Gatu Johnson et al., PRE 2018]. In the present work, a 3D xRAGE simulation including the stalk mount has been brought to bear on the data from that experiment. Furthermore, comprehensive comparisons between simulated and measured observables are made. Good agreement between simulated and measured x-ray image-inferred shell trajectories, bang times and neutron emission widths are seen, showing that the hydrodynamics are well captured in the simulation. Asymmetries seen in simulated and measured time-resolved and time-integrated x-ray images and areal densities also compare well, showing impact of both stalk and mode 2. On the other hand, important differences in measured and simulated neutron emission histories, yield, and ion temperature (T ion ) asymmetries are seen, suggesting that the simulation is overestimating shock yield. Overall, the results clearly demonstrate the importance of considering all asymmetry sources when interpreting measured signatures of asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Design of Cylindrical Implosion Experiments to Demonstrate Scale-Invariant Rayleigh-Taylor Instability Growth

Radiation-hydrodynamics simulations are used to design laser-driven cylindrical implosion experiments to directly measure hydrodynamic instability growth in convergent geometry. Designs for two different size targets, varying in radial dimension by a factor of three, are presented. A set of beam pointings and powers are identified for each scale design that result in a nearly axially uniform implosion of an embedded marker layer. The implosion trajectories are shown to be scale-invariant between designs, with nearly identical scaled acceleration profiles. Finally, linear theory and radiation-hydrodynamics simulations predict that Rayleigh-Taylor instability growth of an azimuthal perturbation, machined on the inner surface of the embedded marker, is also scale-invariant between designs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Technique for an improved and time-resolved carbon areal density measurement with a single channel diagnostic for inertial confinement fusion

Measurement of the areal density and velocity of the carbon ablator shell during peak burn in inertial confinement fusion give powerful information on the state of the ablator and where in the trajectory of implosion it reaches peak burn. Detailed comparison of the absolute densities and velocities of the carbon in implosions has been prevented by the limited ability to resolve shot-to-shot variation within a shot series or within a campaign. Here, a new approach using a single, ultra fast (~10 ps) gamma ray channel can massively reduce uncertainties and will provide insights on improvements to target and drive variables. Small improvements in these experimental design parameters may result in much greater yields.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Effects of spatial variability in vegetation phenology, climate, landcover, biodiversity, topography, and soil property on soil respiration across a coastal ecosystem

Coastal terrestrial-aquatic interfaces (TAIs) are crucial contributors to global biogeochemical cycles and carbon exchange. A systematic evaluation of the interaction between coastal catchment properties and carbon dioxide (CO2) emission by soil respiration is significant for assessing carbon dynamics and predicting the future trajectory of atmospheric CO2 concentrations in coastal TAIs. The soil CO2 efflux in these transition zones is however poorly understood due to the high spatiotemporal dynamics of TAIs, as various sub-ecosystems in this region are compressed and expanded by complex influences of tides, changes in river levels, climate, and land use. We focus on the Chesapeake Bay region to (i) investigate the spatial heterogeneity of the coastal ecosystem and identify spatial zones with similar environmental characteristics based on the spatial data layers, including vegetation index (kNDVI), climate, landcover, diversity, topography, soil property, and relative tidal elevation; (ii) understand the primary driving factors affecting soil respiration within sub-ecosystems of the coastal ecosystem. Specifically, we employed hierarchical clustering analysis to identify spatial regions with distinct environmental characteristics, followed by the determination of main driving factors using Random Forest regression and SHapley Additive exPlanations. Maximum and minimum temperature are the main drivers common to all sub-ecosystems, while each region also has additional unique major drivers that differentiate them from one another. Precipitation exerts an influence on vegetated lands, while soil pH value holds importance specifically in forested lands. In croplands characterized by high clay content and low sand content, the significant role is attributed to bulk density. Wetlands demonstrate the importance of both elevation and sand content, with clay content being more relevant in non-inundated wetlands than in inundated wetlands. The topographic wetness index significantly contributes to the mixed vegetation areas, including shrub, grass, pasture, and forest. Additionally, our research reveals that dense vegetation land covers and urban/developed areas exhibit distinct soil property drivers. Overall, there is no one-size-fits-all approach to modeling carbon fluxes in coastal TAIs, and our study highlights the importance of further research and monitoring practices to improve our understanding of carbon dynamics and promote the sustainable management of coastal TAIs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

The role of source geometry and atmospheric propagation in global bolide infrasound detectability

Global infrasound monitoring provides a persistent means of detecting energetic bolide atmospheric entries, complementing optical observations and extending coverage over remote regions. We present a global assessment of the physical factors governing bolide infrasound detectability by correlating 623 bolide events reported by the Center for Near-Earth Object Studies between 2007 and 2025 with waveform data from the International Monitoring System. We identify 311 events with confirmed infrasound detections, corresponding to a detection rate of approximately 50%, substantially higher than inferred from earlier surveys, reflecting both the maturation of the global infrasound network and advances in automated, multi-frequency array processing. Analysis of flight parameters shows that infrasound detectability is selective rather than uniform across the bolide population. Detected events are preferentially associated with steeper entry angles and lower-altitude energy deposition, while shallow, high-altitude trajectories are less consistently observed. Very high-energy events remain detectable regardless of geometry, but for the more common lower-energy regime, observability depends on specific combinations of entry parameters and propagation conditions. This geometric dependence persists across comparable energy ranges and atmospheric conditions, indicating that entry angle exerts a primary control on detectability, with energy and propagation acting as secondary modulating factors. Furthermore, these results provide new physical constraints on bolide-atmosphere interactions and improve interpretation of global infrasound observations for planetary defense and atmospheric-entry studies.

Bolides↗
DOE OSTI2025

Numerical analysis of air dehumidification through electrospray-enhanced vortical flow cyclone separator

The dehumidification of air is a highly energy-intensive process, typically requiring vapor compression systems to cool the air significantly below its dew point temperature. Here, this study numerically investigates air dehumidification achieved by electrostatic and dielectrophoresis forces in a converging duct and centrifugal principles in cyclone separators. An electrospray was employed to energize water droplets, which then captured water vapor via dielectrophoresis. A novelty of the present study is the inclusion of a cyclone separator positioned downstream of the electrospray. This combination increased the dehumidification efficiency and controlled the flow of the condensate. By integrating electrospray technology and applying an electric field, this study examines the impact of particle separation efficiency, pressure drops, and cyclone performance on air dehumidification. The control parameters for air dehumidification were investigated. Three cyclones were modelled for airflow rates of 5 cubic feet per minute [cfm] (0.0024 m 3 /s), 38 cfm (0.0179 m 3 /s), and 200 cfm (0.094 m 3 /s). The results indicated air dehumidification was up to 5.7 % under environmental conditions, similar to those used to validate the CFD model. The degree of air dehumidification was primarily influenced by the electric charge acquired by the water droplets, their flight trajectory, and the vorticity intensity within the cyclone. Air dehumidification could increase by up to 13 % by increasing the number of droplets injected into the air stream, the flight time, and the electrical charge. The breakdown electric field strength threshold near the electrospray was a limiting factor for further enhancing dehumidification performance. This study examines the optimization of an electro-assisted dehumidification system that includes cyclone-based droplet separation through numerical analysis. The results provide insights into the design of energy-efficient air treatment technologies.

42 ENGINEERING↗
DOE OSTI2022

Locating fast-varying line disturbances with the frequency mismatch

In an attempt to provide an efficient method for line disturbance identification in complex networks of diffusively coupled agents, we recently proposed to leverage the frequency mismatch. The frequency mismatch filters out the intricate combination of interactions induced by the network structure and quantifies to what extent the trajectory of each agent is affected by the disturbance. In this previous work, we provided some analytical evidence of its efficiency when the perturbation is assumed to be slow. Here we claim that the frequency mismatch performs actually well for most disturbance regimes. This is shown through a series of simulations and is backed up by an analytical argument. Therefore, we argue that the frequency mismatch is an efficient and elegant tool for line disturbance location in complex networks of diffusively coupled agents.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Continuous-time look-ahead flexible ramp scheduling in real-time operation

This paper develops a novel continuous-time look-ahead optimization model for co-optimizing balancing energy and flexible ramp products in real-time power systems operation. In addition, a continuous-time day-ahead operation model is developed to commit and schedule an adequate subset of generating units for providing balancing energy and ramping requirements in real-time operation. The up and down flexible ramp products are defined as continuous-time decision trajectories that continuously supply the real-time up and down ramping requirements of the system. The deliverability of flexible ramp products is secured through reserving the respective power requirements in generation capacity constraints of generating units. A function space solution method is proposed to reduce the dimensionality of day-ahead and real-time models through projecting them into finite-dimensional function spaces spanned by Bernstein polynomials, where the real-time operation decisions are modeled with a higher accuracy level than day-ahead decisions. Simulations are conducted on the IEEE-RTS system using CAISO load data, and the numerical results showcase the effectiveness of the proposed model in supplying the energy and flexible ramp requirements and avoiding ramping scarcity in the system, while reducing the operation costs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Temperature statistics in a radiatively heated particle-laden turbulent square duct flow

Radiation absorption by preferentially concentrated particles in a turbulent square duct flow is studied experimentally. The particle-laden flow is exposed to near-infrared radiation, and the gas phase temperature statistics are measured along the wall bisector of the duct. It is found that the instantaneous temperature fluctuations are comparable to the overall mean temperature rise. The temperature statistics at the duct centerline and near the wall are qualitatively different. The former reflects preferential concentration in isotropic flows while the latter displays evidence of particle clustering into streamwise elongated streaks. Comparison of the experimental data to a simplified heat transfer model suggests that the Lagrangian evolution of particle clusters and voids, and turbulent mixing in the vicinity of particle clusters, are important. This work was motivated by particle solar receiver technology, but the findings are also relevant to systems where there is localized heat release or mass transfer from disperse particles or droplets. It shows that obtaining Lagrangian histories of particle trajectories is an important next step towards understanding thermal transport phenomena in particle-laden turbulent flows.

42 ENGINEERING↗
DOE OSTI2021

Experimental study of sprays with isooctane, hexane, ethanol and their binary mixtures under different flash boiling intensities

Here, this work investigated the behaviors of single-component and multi-component sprays under different flash boiling intensities. Isooctane, hexane, ethanol, isooctane-hexane, and isooctane-ethanol mixtures were used as test fuels. The pressure ratio of the ambient pressure over saturation pressure (R p ) was used as the indicator of the superheated degree. Macroscopic characteristics of the sprays under 10 different R p (1.1 to 0.14) were investigated using Diffused Backlight Imaging (DBI). Droplet size, which revealed the spray microscopic characteristics, was also investigated using Phase Doppler Anemometry (PDA). The results showed that ambient pressure reduction led to higher flash boiling intensity, significant spray morphology changes, and smaller droplet size, regardless of the fuel type and blend ratios. Spray collapse was observed for all sprays once R p dropped below 0.3. Spray collapse dramatically increased the liquid penetration length, expanded the near-field spray angle, and reduced the far-field angle for sprays with all the fuels. Fuel properties, such as the latent heat of vaporization, had significant impacts on droplet size distributions, thus altering the liquid penetration length and spray angles. Under similar R p , sprays with high ethanol content generally had larger droplet sizes compared with other fuel sprays due to the high latent heat of vaporization. The larger droplets caused the high ethanol content sprays to be more resistant to the plume trajectory change, resulting in a smaller far-field spray angle variations with changing R p . At R p of 0.28, sprays with the similar latent heat of vaporization have similar droplet sizes, and droplet sizes at different radial locations were also similar to each other.

42 ENGINEERING↗
DOE OSTI2021

3D optical diagnostics for explosively driven deformation and fragmentation

High-speed, optical imaging diagnostics are presented for three-dimensional (3D) quantification of explosively driven metal fragmentation. At early times after detonation, Digital Image Correlation (DIC) provides non-contact measures of 3D case velocities, strains, and strain rates, while a proposed stereo imaging configuration quantifies in-flight fragment masses and velocities at later times. Experiments are performed using commercially obtained RP-80 detonators from Teledyne RISI, which are shown to create a reproducible fragment field at the benchtop scale. DIC measurements are compared with 3D simulations, which have been ‘leveled’ to match the spatial resolution of DIC. Overall, results demonstrate improved ability to identify predicted quantities-of-interest that fall outside of measurement uncertainty and shot-to-shot variability. Similarly, video measures of fragment trajectories and masses allow rapid experimental repetition and provide correlated fragment size-velocity measurements. Measured and simulated fragment mass distributions are shown to agree within confidence bounds, while some statistically meaningful differences are observed between the measured and predicted conditionally averaged fragment velocities. Together these techniques demonstrate new opportunities to improve future model validation.

42 ENGINEERING↗
DOE OSTI2020

A Physically-Based Model for Thermo-Oxidative and Hydrolytic Aging of Elastomers

A computationally efficient model is proposed to capture the loss of mechanical performance due to chemical aging that are formed as the competition of chain scission and cross-link formation/dissolution, such as thermo-oxidative aging or hydrolytic aging. The model should be considered an extension of our recent models [1, 2, 3] which further simplifies the matrix behavior based on the assumption of independence of environmental and mechanical damage. The model uses this assumption to reduce the necessary material parameters needed to model constitutive and inelastic behavior of elastomers during aging. To this end, the model can provide accurate predictions of the material performance with the significantly fewer number of fitting parameters. The model is relevant for all decay mechanisms formed by the occurrence of two simultaneous micro-mechanisms; (i) formation/reduction of the cross-links, and (ii) chain scission, both of which are present in thermo-oxidation and hydrolytic aging. Assuming the alteration of the chain density along the aging trajectory is identical to the peroxide cross-link density for thermo-oxidation, and the change of the average molecular weight for hydrolysis, the strain energy of polymer matrix can be rewritten as a function of deformation, deformation history, storage time and aging temperature. Next, the modified network alteration model is formulated for implementation into Finite Element (FE) simulations. The model is built on the presumption of homogeneous and consistent oxygen/water absorption and thus is mainly relevant for relatively thin samples exposed to environmental loads for a long time. The proposed model includes only six physically inspired material parameters. Therefore, while it is computationally efficient, it shows good agreement with own experimental data, which performed on various range of accelerated aging temperatures and times. With respect to its computational efficiency, simplicity, accuracy, and interpret-ability, the model is the right choice for advanced implementations in FE programs.

42 ENGINEERING↗
DOE OSTI2022

Experimentally validated high-fidelity simulations of a liquid jet in supersonic crossflow

Here, utilizing recent advancements in computational schemes for compressible, multiphase flows, this work features a parametric study of a pure liquid jet in supersonic crossflow that involves simulating the atomization process for four values of momentum-flux ratio. These simulations are validated against experimental results measured with high-speed X-ray imaging, which confirm the accuracy of the numerical approach. Also, the effect of numerical resolution on some flow behavior is investigated, revealing convergence of the jet shape and surface instability wavelength. Analysis of the resulting sprays includes statistical descriptions of the liquid distribution, liquid structures created through breakup, interfacial instabilities, and dominant flow features. As the flowrate increases, the spray penetrates further, it becomes more disperse, and less liquid impacts the wall, but the droplet size distribution changes little. The wavelength of instabilities on the windward side of the jet diverges from measured trends in subsonic crossflows. In a visualization of the time-averaged flow, counter-rotating vortices are observed along the jet core and in the wake, affecting the process of primary atomization and early droplet trajectories.

42 ENGINEERING↗
DOE OSTI2024

The derivation of CRSS in pure Ti and Ti-Al alloys

The work focuses on the determination of the critical resolved shear stress (CRSS) in titanium (Ti) and titanium-aluminum (Ti-Al) alloys, influenced by an array of factors such as non-symmetric fault energies and minimum energy paths, dislocation core-widths, short-range order (SRO) effects which alter the local atomic environment, and tension-compression (T-C) asymmetry affected by intermittent slip motion. To address these multifaceted complexities, an advanced theory has been developed, offering an in-depth understanding of the mechanisms underlying slip behavior. The active slip systems in these materials are basal, prismatic, and pyramidal planes, with the latter involving both $\langle$a$\rangle$ and $\langle$c + a$\rangle$ dislocations. Each slip system is characterized by distinct Wigner-Seitz cell configurations for misfit energy calculations, varying partial dislocation separation distances, and unique dislocation trajectories—all critical to precise CRSS calculations. The theoretical CRSS results were validated against a comprehensive range of experimental data, demonstrating a strong agreement and underscoring the model's efficacy.

Critical stress↗