Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density estimation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Hydraulic Properties Data Package for the E-Area Soils, Cementitious Materials, and Waste Zones–Update

Hydraulic property estimates for the soils, cementitious materials, and waste zones associated with the E-Area low-level radioactive waste disposal units have been provided to support the Performance Assessments (PA) for the E-Area Low-Level Waste Facility (LLWF). Nominal or “best estimate” hydraulic property values for use in the deterministic modeling are provided along with representations of the hydraulic property value uncertainty for use in sensitivity and uncertainty modeling. The hydraulic properties provided for each of the E-Area materials include: porosity ( η ), dry bulk density ( ρ b ), particle density (ρ p ), saturated hydraulic conductivity ( K sat ), characteristic curves (suction head, saturation, and relative permeability), and effective diffusion coefficient ( D e ). A representation of the uncertainty associated with each property, except for the characteristic curves, is provided for each material, except for the E-Area waste zones. These nominal parameter values for each of the E-Area soils, cementitious materials, and waste zones, and, where indicated, parameter uncertainty representations, are based upon the following in order of priority: Site-specific field data, Site-specific laboratory data, Similarity to material with site-specific field or laboratory data, and Literature data. Additionally, a methodology to represent long-term concrete degradation is provided. Finally, because much of the nominal hydraulic property values and uncertainty representations for the E-Area soils, cementitious materials, and waste zones are based on similarity to other materials or literature data, a methodology to prioritize additional work to better define these values and representations is outlined. Prioritization should be based on the importance of the material and/or property to the results of deterministic, sensitivity, and uncertainty modeling. This prioritization should be established through a process of sensitivity modeling. This report supersedes the previous hydraulic properties data package. Revision 1 contains clarifications, updates and corrections to text, figures and tables in a number of places in the report along with some minor editing.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Uncertainty quantification for equations of state: copper as an example

Equations of state are essential for providing a fundamental description of materials properties in thermodynamic equilibrium and are used to provide closure relations for hydrodynamics simulations. Generally, equations of state rely on simple physics-based parameterized materials models to inform on the free energy of a material through out a given thermodynamic state space. Historically the parameters of these models have been tuned by hand to fit various experimental data. However, modern optimization and uncertainty quantification techniques allow us to quickly test thousands of parameter combinations and obtain meaningful uncertainty estimates on the parameters, opening opportunities for assessing systematic uncertainties in experiments, assessing model adequacy, and more. In this report, we use Bayesian inference to fit the solid (fcc) equation of state of copper. We focus on fitting five different experimental datasets, including the isobaric density, isobaric heat capacity, room temperature isotherm, principal isentrope, and principal Hugoniot. We fit all five data types simultaneously, and then explore the extent to which combinations of 2 subsets of the 5 datasets can constrain the EOS parameters, as compared to the fit to all 5. This information is useful for investigating the extent to which different datasets can con strain EOS models and thereby help guide experimental investigations in order to best constrain the EOS. We also discuss ways that the methodologies can be used to investigate systematic discrepancies between experiments, as well as how the methods can be used to assess model uncertainty. The framework we develop is general, in that it can be used with a variety of optimization or uncertainty quantification techniques and with a variety of data sources, including both experimental and ab-inito data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Irradiated Microstructure of Zorita Materials

Reactor baffle plate materials were harvested from the decommissioned Zorita nuclear power reactor after 38 years of service and were characterized with transmission electron microscope. The samples were Type 304 stainless steels and had been irradiated to 0.05, 16.9, ad 48.6 dpa. Faulted dislocation loops, voids, and irradiation-induced precipitates were observed in this material. While dislocation loops were the dominant irradiation defects for all samples, small voids and precipitates were also observed at 16.9 and 48.6 dpa, but not at 0.05 dpa. The size of dislocation loops increased with irradiation dose and saturated below 16.9 dpa. The loop density continued to increase between 16.9 and 48.6 dpa. For voids, both the density and mean size were small, and the volume swelling was negligibly small even at 48.6 dpa. In addition, carbides or G-phase precipitates were detected at 16.9 and 48.6 dpa, and a complex dose dependence was observed. A range of irradiation hardening was estimated at each dose, and the resulting yield strengths were consistent with the previous tensile results measured from the Halden and BOR-60 specimens.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

New Thermonuclear Rate of 7 Li(d,n)2 4 He Relevant to the Cosmological Lithium Problem

Accurate 7 Li(d,n)2 4 He thermonuclear reaction rates are crucial for precise prediction of the primordial abundances of lithium and beryllium and to probe the mysteries beyond fundamental physics and the standard cosmological model. However, uncertainties still exist in current reaction rates of 7 Li(d,n)2 4 He widely used in big bang nucleosynthesis (BBN) simulations. In this work, we reevaluate the 7 Li(d,n)2 4 He reaction rate using the latest data on the three near-threshold 9 Be excited states from experimental measurements. We present for the first time uncertainties that are directly constrained by experiments. Additionally, we take into account for the first time the contribution from the subthreshold resonance at 16.671 MeV of 9 Be. We obtain a 7 Li(d,n)2 4 He rate that is overall smaller than the previous estimation by about a factor of 60 at the typical temperature of the onset of primordial nucleosynthesis. We implemented our new rate in BBN calculations, and we show that the new rates have a very limited impact on the final light element abundances in uniform density models. Typical abundance variations are in the order of 0.002%. For nonuniform density BBN models, the predicted 7 Li production can be increased by 10% and the primordial production of light nuclides with mass number A > 7 can be increased by about 40%. Our results confirm that the cosmological lithium problem remains a long-standing unresolved puzzle from the standpoint of nuclear physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Design and Development of Novel Equiatomic Refractory Multi Principal elemental Alloys Based on MoNbTi System for Use in Irradiation Environments

Multi Principal Elemental Alloys (MPEA) have emerged as promising materials for next-generation nuclear reactors due to their exceptional irradiation resistance. Eight equiatomic MPEA based on the MoNbTi system, comprising of elements with low thermal neutron absorption cross-sections were explored using a combined approach employing empirical parameter estimations, and CALPHAD simulations by which the phases and elemental segregation observed in all the alloys in the as-cast state were predicted. Solution heat treatment at 1500°C transformed five alloys into single-phase matrix materials, enhancing homogeneity and reducing hardness. The densities of the alloys ranged between 6.47 to 7.68 g/cm3, hardness between 472 and 656 VHN, Young’s modulus between 142 GPa to 169 GPa, shear modulus between 54 GPa and 62 GPa, bulk modulus between 117 GPa to 194 GPa and Poisson’s ratio between 0.3 to 0.35. The in-situ high temperature Xray diffraction results, differential scanning calorimetry and dilatometry results up to 1000°C suggested the high temperature phase stability of the MPEA. Subsequent ageing heat treatment at 800 and 1000 oC for 96 hours revealed significant secondary phase precipitation in MoNbTiZr, MoNbTiZrV, and MoNbTiCrA. Oxidation studies at 800 oC for 24 hours in air revealed superior oxidation resistance and cubic rate law dependence in Cr containing MPEA especially in MoNbTiCrAl, while severe mass gain resulting in total disintegration and exfoliation in MoNbTiZr and MoNbTiZrV. This comprehensive study underscores the potential of novel MPEA as promising materials for advanced nuclear reactor applications, shedding light on their microstructural control, mechanical properties, thermal stability, and oxidation resistance.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Recent progress in the design of the K-DEMO divertor

The preliminary conceptual design of the Korean fusion demonstration reactor (K-DEMO) with a major radius of 6.8 m and the fusion power of 2200 MW has been studied since 2012. The overall configuration of the K-DEMO divetor system based on the ITER-like water-cooled tungsten technology is a double-null type symmetric divertor subdivided into 32 toroidal modules for the vertical maintenance. A detached divertor scenario with impurity seeding was considered as the primary approach for the power exhaust to reduce the peak heat flux lower than the engineering limit of 10 MW/m 2 . The power exhaust performance at the scrap off layer was estimated by using UEDGE-2D code, a two-dimensional fluid transport code for collisional edge plasma and neutral species like N, Ne, and Ar. Particle and heat flux on inboard and outboard divertor targets were calculated for the detached cased depending on parameters such as the impurity seeding rate, pumping rate, and the pedestal density. On the other hand, a magnetic solution like X-divertor, snowflake divertor, and super X-divertor to expand the plasma wet area was considered for K-DEMO since the detached divertor increasing a radiation fraction by impurity seeding might be able to be unstable. However, the extremely high current of poloidal coils was required more than the engineering limit, 20 MA, to form magnetic field lines for the alternative divertors. Based on the physical calculation of the edge plasma, engineering analyses were carried out to find out the thermal and structural reliability. The thermo-hydraulic analysis confirmed thermal stability, whether all comprising materials are operating within their allowable temperature windows when the case of the peak heat flux is set to 10 MW/m 2 on the outboard divertor target. The response surface optimization method derived two optimal design candidates employing two kinds of heat sink materials, respectively: the reduced activation ferritic martensitic (RAFM) steel and CuCrZr alloy. Additionally, the drawbacks and merits of the two materials were definite. The optimal design with applying RAFM steel was vulnerable to withstand thermal and mechanical loads since low thermal conductivity caused too thin thickness of the heat sink. On the other hand, the CuCrZr alloy has critical drawbacks in terms of activation and radioactive waste despite its high thermal conductivity. Meanwhile, preliminary electromagnetic (EM) analysis was carried out to estimate the EM loads caused by the abnormal behaviors of plasma since EM loads are one of the most critical external loads for designing a DEMO divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Combining Computational Modeling with Reaction Kinetics Experiments for Elucidating the In Situ Nature of the Active Site in Catalysis

Microkinetic modeling based on density functional theory (DFT) derived energetics is important for addressing fundamental questions in catalysis. The quantitative fidelity of microkinetic models (MKMs), however, is often insufficient to conclusively infer the mechanistic details of a specific catalytic system. Here, this can be attributed to a number of factors such as an incorrect model of the active site for which DFT calculations are performed, deficiencies in the hypothesized reaction mechanism, inadequate consideration of the surface environment under reaction conditions, and intrinsic errors in the DFT exchange-correlation functional. Despite these limitations, we aim at developing a rigorous understanding of the reaction mechanism and of the nature of the active site for heterogeneous catalytic chemistries under reaction conditions. By achieving parity between experimental and modeling outcomes through robust parameter estimation and by ensuring coverage-consistency between DFT calculations and MKM predictions, it is possible to systematically refine the mechanistic model and, thereby, our understanding of the catalytic active site in situ.

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

Extensions of a classical mechanics “piston-model” for understanding the impact of asymmetry on ICF implosions: The cases of mode 2, mode 2/1 coupling, time-dependent asymmetry, and the relationship to coast-time

As long suspected, low mode asymmetry in inertially confined fusion (ICF) implosions has been implicated as a performance limiting factor [Casey et al., “Evidence of three-dimensional asymmetries seeded by high-density carbon-ablator nonuniformity in experiments at the national ignition facility,” Phys. Rev. Lett. 126, 025002 (2021)]. Recently a non-linear, but solvable, theory [Hurricane et al., “An analytic asymmetric-piston model for the impact of mode-1 shell asymmetry on ICF implosions,” Phys. Plasmas 27, 062704 (2020)] based upon the simple picture of a pair of asymmetric pistons has generated new insights and provided some practical formulas for estimating the degradation of an implosion due to mode-1 asymmetry and demonstrated a previously unrecognized connection between measured hot-spot drift velocity, nuclear down-scatter ratio asymmetry, and the concept of residual kinetic energy (RKE). Asymmetry of the implosion “shell,” as opposed to asymmetry of the hot-spot, was key to the classical mechanics model because the majority of the kinetic energy in an implosion is carried by the shell. Herein, the two-piston model is extended to a six-piston model in order to capture mode-2 asymmetry and coupling between mode-1 and mode-2. A key result of this new six-piston model is that the weighted harmonic mean of shell areal density is the fundamental quantity that determines the RKE and performance degradations for a three-dimensional implosion. Agreement is found between the scalings coming from the theory and ICF implosion data from the National Ignition Facility and to large ensembles of detailed simulations. The connection between the piston model's dependence upon the radius of peak velocity and coast-time is also highlighted in this paper. Finally, by extending the two-piston model to include time-dependent “swing,” it is shown in the Appendix that the shell asymmetry at the time of stagnation dominates the solution for RKE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Development of Real-Time High-Density Pulsar Data Transmission and Processing for Grid Synchronization

Taking advantage of the extreme stability of the pulsar period, it can serve as the timing source for grid synchronization to compensate for the timing drift instigated by the loss of GPS signal. Nevertheless, the real-time transmission and processing of the pulsar data suffer from its high-frequency data rate, varying from megahertz to gigahertz, resulting in reduced computing speed and increased time delay. To mitigate this issue, the hardware and software frameworks are implemented for the high-density pulsar data transmission and processing for grid synchronization in this research. Initially, the high-density pulsar data is transferred using open-source software. The complementary duty cycle timing module is designed to coordinate the operation of the dual-channel high-speed interface and software. Subsequently, the multiple-threading is applied to the receiving, parsing, and splicing pulsar data. Next, the pulsar signal extraction method is implemented based on the polyphase filterbank and time of arrival estimation. Ultimately, real-time performance verification experiments are carried out for different components under two hardware platforms. Finally, the results demonstrate that only 0.482 s is required for processing 4 Gigabyte data through multiple-threading, which is 3.8 times faster than the single thread. The pulsar signal extraction can also be executed within 707 ms for 4.8 seconds of data, thereby indicating that real-time requirements can be met.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

TRINIDI (Time-of-Flight Resonance Imaging with Neutrons for Isotopic Density Inference)

This software is an open-source Python library that provides tools for processing hyperspectral neutron time-of-flight radiography data. This type of data allows material decomposed reconstructions to be generated with the use of material characteristic spectral responses and the algorithms provided in this code library. The software library will contain tools for pre-processing the neutron measurement data, estimating measurement system parameters, reconstructing material decomposed radiographs, and computing material decomposed computed tomography (CT). Furthermore, it will have capability to generate and process simulated neutron time-of-flight data with the goal of benchmarking and demonstrating the tools that are provided. The software will include thorough documentation and application examples.

Balke, Thilo↗
DOE OSTI2022

Large-Area (Ag,Cu)(In,Ga)Se 2 Thin-Film Solar Cells with Increased Bandgap and Reduced Voltage Losses Realized with Bulk Defect Reduction and Front-Grading of the Absorber Bandgap

The 1.24 eV bandgap, 18.8% power conversion efficiency Ag-alloyed chalcopyrite (Ag,Cu)(In,Ga)Se 2 (ACIGS) solar cells are characterized to relate voltage and efficiency improvements to electro-optical (EO) characteristics. Shockley-Read-Hall recombination center defect density, identified and characterized through deep level transient spectroscopy and time-resolved photoluminescence (TRPL), is reduced through potassium and copper treatment optimization. Concomitantly, longer minority carrier lifetimes are achieved, which increases open-circuit voltage (V OC ). Near-conduction band defects associated in earlier studies with light-induced current instability are also mitigated. Analysis of charge-carrier dynamics after single- and two-photon excitation is used to separate recombination at the front interface and in the absorber bulk. From TRPL decay simulations, the authors estimate ranges of key solar cell material characteristics: bulk carrier lifetime T bulk = 110-210 ns, charge-carrier mobility u = 110-160 cm 2 V -1 s -1 , and front interface recombination velocity S front = 700-1050 cm s -1 . This lowest-reported S front for ACIGS absorbers originates from the notched conduction band grading, which also makes the impact of the back interface recombination negligible. It is suggested in the results that solar cell performance enhancements can be made most readily with two distinct strategies: improving device architecture and reducing semiconductor defect densities. Using these approaches, power conversion efficiency in large-area solar cells is improved by 1.1% absolute.

14 SOLAR ENERGY↗
DOE OSTI2021

Bayesian automated weighting of aggregated DFT, MD, and experimental data for candidate thermodynamic models of aluminum with uncertainty quantification

Atomic-scale modeling methods such as density functional theory (DFT) and molecular dynamics (MD) can predict the thermodynamic properties of materials at a lower cost than experimental measurements. However, their regular usage in thermodynamic model construction is hampered by the lack of quantitative agreement with experimental measurements and the lack of uncertainty estimates on the data. To make regular usage of this atomistic simulation data, it is important to assess whether the atomistic simulation datasets, by themselves or in combination with experimental measurements, result in the same physics-informed models best supported by experimental measurements alone. Here, models of aluminum thermodynamic properties are discussed using three data sources: atomistic calculations (DFT and MD), experiments, and a combination of atomistic calculations and experiments. The study shows that, after ensuring self-consistency in predicting key invariant points, both experimental measurements and atomistic calculations can significantly contribute to an optimal model.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Cross-machine comparison of runaway electron generation during tokamak start-up for extrapolation to ITER

A cross-machine comparison of global parameters that determine the runaway electron (RE) generation and loss process during tokamak start-up was carried out with the aim to extrapolate these to ITER. The study found that all considered discharges, also those that do not show signs of RE, are non-thermal at the start, i.e. have a streaming parameter larger than 0.1. During the current ramp-up the electric field, E, remains above the critical value, E c , that allows RE in the plasma. The distinction to be made is not if RE can form but, if sufficient RE can form fast enough such that they are detected or start to dominate the dynamics of the tokamak discharge. The dynamics of the value of E, density and temperature during tokamak are key to the formation of RE. It was found that larger devices operate with E closer to E c , due to their higher temperatures, hence the RE generation is relatively slower. The slower time scales for the formation of RE, estimated to be of the order of 100s of ms in ITER simplifies the development of avoidance schemes. The RE confinement time is also an important determinant of the entire process and is found to increase with the device size. The study also revealed that drift orbit losses, a mechanism often attributed as the main RE loss mechanism during the early tokamak discharge, are actually more difficult to achieve. RE losses might be more likely attributed to RE diffusion due to magnetic turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Energy Savings Potential of SSL in Agricultural Applications

A report on the energy-savings potential of solid-state lighting in agricultural applications, which updates a prior edition that focused on horticultural lighting only, with an expanded focus that includes animal lighting. The analysis considered three indoor horticultural lighting categories: lighting supplemented greenhouses, high intensity sole-source farms, and sole-source indoor vertical farms typically growing leafy greens. To characterize animal lighting, the analysis investigated the three largest markets of livestock in the U.S. that use lighting – poultry, swine, and dairy farming. To estimate energy consumption of agricultural lighting, the analysis utilized interviews with lighting manufacturers, growers, utility companies, academic professionals, and industry consultants, as well as various data sources. Each of these resources enabled the determination of total market size and illuminated area, typical lighting configuration, lighting power density for LED and incumbent technologies, operating hours, and the installed lighting technology mix. Among the major findings: 1. If all indoor horticultural lighting today was converted to LED technology, annual horticultural lighting consumption would be reduced to 6.3 TWh of site electricity, or 60 tBtu of source energy, which represents lighting energy savings of 34% or $350 million. 2. If all animal lighting today was converted to LED technology, annual lighting energy consumption for animal production would be reduced to 2.8 TWh site electricity, or 26 tBtu of source energy, which represents lighting energy savings of 25% or $96 million.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

On the Diversity of Asymmetries in Gapped Protoplanetary Disks

Protoplanetary disks with large inner dust cavities are thought to host massive planetary or substellar companions. These disks show asymmetries and rings in the millimeter continuum caused by dust trapping in pressure bumps and potentially vortices or horseshoes. The origin of the asymmetries and their diversity remains unclear. We present a comprehensive study of 16 disks for which the gas surface density profile has been constrained by CO isotopologue data. First, we compare the azimuthal extents of the dust continuum profiles with the local gas surface density in each disk and find that the asymmetries correspond to higher Stokes numbers or low gas surface density. We discuss which asymmetric structures can be explained by a horseshoe, a vortex, or spiral density waves. Second, we reassess the gas gap radii from the {sup 13}CO maps, which are about a factor of 2 smaller than the dust ring radii, suggesting that the companions in these disks are in the brown dwarf (∼15–50 M {sub Jup}) or super-Jovian (∼3–15 M {sub Jup}) mass regime on eccentric orbits. This is consistent with the estimates from contrast curves on companion mass limits. These curves rule out (sub)stellar companions (q > 0.05) for the majority of the sample at the gap location, but it remains possible at even smaller radii. Third, we find that spiral arms in scattered-light images are primarily detected around high-luminosity stars with disks with wide gaps, which can be understood by the dependence of the spiral arm pitch angle on disk temperature and companion mass.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Generating Buoyancy-Driven Convection in Membrane Distillation

Membrane distillation (MD) is a thermally-driven desalination process that can treat hypersaline brines. Considerable MD literature has focused on mitigating temperature and concentration polarization. This literature largely neglects that temperature and concentration polarization increase the feed density near the membrane. With gravity properly oriented, this increase in density could trigger buoyancy-driven convection and increase permeate production. Convection could also be strengthened by heating the feed channel wall opposite the membrane. To investigate that possibility, we perform a series of experiments using a plate-and-frame direct contact MD system with an active membrane area of 300 cm2 and a feed channel wall heated using a resistive heater. The experiments measure the average transmembrane permeate flux for two gravitational orientations, feed Reynolds numbers between 128 and 1128, and wall heat fluxes up to 12 kW/m2. The results confirm that with gravity properly oriented, wall-heating can trigger buoyancy-driven convection for a wide range of feed Reynolds numbers, and increase permeate production between roughly 20 and 130 %. We estimate, however, that at high Reynolds numbers (Re > 800), more than 70 % of the wall heat is carried out of the MD system by the feed flow, without contributing to permeate production. This suggests the need for longer membranes and heat recovery steps in any future practical implementation.

buoyancy-driven convection↗
DOE OSTI2021

Evaluation of Phase Change Plaster/Paste Composites for Building Envelopes

Thermal energy storage (TES) materials increase a building’s thermal mass and provide the flexibility to shave and time-shift energy use. These materials are needed to enable grid-interactive efficient buildings, allowing building occupants to lower costs by leveraging free ambient heating/cooling as well as shifting consumption to cheaper electricity periods. Phase change materials (PCMs) that store heat through a solid-liquid/solid-solid transformation are of particular interest for TES because of their ability to store large amounts of energy across a constant temperature or very small temperature range. Furthermore, previous studies on PCM building envelopes evaluated products that may only be incorporated as part of new construction and are not easily applicable to existing, energy-intensive buildings. In this work, we leverage laboratory PCM characterization and whole-building energy simulation to study the energy and cost savings of PCM plaster/paste coatings, which can be applied post construction to both new builds and retrofits of existing buildings (about 40% of residential buildings in the United States fall into the retrofit category). We evaluated the energy impact of PCM plaster/paste composites for a single-family residential building in the United States, including layer design (location and thickness), composite properties (phase change temperature, energy density, and shape of the enthalpy curve), and U.S. climate zone. The enthalpy curve shape of PCM shows that phase change temperature and energy density are not the only two parameters that must be considered for PCM selection. Maximizing the energy density amid the effective TES zone is key to boosting energy savings. We also studied PCM price targets with the aim of achieving reasonable payback periods based on both block and time-of-use electricity rates. Moreover, we developed an empirical equation to provide a rough estimate of the energy savings potential of incorporating PCM without the need to perform detailed simulations. These results will help guide the selection and design of PCM plaster/paste composite elements for residential building envelopes.

ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATION↗