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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Accurately simulating core-collapse self-interacting dark matter halos

The properties of satellite halos provide a promising probe for dark matter (DM) physics. Observations have motivated current efforts to explain surprisingly compact DM halos. If DM is not collisionless, but has strong self-interactions, halos can undergo gravothermal collapse, leading to higher densities in the central region of the halo. However, it is challenging to model this collapse phase from first principles. To improve on this, we sought to better understand the numerical challenges and convergence properties of self-interacting dark matter (SIDM) N-body simulations in the collapse phase. Especially, our aim was to better understand the evolution of satellite halos. To do so, we ran SIDM N-body simulations of a low-mass halo in isolation and within an external gravitational potential. The simulation set-up was motivated by the perturber of the stellar stream GD-1. We find that the halo evolution is very sensitive to energy conservation errors, and a SIDM kernel size that is too large can artificially speed up the collapse. Moreover, we demonstrate that the King model can describe the density profile at small radii for the late stages that we have simulated. Furthermore, for our most highly resolved simulation (N = 5 × 10 7 ) we have made the data public. It can serve as a benchmark. Overall, we find that the current numerical methods do not suffer from convergence problems in the late collapse phase and provide guidance on how to choose numerical parameters, for example that the energy conservation error is better kept well below 1%. This allows simulations to be run of halos that become concentrated enough to explain observations of GD-1-like stellar streams or strong gravitational lensing systems.

dark matter↗

NASA NTRS1990

Numerical Formulation for the Prediction of Solid/Liquid Change of a Binary Alloy

A computational model is presented for the prediction of solid/liquid phase change energy transport including the influence of free convection fluid flow in the liquid phase region. The computational model considers the velocity components of all non-liquid phase change material control volumes to be zero but fully solves the coupled mass-momentum problem within the liquid region. The thermal energy model includes the entire domain and uses an enthalpy like model and a recently developed method for handling the phase change interface nonlinearity. Convergence studies are performed and comparisons made with experimental data for two different problem specifications. The convergence studies indicate that grid independence was achieved and the comparison with experimental data indicates excellent quantitative prediction of the melt fraction evolution. Qualitative data is also provided in the form of velocity vector diagrams and isotherm plots for selected times in the evolution of both problems. The computational costs incurred are quite low by comparison with previous efforts on solving these problems.

Schneider, G. E.↗

DOE OSTI2024

Virtual Tokamak for Test and Development of Plasma Control Applied to NSTX-U

Safe and efficient operation of tokamak experiments depends on shot preparation prior to the experiment to maximize performance and avoid issues that are predictable consequences of known physics. The General Atomics TokSys toolbox is used to test and develop the plasma control system (PCS) on numerous tokamaks around the world. Tokamaks that use a version of the DIII-D PCS can connect it to a TokSys simulation and control a virtual version of the tokamak with the real PCS. Recent upgrades to the simulation include more detailed profile and transport modeling to support design and implementation of optimized scenarios. The simulation combines a module called “Profiles” that simulates the 1D profiles of density, pressure, and current and a module called “GSevolve'' that evolves the 2D Grad-Shafranov equilibrium. The “Profiles'' module primarily uses simplified models for fueling, heating, and current. The TRANSP code provides precomputed diffusion coefficients and heating/current drive profiles from radio-frequency and neutral beam sources for the scenarios being simulated. The fidelity of a simulation can be enhanced by analyzing it with TRANSP and making corrections to the precomputed values and iterating until convergence. Simulations of the National Spherical Torus Upgrade (NSTX-U) are presented. It is shown that the simulation reproduces the experimental profile evolution. Here, the simulation has also been upgraded to reproduce the experimental outcome at points of bifurcation. This is demonstrated with a simulation of vertical displacement in NSTX-U.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1986

The kinetic and available potential energy budget of a winter extratropical cyclone system

The energy budget of an extratropical cyclone system which traversed North America and intensified through the period January 9-11, 1975 is presented. The objectives of the study are: (1) to document the complete energy budget of a significant winter cyclone event, and (2) to comment on the significance of latent heat release (LHR) in the cyclone's evolution. Results reveal an overall increase in both kinetic (K) and available potential energy (A). K increases are accounted for by boundary flux convergence of K, while A increases are due to generation by LHR and K to A conversion. In addition, the general A increase is accompanied by a 24 h oscillation that is explained largely by the flux quantity in the A budget equation and is correlated with a similar fluctuation in the K to A conversion. LHR does not appear to be critical in the development of this cyclone system. Rather, LHR acts to increase the intensity of the event. It is hypothesized that the direct influence that LHR had on the deepening cyclone's reduced mass was augmented by an indirect influence, in which pre-existing dry dynamical forcing was enhanced by diabatic heating, thus leading to accelerated cyclone development at a later time.

Smith, P. J.↗

NASA NTRS2008

Formation and Evolution of Lakshmi Planum (V-7), Venus: Assessment of Models using Observations from Geological Mapping

Lakshmi Planum is a high-standing plateau (3.5-4.5 km above MPR) surrounded by the highest mountain ranges on Venus. Lakshmi represents a unique type of elevated region different from dome-shaped and rifted rises and tessera-bearing crustal plateaus. The unique characteristics of Lakshmi suggest that it formed by an unusual combination of processes and played an important role in Venus geologic history. Lakshmi was studied with Venera-15/16 and Magellan data, resulting in two classes of models, divergent and convergent, to explain its unusual topographic and morphologic characteristics. Divergent models explain Lakshmi as a site of mantle upwelling due to rising and subsequent collapse of a mantle diapir; such models explain emplacement of a lava plateau inside Lakshmi and, in some circumstances, formation of the mountain ranges. The convergent models consider Lakshmi as a locus of mantle downwelling, convergence, underthrusting, and possible subduction. Key features in these models are the mountain ranges, high topography of Lakshmi interior, and the large volcanic centers in the plateau center. These divergent and convergent models entail principally different mechanisms of formation and suggest different geodynamic regimes on Venus. Almost all models make either explicit or implicit predictions about the type and sequence of major events during formation and evolution of Lakshmi and thus detailed geological mapping can be used to test them. Here we present the results of such geological mapping (the V-7 quadrangle, 50-75degN, 300-360degE; scale 1:5M) that allows testing the proposed models for Lakshmi.

Ivanov, M. A.↗

NASA NTRS1989

Parameter identification for an abstract Cauchy problem by quasilinearization

A parameter identification problem is considered in the context of a linear abstract Cauchy problem with a parameter-dependent evolution operator. Conditions are investigated under which the gradient of the state with respect to a parameter possesses smoothness properties which lead to local convergence of an estimation algorithm based on quasi-linearization. Numerical results are presented concerning estimation of unknown parameters in delay-differential equations.

Brewer, Dennis W.↗

NASA NTRS1993

Parameter identification for an abstract Cauchy problem by quasilinearization

Brewer, Dennis W.↗

NASA NTRS1992

Evolution of solar wind fluctuations and the influence of turbulent 'mixing'

We present various numerical and analytical solutions for the transport of solar wind turbulence. The model used takes into account the effects of convection, expansion, and wave propagation, as well as the recently illuminated effects of (non-WKB) 'mixing' terms. The radial evolution of the fluctuating kinetic energy, magnetic energy, and normalized cross helicity is computed, and, it is demonstrated that in appropriate limits the solutions converge to the WKB forms. In the general case, solutions which differ substantially from those predicted by WKB theory are obtained. The degree of turbulent 'mixing' shows considerable dependence on the nature of the turbulence, giving rise to varying levels, at 1 AU, of the ratio of 'inward' and 'outward' fluctuation energies and the ratio of kinetic to magnetic energies in the fluctuations. The transport properties described here may provide at least a partial explanation for the observed mixing of cross helicities with increasing heliocentric distance in the solar wind.

Oughton, S.↗

NASA NTRS1999

Dynamical Signatures of Living Systems

One of the main challenges in modeling living systems is to distinguish a random walk of physical origin (for instance, Brownian motions) from those of biological origin and that will constitute the starting point of the proposed approach. As conjectured, the biological random walk must be nonlinear. Indeed, any stochastic Markov process can be described by linear Fokker-Planck equation (or its discretized version), only that type of process has been observed in the inanimate world. However, all such processes always converge to a stable (ergodic or periodic) state, i.e., to the states of a lower complexity and high entropy. At the same time, the evolution of living systems directed toward a higher level of complexity if complexity is associated with a number of structural variations. The simplest way to mimic such a tendency is to incorporate a nonlinearity into the random walk; then the probability evolution will attain the features of diffusion equation: the formation and dissipation of shock waves initiated by small shallow wave disturbances. As a result, the evolution never "dies:" it produces new different configurations which are accompanied by an increase or decrease of entropy (the decrease takes place during formation of shock waves, the increase-during their dissipation). In other words, the evolution can be directed "against the second law of thermodynamics" by forming patterns outside of equilibrium in the probability space. Due to that, a specie is not locked up in a certain pattern of behavior: it still can perform a variety of motions, and only the statistics of these motions is constrained by this pattern. It should be emphasized that such a "twist" is based upon the concept of reflection, i.e., the existence of the self-image (adopted from psychology). The model consists of a generator of stochastic processes which represents the motor dynamics in the form of nonlinear random walks, and a simulator of the nonlinear version of the diffusion equation which represents the mental dynamics. It has been demonstrated that coupled mental-motor dynamics can simulate emerging self-organization, prey-predator games, collaboration and competition, "collective brain," etc.

Zak, M.↗

NASA NTRS2005

WRF Simulation of the Genesis of Hurricane Javier (2004) in the Eastern Pacific

The Eastern Pacific has the highest frequency of genesis events per unit area of any region worldwide (Elsberry et al 1987). African easterly waves, mesoscale convective systems (MCSs), and topographic effects are thought to play roles in the genesis of tropical cyclones there (Frank and Clark 1980, Velasco and Fritsch 1987, Zehnder 1991, Zehnder and Gall 1991; Farfan and Zehnder 1997). Mozer and Zehnder (1996), using dry, idealized simulations of flow past a large-scale three-dimensional mountain range comparable to the Sierra Madre Mountains of Mexico, showed that upstream flow blocking led to diversion of the flow primarily to the south of the mountains. This flow diversion led to the formation of a low-level, barotropically unstable jet (at a location comparable to the Isthmus of Tehuantepec) and the continuous formation of synoptic-scale vorticity maxima, which they suggested may play a role in tropical cyclogenesis. Farfan and Zehnder (1 997) examined the synoptic-scale circulations that led to the formation of Hurricane Guillermo (1991). Using numerical simulations, they found that flow blocking led to the formation of a low-level easterly jet south of the mountains of Central America and a northeasterly (gap flow) jet over the Gulf of Tehuantepec, which combined with the flow associated with the Intertropical Convergence Zone (ITCZ) to produce a closed cyclonic circulation in the location of Guillermo s formation. As will be discussed in this paper, the evolution of the flow field that was associated with the genesis of Hurricane Javier was similar to that described in Farfan and Zehnder (1997), with well-defined topographic flow features. Here, using a high- resolution simulation with the WRF model, we investigate whether these topographically induced flows played a significant role in the genesis of Javier.

Braun, Scott A.↗

DOE OSTI2026

JIMWLK on a quantum computer

We propose a method for solving the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) evolution equation on quantum computers. Our approach exploits the reformulation of the JIMWLK equation as a Lindblad master equation governing the rapidity evolution of the hadronic density matrix, as established in prior work. To render the problem tractable for quantum simulation, we introduce several approximations: the two-dimensional transverse plane is reduced to a one-dimensional radial lattice by assuming azimuthal symmetry of the jump operators; the gauge group is restricted to SU(2); and the infinite Wilson lines of the JIMWLK equation are replaced by finite Wilson links along the light-cone direction. The resulting bosonic Hilbert space is truncated using the electric field basis familiar from Hamiltonian lattice gauge theory, with states restricted to angular momenta 𝑗 ≤ 𝑗 max . We derive the matrix elements of the JIMWLK Lindblad jump operators in this basis. As a benchmark, we demonstrate rapid convergence of the fundamental dipole expectation value with 𝑗 max for both pure and mixed Gaussian initial density matrices. For the simplest truncation, 𝑗 max =1/2, we implement the Lindblad evolution using a quantum simulation algorithm verified with the Qiskit statevector simulator by decomposing the non-unitary evolution operator into a linear combination of unitaries. This work establishes a concrete pathway toward quantum simulation of high-energy quantum chromodynamics evolution equations, with direct relevance to the physics program of the Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1987

The effect of along-front temperature variation in a two-dimensional primitive equation model of a surface frontogenesis

The effects of along-front temperature variations on the structure and evolution of surface fronts are investigated using a two-dimensional frontogenesis model. The formulation and characteristics of the confluence-horizontal shear model are described. The structures of three frontal zones forced by confluence and characterized by zero, negative, and positive along-front potential temperature variations are analyzed; particular attention is given to the ageostrophic vertical circulation and their forcings. The temporal evolution of the model frontal zones is examined in terms of the low-level maxima of vorticity, cross-front potential temperature gradient, and convergence. The kinematic and dynamic processes influencing the frontogenesis are discussed. It is noted that the two-dimensional model is applicable to the study of the dynamical properties of frontogenesis.

Keyser, Daniel↗

DOE OSTI2026

Novel additive manufacturing for plasma facing materials ‐ creating a research pathway for minority students

This project addresses two critical and intertwined challenges in fusion energy, namely the shortage of a broadly trained scientific workforce and the lack of scalable manufacturing solutions for plasma-facing components (PFCs). Through a collaboration among Florida International University (FIU), Miami Dade College (MDC), and Purdue University, the project established structured, reproducible educational and research pathways that recruit and advance students from institutions historically outside the fusion energy enterprise, building the human capital that this field urgently needs. The project integrates the complementary research strengths of FIU and Purdue to investigate flash sintering as a transformative processing route for tungsten-based PFCs. Unlike conventional sintering approaches, flash sintering offers rapid densification at significantly reduced thermal budgets, making it a compelling candidate for fabricating complex tungsten geometries that must withstand extreme plasma-facing environments. Systematic experimental and modeling efforts will elucidate the fundamental mechanisms governing microstructure evolution, grain boundary chemistry, and thermomechanical response during flash sintering — knowledge that is presently lacking but essential for translating this technology into reliable manufacturing practice. The convergence of workforce development and cutting-edge manufacturing research positions this project to deliver measurable, durable impact: a pipeline of fusion-ready researchers cultivated through expanded institutional partnerships, and a validated materials processing framework that accelerates domestic readiness for next-generation fusion reactor construction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1989

Gradient methods for identification of distributed parameter systems

Parameter-dependence properties are developed in the context of a linear abstract Cauchy problem governed by a parameter-dependent operator. Of particular interest are problems in which the parameter induces an unbounded perturbation of the evolution operator. Conditions are stated under which the derivative of the state with respect to the parameter possesses certain smoothness properties. These properties lead to local convergence results for a gradient-based estimation algorithm based on quasi-linearization. Numerical results concerning a delay-differential equation indicate the effect of a forcing function on parameter sensitivity.

Brewer, Dennis W.↗

NASA NTRS1986

Meso-beta scale perturbations of the wind field by thunderstorm cells

Data from the high density storm-scale rawinsonde network of the Severe Environmental Storms and Mesoscale Experiment revealed temporal and spatial changes in the divergence fields of the troposphere in response to severe storm evolution on May 2, 1979; these changes were detectable on the meso-beta scale. This unique set of data was subsequently used to study the evolution of the wind, divergence and vertical motion fields in the presence of intense convection. Mid- and upper-tropospheric divergence was superimposed over low-level convergence. The divergence, which has a maximum value of .0004/s, occurred 75 to 100 km upwind as well as over the tornadic cells. To the south of the storm cells, the kinematic pattern was in reverse, upper level convergence was superimposed over low-level divergence. A vertical motion doublet was found to ascend over the squall line and descend about 70 km south of the squall line. It is suggested that the following effects are accountable for the nature of the kinematic fields: (1) blocking of tropospheric environmental flow by the storm cells, (2) anvil outflows, particularly from the tornadic cells, and (3) divergence from the exit region of the jet stream.

Ulanski, S. L.↗

NASA NTRS2021

Assimilation of DAWN Doppler Wind Lidar Data During the 2017 Convective Processes Experiment (CPEX): Impact on Precipitation and Flow Structure

An improved representation of 3-D air motion and precipitation structure through forecast models and assimilation of observations is vital for improvements in weather forecasting capabilities. However, there are few independent data to properly validate a model forecast of precipitation structure when the underlying dynamics are evolving on short convective timescales. Using data from the JPL Ku/Ka-band Airborne Precipitation Radar (APR-2) and the 2 μmDoppler Aerosol Wind (DAWN) lidar collected during the2017 Convective Processes Experiment (CPEX), the NASA Unified Weather Research and Forecasting (WRF) Ensemble Data Assimilation System (EDAS) modeling system was used to quantify the impact of high-resolution sparsely sampled DAWN measurements on the analyzed variables and on the forecast when the DAWN winds were assimilated. Over-all, the assimilation of the DAWN wind profiles had a discernible impact on the wind field as well as the evolution and timing of the 3-D precipitation structure. Analysis of individual variables revealed that the assimilation of the DAWN winds resulted in important and coherent modifications of the environment. It led to an increase in the near-surface convergence, temperature, and water vapor, creating more favorable conditions for the development of convection exactly where it was observed (but not present in the control run). Comparison to APR-2 and observations by the Global Precipitation Measurement (GPM) satellite shows a much-improved forecast after the assimilation of the DAWN winds – development of precipitation where there was none, more organized precipitation where there was some, and a much more intense and organized cold pool, similar to the analysis of the dropsonde data. The onset of the vertical evolution of the precipitation showed similar radar-derived cloud-top heights, but delayed in time. While this investigation was limited to a single CPEX flight date, the investigation design is appropriate for further investigation of the impact of airborne Doppler wind lidar observations upon short-term convective precipitation forecasts

DAWN↗

NASA NTRS1994

Auroral plasma transport processes in the presence of kV potential structures

We have simulated plasma transport processes in the presence of a quasi-two-dimensional current filament, that generated kV potential structure in the auroral region. The simulation consists of a set of one-dimensional flux tube simulations with different imposed time-dependent, field-aligned currents. The model uses the 16 moment system of equations and simultaneously solves coupled continuity and momentum equations and equations describing the transport along the magnetic field lines of parallel and perpendicular thermal energy and heat flows for each species. The lower end of the simulation is at an altitude of 800 km, in the collisional topside ionosphere, while the upper end is at 10 R(sub E) in the magnetosphere. The plasma consists of hot electrons and protons of magnetospheric origin and low-energy electrons, protons, and oxygen ions of ionospheric origin. The dynamical interaction of the individual current filaments with ionospheric and magnetospheric plasma generates a potential structure in the horizontal direction and kilovolt field-aligned potential drops along the field lines. The side-by-side display exhibits the evolution of the implied potential structure in the horizontial direction. In the presence of this potential structure and parallel electric field ionospheric plasma density is depleted and velocity is reduced, while density enhancement and increased velocity is observed in magnetospheric plasma. The ionospheric and magnetospheric electron temperatures increase below 2 R(sub E) due to magnetic mirror force on converging geomagnetic field lines. The primary cross-field motion produced by the horizontal E field (E x B drift) is perpendicular to both of the significant spatial directions and is thus ignorable in this geometry. The effects of other cross-field drift processes are discussed. The simulation thus provides insight into the dynamical evolution of two-dimensional potential structures driven by an imposed finite width, field-aligned current profile.

Ganguli, Supriya B.↗

NASA NTRS1992

Evolution of attached and detached slabs and their associated mantle dynamics

Over the two years of the NASA grant, this project has produced a significant amount of research results related to the plate subduction process and the surface crustal deformation at convergent boundaries (i.e., above subduction zones). While some research objectives are completely accomplished, other research tasks remain active and continue to be investigated at present. A steady state analytic thermal model for subducting slabs was used to examine the torques acting on a descending slab. It is found that gravitational torque vanishes when a slab is dipping either vertically or horizontally, unlike previous studies indicating that the magnitude of gravitational torque decreases as dip angle increases. Subsequently, a new time-dependent, analytic thermal model for a subducting slab was developed. The new model enables us to study transient phenomena associated with plate subduction analytically. On the basis of this model, the nature of slab dip angles was evaluated. Slab dip angles are found to be transient features. As they penetrate into the mantle and increase their lengths, the associated gravitational torque also increases resulting in a downward pulling of the slab to the steeper dip angle. This is especially true once a slab penetrates the olivine-spinel phase boundary at about 400 km depth. However, if the phase transformation does not follow the equilibrium condition, the gravitational torque may have a different behavior. This problem was investigated. Except for fast descending slabs, non-equilibrium phase transformation can only slow down the transient increase of slab dip angles discussed earlier. Its effect is not sufficiently strong to reverse the downward pulling for most of the slabs. However, when slabs subducting at 10 cm/yr or faster, a sufficient amount of metastable olivine can exist beneath 400 km. Because of its low density compared with the surrounding spinel, an upward buoyancy is produced resulting in an upward bending of the slab and possibly an upward rotation of the slab such that smaller dip angles are formed. Seismic studies of the Japanese Slab seem to support this interpretation. The development of oroclinal geometries at convergent boundaries was also examined to study plate obduction which is an important ingredient to the initiation of plate subduction. Although the study suggests that surface features are better modeled by block models, the large scale deformation can be adequately studied by viscous models. Such a model is now under development to complete our original objective to study the initiation of plate subduction. Finally, a three-dimensional, finite element, spherical convective model is developed to study dynamic plate subductions. The model development is now complete and it is being tested to ensure its proper operation. The model is able to generate convection results with a viscosity contrast of about 100. Our research continues to push the viscosity contrast to a level that is appropriate for a subducting slab.

Hsui, Albert T.↗