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 703 records · Page 39

Efficient Parametric Uncertainty Analysis of an Earth Entry Vehicle Concept Using Least Angle Regression

The objective of this work was to outline and apply an efficient and accurate parametric un-certainty propagation approach to the analysis of convective heating on an Earth entry vehicle concept. The described approach was based on Least Angle Regression used to solve a sparse and underdetermined linear system in the point-collocation non-intrusive polynomial chaos surrogate method. This approach involved an iterative process to computing the non-zero terms of the underlying polynomial chaos model using only enough samples to converge uncertainty interval predictions and Sobol index values based global nonlinear sensitivity estimates. The Earth entry vehicle was analyzed at three points along a representative trajectory for a Mars return mission. 329 sources of uncertainty were identified in the computational fluid dynamics model used to predict the forebody convective heating. These included uncertainty in flow field chemical rates, collision integrals, heats of formation, surface finite rate char model reaction rates, wall roughness height, and the turbulent Schmidt number. Results from this study showed that convective heating uncertainty as high as 50% of the nominal was predicted with only about 50 evaluations of the computational model. This was far fewer than would be required for a sampling-based approach or a full basis polynomial chaos model, which would have required over 50,000 samples. Additionally, results showed that over 90% of the total convective heating uncertainty was due to uncertainty in the N2catalytic rate on the surface, while the remainder of the uncertainty was attributed to the turbulent Schmidt number and the wall roughness uncertainties.

Thomas K West IV↗

Quantification of Numerical Uncertainty via Nonlinear Dynamical Approach

Motivations (Ensure a Higher Level of Confidence in the Predictability & Reliability of Numerical Simulation for Multiscale Complex Nonlinear Fluid Problems) - The last two decades have been an era when computation is ahead of analysis & when very large scale practical computations are increasingly used in poorly understood multiscale complex nonlinear physical problems & non-traditional fields (Especially when computations offer the ONLY way of generating this type of data limited simulations). - At present some of the numerical uncertainties can be explained and minimized by traditional numerical analysis and standard CFD practices. However, such practices, usually based on linearized analysis, MIGHT NOT be sufficient for strongly nonlinear and/or stiff problems. - We need a good understanding of the nonlinear behavior of numerical schemes being used as an integral part of code verification, validation and certification.

HEC↗

Simulation of convective transport during frequency chirping of a TAE using the MEGA code

Here, we present a procedure to examine energetic particle phase-space during long range frequency chirping phenomena in tokamak plasmas. To apply the proposed method, we have performed self-consistent simulations using the MEGA code and analyzed the simulation data. We demonstrate a traveling wave in phase-space and that there exist specific slices of phase-space on which the resonant particles lie throughout the wave evolution. For non-linear evolution of an n = 6 toroidicity-induced Alfvén eigenmode (TAE), our results reveal the formation of coherent phase-space structures (holes/clumps) after coarse-graining of the distribution function. These structures cause a convective transport in phase-space which implies a radial drift of the resonant particles. We also demonstrate that the rate of frequency chirping increases with the TAE damping rate. Our observations of the TAE behavior and the corresponding phase-space dynamics are consistent with the Berk–Breizman theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Atomic Simulation of Complex DNA DSBs and the Interactions with the Ku70/80 Heterodimer

DNA double strand breaks (DSBs) induced by ionizing radiation (IR) usually contain modified bases such as 8-oxo-7,8-dihydroguanine (8-oxoG) and thymine glycol, apurinic/apyrimidinic (AP) sites, 2-deoxyribonolactone, or single-strand breaks (SSBs). The presence of such lesions in close proximity to the DSB terminus makes the DNA nicks more difficult to repair and rejoin than endogenously induced simple DSBs, and as such a major determinant of the biological effects of high linear energy transfer (LET) radiation as encountered in space travel. In this study we conducted molecular dynamics simulations on a series of DNA duplexes with various complex lesions of 8-oxoG and AP sites, in an effort to investigate the effects of such lesions to the structural integrity and stability of DNA after insulted by IR. We also simulated the interaction of such complex DSBs with the Ku70/80 heterodimer, the first protein in mammalian cells to embark the non-homologous end joining (NHEJ) DNA repair pathway. The results indicate, compared to DNA with simple DSBs, the complex lesions can enhance the hydrogen bonds opening rate at the DNA terminus, and increase the mobility of the whole duplex, thus they present more deleterious effects to the genome integrity if not captured and repaired promptly in cells. Simulations also demonstrate the binding of Ku drastically reduces structural disruption and flexibility caused by the complex lesions, and the interactions of Ku with complex DSBs have a different potential energy landscape from the bound structure with simple DSB. In all complex DSBs systems, the binding of DSB terminus with Ku70 is softened while the binding of the middle duplex with Ku80 is tightened. This energy shift may help the Ku protein to secure at the DSB terminus for a longer time, so that other end processing factors or repair pathways can proceed at the lesions before NHEJ repair process starts. These atomic simulations may provide valuable new insight into the selective action of repair proteins on damaged DNA.

Hu, Shaowen↗

Simultaneous quantification of poly-dispersed anionic, amphoteric and nonionic surfactants in simulated wastewater samples using C18 high-performance liquid chromatography-quadrupole ion-trap mass spectrometry

This paper describes the development of a guantitative method for direct and simultaneous determination of three frequently encountered surfactants, amphoteric (cocoamphoacetate, CAA), anionic (sodium laureth sulfate, SLES), and nonionic (alcohol ethoxylate, AE) using a reversed-phase C18 HPLC coupled with an ESI ion-trap mass spectrometer (MS). Chemical composition, ionization characteristics and fragmentation pathways of the surfactants are presented. Positive ESI was effective for all three surfactants in agueous methanol buffered with ammonium acetate. The method enables rapid determinations in small sample volumes containing inorganic salts (up to 3.5 g L(-1)) and multiple classes of surfactants with high specificity by applying surfactant specific tandem mass spectrometric strategies. It has dynamic linear ranges of 2-60, 1.5-40, 0.8-56 mg L(-1) with R2 egual or greater than 0.999, 0.98 and 0.999 (10 microL injection) for CAA, SLES, and AE, respectively.

Chromatography, High Pressure Liquid/methods↗

Atomistic insights into the synergistic effects of tensile strain and hydroxyl group on increasing the thermal conductivity of monohydric alcohols as latent heat storage materials

The low thermal conductivity of linear monohydric alcohols significantly limits their performance when used as latent heat storage materials. It was found that imposing certain mechanical strains on nanomaterials with highly ordered structures can lead to appreciable thermal conductivity enhancement. To explore the potential of enhancing the intrinsic thermal conductivity of monohydric alcohols, a simplified ideal crystal model was employed for molecular arrangement in this work. The variation of atomistic heat transfer as a function of tensile strains was exploited by non-equilibrium molecular dynamics simulations. Furthermore, the influence of stretching on the vibration density of states and the corresponding molecular morphology were revealed. The results suggested that the thermal conductivity of monohydric alcohols with such simplified ideal crystalline structures can increase by about 170 % at the strain of 0.1. The distance between neighboring interfaces was found to determine the energy of inter-molecular interactions at the crystal grain boundaries. The hydrogen bonds energy has not reached the critical value when increasing tensile strains from 0 to 0.1. In contrast, the strongest van der Walls interactions for alkanes were obtained at the strain of 0.06 where the interfacial distances of alkanes have achieved their equilibrium distances. In addition, the interfacial heat transfer coefficient of ideal crystal monohydric alcohols at hydroxyl interfaces was found to be threefold than that between methyl interfaces, in consistent with the strength of inter-molecular interactions between functional groups at the interfaces.

36 MATERIALS SCIENCE↗

GPAW: An open Python package for electronic structure calculations

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe–Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn–Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Numerical Investigation of a Central Fuel Property Hypothesis under Boosted Spark-ignition Conditions

In the present work, a central fuel property hypothesis (CFPH), which states that fuel properties are sufficient to provide an indication of a fuel's performance irrespective of its chemical composition, was numerically investigated In particular, the objective of the study was to determine whether Research Octane Number (RON) and Motor Octane Number (MON), as fuel properties, are sufficient to describe a fuel's knock-limited performance under boosted spark-ignition (SI) conditions within the framework of CFPH. To this end, four TPRF-bioblendstock surrogates having different compositions but matched RON (= 98) and MON (= 90), were first generated using a non-linear regression model based on artificial neural network (ANN). Three unconventional bioblendstocks were included in the analysis: Di-isobutylene (DIB), Isobutanol and Anisole. Skeletal reaction mechanisms were generated for the TPRF-DIB, TPRF-isobutanol and TPRF-anisole blends from a detailed kinetic mechanism. Thereafter, numerical simulations were performed for the fuel surrogates using the skeletal mechanisms and a virtual cooperative fuel research (CFR) engine model, under a representative boosted operating condition. In the computational fluid dynamics (CFD) model, the G-equation approach was employed to track the turbulent flame front and the well-stirred reactor model combined with multi-zone binning strategy was used to capture auto-ignition in the end-gas. In addition, laminar flame speed was tabulated for each blend as a function of pressure, temperature and equivalence ratio a priori, and the lookup tables were used to prescribe laminar flame speed as an input to the G-equation model. Parametric spark timing sweeps were performed for each fuel blend to determine the corresponding knock-limited spark advance (KLSA) and 50% burn point (CA50) at the respective KLSA timing. It was observed that despite same RON, MON and engine operating conditions, the TPRF-Anisole blend exhibited markedly different knock-limited performance from the other three blends. This deviation from the octane index (DI) expectation was shown to be caused by differences in laminar flame speed (LFS). However; it was found that relatively large fuel-specific differences in LFS (> 20%) would have to be present to cause any appreciable deviation from the 01 framework. Otherwise, RON and MON would still be robust enough to predict a fuel's knock-limited performance.

Pal, Pinaki↗

On-orbit parametric identification methodology

On-orbit system identification (ID) of large space systems is essential for various reasons. For example, the complex composite structure of such systems cannot be ground-tested; their structural dynamic characteristics must be known accurately in order to accomplish active control. Furthermore, such capability can be used to characterize/identify various disturbances. The identification process is consisted of four principal elements: (1) modeling, (2) the estimation algorithm, (3) input system, and (4) measurement system. These elements are highly correlated and all togerher determine the success of the identification problem. Accurate modeling of large space systems is the most important element of the identification process. Large flexible structures are non-linear and infinite dimensional systems with highly coupled parameters and low frequency packed modes. In addition, these systems are subject to stochastic and time-varying disturbances, they have structural parameters which can vary due to on-orbit assembly deployment, and operations. These systems are generally; however, represented by constant coefficient, finite order differential equations. The non-linearities, coupling and noise effects are also often neglected. Moreover, identification experiment designs which lead to highly complex optimization problems usually require the simultaneous choice of ID algorithm, sensor, and actuator type and placement. On-orbit bandwidth and power restrictions on excitation, limited data window, and restrictions on sensor/actuator type, placement and number, has led to practical questions of implementations.

Hadaegh, Fred Y.↗

Linear operating region in the ozone dial photon counting system

Ozone is a relatively unstable molecule found in Earth's atmosphere. An ozone molecule is made up of three atoms of oxygen. Depending on where ozone resides, it can protect or harm life on Earth. High in the atmosphere, about 15 miles up, ozone acts as a shield to protect Earth's surface from the sun's harmful ultraviolet radiation. Without this shield, we would be more susceptible to skin cancer, cataracts, and impaired immune systems. Closer to Earth, in the air we breathe, ozone is a harmful pollutant that causes damage to lung tissue and plants. Since the early 1980's, airborne lidar systems have been used for making measurements of ozone. The differential absorption lidar (DIAL) technique is used in the remote measurement of O3. This system allows the O3 to be measured as function of the range in the atmosphere. Two frequency-doubled Nd:YAG lasers are used to pump tunable dye lasers. The lasers are operating at 289 nm for the DIAL on-line wavelength of O3, and the other one is operated at 300 nm for the off-line wavelength. The DIAL wavelengths are produced in sequential laser pulses with a time separation of 300 micro s. The backscattered laser energy is collected by telescopes and measured using photon counting systems. The photon counting system measures the light signal by making use of the photon nature of light. The output pulse from the Photo-Multiplier Tube (PE), caused by a photon striking the PMT photo-cathode, is amplified and passed to a pulse height discriminator. The peak value of the pulse is compared to a reference voltage (discrimination level). If the pulse amplitude exceeds the discrimination level, the discriminator generates a standard pulse which is counted by the digital counter. Non-linearity in the system is caused by the overlapping of pulses and the finite response time of the electronics. At low count rates one expects the system to register one event for each output pulse from the PMT corresponding to a photon incident upon the photocathode, however, at higher rates the limitations of the discrimination/counting system will cause the observed count rate to be non-linear with respect to the true count rate. Depending on the pulse height distribution and the discriminator level, the overlapping of pulses (pulse pile-up) can cause count loss or even an additional apparent count gain as the signal levels increase. Characterization of the system, including the pulse height distribution, the signal to noise ratio, and the effect of the discriminator threshold level, is critical in maximizing the linear operating region of the system, thus greatly increasing the useful dynamic range of the system.

Andrawis, Madeleine↗

Orion MPCV Service Module Avionics Ring Pallet Testing, Correlation, and Analysis

The NASA Orion Multi-Purpose Crew Vehicle (MPCV) is being designed to replace the Space Shuttle as the main manned spacecraft for the agency. Based on the predicted environments in the Service Module avionics ring, an isolation system was deemed necessary to protect the avionics packages carried by the spacecraft. Impact, sinusoidal, and random vibration testing were conducted on a prototype Orion Service Module avionics pallet in March 2010 at the NASA Glenn Research Center Structural Dynamics Laboratory (SDL). The pallet design utilized wire rope isolators to reduce the vibration levels seen by the avionics packages. The current pallet design utilizes the same wire rope isolators (M6-120-10) that were tested in March 2010. In an effort to save cost and schedule, the Finite Element Models of the prototype pallet tested in March 2010 were correlated. Frequency Response Function (FRF) comparisons, mode shape and frequency were all part of the correlation process. The non-linear behavior and the modeling the wire rope isolators proved to be the most difficult part of the correlation process. The correlated models of the wire rope isolators were taken from the prototype design and integrated into the current design for future frequency response analysis and component environment specification.

Staab, Lucas↗

STEREO Observations of Interplanetary Coronal Mass Ejections in 2007-2016

We have conducted a survey of 341 interplanetary coronal mass ejections (ICMEs) using STEREO A/B data, analyzing their properties while extending a Level 3 product through 2016. Among the 192 ICMEs with distinguishable sheath region and magnetic obstacle, the magnetic field maxima in the two regions are comparable, and the dynamic pressure peaks mostly in the sheath. The north/south direction of the magnetic field does not present any clear relationship between the sheath region and the magnetic obstacle. About 71% of ICMEs are expanding at 1 au, and their expansion speed varies roughly linearly with their maximum speed except for ICMEs faster than 700 km/s. The total pressure generally peaks near the middle of the well-defined magnetic cloud (MC) passage, while it often declines along with the non-MC ICME passage, consistent with our previous interpretation concerning the effects of sampling geometry on what is observed. The hourly average iron charge state reaches above 12+ approx. 31% of the time for MCs, approx. 16% of the time for non-MC ICMEs, and approx. 1% of the time for non- ICME solar wind. In four ICMEs abrupt deviations of the magnetic field from the nominal field rotations occur in the magnetic obstacles, coincident with a brief drop or increase in field strength-features could be related to the interaction with dust. In comparison with the similar phases of solar cycle 23, the STEREO ICMEs in this cycle occur less often and are generally weaker and slower, although their field and pressure compressions weaken less than the background solar wind.

Jian, L. K.↗

Detection of the large-scale tidal field with galaxy multiplet alignment in the DESI Y1 spectroscopic survey

We explore correlations between the orientations of small galaxy groups, or ‘multiplets’, and the large-scale gravitational tidal field. Using data from the Dark Energy Spectroscopic Instrument (DESI) Y1 survey, we detect the intrinsic alignment (IA) of multiplets to the galaxy-traced matter field out to separations of $100\,h^{-1}$ Mpc. Unlike traditional IA measurements of individual galaxies, this estimator is not limited by imaging of galaxy shapes and allows for direct IA detection beyond redshift $z=1$. Multiplet alignment is a form of higher order clustering, for which the scale-dependence traces the underlying tidal field and amplitude is a result of small-scale ($\lt 1h^{-1}$ Mpc) dynamics. Within samples of bright galaxies, luminous red galaxies (LRG) and emission-line galaxies, we find similar scale-dependence regardless of intrinsic luminosity or colour. This is promising for measuring tidal alignment in galaxy samples that typically display no IA. DESI’s LRG mock galaxy catalogues created from the A BACUS S UMMIT N -body simulations produce a similar alignment signal, though with a 33 per cent lower amplitude at all scales. An analytic model using a non-linear power spectrum (NLA) only matches the signal down to 20 $h^{-1}$ Mpc. Our detection demonstrates that galaxy clustering in the non-linear regime of structure formation preserves an interpretable memory of the large-scale tidal field. Multiplet alignment complements traditional two-point measurements by retaining directional information imprinted by tidal forces, and contains additional line-of-sight information compared to weak lensing. This is a more effective estimator than the alignment of individual galaxies in dense, blue, or faint galaxy samples.

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical Investigation of a Central Fuel Property Hypothesis Under Boosted Spark-Ignition Conditions

In the present work, a central fuel property hypothesis (CFPH), which states that fuel properties are sufficient to provide an indication of a fuel’s performance irrespective of its chemical composition, was numerically investigated. In particular, the objective of the study was to determine whether Research Octane Number (RON) and Motor Octane Number (MON), as fuel properties, are sufficient to describe a fuel’s knock-limited performance under boosted spark-ignition (SI) conditions within the framework of CFPH. To this end, four TPRF-bioblendstock surrogates having different compositions but matched RON (=98) and MON (=90), were first generated using a non-linear regression model based on artificial neural network (ANN). Additionally, three unconventional bioblendstocks were included in the analysis: di-isobutylene (DIB), isobutanol, and Anisole. Skeletal reaction mechanisms were generated for the TPRF-DIB, TPRF-isobutanol, and TPRF-anisole blends from a detailed kinetic mechanism. Thereafter, numerical simulations were performed for the fuel surrogates using the skeletal mechanisms and a virtual cooperative fuel research (CFR) engine model, under a representative boosted operating condition. In the computational fluid dynamics (CFD) model, the G-equation approach was employed to track the turbulent flame front and the well-stirred reactor model combined with the multi-zone binning strategy was used to capture auto-ignition in the end-gas. In addition, laminar flame speed (LFS) was tabulated for each blend as a function of pressure, temperature, and equivalence ratio a priori, and the lookup tables were used to prescribe laminar flame speed as an input to the G-equation model. Parametric spark timing sweeps were performed for each fuel blend to determine the corresponding knock-limited spark advance (KLSA) and 50% burn point (CA50) at the respective KLSA timing. It was observed that despite same RON, MON, and engine operating conditions, the TPRF-anisole blend exhibited markedly different knock-limited performance from the other three blends. This deviation from the octane index (OI) expectation was shown to be caused by differences in laminar flame speed. However, it was found that relatively large fuel-specific differences in LFS (>20%) would have to be present to cause any appreciable deviation from the OI framework. Otherwise, RON and MON would still be robust enough to predict a fuel’s knock-limited performance.

42 ENGINEERING↗

A nonrecursive order N preconditioned conjugate gradient: Range space formulation of MDOF dynamics

While excellent progress has been made in deriving algorithms that are efficient for certain combinations of system topologies and concurrent multiprocessing hardware, several issues must be resolved to incorporate transient simulation in the control design process for large space structures. Specifically, strategies must be developed that are applicable to systems with numerous degrees of freedom. In addition, the algorithms must have a growth potential in that they must also be amenable to implementation on forthcoming parallel system architectures. For mechanical system simulation, this fact implies that algorithms are required that induce parallelism on a fine scale, suitable for the emerging class of highly parallel processors; and transient simulation methods must be automatically load balancing for a wider collection of system topologies and hardware configurations. These problems are addressed by employing a combination range space/preconditioned conjugate gradient formulation of multi-degree-of-freedom dynamics. The method described has several advantages. In a sequential computing environment, the method has the features that: by employing regular ordering of the system connectivity graph, an extremely efficient preconditioner can be derived from the 'range space metric', as opposed to the system coefficient matrix; because of the effectiveness of the preconditioner, preliminary studies indicate that the method can achieve performance rates that depend linearly upon the number of substructures, hence the title 'Order N'; and the method is non-assembling. Furthermore, the approach is promising as a potential parallel processing algorithm in that the method exhibits a fine parallel granularity suitable for a wide collection of combinations of physical system topologies/computer architectures; and the method is easily load balanced among processors, and does not rely upon system topology to induce parallelism.

Kurdila, Andrew J.↗

Solar System Dynamics

In these 18 years, the research has touched every major dynamical problem in the solar system, including: the effect of chaotic zones on the distribution of asteroids, the delivery of meteorites along chaotic pathways, the chaotic motion of Pluto, the chaotic motion of the outer planets and that of the whole solar system, the delivery of short period comets from the Kuiper belt, the tidal evolution of the Uranian arid Galilean satellites, the chaotic tumbling of Hyperion and other irregular satellites, the large chaotic variations of the obliquity of Mars, the evolution of the Earth-Moon system, and the resonant core- mantle dynamics of Earth and Venus. It has introduced new analytical and numerical tools that are in widespread use. Today, nearly every long-term integration of our solar system, its subsystems, and other solar systems uses algorithms that was invented. This research has all been primarily Supported by this sequence of PGG NASA grants. During this period published major investigations of tidal evolution of the Earth-Moon system and of the passage of the Earth and Venus through non-linear core-mantle resonances were completed. It has published a major innovation in symplectic algorithms: the symplectic corrector. A paper was completed on non-perturbative hydrostatic equilibrium.

Wisdom, Jack↗

Gender- and age-related differences in heart rate dynamics: are women more complex than men?

OBJECTIVES. This study aimed to quantify the complex dynamics of beat-to-beat sinus rhythm heart rate fluctuations and to determine their differences as a function of gender and age. BACKGROUND. Recently, measures of heart rate variability and the nonlinear "complexity" of heart rate dynamics have been used as indicators of cardiovascular health. Because women have lower cardiovascular risk and greater longevity than men, we postulated that there are important gender-related differences in beat-to-beat heart rate dynamics. METHODS. We analyzed heart rate dynamics during 8-min segments of continuous electrocardiographic recording in healthy young (20 to 39 years old), middle-aged (40 to 64 years old) and elderly (65 to 90 years old) men (n = 40) and women (n = 27) while they performed spontaneous and metronomic (15 breaths/min) breathing. Relatively high (0.15 to 0.40 Hz) and low (0.01 to 0.15 Hz) frequency components of heart rate variability were computed using spectral analysis. The overall "complexity" of each heart rate time series was quantified by its approximate entropy, a measure of regularity derived from nonlinear dynamics ("chaos" theory). RESULTS. Mean heart rate did not differ between the age groups or genders. High frequency heart rate power and the high/low frequency power ratio decreased with age in both men and women (p < 0.05). The high/low frequency power ratio during spontaneous and metronomic breathing was greater in women than men (p < 0.05). Heart rate approximate entropy decreased with age and was higher in women than men (p < 0.05). CONCLUSIONS. High frequency heart rate spectral power (associated with parasympathetic activity) and the overall complexity of heart rate dynamics are higher in women than men. These complementary findings indicate the need to account for gender-as well as age-related differences in heart rate dynamics. Whether these gender differences are related to lower cardiovascular disease risk and greater longevity in women requires further study.

NASA Discipline Cardiopulmonary↗

Large Eddy Simulation (LES) of Particle-Laden Temporal Mixing Layers

High-fidelity models of plume-regolith interaction are difficult to develop because of the widely disparate flow conditions that exist in this process. The gas in the core of a rocket plume can often be modeled as a time-dependent, high-temperature, turbulent, reacting continuum flow. However, due to the vacuum conditions on the lunar surface, the mean molecular path in the outer parts of the plume is too long for the continuum assumption to remain valid. Molecular methods are better suited to model this region of the flow. Finally, granular and multiphase flow models must be employed to describe the dust and debris that are displaced from the surface, as well as how a crater is formed in the regolith. At present, standard commercial CFD (computational fluid dynamics) software is not capable of coupling each of these flow regimes to provide an accurate representation of this flow process, necessitating the development of custom software. This software solves the fluid-flow-governing equations in an Eulerian framework, coupled with the particle transport equations that are solved in a Lagrangian framework. It uses a fourth-order explicit Runge-Kutta scheme for temporal integration, an eighth-order central finite differencing scheme for spatial discretization. The non-linear terms in the governing equations are recast in cubic skew symmetric form to reduce aliasing error. The second derivative viscous terms are computed using eighth-order narrow stencils that provide better diffusion for the highest resolved wave numbers. A fourth-order Lagrange interpolation procedure is used to obtain gas-phase variable values at the particle locations.

Bellan, Josette↗