Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time dependent behavior”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Solvent absorption rate of perfluorosulphonic acid membranes towards understanding direct coating processes

Here we present a method for measurement of the rate of solvent absorption by perfluorosulfonic acid (PFSA) membranes using a force tensiometer. The method presented here can be used as a tool to understand solvent absorption and should provide a rationale for designing catalyst inks for direct coating processes since it is necessary to understand how the absorption rate – on a time scale of seconds to minutes – compares to the time scales of the coating process in order to minimize membrane swelling. This method allows for rapid screening of the absorption and swelling behavior of different solvents and mixtures. Using this method, the absorption of water/1-propanol mixtures were measured for three thicknesses of Nafion PFSA membrane – Nafion 1135, Nafion 115, and Nafion 117. We find that the absorption rate is dependent on the ratio of the solvents as well as the thickness of the membrane. The analysis indicates that the highest absorption rate occurs when the mass percentage of 1-propanol in the mixture is 50%, whereas the lowest rates are for pure water and 1-propanol. In addition, membrane distortion also occurs most quickly for the 50% 1-propanol mixture. Our results suggest that catalyst inks that are highly rich (=90%) in water or 1-propanol are likely best for direct coating as such formulations will minimize absorption and swelling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Self-consistent hopping theory of activated relaxation and diffusion of dilute penetrants in dense crosslinked polymer networks

We generalize a microscopic statistical mechanical theory of the activated dynamics of dilute spherical penetrants in glass-forming liquids to study the influence of crosslinking in polymer networks on the penetrant relaxation time and diffusivity over a wide range of temperature and crosslink fraction (f n ). Our calculations are relevant to recent experimental studies of a nm-sized molecule diffusing in poly-(n-butyl methacrylate) networks. The theory predicts the penetrant relaxation time increases exponentially with the glass transition temperature, T g (f n ), which grows roughly linearly with the square root of fn due to the coupling of local hopping to longer-range collective elasticity. Moreover, T g is also found to be proportional to a geometric confinement parameter defined as the ratio of the penetrant diameter to the mean network mesh size. The decoupling ratio of the penetrant and Kuhn segment alpha times displays a complex non-monotonic dependence on f n and temperature that is well collapsed based on the variable T g (f n )/T. A model for the penetrant diffusion constant that combines activated relaxation and entropic mesh confinement is proposed, which results in a significantly stronger suppression of mass transport with degree of effective supercooling than predicted for the penetrant alpha time. Furthermore, this behavior corresponds to a new network-based type of “decoupling” of diffusion and relaxation. In contrast to the diffusion of larger nanoparticles in high temperature rubbery networks, our analysis in the supercooled regime suggests that for the penetrants studied the mesh confinement effects are of secondary importance relative to the consequences of crosslink-induced slowing down of activated hopping of glassy physics origin.

Chemistry↗
DOE OSTI2020

Melting temperature measurement of refractory oxide ceramics as a function of oxygen fugacity using containerless methods

An aerodynamic conical nozzle levitator (CNL) has been used to measure the melting temperatures of refractory ceramics. The well-established method of cooling traces has been adapted to allow measurements for air-sensitive samples using a controlled atmosphere, where the oxygen partial pressure can be adjusted and monitored in real time. In this study, we explain the details of a system that we developed for use at the Spallation Neutron Source (Oak Ridge National Laboratory). Measurements of the melting behavior of (Fe,Cr,Al) 3 O 4 spinels were made between log (PO 2 (atm)) of -0.7 and -5. The melting temperature of hercynite is found to be oxygen potential dependent, ranging between 1692 and 1767°C at log (PO 2 ) values from -4.42 to -0.70.

36 MATERIALS SCIENCE↗
NASA NTRS1983

Analytical and Experimental Evaluation of the Heat Transfer Distribution over the Surfaces of Turbine Vanes

Three airfoil data sets were selected for use in evaluating currently available analytical models for predicting airfoil surface heat transfer distributions in a 2-D flow field. Two additional airfoils, representative of highly loaded, low solidity airfoils currently being designed, were selected for cascade testing at simulated engine conditions. Some 2-D analytical methods were examined and a version of the STAN5 boundary layer code was chosen for modification. The final form of the method utilized a time dependent, transonic inviscid cascade code coupled to a modified version of the STAN5 boundary layer code featuring zero order turbulence modeling. The boundary layer code is structured to accommodate a full spectrum of empirical correlations addressing the coupled influences of pressure gradient, airfoil curvature, and free-stream turbulence on airfoil surface heat transfer distribution and boundary layer transitional behavior. Comparison of pedictions made with the model to the data base indicates a significant improvement in predictive capability.

Hylton, L. D.↗
NASA NTRS1988

On nonsteady accretion in stellar wind-fed X-ray sources

The paper reports on the results of hydrodynamical simulations of nonaxisymmetric gas flow past a finite-sized gravitating object. The asymmetry in the flow is assumed to be due to a transverse density gradient. The solutions are found to be highly time-dependent and exhibit phases in which a disk forms. The direction of circulation of the disk alternates between the direct and counter sense on time scales of approximately 1 hr. While the disk is present, the specific angular momentum is high, and the mass accretion rate, low. The angular momentum capture oscillates about a mean which secularly increases. Such temporal behavior may be relevant to those X-ray pulsing systems (e.g., Vela X-1) which exhibit fluctuations in the X-ray luminosity and pulse period on time scales much shorter than the orbital period of the binary system.

Taam, Ronald E.↗
NASA NTRSDate not supplied

Isolator Shock Train Response to Dual-Mode Scramjet Throttling-Like Forcing

Typical dual-mode scramjet designs employ an isolator component situated between the inlet and combustor. The isolator fluidically buffers the inlet from the combustor via the formation of a shock train, an interconnected series of shock-boundary layer interactions. The isolator is required during ramjet-engine operation and at the lower edges of the scramjet operational envelope. In an operational engine, transient behavior imposes fluidic dynamics due to attitude adjustments, combustor throttling, aberrant resonances within the combustor, etc. To accommodate these transients, the shock train adjusts on both temporal and spatial scales. Sufficiently large downstream pressure increases in the combustor play a critical role during unstart, whereby strong compression waves propagate upstream in the isolator to ultimately force the shock train out the inlet, if it cannot stabilize to the pressure increase. This critically reduces the engine mass flow rate, potentially causing a catastrophic reduction in thrust and loss of the vehicle. A nearly instantaneous unstart can, in principle, be predicted by a ratio of the normal shock pressure jump in an inviscid isolator. However, lower-magnitude perturbations, below this bound, also induce unstarts. Understanding the transient behavior inherent in unsteady isolator operation is critical to reducing safety margin in isolator length and refining engine control algorithms to improve engine performance. The Isolator Dynamics Research Lab (IDRL) is a cold-flow, direct-connect isolator facility ideal for assessing the fundamental fluid dynamics of shock trains. The internal flow of a scramjet isolator is modeled using a 2D, symmetric Mach 2.5 nozzle followed by a 695 mm long, 25.4x50.8 mm rectangular isolator test section. A blockage cone on a high-speed ball screw actuator downstream of the isolator is pulsed to simulate the throttling behavior in an operational combustor. Wall flush high-speed static pressure transducers detect the streamwise shock train within the isolator. Time delays between the onset of the throttling and the response motion of the shock train are analyzed to assess the controlling scales and parameters. From this analysis, the response time of the isolator shock train system is found to depend primarily on the isolator back pressure time-rate-of-change but quickly saturates to a maximum propagation speed along the isolator. These results inform dual-mode scramjet engine designers of the critical parameters and limits of an unstart resistant engine.

scramjet↗
NASA NTRS1984

Response of small amplitude slosh to time-dependent effective gravity

Small amplitude, linearized wave theory is used to derive a differential equation for the sloshing wave height under conditions in which a sloshing liquid drains out of a tank of circular but variable cross sectional area, allowing for changes of effective gravity with time. A few illustrative examples are given in which analytical results can be displayed and interpreted. Numerical integrations are presented for various examples of interest with respect to slosh behavior in the bottom dome of the Space Shuttle External Tank during thrust termination, particularly as this bears on scavenging of cryogens.

Traugott, S. C.↗
DOE OSTI2022

Computer Science Research Needs for Parallel Discrete Event Simulation (PDES)

Historically, scientific computing efforts have demonstrated the clear need for, and effective use of, supercomputing with traditional time-stepped simulations. Nevertheless, there are several areas in the mission spaces of the U.S. Department of Energy and other agencies waiting to tap advanced computing research using a different, discrete event style of modeling, simulation, and analysis. These span a wide spectrum of applications including energy grid resilience, urban planning and policy, transportation science, building technologies, emergency response and planning, environmental impact analysis, computational epidemiology, Internet communications, cyber security, and cyber-physical systems, to name only a few. Even within traditional scientific applications, the role of discrete event modes of execution is increasing in the form of new event-based mathematical solvers such as quantized state integration methods and discrete-continuous hybrid system solvers. Co-design of advanced supercomputing hardware systems is another area that exploits discrete event simulation at its core for effective analyses. Complex systems, entity behaviors and interconnections play a significant role in all these applications, which are mapped to large-scale models with discrete event formulations. To make advancements in all the aforementioned scientific areas, many technical aspects need to be more thoroughly studied and deeply understood in parallel discrete event simulation (PDES). The unique dynamics inherent in a discrete event modeling approach, by their very nature, intersect and influence the entire stack of the computing system, including (a) the unique nature of the instruction sets exercised in PDES workloads without a predominance of high-precision floating point operations, (b) virtual time-constrained multi-threaded execution of many logical processes per processor, (c) extremely variable and difficult to predict network traffic characteristics, (d) interfaces and inter-dependencies with machine learning and artificial intelligence codes at higher software layers, and (e) highly challenging load balancing needs, especially in effectively accounting for accelerated/extremely heterogeneous computing in current and future high-performance computing systems. Efficient and accurate parallel execution of PDES workloads is also dominated by challenges in dealing with their asynchronous concurrency fundamentally present at the model level. Conservative synchronization, optimistic/speculative synchronization, and their hybrid schemes open new questions in fundamental computer science with respect to reversibility of computation and prediction (lookahead) of behaviors inherent within model codes. On the implementation front, there are relatively few scalable, general-purpose parallel discrete event simulators in the world, and even fewer have been studied on emerging hardware platforms. To enable scientific advances using PDES, the research needs in computer science must also be pursued and met in the intersection of the algorithmic and hardware-aware aspects of scalable PDES engines. This report is aimed at capturing a computer science-oriented view of this important area of research in PDES, presenting a sample of important applications with their inherent discrete event technology elements. Needs are outlined in core areas of parallel discrete event research as well as cross-cutting directions in computer science research that positively impact scientific advancements across several important application areas. A selection of priority research opportunities in advanced computing for PDES is identified to serve as reference for key research topics and their order of importance for scientific advancements.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1999

Kinetics of template-directed pyrophosphate-linked dideoxyguanylate synthesis as a function of 2-MeImpdG and poly(C) concentration: insights into the mechanism

Aqueous solutions of deoxyguanosine 5'-monophosphate 2-methylimidazolide, 2-MeImpdG, yield primarily deoxyguanosine 5'-monophosphate, 5'dGMP, and pyrophosphate-linked dideoxyguanylate, dG5'ppdG, abbreviated G2p (see Chart 1). The initial rate of G2p formation, d[G2p]/dt in M h-1, determined at 23 degrees C, pH 7.8, 1.0 M NaCl and 0.2 M Mg2+ by timed high-performance liquid chromatography (HPLC) analysis, exhibits a second-order dependence on 2-MeImpdG concentration, [G]o, indicating a bimolecular mechanism of dimerization in the range 0.02 M < or = [G]o < or = 0.09 M. In the presence of polycytidylate, poly(C), G2p synthesis is accelerated and oligodeoxyguanylate products are formed by incorporation of 2-MeImpdG molecules. The kinetics of G2p formation as a function of both monomer and polymer concentration, expressed in C equivalents, were also determined under the above conditions and exhibited a complex behavior. Specifically, at a constant [poly(C)], values of d[G2p]/dt typically increased with [G]o with a parabolic upward curvature. At a constant [G]o, values of d[G2p]/dt increase with [poly(C)], but level off at the higher poly(C) concentrations. As [G]o increases this saturation occurs at a higher poly(C) concentration, a result opposite to expectation for a simple complexation of two reacting monomers with the catalyst prior to reaction. Nevertheless, these results are shown to be quantitatively consistent with a template-directed (TD) mechanism of dimerization where poly(C) acts as the template to bind 2-MeImpdG in a cooperative manner and lead, for the first time, to the formulation of principles that govern template-directed chemistry. Analysis of the kinetic data via a proposed TD cooperative model provides association constants for the affinity between polymer and monomer and the intrinsic reactivity of 2-MeImpdG toward pyrophosphate synthesis. To the best of our knowledge, poly(C)/2-MeImpdG is the first system that could serve as a textbook example of a TD reaction under conditions such that the template is fully saturated by monomers and under conditions that it is not.

NASA Discipline Exobiology↗
DOE OSTI2026

Temperature‐Dependent Crystallization in Two‐Step Perovskite Deposition Revealed by In Situ GIWAXS and Machine Learning‐Guided Analysis

The performance and stability of perovskite solar cells are strongly governed by the crystallization behavior of their active layer. In two-step sequential deposition, early-stage film formation plays a decisive role in determining final phase purity and device quality. Guided by a data-driven analysis of nearly 39 000 devices in the FAIR perovskite database, we identified solvent-mediated quenching and thermal processing as key variables affecting power conversion efficiency (PCE), particularly in two-step fabrication. Here, to investigate these effects in real time, we designed and implemented a custom-built, temperature-controlled spin-coating system, enabling precise thermal modulation during precursor deposition. Using this platform, we performed in situ GIWAXS measurements to study the crystallization dynamics of FA 0.5 MA 0.5 PbI 3 films over a temperature range of 30°C–90°C. Our results reveal a non-monotonic relationship between spin-coating temperature and α-phase formation, governed by the interplay between precursor interdiffusion, PbI 2 crystallinity, and δ-phase suppression. The custom thermal control enabled us to isolate and quantify these competing effects during the earliest stages of film formation, providing mechanistic insight into how spin-coating temperature governs both phase purity and kinetic pathways in two-step perovskite systems. Temperature-dependent SEM and photovoltaic device measurements further demonstrate that early-stage crystallization pathways directly translate into differences in morphology, charge-transport continuity, and device performance. These findings inform targeted strategies for optimizing deposition protocols to balance rapid nucleation, phase stability, and device performance.

Saadawy, Ahmed [King Fahd University of Petroleum ↗
NASA NTRS1985

Viscoelastic behavior and lifetime (durability) predictions

A method for lifetime or durability predictions for laminated fiber reinforced plastics is given. The procedure is similar to but not the same as the well known time-temperature-superposition principle for polymers. The method is better described as an analytical adaptation of time-stress-super-position methods. The analytical constitutive modeling is based upon a nonlinear viscoelastic constitutive model developed by Schapery. Time dependent failure models are discussed and are related to the constitutive models. Finally, results of an incremental lamination analysis using the constitutive and failure model are compared to experimental results. Favorable results between theory and predictions are presented using data from creep tests of about two months duration.

Brinson, R. F.↗
DOE OSTI2020

The Development of a Generalized Riser Flow Regime Map Based Upon Higher Moment and Chaotic Statistics Using Electrical Capacitance Volume Tomography (ECVT)

Dynamic analyses have been applied to the temporal signals from an Electro Capacitance Volume Tomography instrument located near mid-height on the riser of an industrial-scale cold-flow circulating fluidized bed to characterize gas-solids flow behavior in the riser. Twelve capacitance electrodes surround the cylindrical riser over a height of 1.3 m. The instrument used a neural network deconvolution algorithm to determine the spatially resolved solids fraction recorded at 52 Hz. Experiments were carried out over a range of gas and solids flows in the transport regime using a Geldart Group B bed material, high density polyethylene with mean particle size of 880 μm. The radial solids distribution was found to vary from one-time step to the next between profiles typical of laminar and turbulent flow. The duration of time spent in each of these flow profiles depended upon the operating regime – dilute, core-annular, or fast fluidized bed. The chaotic structure of the temporal data was characterized using the three conventional approaches: the first 4 moments from the distribution of signal in time, system memory parameters from the autocorrelation function and the Hurst exponent, and analysis of the correlationentropy and correlation dimension of the attractor. These signal analysis techniques were used to clearly distinguish differences between different transport operating regimes. Specifically, it was experimentally observed that a riser transitions from core annular flow profile to dilute and dense regimes via increasing the frequency of short term transients to either dilute or dense flow profiles, respectively. A regime map was generated based upon these dynamics using solids flux and gas velocity axes. Fast fluidized, core annular, and dilute each exhibited different degree of dynamic characteristics typical of fluid dominated or particle compromising behavior. It should be noted that the magnitude for the different statistics was in the same range regardless of the regime, it was the radial profile for the statistic that changed and subsequently identified that there was a change in the regime. Finally, a reduced regime map was developed consisting of plotting the gas velocity normalized by the upper transport velocity versus the solids flux normalized by the saturation carrying capacity. The use of this reduced plot allowed the data from widely different conditions to be plotted and compared on the same<p>graph. Note that in many instances, some of the statistics identified the operating point as being in one regime while others indicated that it was in another indicating a transition region between dilute or core annular regimes and between the core annular and fast fluidization regimes. This now provides a tool that can be used to optimize process performance, identify changes in operating states, or replicate process dynamics during process scaling or changing operating parameters. </p>

Breault, Ronald↗
DOE OSTI2025

Effects of Strain and Strain Rate on Dynamic Grain Growth and Subgrain Evolution During Plastic Deformation of an Interstitial-Free Steel at 850 ° C

Here, the effects of strain and strain rate on dynamic grain growth (DGG) and subgrain evolution are reported for an interstitial-free steel deformed at 850 ° C. Microstructures produced during tension tests at true-strain rates of 10 -4 and to 10 -3 s -1 true strains ranging from 0.02 to 0.2 were preserved following deformation. These were characterized using electron backscatter diffraction (EBSD), including the application of spherical harmonic transform indexing to produce high-angular-resolution EBSD (HR-EBSD) data. HR-EBSD data resolved the small misorientation angles of subgrain boundaries while imaging much larger data fields than possible with previously available techniques. The resulting data confirmed that steady-state flow stress is inversely proportional to the average subgrain size and that subgrain boundary misorientation angle increases with strain. The following new observations are reported. The rate of DGG increased with respect to time but decreased with respect to strain as strain rate increased. This behavior is rationalized through a simple model using separate rate parameters for the effects of time and strain. Subgrain size was not constant during steady-state deformation, but decreased slowly with increasing strain. Subgrain size distributions and subgrain boundary misorientation angle distributions were measured, and both remained approximately log-normal during steady-state deformation. Subgrain evolution demonstrated no dependence on parent grain size, crystallographic orientation, or Taylor factor. These new data suggest that steady-state flow stress is more likely controlled by the dislocation density internal to subgrains than by the spacing between subgrain boundaries.

dynamic grain growth↗
DOE OSTI2021

Adoption of delivery services in light of the COVID pandemic: Who and how long?

A significant growth in demand for online shopping in light of the Coronavirus Disease (COVID) crisis has received attention from transportation practitioners, policy-makers, and researchers. However, an important question arises in this increase in online shopping and resulting deliveries: How long will this last? Very little is known whether this popularity would last a long time. Here, to address this question, the authors conducted a survey of 915 individuals residing in the U.S. and classified them into the four distinctive consumer types (i.e., the prior adopter, temporary adopter and permanent new adopter, and non-adopter) depending on their usage of delivery services before, during, and after (expected) the COVID crisis. This research aims to gain behavioral insight by exploring the differences between the four consumer types and investigating factors affecting the initial adoption and continuance intention of using delivery services. The descriptive analysis revealed that there are clear differences not only between the four types of consumers but also between the four product types (i.e., grocery, food, home goods, and other packages) considered in the survey. The models found that factors affecting the initial adoption and continuance intention are different from the previous studies conducted before the COVID pandemic. Implications for planning and policymaking are also discussed.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Synthesis and Physical Properties of Manganese Chromium Nitride Thin Films Grown via Molecular Beam Epitaxy

Here, we demonstrate the ability to grow epitaxial thin films of manganese–chromium nitride (Mn x Cr 1–x N) across the entire composition range on magnesium oxide (MgO) (1 0 0) surfaces by molecular beam epitaxy. By independently controlling the ratio of the atomic fluxes of manganese and chromium, the compositional dependence of structural and electrical transport properties is determined. Soft X-ray absorption shows that the valence states of Cr( 3+ ) and Mn( 3+ ) remain identical among these MnCrN epitaxial films. The metallic transport behavior of binary manganese nitride (MnN) and chromium nitride (CrN) is suppressed in the ternary alloys. Changes in room temperature electrical resistivity of up to 100 times and 5 orders of magnitude at cryogenic temperatures are observed. This work expands the variety of nitride compounds available for epitaxial integration and could provide additional flexibility for the development of metal–semiconductor–insulator junctions for spintronic devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Elucidating size effects on the yield strength of single-crystal Cu via the Richtmyer–Meshkov instability

Capturing the dynamic response of a material under high strain-rate deformation often demands challenging and time consuming experimental effort. While shock hydrodynamic simulation methods can aid in this area, a priori characterizations of the material strength under shock loading and spall failure are needed in order to parameterize constitutive models needed for these computational tools. Moreover, parameterizations of strain-rate-dependent strength models are needed to capture the full suite of Richtmyer–Meshkov instability (RMI) behavior of shock compressed metals, creating an unrealistic demand for these training data solely on experiments. Herein, we sweep a large range of geometric, crystallographic, and shock conditions within molecular dynamics (MD) simulations and demonstrate the breadth of RMI in Cu that can be captured from the atomic scale. In this work, yield strength measurements from jetted and arrested material from a sinusoidal surface perturbation were quantified as Y RMI = 0.787 ± 0.374 GPa, higher than strain-rate-independent models used in experimentally matched hydrodynamic simulations. Defect-free, single-crystal Cu samples used in MD will overestimate Y RMI , but the drastic scale difference between experiment and MD is highlighted by high confidence neighborhood clustering predictions of RMI characterizations, yielding incorrect classifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Superconductivity of incoherent electrons in the Yukawa Sachdev-Ye-Kitaev model

We study a model of N fermions in a quantum dot, coupled to M bosons by a disorder-induced complex Yukawa coupling [Yukawa Sachdev-Ye-Kitaev (SYK) model], in order to explore the interplay between non-Fermi liquid and superconductivity in a strongly coupled, (quantum-)critical environment. We analyze the phase diagram of the model for an arbitrary complex interaction and arbitrary ratio of N/M, with special focus on the two regimes of non-Fermi-liquid behavior: an SYK-like behavior with a power-law frequency dependence of the fermionic self-energy and an impuritylike behavior with frequency independent self-energy. We show that the crossover between the two can be reached by varying either the strength of the fermion-boson coupling or the ratio M/N. We next argue that in both regimes the system is unstable to superconductivity if the strength of time-reversal-symmetry-breaking disorder is below a certain threshold. We show how the corresponding onset temperatures vary between the two regimes. We argue that the superconducting state is highly unconventional with an infinite set of minima of the condensation energy at T = 0, corresponding to topologically different gap functions. We discuss in detail similarities and differences between this model and the model of dispersion-full fermions tuned to a metallic quantum-critical point, with an effective singular dynamical interaction V(Ω) ∝ 1/|Ω| γ (the γ model).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Angle-resolved polarized Raman spectroscopy study of phosphorene nanoribbons

We present a systematic angle-resolved polarized Raman spectroscopy (ARPRS) study of black phosphorus (BP) nanostructures formed via electrochemical sodium- and lithium-ion intercalation. Sodium intercalation leads to bundles of densely packed, highly uniform phosphorene nanoribbons (PNRs) separated by parallel amorphous channels, whereas lithium intercalation results in shorter, irregular nanoribbon-like segments with lower aspect ratios. In both cases, six additional Raman peaks (P1–P6) appear alongside the three primary Raman-active modes of BP (A$^{1}_{g}$, B 2g , and A$^{2}_{g}$). These peaks are attributed to the amorphous regions, as confirmed by their isotropic angular dependence in ARPRS measurements. The three BP modes show pronounced angular variations that differ significantly between the two intercalated samples. In sodium-intercalated BP, A$^{1}_{g}$ and A$^{2}_{g}$ modes retain a dumbbell-like angular dependence under parallel polarization with enhanced anisotropy and reduced symmetry under crossed polarization. At the same time, the B 2g mode transitions from four-lobed (cloverleaf) polar plot to a butterfly-like one. In contrast, lithium-intercalated BP exhibits weaker anisotropy and less distinct angular polar plots for all three modes. These differences reflect the sensitivity of phonon behavior to underlying nanostructure morphology. The vibrational frequencies density of states (FDOS) calculations attribute the B 2g mode transformation to phonon band folding and mode mixing in PNRs. This study demonstrates the power of ARPRS in probing phonon-structure relationships and highlights the influence of edge geometry and quantum confinement on phonon dispersion in PNRs.

77 NANOSCIENCE AND NANOTECHNOLOGY↗