Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LOW SPEED STABILITY”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Propagation of a Stress Pulse in a Heterogeneous Elastic Bar

The propagation of a wave pulse due to low-speed impact on a one-dimensional, heterogeneous bar is studied here. Due to the dispersive character of the medium, the pulse attenuates as it propagates. This attenuation is studied over propagation distances that are much longer than the size of the microstructure. A homogenized peridynamic material model can be calibrated to reproduce the attenuation and spreading of the wave. The calibration consists of matching the dispersion curve for the heterogeneous material near the limit of long wavelengths. It is demonstrated that the peridynamic method reproduces the attenuation of wave pulses predicted by an exact microstructural model over large propagation distances.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Behavior and mechanisms of Doppler wind lidar error in varying stability regimes

Abstract. Wind lidars are widespread and important tools in atmospheric observations. An intrinsic part of lidar measurement error is due to atmospheric variability in the remote-sensing scan volume. This study describes and quantifies the distribution of measurement error due to turbulence in varying atmospheric stability. While the lidar error model is general, we demonstrate the approach using large ensembles of virtual WindCube V2 lidar performing a profiling Doppler-beam-swinging scan in quasi-stationary large-eddy simulations (LESs) of convective and stable boundary layers. Error trends vary with the stability regime, time averaging of results, and observation height. A systematic analysis of the observation error explains dominant mechanisms and supports the findings of the empirical results. Treating the error under a random variable framework allows for informed predictions about the effect of different configurations or conditions on lidar performance. Convective conditions are most prone to large errors (up to 1.5 m s−1 in 1 Hz wind speed in strong convection), driven by the large vertical velocity variances in convective conditions and the high elevation angle of the scanning beams (62∘). Range-gate weighting induces a negative bias into the horizontal wind speeds near the surface shear layer (−0.2 m s−1 in the stable test case). Errors in the horizontal wind speed and direction computed from the wind components are sensitive to the background wind speed but have negligible dependence on the relative orientation of the instrument. Especially during low winds and in the presence of large errors in the horizontal velocity estimates, the reported wind speed is subject to a systematic positive bias (up to 0.4 m s−1 in 1 Hz measurements in strong convection). Vector time-averaged measurements can improve the behavior of the error distributions (reducing the 10 min wind speed error standard deviation to <0.3 m s−1 and the bias to <0.1 m s−1 in strong convection) with a predictable effectiveness related to the number of decorrelated samples in the time window. Hybrid schemes weighting the 10 min scalar- and vector-averaged lidar measurements are shown to be effective at reducing the wind speed biases compared to cup measurements in most of the simulated conditions, with time averages longer than 10 min recommended for best use in some unstable conditions. The approach in decomposing the error mechanisms with the help of the LES flow field could be extended to more complex measurement scenarios and scans.

17 WIND ENERGY↗
DOE OSTI2025

Experimental Comparison of Elastomeric Materials for Hydraulic Seal Durability under Reciprocating Conditions

Wave Energy Converters (WECs) rely on hydraulic Power Take-Off (PTO) systems to transform the kinetic energy of ocean waves into electricity, where hydraulic seals play a critical role in preventing fluid leakage and ensuring watertight operation under harsh marine conditions. However, seal degradation due to wear, fatigue, and corrosion remains a key challenge limiting system reliability and increasing maintenance costs. In this study, the tribological performance of elastomeric materials for hydraulic sealing applications that are based on commercial polyurethane and custom-compounded Ethylene Propylene Diene Monomer (EPDM) is investigated. Using a CSM tribometer, reciprocating wear tests were conducted in accordance with ASTM G133 under dry sliding conditions, with normal loads ranging from 3 to 10 N, and sliding speeds of 10 to 30 mm/s. Coefficient of friction (COF), specific wear rate, and worn surface morphology are evaluated to understand how load, velocity, and material composition affect seal durability. Improved wear stability and high Shore A hardness (93–94) were demonstrated by the commercial thermoplastic polyurethane (Duralast 4203 and 4758). EPDM-based compounds showed different filler-dependent behaviors, formulations filled with carbon black (E3) had the lowest wear rate (1.45 × 10?4 mm³ N?¹ m?¹) and increased load-bearing capacity, while samples modified with silica and plasticizer (E4, E5) had better conformability and damping but less hardness. Compounded EPDMs showed wider deformation zones because of their lower modulus and filler-matrix interfacial fluctuations, while commercial grades exhibited shallower wear scars (= 380 µm), according to optical profilometry. The Archard model-based lifetime estimation showed that, under 10 N and 30 mm/s, operational endurance ranged from 3.1 to 6.2 years. Because of its low wear rate and steady interfacial behavior, E3 had the longest anticipated life. For medium-load hydraulic sealing applications, these results demonstrate the possibility of customized EPDM formulations as affordable substitutes for commercial polyurethanes. They also offer a quantitative foundation for material selection and lifespan prediction.

Durability↗
DOE OSTI2025

Gradient flow based phase-field modeling using separable neural networks

Allen–Cahn equation is a reaction–diffusion equation and is widely used for modeling phase separation. Machine learning methods for solving the Allen–Cahn equation in its strong form suffer from inaccuracies in collocation techniques, errors in computing higher-order spatial derivatives, and the large system size required by the space–time approach. To overcome these challenges, we propose solving the gradient flow of the Ginzburg–Landau free energy functional, which is equivalent to the Allen–Cahn equation, thereby avoiding the second-order spatial derivatives associated with the Allen–Cahn equation. A minimizing movement scheme is employed to solve the gradient flow problem, eliminating the complexities of a space–time approach. We utilize a separable neural network that efficiently represents the phase field through low-rank tensor decomposition. As we use the minimizing movement scheme to numerically solve the gradient flow problem, we thus, refer to the proposed method as the Separable Deep Minimizing Movement (SDMM) method. The evaluation of the functional in the minimizing movement scheme using the Gauss quadrature technique bypasses the inaccuracies associated with collocation techniques traditionally used to solve partial differential equations. A hyperbolic tangent transformation is introduced on the phase field prior to the evaluation of the functional to ensure that it remains strictly bounded within the values of the two phases. For this transformation, theoretical guarantee for energy stability of the minimizing movement scheme is established. Our results suggest that this transformation helps to improve the accuracy and efficiency significantly. The proposed method resolves the challenges faced by state-of-the-art machine learning techniques, outperforming them in both accuracy and efficiency. It is also the first machine learning method to achieve an order of magnitude speed improvement over the finite element method. In addition to its formulation and computational implementation, several case studies illustrate the applicability of the proposed method.

42 ENGINEERING↗
DOE OSTI2024

Arctic Ocean Hydroacoustics

As the Arctic warms and loses its perennial ice cover, it is becoming more attractive for a variety of human uses. Hydroacoustic monitoring of this activity will grow in importance over the coming years and decades. Changes to the physical environment affect acoustic propagation and noise, with ramifications for our ability to detect and locate events and activities of interest. In this report, we use two long-term data sets from the Beaufort Sea, western Arctic Ocean, to determine how acoustic propagation conditions and seismic source detections are impacted by the changing environment: 1) oceanographic observations from ice-tethered profilers, and 2) passive acoustic recordings from a hydrophone. We find that changes to Beaufort Sea thermohaline stratification is stabilizing a subsurface duct, leading to more focused acoustic energy arrivals. Detections of catalogued submarine earthquakes show geographic differences in signal strength between seismic and acoustic waves. Signal strength increases with earthquake magnitude, but relationships to other source and path factors are less clear. Ambient noise also has clear seasonal patterns in the Arctic, with relatively low noise in the spring, higher noise near 1 Hz in summer, and higher noise near 10 Hz in winter. Climate change is expected to modify these seasonal noise patterns, impacting event detection. Future work will further investigate the mechanisms of ice effects on sound and couple acoustic modeling to an Earth System Model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Reducing Mass of Steel Auto Bodies using Thin Advanced High Strength Steel with Carbon-Fiber Reinforced Epoxy

Diversitak, a company based in Detroit, MI, has developed a proprietary, low specific gravity, carbon fiber-reinforced epoxy (CFRE) under U.S. patent number 9,963,58832. Preliminary testing on this new material conducted in collaboration with ArcelorMittal Steel Company proved out the CFRE concept. A thin layer of this CFRE was applied to a stamped sheet of steel with residual stamping oils from a mill, in a time corresponding to automotive processing (e.g., ~15 seconds), and processed following automotive e-coat procedures (phosphating + 175–200°C heating), to complete the curing. No problems with adherence or performance were noted. While the CFRE does add weight to a thin gauge steel panel, it weighs much less than what is displaced by using thicker conventional mild steel gauges. The application of the coating showed a significant increased dent resistance, oil canning resistance, and part stiffness.This current two-year project was designed to mature this new technology to near manufacturing readiness to reduce the weight of a vehicle and lower the cost of weight reduction. The process involves the use of thinner gauge steels than are currently used. The collaborative development team included two industrial manufacturers: Diversitak and ArcelorMittal Steel Company; LightMAT; and two National Laboratories: Oak Ridge National Laboratory (ORNL) and Idaho National Laboratory (INL). The team developed a new manufacturing process to reduce the weight of a vehicle and lower the cost of weight reduction, as well as a better understanding of how to apply the coating so that it will perform to a high standard in-service. The team also performed an in-depth study to determine the long-term durability of the materials manufactured using this technology and well-known automotive industry standard tests.The overall process involved stiffening the thinner gauge steel by applying the CFRE on only one side. To accomplish this goal, the optimal reinforcement fiber length and fiber concentration was first determined. This was followed by measuring the coefficient of thermal expansion (CTE) in all three directions, so it could be fed into manufacturing models and methods for rapidly and inexpensively applying the coating. This was followed by panel level evaluations of the coating and steel combination, and then by full part demonstration of the technology on door panels. The final step was corrosion testing of the parts.ArcelorMittal characterized the advanced high strength steel (AHSS) (e.g., metallurgy-heat treatment for required AHSS properties as a function of the sheet thickness, state of internal stress) and quantified CFRE adhesion to the steel as a function of sheet preparation (e.g., rolling and stamping).ORNL optimized the fiber length, fiber concentration, and coating thickness for best vehicle function and performance at the least cost. Along with the suppliers, ORNL developed a durable CFRE application process (e.g., gun material, design, robotic dispensing process) and identified the adhesion stability of the CFRE during process holding. An approach to ensure that application/curing timing conforms to conventional assembly line speed and plant cycle times was determined. ORNL also determined the CTE of the material in all three directions and performed material scanning electron microscopy (SEM) analyses.INL characterized the corrosion properties of the steel panels coated with CFRE. The panels were investigated for corrosion resistance and stability as replacement materials used in automotive body panels to reduce mass. The coupons tested at INL were supplied by Diversitak after an optimized CFRE formulation was achieved in the already coated form for corrosion testing.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

An Assessment of Deploying Advanced Pumped Storage Hydropower Technology in U.S. Electricity Markets

Pumped storage hydropower (PSH) is a type of hydropower technology where energy can be stored and generated by moving water between two reservoirs of differing elevations. In addition to providing 97% of the total utility-scale hydropower storage in the United States, PSH plants have operational characteristics - such as high ramp rates and the ability to provide reserves - that contribute to greater flexibility and reliability of the power grid. New PSH technologies can provide additional flexibility beyond existing, fixed-speed units. With the emergence of high levels of variable renewable energy resources (e.g., wind and solar), energy storage is expected to be crucial to the reliability and reliance of the power grid in a low-carbon future. At diurnal and longer durations, PSH plants have some of the lowest costs per unit of energy, have been proven to be reliable and efficient, are not cycle-limited, and typically have long lives, often exceeding 50 years. The work presented here focuses on a new generation of PSH: ternary PSH and quaternary PSH (together referred to as T/Q-PSH). Given recent experience in Europe, grid operators in the United States and elsewhere are increasingly focusing their attention on T/Q-PSH as a proven, financeable technology that can offer utility-scale, long-duration, fast-acting energy storage capabilities and grid services. T/Q-PSH differs from conventional PSH in that it can provide fast-response ancillary services during both generating and pumping operations. This capability together with fast mode switching times are key attributes that make T/Q-PSH technologies attractive for managing and stabilizing electricity systems with high amounts of variable renewable energy.

ENERGY PLANNING, POLICY, AND ECONOMY,HYDRO ENERGY↗
DOE OSTI2025

Intrinsic property of defective β-Ga 2 O 3 to self-heal under ionizing irradiation

Damage evolution and phase stability in defective β-Ga 2 O 3 and an irradiation-converted γ-Ga 2 O 3 layer have been studied under ionizing irradiation at 300 K. By exploring athermal nonequilibrium processes in β-Ga 2 O 3 , we succeed in identifying a self-healing mechanism that enables recovery pre-existing damage, characterized by a recovery cross-section of ~0.17 nm 2 . Remarkably, this study further demonstrates that the crystallinity of the irradiation-converted γ-Ga 2 O 3 layer improves under ionizing irradiation. More importantly, X-ray diffraction analysis reveals that the highly-strained 𝛾 -phase transforms into a highly-crystalline structure without film disintegration, contrasting to that reported for isochronal annealing at 1000 K. The inelastic thermal spike calculations provide insights into the important effects of energy transfer to electrons in reordering the local atomic arrangement of both defective β- Ga 2 O 3 and 𝛾-Ga 2 O 3 . This behavior suggests a pathway for low-temperature crystallization, offering a promising strategy for fabricating ultrahigh-speed non-volatile memory devices.

36 MATERIALS SCIENCE↗
DOE OSTI2024

DNS of Flame Stabilization Dynamics of a Swirl-Stabilized Spray Burner Using Sustainable Aviation Fuels

Concern for emission reduction has motivated the development of new cost-effective alternative sustainable aviation fuels (SAFs). Drop-in SAFs with blends of the certified and the alternative fuels are beneficial because they do not require engine modifications for use in current aviation engines. Combustion characteristics are of particular concern when comparing a new fuel with Jet-A. Lean blow-off and flame stabilization have been identified to be governed by the fuel propensity to autoignition, i.e. derived cetane number. However, fundamental investigation is required to identify the impact of derived cetane number on the combustion mode of flame stabilization in a realistic combustor. The combustion mode has also a direct impact on turbulent combustion modeling closures. DNS simulations are performed in the low-Mach solver of the Pele Suite called PeleLMeX. Lagrangian multi-phase modeling is used to capture the liquid spray injection of Jet-A (reference fuel) and C1 as a representative of a low cetane number SAF. Adaptive Mesh Refinement (AMR) is used to enable a more efficient simulation of a more realistic domain size and embedded boundary treatment is used to model a bluff-body geometry. Local extinction and edge flame propagation were observed for both fuels. The edge flame propagation mode was quantified in terms of a Damkohler number defined as the ratio between progress variable reaction rate and its diffusive flux. Initial analysis suggests that a mixed mode combustion occurs for the edge flame propagation, with flame propagation assisted by ignition and autoignition co-existing for both Jet-A and C1 flames. The analysis of the extinction region shows a larger progress variable for Jet-A which can further increase the local displacement speed and present a faster reignition of the stoichiometric mixture fraction.

adaptive mesh refinement↗
DOE OSTI2023

Large eddy simulation of atmospheric boundary layer flow over complex terrain in comparison with RANS simulation and on-site measurements under neutral stability condition

Large eddy simulation (LES) of the atmospheric boundary layer (ABL) flow over complex terrain is presented with a validation using meteorological tower (met-tower) data through an improved neutral stability sampling approach. The proposed stability sampling procedure includes a condition based on the most-likely occurrence time-periods of the neutral ABL and reduces the variabilities of the conditional wind statistics calculated at the met-towers in comparison to our previous work. The ABL flow simulations are carried out over a potential wind site with a prominent hill based using the OpenFOAM-based simulator for on/off-shore wind farm applications by applying the Lagrangian-averaged scale-invariant dynamic sub-grid scale turbulence model. A low-dissipative scale-selective discretization scheme for the non-linear convection term in the LES governing equation is adopted implicitly to ensure both the second-order accuracy and bounded solution. The LES inflow is generated through a precursor method with a “tiling” approach based on the flow driving parameters obtained from a corresponding Reynolds-averaged Navier–Stokes (RANS) simulation. Overall, the averaged wind velocity profiles predicted by the LES approach at all met-tower locations show a similar tendency as the RANS results, which are also in reasonable agreement with the met-tower data. An obvious difference in wind speed standard deviation profiles is seen between LES and RANS, especially at regions downstream of the hill edge, where the LES shows under-predicted results at the highest measurement levels in comparison to the tower data. The computational costs of the LES are found to be about 20 times higher than the RANS simulations.

Energy & Fuels↗
DOE OSTI2025

Stable low diffusion flux splitting schemes on unstructured meshes

Shock instabilities are shown to manifest in modern low-diffusion flux-vector splitting (FVS) schemes when used on unstructured meshes, or situations where shocks do not align with the mesh lines. These instabilities occur irrespective of the Mach number of the shock. Three types of dissipative mechanisms that suppress these instabilities are presented. These mechanisms are carefully designed in order to affect only problematic regions of the flux-splittings. The AUSM + and LDFSS schemes are stabilized using the proposed modifications. It is shown that the added dissipation improves the shock behavior of AUSM and LDFSS on unstructured meshes. It is also shown that the AUSM + -up scheme is prone to the “carbuncle” instability, a specific type of shock instability, when used on unstructured meshes. The modifications proposed in this work do not lead to carbuncle instabilities for the problems considered here. Furthermore, the modified schemes are shown to satisfy certain properties that are crucial for accurate shear layer computations, such as stationary contact preservation. Using benchmark problems, it is demonstrated that despite the diffusion added for stabilization, these schemes are not overly diffusive. Furthermore, due to these advantages, the modified FVS schemes presented here are promising candidates for high-speed compressible flow computations on unstructured meshes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Cause and impact of low-frequency chirping modes in DIII-D hybrid discharges

Significant variations in MHD activity and fast-ion transport are observed in the DIIID high-beta, steady-state hybrid discharges with a mixture of Electron Cyclotron (EC) waves and Neutral Beam Injection (NBI). When Electron Cyclotron Heating (ECH) or Current Drive (ECCD) is applied, Alfven Eigenmodes (AEs) are usually suppressed and replaced by low-frequency bursting modes. The analysis of a recently compiled database of hybrid discharges suggests that the change of the fast-ion pressure especially the perpendicular pressure is the main factor responsible for the instability transition. The lower ratio of fast-ion injection speed v inj to Alfven speed v alfven and slight drop of safety factor q min during ECCD may also facilitate the transition. The database shows that AEs mainly occur when the fast-ion fraction P f /P total is less than 0.53 and v inj /v alfven is greater than 0.50, while low-frequency bursting modes appear in the opposite regime. In this research, P f and P total are the central fast-ion pressure from classical prediction and total plasma pressure, respectively. The reason why the instability transition correlates with P f /P total and v inj /v alfven is that they can significantly modify the drive of low-frequency bursting modes and AEs. The explanation is supported by the observation that low-frequency bursting modes are rarely seen in the hybrids with NBI only, with EC waves and counter-NBI, or with high plasma density. A careful check of the low-frequency bursting modes suggests that they are mainly chirping (neoclassical) tearing modes [referred to as chirping (N)TM], i.e. the mode frequency firstly jumps up from the steady (N)TM frequency, then chirps down, and finally returns to the steady (N)TM frequency. Sometimes, the (N)TMs are fully stabilized and replaced with pure fishbones. Lastly, the resonance condition calculation and ’Kick’ model simulations suggest that (N)TMs and fishbones can interact through modification of the fast ion distribution in phase space, which influences the drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Accurate machine-learning predictions of coercivity in high-performance permanent magnets

Increased demand for high-performance permanent magnets in the electric vehicle and wind-turbine industries has prompted the search for cost-effective alternatives. Discovering magnetic materials with the desired intrinsic and extrinsic permanent magnet properties presents a significant challenge to researchers because of issues with the global supply of rare-earth elements, material stability, and a low maximum magnetic energy product BH max . While first-principles density functional theory (DFT) predicts materials’ magnetic moments, magnetocrystalline anisotropy constants, and exchange interactions, it cannot compute extrinsic properties such as coercivity (H c ). Although it is possible to calculate H c theoretically with micromagnetic simulations, the predicted value is larger than the experiment by almost an order of magnitude due to the Brown paradox. To circumvent these issues, we employ machine-learning (ML) methods on an extensive database obtained from experiments, DFT calculations, and micromagnetic modeling. The use of a large experimental dataset enables realistic H c predictions for materials such as Ce-doped Nd 2 ⁢Fe 14 ⁢B, comparing favorably against micromagnetically simulated coercivities. Remarkably, our ML model accurately identifies uniaxial magneto-crystalline anisotropy as the primary contributor to H c . With DFT calculations, we predict the Nd-site-dependent magnetic anisotropy behavior in Nd 2 ⁢Fe 14 ⁢B, confirming that Nd 4⁢g sites mainly contribute to uniaxial magnetocrystalline anisotropy, and also calculate the Curie temperature (T c ). Finally, both calculated results are in good agreement with the experiments. The coupled experimental dataset and ML modeling with DFT input predict H c with far greater accuracy and speed than was previously possible using micromagnetic modeling. Further, we reverse engineer the grain-boundary and intergrain exchange coupling with micromagnetic simulations by employing the ML predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Combustion synthesis of nanoscale magnesium borides with improved hydrogen uptake and release

The overarching goal of the reported project was to develop reversible metal hydrides that can be dehydrogenated and recharged at relatively low temperatures and pressures. The project complemented the research on reversible metal hydrides for hydrogen storage conducted by the Hydrogen Materials – Advanced Research Consortium (HyMARC). The students involved were trained at both the University of Texas at El Paso (UTEP) and Sandia National Laboratories (SNL). The research objectives were to fabricate nanoscale magnesium boride powders with different B/Mg atomic ratios and to investigate the effects of B/Mg ratio and particle size on the hydrogenation properties of the obtained materials. Magnesium diboride (MgB 2 ) and tetraboride (MgB 4 ) were synthesized by using combustion synthesis and by heating in a tube furnace. The latter method resulted in materials with significantly reduced impurity levels. An effective method for the removal of magnesium oxide impurities was identified. High-energy ball milling produced submicron MgB 2 and MgB 4 powders. However, hydrogenation experiments at pressures up to 700 bar have not shown any uptake of hydrogen. Thermogravimetric analysis of the decomposition of the obtained MgB 2 and MgB 4 demonstrated higher thermodynamic stability of the latter. Thus, the boron-rich nature and phase stability of magnesium borides with high B/Mg ratios make them poor candidates for direct hydrogenation. However, the synthesized MgB 4 is a promising energetic additive to solid fuels for high-speed air-breathing propulsion as it offers a lower oxidation onset temperature and greater energy density than MgB 2 .

08 HYDROGEN↗
DOE OSTI2021

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Advanced Perovskite Solar Cells and Modules

The “Advanced perovskite Cells and Modules” research project was the final agreement focused on enhancing perovskite solar cell (PSC) technologies funded by the US Department of Energy's Solar Energy Technologies Office. The project was designed to address three crucial areas in PSC development: stability, manufacturability, and efficiency. The project was then structured around three main tasks, each targeting one of these strategic goals. The team of experienced researchers in these materials worked collaboratively to address the targets outlined in the technical work plan. building on existing PSC research while also exploring promising new concepts arising in the field. An overview of each primary task is summarized below: Task 1 Stability: This first task, aims to identify material characteristics and metrics that can help predict the primary degradation mechanisms impacting PSC stability. This involved developing specific device tests based on hypotheses regarding mechanisms impacting stability, including fast failure procedures to speed up PSC development and improvement. Various strategies to enhance stability, like incorporating additives, post-treatments, novel contact materials etc. were developed using this fast feedback approach. The relationships between indoor and outdoor stresses were also validated. Task 2 Manufacturability: This second task, focused on creating a scalable production process for PSCs. Initially the objective is to establish a best-known method for a 182 cm2 minimodule. However, given resource limitations, these metrics were modified to focus on the other goal of outlined in the TWP. Specifically, this task worked to demonstrate the transferability of this best-known method to another research institution. Work scope in this area was expanded to material purity and understanding of reagent/process relationships. Examination of other difficulties in PSC production and potential solutions for large-scale production were also evaluated. Given challenges observed in process transfer, work to develop data infrastructure and recording tools for processing of material and devices was then also prioritized in this task. Task 3 Efficiency: This task was focused on improvements to PCE, while still considering Task 1 and Task 2 goal. The efforts targeted a PCE greater than 22% with a T95 exceeding 1000 hours at 25°C in a nitrogen environment for lab-scale devices (approximately 0.1 cm2 devices) across a range of solar-relevant perovskite compositions, including wide-gap (around 1.7 eV) and low-gap (around 1.3 eV) materials, using standard metal contacts. This work then provides a foundation for MHP-based tandem efforts undertaken in other projects and the All-MHP tandem efforts outlined in this projects TWP. Work in this project emphasized disseminating its findings through peer-reviewed publications (PRP), conference presentations, and industrial collaborations. Significant products were produced in all these areas, over 53 peer reviewed publications, 32 conference presentations and industrial investment based on NLR assistance on precompetitive challenges. The team also developed significant intellectual property and awards for their technical excellence, innovations and leadership. The team also leveraged traditional and social media platforms to engage with stakeholders and the public.

14 SOLAR ENERGY↗
DOE OSTI2023

Study of Microphonic Effects on the C100 Cryomodule for High Energy Electron Beam Accelerators

The Continuous Electron Beam Accelerator Facility (CEBAF) at Thomas Jefferson National Laboratory (JLab) is a particle accelerator which can accelerate an electron beam to relativistic speeds and apply the beam onto target samples. The C100 superconducting radio frequency (SRF) cavity is the primary accelerating structure of the C100 cryomodule, one of the many cryomodules which compose the CEBAF linear accelerator. SRF cavities are particularly sensitive to internal and external vibrations that can result in a phenomenon called microphonics which degrade the operational stability of a cryomodule. The purpose of this thesis is to investigate the significance of mechanical disturbances on the electromagnetic resonant frequency of a C100 SRF cavity. Knowledge of the mechanical resonance of the cavities and cryomodule sheds light into how these disturbances are most easily realized as deformation which causes radio frequency (RF) detuning. Three studies were conducted: the development and hammer test calibration of a Finite Element Analysis (FEA) model of a C100 cavity, the development and hammer test calibration of an FEA model of a C100 cavity string, and the hammer test of the C100-10R cryomodule at the Cryomodule Test Facility (CMTF). The cavity FEA model was found to accurately predict two modes found in two real cavities in a simply supported configuration. The cavity string FEA model leveraged the calibrated cavity FEA model but was not found to accurately predict the modal behavior of a real cavity string. Even so, the modal behavior of the cavity string inside the C100-10R cryomodule was captured during a hammer test while it was partially assembled. Finally, the C100-10R cryomodule was placed in the CMTF to study RF detuning. The RF detuning spectra during hammer hits and background noise was captured. The results of the hammer testing indicate two strong peaks at low frequencies (9-10 Hz and 22-23 Hz). These two frequencies were found to be nearly coincident to four instances of mechanical resonance found during the hammer testing done on the partially-assembled C100-10R. Because of this, these two modes are believed to contribute to RF detuning of the cryomodule. This test event also included the testing of the effectiveness of a configuration of BNNT canisters designed to act as dampers. While these tests show promising results, the lurking variables render these tests somewhat inconclusive.

Hull, Caleb James↗
DOE OSTI2022

Experimental and analytical study of the hydrodynamic and single and two-phase convective heat transfer performance of flexible PDMS microchannels with micropillar arrays

Various copper and silicon based thermal management systems are used in the cooling of electronics. However, the rigid nature of these materials along with their high thermal and electrical conductivity pose a difficulty in developing direct contact embedded flexible cooling systems that can offer robust cooling performance. The low density, thermal stability, chemical inertness, and electrical insulation of Polydimethylsiloxane (PDMS) make it an ideal material to develop lightweight direct contact thermal management systems for electronics. Its ease of fabrication with tunable flexibility provides the opportunity to go beyond traditional electronics and develop advanced active and passive thermal management systems for a wide–range of applications in foldable and wearable electronics, liquid cooling garments, microgravity, and electric motors. In this study, a flexible PDMS based microchannel with micropillar arrays, which enhance the thermal performance of the device through capillary-assisted flow, has been developed. The hydrodynamic and convection heat transfer performance of three PDMS wick pillar geometries, ranging from a porosity of 0.8–0.91, are investigated and compared under single-phase and two-phase conditions. Dielectric coolant FC-3283 is employed and permeability measurements are made for mass fluxes ranging from 53 kg/m 2 s to 369 kg/m 2 s. Given its conformability, the device demonstrates a deviation from Darcy’s Law, within the laminar regime, with an increasing permeability with mass flux at the rate of ~0.5–0.8 Darcy/(kg/m 2 s). A semi-analytical model has been developed and reported to quantify the conformability of the device. The heat transfer performance is experimentally evaluated using the same dielectric fluid for mass flux ranging from 105 kg/m 2 s to 420 kg/m 2 s with heat fluxes ranging from 1.5 W/cm 2 to 16 W/cm 2 . Heat transfer coefficients of up to 7000 W/m 2 K are observed, which are comparable to copper and silicon microchannels. The effect of porosity on the single phase thermal performance has been evaluated against the pumping power to provide a basis for thermal management system design. Finally, high-speed imaging is performed to study the two-phase flow characteristics to provide insight into the vapor formation and removal.

42 ENGINEERING↗