Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient based optimization”

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.

970 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Design and Analysis of an X-Ray Mirror Assembly Using the Meta-Shell Approach

Lightweight and high resolution optics are needed for future space-based x-ray telescopes to achieve advances in high-energy astrophysics. Past missions such as Chandra and XMM-Newton have achieved excellent angular resolution using a full shell mirror approach. Other missions such as Suzaku and NuSTAR have achieved lightweight mirrors using a segmented approach. This paper describes a new approach, called meta-shells, which combines the fabrication advantages of segmented optics with the alignment advantages of full shell optics. Meta-shells are built by layering overlapping mirror segments onto a central structural shell. The resulting optic has the stiffness and rotational symmetry of a full shell, but with an order of magnitude greater collecting area. Several meta-shells so constructed can be integrated into a large x-ray mirror assembly by proven methods used for Chandra and XMM-Newton. The mirror segments are mounted to the meta-shell using a novel four point semi-kinematic mount. The four point mount deterministically locates the segment in its most performance sensitive degrees of freedom. Extensive analysis has been performed to demonstrate the feasibility of the four point mount and meta-shell approach. A mathematical model of a meta-shell constructed with mirror segments bonded at four points and subject to launch loads has been developed to determine the optimal design parameters, namely bond size, mirror segment span, and number of layers per meta-shell. The parameters of an example 1.3 m diameter mirror assembly are given including the predicted effective area. To verify the mathematical model and support opto-mechanical analysis, a detailed finite element model of a meta-shell was created. Finite element analysis predicts low gravity distortion and low sensitivity to thermal gradients.

module↗
NASA NTRS1996

Adaptive grid methods for RLV environment assessment and nozzle analysis

Rapid access to highly accurate data about complex configurations is needed for multi-disciplinary optimization and design. In order to efficiently meet these requirements a closer coupling between the analysis algorithms and the discretization process is needed. In some cases, such as free surface, temporally varying geometries, and fluid structure interaction, the need is unavoidable. In other cases the need is to rapidly generate and modify high quality grids. Techniques such as unstructured and/or solution-adaptive methods can be used to speed the grid generation process and to automatically cluster mesh points in regions of interest. Global features of the flow can be significantly affected by isolated regions of inadequately resolved flow. These regions may not exhibit high gradients and can be difficult to detect. Thus excessive resolution in certain regions does not necessarily increase the accuracy of the overall solution. Several approaches have been employed for both structured and unstructured grid adaption. The most widely used involve grid point redistribution, local grid point enrichment/derefinement or local modification of the actual flow solver. However, the success of any one of these methods ultimately depends on the feature detection algorithm used to determine solution domain regions which require a fine mesh for their accurate representation. Typically, weight functions are constructed to mimic the local truncation error and may require substantial user input. Most problems of engineering interest involve multi-block grids and widely disparate length scales. Hence, it is desirable that the adaptive grid feature detection algorithm be developed to recognize flow structures of different type as well as differing intensity, and adequately address scaling and normalization across blocks. These weight functions can then be used to construct blending functions for algebraic redistribution, interpolation functions for unstructured grid generation, forcing functions to attract/repel points in an elliptic system, or to trigger local refinement, based upon application of an equidistribution principle. The popularity of solution-adaptive techniques is growing in tandem with unstructured methods. The difficultly of precisely controlling mesh densities and orientations with current unstructured grid generation systems has driven the use of solution-adaptive meshing. Use of derivatives of density or pressure are widely used for construction of such weight functions, and have been proven very successful for inviscid flows with shocks. However, less success has been realized for flowfields with viscous layers, vortices or shocks of disparate strength. It is difficult to maintain the appropriate mesh point spacing in the various regions which require a fine spacing for adequate resolution. Mesh points often migrate from important regions due to refinement of dominant features. An example of this is the well know tendency of adaptive methods to increase the resolution of shocks in the flowfield around airfoils, but in the incorrect location due to inadequate resolution of the stagnation region. This problem has been the motivation for this research.

Thornburg, Hugh J.↗
DOE OSTI2023

Quantitative proton radiography and shadowgraphy for arbitrary intensities

Charged-particle radiography and shadowgraphy data can be directly inverted to obtain a line-integrated transverse Lorentz force or a line-integrated transverse refractive index gradient if intensity modulations due to scattering and absorption are negligible, and angular deflections are small. We develop a new direct-inversion algorithm based on plasma physics and compare it to a new Monge–Ampère code and an existing power diagram code. The measured or source intensity is represented by electrons subject to drag, and the other intensity by fixed ions. The decrease in kinetic plus electrostatic energy determines convergence. The displacement of the electrons from their initial to their equilibrium positions determines the line-integrated force or refractive index gradient. We have implemented two approaches: PIC (particle in cell) and Lagrangian fluid, in 1-D and 2-D. The PIC code works for arbitrary intensities, can work efficiently in parallel, and can make use of existing codes. The Lagrangian code requires less memory and is faster than the PIC code without massively parallel processing, but fails in 2-D for large intensity modulations. The Monge–Ampère code is by far the fastest in 2-D, without massively parallel processing, but fails for intensities with large voids, high contrast ratios and large deflections across the boundaries, and could not obtain the degree of convergence possible with the PIC code. As a result, the power diagram code was by far the slowest and most memory intensive, and failed for large peaks in the measured intensity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Correlating Surface Processing of Nb Superconducting RF Cavities with the Evolution of Surface Electronic States

Superconducting-radio frequency (RF) cavities provide an efficient way to accelerate particle beams with extremely high acceleration gradients while generating very small power dissipation. The few nanometers of the surface play a critical role in defining the RF performance of superconducting Nb based cavities. Over the past two decades, several pioneering surface treatment and processing methods have emerged, enabling remarkable improvements in superconducting cavity performance by simultaneously achieving high quality factors with increasing maximum acceleration gradients. These processing approaches include chemical polishing, distinct multi-step thermal treatments under ultra-high vacuum (UHV) conditions over low to high temperature regimes, as well as high-temperature treatments under controlled nitrogen atmospheres. Beyond their macroscopic impact on RF performance, these methods produce distinct surface oxide configurations characterized by different valence states, oxide thicknesses, chemical uniformity, and oxygen concentration profiles extending into the near-surface bulk of niobium. In this work, we are trying to understand how the surface-processing methods and the resulting oxide/oxygen profiles affect the electronic structure of surface and the mechanism of superconductivity. Using a combination of X-ray photoemission and X-ray absorption spectroscopies, we investigate how the valence-band structure and the electronic density of states (DoS) near the Fermi level evolve under different surface treatments. By employing tunable photon energies across multiple elemental absorption edges, we use resonant photoemission to disentangle and identify the elemental contributions to specific valence-band features. Our observations show that different surface-processing methods lead to distinct temperature evolutions of the DoS and valence-band states near the Fermi level. Our results suggest variations in different Nb-O orbital hybridizations in distinct processes and point towards the possibility of different underlying mechanisms of superconductivity governed by surface chemistry and oxide configuration. We also correlate these distinct superconducting mechanisms with RF cavity performance, specifically focusing on measured surface resistance, the nature of the Q-slope, and quench fields observed in superconducting RF measurements. These results will enable us to identify the potential limiting factors and relevant controllable parameters that can be further optimized to improve the performance of superconducting RF cavities.

Tripathi, Malvika [Fermilab] (ORCID:00000001989251↗
NASA NTRS2007

A Microfabricated Involute-Foil Regenerator for Stirling Engines

A segmented involute-foil regenerator has been designed, microfabricated and tested in an oscillating-flow rig with excellent results. During the Phase I effort, several approximations of parallel-plate regenerator geometry were chosen as potential candidates for a new microfabrication concept. Potential manufacturers and processes were surveyed. The selected concept consisted of stacked segmented-involute-foil disks (or annular portions of disks), originally to be microfabricated from stainless-steel via the LiGA (lithography, electroplating, and molding) process and EDM (electric discharge machining). During Phase II, re-planning of the effort led to test plans based on nickel disks, microfabricated via the LiGA process, only. A stack of nickel segmented-involute-foil disks was tested in an oscillating-flow test rig. These test results yielded a performance figure of merit (roughly the ratio of heat transfer to pressure drop) of about twice that of the 90% random fiber currently used in small ~ 100 W Stirling space-power convertors in the Reynolds Number range of interest (50-100). A Phase III effort is now underway to fabricate and test a segmented-involute-foil regenerator in a Stirling convertor. Though funding limitations prevent optimization of the Stirling engine geometry for use with this regenerator, the Sage computer code will be used to help evaluate the engine test results. Previous Sage Stirling model projections have indicated that a segmented-involute-foil regenerator is capable of improving the performance of an optimized involute-foil engine by 6-9%; it is also anticipated that such involute-foil geometries will be more reliable and easier to manufacture with tight-tolerance characteristics, than random-fiber or wire-screen regenerators. Beyond the near-term Phase III regenerator fabrication and engine testing, other goals are (1) fabrication from a material suitable for high temperature Stirling operation (up to 850 C for current engines; up to 1200 C for a potential engine-cooler for a Venus mission), and (2) reduction of the cost of the fabrication process to make it more suitable for terrestrial applications of segmented involute foils. Past attempts have been made to use wrapped foils to approximate the large theoretical figures of merit projected for parallel plates. Such metal wrapped foils have never proved very successful, apparently due to the difficulties of fabricating wrapped-foils with uniform gaps and maintaining the gaps under the stress of time-varying temperature gradients during start-up and shut-down, and relatively-steady temperature gradients during normal operation. In contrast, stacks of involute-foil disks, with each disk consisting of multiple involute-foil segments held between concentric circular ribs, have relatively robust structures. The oscillating-flow rig tests of the segmented-involute-foil regenerator have demonstrated a shift in regenerator performance strongly in the direction of the theoretical performance of ideal parallel-plate regenerators.

Tew, Roy↗
DOE OSTI2021

Development of High-efficiency and Cost-effective Forged Ingot Niobium Technology for Science Frontiers and Accelerator Applications

Development of Forged Ingot Niobium Technology: Worlds science frontier programs and SRF accelerator applications demand high performance and cost-effective SRF accelerator technology [1-8]. Fine-grain (FG) and Large-grain (LG) ingot niobium technologies have been very well developed and implemented in all the present-day accelerator projects. However, forged ingot niobium technology which is the focus of this development proposal will be much more cost-effective and expected to have several technical advantages. FG niobium sheet production is very complex involving more than ten processing steps making them prone to contamination. As a result, they are very expensive to produce and require stringent QA procedure to be ready for SRF cavity production. The accelerating cavity process steps are also numerous and require strict procedures in order to achieve high accelerating gradients and quality factors needed for science frontier programs. LG niobium disc production, directly sliced from the ingot, is relatively simple and straight forward to keep surface cleanness. The disc production cost is significantly low compared to FG niobium sheet production. However, there are (some) draw backs due to non-homogeneity of the grain boundary distribution, resulting in non-uniform mechanical properties and complex cavity fabrication, although the LG cavities achieve the expected high-gradient performance goals with lower cost. Medium-grain (MG) niobium disk production may be realized with a new approach/process, the disc directly sliced from the forged ingot, involves a simpler process steps contributing major production cost reduction [9]. These discs are expected to be superior as they tend to be homogenous with uniform sub millimeter grains and mechanical properties. We are eagerly looking forward to developing the forged ingot niobium SRF accelerator technology for the benefit of the world-wide science frontier programs, green energy subcritical nuclear energy systems and a wide variety of industrial applications including the production of radio isotopes and nuclear transmutation applications. Measurement of thermal characteristics of the forged ingot niobium: Measurement of the thermal diffusivity, D, of superconducting MG niobium is important to be understood in comparison with FG sheet and LG disc, as well as that of advanced composite material of Nb3Sn film sputtered on forged ingot Nb, which will be determined using transient pump-probe thermo-modulation [10,11]. Transient thermo-modulation is based on using an ultrafast laser pulse to heat the superconducting materials by a few K, then a synchronized laser pulse probes the reflectance of the heated material. For thin films on a substrate (e.g, 100-nm Nb3Sn on Nb), it takes <100 ps for the heat to reach the substrate by diffusion. Therefore, an ultrafast method is needed to probe D of the studied material. D will be measured from room temperature to liquid helium temperatures. Fundamental Research: The increase of rf loss in SRF cavities is related to trapping of residual magnetic field during the cavity cool-down. The study suggested that the micro structure plays the role in flux trapping sensitivity. The research includes the flux trapping and expulsion study as forged ingot niobium goes through several mechanical deformations, crystallization leading to the optimal SRF cavity performance. Furthermore, the cavity made from FG, MG, and LG sheet/disc will be compared to understand the optimal re-crystallization temperature that cavity needed to be heat treated which minimize the flux trapping and increase SRF cavity performance. References: 1. S. Belomestnykh, Overview of recent SRF developments for ERLs, presented at the SRF 2015, Whistler, BC, Canada (2015). 2. A. Yamamoto, M. Yamanaka and G. Myneni, Ingot Nb based SRF Technology for the International Linear Collider, in Science and Technology of Ingot Niobium for Superconducting Radio Frequency Applications, AIP Conf. Proc. 1687, 030005-1 ? 03005-6, 2015 3. Report of the Workshop on Energy and Environmental Applications of Accelerators, DOE Workshop Report. https://science.osti.gov/-/media/hep/pdf/Reports/2020/CASM_WorkshopReport.pdf? la=en&hash=AEB0B318ED0436B1C5FF4EE0FDD6DEB84C2F15B2 4. G. Ciovati, et al., Design of a cw, low-energy, high-power superconducting linac for environmental applications, Phys. Rev. Accel. Beams 21, 091601 (2018). 5. P. Dhakal, et al., Effect of high temperature heat treatments on the quality factor of a large-grain superconducting radio-frequency niobium cavity, Phys. Rev. ST Accel. Beams 16, 042001 (2013). 6. P. Kneisel et al., Review of ingot niobium as a material for superconducting radio frequency accelerating cavities, Nuclear Instruments and Methods in Physics Research A 774, 133 (2015). 7. G. Ciovati, P. Dhakal, and G. R, Myneni, Superconducting radio-frequency cavities made from medium and low-purity niobium ingots, Supercond. Sci. Technol. 29, 064002 (2016). 8. M. Drury et al., commissioning of the prototype C75 cavities in a CEBAF cryomodule, in Proc. of IPAC 2018, Vancouver, BC, Canada (2018) 9. Feasibility of forged-ingot niobium disc and SRF cavity fabrication technology recently demonstrated in cooperation of ATI, BSCE, and KEK (2020). 10. Hani E. Elsayed-Ali, ?Measurements of heat transport in thin films by ultrafast laser-based techniques,? 3rd International Conference on Thermal Issues in Emerging Technologies Theory and Applications, Cairo, Egypt, pp. 347 ? 350 (2010). DOI: 10.1109/ThETA17616.2010 11. W. M. G. Ibrahim, H. E. Elsayed-Ali, M. Schinn, and C. A. Bonner, Jr., ?Ultrafast investigation of electron dynamics in multilayer metals,? Int. J. Heat and Mass Transfer, 47(10?11), 2261?2268 (2004).

Myneni, G.↗
NASA NTRS2007

A Microfabricated Involute-Foil Regenerator for Stirling Engines

A segmented involute-foil regenerator has been designed, microfabricated and tested in an oscillating-flow rig with excellent results. During the Phase I effort, several approximations of parallel-plate regenerator geometry were chosen as potential candidates for a new microfabrication concept. Potential manufacturers and processes were surveyed. The selected concept consisted of stacked segmented-involute-foil disks (or annular portions of disks), originally to be microfabricated from stainless-steel via the LiGA (lithography, electroplating, and molding) process and EDM. During Phase II, re-planning of the effort led to test plans based on nickel disks, microfabricated via the LiGA process, only. A stack of nickel segmented-involute-foil disks was tested in an oscillating-flow test rig. These test results yielded a performance figure of merit (roughly the ratio of heat transfer to pressure drop) of about twice that of the 90 percent random fiber currently used in small approx.100 W Stirling space-power convertors-in the Reynolds Number range of interest (50 to 100). A Phase III effort is now underway to fabricate and test a segmented-involute-foil regenerator in a Stirling convertor. Though funding limitations prevent optimization of the Stirling engine geometry for use with this regenerator, the Sage computer code will be used to help evaluate the engine test results. Previous Sage Stirling model projections have indicated that a segmented-involute-foil regenerator is capable of improving the performance of an optimized involute-foil engine by 6 to 9 percent; it is also anticipated that such involute-foil geometries will be more reliable and easier to manufacture with tight-tolerance characteristics, than random-fiber or wire-screen regenerators. Beyond the near-term Phase III regenerator fabrication and engine testing, other goals are (1) fabrication from a material suitable for high temperature Stirling operation (up to 850 C for current engines; up to 1200 C for a potential engine-cooler for a Venus mission), and (2) reduction of the cost of the fabrication process to make it more suitable for terrestrial applications of segmented involute foils. Past attempts have been made to use wrapped foils to approximate the large theoretical figures of merit projected for parallel plates. Such metal wrapped foils have never proved very successful, apparently due to the difficulties of fabricating wrapped-foils with uniform gaps and maintaining the gaps under the stress of time-varying temperature gradients during start-up and shut-down, and relatively-steady temperature gradients during normal operation. In contrast, stacks of involute-foil disks, with each disk consisting of multiple involute-foil segments held between concentric circular ribs, have relatively robust structures. The oscillating-flow rig tests of the segmented-involute-foil regenerator have demonstrated a shift in regenerator performance strongly in the direction of the theoretical performance of ideal parallel-plate regenerators.

Tew, Roy↗
DOE OSTI2023

Evaluation of Thermal Neutron Scattering Cross Section of Uranium Silicide with Ab Initio Lattice Dynamics

Uranium silicide (U 3 Si 2 ) is a candidate material for the high-density nuclear fuel in commercial light water reactors [1], [2]. Its higher uranium density, 11.3 g-U/cm3, compared to that of uranium dioxide (UO 2 ), 9.7 g-U/cm3, can improve the performance of a nuclear reactor while using low enriched uranium (LEU) and diversify the choice of cladding materials [1]–[3]. It also has a higher thermal conductivity than UO 2 , which can reduce the thermal stress on the material caused by a temperature gradient across the fuel pellet and provide a larger margin for some postulated accidents [1], [2], [4]. Furthermore, compared to U3Si, another high-density fuel candidate, it has better resistance to in-pile swelling due to less irradiation-induced rapid amorphization [1], [3]. Corresponding to its importance in nuclear engineering, many previous studies have reported the properties of U3Si2. Experiments showed that U 3 Si 2 is a paramagnetic (PM) metal, where a slight linear increase in magnetic susceptibility was measured with increasing temperature [5], [6]. In addition, thermodynamic quantities such as thermal expansion coefficient, heat capacity, and thermal conductivity were experimentally determined over a wide temperature range [1], [7], [8]. In several computational studies, ab initio atomistic simulations based on density functional theory (DFT) were performed to calculate various properties including elastic constants, electronic density of states (DOS), and phonon dispersion curves [9]–[12]. Nevertheless, thermal neutron scattering cross sections, which are critical to the prediction of the parameters in reactor physics that are ultimately related to reactor criticality, have not yet been evaluated for U3Si2. The scattering cross section can be calculated from the phonon DOS, or the energy spectrum of lattice vibrations, of the crystalline system [13], [14]. However, there is also no experimental data available for the phonon DOS of U 3 Si 2 . While some computational studies reported the phonon DOS and/or dispersion curves from ab initio simulations [9]–[12], the accuracy cannot be guaranteed because it is unclear whether the spin-polarization behavior of PM U 3 Si 2 was properly described. In the present study, the thermal neutron scattering cross section for U 3 Si 2 is evaluated for the first time by calculating the phonon DOS for U3Si2 from ab initio lattice dynamics (AILD) simulations based on DFT. First, U 3 Si 2 is modeled based on the experimental structure, and AILD simulations are performed on the modeled U3Si2 to optimize the structure. Next, AILD simulations are performed for supercells with atomic displacement to calculate Hellmann-Feynman forces. Based on the calculated forces, partial phonon DOSs for U and Si are obtained, and the thermal neutron scattering law (TSL) for U 3 Si 2 is finally evaluated. To verify the accuracy of the calculations in the present study, the calculation results are compared with experimental data on the structure and heat capacity of U3Si2 [1], [7], [8], [15].

Geometry Optimization↗
NASA NTRS1996

Modeling Orifice Pulse Tube Coolers

We have developed a calculational model that treats all the components of an orifice pulse tube cooler. We base our analysis on 1-dimensional thermodynamic equations for the regenerator and we assume that all mass flows, pressure oscillations and temperature oscillations are small and sinusoidal. Non-linear pressure drop effects are included in the regenerator to account for finite pressure amplitude effects. The resulting mass flows and pressures are matched at the boundaries with the other components of the cooler: compressor, aftercooler, cold heat exchanger, pulse tube, hot heat exchanger, orifice and reservoir. The results of the calculation are oscillating pressures, mass flows and enthalpy flows in the main components of the cooler. By comparing with the calculations of other available models, we show that our model is very similar to REGEN 3 from NIST and DeltaE from Los Alamos National Lab. Our model is much easier to use than other available models because of its simple graphical interface and the fact that no guesses are required for the operating pressures or mass flows. In addition, the model only requires a few minutes of running time allowing many parameters to be optimized in a reasonable time. A version of the model is available for use over the World Wide Web at http://irtek.arc.nasa.gov. Future enhancements include adding a bypass orifice and including second order terms in steady mass streaming and steady heat transfer. A two-dimensional anelastic approximation of the fluid equations will be used as the basis for the latter analysis. Preliminary results are given in dimensionless numbers appropriate for oscillating compressible flows. The model shows how transverse heat transfer reduces enthalpy flow, particularly for small pulse tubes. The model also clearly shows mass recirculation in the open tube on the order of the tube length. They result from the higher order Reynolds stresses. An interesting result of the linearized approach is that the steady mass streaming does not affect the enthalpy flow at second order. The major effect of recirculating mass streaming is to increase transverse temperature gradients, which leads to higher entropy production and reduced efficiency.

Roach, Kittel P.↗
DOE OSTI2024

Evaluation and development of flow condensation correlations using the data from low GWP refrigerants in an axial micro-fin aluminum tube

To mitigate global warming, the world is transitioning to refrigerants with low global warming potential (GWP). Supporting this shift requires a model that can accurately predict the heat transfer and pressure drop of new refrigerants, crucial for designing efficient heat exchangers. Existing models, however, are largely based on currently deployed refrigerants and primarily developed for unexpanded micro-fin tubes with spiral angles of 6° to 30°. Their applicability to new refrigerants, especially in expanded micro-fin tubes, is uncertain. This study assesses the performance of four well-known condensation models for six emerging refrigerants—R-32, R-454B, R-454C, R-455A, R-1234yf, and R-1234ze(E)—against experimental data. Initially, the Han and Lee (2005) model shows the best prediction accuracy with a mean absolute deviation (MAD) of 22.1 %. To enhance the accuracy of heat transfer models for new refrigerants and geometries with large temperature glides, two approaches are proposed. Here, the first approach applies a simple correction factor, reducing the MAD of the Cavallini et al. (2009) model from 68.2 % to 15.4 %. The second approach uses the variable metric method for minimization, fitting new constants to the data. This optimization results in the Kedzierski and Goncalves (1997) model achieving the highest accuracy, with a MAD of 13.1 %. For pressure drop models, the Cavallini et al. (1997) model is the most accurate with a MAD of 6.4 %, followed by the Haraguchi et al. (1993) model with a MAD of 9.4 %. Due to its simplicity, the Haraguchi et al. (1993) model is a practical option for predicting frictional pressure drop.

Axial micro-fin tubes↗
DOE OSTI2023

An adaptive synchronous extraction (ASE) method for estimating intensity and footprint of surface urban heat islands: A case study of 254 North American cities

The urban heat island (UHI) effect has attracted great attention due to its potential impacts on rapidly growing urban areas. Using remotely sensed estimates of land surface temperature (LST), a large number of studies have focused on the surface UHI (SUHI) effect, which can be characterized by its two fundamental properties: intensity and footprint. The SUHI intensity reflects the LST difference between the urban area and the background reference area (BRA), and the SUHI footprint indicates the spatial extent influenced by the heat island. Currently, numerous methods have been developed to estimate the SUHI intensity and footprint, but are still greatly challenged by three main issues. Namely, the discrepancy in BRA selection criterion brings great uncertainty to the estimated SUHI intensity, the estimation of SUHI footprint is largely constrained by the predefined models, and the quantification of SUHI effect is potentially influenced by several confounding factors. Here, we proposed an adaptive synchronous extraction (ASE) method, which is capable of adaptively selecting the most optimal BRA while removing the influence of confounding factors, and achieving synchronous estimation of SUHI intensity and footprint. We applied the ASE method to 254 North American cities and conducted an in-depth comparative analysis to discuss its applicability and benefits. The main results include: (1) The ASE method avoids the limitations of existing methods in BRA selection and model presetting, and shows resilience to parameter variations. This makes the ASE method highly applicable to quantify the SUHI intensity and footprint in cities with various thermal characteristics. (2) The ASE method can better highlight the spatial, seasonal and day-night contrasts in the estimated SUHI intensity. This superiority is particularly evident when comparing it to methods based on the equal-area buffer or the simplified urban-extent algorithm. (3) Confounding factors pose non-negligible impacts on the quantification of the SUHI effect. Typically, ignoring the influence of topographic relief or missing LST data can lead to an overall overestimation of the SUHI intensity, while not removing surrounding urban areas will cause some underestimation of the SUHI intensity. In conclusion, overall, the proposed ASE method provides a new generalizable tool for quantifying the SUHI effect, which has great potentials for future studies and urban climate assessments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Solid State Quantum Refrigeration Superconducting, Absorption and Measurement Based (Final Technical Report)

During this DOE grant, DE-SC0017890, in place for the past six years, all proposed research was carried out and published in peer-reviewed papers, as well as other projects that emerged during the research. In that effort the research team accomplished all proposed research, as well as many closely related research projects discovered and conceived of during the grant. These works included “Efficient Quantum Measurement Engines”, a work published in Physical Review Letters, giving a theory of quantum measurement-based engines, which uses quantum measurement as a resource. These engines are designed to efficiently convert energy from the stochastic quantum measurement process into useful work. Further publications include “Experimental Realization of a Quantum Dot Energy Harvester”, a joint theory and experimental work in collaboration with the group of Charles Smith in Cambridge, UK, as well as long time theoretical collaborators, Rafael Sánchez and Björn Sothmann. This work, featured as an Editor’s Suggestion in Physical Review Letters, realized an earlier theoretical proposal of ours, whereby two resonant tunneling quantum dots are connected to a central electronic cavity that is heated by a hot energy source. We also published “Superconducting Quantum Refrigerator: Breaking and Rejoining Cooper Pairs with Magnetic Field Cycles” a work done in collaboration with experimentalist Francesco Giazotto from ENS Pisa, Italy, which also resulted in a patent. This paper, published in Phys. Rev. Applied, advanced the concept of a cyclic fridge based on the normal/superconducting phase transition together with layered materials separated by tunnel junctions. We also completed the proposed research on a heat transistor, publishing “Thermal transistor and thermometer based on Coulomb-coupled conductors”, carried out as a collaboration between my group and theorists Splettstoesser (Lund U., Sweden), Sothmann (U. Duisburg-Essen, Germany), and Sánchez (U. Autónoma de Madrid, Spain). We carried out an analysis of a quantum coupled to a quantum point contact as a sensitive thermometer and heat transistor. We found the optimal statistical estimator for the temperature and compared it with experiments on the same type of devices. We also investigated autonomous quantum absorption refrigerators using quantum dots to cool by using a very hot thermal reservoir to drive heat between two other reservoirs. In the article “Quantifying the quantum heat contribution from a driven superconducting circuit”, we demonstrated that for a driven superconducting circuit, we showed heat flow provided by a hot source to the qubit can be switched on and off by varying external parameters, the frequency and the intensity of the driving. In the work “Stochastic thermodynamic cycles of a mesoscopic thermoelectric engine”, we reconsidered the autonomous thermoelectric heat engine in terms of underlying cycles. Rather than periodic behavior, the cycles were stochastic in nature. Nevertheless, by undertaking a graph theoretical analysis of the elementary transport processed, great quantitative and qualitative insight could be found. We also considered the quantum measurement process and showed that a quantum version of Maxwell’s demon could be related to the work extraction of a quantum system, closely related to arrow-of-time measures for quantum measurement, as described in our article “Thermodynamics of quantum measurement and Maxwell's demon's arrow of time”. This work was selected in Phys. Rev. A as an Editor’s Suggestion. A recent preprint titled “Cyclic Superconducting Quantum Refrigerators Using Guided Fluxon Propagation” accomplished an important piece of this grant: to propose a new kind of quantum refrigerator using the dynamics of fluxons in a type II superconductor. This invention envisioned a race-track type geometry where fluxons are confined. By applying a gradient of magnetic field together with electric current in a Corbino geometry, the circulating fluxons can actively cool a cold reservoir, realizing a new type of cyclic superconducting refrigerator. We also investigated the possibility of thermal control from different points of view. The application of quantum measurement to the system gives a new kind of control on the system of interest – we have pioneered this approach and shown that measurement can boost the thermal power of quantum engines as described in “Continuous measurement boosted adiabatic quantum thermal machines”. The ability to have heat flows on demand is an outstanding challenge, and we have provided new solutions to this problem in Thermal control across a chain of electronic nanocavities” for a chain of electron cavities using gating voltage control. The control methods using qubit/qubit coupling to create absorption fridges at their most fundamental level have also been developed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Joint Inversion of Surface Electrical Resistivity Tomography and Seismic Refraction Data between the 200 Areas

Geologic stratigraphy on the Hanford Site influences groundwater and contaminant migration through the aquifer system and the vadose zone. The current geologic framework model (GFM) relies heavily on a sparse distribution of borehole data in some locations to map geologic contacts and hydrologic properties in the subsurface. Non-invasive geophysical methods such as electrical resistivity tomography (ERT), transient electromagnetic surveying, and seismic imaging are being used at Hanford to map subsurface structure in areas with limited well observations. This is to develop and mature the capability of geophysical methods to aid in GFM refinement, to identify regions of subsurface complexity, and for optimal well siting. A joint inversion of co-located seismic refraction and ERT data was carried out for data collected on a ~2.3-km profile between the 200 Areas on the Hanford Site. While ERT and seismic refraction images have sensitivity to overlapping physical properties (porosity, moisture content, lithology), the resolution and physics used to acquire each of these datasets are different and therefore information can be different or mutually complementary. Performing a joint inversion provides a reasonable option for a coherent, coupled interpretation for mutually complementary datasets. Between the 200 Areas, there are few boreholes to interpret the geologic framework model, and these data sets were obtained to provide a first line of evidence toward identifying stratigraphic structure. The seismic refraction and ERT data were independently inverted during fiscal year 2022 and broadly showed a two-layer structure with a trough-like feature that is ~1 km wide and upwards of 150 m deep. The depth of the trough feature was greater in the ERT image compared to the seismic image, which indicated a maximum depth of approximately 110 m. The objective of the joint inversion described in this report was to invert the seismic refraction and ERT data together while constraining the ERT image to be structurally similar to the seismic refraction image. The approach was applied using the geophysical inverse modeling program E4D, which has the capability to invert first-arrival times from seismic refraction data and ERT resistances using a “cross-gradient” constraint. The application of cross-gradient constraints with different weights produces ERT models that show a high degree of similarity within the upper 100 m (above ~120 m elevation). None of the ERT models show an improved structural similarity to the seismic result; therefore, it is recommended that further attempts to jointly interpret these models focus on petrophysics and image resolution. Petrophysical measurements of core samples would improve knowledge of what drives the ERT response in this region and, along with downhole geophysical measurements, could be used to “ground truth” the surface-based geophysical results. Image resolution studies would provide insight into which regions of the inverted images are reliable and which regions are poorly constrained.

58 GEOSCIENCES↗
DOE OSTI2024

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗
NASA NTRS2003

Boundary-Layer Separation Control under Low-Pressure Turbine Airfoil Conditions using Glow-Discharge Plasma Actuators

Modem low-pressure turbines, in general, utilize highly loaded airfoils in an effort to improve efficiency and to lower the number of airfoils needed. Typically, the airfoil boundary layers are turbulent and fully attached at takeoff conditions, whereas a substantial fraction of the boundary layers on the airfoils may be transitional at cruise conditions due to the change of density with altitude. The strong adverse pressure gradients on the suction side of these airfoils can lead to boundary-layer separation at the latter low Reynolds number conditions. Large separation bubbles, particularly those which fail to reattach, cause a significant degradation of engine efficiency. A component efficiency drop of the order 2% may occur between takeoff and cruise conditions for large commercial transport engines and could be as large as 7% for smaller engines at higher altitude. An efficient means of of separation elimination/reduction is, therefore, crucial to improved turbine design. Because the large change in the Reynolds number from takeoff to cruise leads to a distinct change in the airfoil flow physics, a separation control strategy intended for cruise conditions will need to be carefully constructed so as to incur minimum impact/penalty at takeoff. A complicating factor, but also a potential advantage in the quest for an efficient strategy, is the intricate interplay between separation and transition for the situation at hand. Volino gives a comprehensive discussion of several recent studies on transition and separation under low-pressure-turbine conditions, among them one in the present facility. Transition may begin before or after separation, depending on the Reynolds number and other flow conditions. If the transition occurs early in the boundary layer then separation may be reduced or completely eliminated. Transition in the shear layer of a separation bubble can lead to rapid reattachment. This suggests using control mechanisms to trigger and enhance early transition. Gad-el-Hak provides a review of various techniques for flow control in general and Volino discusses recent studies on separation control under low-pressure-turbine conditions utilizing passive as well as active devices. As pointed out by Volino, passive devices optimized for separation control at low Reynolds numbers tend to increase losses at high Reynolds numbers, Active devices have the attractive feature that they can be utilized only in operational regimes where they are needed and when turned off would not affect the flow. The focus in the present paper is an experimental Separation is induced on a flat plate installed in a closed-circuit wind tunnel by a shaped insert on the opposite wall. The flow conditions represent flow over the suction surface of a modem low-pressure-turbine airfoil ('Pak-B'). The Reynolds number, based on wetted plate length and nominal exit velocity, is varied from 50,000 to 300,000, covering cruise to takeoff conditions. Low (0.2%) and high (2.5%) Gee-stream turbulence intensities are set using passive grids. A spanwise-oriented phased-plasma-array actuator, fabricated on a printed circuit board, is surface- flush-mounted upstream of the separation point and can provide forcing in a wide frequency range. Static surface pressure measurements and hot-wire anemometry of the base and controlled flows are performed and indicate that the glow-discharge plasma actuator is an effective device for separation control. of active separation control using glow discharge plasma actuators.

Hultgren, Lennart S.↗
DOE OSTI2020

Distributed Minimization of the Power Generation Cost in Prosumer-Based Distribution Networks: Preprint

Traditionally, electrical power was generated in big power plants. The cost of producing energy was related to the cost of fuel, e.g., carbon or gas, and by the cost of maintaining the power plants. With the advent of distributed energy resources, power can be produced directly at the edge of the electrical network by a new type of agents: the prosumers. Prosumers are entities that both consume and generate power, e.g., by means of photovoltaic panels. The cost of the power produced by prosumers is no longer related to fuel consumption since energy coming from distributed generators is essentially free. Rather, the cost is related to the remuneration that is due to the prosumers for the services they provide. The proposed control strategy minimizes the active power generation cost in the aforementioned scenario. The control scheme requires that the prosumers measure their voltage and then adjust the amount of injected power, according to a continuous time feedback control law that is indeed a projected gradient descent strategy. Simulations are provided in order to illustrate the algorithm behavior.

distributed control↗