Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equation recovery”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Global Mass Flux Solutions from GRACE: A Comparison of Parameter Estimation Strategies - Mass Concentrations Versus Stokes Coefficients

The differences between mass concentration (mas con) parameters and standard Stokes coefficient parameters in the recovery of gravity infonnation from gravity recovery and climate experiment (GRACE) intersatellite K-band range rate data are investigated. First, mascons are decomposed into their Stokes coefficient representations to gauge the range of solutions available using each of the two types of parameters. Next, a direct comparison is made between two time series of unconstrained gravity solutions, one based on a set of global equal area mascon parameters (equivalent to 4deg x 4deg at the equator), and the other based on standard Stokes coefficients with each time series using the same fundamental processing of the GRACE tracking data. It is shown that in unconstrained solutions, the type of gravity parameter being estimated does not qualitatively affect the estimated gravity field. It is also shown that many of the differences in mass flux derivations from GRACE gravity solutions arise from the type of smoothing being used and that the type of smoothing that can be embedded in mas con solutions has distinct advantages over postsolution smoothing. Finally, a 1 year time series based on global 2deg equal area mascons estimated every 10 days is presented.

Rowlands, D. D.↗
DOE OSTI2024

An integrated approach to optimizing concentration shock wave electrodialysis using 2D multicell simulation and response surface models

Shock wave electrodialysis (SWED) is a highly promising technique for energy-efficient ion separation in the context of a circular economy. This paper presents a approach way of modeling and improving SWED using a two-dimensional multicell model combined with the COMSOL program and response surface methodology. The model integrates the Nernst-Planck equation, Darcy's law, and first-order electroosmosis to examine the local concentration, flux of ionic species, distribution of current, and velocity of flow in SWED cells under various operating conditions. We first illustrate the clear depiction of concentration, velocity, and electric potential distribution through contours which aids in identifying optimal operating conditions and designing scalable SWED systems. The results emphasize the significance of surface charge density and voltage in influencing the features of shock waves for obtaining effective ion separation while optimizing energy consumption and improving current efficiency by controlling the retention time of feed flow. Here, this study defines two crucial characteristics of shock waves, namely the length of the flat depletion zone of a fully developed shock wave (shock wave height) and the distance of shock wave propagation (shock wave length). These properties significantly impact separation performance, as determined by the simulation results. Additionally, the response surface methodology is incorporated with the COMSOL models to develop predictive models and graph responses, enabling a more comprehensive understanding of the interactions between parameters and performance indicators, such as removal ratio, energy consumption, and water recovery. Finally, this work suggests design tactics for expanding SWED processes and outlines potential areas for further research. This research provides valuable insights into the prospective applications, design optimization, and scalability of SWED in the field of electrokinetic separation technologies for green chemistry and a circular economy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Lyapunov-Based Sensor Failure Detection And Recovery For The Reverse Water Gas Shift Process

Livingstone, a model-based AI software system, is planned for use in the autonomous fault diagnosis, reconfiguration, and control of the oxygen-producing reverse water gas shift (RWGS) process test-bed located in the Applied Chemistry Laboratory at KSC. In this report the RWGS process is first briefly described and an overview of Livingstone is given. Next, a Lyapunov-based approach for detecting and recovering from sensor failures, differing significantly from that used by Livingstone, is presented. In this new method, models used are in terms of the defining differential equations of system components, thus differing from the qualitative, static models used by Livingstone. An easily computed scalar inequality constraint, expressed in terms of sensed system variables, is used to determine the existence of sensor failures. In the event of sensor failure, an observer/estimator is used for determining which sensors have failed. The theory underlying the new approach is developed. Finally, a recommendation is made to use the Lyapunov-based approach to complement the capability of Livingstone and to use this combination in the RWGS process.

Haralambous, Michael G.↗
NASA NTRS2002

LYAPUNOV-Based Sensor Failure Detection and Recovery for the Reverse Water Gas Shift Process

Livingstone, a model-based AI software system, is planned for use in the autonomous fault diagnosis, reconfiguration, and control of the oxygen-producing reverse water gas shift (RWGS) process test-bed located in the Applied Chemistry Laboratory at KSC. In this report the RWGS process is first briefly described and an overview of Livingstone is given. Next, a Lyapunov-based approach for detecting and recovering from sensor failures, differing significantly from that used by Livingstone, is presented. In this new method, models used are in t e m of the defining differential equations of system components, thus differing from the qualitative, static models used by Livingstone. An easily computed scalar inequality constraint, expressed in terms of sensed system variables, is used to determine the existence of sensor failures. In the event of sensor failure, an observer/estimator is used for determining which sensors have failed. The theory underlying the new approach is developed. Finally, a recommendation is made to use the Lyapunov-based approach to complement the capability of Livingstone and to use this combination in the RWGS process.

Haralambous, Michael G.↗
NASA NTRS1988

Solar-cycle modulation of galactic cosmic rays

This paper describes a numerical solution of the spherically-symmetric, time-dependent, cosmic-ray modulation equation, compares it to data taken near the ecliptic, and probes the physical mechanisms responsible for the solar-cycle variations in the cosmic-ray intensity. With a physically defensible diffusion coefficient and other reasonable parameters, the spectra of protons and electrons can be simulated simultaneously over the cycle, the radial intensity gradients for particles greater than about 100 MeV, the overall intensity variation during the cycle, and the time lag in the recovery of low-energy particles behind high-energy ones, known commonly as the 'hysteresis.' The results suggest that cosmic-ray variations near the ecliptic are dominated by turbulent scattering regions in the heliosphere, overwhelming any effects of gradient and curvature drifts.

Perko, John S.↗
NASA NTRS2006

Self-Consistent Ring Current Modeling with Propagating Electromagnetic Ion Cyclotron Waves in the Presence of Heavy Ions

The self-consistent treatment of the RC ion dynamics and EMIC waves, which are thought to exert important influences on the ion dynamical evolution, is an important missing element in our understanding of the storm-and recovery-time ring current evolution. Under certain conditions, relativistic electrons, with energies 21 MeV, can be removed from the outer radiation belt by EMIC wave scattering during a magnetic storm. That is why the modeling of EMIC waves is critical and timely issue in magnetospheric physics. To describe the RC evolution itself this study uses the ring current-atmosphere interaction model (RAM). RAM solves the gyration and bounce-averaged Boltzmann-Landau equation inside of geosynchronous orbit. Originally developed at the University of Michigan, there are now several branches of this model currently in use as describe by Liemohn namely those at NASA Goddard Space Flight Center This study will generalize the self-consistent theoretical description of RC ions and EMIC waves that has been developed by Khazanov and include the heavy ions and propagation effects of EMIC waves in the global dynamic of self-consistent RC - EMIC waves coupling. The results of our newly developed model that will be presented at GEM meeting, focusing mainly on the dynamic of EMIC waves and comparison of these results with the previous global RC modeling studies devoted to EMIC waves formation. We also discuss RC ion precipitations and wave induced thermal electron fluxes into the ionosphere.

Khazanov, George V.↗
NASA NTRS1982

Dynamical features of moving double layers

Numerical simulations of the dynamics of double layers that form in response to an applied potential drop across a plasma are carried out for plasmas of lengths much greater than previously considered. The evolution of finite-size plasmas of lengths 100, 200 and 400 Debye lengths at the low-potential end of the simulation region was followed by the solution of the Vlasov and Poisson equations. Results show that the double layer moves towards the high-potential side of the layer, at a speed directly dependent on the length of the plasma. The moving double layer is accompanied by a moving density front which is similar to an ion acoustic shock created by the expansion of a high-density plasma into a low-density plasma, as well as plasma heating and evacuation, electron and ion current interruptions and recovery. Double layer formation, motion, and accompanying phenomena are found to repeat periodically. Results are in good agreement with laboratory experiments, and may be used to explain certain phenomena observed in auroral plasmas.

Singh, N.↗
DOE OSTI2023

Hydrothermal but Not Mechanical Pretreatment of Wastewater Algae Enhanced Anaerobic Digestion Energy Balance due to Improved Biomass Disintegration and Methane Production Kinetics

This study used pilot-scale high-rate algae ponds to assess algal–bacteria biomass productivity and wastewater nutrient removal as well as the impact of mechanical and hydrothermal pretreatments on biomass disintegration, methane production kinetics, and anaerobic digestion (AD) energy balance. Mechanical pretreatment had a minor effect on biomass disintegration and methane production. By contrast, hydrothermal pretreatment significantly reduced particle size and increased the solubilized organic matter content by 3.5 times. The methane yield and production rate increased by 20–55% and 20–85%, respectively, with the highest values achieved after pretreatment at 121 °C for 60 min. While the 1st-order and pseudo-1st-order reaction equation models fitted methane production from untreated biomass best (R2 > 0.993), the modified Gompertz sigmoidal-type model provided a superior fit for hydrothermally pretreated algae (R2 ≥ 0.99). The AD energy balance revealed that hydrothermal pretreatment improved the total energy output by 25–40%, with the highest values for volume-specific and mass-specific total energy outputs reaching 0.23 kW per digester m3 and 2.3 MW per ton of biomass volatile solids. Additionally, net energy recovery (energy output per biomass HHV) increased from 20% for untreated algae to 32–34% for hydrothermally pretreated algae, resulting in net energy ratio and net energy efficiency of 2.14 and 68%, respectively.

09 BIOMASS FUELS↗
NASA NTRS2016

Computational Study of Flow Interactions in Coaxial Rotors

Although the first idea of coaxial rotors appeared more than 150 years ago, most helicopters have used single main-rotor/tail-rotor combination. Since reactive moments of coaxial rotors are canceled by contra-rotation, no tail rotor is required to counter the torque generated by the main rotor. Unlike the single main rotor design that distributes power to both main and tail rotors, all of the power for coaxial rotors is used for vertical thrust. Thus, no power is wasted for anti-torque or directional control. The saved power helps coaxial rotors reach a higher hover ceiling than single rotor helicopters. Another advantage of coaxial rotors is that the overall rotor diameter can be reduced for a given vehicle gross weight because each rotor provides a maximum contribution to vertical thrust to overcome vehicle weight. However, increased mechanical complexity of the hub has been one of the challenges for manufacturing coaxial rotorcraft. Only the Kamov Design Bureau of Russia had been notably successful in production of coaxial helicopters until Sikorsky built X2, an experimental compound helicopter. Recent developments in unmanned aircraft systems and high-speed rotorcraft have renewed interest in the coaxial configuration. Multi-rotors are frequently used for small electric unmanned rotorcraft partly due to mechanical simplicity. The use of multiple motors provides redundancy as well as cost-efficiency. The multi-rotor concept has rarely been used until recently because of its inherent stability and control problems. However, advances in inexpensive electronic flight control systems have opened the floodgates for small drones using multirotors. Coaxial rotors have started to appear in some multi-rotor configurations. Small coaxial rotors have often been designed using a hundred year old approach that is "sketch, build, fly, and iterate." In that approach, there is no systematic way to explore trade-offs or determine logical next steps. It is neither possible to account for multiple real-world constraints up front in design nor possible to know what performance is possible with a given design. Since unmanned vehicles are sized and optimized for the particular mission, a modern low-fidelity conceptual design and sizing tool that has been used for the design of large helicopters can be used for design of small coaxial rotorcraft. However, unlike most helicopters with single main rotor, the interactions between the upper and lower rotors emerge as an important factor to consider in design because an increase in performance of a multi-rotor system is not proportional to the number of rotors. Interference losses and differences in thrusts between the upper and lower rotors were investigated by theoretical methods as well as a computational fluid dynamics (CFD) method using the Reynolds-Averaged Navier-Stokes (RANS) equations. In this work, hybrid turbulence models are used to investigate the physics of interactions between coaxial rotors and a fuselage that are not well understood. Present study covers not only small-scale drones but also large-scale coaxial rotors for heavy-lifting missions. Considering the recently proposed FAA drone rules that require the flight only in visual line-of-sight, a large multirotor might be used as an airborne carrier for launch and recovery of unmanned aircraft systems with a human operator onboard. For applications to civil operations, their aerodynamic performance and noise levels need to be assessed. Noise is one of the largest limiting factors to rotorcraft operations in urban area. Since the high-frequency noise of multi-rotors may increase the annoyance, noise may turn out to be a key issue that must be addressed for market acceptability. One of the objectives of the present work is to study the effects of inter-rotor spacing and collectives on the performance, efficiency, and acoustics of coaxial rotor systems.

Study↗
DOE OSTI2024

HL-2A's ELM cycle simulations by integrating BOUT++'s drift MHD and transport code

A new integrating model has been developed to couple tokamak edge multiscale magnetohydrodynamic (MHD) events and transport simulations, such as edge-localized mode (ELM) cycles. As a proof of principle, we first start from a set of three-field two-fluid model equations, which includes the pressure, current, and vorticity. Here, the equations are separated into the slowly evolving part of the axisymmetric component by taking a time average of the axisymmetric component. The time-averaged fluxes, which are quadratic in fluctuating quantities, act as driven terms for the time-averaged axisymmetric quantities that determine the plasma transport, and therefore the large-scale evolution of the plasma profiles. Then the HL-2A's ELM cycles are simulated using the model. Good agreements of ELM size and pedestal recovery time have been achieved for the solutions obtained from the coupled simulation compared with experiment. For one ELM cycle simulation, the coupled code can achieve a speedup of a factor of up to 30 over standalone code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1973

The influence of laser ranging on selenodetic control.

A mathematical model for performing an adjustment, using laser distances to improve coordinates of points on the earth and on the moon, as well as the orientation of these two bodies in space, is presented. The observation equations are given, and the orientation of the earth and the moon are defined in terms of three Eulerian angles. Numerical experiments were performed to investigate the expected accuracies of parameters in the adjustment model under differing conditions. The results indicate that the statistics both for the coordinates of points on earth and on the moon, and for the orientation parameters of the two bodies are favorable. However, very high correlations exist between some of the parameters; therefore, the recovery of these parameters from a simultaneous adjustment is questionable.

Mueller, I. I.↗
DOE OSTI2024

Structure–Composition Relationships for Mg–Ni and Mg–Fe Olivine

Olivine is a dynamic and important mineral in the crust and mantle with relevance to processes important to climate change technology, such as geologic carbon storage and critical mineral recovery. In this work, we critically evaluated and compiled a new database of olivine diffraction data, lattice parameters, and composition to enable rapid Ni-Mg-Fe olivine composition determination. A compilation of olivine X-ray diffraction data and chemical compositions from both the literature and the International Centre for Diffraction Data (ICDD) powder database was assembled to plot both the forsterite-fayalite and forsterite-liebenbergite solid solution lines. Here we present an expanded dataset to delineate equations and relationships used for quantifying the correlations between olivine lattice parameters and chemical compositions in Mg 2 SiO 4 -Fe 2 SiO 4 (forsterite-fayalite) and Mg 2 SiO 4 -Ni 2 SiO 4 (forsterite-liebenbergite) olivine solid solution series.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1986

A perturbation approach to cosmic ray transients in interplanetary space

A perturbation approach is used to model the linear response of the cosmic ray distribution as a function of perturbations in the solar wind transport parameters. The analytical technique permits examination of the effects of the different solar wind parameters, i.e., velocity, drift velocity and the diffusion tensor, which are additive. Cosmic ray changes in the solar wind are neglected. Variations in the diffusion coefficient and the convection-diffusion equation are applied to describing Forbush decreases and the 11-yr solar cycle variations. The treatment is limited to energies above 100 MeV, yet is considered valid enough to have identified a hysteresis effect in the 11-yr variation, wherein high-energy particles lead low-energy particles to a magnitude that increases with increasing heliocentric distance within the inner heliosphere. The calculations also indicate that the Forbush decreases are a cumulative effect of precipitous precursors accompanied by slow recoveries from travelling solar wind perturbations.

Chih, P. P.↗
NASA NTRS1977

An approximate calculation of the strong interaction on a transonic airfoil

This paper describes an interactive computational scheme that calculates the transonic flow field over an airfoil. The inviscid flow, the turbulent boundary layer in the weak interaction region, and the elliptic flow in the strong interaction region are computed separately, and the solutions are matched iteratively in a direct, second-order interaction scheme. The pressure distribution, obtained from a transonic, small-perturbation solution, is prescribed on the boundary layer and the Navier-Stokes equations and the ensuing displacement thicknesses are used to modify the inviscid flow. The best overall results are obtained by using a turbulence model that eliminates the eddy viscosity in the strong interaction region. Computed pressure distributions are compared with experimental data and with Navier-Stokes solutions. In weak interactions, the agreement with experimental data is very good, as it is also in strong interactions except for the pressures in the separated-flow region. There the agreement is qualitative only, and the pressure recovery is too high. The results in both cases are superior or equal to those of the Navier-Stokes solutions, and the present method is significantly faster.

Seginer, A.↗
DOE OSTI2023

TOUGH3/EOS7MG (EOS7MG) v1.1

EOS7MG is a new TOUGH3 module for modeling water, brine, and multiple gas mixtures in the subsurface multiphase flow system. The module allows mixture of as more as three gas components selecting from a list of predefined gases (CH4, C2H6, C3H8, H2S, CO2, N2, NO2, H2, and AIR). EOS7MG was developed based on EOS7C and inherited the general modeling capabilities of the TOUGH3 platform. It uses cubic equation of states for gas density, viscosity, and enthalpy calculation. The module provides the selection options for Peng-Robinson, Redlich-Kwong, or Soave-Redlich-Kwong equations of state. The subsurface flow and transport of gas and aqueous phase mixtures over a wide range of pressures and temperatures is modeled with a multiphase Darcy's Law and an accurate chemical equilibrium approach as originally implemented in EOS7C is used for calculation of the solubilities of different mass components. EOS7MG was designed to allow flexible number of components based on the number of gases included in the model. The potential application of this module may include subsurface gas storage, energy storage, CCS, enhanced gas recovery, and more.

Zhang, Keni↗
NASA NTRS1972

Liquid-phase mixing of bipropellant doublets

Experimental results of unlike doublet mixing are correlated with an analytically derived equation predicting fluid cavitation. The correlation relates the minimum orifice pressure drop required to initiate cavitation, with the system back pressure, cold flow simulant vapor pressure, and the orifice flow discharge and contraction coefficients. Stream flow instabilities are also visually correlated with the onset of cavitation and orifice discharge coefficient measurements. The influence of cavitation on the characteristic phenomenon of hydraulic flip is observed for both circular and noncircular shaped orifices. For certain intermediate orifice lengths, some noncircular shapes are shown to produce more fully developed flows (shorter recovery lengths) and therefore a more cohesive jet, which in turn yields slightly higher cold flow mixing uniformities than circular shaped orifices of equal absolute length. The particular noncircular shaped elements evaluated are shown to be more sensitive to liquid stream misimpingement than the corresponding circular orifices.

Hoehn, F. W.↗
NASA NTRS1975

Behaviour of nitric oxide trails deposited in the mesosphere and stratosphere

The transient behavior of a nitric oxide trail deposited at approximately 60 km altitude is studied by the solution of the appropriate multidimensional diffusion equation which includes terms representing the effects of wind shear. Similar analysis is then carried out for the situation in the stratosphere. Trail behavior is found to be relatively independent of altitude and background ozone, but strongly dependent on the magnitude of eddy diffusity and the initial nitric oxide concentration. The nitric oxide trail reacts with ambient ozone to form nitrogen dioxide. For a trail 100 m initial radius, an ozone hole will form to a maximum size in 4 to 6 hours and then decay. The overall recovery time of the atmosphere following the creation of the trail is less than 12 hours.

Eberstein, I. J.↗
NASA NTRS1992

Derivation and application of hard deadlines for real-time control systems

The computation-time delay in the feedback controller of a real-time control system may cause failure to update the control input during one or more sampling periods. If this delay exceeds a certain limit called a hard deadline, either the necessary conditions for system stability are violated or the system leaves the allowed state-space. In such a case a dynamic failure is said to occur to the system. A method for calculating the hard deadlines in linear time-invariant control systems by considering system stability and the allowed state-space is presented. To derive necessary conditions for (asymptotic) system stability, the state difference equation is modified based on an assumed maximum delay and the probability distribution of delays whose magnitudes are less than, or equal to, the assumed maximum delay. Moreover, the allowed state-space - which is derived from given input and state constraints - is used to calculate the hard deadline as a function of time and the system state. A one-shot delay model in which a single event causes a dynamic failure is also considered. The knowledge of hard deadline is then applied to the design of error recovery in a triple modular redundant (TMR) controller computer.

Shin, Kang G.↗