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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Nonequilibrium Physics of Multiphase Flow in Porous Media: Wettability and Disorder

Carbon capture and geologic storage, dissociation of methane hydrates in permafrost, enhanced oil recovery, and water dropout in low-temperature fuel cells, all have something in common: two or more fluids flow simultaneously through a porous medium; and the displacement of one fluid by another is often unstable (either due to gravity or viscous forces). Yet, our ability to model multiphase flow in porous media has remained a challenge. The traditional equations are unable to predict, explain, or even reproduce, the formation of the complex patterns observed in experiments. The overarching goal of the project is to develop new physical understanding of the role of wettability and disorder in multiphase flow through permeable media, and to develop new mathematical and computational models at the pore scale and at the continuum scale.

58 GEOSCIENCES↗
NASA NTRS1977

Observations of ion cyclotron waves within the plasmasphere by Hawkeye 1

A survey of the plasma wave data from the Hawkeye 1 spacecraft has been performed in search of ion cyclotron waves associated with the scattering and loss of ring current ions within and near the plasmapause. During an 18-month period, encompassing about 270 orbits, a total of five events have been found with clearly detectable electric and magnetic fields at frequencies below the proton gyrofrequency. Comparisons of the electric and magnetic field amplitudes for these events provide strong evidence that these waves are ion cyclotron waves. All five events occurred during recovery phases of magnetic storms inside or very close to the plasmapause boundary. The results of this survey confirm and are consistent with the earlier identification of ion cyclotron waves by the Explorer 45 satellites. The Hawkeye 1 observations show that ion cyclotron waves of substantial amplitude occur at magnetic latitudes well away (about 28 deg) from the magnetic equator.

Kintner, P. M.↗
NASA NTRS1988

Vehicle for Space Transfer and Recovery (VSTAR), volume 2: Substantiating analyses and data

The Vehicle Space Transfer and Recovery (VSTAR) system is designed as a manned orbital transfer vehicle (MOTV) with the primary mission of Satellite Launch and Repair (SLR). Reference materials, calculations and trade studies used in the analysis and selection of VSTAR components. Each major VSTAR system is examined separately. Simple graphs and tables are used to make qualitative comparisons of various VSTAR component candidates. Equations and/or calculations used for a particular analysis are also included where applicable.

Source record↗
NASA NTRS1994

On the removal of boundary errors caused by Runge-Kutta integration of non-linear partial differential equations

It has been previously shown that the temporal integration of hyperbolic partial differential equations (PDE's) may, because of boundary conditions, lead to deterioration of accuracy of the solution. A procedure for removal of this error in the linear case has been established previously. In the present paper we consider hyperbolic (PDE's) (linear and non-linear) whose boundary treatment is done via the SAT-procedure. A methodology is present for recovery of the full order of accuracy, and has been applied to the case of a 4th order explicit finite difference scheme.

Abarbanel, Saul↗
NASA NTRS1983

A theory of rotating stall of multistage axial compressors. III - Limit cycles

A theory of rotating stall, based on single parameters for blade-passage lag and external-flow lag and a given compressor characteristic yields limit cycles in velocity space. These limit cycles are governed by Lienard's equation with the characteristic playing the role of nonlinear damping function. Cyclic integrals of the solution determine stall propagation speed and the effect of rotating stall on average performance. Solution with various line-segment characteristics and various throttle settings are found and discussed. There is generally a limiting flow coefficient beyond which no solution is possible; this probably represents stall recovery. This recovery point is independent of internal compressor lag, but does depend on external lags and on the height-to-width ratio of the diagram. Tall diagrams and small external lags (inlet and diffusor) favor recovery. Suggestions for future theoretical and experimental research are discussed.

Moore, F. K.↗
DOE OSTI2024

A numerical method for simulating variable density flows in membrane desalination systems

Here, we present a novel method for simulating unsteady, variable density, fluid flows in membrane desalination systems. By assuming the density varies only with concentration and temperature, the scheme decouples the solution of the governing equations into two sequential blocks. The first solves the governing equations for the temperature and concentration fields, which are used to compute all thermophysical properties. The second block solves the conservation of mass and momentum equations for the velocity and pressure. We show that this is computationally more efficient than schemes that iterate over the full coupled equations in one block. We verify that the method achieves second-order spatial-temporal accuracy, and we use the method to investigate buoyancy-driven convection in a desalination process called vacuum membrane distillation. Specifically, we show that with gravity properly oriented, variations in temperature and concentration can trigger a double-diffusive instability that enhances mixing and improves water recovery. We also show that the instability can be strengthened by providing external heating.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A Dynamic Trajectory Fit to Multisensor Fireball Observations

Meteorites with known orbital origins are key to our understanding of solar system formation and the source of life on Earth. Fireball networks have been developed globally in a unified effort to record and ultimately retrieve these cosmic samples. However, the accuracy of the determined orbit and the likelihood of meteorite recovery depend directly on the accuracy of the chosen meteoroid triangulation method. There are three leading techniques for meteoroid triangulation discussed in the literature: the method of planes, the straight-line least-squares method, and the multiparameter fit method. Here we describe an alternative method to meteoroid triangulation, called the dynamic trajectory fit. This approach uses the meteoroid’s 3D dynamic equations of motion to fit a realistic trajectory directly to multisensor line-of-sight observations. This method has the ability to resolve fragmentation events, fit systematic observatory timing offsets, and determine mass estimates of the meteoroid along its observable trajectory. Through a comprehensive Monte Carlo analysis of over 100,000 trajectory simulations, we find this new method to more accurately estimate meteoroid trajectories of slow entry events (<25 km s{sup −1}) and events observed from low convergence angles (<10°) compared to existing meteoroid triangulation techniques. Additionally, we triangulate an observed fireball event with visible fragmentation using the various triangulation methods to show that the proposed dynamic trajectory fit implementing fragmentation to best match the captured multisensor line-of-sight data.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1976

A numerical solution of the flow field over a transonic airfoil including strong-shock-induced flow separation

An approximate numerical method has been developed for computing the entire transonic flow field over an airfoil, including the boundary layer, the shock-wave/boundary-layer interaction, and the shock-induced-separation. The inviscid flow, the turbulent boundary layer in the weak interaction zone, and the strong interaction elliptic region are solved separately. The inviscid flow is represented by the solution of the transonic small-disturbance equation, while the boundary layer is solved by the generalized Galerkin method, and the strong interaction region (including the wake) is computed from the time-averaged Navier-Stokes equations. A second-order interaction is computed iteratively using the inviscid pressure distribution as a boundary condition for the boundary layer and Navier-Stokes codes, and the viscous displacement thickness to modify the airfoil geometry for the inviscid computation. The computed pressure distributions and shock-wave profiles are compared with experimental data. Agreement is very good except that the pressure recovery on the trailing edge is greater than the experimental data.

Seginer, A.↗
NASA NTRS1985

Accelerated viscoelastic characterization of T300-5208 graphite-epoxy laminates

A viscoelastic response scheme for the accelerated characterization of polymer-based composite laminates in applied to T300/5208 graphite/epoxy. The response of uni-directional specimens is modeled. The transient component of the viscoelastic creep compliance is assumed to follow a power law approximation. A recursive relationship is developed, based upon the Schapery single-integral equation, which allows approximation of a continuous time-varying uniaxial load using discrete steps in stress. The viscoelastic response of T300/5208 to transverse normal and shear stresses is determined unsing 90 deg and 10 deg off-axis tensile specimens. In each case the seven viscoelastic material parameters required in the analysis are determined experimentally using short-term creep and creep recovery tests. It is shown that an accurate measure of the power law exponent is crucial for accurate long-term prediction. A short term test cycle selection procedure is proposed, which should provide useful guidelines for the evaluation of other viscoelastic materials.

Tuttle, M. E.↗
NASA NTRS1995

Plasmasphere Modeling with Ring Current Heating

Coulomb collisions between ring current ions and the thermal plasma in the plasmasphere will heat the plasmaspheric electrons and ions. During a storm such heating would lead to significant changes in the temperature and density of the thermal plasma. This was modeled using a time- dependent, one-stream hydrodynamic model for plasmaspheric flows, in which the model flux tube is connected to the ionosphere. The model simultaneously solves the coupled continuity, momentum, and energy equations of a two-ion (H(+) and O(+) quasineutral, currentless plasma. Heating rates due to collisions with ring current ions were calculated along the field line using a kinetic ring current model. First, diurnally reproducible results were found assuming only photoelectron heating of the thermal electrons. Then results were found with heating of the H(+) ions by the ring current during the recovery phase of a magnetic storm.

Guiter, S. M.↗
NASA NTRS1973

Recovery of sounding rocket payloads by center-of-gravity position control.

A computerized procedure has been developed to predict dynamic pressure at parachute deployment altitude. Inputs to the program are payload geometry, center-of-gravity position, velocity, and altitude at start of reentry. Payload is assumed to be trimmed, and potential and cross-flow drag theory is used to compute trim angle-of-attack and trim drag coefficient as a function of free-stream Mach number. The computations involve only simple geometric quantities and closed-form equations. Payload flight path is assumed to be vertical, and equations of motion of the payload are integrated from start of reentry to chute deployment altitude.

Mcgarvey, J. F.↗
NASA NTRS1988

Anisotropic constitutive modeling for nickel-base single crystal superalloys

An anisotropic constitutive model was developed based on crystallographic slip theory for nickel base single crystal superalloys. The constitutive equations developed utilizes drag stress and back stress state variables to model the local inelastic flow. Specially designed experiments were conducted to evaluate the existence of back stress in single crystal superalloy Rene N4 at 982 C. The results suggest that: (1) the back stress is orientation dependent; and (2) the back stress state variable is required for the current model to predict material anelastic recovery behavior. The model was evaluated for its predictive capability on single crystal material behavior including orientation dependent stress-strain response, tension/compression asymmetry, strain rate sensitivity, anelastic recovery behavior, cyclic hardening and softening, stress relaxation, creep and associated crystal lattice rotation. Limitation and future development needs are discussed.

Sheh, Michael Y.↗
DOE OSTI2023

Experimental and theoretical investigations of rotating algae biofilm reactors (RABRs): Areal productivity, nutrient recovery, and energy efficiency

Microalgae biofilms have been demonstrated to recover nutrients from wastewater and serve as biomass feedstock for bioproducts. However, there is a need to develop a platform to quantitatively describe microalgae biofilm production, which can provide guidance and insights for improving biomass areal productivity and nutrient uptake efficiency. Here, this paper proposes a unified experimental and theoretical framework to investigate algae biofilm growth on a rotating algae biofilm reactor (RABR). Experimental laboratory setups are used to conduct controlled experiments on testing environmental and operational factors for RABRs. We propose a differential–integral equation‐based mathematical model for microalgae biofilm cultivation guided by laboratory experimental findings. The predictive mathematical model development is coordinated with laboratory experiments of biofilm areal productivity associated with ammonia and inorganic phosphorus uptake by RABRs. The unified experimental and theoretical tool is used to investigate the effects of RABR rotating velocity, duty cycle (DC), and light intensity on algae biofilm growth, areal productivity, nutrient uptake efficiency, and energy efficiency in wastewater treatment. Our framework indicates that maintaining a reasonable light intensity range improves biomass areal productivity and nutrient uptake efficiency. Our framework also indicates that faster RABR rotation benefits biomass areal productivity. However, maximizing the nutrient uptake efficiency requires a reasonably low RABR rotating speed. Energy efficiency is strongly correlated with RABR rotating speed and DC.

09 BIOMASS FUELS↗
DOE OSTI2020

A Novel Model of Thyristor in Reverse Recovery Process for Current Interruption Test of HVDC Circuit Breakers

The high voltage thyristor is widely used in high voltage DC (HVDC) grids due to its large current capability. Recently, the high voltage thyristor is implemented to produce a high current in the HVDC circuit breaker (CB) testing process. Unfortunately, the reverse recovery process (RRP) of a thyristor could induce significant turn-off overvoltage and power losses, jeopardizing the safety of the devices. The existing models of a thyristor cannot describe such hard-switching transient accurately. Therefore, this paper divides this process into four stages based on the physical mechanism and proposes a novel cosine exponential (CE), simulating the voltage and current of the thyristor with analytical equations, which are derived from the equivalent circuit. Moreover, the extraction procedure of the proposed model is illustrated in this paper. Finally, we develop an HVDC CB test platform to validate the proposed model. Experimental results show that the proposed CE model can accurately predict the overvoltage with a relative error of less than 8% at various current varying rates.

Dongye, Zhonghao↗
NASA NTRS1973

Parameter estimation supplement to the Mission Analysis Evaluation and Space Trajectory Operations program (MAESTRO)

This Parameter Estimation Supplement describes the PEST computer program and gives instructions for its use in determination of lunar gravitation field coefficients. PEST was developed for use in the RAE-B lunar orbiting mission as a means of lunar field recovery. The observations processed by PEST are short-arc osculating orbital elements. These observations are the end product of an orbit determination process obtained with another program. PEST's end product it a set of harmonic coefficients to be used in long-term prediction of the lunar orbit. PEST employs some novel techniques in its estimation process, notably a square batch estimator and linear variational equations in the orbital elements (both osculating and mean) for measurement sensitivities. The program's capabilities are described, and operating instructions and input/output examples are given. PEST utilizes MAESTRO routines for its trajectory propagation. PEST's program structure and subroutines which are not common to MAESTRO are described. Some of the theoretical background information for the estimation process, and a derivation of linear variational equations for the Method 7 elements are included.

Bjorkman, W. S.↗
DOE OSTI2020

The Sylvester Graphical Lasso (SyGlasso)

This paper introduces the Sylvester graphical lasso (SyGlasso) that captures multiway dependencies present in tensor-valued data. The model is based on the Sylvester equation that defines a generative model. The proposed model complements the tensor graphical lasso that imposes a Kronecker sum model for the inverse covariance matrix by providing an alternative Kronecker sum model that is generative and interpretable. A nodewise regression approach is adopted for estimating the conditional independence relationships among variables. The statistical convergence of the method is established, and empirical studies are provided to demonstrate the recovery of meaningful conditional dependency graphs. Here we apply the SyGlasso to an electroencephalography (EEG) study to compare the brain connectivity of alcoholic and nonalcoholic subjects. Here, we demonstrate that our model can simultaneously estimate both the brain connectivity and its temporal dependencies.

Wang, Yu↗
NASA NTRS2001

A Self-Consistent Model of the Interacting Ring Current Ions and Electromagnetic ICWs. Initial Results: Waves and Precipitation Fluxes

Initial results from the new developed model of the interacting ring current ions and ion cyclotron waves are presented. The model described by the system of two bound kinetic equations: one equation describes the ring current ion dynamics, and another one gives wave evolution. Such system gives a self-consistent description of the ring current ions and ion cyclotron waves in a quasilinear approach. Calculating ion-wave relationships, on a global scale under non steady-state conditions during May 2-5, 1998 storm, we presented the data at three time cuts around initial, main, and late recovery phases of May 4, 1998 storm phase. The structure and dynamics of the ring current proton precipitating flux regions and the wave active ones are discussed in detail.

Khazanov, G. V.↗
NASA NTRS1975

Hot-cold plasma interactions in the earth's magnetosphere

The paper gives a synopsis of observations and results obtained from the Explorer-45 satellite which pertain directly to the interaction of the magnetospheric hot and cold plasma populations. The specific case discussed is the interaction of the hot ring current plasma with the cold plasmaspheric plasma in the evening to early morning local time sector during magnetic storm recovery phase. It was found that above the plasmapause region, the hot ring current plasma is stable with negligible losses due to pitch angle diffusion; the hot ring current plasma enters a moderate pitch angle regime in the plasmapause region, whereby the addition of cold plasma destabilizes the hot plasma. Analysis of the energy, spatial, and temporal dependence of the above destabilization along with the ion-cyclotron resonant energy equation, and comparison of this analysis with an in situ estimate of the plasma density strongly indicate that the mechanism responsible for destabilization of the hot plasma is the amplification of ion-cyclotron waves due to the interaction of the cold plasmaspheric plasma with the hot ring current plasma in a manner similar to that discussed by Cornwall et al. (1970).

Williams, D. J.↗