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 199 records · Page 11

Oblique MHD cosmic-ray modified shocks: Two-fluid numerical simulations

We present the first results of time dependent, two-fluid, cosmic-ray (CR) modified, MHD shock simulations. The calculations were carried out with a new numerical code for 1-D ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they appear to have converged to dynamical steady states. We also present results from an initial time dependent simulation which extend the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescale than parallel shocks.

Frank, Adam↗

Oblique magnetohydrodynamic cosmic-ray-modified shocks: Two-fluid numerical simulations

We present the first results of time-dependent two-fluid cosmic-ray (CR) modified magnetohydrodynamic (MHD) shock simulations. The calculations were carried out with a new numerical code for one-dimensional ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks, including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they apppear to have converged to dynamical steady states. We also present results from an initial time-dependent simulation which extends the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescales than parallel shocks.

Frank, Adam↗

Model studies on the role of moist convection as a mechanism for interaction between the mesoscales

A three year research effort is described which had as its goal the development of techniques to improve the numerical prediction of cumulus convection on the meso-beta and meso-gamma scales. Two MESO models are used, the MASS (mesoscale) and TASS (cloud scale) models. The primary meteorological situation studied is the 28-29 Jun. 1986 Cooperative Huntsville Meteorological Experiment (COHMEX) study area on a day with relatively weak large scale forcing. The problem of determining where and when convection should be initiated is considered to be a major problem of current approaches. Assimilation of moisture data from satellite, radar, and surface data is shown to significantly improve mesoscale simulations. The TASS model is shown to reproduce some observed mesoscale features when initialized with 3-D observational data. Convection evolution studies center on comparison of the Kuo and Fritsch-Chappell cumulus parameterization schemes to each other, and to cloud model results. The Fritsch-Chappell scheme is found to be superior at about 30 km resolution, while the Kuo scheme does surprisingly well in simulating convection down to 10 km in cases where convergence features are well-resolved by the model grid. Results from MASS-TASS interaction experiments are presented and discussed. A discussion of the future of convective simulation is given, with the conclusion that significant progress is possible on several fronts in the next few years.

Waight, Kenneth T., III↗

Designing a Unique Single Point Cross Over Method

The idea behind genetic algorithms is to extract optimization strategies nature uses successfully - known as Darwinian Evolution - and transform them for application in mathematical optimization theory to find the global optimum in a defined phase space. One could imagine a population of individual 'explorers' sent into the optimization phase-space. Each explorer is defined by its genes, what means, its position inside the phase-space is coded in his genes. Every explorer has the duty to find a value of the quality of his position in the phase space. (Consider the phase-space being a number of variables in some technological process, the value of quality of any position in the phase space - in other words: any set of the variables - can be expressed by the yield of the desired chemical product.) Then the struggle of 'life' begins. The three fundamental principles are selection, mating/crossover, and mutation. Only explorers (= genes) sitting on the best places will reproduce and create a new population. This is performed in the second step (mating/crossover). The 'hope' behind this part of the algorithm is, that 'good' sections of two parents will be recombined to yet better fitting children. In fact, many of the created children will not be successful (as in biological evolution), but a few children will indeed fulfill this hope. These good sections are named in some publications as building blocks. Now there appears a problem. Repeating these steps, no new area would be explored. The two former steps would only exploit the already known regions in the phase space, which could lead to premature convergence of the algorithm with the consequence of missing the global optimum by exploiting some local optimum. The third step, mutation, ensures the necessary accidental effects. One can imagine the new population being mixed up a little bit to bring some new information into this set of genes. Whereas in biology a gene is described as a macro-molecule with four different bases to code the genetic information, a gene in genetic algorithms is usually defined as a bitstring (a sequence of b 1's and 0's).

Wilson, Richard Phillip↗

Do convergent developmental mechanisms underlie convergent phenotypes?

Convergence is a pervasive evolutionary process, affecting many aspects of phenotype and even genotype. Relatively little is known about convergence in developmental processes, however, nor about the degree to which convergence in development underlies convergence in anatomy. A switch in the ecology of sea urchins from feeding to nonfeeding larvae illustrates how convergence in development can be associated with convergence in anatomy. Comparisons to more distantly related taxa, however, suggest that this association may be limited to relatively close phylogenetic comparisons. Similarities in gene expression during development provide another window into the association between convergence in developmental processes and convergence in anatomy. Several well-studied transcription factors exhibit likely cases of convergent gene expression in distantly related animal phyla. Convergence in regulatory gene expression domains is probably more common than generally acknowledged, and can arise for several different reasons. Copyright 2002 S. Karger AG, Basel.

NASA Discipline Evolutionary Biology↗

First-Principles Simulation of Beam-Induced Processes Underlying Atomic Manipulation in Electron Microscopes

The development of experimental methods and apparatuses capable of promoting atomically precise material manipulations holds great promise for realizing the ultimate limit of feature miniaturization in materials and devices. The ability to modify materials atom by atom is anticipated to usher in new technologies in areas as diverse as separation science, medicine, and quantum information science. Historically, scanning probe-based techniques have been the most prominent approaches in this space. However, these methods are best suited for the manipulation of surface-exposed regions of materials, as the strong perturbations required for bond scission are delivered most effectively to atoms in the near-proximity to the scanning probe. In contrast, convergent electron beams with energies tuned slightly below the threshold for inducing irreversible knock-on damage have recently been employed (within scanning transmission electron microscopy) to promote atomic-scale bond rearrangements in various beam-stable solids. Currently, however, the efficiency and selectivity of beam-induced atomic manipulation processes with focused electron beams are such that long irradiation times are required to induce a desired atomic rearrangement. With a better understanding of the underlying physics dictating the outcome of a given irradiation event, methods can be devised to improve the efficiency of these techniques so that their promise can be fully realized through widespread adoption.To this end, this Account details our recent efforts to develop and apply tractable first-principles simulation approaches for studying the response of materials to electric beam-like external electric potentials applied in real space. We briefly review the concepts and capabilities in the area of atomically precise materials manipulation and review the early demonstrations of accomplishments in this area, focusing on studies using scanned convergent electron beam probes in particular. We expound upon the depth of the challenge and identify critical shortcomings of theoretical methods that have previously been employed in the simulation of beam-induced processes. We then describe the computational methods that we have generalized from the concepts and tools most commonly applied to the study of molecular photochemistry and how our adaptations of these methods can be employed to capture the relevant dynamical phenomena for beam-induced processes ranging from the initial electron scattering to the ensuing multistate reactions. Here, we contextualize these methods within the current state of the art in this area, which has historically focused primarily on the simulation of inelastic image formation in the electron microscope for the purpose of interpreting the results of quantitative electron microscopy experiments. We demonstrate that the spatial distribution of state-specific excitation rates due to the presence of an external (probe) electric charge is inhomogeneous, such that irradiation at particular locations in materials can favor specific electronic transitions (and disallow others). In addition to the potential for excited-state reaction pathways to be accessed through the initial inelastic scattering of the tightly focused electron beam from the targeted atoms, we also identify favorable conditions for the electronically nonadiabatic evolution of the highly vibrationally excited system to open complex multistate reaction pathways. Implications of the early results for understanding the mechanisms and potential routes to improved efficiency and selectivity in beam-induced reactions are discussed. We conclude with a summary of the current state of theory and modeling capabilities in this area and provide our perspective on future directions for theoretical and experimental developments that we view as crucial to advancing the use of convergent electron beams in mode-specific, atomically precise platforms for direct-write materials modifications.

36 MATERIALS SCIENCE↗

Numerical studies of identification in nonlinear distributed parameter systems

An abstract approximation framework and convergence theory for the identification of first and second order nonlinear distributed parameter systems developed previously by the authors and reported on in detail elsewhere are summarized and discussed. The theory is based upon results for systems whose dynamics can be described by monotone operators in Hilbert space and an abstract approximation theorem for the resulting nonlinear evolution system. The application of the theory together with numerical evidence demonstrating the feasibility of the general approach are discussed in the context of the identification of a first order quasi-linear parabolic model for one dimensional heat conduction/mass transport and the identification of a nonlinear dissipation mechanism (i.e., damping) in a second order one dimensional wave equation. Computational and implementational considerations, in particular, with regard to supercomputing, are addressed.

Banks, H. T.↗

Supercomputing Aspects for Simulating Incompressible Flow

The primary objective of this research is to support the design of liquid rocket systems for the Advanced Space Transportation System. Since the space launch systems in the near future are likely to rely on liquid rocket engines, increasing the efficiency and reliability of the engine components is an important task. One of the major problems in the liquid rocket engine is to understand fluid dynamics of fuel and oxidizer flows from the fuel tank to plume. Understanding the flow through the entire turbo-pump geometry through numerical simulation will be of significant value toward design. One of the milestones of this effort is to develop, apply and demonstrate the capability and accuracy of 3D CFD methods as efficient design analysis tools on high performance computer platforms. The development of the Message Passage Interface (MPI) and Multi Level Parallel (MLP) versions of the INS3D code is currently underway. The serial version of INS3D code is a multidimensional incompressible Navier-Stokes solver based on overset grid technology, INS3D-MPI is based on the explicit massage-passing interface across processors and is primarily suited for distributed memory systems. INS3D-MLP is based on multi-level parallel method and is suitable for distributed-shared memory systems. For the entire turbo-pump simulations, moving boundary capability and efficient time-accurate integration methods are built in the flow solver, To handle the geometric complexity and moving boundary problems, an overset grid scheme is incorporated with the solver so that new connectivity data will be obtained at each time step. The Chimera overlapped grid scheme allows subdomains move relative to each other, and provides a great flexibility when the boundary movement creates large displacements. Two numerical procedures, one based on artificial compressibility method and the other pressure projection method, are outlined for obtaining time-accurate solutions of the incompressible Navier-Stokes equations. The performance of the two methods is compared by obtaining unsteady solutions for the evolution of twin vortices behind a flat plate. Calculated results are compared with experimental and other numerical results. For an unsteady flow, which requires small physical time step, the pressure projection method was found to be computationally efficient since it does not require any subiteration procedure. It was observed that the artificial compressibility method requires a fast convergence scheme at each physical time step in order to satisfy the incompressibility condition. This was obtained by using a GMRES-ILU(0) solver in present computations. When a line-relaxation scheme was used, the time accuracy was degraded and time-accurate computations became very expensive.

Kwak, Dochan↗

Time-Accurate Solutions of Incompressible Navier-Stokes Equations for Potential Turbopump Applications

Two numerical procedures, one based on artificial compressibility method and the other pressure projection method, are outlined for obtaining time-accurate solutions of the incompressible Navier-Stokes equations. The performance of the two method are compared by obtaining unsteady solutions for the evolution of twin vortices behind a at plate. Calculated results are compared with experimental and other numerical results. For an un- steady ow which requires small physical time step, pressure projection method was found to be computationally efficient since it does not require any subiterations procedure. It was observed that the artificial compressibility method requires a fast convergence scheme at each physical time step in order to satisfy incompressibility condition. This was obtained by using a GMRES-ILU(0) solver in our computations. When a line-relaxation scheme was used, the time accuracy was degraded and time-accurate computations became very expensive.

Kiris, Cetin↗

2D kinetic-ion simulations of inverted corona fusion targets

Laser-driven “inverted corona” fusion targets have attracted interest as a low-convergence neutron source and platform for studying kinetic physics. The scheme consists of a hollow or gas-filled spherical shell made of deuterated plastic. The shell has one or more laser entrance holes (LEH), resembling a spherical hohlraum. The laser passes through the LEH’s and illuminates the interior surface of the shell, ablating a plasma that travels inward towards the target center. Long ion mean free paths in the converging plasma can lead to significant interpenetration, atomic mix, and other kinetic effects. Here, in this work we report on numerical simulations of inverted corona targets using the kinetic-ion, fluid–electron hybrid particle-in-cell (PIC) approach in 2D RZ geometry. 2D simulations suggest that shape effects do not have a significant impact on plasma evolution and observed yield trends are primarily the result of 1D kinetic mix mechanisms. Simulations are also compared against available experimental data recorded at the OMEGA laser facility. In particular, synthetic x-ray emission images show good qualitative agreement with experimental results, albeit with an apparent timing discrepancy for the two-sided vacuum target. More generally, we demonstrate the potential of hybrid-PIC simulations for full-system modeling and experimental design, including collisional absorption of laser energy, plasma evolution, mix, and fusion burn.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tectonics and volcanism of Eastern Aphrodite Terra: No subduction, no spreading

Eastern Aphrodite Terra is approximately equal in size to the western North American Cordillera, from Mexico to Alaska. Its size and unique landforms make it an important area for understanding the tectonics of Venus, yet models for its formation are diametrically opposed. This region is part of the Equatorial Highlands, which was proposed as a region of lithospheric thinning, isostatic uplift, and attendant volcanism. Eastern Aphrodite Terra is dominated by circular structures within which deformation and volcanism are intimately related. These structures are marked by radial and concentric fractures, and volcanic flows that emanate from a central vent, as well as from concentric fracture sets. Cross-cutting relations between flows and concentric fracture sets indicate that outer concentric fracture sets are younger than inner fracture sets. The circular structures are joined by regional northeast- to east-trending fractures that dominantly postdate formation of the circular structures. We propose that the circular structures 'grow' outward with time. Although these structures probably represent addition of crust to the lithosphere, they do not represent significant lithospheric spreading or convergence, and the region does not mark the boundary between two distinct tectonic plates. This region is not easily explained by analogy with either terrestrial midocean rifts or subduction zones. It is perhaps best explained by upwelling of magma diapirs that blister the surface, but do not cause significant lithospheric spreading. Further study of the structural and volcanic evolution of this region using Magellan altimetry and SAR data should lead to better understanding of the tectonic evolution of this region.

Hansen, Vicki L.↗

MRSR: Rationale for a Mars Rover/Sample Return mission

The Solar System Exploration Committee of the NASA Advisory Council has recommended that a Mars Rover/Sample Return mission be launched before the year 2000. The recommendation is consistent with the science objectives as outlined by the National Academy of Sciences committees on Planetary and Lunar Exploration, and Planetary Biology and Chemical Evolution. Interest has also focused on Mars Rover/Sample Return (MRSR) missions, because of their crucial role as precursors for human exploration. As a result of this consensus among the advisory groups, a study of an MRSR mission began early in 1987. The study has the following goals: (1) to assess the technical feasibility of the mission; (2) to converge on two or three options for the general architecture of the mission; (3) to determine what new technologies need to be developed in order to implement the mission; (4) to define the different options sufficiently well that preliminary cost estimates can be made; and (5) to better define the science requirements. This chapter briefly describes Mars Rover/Sample Return missions that were examined in the late 1980s. These missions generally include a large (1000 kg) rover and return of over 5 kg of sample.

Carr, Michael H.↗

Exploration of Rayleigh-Taylor seeding mechanisms in laser- and pulsed-power-driven inertial fusion

The ubiquitous Rayleigh-Taylor instability (RTI) has been observed in fluids and plasmas in a wide range of parameter regimes from astrophysical to laboratory plasmas. In high-energy-density (HED) laboratory plasmas, such as laser-driven inertial confinement fusion (ICF) plasmas, the RTI can have a detrimental impact on achieving fusion ignition through the generation of hydrodynamic mix. Hence, an understanding of the seeding mechanisms that produce the RTI and identifying ways to mitigate hydrodynamic mix is of critical importance. Appropriately aligned magnetic fields have been known to stabilize short-wavelength RTI. A number of numerical and experimental studies have demonstrated the benefits of using imposed magnetic fields in laser-driven ICF to achieve higher ion temeratures, higher neutron yields, and a relaxation of the ignition criteria. This work addresses seeding mechanisms in laser-driven implosions to understand critical early-stage physics that ultimately leads to substantial growth of the RTI along with mechanisms for mitigation of this growth. Surface perturbations due to machining tolerances and single-feature seeds (for example due to fill tubes) can produce substantial RTI growth in the ignition-relevant high-convergence ratio targets for laser- driven ICF implosions. RTI growth from single-feature seeding in laser-based implosions has the potential to be mitigated through appropriately aligned externally applied magnetic fields and the goal of this work has been to quantify that numerically and experimentally. Resistive-magnetohydrodynamic (MHD) simu- lations are used to study the seeding and evolution of the RTI leveraging previous support from the DOE HEDLP program. Furthermore, experimental data has been obtained and applied towards code validation of unmagnetized and magnetized evolution of single-feature seeded RTI growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Implicit solvers for unstructured meshes

Implicit methods for unstructured mesh computations are developed and tested. The approximate system which arises from the Newton-linearization of the nonlinear evolution operator is solved by using the preconditioned generalized minimum residual technique. These different preconditioners are investigated: the incomplete LU factorization (ILU), block diagonal factorization, and the symmetric successive over-relaxation (SSOR). The preconditioners have been optimized to have good vectorization properties. The various methods are compared over a wide range of problems. Ordering of the unknowns, which affects the convergence of these sparse matrix iterative methods, is also investigated. Results are presented for inviscid and turbulent viscous calculations on single and multielement airfoil configurations using globally and adaptively generated meshes.

Venkatakrishnan, V.↗

Diagnostic budgets of analyzed and modelled tropical plumes

Blackwell et al. successfully simulated tropical plumes in a global barotropic model valid at 200 mb. The plume evolved in response to strong equatorial convergence which simulated a surge in the Walker Circulation. The defining characteristics of simulated plumes are: a subtropical jet with southerlies emanating from the deep tropics; a tropical/mid-latitude trough to the west; a convergence/divergence dipole straddling the trough; and strong cross contour flow at the tropical base of the jet. Diagnostic budgets of vorticity, divergence, and kinetic energy are calculated to explain the evolution of the modelled plumes. Budgets describe the unforced (basic) state, forced plumes, forced cases with no plumes, and ECMWF analyzed plumes.

Mcguirk, James P.↗

The Triggering Mechanism of Quiet-Region Coronal Jet Eruptions: Flux Cancelation

Coronal jets are frequent transient features on the Sun, observed in EUV and X-ray emissions. They occur in active regions, quiet Sun and coronal holes, and appear as a bright spire with base brightenings. Recent studies show that many coronal jets are driven by the eruption of a minifilament. Here we investigate the magnetic cause of jet-driving minifilament eruptions. We study ten randomly-found on-disk quiet-region coronal jets using SDO/AIA intensity images and SDO/HMI magnetograms. For all ten events, we track the evolution of photospheric magnetic flux in the jet-base region in EUV images and find that (a) a cool (transition-region temperature) minifilament is present prior to each jet eruption; (b) the pre-eruption minifilament resides above the polarity-inversion line between majority-polarity and minority-polarity magnetic flux patches; (c) the opposite-polarity flux patches converge and cancel with each other; (d) the cancelation between the majority-polarity and minority-polarity flux patches eventually destabilizes the field holding the minifilament to erupt outwards; (e) the envelope of the erupting field barges into ambient oppositely-directed far-reaching field and undergoes external reconnection (interchange reconnection); (f) the external reconnection opens the envelope field and the minifilament field inside, allowing reconnected-heated hot material and cool minifilament material to escape along the far-reaching field, producing the jet spire. In summary, we found that each of our ten jets resulted from a minifilament eruption following flux cancelation at the magnetic neutral line under the pre-eruption minifilament. These observations show that flux cancelation is usually the trigger of quiet-region coronal jet eruptions.

coronal jet↗

Saturation and Recurrence of Quantum Complexity in Random Local Quantum Dynamics

Quantum complexity is a measure of the minimal number of elementary operations required to approximately prepare a given state or unitary channel. Recently, this concept has found applications beyond quantum computing—in studying the dynamics of quantum many-body systems and the long-time properties of anti–de Sitter black holes. In this context, Brown and Susskind [] conjectured that the complexity of a chaotic quantum system grows linearly in time up to times exponential in the system size, saturating at a maximal value, and remaining maximally complex until undergoing recurrences at doubly exponential times. In this work, we prove the saturation and recurrence of complexity in two models of chaotic time evolutions based on (i) random local quantum circuits and (ii) stochastic local Hamiltonian evolution. Our results advance an understanding of the long-time behavior of chaotic quantum systems and could shed light on the physics of black-hole interiors. From a technical perspective, our results are based on establishing new quantitative connections between the Haar measure and high-degree approximate designs, as well as the fact that random quantum circuits of sufficiently high depth converge to approximate designs. Published by the American Physical Society 2024

Oszmaniec, Michał (ORCID:0000000249466835)↗

Rapid temporal evolution of radiation from non-thermal electrons in solar flares

Solutions of the time dependent Fokker-Planck equation was found for accelerated electrons undergoing Coulomb collisions in a magnetized, fully ionized plasma. An exact solution was found for arbitrary pitch angle and energy distribution in a uniform background plasma. Then, for an inhomogeneous plasma, a solution was found for particles with small pitch angles. These solutions were used to calculate the temporal evolution of bremsstrahlung x-rays from short bursts of nonthermal electron beams, and these spectra were compared with observed high time resolution spectra of short timescale solar hard x-ray bursts. It is shown that the observed softening in time of the spectra rules out a homogeneous background and therefore the possibility of electrons being confined to the corona either because of converging magnetic field or high densities. The inhomogeneous solution was also applied to a model with constant coronal density and exponentially rising chromospheric density. The spectra are shown to be consistent with that produced by a collimated beam of electrons accelerated in the corona with certain given conditions. These conditions could be violated if large pitch angle electrons are present.

Lu, Edward T.↗