Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “method of multipliers”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Development of an apparatus for obtaining molecular beams in the energy range from 2 to 200 eV

The formation and detection of molecular beams obtained by charge exchange from a low-energy ion source is discussed. Dispersion in energy of the ion source was measured and problems concerning detection of neutral beams were studied. Various methods were used, specifically secondary electron emissivity of a metallic surface and ionization of a gas target with a low ionization voltage. The intensities of neutral beams as low as 10 eV are measured by a tubular electron multiplier and a lock-in amplifier.

Clapier, R.↗
NASA NTRS1985

Formulation and application of Russell's method

It is shown that the numerical technique of Russell's momentum approach can be derived by using Hamilton's principle and Vance's numerical scheme. It results in a set of first order differnce equations for solving the angular velocities. The numerical examples show that the method is reliable. The algorithm is modified next to perform the analysis of N-body systems with closed loop topology. To increase the formulation flexibility, the equations of motion are represented by using Cartesian coordinates and Lagrange multipliers. The algorithm consists of two parts, Vance's scheme and an unconstrained minimization. The Vance's scheme is used to find the angular velocities, and the unconstrained minimization is applied to provide the correct angular displacements. The proposed scheme is further extended to find the design sensitivity of an N-body system with closed loop configuration, and to carry out the design optimization as well. The numerical example of a small-scaled mechanical system is presented to verify the proposed formulation.

Hou, J. W.↗
NASA NTRS1973

Advanced control concepts

The problems of excess control devices and insufficient trim control capability on shuttle ascent vehicles were investigated. The trim problem is solved at all time points of interest using Lagrangian multipliers and a Simplex based iterative algorithm developed as a result of the study. This algorithm has the capability to solve any bounded linear problem with physically realizable constraints, and to minimize any piecewise differentiable cost function. Both solution methods also automatically distribute the command torques to the control devices. It is shown that trim requirements are unrealizable if only the orbiter engines and the aerodynamic surfaces are used.

Sharp, J. B.↗
NASA NTRS1987

A static data flow simulation study at Ames Research Center

Demands in computational power, particularly in the area of computational fluid dynamics (CFD), led NASA Ames Research Center to study advanced computer architectures. One architecture being studied is the static data flow architecture based on research done by Jack B. Dennis at MIT. To improve understanding of this architecture, a static data flow simulator, written in Pascal, has been implemented for use on a Cray X-MP/48. A matrix multiply and a two-dimensional fast Fourier transform (FFT), two algorithms used in CFD work at Ames, have been run on the simulator. Execution times can vary by a factor of more than 2 depending on the partitioning method used to assign instructions to processing elements. Service time for matching tokens has proved to be a major bottleneck. Loop control and array address calculation overhead can double the execution time. The best sustained MFLOPS rates were less than 50% of the maximum capability of the machine.

Barszcz, Eric↗
DOE OSTI2023

A Novel Partitioned Approach for Reduced Order Model—Finite Element Model (ROM-FEM) and ROM-ROM Coupling

Partitioned methods allow one to build a simulation capability for coupled problems by reusing existing single-component codes. In so doing, partitioned methods can shorten code development and validation times for multiphysics and multiscale applications. In this work, we consider a scenario in which one or more of the “codes” being coupled are projection-based reduced order models (ROMs), introduced to lower the computational cost associated with a particular component. We simulate this scenario by considering a model interface problem that is discretized independently on two non-overlapping subdomains. Here we then formulate a partitioned scheme for this problem that allows the coupling between a ROM “code” for one of the subdomains with a finite element model (FEM) or ROM “code” for the other subdomain. The ROM “codes” are constructed by performing proper orthogonal decomposition (POD) on a snapshot ensemble to obtain a low-dimensional reduced order basis, followed by a Galerkin projection onto this basis. The ROM and/or FEM “codes” on each subdomain are then coupled using a Lagrange multiplier representing the interface flux. To partition the resulting monolithic problem, we first eliminate the flux through a dual Schur complement. Application of an explicit time integration scheme to the transformed monolithic problem decouples the subdomain equations, allowing their independent solution for the next time step. We show numerical results that demonstrate the proposed method’s efficacy in achieving both ROM-FEM and ROM-ROM coupling.

de Castro, Amy↗
DOE OSTI2023

Natural Abundance Isotope Ratio Measurements of Organic Molecules Using 21 T FTICR MS

Subtle variations in stable isotope ratios at natural abundance are challenging to measure but can yield critical insights into biological, physical, and geochemical processes. Well established methods, particularly multi-collector, gas-source or plasma isotope ratio mass spectrometry (IRMS) are the gold standard for stable isotope measurement but inherent limitations in these approaches make them ill-suited to determining site-specific and multiply substituted isotopic abundances of all but a few compounds, or to characterizing mixtures or larger intact molecules. Additionally, Fourier transform mass spectrometry (FTMS), namely Orbitrap mass spectrometry, has recently demonstrated the ability to measure natural abundance isotope ratios with chemically informative accuracy and precision. Here, we report the first use of Fourier transform ion cyclotron resonance mass spectrometry (FTICR MS) for the accurate (<1‰) and precise (<1‰ standard error) simultaneous determination of $δ$ 13 C and $δ$ 15 N in caffeine isotopologues and provide a discussion of the critical instrumental parameters necessary to make such measurements. Finally, we report the ability to make these measurements with online liquid chromatography, expanding the ability of this technique to explore mixtures in the future.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

An analysis of the hybrid internal mass modeling approach in EnergyPlus

Accurate simulation of building system dynamics is particularly important for understanding building energy flexibility. Among all dynamics in a building, a zone temperature’s variation is especially important, as it significantly affects a building’s electricity load profile when its heating, ventilating, and air conditioning (HVAC) system is controlled with on/off cycles or setpoint reset strategies. To accurately simulate a zone temperature’s dynamics, internal mass needs to be modeled carefully. In this paper, we compare the two internal mass modeling approaches provided by EnergyPlus, i.e., internal mass object and zone air capacitance multiplier, to better understand their impacts on zone temperature simulation. Real building zone temperature dynamic data from small- and medium-sized office buildings are analyzed and compared with simulated data. In particular, we illustrate the effectiveness of a hybrid method of the two EnergyPlus modeling approaches, which yields more realistic zone temperature dynamics, especially when the zone is conditioned with heat pump systems with on/off cycling.

Chen, Zhelun↗
NASA NTRS2008

Shaping Diffraction-Grating Grooves to Optimize Efficiency

A method of shaping diffraction-grating grooves to optimize the spectral efficiency, spectral range, and image quality of a spectral imaging instrument is under development. The method is based on the use of an advanced design algorithm to determine the possibly complex shape of grooves needed to obtain a desired efficiency-versus-wavelength response (see figure). Then electron- beam fabrication techniques are used to realize the required groove shape. The method could be used, for example, to make the spectral efficiency of the grating in a given wavelength range proportional to the inverse of the spectral efficiency of a photodetector array so that the overall spectral efficiency of the combination of the grating and the photodetector array would be flat. The method has thus far been applied to one-dimensional gratings only, but in principle, it is also applicable to two-dimensional gratings. The algorithm involves calculations in the spatial-frequency domain. The spatial-frequency spectrum of a grating is represented as a diffraction-order spectral-peak-width function multiplied by an efficiency function for a single grating groove. This representation affords computational efficiency and accuracy by making it possible to consider only the response from one grating groove (one period of the grating), instead of from the whole grating area, in determining the response from the entire grating. This combination of efficiency and accuracy is crucial for future extensions of the algorithm to two-dimensional designs and to designs in which polarization must also be taken into account. The algorithm begins with the definition of target values of relative efficiency that represent the desired spectral response of the grating in certain spectral frequencies calculated from the diffraction order and wavelength. The grating period is divided into a number of cells - typically, 100. The phase contribution from each cell is determined from the phase of the incident electromagnetic wave and the height of the grating surface in the cell. The total contribution from all cells to each target value is then calculated. Then a method known to specialists as the optimum-rotation-angle method is used to adjust the height of each cell so that the total response from all cells is optimized. The computation is iterative and continues until the desired response is obtained. In the event that the desired response is unphysical, the algorithm nevertheless strives to generate a grating-grove profile for which the response approximates the desired one as closely as possible.

Backlund, John↗
NASA NTRS2018

Discovery of Activities via Statistical Clustering of Fixation Patterns

Human behavior often consists of a series of distinct activities, each characterized by a unique pattern of interaction with the visual environment. This is true even in a restricted domain, such as a piloting an aircraft, where activities with distinct visual signatures might be things like communicating, navigating, and monitoring. We propose a novel analysis method for gaze-tracking data, to perform blind discovery of these hypothetical activities. The method is in some respects similar to recurrence analysis, but here we compare not individual fixations, but groups of fixations aggregated over a fixed time interval. The duration of this interval is a parameter that we will refer to as delta. We assume that the environment has been divided into a set of N different areas-of-interest (AOIs). For a given interval of time of duration delta, we compute the proportion of time spent fixating each AOI, resulting in an N-dimensional vector. These proportions can be converted to integer counts by multiplying by delta divided by the average fixation duration (another parameter that we fix at 280 milliseconds). We compare different intervals by computing the chi-square statistic. The p-value associated with the statistic is the likelihood of observing the data under the hypothesis that the data in the two intervals were generated by a single process with a single set of probabilities governing the fixation of each AOI. The method has been applied to approximately 100 hours of eye movement data collected from pilots in a high-fidelity B747 flight simulator, and the results have been compared to synthetic data in which the each activity is represented as first-order Markov process with random probabilities assigned to the AOIs. Randomly-generated synthetic activities can require thousands of fixations to be discriminated with statistical significance, while the human data can be clustered using averaging windows of some 10's of seconds, suggesting that the actual activities are much more narrowly focused than random Markov models.

activity analysis↗
DOE OSTI2025

A Measurement of the Neutron Electromagnetic Form Factor Ratio from a Rosenbluth Technique with Simultaneous Detection of Neutrons and Protons

The internal structure of protons and neutrons provides insight into both the dynamical behavior of the constitute quarks and gluons, and emergent properties of the nucleons (such as mass, spin, and electromagnetic distributions). Elastic electron-nucleon scattering can probe the elastic electromagnetic form factors of the nucleon. The electric and magnetic form factors, respectively, encode information about the internal charge and magnetization distributions within the nucleon. Precision data for these form factors, over a broad range of the four-momentum transfer squared, Q^2, can benchmark theoretical models describing the strong interaction of nuclear physics. The Super BigBite Spectrometer (SBS) program in Hall A at Jefferson Lab, is a series of high-precision experiments which seek to significantly extend the Q^2 reach of previous data for the nucleon electromagnetic form factors. The first two experiments of this program are known as G_M^n and the neutron Two Photon Exchange (nTPE) and the data were collected from October 2021 to February 2022. Both experiments were conducted with the simultaneous measurement of D(e,e'n) and D(e,e'p) reactions for quasi-elastic electron-deuteron scattering. The scattered electrons were detected in the BigBite Spectrometer, which features multiple large-acceptance Gas Electron Multiplier (GEM) detectors. The Super BigBite Spectrometer provided simultaneous detection of scattered nucleons, and utilized a large acceptance dipole magnet and Hadron Calorimeter (HCal). The G_M^n experiment provides precision measurements of the neutron magnetic form factor, via the ratio method, over a Q^2 range of 3.0 to 13.5 (GeV/c)2. From this data analysis, preliminary values for G_M^n/µ_n G_D are extracted. For Q^2=4.48 (GeV/c)2 we find G_M^n/µ_n G_D=0.9546±0.0132 and for Q^2=4.476 (GeV/c)2 we find G_M^n/µ_n G_D=0.9563±0.0110. These preliminary G_M^n/µ_n G_D values are more precise than existing world data in this Q^2 regime and are consistent with the most recent parameterization of the G_M^n/µ_n G_D world data. The nTPE experiment provides a first measurement of the neutron Rosenbluth Slope and seeks to quantify the two-photon exchange(TPE) contribution to elastic electron-neutron scattering at a fixed Q^2=4.5 (GeV/c)2 with two different beam energies and scattering angle values. For data of the proton form factor ratio, µ_p G_E^p/G_M^p, significant discrepancies exist between values obtained from Rosenbluth Separation and polarization transfer measurement, particularly at large Q^2, and TPE contributions are thought to resolve this discrepancy. The impacts of TPE contributions have not yet been experimentally established for the neutron. From the data analysis presented in this dissertation, a preliminary result for the neutron Rosenbluth Slope is found as S^n=(G_E^n )^2/t_n (G_M^n )^2=0.0916±0.0476 for Q^2=4.48 (GeV/c)2. This value of the neutron Rosenbluth Slope is consistent with the world data extrapolation and the absence of large TPE corrections.

Wertz, Ezekiel [Thomas Jefferson National Accelera↗
DOE OSTI2023

Self-Powered Autonomous Sensing System for Arctic Ocean using a Frequency-multiplied Cylindrical Triboelectric Nanogenerator

An autonomous sensing system for collecting environmental observations in the Arctic region is critical for the estimation and prediction of climate change. Ocean waves are a great source of energy for these sensing systems but there has been limited research done on small-scale energy harvesting applications in the Arctic Ocean (subsea or surface). The available wave energy in the Arctic Ocean is lower than the typical ocean wave energy due to the low wave frequency, height, and operating months. Although the available wave energy depends on the specific location in the Arctic Ocean, the wave height decreases everywhere during the winter (Jan-March) [1]. Our target location in this work is the Beaufort and Chukchi Seas, which are ice covered for about 195 days per year leaving only a few months (June-November) for wave energy harvesting [2]. During these few months, the average wave frequency is 0.2 Hz and the most common wave frequency is around 0.15 Hz. Many energy harvesting methods are unsuitable for use in the Arctic Ocean because of to the cold temperature, low wave frequency, and low wave height. Triboelectric nanogenerators (TENG) are one of the few energy harvesting methods that excel in these conditions. Zhong et al. [3] designed a stacked pendulum-structured TENG for low-frequency ocean wave energy and reported the device generated a peak power density of 11.2 W/m3 under the low wave frequency of 0.2 Hz. The autonomous sensing system we developed for the Arctic Ocean, the Arctic-TENG, is based on a frequency-multiplied cylindrical TENG (FMC-TENG) which is an optimized TENG configuration for Arctic Ocean conditions due to the high power density under low-frequency wave conditions [4]. Figure 1 shows an FMC-TENG with multiple pairs of free-standing triboelectric-layer mode materials (Aluminum and FEP). The mass and magnet attached to the rotor store gravitational potential energy which is released as kinetic energy when the potential energy overcomes the repulsive magnetic force generated by the opposing magnet attached to the stator. This force unbalance triggers a sudden rotation and swinging motion of the mass which increases the angular velocity of the system and therefore enhances the output power of the TENG system. The main components of the Arctic-TENG are a 3D-printed rotor and stator, electric and dielectric material adhesive tape, and bearings. All these materials were cold-soaked and tested at -40 °C in a chest freezer. Both the 3D-printed parts and adhesives were confirmed to have a minimal effect from the cold temperature. Multiple bearings were tested and the starting torque of each one was compared both at room temperature and -40 °C. The bearing with the lowest starting torque at -40 °C was selected for use in the Arctic-TENG. The Arctic-TENG was tested using an out-of-water motor-driven wave simulator (Figure 2). The wave simulator allows for controlled testing at a wave height of 0.2 m and frequencies between 0.1 Hz to 0.5 Hz. This wave simulator was used for both room temperature and -40 °C testing. The Arctic-TENG generated significantly more power at -40 °C compared to room temperature for each frequency tested. The system stored energy in a supercapacitor via a power management circuit, and the amount of energy stored per day was calculated at different wave frequencies. Based on the conditions at the proposed deployment location (days of non-ice-covered ocean and wave frequency), the amount of stored energy per year from Arctic-TENG was calculated to be enough energy for two transmissions every day. Durability testing was completed on the Arctic-TENG at room temperature to determine the lifetime. The rotor of the Arctic-TENG was attached to a DC motor and spun continuously for several millions of cycles without degradation of the electrical output, demonstrating the feasibility for long-term operation.

Jung, Hyunjun↗
NASA NTRS2019

Discovery of Activities via Statistical Clustering of Fixation Patterns

Human behavior often consists of a series of distinct activities, each characterized by a unique signature of visual behavior. This is true even in a restricted domain, such as piloting an aircraft, where patterns of visual signatures might represent activities like communicating, navigating, and monitoring. We propose a novel analysis method for gaze-tracking data, to perform blind discovery of these activities based on their behavioral signatures. The method is in some respects similar to recurrence analysis, but here we compare not individual fixations, but groups of fixations aggregated over a fixed time interval. The duration of this interval is a parameter that we will refer to as τ. We assume that the environment has been divided into a set of N different areas-of-interest (AOIs). For a given interval of time of duration τ, we compute the proportion of time spent fixating each AOI, resulting in an N-dimensional vector. These proportions can be converted to counts by multiplying by τ divided by the average fixation duration (another parameter that we fix at 280 milliseconds). We compare different intervals by computing the chi-square statistic. The p-value associated with the statistic is the likelihood of observing the data under the hypothesis that the data in the two intervals were generated by a single process with a single set of probabilities governing the fixation of each AOI. We have investigated the method using a set of 10 synthetic "activities," that sample 4 AOIs. Four of these activities visit 3 of the 4 AOIs, with equal probability; as there are four different ways to leave-one- out, there are four such activities. Similarly, there are six different activities that leave-two-out. Sequences of simulated behavior were generated by running each activity for 40 seconds, in sequence, for a total of 6.7 minutes. The figure to the right shows the matrix of chi-square statistics, using a value of 2.8 seconds for τ, corresponding to 10 fixations. Low values (dark) indicate poor evidence for activity differences, while high values (bright) indicate strong evidence. The dark squares along the main diagonal each correspond to the forty second intervals in which the activity was held constant; the 4x4 block at the lower left corresponds to the four leave-one-out activities, while the 6x6 block in the upper right corresponds to the leave-two-out activities. (The anti-diagonal pattern of white squares indicates those activity pairs that share no AOIs.) The chi-square values can be binarized by choosing a particular significance level; we are interested in grouping bins that represent the same activity, effectively accepting the null hypothesis. Therefore, we may adopt a relatively lax criterion; for example, choosing a p-value of 0.2 means that two behaviors that have only a 1-in-5 chance of being produced by a single activity might nevertheless be clustered together. We have explored several methods to perform clustering on the data and solving for the activity probabilities. Greedy methods begin by selecting the time bin that is similar to the most (or least) other bins, and then forming a cluster from it and all other non-discriminable bins. These methods show mediocre performance, as they do not take into account temporal contiguity. Preliminary results indicate that methods that "grow" clusters in time from seed points perform better.

activity analysis↗
NASA NTRS1993

A Fractional-Step Method Of Computing Incompressible Flow

Method of computing time-dependent flow of incompressible, viscous fluid involves numerical solution of Navier-Stokes equations on two- or three-dimensional computational grid based on generalized curvilinear coordinates. Equations of method derived in primitive-variable formulation. Dependent variables are pressure at center of each cell of computational grid and volume fluxes across faces of each cell. Volume fluxes replace Cartesian components of velocity; these fluxes correspond to contravariant components of velocity multiplied by volume of computational cell, in staggered grid. Choice of dependent variables enables simple extension of previously developed staggered-grid approach to generalized curvilinear coordinates and facilitates enforcement of conservation of mass.

Kwak, Dochan↗
NASA NTRS1953

Effect of variable viscosity and thermal conductivity on high-speed slip flow between concentric cylinders

The differential equations of slip flow, including the Burnett terms, were first solved by Schamberg assuming that the coefficients of viscosity and heat conduction of the gas were constants. The problem is solved herein for variable coefficients of viscosity and thermal conductivity by applying a transformation leading to an iteration method. The method, starting with the solution for constant coefficients, enables one to approximate the solution for variable coefficients very closely after one or two steps. Satisfactory results are shown to follow from Schamberg's solution by using his values of the constant coefficients multiplied by a constant factor 'N', leading to what are denoted as the effective coefficients of viscosity and thermal conductivity.

RESEARCH TECHNIQUE, AERODYNAMICS↗
DOE OSTI2023

Construction and Commissioning of Gas Electron Multiplier (GEM) Detectors in Advanced Assembly Design for Low-Energy Applications at High Rates and Analysis of GEM Data From the Muse Experiment at PSI

The search beyond the Standard Model explores dark matter and a potential fifth force. DarkLight@ARIEL investigates the dark photon as a mediator between ordinary and dark matter. The experiment aims to measure the process e?Ta ? e?TaX ? e?Ta(e?e+), by detecting a charged lepton pair in the final state. The spectrometers will be instrumented with Gas Electron Multiplier (GEM) detectors with minimal material budget for tracking. A novel GEM construction technique is employed for fabricating 25 cm x 40 cm GEM detectors, where all layers are mechanically stretched and assembled within a double frame. The dissertation outlines physics motivation, methods, experimental setup, and the role of GEM detectors in DarkLight. A comprehensive discussion on GEM detectors, focusing on the "NS2" technique, including design, parts, assembly, testing, and performance evaluation will be discussed. The proton has garnered attention due to discrepancies in measuring its charge radius using muonic hydrogen and electron-based methods. The 2010 muonic hydrogen measurement of the proton charge radius Rp = 0.84184(67) fm, showed a significant 7? discrepancy compared to the previously known value of Rp = 0.8775(51) fm, giving rise to the proton radius puzzle. Through simultaneous measurements of e ? p and ? ? p elastic scattering, MUSE facilitates a precise and direct comparison of the proton radius. The physics background of the proton radius puzzle, measurement techniques, the approach adopted by the MUSE experiment, its significance and the experimental setup will be discussed. Furthermore, the role of GEM detectors, along with an in-depth analysis of their data and efficiency will be discussed and evaluated.

Mohammed Prem Nazeer, Sahara↗
DOE OSTI2022

Implement and Test 3D Mortar Contact in BISON

We leverage the extension of the generation of mortar segment meshes to three dimensions in MOOSE’s framework to extend thermomechanical modeling capabilities to problems with three dimensions. A modular approach to gap heat transfer physics using the mortar finite element method was created and documented, mechanical contact was extended to three dimensions—including frictional behavior, performance and ease of use were improved, and steps towards scalability of solid mechanics problems involving contact were taken. Many of these new developments are demonstrated in the simulation of 3D light-water reactor (LWR) problems, where the thermomechanical interface problem is solved using the mortar finite element method. Usage of the mortar framework has improved convergence in 2D problems and has enabled employing friction in 3D problems, of which we show results of a short, local stack of 3D pellets. Consequently, the benefits of mortar in terms of solution convergence and quality are extended to three dimensions. Section 2 discusses fundamental developments that enabled the simulation of practical mortar problems in three dimensions and other general improvements, including the reduction of the derivative container size, the modification of dual basis computations when edge dropping (lack of secondary element projection) takes place, the improvement of conditioning when employing the VCP in-edge dropping conditions, and code usability and quality improvements. These latter code enhancements include the migration of tests using “old” mortar contact constraints to using dual mortar with a semi-smooth Newton solution strategy and the reuse of lower dimensional domains for straightforwardly setting up a mortar thermomechanical LWR problem, i.e. the MOOSE action is employed for mechanical contact and the thermal LWR action is employed to capture the gas conductance, contact, and radiation components of gap heat transfer physics. Independently of the mortar LWR thermal action, we developed a modular approach to gap heat transfer that resides in MOOSE and can be leveraged, e.g., in metallic fuel problems. This approach, whose code design based on MOOSE’s user objects to model specific physics was proposed by the maintenance activity, is detailed in Section 3. Based on the dual mortar finite element method, the frictional contact constraints were extended to three dimensions. A block sheared in two directions in and out of contact with a rigid plane is employed in Section 4 to show the way the approach handles changes in frictional states (e.g. stick to slip) within a competitive number of Newton iterations. Equations and numerical results on the use of the VCP with Cartesian Lagrange multipliers, whose combination enables their direct condensation, are described in Section 5.3. Two-dimensional and three-dimensional BISON LWR simulations are discussed in Section 6. Particularly, a stack of five eccentric pellets with a surface defect is simulated and the effect of pellet-cladding friction is assessed. Finally, conclusions are outlined in Section 7.

42 ENGINEERING↗
NASA NTRS2010

Radio-Frequency Tank Eigenmode Sensor for Propellant Quantity Gauging

Although there are several methods for determining liquid level in a tank, there are no proven methods to quickly gauge the amount of propellant in a tank while it is in low gravity or under low-settling thrust conditions where propellant sloshing is an issue. Having the ability to quickly and accurately gauge propellant tanks in low-gravity is an enabling technology that would allow a spacecraft crew or mission control to always know the amount of propellant onboard, thus increasing the chances for a successful mission. The Radio Frequency Mass Gauge (RFMG) technique measures the electromagnetic eigenmodes, or natural resonant frequencies, of a tank containing a dielectric fluid. The essential hardware components consist of an RF network analyzer that measures the reflected power from an antenna probe mounted internal to the tank. At a resonant frequency, there is a drop in the reflected power, and these inverted peaks in the reflected power spectrum are identified as the tank eigenmode frequencies using a peak-detection software algorithm. This information is passed to a pattern-matching algorithm, which compares the measured eigenmode frequencies with a database of simulated eigenmode frequencies at various fill levels. A best match between the simulated and measured frequency values occurs at some fill level, which is then reported as the gauged fill level. The database of simulated eigenmode frequencies is created by using RF simulation software to calculate the tank eigenmodes at various fill levels. The input to the simulations consists of a fairly high-fidelity tank model with proper dimensions and including internal tank hardware, the dielectric properties of the fluid, and a defined liquid/vapor interface. Because of small discrepancies between the model and actual hardware, the measured empty tank spectra and simulations are used to create a set of correction factors for each mode (typically in the range of 0.999 1.001), which effectively accounts for the small discrepancies. These correction factors are multiplied to the modes at all fill levels. By comparing several measured modes with the simulations, it is possible to accurately gauge the amount of propellant in the tank. An advantage of the RFMG approach of applying computer simulations and a pattern-matching algorithm is that the Although there are several methods for determining liquid level in a tank, there are no proven methods to quickly gauge the amount of propellant in a tank while it is in low gravity or under low-settling thrust conditions where propellant sloshing is an issue. Having the ability to quickly and accurately gauge propellant tanks in low-gravity is an enabling technology that would allow a spacecraft crew or mission control to always know the amount of propellant onboard, thus increasing the chances for a successful mission. The Radio Frequency Mass Gauge (RFMG) technique measures the electromagnetic eigenmodes, or natural resonant frequencies, of a tank containing a dielectric fluid. The essential hardware components consist of an RF network analyzer that measures the reflected power from an antenna probe mounted internal to the tank. At a resonant frequency, there is a drop in the reflected power, and these inverted peaks in the reflected power spectrum are identified as the tank eigenmode frequencies using a peak-detection software algorithm. This information is passed to a pattern-matching algorithm, which compares the measured eigenmode frequencies with a database of simulated eigenmode frequencies at various fill levels. A best match between the simulated and measured frequency values occurs at some fill level, which is then reported as the gauged fill level. The database of simulated eigenmode frequencies is created by using RF simulation software to calculate the tank eigenmodes at various fill levels. The input to the simulations consists of a fairly high-fidelity tank model with proper dimensions and including internal tank harare, the dielectric properties of the fluid, and a defined liquid/vapor interface. Because of small discrepancies between the model and actual hardware, the measured empty tank spectra and simulations are used to create a set of correction factors for each mode (typically in the range of 0.999 1.001), which effectively accounts for the small discrepancies. These correction factors are multiplied to the modes at all fill levels. By comparing several measured modes with the simulations, it is possible to accurately gauge the amount of propellant in the tank. An advantage of the RFMG approach of applying computer simulations and a pattern-matching algorithm is that the

Zimmerli, Gregory A.↗
NASA NTRS1984

Evolution of assumed stress hybrid finite element

Early versions of the assumed stress hybrid finite elements were based on the a priori satisifaction of stress equilibrium conditions. In the new version such conditions are relaxed but are introduced through additional internal displacement functions as Lagrange multipliers. A rational procedure is to choose the displacement terms such that the resulting strains are now of complete polynomials up to the same degree as that of the assumed stresses. Several example problems indicate that optimal element properties are resulted by this method.

Pian, T. H. H.↗