Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical approximation & analysis”

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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Radial alignment simulation of solar wind streams observed by Pioneers 10 and 11 in 1974

A particularly favorable radial lineup between spacecraft in Earth orbit and Pioneers 10 and 11 (near the ecliptic plane at approximately 5.5 AU and approximately 4.5 AU, respectively) occurred in mid-1974, when the solar corona was in a stable and well-defined warped-dipole configuration. The radial alignment study reported here differs from previous applications of the technique in two respects: (1) It is the first time a two-dimensional (2-D) MHD model has been tested over such a lengthy propagation interval; the 2-D capability is crucial for treatment of the nonradial shearing motions occurring across the stream interface. (2) The three-dimensional (3-D) structure observed in the white light corona at that time is related to the systematic patterns of nonradial flow deflections appearing at the Pioneer corotating interaction region (CIR) fronts. Comparison of predicted and observed flows for pairs of streams in two successive rotations reveals that when the parent coronal hole projects far across the solar equator in a predominantly north-south orientation (i.e., nearest the 2-D idealization), the mapping is accurate down to details of the flow structures. But where the spacecraft tracks along a latitudinal boundary of a hole or the associated stream front is inclined at a shallow angle to the equator, the numerical projections deviate systematically from the observations. Among the sources of error are 3-D dynamical interactions neglected in the model, differential rotation effects, and slow temporal evolution of the coronal structures. A characteristic pattern of north-south and east-west deflections indicative of the 3-D geometry appears across the CIRs, but these patterns evidently reflect primarily the local, as opposed to global, orientation of the 3-D stream fronts. Such patterns appear common in CIRs observed by Pioneer during this period. These findings thus hold special relevance for the analysis of Ulysses observations, since the present coronal configuration is similar to that of 1974.

Pizzo, V. J.↗
NASA NTRS2002

Formation and Migration of Trans-Neptunian Objects

We consider that trans-Neptunian objects (TNOs) with diameter d greater than 100 kilometers moving now in not very eccentric orbits could be formed directly by the compression of large rarefied dust condensations (with a greater than 30 AU), but not by the accretion of smaller solid planetesimals. Probably, some planetesimals with d approximately 100- 1000 kilometers in the feeding zone of the giant planets and even large main-belt asteroids also could be formed directly by such compression. Some smaller objects (TNOs, planetesimals, asteroids) could be debris of larger objects, and other such objects could be formed directly by compression of condensations. A small portion of planetesimals from the feeding zone of the giant planets that entered into the trans-Neptunian region could left in eccentrical orbits beyond Neptune and became so called "scattered objects". Our computer runs in which gravitational interactions of bodies were taken into account with the use of the spheres method, showed that the embryos of Uranus and Neptune could increase their semimajor axes from less than 10 AU to their present values, moving permanently in orbits with small eccentricities, due to gravitational interactions with the migrating planetesimals. Later on, Thommes et al. considered a similar model using direct numerical integration. The total amount of water delivered to the Earth during the formation of the giant, planets was about the mass of water in the Earth oceans. The end of such bombardment could be caused mainly by the planetesimals that became scattered objects. For the present mass of the trans-Neptunian belt, the collisional lifetime of 1-kilometer TNO is about the age of the Solar system, but only a small portion (less than 1%) of 100-km TNOs could be destroyed during this age. The probability of destruction of a typical TNO (with 30 less than a,less than 50 AU) by scattered objects can be of the same order of magnitude as that by typical TNOs. TNOs could be even more often destroyed during planet formation than during last 4 Gyr. The analysis of the results of the orbital evolution of Jupiter-crossing objects (JCOs) and TNOs showed that, in principle, the trans-Neptunian belt can provide up to 100% of Earth-crossing objects, but, of course, some of them came from the main asteroid belt. Most of the collisions of former JCOs with the Earth were from orbits with aphelia inside Jupiter's orbit and belonged to a small portion of objects which moved in such orbits for a long time. About 1 of 300 JCOs collided with the Sun. The ratio of the total mass of icy planetesimals that migrated from the feeding zone of the giant planets and collided with the planet to the mass of this planet was greater (by a factor of 3 in our runs) for Mars than that for Earth and Venus.

Ipatov, S. I.↗
DOE OSTI2025

Renormalization of states and quasiparticles in many-body downfolding

We explore the principles of many-body Hamiltonian complexity reduction via downfolding on an effective low-dimensional representation. We show that the renormalization factor provides a unique measure of the quality of the compression as it directly represents the projection between the approximate stationary state of the many-body Hamiltonian and the full many-body wavefunction. Hence, the renormalization factor is a measure of fidelity between the effective (reduced-rank) description and the full many-body treatment for arbitrary (i.e., ground and excited) states. When the entire problem is mapped on a system of interacting quasiparticles [Romanova et al., npj Comput. Mater. 9, 126 (2023)], the effective Hamiltonians can faithfully reproduce the physics only when a clear energy scale separation exists between the subsystems and their environment. We also demonstrate that it is necessary to include quasiparticle renormalization at distinct energy scales, capturing the distinct interaction between subsystems and their surrounding environments. Numerical results from simple, exactly solvable models highlight the limitations and strengths of this approach, particularly for ground and low-lying excited states. This work lays the groundwork for applying dynamical downfolding techniques to problems concerned with (quantum) interfaces.

Green-functions technique↗
NASA NTRS2003

Analysis of human postural responses to recoverable falls

We studied the kinematics and kinetics of human postural responses to "recoverable falls." To induce brief falling we used a Hold and Release (H&R) paradigm. Standing subjects actively resisted a force applied to their sternum. When this force was quickly released they were suddenly off balance. For a brief period, approximately 125 ms, until restoring forces were generated to shift the center of foot pressure in front of the center of mass, the body was in a forward fall acted on by gravity and ground support forces. We were able to describe the whole-body postural behavior following release using a multilink inverted pendulum model in a regime of "small oscillations." A three-segment model incorporating upper body, upper leg, and lower leg, with active stiffness and damping at the joints was fully adequate to fit the kinematic data from all conditions. The significance of our findings is that in situations involving recoverable falls or loss of balance the earliest responses are likely dependent on actively-tuned, reflexive mechanisms yielding stiffness and damping modulation of the joints. We demonstrate that haptic cues from index fingertip contact with a stationary surface lead to a significantly smaller angular displacement of the torso and a more rapid recovery of balance. Our H&R paradigm and associated model provide a quantifiable approach to studying recovery from potential falling in normal and clinical subjects.

NASA Discipline Space Human Factors↗
DOE OSTI2024

A chemical kinetic analysis of knock propensity of methanol-to-gasoline fuel

Production of low carbon gasoline-like fuels such as methanol-to-gasoline (MTG) is a promising approach to achieve rapid greenhouse gas emission reduction of the transportation sector. Despite the fact that gasoline that meets the ASTM D4814 standard for automotive spark-ignition engine fuel can be readily produced from these processes, it is unclear how the composition of MTG may affect engine performance and emissions. Here, in this paper, a surrogate for an MTG is used to numerically study the effects of gasoline composition on knock propensity and on the sensitivity of knock to thermal and fuel stratification, to oxygen dilution and to nitric oxide from exhaust gas recirculation of residual gases. Simulations were performed in ANSYS CHEMKIN-PRO using a comprehensive chemical kinetic mechanism for gasoline surrogates, and results of the MTG surrogate were compared against those of a petroleum-based regular E10 gasoline, termed PACE-20. A premium-grade MTG fuel was also formulated by adding ethanol to the MTG surrogate, and results were compared against those of four premium-grade, gasoline-like fuels representative of future alternative gasoline formulations. Surrogates and mechanism were evaluated by comparison against experimental engine data, and the model showed high accuracy at stoichiometric conditions (mean absolute error of ignition timing equal to 1.46 crank angle degrees) but larger deviations at lean conditions (mean absolute error of ignition timing equal to 5.52 crank angle degrees). Despite the fact that the MTG surrogate has a RON 1.1 units higher than that of PACE-20, it may show higher knock propensity at medium temperature conditions due to a less intense NTC behavior. MTG autoignition was more temperature- and equivalence ratio-sensitive than that of PACE20, suggesting that MTG can benefit more from naturally-occurring thermal stratification or from induced fuel stratification of the end gas to mitigate knock intensity. The sensitivity of autoignition reactivity to oxygen dilution and to NO concentration was higher for MTG than for regular gasoline at medium loads, but the opposite trend was observed at high loads due to the effect of pressure on the low-temperature chemistry of regular gasoline. Approximately 14 % vol ethanol content was required to upgrade the octane rating of MTG from regular grade to premium grade. Adding 13.6 % vol ethanol made the fuel autoignition less sensitive to both oxygen dilution and NO content (ignition time varies approx. 17 % and 50 % less with oxygen dilution and NO addition, respectively, when adding ethanol at high engine loads).

02 PETROLEUM↗
NASA NTRS2002

Digital Distortion Caused by Traveling- Wave-Tube Amplifiers Simulated

Future NASA missions demand increased data rates in satellite communications for near real-time transmission of large volumes of remote data. Increased data rates necessitate higher order digital modulation schemes and larger system bandwidth, which place stricter requirements on the allowable distortion caused by the high-power amplifier, or the traveling-wave-tube amplifier (TWTA). In particular, intersymbol interference caused by the TWTA becomes a major consideration for accurate data detection at the receiver. Experimentally investigating the effects of the physical TWTA on intersymbol interference would be prohibitively expensive, as it would require manufacturing numerous amplifiers in addition to acquiring the required digital hardware. Thus, an accurate computational model is essential to predict the effects of the TWTA on system-level performance when a communication system is being designed with adequate digital integrity for high data rates. A fully three-dimensional, time-dependent, TWT interaction model has been developed using the electromagnetic particle-in-cell code MAFIA (Solution of Maxwell's equations by the Finite-Integration-Algorithm). It comprehensively takes into account the effects of frequency-dependent AM (amplitude modulation)/AM and AM/PM (phase modulation) conversion, gain and phase ripple due to reflections, drive-induced oscillations, harmonic generation, intermodulation products, and backward waves. This physics-based TWT model can be used to give a direct description of the effects of the nonlinear TWT on the operational signal as a function of the physical device. Users can define arbitrary excitation functions so that higher order modulated digital signals can be used as input and that computations can directly correlate intersymbol interference with TWT parameters. Standard practice involves using communication-system-level software packages, such as SPW, to predict if adequate signal detection will be achieved. These models use a nonlinear, black-box model to represent the TWTA. The models vary in complexity, but most make several assumptions regarding the operation of the high-power amplifier. When the MAFIA TWT interaction model was used, these assumptions were found to be in significant error. In addition, digital signal performance, including intersymbol interference, was compared using direct data input into the MAFIA model and using the system-level analysis tool SPW for several higher order modulation schemes. Results show significant differences in predicted degradation between SPW and MAFIA simulations, demonstrating the significance of the TWTA approximations made in the SPW model on digital signal performance. For example, a comparison of the SPW and MAFIA output constellation diagrams for a 16-ary quadrature amplitude modulation (16-QAM) signal (data shown only for second and fourth quadrants) is shown. The upper-bound degradation was calculated from the corresponding eye diagrams. In comparison to SPW simulations, the MAFIA data resulted in a 3.6-dB larger degradation.

Kory, Carol L.↗
DOE OSTI2020

Pass-efficient methods for compression of high-dimensional turbulent flow data

The future of high-performance computing, specifically on future Exascale computers, will presumably see memory capacity and bandwidth fail to keep pace with data generated, for instance, from massively parallel partial differential equation (PDE) systems. Current strategies proposed to address this bottleneck entail the omission of large fractions of data, as well as the incorporation of in situ compression algorithms to avoid overuse of memory. To ensure that post-processing operations are successful, this must be done in a way that a sufficiently accurate representation of the solution is stored. Moreover, in situations where the input/output system becomes a bottleneck in analysis, visualization, etc., or the execution of the PDE solver is expensive, the number of passes made over the data must be minimized. In the interest of addressing this problem, this work focuses on the utility of pass-efficient, parallelizable, low-rank, matrix decomposition methods in compressing high-dimensional simulation data from turbulent flows. Additionally, a particular emphasis is placed on using coarse representation of the data – compatible with the PDE discretization grid – to accelerate the construction of the low-rank factorization. This includes the presentation of a novel single-pass matrix decomposition algorithm for computing the so-called interpolative decomposition. The methods are described extensively and numerical experiments on two turbulent channel flow data are performed. In the first (unladen) channel flow case, compression factors exceeding 400 are achieved while maintaining accuracy with respect to first- and second-order flow statistics. In the particle-laden case, compression factors of 100 are achieved and the compressed data is used to recover particle velocities. These results show that these compression methods can enable efficient computation of various quantities of interest in both the carrier and disperse phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

Venus small volcano classification and description

The high resolution and global coverage of the Magellan radar image data set allows detailed study of the smallest volcanoes on the planet. A modified classification scheme for volcanoes less than 20 km in diameter is shown and described. It is based on observations of all members of the 556 significant clusters or fields of small volcanoes located and described by this author during data collection for the Magellan Volcanic and Magmatic Feature Catalog. This global study of approximately 10 exp 4 volcanoes provides new information for refining small volcano classification based on individual characteristics. Total number of these volcanoes was estimated to be 10 exp 5 to 10 exp 6 planetwide based on pre-Magellan analysis of Venera 15/16, and during preparation of the global catalog, small volcanoes were identified individually or in clusters in every C1-MIDR mosaic of the Magellan data set. Basal diameter (based on 1000 measured edifices) generally ranges from 2 to 12 km with a mode of 34 km, and follows an exponential distribution similar to the size frequency distribution of seamounts as measured from GLORIA sonar images. This is a typical distribution for most size-limited natural phenomena unlike impact craters which follow a power law distribution and continue to infinitely increase in number with decreasing size. Using an exponential distribution calculated from measured small volcanoes selected globally at random, we can calculate total number possible given a minimum size. The paucity of edifice diameters less than 2 km may be due to inability to identify very small volcanic edifices in this data set; however, summit pits are recognizable at smaller diameters, and 2 km may represent a significant minimum diameter related to style of volcanic eruption. Guest, et al, discussed four general types of small volcanic edifices on Venus: (1) small lava shields; (2) small volcanic cones; (3) small volcanic domes; and (4) scalloped margin domes ('ticks'). Steep-sided domes or 'pancake domes', larger than 20 km in diameter, were included with the small volcanic domes. For the purposes of this study, only volcanic edifices less than 20 km in diameter are discussed. This forms a convenient cutoff since most of the steep-sided domes ('pancake domes') and scalloped margin domes ('ticks') are 20 to 100 km in diameter, are much less numerous globally than are the smaller diameter volcanic edifices (2 to 3 orders of magnitude lower in total global number), and do not commonly occur in large clusters or fields of large numbers of edifices.

Aubele, J. C.↗
NASA NTRS1998

Potential Mars Surveyor 2001 Landing Sites: Low-Elevation Cratered "Highlands" in Central and Eastern Sinus Meridiani and Near Amenthes Fossae

The main scientific goal for the Mars Surveyor Program 2001 (MSP 01) landed mission is to collect and characterize 91 rock and 13 soil core samples using an integrated instrument suite onboard the Athena Rover. If possible, these samples will be retrieved and returned to Earth via a MSP 05 or MSP 07 mission. Preliminary engineering constraints for the MSP 01 landing site call for a location that lies between 15 S and 30 N, and below about 2 km elevation (based on Viking-era topography). Desirable landing sites for MSP 01 are to be located "in the ancient highlands where the environmental conditions may have been favorable to the preservation of evidence of possible prebiotic or biotic processes including the emergence (and, potentially, the persistence) of life". We interpret this to mean that the desirable sites include those that have evidence of aqueous sediments that might have been deposited during the Noachian and/or Hesperian Epochs of Mars' history. In addition to the search for subaqueous sedimentary deposits, we took into consideration the fact that the rover, Athena, will need to be able to access these materials. Thus, a site where aeolian deflation has occurred might be desirable because it might expose, in situ, layered sedimentary deposits. Deflated areas, of course, might include potential landing hazards in the form of meterscale buttes and mesas (e.g., Christmas Lake Valley, OR), thus careful study of such sites with high resolution images will be required before a decision is made to land. We have been examining three regions that have potential to be considered for MSP 01 landing sites. This work is based on Viking (VIS, IRTM) and Phobos 2 (Termoskan) observations and should be regarded as preliminary because we believe that the final site selection should also be based upon analysis of Mars Global Surveyor observations that help constrain mineralogy (TES) and local geomorphology (MOC, MOLA). (1) Eastern Sinus Meridiani Region (proposed by K. S. Edgett) Sinus Meridiani is a persistent low-albedo (< 0. 16) region on the martian equator that has been recognized for about 400 years. All of Sinus Meridiani is below the 2 km elevation constraint for MSP 01. The region includes cratered highlands, valley networks, aeolian dunes, and possible aqueous sedimentary deposits. The landing site study region is located in the eastern portion of Sinus Meridiani. It is bounded by latitudes 10 S to 2 N, longitudes 355 W to 345 W, and the elevations are mostly 1-2 km. The center of this area contains a medium-albedo (0.19-0.21), relatively smooth-surfaced deposit that was suggested by Rice to be a lacustrine deposit. Several potential landing areas can be suggested within this region. Based on Viking and Phobos 2 data, the favored sites so far are centered at 7.6 S, 346.9 W and 0.8 S, 349 W. Thermal inertias are 3.2-7.0 x 10(exp -3) cal /sq cm s(exp -0.5)/K; and rock abundances are around 2-6%. The site at 7.6 S, 346.9 W is at the southern end of the smooth, medium-albedo unit that might have a lacustrine origin. At this location, numerous channels appear to have drained toward the smooth unit. Viking images from orbit 747A show this area at about 15 m/pixel ground resolution. The images reveal that aeolian deflation has occurred along the deposit's margins. Bright (i.e., albedo >= 0.21) aeolian dunes are present on the channel floors and in some of the depressions on the smooth unit. The bright, apparently active dunes might consist of material (perhaps lakedeposited sands) that has been eroded from the smooth unit. The site at approximately 0.8 S, 349 W is selected because it offers an opportunity to solve a long-standing puzzle about Mars remote sensing. There are three main "color" units on Mars: "dark red, dark gray, and bright red". This landing site would allow the Athena rover an opportunity to investigate all three materials within close proximity (the best place on Mars to do so). There are no high resolution (better than 100 m/pixel) Viking or Mariner images of this site. (2) Central Sinus Meridiani Region (proposed by K. S. Edgett and T. J. Parker) Central Sinus Meridiani is characterized by two types of surfaces [4]. One is like typical martian cratered highlands elsewhere- there are old valley networks and old impact craters. The other is relatively smooth and flat. These two units are in contact around 3.1 S between 5 W and 4 E longitudes. Valley networks- including one at 6'S, 358 W that rivals the Grand Canyon of Arizona- once drained toward the smooth unit. Edgett and Parker [4] proposed that the smooth unit might consist of sediments laid down in a large Noachian-aged sea/ocean that would have covered much of the northern hemisphere. Schultz and Lutz [I I I suggested that it is a paleopolar layered deposit. Regardless, the smooth unit where it contacts the cratered terrain would make an excellent site for Athena rover to investigate. The site is best seen in Viking high resolution images from orbits 408B (about 30 m/pixel) and 746A (about 12 m/pixel). These images suggest that aeolian deflation has occurred along the margin of the smooth unit, and this deflation has exposed horizontal larrs of material. The elevation is about 0.5 km; thermal inertias are 6.5-8.0 x 10(exp -3) cal /sq cm s(exp -0.5) / K; rock abundances are 2-4%; and the surface is probably sandy with dark drifts and ripples but almost no actual dunes. We suggest a landing around 3.2 S, 3.0 W would test the aqueous sediment hypothesis and provide a potentially smooth surface on which to land. (3) Amenthes Fossae Region (proposed by S. N. Huntwork and K. S. Edgett) The Amenthes Fossae are a series of graben/fissures that are circumferential to the southeast side of Isidis Planitia. These fissures cross a variety of ancient, heavily cratered Noachian terrain and younger, Hesperian and Amazonian terrain. We focused our search on a region 0-15 N, 250 - 270 W. Elevations are -0.5 to 2 km. Depending upon whether Isidis Planitia was ever a water-rich environment, this region might have been influenced by aqueous sedimentation. Valley networks are common, and they drained toward the north and northwest. We focused our work on a set of Viking orbiter high-resolution images, 719A 1-48. These have resolutions 16-24 rri/pixel. We examined images 20-23, centered at 2 N, 258 W. This site, on the plains just southwest of a 42 km-diameter crater, includes a valley network channel, a relatively young crater ejecta deposit, a few buttes composed of presumably ancient, Noachian bedrock, and a "plains" unit. The plains might make an ideal landing surface, except for the presence of some fine-scale ridges (oriented approximately N-S). The ridges are probably yardangs, thus this site offers a place where aeolian deflation has probably exposed some of the layered rock units that comprise the "plains". Thermal inertias are 7.9-8.3 x 10(exp -3) cal/sq cm s(exp -0.5) /K and rock abundances are 10-15%. A rover traverse might include the opportunity to go down to the floor (and sample along the walls) of the valley network channel at 2 N 258.1 W.

Edgett, K. S.↗
DOE OSTI2020

Electro-osmo-dialysis through nanoporous layers physically conjugated to micro-perforated ion-exchange membranes: Highly selective accumulation of trace coions

This study explores theoretically and numerically the phenomenon of counter-flow accumulation of trace coions during electro-osmo-dialysis through composite membranes containing nanoporous and perforated ion-exchange layers. The combination of charged nanopores and perforations allows electroosmosis, and ion concentration polarization appears near the ion-exchange layer. Accumulation of trace coions occurs because their advective and electromigration flow velocity components have opposite signs and are locally compensating in depletion regions where electric fields are high. This local compensation of flow components can give rise to high accumulation factors (>1000 in some cases) for trace coions at some distance from the interface between the nanoporous and ion-exchange layers. Moreover, the extent of this accumulation depends on the trace-coion diffusion coefficient. The resulting selectivity between different coions might serve as a basis for environmentally friendly separation of ionic species such as synthetic peptides or rare-earth elements. For example, a periodic process might include accumulation followed by species release into smaller volumes of solution. For both rare-earth ions and peptide-like species, the accumulation selectivity between ions with small differences in mobilities depends on system parameters such as nanoporous-layer thickness, distance between the perforations, and the zeta potential of the nanopore surface. Here, the use of optimal values of these membrane parameters may enable tuning of the electro-osmo-dialysis system for separation of specific ions. Modelling reveals that despite inhomogeneous patterns of fluid and ion flows close to the nanoporous/perforated ion-exchange interface, at relatively short distances from this interface (less than half the distance between the perforations) all the flows become practically one-dimensional (1D). Similarly, the concentration profiles of trace coions also become 1D. The flow one-dimensionality at short distances from the current-polarized interface enables a simple 1D approximation that is useful for an approximate analysis of trends in this multi-parameter system.

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

Preliminary Planar Formation: Flight Dynamics Near Sun-Earth L2 Point

NASA's Goddard Space Flight Center is planning a series of missions in the vicinity of the Sun-Earth L2 libration point. Some of these projects will involve a distributed space system of telescope spacecraft acting together as a single telescope for high-resolution. The individual telescopes will be configured in a plane, surrounding a hub, where the telescope plane can be aimed toward various astronomical targets of interest. In preparation for these missions, it is necessary to develop an improved understanding of the dynamical behavior of objects in a planar configuration near L2. The classical circular restricted three body problem is taken as the basis for the analysis. At first order, the motion of such a telescope relative to the hub is described by a system of linear second order differential equations. These equations are identical to the circular restricted problem's linear equations describing the hub motion about L2. Therefore, the fundamental frequencies, both parallel to and normal to the ecliptic plane, are the same for the relative telescope motion as for the hub motion. To maintain the telescope plane for the duration necessary for the planned observations, a halo-type orbit of the telescopes about the hub is investigated. By using a halo orbit, the individual telescopes remain in approximately the same plane over the observation duration. For such an orbit, the fundamental periods parallel to and normal to the ecliptic plane are forced to be the same by careful selection of the initial conditions in order to adjust the higher order forces. The relative amplitudes of the resulting oscillations are associated with the orientation of the telescope plane relative to the ecliptic. As in the circular restricted problem, initial conditions for the linearized equations must be selected so as not to excite the convergent or divergent linear modes. In a higher order analysis, the telescope relative motion equations include the effects of the position of the hub relative to L2. In this paper, the differential equations are developed through second order in the distance of the hub from the libration point. A modified Lindstedt-Poincad perturbation method is employed to construct the solution of these differential equations through that same order of magnitude. In the course of the solution process, relationships are determined between the initial conditions of the telescopes, selected in order to avoid resonance excitation. As the differential equations include the hub position, it is necessary to simultaneously develop the solution for the hub. As has been done in past analyses of the circular restricted problem, the hub position is written in a power series formulation in terms of its distance from L2. Then, in order to be included in the telescope equations, the hub solution is cast in terms of the nonlinear frequency of the relative telescope motion. In the course of the analysis, it is determined that the hub should also maintain a halo orbit - about L2. Additionally, relationships are formed between the initial conditions of the telescopes and the hub. These relationships may be used to associate sets of initial conditions with particular orientations of the telescope plane. The accuracy of the analytical solution is verified through various simulations and comparison to numerical integration of the differential equations. The results of the simulations are presented, along with a graphical representation of the relationships between the initial conditions of the telescopes and hub.

Segerman, Alan M.↗
NASA NTRS2003

Small Particle Response to Fluid Motion Using Tethered Particles to Simulate Microgravity

This paper reports on ground based work conducted to support the Spaceflight Definition project SHIVA (Spaceflight Holography Investigation in a Virtual Apparatus). SHIVA will advance our understanding of the movement of a particle in a fluid. Gravity usually dominates the equations of motion, but in microgravity as well as on earth other terms can become important. Through an innovative application of fractional differential equations, two members of our team produced the first analytical solution of a fundamental equation of motion, which had only been solved numerically or by approximation before. The general solution predicts that the usually neglected history term becomes important in particle response to a sinusoidal fluid movement when the characteristic viscous time is in the same order as the fluid oscillation period and peaks when the two times are equal. In this case three force terms, the Stokes drag, the added mass, and the history drag must all be included in predicting particle movement. We have developed diagnostic recording methods using holography to save all of the particle field data, allowing the experiment to essentially be transferred from space back to earth in what we call the virtual apparatus for on-earth microgravity experimentation. We can quantify precisely the three-dimensional motion of sets of particles, allowing us to test and apply the new analytical solutions. We are examining the response of particles up to 2 mm radius to fluid oscillation at frequencies up to 80 Hz with amplitudes up to 200 microns. Ground studies to support the flight development program have employed various schemes to simulate microgravity. One of the most reliable and meaningful methods uses spheres tethered to a fine hair suspended in the fluid. We have also investigated particles with nearly neutral buoyancy. Recordings are made at the peak amplitudes of vibration of the cell providing a measure of the ratio of fluid to particle amplitude. The experiment requires precise location of the particle to within microns during recording, and techniques for achieving this are one of the project challenges. Focused microscopic images and diffraction patterns are used. To make the experiment more versatile, the spaceflight system will record holograms both on film and electronically. A cross correlation procedure enables sub pixel accuracies for electronic recordings, partially accommodating the lower spatial resolution of CCDs. The electronic holograms can be down linked providing real time data. Results of the ground experiments, the flight experiment design, and data analysis procedures are reported.

Trolinger, James↗
DOE OSTI2023

Shock-induced transformation of nitinol shape memory alloy: Effect of stress state on transformation

Due to its numerous practical applications and intriguing phase transformation behavior, shape memory alloys (SMAs) have garnered significant research and development interests. In the past, most studies on the mechanical behavior of SMAs have been conducted under uniaxial stress loadings. Limited research on SMAs under shock loading has not provided conclusive results regarding their transformation behavior and transformation stress under such loading. Additionally, there is a lack of comprehensive understanding regarding the effects of different stress states on transformation behavior. The main objectives of this study are to address these issues. To achieve these objectives, a series of shock wave experiments were designed and conducted. Additionally, quasi-static and dynamic uniaxial stress experiments were carried out to establish a baseline for comparison. The results revealed that the transformation stress under dynamic uniaxial strain shock loading was approximately 1.92 GPa in contrast to 0.5 GPa (quasi-static) to 0.8 GPa (dynamic) observed in uniaxial stress loading. The transformation behavior exhibited noticeable rate sensitivity for both types of loading. There appeared to be a critical strain rate above which the austenite phase was driven to a metastable state. This estimated critical axial strain rate along the loading direction was approximately 2 × 10 3 /s–4 × 10 3 /s for uniaxial stress loading and approximately 2 × 10 6 /s for uniaxial strain loading. The apparent high transformation stress for uniaxial strain loading can likely be attributed to a combination of high-pressure confinement and high strain rate. Furthermore, determining their relative contributions remains an open issue.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Advances in constraining intrinsic alignment models with hydrodynamic simulations

We use galaxies from the illustristng, massiveblack-ii, and illustris-1 hydrodynamic simulations to investigate the behaviour of large scale galaxy intrinsic alignments. Our analysis spans four redshift slices over the approximate range of contemporary lensing surveys z = 0–1. We construct comparable weighted samples from the three simulations, which we then analyse using an alignment model that includes both linear and quadratic alignment contributions. Our data vector includes galaxy–galaxy, galaxy–shape, and shape–shape projected correlations, with the joint covariance matrix estimated analytically. In all of the simulations, we report non-zero IAs at the level of several σ. For a fixed lower mass threshold, we find a relatively strong redshift dependence in all three simulations, with the linear IA amplitude increasing by a factor of ~2 between redshifts z = 0 and z = 1. We report no significant evidence for non-zero values of the tidal torquing amplitude, A 2 , in TNG, above statistical uncertainties, although MBII favours a moderately negative A 2 ~ –2. Examining the properties of the TATT model as a function of colour, luminosity and galaxy type (satellite or central), our findings are consistent with the most recent measurements on real data. We also outline a novel method for constraining the TATT model parameters directly from the pixelized tidal field, alongside a proof-of-concept exercise using TNG. This technique is shown to be promising, although comparison with previous results obtained via other methods is non-trivial.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

An infrasound source analysis of the OSIRIS-REx sample return capsule hypersonic re-entry

The OSIRIS-REx sample return capsule's hypersonic re-entry into the atmosphere is a rare opportunity to test a variety of sonic boom source models since the projectile dimensions are well characterized. While the as-flown flight path is unknown, the predicted flight path enables a rough approximation of the source Mach number and location. Six infrasound microphones deployed in the boom carpet along the predicted flight path recorded impulsive signals from the OSIRIS-REx re-entry. Using a suite of atmosphere profiles and the geometric acoustics approximation, we estimate locations with uncertainty estimates along the flight path from which the signals were emitted. Acoustic overpressure and signal duration predictions from Whitham's far field theory, Carlson's simplified sonic boom prediction method, and a drag-dominated hypersonic model are analyzed with uncertainty estimates from the location estimate. While the Carlson simplified sonic boom prediction method could be accurate, our preference is for the drag-dominated source model. Using this source model with an inviscid Burgers's equation solver for propagation, we obtained an excellent match to the recorded data. In conclusion, these results will help better inform future sample return capsule re-entry observation campaigns as well as contribute to a better understanding of high altitude infrasonic sources.

58 GEOSCIENCES↗
DOE OSTI2021

Analysis and prediction of reaction kinetics using the degree of rate control

“Degree of rate control” (DRC) analysis provides a quantitative approach for analysing the kinetics of multi-step reaction mechanisms that has been widely applied to both heterogeneous and homogeneous catalysis research, as well as electrocatalysis. The DRC of any given transition state or intermediate is defined as a partial derivative such that it approximately equals the fractional increase in net rate to the product of interest per differential decrease in its standard-state free energy for that species (÷RT), holding constant the standard-state free energies of all other transition states and intermediates. Even very complex mechanisms usually have only a few species with non-zero DRCs and are thus the species whose interactions with the catalyst most strongly affect the net rate. These key DRC values thus offer a simple and intuitive route to optimize catalyst materials, especially with the assistance of computational methods like density functional theory (DFT). These high-DRC species are also the species whose energetics must be most accurately measured or calculated to achieve an accurate kinetic model for any reaction mechanism. In simple cases with a single “rate-determining step”, the DRC for its transition state (TS) is + 1. Catalyst-bound intermediates, on the other hand, often have negative DRCs equal to a small integer times their fractional occupancy of catalyst sites. The apparent activation energy equals a weighted average of the standard-state enthalpies (relative to reactants) of all the species (intermediates, transition states and products) in the reaction mechanism, each weighted by its DRC (+RT). It has been shown that the apparent transfer coefficient in electrocatalysis, an inverted form of the Tafel slope, is a weighted average of the number of electrons transferred to generate each intermediate or product species in the mechanism, weighted again by the DRC. Quantitative analysis of kinetic isotope effects (KIEs, or the ratio of net rates for different reactant isotopes) in complex mechanisms has shown that the logarithm of the KIE equals the weighted average over all species in the mechanism of the difference between the two isotopes in their standard-state free energies (÷RT), again weighted by the DRC. The reaction orders with respect to fluid-phase concentrations of reactants, products and intermediates have also been proven to be directly related to DRCs. Thus, there are numerous experimental observables which equate to short linear combinations of DRCs, so that combinations of experimental measurements might provide access to DRC values. Since its invention for transition states in 1994 and its generalization to include intermediates in 2009, the DRC has thus far mainly been calculated for microkinetic models based either entirely on DFT or on DFT with the key energies (i.e., those for high-DRC species) being fine-tuned to match experiments. The relationships summarized in this work provide new opportunities for using experiments earlier in the development and optimization of microkinetic models that require input from computational catalysis.

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

Status of the NASA GMAO Observing System Simulation Experiment

An Observing System Simulation Experiment (OSSE) is a pure modeling study used when actual observations are too expensive or difficult to obtain. OSSEs are valuable tools for determining the potential impact of new observing systems on numerical weather forecasts and for evaluation of data assimilation systems (DAS). An OSSE has been developed at the NASA Global Modeling and Assimilation Office (GMAO, Errico et al 2013). The GMAO OSSE uses a 13-month integration of the European Centre for Medium- Range Weather Forecasts 2005 operational model at T511/L91 resolution for the Nature Run (NR). Synthetic observations have been updated so that they are based on real observations during the summer of 2013. The emulated observation types include AMSU-A, MHS, IASI, AIRS, and HIRS4 radiance data, GPS-RO, and conventional types including aircraft, rawinsonde, profiler, surface, and satellite winds. The synthetic satellite wind observations are colocated with the NR cloud fields, and the rawinsondes are advected during ascent using the NR wind fields. Data counts for the synthetic observations are matched as closely as possible to real data counts, as shown in Figure 2. Errors are added to the synthetic observations to emulate representativeness and instrument errors. The synthetic errors are calibrated so that the statistics of observation innovation and analysis increments in the OSSE are similar to the same statistics for assimilation of real observations, in an iterative method described by Errico et al (2013). The standard deviations of observation minus forecast (xo-H(xb)) are compared for the OSSE and real data in Figure 3. The synthetic errors include both random, uncorrelated errors, and an additional correlated error component for some observational types. Vertically correlated errors are included for conventional sounding data and GPS-RO, and channel correlated errors are introduced to AIRS and IASI (Figure 4). HIRS, AMSU-A, and MHS have a component of horizontally correlated error. The forecast model used by the GMAO OSSE is the Goddard Earth Observing System Model, Version 5 (GEOS-5) with Gridpoint Statistical Interpolation (GSI) DAS. The model version has been updated to v. 5.13.3, corresponding to the current operational model. Forecasts are run on a cube-sphere grid with 180 points along each edge of the cube (approximately 0.5 degree horizontal resolution) with 72 vertical levels. The DAS is cycled at 6-hour intervals, with 240 hour forecasts launched daily at 0000 UTC. Evaluation of the forecasting skill for July and August is currently underway. Prior versions of the GMAO OSSE have been found to have greater forecasting skill than real world forecasts. It is anticipated that similar forecast skill will be found in the updated OSSE.

OSSE↗
DOE OSTI2021

Hydro-chemo-mechanical phase field formulation for corrosion induced cracking in reinforced concrete

Steel corrosion induced cracking in reinforced concrete structures is often caused by volumetric straining at the steel-concrete interface. This volumetric strain is mainly controlled by the rate of generation of corrosion product, which depends on the diffusion of Cl{sup −} and dissolved O{sub 2} as well as the chemical reaction at the steel-concrete interface. While many numerical modeling techniques have been proposed to tackle this complex phenomenon, most have been limited to using staggered solution schemes, where the diffusion, chemical reaction, mechanical strain, and crack fields are solved for separately. Unfortunately, such approaches introduce numerical errors that can lead to unrealistic predictions. Therefore, a coupled solution, is needed to remove the formulation induced errors and understand this complex phenomenon. To provide a unifying solution framework for corrosion induced cracking in reinforced concrete, a phase field formulation is presented and implemented with FEM algorithms. The governing equations are derived from the system's free energy and Fick's law of diffusion. The proposed formulation was first validated in 1D by comparing the modeling results with the semi-analytical solution approximated via Laurent series to avoid singularities in the complex domain. 2D plane strain solutions to the embedded rebar undergoing corrosion induced cracking with pre-cracks are also presented and compared with the approximated solutions from linear elastic fracture mechanics theory (LEFM). The modeling of crack propagation enabled by the proposed formulation is also presented along with parametric studies to reveal the roles of the fracture toughness and permeability of the concrete on the corrosion induced cracking. A case study was also conducted on H-pile steel with concrete jackets to demonstrate the feasibility of modeling corrosion induced fracture in reinforced concrete structures with complex steel-concrete interfaces.

36 MATERIALS SCIENCE↗