Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational approach”

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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Apollo 16 core 60013/14 as a product of path 1 and path 2 regolith evolution processes

We have made a petrographic study of 12 samples along the 62 cm profile of Apollo 16 core 60013/14 taken at Station 10 near the lunar module. The basic core description and maturity variations have been discussed earlier. This core is important because it enables a stratigraphic comparison with nearby cores 60001-7 and 60009/10. These cores form a triangle with sides 35-40 m. FMR data show a maturity-match between 60009/10 and 60013/14. Unfortunately, segments of the core 60001-7 are too disturbed for any cm scale comparison. We present new data on the grain size distribution and the modal petrography for six samples from six levels of 60013 using the same approach and rationale used for 60014. We have combined the new data with our previous data on 60014 to plot the variations of grain size distributions and grain type abundances, in the 500-1000 m and 90-150 micrometers size fractions, with depth. Agglutinate abundances match the I sub s/FeO profile and both show a significant break at a depth of about 45 cm, with lower maturity material below that point. Grain size distributions show that for all analyzed samples, the mean grain size is below 90 micrometers. The two lowest samples are somewhat more coarse-grained than the rest of the core. The break at about 45 cm is the most prominent feature of the core and was described during the initial dissection as a pronounced change in color from dark to light. It now seems clear that this break is primarily one of maturity rather than of bulk chemical composition.

Basu, A.↗
NASA NTRS1974

A unified approach for the application of general perturbation theories to the artificial satellite problem

An Encke-type method is developed as well as a variation of parameters method, both of which use a non-Keplerian reference orbit. The regularized time is used in the numerical integration and the optimum value of n is found for each type of orbit investigated. Accuracy and computation time comparisons are made with a classical Cowell method. It should be noted that the element formulations developed and tested were found to be exceptionally numerically stable in the sense that it was possible to achieve numerical consistency order of very long integration periods - a property not available with the Cowell formulation - and therefore these formulations may be helpful for high precision calculations.

Alfriend, K. T.↗
DOE OSTI2023

Quasi-Steady Simulation of Glaze Ice Accretion and Heat Transfer in the Supercooled Large Droplet Regime in Atmospheric Turbulent Air Flow

Abstract To avoid the use of a computationally intensive full unsteady simulation while providing an accurate solution, a quasi-steady simulation has been performed to study glaze and mixed ice accretion in the supercooled large droplet (SLD) regime in turbulent flow of atmospheric air. We have attempted to find the minimum time-step necessary to adequately simulate the icing process on an airfoil surface. Based on node displacement, a mesh morphing scheme has been adopted in the computations to account for the moving boundaries that are caused by the continuous ice buildup. We have modeled ice accretion on an airfoil surface for a time period of 232 s using several time steps. At each time-step we solved the steady-state conservation equations for the air and droplet phases and then used the results as initial conditions for the next time position. The magnitude of the time-step ranged from a least accurate two-shot simulation (where the time-step is 116 s) to a most accurate 2320-shot simulation (where the time-step is 0.1 s). In-between these two extreme time steps, we have performed a three-shot simulation (where the time-step is 77.33 s), a four-shot simulation (where the time-step is 58 s), a six-shot simulation (where the time-step is 38.67 s), a 46-shot simulation (where the time-step is 5 s), and a 232-shot simulation (where the time-step is 1 s). We have done so to find out the degree of accuracy (or inaccuracy) of the multishot simulation approach and to find out the appropriate time-step needed for a successful and valid quasi-steady simulation. A valid quasi-steady simulation needs to use a time-step that is small enough to reproduce the full time-dependent solution within a very small error band. We have found that both the 1 s and the 0.1 s time steps produce virtually identical results. This is the primary litmus test that proves the validity of the quasi-steady-state assumption. The results adopted in this paper are thus all based on the more conservative 0.1 s time-step. In the process of performing the simulations, remeshing was required in order to maintain the grid density in zones of high curvature to be able to capture the full physics in those zones. After successfully modeling glaze ice accretion over the airfoil surface using the 0.1 s quasi-steady simulation approach, the effects of supercooled large droplets (SLDs) impacting the surface have been examined and presented in terms of the variation of the local collection efficiency, the water film thickness, and the heat transfer rate. Examination of the variation of the angle that the ice horn makes with the airfoil chord line demonstrated a 20% improvement in angle prediction when the time-step is reduced from 116 s to 0.1 s. The analysis also reveals a 12% or 8.5% increase in the maximum collection efficiency, βmax, depending on whether the value of the liquid water content (LWC) has been doubled from 0.5 g/m3 to 1 g/m3, or the value of the freestream velocity has been doubled from 75 m/s to 150 m/s, respectively. Because of the need to monitor the local collection efficiency and convective heat fluxes at each shot (regardless of the number of shots employed), the approach adopted here was found to be effective in successively and successfully reproducing the curvature of the glaze ice horn.

Engineering↗
DOE OSTI2021

Improving Resilience of Bus Bunching Holding Strategy through a Rolling Horizon Approach

Providing public transportation with quality service is critical to attracting more passengers to the system. However, high-demand routes are prone to the so-called bus bunching -- a tendency of buses to group as a consequence of variations in travel times and demands. Bus holding is applied to overcome this effect. In this study, we present a novel method for bus holding in which the control law is based only on the buses' position using a computationally efficient rolling horizon approach. The method uses similar inputs as linear control approaches while not increasing significantly the computational time. However, the method overcomes key a weakness of the linear control approach thanks to the explicit constraint handling that always ensures the control action effectiveness. Simulation experiments in a validation case and a model-specific for a bus rapid transit line in Curitiba, Brazil showed a reduced holding time and improved resilience, delivering more than 20% reduction in delay time accounting for the on-board and station delays.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Extrapolation of the Rainflow-Counted Load Ranges for Fatigue Assessment of the Wind Turbine's Blades

Wind turbine design standards recommend the use of statistical modeling coupled with extrapolation of the short-term load data to long-term periods for fatigue reliability assessment. However, statistical error and computational expense can limit the accuracy of such approaches. In the case of wind turbine blades, the errors are more significant because of the high material fatigue exponent that makes the damage estimations more sensitive to variations. In addition, due to different excitation sources, the flapwise load range histogram is not unimodal, and thus its statistical modeling is complex. In the present work, we provide three methods for statistical modeling of the flapwise bending moment ranges including a novel approach based on frequency-based separation of the modes. The first two methods are simplified approaches for modeling the most crucial load ranges using unimodal distributions and the third method involves multimodal distribution fitting. The research is based on 3600 10-minute aeroelastic simulations of DTU 10MW case study wind turbine from which a benchmark damage equivalent load (DEL) is calculated. The DEL calculated by each of the three proposed methods is compared to this reference. The results show that the conventional approach based on using 6 seeds as well as using mixture models fitted on the limited data lead to under-conservative results with errors up to 23%. On the other hand, the simplified unimodal approaches provided in this work can provide conservative estimations of the fatigue damage with mean values 5% and 12% higher than the benchmark. However, the variability of the DEL estimates is higher when using unimodal extrapolation of the load ranges, and the data can be conservative by 17.5%. The proposed unimodal fits suggested for modeling and extrapolation of the blade's load ranges provide less errors relatively and most importantly conservative DEL estimations while maintaining computational efficiency.

blade fatigue↗
DOE OSTI2023

Adaptive cold-load pickup considerations in 2-stage microgrid unit commitment for enhancing microgrid resilience

In an extended main grid outage spanning multiple days, load shedding serves as a critical mechanism for islanded microgrids to maintain essential power and energy reserves that are indispensable for fulfilling reliability and resiliency mandates. However, using load shedding for such purposes leads to increasing occurrence of cold load pickup (CLPU) events. Here, this study presents an innovative adaptive CLPU model that introduces a method for determining and incorporating parameters related to CLPU power and energy requirements into a two-stage microgrid unit commitment (MGUC) algorithm. In contrast to the traditional fixed-CLPU-curve approach, this model calculates CLPU duration, power, and energy demands by considering outage durations and ambient temperature variations within the MGUC process. By integrating the adaptive CLPU model into the MGUC problem formulation, it allows for the optimal allocation of energy resources throughout the entire scheduling horizon to fulfill the CLPU requirements when scheduling multiple CLPU events. The performance of the enhanced MGUC algorithm considering CLPU needs is assessed using actual load and photovoltaic (PV) data. Simulation results demonstrate significant improvements in dispatch optimality evaluated by the amount of load served, customer comfort, energy storage operation, and adherence to energy schedules. These enhancements collectively contribute to reliable and resilient microgrid operation.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

The inclusion of uncertainty in circularity transition modeling: A case study on wind turbine blade end-of-life management

The transition to a more circular economy (CE) is complex and hard to predict, including in sustainable energy technologies. Many sources of uncertainty make it challenging to model CE scenarios and their potential benefits. As an example, the high variability in costs and revenues of different recycling options makes future wind turbine blade recycling highly uncertain. To better understand this challenge, the circular economy, life cycle assessment and visualization (CELAVI) framework - a discrete event simulation and life cycle assessment framework - is modified to incorporate uncertainty analysis capabilities. Moreover, a 3-step procedure that covers different aspects of uncertainty in CE studies and includes a Monte-Carlo analysis is proposed. The procedure is tested in a case study on wind turbine blade recycling using CELAVI. Results highlight that grinding and landfilling costs are the most influential parameters for wind power circularity. The model's output coefficients of variation (when input parameter uncertainties are propagated) are between 92% and 384% depending on the indicator. The approach developed in this study may help researchers and decision-makers who study circularity prioritize their data collection effort. Finally, our method contributes to a mounting yet critical body of research: the measurement of uncertainty in circularity transitions.

17 WIND ENERGY↗
DOE OSTI2025

Towards large-scale quantum optimization solvers with few qubits

Quantum computers hold the promise of more efficient combinatorial optimization solvers, which could be game-changing for a broad range of applications. However, a bottleneck for materializing such advantages is that, in order to challenge classical algorithms in practice, mainstream approaches require a number of qubits prohibitively large for near-term hardware. Here we introduce a variational solver for MaxCut problems over $m={{\mathcal{O}}}({n}^{k})$ binary variables using only n qubits, with tunable k > 1. The number of parameters and circuit depth display mild linear and sublinear scalings in m , respectively. Moreover, we analytically prove that the specific qubit-efficient encoding brings in a super-polynomial mitigation of barren plateaus as a built-in feature. Altogether, this leads to high quantum-solver performances. For instance, for m = 7000, numerical simulations produce solutions competitive in quality with state-of-the-art classical solvers. In turn, for m = 2000, experiments with n = 17 trapped-ion qubits feature MaxCut approximation ratios estimated to be beyond the hardness threshold 0.941. Our findings offer an interesting heuristics for quantum-inspired solvers as well as a promising route towards solving commercially-relevant problems on near-term quantum devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Development of MOSCATO: A CFD-Level Electrochemistry and Corrosion Simulator for Molten Salt Systems

For both coolant and fueled variants of molten salt reactors (MSRs), the corrosion of structural materials is a significant challenge. The corrosion stems from chemical and electrochemical reactions initiated by fissile material, fission products, and impurities in the salt. Lower-fidelity models rely on empirical correlations for mass transfer, simplified lumped temperature profiles, and similar assumptions. They do not capture detailed spatial variations in complex geometries, creating the need for high-fidelity modeling to bridge this gap.As we approach the demonstration and possible deployment of MSRs in this decade, the development of a high-fidelity, high-performance simulator becomes imperative. To simulate the complex electrochemical environment and corrosion within molten salt systems, we have developed the Molten Salt Chemistry And TranspOrt (MOSCATO) code. This endeavor is comprised of three essential components. First, mass transfer equations are coupled with the Navier-Stokes equations in order to account for the transport of species in the salt. Second, the diffusion of alloy constituents, such as Cr, Fe, Ni, etc. is simulated within the structural metals. Third, the alloy and salt domains are coupled to account for the heterogeneous chemical and electrochemical reactions that occur at the salt-alloy interface.MOSCATO manages all three components within the framework of the highly scalable, open-source spectral element method computational fluid dynamics code Nek5000/NekRS. This integration enables MOSCATO to harness the immense computational power of modern high-performance computing resources, ensuring both high fidelity and computational speed.In addition to code development, we have initiated a comprehensive verification and validation campaign, utilizing data from diverse sources. First, MOSCATO's electrochemical solver was verified with reference numerical data. Then validation occurred against experiments: one of a thermal galvanic cell and the other for corrosion in flowing molten salt of FLiNaK (LiF-NaF-KF). This campaign verified and validated MOSCATO as a reliable tool for simulating electrochemical environments and corrosion in molten salt systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Dual-cone variational calculation of the two-electron reduced density matrix

The computation of strongly correlated quantum systems is challenging because of its potentially exponential scaling in the number of electron configurations. Variational calculation of the two-electron reduced density matrix (2-RDM) without the many-electron wave function exploits the pairwise nature of the electronic Coulomb interaction to compute a lower bound on the ground-state energy with polynomial computational scaling. Recently, a dual-cone formulation of the variational 2-RDM calculation was shown to generate the ground-state energy, albeit not the 2-RDM, at a substantially reduced computational cost, especially for higher $\textit{N}$-representability conditions such as the T2 constraint. Here we generalize the dual-cone variational 2-RDM method to compute not only the ground-state energy but also the 2-RDM. The central result is that we can compute the 2-RDM from a generalization of the Hellmann-Feynman theorem. Specifically, we prove that in the Lagrangian formulation of the dual-cone optimization the 2-RDM is the Lagrange multiplier. Further, we apply the method to computing the energies and properties of strongly correlated electrons—including atomic charges, electron densities, dipole moments, and orbital occupations—in an illustrative hydrogen chain and the nitrogen-fixation catalyst FeMoco. The dual variational computation of the 2-RDM with T2 or higher N -representability conditions provides a polynomially scaling approach to strongly correlated molecules and materials with significant applications in atomic and molecular and condensed-matter chemistry and physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Bottomonium suppression in 5.02 and 8.16 TeV 𝑝-Pb collisions

Here, we compute the suppression of ϒ⁡(1⁢𝑆), ϒ⁡(2⁢𝑆), and ϒ⁡(3⁢𝑆) states in 𝑝-Pb collisions relative to 𝑝⁢𝑝 collisions, including nuclear parton distribution function (nPDF) effects, coherent energy loss, momentum broadening, and final-state interactions in the quark-gluon plasma. We employ the EPPS21 nPDFs and calculate the uncertainty resulting from variation over the associated error sets. To compute coherent energy loss and momentum broadening, we follow the approach of Arleo, Peigne, and collaborators. The 3+1⁢D viscous hydrodynamical background evolution of the quark-gluon plasma is generated by anisotropic hydrodynamics. The in-medium suppression of bottomonium in the quark-gluon plasma is computed using a next-to-leading-order open quantum system framework formulated within potential nonrelativistic quantum chromodynamics. We find that inclusion of all these effects provides a reasonable description of experimental data from the ALICE, ATLAS, CMS, and LHCb Collaborations for the suppression of ϒ⁡(1⁢𝑆), ϒ⁡(2⁢𝑆), and ϒ⁡(3⁢𝑆) as a function of both transverse momentum and rapidity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A Hybrid Fuel Cell and Battery Storage Power Management for Grid-Interactive EV Charging Station

With the increasing adoption of renewable energy sources in grid-interactive Electric Vehicle (EV) charging stations, the role of energy storage systems has become critical. While large energy storage systems have mitigated the intermittency of renewable energy, integrating multi-source energy management with prioritized charging can further enhance the reliability of charging stations (CS). This paper presents a decentralized energy management (DEM) approach combining battery energy storage (BES) and fuel cell (FC) systems using a rule-based line resistance correction droop (LRCD) control technique. The proposed droop control dynamically adjusts the gain to balance the state-of-charge (SoC) of the BES, enhancing power support longevity and improving battery life under varying capacity conditions by reducing current stress. Additionally, the paper addresses the challenges of using fuel cells in linear regions to optimize efficiency and manage various charging scenarios. The CS integrates unity power factor grid interaction, and power support for auxiliary loads, maintaining harmonic distortion within 5% during grid islanding. The approach evaluates DC bus voltage regulation under various scenarios of PV array power fluctuations and dynamic load variations, in both grid-connected and standalone operations. In conclusion, the proposed control strategy is validated on a laboratory prototype through various dynamic load variation and grid islanding scenarios.

Khalid, Mohd [Oak Ridge National Laboratory (ORNL)↗
DOE OSTI2022

MetaPop: a pipeline for macro- and microdiversity analyses and visualization of microbial and viral metagenome-derived populations

Abstract Background Microbes and their viruses are hidden engines driving Earth’s ecosystems from the oceans and soils to humans and bioreactors. Though gene marker approaches can now be complemented by genome-resolved studies of inter-(macrodiversity) and intra-(microdiversity) population variation, analytical tools to do so remain scattered or under-developed. Results Here, we introduce MetaPop, an open-source bioinformatic pipeline that provides a single interface to analyze and visualize microbial and viral community metagenomes at both the macro - and microdiversity levels. Macrodiversity estimates include population abundances and α- and β-diversity. Microdiversity calculations include identification of single nucleotide polymorphisms, novel codon-constrained linkage of SNPs, nucleotide diversity ( π and θ ), and selective pressures (pN/pS and Tajima’s D ) within and fixation indices ( F ST ) between populations. MetaPop will also identify genes with distinct codon usage. Following rigorous validation, we applied MetaPop to the gut viromes of autistic children that underwent fecal microbiota transfers and their neurotypical peers. The macrodiversity results confirmed our prior findings for viral populations (microbial shotgun metagenomes were not available) that diversity did not significantly differ between autistic and neurotypical children. However, by also quantifying microdiversity, MetaPop revealed lower average viral nucleotide diversity ( π ) in autistic children. Analysis of the percentage of genomes detected under positive selection was also lower among autistic children, suggesting that higher viral π in neurotypical children may be beneficial because it allows populations to better “bet hedge” in changing environments. Further, comparisons of microdiversity pre- and post-FMT in autistic children revealed that the delivery FMT method (oral versus rectal) may influence viral activity and engraftment of microdiverse viral populations, with children who received their FMT rectally having higher microdiversity post-FMT. Overall, these results show that analyses at the macro level alone can miss important biological differences. Conclusions These findings suggest that standardized population and genetic variation analyses will be invaluable for maximizing biological inference, and MetaPop provides a convenient tool package to explore the dual impact of macro - and microdiversity across microbial communities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Interactions of Microhabitat and Time Control Grassland Bacterial and Fungal Composition

Dryland grasslands are vast and globally important and, as in all terrestrial ecosystems, soil microbial communities play fundamental roles in regulating dryland ecosystem function. A typical characteristic of drylands is the spatial mosaic of vascular plant cover surrounded by interspace soils, where biological soil crusts (biocrusts)—a complex community of organisms including bacteria, fungi, algae, mosses, and lichens—are common. The implications of this heterogeneity, where plants and biocrust cover co-occur, are often explored in the context of soil fertility and hydrology, but rarely has the impact of these multiple microhabitat types been simultaneously explored to determine the influence on bacterial and fungal communities, key biological players in these ecosystems. Further, our understanding of the temporal dynamics of bacterial and fungal communities in grasslands, and of how these dynamics depend on the microhabitat within the ecosystem, is notably poor. Here we used a temporally and spatially explicit approach to assess bacterial and fungal communities in a grassland on the Colorado Plateau, and to link variation in these communities to edaphic characteristics. We found that microhabitat (e.g., vascular plant rhizosphere, biocrust, and below biocrust) was the strongest driver of differences in bacterial and fungal community richness, diversity, and composition. Microhabitat type also significantly mediated the impact of temporal change in shaping community composition. Taken together, 29% of the variation in bacterial community composition could be explained by microhabitat, date, and microhabitat-by-date interactions, while only 11% of the variation in fungal community composition could be explained by the same factors, suggesting important differences in community assembly processes. Soil microbial communities dictate myriad critical ecosystem functions, thus understanding the factors that control their compostition is crucial to considering and forecasting how terrestrial ecosystems work. Overall, this case study provides insights for future studies on the spatial and temporal dynamics of bacterial and fungal communities in dryland grasslands.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Thermomechanical Processing for Improved Mechanical Properties of HT9 Steels

Thermomechanical processing (TMP) of ferritic–martensitic (FM) steels, such as HT9 (Fe–12Cr–1MoWV) steels, involves normalizing, quenching, and tempering to create a microstructure of fine ferritic/martensitic laths with carbide precipitates. HT9 steels are used in fast reactor core components due to their high-temperature strength and resistance to irradiation damage. However, traditional TMP methods for these steels often result in performance limitations under irradiation, including embrittlement at low temperatures (<~430 °C), insufficient strength and toughness at higher temperatures (>500 °C), and void swelling after high-dose irradiation (>200 dpa). This research aimed to enhance both fracture toughness and strength at high temperatures by creating a quenched and tempered martensitic structure with ultrafine laths and precipitates through rapid quenching and unconventional tempering. Mechanical testing revealed significant variations in strength and fracture toughness depending on the processing route, particularly the tempering conditions. Tailored TMP approaches, combining rapid quenching with limited tempering, elevated strength to levels comparable to nano-oxide strengthened ferritic alloys while preserving fracture toughness. For optimal properties in high-Cr steels for future reactor applications, this study recommends a modified tempering treatment, i.e., post-quench annealing at 500 °C or 600 °C for 1 h, possibly followed by a brief tempering at a slightly higher temperature.

36 MATERIALS SCIENCE↗
NASA NTRS1972

The detection and mapping of water vapor in the Martian atmosphere.

The objectives of the water vapor investigations, to be carried out during the Viking missions, are described in the light of our current knowledge of Martian atmospheric conditions and the diurnal and seasonal variations of the water abundance. A discussion is given of the relative merits of the different experimental approaches which can be adopted for this purpose, followed by a summary of the detection limits achievable in the available spectral regions. These considerations led to the choice of a spectrometer operating in the 1.4-micron region, with a spectral resolution of approximately one wavenumber and a surface spatial resolution at periapsis of 3 by 24 km. The design of the spectrometer and the results of laboratory simulation tests which have been carried out to verify the expected instrumental performance are described.

Farmer, C. B.↗
NASA NTRS1972

Stability of orthotropic stiffened composite plates.

It is pointed out that the buckling problem of a plate with restrictions resulting from boundary conditions is suited for an application of Lagrangian multipliers. An approach used by Budiansky and Hu (1946) is extended to the study of orthotropic stiffened plates. The variations of the total energy with respect to the three deformation components are considered. The equilibrium equations and the boundary conditions are obtained. It is shown that the deflection functions must satisfy the equilibrium equations. The Lagrangian multipliers are then used to solve stability problems on the basis of the total energy involved, the equilibrium equations, and the boundary conditions. The theoretical analysis is supplemented by an experimental investigation.

Chiu, K. D.↗
NASA NTRS1974

A study of light airplane pilot landing performance

The results of a total of 289 landings performed with a low-wing aircraft by a group of private pilots on long and short runways have been analyzed to determine the landing performance of these pilots. Both the long- and short-runway landings show significant variation from straight or uniformly curved trajectories, denoting considerable 'jockeying' on the part of the pilot during approach. The long-runway approach speeds were variable and higher than recommended. Considerable float followed by touchdown at speeds well above stall were noted. The short-runway landings were quite similar, except that approach speeds were slightly lower as a result of the use of two-thirds to full flaps, and touchdown occurred closer to the threshold.

Obryan, T. C.↗