Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stability regime”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

New Perspective On Slosh Dynamics in High-Gravity Regimes for Lunar Missions

Slosh dynamics particularly in high gravity regime is of high interest to the NASA Human Landing Systems program, which aims to return humans to the Moon. This paper reexamines the commonly accepted equivalence between pendulum and mass-spring-damper mechanical analogies, demonstrating the differences through simulations and their impact on design and analysis. This study re-derives equations of motion for these mechanical models including mass-spring-damper (MSD) and pendulum by relaxing simplifying assumptions that had been typically applied to large launch vehicles and may be violated for lunar and planetary landing applications. The newly derived equations uncovers potential modeling gaps when applying the pendulum and MSD models, and the need for additional considerations in their application. Finally, the newly derived equations are validated using MathWorks® Simscape™ Multibody™ dynamics toolbox simulation.

Slosh↗
DOE OSTI2024

Reexpansion of charged nanoparticle assemblies in concentrated electrolytes

Electrostatic forces in solutions are highly relevant to a variety of fields, ranging from electrochemical energy storage to biology. However, their manifestation in concentrated electrolytes is not fully understood, as exemplified by counterintuitive observations of colloidal stability and long-ranged repulsions in molten salts. Highly charged biomolecules, such as DNA, respond sensitively to ions in dilute solutions. Here, we use non-base-pairing DNA-coated nanoparticles (DNA-NP) to analyze electrostatic interactions in concentrated salt solutions. Despite their negative charge, these conjugates form colloidal crystals in solutions of sufficient divalent cation concentration. We utilize small-angle X-ray scattering (SAXS) to study such DNA-NP assemblies across the full accessible concentration ranges of aqueous CaCl 2 , MgCl 2 , and SrCl 2 solutions. SAXS shows that the crystallinity and phases of the assembled structures vary with cation type. For all tested salts, the aggregates contract with added ions at low salinities and then begin expanding above a cation-dependent threshold salt concentration. Wide-angle X-ray scattering (WAXS) reveals enhanced positional correlations between ions in the solution at high salt concentrations. Complementary molecular dynamics simulations show that these ion–ion interactions reduce the favorability of dense ion configurations within the DNA brushes below that of the bulk solution. Measurements in solutions with lowered permittivity demonstrate a simultaneous increase in ion coupling and decrease in the concentration at which aggregate expansion begins, thus confirming the connection between these phenomena. Our work demonstrates that interactions between charged objects continue to evolve considerably into the high-concentration regime, where classical theories project electrostatics to be of negligible consequence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRSDate not supplied

New Perspective on Slosh Dynamics in High-Gravity Regimes for Lunar Missions

Slosh dynamics particularly in high gravity regime is of high interest to the NASA Human Landing Systems program, which aims to return humans to the Moon. This paper reexamines the commonly accepted equivalence between pendulum and massspring-damper mechanical analogies, demonstrating the differences through simulations and their impact on design and analysis. This study re-derives equations of motion for these mechanical models including mass-spring-damper (MSD) and pendulum by relaxing simplifying assumptions that had been typically applied to large launch vehicles and may be violated for lunar and planetary landing applications. The newly-derived equations uncover potential modeling gaps when applying the pendulum and MSD models, and the need for additional considerations in their application. Finally, the newly-derived equations are validated using MathWorks® Simscape™ Multibody™ dynamics toolbox simulation.

Pendulum↗
DOE OSTI2024

A large deformation multiphase continuum mechanics model for shock loading of soft porous materials

A large deformation, coupled finite-element (FE) model is developed to simulate the multiphase response of soft porous materials subjected to high strain-rate loading. The approach is based on the theory of porous media (TPM) at large deformations. Simplifications to the one-dimensional regime studied in the numerical simulations follow. An overview of several different time integration schemes is presented for the purpose of solving the nonlinear dynamic coupled balance of momenta (mixture and fluid) and balance of mass of the mixture equations. Numerical examples are presented for (i) verification against closed-form analytical solutions assuming small loads, (ii) demonstrating large deformation effects at high strain-rate, and (iii) showing differences in deformations between a single-phase elastodynamics model with occluded compressible pore fluid and a multiphase poroelastodynamics model at high strain-rate. The multiphase model shows that the relative motion of the pore fluid significantly dampens the deformation response of the solid skeleton as compared to the single-phase model, and makes it possible to extract quantitative values for the stresses of the different constituents, thereby allowing one to form preliminary conclusions about the onset of damage in the solid skeleton. The novelty of the current work is developing a multiphase, large deformation, mixture theory numerical model for high strain-rate loading of soft porous materials. It was discovered that explicit, adaptive time-stepping Runge–Kutta schemes offer high accuracy at relatively low cost when compared to traditional implicit or explicit central difference time-stepping schemes for shock-like loadings. Here, shock viscosity is added to the mixture momentum balance equation to regularize the shock front, and a stabilization term is added to the mixture mass balance equation to stabilize equal order interpolation finite elements for the coupled finite element solution of multiphase materials.

Engineering↗
DOE OSTI2024

Invasive Plant Species Management Plan for Los Alamos National Laboratory (Rev.3)

Native species are plants and animals that continually occupy a natural range without direct or indirect introduction and/or care by humans. They are adapted to the environmental conditions and processes of the ecosystem in which they reside. A species introduced into a novel ecosystem can either exploit that ecosystem and thrive or be unable to survive in that ecosystem (Hobbs et al. 2006). Alien or non-native species are species that are intentionally or accidentally introduced into a novel ecosystem and are capable of living and propagating within the physical parameters of that ecosystem. Invasive species is a species that is non-native (or alien) to the ecosystem under consideration and a noxious species are those whose introduction causes or is likely to cause economic or environmental harm. The term invasive species is applicable to plants and animals alike; however, this invasive species management plan currently focuses on invasive vegetation. Invasive plant species are usually capable of rapid colonization of disturbed ground, such as after a change in wildfire regime intensity and frequency (Reilly et al. 2020) or anthropogenic ground disturbance (Burke and Grime 1996; Hobbs and Huenneke 1992). Climate change facilitates the spread and establishment of many alien species and creates new opportunities for them to become invasive (Turbelin and Catford 2021). As climate change impacts increase in the coming decades, we may see in an increase in invasive species establishment. See a list of definitions of terms pertinent to this document in Appendix A: Definitions. Los Alamos National Laboratory (LANL) hosts populations of non-native and invasive species all typical of the northern New Mexico region (Martin 2004; NMDA 2020). By implementing an invasive species management plan, LANL will have readily accessible management strategies for invasive species that are found on-site. The benefits of managing invasive species include a decrease in wildland fire risks, an increase in soil productivity, an increase in (productive or beneficial) wildlife habitat, an increase in water quantity and quality, and the restoration of impacted areas (Burke and Grime 1996; Hobbs and Huenneke 1992; D’Antonio and Hobbie 2005; MacDougall et al. 2013; Reilly et al. 2020). The aforementioned benefits from invasive species management directly enable the LANL mission by ensuring compliance requirements are met and site-wide programs, such as the Vulnerability Assessment and Resilience Plan, are supported for a mutually beneficial outcome. An example for the LANL site specifically, controlling non-native annual plant species, for example, could reduce the costs associated with stabilizing soils during the stormwater pollution prevention compliance process. Roadway and utility right-of-way areas are another place where an integrated vegetation management strategy would promote low growing perennial plants in a way that is mutually beneficial to habitat and the institution through lowered maintenance costs. There is also an economic benefit to managing invasive species. In one nationwide study, invasive plants had an estimated impact cost of $190.45 billion (Fantel-Lepczyk et al. 2022). Investing in preventative measures and surveillance could help to offset future control and management costs of invasive species that have the potential to become established. The State of New Mexico has developed plant species lists and recommendations through the New Mexico Department of Agriculture’s Noxious Weeds Management Act, Article 7D (NM Statute § 76-7D-4 2021); however, effective invasive species management must rely on local knowledge of the site and region. Los Alamos County (LAC) has already compiled a target invasive plant species list and species-specific management objectives (Martin 2004). Therefore, the New Mexico Noxious Weed List and the LAC invasive plant species list, as well as management objectives from those documents, are integrated into LANL’s invasive plant species management plan.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Applicability of semiclassical theories in the strong-field plasma regime

For many purposes, classical plasma dynamics models can work surprisingly well, even for strong electromagnetic fields, approaching the Schwinger critical fields, and high frequencies, approaching the Compton frequency. However, the applicability of classical models tends to depend rather sensitively on the details of the problem. In the present paper, we study the specific case of plasma oscillations to draw a line between the classical and quantum relativistic regimes. Here, due to the field geometry of study, mechanisms like radiation reaction and Breit-Wheeler pair production, which tend to be important for electromagnetic fields, are rather effectively suppressed. Moreover, we find that the polarization current due to the electron spin is generally negligible for frequencies below the Compton frequency, compared with the free current, whose magnitude is well-approximated by the classical Vlasov theory. However, we show that pair creation due to the Schwinger mechanism can sometimes be important for surprisingly modest field strengths, of the order of 10% of the critical field or even smaller. A rough guideline for when the classical Vlasov theory can be applied is given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Convergence acceleration of implicit schemes in the presence of high aspect ratio grid cells

The performance of Navier-Stokes codes are influenced by several phenomena. For example, the robustness of the code may be compromised by the lack of grid resolution, by a need for more precise initial conditions or because all or part of the flowfield lies outside the flow regime in which the algorithm converges efficiently. A primary example of the latter effect is the presence of extended low Mach number and/or low Reynolds number regions which cause convergence deterioration of time marching algorithms. Recent research into this problem by several workers including the present authors has largely negated this difficulty through the introduction of time-derivative preconditioning. In the present paper, we employ the preconditioned algorithm to address convergence difficulties arising from sensitivity to grid stretching and high aspect ratio grid cells. Strong grid stretching is particularly characteristic of turbulent flow calculations where the grid must be refined very tightly in the dimension normal to the wall, without a similar refinement in the tangential direction. High aspect ratio grid cells also arise in problems that involve high aspect ratio domains such as combustor coolant channels. In both situations, the high aspect ratio cells can lead to extreme deterioration in convergence. It is the purpose of the present paper to address the reasons for this adverse response to grid stretching and to suggest methods for enhancing convergence under such circumstances. Numerical algorithms typically possess a maximum allowable or optimum value for the time step size, expressed in non-dimensional terms as a CFL number or vonNeumann number (VNN). In the presence of high aspect ratio cells, the smallest dimension of the grid cell controls the time step size causing it to be extremely small, which in turn results in the deterioration of convergence behavior. For explicit schemes, this time step limitation cannot be exceeded without violating stability restrictions of the scheme. On the other hand, for implicit schemes, which are typically unconditionally stable, there appears to be room for improvement through careful tailoring of the time step definition based on results of linear stability analyses. In the present paper, we focus on the central-differenced alternating direction implicit (ADI) scheme. The understanding garnered from this analyses can then be applied to other implicit schemes. In order to systematically study the effects of aspect ratio and the methods of mitigating the associated problems, we use a two pronged approach. We use stability analyses as a tool for predicting numerical convergence behavior and numerical experiments on simple model problems to verify predicted trends. Based on these analyses, we determine that efficient convergence may be obtained at all aspect ratios by getting a combination of things right. Primary among these are the proper definition of the time step size, proper selection of viscous preconditioner and the precise treatment of boundary conditions. These algorithmic improvements are then applied to a variety of test cases to demonstrate uniform convergence at all aspect ratios.

Buelow, B. E. O.↗
DOE OSTI2023

Asymptotic vacuum solution at tokamak X-point tip

In the H-mode regime of diverted tokamaks, the presence of strong pressure gradients in the pedestal gives rise to a sizable bootstrap current, together with the Ohmic and Pfirsch–Schlueter currents, close to the separatrix. For such equilibria, the presence of finite current density close to the separatrix requires the reexamination of equilibrium properties. It is almost universally assumed that the two branches of the separatrix (the stable and unstable manifolds) are straight as they cross at the X-point. However, the opposite angles of the plasma-filled segment and vacuum one cannot be equal if the current density does not vanish at the separatrix on the plasma side. We solve this difficulty by chipping off a thin layer of plasma edge so that the sharp corner of the plasma-filled segment becomes a hyperbola. Using the conformal transformation, we found that in the assumption of a hyperbolic boundary, the X point moves beyond the plasma boundary to fall in the vacuum region. An acute angle of the plasma-filled segment leads to an obtuse opposite angle of vacuum segment and vice versa. In the case of an acute angle of the plasma-filled segment, the new X point shifts inside the X point formed by the asymptotes of a hyperbolic boundary; in the case of an obtuse angle of the plasma-filled segment, the new X point shifts outside the X point formed by the asymptotes of a hyperbolic plasma boundary. Furthermore, the results are important for understanding the X point features, which affect the tokamak edge stability and transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2022

Linear Disturbance Amplification Over Blunted Flat Plates in High-Speed Flows

Modal and nonmodal instability characteristics of cylindrically blunted flat plates with varying leading edge radii are described for Mach 4 and Mach 6 freestream conditions. The selection of leading edge radii and freestream parameters is informed by experimental conditions. The investigation of this 2D problem provides a slow entropy layer swallowing which allows for an isolated development of perturbations seeded upstream within different wall-normal regions of the flow. At both Mach numbers, a decrease in modal instability amplification was seen as the leading edge radius was increased. Nonmodal analysis reveals amplifying perturbations in the boundary layer as well as in the entropy layer. The medium bluntness regime exhibits the strongest amplification of nonmodal disturbances that is nonmonotonic in character. Small amplitude boundary forcing at the plate surface or volumetric forcing at various wall-normal heights was used to account for receptivity effects. While wall forcing effectively induced modal instabilities, only an actuation above the boundary layer captured disturbances that amplify within the entropy layer. The optimal nonmodal theory’s entropy-layer disturbance evolution exhibited outstanding agreement with controlled forcing, including receptivity effects. The evolution of entropy-layer disturbances from the optimal nonmodal theory showed excellent agreement with the results of receptivity to controlled forcing. Therefore, the nonmodal optimal growth analysis may provide a useful as well as efficient technique to identify the complete disturbance spectrum in blunt hypersonic configurations, where both modal and nonmodal disturbances can amplify in the boundary-layer and entropy-layer regions.

Boundary layer transition↗
NASA NTRS2016

Cislunar Near Rectilinear Halo Orbit for Human Space Exploration

In order to conduct sustained human exploration beyond Low Earth Orbit (LEO), spacecraft systems are designed to operate in a series of missions of increasing complexity. Regardless of the destination, Moon, Mars, asteroids or beyond, there is a substantial set of common objectives that must be met. Many orbit characterization studies have endeavored to evaluate the potential locations in cislunar space that are favorable for meeting common human exploration objectives in a stepwise approach. Multiple studies, by both NASA and other international space agencies, have indicated that Earth-­‐moon libration point orbits are attractive candidates for staging operations in the proving ground and beyond. In particular, the Near Rectilinear Orbit (NRO) has been demonstrated to meet multi-­‐mission and multi-­‐destination architectural constraints. However, a human mission to a selected NRO presents a variety of new challenges for mission planning. While a growing number of robotic missions have completed successful operations to various specific libration point orbits, human missions have never been conducted to orbits of this class. Human missions have unique challenges that differ significantly from robotic missions, including a lower tolerance for mission risk and additional operational constraints that are associated only with human spacecraft. In addition, neither robotic nor human missions have been operated in the NRO regime specifically, and NROs exhibit dynamical characteristics that can differ significantly as compared to other halo orbits. Finally, multi-­‐body orbits, such as libration point orbits, are identified to exist in a simplified orbit model known as the Circular Restricted Three Body Problem (CRTBP) and must then be re-­‐solved in the full ephemeris model. As a result, the behavior of multi-­‐body orbits cannot be effectively characterized within the classical two-­‐body orbit dynamics framework more familiar to the human spaceflight community. In fact, a given NRO is not identified by a set of Keplerian orbit parameters, and a valid epoch specific state vector must be first obtained from a multi-body dynamical model. In this paper, the significant performance and operational challenges of conducting human missions to the NRO are evaluated. First, a systematic process for generating full ephemeris based ballistic NROs of various families is outlined to demonstrate the relative ease in which a multi-­‐revolution orbit can be found for any epoch and for various orbit geometries. In the Earth-­‐Moon system, NROs, which are halo orbits with close passage over a lunar pole, can exist with respect to libration point 1 (L1) or libration point 2 (L2) and are either from a North or South family orbit class with respect to the ecliptic. Second, the ability to maintain the orbit over the lifetime of a habitat mission by applying a reliable station-keeping strategy is investigated. The NRO, while similar to the quasi-­‐halo orbits that the Artemis mission flew, requires an updated station keeping strategy. This is due to several dynamical differences such as the increased relative stability of the NRO compared to other halo orbits and the close passage over the lunar surface as shown in Figure 1. Multiple station-keeping strategies are being investigated to ensure a human spacecraft remains on a predictable path. As the NRO is not described in simple two-­‐body parameters, analysis must determine the best strategy for targeting a reference NRO as well as how closely a future state should be constrained. In addition, costs will be minimized by determining maneuver directionality based on an identified pattern in the optimal station-keeping solutions or an analytically derived relationship. The candidate station-keeping algorithm must be stable and robust to environmental and vehicle uncertainties as well to navigation estimation and flight control execution errors. To that end, navigation accuracies, the impact on the station-keeping execution errors as well as other vehicle uncertainties need to be assessed. Starting with Orion, current navigation accuracies are evaluated and then navigation requirements are derived assuming a desired station-keeping propellant budget. Third, the performance requirements to and from the NRO are evaluated. Important parameters for developing expected propellant costs include epoch of operation, size and type of NRO, Earth departure and return constraints, as well as abort or early-­‐return capability. Finally, rendezvous and proximity operations are vital aspects of multi-­‐mission human exploration endeavors. The ability to conduct rendezvous and the associated propellant costs are assessed as well as the impacts of various profile assumptions including the location within the NRO the rendezvous is performed. The results of these studies will influence plans for international cooperation on both nearer term proving ground missions and beyond.

Whitley, Ryan↗
DOE OSTI2022

Parametric model-order-reduction development for unsteady convection

A time-averaged error indicator with POD- h Greedy is developed to drive parametric model order reduction (pMOR) for 2D unsteady natural convection in a high-aspect ratio slot parameterized with the Prandtl number, Rayleigh number, and slot angle with respect to the gravity. The error indicator is extended to accommodate the energy equation and Leray regularization. Despite being two-dimensional and laminar, the target flow regime presents several challenges: 1) there is a bifurcation in the angle parameter space; 2) the solution can be multivalued, even at steady state; and 3) the solution exhibits spatio-temporal chaos at several points in the parameter space. The authors explore several reduced-order models (ROMs) and demonstrate that Leray-regularized Galerkin ROMs provide a robust solution approach for this class of flows. They further demonstrate that error-indicated pMOR can efficiently predict several QOIs, such as mean flow, mean Nusselt number and mean turbulent kinetic energy, even in the presence of a bifurcation. Finally, they show that spatio-temporal chaos can lead to lack of reproducibility in both the full-order model and the reduced-order model and that the variance in the full-order model provides a lower bound on the pMOR error in these cases.

leray regularization↗
NASA NTRS1993

Long-term microparticle flux variability indicated by comparison of Interplanetary Dust Experiment (IDE) timed impacts for LDEF's first year in orbit with impact data for the entire 5.77-year orbital lifetime

The electronic sensors of the Interplanetary Dust Experiment (IDE) recorded precise impact times and approximate directions for submicron to approximately 100 micron size particles on all six primary sides of the spacecraft for the first 346 days of the LDEF orbital mission. Previously-reported analyses of the timed impact data have established their spatio-temporal features, including the demonstration that a preponderance of the particles in this regime are orbital debris and that a large fraction of the debris particles are encountered in megameter-size clouds. Short-term fluxes within such clouds can rise several orders of magnitude above the long-term average. These unexpectedly large short-term variations in debris flux raise the question of how representative an indication of the multi-year average flux is given by the nearly one year of timed data. One of the goals of the IDE was to conduct an optical survey of impact sites on detectors that remained active during the entire LDEF mission, to obtain full-mission fluxes. We present here the comparisons and contrasts among the new IDE optical survey impact data, the IDE first-year timed impact data, and impact data from other LDEF micrometeoroid and debris experiments. The following observations are reported: (1) the 5.77 year long-term integrated microparticle impact fluxes recorded by IDE detectors matched the integrated impact fluxes measured by other LDEF investigators for the same period; (2) IDE integrated microparticle impact fluxes varied by factors from 0.5 to 8.3 for LDEF days 1-346, 347-2106 and 1-2106 (5.77 years) on rows 3 (trailing edge, or West), 6 (South side), 12 (North side), and the Earth and Space ends; and (3) IDE integrated microparticle impact fluxes varied less than 3 percent for LDEF days 1-346, 347-2106 and 1-2106 (5.77 years) on row 9 (leading edge, or East). These results give further evidence of the accuracy and internal consistency of the recorded IDE impact data. This leads to the further conclusion that the utility of long-term ratios for impacts on various sides of a stabilized satellite in low Earth orbit (LEO) is extremely limited. These observations and their consequences highlight the need for continuous, real time monitoring of the dynamic microparticle environment in LEO.

Simon, Charles G.↗
DOE OSTI2024

Investigation of the Impact of Flow on MHD Perturbations with the NIMROD Code (Final Report)

This work will enable improved understanding of tokamak edge flows and impurity-species dynamics and their impact on MHD stability as it relates to ELMs, RMPs and QH-mode by leveraging computational developments of the NIMROD code. Experimentally, it is established that the tokamak flows and impurities species can have a substantial effect on MHD stability. This is particularly true for the edge plasma where, for example, large flow shear is correlated with the occurrence of the Quiescent H-mode (QH- mode) state as opposed to operation with edge localized modes (ELMs). Understanding plasma flows is challenging because of the many physical effects that come into play: particle orbits loss, neutrals dynamics, and interaction between multiple ion species with decoupled motion. We will develop a model that incorporates multiple ion species into the MHD framework that includes self-consistent magnetic-field evolution. With this model we will understand how the presence of impurities impacts the tokamak edge-pedestal flows. Finally, we will quantify the impact of the new multiple species modeling on 3D QH-mode simulations. This work will have broader impacts on all MHD studies when multiple species are present (as is the case in modern tokamaks), as well as contributing to the edge and transport communities who also are interested in the studies of multiple-ion-species flows in the edge pedestal region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2010

Possible Influences of Air Pollution, Dust and Sandstorms on the Indian Monsoon

In Asian monsoon countries, such as China and India, human health and safety problems caused by air pollution are becoming increasingly serious, due to the increased loading of atmospheric pollutants from waste gas emissions and from rising energy demand associated with the rapid pace of industrialization and modernization. Meanwhile, uneven distribution of monsoon rain associated with flash floods or prolonged drought, has caused major loss of human life and damage to crops and.property with devastating societal impacts. Historically, air-pollution and monsoons research are treated as separate problems. However recent studies have suggested that the two problems may be intrinsically linked and need to be studied jointly. Fundamentally, aerosols can affect precipitation through radiative effects cif suspended particles in the atmosphere (direct effect) and/or by interfering and changing: the cloud and precipitation formation processes (indirect effect). Based on their optical properties, aerosols can be classified into two types.: those that absorb solar radiation, and those that do not. Both types of aerosols scatter sunlight and reduce the amount of solar radiation from reaching the Earth's surface, causing it to cool. The surface cooling increases atmospheric stability and reduces convection potential, Absorbing aerosols, however, in addition to cooling the surface, can heat the atmosphere. The heating of the atmosphere may reduce the amount of low clouds by increased evaporation in cloud drops. The heating, however, may induce rising motion, enhance low-level moisture, convergence and, hence, increases rainfall, The latent heating from enhanced rainfall may excite feedback processes in the large-scale circulation, further amplify.the initial response to aerosol heating and producing more rain. Additionally, aerosols can increase the concentration of cloud condensation nuclei (CCN), increase cloud amount and decrease coalescence and collision rates, leading to reduced precipitation. However, in the presence of increasing moist and warm air, the reduced coalescence/collision may lead to supercooled drops at higher altitudes where ice precipitation falls and melts. The latent heat release from freezing aloft and melting below implies greater upward heat transport in polluted clouds and invigorate deep convection. In this way, aerosols may lead to increased local convection. Hence, depending on the ambient large-scale conditions and dynamical feedback processes, aerosols' effect on precipitation can be positive, negative or mixed. In the Asian monsoon and adjacent regions, the aerosol forcing and responses of the water cycle are even more complex, Both direct and indirect effects may take place locally and simultaneously, interacting with each other. in addition to local effects, monsoon rainfall may be affected by aerosols transported from other regions and intensified through large-scale circulation and moisture feedback. Thus, dust transported by the large-scale circulation from the adjacent deserts to northern India may affect rainfall over the Bay of Bengal; sulphate and black carbon front industrial pollution in central, southern China and northern India may affect the rainfall regime over the Korean peninsula and Japan; organic and black carbon front biomass burning from Indo-China may modulate the pre-monsoon rainfall regime over southern China and coastal regions, contributing to variability in differential heating and cooling of the atmosphere and to the land-sea thermal contrast. During the pre-monsoon season and monsoon breaks, it has been suggested that radiative forcing by absorbing aerosols have nearly the same order of magnitude as the forcing due to latent heating from convection and surface fluxes. The magnitude of the total aerosol radiative cooling due to sulphates and soot is of the order of 20-40 W/m2 over the Asian monsoon land region in the pre-monsoon season, compared to about 1-2 W/m2 for global warng. However, the combined forcing at the surface and in the atmosphere, including all species. if aerosols, and details of aerosol mixing, and impacts on the energy and water cycles in the monsoon land regions, are not well known.

Lau, William K. M.↗
NASA NTRS1999

Evolution of Local Microstructures: Spatial Instabilities of Coarsening Clusters

This work examines the diffusional growth of discrete phase particles dispersed within a matrix. Engineering materials are microstructurally heterogeneous, and the details of the microstructure determine how well that material performs in a given application. Critical to the development of designing multiphase microstructures with long-term stability is the process of Ostwald ripening. Ripening, or phase coarsening, is a diffusion-limited process which arises in polydisperse multiphase materials. Growth and dissolution occur because fluxes of solute, driven by chemical potential gradients at the interfaces of the dispersed phase material, depend on particle size. The kinetics of these processes are "competitive," dictating that larger particles grow at the expense of smaller ones, overall leading to an increase of the average particle size. The classical treatment of phase coarsening was done by Todes, Lifshitz, and Slyozov, (TLS) in the limit of zero volume fraction, V(sub v), of the dispersed phase. Since the publication of TLS theory there have been numerous investigations, many of which sought to describe the kinetic scaling behavior over a 0 range of volume fractions. Some studies in the literature report that the relative increase in coarsening rate at low (but not zero) volume fractions compared to that predicted by TLS is proportional to v(sub v)(exp 1/2), whereas others suggcest V(sub v)(exp 1/3). This issue has been resolved recently by simulation studies at low volume fractions in three dimensions by members of the Rensselaer/MSFC team. Our studies of ripening behavior using large-scale numerical simulations suggest that although there are different circumstances which can lead to either scaling law, the most important length scale at low volume fractions is the diffusional analog of the Debye screening length. The numerical simulations we employed exploit the use of a recently developed "snapshot" technique, and identifies the nature of the coarsening dynamics at various volume fractions. Preliminary results of numerical and experimental investigations, focused on the growth of finite particle clusters, provide important insight into the nature of the transition between the two scaling regimes. The companion microgravity experiment centers on the growth within finite particle clusters, and follows the temporal dynamics driving microstructural evolution, using holography.

Frazier, Donald O.↗
NASA NTRS1997

The Stability of Radiatively Cooling Jets I. Linear Analysis

The results of a spatial stability analysis of a two-dimensional slab jet, in which optically thin radiative cooling is dynamically important, are presented. We study both magnetized and unmagnetized jets at external Mach numbers of 5 and 20. We model the cooling rate by using two different cooling curves: one appropriate to interstellar gas, and the other to photoionized gas of reduced metallicity. Thus, our results will be applicable to both protostellar (Herbig-Haro) jets and optical jets from active galactic nuclei. We present analytical solutions to the dispersion relations in useful limits and solve the dispersion relations numerically over a broad range of perturbation frequencies. We find that the growth rates and wavelengths of the unstable Kelvin-Helmholtz (K-H) modes are significantly different from the adiabatic limit, and that the form of the cooling function strongly affects the results. In particular, if the cooling curve is a steep function of temperature in the neighborhood of the equilibrium state, then the growth of K-H modes is reduced relative to the adiabatic jet. On the other hand, if the cooling curve is a shallow function of temperature, then the growth of K-H modes can be enhanced relative to the adiabatic jet by the increase in cooling relative to heating in overdense regions. Inclusion of a dynamically important magnetic field does not strongly modify the important differences between an adiabatic jet and a cooling jet, provided the jet is highly supermagnetosonic and not magnetic pressure-dominated. In the latter case, the unstable modes behave more like the transmagnetosonic magnetic pressure-dominated adiabatic limit. We also plot fluid displacement surfaces associated with the various waves in a cooling jet in order to predict the structures that might arise in the nonlinear regime. This analysis predicts that low-frequency surface waves and the lowest order body modes will be the most effective at producing observable features in the jet.

Hardee, Philip E.↗
NASA NTRS1995

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗