Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Thermal fluid physics”

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 361 records · Page 20

Internally-cooled compressor diaphragm

An internally-cooled diaphragm for an internally-cooled compressor is provided. The internally-cooled diaphragm may include an annular body configured to cool a process fluid flowing through a fluid pathway of the internally-cooled compressor. The annular body may define a return channel of the fluid pathway, and a cooling pathway in thermal communication with the fluid pathway. The return channel may be configured to at least partially diffuse and de-swirl the process fluid flowing therethrough, and the cooling pathway may be configured to receive a coolant to absorb heat from the process fluid flowing through the return channel.

42 ENGINEERING↗

Development of an Advanced Hydrogen Energy Storage System using Aerogel in a Cryogenic Flux Capacitor (CFC)

The Cryogenic Flux Capacitor (CFC) is a cold, dense fluid storage core with integrated design features that afford the designer flexibility and provide new possibilities for the storage and discharge of energy. The stored energy, in this case, is represented by hydrogen physically bonded within the nanoscale pores within the aerogel composite blanket material, and the process of bonding or debonding is governed by principles of physical adsorption (physisorption) and thermodynamics. The large surface area afforded by the nanoporous aerogel (~1,000 m2/g) allows for storage densities close to, or in some cases exceeding, that of normal boiling point liquids. Its performance easily exceeds what can be achieved via ambient temperature, high-pressure gas storage for an equivalent volume. CFC storage is predicted to be easily scalable, is constructed from readily available commercial materials, lends itself to a range of pressure applications, and is geometry insensitive. They can also be used in modular designs, affording even more flexibility for potential deployment. In addition to the aerogel adsorbent within the blanket material, a CFC includes thermally-conductive membrane layered with the aerogel, which acts as a large-area, quick-response thermal management system. The system conducts heat throughout the volume to discharge the unit quickly. This same thermal management system can also be connected to a refrigeration system, or cold fluid such as liquid nitrogen (LN2), to facilitate the charging up of the CFC. The charging and discharging can be performed as fast or as slow as desired by controlling the cooling or heating supplied to the CFC.

08 HYDROGEN↗

A numerical method for predicting hypersonic flowfields

The flow about a body traveling at hypersonic speed is energetic enough to cause the atmospheric gases to chemically react and reach states in thermal nonequilibrium. The prediction of hypersonic flowfields requires a numerical method capable of solving the conservation equations of fluid flow, the chemical rate equations for specie formation and dissociation, and the transfer of energy relations between translational and vibrational temperature states. Because the number of equations to be solved is large, the numerical method should also be as efficient as possible. The proposed paper presents a fully implicit method that fully couples the solution of the fluid flow equations with the gas physics and chemistry relations. The method flux splits the inviscid flow terms, central differences of the viscous terms, preserves element conservation in the strong chemistry source terms, and solves the resulting block matrix equation by Gauss Seidel line relaxation.

Maccormack, Robert W.↗

Simulation and flow physics of a shocked and reshocked high-energy-density mixing layer

This work describes a computational investigation of multimode instability growth and multimaterial mixing induced by multiple shock waves in a high-energy-density (HED) environment, where pressures exceed 1 Mbar. The simulations are based on a series of experiments performed at the National Ignition Facility (NIF) and designed as an HED analogue of non-HED shock-tube studies of the Richtmyer–Meshkov instability and turbulent mixing. A three-dimensional computational modelling framework is presented. It treats many complications absent from canonical non-HED shock-tube flows, including distinct ion and free-electron internal energies, non-ideal equations of state, radiation transport and plasma-state mass diffusivities, viscosities and thermal conductivities. The simulations are tuned to the available NIF data, and traditional statistical quantities of turbulence are analysed. Integrated measures of turbulent kinetic energy and enstrophy both increase by over an order of magnitude due to reshock. Large contributions to enstrophy production during reshock are seen from both the baroclinic source and enstrophy–dilatation terms, highlighting the significance of fluid compressibility in the HED regime. Dimensional analysis reveals that Reynolds numbers and diffusive Péclet numbers in the HED flow are similar to those in a canonical non-HED analogue, but conductive Péclet numbers are much smaller in the HED flow due to efficient thermal conduction by free electrons. It is shown that the mechanism of electron thermal conduction significantly softens local spanwise gradients of both temperature and density, which causes a minor but non-negligible decrease in enstrophy production and small-scale mixing relative to a flow without this mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coupled Monte Carlo Transport and Conjugate Heat Transfer for Wire-Wrapped Bundles Within the MOOSE Framework

Cardinal is an open-source application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics (CFD) to the Multiphysics Object-Oriented Simulation Environment (MOOSE), closing neutronics and thermal-fluid gaps in conducting high-resolution multiscale and multiphysics analyses of nuclear systems. We first provide a brief introduction to Cardinal's software design, data mapping, and coupling strategy to highlight our approach to overcoming common challenges in high-fidelity multiphysics simulations. Here we then present two Cardinal simulations for hexagonal pin bundles. The first is a validation of Cardinal's conjugate heat transfer coupling of NekRS's Reynolds-Averaged Navier Stokes model with MOOSE's heat conduction physics for a bare seven-pin Freon-12 bundle flow experiment. Predictions for pin surface temperatures under three different heating modes agree reasonably well with experimental data and similar CFD modeling from the literature. The second simulation is a multiphysics coupling of OpenMC, NekRS, and BISON for a reduced-scale, seven-pin wire-wrapped version of an Advanced Burner Reactor bundle. Wire wraps are approximated using a momentum source model, and coupled predictions are provided for velocity, temperature, and power distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Thermal Hydraulics Analysis of a Divertor Monoblock Using SALAMANDER

Divertors are critical components in magnetic confinement fusion devices. One of the Divertor's crictial roles is absorbing the highest heat flux from the plasma. The Divertor accomplishes this through the use of divertor monoblocks and cooling channels. This work presents a thermal hydraulic analysis of a divertor monoblock using the Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering and Research (SALAMANDER). SALAMANDER is an open source tool developed at Idaho National Laboratory for conducting multiphysics and multscale analysis of Plasma Facing Components (PFCs). For the monoblock, solid heat conduction is modeled in the block with a constant heat flux of 1E7 W/m^2 applied to the top of the block to represent peak heat flux from the plasma. The cooling channel (which is pressurized water at a Reynolds number of ~1,000,000) is modeled using a k-epsilon turbulence model with standard wall functions. The heat transfer between the solid and fluid domain is modeled using the Dittus-Boelter correlation for the convective heat transfer coefficient. This work aims to show the importance of a multiphysics modeling approach when performing simulations on PFCs.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Center of Excellence Collaboration Projects: Second Wave

The NEAMS program aims to develop an integrated multi-physics simulation capability “pellet-to-plant” for the design and analysis of future generations of nuclear power plants. In particular, the Reactor Product Line code suite's multi-resolution hierarchy is being designed to ultimately span the full range of length and time scales present in relevant reactor design and safety analyses, as well as scale from desktop to petaflop computing platforms. In particular the NEAMS program is supporting the development of novel thermal-hydraulic codes. The Center of Excellence for Thermal Fluids Application in Nuclear Energy, launched in 2018, has as its key goals to serve as a front door to industry. The Center of Excellence for Thermal Fluids Applications in Nuclear Energy has recently launched a program to start collaborative efforts between the laboratories and industry with the objective of stimulating cooperation and increasing adoption of T/H tools developed under NEAMS by the industry-at-large. In particular, two industry partners agreed to participate in a second wave of short-term collaborations aimed at demonstrating the value of NEAMS tools to their designs. These are Framatome and TerraPower who proposed projects related to LWR (Light Water Reactor) and MSFR (Molten Salt Fast Reactor) designs respectively. In this report we present the results of these two collaborations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development for Integrated System-Level Analysis Capabilities in SAM for Molten Salt Reactors

In recent years, there has been renewed interest in Molten Salt Reactors (MSRs) for their potential advantages compared to reactors that rely on solid fuel. In response to such interest, many methods and codes have been developed to capture the unique features of MSRs. Among them, the System Analysis Module (SAM) is a modern system analysis tool that provides fast-running, modest-fidelity, whole-plant transient analysis capabilities, essential for fast-turnaround design scoping and engineering analyses of advanced reactor concepts. For liquid-fuel MSRs, the complex physics and chemistry involved in MSR operation—such as reactor kinetics, fluid flow, heat transfer, and salt composition dynamics—pose significant challenges for system-level modeling. Specific modeling capabilities are needed for system-level transient simulation. This paper presents recent advancements in SAM capability enhancements for system-level modeling of MSRs, focusing on improved simulation fidelity, computational efficiency, and multi-physics integration. Key enhancements include the development of species transport, Delayed Neutron Precursor (DNP) drift, modified Point Kinetics Equations (PKE), decay heat modeling, key fission product behavior, salt corrosion, and thermal-hydraulic coupling, as well as code robustness and performance enhancements for MSR applications. The code enhancement allows for better predictive accuracy in safety analysis, transient behavior, and operational optimization, thus supporting the design and licensing of next-generation MSRs. Results from case studies are presented to demonstrate the benefits of these enhancements in accurately capturing key reactor transient behaviors.

Hu, Rui (ORCID:0000000237712920)↗

Validation study of RWM stability in DIII-D high- β N plasmas

The n = 1 (n is the toroidal mode number) resistive wall mode (RWM) stability is numerically investigated for two DIII-D high-β N discharges 176440 and 172461, utilizing the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu et al 2008 Phys. Plasmas 15 112503) codes. Systematic validation efforts are attempted, for the first time, for discharges with very slow or vanishing toroidal flow for a large fraction of the plasma volume. While gaining physics insights in accessing stable operation regime at β N exceeding the Troyon no-wall limit in these slow-rotation experiments, the predictive capability of fluid and non-perturbative magnetohydrodynamic-kinetic hybrid models for the RWM is further confirmed. The MARS-F fluid model, with a strong but numerically tunable viscosity mimicking ion Landau damping of parallel sound waves, finds complete stabilization of the n = 1 RWM in the considered DIII-D plasmas under the experimental flow conditions. Similarly, either full stabilization (for discharge 176440) or marginal stability (for discharge 172461) of the mode is computed by the MARS-K hybrid model, which is first-principle based without free model parameters. In particular, all drift kinetic resonances, including those of thermal and energetic particles, are found to synergistically act to marginally stabilize the RWM in discharge 172461. These MARS-F/K modeling results explain the experimentally observed stable operational regime in DIII-D, as far as the RWM stability is concerned. Extensive numerical sensitivity studies, with respect to the plasma toroidal flow speed as well as the radial location of the resistive wall, are also carried out to further support the validation study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neutron transport methods for multiphysics heterogeneous reactor core simulation in Griffin

Griffin is a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). This work discloses the methods, algorithms, and implementation for simulating heterogeneous reactor dynamics models. Griffin utilizes a discontinuous finite-element method with discrete ordinates (DFEM-S ) to discretize the field variable of the multigroup neutron transport equation. Multiphysics feedback is handled using two-step tabulated cross-section methodology. Feedback quantities are evaluated using the MOOSE-MultiApp system to couple various engineering phenomena, such as heat conduction and thermal fluids. The multiphysics DFEM-S system is solved using fixed-point iteration with a fully asynchronous parallel sweeper, unstructured coarse-mesh finite difference acceleration, and a multi-timescale improved quasi-static method scheme. The implementation is applied to a multiphysics microreactor model, with two transients: one initiated by a single heat-pipe failure and another by control drum rotation. Importantly, these examples demonstrate the ability of Griffin to tractably solve the neutron transport equation considering seven independent variables and feedback.

97 MATHEMATICS AND COMPUTING↗

A theory of local and global processes which affect solar wind electrons. 1: The origin of typical 1 AU velocity distribution functions: Steady state theory

A detailed first principle kinetic theory for electrons which is neither a classical fluid treatment nor an exospheric calculation is presented. This theory illustrates the global and local properties of the solar wind expansion that shape the observed features of the electron distribution function, such as its bifurcation, its skewness and the differential temperatures of the thermal and suprathermal subpopulations. Coulomb collisions are substantial mediators of the interplanetary electron velocity distribution function and they place a zone for a bifurcation of the electron distribution function deep in the corona. The local cause and effect precept which permeates the physics of denser media is modified for electrons in the solar wind. The local form of transport laws and equations of state which apply to collision dominated plasmas are replaced with global relations that explicitly depend on the relative position of the observer to the boundaries of the system.

Scudder, J. D.↗

Data Base for CFD Validation

The flow behind the shock wave formed around objects which fly at hypervelocity behaves differently from that of a perfect gas. Molecules become vibrationally excited, dissociated, and ionized. The hot gas may emit or absorb radiation. When the atoms produced by dissociation reach the wall surface, chemical reactions, including recombination, may occur. The thermochemical phenomena of vibration, dissociation, ionization, surface chemical reaction, and radiation are referred to commonly as high-temperature real-gas phenomena. The phenomena cause changes in the dynamic behavior of the flow and the surface pressure and heat transfer distribution around the object. The character of a real gas is described by the internal degrees of freedom and state of constituent molecules; nitrogen and oxygen for air. The internal energy states, rotation, vibration and electronic, of the molecules are excited and, in the limit, the molecular bonds are exceeded and the gas dissociated into atomic and, possibly, ionic constituents. The process of energy transfer causing excitation, dissociation and recombination is a rate process controlled by particle collisions. Binary, two-body, collisions are sufficient to cause internal excitation, dissociation and ionization while three-body collisions are required to recombine the particles into molecular constituents. If the rates of energy transfer are fast with respect to the local fluid dynamic time scale the gas is in, or nearly in, equilibrium. If the energy transfer rates are very slow the gas can be described as frozen. In all other instances, wherein any of the energy exchange rates are comparable to the local fluid time scale, the gas will be thermally or chemically reacting and out of equilibrium. Real gas thermochemical nonequilibrium processes are important in the determination of aerodynamic heating; both convective (including wall catalytic effects) and radiative heating. To illustrate this we consider the hypervelocity flow over a bluff body typical of an atmospheric entry vehicle or an aerospace transfer vehicle (ASTV). The qualitative aspects of a hypersonic flow field over a bluff body are discussed in two parts, forebody and afterbody, with attention to which particular physical effects must be included in an analysis. This will indicate what type of numerical modeling will be adequate in each region of the flow. A bluff forebody flow field is dominated by the presence of the strong bow shock wave and the consequent heating, and chemical reaction of the gas. At high altitude hypersonic flight conditions the thermal excitation and chemical reaction of the gas occur slowly enough that a significant portion of the flow field is in a state of thermochemical nonequilibrium. A second important effect is the presence of the thick boundary layer along the forebody surface. In this region there are large thermal and chemical species gradients due to the interaction of the gas with the wall. Also at high altitudes the shock wave and the boundary layer may become so thick that they merge; in this case the entire shock layer is dominated by viscous effects.

Deiwert, George S.↗

Growth and Morphology of Supercritical Fluids, a Fluid Physics Experiment Conducted on Mir, Complete

The Growth and Morphology of Supercritical Fluids (GMSF) is an international experiment facilitated by the NASA Glenn Research Center and under the guidance of U.S. principal investor Professor Hegseth of the University of New Orleans and three French coinvestigators: Daniel Beysens, Yves Garrabos, and Carole Chabot. The GMSF experiments were concluded in early 1999 on the Russian space station Mir. The experiments spanned the three science themes of near-critical phase separation rates, interface dynamics in near-critical boiling, and measurement of the spectrum of density fluctuation length scales very close to the critical point. The fluids used were pure CO2 or SF6. Three of the five thermostats used could adjust the sample volume with the scheduled crew time. Such a volume adjustment enabled variable sample densities around the critical density as well as pressure steps (as distinct from the usual temperature steps) applied to the sample. The French-built ALICE II facility was used for these experiments. It allows tightly thermostated (left photograph) samples (right photograph) to be controlled and viewed/measured. Its diagnostics include interferometry, shadowgraph, high-speed pressure measurements, and microscopy. Data were logged on DAT tapes, and PCMCIA cards and were returned to Earth only after the mission was over. The ground-breaking near critical boiling experiment has yielded the most results with a paper published in Physical Review Letters (ref. 1). The boiling work also received press in Science Magazine (ref. 2). This work showed that, in very compressible near-critical two-phase pure fluids, a vapor bubble was induced to temporarily overheat during a rapid heating of the sample wall. The temperature rise in the vapor was 23-percent higher than the rise in the driving container wall. The effect is due to adiabatic compression of the vapor bubble by the rapid expansion of fluid near the boundary during heatup. Thermal diffusivity is low near the critical point, so getting heat out of the compressed bubble is observably slow. This gives the appearance of a backward heat flow, or heat flow from a cold surface to a warm fluid.

Wilkinson, R. Allen↗

Advances in Design Capabilities for Planetary Missions from the NASA Entry Systems Modeling and Instrumentation Portfolio

The Entry Systems Modeling project (ESM) is supported by both the NASA Space Technology and the Science Mission Directorates and focuses on developing simulation tools and validated models for characterizing the performance of entry systems tailored to planetary destinations across the Solar System. ESM is organized into six technical capability areas that together address all relevant factors related to spacecraft entry, as well as some aspects of descent: Thermal Protection System (TPS) Materials; Aerothermodynamics; Entry & Descent Vehicle Dynamics; Guidance, Navigation, and Control; Vehicle Systems Analysis; and Advanced Tools and Numerical Methods. Development within the capability areas is undertaken explicitly with a focus on transition and infusion to science missions, human exploration missions, and commercial space activities. The present talk details developments that specifically impact science missions, including simulation tool capabilities that aid in mission design and model development to understand entry system performance at a given destination. Examples of the successful infusion and transition of such project outcomes to science missions also are provided. Several simulation tool development efforts within ESM have resulted in new design capabilities for missions. One such outcome is improved toolsets for mission trajectory and concept of operations design. Specifically, an initiative to couple a leading tool for entry, ascent/descent, and orbital trajectory optimization (Program to Optimize Simulated Trajectories II or POST2) to those used within the Agency for interplanetary trajectory optimization (Copernicus and Monte) has made substantial progress, with the outcomes to date promising to allow efficient trajectory optimization across mission phases. Additionally, toolchains for the evaluation of vehicle performance during entry and descent have been developed that allow assessment of multi-dimensional aeroheating on detailed vehicle geometries, characterization of deployment and inflation of parachutes, and assessment of vehicle dynamic stability during descent. These capabilities are achieved by coupling diverse sets of physics together – material response, computational fluid dynamics, radiation, and vehicle dynamics – to suitably describe complex entry and descent phenomena. Several model development and validation efforts for specific destinations and entry regimes also are underway within the ESM project. For instance, new experimental capabilities to validate radiation models at low densities/high altitudes recently have been established with project support, specifically the Low-Density Shock Tube (LDST) at the NASA Ames Research Center Electric Arc Shock Tube (EAST) facility. The LDST is being leveraged to develop improved models of shock layer kinetics and radiation in Titan atmospheres, while future studies will be conducted in the LDST and the existing high velocity shock tube to provide validation data for radiation models of Venus, Ice Giants, and Mars atmospheres. Models describing the aerothermal and thermo-structural performance of Thermal Protection System (TPS) materials has been another focus, with multiscale modeling activities on-going for the two leading TPS materials applicable to a range of entry conditions and science missions: the Phenolic-Impregnated Carbon Ablator (PICA) and woven materials like 3D Mid-Density Carbon Phenolic (3MDCP). A continual effort is made to infuse and transition outcomes from ESM simulation tool and model development activities into relevant science missions. Significant progress has been made on this front, with missions such as Dragonfly, DAVINCI, and Mars Missions benefitting from project outcomes. The groundwork also is being laid to provide insights into forward looking missions to Gas/Ice Giants as well as for potential sample returns.

Justin Haskins↗

A Code-Agnostic Driver Application for Coupled Neutronics and Thermal-Hydraulic Simulations

While the literature has numerous examples of Monte Carlo and computational fluid dynamics (CFD) coupling, most are hard-wired codes intended primarily for research rather than as standalone, general-purpose applications. In this work, we describe an open source application, ENRICO, that enables coupled neutronic and thermal-hydraulic simulations between multiple codes that can be chosen at runtime (as opposed to a coupling between two specific codes). The application has been designed such that the control flow logic, domain mapping, nonlinear fixed-point iteration, solution transfers, and convergence checks are all agnostic to the underlying physics solvers used. Special emphasis has also been placed on enabling efficient execution on distributed-memory computing environments. The transfer of solution fields between solvers is performed in memory rather than through filesystem I/O. Additionally, solvers can be configured to run on overlapping or disjoint sets of processes. To date, coupling with the OpenMC and Shift Monte Carlo codes, the Nek5000 CFD code, and a simplified heat diffusion and subchannel solver has been implemented in ENRICO. We present results for coupled simulations of a single light-water reactor fuel assembly based on the NuScale reactor using various combinations of the physics solvers. For this problem, the coupled simulations are shown to converge in about four Picard iterations. A comparison of the heat source and temperature distributions computed by ENRICO using OpenMC coupled with Nek5000 and Shift coupled with Nek5000 illustrates remarkable agreement between the codes.

42 ENGINEERING↗

A hierarchy for modeling high speed propulsion systems

General research efforts on reduced order propulsion models for control systems design are overviewed. Methods for modeling high speed propulsion systems are discussed including internal flow propulsion systems that do not contain rotating machinery such as inlets, ramjets, and scramjets. The discussion is separated into four sections: (1) computational fluid dynamics model for the entire nonlinear system or high order nonlinear models; (2) high order linearized model derived from fundamental physics; (3) low order linear models obtained from other high order models; and (4) low order nonlinear models. Included are special considerations on any relevant control system designs. The methods discussed are for the quasi-one dimensional Euler equations of gasdynamic flow. The essential nonlinear features represented are large amplitude nonlinear waves, moving normal shocks, hammershocks, subsonic combustion via heat addition, temperature dependent gases, detonation, and thermal choking.

Hartley, Tom T.↗

Thermal and Mechanical Microspacecraft Technologies for X-2000 Future Deliveries

Thermal and mechanical technologies are an important part of the X-2000 Future Delivery (X-2000 FD) microspacecraft. A wide range of future space missions are expected to utilize the technologies and the architecture developed by the X-2000 FD. These technologies, besides being small in physical size, make the tiny spacecraft robust and flexible. The X2000 FD architecture is designed to be highly reliable and suitable for a wide range of missions such as planetary landers/orbiters/flybys, earth orbiters, cometary flybys/landers/sample returns, etc. One of the key ideas used in the development of these technologies and architecture is that several functions be in included in each of the thermal and mechanical elements. One of the thermal architecture being explored for the X-2000 FD microspacecraft is integrated thermal energy management of the complete spacecraft using a fluid loop. The robustness and the simplicity of the loop and the flexibility with which it can be integrated in the spacecraft have made it attractive for applications to X-2000 FD. Some of the thermal technologies to be developed as a part of this architecture are passive and active cooling loops, electrically variable emittance surfaces, miniature thermal switches, and specific high density electronic cooling technologies. In the mechanical area, multifunction architecture for the structural elements will be developed. The multifunction aspect is expected to substantially reduce the mass and volume of the spacecraft. Some of the technologies that will be developed are composite material panels incorporating electronics, cabling, and thermal elements in them. The paper to be presented at the 1999 conference, will describe the progress made so far in the microspacecraft thermal and mechanical technologies and approaches for the X2000 Future Deliveries microspacecraft.

Birur, Gaj↗

Superfluid effective field theory for dark matter direct detection

We develop an effective field theory (EFT) framework for superfluid 4 He to model the interactions among quasiparticles, helium atoms and probe particles. Our effective field theory approach brings together symmetry arguments and power-counting and matches to classical fluid dynamics. We then present the decay and scattering rates for the relevant processes involving quasiparticles and helium atoms. The presented EFT framework and results can be used to understand the dynamics of thermalization in the superfluid, and can be further applied to sub-GeV dark matter direct detection with superfluid 4 He.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗