Search NASASearch

SEARCH · Search NASA

Results for “STEP FUNCTION”

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 127 records · Page 7

Electrodeposited CulnSe2 Thin Film Junctions

We have investigated thin films and junctions based on copper indium diselenide (CIS) which have been grown by electrochemical deposition. CIS is a leading candidate for use in polycrystalline thin film photovoltaic solar cells. Electrodeposition is a cost-effective method for producing thin-film CIS. We have produced both p and n type CIS thin films from the same aqueous solution by simply varying the deposition potential. A CIS pn junction was deposited using a step-function potential. Stoichiometry of the single layer films was determined by energy dispersive spectroscopy. Carrier densities of these films increased with deviation from stoichiometry, as determined by the capacitance versus voltage dependence of Schottky contacts. Optical bandgaps for the single layer films as determined by transmission spectroscopy were also found to increase with deviation from stoichiometry. Rectifying current versus voltage characteristics were demonstrated for the Schottky barriers and for the pn junction.

Raffaelle, R. P.

MOLA Topography of the Crustal Dichotomy Boundary Zone, Mars

Mars Orbiter Laser Altimeter (MOLA) profiles frequently cross the crustal dichotomy boundary where the transition zone (TZ) between cratered highland terrain (CT) and lowland smooth plains (SP) is marked by mesas and knobby terrain. The detailed topographic character of the boundary zone is longitudinally variable, as is the geomorphology of the TZ. Some portions of the boundary are associated with an outer ring of the Utopia impact basin; MOLA topography is consistent with this. The regional character of the boundary topography is a 2-4 km step function from nearly flat SP to almost as flat CT. This rise has a regional slope of 1-2 degrees, 50-100 times that of the Cr and SP away from TZ, which suggests a significant change in crustal properties (thickness, composition or both) across the TZ. The overall topography is very similar to that at some passive continent-oceanic crustal margins on the Earth, with the seafloor allowed to adjust upward after removal of the overlying water. A possible temporal constraint on the CT/SP elevation difference comes from two MOLA profiles which pass through two large (150 km diameter) craters located at the boundary in Aeolis. The N and S rims of the more degraded crater are at the same elevation; north of the N rim the topography drops by greater than 2 km to the floor of the TZ. This crater predates the elevation offset between CT and TZ floor. The better preserved crater (Gale) has a N rim 2 km lower than its S rim, and appears to have been emplaced on a pre-existing regional slope of about I degree. Gale probably post- dates the elevation difference between CT and TZ floor. Based on the stratigraphy of the units in which these craters are found, the elevation difference appears to have been in place in the Mid to Late Noachian.

Frey, Herbert V.

Co-Seismic Mass Dislocation and Its Effect on Earth's Rotation and Gravity

Mantle processes often involve large-scale mass transport, ranging from mantle convection, tectonic motions, glacial isostatic adjustment, to tides, atmospheric and oceanic loadings, volcanism and seismicity. On very short time scale of less than an hour, co-seismic event, apart from the "shaking" that is the earthquake, leaves behind permanent (step-function-like) dislocations in the crust and mantle. This redistribution of mass changes the Earth's inertia tensor (and hence Earth's rotation in both length-of-day and polar motion), and the gravity field (in terms of spherical harmonic Stokes coefficients). The question is whether these effects are large enough to be of any significance. In this paper we report updated calculation results. The calculation uses the normal mode summation scheme, applied to 15,814 major earthquakes that occurred during 1976-1998, according to source mechanism solutions given by the Harvard Central Moment Tensor catalog. Compared to the truly large ones earlier in the century, the earthquakes we study are individually all too small to have left any discernible signature in geodetic records of Earth rotation or global gravity field. However, their collective effects continue to exhibit an extremely strong statistical tendencies. For example, earthquakes conspire to decrease J(sub 2) and J(sub 22) while shortening LOD, resulting in a rounder and more compact Earth. Strong tendency is also seen in the earthquakes trying to "nudge" the Earth rotation pole towards about 140 degree E, roughly opposite to the observed polar drift direction. The geophysical significance and implications will be further studied.

Chao, Benjamin F.

Co-Seismic Mass Dislocation and its Effect on Earth's Rotation and Gravity

Mantle processes often involve large-scale mass transport, ranging from mantle convection, tectonic motions, glacial isostatic adjustment, to tides, atmospheric and oceanic loadings, volcanism and seismicity. On very short time scale of less than an hour, co-seismic event, apart from the shaking that is the earthquake, leaves behind permanent (step-function-like) dislocations in the crust and mantle. This redistribution of mass changes the Earth's inertia tensor (and hence Earth's rotation in both length-of-day and polar motion), and the gravity field (in terms of spherical harmonic Stokes coefficients). The question is whether these effects are large enough to be of any significance. In this paper we report updated calculation results based on Chao & Gross (1987). The calculation uses the normal mode summation scheme, applied to nearly twenty thousand major earthquakes that occurred during 1976-2002, according to source mechanism solutions given by the Harvard Central Moment Tensor catalog. Compared to the truly large ones earlier in the century, the earthquakes we study are individually all too small to have left any discernible signature in geodetic records of Earth rotation or global gravity field. However, their collective effects continue to exhibit an extremely strong statistical tendencies. For example, earthquakes conspire to decrease J2 and J22 while shortening LOD, resulting in a rounder and more compact Earth. Strong tendency is also seen in the earthquakes trying to nudge the Earth rotation pole towards approximately 140 degrees E, roughly opposite to the observed polar drift direction. The geophysical significance and implications will be further studied.

Chao, B. F.

Finite Element and Plate Theory Modeling of Acoustic Emission Waveforms

A comparison was made between two approaches to predict acoustic emission waveforms in thin plates. A normal mode solution method for Mindlin plate theory was used to predict the response of the flexural plate mode to a point source, step-function load, applied on the plate surface. The second approach used a dynamic finite element method to model the problem using equations of motion based on exact linear elasticity. Calculations were made using properties for both isotropic (aluminum) and anisotropic (unidirectional graphite/epoxy composite) materials. For simulations of anisotropic plates, propagation along multiple directions was evaluated. In general, agreement between the two theoretical approaches was good. Discrepancies in the waveforms at longer times were caused by differences in reflections from the lateral plate boundaries. These differences resulted from the fact that the two methods used different boundary conditions. At shorter times in the signals, before reflections, the slight discrepancies in the waveforms were attributed to limitations of Mindlin plate theory, which is an approximate plate theory. The advantages of the finite element method are that it used the exact linear elasticity solutions, and that it can be used to model real source conditions and complicated, finite specimen geometries as well as thick plates. These advantages come at a cost of increased computational difficulty, requiring lengthy calculations on workstations or supercomputers. The Mindlin plate theory solutions, meanwhile, can be quickly generated on personal computers. Specimens with finite geometry can also be modeled. However, only limited simple geometries such as circular or rectangular plates can easily be accommodated with the normal mode solution technique. Likewise, very limited source configurations can be modeled and plate theory is applicable only to thin plates.

Prosser, W. H.

Some Useful Innovations with Trasys and Sinda-85

Several innovative methods have been used to allow more efficient and accurate thermal analysis using SINDA-85 and TRASYS, including model integration and reduction, planetary surface calculations, and model animation. Integration with other modeling and analysis codes allows an analyst to import a geometry from a solid modeling or computer-aided design (CAD) software package, rather than building the geometry "by hand." This is more efficient as well as potentially more accurate. However, the use of solid modeling software often generates large analytical models. The problem of reducing large models has been elegantly solved using the response of the transient derivative to a forcing step function. The thermal analysis of a lunar rover implemented two unusual features of the TRASYS/SINDA system. A little-known TRASYS routine SURFP calculates the solar heating of a rover on the lunar surface for several different rover positions and orientations. This is used not only to determine the rover temperatures, but also to automatically determine the power generated by the solar arrays. The animation of transient thermal results is an effective tool, especially in a vivid case such as the 14-day progress of the sun over the lunar rover. An animated color map on the solid model displays the progression of temperatures.

Amundsen, Ruth M.

Co-Seismic Mass Displacement and its Effect on Earth's Rotation and Gravity

Mantle processes often involve large-scale mass transport, ranging from mantle convection, tectonic motions, glacial isostatic adjustment, to tides, atmospheric and oceanic loadings, volcanism and seismicity. On very short time scale of less than an hour, co-seismic event, apart from the "shaking" that is the earthquake, leaves behind permanent (step-function-like) displacements in the crust and mantle. This redistribution of mass changes the Earth's inertia tensor (and hence Earth's rotation in both length-of-day and polar motion), and the gravity field. The question is whether these effects are large enough to be of any significance. In this paper we report updated calculation results based on Chao & Gross. The calculation uses the normal mode summation scheme, applied to over twenty thousand major earthquakes that occurred during 1976-2002, according to source mechanism solutions given by the Harvard Centroid Moment Tensor catalog. Compared to the truly large ones earlier in the century, the earthquakes we study are individually all too small to have left any discernible signature in geodetic records of Earth rotation or global gravity field. However, their collective effects continue to exhibit an extremely strong statistical tendencies, conspiring to decrease J2 and J22 while shortening LOD, resulting in a rounder and more compact Earth. Strong tendency is also seen in the earthquakes trying to "nudge" the Earth rotation pole towards approx. 140 deg.E, roughly opposite to the observed polar drift direction. Currently, the Gravity Recovery And Climate Experiment (GRACE) is measuring the time-variable gravity to high degree and order with unprecedented accuracy. Our results show that great earthquakes such as the 1960 Chilean or 1964 Alaskan events cause gravitational field changes that are large enough to be detected by GRACE.

Chao, B. F.

Two-stage model of radon-induced malignant lung tumors in rats: effects of cell killing

A two-stage stochastic model of carcinogenesis is used to analyze lung tumor incidence in 3750 rats exposed to varying regimens of radon carried on a constant-concentration uranium ore dust aerosol. New to this analysis is the parameterization of the model such that cell killing by the alpha particles could be included. The model contains parameters characterizing the rate of the first mutation, the net proliferation rate of initiated cells, the ratio of the rates of cell loss (cell killing plus differentiation) and cell division, and the lag time between the appearance of the first malignant cell and the tumor. Data analysis was by standard maximum likelihood estimation techniques. Results indicate that the rate of the first mutation is dependent on radon and consistent with in vitro rates measured experimentally, and that the rate of the second mutation is not dependent on radon. An initial sharp rise in the net proliferation rate of initiated cell was found with increasing exposure rate (denoted model I), which leads to an unrealistically high cell-killing coefficient. A second model (model II) was studied, in which the initial rise was attributed to promotion via a step function, implying that it is due not to radon but to the uranium ore dust. This model resulted in values for the cell-killing coefficient consistent with those found for in vitro cells. An "inverse dose-rate" effect is seen, i.e. an increase in the lifetime probability of tumor with a decrease in exposure rate. This is attributed in large part to promotion of intermediate lesions. Since model II is preferable on biological grounds (it yields a plausible cell-killing coefficient), such as uranium ore dust. This analysis presents evidence that a two-stage model describes the data adequately and generates hypotheses regarding the mechanism of radon-induced carcinogenesis.

NASA Program Radiation Health

Implications of Dynamic Pressure Transducer Mounting Variations on Measurements in Pyrotechnic Test Apparatus

Accurate dynamic pressure measurements are often difficult to make within small pyrotechnic devices, and transducer mounting difficulties can cause data anomalies that lead to erroneous conclusions. Delayed initial pressure response followed by data ringing has been observed when using miniaturized pressure transducer mounting adapters required to interface transducers to small test chambers. This delayed pressure response and ringing, combined with a high data acquisition rate, has complicated data analysis. This paper compares the output signal characteristics from different pressure transducer mounting options, where the passage distance from the transducer face to the pyrotechnic chamber is varied in length and diameter. By analyzing the data and understating the associated system dynamics, a more realistic understanding of the actual dynamic pressure variations is achieved. Three pressure transducer mounting configurations (elongated, standard, and face/flush mount) were simultaneously tested using NASA standard initiators in closed volume pressure bombs. This paper also presents results of these pressure transducer mounting configurations as a result of a larger NASA Engineering and Safety Center pyrovalve test project. Results from these tests indicate the improved performance of using face/flush mounted pressure transducers in this application. This type of mounting improved initial pressure measurement response time by approximately 19 s over standard adapter mounting, eliminating most of the lag time; provided a near step-function type initial pressure increase; and greatly reduced data ringing in high data acquisition rate systems. The paper goes on to discuss other issues associated with the firing and instrumentation that are important for the tester to understand.

Dibbern, Andreas

A Test of Thick-Target Nonuniform Ionization as an Explanation for Breaks in Solar Flare Hard X-Ray Spectra

Solar nonthermal hard X-ray (HXR) flare spectra often cannot be fitted by a single power law, but rather require a downward break in the photon spectrum. A possible explanation for this spectral break is nonuniform ionization in the emission region. We have developed a computer code to calculate the photon spectrum from electrons with a power-law distribution injected into a thick-target in which the ionization decreases linearly from 100% to zero. We use the bremsstrahlung cross-section from Haug (1997), which closely approximates the full relativistic Bethe-Heitler cross-section, and compare photon spectra computed from this model with those obtained by Kontar, Brown and McArthur (2002), who used a step-function ionization model and the Kramers approximation to the cross-section. We find that for HXR spectra from a target with nonuniform ionization, the difference (Delta-gamma) between the power-law indexes above and below the break has an upper limit between approx.0.2 and 0.7 that depends on the power-law index delta of the injected electron distribution. A broken power-law spectrum with a. higher value of Delta-gamma cannot result from nonuniform ionization alone. The model is applied to spectra obtained around the peak times of 20 flares observed by the Ramaty High Energy Solar Spectroscopic Imager (RHESSI from 2002 to 2004 to determine whether thick-target nonuniform ionization can explain the measured spectral breaks. A Monte Carlo method is used to determine the uncertainties of the best-fit parameters, especially on Delta-gamma. We find that 15 of the 20 flare spectra require a downward spectral break and that at least 6 of these could not be explained by nonuniform ionization alone because they had values of Delta-gamma with less than a 2.5% probability of being consistent with the computed upper limits from the model. The remaining 9 flare spectra, based on this criterion, are consistent with the nonuniform ionization model.

Holman, gordon

Shuttle Risk Progression: Use of the Shuttle Probabilistic Risk Assessment (PRA) to Show Reliability Growth

It is important to the Space Shuttle Program (SSP), as well as future manned spaceflight programs, to understand the early mission risk and progression of risk as the program gains insights into the integrated vehicle through flight. The risk progression is important to the SSP as part of the documentation of lessons learned. The risk progression is important to future programs to understand reliability growth and the first flight risk. This analysis uses the knowledge gained from 30 years of operational flights and the current Shuttle PRA to calculate the risk of Loss of Crew and Vehicle (LOCV) at significant milestones beginning with the first flight. Key flights were evaluated based upon historical events and significant re-designs. The results indicated that the Shuttle risk tends to follow a step function as opposed to following a traditional reliability growth pattern where risk exponentially improves with each flight. In addition, it shows that risk can increase due to trading safety margin for increased performance or due to external events. Due to the risk drivers not being addressed, the risk did not improve appreciably during the first 25 flights. It was only after significant events occurred such as Challenger and Columbia, where the risk drivers were apparent, that risk was significantly improved. In addition, this paper will show that the SSP has reduced the risk of LOCV by almost an order of magnitude. It is easy to look back afte r 30 years and point to risks that are now obvious, however; the key is to use this knowledge to benefit other programs which are in their infancy stages. One lesson learned from the SSP is understanding risk drivers are essential in order to considerably reduce risk. This will enable the new program to focus time and resources on identifying and reducing the significant risks. A comprehensive PRA, similar to that of the Shuttle PRA, is an effective tool quantifying risk drivers if support from all of the stakeholders is given.

Hamlin, Teri L.

Big Impacts and Transient Oceans on Titan

We have studied the thermal consequences of very big impacts on Titan [1]. Titan's thick atmosphere and volatile-rich surface cause it to respond to big impacts in a somewhat Earth-like manner. Here we construct a simple globally-averaged model that tracks the flow of energy through the environment in the weeks, years, and millenia after a big comet strikes Titan. The model Titan is endowed with 1.4 bars of N2 and 0.07 bars of CH4, methane lakes, a water ice crust, and enough methane underground to saturate the regolith to the surface. We assume that half of the impact energy is immediately available to the atmosphere and surface while the other half is buried at the site of the crater and is unavailable on time scales of interest. The atmosphere and surface are treated as isothermal. We make the simplifying assumptions that the crust is everywhere as methane saturated as it was at the Huygens landing site, that the concentration of methane in the regolith is the same as it is at the surface, and that the crust is made of water ice. Heat flow into and out of the crust is approximated by step-functions. If the impact is great enough, ice melts. The meltwater oceans cool to the atmosphere conductively through an ice lid while at the base melting their way into the interior, driven down in part through Rayleigh-Taylor instabilities between the dense water and the warm ice. Topography, CO2, and hydrocarbons other than methane are ignored. Methane and ethane clathrate hydrates are discussed quantitatively but not fully incorporated into the model.

methane

Enhancements of Bayesian Blocks; Application to Large Light Curve Databases

Bayesian Blocks are optimal piecewise linear representations (step function fits) of light-curves. The simple algorithm implementing this idea, using dynamic programming, has been extended to include more data modes and fitness metrics, multivariate analysis, and data on the circle (Studies in Astronomical Time Series Analysis. VI. Bayesian Block Representations, Scargle, Norris, Jackson and Chiang 2013, ApJ, 764, 167), as well as new results on background subtraction and refinement of the procedure for precise timing of transient events in sparse data. Example demonstrations will include exploratory analysis of the Kepler light curve archive in a search for "star-tickling" signals from extraterrestrial civilizations. (The Cepheid Galactic Internet, Learned, Kudritzki, Pakvasa1, and Zee, 2008, arXiv: 0809.0339; Walkowicz et al., in progress).

Kepler light curve archive

Calibrated Noise Measurements with Induced Receiver Gain Fluctuations

The lack of well-developed techniques for modeling changing statistical moments in our observations has stymied the application of stochastic process theory in science and engineering. These limitations were encountered when modeling the performance of radiometer calibration architectures and algorithms in the presence of non stationary receiver fluctuations. Analyses of measured signals have traditionally been limited to a single measurement series. Whereas in a radiometer that samples a set of noise references, the data collection can be treated as an ensemble set of measurements of the receiver state. Noise Assisted Data Analysis is a growing field of study with significant potential for aiding the understanding and modeling of non stationary processes. Typically, NADA entails adding noise to a signal to produce an ensemble set on which statistical analysis is performed. Alternatively as in radiometric measurements, mixing a signal with calibrated noise provides, through the calibration process, the means to detect deviations from the stationary assumption and thereby a measurement tool to characterize the signal's non stationary properties. Data sets comprised of calibrated noise measurements have been limited to those collected with naturally occurring fluctuations in the radiometer receiver. To examine the application of NADA using calibrated noise, a Receiver Gain Modulation Circuit (RGMC) was designed and built to modulate the gain of a radiometer receiver using an external signal. In 2010, an RGMC was installed and operated at the National Institute of Standards and Techniques (NIST) using their Noise Figure Radiometer (NFRad) and national standard noise references. The data collected is the first known set of calibrated noise measurements from a receiver with an externally modulated gain. As an initial step, sinusoidal and step-function signals were used to modulate the receiver gain, to evaluate the circuit characteristics and to study the performance of a variety of calibration algorithms. The receiver noise temperature and time-bandwidth product of the NFRad are calculated from the data. Statistical analysis using temporal-dependent calibration algorithms reveals that the natural occurring fluctuations in the receiver are stationary over long intervals (100s of seconds); however the receiver exhibits local non stationarity over the interval over which one set of reference measurements are collected. A variety of calibration algorithms have been applied to the data to assess algorithms' performance with the gain fluctuation signals. This presentation will describe the RGMC, experiment design and a comparative analysis of calibration algorithms.

Racette, Paul

Data, Meet Compute: NASA's Cumulus Ingest Architecture

NASA's Earth Observing System Data and Information System (EOSDIS) houses nearly 30PBs of critical Earth Science data and with upcoming missions is expected to balloon to between 200PBs-300PBs over the next seven years. In addition to the massive increase in data collected, researchers and application developers want more and faster access - enabling complex visualizations, long time-series analysis, and cross dataset research without needing to copy and manage massive amounts of data locally. NASA has looked to the cloud to address these needs, building its Cumulus system to manage the ingest of diverse data in a wide variety of formats into the cloud. In this talk, we look at what Cumulus is from a high level and then take a deep dive into how it manages complexity and versioning associated with multiple AWS Lambda and ECS microservices communicating through AWS Step Functions across several disparate installations

EOSDIS

Simulation of Stochastic Mud-Crack Damage Formation in an Environmental Barrier Coating

The integrated Finite Element Analysis–Micromechanics Analysis Code/Ceramics Analysis and Reliability Evaluation of Structures (FEAMAC/CARES) program was used to simulate the formation of mudflat-cracks from thermomechanical loading on a multi-layered Environmental Barrier Coating (EBC) system deposited on a ceramic substrate. FEAMAC/CARES combines MAC/GMC multiscale composite micromechanics code with CARES/Life probabilistic multiaxial failure criteria code and Abaqusfinite element analysis. In this work, step function elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio finite elements enabled the depiction of crack boundaries and formation of mudflat patterned damage. Demonstrated examples include finite element models of button–sized disk–shaped 3-D specimen, and a 2-D model of through-the-thickness cross-section. All models were subjected to a progressive cool down from 1300oC to room temperature. Mudflat crack damage in the coating system resulted from the buildup of residual tensile stresses between the individual material constituents from thermal expansion mismatch. A 2-parameter Weibull distribution characterized the coating layer stochastic strength response and the effect of the Weibull modulus on the formation of damage was studied here.

residual tensile stress

Weibull-Based Stochastic Simulation of Mud-Crack Damage Formation in an Environmental Barrier Coating

The integrated Finite Element Analysis Micromechanics Analysis Code/Ceramics Analysis and Reliability Evaluation of Structures (FEAMAC/CARES) program was used to simulate the formation of mudflat-cracks from thermomechanical loading on a multi-layered Environmental Barrier Coating (EBC) system deposited on a ceramic substrate. FEAMAC/CARES combines MAC/GMC (Micromechanics Analysis Code/Generalized Method of Cells) multiscale composite micromechanics code with CARES/Life probabilistic multiaxial failure criteria code and Abaqusfinite element analysis. In this work, step function elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio finite elements enabled the depiction of crack boundaries and formation of mudflat patterned damage. Demonstrated examples include finite element models of button-sized disk-shaped 3-D specimen, and a 2-D model of through-the-thickness cross-section. All models were subjected to a progressive cool down from 1300 degrees Centigrade to room temperature. Mudflat crack damage in the coating system resulted from the buildup of residual tensile stresses between the individual material constituents from thermal expansion mismatch. A 2-parameter Weibull distribution characterized the coating layer stochastic strength response and the effect of the Weibull modulus on the formation of damage was studied here.

Residual Tensile Stress

Analysis of Oxygen Liquefaction with Transient Flow Rates for ISRU Systems

Human occupation of lunar and Martian surfaces requires in-situ resource utilization (ISRU) to create a sustainable environment with the limited resources available in space. Fuel and oxidizer generation is essential for developing a refueling capability for spacecraft, significantly reducing the propellant mass required for landing. Liquefaction systems are a key step in producing and storing cryogenic liquids such as oxygen, hydrogen, and methane. All those fluids are critical to propulsion, life support, and other spacecraft systems. NASA has previously confirmed the presence of water on the Moon and an electrolysis process can be used to separate the oxygen from the hydrogen molecules. Oxygen has also been found within lunar regolith and can be separated through ISRU processes. The extracted oxygen gas can then be liquefied inside a storage tank with tube-on-tank heat exchangers for future use as an oxidizer for propellants. While steady state liquefaction performance is easy to model, power limitations and cyclical environments caused by the change in Sun exposure between day and night periods potentially call for transient operations. The following analysis was performed using the Thermal Desktop software to investigate how various transient gaseous oxygen (GOX) flow rates impact the production and storage of liquid oxygen (LOX) in a 1-g liquefaction system. Variations of sinusoidal and exponential step functions were selected to model potential GOX flow rates that could be experienced by ISRU systems on lunar or Martian surfaces. This early computational analysis provides some insight on how transient operations impact oxygen liquefaction systems and helps illustrate operational questions to explore experimentally throughout the design process.

Cryogenic Fluid Management