Search NASA⌕ Search

SEARCH · Search NASA

Results for “Working Fluid demonstrates”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Windowless Optical Cell Flow System

One of the goals for advanced integrated nuclear safeguards is on-line optical analysis of molten salt in molten salt reactors and electrochemical fuel reprocessing facilities. While progress has been made in the development of optical techniques that may be suitable for on-line analysis of molten salts, significant technology gaps exist with regard to integrating these analytical techniques into industrial-scale processes. One key technological hurdle is the short service life of optical window materials, which are subject to clouding and scaling during extended exposure to molten salts. To address this issue, Argonne is developing a gravity flow cell with open-apertures to support on-line optical analysis without the need for salt-wetted windows. The flow cell uses gravity flow and fluid dynamics principles to enable molten salt flow through a pipe that has holes in its wall through which salt does not flow. In FY21, an open-aperture gravity flow cell with windowless optical access to flowing salt was successfully demonstrated. Future work should include the integration of the flow cell into a sampling loop, stakeholder outreach and collaboration, integration of multiple optical analysis techniques, evaluation of the integrated methods, and analyses to determine how the various tools might fit into an integrated safeguards monitoring system of unattended near real time monitoring tools.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1990

Pool boiler heat transport system for a 25 kWe advanced Stirling conversion system

Experiments to determine alkali metal/enhanced surface combinations that have stable boiling at the temperatures and heat fluxes that occur in the Stirling engine are reported. Two enhanced surfaces and two alkali metal working fluids were evaluated. The enhanced surfaces were an EDM hole covered surface and a sintered-powder-metal porous layer surface. The working fluids tested were potassium and eutectic sodium-potasium alloy (NaK), both with and without undissolved noncondensible gas. Noncondensible gas (He and Xe) was added to the system to provide gas in the nucleation sites, preventing quenching of the sites. The experiments demonstrated the potential of an alkali metal pool boiler heat transport system for use in a solar-powered Stirling engine. The most favorable fluid/surface combination tested was NaK boiling on a -100 +140 mesh 304L stainless steel sintered porous layer with no undissolved noncondensible gas. This combination provided stable, high-performance boiling at the operating temperature of 700 C. Heat fluxes into the system ranged from 10 to 50 W/sq cm. The transition from free convection to nucleate boiling occurred at temperatures near 540 C. Based on these experiments, a pool boiler was designed for a full-scale 25-kWe Stirling system.

Anderson, W. G.↗
NASA NTRS2008

Development of the Aquarius Antenna Deployment Mechanisms and Spring/Damper Actuator

The Aquarius Instrument s large radar reflector dish needed to be stowed for launch, and then deployed on-orbit. The Deployment Subsystem consisted of a cantilevered boom structure and two single-axis hinge mechanisms to accurately deploy and position the reflector dish relative to the radar feed horns. The cantilevered design demanded high stiffness and accuracy from the deployment mechanism at the root of the boom. A preload-generating end-of-travel latch was also required. To largely eliminate the need for control systems, each deployment mechanism was actuated by a passive spring motor with viscous-fluid damping. Tough requirements and adaptation of a heritage actuator to the new application resulted in numerous challenges. Fabrication, assembly, and testing encountered additional problems, though ultimately the system was demonstrated very successfully. This paper revisits the development to highlight which design concepts worked and the many important lessons learned.

Johnson, Joel A.↗
NASA NTRS2003

Microscale Thermal-Transpiration Gas Pump

A recent addition to the growing class of microelectromechanical systems (MEMS) is a single stage of a Knudsen compressor. This device was fabricated and tested to demonstrate the feasibility of Knudsen compressors as miniature vacuum pumps for future portable scientific instruments. The attributes of Knudsen compressors that make them attractive as miniature vacuum pumps are that they contain no moving parts and operate without need for lubricants or working fluids.

Vargo, Stephen↗
NASA NTRS2015

Pressure Loss Predictions of the Reactor Simulator Subsystem at NASA GRC

Testing of the Fission Power System (FPS) Technology Demonstration Unit (TDU) is being conducted at NASA GRC. The TDU consists of three subsystems: the Reactor Simulator (RxSim), the Stirling Power Conversion Unit (PCU), and the Heat Exchanger Manifold (HXM). An Annular Linear Induction Pump (ALIP) is used to drive the working fluid. A preliminary version of the TDU system (which excludes the PCU for now), is referred to as the RxSim subsystem and was used to conduct flow tests in Vacuum Facility 6 (VF 6). In parallel, a computational model of the RxSim subsystem was created based on the CAD model and was used to predict loop pressure losses over a range of mass flows. This was done to assess the ability of the pump to meet the design intent mass flow demand. Measured data indicates that the pump can produce 2.333 kg/sec of flow, which is enough to supply the RxSim subsystem with a nominal flow of 1.75 kg/sec. Computational predictions indicated that the pump could provide 2.157 kg/sec (using the Spalart-Allmaras turbulence model), and 2.223 kg/sec (using the k- turbulence model). The computational error of the predictions for the available mass flow is -0.176 kg/sec (with the S-A turbulence model) and -0.110 kg/sec (with the k-epsilon turbulence model) when compared to measured data.

Turbulence Models↗
DOE OSTI2025

High fidelity multiphysics tightly coupled model for a lead cooled fast reactor concept and application to statistical calculation of hot channel factors

A tightly coupled multiphysics code system is established using the MOOSE framework for hot channel factor (HCF) evaluation on a Lead Fast Reactor (LFR) concept. The coupled system is driven by the Griffin multiphysics coupling capability under which the MOOSE Heat Transfer module and NekRS computational fluid dynamics solver are coupled for conjugate heat transfer using the Cardinal application. The coupled capability is demonstrated on an LFR assembly model based on materials and geometry of a prototypical lead-cooled fast reactor design by Westinghouse Electric Company, LLC. Moreover, the work integrates the Multiphysics Object Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM) to perform calculations for statistical analysis of HCF. Furthermore, the coupling strategy and workflow demonstrated in this paper is not only useful for predicting accurate hot channel factors for different kinds of advanced reactors but also for other engineering applications such as control rod worth assessment, generation of high-fidelity database for Artificial intelligence (AI)/machine learning (ML) training, design optimization and multi-resolution modeling.

Cardinal↗
NASA NTRS1986

Numerical modeling of physical vapor transport in a vertical cylindrical ampoule, with and without gravity

Numerical modeling has been performed of the fluid dynamics in a prototypical physical vapor transport crystal growing situation. Cases with and without gravity have been computed. Dependence of the flows upon the dimensionless parameters aspect ratio and Peclet, Rayleigh, and Schmidt numbers is demonstrated to a greater extent than in previous works. Most notably, it is shown that the effects of thermally-induced buoyant convection upon the mass flux on the growth interface crucially depend upon the temperature boundary conditions on the sidewall (e.g., whether adiabatic or of a fixed profile, and in the latter case the results depend upon the shape of the profile assumed).

Miller, T. L.↗
NASA NTRS2002

Magnetohydrodynamic Augmented Propulsion Experiment

A fundamental obstacle to routine space access is the specific energy limitations associated with chemical fuels. In the case of vertical take-off, the high thrust needed for vertical liftoff and acceleration to orbit translates into power levels in the 10 GW range. Furthermore, useful payload mass fractions are possible only if the exhaust particle energy (i.e., exhaust velocity) is much greater than that available with traditional chemical propulsion. The electronic binding energy released by the best chemical reactions (e.g., LOX/LH2 for example, is less than 2 eV per product molecule (approx. 1.8 eV per H2O molecule), which translates into particle velocities less than 5 km/s. Useful payload fractions, however, will require exhaust velocities exceeding 15 km/s (i.e., particle energies greater than 20 eV). As an added challenge, the envisioned hypothetical RLV (reusable launch vehicle) should accomplish these amazing performance feats while providing relatively low acceleration levels to orbit (2-3g maximum). From such fundamental considerations, it is painfully obvious that planned and current RLV solutions based on chemical fuels alone represent only a temporary solution and can only result in minor gains, at best. What is truly needed is a revolutionary approach that will dramatically reduce the amount of fuel and size of the launch vehicle. This implies the need for new compact high-power energy sources as well as advanced accelerator technologies for increasing engine exhaust velocity. Electromagnetic acceleration techniques are of immense interest since they can be used to circumvent the thermal limits associated with conventional propulsion systems. This paper describes the Magnetohydrodynamic Augmented Propulsion Experiment (MAPX) being undertaken at NASA Marshall Space Flight Center (MSFC). In this experiment, a 1-MW arc heater is being used as a feeder for a 1-MW magnetohydrodynamic (MHD) accelerator. The purpose of the experiment is to demonstrate that an MHD accelerator can be an effective augmentation system for increasing engine exhaust velocity. More specifically, the experiment is intended to show that electromagnetic effects are effective at producing flow acceleration whereas electrothermal effects do not cause unacceptable heating of the working fluid. The MHD accelerator was designed as an externally diagonalized segmented Faraday channel, which will be inserted into an existing 2-tesla electromagnet. This allows the external power to be connected through two terminals thereby minimizing the complexity and cost associated with powering each segment independently. The design of the accelerator and other components in the flow path has been completed and fabrication activities are underway. This paper provides a full description of MAPX including performance analysis, design, and test plans, and current status.

Litchford, Ron J.↗
NASA NTRS1988

Overview of heat transfer and fluid flow problem areas encountered in Stirling engine modeling

NASA Lewis Research Center has been managing Stirling engine development programs for over a decade. In addition to contractual programs, this work has included in-house engine testing and development of engine computer models. Attempts to validate Stirling engine computer models with test data have demonstrated that engine thermodynamic losses need better characterization. Various Stirling engine thermodynamic losses and efforts that are underway to characterize these losses are discussed.

Tew, Roy C., Jr.↗
DOE OSTI2021

Renormalized site density functional theory

Site density functional theory (SDFT) provides a rigorous framework for statistical mechanics analysis of inhomogeneous molecular liquids. The key defining feature of these systems is the presence of two very distinct interactions scales (intra- and inter-molecular), and as such proper description of both effects is critical to the accuracy of the calculations. Current SDFT applications utilize the same approximation scheme for both interaction motifs, which negatively impacts the results. Dual space methodology, used in this work, alleviates this issue by providing the flexibility of evaluating part of the interactions in traditional field representation. For molecular liquid this translates into retaining density representation for inter-molecular interactions but describing stiff intra-molecular remainder by more appropriate conventional field based methods. This opens the way to decouple analysis of the two interactions scales - the idea which is developed further in this work for the case of homogeneous reference approximation of inter-molecular interactions. We demonstrate that by defining new collective variables, the behavior of the original molecular liquid system at the inter-molecular level can be transformed to resemble that of an effective simple fluid mixture. The latter is linked to intra-molecular scale through renormalized interaction parameters. We illustrate this renormalization procedure for several types of diatomic liquids, showing that this approach cures many of the shortcomings of existing SDFT methods.

Chuev, Gennady N.↗
NASA NTRS1993

The Influence of Swirl Brakes and a Tip Discharge Orifice on the Rotordynamic Forces Generated by Discharge-to-Suction Leakage Flows in Shrouded Centrifugal Pumps

Recent experiments conducted in the Rotor Force Test Facility at the California Institute of Technology have examined the effects of a tip leakage restriction and swirl brakes on the rotordynamic forces due to leakage flows on an impeller undergoing a prescribed circular whirl. The experiments simulate the leakage flow conditions and geometry of the Alternate Turbopump Design (ATD) of the Space Shuttle High Pressure Oxygen Turbopump and are critical to evaluating the pump's rotordynamic instability problems. Previous experimental and analytical results have shown that discharge-to-suction leakage flows in the annulus of a shrouded centrifugal pump contribute substantially to the fluid induced rotordynamic forces. Also, previous experiments have shown that leakage inlet (pump discharge) swirl can increase the cross-coupled stiffness coefficient and hence increase the range of positive whirl for which the tangential force is destabilizing. In recent experimental work, the present authors demonstrated that when the swirl velocity within the leakage path is reduced by the introduction of ribs or swirl brakes, then a substantial decrease in both the destabilizing normal and tangential forces could be achieved. Motivation for the present research is that previous experiments have shown that restrictions such as wear rings or orifices at pump inlets affect the leakage forces. Recent pump designs such as the Space Shuttle Alternate Turbopump Design (ATD) utilize tip orifices at discharge for the purpose of establishing axial thrust balance. The ATD has experienced rotordynamic instability problems and one may surmise that these tip discharge orifices may also have an important effect on the normal and tangential forces in the plane of impeller rotation. The present study determines if such tip leakage restrictions contribute to undesirable rotordynamic forces. Additional motivation for the present study is that the widening of the leakage path annular clearance and the installation of swirl brakes in the ATD has been proposed to solve its instability problems. The present study assesses the effect of such a design modification on the rotordynamic forces. The experimental apparatus consists of a solid or dummy impeller, a housing instrumented for pressure measurements, a rotating dynamometer and an eccentric whirl mechanism. The solid impeller is used so that leakage flow contributions to the forces are measured, but the main throughflow contributions are not experienced. The inner surface of the housing has been modified to accommodate meridional ribs or swirl brakes within the leakage annulus. In addition, the housing has been modified to accommodate a discharge orifice that qualitatively simulates one side of the balance piston orifice of the Space Shuttle ATD. Results indicate the detrimental effects of a discharge orifice and the beneficial effects of brakes. Plots of the tangential and normal forces versus whirl ratio show a substantial increase in these forces along with destabilizing resonances at some positive whirl ratios when a discharge orifice is added. When brakes are added, some of the detrimental effects of the orifice are reduced. For the tangential force, a plot versus whirl ratio shows a significant reduction and a destabilizing resonance appears to be eliminated. For the normal force, although the overall force is not reduced, again a destabilizing resonance appears to be eliminated.

Sivo, Joseph M.↗
NASA NTRS2011

Post-Test Analysis of a 10-Year Sodium Heat Pipe Life Test

High-temperature heat pipes are being evaluated for use in energy conversion applications such as fuel cells, gas turbine re-combustors, Stirling cycle heat sources; and with the resurgence of space nuclear power both as reactor heat removal elements and as radiator elements. Long operating life and reliable performance are critical requirements for these applications. Accordingly, long-term materials compatibility is being evaluated through the use of high-temperature life test heat pipes. Thermacore, Inc., has carried out a sodium heat pipe 10-year life test to establish long-term operating reliability. Sodium heat pipes have demonstrated favorable materials compatibility and heat transport characteristics at high operating temperatures in air over long time periods. A representative one-tenth segment Stirling Space Power Converter heat pipe with an Inconel 718 envelope and a stainless steel screen wick has operated for over 87,000 hr (10 years) at nearly 700 C. These life test results have demonstrated the potential for high-temperature heat pipes to serve as reliable energy conversion system components for power applications that require long operating lifetime with high reliability. Detailed design specifications, operating history, and post-test analysis of the heat pipe and sodium working fluid are described. Lessons learned and future life test plans are also discussed.

Rosenfeld, John H.↗
NASA NTRS2012

Ten Year Operating Test Results and Post-Test Analysis of a 1/10 Segment Stirling Sodium Heat Pipe, Phase III

High-temperature heat pipes are being evaluated for use in energy conversion applications such as fuel cells, gas turbine re-combustors, Stirling cycle heat sources; and with the resurgence of space nuclear power both as reactor heat removal elements and as radiator elements. Long operating life and reliable performance are critical requirements for these applications. Accordingly, long-term materials compatibility is being evaluated through the use of high-temperature life test heat pipes. Thermacore, Inc., has carried out a sodium heat pipe 10-year life test to establish long-term operating reliability. Sodium heat pipes have demonstrated favorable materials compatibility and heat transport characteristics at high operating temperatures in air over long time periods. A representative one-tenth segment Stirling Space Power Converter heat pipe with an Inconel 718 envelope and a stainless steel screen wick has operated for over 87,000 hr (10 yr) at nearly 700 C. These life test results have demonstrated the potential for high-temperature heat pipes to serve as reliable energy conversion system components for power applications that require long operating lifetime with high reliability. Detailed design specifications, operating history, and post-test analysis of the heat pipe and sodium working fluid are described.

Rosenfeld, John, H↗
NASA NTRS2016

Pressure Loss Predictions of the Reactor Simulator Subsystem at NASA Glenn Research Center

Testing of the Fission Power System (FPS) Technology Demonstration Unit (TDU) is being conducted at NASA Glenn Research Center. The TDU consists of three subsystems: the reactor simulator (RxSim), the Stirling Power Conversion Unit (PCU), and the heat exchanger manifold (HXM). An annular linear induction pump (ALIP) is used to drive the working fluid. A preliminary version of the TDU system (which excludes the PCU for now) is referred to as the "RxSim subsystem" and was used to conduct flow tests in Vacuum Facility 6 (VF 6). In parallel, a computational model of the RxSim subsystem was created based on the computer-aided-design (CAD) model and was used to predict loop pressure losses over a range of mass flows. This was done to assess the ability of the pump to meet the design intent mass flow demand. Measured data indicates that the pump can produce 2.333 kg/sec of flow, which is enough to supply the RxSim subsystem with a nominal flow of 1.75 kg/sec. Computational predictions indicated that the pump could provide 2.157 kg/sec (using the Spalart-Allmaras (S‒A) turbulence model) and 2.223 kg/sec (using the k- turbulence model). The computational error of the predictions for the available mass flow is ‒0.176 kg/sec (with the S-A turbulence model) and -0.110 kg/sec (with the k- turbulence model) when compared to measured data.

Mass Flow↗
NASA NTRS2026

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗
DOE OSTI2021

Progress on Optimizing Wind Farms and Rotor Designs Using Adjoints

Modern wind plants are increasingly tasked with multiple performance objectives. In addition to designing plants that maximize power output and minimize the levelized cost of energy (LCOE), the design and operation of wind plants is increasingly influenced by challenges regarding grid integration of variable generation renewables. This places a growing emphasis on making wind plants more controllable and predictable. WindSE is a Reynolds-averaged Navier-Stokes (RANS) model designed around analytical gradient and adjoint methods, with the ability to capture terrain-induced effects, as shown in Figure 1. The recent addition of an unsteady solver with an actuator line method (ALM) and ongoing work to enable massively parallel optimizations gives it a unique niche to explore coupled plant-level controls and design problems. This code is an open source python package built on the FEniCS framework that utilizes fast, parallel PETSc solvers to model fluid flow throughout wind-farm scale domains. Two recent studies performed using WindSE demonstrate the capability to optimize under a wide variety of flow conditions and objective functions. In the first, we present an optimization focused on modifying the layout of a wind farm with a fixed number of turbines for maximum total power output [1]. This study highlights the ability to quickly perform simulations using the steady Navier-Stokes solver combined with rotors represented as actuator disks while also stressing the importance of capturing terrain-induced effects. Gradient-based optimization using the RANS equations is viable due to the inclusion of efficiently computed adjoint derivatives. We interpret the physical results of the optimal layout and also discuss the computational cost of scaling to larger problems. In the second study, we present the capabilities of the unsteady Navier-Stokes solver, where rotor-blade profiles represented by actuator lines are optimized to enhance wake steering effects and overall power production [2]. We quantify the wind plant performance gains obtained from this type of simultaneous control co-design optimization as compared to optimizing the blade design and yaw independently. Figure 2 shows the differences between a baseline two-turbine system and an optimized system where we fine-tune the blade chord profile. Results and challenges from each study are quickly summarized and used to motivate the current development efforts within WindSE. Current and future work is focused on enabling higher-resolution studies with more degrees of freedom through parallelization of both the simulation and optimization algorithms. We present benchmarking results to show that WindSE performs well in both weak- and strong-scaling tests and further demonstrate that the optimizer obtains the same convergence rates in both shared- and distributed-memory environments. Using larger wind farms, we can study deep-array effects within an optimization context, allowing the use of objective functions that have been previously unstudied. As an example, we present ongoing work on a blockage metric which characterizes the loss of available kinetic energy due to wake effects from multiple upstream turbines.

adjoint optimization↗
DOE OSTI2021

High-temperature high-pressure microfluidic system for rapid screening of supercritical CO 2 foaming agents

CO 2 foam helps to increase the viscosity of CO 2 flood fluid and thus improve the process efficiency of the anthropogenic greenhouse gas’s subsurface utilization and sequestration. Successful CO 2 foam formation mandates the development of high-performance chemicals at close to reservoir conditions, which in turn requires extensive laboratory tests and evaluations. This work demonstrates the utilization of a microfluidic reservoir analogue for rapid evaluation and screening of commercial surfactants (i.e., Cocamidopropyl Hydroxysultaine, Lauramidopropyl Betaine, Tallow Amine Ethoxylate, N,N,N' Trimethyl-N'-Tallow-1,3-diaminopropane, and Sodium Alpha Olefin Sulfonate) based on their performance to produce supercritical CO 2 foam at high salinity, temperature, and pressure conditions. The microfluidic analogue was designed to represent the pore sizes of the geologic reservoir rock and to operate at 100 °C and 13.8 MPa. Values of the pressure drop across the microfluidic analogue during flow of the CO 2 foam through its pore network was used to evaluate the strength of the generated foam and utilized only milliliters of liquid. The transparent microfluidic pore network allows in-situ quantitative visualization of CO 2 foam to calculate its half-life under static conditions while observing if there is any damage to the pore network due to precipitation and blockage. The microfluidic mobility reduction results agree with those of foam loop rheometer measurements, however, the microfluidic approach provided more accurate foam stability data to differentiate the foaming agent as compared with conventional balk testing. The results obtained here supports the utility of microfluidic systems for rapid screening of chemicals for carbon sequestration or enhanced oil recovery operations.

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

Machine learning of consistent thermodynamic models using automatic differentiation

In this study, we propose a data-driven method to describe consistent equations of state (EOS) for arbitrary systems. Complex EOS are traditionally obtained by fitting suitable analytical expressions to thermophysical data. A key aspect of EOS is that the relationships between state variables are given by derivatives of the system free energy. In this work, we model the free energy with an artificial neural network and utilize automatic differentiation to directly learn the derivatives of the free energy. We demonstrate this approach on two different systems, the analytic van der Waals EOS and published data for the Lennard-Jones fluid, and we show that it is advantageous over direct learning of thermodynamic properties (i.e., not as derivatives of the free energy but as independent properties), in terms of both accuracy and the exact preservation of the Maxwell relations. Furthermore, the method implicitly provides the free energy of a system without explicit integration.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗