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 559 records · Page 31

Hyper-fidelity depletion coupled with discrete pebble motion in pebble bed reactors

Pebble bed reactors have raised new interest during the past decade due to their attractive characteristics. Therefore, accurate simulations must be performed to better understand these systems and ensure optimal and safe designs. Most current methods use lower fidelity approaches with representative unit-cells or macro-zones with uniform fluxes, which have accuracy and flexibility limitations. A novel hyper-fidelity method for pebble bed reactors depletion is presented and internally couples Serpent 2 and a pseudo-motion routine. Pseudo-motion is applied handling vertical shifts of compositions in a static pebble bed, random reinsertion of used pebbles, insertion of fresh pebbles and used pebbles discarding. Associated with individual depletion to correctly determine the flux spectrum and composition in each pebble, this hyper-fidelity approach paves the way towards more accurate depletion calculation in pebble bed reactors. Using this method, a demonstration is completed on a small-scale reactor. In this application, the core reaches equilibrium, and the following data is extracted: core-wise parameters evolution, pebble-wise spatial and statistical distribution. Discarded pebbles are analyzed, and relevant information is shown. This work proves the feasibility of hyper fidelity depletion with Serpent 2, and the range of use for this method: reactor design and analysis for equilibrium and slow transients, lower fidelity methods validation and feeding fuel performance, thermal-hydraulics, or waste management models. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Effect of Unsaturated Flow Modes on Partitioning Dynamics of Gravity-Driven Flow at a Simple Fracture Intersection: Laboratory Study and Three-Dimensional Smoothed Particle Hydrodynamics Simulations

In this work, we study unsaturated gravity-driven free-surface flows on a synthetic surface intersected by a horizontal fracture and the importance of flow modes (droplet; rivulet) on the partitioning behavior at the fracture intersection. We present laboratory experiments, three-dimensional smoothed particle hydrodynamics (SPH) simulations using a heavily parallelized code, and an analytical solution. The flow-rate-dependent mode switching from droplets to rivulets is reproduced by the SPH model, and the transition ranges agree with the laboratory experiments. We show that flow modes heavily influence the bypass behavior. Flows favoring the formation of droplets exhibit a much stronger bypass capacity compared to rivulet flows, where nearly the whole fluid mass is initially stored within the horizontal fracture. This behavior is demonstrated for a multi-inlet laboratory setup where the inlet- specific flow rate is chosen so that either a droplet or rivulet flow persists. The effect of fluid buffering within the horizontal fracture is presented in terms of dimensionless fracture in- flow so that characteristic scaling regimes can be recovered. For both cases (rivulets and droplets), the flow within the horizontal fracture transitions into a Washburn regime until a critical threshold is reached and the bypass efficiency increases. For rivulet flows, the initial filling of the horizontal fracture is described by classical plug flow. Meanwhile, for droplet flows, a size-dependent partitioning behavior is observed, and the filling of the fracture takes longer. For the case of rivulet flow, we provide an analytical solution that demonstrates the existence of classical Washburn flow within the horizontal fracture.

unsaturated flow, fractured media, smoothed partic↗

Implementation and Synthesis of Math Library Functions

Achieving speed and accuracy for math library functions like exp, sin, and log is difficult. This is because low-level implementation languages like C do not help math library developers catch mathematical errors, build implementations incrementally, or separate high-level and low-level decision making. This ultimately puts development of such functions out of reach for all but the most experienced experts. To address this, we introduce MegaLibm, a domain-specific language for implementing, testing, and tuning math library implementations. MegaLibm is safe, modular, and tunable. Implementations in MegaLibm can automatically detect mathematical mistakes like sign flips via semantic wellformedness checks, and components like range reductions can be implemented in a modular, composable way, simplifying implementations. Once the high-level algorithm is done, tuning parameters like working precisions and evaluation schemes can be adjusted through orthogonal tuning parameters to achieve the desired speed and accuracy. MegaLibm also enables math library developers to work interactively, compiling, testing, and tuning their implementations and invoking tools like Sollya and type-directed synthesis to complete components and synthesize entire implementations. MegaLibm can express 8 state-of-the-art math library implementations with comparable speed and accuracy to the original C code, and can synthesize 5 variations and 3 from-scratch implementations with minimal guidance.

97 MATHEMATICS AND COMPUTING↗

Nonlinear modeling of ELM mitigation with RMP on HL-2A

Abstract Nonlinear modeling of mitigation of the edge localized mode (ELM) with resonant magnetic perturbation (RMP) is performed for the HL-2A tokamak, utilizing the three-dimensional (3D) magnetohydrodynamic code JOREK. Based on the 3D equilibrium established after application of the n = 1 ( n is the toroidal mode number) RMP at 4.9 kAt coil current with odd parity, ELM mitigation is successfully simulated consistent with the experimental result. Nonlinear simulations show strong mode coupling among toroidal Fourier harmonics, allowing redistribution of the magnetic energy such that the most unstable toroidal mode saturates at a lower level. This magnetic energy cascade offers an explanation of the RMP-induced ELM mitigation achieved in HL-2A. Detailed examination of the simulation results shows persistent resonant field screening even during the ELM mitigation phase. Finite plasma resistivity however does enable partial penetration of the resonant field thus modifying the edge magnetic topology and characteristics of the edge transport. Plasma radial profiles undergo pronounced changes around the pedestal region, when the magnetic energy of the most unstable toroidal mode reaches the maximum value. Systematic scans of the applied RMP coil current with the JOREK simulations find a threshold value of around 4.5 kAt required for achieving the ELM mitigation on HL-2A.

Physics↗

Ab initio Bogoliubov coupled cluster theory for open-shell nuclei

Background: Ab initio many-body methods have been developed over the past 10 yr to address closed-shell nuclei up to mass A≈130 on the basis of realistic two- and three-nucleon interactions. A current frontier relates to the extension of those many-body methods to the description of open-shell nuclei. Several routes to address open-shell nuclei are currently under investigation, including ideas that exploit spontaneous symmetry breaking. Purpose: Singly open-shell nuclei can be efficiently described via the sole breaking of U(1) gauge symmetry associated with particle-number conservation as a way to account for their superfluid character. While this route was recently followed within the framework of self-consistent Green's function theory, the goal of the present work is to formulate a similar extension within the framework of coupled cluster theory. Methods: We formulate and apply Bogoliubov coupled cluster (BCC) theory, which consists of representing the exact ground-state wave function of the system as the exponential of a quasiparticle excitation cluster operator acting on a Bogoliubov reference state. Equations for the ground-state energy and the cluster amplitudes are derived at the singles and doubles level (BCCSD) both algebraically and diagrammatically. The formalism includes three-nucleon forces at the normal-ordered two-body level. The first BCC code is implemented in m scheme, which will permit the treatment of doubly open-shell nuclei via the further breaking of SU(2) symmetry associated with angular momentum conservation. Results: Proof-of-principle calculations in an N max =6 spherical harmonic oscillator basis for 16,18 O and 18 Ne in the BCCD approximation are in good agreement with standard coupled cluster results with the same chiral two-nucleon interaction, while 20 O and 20 Mg display underbinding relative to experiment. The breaking of U(1) symmetry, monitored by computing the variance associated with the particle-number operator, is relatively constant for all five nuclei, in both the Hartree-Fock-Bogoliubov and BCCD approximations. Conclusions: The newly developed many-body formalism increases the potential span of ab initio calculations based on single-reference coupled cluster techniques tremendously, i.e., potentially to reach several hundred additional midmass nuclei. The new formalism offers a wealth of potential applications and further extensions dedicated to the description of ground and excited states of open-shell nuclei. Short-term goals include the implementation of three-nucleon forces at the normal-ordered two-body level. Midterm extensions include the approximate treatment of triples corrections and the development of the equation-of-motion methodology to treat both excited states and odd nuclei. Long-term extensions include exact restoration of U(1) and SU(2) symmetries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Gravitational wave signatures of ab initio two-dimensional core collapse supernova explosion models for 12 – 25 M ⊙ stars

Here, we present the gravitational waveforms computed in ab initio two-dimensional core collapse supernova models evolved with the chimera code for progenitor masses between 12 and 25 M. For all models employ multifrequency neutrino transport in the ray-by-ray approximation, state-of-the-art weak interaction physics, relativistic transport corrections such as the gravitational redshift of neutrinos, two-dimensional hydrodynamics with the commensurate relativistic corrections, Newtonian self-gravity with a general-relativistic monopole correction, and the Lattimer-Swesty equation of state with 220 MeV compressibility, and begin with the most recent Woosley-Heger nonrotating progenitors in this mass range. All of our models exhibit robust explosions. Moreover, our waveforms capture all stages of supernova development: 1) a relatively short and weak prompt signal, 2) a quiescent stage, 3) a strong signal due to convection and standing accretion shock instability activity, 4) termination of active accretion onto the proto-neutron star, and 5) a slowly increasing tail that reaches a saturation value. Fourier decomposition shows that the gravitational wave signals we predict should be observable by AdvLIGO for Galactic events across the range of progenitors considered here. The fundamental limitation of these models is in their imposition of axisymmetry. Further progress will require counterpart three-dimensional models, which are underway.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of the nonlinear behavior of shear-Alfvén modes in tokamaks based on Hamiltonian mapping techniques

We present a series of numerical simulation experiments set up to illustrate the fundamental physics processes underlying the nonlinear dynamics of Alfvénic modes resonantly excited by energetic particles in tokamak plasmas and of the ensuing energetic particle transports. These phenomena are investigated by following the evolution of a test particle population in the electromagnetic fields computed in self-consistent MHD-particle simulation performed by the HMGC code. Hamiltonian mapping techniques are used to extract and illustrate several features of wave-particle dynamics. The universal structure of resonant particle phase space near an isolated resonance is recovered and analyzed, showing that bounded orbits and untrapped trajectories, divided by the instantaneous separatrix, form phase space zonal structures, whose characteristic non-adiabatic evolution time is the same as the nonlinear time of the underlying fluctuations. Bounded orbits correspond to a net outward resonant particle flux, which produces a flattening and/or gradient inversion of the fast ion density profile around the peak of the linear wave-particle resonance. The connection of this phenomenon to the mode saturation is analyzed with reference to two different cases: a Toroidal Alfvén eigenmode in a low shear magnetic equilibrium and a weakly unstable energetic particle mode for stronger magnetic shear. It is shown that, in the former case, saturation is reached because of radial decoupling (resonant particle redistribution matching the mode radial width) and is characterized by a weak dependence of the mode amplitude on the growth rate. In the latter case, saturation is due to resonance detuning (resonant particle redistribution matching the resonance width) with a stronger dependence of the mode amplitude on the growth rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Toroidal modeling of plasma flow damping and density pump-out by RMP during ELM mitigation in HL-2A

Reduction of both the plasma density and toroidal flow speed, due to application of the predominantly n = 1 (n is the toroidal mode number) resonant magnetic perturbation (RMP) for controlling the edge localized mode in the HL-2A tokamak, is numerically investigated utilizing the quasi-linear initial-value code MARS-Q (Liu et al 2013 Phys. Plasmas 20 042503). Simulation results reveal that the neoclassical toroidal viscosity (NTV) due to three dimensional fields plays the key role in modifying the plasma momentum and particle transport in the HL-2A discharge. By comparing the modeling results with the measured density pump-out in the experiment, the electron NTV particle flux model, in combination with the free-boundary condition for the axisymmetric change of the density at the plasma edge, is found to yield the best agreement in terms of both the pump-out level and the overall time scale. Further sensitivity studies show that the simulated density pump-out level is reasonably robust against variations in the model assumptions, including the particle diffusion model and the non-ambipolar versus ambipolar NTV particle flux. The latter however affects the time scale for reaching the steady state solution. Finally, it is found that the plasma edge-peeling response, the NTV torque, as well as the plasma momentum and particle transport, all are sensitive to the toroidal phase difference between the upper and lower rows of the RMP coil currents in HL-2A, with the 30 ∘ coil phasing producing the minimal side effects on the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimal production of Phanerochaete chrysosporium manganese peroxidases and Trametes sp. C30 laccase hybrid Lac131 in Aspergillus niger for lignin bioconversion

Background Incorporating the production of related ligninolytic enzymes into industrial filamentous fungus Aspergillus niger will enhance the bioconversion of lignocelluloses to various chemical products. Results In this study, transgenic expression of Phanerochaete chrysosporium manganese peroxidases (mnps) and Trametes sp. C30 laccase hybrid Lac131 (lac131) were examined and optimized in A. niger 11414 prtT∆ strain. Five mnps (mnp1, mnp2, mnp3, mnp4, and mnp5) and lac131 genes were expressed separately or in combination. The transgenic strain containing the entire mnp2 genomic coding sequence (gmnp2) exhibited the highest mnP activity among the five mnp over-expression strains in the modified minimal medium (mMM) with addition of 5 g/L bovine hemoglobin (bHg). We examined the effects of hemin and bHg on mnP production in the gmnp2 strain cultures and found that at least 1 g/L bHg was required, while hemin was not. Culture conditions for mnP production were further optimized for the gmnp2 strain and the highest mnP activities were detected in the cultures grown at 25 °C and 200 rpm with an initial pH of 4.5. Effects of soy protein, skim milk, and bovine serum albumin on mnP production were investigated; 5 g/L of soy proteins or skim milk had comparable effects to 2.5 g/L bHg, while cultures with bovine serum albumin had diminished mnP activity. Disruption of both prtT and vsm1 substantially augmented the mnP production and its activity reached 575 U/L. Trametes sp. C30 laccase hybrid lac131 was strongly expressed in either A. niger gmnp2 (1975 U/L) or 11414prtT∆ (3895 U/L) strain. Both mnP and laccase in the culture supernatants effectively decolorized selected phenolic compounds (dyes) and cleaved tagged model lignin dimers. Conclusion The mnP was successfully produced in A. niger by optimizing the culture conditions and host strain. Co-expression of all four mnp genes in the same expression host by multiplex CRISPR will lead to the mnP production reaching levels comparable to P. chrysosporium, while only requiring 36 h at 25 °C. The Lac131 activity in transgenic A. niger strain is 4- to 7-times higher than that in previous studies. Co-production of mnP and laccase in A. niger will enhance the lignin bioconversion efficiency.

Aspergillus niger↗

SU-E-J-51: Dose Response of Common Solid State Detectors in Homogeneous Transverse and Longitudinal Magnetic Fields

Purpose: Solid state radiation detectors are often used for dose profiles and percent depth dose measurements. The dose response of selected solid state detectors is evaluated in varying transverse and longitudinal magnetic fields for eventual use in MR-Linac devices. Methods: A PTW 60003 and IBA PFD detector were modeled in the Monte Carlo code PENELOPE, incorporating a magnetic field which was varied in strength and oriented both transversely and longitudinally with respect to the incident photon beam. The detectors' long axis was in turn oriented either parallel or perpendicular to the photon beam. Dose to the active volume of each detector was scored, and its ratio to dose with zero magnetic field strength (dose response) was determined. Accuracy of the simulations was evaluated by measurements using both chambers taken at low field with a small electromagnet. Simulations were also performed in a water phantom to compare to the in air results. Results: Significant dose response was found in transverse field geometries, nearing 20% at 1.5T. The response is highly dependent on relative orientations to the magnetic field and photon beam, and on detector composition. Low field measurements confirm these results. In the presence of longitudinal magnetic fields, the detectors exhibit little dose response, reaching 0.5–1% at 1.5T regardless of detector orientation. Water tank simulations compared well to the in air simulations when not at the beam periphery, where in transverse magnetic fields only, the water tank simulations differed from the in air results. Conclusion: Transverse magnetic fields can cause large deviations in dose response, and are highly position orientation dependent. Comparatively, longitudinal magnetic fields exhibit little to no dose response in each detector as a function of magnetic field strength. Water tank simulations show longitudinal fields are generally easier to work with, but each detector must be evaluated separately.

07 ISOTOPE AND RADIATION SOURCES↗

Effect of Impurity Atoms on Thermonuclear Yield

The presence of impurity atoms in an inertial confinement fusion (ICF) implosion can greatly affect the burning of the thermonuclear fuel. This is a serious concern for attaining ignition as high atomic number atoms, such as carbon, can mix into the hot spot and reach the deuterium–tritium (DT) fuel, absorb energy from the implosion, and radiate away that energy. The temperature of the DT fuel drops and reduces the reactivity of the fuel resulting in lower thermonuclear yield. Fuel impurities can also be a useful diagnostic as when higher-Z atoms, such as Ar, Kr, and Xe, are purposefully mixed into the fuel. From the spectra emitted by these atoms, the temperature and density of the fuel can be deduced. The long history of using such tracer atoms by the national ICF program shows the great advances in experimental technique, diagnostic instruments, and atomic-physics modeling over the past 50 years. A particular set of experiments, the “High-Z” campaign, is reviewed in depth as this was the most systematic study of tracer element composition and quantity ever reported. Difficulties in interpreting the results of those experiments due to the limited code capabilities at that time are discussed. A main concern for the High-Z campaign was the then current state of non-Local Thermodynamic Equilibrium (nLTE) atomic physics models. The present state of nLTE modeling and its implementation in Los Alamos’s multi-physics codes is discussed in this report. To further the understanding of the complex interaction of high-Z atoms in the gas, a survey of relevant existing experimental platforms at the National Ignition Facility (NIF) and the Omega Laser Facility are presented along with the current state of pertinent diagnostics at NIF and Omega. This review concludes with a discussion of the trade-off available between temperature and volume of the imploded ICF capsule with the goal of maximizing the total yield.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

COSMIC-RAY POSITRONS FROM MILLISECOND PULSARS

Observations by the Fermi Large Area Telescope of γ-ray millisecond pulsar (MSP) light curves imply copious pair production in their magnetospheres, and not exclusively in those of younger pulsars. Such pair cascades may be a primary source of Galactic electrons and positrons, contributing to the observed enhancement in positron flux above ∼10 GeV. Fermi has also uncovered many new MSPs, impacting Galactic stellar population models. We investigate the contribution of Galactic MSPs to the flux of terrestrial cosmic-ray electrons and positrons. Our population synthesis code predicts the source properties of present-day MSPs. We simulate their pair spectra invoking an offset-dipole magnetic field. We also consider positrons and electrons that have been further accelerated to energies of several TeV by strong intrabinary shocks in black widow (BW) and redback (RB) systems. Since MSPs are not surrounded by pulsar wind nebulae or supernova shells, we assume that the pairs freely escape and undergo losses only in the intergalactic medium. We compute the transported pair spectra at Earth, following their diffusion and energy loss through the Galaxy. The predicted particle flux increases for non-zero offsets of the magnetic polar caps. Pair cascades from the magnetospheres of MSPs are only modest contributors around a few tens of GeV to the lepton fluxes measured by the Alpha Magnetic Spectrometer, PAMELA, and Fermi, after which this component cuts off. The contribution by BWs and RBs may, however, reach levels of a few tens of percent at tens of TeV, depending on model parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-T-562: Scanned Percent Depth Dose Curve Discrepancy for Photon Beams with Physical Wedge in Place (Varian IX) Using Different Sensitive Volume Ion Chambers

Purpose: To investigate and report the discrepancy of scanned percent depth dose (PDD) for photon beams with physical wedge in place when using ion chambers with different sensitive volumes. Methods/Materials: PDD curves of open fields and physical wedged fields (15, 30, 45, and 60 degree wedge) were scanned for photon beams (6MV and 10MV, Varian iX) with field size of 5x5 and 10x10 cm using three common scanning chambers with different sensitive volumes - PTW30013 (0.6cm3), PTW23323 (0.1cm3) and Exradin A16 (0.007cm3). The scanning system software used was OmniPro version 6.2, and the scanning water tank was the Scanditronix Wellhoffer RFA 300.The PDD curves from the three chambers were compared. Results: Scanned PDD curves of the same energy beams for open fields were almost identical between three chambers, but the wedged fields showed non-trivial differences. The largest differences were observed between chamber PTW30013 and Exradin A16. The differences increased as physical wedge angle increased. The differences also increased with depth, and were more pronounced for 6MV beam. Similar patterns were shown for both 5x5 and 10x10 cm field sizes. For open fields, all PDD values agreed with each other within 1% at 10cm depth and within 1.62% at 20 cm depth. For wedged fields, the difference of PDD values between PTW30013 and A16 reached 4.09% at 10cm depth, and 5.97% at 20 cm depth for 6MV with 60 degree physical wedge. Conclusion: We observed a significant difference in scanned PDD curves of photon beams with physical wedge in place obtained when using different sensitive volume ion chambers. The PDD curves scanned with the smallest sensitive volume ion chamber showed significant difference from larger chamber results, beyond 10cm depth. We believe this to be caused by varying response to beam hardening by the wedges.

07 ISOTOPE AND RADIATION SOURCES↗

Bay Area Regional Energy: Network Integrated Commercial Retrofits (BRICR) Project. Final Report

The BRICR project applied large-scale building energy modeling concepts with the aim of reducing the cost of energy efficiency targeting, design, and project development, and measurement of energy savings for energy efficiency programs implemented by local governments that serve small and medium commercial buildings (SMB). The project leveraged the services and resources of existing local government energy programs serving disadvantaged and hard-to-reach SMB customers. In contrast to programs run by utilities, local government programs generally do not have direct access to energy billing records for an entire class of customers in a geographic area, which prior research demonstrated useful for large-scale building energy model baseline development and calibration. , However, local governments are rich in public records that offer important clues about physical attributes and uses that, along with behavior, determine energy use. Relying only on public records, BRICR demonstrated development of credible baseline energy models for 3,792 office, retail, and hotel buildings. Publicly disclosed annual energy use data from a local energy benchmarking program and anonymized data from the Building Performance Database, the nation’s largest dataset about energy-related characteristics of buildings, were utilized to validate and calibrate energy models via an innovative method comparing distributions of energy intensity by fuel type for portfolios of buildings of similar size, vintage, and use. Portfolio calibration does not provide certainty that an energy model fits an individual building; the method is useful when billing data is not accessible – a common situation for researchers, energy service providers and ESCOs, local governments, and any party other than a utility. A software component was developed, the BRICR gem, which automates simulation when relevant data is added or edited by the user to a file saved in the standardized BuildingSync XML schema for energy audit data. The component was demonstrated as a simplified means to generate a mass of energy models corresponding to public records containing basic attributes such as building scale, location, use, year built, and aspect ratio in combination with building energy code prototype data corresponding to use and vintage. The component was also demonstrated as a simplified means to automate energy simulation when attributes are revised; the intention was to enable iterative improvement of the baseline model and energy savings estimates for common energy conservation measures as users revise relevant attributes based on their observations. In the context of institutional change and uncertainty for the participating local government energy programs, 13 whole building retrofits were completed. Impacts were measured by applying the CalTRACK2.0 methods to standardize measurement of normalized metered energy consumption. The GRIDMeter methods of stratified sampling and individual load shape analysis were applied to adjust for impacts of the effect of COVID-19 on retrofitted buildings in the context of all local buildings of similar size and use. Excluding impacts of the pandemic, retrofitted buildings demonstrated between 1.6% and 25.1% reduction in energy use. The project contributed use cases and feedback that helped inform evolution of the software tools and data formats that were combined for the first time in the BRICR project.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

A finite element procedure for radio-frequency sheath–plasma interactions based on a sheath impedance model

A finite element code that solves self-consistent radio-frequency (RF) sheath-plasma interaction problems is improved by incorporating a generalized sheath boundary condition in the macroscopic solution scheme. This sheath boundary condition makes use of a complex sheath impedance including both the sheath capacitance and resistance, which enables evaluation of not only the RF voltage across the sheath but also the power dissipation in the sheath. The newly developed finite element procedure is applied to cases where the background magnetic field is perpendicular to the sheath surface in one- and two-dimensional domains filled by uniform low- and high-density plasmas. The numerical results are compared with those obtained by employing the previous capacitive sheath model at a typical frequency for ion cyclotron heating used in fusion experiments. It is shown that for sheaths on the order of 100 V in a high-density plasma, localized RF power deposition can reach a level which causes material damage. It is also shown that the sheath-plasma wave resonances predicted by the capacitive sheath model do not occur when parameters are such that the generalized sheath impedance model substantially modifies the capacitive character of the sheath. Here, possible explanations for the difference in the maximum RF sheath voltage depending on the plasma density are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical investigation of laser-driven shock interaction with a deformable particle

A laser-driven shock propagating through an isolated particle embedded in a plastic (CH) target was studied using the radiation-hydrodynamic code FLASH. Preliminary simulations using IONMIX equations of state (EOS) showed significant differences in the shock Hugoniot of aluminum compared to experimental data in the low-pressure regime [O(10) GPa], resulting in higher streamwise compression and deformation of an aluminum particle. Hence, a simple modification to the ideal gas EOS was developed and employed to describe the target materials and examine the particle dynamics. The evolution of the pressure field demonstrated a complex wave interaction, resulting in a highly unsteady particle drag which featured two drag minima due to shock focusing at the rear end of the particle and rarefaction stretching due to laser shut-off. Although ~30% lateral expansion and ~25% streamwise compression were observed, the aluminum particle maintained considerable integrity without significant distortion. Additional simulations examined the particle response for a range of particle densities, sizes, and acoustic impedances. The data revealed that lighter particles such as aluminum gained significant momentum, reaching up to ~96% of the shocked CH's speed, compared to ~29% for the heavier tungsten particles. Despite the differences seen in the early stage of shock interaction, particles with varying acoustic impedances ultimately reached the same peak velocity. This identified particle-to-host density ratio is an important factor in determining the inviscid terminal velocity of the particle. In addition, the modified EOS model presented in this study could be used to approximate solid materials in hydrocodes that lack material strength models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

TEAMER: Pitching Foil Crossflow Turbine Efficiency Data

This dataset documents the efficiency testing of a pitching foil crossflow turbine, conducted at the University of New Hampshire's (UNH) Chase Ocean Engineering Laboratory tow tank facility. The tests explored various pitch phases and amplitudes, ranging from 0 to 18 degrees, across different flow speeds and turbine RPMs. Specifically, testing was performed at pitch amplitudes of 0,3,6,9,12,15 and 18 degrees. The results indicate that turbine efficiency improves with increased pitch amplitude, reaching optimal performance at 12 degrees. Modifications were made to the test frame by UNH to enhance the stiffness of the support structure. The testing specifically mapped the efficiency of a single quadrant of the Ocean Renewable Power Company (ORPC) Autonomous Turbine Generator Unit (ATGU), covering the full range of allowable pitching amplitudes and phases. The dataset provides both raw and processed data, including drag and torque measurements, and contains Python scripts used for data processing and visualization. This research was supported by funding from TEAMER RFTS 1 (Request for Technical Support).

16 TIDAL AND WAVE POWER↗

Implementation and Comparison of Different Algorithms for on the fly Doppler Broadening in RMC Code

Accurate cross sections are of great importance to Monte Carlo particle transport calculation codes which use pointwise cross sections. It is well known that cross sections have much dependence on temperature on account of Doppler effect which is particularly vital for resonance absorption cross sections. In general, pre-prepared continuous-energy libraries are used in Monte Carlo codes, which is feasible in many applications. However, in multi-physics calculations where neutronics and thermal hydraulics are coupled, cross sections at a fine interval of temperatures are desired, which brings a great challenge to computer memory. Consequently, different methods about on the fly Doppler broadening for cross sections are proposed, which can generate cross sections at the specific temperature points during neutron simulations. These methods can be roughly classified into two categories: one is to calculate the required cross sections by different algorithms, such as interpolation using a high order function expansion, multipole representation, and different integration algorithms based on Doppler broadening equation. The other, for instance, target motion sampling (TMS), does not calculate the Doppler broadening cross sections but takes into account Doppler effect by rejection technique. In this paper, several on the fly Doppler broadening algorithms based on Doppler broadening equation are implemented in the Reactor Monte Carlo code RMC, and then the accuracy and efficiency are discussed and compared. In this paper, several algorithms for on the fly Doppler broadening based on Doppler broadening equation have been implemented in RMC and the focus falls on the comparison of the Gauss Legendre quadrature and Gauss Hermite quadrature. Gauss Legendre quadrature is more accurate than Gauss Hermite quadrature but also more time consuming. Base temperature has a significant impact on Gauss Hermite quadrature while increasing base temperature to room temperature (300 K) can give satisfactory results. In addition, the proportion of overhead of Doppler broadening treatment in total computational time tends to decrease with an increase number of nuclides in fuel. For depletion calculation, both algorithms are in good agreement with reference and the relative errors of K{sub inf} are within 5 times of standard deviation. As for total computational time, the ratio between Gauss Legendre quadrature or Gauss Hermite quadrature and reference decreases and reaches 9.6 or 3.4 at last, respectively. Four important isotopes are chosen to compare the atom density and all relative errors are within 0.25%. After comparisons of different algorithms and different base temperatures, it can be found that Gauss Hermite quadrature broadening from 300 K can perform on the fly cross section Doppler broadening with least total computational time and acceptable precision. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗