Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “critical points”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Atomistic Probing of Defect-Engineered 2H-MoTe 2 Monolayers

Point defects dictate various physical, chemical, and optoelectronic properties of two-dimensional (2D) materials, and therefore, a rudimentary understanding of the formation and spatial distribution of point defects is a key to advancement in 2D material-based nanotechnology. In this work, we performed the demonstration to directly probe the point defects in 2H-MoTe 2 monolayers that are tactically exposed to (i) 200 °C-vacuum-annealing and (ii) 532 nm-laser-illumination; and accordingly, we utilize a deep learning algorithm to classify and quantify the generated point defects. We discovered that tellurium-related defects are mainly generated in both 2H-MoTe 2 samples; but interestingly, 200 °C-vacuum-annealing and 532 nm-laser-illumination modulate a strong n-type and strong p-type 2H-MoTe 2 , respectively. While 200 °C-vacuum-annealing generates tellurium vacancies or tellurium adatoms, 532 nm-laser-illumination prompts oxygen atoms to be adsorbed/chemisorbed at tellurium vacancies, giving rise to the p-type characteristic. This work significantly advances the current understanding of point defect engineering in 2H-MoTe 2 monolayers and other 2D materials, which is critical for developing nanoscale devices with desired functionality.

2H-MoTe 2↗
DOE OSTI2020

Superuniversality from disorder at two-dimensional topological phase transitions

We investigate the effects of quenched randomness on topological quantum phase transitions in strongly interacting two-dimensional systems. We focus first on transitions driven by the condensation of a subset of fractionalized quasiparticles (“anyons”) identified with “electric charge” excitations of a phase with intrinsic topological order. All other anyons have nontrivial mutual statistics with the condensed subset and hence become confined at the anyon condensation transition. Furthermore, using a combination of microscopically exact duality transformations and asymptotically exact real-space renormalization group techniques applied to these two-dimensional disordered gauge theories, we argue that the resulting critical scaling behavior is “superuniversal” across a wide range of such condensation transitions and is controlled by the same infinite-randomness fixed point as that of the 2D random transverse-field Ising model. We validate this claim using large-scale quantum Monte Carlo simulations that allow us to extract zero-temperature critical exponents and correlation functions in (2+1)D disordered interacting systems. We discuss generalizations of these results to a large class of ground-state and excited-state topological transitions in systems with intrinsic topological order as well as those where topological order is either protected or enriched by global symmetries. When the underlying topological order and the symmetry group are Abelian, our results provide prototypes for topological phase transitions between distinct many-body localized phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Interactions of clathrate hydrate promoters sodium dodecyl sulfate and tetrahydrofuran investigated using 1H diffusion nuclear magnetic resonance at hydrate-forming conditions

Thermodynamic hydrate promoters and kinetic hydrate promoters can be used to reduce the P–T conditions for clathrate hydrate synthesis to decrease the nucleation induction time while increasing growth rates. Two commonly used promoters for hydrate research are tetrahydrofuran (THF) and sodium dodecyl sulfate (SDS), which can increase the overall hydrate promotion when used in tandem as compared to individually. There are several molecular theories regarding how SDS promotes hydrate growth. This study explores the micellular theory, for which hydrate formation depends on surfactant aggregates (micelles) at a critical micelle concentration (CMC) to increase the interfacial surface area. The micellular theory is the most investigated and criticized surfactant hydrate promotion theory. To address questions related to micellar behavior, this study investigates the intermolecular behavior between SDS and THF for the identification of micelles at hydrate-forming conditions. The systems explored contained THF at 3 and 5 wt. % with varying concentrations of SDS below and above the CMC. Several methods including a qualitative visual method, conductivity, interfacial tensiometry, 13C Liquid-state Nuclear Magnetic Resonance (NMR) spectroscopy, and 1H diffusion NMR spectroscopy were evaluated at temperatures below the Krafft point of SDS and above 0 °C. The presence of THF at low concentrations decreased the critical temperature for the formation of SDS micelles, where SDS is solubilized in THF/water solution at hydrate-forming temperatures without precipitation. The CMC of SDS was decreased significantly even at hydrate-forming conditions. Mixed surfactant–cosolvent micellular behavior of SDS in the presence of low concentrations of THF was confirmed at hydrate-forming conditions above 0 °C.

Chemistry↗
DOE OSTI2022

A Note on the Probability of Initiation Problem

The purpose of this note is to draw attention to Pál for his derivation of the Pál-Bell equation which is used by LLNL and LANL to model Probability of Initiation (POI) problems. Both laboratories credit the equation to Bell and Lee. Although the equation formulated by Bell and Lee and the equation derived by Pál appear to be the same, the equation by Bell and Lee is flawed. Scattering, which was absent in Bell’s original formulation of a POI problem, was later introduced by Bell and Lee by replacing a fission term with a scattering term. Such a replacement, however, leads to an incorrect average number $\overline{v}$ of neutrons that is produced in a fission event. By approaching a POI problem from probabilistic point of view, Pál formulates the problem as a branching process that i) describes the growth of a population of neutrons in a fissile medium, and ii) predicts $\overline{v}$ correctly. There are five reasons for writing this report. The first is to clarify the failing in Bell and Lee’s derivation. The second is to establish the origin of the equation used by the laboratories. The third is to elucidate the fixed point technique formulated by Mokhtar-Kharroubi to solve the Pál-Bell equation. The fourth is to derive the α eigenvalue that is associated with a super-critical situation. The fifth is to describe in the Appendix a second order positivity preserving technique that accelerates the convergence of the first order fixed point method developed by Mokhtar-Kharroubi.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Report on Year-4 of Water NSTF Matrix Testing: Facility Maintenance and Accident Testing

Under support from the Department of Energy (DOE) and the Office of Advanced Reactor Technologies (ART), a large-scale test facility has been constructed at Argonne National Laboratory to generate NQA-1 qualified validation data for passive decay heat removal systems in advanced reactors. The Natural convection Shutdown heat removal Test Facility (NSTF) reflects key features of a ½ scale, water-based, Reactor Cavity Cooling System (RCCS) and is intended to study the behavior, bound performance, and ultimately guide design decisions for passive decay heat removal systems for advanced reactors. In addition to the experimental activities detailed in this report, a supportive computational modeling effort is on-going which has been demonstrated to significantly strengthen the experimental program while also improving accuracy of the computer models. Together these create a mutually beneficial relationship integral to meeting the overall program objective of examining the heat removal performance of the RCCS concept. This report serves as a summary of maintenance and experimental activities during the program’s fourth year of water-based operation. A planned six-month maintenance period began in August 2021, during which major inspections, repairs, cleaning, and installation of new instrumentation and data acquisition hardware were conducted. Most significantly, two heaters that faulted during Year-3 were repaired, allowing the facility to resume use of the full heated section area and full range of available electric power. The remainder of the year consisted of eight months of test operations, during which the facility logged 211 hours of active heating across one bake-out (following the maintenance period) and seven matrix tests; five classifieds as Accepted per NQA-1, one as Trending, and one as Failed. Testing began by performing two repeat cases to confirm expected facility response and behavior during both single- and two-phase flow conditions, ensuring no changes were introduced during the maintenance period that might have altered the thermal-hydraulic characteristics of the facility. In continuation of the power parametric series initiated in previous years, a high-power test case was then performed examining heat removal performance at a decay heat load equivalent to 2.4 MWt, full-scale, a level exceeding maximum design targets. Additional testing then introduced various blockage scenarios along the network piping, examining the effects of partial and complete blockages of the flow paths on the system behavior and heat removal performance. A study of static boiling tests directed at understanding the geysering two-phase instability was also conducted. The loop was filled only to the bottom of the tank outlet, creating an open loop configuration that prevents any natural circulation flow from occurring, and the heaters were powered on until the facility reached saturation conditions. Following, a series of quasi-steady-state conditions were introduced by adjusting the inventory level in the adiabatic chimney piping at decreasingly lower elevations above the heated region. A strong correlation of geysering characteristics and loop level was observed, with flow and temperature excursions decreasing in intensity, but increasing in frequency, as the fill were reduced to lower elevations along the chimney piping. Once the level fell very low in the chimney, at points near the top of the heated section, the system reached a stable state of continuous boiling without any occurrence of geysering eruptions. A final significant testing accomplishment this year was successful completion of an “accident scenario” test, whose operating conditions were based on a prototypic decay heat curve provided by Framatome and scaled for the NSTF. This test began by establishing steady-state, single-phase “normal operation” conditions, before simulating an accident trip where the availability of active cooling systems was lost. Loop temperatures gradually increased until reaching saturation and subsequent two-phase boiling flow. Over the course of an extended operational period along the defined decay heat curve, steam boil-off caused gradual but continued depletion of liquid inventory until reaching a critically low level causing flow stagnation and cessation of natural circulation heat removal. At this point, after nearly 72 hours of continuous operation, a cold refill was performed to replenish the system inventory and allow the facility to re-establish closed loop natural circulation flow and return to a safe operational state.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Advances in physics and applications of 3D magnetic perturbations on the J-TEXT tokamak

In the last two years, three major technical improvements have been made on J-TEXT in supporting of the expanded operation regions and diagnostic capabilities. (1) The successful commission of the 105 GHz/500 kW/1 s electron cyclotron resonance heating (ECRH) system increasing the core electron temperature from 0.9 keV up to around 1.5 keV. (2) The poloidal divertor configuration with an X-point in the high-field side has been achieved. In particular, the 400 kW electron cyclotron wave has also been successfully injected into the diverted plasma. (3) A 256-channel electron cyclotron emission imaging diagnostic system and two sets of four-channel Doppler backscattering diagnostics have been successfully developed on J-TEXT, allowing detailed measurement of the electron temperature and density fluctuations for turbulence and MHD research. The locked mode (LM), especially the 2/1 LM, is one of the biggest threats to the plasma operation. Both the thresholds of 2/1 and 3/1 LM are observed to vary non-monotonically on electron density. The electrode biasing was applied successfully to unlock the LM from either a rotating or static resonant magnetic perturbation (RMP) field. In the presence of 2/1 LM, three kinds of standing wave (SW) structures have been observed to share a similar connection to the island structure, i.e. the nodes of the SWs locate around the O- or X-points of the 2/1 island. The control and mitigation of disruption is essential to the safe operation of ITER, and it has been systematically studied by applying a RMP field, massive gas injection (MGI) and shattered pellet injection on J-TEXT. When the RMP-induced 2/1 LM is larger than a critical width, the MGI shutdown process can be significantly influenced. If the phase difference between the O-point of LM and the MGI valve is +90° (or –90°), the penetration depth and the assimilation of impurities can be enhanced (or suppressed) during the pre-thermal quench (TQ) phase and result in a faster (or slower) TQ. A secondary MGI can also suppress the runaway electron (RE) generation, if the additional high-Z impurity gas arrives at the plasma edge before TQ. When the secondary MGI has been applied after the formation of the RE current plateau, the RE current can be dissipated, and the dissipation rate increases with the injected impurity quantity but saturates with a maximum of 28 MA s –1 . The non-local transport is experimentally observed in the ion transport channel. The electron thermal diffusivity significantly increases with the ECRH power. Theoretical work shows that significant intrinsic current can be driven by electromagnetic turbulence, and the robust formation mechanism of the E × B staircase is identified from the Hasegawa–Wakatani system.

divertor↗
DOE OSTI2023

Estimating source-sink distributions and fluxes of reactive nitrogen and sulfur within a mixed forest canopy

The vertical source-sink distribution of air pollutants within and above forested canopies is necessary for describing the biological, physical, and chemical processes influencing the soil-vegetation-atmosphere exchange. Here, this study implemented inverse modeling methods to estimate the source-sink and flux profiles of reactive nitrogen (N) and sulfur (S) compounds from measurements of the mean concentration profiles of ammonia (NH 3 ), nitric acid (HNO 3 ), sulfur dioxide (SO 2 ), and particulate ammonium (NH 4 + ), nitrate (NO 3 − ), and sulfate (SO 4 2− ) at a forest site in the southern Appalachian Mountains. Three inverse approaches utilizing different approximations to scalar transport within the canopy were developed and evaluated against sensible heat flux measurements. The Eulerian model (EUL), which incorporates vertical velocity skewness, performed well in reproducing the turbulent heat fluxes and was subsequently used to calculate the chemical source-sink and flux profiles. Above-canopy fluxes of NH 3 were downward, indicating that the forest was a net sink of NH 3 . The soil/litter layer was both a source and a sink for NH 3 but the exchange rate at the forest floor was small. Fluxes of HNO 3 , SO 2 , NO 3 − , NH 4 + , and SO 4 2- were uni-directional (deposition only) between the air and the canopy/ground and increased monotonically from the forest floor to the canopy top. Crown foliage dominated the uptake of reactive N and S during the growing season, accounting for 80–90% of the total canopy-scale flux. Fluxes and canopy-ground partitioning estimated using the resistance-based Surface Tiled Aerosol and Gas Exchange (STAGE) model were generally comparable to EUL. The comparison highlights the need for improved parameterizations of litter exchange and NH 3 compensation points in resistance models for forest ecosystems. The findings here benefit the application of critical loads in forest ecosystems and guide further development of resistance-based exchange models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Hydrothermal liquefaction of wastewater-grown algae to produce synthetic aviation fuel: A combined experimental study and techno-economic assessment

Large-scale algae farms may someday become a consistent source of biomass feedstock for biofuels. Near-term supplies of algal biomass are available at certain water resource recovery facilities as algae cultivation is used as a method for nutrient recovery from specific effluent streams. Algae grown as a service shifts the value to the service rather than its sole use as a feedstock, which could enable the provision of algal biomass at low to no cost to biofuel producers. Hydrothermal liquefaction (HTL) can readily upgrade wet feedstock slurries, such as algae, to produce a carbon-enriched biocrude. The HTL biocrude can be hydrotreated and distilled, producing a variety of distillate fuels, including synthetic aviation fuel (SAF). We present a pathway, showing the experimental production of SAF from wastewater-grown algae via HTL, along with a techno-economic assessment to identify opportunities for process improvements. Critical quality attributes of the SAF, such as density, viscosity, surface tension, and freeze point, were estimated within the expected fuel experience ranges when compared against petroleum jet fuel. The average minimum fuel selling price of fuels from wastewater-grown algae for breakeven economics was $\$9.04$ per gasoline gallon equivalent (GGE). The sale of co-products such as struvite fertilizers and cement additives can add revenue to reduce the net cost. Ultimately, the selling price is influenced by the scale of the HTL processing facility. Adjusting estimations in the process scale, algae yield, and capital cost estimation can lower the price to $\$6.51$/GGE or raise it to $13.07/GGE.

Biofuels↗
DOE OSTI2023

Dislocation loop bias and void swelling in irradiated α-iron from mesoscale and atomistic simulations

Abstract Dislocation loops are ubiquitous in irradiated materials, and dislocation loop bias plays a critical role in void swelling. However, due to complicated interactions between dislocation loops and point defects, it is challenging to evaluate the bias factors of dislocation loops. Here, we determine the bias of sessile < 100 > loops in α-iron using a recently developed atomistic approach based on the lifetime of point defects. We establish a mechanistic understanding of the loop interaction based on the diffusion tendency of point defects near the loop core region. Mobile self-interstitial atoms tend to be absorbed from the edge of the loop, and a trapping region perpendicular to the habit plane of the loop exists. The dislocation loop bias is found to be substantially lower than those of straight dislocations in α-iron and should be included in swelling rate estimates. With the obtained sink strength and bias values, agreement is achieved with experimental results for both absolute values and temperature dependence.

Yu, Ziang↗
DOE OSTI2022

Solving quasi-free and quadratic Lindblad master equations for open fermionic and bosonic systems

The dynamics of Markovian open quantum systems are described by Lindblad master equations. For fermionic and bosonic systems that are quasi-free, i.e. with Hamiltonians that are quadratic in the ladder operators and Lindblad operators that are linear in the ladder operators, we derive the equation of motion for the covariance matrix. This determines the evolution of Gaussian initial states and the steady states, which are also Gaussian. Using ladder super-operators (a.k.a. third quantization), here we show how the Liouvillian can be transformed to a many-body Jordan normal form which also reveals the full many-body spectrum. Extending previous work by Prosen and Seligman, we treat fermionic and bosonic systems on equal footing with Majorana operators, shorten and complete some derivations, also address the odd-parity sector for fermions, give a criterion for the existence of bosonic steady states, cover non-diagonalizable Liouvillians also for bosons, and include quadratic systems. In extension of the quasi-free open systems, quadratic open systems comprise additional Hermitian Lindblad operators that are quadratic in the ladder operators. While Gaussian states may then evolve into non-Gaussian states, the Liouvillian can still be transformed to a useful block-triangular form, and the equations of motion for k-point Green’s functions form a closed hierarchy. Based on this formalism, results on criticality and dissipative phase transitions in such models are discussed in a companion paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Stasis in an expanding universe: A recipe for stable mixed-component cosmological eras

One signature of an expanding universe is the time variation of the cosmological abundances of its different components. For example, a radiation-dominated universe inevitably gives way to a matter-dominated universe, and critical moments such as matter-radiation equality are fleeting. In this paper, we point out that this lore is not always correct and that it is possible to obtain a form of “stasis” in which the relative cosmological abundances $Ω_i$ of the different components remain unchanged over extended cosmological epochs, even as the universe expands. Moreover, here we demonstrate that such situations are not fine-tuned but are actually global attractors within certain cosmological frameworks, with the universe naturally evolving toward such long-lasting periods of stasis for a wide variety of initial conditions. The existence of this kind of stasis therefore gives rise to a host of new theoretical possibilities across the entire cosmological timeline, ranging from potential implications for primordial density perturbations, dark-matter production, and structure formation all the way to early reheating, early matter-dominated eras, and even the age of the Universe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Isomeric yield ratios of fission products: A missing piece in reactor antineutrino summation calculations

The calculation of the spectrum of antineutrinos ($\bar{v}_e$) from a reactor is a complicated problem requiring several nuclear data and physics inputs. Many of these have been investigated thoroughly to improve calculations and properly account for uncertainties. One input which has heretofore escaped consideration is the fission-yield distribution between ground and isomeric states. Here, in this work, we explore the effect of incorporating newly evaluated isomeric yield ratios (IYR) for 43 fission products into summation calculations and identify the disproportionate importance of certain isotopes, particularly at higher energies. Our analysis indicates that updated IYRs contribute to a significant increase in the $\bar{v}_e$ spectrum around and above 7 MeV, with increases of more than 50% at higher energies. Through a detailed sensitivity study, we highlight a number of isotopes for which the IYR has a substantial effect on the $\bar{v}_e$ spectrum. This work stresses the critical role of isomeric yields in calculations of reactor $\bar{v}_e$ spectra and points to the necessity for their accurate experimental determination, especially for key fission products, in order to refine our understanding and address the observed anomalies between measured and calculated $\bar{v}_e$ spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Towards Automated Reasoning Chains for Verification of LLM-Generated Scientific Code

With the rise of Large Language Model (LLM) generated code, including in domains like scientific computing, ensuring not only syntactical, but also mathematical correctness, has become a critical task. Traditional formal methods approaches often struggle with the ambiguity of floating-point code, and full symbolic execution is extremely costly and limited. We propose a chain-of-reasoning approach that iteratively lifts basic semantics from code into the SPIRAL system and then establishes numerical equivalency to the desired mathematical operation. Here, we leverage the ample mathematical knowledge already formalized in SPIRAL to enable the system to recognize not just different implementations of the same algorithm but fully separate approaches to solving the given problem. The chain establishes tight error bounds on the output of given code with respect to the true continuous solution it approximates, quantifying all sources of error. We demonstrate this approach by establishing the correctness of a pseudospectral solver for a simple 1-dimensional Poisson problem.

Oschatz, Quentin [Carnegie Mellon University,Pitts↗
DOE OSTI2021

Cooperative Inverters to Overcome PLL Malfunctions

This paper presents a cooperative strategy between smart inverters to overcome instability and undesired trips of inverters caused by Phase-Locked Loop (PLL) malfunctions. For the grid-interactive inverter, PLL is typically required to measure the phase angle of the AC voltage at the Point of Common Coupling (PCC) of the inverter. Accurate phase-angle measurement is critical for grid-interactive inverters to stay synchronized with the power grid. If the inverters are interconnected through a communication link, they can support each other by sharing their grid information, such as phase-angles at their PCCs. The cooperative strategy utilizes the phase-angles and voltages at PCCs in a neighborhood to identify any possible PLL malfunctions or measurement inaccuracies and then enable the inverters to make corrective decisions. The efficacy of the cooperative strategy has been verified through simulation results in MATLAB/Simulink environment in this paper.

Sadeque, Fahmid↗
DOE OSTI2025

Isomeric Yield Ratios of fission products: a missing piece in reactor antineutrino summation calculations

The calculation of the spectrum of antineutrinos ( $\overline{V}$ e ) from a reactor is a complicated problem requiring several nuclear data and physics inputs. Many of these have been investigated thoroughly to improve calculations and properly account for uncertainties. One input which has heretofore escaped consideration is the fission yield distribution between ground and isomeric states. In this work, we explore the effect of incorporating newly evaluated isomeric yield ratios (IYR) for 43 fission products into summation calculations and identify the disproportionate importance of certain isotopes, particularly at higher energies. Our analysis indicates that updated IYRs contribute to a significant increase in the $\overline{V}$ e spectrum around and above 7 MeV, with increases of more than 50% at higher energies. Through a detailed sensitivity study, we highlight a number of isotopes for which the IYR has a substantial effect on the $\overline{V}$ e spectrum. This work stresses the critical role of isomeric yields in calculations of reactor $\overline{V}$ e spectra and points to the necessity for their accurate experimental determination, especially for key fission products, in order to refine our understanding and address the observed anomalies between measured and calculated $\overline{V}$ e spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Understanding Subsidies to Achieve Diesel Powertrain Financial Parity for Heavy-Duty Fuel Cell Electric Vehicles

The development of a long-term sustainable hydrogen energy economy for commercial vehicle transportation will need to overcome key critical technical and logistics considerations in the near term. As compared to zero-emission powertrains, fossil-fuel-based powertrains provide mission flexibility and high uptime at a comparatively low total cost of ownership (TCO). While the incumbent carbon-intensive powertrains suffer from poor efficiency and are not sustainable to support global climate change initiatives in transportation decarbonization, techno-economic challenges continue to create complex barriers to the large-scale displacement of these with highly electrified powertrains architectures. Here, this article specifically addresses opportunities that well-targeted subsidies would afford in achieving fuel cell electric powertrain financial parity with diesel powertrains in heavy-duty trucks (HDTs). It shows that a rigorous assessment of system TCO factors will be critical in designing a viable subsidy strategy. Based on powertrain component and fuel price points, subsidy variations (for both capital expenses and operating fuel expenses) are assessed and quantified to determine where and how much subsidies should be leveraged to accelerate achieving diesel system financial parity. It shows that subsidies on hydrogen fuel expenses (and associated constituents) will play an important role in achieving this goal while subsidies on capital expenses (the more typical subsidy formulation) may not always be the preferred subsidy strategy. The results provide subsidy agencies with clear guidance and a database approach to characterize and implement viable HDT hydrogen powertrain system subsidy strategies. The approach may be readily tailored for other applications and system architectures.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Exploring ship track spreading rates with a physics-informed Langevin particle parameterization

Abstract. The rate at which aerosols spread from a point source injection, such as from a ship or other stationary pollution source, is critical for accurately representing subgrid plume spreading in a climate model. Such climate model results will guide future decisions regarding the feasibility and application of large-scale intentional marine cloud brightening (MCB). Prior modeling studies have shown that the rate at which ship plumes spread may be strongly dependent on meteorological conditions, such as precipitating versus non-precipitating boundary layers and shear. In this study, we apply a Lagrangian particle model (PM-ABL v1.0), governed by a Langevin stochastic differential equation, to create a simplified framework for predicting the rate of spreading from a ship-injected aerosol plume in sheared, precipitating, and non-precipitating boundary layers. The velocity and position of each stochastic particle is predicted with the acceleration of each particle being driven by the turbulent kinetic energy, dissipation rate, momentum variance, and mean wind. These inputs to the stochastic particle velocity equation are derived from high-fidelity large-eddy simulations (LES) equipped with a prognostic aerosol–cloud microphysics scheme (UW-SAM) to simulate an aerosol injection from a ship into a cloud-topped marine boundary layer. The resulting spreading rate from the reduced-order stochastic model is then compared to the spreading rate in the LES. The stochastic particle velocity representation is shown to reasonably reproduce spreading rates in sheared, precipitating, and non-precipitating cases using domain-averaged turbulent statistics from the LES.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Review of Wave Energy Converter Power Take-Off Systems, Testing Practices, and Evaluation Metrics: Preprint

While the field of wave energy has been the subject of numerical simulation, scale model testing, and precommercial project testing for decades, wave energy technologies remain in the early stages of development and must continuing proving themselves as a promising modern renewable energy field. One of the difficulties that wave energy systems have been struggling to overcome is the design of highly efficient energy conversion systems that can convert the mechanical power, derived from the oscillation of wave activated bodies, into another useful product. Often the power take-off (PTO) is defined as the single unit responsible for converting mechanical power into another usable form such as electricity, pressurized fluid, compressed air, and others. The PTO, and the entire power conversion chain (PCC), is of great importance as it affects not only how efficient wave power is converted into electricity, but also contributes to the mass, size, structural dynamics, and levelized cost of energy (LCOE) of the wave energy converter (WEC). Unlike wind and solar, there is no industrial standard device, or devices, for wave energy conversion and this diversity is transferred to the PTO system. The majority of current WEC PTO systems incorporate a mechanical or hydraulic drive train, power generator, and an electrical control system. The challenge of WEC PTO designs is designing a mechanical-to-electrical component that can efficiently convert irregular, bi-directional, low frequency and low alternating velocity wave motions. While gross average power levels can be predicted in advance, the variable wave elevation input has to be converted into smooth electrical output and hence usually necessitates some type of energy storage system, such as battery storage, accumulator super capacitors, etc., or other means of compensation such as an array of devices. One of the primary challenges for wave energy converter systems is the fluctuating nature of wave resources, which require WEC components to be designed to handle loads (i.e. torques, forces, and powers) that are many times greater than the average load. This approach requires a much greater PTO capacity than the average power output and indicates a higher cost. In addition, supporting mechanical coupling and or gearing can be added to the PCC to help alleviate the difficulties with transmission and control of fluctuating large loads with low frequencies (indicative of wave forcing) into smaller loads at higher frequencies (optimum for conventional electrical machine design) can quickly increase the complexity of the PCC which could result in a greater number of failure modes and increased maintenance costs. All of the previous points demonstrate how the PTO influences WEC dynamics, reliability, performance and cost which are critical design factors. This paper further explores these topics by providing a review of the state-of-the-art PTO systems currently under development, how these novel PTO systems are tested and derisked prior to precommercial deployment, and the evaluation metrics historically used to differentiate between PTO designs and how they can be improved to support control co-design focused development of wave energy systems.

laboratory testing↗