Search NASA⌕ Search

SEARCH · Search NASA

Results for “mixed methods”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Analysis of Modeling Assumptions used in Production Cost Models for Renewable Integration Studies

Renewable energy integration studies have been published for many different regions exploring the question of how higher penetration of renewable energy will impact the electric grid. These studies each make assumptions about the systems they are analyzing; however the effect of many of these assumptions has not been yet been examined and published. In this paper we analyze the impact of modeling assumptions in renewable integration studies, including the optimization method used (linear or mixed-integer programming) and the temporal resolution of the dispatch stage (hourly or sub-hourly). We analyze each of these assumptions on a large and a small system and determine the impact of each assumption on key metrics including the total production cost, curtailment of renewables, CO2 emissions, and generator starts and ramps. Additionally, we identified the impact on these metrics if a four-hour ahead commitment step is included before the dispatch step and the impact of retiring generators to reduce the degree to which the system is overbuilt. We find that the largest effect of these assumptions is at the unit level on starts and ramps, particularly for the temporal resolution, and saw a smaller impact at the aggregate level on system costs and emissions. For each fossil fuel generator type we measured the average capacity started, average run-time per start, and average number of ramps. Linear programming results saw up to a 20% difference in number of starts and average run time of traditional generators, and up to a 4% difference in the number of ramps, when compared to mixed-integer programming. Utilizing hourly dispatch instead of sub-hourly dispatch saw no difference in coal or gas CC units for either start metric, while gas CT units had a 5% increase in the number of starts and 2% increase in the average on-time per start. The number of ramps decreased up to 44%. The smallest effect seen was on the CO2 emissions and total production cost, with a 0.8% and 0.9% reduction respectively when using linear programming compared to mixed-integer programming and 0.07% and 0.6% reduction, respectively, in the hourly dispatch compared to sub-hourly dispatch.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

Optimizing district energy systems by integrating Borehole Thermal Energy Storage Using a Mixed-Integer Linear Programming g-function framework with a Multi-Timescale Rolling Horizon method

Shallow geothermal has gained increasing attention in recent years; however, a reliable framework for its accurate incorporation into large-scale energy system optimization remains lacking. This study proposes a Mixed-Integer Linear Programming (MILP) framework combined with the g-function approach to integrate Borehole Thermal Energy Storage (BTES) technology into energy system optimization. Validation against a Modelica-based reservoir network simulation demonstrates that the proposed framework effectively captures the ground thermal response under varying energy loads and accurately estimates the borefield energy supply. To enhance scalability, a Rolling Horizon with Multi-Timescale (RH-MTS) method is further introduced, reducing computational time by 73 % for the 1-year optimization model with only minor loss of optimality. The framework is demonstrated through the case study of the UC Berkeley campus. Results indicate that BTES is a cost-effective and low-carbon solution: two borefields comprising 382 boreholes can meet 8.0 % and 6.6 % of the total campus heating and cooling demand, respectively, at an average energy rate of 0.70–0.77 USD/kWh and carbon intensity of 0.54 kg-CO2/kWh. Short-term analysis reveals a 35%–65% decline in BTES energy flow after 3–6 months of continuous heating/cooling operation, while long-term simulation shows that annual energy production of BTES can vary by up to 12.0 % after four years before stabilizing. Overall, this study develops a novel optimization framework that couples physics-based g-function method with MILP optimization framework, thereby advancing methodological development for shallow-geothermal integration and providing actionable guidance for BTES deployment in district-energy systems.

Yang, Jiahui↗
DOE OSTI2017

Multi-tracer investigation of river and groundwater interactions: a case study in Nalenggele River basin, northwest China

Environmental tracers (such as major ions, stable and radiogenic isotopes, and heat) monitored in natural waters provide valuable information for understanding the processes of river–groundwater interactions in arid areas. An integrated framework is presented for interpreting multi-tracer data (major ions, stable isotopes ({sup 2}H, {sup 18}O), the radioactive isotope {sup 222}Rn, and heat) for delineating the river–groundwater interactions in Nalenggele River basin, northwest China. Qualitative and quantitative analyses were undertaken to estimate the bidirectional water exchange associated with small-scale interactions between groundwater and surface water. Along the river stretch, groundwater and river water exchange readily. From the high mountain zone to the alluvial fan, groundwater discharge to the river is detected by tracer methods and end-member mixing models, but the river has also been identified as a losing river using discharge measurements, i.e. discharge is bidirectional. On the delta-front of the alluvial fan and in the alluvial plain, in the downstream area, the characteristics of total dissolved solids values, {sup 222}Rn concentrations and δ{sup 18}O values in the surface water, and patterns derived from a heat-tracing method, indicate that groundwater discharges into the river. With the environmental tracers, the processes of river–groundwater interaction have been identified in detail for better understanding of overall hydrogeological processes and of the impacts on water allocation policies.

58 GEOSCIENCES↗
DOE OSTI2012

Co{sub 3}O{sub 4}-Co{sub 2}ZnO{sub 4} spinels: The case for a solid solution

In prior first-principles theoretical work we predicted a complete solid solution in the Co{sub 3}O{sub 4}-Co{sub 2}ZnO{sub 4} system, with a negligibly small mixing enthalpy. In this work we tested this prediction on bulk, large-grained specimens across the Co{sub 3}O{sub 4}-Co{sub 2}ZnO{sub 4} join, combining oxide melt solution calorimetry, differential scanning calorimetry, precise lattice parameter measurements, anomalous X-ray and neutron diffraction, and in situ electrical measurements. The calorimetric results confirm the presence of a solid solution at high temperatures, but with a large enthalpy of mixing that exceeds the predicted value. Because Co{sub 3}O{sub 4} and Co{sub 2}ZnO{sub 4} have essentially identical lattice parameters, this energetic destabilization must arise from factors other than the strain energy resulting from size mismatch. Changes in Co{sup 3+} spin states vs. temperature and zinc content are proposed to account for the positive excess enthalpy, and may also provide additional entropy to stabilize the solid solution at high temperature. - Graphical abstract: The first-principles theoretical prediction of a complete solid solution in Co{sub 3}O{sub 4}-Co{sub 2}ZnO{sub 4} spinels has been confirmed experimentally by calorimetric, structural, and electrical measurements. The large mixing enthalpy observed exceeds the predicted value and cannot be accounted for by strain or cation distribution contributions alone, since such effects are small in this system; instead, changes in the spin state of Co{sup 3+} are proposed to account for the mixing thermodynamics. Highlights: Black-Right-Pointing-Pointer Experiments confirm theory prediction of a solid solution in Co{sub 3}O{sub 4}-Co{sub 2}ZnO{sub 4} spinels. Black-Right-Pointing-Pointer Thermodynamics of mixing studied by calorimetric, electrical, structural methods. Black-Right-Pointing-Pointer Cation distributions found to be mostly normal for various compositions. Black-Right-Pointing-Pointer Mixing enthalpy is larger than prediction and inconsistent with strain effects. Black-Right-Pointing-Pointer Changes in Co{sup 3+} spin state proposed as source of large mixing enthalpy/entropy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Virtual Design and Analysis of the Balancing Element of an External Gear Machine Considering Cavitation and Mixed Lubrication Effects

This report presents a simulation based approach for sizing the axial balance elements of a pressure compensated external gear machine. The reference hydraulic unit must be able to guarantee operation in two quadrants, as either a pump or a motor. It also needs to operate in a wide operating range of speed and pressure, as it is to be used as primary hydraulic unit for an electro-hydraulic actuator (EHA). The design procedure builds upon past work at the author’s center, but it extends it to the case of multiple quadrant units. Also, a unique method for modeling of both mixed lubrication and cavitation that might occur in the lubricating film is used. After describing the design procedure, the paper discusses the features of the balancing of a reference unit, along with the effects of both cavitation and mixed lubrication.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2013

Foam patterns

A method of creating a foam pattern comprises mixing a polyol component and an isocyanate component to form a liquid mixture. The method further comprises placing a temporary core having a shape corresponding to a desired internal feature in a cavity of a mold and inserting the mixture into the cavity of the mold so that the mixture surrounds a portion of the temporary core. The method optionally further comprises using supporting pins made of foam to support the core in the mold cavity, with such pins becoming integral part of the pattern material simplifying subsequent processing. The method further comprises waiting for a predetermined time sufficient for a reaction from the mixture to form a foam pattern structure corresponding to the cavity of the mold, wherein the foam pattern structure encloses a portion of the temporary core and removing the temporary core from the pattern independent of chemical leaching.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Near Infrared Spectroscopy Calibration for Wood Chemistry: Which Chemometric Technique Is Best for Prediction and Interpretation?

This paper addresses the precision in factor loadings during partial least squares (PLS) and principal components regression (PCR) of wood chemistry content from near infrared reflectance (NIR) spectra. The precision of the loadings is considered important because these estimates are often utilized to interpret chemometric models or selection of meaningful wavenumbers. Standard laboratory chemistry methods were employed on a mixed genus/species hardwood sample set. PLS and PCR, before and after 1st derivative pretreatment, was utilized for model building and loadings investigation. As demonstrated by others, PLS was found to provide better predictive diagnostics. However, PCR exhibited a more precise estimate of loading peaks which makes PCR better for interpretation. Application of the 1st derivative appeared to assist in improving both PCR and PLS loading precision, but due to the small sample size, the two chemometric methods could not be compared statistically. This work is important because to date most research works have committed to PLS because it yields better predictive performance. But this research suggests there is a tradeoff between better prediction and model interpretation. Future work is needed to compare PLS and PCR for a suite of spectral pretreatment techniques.

09 BIOMASS FUELS↗
DOE OSTI2014

A DNS evaluation of mixing models for transported PDF modelling of turbulent nonpremixed flames

Transported probability density function (TPDF) approaches are well suited to modelling turbulent, reacting, variable density flows. One of the main challenges to the successful deployment of TPDF methods is accurately modelling the unclosed molecular mixing term. This study examines three of the most widely used mixing models: the Interaction by Exchange with the Mean (IEM), Modified Curl (MC) and Euclidean Minimum Spanning Tree (EMST) models. Direct numerical simulation (DNS) data-sets were used to provide both initial conditions and inputs needed over the course of the runs, including the mean flow velocities, mixing frequency, and the turbulent diffusion coefficient. The same chemical mechanism and thermodynamic properties were used, allowing the study to focus on the mixing model. The simulation scenario was a one-dimensional, nonpremixed, turbulent jet flame burning either a syngas or ethylene fuel stream that featured extinction and reignition. This test scenario was selected because extinction and reignition phenomena are sensitive to the mixing model. Three DNS cases were considered for both the syngas and ethylene cases with a parametric variation of Reynolds and Damköhler numbers, respectively. Extinction events became more popular with increasing Reynolds number in the syngas cases and with decreasing Damköhler number in the ethylene cases. The model was first tested with the mixing frequency defined from the dissipation rate and variance of mixture fraction. With this definition, for the syngas cases this study finds that the TPDF method is successful at predicting flame extinction and reignition using all three mixing models for the relatively lower and intermediate Reynolds number cases, but that all models under-predict reignition in the relatively higher Reynolds number case. In the ethylene fuelled cases, only the EMST mixing model correctly predicts the reignition event for the two higher Damköhler number cases, yet, in the lowest Damköhler number case the EMST model over-predicts reignition and the IEM and MC models under-predict it. Mixing frequency was then modelled based on the turbulence frequency and a model constant $C_{\phi}$, the ratio of scalar to mechanical mixing rates. The DNS cases were reexamined with this definition and the findings suggested that the optimal value for $C_{\phi}$ is mixing model and case dependent. In particular, it was found in the ethylene case considered that reignition could be achieved with the IEM and MC models by adjusting the value of $C_{\phi}$.

42 ENGINEERING↗
DOE OSTI2024

Systematic Improvement of Quantum Monte Carlo Calculations in Transition Metal Oxides: sCI-Driven Wavefunction Optimization for Reliable Band Gap Prediction

Accurate determination of the electronic properties of correlated oxides remains a significant challenge for computational theory. Traditional Hubbard-corrected density functional theory (DFT+U) frequently encounters limitations in precisely capturing electron correlation, particularly in predicting band gaps. We introduce a systematic methodology to enhance the accuracy of diffusion Monte Carlo (DMC) simulations for both ground and excited states, focusing on LiCoO 2 as a case study. By employing a selected configuration interaction (sCI) approach, we demonstrate the capability to optimize wavefunctions beyond the constraints of single-reference DFT+U trial wavefunctions. Here, we show that the sCI framework enables accurate prediction of band gaps in LiCoO 2 , closely aligning with experimental values and substantially improving traditional computational methods. The study uncovers a nuanced mixed state of t 2g and e g orbitals at the band edges that is not captured by conventional single-reference methods, further elucidating the limitations of PBE+U in describing d-d excitations. Our findings advocate for the adoption of beyond-DFT methodologies, such as sCI, to capture the essential physics of excited-state wavefunctions in strongly correlated materials. The improved accuracy in band gap predictions and the ability to generate more reliable trial wavefunctions for DMC calculations underscore the potential of this approach for broader applications in the study of correlated oxides. This work not only provides a pathway for more accurate simulations of electronic structures in complex materials but also suggests a framework for future investigations of the excited states of other challenging systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Similarity of the turbulent kinetic energy dissipation rate distribution in the upper mixed layer of the tropical Indian Ocean

Turbulence within the upper ocean mixed layer plays a key role in various physical, biological, and chemical processes. Between September and November 2011, a dataset of 570 vertical profiles of the turbulent kinetic energy (TKE) dissipation rate, as well as conventional hydrological and meteorological data, were collected in the upper layer of the tropical Indian Ocean. These data were used to statistically analyze the vertical distribution of the TKE dissipation rate in the mixed layer. The arithmetic-mean method made the statistical TKE dissipation rate profile more scattered than the median and geometric-mean methods. The statistical TKE dissipation rate were respectively scaled by the surface buoyancy flux and the TKE dissipation rate at the mixed-layer base. It was found that the TKE dissipation rate scaled by that at the mixed-layer base exhibited better similarity characteristics than that scaled by the surface buoyancy flux, whether the stability parameter D/|L MO | was greater or less than 10, indicating that the TKE dissipation rate at the mixed-layer base is a better characteristic scaling parameter for reflecting the intrinsic structure of the TKE dissipation rate in the mixed layer, where D and L MO are respectively the mixed-layer thickness and the Monin-Obukhov length scale. The parameterization of the TKE dissipation rate at the mixed-layer base on the shear-driven dissipation rate and the surface buoyancy flux was further explored. It was found that the TKE dissipation rate at the mixed-layer base could be well fitted by a linear combination of three terms: the wind-shear-driven dissipation rate, the surface buoyancy flux, and a simple nonlinear coupling term of these two.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1991

Thermal postbuckling of thin-walled composite stiffeners

A study is made of the thermal postbuckling response of composite stiffeners subjected to prescribed edge displacement and a temperature rise. The flanges and web of the stiffeners are modeled by using two-dimensional plate finite elements. A mixed formulation is used with the fundamental unknowns consisting of the generalized displacements and the stress resultants of the plate. A reduction method is used in conjunction with mixed finite element models for determining the postbuckling response of the stiffeners. Sensitivity derivatives are evaluated and used to study the effects of variations in the different lamination and material parameters of the stiffeners on their postbuckling response characteristics. Numerical studies are presented for anisotropic stiffeners with Zee and channel sections.

Noor, Ahmed K.↗
DOE OSTI2012

An efficient modeling method for thermal stratification simulation in a BWR suppression pool

The suppression pool in a BWR plant not only is the major heat sink within the containment system, but also provides major emergency cooling water for the reactor core. In several accident scenarios, such as LOCA and extended station blackout, thermal stratification tends to form in the pool after the initial rapid venting stage. Accurately predicting the pool stratification phenomenon is important because it affects the peak containment pressure; and the pool temperature distribution also affects the NPSHa (Available Net Positive Suction Head) and therefore the performance of the pump which draws cooling water back to the core. Current safety analysis codes use 0-D lumped parameter methods to calculate the energy and mass balance in the pool and therefore have large uncertainty in prediction of scenarios in which stratification and mixing are important. While 3-D CFD methods can be used to analyze realistic 3D configurations, these methods normally require very fine grid resolution to resolve thin substructures such as jets and wall boundaries, therefore long simulation time. For mixing in stably stratified large enclosures, the BMIX++ code has been developed to implement a highly efficient analysis method for stratification where the ambient fluid volume is represented by 1-D transient partial differential equations and substructures such as free or wall jets are modeled with 1-D integral models. This allows very large reductions in computational effort compared to 3-D CFD modeling. The POOLEX experiments at Finland, which was designed to study phenomena relevant to Nordic design BWR suppression pool including thermal stratification and mixing, are used for validation. GOTHIC lumped parameter models are used to obtain boundary conditions for BMIX++ code and CFD simulations. Comparison between the BMIX++, GOTHIC, and CFD calculations against the POOLEX experimental data is discussed in detail.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Pushered single shell implosions for mix and radiation trapping studies using high-Z layers on National Ignition Facility

Pushered Single Shells (PSSs) are another approach to Inertial Confinement Fusion implosions that employ high-Z materials in the innermost capsule layer (pusher) as a means to enhance radiation trapping and lower core ignition requirements. However, adding high-Z materials can also increase losses due to mix, provide extra tamping, and make the capsule emission opaque to x-ray diagnostics. The first PSS implosions performed on the National Ignition Facility use plastic ablators with a germanium (Ge) dopant as a high-Z surrogate in the pusher to isolate the effects of high-Z mix and radiation trapping without changing tamping. Using a 2-shock laser pulse, the PSS implosions are designed and symmetrized to reach 3.7 keV core temperatures. A low concentration (2.8%) Ge dopant is added to the innermost layer, and the resulting effects on mix and x-ray opacity are observed. The method of separated reactants is used to infer information about mixing between the deuterated plastic pusher and the capsule fill gas (25% tritium) from the resulting nuclear DT reactions. Radiation transport is studied via capsule emission x-ray spectroscopy and imaging. Both nuclear and x-ray data corroborate the hypothesis that the addition of Ge strongly affects the mix region through radiation losses but has a minimal effect on the core and the warm, unmixed regions. Simulations using diffusive and turbulent mix models agree qualitatively with data, but quantitative agreement may require hybrid mix models that can model the transitional regime between turbulence and diffusion. Simulations matching the observables show increased core radiation trapping when Ge is added.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Transient cycling of nitrogen, organic carbon and oxygen within the free-flowing Columbia River corridor: Linking exposure time dependent biogeochemical reactions to river stage fluctuations (Final project report)

The objective of this project was to develop next-generation techniques for representing the transformations of complex reaction networks like those impacting transient river-corridors, then apply them to develop upscaling tools. The project contributed significantly to the development of Lagrangian “mass transfer particle tracking” (MTPT) tools that are the first numerical framework that explicitly separates mixing and spreading processes; a crucial distinction for accurate representations of reaction rates across scales. The MTPT approach uses an operator splitting scheme where physical transport processes are simulated using classical random walk methods and each particle is treated as a “container” that carries with it any number of chemical species. Mixing between particles is simulated using a colocation-probability based mass transfer kernel and reactions are evaluated on each particle after mixing. The method is stable, accurate, and also allows the explicit modeling of residence time distributions of the individual containers of mass. The MTPT scheme allowed us to explore several research questions related to reactive transport in river corridors and complex reaction networks. We found that it is possible to quantify the exposure time distributions (ETDs) of reactants and that this does suggest some pathways to upscaling. However, comparable approximations can be obtained using simplified (reduced dimensionality) MTPT simulations in less time with similar uncertainty, so development of extensive ETD-based methods was not productive. The major findings are that i) the MTPT schemes are robust and highly accurate across a range of arbitrarily complex reactions, ii) the methods can be efficiently parallelized and the parallel performance characteristics of MTPT are predictable, and iii) the MTPT tools allow simultaneous tracking of residence time. The development of these research tools into efficient software packages continues. The methods are already available to the community because open-source, working examples have been included with all publications.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1972

Initial development of a hypersonic free mixing layer

A preliminary experimental investigation to establish some of the characteristics and further the understanding of the initial development of a turbulent free mixing layer for hypersonic speeds has been conducted. Mean profile data at about 6 inches downstream of the exit of a hypersonic nozzle have been obtained in nitrogen for a nominal Mach number of 19.5, total temperature of about 1670 K and Reynolds number range from about 50,000 to 110,000 per foot and have been compared with profiles upstream of the nozzle exit. Static pressure varied across the shear layer for the present tests. The outer 80 percent of the high-velocity portion of the free shear layer can be calculated by a rotational method of characteristics. However, turbulent mixing is evidently important in the low-velocity region, and effects of eddy viscosity and eddy conductivity should be included in a theoretical analysis.

Harvey, W. D.↗
NASA NTRS1989

Advances in ambient temperature secondary lithium cells

The Jet Propulsion Laboratory is involved in a Research and Development program sponsored by NASA/OAST on the development of ambient temperature secondary lithium cells for future space applications. Some of the projected applications are planetary spacecraft, planetary rovers, and astronaut equipment. The main objective is to develop secondary lithium cells with greater than 100 Wh/kg specific energy while delivering 1000 cycles at 50 percent Depth of Discharge (DOD). To realize these ambitious goals, the work was initially focused on several important basic issues related to the cell chemistry, selection of cathode materials and electrolytes, and component development. The performance potential of Li-TiS2, Li-MoS3, Li-V6O13 and Li-NbSe3 electrochemical systems was examined. Among these four, the Li-TiS2 system was found to be the most promising system in terms of realizable specific energy and cycle life. Some of the major advancements made so far in the development of Li-TiS2 cells are in the areas of cathode processing technology, mixed solvent electrolytes, and cell assembly. Methods were developed for the fabrication of large size high performance TiS2 cathodes. Among the various electrolytes examined, 1.5M LiAsF6/EC + 2-MeTHF mixed solvent electrolyte was found to be more stable towards lithium. Experimental cells activated with this electrolyte exhibited more than 300 cycles at 100 percent Depth of Discharge. Work is in progress in other areas such as selection of lithium alloys as candidate anode materials, optimization of cell design, and development of 5 Ah cells. The advances made at the Jet Propulsion Laboratory on the development of secondary lithium cells are summarized.

Subbarao, S.↗