Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “non linear dynamics”

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 523 records · Page 29

The adequate stimulus for avian short latency vestibular responses to linear translation

Transient linear acceleration stimuli have been shown to elicit eighth nerve vestibular compound action potentials in birds and mammals. The present study was undertaken to better define the nature of the adequate stimulus for neurons generating the response in the chicken (Gallus domesticus). In particular, the study evaluated the question of whether the neurons studied are most sensitive to the maximum level of linear acceleration achieved or to the rate of change in acceleration (da/dt, or jerk). To do this, vestibular response thresholds were measured as a function of stimulus onset slope. Traditional computer signal averaging was used to record responses to pulsed linear acceleration stimuli. Stimulus onset slope was systematically varied. Acceleration thresholds decreased with increasing stimulus onset slope (decreasing stimulus rise time). When stimuli were expressed in units of jerk (g/ms), thresholds were virtually constant for all stimulus rise times. Moreover, stimuli having identical jerk magnitudes but widely varying peak acceleration levels produced virtually identical responses. Vestibular response thresholds, latencies and amplitudes appear to be determined strictly by stimulus jerk magnitudes. Stimulus attributes such as peak acceleration or rise time alone do not provide sufficient information to predict response parameter quantities. Indeed, the major response parameters were shown to be virtually independent of peak acceleration levels or rise time when these stimulus features were isolated and considered separately. It is concluded that the neurons generating short latency vestibular evoked potentials do so as "jerk encoders" in the chicken. Primary afferents classified as "irregular", and which traditionally fall into the broad category of "dynamic" or "phasic" neurons, would seem to be the most likely candidates for the neural generators of short latency vestibular compound action potentials.

NASA Discipline Neuroscience↗

Extending the Operational Envelope of a Turbofan Engine Simulation into the Sub-Idle Region

In many non-linear gas turbine simulations, operation in the sub-idle region can lead to model instability. This paper lays out a method for extending the operational envelope of a map based gas turbine simulation to include the sub-idle region. This method develops a multi-simulation solution where the baseline component maps are extrapolated below the idle level and an alternate model is developed to serve as a safety net when the baseline model becomes unstable or unreliable. Sub-idle model development takes place in two distinct operational areas, windmilling/shutdown and purge/cranking/startup. These models are based on derived steady state operating points with transient values extrapolated between initial (known) and final (assumed) states. Model transitioning logic is developed to predict baseline model sub-idle instability, and transition smoothly and stably to the backup sub-idle model. Results from the simulation show a realistic approximation of sub-idle behavior as compared to generic sub-idle engine performance that allows the engine to operate continuously and stably from shutdown to full power.

gas path dynamics↗

Stimulus configuration determines the detectability of motion signals in noise

We measured the detectability of moving signal dots in dynamic noise to determine whether local motion signals are preferentially combined along an axis parallel to the direction of motion. Observers were asked to detect a signal composed of three dots moving in a linear trajectory among dynamic noise dots. The signal dots were collinear and equally spaced in a configuration that was either parallel to or perpendicular to their trajectory. The probability of detecting the signal was measured as a function of noise density, over a range of signal dot spacings from 0.5 degrees to 5.0 degrees. At any given noise density, the signal in the parallel configuration was more detectable than that in the perpendicular configuration. Our four observers could tolerate 1.5-2.5 times more noise in the parallel configuration. This improvement is not due merely to temporal summation between consecutive dots in the parallel trajectory. Temporal summation functions measured on our observers indicate that the benefit from spatial coincidence of the dots lasts for no more than 50 ms, whereas the increased detectability of the parallel configuration is observed up to the largest temporal separations tested (210 ms). These results demonstrate that dots arranged parallel to the signal trajectory are more easily detected than those arranged perpendicularly. Moreover, this enhancement points to the existence of visual mechanisms that preferentially organize motion information parallel to the direction of motion.

Non-NASA Center↗

M.I.T./Canadian vestibular experiments on the Spacelab-1 mission: 6. Vestibular reactions to lateral acceleration following ten days of weightlessness

Tests of otolith function were performed pre-flight and post-flight on the science crew of the first Spacelab Mission with a rail-mounted linear acceleration sled. Four tests were performed using horizontal lateral (y-axis) acceleration: perception of linear motion, a closed loop nulling task, dynamic ocular torsion, and lateral eye deviations. The motion perception test measured the time to detect the onset and direction of near threshold accelerations. Post-flight measures of threshold and velocity constant obtained during the days immediately following the mission showed no consistent pattern of change among the four crewmen compared to their pre-flight baseline other than an increased variability of response. In the closed loop nulling task, crewmen controlled the motion of the sled and attempted to null a computer-generated random disturbance motion. When performed in the light, no difference in ability was noted between pre-flight and post-flight. In the dark, however, two of the four crewmen exhibited somewhat enhanced performance post-flight. Dynamic ocular torsion was measured in response to sinusoidal lateral acceleration which produces a gravitionertial stimulus equivalent to lateral head tilt without rotational movement of the head. Results available for two crewmen suggest a decreased amplitude of sinusoidal ocular torsion when measured on the day of landing (R+0) and an increasing amplitude when measured during the week following the mission.

short duration↗

Extending the Operational Envelope of a Turbofan Engine Simulation into the Sub-Idle Region

In many non-linear gas turbine simulations, operation in the sub-idle region can lead to model instability. This paper lays out a method for extending the operational envelope of a map based gas turbine simulation to include the sub-idle region. This method develops a multi-simulation solution where the baseline component maps are extrapolated below the idle level and an alternate model is developed to serve as a safety net when the baseline model becomes unstable or unreliable. Sub-idle model development takes place in two distinct operational areas, windmilling/shutdown and purge/cranking/ startup. These models are based on derived steady state operating points with transient values extrapolated between initial (known) and final (assumed) states. Model transitioning logic is developed to predict baseline model sub-idle instability, and transition smoothly and stably to the backup sub-idle model. Results from the simulation show a realistic approximation of sub-idle behavior as compared to generic sub-idle engine performance that allows the engine to operate continuously and stably from shutdown to full power.

aerothermodynamics↗

Assessment of CFD-based Response Surface Model for Ares I Supersonic Ascent Aerodynamics

The Ascent Force and Moment Aerodynamic (AFMA) Databases (DBs) for the Ares I Crew Launch Vehicle (CLV) were typically based on wind tunnel (WT) data, with increments provided by computational fluid dynamics (CFD) simulations for aspects of the vehicle that could not be tested in the WT tests. During the Design Analysis Cycle 3 analysis for the outer mold line (OML) geometry designated A106, a major tunnel mishap delayed the WT test for supersonic Mach numbers (M) greater than 1.6 in the Unitary Plan Wind Tunnel at NASA Langley Research Center, and the test delay pushed the final delivery of the A106 AFMA DB back by several months. The aero team developed an interim database based entirely on the already completed CFD simulations to mitigate the impact of the delay. This CFD-based database used a response surface methodology based on radial basis functions to predict the aerodynamic coefficients for M > 1.6 based on only the CFD data from both WT and flight Reynolds number conditions. The aero team used extensive knowledge of the previous AFMA DB for the A103 OML to guide the development of the CFD-based A106 AFMA DB. This report details the development of the CFD-based A106 Supersonic AFMA DB, constructs a prediction of the database uncertainty using data available at the time of development, and assesses the overall quality of the CFD-based DB both qualitatively and quantitatively. This assessment confirms that a reasonable aerodynamic database can be constructed for launch vehicles at supersonic conditions using only CFD data if sufficient knowledge of the physics and expected behavior is available. This report also demonstrates the applicability of non-parametric response surface modeling using radial basis functions for development of aerodynamic databases that exhibit both linear and non-linear behavior throughout a large data space.

Hanke, Jeremy L.↗

A dynamical model for reflex activated head movements in the horizontal plane

We present a controls systems model of horizontal-plane head movements during perturbations of the trunk, which for the first time interfaces a model of the human head with neural feedback controllers representing the vestibulocollic (VCR) and the cervicocollic (CCR) reflexes. This model is homeomorphic such that model structure and parameters are drawn directly from anthropomorphic, biomechanical and physiological studies. Using control theory we analyzed the system model in the time and frequency domains, simulating neck movement responses to input perturbations of the trunk. Without reflex control, the head and neck system produced a second-order underdamped response with a 5.2 dB resonant peak at 2.1 Hz. Adding the CCR component to the system dampened the response by approximately 7%. Adding the VCR component dampened head oscillations by 75%. The VCR also improved low-frequency compensation by increasing the gain and phase lag, creating a phase minimum at 0.1 Hz and a phase peak at 1.1 Hz. Combining all three components (mechanics, VCR and CCR) linearly in the head and neck system reduced the amplitude of the resonant peak to 1.1 dB and increased the resonant frequency to 2.9 Hz. The closed loop results closely fit human data, and explain quantitatively the characteristic phase peak often observed.

NASA Discipline Neuroscience↗

Comparative Study on High-Order Positivity-preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order WENO schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more difficult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium flow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

WENO↗

Function of human eccrine sweat glands during dynamic exercise and passive heat stress

The purpose of this study was to identify the pattern of change in the density of activated sweat glands (ASG) and sweat output per gland (SGO) during dynamic constant-workload exercise and passive heat stress. Eight male subjects (22.8 +/- 0.9 yr) exercised at a constant workload (117.5 +/- 4.8 W) and were also passively heated by lower-leg immersion into hot water of 42 degrees C under an ambient temperature of 25 degrees C and relative humidity of 50%. Esophageal temperature, mean skin temperature, sweating rate (SR), and heart rate were measured continuously during both trials. The number of ASG was determined every 4 min after the onset of sweating, whereas SGO was calculated by dividing SR by ASG. During both exercise and passive heating, SR increased abruptly during the first 8 min after onset of sweating, followed by a slower increase. Similarly for both protocols, the number of ASG increased rapidly during the first 8 min after the onset of sweating and then ceased to increase further (P > 0.05). Conversely, SGO increased linearly throughout both perturbations. Our results suggest that changes in forearm sweating rate rely on both ASG and SGO during the initial period of exercise and passive heating, whereas further increases in SR are dependent on increases in SGO.

NASA Discipline Cardiopulmonary↗

Vestibular adaptation to space in monkeys

Otolith-induced eye movements of rhesus monkeys were studied before and after the 1989 COSMOS 2044 and the 1992 to 1993 COSMOS 2229 flights. Two animals flew in each mission for approximately 2 weeks. After flight, spatial orientation of the angular vestibulo-ocular reflex was altered. In one animal the time constant of postrotatory nystagmus, which had been shortened by head tilts with regard to gravity before flight, was unaffected by the same head tilts after flight. In another animal, eye velocity, which tended to align with a gravitational axis before flight, moved toward a body axis after flight. This shift of orientation disappeared by 7 days after landing. After flight, the magnitude of compensatory ocular counter-rolling was reduced by about 70% in both dynamic and static tilts. Modulation in vergence in response to naso-occipital linear acceleration during off-vertical axis rotation was reduced by more than 50%. These changes persisted for 11 days after recovery. An up and down asymmetry of vertical nystagmus was diminished for 7 days. Gains of the semicircular canal-induced horizontal and vertical angular vestibulo-ocular reflexes were unaffected in both flights, but the gain of the roll angular vestibulo-ocular reflex was decreased. These data indicate that there are short- and long-term changes in otolith-induced eye movements after adaptation to microgravity. These experiments also demonstrate the unique value of the monkey as a model for studying effects of vestibular adaptation in space. Eye movements can be measured in three dimensions in response to controlled vestibular and visual stimulation, and the results are directly applicable to human beings. Studies in monkeys to determine how otolith afferent input and central processing is altered by adaptation to microgravity should be an essential component of future space-related research.

Cosmos 2229 Project↗

A model for the characterization of the spatial properties in vestibular neurons

Quantitative study of the static and dynamic response properties of some otolith-sensitive neurons has been difficult in the past partly because their responses to different linear acceleration vectors exhibited no "null" plane and a dependence of phase on stimulus orientation. The theoretical formulation of the response ellipse provides a quantitative way to estimate the spatio-temporal properties of such neurons. Its semi-major axis gives the direction of the polarization vector (i.e., direction of maximal sensitivity) and it estimates the neuronal response for stimulation along that direction. In addition, the semi-minor axis of the ellipse provides an estimate of the neuron's maximal sensitivity in the "null" plane. In this paper, extracellular recordings from otolith-sensitive vestibular nuclei neurons in decerebrate rats were used to demonstrate the practical application of the method. The experimentally observed gain and phase dependence on the orientation angle of the acceleration vector in a head-horizontal plane was described and satisfactorily fit by the response ellipse model. In addition, the model satisfactorily fits neuronal responses in three-dimensions and unequivocally demonstrates that the response ellipse formulation is the general approach to describe quantitatively the spatial properties of vestibular neurons.

Non-NASA Center↗

Feasibility of Varying Geo-Fence Around an Unmanned Aircraft Operation Based on Vehicle Performance and Wind

Managing trajectory separation is critical to ensuring accessibility, efficiency, and safety in the unmanned airspace. The notion of geo-fences is an emerging concept, where distance buffers enclose individual trajectories and areas of operation in order to manage the airspace. Currently, the Air Traffic Management system for commercial travel defines static distance buffers around the aircraft; however, commercial UASs are envisioned to operate in significantly closer proximity to other UAS requiring a geo-fence for spacing operations. The geo-fence size can be determined based on vehicle performance characteristics, state of the airspace, weather, and other unforeseen events such as emergency or disaster response. Calculation of the geo-fence size could be determined as part of pre-flight planning and during real-time operations. A largely non-homogeneous fleet of UASs will be operating in low altitude and will likely be commercially developed. Due to intellectual property concerns, the operators may not provide detailed specifications of the control system to UTM. In addition, the huge variety of UAS makes modeling each control system prohibitive and flight data for these vehicles may not exist. Therefore, a generalized, simple geo-fence sizing algorithm must be developed such that it does not rely on detailed knowledge of the vehicle control system, accounts for the presence of urban winds, and is sufficiently accurate. In this work, two simple models are investigated to determine its feasibility as an adequate means for calculating the geo-fence size. The vehicle data used in this work are provided by UAS manufactures who have partnered with NASA's UTM project and some publicly available websites. The first model utilizes wind data processed from the NOAA HRRR (Hourly Rapid Refresh) product and Sonar Annemometer data provided by San Jose State. The second model utilizes OpenFOAM which is a CFD code used to generate a wind field for flow around a single building. The key vehicle performance parameters can include UAS response time to disturbances, command to actuation latency, control system rate limits, time to recovery to desired path, and aerodynamics. It was found that the first model provides an initial understanding of geo-fence sizing, but does not provide enough accuracy to provide UTM with an efficient means of scheduling vehicles. The results of the second model reveal that modeling UAS controls systems with a linearized plant and gain scheduled PID controller does not allow capture the UAS flight dynamics within a significant envelope of the wind disturbances.

uncertainty↗

Instrumentation for investigation of the depth-dose distribution by the Liulin-5 instrument of a human phantom on the Russian segment of ISS for estimation of the radiation risk during long term space flights

Described is the Liulin-5 experiment and instrumentation, developed for investigation of the space radiation doses depth distribution in a human phantom on the Russian Segment of the International Space Station (ISS). Liulin-5 experiment is a part of the international project MATROSHKA-R on ISS. The experiment MATROSHKA-R is aimed to study the depth dose distribution at the sites of critical organs of the human body, using models of human body-anthropomorphic and spherical tissue-equivalent phantoms. The aim of Liulin-5 experiment is long term (4-5 years) investigation of the radiation environment dynamics inside the spherical tissue-equivalent phantom, mounted in different places of the Russian Segment of ISS. Energy deposition spectra, linear energy transfer spectra, flux and dose rates for protons and the biologically-relevant heavy ion components of the galactic cosmic radiation will be measured simultaneously with near real time resolution at different depths of the phantom by a telescope of silicon detectors. Data obtained together with data from other active and passive dosimeters will be used to estimate the radiation risk to the crewmembers, verify the models of radiation environment in low Earth orbit, validate body transport model and correlate organ level dose to skin dose. Presented are the test results of the prototype unit. The spherical phantom will be flown on the ISS in 2004 year and Liulin-5 experiment is planned for 2005 year. c2004 COSPAR. Published by Elsevier Ltd. All rights reserved.

Cosmic Radiation↗

Dynamic posture analysis of Spacelab-1 crew members

Dynamic posture testing was conducted on the science crew of the Spacelab-1 mission on a single axis linear motion platform. Tests took place in pre- and post-flight sessions lasting approximately 20 min each. The pre-flight tests were widely spaced over the several months prior to the mission while the post-flight tests were conducted over the first, second, fourth, and sixth days after landing. Two of the crew members were also tested on the day of landing. Consistent with previous postural testing conducted on flight crews, these crew members were able to complete simple postural tasks to an acceptable level even in the first few hours after landing. Our tests were designed to induce dynamic postural responses using a variety of stimuli and from these responses, evaluate subtle changes in the postural control system which had occurred over the duration of the flight. Periodic sampling post-flight allowed us to observe the time course of readaptation to terrestrial life. Our observations of hip and shoulder position, when subjected to careful analysis, indicated modification of the postural response from pre- to post-flight and that demonstrable adjustments in the dynamic control of their postural systems were taking place in the first few days after flight. For transient stimuli where the platform on which they were asked to stand quickly moved a few centimeters fore or aft then stopped, ballistic or open loop 'programs' would closely characterize the response. During these responses the desired target position was not always achieved and of equal importance not always properly corrected some 15 seconds after the platform ceased to move. The persistent observation was that the subjects had a much stronger dependence on visual stabilization post-flight than pre-flight. This was best illustrated by a slow or only partial recovery to an upward posture after a transient base-of-support movement with eyes open. Postural responses to persistent wideband pseudorandom base-of-support translation were modeled as time invarient linear systems arrived at by Kalman adaptive filter techniques. Derived model parameters such as damping factor and fundamental frequency of the closed loop system showed significant modification between pre- and post-flight. This phenomenon is best characterized by movement of the poles toward increasing stability. While pre-flight data tended to show shoulders and hips moving in phase with each other, post-flight data showed a more disjoint behavior.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗

Short latency compound action potentials from mammalian gravity receptor organs

Gravity receptor function was characterized in four mammalian species using far-field vestibular evoked potentials (VsEPs). VsEPs are compound action potentials of the vestibular nerve and central relays that are elicited by linear acceleration ramps applied to the cranium. Rats, mice, guinea pigs, and gerbils were studied. In all species, response onset occurred within 1.5 ms of the stimulus onset. Responses persisted during intense (116 dBSPL) wide-band (50 to 50 inverted question mark omitted inverted question mark000 Hz) forward masking, whereas auditory responses to intense clicks (112 dBpeSPL) were eliminated under the same conditions. VsEPs remained after cochlear extirpation but were eliminated following bilateral labyrinthectomy. Responses included a series of positive and negative peaks that occurred within 8 ms of stimulus onset (range of means at +6 dBre: 1.0 g/ms: P1=908 to 1062 micros, N1=1342 to 1475 micros, P2=1632 to 1952 micros, N2=2038 to 2387 micros). Mean response amplitudes at +6 dBre: 1.0 g/ms ranged from 0.14 to 0.99 microV. VsEP input/output functions revealed latency slopes that varied across peaks and species ranging from -19 to -51 micros/dB. Amplitude-intensity slopes also varied ranging from 0.04 to 0.08 microV/dB for rats and mice. Latency values were comparable to those of birds although amplitudes were substantially smaller in mammals. VsEP threshold values were considerably higher in mammals compared to birds and ranged from -8.1 to -10.5 dBre 1.0 g/ms across species. These results support the hypothesis that mammalian gravity receptors are less sensitive to dynamic stimuli than are those of birds.

Non-NASA Center↗

The dynamic contributions of the otolith organs to human ocular torsion

We measured human ocular torsion (OT) monocularly (using video) and binocularly (using search coils) while sinusoidally accelerating (0.7 g) five human subjects along an earth-horizontal axis at five frequencies (0.35, 0.4, 0.5, 0.75, and 1.0 Hz). The compensatory nature of OT was investigated by changing the relative orientation of the dynamic (linear acceleration) and static (gravitational) cues. Four subject orientations were investigated: (1) Y-upright-acceleration along the interaural (y) axis while upright; (2) Y-supine-acceleration along the y-axis while supine; (3) Z-RED-acceleration along the dorsoventral (z) axis with right ear down; (4) Z-supine-acceleration along the z-axis while supine. Linear acceleration in the Y-upright, Y-supine and Z-RED orientations elicited conjugate OT. The smaller response in the Z-supine orientation appeared disconjugate. The amplitude of the response decreased and the phase lag increased with increasing frequency for each orientation. This frequency dependence does not match the frequency response of the regular or irregular afferent otolith neurons; therefore the response dynamics cannot be explained by simple peripheral mechanisms. The Y-upright responses were larger than the Y-supine responses (P < 0.05). This difference indicates that OT must be more complicated than a simple low-pass filtered response to interaural shear force, since the dynamic shear force along the interaural axis was identical in these two orientations. The Y-supine responses were, in turn, larger than the Z-RED responses (P < 0.01). Interestingly, the vector sum of the Y-supine responses plus Z-RED responses was not significantly different (P = 0.99) from the Y-upright responses. This suggests that, in this frequency range, the conjugate OT response during Y-upright stimulation might be composed of two components: (1) a response to shear force along the y-axis (as in Y-supine stimulation), and (2) a response to roll tilt of gravitoinertial force (as in Z-RED stimulation).

NASA Discipline Neuroscience↗

Chemical Thermodynamics and the Mathematical Integration of Reaction Kinetics

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

STMD↗

Simultaneous concentration and purification through gradient deformation chromatography

Mobile-phase additives, commonly used to modulate absorbate retention in gradient elution chromatography, are usually assumed to be either linearly retained or unretained. Previous theoretical work from our laboratory has shown that these modulators, such as salts in ion-exchange and hydrophobic interaction chromatography and organic modifiers in reversed-phase chromatography, can absorb nonlinearly, giving rise to gradient deformation. Consequently, adsorbate peaks that elute in the vicinity of the head of the deformed gradient may exhibit unusual shapes, form shoulders, and/or be concentrated. These effects for a reversed-phase sorbent with aqueous acetonitrile (ACN) as the modulator are verified experimentally. Gradient deformation is demonstrated experimentally and agrees with simulations based on ACN isotherm parameters that are independently determined from batch equilibrium studies using the layer model. Unusual absorbate peak shapes were found experimentally for single-component injections of phenylalanine, similar to those calculated by the simulations. A binary mixture of tryptophan and phenylalanine is used to demonstrate simultaneous concentration and separation, again in agreement with simulations. The possibility of gradient deformation in ion-exchange and hydrophobic interaction chromatography is discussed.

NASA Discipline Number 93-10↗