Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Continuous symmetries”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Self Induced Buoyant Blow Off in Upward Flame Spread on Thin Solid Fuels

Upward flame spread experiments were conducted on a thin fabric cloth consisting of 75% cotton and 25% fiberglass. The sample is sandwiched symmetrically with stainless steel plates with the exposed width varying between 2 to 8.8 cm from test to test and >1.5m tall. The bottom edge was ignited resulting in a symmetric two sided flame. For the narrower samples (. 5cm), two sided flame growth would proceed until reaching some limiting value (15-30 cm depending on sample width). Fluctuation or instability of the flame base on one side would initially become visible and then the flame base would retreat downstream and cause extinguishment on one side. Detailed examination of the still images shows that the fuel continues to vaporize from the extinguished side due to the thermally thin nature of the fuel. But, due to the remaining inert fiberglass mesh, which acts as a flashback arrestor, the extinguished side was not able to be reignited by the remaining flame. The remaining flame would then shrink in length due to the reduced heat transfer to the solid to a shorter length. The one-sided flame will spread stably with a constant speed and a constant flame length to the end of the sample. A constant length flame implies that the pyrolysis front and the burnt out fronts move at the same speed. For the wider samples (. 7cm), no one-sided extinction is observed. Two-sided flames spread all the way to the top of the sample. For these wider widths, the flames are still growing and have not reached their limiting length if it exists. Care was taken to minimize the amount of non-symmetries in the experimental configuration. Repeated tests show that blow-off can occur on either side of the sample. The flame growth is observed to be very symmetric during the growth phase and grew to significant length (>10cm) before extinction of the flame on one side. Our proposed explanation of this unusual phenomenon (i.e. stronger two ]sided flame cannot exist but weaker one-sided flame can) is as follows: The observed one-sided extinction is a blow- off induced by buoyant entrainment. It is known that the flammable diffusion flame regime is bounded by quenching and blow ]off limits when varying incoming air velocity. The narrowest samples tested (between 2 and 5 cm) begin within the flammable range, but as the flame grows, the buoyancy driven air velocity increases at the neighborhood of the flame base. The initially stable flame crosses the extinguishment boundary resulting in a flame blow-off. When one-side of the flame extinguishes, the remaining side shrinks due to the reduced heat transfer to the solid. This reduces the induced velocity and the flame becomes stable. It is proposed that this may have implications to upward flame growth beyond this experiment.

Johnston, Michael C.↗
DOE OSTI2023

Status of CHIP-TRAP: The Central Michigan University High-Precision Penning Trap

Precise and accurate atomic mass data provide crucial information for applications in a wide range of fields in physics and beyond, including astrophysics, nuclear structure, particle and neutrino physics, fundamental symmetries, chemistry, and metrology. The most precise atomic mass measurements are performed on charged particles confined in a Penning trap. Here, we describe the development, status, and outlook of CHIP-TRAP: the Central Michigan University high-precision Penning trap. CHIP-TRAP aims to perform ultra-high precision (~1 part in 10 11 fractional precision) mass measurements on stable and long-lived isotopes produced with external ion sources and transported to the Penning traps. Along the way, ions of a particular m/q are selected with a multi-reflection time-of-flight mass separator (MR-TOF-MS), with further filtering performed in a cylindrical capture trap before the ions are transported to a pair of hyperbolic measurement traps. In this paper, we report on the design and status of CHIP-TRAP and present results from the commissioning of the ion sources, MR-TOF-MS, and capture trap. We also provide an outlook on the continued development and commissioning of CHIP-TRAP.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Fluctuating hydrodynamics of chiral active fluids

Active materials are characterized by continuous injection of energy at the microscopic level and typically cannot be adequately described by equilibrium thermodynamics. Here we study a class of active fluids in which equilibrium-like properties emerge when fluctuating and activated degrees of freedom are statistically decoupled, such that their mutual information is negligible. We analyse three paradigmatic systems: chiral active fluids composed of spinning frictional particles that are free to translate, oscillating granular gases and active Brownian rollers. In all of these systems, a single effective temperature generated by activity parameterizes both the equation of state and the emergent Boltzmann statistics. The same effective temperature, renormalized by velocity correlations, relates viscosities to steady-state stress fluctuations via a Green–Kubo relation. To rationalize these observations, we develop a theory for the fluctuating hydrodynamics of these non-equilibrium fluids and validate it through large-scale molecular dynamics simulations. Our work sheds light on the microscopic origin of odd viscosities and stress fluctuations characteristic of parity-violating fluids, in which mirror symmetry and detailed balance are broken.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2005

Video-Camera-Based Position-Measuring System

A prototype optoelectronic system measures the three-dimensional relative coordinates of objects of interest or of targets affixed to objects of interest in a workspace. The system includes a charge-coupled-device video camera mounted in a known position and orientation in the workspace, a frame grabber, and a personal computer running image-data-processing software. Relative to conventional optical surveying equipment, this system can be built and operated at much lower cost; however, it is less accurate. It is also much easier to operate than are conventional instrumentation systems. In addition, there is no need to establish a coordinate system through cooperative action by a team of surveyors. The system operates in real time at around 30 frames per second (limited mostly by the frame rate of the camera). It continuously tracks targets as long as they remain in the field of the camera. In this respect, it emulates more expensive, elaborate laser tracking equipment that costs of the order of 100 times as much. Unlike laser tracking equipment, this system does not pose a hazard of laser exposure. Images acquired by the camera are digitized and processed to extract all valid targets in the field of view. The three-dimensional coordinates (x, y, and z) of each target are computed from the pixel coordinates of the targets in the images to accuracy of the order of millimeters over distances of the orders of meters. The system was originally intended specifically for real-time position measurement of payload transfers from payload canisters into the payload bay of the Space Shuttle Orbiters (see Figure 1). The system may be easily adapted to other applications that involve similar coordinate-measuring requirements. Examples of such applications include manufacturing, construction, preliminary approximate land surveying, and aerial surveying. For some applications with rectangular symmetry, it is feasible and desirable to attach a target composed of black and white squares to an object of interest (see Figure 2). For other situations, where circular symmetry is more desirable, circular targets also can be created. Such a target can readily be generated and modified by use of commercially available software and printed by use of a standard office printer. All three relative coordinates (x, y, and z) of each target can be determined by processing the video image of the target. Because of the unique design of corresponding image-processing filters and targets, the vision-based position- measurement system is extremely robust and tolerant of widely varying fields of view, lighting conditions, and varying background imagery.

Lane, John↗
DOE OSTI2020

A micromirror array with annular partitioning for high-speed random-access axial focusing

Dynamic axial focusing functionality has recently experienced widespread incorporation in microscopy, augmented/virtual reality (AR/VR), adaptive optics and material processing. However, the limitations of existing varifocal tools continue to beset the performance capabilities and operating overhead of the optical systems that mobilize such functionality. The varifocal tools that are the least burdensome to operate (e.g. liquid crystal, elastomeric or optofluidic lenses) suffer from low (≈100 Hz) refresh rates. Conversely, the fastest devices sacrifice either critical capabilities such as their dwelling capacity (e.g. acoustic gradient lenses or monolithic micromechanical mirrors) or low operating overhead (e.g. deformable mirrors). Here, we present a general-purpose random-access axial focusing device that bridges these previously conflicting features of high speed, dwelling capacity and lightweight drive by employing low-rigidity micromirrors that exploit the robustness of defocusing phase profiles. Geometrically, the device consists of an 8.2 mm diameter array of piston-motion and 48-μm-pitch micromirror pixels that provide 2π phase shifting for wavelengths shorter than 1100 nm with 10–90% settling in 64.8 μs (i.e., 15.44 kHz refresh rate). The pixels are electrically partitioned into 32 rings for a driving scheme that enables phase-wrapped operation with circular symmetry and requires <30 V per channel. Optical experiments demonstrated the array’s wide focusing range with a measured ability to target 29 distinct resolvable depth planes. Overall, the features of the proposed array offer the potential for compact, straightforward methods of tackling bottlenecked applications, including high-throughput single-cell targeting in neurobiology and the delivery of dense 3D visual information in AR/VR.

42 ENGINEERING↗
NASA NTRS1999

Surface Tension Driven Convection Experiment-2 (STDCE-2)

Thermocapillary flows are known to become oscillatory (time-periodic), but how and when they become oscillatory in containers of unit-order aspect ratio are not yet fully understood. The present work is a part of our continuous effort to obtain a better understanding of the phenomenon. Thermocapillary flow experiments in normal gravity are limited to a narrow parametric range in order to minimize gravity and buoyancy effects, which is an important reason for our lack of full understanding of the oscillation phenomenon. One important unanswered question is what role, if any, free surface deformation plays in the oscillation mechanism. For that reason we performed thermocapillary flow experiments, called the Surface Tension Driven Convection Experiment-2 (STDCE-2), aboard the USML-2 Spacelab in 1995. The main objectives of the experiments were to investigate oscillatory thermocapillary flows in microgravity and to clarify the importance of free surface deformation in such flows. Steady and oscillatory thermocapillary flows were generated in cylindrical containers by employing two heating modes. A CO2 laser with adjustable power and beam diameter was used in the Constant Flux (CF) configuration to heat the free surface. The other configuration investigated in STDCE-2 was the Constant Temperature (CT) configuration in which a submerged cylindrical cartridge heater placed at the symmetry (axial) axis of the test container heated the fluid. Both heating modes cause non-uniform temperature distributions on the free surface, which generates thermocapillary flow. The flow field was investigated by flow visualization, and the temperature field was measured by thermistors and an infrared imager. The free surface shape and motion were measured by a Ronchi system. The hardware performed well and we were able to conduct more tests than originally planned. From the successful experiments a large amount of data was acquired. The analysis of the data is now nearly complete. Some important results are presented and discussed herein.

Masud, J.↗
DOE OSTI2021

Near-infrared C-term MCD spectroscopy of octahedral uranium (V) complexes

C-term magnetic circular dichroism (MCD) spectroscopy is a powerful method for probing d–d and f–f transitions in paramagnetic metal complexes. However, this technique remains underdeveloped both experimentally and theoretically for studies of U(V) complexes of Oh symmetry, which have been of longstanding interest for probing electronic structure, bonding, and covalency in 5f systems. In this study, C-term NIR MCD of the Laporte forbidden f–f transitions of [UCl 6 ] - and [UF 6 ] - are reported, demonstrating the significant fine structure resolution possible with this technique including for the low energy Γ 7 → Γ 8 transitions in [UF 6 ] - . The experimental NIR MCD studies were further extended to [U(OC 6 F 5 ) 6 ] - , [U(CH 2 SiMe 3 ) 6 ] - , and [U(NC(tBu)(Ph)) 6 ] - to evaluate the effects of ligand-type on the f–f MCD fine structure features. Theoretical calculations were conducted to determine the Laporte forbidden f–f transitions and their MCD intensity experimentally observed in the NIR spectra of the U(V) hexahalide complexes, via the inclusion of vibronic coupling, to better understand the underlying spectral fine structure features for these complexes. These spectra and simulations provide an important platform for the application of MCD spectroscopy to this widely studied class of U(V) complexes and identify areas for continued theoretical development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2014

A Multiscale Progressive Failure Modeling Methodology for Composites that Includes Fiber Strength Stochastics

A multiscale modeling methodology was developed for continuous fiber composites that incorporates a statistical distribution of fiber strengths into coupled multiscale micromechanics/finite element (FE) analyses. A modified two-parameter Weibull cumulative distribution function, which accounts for the effect of fiber length on the probability of failure, was used to characterize the statistical distribution of fiber strengths. A parametric study using the NASA Micromechanics Analysis Code with the Generalized Method of Cells (MAC/GMC) was performed to assess the effect of variable fiber strengths on local composite failure within a repeating unit cell (RUC) and subsequent global failure. The NASA code FEAMAC and the ABAQUS finite element solver were used to analyze the progressive failure of a unidirectional SCS-6/TIMETAL 21S metal matrix composite tensile dogbone specimen at 650 degC. Multiscale progressive failure analyses were performed to quantify the effect of spatially varying fiber strengths on the RUC-averaged and global stress-strain responses and failure. The ultimate composite strengths and distribution of failure locations (predominately within the gage section) reasonably matched the experimentally observed failure behavior. The predicted composite failure behavior suggests that use of macroscale models that exploit global geometric symmetries are inappropriate for cases where the actual distribution of local fiber strengths displays no such symmetries. This issue has not received much attention in the literature. Moreover, the model discretization at a specific length scale can have a profound effect on the computational costs associated with multiscale simulations.models that yield accurate yet tractable results.

Generalized Method of Cells↗
DOE OSTI2026

Constraining hydrodynamic models of inertial confinement fusion implosions using capsule surrogate experiments

We conduct capsule surrogate experiments at the National Ignition Facility to calibrate radiation hydrodynamic simulations to infer hydrodynamic conditions that are not observable in indirect drive ignition implosions. We tune the simulations by applying laser power and cross beam energy transfer (CBET) saturation multipliers to match the observables from capsule surrogate experiments. Shock timing, velocity, and symmetry are measured in liquid D 2 filled Keyhole capsule surrogate experiments and implosion trajectory, stagnation time, and shape time history are measured in in-flight 2D backlit x-ray radiography experiments (“2DConA”) of D 2 gas filled capsule implosions. Calibrated simulations suggest that the N210808 ignition implosion (fusion target gain = 0.7) had a shell mass remaining at stagnation of less than the nominal %5 (3.8%) and resulted in less confinement. For N221204, the shell was made 5.75 μm thicker to trade implosion velocity for increased confinement and resulted in a target gain = 1.5 with a shell mass remaining of 5.7%. Furthermore, a single adjusted model can reproduce all shock timing data as changes are made to shell thickness (79–85 μm) and laser wavelength separation (1.8–4.0 Å). However, for the 2DConA implosions, a 5% variation in the peak power laser multipliers and a 30% variation in late-time CBET between experiments are needed to match the observed stagnation times, in-flight $P_2$ shape, and hot-spot $P_2$ shape. While progress is being made to improve the models in simulations using focused experiments, capsule surrogate experiments will continue to be needed to optimize future ignition designs.

Lasers↗
DOE OSTI2025

Monomer-dimer tensor-network basis for qubit-regularized lattice gauge theories

Traditional SU⁡(𝑁) lattice gauge theories (LGTs) can be formulated using an orthonormal basis constructed from the irreducible representations (irreps) 𝑉 𝜆 of the SU⁡(𝑁) gauge symmetry. On a lattice, the elements of this basis are tensor networks comprising dimer tensors on the links labeled by a set of irreps {𝜆 ℓ } and monomer tensors on sites labeled by {𝜆 𝑠 }. These tensors naturally define a local site Hilbert space, ℋ$^𝑔_𝑠$, on which gauge transformations act. Gauss’s law introduces an additional index 𝛼 𝑠 =1,2,…,𝒟⁡(ℋ$^𝑔_𝑠$) that labels an orthonormal basis of the gauge-invariant subspace of ℋ$^𝑔_𝑠$. This monomer-dimer tensor-network (MDTN) basis, |{𝜆 𝑠 },{𝜆 ℓ },{𝛼 𝑠 }⟩, of the physical Hilbert space enables the construction of new qubit-regularized SU⁡(𝑁) gauge theories that are free of sign problems while preserving key features of traditional LGTs. Here, we investigate finite-temperature confinement-deconfinement transitions in a simple qubit-regularized SU(2) and SU(3) gauge theory in 𝑑 =2 and 𝑑 =3 spatial dimensions, formulated using the MDTN basis, and show that they reproduce the universal results of traditional LGTs at these transitions. Additionally, in 𝑑 =1, we demonstrate using a plaquette chain that the string tension at zero temperature can be continuously tuned to zero by adjusting a model parameter that plays the role of the gauge coupling in traditional LGTs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2016

Spiral Lamber's Problem with Generalized Logarithmic Spirals

Lambert’s problem subject to a continuous acceleration is solved using the family of generalized logarithmic spirals. Thanks to the existence of two first integrals related to the energy and angular momentum surprising analogies with the Keplerian case are found. A minimum-energy spiral transfer exists. Increasing the value of the constant of the generalized energy yields pairs of conjugate spiral trajectories. The properties of such spirals are strongly connected with the properties of conjugate Keplerian orbits. When the generalized constant of the energy reaches a critical value the two solutions degenerate into a pair of parabolic spirals, one of which connects the two vectors through infinity. From that point the spiral transfers become hyperbolic. Generalized logarithmic spirals admit closed-form solutions to all the required magnitudes including the time of flight, providing a deep insight into the dynamics of the problem. In addition, the maximum acceleration along the transfer is found analytically so the solutions that violate the design constraints on the maximum thrust acceleration can be rejected without any further computations. When the time of flight is fixed there is still a degree of freedom in the solution, related to a control parameter. Resonant transfers appear naturally thanks to the symmetry properties of the generalized logarithmic spirals. The problem of designing a low-thrust transfer between two bodies can be reduced to solving the corresponding spiral Lambert’s problem. In order to show the versatility of the method it is applied to the design of an asteroid tour and to explore launch opportunities to Mars.\

Roa, Javier↗
DOE OSTI2021

A Machine Learning Approach for the Prediction of Formability and Thermodynamic Stability of Single and Double Perovskite Oxides

Perovskite oxides continue to attract huge interest due to their fascinating and wide-ranging properties for diverse applications. The tunability of these properties may be further enhanced by increasing their compositional complexity via double perovskite-ordered configurations containing multiple cations. In this work, we focus on an exhaustive chemical space of single and double oxide perovskites and optimally explore this space to identify novel compositions that are likely to form stable compounds. Critically, we examine the relationship between formability, the practical ability to synthesize a compound, and stability, the thermodynamic preference to form the structure. Our formability and stability training data sets were enumerated from the available experimental literature and in-house density functional theory computations and contained 1505 and 3469 examples, respectively, representing state-of-the-art in the current open literature in perovskite and double perovskite compounds. Subsequently, cross-validated and highly accurate machine learning classification models are built using these training data sets and employed to screen for novel stable oxide perovskites. The study identifies (1) atomic features relevant to prediction of formability and stability in perovskite and double perovskite compounds, (2) the importance of including energy contributions due to local structural relaxations going beyond the high symmetry perovskite phase, and (3) 437,828 double perovskite compounds that are likely to be stable and 891,188 compounds that are likely to be formable. From the intersection of this large chemical space of formable and stable oxide perovskites, 414 compositions are identified as the most promising candidates for future experimental synthesis of novel oxide perovskites. The developed models may be generalized and have implications beyond perovskite discovery if applied to other families of compounds.

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

Analyticity of replica correlators and modular ETH

We study the two point correlation function of a local operator on an n -sheeted replica manifold corresponding to the half-space in the vacuum state of a conformal field theory. In analogy with the inverse Laplace transform, we define the Renyi transform of this correlation function, which is a function of one complex variable w, dual to the Renyi parameter n . Inspired by the inversion formula of Caron-Huot, we argue that if the Renyi transform ƒ(w) has bounded behavior at infinity in the complex w plane, the discontinuity of the Renyi transform disc ƒ(w) provides the unique analytic continuation in n of the original replica correlation function. We check our formula by explicitly calculating the Renyi transform of a particular replica correlator in a large N holographic CFT d in dimensions d > 2. We also discover that the discontinuity of the Renyi transform is related to the matrix element of local operators between two distinct eigenstates of the modular Hamiltonian. We calculate the Renyi transform in 2 d conformal field theories, and use it to extract the off-diagonal elements of (modular) ETH. We argue that in 2 d , this is equivalent to the off-diagonal OPE coefficients of a CFT and show that our technique exactly reproduces recent results in the literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2010

The OMPS Limb Profiler Instrument: Two-Dimensional Retrieval Algorithm

The upcoming Ozone Mapper and Profiler Suite (OMPS), which will be launched on the NPOESS Preparatory Project (NPP) platform in early 2011, will continue monitoring the global distribution of the Earth's middle atmosphere ozone and aerosol. OMPS is composed of three instruments, namely the Total Column Mapper (heritage: TOMS, OMI), the Nadir Profiler (heritage: SBUV) and the Limb Profiler (heritage: SOLSE/LORE, OSIRIS, SCIAMACHY, SAGE III). The ultimate goal of the mission is to better understand and quantify the rate of stratospheric ozone recovery. The focus of the paper will be on the Limb Profiler (LP) instrument. The LP instrument will measure the Earth's limb radiance (which is due to the scattering of solar photons by air molecules, aerosol and Earth surface) in the ultra-violet (UV), visible and near infrared, from 285 to 1000 nm. The LP simultaneously images the whole vertical extent of the Earth's limb through three vertical slits, each covering a vertical tangent height range of 100 km and each horizontally spaced by 250 km in the cross-track direction. Measurements are made every 19 seconds along the orbit track, which corresponds to a distance of about 150km. Several data analysis tools are presently being constructed and tested to retrieve ozone and aerosol vertical distribution from limb radiance measurements. The primary NASA algorithm is based on earlier algorithms developed for the SOLSE/LORE and SAGE III limb scatter missions. All the existing retrieval algorithms rely on a spherical symmetry assumption for the atmosphere structure. While this assumption is reasonable in most of the stratosphere, it is no longer valid in regions of prime scientific interest, such as polar vortex and UTLS regions. The paper will describe a two-dimensional retrieval algorithm whereby the ozone distribution is simultaneously retrieved vertically and horizontally for a whole orbit. The retrieval code relies on (1) a forward 2D Radiative Transfer code (to model limb radiances within a non-uniform atmosphere and evaluate 2D analytical partial derivatives) and (2) an optimal estimator inversion routine. The algorithm uses the typically sparse nature of the kernel matrices as well as fast matrix inversion techniques to allow for fast inversion of limb data with efficient memory management (as was done for MIPAS data processing). While the method has so far only been developed in the context of Single Scatter, the paper will show how the CPU intensive Multiple Scatter modeling can be implemented using parallel CPU processing. Initial results will be presented in terms of retrieved ozone profiles and code performance.

Rault, Didier F.↗
DOE OSTI2024

Graviton topology

Over the past three decades, it has been shown that discrete and continuous media can support topologically nontrivial waves. Recently, it was shown that the same is true of the vacuum, in particular, right (R) and left (L) circularly polarized photons are topologically nontrivial. Here, we study the topology of another class of massless particles, namely gravitons. We show that the collection of all gravitons forms a topologically trivial vector bundle over the lightcone, allowing us to construct a globally smooth basis for gravitons. The graviton bundle also has a natural geometric splitting into two topologically nontrivial subbundles, consisting of the R and L gravitons. The R and L gravitons are unitary irreducible bundle representations of the Poincaré group, and are thus elementary particles; their topology is characterized by the Chern numbers ∓4. This nontrivial topology obstructs the splitting of graviton angular momentum into spin and orbital angular momentum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2008

Configuration Effects on Acoustic Performance of a Duct Liner

Continued success in aircraft engine noise reduction necessitates ever more complete understanding of the effect that flow path geometry has on sound propagation in the engine. The Curved Duct Test Rig (CDTR) has been developed at NASA Langley Research Center to investigate sound propagation through a duct of comparable size (approximately the gap of GE90) and physical characteristics to the aft bypass duct of typical aircraft engines. The liner test section is designed to mimic the outer/inner walls of an engine exhaust bypass duct that has been unrolled circumferentially. Experiments to investigate the effect of curvature along the flow path on the acoustic performance of a test liner are performed in the CDTR and reported in this paper. Flow paths investigated include both straight and curved with offsets from the inlet to the discharge plane of and 1 duct width, respectively. The test liners are installed on the side walls of the liner test section. The liner samples are perforate over honeycomb core, which design is typical of liners installed in aircraft nacelles. In addition to fully treated side walls, combinations of treated and acoustically rigid walls are investigated. While curvature in the hard wall duct is found not to reduce the incident sound significantly, it does cause mode scattering. It is found that asymmetry of liner treatment causes scattering of the incident mode into less attenuated modes, which degrades the overall liner attenuation. It is also found that symmetry of liner treatment enhances liner performance by eliminating scattering into less attenuated modes. Comparisons of measured liner attenuation with numerical results predicted by an analytic model based on the parabolic approximation (CDUCT-LaRC) have also been made and are reported in this paper. The effect of curvature in the rigid wall configuration estimated by CDUCT-LaRC is similar to the observed results, and the mode scattering seen in the measurements also occurs in the analytic model results. The analytic model and experiment show similar differences of overall attenuation between one wall treated and both walls treated.

Gerhold, Carl H.↗
DOE OSTI2024

Proton and neutron density distributions at supranormal density in low- and medium-energy heavy-ion collisions. II. Central Pb + Pb collisions

This paper represents a continuation of our investigation into the limitations on total, proton, and neutron particle number densities, as well as the asymmetry of proton and neutron density distributions achievable in central heavy-ion collisions. We explore these aspects at low and medium energies within the framework of the Boltzmann-Uhlenbeck-Uehling (pBUU) transport and time-dependent Hartree-Fock (TDHF) models. Previous studies, focusing on symmetric and asymmetric collisions of Ca and Sn nuclei [1], and initial results on Pb-nuclei collisions [2], emphasized the role of the Coulomb interaction in these events. Our findings indicated that: (i) the highest total densities predicted at 𝐸 beam = 800 MeV/nucleon were on the order of ≈ 2.5⁢𝜌 0 (𝜌 0 = 0.16 fm −3 ), (ii) the proton-neutron asymmetry for maximal densities, 𝛿=(𝜌$^{max}_{n}$− 𝜌$^{max}_{p}$)/(𝜌$^{max}_{n}$+𝜌$^{max}_{p}$), did not generally exceed the asymmetry in the initial state of the collision at all beam energies and tended to decrease during the reaction, and (iii) a significant portion of this asymmetry had its microscopic origin in Coulomb forces, masking the pure nuclear contribution. These new findings, particularly relevant in the astrophysical context, are further examined in this work, focusing on the heaviest target-projectile combination 212,208 Pb accessible in an experiment. We introduce the SkT3 Skyrme force model, not previously used for the Pb system, and compare it to the SV-bas and SV-sym34 models to explore the symmetry-energy dependence of the results. Contour plots of nucleonic densities are presented, contrasting the time evolution of the density distributions in low (TDHF) and high (pBUU) models. We also present the evolution of normalized maximal proton, neutron, and total nucleon number density with increasing beam energy in the full pBUU model and the Vlasov approximation, aiming to explore the impact of correlations in the reaction. The time evolution of the proton and neutron density distributions in the plane transverse to the beam direction is illustrated at both low and high beam energy. In conclusion, our detailed examination of the Pb system in this work provides further essential evidence that the aforementioned findings (i)–(iii) are only weakly dependent on system size and a symmetry-energy model, and thus, they are of more general importance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Non-linear MHD modelling of edge localized modes suppression by resonant magnetic perturbations in ITER

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗