Search NASA⌕ Search

SEARCH · Search NASA

Results for “temperature algorithm”

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 577 records · Page 32

A Machine Learning–Based Approach to Time-series Wave Identification in the Solar Wind

The Wind spacecraft has yielded several decades of high-resolution magnetic field data, a large fraction of which displays small-scale structures. In particular, the solar wind is full of wavelike fluctuations that appear in both the field magnitude and its components. The nature of these fluctuations can be tied to the properties of other structures in the solar wind, such as shocks, that have implications for the time evolution of the solar wind. As such, having a large collection of wave events would facilitate further study of the effects that these fluctuations have on solar wind evolution. Given the large volume of magnetic field data available, machine learning is the most practical approach to classifying the myriad small-scale structures observed. To this end, a subset of Wind data is labeled and used as a training set for a multibranch 1D convolutional neural network aimed at classifying circularly polarized wave modes. Using this algorithm, a preliminary statistical study of 1 yr of data is performed, yielding about 300,000 wave intervals out of about 5,000,000 solar wind intervals. The wave intervals come about more often in the fast solar wind and at higher temperatures, and the number of waves per day is highly periodic. This machine learning–based approach to wave detection has the potential to be a powerful, inexpensive way to catalog waves throughout decades of spacecraft data.

79 ASTRONOMY AND ASTROPHYSICS↗

Enhancing Thomson scattering polychromator performance with multi-pass spectral filters

In photon-deficient, noncollective Thomson scattering diagnostics, filter polychromators are typically employed in the spectral analysis of Thomson-scattered signals to achieve acceptable signal-to-noise performance. Currently, the most common polychromator filter configuration employs a set of single-passband optical filters that define individual spectral channels. Here, we introduce a new spectral analysis method for Thomson scattering based on spectral filters with multiple passbands, referred to as Thomson scattering spectral multiplexing. Implementing multi-bandpass spectral filters on polychromators increases the achievable range of electron temperature measurement for a given number of filters employed. In addition, Thomson scattering spectral multiplexing reduces systematic measurement uncertainty, with fewer required spectral channels, thereby decreasing light loss from reduced optical element interactions. A multi-bandpass filter set, optimized by a genetic algorithm, has been successfully installed and tested on the Helically Symmetric eXperiment (HSX), demonstrating the benefits of the Thomson scattering spectral multiplexing method.

Instruments & Instrumentation↗

Resilience Assessment Framework For Electric Distribution Systems Performance Under Extreme Conditions

The devastating impact of extreme weather-related events is increasingly evident on power grids, especially on distribution grids. The severity of their potential impact calls for 1) developing a suitable resilience assessment framework to capture the system performance and 2) assessing relevant mitigative strategies to lessen the impact of such events. This paper proposes a framework to identify grid vulnerabilities using the energy-at-risk concept to select, disconnect and isolate grid portions due to a resilience event. The proposed framework mainly consists of two steps; i) processing the utility's available infrastructure, i.e., a network model, possible switching combinations, and outage information for those combinations, as a graph-based database, and ii) implementing a novel optimal switching algorithm leveraging database and grid simulated metrics. These switching actions are generated to implement load curtailment in a rolling manner during anticipated grid scarcity conditions. In this study, a test case is created using two taxonomy feeders and is simulated against an extreme temperature event, e.g., a long, relatively cold, and prolonged freeze peak, thereby creating stress on the grid. It is demonstrated that the proposed framework allows utilities to predict the energy-at-risk during such resilience events and design suitable outage management strategies.

Poudel, Shiva↗

Optimal Estimation Retrievals and Their Uncertainties: What Every Atmospheric Scientist Should Know

Remote sensing instruments are heavily used to provide observations for both the operational and research communities. These sensors do not provide direct observations of the desired atmospheric variables, but instead, retrieval algorithms are necessary to convert the indirect observations into the variable of interest. It is critical to be aware of the underlying assumptions made by many retrieval algorithms, including that the retrieval problem is often ill posed and that there are various sources of uncertainty that need to be treated properly. In short, the retrieval challenge is to invert a set of noisy observations to obtain estimates of atmospheric quantities. The problem is often complicated by imperfect forward models, by imperfect prior knowledge, and by the existence of nonunique solutions. Optimal estimation (OE) is a widely used physical retrieval method that combines measurements, prior information, and the corresponding uncertainties based on Bayes’s theorem to find an optimal solution for the atmospheric state. Furthermore, OE also allows the relative contributions of the different sources of error to the uncertainty in the final retrieved atmospheric state to be understood. Here, we provide a novel Python library to illustrate the use of OE for inverse problems in the atmospheric sciences. We introduce two example problems: how to retrieve drop size distribution parameters from radar observations and how to retrieve the temperature profile from ground-based microwave sensors. Using these examples, we discuss common pitfalls, how the various error sources impact the retrieval, and how the quality of the retrieval results can be quantified.

54 ENVIRONMENTAL SCIENCES↗

Real-time chiral dynamics at finite temperature from quantum simulation

In this study, we explore the real-time dynamics of the chiral magnetic effect (CME) at a finite temperature in the (1+1)-dimensional QED, the massive Schwinger model. By introducing a chiral chemical potential μ 5 through a quench process, we drive the system out of equilibrium and analyze the induced vector currents and their evolution over time. The Hamiltonian is modified to include the time-dependent chiral chemical potential, thus allowing the investigation of the CME within a quantum computing framework. We employ the quantum imaginary time evolution (QITE) algorithm to study the thermal states, and utilize the Suzuki-Trotter decomposition for the real-time evolution. This study provides insights into the quantum simulation capabilities for modeling the CME and offers a pathway for studying chiral dynamics in low-dimensional quantum field theories.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Composited structure of non-precipitating shallow cumulus clouds

The normalized distributions of thermodynamic and dynamical variables both within and outside shallow clouds are investigated through a composite algorithm using large-eddy simulations of oceanic and continental cases. The normalized magnitude is maximum near the cloud centre and decreases outwards. While relative humidity (RH) and cloud liquid water ( q l ) decrease smoothly to match the environment, the vertical velocity, virtual potential temperature (θ v ), and potential temperature (θ) perturbations have more complicated behaviour towards the cloud boundary. Below the inversion layer, θ' v becomes negative before the vertical velocity has turned from an updraft to a subsiding shell outside the cloud, indicating the presence of a transition zone where the updraft is negatively buoyant. Due to the downdraft outside the cloud and enhanced horizontal turbulent mixing across the edge, the normalized turbulent kinetic energy (TKE) and horizontal turbulent kinetic energy (HTKE) decrease more slowly from the cloud centre outwards than the thermodynamic variables. The distributions all present asymmetric structures in response to the vertical wind shear, with more negatively buoyant air, stronger downdrafts, and larger TKE on the downshear side. We discuss several implications of the distributions for theoretical models and parameterizations. Positive buoyancy near the cloud base is mostly due to the virtual effect of water vapour, emphasizing the role of moisture in triggering. The mean vertical velocity is found to be approximately half the maximum vertical velocity within each cloud, providing a constraint to achieve possible power-law distributions for some models. Finally, the normalized distributions for different variables are used to estimate the vertical heat and moisture fluxes within clouds. The results suggest that distributions near the cloud edge and variability of maximum perturbations need careful treatment. The fluxes are underestimated in the inversion layer because cloud-top downdrafts cannot be captured well.

54 ENVIRONMENTAL SCIENCES↗

Performance investigation of a high-field active magnetic regenerator

Regenerative magnetic cycles are of interest for small-scale, high-efficiency cryogen liquefiers; however, commercially relevant performance has yet to be demonstrated. To develop improved engineering prototypes, an efficient modeling tool is required to screen the multi-parameter design space. In this work, we describe an Active Magnetic Regenerative Refrigerator (AMRR) prototype using a high-field superconducting magnet that produces a 100 K temperature span. Using the experimental data, a semi-analytic AMR element model is validated and enhanced system performance is simulated using liquid propane as a heat transfer fluid. In addition, the regenerator composition and fluid flow are simultaneously optimized using a differential evolution algorithm. Simulation results indicate that a natural gas liquefier with a 160 K span and a second-law efficiency exceeding 20% is achievable.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Physics-based adaptivity of a spectral method for the Vlasov–Poisson equations based on the asymmetrically-weighted Hermite expansion in velocity space

We propose a spectral method for the 1D-1V Vlasov–Poisson system where the discretization in velocity space is based on asymmetrically-weighted Hermite functions, dynamically adapted via a scaling α and shifting u of the velocity variable. Specifically, at each time instant an adaptivity criterion selects new values of α and u based on the numerical solution of the discrete Vlasov–Poisson system obtained at that time step. Once the new values of the Hermite parameters α and u are fixed, the Hermite expansion is updated and the discrete system is further evolved for the next time step. The procedure is applied iteratively over the desired temporal interval. The key aspects of the adaptive algorithm are: the map between approximation spaces associated with different values of the Hermite parameters that preserves total mass, momentum and energy; and the adaptivity criterion to update α and u based on physics considerations relating the Hermite parameters to the average velocity and temperature of each plasma species. For the discretization of the spatial coordinate, we rely on Fourier functions and use the implicit midpoint rule for time stepping. The resulting numerical method possesses intrinsically the property of fluid-kinetic coupling, where the low-order terms of the expansion are akin to the fluid moments of a macroscopic description of the plasma, while kinetic physics is retained by adding more spectral terms. Moreover, the scheme features conservation of total mass, momentum and energy associated in the discrete, for periodic boundary conditions. A set of numerical experiments confirms that the adaptive method outperforms the non-adaptive one in terms of accuracy and stability of the numerical solution.

97 MATHEMATICS AND COMPUTING↗

Novel ammonia recovery from anaerobic digestion by integrating biogas stripping and gypsum absorption

Recovering ammonia from anaerobic digestion (AD) effluent offers the opportunity to produce a renewablefertilizer and reduce environmental impact. A novel process developed for such an application was evaluatedduring this study. The feature of the process is integrating biogas stripping and gypsum (CaSO 4 •2H 2 O) absorption. The result revealed that the lower CO 2 content and the absorption temperature contributed the higherefficiency of ammonia recovery. The gypsum solution recovered around 95% of the stripped ammonia when thestripping gas contained 10% CO 2 and it rose to 100% when the gas contained more than 30% CO 2 . Furthermore,a statistical model and genetic algorithm were applied to estimate the total ammonia recovery and optimize theoperating parameters. The experimental and modeling results indicated that integrating biogas stripping withgypsum absorption can effectively harvest ammonia nitrogen and reduce CO 2 emissions while producing(NH 4 ) 2 SO 4 fertilizer from AD effluent.

59 BASIC BIOLOGICAL SCIENCES↗

Dataset from: Improved multispecies Dougherty collisions

The Dougherty model Fokker-Planck operator is extended to describe nonlinear full-f collisions between multiple species in plasmas. Simple relations for cross-species primitive moments are developed which obey conservation laws, and reproduce familiar velocity and temperature relaxation rates. This treatment of multispecies Dougherty collisions, valid for arbitrary mass ratios, avoids unphysical temperatures and satisfies the H-theorem unlike an analogous Bhatnagar-Gross-Krook operator. Formulas for both a Cartesian velocity-space and a gyroaveraged operator are provided for use in Vlasov as well as long-wavelength gyrokinetic models. We present an algorithm for the discontinuous Galerkin discretization of this operator, and provide results from relaxation and Landau damping benchmarks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Rapid Quantum Ground State Preparation via Dissipative Dynamics

Inspired by natural cooling processes, dissipation has become a promising approach for preparing low-energy states of quantum systems. However, the potential of dissipative protocols remains unclear beyond certain commuting Hamiltonians. This work provides significant analytical and numerical insights into the power of dissipation for preparing the ground state of noncommuting Hamiltonians. For quasi-free dissipative dynamics, including certain 1D spin systems with boundary dissipation, our results reveal a new connection between the mixing time in trace distance and the spectral properties of a non-Hermitian Hamiltonian, leading to an explicit and sharp bound on the mixing time that scales polynomially with system size. For more general spin systems, we develop a tensor network-based algorithm for constructing the Lindblad jump operator and for simulating the dynamics. Using this algorithm, we demonstrate numerically that dissipative ground state preparation protocols can achieve rapid mixing for certain 1D local Hamiltonians under bulk dissipation, with a mixing time that scales logarithmically with the system size. We then prove the rapid mixing result for certain weakly interacting spin and fermionic systems in arbitrary dimensions, extending recent results for high-temperature quantum Gibbs samplers to the zero-temperature regime. Together, these results show that dissipation can be a powerful tool for ground state preparation, with potential applications across condensed matter physics, quantum materials science, and beyond.

decoherence↗

A novel triple power cycle featuring a gas turbine cycle with supercritical carbon dioxide and organic Rankine cycles: Thermoeconomic analysis and optimization

In this study, a novel triple power cycle is proposed where waste heat from a gas turbine cycle is utilized to drive a supercritical carbon dioxide (s-CO 2 ) recompression cycle and a recuperative organic Rankine cycle (ORC) in sequence. A detailed thermoeconomic model is developed and implemented in MATLAB to evaluate the performance of the proposed cycle under different operating conditions. Optimization using a particle swarm optimization (PSO) algorithm is performed to minimize the levelized cost of electricity (LCOE) and determine the optimum design conditions of the cycle. The optimization results show that for a 100 MW cycle, the overall thermal efficiency and LCOE are 0.521 and $\$52.819$/MWh, respectively. The turbine inlet temperature of the gas turbine and s-CO 2 cycles are found as the most influential parameters on the thermoeconomic performance of the triple cycle. The proposed triple cycle shows an excellent waste energy recovery potential and superiority over a very thermodynamically efficient cycle from the literature, which included a gas turbine topping cycle and a complex cascade s-CO 2 power cycle used a bottoming cycle. As a result, the proposed triple cycle shows up to 0.9% points higher efficiency while it has fewer heat exchangers and turbomachinery than the cycle from the literature.

30 DIRECT ENERGY CONVERSION↗

Resolve the Clear-Sky Continuous Diurnal Cycle of High-Resolution ECOSTRESS Evapotranspiration and Land Surface Temperature

Although numerous evapotranspiration (ET) estimates have been developed from remote sensing measurements, a diurnally resolved product with high spatial resolution (<100 m) is still lacking, which is critically needed for agricultural and ecological monitoring. The ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) provides, for the first time, ET at 70 m resolution with diurnal samplings, enabling a step toward filling this data gap. This study developed the first framework to resolve the full diurnal cycle of ET and Land Surface Temperature (LST) from sporadic ECOSTRESS measurements. We first constructed 70 m diurnal LST utilizing a diurnal temperature cycle (DTC) model that fused ECOSTRESS and Geostationary Operational Environmental Satellite (GOES) LST. Next, we derived 70 m diurnal ET from the diurnally resolved LST, along with ancillary meteorological and reflectance data sets, using the Priestley-Taylor Jet Propulsion Laboratory (PT-JPL) algorithm. This framework was tested in two regions: A relatively homogeneous semiarid grassland landscape (centered at the US-Seg Ameriflux site) and a heterogeneous humid cropland landscape (centered at US-ARM). Our diurnally resolved LST and ET successfully reproduced the spatial patterns in the native ECOSTRESS measurements during overpasses for both landscapes (LST: r > 0.97, RMSE ~2°C; ET: r > 0.99, and RMSE = 1–7 W m -2 ). The constructed time series also well captured the temporal variability of the in situ measurements (especially under clear-sky conditions) at both US-Seg and US-ARM, albeit performance varying with availability and quality of ECOSTRESS samples. This study sets the stage for testing and applying our framework to broader climates, biomes, and landscapes toward eventually generating diurnally resolved 70 m global operational LST/ET products to enhance ecological/agricultural applications.

54 ENVIRONMENTAL SCIENCES↗

Tree-based solvers for adaptive mesh refinement code $\scriptsize{FLASH}$ – IV. An X-ray radiation scheme to couple discrete and diffuse X-ray emission sources to the thermochemistry of the interstellar medium

X-ray radiation, in particular radiation between 0.1 and 10 keV, is evident from both point-like sources, such as compact objects and T-Tauri young stellar objects, and extended emission from hot, cooling gas, such as in supernova remnants. The X-ray radiation is absorbed by nearby gas, providing a source of both heating and ionization. While protoplanetary chemistry models now often include X-ray emission from the central young stellar object, simulations of star-forming regions have yet to include X-ray emission coupled to the chemo-dynamical evolution of the gas. We present an extension of the $\scriptsize{TREERAY}$ reverse ray trace algorithm implemented in the flash magnetohydrodynamic code which enables the inclusion of X-ray radiation from 0.1 keV < E γ < 100 keV, dubbed $\scriptsize{XRAYTHESPOT}$. $\scriptsize{XRAYTHESPOT}$ allows for the use of an arbitrary number of bins, minimum and maximum energies, and both temperature-independent and temperature-dependent user-defined cross-sections, along with the ability to include both point and extended diffuse emission and is coupled to the thermochemical evolution. We demonstrate the method with several multibin benchmarks testing the radiation transfer solution and coupling to the thermochemistry. Finally, we show two example star formation science cases for this module: X-ray emission from protostellar accretion irradiating an accretion disc and simulations of molecular clouds with active chemistry, radiation pressure, and protostellar radiation feedback from infrared to X-ray radiation.

79 ASTRONOMY AND ASTROPHYSICS↗

A clustering-based approach to ocean model–data comparison around Antarctica

The Antarctic Continental Shelf seas (ACSS) are a critical, rapidly changing element of the Earth system. Analyses of global-scale general circulation model (GCM) simulations, including those available through the Coupled Model Intercomparison Project, Phase 6 (CMIP6), can help reveal the origins of observed changes and predict the future evolution of the ACSS. However, an evaluation of ACSS hydrography in GCMs is vital: previous CMIP ensembles exhibit substantial mean-state biases (reflecting, for example, misplaced water masses) with a wide inter-model spread. Because the ACSS are also a sparely sampled region, grid-point-based model assessments are of limited value. Our goal is to demonstrate the utility of clustering tools for identifying hydrographic regimes that are common to different source fields (model or data), while allowing for biases in other metrics (e.g., water mass core properties) and shifts in region boundaries. We apply K-means clustering to hydrographic metrics based on the stratification from one GCM (Community Earth System Model version 2; CESM2) and one observation-based product (World Ocean Atlas 2018; WOA), focusing on the Amundsen, Bellingshausen and Ross seas. When applied to WOA temperature and salinity profiles, clustering identifies “primary” and “mixed” regimes that have physically interpretable bases. For example, meltwater-freshened coastal currents in the Amundsen Sea and a region of high-salinity shelf water formation in the southwestern Ross Sea emerge naturally from the algorithm. Both regions also exhibit clearly differentiated inner- and outer-shelf regimes. The same analysis applied to CESM2 demonstrates that, although mean-state model biases in water mass T–S characteristics can be substantial, using a clustering approach highlights that the relative differences between regimes and the locations where each regime dominates are well represented in the model. CESM2 is generally fresher and warmer than WOA and has a limited fresh-water-enriched coastal regimes. Given the sparsity of observations of the ACSS, this technique is a promising tool for the evaluation of a larger model ensemble (e.g., CMIP6) on a circum-Antarctic basis.

54 ENVIRONMENTAL SCIENCES↗

Simulated meteorological impacts of offshore wind turbines and sensitivity to the amount of added turbulence kinetic energy

Offshore wind energy projects are currently in development off the east coast of the United States and may influence the local meteorology of the region. Wind power production and other commercial uses in this area are related to atmospheric conditions, and so it is important to understand how future wind plants may change the local meteorology. In the absence of measurements of potential wind plant impacts on meteorology, simulations offer the next-best possible insight into wake effects on boundary layer height, temperature, fluxes, and wind speeds. However, simulation tools that capture these effects offer multiple options for representing the amount of turbine-added turbulence that may impact assessments of micrometeorological effects. To explore this sensitivity, we compare 1 year of simulations from the Weather Research and Forecasting (WRF) model with and without wind plants incorporated, focusing on the lease area south of Massachusetts and Rhode Island. The simulations with wind plants are repeated to include both the maximum and minimum amounts of added turbulence to provide bounds on the potential impacts. We assess changes in wind speeds, 2 m temperature, surface heat flux, turbulence kinetic energy (TKE), and boundary layer height during different stability classifications and ambient wind speeds over the entire year and compare results for the degree of added turbulence in the wind plant simulations. Because the wake behavior may be a function of boundary layer stability, in this paper, we also present a machine learning algorithm to quantify the area and distance of the wake generated by the wind plant. This analysis enables us to identify the relationship between wake extent and boundary layer height. We find that hub-height wind speed is reduced within and downwind of the wind plant, with the strongest impacts occurring during stable conditions and faster wind speeds in region 3 of the turbine power curve, although impacts lessen as wind speeds increase past 15 m s−1. In contrast, wind speeds near the surface decrease when no turbine-added turbulence is included but can increase for stably stratified conditions when 100 % of possible TKE is included in the simulations. TKE increases at hub height in the simulations with added TKE for all stability classes, suggesting that atmospheric stability does not immediately modify the TKE generated by turbines. Negligible changes in hub-height TKE manifest in the simulations without the added TKE. At the surface, TKE increases in the simulations with maximum added turbulence only for unstable conditions. In the no-added-turbulence simulations, surface TKE decreases slightly in neutral and unstable simulations. Differences in 2 m temperatures and surface heat fluxes are small but vary considerably with atmospheric stability and the amount of added TKE. Boundary layer heights increase within the wind plant when turbine-added turbulence is included and decrease slightly downwind during stable conditions. In contrast, with no added turbulence, the boundary layer height is in general reduced in stable conditions with wind speeds less than 15 m s −1 and slightly increased in neutral conditions. Finally, shallower upwind boundary layer heights tend to correlate with larger wake areas and distances, though other factors likely also play a role in determining the extent of the wind plant wake. These simulation-based results provide a bound for micrometeorological impacts of wind plant wakes: simulations that couple the atmosphere to the ocean may reduce these impacts, and we await observational verification.

17 WIND ENERGY↗

High-order fully implicit solver for all-speed fluid dynamics: AUSM ride from nearly incompressible variable-density flows to shock dynamics

We introduce a novel Newton–Krylov (NK)-based fully implicit algorithm for solving fluid flows in a wide range of flow conditions—from variable density nearly incompressible to supersonic shock dynamics. The key enabling feature of our all-speed solver is the ability to efficiently solve conservation laws by choosing a set of independent variables that produce a well-conditioned Jacobian matrix for the linear iterations of the global nonlinear iterative solver. In particular, instead of choosing to discretize the conservative variables (density, momentum, total energy), which is traditionally used in Eulerian high-speed compressible fluid dynamics, we demonstrate superior performance by discretizing the primitive variables—pressure–velocity–temperature in the very low-Mach flow limits or density–velocity–temperature/entropy in the shock dynamics range. Moreover, our method allows us to avoid direct inversion of the mass matrix in discrete time derivatives, which is usually an additional source for stiffness, especially pronounced when going to very high-order schemes with non-orthogonal basis functions. Here, we show robust solutions obtained for discontinuous finite element discretization up to seventh-order accuracy. Another important aspect of the solution algorithm is the Advection Upstream Splitting Method (AUSM), adopted to compute numerical fluxes within our reconstructed discontinuous Galerkin (rDG) spatial discretization scheme. The use of the low-Mach modification of the hyperbolic flux operator is found to be necessary for enabling robust simulations of very stiff liquids and metals for Mach numbers below M=10 -5 , which is well known to be very computationally challenging for compressible solvers. We demonstrate that our fully implicit rDG-NK solver with the AUSM + -up flux treatment produces efficient and high-resolution numerical solutions at all speeds, ranging from vanishing Mach numbers to transonic and supersonic, without substantial modifications of the solution procedures. (At high speed, we add limiting and use a simpler preconditioning of the Krylov solver.) Numerical examples include nearly incompressible constant-property flow past a backward-facing step with heat transfer, low-Mach variable-property channel flow of water at supercritical state, phase change and melt pool dynamics for laser spot welding and selective laser melting in additive manufacturing, and Mach 3 flow in a wind tunnel with a step.

97 MATHEMATICS AND COMPUTING↗