Search NASA⌕ Search

SEARCH · Search NASA

Results for “convergent evolution”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Forced Reconnection in the Near Magnetotail: Onset and Energy Conversion in PIC and MHD Simulations

Using two-dimensional particle-in-cell (PIC) together with magnetohydrodynamic (MHD) Q1 simulations of magnetotail dynamics, we investigate the evolution toward onset of reconnection and the subsequent energy transfer and conversion. In either case, reconnection onset is preceded by a driven phase, during which magnetic flux is added to the tail at the high-latitude boundaries, followed by a relaxation phase, during which the configuration continues to respond to the driving. The boundary deformation leads to the formation of thin embedded current sheets, which are bifurcated in the near tail, converging to a single sheet farther out in the MHD simulations. The thin current sheets in the PIC simulation are carried by electrons and are associated with a strong perpendicular electrostatic field, which may provide a connection to parallel potentials and auroral arcs and an ionospheric signal even prior to the onset of reconnection. The PIC simulation very well satisfies integral entropy conservation (intrinsic to ideal MHD) during this phase, supporting ideal ballooning stability. Eventually, the current intensification leads to the onset of reconnection, the formation and ejection of a plasmoid, and a collapse of the inner tail. The earthward flow shows the characteristics of a dipolarization front: enhancement of Bz, associated with a thin vertical electron current sheet in the PIC simulation. Both MHD and PIC simulations show a dominance of energy conversion from incoming Poynting flux to outgoing enthalpy flux, resulting in heating of the inner tail. Localized Joule dissipation plays only a minor role.

two-dimensional↗
NASA NTRS2020

Entropy Stable h/p-Nonconforming Discretization with the Summation-by-Parts Property for the Compressible Euler and Navier–Stokes Equations

In this paper, we extend the entropy conservative/stable algorithms presented by Del Rey Fernandez and coauthors for the compressible Euler and Navier-Stokes equations on nonconforming p-refined/coarsened curvilinear grids to h/p refinement/coarsening. The main difficulty in developing nonconforming algorithms is the construction of appropriate coupling procedures across nonconforming interfaces. Here, we utilize a computationally simple and efficient approach based upon using decoupled interpolation operators. The resulting scheme is entropy conservative/stable and elementwise conservative. Numerical simulations of the isentropic vortex and viscous shock propagation con firm the entropy conservation/stability and accuracy properties of the method (achieving ~ p + 1 convergence), which are comparable to those of the original conforming scheme. Simulations of the Taylor{Green vortex at R(e) = 1,600 and turbulent flow past a sphere at R(e(infinity)) = 2,000 show the robustness and stability properties of the overall spatial discretization for unstructured grids. Finally, to demonstrate the entropy conservation property of a fully-discrete explicit entropy stable algorithm with h=p refinement/coarsening, we present the time evolution of the entropy function obtained by simulating the propagation of the isentropic vortex using a relaxation Runge-Kutta scheme.

Nonconforming interfaces↗
DOE OSTI2024

Non-locality of mean scalar transport in two-dimensional Rayleigh–Taylor instability using the macroscopic forcing method

The importance of non-locality of mean scalar transport in two-dimensional Rayleigh–Taylor Instability (RTI) is investigated. The macroscopic forcing method is utilized to measure spatio-temporal moments of the eddy diffusivity kernel representing passive scalar transport in the ensemble averaged fields. Presented in this work are several studies assessing the importance of the higher-order moments of the eddy diffusivity, which contain information about non-locality, in models for RTI. First, it is demonstrated through a comparison of leading-order models that a purely local eddy diffusivity is insufficient to capture the mean field evolution of the mass fraction in RTI. Therefore, higher-order moments of the eddy diffusivity operator are not negligible. Models are then constructed by utilizing the measured higher-order moments. It is demonstrated that an explicit operator based on the Kramers–Moyal expansion of the eddy diffusivity kernel is insufficient. An implicit operator construction that matches the measured moments is shown to offer improvements relative to the local model in a converging fashion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1986

Desingularization of periodic vortex sheet roll-up

An analytical approach is used in an attempt to model the evolution of a vortex sheet past the critical time by means of a desingularization method. Evolution of the sheet, which is embedded in a two-dimensional flow, is described by approximating the total circulation between a fixed material point and an imaginary point on the curve of the flow. A linear dispersion relation is defined which shows that the short wavelength modes of the sheet are not unstable and therefore do adversely effect the computations as artifacts of round-off errors. The desingularization approach is demonstrated to converge beyond the critical time for the vortex sheet. Application of the technique for the study of the vortex sheet shed from an elliptically loaded wing is indicated.

Krasny, R.↗
NASA NTRS1989

The role of diabatic heating in maintaining the upper-tropospheric baroclinic zone in the South Pacific

The four-dimensional structure of the region in the South Pacific convergence zone (SPCZ) during January 10-18, 1979 is studied using a modified set of ECMWF FGGE level III-b analyses. The effects of kinematic and thermodynamic processes on maintaining the upper baroclinic region corresponding with the subtropical jet are analyzed. The role of adiabatic and diabatic processes in the maintenance of baroclinicity in the area of the SPCZ is examined using the frontogenetical function. It is observed that the processes affecting the evolution of baroclinicity in the SPCZ region differ from those in the middle latitude; possible reasons for these differences are proposed.

Robertson, Franklin R.↗
DOE OSTI2024

Exploring the effect of ELM and code-coupling frequencies on plasma and material modeling of dynamic recycling in divertors

Abstract Integrated modeling of plasma-surface interactions provides a comprehensive and self-consistent description of the system, moving the field closer to developing predictive and design capabilities for plasma facing components. One such workflow, including descriptions for the scrape-off-layer plasma, ion-surface interactions and the sub-surface evolution, was previously used to address steady-state scenarios and has recently been extended to incorporate time-dependence and two-way information flow. The new model can address dynamic recycling in transient scenarios, such as the application presented in this paper: the evolution of W samples pre-damaged by helium and exposed to ELMy H-mode plasmas in the DIII-D DiMES. A first set of simulations explored the effect of ELM frequency. This study was discussed in detail in this conference’s proceedings and is summarized here. The 2nd set of simulations, which is the focus of this paper, explores the effect of code-coupling frequency. These simulations include initial SOLPS solutions converged to the inter-ELM state, ion impact energy ( E in ) and angles ( A in ) calculated by hPIC2, and an improved heat transfer description in Xolotl. The model predicts increases in particle fluxes and decreases in heat fluxes by 10%–20% with the coupling time-step. Compared with the first set of simulations, the less shallow impact angle leads to smaller reflection rates and significant D implantation. The higher fraction of implanted flux (and deeper), in particular during ELMs, increases the accumulated D content in the W near-surface region. Future expansion of the workflow includes coupling to hPIC2 and GITR to ensure accurate descriptions of E in and A in , and W impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1989

Determination of solar flare accelerated ion angular distributions from SMM gamma ray and neutron measurements and determination of the He-3/H ratio in the solar photosphere from SMM gamma ray measurements

Comparisons of Solar Maximum Mission (SMM) observations of gamma-ray line and neutron emission with theoretical calculation of their expected production by flare accelerated ion interactions in the solar atmosphere have led to significant advances in the understanding of solar flare particle acceleration and interaction, as well as the flare process itself. These comparisons have enabled the determination of, not only the total number and energy spectrum of accelerated ions trapped at the sun, but also the ion angular distribution as they interact in the solar atmosphere. The Monte Carlo program was modified to include in the calculations of ion trajectories the effects of both mirroring in converging magnetic fields and of pitch angle scattering. Comparing the results of these calculations with the SMM observations, not only the angular distribution of the interacting ions can be determined, but also the initial angular distribution of the ions at acceleration. The reliable determination of the solar photospheric He-3 abundance is of great importance for understanding nucleosynthesis in the early universe and its implications for cosmology, as well as for the study of the evolution of the sun. It is also essential for the determinations of the spectrum and total number of flare accelerated ions from the SMM/GRS gamma-ray line measurements. Systematic Monte Carlo calculations of the time dependence were made as a function of the He-3 abundance and other variables. A new series of calculations were compared for the time-dependent flux of 2.223 MeV neutron capture line emission and the ratio of the time-integrated flux in the 2.223 MeV line to that in the 4.1 to 6.4 MeV nuclear deexcitation band.

Lingenfelter, Richard E.↗
NASA NTRS1993

Thermodynamic balance of three-dimensional stratospheric winds derived from a data assimilation procedure

The NASA/Goddard three-dimensional chemistry and transport model is driven by winds from a stratospheric data assimilation system. Synoptic- and planetary-scale patterns, apparent in satellite observations of trace constituents, are successfully reproduced for seasonal integrations. As model integrations proceed, however, the quality of simulations decreases, and systematic differences between calculation and measurement appear. The differences are explained by examining the zonal-mean residual circulation. The vertical velocity w-bar (sup star) is calculated two ways: (1) from the diabatic heating rates and temperature tendency and (2) from the Eulerian vertical velocity and the horizontal eddy heat flux convergence. The results from these calculations differ substantially. Periodic insertion of observational data during the assimilation process continually shocks the general circulation model and produces these differences, which leads to an overestimate of the mean vertical heat and cconstituent transport. Such differences are expected to be general to all data assimilation products. This interpretation is corroborated by two-dimensional (2D) model calculations. When w-bar(sup star) is calculated from (2), the 2D ozone evolution is unrealistic and qualitatively similar to the 3D model simulation. The 2D ozone evolution is reasonable when w-bar (sup star) from (1).

Weaver, Clark J.↗
NASA NTRS1980

Theoretical studies of massive stars. II - Evolution of a 15 solar-mass star from carbon shell burning to iron core collapse

The evolution of a Population I star of 15 solar masses is described from the carbon shell burning stage to the formation and collapse of an iron core. An unusual aspect of the evolution is that neon ignition occurs off-center and neon burning propagates inward by a series of shell flashes. The extent of the core burning is generally smaller than the Chandrasekhar mass, so that most of the nuclear energy generation occurs in shell sources. Because of degeneracy and the influence of rapid convective mixing, these shell sources are unstable and the core goes through large excursions in temperature and density. The small core also causes the shell sources to converge into a narrow mass region slightly above the Chandrasekhar mass. Thus, the final nucleosynthesis yields are generally small, with silicon being most strongly enhanced with respect to solar system abundances.

Sparks, W. M.↗
NASA NTRS2009

Steady-State ALPS for Real-Valued Problems

The two objectives of this paper are to describe a steady-state version of the Age-Layered Population Structure (ALPS) Evolutionary Algorithm (EA) and to compare it against other GAs on real-valued problems. Motivation for this work comes from our previous success in demonstrating that a generational version of ALPS greatly improves search performance on a Genetic Programming problem. In making steady-state ALPS some modifications were made to the method for calculating age and the method for moving individuals up layers. To demonstrate that ALPS works well on real-valued problems we compare it against CMA-ES and Differential Evolution (DE) on five challenging, real-valued functions and on one real-world problem. While CMA-ES and DE outperform ALPS on the two unimodal test functions, ALPS is much better on the three multimodal test problems and on the real-world problem. Further examination shows that, unlike the other GAs, ALPS maintains a genotypically diverse population throughout the entire search process. These findings strongly suggest that the ALPS paradigm is better able to avoid premature convergence then the other GAs.

Hornby, Gregory S.↗
NASA NTRS1993

On the 'delta-equations' for vortex sheet evolution

We use a set of equations, sometimes referred to as the 'delta-equations', to approximate the two-dimensional inviscid motion of an initially circular vortex sheet released from rest in a cross-flow. We present numerical solutions of these equations for the case with delta-square = 0 (for which the equations are exact) and for delta-square greater than 0. For small values of the smoothing parameter delta, a spectral filter must be used to eliminate spurious instabilities due to round-off error. Two singularities appear simultaneously in the vortex sheet when delta-square = 0 at a critical time t(c). After t(c), the solutions do not converge as the computational mesh is refined. With delta-square greater than 0, converged solutions were found for all values of delta-square when t is less than t(c), and for all but the two smallest values of delta-square used when t is greater than t(c). Our results show that, when delta-square is greater than 0, the vortex sheet deforms into two doubly branched spirals some time after t(c). The limiting solution as delta approaching 0 clearly exists and equals the delta = 0 solution when t is less than t(c).

Rottman, James W.↗
NASA NTRS2000

Coupled Model Simulation of Snowfall Events Over the Black Hills

Although many long-term simulations of snow accumulation and oblation have been made using stand-alone land surface models and surface models coupled with GCMs, less research has focused on short-term event simulations. Actually, accurate event simulations of snow-related processes are the basis for successful long-term simulation. Three advantages of event simulations of snowfall and snow melting are availability of: (1) intensive observation data from field experiments for validation; (2) more physically-realistic precipitation schemes for use in atmospheric models to simulate snowfall; and (3) a more detailed analysis of the snow melting processes. In addition to the complexities of snow related processes themselves, terrain-induced effects on snowfall/snow melting make simulations of snow events more difficult. Climatological observations indicate that terrain features such as the Black Hills of South Dakota and Wyoming can exert important effects on snow accumulation and snow oblation processes. One of the primary effects is that the orography causes forced uplift of airflow and causes atmospheric waves to form both upwind and downwind of it. Airflow often splits around the obstacle, converging on the lee side. This convergence may lead to precipitation enhancement. It also provides an elevated heat and moisture source that enhances atmospheric instability. During the period of April 5-May 5, 1999, the Upper Missouri River Basin Pilot Project (UMRBPP) made intensive observations on precipitation events occurring in the Black Hills. Two moderate snowfall events were captured during the period. The resulting high temporal and spatial resolution data provides opportunities to investigate terrain effects on snowfall amount, distribution, and melting. Successful simulation of snowfall amount, distribution, and evolution using atmospheric models is important to subsequent modeling of snow melting using snow sub-models in land surface schemes. In this paper, a coupled model system, consisting of an atmosphere model (ARPS) and a land-surface model (revised NCAR LSM), is used to simulate one of these UMRBPP cases.

Wang, Jianzhong↗
DOE OSTI2024

Bipartite mutual information in classical many-body dynamics

Information theoretic measures have helped to sharpen our understanding of many-body quantum states. As perhaps the most well-known example, the entanglement entropy (or more generally, the bipartite mutual information) has become a powerful tool for characterizing the dynamical growth of quantum correlations. By contrast, although computable, the bipartite mutual information (MI) is almost never explored in classical many particle systems; this owes in part to the fact that computing the MI requires keeping track of the evolution of the full probability distribution, a feat which is rarely done (or thought to be needed) in classical many-body simulations. Here, we utilize the MI to analyze the spreading of information in 1D elementary cellular automata (CA). Broadly speaking, we find that the behavior of the MI in these dynamical systems exhibits a few different types of scaling that roughly correspond to known CA universality classes. Of particular note is that we observe a set of automata for which the MI converges parametrically slowly to its thermodynamic value. We develop a microscopic understanding of this behavior by analyzing a two-species model of annihilating particles moving in opposite directions. Furthermore, our work suggests the possibility that information theoretic tools such as the MI might enable a more fine-grained characterization of classical many-body states and dynamics.

Cellular automata↗
NASA NTRS2000

The Influence of Soil Moisture, Coastline Curvature, and Land-Breeze Circulations on Sea-Breeze Initiated Precipitation

Idealized numerical simulations are performed with a coupled atmosphere/land-surface model to identify the roles of initial soil moisture, coastline curvature, and land breeze circulations on sea breeze initiated precipitation. Data collected on 27 July 1991 during the Convection and Precipitation Electrification Experiment (CAPE) in central Florida are used. The 3D Goddard Cumulus Ensemble (GCE) cloud resolving model is coupled with the Goddard Parameterization for Land-Atmosphere-Cloud Exchange (PLACE) land surface model, thus providing a tool to simulate more realistically land-surface/atmosphere interaction and convective initiation. Eight simulations are conducted with either straight or curved coast-lines, initially homogeneous soil moisture or initially variable soil moisture, and initially homogeneous horizontal winds or initially variable horizontal winds (land breezes). All model simulations capture the diurnal evolution and general distribution of sea-breeze initiated precipitation over central Florida. The distribution of initial soil moisture influences the timing, intensity and location of subsequent precipitation. Soil moisture acts as a moisture source for the atmosphere, increases the connectively available potential energy, and thus preferentially focuses heavy precipitation over existing wet soil. Strong soil moisture-induced mesoscale circulations are not evident in these simulations. Coastline curvature has a major impact on the timing and location of precipitation. Earlier low-level convergence occurs inland of convex coastlines, and subsequent precipitation occurs earlier in simulations with curved coastlines. The presence of initial land breezes alone has little impact on subsequent precipitation. however, simulations with both coastline curvature and initial land breezes produce significantly larger peak rain rates due to nonlinear interactions.

Baker, David R.↗
DOE OSTI2024

Environmental controls on isolated convection during the Amazonian wet season

The Amazon rainforest is a vital component of the global climate system, influencing the hydrological cycle and tropical circulation. However, understanding and modeling the evolution of convection in this region remain a scientific challenge. Here, we assess the environmental conditions associated with shallow, congestus, and isolated deep convection days during the wet season (December to April), employing measurements from the Green Ocean Amazon 2014–2015 (GoAmazon2014/5) experiment and large-scale wind fields from the constrained variational analysis. Composites of deep days show moister than average conditions below 3 km early in the morning. Analyzing the water budget at the surface through observations only, we estimated the water vapor convergence term as a residual of the water balance closure. Convergence remains nearly zero during the deep days until early afternoon (13:00 LST), when it becomes a dominant factor in the water budget. At 14:00 LST, the deep days experience a robust upward large-scale vertical velocity, especially above 4 km, which supports the shallow-to-deep convective transition occurring around 16:00–17:00 LST. In contrast, shallow and congestus days exhibit drier pre-convective conditions, along with diurnal water vapor divergence and large-scale subsidence that extend from the surface to the lower free troposphere. Moreover, afternoon precipitation exhibits the strongest linear correlation (0.6) with large-scale vertical velocity, nearly double the magnitude observed for other environmental factors, even moisture, at different levels and periods of the day. Precipitation also exhibits a moderate increase with low-level wind shear, while upper-level shear has a relatively minor negative impact on convection.

Environmental Sciences & Ecology↗
NASA NTRS1992

Essentially nonoscillatory postprocessing filtering methods

High order accurate centered flux approximations used in the computation of numerical solutions to nonlinear partial differential equations produce large oscillations in regions of sharp transitions. Here, we present a new class of filtering methods denoted by Essentially Nonoscillatory Least Squares (ENOLS), which constructs an upgraded filtered solution that is close to the physically correct weak solution of the original evolution equation. Our method relies on the evaluation of a least squares polynomial approximation to oscillatory data using a set of points which is determined via the ENO network. Numerical results are given in one and two space dimensions for both scalar and systems of hyperbolic conservation laws. Computational running time, efficiency, and robustness of method are illustrated in various examples such as Riemann initial data for both Burgers' and Euler's equations of gas dynamics. In all standard cases, the filtered solution appears to converge numerically to the correct solution of the original problem. Some interesting results based on nonstandard central difference schemes, which exactly preserve entropy, and have been recently shown generally not to be weakly convergent to a solution of the conservation law, are also obtained using our filters.

Lafon, F.↗
NASA NTRS2005

The Cycle of Cosmic Dust

The presence of dust in the interstellar medium of galaxies dramatically affects their spectral appearance, and the interpretation of their underlying physical properties. Consequently, the diagnostic of galaxy spectra depends crucially on our understanding of the nature and properties of these dust particles. Models of interstellar dust particles must be able to reproduce a basic set of observational constraints, including the general interstellar extinction and infrared emission observed in the diffuse interstellar medium (ISM). Recent analysis of the solar spectrum have resulted in a convergence between the solar and B star abundances. This development, and the steadiness in the solar abundance determination of the primary refractory elements Mg, Si, and Fe, strongly suggest that any viable dust model must also obey interstellar abundances constraints. I will present fifteen dust models that differ in composition and size distribution and that simultaneously satisfy the local extinction, infrared emission, and abundances constraints have been shown to exist. This multitude of viable dust models provides us with an increased flexibility in understanding dust evolution and the many variations in dust properties in different phases of the ISM and stellar environments.

Dwek, Eliahu↗
NASA NTRS2002

Photoevaporation of Clumps in Photodissociation Regions

We present the results of an investigation of the effects of Far Ultraviolet (FUV) radiation (6.0eV < hv < 13.6eV) from hot early type OB stars on clumps in star-forming molecular clouds. Clumps in FUV-illuminated regions (or photodissociation regions or PDRs) undergo external heating and photodissociation as they are exposed to the FUV field, resulting in a loss of cold, molecular lump mass as it is converted to warm atomic gas. The heating, if rapid, creates strong photoevaporative mass flows off the clump surfaces, and drives shocks into the clumps, compressing them to high densities. The clumps lose mass on relatively short timescales. The evolution of an individual clump is found to be sensitive to three dimensionless parameters: Nc0, the ratio of the initial column density of the clump to the column N(0) approx. 10(exp 21) cm(exp -2) of a warm FUV-heated surface region; upsilon, the ratio of the sound speed in the heated surface to that in the cold clump material: and t(FUV)t(c), the ratio of the "turn-on time" t(FUV) of the heating flux on a clump to its initial sound crossing-time t(c). The evolution also depends on whether a confining interclump medium exists, or whether the interclump region has negligible pressure, as is the case for turbulence-generated clumps. In this paper, we use spherical 1-D numerical hydrodynamic models as well as approximate analytical models to study the dependence of clump photoevaporation on the physical parameters of the clump, and to derive the dynamical evolution, mass loss rates and photoevaporative timescales of a clump for a variety of astrophysical situations. Turbulent clumps evolve so that their column densities are equal to a critical value determined by the local FUV field, and typically have short photo evaporation timescales, approx. 10(exp 4-5) years for a 1 M(solar mass) clump in a typical star-forming region (Nc0 = 10, upsilon = 10). Clumps with insufficient magnetic pressure support, and in strong FUV fields may be driven to collapse by the compressional effect of converging shock waves. We also estimate the rocket effect on photoevaporating clumps and find that it is significant only for the smallest clumps, with sizes much less than the extent of the PDR itself. Clumps that are confined by all interclump medium may either get completely photoevaporated, or may preserve a shielded core with a warm, dissociated, protective shell that, absorbs the incident FUV flux. We compare our results with observations of some well studied PDRs: the Orion Bar, M17SW NGC 2023 and the Rosette Nebula. The data are consistent with both interpretations of clump origin. turbulence and pressure confinement, with a slight indication for favouring the turbulent model for clumps over pressure-confined

Gorti, Uma↗