SEARCH · Engineering Papers
Results for “parabolic”
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.
Design and characterization of a single-channel microwave interferometer for the Helicon Physics Prototype eXperiment
A heterodyne interferometer for measuring the plasma electron density in Helicon Physics Prototype eXperiment (HPPX) has been developed and characterized in the laboratory. The heterodyne-type interferometer uses a probing beam of 85 GHz and a reference beam of 85.7 GHz by using two Gunn oscillators for facilitating frequency modulation. The cutoff density for this frequency is 8.96 x 10 19 m –3 that is much higher than the nominal electron density of the HPPX plasmas. The maximum deflection angle due to the refraction effect is estimated to be 0.63° assuming a parabolic electron density profile, and the impact on the measurement is negligible for the plasmas with a size of 100 mm. Two identical HDPE lenses act as the vacuum windows, and the lenses working together with the transmitter and receiver horns, form a beam with the beam waist of about 20 mm located at the plasma center. The characterization results of the lens and the critical microwave components show good performance. The plasma electron density measurement results are in line with the expectation. Furthermore, the Signal to Noise Ratio (SNR) of the measurement is defined as SNR = 20lg($\bar{n}_{e}$/S), where $\bar{n}_{e}$ and S are respectively the average electron density and the standard deviation of the signal, and the SNR is 80 dB. Finally, this demonstrates the system's capability of providing a clear and detailed density resolution.
Alluvial fan morphology: A self-similar free boundary problem description
In this work, we examine approximate geometrically self-similar solutions to a parabolic free boundary value problem applied to alluvial fan surface morphology and growth. Alluvial fans are fan- or cone-shaped sedimentary deposits caused by the rapid deposition of sediment from a canyon discharging onto a flatter plain. Longitudinal, topographic profiles of fans can be readily described by a seemingly time independent dimensionless profile (DeChant et al., 1999). However, because an alluvial fan can be expected to grow over time, it is clear that this “steady” profile is certainly time dependent and can be described using a space-time self-similar solution. In an experimental and theory-based study, Guerit et al. (2014) developed a self-similar (or as they describe it a self-affine) linear solution based upon an approximate first order small parameter expansion solution for a 1-d homogeneous nonlinear diffusion equation. Direct substitution of this result into a linear diffusion equation suggests that this first order expression may not fully satisfy the associated governing equation. In contrast, we develop a more complete solution based upon a modeled approximation for the axi-symmetric formulation such that the associated temporal behavior is consistent with a 1/3 time power-law as described by Reitz and Jerolmack (2014). The resulting expression is an exact solution to a linear heat equation. Furthermore, we emphasize that a small parameter is not inherent to the resulting profile result and is not included in our model development. Though developed using rather different approaches, the formal solution developed here is in good agreement with the simple polynomial described by DeChant et al. (1999) suggesting that this self-similar solution is a suitable time dependent representation of alluvial fan longitudinal profile form and improves on earlier work.
Multi-arrival infrasound from meteoroids: Fragmentation signatures versus propagation effects in a fine-scale layered atmosphere
Infrasonic signatures of meteoroid fragmentation are frequently ambiguous: do multiple arrivals signify a complex breakup or merely the distorting effects of a layered atmosphere? Resolving this ambiguity is critical for accurate energy estimates and source reconstruction. In this study, we address this challenge by analyzing a unique regional dataset of well-constrained meteoroid events observed by the Southern Ontario Meteor Network and the co-located Elginfield Infrasound Array. We employ pseudo-differential parabolic equation (PPE) simulations to quantify how fine-scale gravity-wave structures in the stratosphere and lower thermosphere modify acoustic waveforms at ranges <300 km. Our modeling reveals that while fine-scale layering can stretch signals and generate diffuse oscillatory tails, it does not produce discrete, high-amplitude pulse splitting at ranges below ∼140 km. By applying these results to the rare multi-arrival event 20060305, we demonstrate that its distinct double arrival at 100 km range is inconsistent with atmospheric multipathing and provides definitive evidence of separate fragmentation episodes. These findings establish new diagnostic criteria for separating source physics from propagation artifacts, improving the reliability of infrasound as a monitoring tool for natural bolides, space debris re-entries, and catastrophic launch failures.
Gas nitriding behavior of refractory metals and implications for multi-principal element alloy design
Multi-principal element alloys (MPEAs) comprise a large, flexible compositional space that enables tuning of their chemistry, structure, and properties. Here, to facilitate the development of nitriding-based surface-enhancement strategies that harness a broad compositional space, this study examined the gas nitriding behavior of Hf, Mo, Nb, Ta, Ti, and Zr as a function of time, temperature (750 and 1000 °C), and nitriding potential (i.e. ammonia-to-hydrogen ratio). These metals were selected because they have a strong driving force to form nitrides, and appear in many promising refractory MPEA compositions. The nitriding temperatures were selected based on the phase transformation temperature of Ti and Zr, and the nitriding potentials were chosen such that all elements are expected to form nitrides. Mass gain measurements indicate that all six elements follow parabolic kinetics. The microstructure observations and quantitative microchemical analysis show formation of dense and well-adhered compound layers for Mo, Nb, and Ta. Thick diffusion zones appear in Hf, Ta, Ti, and Zr, and diffusion coefficients were fit to the composition profiles. Partial delamination of the compound layer occurred for Ti and Zr. Peak hardness values above 30 GPa are obtained in the dense compound layers, and the solute hardening of the underlying alloy is correlated with the nitrogen content. The results provide insight into the dynamics of nitride compound formation relative to interstitial dissolution of nitrogen, and are discussed in the context of MPEA composition and processing design.
A linearity preserving nodal variation limiting algorithm for continuous Galerkin discretization of ideal MHD equations
A stabilized continuous Galerkin (CG) method for magnetohydrodynamics (MHD) is presented herein. Ideal, compressible inviscid MHD equations are discretized in space on unstructured meshes using piecewise linear or bilinear finite element bases to get a semi-discrete scheme. Stabilization is then introduced to the semi-discrete method in a strategy that follows the algebraic flux correction paradigm. This involves adding some artificial diffusion to the high order, semi-discrete method and mass lumping in the time derivative term. The result is a low order method that provides local extremum diminishing properties for hyperbolic systems. The difference between the low order method and the high order method is scaled element-wise using a limiter and added to the low order scheme. The limiter is solution dependent and computed via an iterative linearity preserving nodal variation limiting strategy. The stabilization also involves an optional consistent background high order dissipation that reduces phase errors. The resulting stabilized scheme is a semi-discrete method that can be applied to inviscid shock MHD problems and may be even extended to resistive and viscous MHD problems. To satisfy the divergence free constraint of the MHD equations, we add parabolic divergence cleaning to the system. Various time integration methods can be used to discretize the scheme in time. We demonstrate the robustness of the scheme by solving several shock MHD problems.
A unified formulation of splitting-based implicit time integration schemes
Splitting-based time integration approaches such as fractional step, alternating direction implicit, operator splitting, and locally one dimensional methods partition the system of interest into components, and solve individual components implicitly in a cost-effective way. Here this work proposes a unified formulation of splitting time integration schemes in the framework of general-structure additive Runge–Kutta (GARK) methods. Specifically, we develop implicit-implicit (IMIM) GARK schemes, provide the order conditions for this class, and explain their application to partitioned systems of ordinary differential equations. We show that classical splitting methods belong to the IMIM GARK family, and therefore can be studied in this unified framework. New IMIM-GARK splitting methods are developed and tested using parabolic systems.
Neural operator prediction of linear instability waves in high-speed boundary layers
Here, we investigate if neural operators can predict the linear evolution of instability waves in high-speed boundary layers. To this end, we extend the design of the DeepOnet to ensure accurate and robust predictions, and also to perform data assimilation. In particular, we train DeepONet to take as inputs an upstream disturbance and a downstream location of interest, and to provide as output the perturbation field downstream in the boundary layer. DeepONet thus approximates the linearized and parabolized Navier-Stokes operator for this flow. For successful application to the high-speed boundary layer problem, we add sample weighting and Fourier input features to the regular DeepONet formulation. Once trained, the DeepOnet can perform fast and accurate predictions of the downstream disturbances within the range of training frequencies (inside the distribution). In addition, we show that DeepONet can solve the inverse problem, where downstream wall measurements are adopted as input, and a trained network can predict the upstream disturbances that led to these observations. This capability, along with the forward predictions, allows us to perform a full data assimilation cycle efficiently: starting from wall-pressure data, we predict the upstream disturbance using the inverse DeepONet and its evolution using the forward DeepONet. Finally, we introduce three new metrics to benchmark the training, evaluation and break-even cost of neural operators.
Osmotic control of the spacing of parallel shear cracks in shale growing subcritically in geologic past
The geological genesis of natural cracks in sedimentary rocks such as shale is a problem that needs to be understood to improve the technology of hydraulic fracturing as well as deep sequestration of harmful fluids. Why are the vertical natural cracks roughly parallel and equidistant, and why is the spacing roughly 10 cm rather than 1 cm or 100 cm? Fracture mechanics of critical cracks cannot answer this question. Neither can the material heterogeneity. The growth of critical parallel cracks is impossible because the relative crack face displacements would immediately localize into one crack, leading to an earthquake. The cracks must have formed, on the tectonic time scale, by a slow growth of subcritical shear cracks governed by the Charles-Evans law. The idea advanced here is that what controls the crack spacing is the balance between the reduction, due to shear dilatancy, of the concentration of ions such as Na + and Cl - in each fracture process zone (PFZ), which decelerates the cracks, and the restoration of ion concentration by diffusion of ions from the space between the cracks into the FPZ. This diffusion of water is driven mainly by the osmotic pressure gradient, which offsets the deceleration and depends strongly on the crack spacing. A simple analytical solution of the steady state is rendered possible by approximating the ion concentration profiles between adjacent cracks by parabolic arcs. Applying this theory to Woodford shale yields the approximate crack spacing of 10 cm, which is realistic. Furthermore, the stability of unlimited parallel mode II frictional crack growth is proven by examining the second variation of the free energy. Water concentration drop in the FPZ due to shear dilatancy and its restoration by water diffusion from the inter-crack space have similar effect, although probably much weaker.
Steam oxidation of chromium corrosion barrier coatings for sic-based accident tolerant fuel cladding
SiC-based materials are currently being considered for accident tolerant fuel (ATF) cladding, however, there are concerns regarding hydrothermal corrosion rates and hermeticity. These issues may be addressed with the utilization of a Cr corrosion barrier coating (CBC) but the interaction of Cr with SiC under loss of coolant accident (LOCA) conditions needs to be investigated. In the present work, Cr coatings were deposited onto chemically vapor deposited (CVD) SiC with either high power impulse magnetron sputtering (HiPIMS), cathodic arc (CA), or a combination of both. Annealing and steam thermogravimetric analysis (TGA) was used to study reaction product formation with and without oxidation. Additionally, after 1200 °C annealing, formation of a Cr x C y carbide layer at the ambient interface and a Cr x Si y C z silicide layer at the SiC interface was observed. After 1200 °C TGA exposure, similar carbide and silicide layers were observed but with a Cr 2 O 3 outer reaction layer and a Si and O rich interfacial layer forming between the carbide and silicide layers. Initial mass gain of the coatings during the 1200 °C TGA exposure was parabolic and were similar to reported rates of chromia formation. Mass gain behavior and observed transitions in chromia microstructure indicate that metallic Cr may be consumed within the first ~1 h at 1200 °C and that subsequent chromia formation occurs at the expense of the underlying Cr carbide layer. Overall, the results show no detrimental impact of Cr coatings on the steam oxidation resistance of SiC up to 4 h at 1200 °C.
Grain growth kinetics of the gamma phase metallic uranium
We report metallic uranium is a leading fuel form for sodium cooled fast reactors as an enabling technology of future nuclear energy systems. Mechanistic understanding of fuel behaviors and kinetics under thermodynamic equilibrium and highly non-equilibrium conditions are essential for evaluating fuel performance. It is important to understand and predict the grain and pore evolutions of metallic fuels under thermal and irradiation conditions. However, very limited data are available on the grain growth kinetics and mechanisms of pure gamma phase uranium. In this paper, the pure gamma uranium pellets with different grain structures were fabricated by combining high-energy ball milling and spark plasma sintering. Isothermal annealing tests were performed to investigate the grain growth behavior of the pure gamma phase uranium with different initial grain sizes. A parabolic relationship in grain growth with time was identified for the submicron-sized (374 nm) sample. In contrast, for the nano-sized (137 nm) sample, the grain growth shows a linear relationship with time. The activation energies of grain growth were determined as 199.5 KJ/mol and 80.6 KJ/mol for nano-sized and submicron-sized grain structures, respectively. For the nano-sized sample, the rate-control step of grain growth is dominated by the triple-junction migration, in which the grain boundary triple junction drags the grain growth, leading to a higher activation energy than the bulk diffusion. The dominating mechanism for the submicron-sized sample is grain boundary diffusion. The mechanistic understanding and critical data obtained on the kinetics of pure uranium phases will be useful to evaluate fuel behavior under thermodynamic equilibrium conditions and develop a high fidelity model to predict fuel performance.
Comparison of unirradiated and irradiated AGR-2 TRISO fuel particle oxidation response
The silicon carbide (SiC) coating in a tristructural isotropic (TRISO) particle acts as a barrier to fission product release during reactor operation and accident scenarios. Oxidation and subsequent failure of the SiC layer during a rare air ingress event is a proposed mechanism for fission product release in a high-temperature gas-cooled reactor (HTGR). Although previous oxidation studies have analyzed unirradiated TRISO particle response, this study compared the oxidation behavior of irradiated and unirradiated TRISO particles from the second Advanced Gas Reactor Fuel Development and Qualification Program irradiation experiment (AGR-2). Particles with exposed SiC were subjected to six varying oxidizing tests in the Furnace for Irradiated TRISO Testing (FITT), examined for failure fraction with the Irradiated Microsphere Gamma Analyzer (IMGA) and characterized with focused ion beam and scanning/transmission electron microscopy techniques to analyze the oxide layer. Uncorrelated unirradiated particle failures throughout the series of exposures suggests that external factors inherent to the experiment increased particle failure sensitivity. However, irradiated particle observations indicated an increased failure response at 400 h 1400 °C in both 2% and 21% O 2 atmospheres above failure associated with external factors. Oxide thickness measurements after 400 h at 1400 °C revealed a greater oxidation rate than predicted by parabolic growth, which was attributed to the increased complexity of the oxide structure at longer exposure times. Altering the atmosphere from 21% to 2% O 2 reduced the average oxide thickness by approximately 12%–14% in both irradiated and unirradiated particles at 400 h 1400 °C. Altogether, the minor variations observed between irradiated and unirradiated particles in this study led to the conclusion that unirradiated TRISO particles can be used to approximate irradiated TRISO oxidation kinetics.
Speciation mapping of the oxidation layer on aged uranium dioxide using scanning transmission x-ray spectromicroscopy
In this study, UO 2 was aged in humid air and prepared as a thin section using a focused ion beam (FIB) instrument. The specimen was measured using synchrotron radiation spectromicroscopy techniques at the Beamline 11.0.2 STXM end station of the Advanced Light Source (ALS). Non-negative matrix factorization (NMF) methods were used to identify and map three component x-ray absorption near-edge structure (XANES) spectra in the oxygen K-edge data, revealing a surface layer of U 4 O 9 with a thickness of 206 ± 21 nm, and the bulk of the sample remaining as UO 2 . Uranium N 4,5 -edge XANES spectromicroscopy supports these results. The diffusion-controlled parabolic rate constant for UO 2 oxidation to U 4 O 9 was calculated from the observed layer thickness and compared to literature values. Complementary transmission electron microscopy (TEM) was used to image the sample and identify the phases present in various regions, confirming the STXM results.
Improved correlations of the fuel/cladding liquid penetration rate with the out-of-pile transient database
Here, utilizing data from the Out-of-Pile Transient Database (OPTD), we developed advanced correlations for the cladding liquid penetration rate in metallic fuels used in sodium-cooled fast reactors (SFRs). Based on the available OPTD data, our research highlights the key role of Pu concentration in influencing the liquid penetration behavior, showing different threshold fuel-cladding interaction temperatures, with burnup playing a substantial role, particularly in fuels with a 19% Pu concentration. The penetration rate exhibits a parabolic, diffusion-limited behavior, thereby facilitating predictions of maximum liquid penetration depth over time. These correlations are significantly more accurate, under a variety of conditions, compared to existing single-parameter models, and enhance our understanding of SFR metallic fuel behavior, promoting safer and more efficient reactor operations.
Assessment of the polygeneration approach in wastewater treatment plants for enhanced energy efficiency and green hydrogen/ammonia production
Wastewater treatment plants (WWTPs) offer opportunities to optimize resource utilization and enhance energy efficiency. Here, this study provides a comprehensive analysis of using the polygeneration approach in WWTPs to reduce grid energy dependence, optimize energy distribution, and utilize surplus energy for hydrogen (H 2 ) and ammonia (NH 3 ) production. Several models were employed, including photovoltaic (PV) cells, parabolic trough collectors (PTCs), steam methane reforming, and polymer electrolyte membranes, to assess the feasibility of this approach. Three scenarios were evaluated and compared: Scenario 1 (Baseline) represents the current situation, Scenario 2 maximizes the Net Present Value (NPV), and Scenario 3 minimizes NH 3 production costs. Real data from As-Samra WWTP in Jordan was used to accurately assess the feasibility of each scenario. The results show that Scenario 2 offers the highest profitability and efficiency, with a NPV of 87.48 million USD and an annual NH 3 production of 15,417 tons, reducing both grid dependency and biogas fuel consumption. Both Scenarios 2 and 3 demonstrate the ability to meet thermal demands efficiently while generating significant revenue from NH 3 production. Scenario 3, in particular, achieves competitive H 2 and NH 3 production costs. Environmentally, Scenario 2 significantly reduces annual greenhouse gas emissions by 12.66 kilotons of CO 2eq , with near-zero carbon intensity for thermal energy due to solar reliance. In conclusion, the polygeneration approach offers a promising pathway for WWTPs to achieve greater sustainability, economic gains, and reduced environmental impact, providing valuable insights for decision-makers.
Minimizing grid energy consumption in wastewater treatment plants: Towards green energy solutions, water sustainability, and cleaner environment
Wastewater treatment plants (WWTPs) consume significant amount of energy to sustain their operation. From this point, the current study aims to enhance the capacity of these facilities to meet their energy needs by integrating renewable energy sources. The study focused on the investigation of two primary solar energy systems in As Samra WWTP in Jordan. The first system combines parabolic trough collectors (PTCs) with thermal energy storage (TES). This system primarily serves to fulfill the thermal energy demands of the plant by reducing the demands from boiler units, which allows more biogas for electricity generation. The second system is a photovoltaic (PV) system with Lithium-Ion batteries, which directly produces electricity that will be used to cover part of the electrical energy demands of plant. To assess the optimal configuration, two distinct scenarios have been formulated and compared to the current case scenario (SC#1). The first scenario focuses on maximizing the net present value (NPV) and minimizing the levelized cost of electricity (LCOE). The second scenario is centred on minimizing the levelized cost of heat (LCOH). The findings indicate that both scenarios succeeded in reducing the reliance on the grid to a value that reach 1 %. Moreover, they both reduced biogas percentage in energy production from 88 % to approximately 65 % through the integration of the PV system. In terms of thermal demand, SC#2 reduced the reliance on biogas boiler units from 100 % to 25 %, while SC#3 achieved an even more impressive reduction to just 8 %. The best LCOE value was attained in SC#2, at 0.0895 USD/kWh, with an NPV of 10.54 million USD. Conversely, SC# 3 yielded an LCOH value of 0.0432 USD/kWh th compared to 0.0534 USD/kWh th USD for SC#2. In conclusion, despite their relatively high capital and operating costs, SC#2 and SC#3 managed to substantially decrease the annual electricity expenditure from approximately 2 million USD to 86,000 USD and 0 USD, respectively.
A polarimetry-based field-deployable non-interruptive mirror soiling detection method
The soiling level of heliostat mirrors in Concentrated Solar Power (CSP) fields is one of the key factors that significantly influences optical efficiency. State-of-the-art methods of monitoring heliostats soiling levels still face various challenges, including slow speed, labor-intensive operations, resolution and accuracy constraints or interruptions to solar field operations. Here, we present a rapid, cost-effective, and non-intrusive method for mirror soiling detection based on polarimetric imaging, referred to as Polarimetric Imaging-based Mirror Soiling (PIMS). The compact PIMS device is designed for integration with unmanned aerial vehicles (UAVs), enabling rapid, large-area assessments of heliostat mirrors for efficient soiling detection. Our method utilizes the correlation between the Degree of Linear Polarization (DoLP) and surface soiling level based on Mie scattering theory and Monte Carlo simulations. Field deployment of the PIMS method requires minimal device installation, and its UAV-based operation allows for soiling detection without interrupting plant activities. The PIMS method holds the potential for mirror soiling detection across various concentrated solar power (CSP) plants and can be further adapted for other types of solar fields, such as parabolic trough systems.
Wind Conditions for Solar Energy Facility Design Based on Field Measurements at Nevada Solar One
This study characterises near-surface wind conditions at the Nevada Solar One (NSO) concentrated solar power plant using two years of measurements from a 15-m mast located 30 m west of the solar arrays. The surrounding terrain is fairly flat within 3 km, but nearby solar photovoltaic (PV) and concentrated solar power (CSP) installations introduce strong directional wind inhomogeneity, which motivates a sector-based analysis that distinguishes undisturbed inflow from waked flow. The analysis uses data from three sonic anemometers to assess whether turbulence in the surface layer can be represented as Gaussian and whether standard surface-layer scaling remains valid below 10 m height. The roughness length and displacement height are estimated for twelve directional sectors under near-neutral wind conditions, and all three velocity components show departures from Gaussian behaviour. Spectral analysis indicates that the turbulence contains low-frequency fluctuations with amplitudes higher than those predicted by standard models such as the Kaimal model or the Simiu and Scanlan model. The study introduces a new turbulence spectral model based on the eddy-surface-layer framework to better represent near-surface wind turbulence. The analysis further quantifies the influence of the surrounding power plants on mean wind profiles, one-point spectra, and vertical coherence for all velocity components. The solar collector arrays located 30 m east of the mast significantly modify the turbulence structure, whereas PV installations situated a few hundred metres away have only a limited influence on the turbulence measured at the mast. These findings highlight the limitations of conventional turbulence models for near-surface wind over solar infrastructure and provide improved parameterisations to support the structural design. The implications for the structural response of parabolic troughs and PV modules to such turbulence remain uncertain and require further investigation.