Results for “numerical continuation” — Page 42

NASA NTRS1992

Magmatic intrusions and hydrothermal systems on Mars

We are continuing our investigation of Martian hydrothermal systems and the formation of fluvial valleys on Mars. Here we present our initial numerical modeling results of hydrothermal systems associated with magmatic intrusions on Mars. To model such hydrothermal systems, we consider single, cylindrical intrusions of 4 km height 2 km below the surface. Our preliminary results of modeling hydrothermal systems associated with magmatic intrusions in the Martian environment suggest that such systems, if associated with intrusions of several 10(exp 2) km(exp 3) or larger, in the presence of permeable, water-rich subsurface should be able to provide adequate discharges of water over the time periods needed to form fluvial valleys. However, it is important to note that the amount of water ultimately available for erosion also depends on the surface and subsurface lithology.

Gulick, V. C.↗

NASA NTRS1993

Interactive access and management for four-dimensional environmental data sets using McIDAS

Significant accomplishments in the past year are presented and include the following: (1) enhancements to VIS-5D; (2) Implementation of the VIS AD System; and (3) numerical modeling applications. Focus of current research and plans for next year in the following areas are briefly discussed: (1) continued development and application of the VIS-AD system; (2) further enhancements to VIS-5D; and (3) plans for modeling applications.

Hibbard, William L.↗

NASA NTRS1994

Galactic oscillations

A stable galaxy, if excited above its ground state, oscillates about that ground state. If it is resonably robust, it can support oscillations of large amplitude. Normal mode oscillations, with surprisingly large amplitudes, have been seen in numerical experiments. Observational evidence shows that real galaxies also oscillate. Galaxies ring like a bell in the experiments, and ringing continues undamped long after initial transients have died out. Their total kinetic energy oscillates with an amplitude as large as 10% of the mean. A fundamental mode dominates. It is homologous expansion/contraction of the entire galaxy (no nodes). Inward or outward velocities due to this mode are sufficiently large in the outer reaches of a galaxy to account for kinematic warps in observed velocity fields. A second spherically symmetrical mode has one node and is important near the center of the galaxy. It may be the driving force behind bulges in spiral galaxies. Two other normal modes have been identified as well. This appears to be the first experimental demonstration of normal mode oscillations within stable galaxy models.

Miller, R. H.↗

NASA NTRS2001

The CE/SE Method: a CFD Framework for the Challenges of the New Millennium

The space-time conservation element and solution element (CE/SE) method, which was originated and is continuously being developed at NASA Glenn Research Center, is a high-resolution, genuinely multidimensional and unstructured-mesh compatible numerical method for solving conservation laws. Since its inception in 1991, the CE/SE method has been used to obtain highly accurate numerical solutions for 1D, 2D and 3D flow problems involving shocks, contact discontinuities, acoustic waves, vortices, shock/acoustic waves/vortices interactions, shock/boundary layers interactions and chemical reactions. Without the aid of preconditioning or other special techniques, it has been applied to both steady and unsteady flows with speeds ranging from Mach number = 0.00288 to 10. In addition, the method has unique features that allow for (i) the use of very simple non-reflecting boundary conditions, and (ii) a unified wall boundary treatment for viscous and inviscid flows. The CE/SE method was developed with the conviction that, with a solid foundation in physics, a robust, coherent and accurate numerical framework can be built without involving overly complex mathematics. As a result, the method was constructed using a set of design principles that facilitate simplicity, robustness and accuracy. The most important among them are: (i) enforcing both local and global flux conservation in space and time, with flux evaluation at an interface being an integral part of the solution procedure and requiring no interpolation or extrapolation; (ii) unifying space and time and treating them as a single entity; and (iii) requiring that a numerical scheme be built from a nondissipative core scheme such that the numerical dissipation can be effectively controlled and, as a result, will not overwhelm the physical dissipation. Part I of the workshop will be devoted to a discussion of these principles along with a description of how the ID, 2D and 3D CE/SE schemes are constructed. In Part II, various applications of the CE/SE method, particularly those involving chemical reactions and acoustics, will be presented. The workshop will be concluded with a sketch of the future research directions.

Chang, Sin-Chung↗

NASA NTRS2003

Conducting the Heat

Heat conduction plays an important role in the efficiency and life span of electronic components. To keep electronic components running efficiently and at a proper temperature, thermal management systems transfer heat generated from the components to thermal surfaces such as heat sinks, heat pipes, radiators, or heat spreaders. Thermal surfaces absorb the heat from the electrical components and dissipate it into the environment, preventing overheating. To ensure the best contact between electrical components and thermal surfaces, thermal interface materials are applied. In addition to having high conductivity, ideal thermal interface materials should be compliant to conform to the components, increasing the surface contact. While many different types of interface materials exist for varying purposes, Energy Science Laboratories, Inc. (ESLI), of San Diego, California, proposed using carbon velvets as thermal interface materials for general aerospace and electronics applications. NASA s Johnson Space Center granted ESLI a Small Business Innovation Research (SBIR) contract to develop thermal interface materials that are lightweight and compliant, and demonstrate high thermal conductance even for nonflat surfaces. Through Phase II SBIR work, ESLI created Vel-Therm for the commercial market. Vel-Therm is a soft, carbon fiber velvet consisting of numerous high thermal conductivity carbon fibers anchored in a thin layer of adhesive. The velvets are fabricated by precision cutting continuous carbon fiber tows and electrostatically flocking the fibers into uncured adhesive, using proprietary techniques.

Source record↗

NASA NTRS1989

Rootcap structure in wild type and in a starchless mutant of Arabidopsis

Rootcaps of the wild type (WT) and of a starchless, gravitropic mutant (TC7) of Arabidopsis thaliana L. were examined by electron microscopy to identify cellular polarities with respect to gravity. In columella cells, nuclei are located proximally, and the nuclear envelope is continuous with endoplasmic reticulum (ER) that is in turn connected to nearby plasmodesmata. Impregnation of ER with osmium ferricyanide revealed numerous contacts between columella plastids and ER in both genotypes. ER is present mostly in the outer regions of the columella protoplast except in older columella cells that are developing into peripheral cells. In vertical roots, only columella cells that are intermediate in development (story 2 cells) have a higher surface density (S) of ER in the distal compared to proximal regions of the cell. The distal but not the proximal S of the ER is constant throughout columella development. Plastids are less sedimented in TC7 columella cells compared to those of the WT. It is hypothesized that plastid contact with the ER plays a role in gravity perception in both genotypes.

NASA Discipline Plant Biology↗

NASA NTRS2004

Large-Scale Alfvenic Impulses on the Sun: How They Are Generated and What We Learn From Them

NASA GSFC The Sun's atmosphere hosts a wide variety of magnetosonic disturbances. These wave modes are detected, almost exclusively, by examining images of the Sun's magnetic atmosphere and looking for propagating distortions. Although none of the Sun's plasma parameters are measured directly, we derive a great deal of information from these observations. In fact, by modeling these propagating disturbances, we may be able to derive the most accurate estimates plasma parameters. From observations absorption, refraction, reflection, and coupling of numerous wave modes, we advance our knowledge of the Sun's magnetic field, temperature, density, and current. The Sun's continuous oscillation, coronal mass ejections, flares, and other dynamic phenomena can produce wave disturbances which are observable from near-Earth space. Several of these disturbances have been traced from the inner corona out into the heliosphere. From the generation of these disturbances, we are able to learn about the phenomena which create them as well as the media through which they re-propagating. The presentation will include a discussion of the generation of Alfvenic disturbances on the Sun, ways we observe these disturbances, and how recent advances in modeling and analysis have brought us closer to determining solar in situ parameters.

Thompson, Barbara↗

NASA NTRS2010

Comet Gas and Dust Dynamics Modeling

This software models the gas and dust dynamics of comet coma (the head region of a comet) in order to support the Microwave Instrument for Rosetta Orbiter (MIRO) project. MIRO will study the evolution of the comet 67P/Churyumov-Gerasimenko's coma system. The instrument will measure surface temperature, gas-production rates and relative abundances, and velocity and excitation temperatures of each species along with their spatial temporal variability. This software will use these measurements to improve the understanding of coma dynamics. The modeling tool solves the equation of motion of a dust particle, the energy balance equation of the dust particle, the continuity equation for the dust and gas flow, and the dust and gas mixture energy equation. By solving these equations numerically, the software calculates the temperature and velocity of gas and dust as a function of time for a given initial gas and dust production rate, and a dust characteristic parameter that measures the ability of a dust particle to adjust its velocity to the local gas velocity. The software is written in a modular manner, thereby allowing the addition of more dynamics equations as needed. All of the numerical algorithms are added in-house and no third-party libraries are used.

Von Allmen, Paul A.↗

NASA NTRS2016

Trajectories for Flyby Sample Return at Saturn's Moons

Ballistic trajectories are computed which would enable a sample return mission to Titan or Enceladus without capturing, descending, or landing. The low-cost mission concept utilizes a free return trajectory that also involves a close flyby of the moon. This work extends the concept, and related trajectory analysis methodology, previously applied to a Europa mission. Specifically, a broad search algorithm is employed to systematically locate potentially feasible itineraries over an entire Saturn period. High-quality approximate solutions are then optimized to be continuous using high-fidelity dynamics. Techniques and software from the Europa analysis, were readily adapted and able to find numerous mission enabling trajectories. A direct mission to Titan is possible with flight time under 16 years and Earth-relative speeds below 11.0 km/sec. The VEEGA option is shown to substantially reduce launch C3, but flight times exceed 21 years. Unfortunately, an Enceladus mission requires a flight time of 25 years or more, and incurs fairly high relative speeds. Nevertheless, an optimized reference mission is computed.

Jones, Drew Ryan↗

NASA NTRSDate not supplied

A Fundamental Study of Bifurcation Acoustic Treatment Effects on Aft-Fan Engine Noise

Increasing air traffic and more stringent aircraft noise regulations continue to expand requirements on aircraft noise levels for conventional and unconventional aircraft configurations. A major component of the overall aircraft noise is the sound associated with the propulsion system mounted in the engine nacelle. Acoustic liners mounted in the aircraft engine nacelles provide a significant portion of the current fan noise reduction. However, they must be further optimized if challenging noise reduction goals are to be achieved. One location within the aft bypass duct that may be an excellent candidate for increased attention is the acoustic treatment on the engine bifurcations (i.e., engine pylon and lower bifurcation). This paper presents the continuation of a fundamental study of the effects of bifurcation treatment on simulated aft fan noise and the validation of numerical tools to predict such effects. Five bifurcation configurations (four treated and one hardwall) were fabricated and tested in the NASA Langley Curved Duct Test Rig. Results show that mode scattering may occur due to both the presence of the bifurcation, as well as variable impedance distributions on the bifurcation surface. Future work will also include optimization of bifurcation treatments for testing in the Curved Duct Test Rig. These initial results are promising and this work provides valuable information for further study and improvement of the performance of bifurcation acoustic treatments.

bifurcation effects↗

NASA NTRSDate not supplied

A Fundamental Study of Bifurcation Acoustic Treatment Effects on Aft-Fan Engine Noise

Increasing air traffic and more stringent aircraft noise regulations continue to expand requirements on aircraft noise levels for conventional and unconventional aircraft configurations. A major component of the overall aircraft noise is the sound associated with the propulsion system mounted in the engine nacelle. Acoustic liners mounted in the aircraft engine nacelles provide a significant portion of the current fan noise reduction. However, they must be further optimized if challenging noise reduction goals are to be achieved. One location within the aft bypass duct that may be an excellent candidate for increased attention is the acoustic treatment on the engine bifurcations (i.e., engine pylon and lower bifurcation). This paper presents the continuation of a fundamental study of the effects of bifurcation treatment on simulated aft fan noise and the validation of numerical tools to predict such effects. Five bifurcation configurations (four treated and one hardwall) were fabricated and tested in the NASA Langley Curved Duct Test Rig. Results show that mode scattering may occur due to both the presence of the bifurcation, as well as variable impedance distributions on the bifurcation surface. Future work will also include optimization of bifurcation treatments for testing in the Curved Duct Test Rig. These initial results are promising and this work provides valuable information for further study and improvement of the performance of bifurcation acoustic treatments.

bifurcation effects↗

NASA NTRSDate not supplied

A Fundamental Study of Bifurcation Acoustic Treatment Effects on Aft-Fan Engine Noise

Increasing air traffic and more stringent aircraft noise regulations continue to expand requirements on aircraft noise levels for conventional and unconventional aircraft configurations. A major component of the overall aircraft noise is the sound associated with the propulsion system mounted in the engine nacelle. Acoustic liners mounted in the aircraft engine nacelles provide a significant portion of the current fan noise reduction. However, they must be further optimized if challenging noise reduction goals are to be achieved. One location within the aft bypass duct that may be an excellent candidate for increased attention is the acoustic treatment on the engine bifurcations (i.e., engine pylon and lower bifurcation). This paper presents the continuation of a fundamental study of the effects of bifurcation treatment on simulated aft fan noise and the validation of numerical tools to predict such effects. Five bifurcation configurations (four treated and one hardwall) were fabricated and tested in the NASA Langley Curved Duct Test Rig. Results show that mode scattering may occur due to both the presence of the bifurcation, as well as variable impedance distributions on the bifurcation surface. Future work will also include optimization of bifurcation treatments for testing in the Curved Duct Test Rig. These initial results are promising and this work provides valuable information for further study and improvement of the performance of bifurcation acoustic treatments.

bifurcation effects↗

DOE OSTI2020

Calculation Of The First Moment Of Energy Using D-T Reactivity Formalisms Under The Maxwell-Boltzmann Distribution--Part II

Nuclear fusion science is an example of a scientific field with a rich history of expert involvement and scientific publications, which together, form an expert-knowledge base. One example of a nuclear fusion formalism is the utilization of published reaction rates from a variety of authors. Investigators for Deuterium- Tritium (D-T) ion fusion can choose from using frequently cited methods: the Bosch and Hal reactivity, thermonuclear reaction rates from Caughlan and Fowler, and the reactivity evaluation from Miley, Towner & Ivich which forms the basis of the Naval Research Lab (NRL) formulary. There are other choices available. Each of the reactivity formulations considered here, are based upon the Maxwell-Boltzmann velocity distribution for D-T fusion ion reactants. Numerical methods for computer codes simulating hot, energetic plasmas, include tabulations of the reactivity, and the first moment of energy. This report continues the step toward building understanding of nuclear fusion reactivity formalisms. It is part of a series of reports with the same goal, [5-10] and is the continuation of the Part I paper for defining the mathematical relationship of the first moment of D-T fusion ion kinetic energy, <$E$>, with the fusion cross-section, fusion reactivity and its derivative with ion-temperature. In Part I, three variants of the first moment <$E$> were analytically developed and explored: 1) constant cross-section, 2) a normalized first moment, and 3) a particular function of the first moment from Brysk. In Part II, attention is given to the definition of <$K$>, originally described as a ratio of moments from Brysk, and its relationship to the first moment definitions from Part I. One measure of the progress made in these documents is the identification that Brysk’s ratio of the second moment to the first moment ratio, <$K$>, does not correspond to his provided solution of the first moment of energy. Another measure of (our) progress from this work is the comparison of first moment variants. That comparison includes confirming the importance of cross sections defined in terms of energy. The analytical relationships we have developed among important physics quantities are useful tools in validation and verification (V&V). For example, we can calculate the kinetic energy as a mean or as a first moment, or as a function of the first moment. These analytically-determined values can be compared directly with numerically-determined values, supplied to the authors, representing <$E$>.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1978

LSSA large area silicon sheet task continuous Czochralski process development

A Czochralski crystal growing furnace was converted to a continuous growth facility by installation of a premelter to provide molten silicon flow into the primary crucible. The basic furnace is operational and several trial crystals were grown in the batch mode. Numerous premelter configurations were tested both in laboratory-scale equipment as well as in the actual furnace. The best arrangement tested to date is a vertical, cylindrical graphite heater containing small fused silicon test tube liner in which the incoming silicon is melted and flows into the primary crucible. Economic modeling of the continuous Czochralski process indicates that for 10 cm diameter crystal, 100 kg furnace runs of four or five crystals each are near-optimal. Costs tend to asymptote at the 100 kg level so little additional cost improvement occurs at larger runs. For these conditions, crystal cost in equivalent wafer area of around $20/sq m exclusive of polysilicon and slicing was obtained.

Rea, S. N.↗

DOE OSTI2022

An extrapolation method for strain ranges and hold times in developing the EPP+SMT creep-fatigue design curves for Alloy 617

Experimental and numerical studies in developing the integrated Elastic–Perfectly Plastic (EPP) plus Simplified Model Test (SMT) design methodology, referred to as the EPP+SMT method, continued in FY2022. This report focuses on the methods for extrapolating the EPP+SMT creep-fatigue (CF) design curves at long hold times and low strain ranges. In this study, the available CF failure data on Alloy 617 at 950°C were analyzed to determine a set of CF failure criteria. At very low strain ranges and long hold times, CF failure data are not accessible by experiments because of the extraordinarily long test durations and the inability of the test machines to accurately control these small strain ranges. A CF experimental approach with the concept of block-strain range CF testing protocol was developed. Tests using this protocol were conducted to generate the needed information for calibrating material parameters of the numerical material models. The Time Fraction based method and Dissipated Energy method were used to extrapolate the CF life curves to low strain ranges and long hold times. Based on the new experimental approach and CF life prediction methods, the CF life curves with various hold times were developed for Alloy 617 at 950°C. In addition, an experiment was designed and is being performed to verify the predicted CF curves at 950°C. The extrapolation procedure will be applied at lower temperatures to complete the development of the EPP+SMT CF design curves for Alloy 617 in F2023.

36 MATERIALS SCIENCE↗

NASA NTRS1985

Numerical solution of the vertical structure equation in the normal mode method

In the present model of multilayered stability stratification, aimed at obtaining the analytic eigensolutions of the vertical structure equation, each layer is characterized by its own static stability value. By requiring continuity of pressure and vertical velocity across each interface level, and by imposing suitable upper and lower boundary conditions, matching eigensolutions are obtained in terms of the Bessel functions. Attention is given to an explicit example of a double-layered stratified atmosphere which demonstrates the mathematical manipulations involved; the resultant vertical structure functions are used to check the accuracy of the numerical solutions by the finite difference and finite element methods.

Sasaki, Y. K.↗

DOE OSTI2025

Numerical Simulation of Electron Magnetohydrodynamics with Landau-quantized Electrons in Magnetar Crusts

Abstract In magnetar crusts, magnetic fields are sufficiently strong to confine electrons into a small to moderate number of quantized Landau levels. This can have a dramatic effect on the crust's thermodynamic properties, generating field-dependent de Haas–van Alphen oscillations. We previously argued that the large-amplitude oscillations of the magnetic susceptibility could enhance the ohmic dissipation of the magnetic field by continuously generating small-scale, rapidly dissipating field features. This could be important to magnetar field evolution and contribute to their observed higher temperatures. To study this, we performed quasi-3D numerical simulations of electron MHD in a representative volume of neutron star crust matter, for the first time including the magnetization and magnetic susceptibility resulting from Landau quantization. We find that the potential enhancement in the ohmic dissipation rate due to this effect can be a factor ∼3 for temperatures of the order of 10 8 K, and ∼4.5 for temperatures of the order of 5 × 10 7 K, depending on the magnetic field configuration. The nonlinear Hall term is crucial to this amplification: without it, the magnetic field decay is only enhanced by a factor ≲2 even at 5 × 10 7 K. These effects generate a high wavenumber plateau in the magnetic energy spectrum associated with the small-scale de Haas–van Alphen oscillations. Our results suggest that this mechanism could help explain the magnetar heating problem, though due to the effect's temperature-dependence, full magneto-thermal evolution simulations in a realistic stellar model are needed to judge whether it is viable explanation.

Rau, Peter B. (ORCID:0000000152209277)↗

NASA NTRS1981

Geometrically derived difference formulae for the numerical integration of trajectory problems

The term 'trajectory problem' is taken to include problems that can arise, for instance, in connection with contour plotting, or in the application of continuation methods, or during phase-plane analysis. Geometrical techniques are used to construct difference methods for these problems to produce in turn explicit and implicit circularly exact formulae. Based on these formulae, a predictor-corrector method is derived which, when compared with a closely related standard method, shows improved performance. It is found that this latter method produces spurious limit cycles, and this behavior is partly analyzed. Finally, a simple variable-step algorithm is constructed and tested.

Mcleod, R. J. Y.↗