Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical phase transitions”

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 109 records · Page 6

Solid-liquid phase transition in argon

Starting from the Lennard-Jones interatomic potential, a modified cell theory has been used to describe the solid-liquid phase transition in argon. The cell-size variations may be evaluated by a self-consistent condition. With the inclusion of cell-size variations, the transition temperature, the solid and liquid densities, and the liquid-phase radial-distribution functions have been calculated. These ab initio results are in satisfactory agreement with molecular-dynamics calculations as well as experimental data on argon.

Tsang, T.↗

Liftoff and Transition Database Generation for Launch Vehicles Using Data-Fusion-Based Modeling

A data fusion technique for merging multiple data sources with differing fidelity and resolution was developed to support the production of aerodynamic line load databases for the Liftoff and Transition (LOT) flight phase of the Space Launch System (SLS). The technique uses a reduced order model based on a high-fidelity line load data set from Computational Fluid Dynamics (CFD) to predict solutions for a much larger solution space. Even higher-fidelity force and moment information (from wind-tunnel tests) is then used to adjust the model. The adjustment uses constrained optimization through the method of Lagrange multipliers in order to minimize the deviation of the line load distribution from the spatially-dense CFD solution, while ensuring that the integrated force and moment values match those observed in physical wind tunnel measurements. Though the wind-tunnel data are operationally-dense (available at many flow conditions), they are spatially coarse (as only the overall forces and moments are available). Conversely, CFD for such complex configurations is expensive, and thus operationally sparse. Data fusion techniques are necessary to make the most efficient use of available information, delivering accurate results within time and resource constraints.

Wignall, T. J.↗

Statistics and dynamics of persistent anomalies

Persistent anomalies with recurrent spatial patterns play an important role in the atmosphere's low-frequency variability. A connection between statistical and dynamical methods of description and prediction of persistent anomalies is established by computing and analyzing the empirical orthogonal functions (EOFs) in a simple deterministic model, on the one hand, and in Southern Hemisphere geopotential heights, on the other. The dynamical model is governed by the fully nonlinear, equivalent-barotropic vorticity equation on the sphere, with simplified forcing, dissipation and topography. The Southern Hemisphere data consist in gridded daily maps of 500 mb heights from June 1972 to July 1983. Two types of persistent anomalies appear in this time series, both having a strong wavenumber-three component; they differ by the value of the constant phase of this wave and by the strength of the wavenumber-one component. The first two EOFs bear a striking resemblance to these two patterns. It is concluded that the dynamical interpretation of EOFs is their pointing from the time mean to the most populated regions of the system's phase space. Pursuing this interpretation, a Markov-chain formulation of transitions from one persistent anomaly regime to another is introduced, and the implications for long-range forecasting are discussed.

Mo, Kingtse C.↗

Access to Mars from Earth-Moon Libration Point Orbits:

This investigation is focused specifically on transfers from Earth-Moon L(sub 1)/L(sub 2) libration point orbits to Mars. Initially, the analysis is based in the circular restricted three-body problem to utilize the framework of the invariant manifolds. Various departure scenarios are compared, including arcs that leverage manifolds associated with the Sun-Earth L(sub 2) orbits as well as non-manifold trajectories. For the manifold options, ballistic transfers from Earth-Moon L(sub 2) libration point orbits to Sun-Earth L(sub 1)/L(sub 2) halo orbits are first computed. This autonomous procedure applies to both departure and arrival between the Earth-Moon and Sun-Earth systems. Departure times in the lunar cycle, amplitudes and types of libration point orbits, manifold selection, and the orientation/location of the surface of section all contribute to produce a variety of options. As the destination planet, the ephemeris position for Mars is employed throughout the analysis. The complete transfer is transitioned to the ephemeris model after the initial design phase. Results for multiple departure/arrival scenarios are compared.

Dynamical Systems↗

A Missed Thrust Framework for Low-Thrust Spiral Trajectories to the NRHO

A framework is developed by which end-to-end optimization of many-revolution low-thrust spiral trajectories can be completed in the presence of missed thrust events. This framework is applied to the Lunar Transit trajectory by which the initial capability of NASA’s Gateway lunar space station will be delivered to a Near Rectilinear Halo Orbit. This low-thrust mission consists of three subphases, each designed according to the specific objectives and dynamical regimes encountered as the mission progresses from a medium Earth insertion orbit to cislunar space. The presented framework accounts for the unique considerations demanded by each mission phase and incorporates appropriate capabilities into a novel mission analysis tool. This methodology enables large scale and reliable analyses of missed thrust events across the end-to-end Lunar Transit to verify the robustness of flight trajectories across the full range of considered launch dates.

missed thrust↗

Computational Fluid Dynamics Methods Used in the Development of the Space Launch System Liftoff and Transition Lineloads Databases

The objective of this paper is to document the reasoning and trade studies that supported the selection of appropriate tools for constructing aerodynamic lineload databases for the Liftoff and Transition phases of flight for launch vehicles. These decisions were made amid the maturation of an evolving workflow for generating databases on variants of the Space Launch System launch vehicle, with most being based on results from brief developmental studies performed in response to specific, unforeseen challenges that were encountered in analyzing a given configuration. This report is intended to provide a summary of the results and the decision-making processes chronologically over the design cycles of various configurations, starting with isolated free-air bodies for the Block 1 Crew, then the Block 1B Crew and Cargo configurations, and most recently the Block 1B Crew configuration in proximity to the launch tower. The results from these analyses led to the selection of the CREATE-AV Kestrel flowsolver for simulating these problems. The need to accurately capture the expected leeward-wake flow field characteristics required the use of Delayed Detached Eddy Simulation (DDES) method, for which the vorticity magnitude was employed as the solution Adaptive Mesh Refinement (AMR) function over the off-body Cartesian grid region. In addition, the Spalart-Allmaras (SA) model is used to account for the flow turbulence effects.

Ratnayake, Nalin A.↗

Critical Slowing Down in Time-to-Extinction: An Example of Critical Phenomena in Ecology

We study a model for two competing species that explicitly accounts for effects due to discreteness, stochasticity and spatial extension of populations. The two species are equally preferred by the environment and do better when surrounded by others of the same species. We observe that the final outcome depends on the initial densities (uniformly distributed in space) of the two species. The observed phase transition is a continuous one and key macroscopic quantities like the correlation length of clusters and the time-to-extinction diverge at a critical point. Away from the critical point, the dynamics can be described by a mean-field approximation. Close to the critical point, however, there is a crossover to power-law behavior because of the gross mismatch between the largest and smallest scales in the system. We have developed a theory based on surface effects, which is in good agreement with the observed behavior. The course-grained reaction-diffusion system obtained from the mean-field dynamics agrees well with the particle system.

Gandhi, Amar↗

Architecture and Dynamics of Kepler's Multi-Transiting Planet Systems.III. Comprehensive Investigation Using All Four Years of Kepler Mission Data

We study the orbital architectures of planetary systems orbiting within ∼1 AU of their stars by analyzing the ensemble of Kepler systems having two or more planet candidates. We use data from the entire Kepler mission, and in many cases we apply improved analysis techniques (e.g., replacing histograms by top-hat Kernel Density Estimators that avoid the loss of information resulting from choosing a particular phase for the bin boundaries) to extend and enhance the studies of Lissauer et al. (2011, ApJS 197, 8) and Fabrycky et al. (2014, ApJ 790, 146).These data show ~ 1700 transiting planet candidates in > 600 multiple-planet systems, far more than were available for our previous two studies. The increased numbers and better information about planetary radii and the properties of stellar hosts made possible by Gaia DR2 allow more statistically-robust analyses of the entire ensemble of Kepler multis as well as independent analyses of subsets of the population. We are thus able to contrast the dynamical configurations of small and large planets, short-period and longer-period planets, and planets orbiting various types of host stars. We reinforce our previous findings that most pairs of planets within the same system are neither in nor near low-order mean motion resonances and that there is a substantial excess of planets having period ratios slightly larger than those of first-order mean-motion resonances. However, neglecting three systems whose planets are locked in 3- body resonances and summing over all first-order mean motion resonances, the deficit of planet pairs with period ratios just narrow of resonance is as large as the excess of planets wide of resonance (within statistical uncertainties), suggesting that overall there is no overall excess of planet pairs in the vicinity of resonance. Other aspects of our study, including estimates of the typical relative inclinations of planetary orbits and their variations as functions of orbital period, planet sizes and stellar properties, are in progress, with results expected to be available for presentation by the time of the conference.

Lissauer, Jack↗

Better, Cheaper, Faster Molecular Dynamics

Recent, revolutionary progress in genomics and structural, molecular and cellular biology has created new opportunities for molecular-level computer simulations of biological systems by providing vast amounts of data that require interpretation. These opportunities are further enhanced by the increasing availability of massively parallel computers. For many problems, the method of choice is classical molecular dynamics (iterative solving of Newton's equations of motion). It focuses on two main objectives. One is to calculate the relative stability of different states of the system. A typical problem that has' such an objective is computer-aided drug design. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), "native" state. Perhaps the best example of such a problem is protein folding. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to "quasi non-ergodicity", whereby a part of phase space is inaccessible on time scales of the simulation. To overcome this difficulty and to extend molecular dynamics to "biological" time scales (millisecond or longer) new physical formulations and new algorithmic developments are required. To be efficient they should account for natural limitations of multi-processor computer architecture. I will present work along these lines done in my group. In particular, I will focus on a new approach to calculating the free energies (stability) of different states and to overcoming "the curse of rare events". I will also discuss algorithmic improvements to multiple time step methods and to the treatment of slowly decaying, log-ranged, electrostatic effects.

Pohorille, Andrew↗

Self-organized criticality in sandpiles - Nature of the critical phenomenon

The scaling behavior of sandpile models is investigated analytically. First, it is shown that sandpile models contain a set of domain walls, referred to as troughs, which bound regions that can experience avalanches. It is further shown that the dynamics of the troughs is governed by a simple set of rules involving birth, death, and coalescence events. A simple trough model is then introduced, and it is proved that the model has a phase transition with the density of the troughs as an order parameter and that, in the thermodynamic limit, the trough density goes to zero at the transition point. Finally, it is shown that the observed scaling behavior is a consequence of finite-size effects.

Carlson, J. M.↗

Experimental Investigation of Flow Orientation Effects on Cryogenic Flow Boiling

Cryogenic flow boiling physics is of paramount importance for future cryogenic space applications employing Cryogenic Fluid Management (CFM) technologies. This experimental study investigates the two-phase flow and heat transfer performance of LN 2 flow boiling across five distinct flow orientations: vertical upflow, vertical downflow, horizontal flow, 45° inclined upflow, and 45° inclined downflow. The study employed the steady- state heating method within a circular heated tube featuring an 8.5-mm inner diameter and a 680-mm heated length. High-speed video recordings were utilized to capture two-phase flow patterns and interfacial behaviors across various flow orientations. The experiments covered a wide range of operating conditions, including mass velocities ranging from 351.80 to 1572.77 kg/m 2 s and inlet pressures from 297.14 to 1032.97 kPa, primarily with near-saturated inlet subcooling. Distinct two-phase flow patterns and regime transitions were identified for each flow orientation. Symmetrical flow patterns were evident in vertical orientations, whereas non-vertical orientations exhibited asymmetric flow stratifications, primarily influenced by the buoyancy force in a terrestrial gravity environment. Bubble dynamic parameters were quantified, and bubble collision and dispersion phenomena were visualized. Analyzing heat transfer performance based on local flow boiling curves and variations in heat transfer coefficients (HTCs), it was observed that vertical upflow demonstrated the most enhanced heat transfer performance, while vertical downflow exhibited the lowest. As mass velocity exceeded 830 kg/ m 2 s, the differences in heat transfer among orientations became less distinct, emphasizing the role of flow inertia in mitigating the influence of flow orientation. A direct comparison of microgravity HTC data from a recent study by the authors against the 1-𝘨 e HTC data indicated enhanced heat transfer performance in microgravity for flow orientation configurations in terrestrial gravity. However, this enhancement progressively diminished as the heat flux increased. Distinct HTC trends were observed as mass velocity increased for different flow orientations and gravity levels.

Cryogens↗

Physics of Hard Spheres Experiment: Significant and Quantitative Findings Made

Direct examination of atomic interactions is difficult. One powerful approach to visualizing atomic interactions is to study near-index-matched colloidal dispersions of microscopic plastic spheres, which can be probed by visible light. Such spheres interact through hydrodynamic and Brownian forces, but they feel no direct force before an infinite repulsion at contact. Through the microgravity flight of the Physics of Hard Spheres Experiment (PHaSE), researchers have sought a more complete understanding of the entropically driven disorder-order transition in hard-sphere colloidal dispersions. The experiment was conceived by Professors Paul M. Chaikin and William B. Russel of Princeton University. Microgravity was required because, on Earth, index-matched colloidal dispersions often cannot be density matched, resulting in significant settling over the crystallization period. This settling makes them a poor model of the equilibrium atomic system, where the effect of gravity is truly negligible. For this purpose, a customized light-scattering instrument was designed, built, and flown by the NASA Glenn Research Center at Lewis Field on the space shuttle (shuttle missions STS 83 and STS 94). This instrument performed both static and dynamic light scattering, with sample oscillation for determining rheological properties. Scattered light from a 532- nm laser was recorded either by a 10-bit charge-coupled discharge (CCD) camera from a concentric screen covering angles of 0 to 60 or by sensitive avalanche photodiode detectors, which convert the photons into binary data from which two correlators compute autocorrelation functions. The sample cell was driven by a direct-current servomotor to allow sinusoidal oscillation for the measurement of rheological properties. Significant microgravity research findings include the observation of beautiful dendritic crystals, the crystallization of a "glassy phase" sample in microgravity that did not crystallize for over 1 year in 1g (Earth's gravity), and the emergence of face-centered-cubic (FCC) crystals late in the coarsening process (as small crystallites lost particles to the slow ripening of large crystallites). Significant quantitative findings from the microgravity experiments have been developed describing complex interactions among crystallites during the growth process, as concentration fields overlap in the surrounding disordered phase. Time-resolved Bragg scattering under microgravity captures one effect of these interactions quite conclusively for the sample at a volume fraction of 0.528. From the earliest time until the sample is almost fully crystalline, the size and overall crystallinity grow monotonically, but the number of crystallites per unit volume (number density) falls. Apparently nucleation is slower than the loss of crystallites because of the transfer of particles from small to large crystals. Thus, coarsening occurs simultaneously with growth, rather than following the completion of nucleation and growth as is generally assumed. In the same sample, an interesting signature appears in the apparent number density of crystallites and the volume fraction within the crystallites shortly before full crystallinity is reached. A brief upturn in both indicates the creation of more domains of the size of the average crystallite simultaneous with the compression of the crystallites. Only the emergence of dendritic arms offers a reasonable explanation. The arms would be "seen" by the light scattering as separate domains whose smaller radii of curvature would compress the interior phase. In fiscal year 1999, numerous papers, a doctoral dissertation, and the PHaSE final report were produced. Although this flight project has been completed, plans are in place for a follow-on colloid experiment by Chaikin and Russel that employs a light microscope within Glenn's Fluids and Combustion Facility on the International Space Station. PHaSE is providing us with a deeper understanding of the nure of phase transitions. The knowledge derived has added to the understanding of condensed matter. In addition, the burgeoning study of the dynamics of colloidal self-assembly may lead to the development of a range of photonic materials that control the desirable properties of light. Thus, applications of ordered colloidal structures include not only ultrastructure ceramics, but also photonic crystals and photothermal nanosecond light-switching devices. Industries dealing with semiconductors, electro-optics, ceramics, and composites stand to benefit from such advancements.

Doherty, Michael P.↗

Photochemical modeling of the Antarctic stratosphere: Observational constraints from the airborne Antarctic ozone experiment and implications for ozone behavior

The rapid decrease in O3 column densities observed during Antarctic spring has been attributed to several chemical mechanisms involving nitrogen, bromine, or chlorine species, to dynamical mechanisms, or to a combination of the above. Chlorine-related theories, in particular, predict greatly elevated concentrations of ClO and OClO and suppressed abundances of NO2 below 22 km. The heterogeneous reactions and phase transitions proposed by these theories could also impact the concentrations of HCl, ClNO3 and HNO3 in this region. Observations of the above species have been carried out from the ground by the National Ozone Expedition (NOZE-I, 1986, and NOZE-II, 1987), and from aircrafts by the Airborne Antarctic Ozone Experiment (AAOE) during the austral spring of 1987. Observations of aerosol concentrations, size distribution and backscattering ratio from AAOE, and of aerosol extinction coefficients from the SAM-II satellite can also be used to deduce the altitude and temporal behavior of surfaces which catalyze heterogeneous mechanisms. All these observations provide important constraints on the photochemical processes suggested for the spring Antarctic stratosphere. Results are presented for the concentrations and time development of key trace gases in the Antarctic stratosphere, utilizing the AER photochemical model. This model includes complete gas-phase photochemistry, as well as heterogeneous reactions. Heterogeneous chemistry is parameterized in terms of surface concentrations of aerosols, collision frequencies between gas molecules and aerosol surfaces, concentrations of HCl/H2O in the frozen particles, and probability of reaction per collision (gamma). Values of gamma are taken from the latest laboratory measurements. The heterogeneous chemistry and phase transitions are assumed to occur between 12 and 22 km. The behavior of trace species at higher altitudes is calculated by the AER 2-D model without heterogeneous chemistry. Calculations are performed for solar illumination conditions typical of 60, 70, and 80 S, from July 15 to October 31.

Rodriguez, Jose M.↗

Surface Instability of Liquid Propellant under Vertical Oscillatory Forcing

Fluid motion in a fuel tank produced during thrust oscillations can circulate sub-cooled hydrogen near the liquid-vapor interface resulting in increased condensation and ullage pressure collapse. The first objective of this study is to validate the capabilities of a Computational Fluid Dynamics (CFD) tool, CFD-ACE+, in modeling the fundamental interface transition physics occurring at the propellant surface. The second objective is to use the tool to assess the effects of thrust oscillations on surface dynamics. Our technical approach is to first verify the CFD code against known theoretical solutions, and then validate against existing experiments for small scale tanks and a range of transition regimes. A 2D axisymmetric, multi-phase model of gases, liquids, and solids is used to verify that CFD-ACE+ is capable of modeling fluid-structure interaction and system resonance in a typical thrust oscillation environment. Then, the 3D mode is studied with an assumed oscillatory body force to simulate the thrust oscillating effect. The study showed that CFD modeling can capture all of the transition physics from solid body motion to standing surface wave and to droplet ejection from liquid-gas interface. Unlike the analytical solutions established during the 1960 s, CFD modeling is not limited to the small amplitude regime. It can extend solutions to the nonlinear regime to determine the amplitude of surface waves after the onset of instability. The present simulation also demonstrated consistent trends from numerical experiments through variation of physical properties from low viscous fluid to high viscous fluids, and through variation of geometry and input forcing functions. A comparison of surface wave patterns under various forcing frequencies and amplitudes showed good agreement with experimental observations. It is concluded that thrust oscillations can cause droplet formation at the interface, which results in increased surface area and enhanced heat transfer between the liquid and gas phases as the ejected droplets travel well into the warmer gas region.

Yang, H. Q.↗

Comparison of dynamic fatigue behavior between SiC whisker-reinforced composite and monolithic silicon nitrides

The dynamic fatigue behavior of 30 vol percent silicon nitride whisker-reinforced composite and monolithic silicon nitrides were determined as a function of temperature from 1100 to 1300 C in ambient air. The fatigue susceptibility parameter, n, decreased from 88.1 to 20.1 for the composite material, and from 50.8 to 40.4 for the monolithic, with increasing temperature from 1100 to 1300 C. A transition in the dynamic fatigue curve occurred for the composite material at a low stressing rate of 2 MPa/min at 1300 C, resulting in a very low value of n equals 5.8. Fractographic analysis showed that glassy phases in the slow crack growth region were more pronounced in the composite compared to the monolithic material, implying that SiC whisker addition promotes the formation of glass rich phases at the grain boundaries, thereby enhancing fatigue. These results indicate that SiC whisker addition to Si3 N4 matrix substantially deteriorates fatigue resistance inherent to the matrix base material for this selected material system.

Choi, Sung R.↗

Multi-Arc Filtering During the Navigation Campaign of the OSIRIS-REx Mission

The Navigation Campaign of the OSIRIS-REx mission consisted of three phases: Approach, Preliminary Survey and Orbital-A. These phases were designed to optimize the initial characterization of Bennu's mass, shape and spin-state to support a safe orbit insertion and a quick transition to landmark-based optical navigation tracking. The standard orbit determination filtering techniques used to navigate the spacecraft were unable to fit data from these three phases simultaneously due to numerical issues associated with the nonlinear dynamics and the long arc length. Consequently, a multi-arc filtering algorithm was implemented in order to combine the information from each of these arcs. Multi-arc solutions for Bennu's spin state and gravity field are presented here.

French, Andrew S.↗

Cluster analysis of multiple planetary flow regimes

A modified cluster analysis method was developed to identify spatial patterns of planetary flow regimes, and to study transitions between them. This method was applied first to a simple deterministic model and second to Northern Hemisphere (NH) 500 mb data. The dynamical model is governed by the fully-nonlinear, equivalent-barotropic vorticity equation on the sphere. Clusters of point in the model's phase space are associated with either a few persistent or with many transient events. Two stationary clusters have patterns similar to unstable stationary model solutions, zonal, or blocked. Transient clusters of wave trains serve as way stations between the stationary ones. For the NH data, cluster analysis was performed in the subspace of the first seven empirical orthogonal functions (EOFs). Stationary clusters are found in the low-frequency band of more than 10 days, and transient clusters in the bandpass frequency window between 2.5 and 6 days. In the low-frequency band three pairs of clusters determine, respectively, EOFs 1, 2, and 3. They exhibit well-known regional features, such as blocking, the Pacific/North American (PNA) pattern and wave trains. Both model and low-pass data show strong bimodality. Clusters in the bandpass window show wave-train patterns in the two jet exit regions. They are related, as in the model, to transitions between stationary clusters.

Mo, Kingtse↗