Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time scale separation”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Local primordial non-Gaussian bias from time evolution

Primordial non-Gaussianity (PNG) is a signature of fundamental physics in the early Universe that is probed by cosmological observations. Here, it is well known that the local type of PNG generates a strong signal in the two-point function of large-scale structure tracers, such as galaxies. This signal, often termed “scale-dependent bias” is a generic feature of modulation of gravitational structure formation by a large-scale mode. It is less well appreciated that the coefficient controlling this signal, b ϕ , is closely connected to the time evolution of the tracer number density. This correspondence between time evolution and local PNG can be simply explained for a universal tracer whose mass function only depends on peak height and, more generally, for nonuniversal tracers in the separate universe picture, which we validate in simulations. We also describe how to recover the bias of tracers subject to a survey selection function and perform a simple demonstration on simulated galaxies. Since the local PNG amplitude in n-point statistics ($f$ NL ) is largely degenerate with the coefficient b ϕ , this proof of concept study demonstrates that Galaxy survey data can allow for more optimal and robust extraction of local PNG information from upcoming surveys.

Sullivan, James M. [University of California, Berk↗
DOE OSTI2024

Electrically Powered High-Salinity Brine Separation Using Dimethyl Ether

Dewatering highly saline aqueous streams, from mining and geothermal leachates to industrial wastewater, is essential for effective resource recovery and safe disposal. Membraneless water extraction (MWE) uses a low-polarity solvent to separate water from concentrated aqueous solutions. In this study, we design a new MWE that uses dimethyl ether (DME) to selectively extract water from high-salinity brines, leveraging the volatility of DME to achieve rapid solvent recovery. By separating water and dissolved salts at a liquid–liquid interface, MWE minimizes the deleterious effects of scaling on vulnerable membrane and heat exchanger surfaces, reducing the need for extensive pretreatment and expensive materials. We begin by developing a computational framework for a multistage counterflow liquid–liquid contactor, which extracts water into DME, coupled with a multistage solvent regenerator that uses vapor compression to efficiently separate the desalinated water from the DME extractant. Excess Gibbs free energy and equation of state frameworks are used to model fluid phase equilibria in water–DME–sodium chloride (NaCl) mixtures, with interaction parameters estimated from experimental data. Incorporating equilibrium calculations into a system-scale computational model, we examine the performance of MWE using DME for the first time. Our analysis demonstrates that MWE can concentrate seawater desalination brine (>1.0 mol NaCl kg –1 ) to zero-liquid discharge salinities, with an energy consumption of under 50 kW h per m 3 of water extracted with a solvent recovery ratio greater than 99.9%. We highlight the importance of staging the vapor compression process to simultaneously minimize energy consumption while enabling brine concentration and product water solvent contamination. Finally, the thermodynamic framework developed here allows for the robust evaluation of new MWE solvents and systems for critical brine concentration and fractional precipitation applications.

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

Nanoscale covariance magnetometry with diamond quantum sensors

Nitrogen vacancy (NV) centers in diamond are atom-scale defects that can be used to sense magnetic fields with high sensitivity and spatial resolution. Typically, the magnetic field is measured by averaging sequential measurements of single NV centers, or by spatial averaging over ensembles of many NV centers, which provides mean values that contain no nonlocal information about the relationship between two points separated in space or time. In this work, we propose and implement a sensing modality whereby two or more NV centers are measured simultaneously, and we extract temporal and spatial correlations in their signals that would otherwise be inaccessible. We demonstrate measurements of correlated applied noise using spin-to-charge readout of two NV centers and implement a spectral reconstruction protocol for disentangling local and nonlocal noise sources.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Constraining dark matter substructure with Gaia wide binaries

ABSTRACT We use a catalogue of stellar binaries with wide separations (up to 1 pc) identified by the Gaia satellite to constrain the presence of extended substructure within the Milky Way galaxy. Heating of the binaries through repeated encounters with substructure results in a characteristic distribution of binary separations, allowing constraints to be placed independent of the formation mechanism of wide binaries. Across a wide range of subhalo density profiles, we show that subhaloes with masses $\gtrsim 65\, \mathrm{ M}_\odot$ and characteristic length scales similar to the separation of these wide binaries cannot make up 100 per cent of the Galaxy’s dark matter. Constraints weaken for subhaloes with larger length scales and are dependent on their density profiles. For such large subhaloes, higher central densities lead to stronger constraints. Subhaloes with density profiles similar to those expected from cold dark matter must be at least ∼5000 times denser than predicted by simulation to be constrained by the wide binary catalogue.

Astronomy & Astrophysics↗
DOE OSTI2024

FORCE Update 2024

The Framework for Optimization of Resources and Economics (FORCE) tool suite is the U.S. Department of Energy’s Nuclear Integrated Energy Systems (IES) Program flagship tool suite for technoeconomic IES analysis of IES. This tool suite is useful for analysis designed to evaluate and improve the technoeconomics of energy production systems, particularly for systems including nuclear technology. In this report, we document the development activity for the FORCE tool suite to extend its capabilities as performed during fiscal year 2024. In addition to reliability and accessibility, capability is one of the three standards guiding the development of the FORCE tool suite and the software codes that are its constituent parts. Extending the capabilities of the FORCE tool suite allows analysis both within the IES program as well as industry, university, and laboratory partners to perform analysis with more accuracy, insight, and impactful narrative. Four areas of capability development were the focus of activity this year: economic parameter uncertainty quantification, multiresolution analysis, components-to-optimization workflow automation, and statespace construction workflows for real-time optimal control. In economic parameter uncertainty quantification, the ability of HERON to capture risk due to scenarios (weather and energy demand uncertainty) was expanded to also include uncertainties in financial parameters such as capital cost or operation and maintenance costs. By including these sources of uncertainty, which are sometimes very large compared with scenario uncertainty, HERON is better able to capture the risk posed by investment in various IES technology. Because of this, analysts can also consider the reduction in risks that can be realized by choice of some technologies. In multiresolution analysis, development activity extended on work completed previously. In fiscal year 2023, methods for decomposing time series signals, such as demand, solar and wind availability, and price profiles, were analyzed and down-selected to those most effective at splitting signals into different resolutions. These resolutions allow considering the influence of different energy demand and supply behaviors across different time scales. For example, energy demand might be divided into seasonal, weekly, and hourly profiles. In fiscal year 2024, this preliminary work was extended and implemented within the Risk Analysis Virtual Environment (RAVEN) risk and uncertainty analysis platform, which is used throughout the FORCE framework. This development of the “multi-resolution time series analysis” (MR-TSA) module in RAVEN allows training synthetic history generators on complex time series. These synthetic history generators can then be used in HERON for generating scenarios that represent possible market and weather scenarios that can be analyzed on different time scales. We envision completing this work in the future, implementing multiresolution dispatch optimization strategies that can make the most beneficial use of these stratified time histories. In components-to-optimization workflow development, workflows for translating user inputs of components into algorithms for algebraic optimization were selected and implemented. Similar algorithms within the Holistic Energy Resource Optimization Network (HERON) were separated from the main code base of HERON and gathered with the components-to-optimization workflows in the new Dispatch Optimization Variable Engine (DOVE) software library. This modularization allows FORCE users to analyze dispatch optimization and energy system duty cycles independently of HERON, which previously was a burdensome task. Additionally, these dispatch optimization algorithms, set up in an independent library, can now be used across all software applications within FORCE, especially including the real-time optimal control software Optimization of Real-time Capacity Allocation (ORCA). Allowing FORCE software to share dispatch optimization algorithms within a single library allows for improved software maintenance and reliability. In statespace characterization workflow development, alternative workflows for optimizing dispatch with additional technical accuracy was the focus, particularly to improve the real-time optimization decision making in ORCA. Using algorithms and workflows initially developed for the Feasible Actuator Range Modifier (FARM), workflows for determining the statespace representation of IES were identified and demonstrated. The resulting dispatch optimization required a more robust optimization algorithm than that originally used in HERON (and moved to DOVE), which required adding an alternate workflow to DOVE that can more accurately match the behavior of physical systems using a partial differential equation representation. In conclusion, capability developments in the FORCE tool suite in fiscal year 2024 have improved the ability of the FORCE tool suite to perform

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Rapid charging made practical in graphite-based lithium batteries: surface-acoustic wave turbulent electrolyte mixing to overcome diffusion limited charging rates [Final Report]

One of the key limits in rapidly recharging a lithium-ion battery is the depletion of lithium ions within the electrolyte adjacent the anode during charging and long diffusion time to overcome this depletion. It also causes dendrite formation, inefficient use of the lithium, and battery degradation over many charge-discharge cycles. Because the liquid electrolyte remains quiescent and unmixed, this depletion layer’s depth rapidly grows to match the anode-cathode separation distance at even modest charge rates. The solution proposed by PI Prof. James Friend and Co-PI Prof. Ping Liu, both from the University of California, San Diego, is to mix the electrolyte and minimize the Li+ ion concentration gradient during charging. Even in the presence of the separator, the charging rate could then be significantly increased. They propose to accomplish this using surface acoustic wave (SAW)-driven acoustic streaming, a technique employing 10-mW fingernail-sized solid-state devices from the telecommunications industry to drive turbulent mixing to submicron length scales in a manner completely compatible with the typical 20700 and 18650 cells used in electric vehicles. 0.1. Turbulent acoustic streaming mixes the electrolyte during charging—even with separator Our proposed 100-MHz SAW device used to recirculate the electrolyte is compatible with lithium- ion battery electrochemistry, as it is made in our lab of single-crystal lithium niobate. Their solution is straightforward to drive from a DC power source alongside the signal provided during battery charging. The device is only needed during charging and does not consume power dur- ing battery discharge. Uniquely, SAW generates extreme accelerations of over 1 billion meters per second squared in the fluid, driving turbulent mixing from centimeter to submicron length scales, even through the porous separator materials that tend to be used in batteries, all while avoiding interfering with the anode’s solid electrolyte interphase layer that forms during use. By employing novel fluid mechanics, the investigators are proposing a new direction for battery research away from direct use of materials science and electrochemistry. The chemistry agnostic solution may be employed in any battery chemistry that makes use of liquid electrolytes, providing a broader transformative benefit to the battery research community. They are to produce a series of prismatic and 20700-cell sized 2 Ah batteries capable of being charged and discharged at least 500 times without more than 20% loss in battery capacity, and to provide analysis tools useful to the battery research and development community for adopting this approach for other battery chemistries and configurations.

25 ENERGY STORAGE↗
DOE OSTI2021

Mass-conserving implicit–explicit methods for coupled compressible Navier–Stokes equations

Earth system models are composed of coupled components that separately model systems such as the global atmosphere, ocean, and land surface. While these components are well developed, coupling them in a single system can be a significant challenge. Computational efficiency, accuracy, and stability are principal concerns. In this study we focus on these issues. In particular, implicit–explicit (IMEX) tight and loose coupling strategies are explored for handling different time scales. For a simplified model for the air–sea interaction problem, we consider coupled compressible Navier–Stokes equations with an interface condition. Under the rigid-lid assumption, horizontal momentum and heat flux are exchanged through the interface. Several numerical experiments are presented to demonstrate the stability of the coupling schemes. Here, we show both numerically and theoretically that our IMEX coupling methods are mass conservative for a coupled compressible Navier–Stokes system with the rigid-lid condition

42 ENGINEERING↗
DOE OSTI2025

Unified and optimal frame choice for generalized parton distributions

Reconstructing the internal three-dimensional quark and gluon structures of hadrons through generalized parton distributions (GPDs) from hard exclusive scattering processes is one of the most challenging tasks in nuclear and particle physics. In this paper, we introduce a new optimized reference frame that, for the first time, enables a unified view of all the reactions sensitive to GPDs and facilitates the interpretation of a variety of phase-space patterns that were previously hardly accessible and interpretable. Similarly to how the heliocentric description advanced our understanding of the solar system and gravitation, our new frame centers around a quasireal state, allows for a consistent separation of physical scales, and reveals a novel quantum interference mechanism. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Mitigating cosmic-ray-like correlated events with a modular quantum processor

Quantum processors based on superconducting qubits are being scaled to larger qubit numbers, enabling the implementation of small-scale quantum error-correction codes. However, catastrophic chip-scale correlated errors have been observed in these processors, attributed to, e.g., cosmic ray impacts, which challenge conventional error-correction codes such as the surface code. These events are characterized by a temporary but pronounced suppression of the qubit-energy relaxation times. Here, in this study, we explore the potential for modular quantum computing architectures to mitigate such correlated energy decay events. We measure cosmic-ray-like events in a quantum processor comprising a motherboard and two flip-chip bonded daughterboard modules, each module containing two superconducting qubits. We monitor the appearance of correlated qubit decay events within a single module and across the physically separated modules. We find that while decay events within one module are strongly correlated (over 85%), events in separate modules only display approximately 2% correlations. We also report coincident decay events in the motherboard and in either of the two daughterboard modules, providing further insight into the nature of these decay events. These results suggest that modular architectures, combined with bespoke errorcorrection codes, offer a promising approach for protecting future quantum processors from chip-scale correlated errors.

Wu, Xuntao [Univ. of Chicago, IL (United States)] ↗
DOE OSTI2021

An integrated numerical model for coupled poro-hydro-mechanics and fracture propagation using embedded meshes

Integrated models for fluid-driven fracture propagation and general multiphase flow in porous media are valuable to the study and engineering of several systems, including hydraulic fracturing, underground disposal of waste, and geohazard mitigation across such applications. This work extends the coupled model multiphase flow and poromechanical model of Ren et al. (2018) to admit fracture propagation (FP). The coupled XFEM-EDFM scheme utilizes a separate fracture mesh that is embedded on a static background mesh. The onset and dynamics of fracture propagation are governed by the equivalent stress intensity factor (SIF) criterion. A domain-integral method (J integral) is applied to compute this information. Additionally, an adaptive time-marching scheme is proposed to rapidly restrict and grow temporal resolution to match the underlying time-scales. The proposed model is verified with analytical solutions, and shows the capability to accurately and adaptively co-simulate fluid transport and deformation as well as the propagation of multiple fractures.

42 ENGINEERING↗
DOE OSTI2023

Measuring Photoexcited Electron and Hole Dynamics in ZnTe and Modeling Excited State Core-Valence Effects in Transient Extreme Ultraviolet Reflection Spectroscopy

Transient extreme ultraviolet (XUV) spectroscopy is becoming a valuable tool for characterizing solar energy materials because it can separate photoexcited electron and hole dynamics with element specificity. We use surface-sensitive femtosecond XUV reflection spectroscopy to separately measure photoexcited electron, hole, and band gap dynamics of ZnTe, a promising photocathode for CO 2 reduction. We develop an ab initio theoretical framework based on density functional theory and the Bethe-Salpeter equation to robustly assign the complex transient XUV spectra to the material's electronic states. Applying this framework, we identify the relaxation pathways and quantify their time scales in photoexcited ZnTe, including subpicosecond hot electron and hole thermalization, surface carrier diffusion, ultrafast band gap renormalization, and evidence of acoustic phonon oscillations.

14 SOLAR ENERGY↗
DOE OSTI2024

Initial energy-momentum to final flow: A general framework for heavy-ion collisions

The evolution of a relativistic heavy-ion collision is typically understood as a process that transmutes the initial geometry of the system into the final momentum distribution of observed hadrons, which can be described via a cumulant expansion of the initial distribution of energy density and is represented at leading order as the well-known eccentricity scaling of anisotropic flow. We extend this framework to include the contribution from initial momentum-space properties, as encoded in other components of the energy-momentum tensor. Furthermore, we confirm the validity of the framework in state-of-the-art hydrodynamic simulations of large and small systems. With this framework, it is possible to separate the effects of early time dynamics from those of final-state evolution, even in the case when the distribution of energy does not fully determine subsequent evolution, as for example, in small systems. Specifically, we answer the question of when and how azimuthal correlations from the initial state survive to the final state. In very small systems such as 𝑝−𝑝, for example, initial momentum degrees of freedom dominate over energy. Thus, even if the system forms a quark-gluon plasma that is well described by hydrodynamics, the usual hydrodynamic picture of the transmutation of initial geometry to final momentum anisotropy is broken. Nevertheless, we show that the hydrodynamic response to the full energy-momentum tensor can be well understood in a similar manner as larger systems. Additionally, this framework elucidates the generic features of the system's evolution that are responsible for the impressive success of hydrodynamic simulations, but which may still hold even in cases where hydrodynamics is not applicable.

Relativistic heavy-ion collisions↗
DOE OSTI2024

Galaxy bispectrum in the spherical Fourier-Bessel basis

The bispectrum, the three-point correlation in Fourier space, is a crucial statistic for studying many effects targeted by the next-generation galaxy surveys, such as primordial non-Gaussianity (PNG) and general relativistic (GR) effects on large scales. In this work we develop a formalism for the bispectrum in the spherical Fourier-Bessel (SFB) basis—a natural basis for computing correlation functions on the curved sky, as it diagonalizes the Laplacian operator in spherical coordinates. Working in the SFB basis allows for line-of-sight effects such as redshift space distortions and GR to be accounted for exactly, i.e., without having to resort to perturbative expansions to go beyond the plane-parallel approximation. Only analytic results for the SFB bispectrum exist in the literature given the intensive computations needed. We numerically calculate the SFB bispectrum for the first time, enabled by a few techniques: We implement a template decomposition of the redshift-space kernel Z 2 into Legendre polynomials, and separately treat the PNG and velocity-divergence terms. We derive an identity to integrate a product of three spherical harmonics connected by a Dirac delta function as a simple sum and use it to investigate the limit of a homogeneous and isotropic Universe. Furthermore, we present a formalism for convolving the signal with separable window functions and use a toy spherically symmetric window to demonstrate the computation and give insights into the properties of the observed bispectrum signal. While our implementation remains computationally challenging, it is a step toward a feasible full extraction of information on large scales via a SFB bispectrum analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Visualization of shocked material instabilities using a fast-framing camera and XFEL four-pulse train

Many questions regarding dynamic materials could be answered by using time-resolved ultra-fast imaging techniques to characterize the physical and chemical behavior of materials in extreme conditions and their evolution on the nanosecond scale. In this work, we perform multi-frame phase-contrast imaging (PCI) of micro-voids in low density polymers under laser-driven shock compression. At the Matter in Extreme Conditions (MEC) Instrument at the Linac Coherent Light Source (LCLS), we used a train of four x-ray free electron laser (XFEL) pulses to probe the evolution of the samples. To visualize the void and shock wave interaction, here, we deployed the Icarus V2 detector to record up to four XFEL pulses, separated by 1-3 nanoseconds. In this work, we image elastic waves interacting with the micro-voids at a pressure of several GPa. Monitoring how the material’s heterogeneities, like micro-voids, dictate its response to a compressive wave is important for benchmarking the performances of inertial confinement fusion energy materials. For the first time in a single sample, we have combined an ultrafast x-ray framing camera and four XFEL pulse train to create an ultrafast movie of micro-void evolution under laser-driven shock compression. Eventually, we hope this technique will resolve the material density as it evolves dynamically under laser shock compression.

fusion↗
DOE OSTI2026

A low-latitude, high-elevation 10Be production-rate calibration from the Rwenzori Mountains, Uganda

Cosmogenic nuclide surface-exposure dating has become a key tool in geomorphology. Robust application of the technique relies on determination of local nuclide production rates using independently dated calibration sites. However, there are few existing calibration sites in certain regions, particularly in low-latitude, high-elevation regions where cosmogenic nuclide production is theoretically more sensitive to changes in Earth's magnetic field and atmospheric conditions. We present a new cosmogenic beryllium-10 (10Be) production-rate calibration from the equatorial Rwenzori Mountains in Uganda. The calibration is based on radiocarbon dating of basal lake sediments (∼21.3 kyr BP) from Lake Mahoma (2900 m asl) and measured 10Be concentrations from boulders on the crest of a moraine which dams the lake. Our results indicate a time-integrated 10Be production rate of 15.83 ± 0.81 at/g/yr at Lake Mahoma. We then use these calibration data and two public online calculators to calculate ages for moraines in a separate Rwenzori valley for which there exists a minimum-limiting radiocarbon age constraint (∼12 kyr BP) on surficial glacial sediments. This enables us to assess the performance of different parameters used to scale production rates spatially and temporally (e.g., scaling frameworks, geomagnetic field reconstructions, atmospheric models) for surface-exposure age calculation in the region. Our results highlight the sensitivity of low-latitude cosmogenic nuclide surface-exposure ages to discrete calculation parameters, such as scaling framework or reconstructions of Earth's magnetic field, and the utility of determining production rates at sites where limiting age data exists for other nearby landforms.

Anderson, N [Dartmouth College]↗
DOE OSTI2021

Bayesian forecasts for dark matter substructure searches with mock pulsar timing data

Dark matter substructure, such as primordial black holes (PBHs) and axion miniclusters, can induce phase shifts in pulsar timing arrays (PTAs) measurements due to gravitational effects. In order to gain a more realistic forecast for the detectability of such models of dark matter with PTAs, we propose a Bayesian inference framework to search for phase shifts generated by PBHs and perform the analysis on mock PTA data. Furthermore, for most PBH masses the constraints on the dark matter abundance agree with previous (frequentist) analyses (without mock data) to $\mathcal{O}$(1) factors. This further motivates a dedicated search for PBHs (and dense small scale structures) in the mass range from 10 -8 M ⊙ to well above 10 2 M ⊙ with the Square Kilometer Array. Moreover, with a more optimistic set of timing parameters, future PTAs are predicted to constrain PBHs down to 10 -11 M ⊙ . Lastly, we discuss the impact of backgrounds, such as Supermassive Black Hole Mergers, on detection prospects, suggesting a future program to separate a dark matter signal from other astrophysical sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

First spectroscopic study of odd-odd Cu 78

Nuclei in the vicinity of 78 Ni are important benchmarks for nuclear structure, which can reveal changes in the shell structure far from stability. Spectroscopy of the odd-odd isotope 78 CU was performed for the first time in an experiment with the EURICA setup at the Radioactive Isotope Beam Factory at RIKEN Nishina Center. Excited states in the neutron-rich isotope were populated following the β decay of 78 Ni produced by in-flight fission and separated by the BigRIPS separator. A level scheme based on the analysis of $\gamma$-$\gamma$ coincidences is presented. Tentative spin and parity assignments were made when possible based on the β -decay feeding intensities and $\gamma$-decay properties of the excited states. Time correlations between β and $\gamma$ decay show clear indications of an isomeric state with a half-life of 3.8(4) ms. Large-scale Monte Carlo shell-model calculations were performed using the A3DA-m interaction and a valence space comprising the full f p shell and the 1⁢g 9/2 and 2⁢d 5/2 orbitals for both protons and neutrons. The comparison of the experimental results with the shell-model calculations allows interpreting the excited states in terms of spin multiplets arising from the proton-neutron interaction. In conclusion, the results provide further insight into the evolution of the proton single-particle orbitals as a function of neutron number, and quantitative information about the proton-neutron interaction outside the doubly magic 78 Ni core.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Will LISA Detect Harmonic Gravitational Waves from Galactic Cosmic String Loops?

Rapid advancement in the observation of cosmic strings has been made in recent years placing increasingly stringent constraints on their properties, with $G\mu\lesssim 10^{-11}$ from Pulsar Timing Array (PTA). Cosmic string loops with low string tension clump in the Galaxy due to slow loop decay and low gravitational recoil, resulting in great enhancement to loop abundance in the Galaxy. With an average separation of down to just a fraction of a kpc, and the total power of gravitational wave (GW) emission dominated by harmonic modes spanning a wide angular scale, resolved loops located in proximity are powerful, persistent, and highly monochromatic sources of GW with a harmonic signature not replicated by any other sources, making them prime targets for direct detection by the upcoming Laser Interferometer Space Antenna (LISA), whose frequency range is well-matched. Unlike detection of bursts where the detection rate scales with loop abundance, the detection rate for harmonic signal is the result of a complex interplay between the strength of GW emission, loop abundance, and other sources of noise, and is most suitably studied through numerical simulations. We develop a robust and flexible framework for simulating loops in the Galaxy for predicting direct detection of harmonic signal from resolved loops by LISA. Our simulation reveals that the most accessible region in the parameter space consists of large loops $\alpha=0.1$ with low tension $10^{-21}\lesssim G\mu\lesssim 10^{-19}$. Direct detection of field theory cosmic strings is unlikely, with the detection probability $p_{\mathrm{det}}\lesssim 2\%$ for a 1-year mission. An extension suggests that direct detection of cosmic superstrings with a low intercommutation probability is very promising. Searching for harmonic GW signal from resolved loops through LISA observations will potentially place physical constraints on string theory.

79 ASTRONOMY AND ASTROPHYSICS↗