Search NASA⌕ Search

SEARCH · Search NASA

Results for “Surface Operations”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A comprehensive conductivity model for drift and micro-tearing modes

The parallel electrical conductivity is a crucial parameter in the study of the linear stability of drift-modes like the micro-tearing mode (MTM). The conductivity enters by closing the electromagnetic tearing layer equations. Recent progress in the understanding of the pedestal suggests that the MTM could play an important role in its structure and evolution. For this reason, we revisit and improve previous model conductivities. This parameter is generally derived from the linearized drift kinetic equation. In the past literature, it has been computed using either simplified collision operators or neglecting the spatial dependence away from the rational surface. Here, a fully consistent expression for the conductivity that would accurately model the pedestal has not been available. By applying a novel variational procedure and with the full Fokker Plank collision operator, including electron–electron collisions, we compute a closed expression for the parallel electrical conductivity in the form of a rational function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Ultrathin Microporous Transport Layers: Implications for Low Catalyst Loadings, Thin Membranes, and High Current Density Operation for Proton Exchange Membrane Electrolysis

Porous transport layers (PTL) and their surface properties have the potential to improve the performance of proton exchange membrane water electrolyzers (PEMWE), which is imperative to reduce feedstock costs and lead to their widespread implementation. This work introduces a novel generation of titanium microporous layers (MPLs) with ultra-low thicknesses of approx. 20 um which reduces raw material costs. They also feature advanced interfacial properties tailored to maximize catalyst utilization at low Ir-loadings. The bulk morphology and surface properties of the hierarchically structured PTLs are assessed by X-ray tomographic microscopy. The low surface roughness of the MPL allows the use of thinner membranes since it minimizes possible deformations in the membrane. Cells containing the MPLs outperformed those containing state-of-the-art commercially available PTL materials by up to 100 mV at 7 A cm-2 in combination with low-loaded catalyst-coated membranes of 0.4 mgIr cm-2. Hydrogen crossover is also reduced, especially at low current densities, leading to a larger turndown ratio which can enable more cost-effective operating strategies. Finally, these rationally designed MPLs also lead to high catalyst utilization by overcoming the naturally occurring high in-plane resistance of low-loaded catalyst layers.

hydrogen crossover↗
DOE OSTI2020

Extended Interfacial Stability through Simple Acid Rinsing in a Li-Rich Oxide Cathode Material

Layered Li-rich Ni, Mn, Co (NMC) oxide cathodes in Li-ion batteries provide high specific capacities (>250 mAh/g) via O-redox at high voltages. However, associated high-voltage interfacial degradation processes require strategies for effective electrode surface passivation. Here, we show that an acidic surface treatment of a Li-rich NMC layered oxide cathode material leads to a substantial suppression of CO 2 and O 2 evolution, ~90% and ~100% respectively, during the first charge up to 4.8 V vs. Li +/0 . CO 2 suppression is related to Li 2 CO 3 removal as well as effective surface passivation against electrolyte degradation. This treatment does not result in any loss of discharge capacity and provides superior long-term cycling and rate performance compared to as-received, untreated materials. We also quantify the extent of lattice oxygen participation in charge compensation (“O-redox”) during Li + removal by a novel ex-situ acid titration. Our results indicate that the peroxo-like species resulting from O-redox originate on the surface at least 300 mV earlier than the activation plateau region around 4.5 V. X-ray photoelectron spectra and Mn -L X-ray absorption spectra of the cathode powders reveal a Li + deficiency and a partial reduction of Mn ions on the surface of the acid-treated material. More interestingly, although the irreversible oxygen evolution is greatly suppressed through the surface treatment, our O K-edge resonant inelastic X-ray scattering shows the lattice O-redox behavior largely sustained. The acidic treatment, therefore, only optimizes the surface of the Li-rich material and almost eliminates the irreversible gas evolution, leading to improved cycling and rate performance. This work therefore presents a simple yet effective approach to passivate cathode surfaces against interfacial instabilities during high-voltage battery operation.

25 ENERGY STORAGE↗
DOE OSTI2023

Kinetically versus thermodynamically controlled factors governing elementary pathways of GaP(111) surface oxidation

GaP and related III-V semiconductors have attracted interest as high-efficiency photoelectrochemical electrodes for solar water splitting. However, their efficacy is linked to the presence, identity, and integrity of native surface oxides, which are structurally and chemically complex and evolve during operation. Using ambient pressure X-ray photoelectron spectroscopy (APXPS) coupled with ab initio simulations, we track key chemical motifs expressed during evolution of GaP(111) surface oxides, their associated reaction kinetics, and their correlations with electronic properties. We identify two distinct thermal regimes corresponding to kinetically and thermodynamically controlled oxidation. Below 600 K, exposure to O2 generates kinetically facile Ga-O-Ga configurations, whereas higher temperatures cause activated oxygen to insert into Ga–P bonds as part of a thermodynamically driven transformation into a complex, heterogeneous 3D network of surface POx (1≤x ≤ 4) groups and Ga2O3 species, the latter of which eventually dominates upon depletion of surface phosphorus. Here our study highlights the critical competition between kinetic and thermodynamic factors during GaP oxidation, yielding insights for fabricating stable III-P-based photoelectrodes with precisely engineered surface properties.

08 HYDROGEN↗
DOE OSTI2021

Copper Dendrites and Surface Engineering for Enhanced CO2 Reduction Research Report Paper

The electrocatalytic reduction of CO 2 (CO 2 R) into hydrocarbon fuels, such as ethylene and ethanol, provides an attractive pathway towards closing the industrial carbon cycle and producing chemicals using renewable electricity. Existing CO 2 R technologies can exhibit C 2+ product selectivity above 74% and operating current densities over 1.3 A/cm 2 , but none can achieve both simultaneously, limiting the commercial viability and scalability of this technology. Low CO 2 concentrations at the catalyst surface and competing side reactions such as the hydrogen evolution reaction (HER) inhibit C 2+ faradaic efficiencies and operational current densities. Additionally, the use of high surface area cathode geometries to increase operational current densities is under-investigated. In this work, we electrodeposit a macroporous copper dendrite film to form high surface area cathodes. These structures are then coated with a combination of hydrophobic/hydrophilic ionomers to achieve high electrochemically active surface areas (ECSA), form regions of intrinsic porosity, and maximize catalytic availability of CO 2 . Using these techniques, we hope to maximize C 2+ faradaic efficiencies (FE) and geometric partial current densities. The treated cathode compound macrostructure (~10μm-sized pores, and ~μm-long dendrites) and concomitant surface area allows increased geometric current densities while the ionomer coating simultaneously inhibits HER and increases hydrogen availability to improve C 2+ selectivity. Our results show that careful engineering of the catalyst-electrolyte interface can enhance CO 2 reduction product selectivity and efficiency.

10 SYNTHETIC FUELS↗
DOE OSTI2020

Machine learning coupled multi-scale modeling for redox flow batteries

The reaction distribution in macro or device-scale has been studied for redox flow batteries. The reaction distribution on electrode pore-scale structure however is not well understood, lacking especially on how the reaction distribution on the pore-scale may impact the overall performance of a flow battery. This study introduces for the first time a framework of a multi-scale model that provides understanding of the relationship between the pore-scale electrode structure reaction and the device-scale electrochemical reaction uniformity within the flow battery. A reduced order model is constructed based on 128 pore-scale simulations, which provide a quantitative relationship between the battery operation conditions (inlet velocity, current density, inlet concentration) and the surface reaction uniformity for the pore-scale sample. The multi-scale framework upscales this pore-scale surface reaction uniformity to device-scale combined uniformity. Based on the multi-scale model, a time-varying optimization of the inlet velocity is established, leading to significant reduction on pump power consumption with targeted surface reaction uniformity. The multi-scale model establishes the critical link between the micro-structure of a flow battery component and its performance at the macro-scale, therefore providing rationale for further operational or material optimization.

flow batteries, machine learning, multi-scale mode↗
DOE OSTI2021

Quantum reactive scattering calculations for the cold and ultracold Li + LiNa → Li 2 + Na reaction

A first-principles based quantum dynamics study of the Li + LiNa(ν = 0, j = 0) → Li 2 (ν', j') + Na reaction is reported for collision energies spanning the ultracold (1 nK) to cold (1 K) regimes. A full-dimensional ab initio potential energy surface for the ground electronic state of Li 2 Na is utilized that includes an accurate treatment of the long-range interactions. The Li + LiNa reaction is barrierless and exoergic and exhibits a deep attractive potential well that supports complex formation. Thus, significant reactivity occurs even for collision temperatures approaching absolute zero. The reactive scattering calculations are based on a numerically exact time-independent quantum dynamics methodology in hyperspherical coordinates. Total and rotationally resolved rate coefficients are reported at 56 collision energies and include all contributing partial waves. Several shape resonances are observed in many of the rotationally resolved rate coefficients and a small resonance feature is also reported in the total rate coefficient near 50 mK. Furthermore, of particular interest, the angular distributions or differential cross sections are reported as a function of both the collision energy and scattering angle. Unique quantum fingerprints (bumps, channels, and ripples) are observed in the angular distributions for each product rotational state due to quantum interference and shape resonance contributions. The Li + LiNa reaction is under active experimental investigation so that these intriguing features could be verified experimentally when sufficient product state resolution becomes feasible for collision energies below 1 K.

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

Architecture Design for Remote Operation of Microreactors: Poster

The nuclear sector is pursuing a number of advanced reactor concepts, one of which is the microreactor, an advanced reactor characterized by a power output of less than 20 MWth. These reactors are intended for use in applications where traditional power solutions are economically or logistically impractical. One key feature required for the successful deployment of microreactors is a remote operation capability, which can dramatically cut staffing expenses by eliminating the need for licensed operators to be present at the site of each reactor and instead concentrate in a centrally located operations center. In moving to a remote operations framework for nuclear reactors, the number of potential attack surfaces for a cyber adversary looking to cause harm or disruption increases. Therefore, a robust cybersecurity architecture is required to mitigate these potential vulnerabilities. This paper explores the functional infrastructure, security, and communication requirements to adapt remote operations to nuclear applications. It then presents a reference architecture for remote operations of microreactors that applies best practices in cybersecurity and remote communications in the context of a digital twin remote operation system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Online Heatmap Generation with Both High and Low Weights

Heatmap is a common geovisualization method that interpolates and visualizes a set of point observations on a map surface. Most of online web mapping libraries implement a one-pass heatmap algorithm using HTML5 canvas or WebGL for efficient heatmap generation. However, such implementation applies additive operations that accumulate the rendering of point weights on the map surface grid, making it inappropriate for visualizations that require the highlighting of both low and high weights. We introduce \textit{hilomap}, an online heatmap algorithm that highlights surface areas where points with both low and high trends are located. An HTML5 Canvas-based reference implementation on OpenLayers is presented and evaluated.

Liu, Yan↗
DOE OSTI2022

Estimation of the Surface Fluxes for Heat and Momentum in Unstable Conditions with Machine Learning and Similarity Approaches for the LAFE Data Set

Abstract Measurements of three flux towers operated during the land atmosphere feedback experiment (LAFE) are used to investigate relationships between surface fluxes and variables of the land–atmosphere system. We study these relations by means of two machine learning (ML) techniques: multilayer perceptrons (MLP) and extreme gradient boosting (XGB). We compare their flux derivation performance with Monin–Obukhov similarity theory (MOST) and a similarity relationship using the bulk Richardson number (BRN). The ML approaches outperform MOST and BRN. Best agreement with the observations is achieved for the friction velocity. For the sensible heat flux and even more so for the latent heat flux, MOST and BRN deviate from the observations while MLP and XGB yield more accurate predictions. Using MOST and BRN for latent heat flux, the root mean square errors (RMSE) are 107 Wm $$^{-2}$$ - 2 and 121 Wm $$^{-2}$$ - 2 , respectively, as well as the intercepts of the regression lines are $$\approx 110$$ ≈ 110 Wm $$^{-2}$$ - 2 . For the ML methods, the RMSEs reduce to 31 Wm $$^{-2}$$ - 2 for MLP and 33 Wm $$^{-2}$$ - 2 for XGB as well as the intercepts to just 4 Wm $$^{-2}$$ - 2 for MLP and $$-1$$ - 1 Wm $$^{-2}$$ - 2 for XGB with slopes of the regression lines close to 1, respectively. These results indicate significant deficiencies of MOST and BRN, particularly for the derivation of the latent heat flux. In fact, in contrast to the established theories, feature importance weighting demonstrates that the ML methods base their improved derivations on net radiation, the incoming and outgoing shortwave radiations, the air temperature gradient, and the available water contents, but not on the water vapor gradient. The results imply that further studies of surface fluxes and other turbulent variables with ML techniques provide great promise for deriving advanced flux parameterizations and their implementation in land–atmosphere system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Tunable quantum-cascade VECSEL operating at 19 THz

We report a terahertz quantum-cascade vertical-external-cavity surface-emitting laser (QC-VECSEL) emitting around 1.9 THz with up to 10% continuous fractional frequency tuning of a single laser mode. The device shows lasing operation in pulsed mode up to 102 K in a high-quality beam, with the maximum output power of 37 mW and slope efficiency of 295 mW/A at 77 K. Challenges for up-scaling the operating wavelength in QC metasurface VECSELs are identified.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Low temperature hydrogen plasma permeation in palladium and its alloys for fuel recycling in fusion systems

Superpermeation of hydrogen isotopes through metal foils is a critical component for efficient fuel recycling in fusion power systems. In that context hydrogen permeation through foils of palladium and its alloys with silver and copper was studied at low temperature (60 - 200°C) under plasma exposure. These alloys differ significantly in both bulk and surface properties, and comparisons can provide mechanistic insights. Here, permeation was observed only during plasma operation, confirming the negligible contribution of molecular hydrogen to the observed flux. As-received foils required surface treatment to achieve top performance. For Pd and Pd75Ag25 an oxidation treatment increased permeation an order of magnitude, but proved unstable as this desirable surface was reduced under hydrogen plasma exposure. In contrast, an Ar plasma cleaning step provided both high and stable flux. As-received FCC phase Pd60Cu40 foils required annealing to transform it into the high permeability BCC phase that delivered top performance. All foils displayed similar temperature dependence with flux declining with temperature, suggesting that the primary rate-limiting step is absorption of superthermal hydrogen. Among these foils the hydrogen flux through PdCu was 3 - 5X greater than that of Pd or PdAg, which were similar. The superiority of BCC PdCu is attributed to its superior hydrogen desorption kinetics. Using PdCu 100% permeation of supplied hydrogen was achieved, and the flux saturated with increasing plasma power at values >10 –2 mol H•m –2 •s– 1 , and under these conditions permeation rates are equivalent with or without the membrane present. The fluxes achieved are the highest reported to date at these conditions, and the results highlight the important roles of both surface and bulk properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Evaluation of Hybrid Perovskite Prototypes After 10-Month Space Flight on the International Space Station

Metal halide perovskites (MHPs) have emerged as a prominent new photovoltaic material combining a very competitive power conversion efficiency that rivals crystalline silicon with the added benefits of tunable properties for multijunction devices fabricated from solution which can yield high specific power. Perovskites have also demonstrated some of the lowest temperature coefficients and highest defect tolerance, which make them excellent candidates for aerospace applications. However, MHPs must demonstrate durability in space which presents different challenges than terrestrial operating environments. To decisively test the viability of perovskites being used in space, a perovskite thin film is positioned in low earth orbit for 10 months on the International Space Station, which was the first long-duration study of an MHP in space. Postflight high-resolution ultrafast spectroscopic characterization and comparison with control samples reveal that the flight sample exhibits superior photo-stability, no irreversible radiation damage, and a suppressed structural phase transition temperature by nearly 65 K, broadening the photovoltaic operational range. Further, significant photo-annealing of surface defects is shown following prolonged light-soaking postflight. Here these results emphasize that methylammonium lead iodide can be packaged adequately for space missions, affirming that space stressors can be managed as theorized.

14 SOLAR ENERGY↗
DOE OSTI2024

Estimation of the Thermal-Fluids and Thermal-Structural Performance of Helium-Cooled Modular Finger-Type Divertors

Over the past decade, our group has investigated the thermal-fluid performance of the helium-cooled modular divertor with multiple jets (HEMJ) and a simplified “flat” design of the HEMJ for long-pulse magnetic fusion energy (MFE) reactors. Experimental studies were performed in a helium (He) loop at the prototypical pressure of 10 MPa, nearly prototypical He temperatures and incident heat fluxes using test sections made from stainless steel and tungsten alloys. Correlations for average Nusselt numbers and pressure loss coefficients were developed from the data and are used to validate computational fluid dynamics (CFD) models. Here, this work presents updated thermal-fluids performance curves based on these correlations that estimate the maximum heat flux that can be accommodated by the plasma-facing surface and coolant pumping power requirements at prototypical operating conditions. Thermal-structural performance curves developed from ITER structural design criteria are introduced, which include protection against ductile and non-ductile failure, ratcheting fatigue, and creep fatigue. The performance design curves for these finger-type divertors demonstrate that the “flat” design, with a significantly less complicated geometry than the HEMJ, has thermal-fluid and thermal-structural performance comparable to the original HEMJ concept, and may be superior in terms of non-ductile failure criteria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Evaluating the Origins of Aerobic Oxidation Catalysis with TAM-3, a MOF with Accessible Co(II) Sites and Large Pores

Metal-organic frameworks (MOFs) are attractive platforms that merge concepts of homogeneous and heterogeneous catalysis. Catalyst design and optimization are enabled by an array of synthetic methods that offer independent control over the local chemical structure of lattice-embedded metal ions (i.e., ligand identity and geometry) and the long-range materials properties (i.e., porosity). Establishing the origin of catalytic activity in MOF-promoted reactions remains a significant challenge: The relative rates of catalyst turnover and substrate diffusion dictate the extent to which interstitial sites are accessible and operational in catalysis. To minimize the contributions of surface sites in catalysis, materials with large pore dimensions are often sought, however, the impact of pore expansion on the origins of catalytic activity is similarly challenging to establish. Here, we describe TAM-3, a Co(II) based MOF with accessible metal sites supported by a facially coordinating tris-tetrazole ligand set. TAM-3 features large channel-like pores (17 × 23 Å) and promotes aerobic C−H oxidation and olefin epoxidation. Using a set of simple kinetics experiments, based on the analysis of kinetic isotope effects and olefin oxidation diastereoselectivities, we demonstrate that despite the large pores, interstitial metal ions do not significantly contribute to the observed substrate oxidation. This study highlights the importance of conducting kinetic experiments to assess the origin of apparent catalytic activity with MOFs and the challenge of harnessing reactive oxidants with microporous catalyst materials.

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

Catalytic Performance and Near-Surface X-ray Characterization of Titanium Hydride Electrodes for the Electrochemical Nitrate Reduction Reaction

The electrochemical nitrate reduction reaction (NO3RR) on titanium introduces significant surface reconstruction and forms titanium hydride (TiH x , 0 < x ≤ 2). With ex situ grazing-incidence X-ray diffraction (GIXRD) and X-ray absorption spectroscopy (XAS), we demonstrated near-surface TiH2 enrichment with increasing NO3RR applied potential and duration. This quantitative relationship facilitated electrochemical treatment of Ti to form TiH 2 /Ti electrodes for use in NO3RR, thereby decoupling hydride formation from NO 3 RR performance. A wide range of NO 3 RR activity and selectivity on TiH2/Ti electrodes between -0.4 and -1.0 VRHE was observed and analyzed with density functional theory (DFT) calculations on TiH 2 (111). Finally, this work underscores the importance of relating NO 3 RR performance with near-surface electrode structure to advance catalyst design and operation.

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

Peering into cloud physics using ultra-fine resolution radar and lidar systems

Cloud microphysical processes, such as droplet activation, condensational growth, and collisional growth, play a central role in the evolution of clouds and precipitation. Accurate representations of these processes in numerical models are challenging partially due to incomplete understanding of them at the process-level arising from limited systematic observations. Most surface-based active remote sensors, including today’s operational cloud radars and lidars, have a resolution on the order of tens of meters. This resolution is insufficient to resolve cloud microphysical processes that manifest at finer (meter and sub-meter) scales. A new set of ultra-high-resolution ground-based radar and lidar systems have been developed to address this observational gap. The newly developed 94-GHz cloud radar has a range resolution down to 2.8 m, or a factor of 10 finer than typical radars, using a large bandwidth and quadratic phase coding techniques. The lidar has a range resolution down to 10 cm, or a factor of 100 finer than typical lidars, using a time-gated time-correlated single photon counting technique. Such high-resolution observations were previously only achievable through in situ aircraft measurements. Even then, aircraft measurements do not permit continuous long-term cloud observation as is possible with ground-based remote sensing instruments. In this study, the first-light cloud observations from the new radar and lidar systems are shown to reveal detailed cloud structures that conventional sensors could only perceive in a bulk sense, thus providing new avenues to investigate cloud microphysical processes and their impact on weather and climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Coupling Subsurface and Above-Surface Models for Optimizing the Design of Borefields and District Heating and Cooling Systems

Accurate dynamic energy simulation is important for the design and sizing of district heating and cooling systems with geothermal heat exchange for seasonal energy storage. Current modeling approaches in building and district energy simulation tools typically consider heat conduction through the ground between boreholes without flowing groundwater. While detailed simulation tools for subsurface heat and mass transfer exist, these fall short in simulating above-surface energy systems. To support the design and operation of such systems, the study developed a coupled model including a software package for building and district energy simulation, and software for detailed heat and mass transfer in the subsurface. For the first, it uses the open-source Modelica Buildings Library, which includes dynamic simulation models for building and district energy and control systems. For the heat and mass transfer in the soil, it uses the TOUGH simulator. The TOUGH family of codes can model heat and multi-phase, multi-component mass transport for a variety of fluid systems, as well as chemical reactions, in fractured porous media. The study validated the coupled modeling approach by comparing the simulation results with one from the g-function based ground response model. It then looked into effects when the water table and the regional groundwater flow are considered in the ground, from the perspective of heat exchange between borehole and ground, and the electrical consumption of the district heating and cooling systems. To access the simulation models, please find the links in the submission: -- For coupled approach validation: see model Buildings.Fluid.Geothermal.Borefields.Examples.BorefieldsWithTough and Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from the "Modelica Building Library" resource, branch issue1495_tough_interface, commit a2667c0. -- For the study of the effect of water table: see model Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO, commit 760de49. -- For the study of the effect of regional groundwater flow: see Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO_3D, commit c2a2d2a. The coupling interface script "GrounResponse.py" can be found from the above links in the folder Buildings/Resources/Python-Sources. Also, the needed files for TOUGH simulation are in the folder Buildings/Resources/Python-Sources/ToughFiles that can be accessed through the above links. A brief description of these files is given below; detailed specifications for the first three files may be found in the TOUGH3 Users Guide (Jung et al., 2018) https://tough.lbl.gov/documentation/tough-manuals/. (1) INCON - initial conditions for each grid block (2) INFILE - main input file with material properties and control parameters (3) MESH - description of the computational grid (4) readsave - Modelica/TOUGH interface program: read the final output of TOUGH simulation after TOUGH time step and prepare for transfer to Modelica for next Modelica time step (5) readsave.inp - input parameters for program readsave (6) writeincon - Modelica/TOUGH interface program: write the output of Modelica after Modelica time step and prepare for transfer to TOUGH as initial conditions for the next TOUGH step (7) writeincon.inp - input parameters for program writeincon

15 GEOTHERMAL ENERGY↗