Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reach code”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Additive Manufacturing: The Rise of a Technology and the Need for Quality Control and Inspection Techniques

Additive manufacturing (AM) technology, also known as three-dimensional (3D) printing, rapid prototyping (RP), or solid free forms (SFF), was developed back in the early 1980s. The technology was quickly commercialized in the late 1980s, and since then, the technology has been widely used in different applications at several industrial sectors. AM offers distinct advantages from traditional, subtractive manufacturing techniques. Traditional subtractive manufacturing techniques, including milling, drilling, and turning, are techniques that start with a large block of stock material, which is then subtracted from until the final desired product is reached, leaving much of the initial block as wasted scrap. In the AM process, the part is constructed by depositing material layer by layer in the Z-direction until the final product is produced, leaving little to no waste. The process starts with a computer-aided design (CAD) model that gets sliced into individual layers using software that generates instructions, known as G-code, which are sent to the AM machine. Here, this process allows complex geometry construction that traditional methods would have diffi-culty accomplishing or fail completely. In some cases, a secondary finishing process is needed such as polishing, sintering, curing, sanding, powder removal, or painting.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Decomposition Algorithms for Scalable Quantum Annealing

The presented python code provides wrapper functions for two decomposition algorithms: One algorithm to decompose Maximum Clique problems and one to decompose Minimum Vertex Cover problems. The functions take as input a networkx graph object, and decompose either problem on the input graph recursively into subproblems such that the optimal solution can be constructed from the optimal solutions of both subproblems. The recursion ends as soon as the subproblems reach a pre-specified size limit by the user, and they can be solved using any method provided in advance by the user as external function. This includes exact solvers or an adiabatic quantum annealer.

Pelofske, Elijah↗
DOE OSTI2017

ULTRAVIOLET ESCAPE FRACTIONS FROM GIANT MOLECULAR CLOUDS DURING EARLY CLUSTER FORMATION

The UV photon escape fraction from molecular clouds is a key parameter for understanding the ionization of the interstellar medium and extragalactic processes such as cosmic reionization. We present the ionizing photon flux and the corresponding photon escape fraction ( f {sub esc}) arising as a consequence of star cluster formation in a turbulent, 10{sup 6} M {sub ⊙} giant molecular cloud, simulated using the code FLASH. We make use of sink particles to represent young, star-forming clusters coupled with a radiative transfer scheme to calculate the emergent UV flux. We find that the ionizing photon flux across the cloud boundary is highly variable in time and space due to the turbulent nature of the intervening gas. The escaping photon fraction remains at ∼5% for the first 2.5 Myr, followed by two pronounced peaks at 3.25 and 3.8 Myr with a maximum f {sub esc} of 30% and 37%, respectively. These peaks are due to the formation of large H ii regions that expand into regions of lower density, some of which reaching the cloud surface. However, these phases are short-lived, and f {sub esc} drops sharply as the H ii regions are quenched by the central cluster passing through high-density material due to the turbulent nature of the cloud. We find an average f {sub esc} of 15% with factor of two variations over 1 Myr timescales. Our results suggest that assuming a single value for f {sub esc} from a molecular cloud is in general a poor approximation, and that the dynamical evolution of the system leads to large temporal variation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Buoyancy-driven air flow within plenum-to-plenum facility down-comer channel

Reliable measurements of thermal and velocity fields were experimentally investigated under different natural circulation intensities in a dual channel facility designed and developed with a representative geometry of prismatic modular reactor core. Experiments were conducted under steady state conditions to obtain statistically stationary time series signals of temperature and velocity. Insignificant radial temperature variation was observed along the down-comer channel implying that the majority of supplied heat is removed through the upper plenum. Velocity results obtained emphasize the delicate nature of natural circulation phenomena in terms of flow destabilization, and recirculation penetration length. Observed air velocity distribution reveals that downward velocity is proportional to the extent of cooling applied to the setup. The calculated values of the dimensionless group Froude number ( Fr ) in the range of 2 to 10 indicate that air is dominated by momentum and is exiting the channel as a jet. Quantification of turbulent intensities implies that flow destabilization is noticeable along the down-comer channel particularly for the case of high cooling intensity. For 5 °C cooling water temperature, turbulent intensity ($v$ rms / $v$) peaks at the mid-channel (z/L = 0.5) reaching a value of 0.5. This value of turbulent intensity is much higher than the inlet (z/L = 0.96) and outlet (z/L = 0.04) turbulent intensities of 0.2 and 0.06 respectively. This destabilization is because of the flow reversal and heat conduction effects through stainless steel flanges connecting both channels. In conclusion, current results provide detailed velocity and temperature variation which can be useful for validating computational fluid dynamics codes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Numerical challenges for energy conservation in N -body simulations of collapsing self-interacting dark matter halos

Dark matter (DM) halos can be subject to gravothermal collapse if the DM is not collisionless, but engaged in strong self-interactions instead. When the scattering is able to efficiently transfer heat from the centre to the outskirts, the central region of the halo collapses and reaches densities much higher than those for collisionless DM. This phenomenon is potentially observable in studies of strong lensing. Current theoretical efforts are motivated by observations of surprisingly dense substructures. However, a comparison with observations requires accurate predictions. One method to obtain such predictions is to use N-body simulations. Collapsed halos are extreme systems that pose severe challenges when applying state-of-the-art codes to model self-interacting dark matter (SIDM). In this work, we investigate the root of such problems, with a focus on energy non-conservation. Moreover, we discuss possible strategies to avoid them. We ran N-body simulations, both with and without SIDM, of an isolated DM-only halo and we adjusted the numerical parameters to check the accuracy of the simulation. We find that not only the numerical scheme for SIDM can lead to energy non-conservation, but also the modelling of gravitational interaction and the time integration are problematic. The main issues we find are: (a) particles changing their time step in a non-time-reversible manner; (b) the asymmetry in the tree-based gravitational force evaluation; and (c) SIDM velocity kicks breaking the time symmetry. Tuning the parameters of the simulation to achieve a high level of accuracy allows us to conserve energy not only at early stages of the evolution, but also later on. However, the cost of the simulations becomes prohibitively large as a result. Some of the problems that make the simulations of the gravothermal collapse phase inaccurate can be overcome by choosing appropriate numerical schemes. However, other issues still pose a challenge. Our findings motivate further works on addressing the challenges in simulating strong DM self-interactions.

dark matter↗
DOE OSTI2020

Stable heat and particle flux detachment with efficient particle exhaust in the island divertor of Wendelstein 7-X

The island divertor concept is an innovative and promising idea to handle heat and particle exhaust in stellarators. At the Wendelstein 7-X (W7-X) stellarator, this divertor concept plays a central role in the device mission to demonstrate reactor relevant plasma confinement for steady-state time scales of up to 30 minutes in the high-performance campaign (OP2) starting in 2022. During the recently concluded first campaign with the inertially cooled island divertor, a large step in the experimental qualification of this divertor concept has been made. In discharges heated with electron cylotron resonance heating of 5-6 MW, central densities in the range of 0.7-1.2 × 10 20 m -3 have been reached in combination with full divertor heat flux detachment. Also, significant neutral gas pressures and neutral compression ratios were shown for the first time in combination with reduced divertor particle flux. The divertor heat loads drop by an order of magnitude from >5 MW m -2 to below 0.5 MW m-2 with increasing density, and substantial compression of neutrals reaching neutral pressure in the sub-divertor volume of >6.0 × 10 -4 mbar was seen. These elevated neutral pressure levels can be obtained and maintained with an up to 80% reduction of the particle fluxes onto the divertor target tiles. This discharge scenario was held stably detached for up to 28 seconds, which is equivalent to several hundred energy confinement times τ E and longer than the time scales for current relaxation. No impurity accumulation was seen at constant Z eff ≈ 1.5 and the stored energy stayed constant at levels of W dia >600 kJ. The level of neutral pressure and compression reached in this scenario extrapolates well to the steady-state particle exhaust requirements for high-performance steady-state operation in OP2, in which the fully actively cooled high-heat-flux divertor will be available. Finally, an overview of this recently discovered divertor regime is given and the status of the physics understanding based on modeling of these regimes with the EMC3-EIRENE code is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Designing a validation experiment for radio frequency condensation

Abstract Theoretical studies have suggested that nonlinear effects can lead to ‘radio frequency (RF) condensation’, where an initially broad current profile can coalesce in islands when they reach sufficient width. In suitable conditions, RF condensation can ‘self-focus’ the driven current to the center of an island, improving stabilization efficiency and reducing control complexity. In unsuitable conditions, the effect can prematurely deplete the RF energy before it reaches the island center, impairing stabilization. It is predicted that the RF condensation effect can significantly impact reactor-scale tokamaks. This paper presents a set of simulations investigating the conditions under which RF condensation might be encountered in present-day tokamaks. For concreteness, the calculations use equilibrium reconstructions for two shots from DIII-D and AUG. The Current Condensation Amid Magnetic Islands (OCCAMI) simulation code has been used for this investigation. The code takes as its input a numerically specified axisymmetric EFIT equilibrium solution, and it perturbatively constructs a 3D field with an island embedded at the appropriate rational surface. In the OCCAMI code, the GENRAY code is used for ray tracing and for calculating the power deposition along a ray trajectory, and GENRAY is coupled self-consistently to a solution of the thermal diffusion equation in the island. The simulation results described in the paper illuminate the conditions required for experimental validation of the theory of RF condensation. The simulations also provide an explanation of why the effect was not noticed in experiments prior to the publication of theoretical papers on the subject.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Magnetically driven coupling in relativistic radiation-mediated shocks

ABSTRACT The radiation drag in photon-rich environments of cosmic explosions can seed kinetic instabilities by inducing velocity spreads between relativistically streaming plasma components. Such microturbulence is likely imprinted on the breakout signals of radiation-mediated shocks. However, large-scale, transverse magnetic fields in the deceleration region of the shock transition can suppress the dominant kinetic instabilities by preventing the development of velocity separations between electron–positron pairs and a heavy ion species. We use a 1D five-fluid radiative transfer code to generate self-consistent profiles of the radiation drag force and plasma composition in the deceleration region. For increasing magnetization, our models predict rapidly growing pair multiplicities and a substantial radiative drag developing self-similarly throughout the deceleration region. We extract the critical magnetization parameter σc, determining the limiting magnetic field strength at which a three-species plasma can develop kinetic instabilities before reaching the isotropized downstream. For a relativistic, single ion plasma drifting with γu = 10 in the upstream of a relativistic radiation-mediated shock, we find the threshold σc ≈ 10−7 for the onset of microturbulence. Suppression of plasma instabilities in the case of multi-ion composition would likely require much higher values of σc. Identifying high-energy signatures of microturbulence in shock breakout signals and combining them with the magnetization limits provided in this work will allow a deeper understanding of the magnetic environment of cosmic explosions like supernovae, gamma-ray bursts, and neutron star binary mergers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

20 years of the CEA/DAM NNSA/DP Agreement

For the past twenty years, there has been a very active, productive International Agreement between France and the United States of America for Cooperation on Fundamental Science supporting Stockpile Stewardship. Under this Agreement the scientists at the nuclear weapons laboratories in both countries have collaborated on many unclassified research projects in areas such as Materials in Extreme Conditions, Nuclear Physics, and Atomic and Plasma Physics. The results of their efforts have not only been published in the open literature but have enhanced the physics base of the computer codes essential to the mission of ensuring that the nuclear stockpiles are safe, secure and effective. Work on these collaborations is extremely important. Not only does such joint work bring more brilliant minds to work on pressing research problems for both countries but also the collaborative effort sharpens the scientific skills of and presents scientific challenges to the scientific and technical staff of the laboratories. As we reach the 20th Anniversary of the formal signing of the Agreement, we feel it is important to thank all of the individuals who have contributed to its continuing success. Our expectations for brilliant, exciting, challenging joint research projects under this Agreement are at an all-time high. We toast the achievement of this milestone and look forward to the research results to come.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Synthetic Spectra of Pair-instability Supernovae in 3D

Pair-instability supernovae (PISNe) may signal the deaths of extremely massive stars in the local Universe or massive primordial stars after the end of the Cosmic Dark Ages. Hydrodynamic simulations of these explosions, performed in 1D, 2D, and 3D geometry, have revealed the strong dependence of mixing in the PISN ejecta on dimensionality. This chemical rearrangement is mainly driven by Rayleigh–Taylor instabilities that start to grow shortly after the collapse of the carbon–oxygen core. We report on the effects of such mixing on the spectroscopic evolution of PISNe by post-processing explosion profiles with the radiation diffusion-equilibrium code SNEC and the implicit Monte Carlo–discrete diffusion Monte Carlo radiation transport code SuperNu . The first 3D radiation transport calculation of a PISN explosion is presented, yielding viewing-angle-dependent synthetic spectra and light curves. Our findings show that while 2D and 3D mixing does not significantly affect the light curves of PISNe, their spectroscopic and color evolution are impacted. Strong features of intermediate-mass elements dominated by silicon, magnesium, and oxygen appear at different phases and reach different intensities depending on the extent of mixing in the silicon/oxygen interface of the PISN ejecta. In contrast, we do not find a significant dependence of PISN light curves and spectra on viewing angle. These results showcase the capabilities of SuperNu to handle 3D radiation transport and highlight the importance of modeling time series of spectra in identifying PISNe with future missions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Synthesis of nano-forsterite powder by making use of natural silica sand

Nano-forsterite powder with natural silica sand and magnesium powder as the raw materials have been succesfully synthesized. The silica sand was purified followed by a coprecipitation process to obtain colloidal silica. The magnesium powder was dissolved in a chloric acid solution to obtain MgCl{sub 2} solution. The nanoforsterite powder was synthesised using a sol-gel method which included the mixing the colloidal silica and the MgCl{sub 2} solution with various aging and filtering processes. The samples were dried at 100 °C using a hot plate and then the dried powders were calcinated at 900 °C for 2 hours. The samples were characetised for their elements and phase compositions using X-ray Flourescence (XRF) and X-ray Diffraction (XRD) methods, respectively. The diffraction data were qualitatively analyzed using Match!2 software and quantitatively using Rietica software. The crystallite size was verified using Transmission Electron Microscopy (TEM). Results of XRD data analysis showed that the forsterite content reached up to 90.5% wt. The TEM average crystallite size was approximately 53(6) nm.

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

Note: Fully integrated 3.2 Gbps quantum random number generator with real-time extraction

We present a real-time and fully integrated quantum random number generator (QRNG) by measuring laser phase fluctuations. The QRNG scheme based on laser phase fluctuations is featured for its capability of generating ultra-high-speed random numbers. However, the speed bottleneck of a practical QRNG lies on the limited speed of randomness extraction. To close the gap between the fast randomness generation and the slow post-processing, we propose a pipeline extraction algorithm based on Toeplitz matrix hashing and implement it in a high-speed field-programmable gate array. Further, all the QRNG components are integrated into a module, including a compact and actively stabilized interferometer, high-speed data acquisition, and real-time data post-processing and transmission. The final generation rate of the QRNG module with real-time extraction can reach 3.2 Gbps.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Resistive wall tearing mode generated finite net electromagnetic torque in a static plasma

The MARS-F code [Y. Q. Liu et al., Phys. Plasmas 7, 3681 (2000)] is applied to numerically investigate the effect of the plasma pressure on the tearing mode stability as well as the tearing mode-induced electromagnetic torque, in the presence of a resistive wall. The tearing mode with a complex eigenvalue, resulted from the favorable averaged curvature effect [A. H. Glasser et al., Phys. Fluids 18, 875 (1975)], leads to a re-distribution of the electromagnetic torque with multiple peaking in the immediate vicinity of the resistive layer. The multiple peaking is often caused by the sound wave resonances. In the presence of a resistive wall surrounding the plasma, a rotating tearing mode can generate a finite net electromagnetic torque acting on the static plasma column. Meanwhile, an equal but opposite torque is generated in the resistive wall, thus conserving the total momentum of the whole plasma-wall system. The direction of the net torque on the plasma is always opposite to the real frequency of the mode, agreeing with the analytic result by Pustovitov [Nucl. Fusion 47, 1583 (2007)]. When the wall time is close to the oscillating time of the tearing mode, the finite net torque reaches its maximum. Without wall or with an ideal wall, no net torque on the static plasma is generated by the tearing mode. However, re-distribution of the torque density in the resistive layer still occurs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS↗
DOE OSTI2018

Integrated circuits based on conjugated polymer monolayer

It is still a great challenge to fabricate conjugated polymer monolayer field-effect transistors (PoM-FETs) due to intricate crystallization and film formation of conjugated polymers. Here we demonstrate PoM-FETs based on a single monolayer of a conjugated polymer. The resulting PoM-FETs are highly reproducible and exhibit charge carrier mobilities reaching 3 cm 2 V -1 s -1 . The high performance is attributed to the strong interactions of the polymer chains present already in solution leading to pronounced edge-on packing and well-defined microstructure in the monolayer. The high reproducibility enables the integration of discrete unipolar PoM-FETs into inverters and ring oscillators. Real logic functionality has been demonstrated by constructing a 15-bit code generator in which hundreds of self-assembled PoM-FETs are addressed simultaneously. Lastly, our results provide the state-of-the-art example of integrated circuits based on a conjugated polymer monolayer, opening prospective pathways for bottom-up organic electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Fuel Behavior Implications of Reactor Design Choices in Pressurized Water SMRs

Small pressurized water reactors (PWRs) can feature boron free operation, natural circulation mode, reduced height assemblies and/or long refueling cycles. This paper attempts to explore core design optimization for each of these design evolutions. In consequence, five core design layouts are developed incorporating boron free operation with continuous control rods insertion, natural circulation with low burnup/low power density design, natural circulation with high burnup/low power density design, forced circulation with standard core power density design, and forced circulation with high power density design. These cores’ performance is compared to a standard 4-loop PWR. The design process aims to improve the fuel cycle cost under safety constraints through core design optimization using CASMO4E/SIMULATE3 reactor physics codes and FRAPCON4.1 fuel performance assessment tool. Core modeling assumes standard 17x17 PWR fuel assemblies loaded with low enriched uranium (LEU) up to 5wt% or LEU+ (i.e., below 10wt% enrichment) pellets with gadolinium oxide (Gd2O3) as the burnable poison. Satisfactory core and fuel performances are obtained for all the designed cores under steady state and considered overpower transients. For low power density operation, long cycle lengths are achieved reaching a 2.5- and a 5-year cycles and peak rod-average burnup is pushed to 83 MWd/kgU. Other cycle lengths are maintained at 18 months. Boron free operation exhibits the ability to achieve longer cycle lengths at the cost of higher peaking factors leading to high local power and fuel temperatures which prevents sizable power uprates and is deemed uneconomical. Fuel assembly height reduction allows coolant velocity retrofit which enables higher core power density without violating structural integrity of the fuel assembly. As a result, a core power density of 123 kW/l is reached where total cladding hoop strain becomes the limiting parameter.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Impact of a midband gravitational wave experiment on detectability of cosmological stochastic gravitational wave backgrounds

In this report we make forecasts for the impact a future “midband” space-based gravitational wave experiment, most sensitive to 10 -2 - 10Hz,could have on potential detections of cosmological stochastic gravitational wave backgrounds (SGWBs). Specific proposed midband experiments considered are TianGo, B-DECIGO, and AEDGE. We propose a combined power-law integrated sensitivity (CPLS) curve combining GW experiments over different frequency bands, which shows the midband improves sensitivity to SGWBs by up to two orders of magnitude at 10 -2 - 10 Hz . We consider GW emission from cosmic strings and phase transitions as benchmark examples of cosmological SGWBs. We explicitly model various astrophysical SGWB sources, most importantly from unresolved black hole mergers. Using Markov Chain Monte Carlo, we demonstrated that midband experiments can, when combined with LIGO A+ and LISA, significantly improve sensitivities to cosmological SGWBs and better separate them from astrophysical SGWBs. In particular, we forecast that a midband experiment improves sensitivity to cosmic string tension G μ by up to a factor of 10, driven by improved component separation from astrophysical sources. For phase transitions, a midband experiment can detect signals peaking at 0.1–1 Hz, which for our fiducial model corresponds to early Universe temperatures of T * ~10 4 – 10 6 GeV , generally beyond the reach of LIGO and LISA. The midband closes an energy gap and better captures characteristic spectral shape information. It thus substantially improves measurement of the properties of phase transitions at lower energies of T * ~ O(10 3 ) GeV , potentially relevant to new physics at the electroweak scale, whereas in this energy range LISA alone will detect an excess but not effectively measure the phase transition parameters. Our modeling code and chains are publicly available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Nonlinear excitation of low-n harmonics in reduced magnetohydrodynamic simulations of edge-localized modes

Nonlinear simulations of the early edge-localized mode (ELM) phase based on a typical type-I ELMy ASDEX Upgrade discharge have been carried out using the reduced MHD code JOREK. The analysis is focused on the evolution of the toroidal Fourier spectrum. It is found that during the nonlinear evolution, linearly subdominant low-n Fourier components, in particular the n = 1, grow to energies comparable with linearly dominant harmonics. A simple model is developed, based on the idea that energy is transferred among the toroidal harmonics via second order nonlinear interaction. The simple model reproduces and explains very well the early nonlinear evolution of the toroidal spectrum in the JOREK simulations. Furthermore, it is shown for the n = 1 harmonic, that its spatial structure changes significantly during the transition from linear to nonlinearly driven growth. The rigidly growing structure of the linearly barely unstable n = 1 reaches far into the plasma core. In contrast, the nonlinearly driven n= 1 has a rigidly growing structure localized at the plasma edge, where the dominant toroidal harmonics driving the n = 1 are maximal and in phase. The presented quadratic coupling model might explain the recent experimental observation of strong low-n components in magnetic measurements [Wenninger et al., “Non-linear magnetic perturbations during edge localized modes in TCV dominated by low n mode components,” Nucl. Fusion (submitted)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗