Search NASA⌕ Search

SEARCH · Search NASA

Results for “First-principles calculations”

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.

613 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Helmholtz Motor: A Novel Electric Machine That Enables High Temperature Superconducting Armatures

Electrified aircraft are being developed to increase the efficiency and reduce the cost of operating subsonic transport aircraft. Achieving a substantial impact necessitates focusing on single- and twin-aisle aircraft which use propulsion systems with >20 MW ratings. For these aircraft, multi-MW superconducting electric machines are being developed due to their high specific power and high efficiency. Cryogenic electric machines are a promising technology for multi megawatt electric aircraft drivetrains. Second generation high temperature superconductors (HTS) are the ideal conductor for the field winding (typically the rotor) of superconducting machines, and they have attractive features for the armature winding (typically the stator) of these machines. However, their use in armature windings has been severely limited due to excessive AC losses in conventional electric machines. This presentation will present the working principle, design, and a performance trade study of a novel electric machine – termed a Helmholtz machine – that enables HTS armatures by significantly reducing AC losses. The working principle of the patent pending Helmholtz motor is to optimize the rotor-produced magnetic field (i.e., the field winding) so that it produces a magnetic field almost purely in the plane of the HTS. This is accomplished using 2 rotors that contain a set of matched Helmholtz coil pairs (or matched permanent magnet Halbach arrays). The armature (stator) HTS coils are positioned between the rotors and oriented to minimize the out of plane magnetic field. The armature HTS coils are also designed to minimize the out of plane component of their self field. The presentation will include the lessons learned in the sizing and design of this motor, including equations for how to match the magnet arrays. A performance trade study will be presented for a partially superconducting, radial flux version of the motor. The trade study involved an optimization of a motor design code based on analytical calculations (electromagnetic, AC loss, thermofluid, mechanical, and motor sizing) that determine the motor’s total efficiency and total specific power. The design code will be described in the presentation, with an emphasis on the calculation of magnetic fields, AC loss, and temperature distribution.

Superconducting electric machines↗
NASA NTRS1985

Nondestructive Evaluation of Composites By Optical Time Domain Reflectometry

Results in the measurement of slowly varying mechanical loading functions on composites using optical time domain reflectometry (OTDR) in imbedded optical fiber during both simulated manufacture and use are reviewed. First, the basic theoretical and experimental principles of the OTDR system are described. Next, the mechanical system of the composite and the imbedded fiber is analyzed. Results of measurement obtained for various loading functions applied to material specimens are then described and system limitations on spatial resolution, strain amplitude sensitivity, and frequency response are explained in terms of the range and resolution limits of the OTDR system.

R O Claus↗
DOE OSTI2026

Multiple shape coexistence near {sup 118}Sn: First 0+ 3 lifetime measurement

The intruder bands in Sn isotopes, built on the 2p-2h excitation across the Z = 50 proton shell gap, are well-known examples of shape coexistence near the neutron mid-shell region. Spectroscopic signatures for shape coexistence include enhanced E0 transitions between the 0+ band heads. However, the underlying shape coexistence and mixing has been unclear because lifetime information for the excited 0+ states was incomplete in 118Sn. We thus present here the first measurement of the 0+ 3 lifetime in 118Sn using the fast-timing technique following thermal-neutron capture. The observed enhancement in ρ2(E0;0+ 3 → 0+ 2 ) of 150(30) milliunits provides compelling indications for multiple shape coexistence in 118Sn. Additionally, three distinct shapes in 116,118,120Sn naturally emerged in theoretical calculations based on the quantum-number-projected generator coordinate method employing a relativistic energy density functional.

Wu, F.↗
DOE OSTI2024

Unravelling ultralow thermal conductivity in perovskite Cs2AgBiBr6: dominant wave-like phonon tunnelling and strong anharmonicity

Abstract Understanding the lattice dynamics and heat transport physics in the lead-free halide double perovskites remains an outstanding challenge due to their lattice dynamical instability and strong anharmonicity. In this work, we investigate the microscopic mechanisms of anharmonic lattice dynamics and thermal transport in lead-free halide double perovskite Cs 2 AgBiBr 6 from first principles. We combine self-consistent phonon calculations with bubble diagram correction and a unified theory of lattice thermal transport that considers both the particle-like phonon propagation and wave-like tunnelling of phonons. An ultra-low thermal conductivity at room temperature (~0.21 Wm −1 K −1 ) is predicted with weak temperature dependence( ~ T −0.34 ), in sharp contrast to the conventional ~T −1 dependence. Particularly, the vibrational properties of Cs 2 AgBiBr 6 are featured by strong anharmonicity and wave-like tunnelling of phonons. Anharmonic phonon renormalization from both the cubic and quartic anharmonicities are found essential in precisely predicting the phase transition temperature in Cs 2 AgBiBr 6 while the negative phonon energy shifts induced by cubic anharmonicity has a significant influence on particle-like phonon propagation. Further, the contribution of the wave-like tunnelling to the total thermal conductivity surpasses that of the particle-like propagation above around 310 K, indicating the breakdown of the phonon gas picture conventionally used in the Peierls-Boltzmann Transport Equation. Importantly, further including four-phonon scatterings is required in achieving the dominance of wave-like tunnelling, as compared to the dominant particle-like propagation channel when considering only three-phonon scatterings. Our work highlights the importance of lattice anharmonicity and wave-like tunnelling of phonons in the thermal transport in lead-free halide double perovskites.

Chemistry↗
DOE OSTI2024

Converting Second‐Order Saddle Points to Transition States: New Principles for the Design of 4π Photoswitches

Abstract Molecular photoswitches have demonstrated potential for storing solar energy at the molecular level, with power densities comparable to commercial batteries and hydroelectric energy storage. However, development of efficient photoswitches is hindered by limitations in cyclability and optical properties of existing materials. We here demonstrate that certain limitations in photoswitches based on electrocyclizations stem from the issue of controlling competition between Woodward‐Hoffmann allowed and forbidden pathways. Our approach moves beyond the traditional view of activation barriers and reveals that second‐order saddle points are crucial in dictating the competition between disrotatory and conrotatory pathways. These insights suggest new opportunities to manipulate the competition between these pathways through geometric constraints, fundamentally altering the connectivity of the potential energy surface. Our study also emphasizes the necessity of multi‐reference methods and the need to conduct higher‐dimensional explorations for competing pathways beyond photoswitch design.

Chemistry↗
DOE OSTI2025

Deriving the Landauer Principle From the Quantum Shannon Entropy

We derive an expression to determine the equilibrium probability distribution of a quantum state in contact with a noisy thermal environment that formally separates contributions from quantum and classical forms of probabilistic uncertainty. A statistical mechanical interpretation of this probability distribution enables us to derive an expression for the minimum free energy costs for arbitrary (reversible or irreversible) quantum state changes. In conclusion, based on this derivation, we demonstrate that–in contrast to classical systems–the free energy required to erase or reset a qubit depends sensitively on both the fidelity of the target state and on the physical properties of the environment, such as the number of quantum bath states, due primarily to the entropic effects of system-bath entanglement.

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

Energy-Optimal Vehicle Longitudinal Motion Control via Pontryagin’s Minimum Principle and Ultra-Local Model

Longitudinal vehicle motion control is essential for enhancing performance and optimizing a vehicle’s energy usage. However, it remains a challenging task due to the nonlinear and uncertain nature of vehicle dynamics, along with varying driving conditions. This paper presents a novel ultra-local optimal control approach based on Pontryagin’s Minimum Principle (PMP) that circumvents the need for detailed system identification by employing an ultra-local model. The control objective is to minimize the total energy consumption under boundary conditions while ensuring smooth traction force generation. The proposed approach is evaluated using a high-fidelity vehicle model in three representative scenarios: (i) nominal driving, (ii) a change in tire road friction coefficient (TRFC) from 0.5 to 0.65 and road slope from 0% to 5% during the maneuver, with target velocity unchanged, and (iii) a change in target velocity from 20 m/s to 0 m/s during the maneuver, while maintaining nominal TRFC and slope conditions. The simulation results demonstrate that the proposed method delivers robust performance, effectively balancing consumption and tracking accuracy in all tested scenarios.

Waleed khan, Muhammad [The University of Texas at ↗
DOE OSTI2025

PAH101: A GW+BSE Dataset of 101 Polycyclic Aromatic Hydrocarbon (PAH) Molecular Crystals

Abstract The excited-state properties of molecular crystals are important for applications in organic electronic devices. TheGWapproximation and Bethe-Salpeter equation (GW+BSE) is the state-of-the-art method for calculating the excited-state properties of crystalline solids with periodic boundary conditions. We present the PAH101 dataset ofGW+BSE calculations for 101 molecular crystals of polycyclic aromatic hydrocarbons (PAHs) with up to ~500 atoms in the unit cell. To the best of our knowledge, this is the firstGW+BSE dataset for molecular crystals. The data records include theGWquasiparticle band structure, the fundamental band gap, the static dielectric constant, the first singlet exciton energy (optical gap), the first triplet exciton energy, the dielectric function, and optical absorption spectra for light polarized along the three lattice vectors. The dataset can be used to (i) discover materials with desired electronic/optical properties, (ii) identify correlations between DFT andGW+BSE quantities, and (iii) train machine learned models to help in materials discovery efforts.

Science & Technology - Other Topics↗
NASA NTRS1977

Calculations of the Moon's Thermal History at Different Concentrations of Radioactive Elements, Taking Into Account Differentiation on Melting

Calculations of the thermal history of the moon were carried out by solving the thermal conductivity equation for the case in which the heat sources are the long-lived radioactive elements Th, U, and K 40 . The concentrations of these elements were adjusted to give four variations of the heat flow: 1.35 x 10 -8 and 0.91 x 10 -6 cal cm -2 s -1 (I and I**), 0.61 x 10 -6 (variant II, the terrestrial mixture of Lyubimova and Starkova), and 0.236 x 10 -6 (variant III, the chondrite model of Urey and MacDonald). In the same calculations, we considered layering of the differentiated material with transport to the surface of the radioactive elements after the temperature of the layer rose to 200 K above the melting temperature, which is considered in five variants that differ in the amount of transported radioactive elements: 100 percent ( n = 1), 80 percent n = 0.8), 40 percent, 20 percent ( n = 0.6, 0.4, 0.2). During fusion the heat capacity and heat conductivity were changed. We considered two variants of an initially cold (273 K) and hot (900 K) Moon. Calculations show that the interior of the Moon was heated to melting during the first 0.7 to 2.3 x 10 9 years. The maximum fusion involved practically the entire Moon to a distance from 15 to 45 km beneath the surface, and started 3.5 to 4.0 x 10 9 years ago (I, I**), or 2.5 to 3.0 x 10 9 years ago (II, III) and continued for 1 to 2 x 10 9 years. Today the Moon is cooling. The current thickness of the solid crust is from 150 to 200 km and the heat flow exceeds the stationary value 1.5 fold. Apparently the most realistic variant is II (terrestrial mixture) for an initially hot Moon, and gives, regardless of the moderate concentration of radioactive elements, a heat flow of 0.9 to 0.95 x 10 -6 cal cm -2 s -1 , which agrees with the radioastronomical measurements of Troitsky and Krotikov and with the Apollo 15 data.

O I Ornatskaya↗
NASA NTRSDate not supplied

Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)

Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.

Bohdan O Wesely↗
NASA NTRS1960

A Coaxial Flow Reactor - A Gaseous Nuclear-Rocket Concept

Aspects of a new gaseous nuclear-rocket concept are considered. The principle of this type of cavity reactor is one of maintaining the separation of two coaxial streams - fuel within propellant - rather than attaining the separation of fuel propellant gases within the cavity. An acceptably low uranium consumption rate is possible if the fuel stream can be maintained at an axial velocity which is 1/100 to 1/1000 of the propellant velocity. An elementary hydrodynamic analysis for laminar incompressible flow indicates that the average uranium core velocity may only reach 2 to 11 times the initial core velocity for a typical reactor. The results of a 1-D criticality analysis for a D2O/(235)U system show that the effect of shrinking the fuel region within the reactor cavity is not serious; the effects of moderator heating and structural materials on critical mass are also shown. Overall rocket performance calculations are made for reasonable operating conditions; these results indicate that thrusts from 25,000 to 50,000 pounds and accelerations from 10(exp -1) to 3 g are possible.

Herbert Weinstein↗
DOE OSTI2024

The first high-redshift cavity power measurements of cool-core galaxy clusters with the International LOFAR Telescope

Radio-mode feedback associated with the active galactic nuclei (AGNs) at the cores of galaxy clusters injects a large amount of energy into the intracluster medium (ICM), offsetting radiative losses through X-ray emission. This mechanism prevents the ICM from rapidly cooling down and fueling extreme starburst activity as it accretes onto the central galaxies, and it is therefore a key ingredient in the evolution of galaxy clusters. However, the influence and mode of feedback at high redshifts (z ∼ 1) remains largely unknown. Low-frequency sub-arcsecond-resolution radio observations taken with the International LOFAR Telescope have demonstrated their ability to assist X-ray observations with constraining the energy output from the AGNs (or “cavity power”) in galaxy clusters, thereby enabling research at higher redshifts than before. In this pilot project, we tested this hybrid method on a high-redshift (0.6 < z < 1.3) sample of 13 galaxy clusters for the first time with the aim of verifying the performance of this method at these redshifts and providing the first estimates of the cavity power associated with the central AGN for a sample of distant clusters. We were able to detect clear radio lobes in three out of 13 galaxy clusters at redshifts of 0.7 < z < 0.9, and we used these detections in combination with ICM pressures surrounding the radio lobes obtained from standard profiles to calculate the corresponding cavity powers of the AGNs. Combining our results with the literature, the current data appear to suggest that the average cavity power peaked at a redshift ofz ∼ 0.4 and slowly decreases toward higher redshifts. However, we require more and tighter constraints on the cavity volume and a better understanding of our observational systematics to confirm any deviation of the cavity power trend from a constant level.

Astronomy & Astrophysics↗
DOE OSTI2026

Role of ion acoustic instability in magnetic reconnection

We report on a first-principles numerical study of magnetic reconnection in plasmas with different initial ion-to-electron temperature ratios. In cases where this ratio is significantly below unity, we observe intense wave activity in the diffusion region, driven by the ion-acoustic instability. Our analysis shows that the dominant macroscopic effect of this instability is to drive substantial ion heating. In contrast to earlier studies reporting significant anomalous resistivity, we find that anomalous contributions due to the ion-acoustic instability are minimal. These results shed light on the dynamical impact of this instability on reconnection processes, offering new insights into the fundamental physics governing collisionless reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

First N 2 Profile for Venus’ Deep Lower Atmosphere

We present the first N 2 profile for Venus’ deep lower atmosphere (<15 km) and a constrained isotopic composition for cloud N 2 [1]. This work directly addresses unresolved questions for Venus using legacy observations. Prior to this work, there were no reported measurements for N 2 abundances at <15 km and the isotopic composition remained unconstrained [2]. The N 2 parameters are critical to understanding the evolution and thermal properties of the atmosphere [3-7]. Our N 2 results were obtained by re-analysis of data acquired in 1978 by the Pioneer Venus Large Probe Neutral Mass Spectrometer (LNMS) [8]. The archived mass spectra from 64.2 to 0.2 km were treated using the analytical procedures specifically developed for the LNMS [9-11]. Judicious peak fitting permitted disambiguation of (A) N 2 + and CO + at 28 u and (B) 14 N 15 N+, 13 CO + , and C 2 H 5 + at 29 u. Quality controls included comparing the (A) LNMS CO + /CO 2 + ratios to the NIST database and literature and (B) fitted counts for CO + to the expected counts of CO + calculated from C 18 O + and 13 CO + using the LNMS 16 O/ 18 O and 12 C/ 13 C ratios (obtained from CO 2 ). Our results show that N 2 is uniformly mixed across the deep lower atmosphere between ~0.2 and 15 km (2.49 ± 0.10 v%). In contrast, N 2 is non-uniformly mixed across the sub-cloud atmosphere and clouds (~15–59 km), where N 2 abundances increase by ~ 2-fold between ~15 km (2.45 ± 0.32 v%) and ~59-51 km (5.21 ± 0.18 v%). Using the cloud data, we also obtained a constrained 15 N/ 14 N ratio (2.93×10 -3 ± 0.13×10 -3 ) and δ 15 N value (-204 ± 35‰). Thus, the LNMS results [1] suggest that (A) the atmosphere is unstable at <15 km, (B) N 2 is not well-mixed >15 km, and (C) the cloud δ 15 N falls between Earth and the solar wind [12, 13]. Comparisons of the N 2 abundances and isotopic compositions for nitrogen, carbon, and oxygen to other Venus measurements will be discussed.

deep lower atmosphere↗
NASA NTRSDate not supplied

Usage-based Lifing of Lithium-Ion Battery with HybridPhysics-Informed Neural Networks

Lithium-ion batteries are commonly used to power unmanned aircraft vehicles (UAVs).The ability to model and forecast the remaining useful life of these batteries enables UAV reliability assurance. Building accurate models for battery state of charge and state of health based on first principles is challenging due to the complex electrochemistry that governs battery operations and computational complexity required to solve them. Therefore, reduced order models are often used due to their ability to capture the overall battery discharge. Un-fortunately, these simplifications lead to residual discrepancy between model predictions and observed data. In this paper, we present a hybrid modeling approach merging reduced-order models and neural networks. In this approach, while most of the input-output relationship is captured by Nernst and Butler-Volmer equations, data-driven kernels reduce the gap between predictions and observations. We validate our approach using data publicly available through the NASA Prognostics Center of Excellence repository. Results showed that our hybrid battery prognosis model can be successfully calibrated, even with a limited number of observations.

Lithium-ion Battery↗
DOE OSTI2026

Working Fluid Characterization and Performance Assessment of Subcritical Organic Rankine Cycles Based on the Lee–Kesler Approach for Energy Recovery

Here, a generalized model using the Lee–Kesler approach based on the corresponding states principle is developed to assess the performance of subcritical Organic Rankine Cycles operating with different working fluids. Each fluid is characterized by five parameters: the acentric factor, critical temperature, critical pressure, molar mass, and the ideal-gas ratio of specific heats at the critical temperature. The model was developed using the compressibility factor modified version of the Benedict–Webb–Rubin equation proposed by Lee and Kesler and the enthalpy and entropy functions to calculate thermodynamic state properties. The model was validated by comparing the results calculated with the model and working fluid thermodynamic properties obtained with the CoolProp database. This comparison was conducted for 91 working fluids, obtaining a relative error below 5% for 88 out of the 91 fluids (∼97%). A generalized parametric study was conducted to determine the influence of the pinch point and each fluid parameter on the performance of Organic Rankine Cycle (ORC) systems. It was found that efficiency increases with critical temperature, ideal-gas ratio of specific heats at the critical temperature, and acentric factor, reaching up to 13%. The developed model enables the evaluation of ORC system performance for existing working fluids. It also allows the formulation and evaluation of new fluids to enhance the performance of the ORC while retrieving energy from any kind of source; and likewise, the methodology can be applied to other power generation cycles.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2014

Influential Factors for Liquid Acquisition Device Screen Selection for Cryogenic Propulsion Systems

This paper presents the influential factors which govern screen selection for liquid acquisition devices (LADs) operating in microgravity conditions for future in-space cryogenic propulsion engines and cryogenic propellant depots. Space flight requirements, which include mass flow rate, acceleration level and direction, and thermal environment, dictate screen selection for a particular mission. The five influential factors include bubble point pressure, flow-through-screen pressure drop, wicking rate, screen compliance, and material compatibility. Governing equations and analytical models for these parameters are developed from first principles. A comprehensive survey of the historical data on coarser LAD meshes over four decades of work is conducted, and liquid hydrogen data for finer Dutch Twill meshes (325 x 2300, 450 x 2750, 510 x 3600) from recently concluded experiments is also presented to validate analytical models. Each of these parameters is measurable from ground based tests, making it facile to predict flight system performance. Therefore analytical models in this paper will be valuable for future LAD designs for both cryogenic and storable propulsion systems. Additionally, analysis will be given on the impact of the factors on liquid hydrogen systems.

Fuel Depot↗