Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coupled process modeling”

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,693 records · Page 94

Rheo-Structural Spectroscopy: Fingerprinting the In Situ Response of Fluids to Arbitrary Flow Fields

The objectives of this project were to develop new sample environments, measurement methodologies and associated modeling tools for characterizing the structural response to arbitrarily complex processing flows using small angle scattering, and to apply these new tools for understanding the fundamental physics governing the structuring of anisotropic particulate and polymeric materials under flow histories and conditions relevant to industrial processing flows. The research resulted in the development and implementation of a new sample environment, the fluidic four roll mill (FFoRM), for in situ small angle neutron and X-ray scattering (SANS/SAXS) measurements. These measurements are capable of generating large data sets that “fingerprint” how a complex fluid responds to a wide range of flow histories involving time variations in deformation type and rate. New modeling tools were developed to extract detailed microstructural information from such data sets, including orientation distribution functions and interparticle correlation functions, as well as reduced-order parametric descriptors of these high-dimensional functions that can be used to readily map, visualize and interpret a fluid’s structural response to its flow history. These new tools were applied to a range of model materials involving elongated particle suspensions in order to provide new insights into the physics of how flow couples with orientational and structural order in complex flows, particularly under non-dilute conditions for which no accurate theories currently exist. Using these investigations, we elucidated a number of new insights into the fundamental phenomena driving such process-structure-property relationships. These findings provide guidance for the further development of rheological models, and ultimately can inform the rational and model-based design of flow processes to achieve optimized orientational ordering that is key to the properties and function of a wide range of energy-relevant materials.

36 MATERIALS SCIENCE↗

Collisional-Radiative Nonequilibrium and Precursor Effects in a Nitrogen Shock Wave

Improvements to a plasma code with a Collisional-Radiative (CR) non-equilibrium model are made, allowing for a more accurate description of the physical processes. The code allows for non-Boltzmann distributions of the electronic excited states by convecting separately each excited state, as a pseudo-specie. Each molecular state has also its own vibrational temperature, while a global rotational temperature is assumed. The free electron temperature is different from those of the excited states, and the electron heat conduction is also included. The CR model also uses a unique coupling between chemistry and vibrational energy (C-V coupling), which is fully coherent, and has the property of establishing thermal equilibrium as well as chemical equilibrium, on its own. We have also included a coupling between electronic excitations and vibrational energy (X-V coupling), which can have a strong influence on the vibrational temperature of some states. The recent improvements include the multi- temperature dependence of the chemical rates for associative ionization, as well as the estimation of the internal energies transferred during this process. Additionally, the distribution of energy into different translational modes (electron and heavy particles) is now correctly modeled. This provides a very rapid heating mechanism for the free electrons, since it is found that the electrons are generated with an average thermal energy of the same order as the heavy particle translational energy. This effect was observed by Gorelov et al in a recent paper, and lead to pronounced peaks in electron temperature immediately behind the shock. We will attempt ro reproduce this phenomenon. The last modification concerns the inclusion of the radiative terms into the calculations, thus enabling us to observe the effect of radiative losses and radiation transport. Preliminary tests have shown that the radiative losses are not negligible, i.e. the shock velocity drops when the radiative emission is included. In addition, we have observed some precursor excitation and ionization phenomena. These may be important at high shock velocities, and this will be investigated in detail. Precursor phenomena were also observed (in air) by Gorelov et al, and our results (for pure nitrogen) will be qualitatively compared with theirs.

Cambier, Jean-Luc↗

Resource-Efficient Chemistry on Quantum Computers with the Variational Quantum Eigensolver and The Double Unitary Coupled-Cluster approach

Applications of quantum simulation algorithms to obtain electronic energies of molecules on noisy intermediate-scale quantum (NISQ) devices require careful consideration of resources describing the complex electron correlation effects. In modeling second-quantized problems, the biggest challenge confronted is that the number of qubits scales linearly with the size of molecular basis. This poses a significant limitation on the size of the basis sets and the number of correlated electrons included in quantum simulations of chemical processes. To address this issue and to enable more realistic simulations on NISQ computers, we employ the double unitary coupled-cluster (DUCC) method to effectively downfold correlation effects into the reduced-size orbital space, commonly referred to as the active space. Using downfolding techniques, we demonstrate that properly constructed effective Hamiltonians can capture the effect of the whole orbital space in small-size active spaces. Combining the downfolding pre-processing technique with the Variational Quantum Eigensolver, we solve for the ground-state energy of H2 and Li2 in the cc-pVTZ basis using the DUCC-reduced active spaces. We compare these results to full configuration-interaction and high-level coupled-cluster reference calculations.

quantum computing, variational quantum solver, cou↗

A coupled hydrodynamic (HEC-RAS 2D) and water quality model (WASP) for simulating flood-induced soil, sediment, and contaminant transport

Increased intensity and frequency of floods raise concerns about the release and transport of contaminated soil and sediment to and from rivers and streams. To model these processes during flooding events, we developed an External Coupler in Python to link the Hydrologic Engineering Center-River Analysis System (HEC-RAS) 2D hydrodynamic model to the Water Quality Analysis Simulation Program (WASP). Accurate data transfer from a hydrodynamic model to a water quality model is critical. Our test results showed the External Coupler successfully linked HEC-RAS and WASP and addressed technical challenges in aggregating flow data and conserving mass during the flood event. We ran the coupled models for a 100-year flood event to calculate flood-induced transport of sediment-associated arsenic in Woodbridge Creek, NJ. Change in surface sediment and arsenic at the end of 48-h flood simulation ranged from a net loss of 13.5 cm to a net gain of 11.6 cm, and 16.2 to 2.9 mg/kg, respectively, per model segment, which demonstrates the capability of the coupled model for simulating sediment and contaminant transport in flood.

54 ENVIRONMENTAL SCIENCES↗

Quantification of uncertainties in coupled material degradation processes - High temperature, fatigue and creep

This paper describes the development of methodology that provides for quantification of uncertainties in lifetime strength of aerospace materials subjected to a number of diverse effects. A probabilistic material degradation model, in the form of a randomized multifactor interaction equation, has been postulated for lifetime strength degradation of structural components of aerospace propulsion systems. The model includes effects that typically reduce lifetime strength and may include temperature, mechanical fatigue, creep and others. The paper also includes the analysis of experimental data from the open literature for Inconel 718. These data are used to provide an initial check for model validity, as well as for calibration of the model's empirical material constants. The model validity check and calibration is carried out for three effects, namely, high temperature, mechanical fatigue and creep.

Boyce, L.↗

On the connection between continental-scale land surface processes and the tropical climate in a coupled ocean-atmosphere-land system

The impact of global tropical climate to perturbations in land surface processes (LSP) are evaluated using perturbations given by different LSP representations of continental-scale in a global climate model that includes atmosphere-ocean interactions. One representation is a simple land scheme, which specifies climatological albedos and soil moisture availability. The other representation is the more comprehensive Simplified Simple Biosphere Model, which allows for interactive soil moisture and vegetation biophysical processes. The results demonstrate that LSP processes such as interactive soil moisture and vegetation biophysical processes have strong impacts on the seasonal mean states and seasonal cycles of global precipitation, clouds, and surface air temperature. The impact is especially significant over the tropical Pacific. To explore the mechanisms for such impact, different LSP representations are confined to selected continental-scale regions where strong interactions of climate-vegetation biophysical processes are present. We find that the largest impact is mainly from LSP perturbations over the tropical African continent. The impact is through anomalous convective heating in tropical Africa due to changes in the surface heat fluxes, which in turn affect basinwide teleconnections in the Pacific through equatorial wave dynamics. The modifications in the equatorial Pacific climate are further enhanced by strong air-sea coupling between surface wind stress and upwelling, as well as effect of ocean memory. Our results further suggest that correct representations of land surface processes, land use change and the associated changes in the deep convection over tropical Africa are crucial to reducing the uncertainty when performing future climate projections under different climate change scenarios.

Ma, Hsi-Yen↗

First-Principles Statistical Mechanics Study of Magnetic Fluctuations and Order–Disorder in the Spinel LiNi 0.5 Mn 1.5 O 4 Cathode

While significant magnetic interactions exist in lithium transition metal oxides, commonly used as Li-ion cathodes, the interplay between magnetic couplings, disorder, and redox processes remains poorly understood. In this work, we focus on the high-voltage spinel LiNi 0.5 Mn 1.5 O 4 (LNMO) cathode as a model system on which to apply a computational framework that uses first principles-based statistical mechanics methods to predict the finite temperature magnetic properties of materials and provide insights into the complex interplay between magnetic and chemical degrees of freedom. Density functional theory calculations on multiple distinct Ni–Mn orderings within the LNMO system, including the ordered ground-state structure (space group P4332), reveal a preference for a ferrimagnetic arrangement of the Ni and Mn sublattices due to strong antiferromagnetic superexchange interactions between neighboring Mn 4+ and Ni 2+ ions and ferromagnetic Mn–Mn and Ni–Ni couplings, as revealed by magnetic cluster expansions. These results are consistent with qualitative predictions using the Goodenough-Kanamori-Anderson rules. Simulations of the finite temperature magnetic properties of LNMO are conducted using Metropolis Monte Carlo. We find that a “semiclassical” Monte Carlo sampling method based on the Heisenberg Hamiltonian accurately predicts experimental magnetic transition temperatures observed in magnetometry measurements. This study highlights the importance of a robust computational toolkit that accurately captures the complex chemomagnetic interactions and predicts finite temperature magnetic behavior to help analyze experimental magnetic and magnetic resonance spectroscopy data acquired ex situ and operando.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-level theoretical characterization of the vinoxy radical ( • CH 2 CHO) + O 2 reaction

Numerous processes in atmospheric and combustion chemistry produce the vinoxy radical ( • CH 2 CHO). To understand the fate of this radical and to provide reliable energies needed for kinetic modeling of such processes, we have examined its reaction with O 2 using highly reliable theoretical methods. Utilizing the focal point approach, the energetics of this reaction and subsequent reactions were obtained using coupled-cluster theory with single, double, and perturbative triple excitations [CCSD(T)] extrapolated to the complete basis set limit. These extrapolated energies were appended with several corrections including a treatment of full triples and connected quadruple excitations, i.e., CCSDT(Q). In addition, this study models the initial vinoxy radical + O 2 reaction for the first time with multireference methods. We predict a barrier for this reaction of approximately 0.4 kcal mol –1 . This result agrees with experimental findings but is in disagreement with previous theoretical studies. The vinoxy radical + O 2 reaction produces a 2-oxoethylperoxy radical which can undergo a number of unimolecular reactions. Abstraction of a β-hydrogen (a 1,4-hydrogen shift) and dissociation back to reactants are predicted to be competitive to each other due to their similar barriers of 21.2 and 22.3 kcal mol –1 , respectively. The minimum-energy β-hydrogen abstraction pathway produces a hydroperoxy radical (QOOH) that eventually decomposes to formaldehyde, CO, and • OH. Two other unimolecular reactions of the peroxy radical are α-hydrogen abstraction (38.7 kcal mol –1 barrier) and HO 2 • elimination (43.5 kcal mol –1 barrier). Furthermore, these pathways lead to glyoxal + • OH and ketene + HO 2 • formation, respectively, but they are expected to be uncompetitive due to their high barriers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling discontinuous dynamic recrystallization containing second phase particles in magnesium alloys utilizing phase field method

Pre-existing second phase particles can induce particle stimulated nucleation (PSN) and pinning effect on grain boundaries during DRX process of magnesium alloys. The interaction among pinning effect, PSN mechanism and strain-induced grain boundary migration (SIBM) nucleation mechanism imposes challenge in quantitative study of microstructure evolution. A phase field model that describes discontinuous dynamic recrystallization process considering second phase particles (PF-DDRX-SP) which mainly distribute along grain boundaries is developed. The second phase order parameter with diffusion interface is introduced to the free energy function. With considering the pinning effect of second phase particles on grain boundaries, the grain boundary energy couples the interaction between second phase and grains. The effect of second phase on dislocation density evolution is considered in the model, and the PSN and SIBM mechanism are both executed by limiting the nucleation position to the grain boundary and the second phase interface. For the validation, the initial topology consistent with initial microstructure is applied to the simulation of DRX kinetics and flow behavior of AZ80 magnesium alloy, and great reliability and accuracy of developed model is indicated. Furthermore, through the model, the microstructure evolution with different particle size and with different volume fraction of round particle during DRX process were simulated. Finally, the results show that particles with higher dispersion reduce the number of DRX nuclei, which indicate that second phase particles suppress the nucleation at grain boundaries because of the pinning effect. More importantly, the proposed model could also quantitatively predict the relationships between the parameters of second phase particles and DRX behaviors, and enable to optimize the initial second phase structure in a uniform grain structure during thermomechanical process.

36 MATERIALS SCIENCE↗

Energy Coupling Mechanisms and Scaling Behavior Associated with Laser Powder Bed Fusion Additive Manufacturing

In situ optical absorptivity measurements are carried out to clarify the physics of the laser-material interactions involved and to validate both finite element and analytical models describing laser powder bed fusion processing. Absorption of laser energy is evaluated directly using precise calorimetry measurements and compared to melt pool depths for common structural metal alloys (Ti-6Al-4V, Inconel 625, and stainless steel 316L) as a function of incident laser power, scan velocity and laser beam diameter. Changes in both absorptivity and melt pool depths for all materials are found to vary strongly across the conduction-keyhole mode threshold. A hydrodynamic finite element model is coupled to a ray-tracing-based absorptivity model, yielding excellent agreement with the experimental results and elucidating additional physics. The experimental results are further analyzed using normalized enthalpy ( β ) and normalized thermal diffusion length ( L* th ) relations and demonstrate that the normalized melt pool depth ( d* = d/a , where d is melt pool depth and a is beam radius) is proportional to βL* th , while absorptivity follows an asymptotic exponential function against βL* th . Expressions for melt pool depth and laser absorptivity across different materials and laser scan systems are derived and thus provide useful tools to accelerate the optimization of laser processing parameters for metal 3D printing processes.

36 MATERIALS SCIENCE↗

Dynamics of the Transition Corona

Magnetic reconnection between the open and closed magnetic fields in the corona is believed to play a crucial role in the corona/heliosphere coupling. At large scale, the exchange of open/closed connectivity is expected to occur in pseudo-streamer (PS) structures. However, there is neither clear observational evidence of how such coupling occurs in PSs, nor evidence for how the magnetic reconnection evolves. Using a newly developed technique, we enhance the off-limb magnetic fine structures observed with the Atmospheric Imaging Assembly and identify a PS-like feature located close to the northern coronal hole. We first identify that the magnetic topology associated with the observation is a PS, null-point (NP) related topology bounded by the open field. By comparing the magnetic field configuration with the extreme ultraviolet (EUV) emission regions, we determined that most of the magnetic flux associated with plasma emission are small loops below the PS basic NP and open field bounding the PS topology. In order to interpret the evolution of the PS, we referred to a three-dimensional MHD interchange reconnection modeling the exchange of connectivity between small closed loops and the open field. The observed PS fine structures follow the dynamics of the magnetic field before and after reconnecting at the NP obtained by the interchange model. Moreover, the pattern of the EUV plasma emission is the same as the shape of the expected plasma emission location derived from the simulation. These morphological and dynamical similarities between the PS observations and the results from the simulation strongly suggest that the evolution of the PS, and in particular the opening/closing of the field, occurs via interchange/slipping reconnection at the basic NP of the PS. Besides identifying the mechanism at work in the large-scale coupling between the open and closed fields, our results highlight that interchange reconnection in PSs is a gradual physical process that differs from the impulsive reconnection of the solar-jet model.

magnetic fields↗

A numerical model of non-equilibrium thermal plasmas. II. Governing equations

Governing equations and the corresponding physical properties of the plasmas are both prerequisites for studying the fundamental processes in a non-equilibrium thermal plasma system numerically. In this paper, a kinetic derivation of the governing equations used for describing the complicated thermo-electro-magneto-hydrodynamic-chemical coupling effects in non-equilibrium thermal plasmas is presented. This derivation, which is achieved using the Chapman-Enskog method, is completely consistent with the theory of the transport properties reported in the previous paper by the same authors. It is shown, based on this self-consistent theory, that the definitions of the specific heat at constant pressure and the reactive thermal conductivity of two-temperature plasmas are not necessary. The governing equations can be reduced to their counterparts under local thermodynamic equilibrium (LTE) and local chemical equilibrium (LCE) conditions. The general method for the determination of the boundary conditions of the solved variables is also discussed briefly. The two papers establish a self-consistent physical-mathematical model that describes the complicated physical and chemical processes in a thermal plasma system for the cases both in LTE or LCE conditions and under non-equilibrium conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Visualization and prediction of supercritical CO 2 distribution in sandstones during drainage: An in situ synchrotron X-ray micro-computed tomography study

Pore-scale distribution of supercritical CO 2 (scCO 2 ) exerts significant control on a variety of key hydrologic as well as geochemical processes, including residual trapping and dissolution. Despite such importance, only a small number of experiments have directly characterized the three-dimensional distribution of scCO 2 in geologic materials during the invasion (drainage) process. Here, we present a study which couples dynamic high-resolution synchrotron X-ray micro-computed tomography imaging of a scCO 2 /brine system at in situ pressure/temperature conditions with quantitative pore-scale modeling to allow direct validation of a pore-scale description of scCO2 distribution. The experiment combines high-speed synchrotron radiography with tomography to characterize the brine saturated sample, the scCO 2 breakthrough process, and the partially saturated state of a sandstone sample from the Domengine Formation, a regionally extensive unit within the Sacramento Basin (California, USA). The availability of a 3D dataset allowed us to examine correlations between grains and pores morphometric parameters and the actual distribution of scCO 2 in the sample, including the examination of the role of small-scale sedimentary structure on CO2 distribution. The segmented scCO 2 /brine volume was also used to validate a simple computational model based on the local thickness concept, able to accurately simulate the distribution of scCO 2 after drainage. The same method was also used to simulate Hg capillary pressure curves with satisfactory results when compared to the measured ones. Finally, this predictive approach, requiring only a tomographic scan of the dry sample, proved to be an effective route for studying processes related to CO 2 invasion structure in geological samples at the pore scale.

3D quantitative morphometric analysis↗

An advanced joint inversion system for CO 2 storage modeling with large date sets for characterization and real-time monitoring-enhancing storage performance and reducing failure risks under uncertainties

As large-scale, commercial storage projects become operational, the problem of utilizing information from diverse sources becomes more critically important. In this project, we developed, tested, and applied an advanced joint data inversion system for CO 2 storage modeling with large data sets for use in site characterization and real-time monitoring. Emphasis was on the development of advanced and efficient computational algorithms for joint inversion of hydro-geophysical data, coupled with state-of-the-art forward process simulations. The developed system consists of (1) inversion tools using characterization data, such as 3D seismic survey (amplitude images), borehole log and core data, as well as hydraulic, tracer and thermal tests before CO 2 injection, (2) joint inversion tools for updating the geologic model with the distribution of rock properties, thus reducing uncertainty, using hydro-geophysical monitoring data, and (3) highly efficient algorithms for directly solving the dense or sparse linear algebra systems derived from the joint inversion. The system combines methods from stochastic analysis, fast linear algebra, and high performance computing. The developed joint inversion tools have been tested through synthetic CO 2 storage examples.

01 COAL, LIGNITE, AND PEAT↗

Precision masses for studies of the astrophysical r process

Half of the elements heavier than iron (Z > 26) are thought to be created through the astrophysical r process, whereby nuclides are produced via a rapid series of nuclear reactions that are postulated to occur in high temperature and neutron density environments such as supernovae or merging neutron stars. The nucleosynthetic path that describes the sequence of reactions through the chart of nuclides strongly depends on the neutron-separation energies of the nuclei. Until recently, however, almost all of these neutron-rich nuclei were not within reach of accelerator facilities, and therefore simulations of the r process had to rely on mass models for input into the calculations. Now, with the advent of facilities such as CARIBU at Argonne National Laboratory, the masses of many nuclides along the r-process path can be determined precisely with Penning trap mass spectrometers coupled to these facilities. More than 70 nuclides have been measured with the Canadian Penning Trap mass spectrometer alone in the past year, which overlap and complement results from other Penning trap mass spectrometers, and first calculations with these new masses suggest the timescale of the r process through the tin isotopes is delayed much more strongly than mass models would suggest.

Clark, Jason↗

Multi-Sensors Observations of Pre-Earthquake Signals. What We Learned from the Great Tohoku Earthquake?

The lessons learned from the Great Tohoku EQ (Japan, 2011) will affect our future observations and an analysis is the main focus of this presentation. Multi-sensors observations and multidisciplinary research is presented in our study of the phenomena preceding major earthquakes Our approach is based on a systematic analysis of several physical and environmental parameters, which been reported by others in connections with earthquake processes: thermal infrared radiation; temperature; concentration of electrons in the ionosphere; radon/ion activities; and atmospheric temperature/humidity [Ouzounov et al, 2011]. We used the Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) model, one of several possible paradigms [Pulinets and Ouzounov, 2011] to interpret our observations. We retrospectively analyzed the temporal and spatial variations of three different physical parameters characterizing the state of the atmosphere, ionosphere the ground surface several days before the March 11, 2011 M9 Tohoku earthquake Namely: (i) Outgoing Long wave Radiation (OLR) measured at the top of the atmosphere; (ii) Anomalous variations of ionospheric parameters revealed by multi-sensors observations; and (iii) The change in the foreshock sequence (rate, space and time); Our results show that on March 8th, 2011 a rapid increase of emitted infrared radiation was observed and an anomaly developed near the epicenter with largest value occurring on March 11 at 07.30 LT. The GPS/TEC data indicate an increase and variation in electron density reaching a maximum value on March 8. Starting from this day in the lower ionosphere there was also observed an abnormal TEC variation over the epicenter. From March 3 to 11 a large increase in electron concentration was recorded at all four Japanese ground-based ionosondes, which returned to normal after the main earthquake. We use the Japanese GPS network stations and method of Radio Tomography to study the spatiotemporal structure of ionospheric perturbations, and to distinguish ionospheric responses to processes of EQ preparation against the effects of other factors. The 2-D snapshots of the electron density over Japan showed abnormal increase over the maximum stress during the night, a few hours before the main shock. Our results from recording atmospheric and ionospheric conditions during the earthquake indicate the presence of anomalies in the atmosphere and ionospheres occurring consistently over regions of maximum stress near the epicenter. Due to their long duration (hours and days) and spatial appearance (only over the Sendai region) these results do not appear to be caused by meteorological or magnetic activity. They reveal the existence of atmospheric and ionospheric phenomena occurring prior to the earthquake, which indicates new evidence of a distinct coupling between the lithosphere and atmosphere/ionosphere. Similar results have been reported before the catastrophic events in Chile (M8.8, 2010), Italy (M6.3, 2009) and Sumatra (M9.3, 2004).

Ouzonounov, D.↗

From cells to laminate: probing and modeling residual stress evolution in thin silicon photovoltaic modules using synchrotron X-ray micro-diffraction experiments and finite element simulations

Fracture of silicon crystalline solar cells has recently been observed in increasing percentages especially in solar photovoltaic (PV) modules involving thinner silicon solar cells (<200 μm). Many failures due to fracture have been reported from the field because of environmental loading (snow, wind, etc.) as well as mishandling of the solar PV modules (during installation, maintenance, etc.). However, a significantly higher number of failures have also been reported during module encapsulation (lamination) indicating high residual stress in the modules and thus more prone to cell cracking. Here in this paper we report through the use of synchrotron X-ray submicron diffraction coupled with physics-based finite element modeling, the complete residual stress evolution in mono-crystalline silicon solar cells during PV module integration process. For the first time, we unravel the reason for the high stress and cracking of silicon cells near soldered inter-connects. Our experiments revealed a significant increase of residual stress in the silicon cell near the solder joint after lamination. Moreover, our finite element simulations show that this increase of stress during lamination is a result of highly localized bending of the cell near the soldered inter-connects. Further, the synchrotron X-ray submicron diffraction has proven to be a very effective way to quantitatively probe mechanical stress in encapsulated silicon solar cells. Thus, this technique has ultimately enabled these findings leading to the enlightening of the role of soldering and encapsulation processes on the cell residual stress. This model can be further used to suggest methodologies that could lead to lower stress in encapsulated silicon solar cells, which are the subjects of our continued investigations.

14 SOLAR ENERGY↗