Search NASA⌕ Search

SEARCH · Search NASA

Results for “Regularly Separable”

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.

239 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Experimental Investigation of Pool Boiling Heat Transfer Enhancement in Microgravity in the Presence of Electric Fields

In boiling high heat fluxes are possible driven by relatively small temperature differences, which make its use increasingly attractive in aerospace applications. The objective of the research is to develop ways to overcome specific problems associated with boiling in the low gravity environment by substituting the buoyancy force with the electric force to enhance bubble removal from the heated surface. Previous studies indicate that in terrestrial applications nucleate boiling heat transfer can be increased by a factor of 50, as compared to values obtained for the same system without electric fields. The goal of our research is to experimentally explore the mechanisms responsible for EHD heat transfer enhancement in boiling in low gravity conditions, by visualizing the temperature distributions in the vicinity of the heated surface and around the bubble during boiling using real-time holographic interferometry (HI) combined with high-speed cinematography. In the first phase of the project the influence of the electric field on a single bubble is investigated. Pool boiling is simulated by injecting a single bubble through a nozzle into the subcooled liquid or into the thermal boundary layer developed along the flat heater surface. Since the exact location of bubble formation is known, the optical equipment can be aligned and focused accurately, which is an essential requirement for precision measurements of bubble shape, size and deformation, as well as the visualization of temperature fields by HI. The size of the bubble and the frequency of bubble departure can be controlled by suitable selection of nozzle diameter and mass flow rate of vapor. In this approach effects due to the presence of the electric field can be separated from effects caused by the temperature gradients in the thermal boundary layer. The influence of the thermal boundary layer can be investigated after activating the heater at a later stage of the research. For the visualization experiments a test cell was developed. All four vertical walls of the test cell are transparent, and they allow transillumination with laser light for visualization experiments by HI. The bottom electrode is a copper cylinder, which is electrically grounded. The copper block is heated with a resistive heater and it is equipped with 6 thermocouples that provide reference temperatures for the measurements with HI. The top electrode is a mesh electrode. Bubbles are injected with a syringe into the test cell through the bottom electrode. The working fluids presently used in the interferometric visualization experiments, water and PF 5052, satisfy requirements regarding thermophysical, optical and electrical properties. A 30kV power supply equipped with a voltmeter allows to apply the electric field to the electrodes during the experiments. The magnitude of the applied voltage can be adjusted either manually or through the LabVIEW data acquisition and control system connected to a PC. Temperatures of the heated block are recorded using type-T thermocouples, whose output is read by a data acquisition system. Images of the bubbles are recorded with 35mm photographic and 16mm high-speed cameras, scanned and analyzed using various software packages. Visualized temperature fields HI allows the visualization of temperature fields in the vicinity of bubbles during boiling in the form of fringes. Typical visualized temperature distributions around the air bubbles injected into the thermal boundary layer in PF5052 are shown. The temperature of the heated surface is 35 C. The temperature difference for a pair of fringes is approximately 0.05 C. The heat flux applied to the bottom surface is moderate, and the fringe patterns are regular. In the image a bubble penetrating the thermal boundary layer is visible. Because of the axial symmetry of the problem, simplified reconstruction techniques can be applied to recover the temperature field. The thermal plume developing above the heated surface for more intensive heating is shown. The temperature distribution in the liquid is clearly 3D, and tomographic techniques have to be applied to recover the temperature distribution in such a physical situation. A sequence of interferometric images showing the temperature distribution around the rising bubble, recorded with a high-speed camera is shown. Again, the temperature distribution is 3D, and a more complex approach to the evaluation, the tomographic reconstruction has to be taken. Measurement of the temperature distribution from the fringe pattern temperature distributions that yield important information regarding heat transfer are determined. Two algorithms that allow the quantitative evaluation of interferometric fringe patterns and the reconstruction of temperature fields during boiling have been developed at the Heat Transfer Laboratory of the Johns Hopkins University. In the first algorithm the bubble is assumed to be axially symmetrical, which significantly reduces the computational effort for quantifying the temperature distribution around the bubble. For this purpose the thermal boundary layer around the bubble is divided into equidistant concentric shells, and the refractive index is assumed to be constant in each of the shells. Since large temperature gradients are expected in the vicinity of the bubble during boiling, the deflection of the light beam cannot be neglected in boiling experiments. Since the exit angle of the light beam is known, this allows to account for the deflections and phase shifts outside the boundary layer (in the bulk fluid and in the windows of the test cell). Three dimensional temperature distributions in the vicinity of the bubble are reconstructed using tomographic techniques. In tomography, the measurement volume is sliced into 2D planes. In the present study these planes are parallel to the heated surface. The objective is to determine the values of the field parameter of interest in form of the field function in these 2D planes. The field parameter is the change of the refractive index of the liquid in the measurement volume caused by temperature changes. By superimposing data for many 2D planes recorded at the same time instant, the 3D temperature distribution in the measurement volume is recovered.

Herman, Cila↗
DOE OSTI2022

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Evaluation of Zero-Net-Rate Pumping Tests

Accurately estimating the distribution of aquifer properties is key to understanding contaminant movement in the subsurface. The distribution of aquifer properties is typically addressed using slug or constant-rate well tests, and the pros and cons of these tests are well known. Slug tests are appealing because they avoid removing contaminated water, but their results are affected by well skin and the small volume of displaced water limits the volume of aquifer that can be evaluated. Constant-rate well tests have the disadvantage of requiring disposal of potentially contaminated water, but they can generate properties that are more representative than slug tests, and they can be used to estimate well efficiency and storativity, which are difficult to characterize using slug tests. Periodic pumping tests are appealing because they have many of the advantages and few of the disadvantages of slug and constant-rate well tests. Periodic pumping tests involve cycling the pumping rate with a regular period and measuring the resulting response in monitoring wells. Some of these tests involve imposing a periodic rate on a constant mean rate. However, other tests involve moving water out and back into the well at rates that are balanced so the net rate after each period is zero. Some zero-net-rate (ZNR) well tests use rates that follow a sinusoidal pattern, whereas others use a square wave pattern that switches between a constant rate of pumping and a constant rate of injection to achieve a zero-net rate. ZNR well tests appear to be a useful compromise between slug and constant rate tests, but methods for conducting and analyzing the results from these tests have received limited evaluation. Periodic pumping tests are typically evaluated through recorded pressure responses to a disturbance source. Recent work has demonstrated that the strain field in the vadose zone shows a response to disturbances in the underlying aquifer. Various instrument designs are available that can record vertical and horizontal strain to very precise resolutions (10-9 ε). Measuring the strain in the vadose zone can present lower costs than measuring pressure in a monitoring well, so using strain could improve the resolution of ZNR tests, as well as make these tests potentially cheaper. A ZNR periodic pumping system was constructed and used to perform tests in the Clemson, South Carolina area. The system was designed to pump and inject with a periodic square-wave, so pumping occurs at a constant rate for half the period and is followed by injection at the same rate for half the period resulting in no net gain or loss of water from the aquifer. The system is designed to generate flow rates from 1 to 3.5 gallons per minute (gpm) with a capacity of 875-gallon per half period. Pressure and strain monitoring points are located in the vicinity of the pumping well. A series of 10 ZNR periodic pumping tests were conducted at the field site, with periods ranging from 2 to 540 minutes. Pressure data were collected and recorded in four monitoring wells in various locations around the pumping well. In addition, vertical and horizontal strain and tilt data were recorded in various locations around the pumping well. Traditional aquifer tests, including slug tests and constant-rate pumping tests were conducted at the field site to provide a baseline of aquifer parameter estimates that are compared the results of ZNR tests. Pressure data measured at monitoring wells during ZNR tests were analyzed utilizing an analytical solution (Streltsova, 1988) that assumes a confined aquifer and estimates hydraulic diffusivity and transmissivity using the time lag and amplitude of the pressure signal. The storativity can be separated out from the hydraulic diffusivity using the transmissivity estimate. The time-lag and amplitude of the pressure signals were estimated using a Fourier transform. The time-lag of the primary period of the pressure increases as a linear function of distance from the pumping well, and it increases as the square root of the period of the ZNR test over a range of periods spanning two orders of magnitudes. The first few harmonics follow similar trends, but the higher frequency harmonics diverge from this trend. Strain and tilt data from various instruments in the vadose and saturated zones are periodic during ZNR tests. The time lags of the strain components increase roughly linearly with distance from the well. Many of the time lags of the strain data from the vadose zone are similar to the time lags observed in the pressure data from a similar distance and pumping period. Theoretical experiments were performed to understand how pressure and strain responds to ZNR periodic pumping. The time-lags of the simulated pressure responses match the field data trends. The results were used to validate the Streltsova (1988) solution for periodic pumping to estimate hydraulic diffusivity. Hydraulic diffusivity was estimated to be 1.8 x 10-2 m2s-1 < Dh < 9.0 x 10-2 m2s-1 using pressure data from all the wells during ZNR tests. Transmissivity was estimated to be 0.8 x 10-4 m2s-1 < T < 4.3 x 10-4 m2s-1. Assuming Storativity = T/Dh gives 1.1 x 10-3 < S < 19 x10-3. Transmissivity is used to estimate the hydraulic conductivity, K, by dividing by the assumed aquifer thickness. Two constant-rate pumping tests were conducted and analyzed using two conceptual models: The Hantush (1961) solution was used to analyze data assuming confined conditions and the Neuman (1974) solution was used to analyze data assuming unconfined conditions. The Hantush (1961) solution estimated transmissivity to be 1.6 x 10-4 m2s-1 < T < 1.8 x 10-4 m2s-1 when using data from all the wells. Storativity was estimated to be 4.4 x 10-3 < S < 8.2 x 10-3. Assuming Dh = T/S gives 0.22 x 10-2 m2s-1 < Dh < 3.6 x 10-2 m2s-1. The Neuman (1974) solution estimated transmissivity to be 0.80 x 10-4 m2s-1 < T < 0.84 x 10-4 m2s-1. Total storativity (S + Sy) was estimated to be 22 x 10-3 < S < 110 x 10-3. Storativity excluding the Sy term was estimated to be 2.2 x 10-3 < S < 2.5 x 10-3. Assuming Dh = T/S gives 0.073 x 10-2 m2s-1 < S < 0.36 x 10-2 m2s-1. Two slug tests were performed on the pumping well, PW-2. The Bouwer and Rice (1976) solution for slug tests was used to analyze the data. This solution works for both confined and unconfined conditions, so both conceptual models were covered through the analysis. Hydraulic conductivity was estimated to be 1.95 x 10-6 ms-1 < K < 2.49 x 10-6 ms-1. Statistical comparison show generally no significant difference in the estimates of Dh, K, or S made using ZNR and conventional pumping tests.. However, the total storage estimated by the Neuman (1974) solution for unconfined settings is larger than that estimate using the Streltsove analysis applied to the ZNR tests. This is likely because the Neuman solution considers delayed yield from storage that the water table, whereas the Streltsova (1988) solution assumes confined conditions. The analysis also indicates that there is a statistically significant difference between parameters estimated with conventional slug tests and those measured with either constant-rate pumping tests or ZNR tests. Hydraulic diffusivity is estimated from the time-lag and the distance from the well. The time-lags of the strain data from the vadose zone are similar to the lags of the pressure data from a similar distance. This indicates that hydraulic diffusivities estimated from the strain data measured in the vadose zone would be similar to the estimates from the pressure data measured in the aquifer. There appear to be sign reversals in some of the strain data that were corrected to estimate the time lag. The dissolved oxygen (CO) concentration was measured during several ZNR tests to evaluate the feasibility of using the procedure to increase contaminant degradation kinetics that are related to DO. The results indicate that DO can be increased during the cycling of ZNR tests in some cases.

Smith-Jones, Austin↗
NASA NTRSDate not supplied

Demand Capacity Balancing at Vertiports for Initial Strategic Conflict Management of Urban Air Mobility Operations

Urban Air Mobility (UAM) is a new transportation concept that enables highly automated, cooperative, passenger or cargo-carrying air transportation services in and around urban areas. To achieve the high level of operational density and complexity desired by the UAM community, an airspace system that allows UAM operators to readily access and operate safely and efficiently in the airspace is needed. This airspace system will require air traffic management designed to reduce the risk of conflicts and loss of separation between UAM flights. In general, strategic conflict management is considered as the first layer of conflict management for safe flight operations to condition the traffic to reduce the need for airborne separation provision, the second layer of conflict management. Demand Capacity Balancing (DCB) is one of the concept components to achieve strategic conflict management. DCB strategically evaluates traffic demand and resource capacities to allow UAM operators to determine when, where and how they operate, while mitigating conflicting needs for airspace and vertiport capacity. DCB can be applied whenever UAM demand exceeds the capacity in airspace or at vertiports. As the UAM ecosystem evolves with advanced technologies and matured operational procedures, more complicated conflict management will likely be needed. In the current UAM ‘Concept of Operation (ConOps) 1.0’ operational stage defined by FAA, however, it will be meaningful to explore the demand capacity balancing at vertiports only, as an initial strategic conflict management approach for UAM operations because vertiport capacity seems to be a bottleneck of UAM traffic. For this research, we developed a demand-capacity imbalance detection and resolution service for UAM. This DCB service identifies the demand from operators and compares the demand to a given capacity at the shared resources (i.e., vertiports) over the upcoming time horizon which is divided into time bins having a constant interval. When a new flight plan is submitted, the algorithm embedded in the DCB service checks the available time bins based on the desired departure time and estimated arrival time at origin and destination vertiports, respectively. If the time bins for the originally desired times are already occupied by other flights (i.e., demand is at or above capacity), the algorithm finds the next available time bins for takeoff and landing and shifts the conflicting departure time to the earliest time that satisfies the capacity constraints at both origin and destination vertiports. The details of the algorithm will be described in the final manuscript. Figure 1 shows that the proposed DCB algorithm works well for a sample traffic scenario. In this example, a total of 144 flights, split between two operators, are planned over 2 hours, traveling 10 routes between five vertiports. In the heatmaps, the horizontal axis shows 12 time bins where each bin represents a 12-minute interval, and the vertical axis shows five vertiports. The number in each cell shows the number of operations, counting both departures and arrivals, at a specific vertiport in each time bin. For the given capacity of 2 operations/vertiport/bin, Figure 1 shows that the original demand sometimes exceeds the capacity, but the modified demand is reduced to the given capacity after resolving demand-capacity imbalances. When UAM flights are operated, it is expected that many practical issues would arise in the federated system architecture with multiple operators. UAM operators may experience a time synchronization issue due to communication delay between operator and vehicle. UAM vehicles would fly at different flight speeds, depending on vehicle models. Actual departure and arrival times can have large variations, compared to the schedule. The lead time from flight plan submission to desired departure time can vary by service type (e.g., regular shuttle service vs. on-demand service). Using the proposed DCB algorithm, we also investigated how the actual flight schedule and DCB performance are affected by these uncertainties such as unsynchronized times between operators, flight speed differences, lead time differences, and departure time errors. The final manuscript will include the background of this research work, the description of the DCB algorithm and its use cases with traffic scenarios. It will also provide the analytical results about the impact of various uncertainties that can occur in actual UAM operations on the DCB at vertiports, in terms of demand distribution changes, number of simultaneous operations, and delay propagation.

Urban Air Mobility↗
NASA NTRS2016

Investigating the Simulink Auto-Coding Process

Model based program design is the most clear and direct way to develop algorithms and programs for interfacing with hardware. While coding "by hand" results in a more tailored product, the ever-growing size and complexity of modern-day applications can cause the project work load to quickly become unreasonable for one programmer. This has generally been addressed by splitting the product into separate modules to allow multiple developers to work in parallel on the same project, however this introduces new potentials for errors in the process. The fluidity, reliability and robustness of the code relies on the abilities of the programmers to communicate their methods to one another; furthermore, multiple programmers invites multiple potentially differing coding styles into the same product, which can cause a loss of readability or even module incompatibility. Fortunately, Mathworks has implemented an auto-coding feature that allows programmers to design their algorithms through the use of models and diagrams in the graphical programming environment Simulink, allowing the designer to visually determine what the hardware is to do. From here, the auto-coding feature handles converting the project into another programming language. This type of approach allows the designer to clearly see how the software will be directing the hardware without the need to try and interpret large amounts of code. In addition, it speeds up the programming process, minimizing the amount of man-hours spent on a single project, thus reducing the chance of human error as well as project turnover time. One such project that has benefited from the auto-coding procedure is Ramses, a portion of the GNC flight software on-board Orion that has been implemented primarily in Simulink. Currently, however, auto-coding Ramses into C++ requires 5 hours of code generation time. This causes issues if the tool ever needs to be debugged, as this code generation will need to occur with each edit to any part of the program; additionally, this is lost time that could be spent testing and analyzing the code. This is one of the more prominent issues with the auto-coding process, and while much information is available with regard to optimizing Simulink designs to produce efficient and reliable C++ code, not much research has been made public on how to reduce the code generation time. It is of interest to develop some insight as to what causes code generation times to be so significant, and determine if there are architecture guidelines or a desirable auto-coding configuration set to assist in streamlining this step of the design process for particular applications. To address the issue at hand, the Simulink coder was studied at a foundational level. For each different component type made available by the software, the features, auto-code generation time, and the format of the generated code were analyzed and documented. Tools were developed and documented to expedite these studies, particularly in the area of automating sequential builds to ensure accurate data was obtained. Next, the Ramses model was examined in an attempt to determine the composition and the types of technologies used in the model. This enabled the development of a model that uses similar technologies, but takes a fraction of the time to auto-code to reduce the turnaround time for experimentation. Lastly, the model was used to run a wide array of experiments and collect data to obtain knowledge about where to search for bottlenecks in the Ramses model. The resulting contributions of the overall effort consist of an experimental model for further investigation into the subject, as well as several automation tools to assist in analyzing the model, and a reference document offering insight to the auto-coding process, including documentation of the tools used in the model analysis, data illustrating some potential problem areas in the auto-coding process, and recommendations on areas or practices in the current Ramses model that should be further investigated. Several skills were required to be built up over the course of the internship project. First and foremost, my Simulink skills have improved drastically, as much of my experience had been modeling electronic circuits as opposed to software models. Furthermore, I am now comfortable working with the Simulink Auto-coder, a tool I had never used until this summer; this tool also tested my critical thinking and C++ knowledge as I had to interpret the C++ code it was generating and attempt to understand how the Simulink model affected the generated code. I had come into the internship with a solid understanding of Matlab code, but had done very little in using it to automate tasks, particularly Simulink tasks; along the same lines, I had rarely used shell script to automate and interface with programs, which I gained a fair amount of experience with this summer, including how to use regular expression. Lastly, soft-skills are an area everyone can continuously improve on; having never worked with NASA engineers, which to me seem to be a completely different breed than what I am used to (commercial electronic engineers), I learned to utilize the wealth of knowledge present at JSC. I wish I had come into the internship knowing exactly how helpful everyone in my branch would be, as I would have picked up on this sooner. I hope that having gained such a strong foundation in Simulink over this summer will open the opportunity to return to work on this project, or potentially other opportunities within the division. The idea of leaving a project I devoted ten weeks to is a hard one to cope with, so having the chance to pick up where I left off sounds appealing; alternatively, I am interested to see if there are any opening in the future that would allow me to work on a project that is more in-line with my research in estimation algorithms. Regardless, this summer has been a milestone in my professional career, and I hope this has started a long-term relationship between JSC and myself. I really enjoy the thought of building on my experience here over future summers while I work to complete my PhD at Missouri University of Science and Technology.

Gualdoni, Matthew J.↗