Search NASA⌕ Search

SEARCH · Search NASA

Results for “core simulation”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

An elastic-plastic analysis of fatigue crack closure in modes I and II

In this paper, an efficient elastic-plastic finite element procedure to analyze crack-closure and, its effects on fatigue crack growth under general spectrum loading, is presented. A hybrid-displacement finite element procedure is used to properly treat the stress and strain singularities near the crack-tip; and crack-growth under cyclic loading is simulated by the translation of certain 'core' elements, near the crack-tip, in which proper stress and strain singularities were embedded. Both pure mode I and mode II types of cyclic loading are considered. In the mode I case, four types of cyclic loading, viz., constant amplitude block loading, high-to-low block loading, low-to-high block, and a single overload in an otherwise constant amplitude block loading are considered; whereas in mode II, only a constant amplitude block loading is considered. Detailed results are presented for the crack-closure and opening-stresses, crack-surface deformation profiles, etc., in each case. Certain observations, based on the present numerical results, concerning various factors that cause crack growth acceleration or retardation under general spectrum loading, are presented and discussed.

Nakagaki, M.↗
NASA NTRS1984

Observations of the Earth's magnetic field from the Space Station: Measurement at high and extremely low altitude using Space Station-controlled free-flyers

Simulation studies established that the main (core), crustal and electrojet components of the Earth's magnetic field can be observed with greater resolution or over a longer time-base than is presently possible by using the capabilities provided by the space station. Two systems are studied. The first, a large lifetime, magnetic monitor would observe the main field and its time variation. The second, a remotely-piloted, magnetic probe would observe the crustal field at low altitude and the electrojet field in situ. The system design and the scientific performance of these systems is assessed. The advantages of the space station are reviewed.

Webster, W., Jr.↗
NASA NTRS1987

Incompressible Navier-Stokes computations of vortical flows over double-delta wings

An implicit flux-difference splitting scheme is used to compute incompressible laminar vortical flows about a thin round-edged double-delta wing with 80 deg and 60 deg leading-edged sweep for the strake and the main wing respectively. The numerical scheme combines approximate factorization in crossflow planes with a symmetric planar Gauss-Seidel relaxation in the remaining spatial direction. It is second-order accurate spatially by applying TVD-like upwind discretization to the inviscid fluxes and central differencing to the viscous shear fluxes. The interaction between the two primary vortices emanating from the apex and kink leading-edges is successfully simulated. The computed trajectories of vortical cores compare well with experimental data.

Hsu, Chung-Hao↗
NASA NTRS1995

Evolution of Photon and Particle Spectra in Compact, Luminous Objects

Physical conditions in the radiating plasma in the cores of radio-strong quasars and active galactic nuclei cannot be derived from observations until the effects of relativistic aberration are understood. This requires determining both the bulk flow speeds and any wave or signal speed in the parsec-scale nuclear jets. In this project we studied several aspects of such waves. We considered constraints on jet deceleration by mass pickup, and found that bolometric luminosities of the active nuclei cannot constrain core jet speeds usefully. We also simulated observations of ballistic, helical trajectories and helical waves moving directly outwards along the jet. We found that ballistic trajectories are not allowed by the data; the helical features seen are very likely to be helical waves. We believe these are waves propagating in the jet plasma. To this end, we studied waves propagating in relativistic pair plasma jets. In particular, we undertook a program whose goal was to determine the nature of waves which can propagate in relativistic pair plasmas, and how such waves propagating in streaming jet plasma would be observed by an external observer. We developed the possibility of using pulsars as test cases for our models; this takes advantage of new technology in pulsar observations, and the similarity of the physical conditions in the pulsar magnetosphere to the dense, relativistic pair plasmas which exist in radio-strong quasars.

Eilek, Jean A.↗
NASA NTRS2001

Quasar Host Galaxies/Neptune Rotation/Galaxy Building Blocks/Hubble Deep Field/Saturn Storm

Computerized animations simulate a quasar erupting in the core of a normal spiral galaxy, the collision of two interacting galaxies, and the evolution of the universe. Hubble Space Telescope (HST) images show six quasars' host galaxies (including spirals, ellipticals, and colliding galaxies) and six clumps of galaxies approximately 11 billion light years away. A false color time lapse movie of Neptune displays the planet's 16-hour rotation, and the evolution of a storm on Saturn is seen though a video of the planet's rotation. A zoom sequence starts with a ground-based image of the constellation Ursa major and ends with the Hubble Deep Field through progressively narrower and deeper views.

Source record↗
NASA NTRS2002

COR1 Engineering Test Unit Measurements at the NCAR/HAO Vacuum Tunnel Facility, October-November 2002

The Engineering Test Unit (ETU) of COR1 was made in two configurations. The first configuration, ETU-1, was for vibration testing, while the second, ETU-2, was for optical testing. This is a report on the optical testing performed on ETU-2 at the NCAR/HAO Vacuum Tunnel Facility during the months of October and November, 2002. This was the same facility used to test the two previous breadboard models. In both configurations, the first two tube sections were complete, with all optical elements aligned. The vibration model ETU-1 had the remaining tube sections attached, with mass models for the remaining optics, for the various mechanisms, and for the focal plane assembly. It was then converted into the optical model ETU-2 by removing tube sections 3 to 5, and mounting the remaining optics on commercial mounts. (The bandpass filter was also installed into tube 2, which had been replaced in ETU-1 by a mass model, so that pre- and post-vibration optical measurements could be made.) Doublet 2 was installed in a Newport LP-2 carrier, and aligned to the other optics in the first two tube sections. The LP-2 adjustment screws were then uralened so that the alignment could be maintained during shipping. Because neither the flight polarizer nor Hollow Core Motor were available, they were simulated by a commercial polarizer and rotational mount, both from Oriel corporation. The Oriel rotational stage was not designed for vacuum use, but it was determined after consultation with the company, and lab testing, that the stage could be used in the moderate vacuum conditions at the NCAR/HAO facility. The shutter and focal plane assembly were simulated with the same camera used for the previous two breadboard tests. The focal plane mask was simulated with a plane of BK7 glass with a mask glued on, using the same procedure as for the Lyot spot on Doublet 1, and mounted in an adjustable LP-2 carrier. Two masks were made, one made to the precise specifications of the optical design, the other slightly bigger to make alignment easier.

Thompson, William↗
NASA NTRS2011

Ndarts

Ndarts software provides algorithms for computing quantities associated with the dynamics of articulated, rigid-link, multibody systems. It is designed as a general-purpose dynamics library that can be used for the modeling of robotic platforms, space vehicles, molecular dynamics, and other such applications. The architecture and algorithms in Ndarts are based on the Spatial Operator Algebra (SOA) theory for computational multibody and robot dynamics developed at JPL. It uses minimal, internal coordinate models. The algorithms are low-order, recursive scatter/ gather algorithms. In comparison with the earlier Darts++ software, this version has a more general and cleaner design needed to support a larger class of computational dynamics needs. It includes a frames infrastructure, allows algorithms to operate on subgraphs of the system, and implements lazy and deferred computation for better efficiency. Dynamics modeling modules such as Ndarts are core building blocks of control and simulation software for space, robotic, mechanism, bio-molecular, and material systems modeling.

Jain, Abhinandan↗
NASA NTRS2013

Time-Accurate Local Time Stepping and High-Order Time CESE Methods for Multi-Dimensional Flows Using Unstructured Meshes

With the wide availability of affordable multiple-core parallel supercomputers, next generation numerical simulations of flow physics are being focused on unsteady computations for problems involving multiple time scales and multiple physics. These simulations require higher solution accuracy than most algorithms and computational fluid dynamics codes currently available. This paper focuses on the developmental effort for high-fidelity multi-dimensional, unstructured-mesh flow solvers using the space-time conservation element, solution element (CESE) framework. Two approaches have been investigated in this research in order to provide high-accuracy, cross-cutting numerical simulations for a variety of flow regimes: 1) time-accurate local time stepping and 2) highorder CESE method. The first approach utilizes consistent numerical formulations in the space-time flux integration to preserve temporal conservation across the cells with different marching time steps. Such approach relieves the stringent time step constraint associated with the smallest time step in the computational domain while preserving temporal accuracy for all the cells. For flows involving multiple scales, both numerical accuracy and efficiency can be significantly enhanced. The second approach extends the current CESE solver to higher-order accuracy. Unlike other existing explicit high-order methods for unstructured meshes, the CESE framework maintains a CFL condition of one for arbitrarily high-order formulations while retaining the same compact stencil as its second-order counterpart. For large-scale unsteady computations, this feature substantially enhances numerical efficiency. Numerical formulations and validations using benchmark problems are discussed in this paper along with realistic examples.

Chang, Chau-Lyan↗
DOE OSTI2024

Numerical challenges in modeling gravothermal collapse in Self-Interacting Dark Matter halos

When dark matter has a large cross section for self scattering, halos can undergo a process known as gravothermal core collapse, where the inner core rapidly increases in density and temperature. To date, several methods have been used to implement Self-Interacting Dark Matter (SIDM) in N-body codes, but there has been no systematic study of these different methods or their accuracy in the core-collapse phase. In this paper, we compare three different numerical implementations of SIDM, including the standard methods from the GIZMO and Arepo codes, by simulating idealized dwarf halos undergoing significant dark matter self interactions (σ/m = 50 cm 2 /g). When simulating these halos, we also vary the mass resolution, time-stepping criteria, and gravitational force-softening scheme. The various SIDM methods lead to distinct differences in a halo's evolution during the core-collapse phase, as each results in spurious scattering rate differences and energy gains/losses. The use of adaptive force softening for gravity can lead to numerical heating that artificially accelerates core collapse, while an insufficiently small simulation time step can cause core evolution to stall or completely reverse. Additionally, particle numbers must be large enough to ensure that the simulated halos are not sensitive to noise in the initial conditions. Even for the highest-resolution simulations tested in this study (10 6 particles per halo), we find that variations of order 10% in collapse time are still present. Here, the results of this work underscore the sensitivity of SIDM modeling on the choice of numerical implementation and motivate a careful study of how these results generalize to halos in a cosmological context.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Black Hole Supernovae, Their Equation of State Dependence, and Ejecta Composition

Abstract Recent literature on core-collapse supernovae suggests that a black hole (BH) can form within ∼1 s of shock revival, while still culminating in a successful supernova. We refer to these as BH supernovae, as they are distinct from other BH formation channels in both timescale and impact on the explosion. We simulate these events self-consistently from core collapse until 20–50 days after collapse using three axisymmetric models of a 60 M ⊙ zero-age main-sequence progenitor star and investigate how the composition of the ejecta is impacted by the BH formation. We employ Skyrme-type equations of state (EOSs) and vary the uncertain nucleonic effective mass, which affects the pressure inside the proto–neutron star through the thermal part of the EOS. This results in different BH formation times and explosion energies at BH formation, yielding final explosion energies between 0.06 and 0.72 × 10 51 erg with 21.8–23.3 M ⊙ of ejecta, of which 0–0.018 M ⊙ is 56 Ni. Compared to expectations from 1D simulations, we find more nuanced EOS dependences of the explosion dynamics, the mass of the BH remnant, and the elemental composition of the ejecta. We investigate why the explosions survive despite the massive overburden and link the shape of the diagnostic energy curve and character of the ejecta evolution to the progenitor structure.

Eggenberger Andersen, Oliver (ORCID:00000002966079↗
NASA NTRS2001

Warm-Core Intensification Through Horizontal Eddy Heat Transports into the Eye

A simulation of Hurricane Bob (1991) using the PSU/NCAR MM5 mesoscale model with a finest mesh spacing of 1.3 km is used to diagnose the heat budget of the hurricane. Heat budget terms, including latent and radiative heating, boundary layer forcing, and advection terms were output directly from the model for a 6-h period with 2-min frequency. Previous studies of warm core formation have emphasized the warming associated with gentle subsidence within the eye. The simulation of Hurricane Bob confirms subsidence warming as a major factor for eye warming, but also shows a significant contribution from horizontal advective terms. When averaged over the area of the eye, subsidence is found to strongly warm the mid-troposphere (2-9 km) while horizontal advection warms the mid to upper troposphere (5-13 km) with about equal magnitude. Partitioning of the horizontal advective terms into azimuthal mean and eddy components shows that the mean radial circulation does not, as expected, generally contribute to this warming, but that it is produced almost entirely by the horizontal eddy transport of heat into the eye. A further breakdown of the eddy components into azimuthal wave numbers 1, 2, and higher indicates that the warming is dominated by wave number 1 asymmetries, with smaller coming from higher wave numbers. Warming by horizontal eddy transport is consistent with idealized modeling of vortex Rossby waves and work is in progress to identify and clarify the role of vortex Rossby waves in warm-core intensification in both the full-physics model and idealized models.

Braun, Scott A.↗
NASA NTRS1986

The TAVERNS emulator: An Ada simulation of the space station data communications network and software development environment

The Space Station DMS (Data Management System) is the onboard component of the Space Station Information System (SSIS) that includes the computers, networks and software that support the various core and payload subsystems of the Space Station. TAVERNS (Test And Validation Environment for Remote Networked Systems) is a distributed approach for development and validation of application software for Space Station. The TAVERNS concept assumes that the different subsystems will be developed by different contractors who may be geographically separated. The TAVERNS Emulator is an Ada simulation of a TAVERNS on the ASD VAX. The software services described in the DMS Test Bed User's Manual are being emulated on the VAX together with simulations of some of the core subsystems and a simulation of the DCN. The TAVERNS Emulator will be accessible remotely from any VAX that can communicate with the ASD VAX.

Howes, Norman R.↗
NASA NTRS2018

Optimization of Elastodynamic Finite Integration Technique on Intel Xeon Phi Knights Landing Processors

This work describes the development and optimization of an implementation of an isotropic elastodynamic finite integration technique (EFIT) code for parallelized computation on Intel Knights Landing (KNL) hardware. EFIT is a numerical approach resulting in standard staggered-grid finite difference equations for the elastodynamic equations of motion to simulate bulk waves is solids. The computationally efficient simulation of elastodynamic wave propagation and interactions in aerospace materials is of high-interest in the fields of nondestructive evaluation (NDE) and structural health monitoring (SHM). Ultrasonic inspection uses an ultrasonic signal, generated at the surface of the material/structure via use of a piezoelectric transducer, to propagate sound waves into the material where it interacts with any existing defects, as well as with structural boundaries and any material inhomogeneity. Reflections from defects and boundaries are then measured by a transducer. Realistic ultrasound simulation tools can significantly aid the development and optimization of inspection techniques and can assist in the interpretation of experimental data. The optimization of an elastodynamics simulation code for the KNL Many Integrated Core processor was performed. The optimization focused on data locality and vectorization. Results show that tiling of the data to exploit the cache behavior and allow for significant utilization of the KNL hardware. The MPI implementation allows for a scalable implementation enabling large problems to be simulated. The model results were validated against theoretical dispersion curves to within 2% of the group velocity, and within 0.5% of the phase velocity of the A0 mode. Aggressive use of tiling, threading, and vectorization techniques allowed for dramatically improved time to solution.

William C. Schneck↗
DOE OSTI2024

Enabling Parallel Performance and Portability of Solid Mechanics Simulations Across CPU and GPU Architectures

Efficiently simulating solid mechanics is vital across various engineering applications. As constitutive models grow more complex and simulations scale up in size, harnessing the capabilities of modern computer architectures has become essential for achieving timely results. This paper presents advancements in running parallel simulations of solid mechanics on multi-core CPUs and GPUs using a single-code implementation. This portability is made possible by the C++ matrix and array (MATAR) library, which interfaces with the C++ Kokkos library, enabling the selection of fine-grained parallelism backends (e.g., CUDA, HIP, OpenMP, pthreads, etc.) at compile time. MATAR simplifies the transition from Fortran to C++ and Kokkos, making it easier to modernize legacy solid mechanics codes. We applied this approach to modernize a suite of constitutive models and to demonstrate substantial performance improvements across different computer architectures. This paper includes comparative performance studies using multi-core CPUs along with AMD and NVIDIA GPUs. Results are presented using a hypoelastic–plastic model, a crystal plasticity model, and the viscoplastic self-consistent generalized material model (VPSC-GMM). The results underscore the potential of using the MATAR library and modern computer architectures to accelerate solid mechanics simulations.

Morgan, Nathaniel (ORCID:0000000276118449)↗
NASA NTRS2003

Terrestrial Planet Formation in Binary Star Systems

Most stars reside in multiple star systems; however, virtually all models of planetary growth have assumed an isolated single star. Numerical simulations of the collapse of molecular cloud cores to form binary stars suggest that disks will form within such systems. Observations indirectly suggest disk material around one or both components within young binary star systems. If planets form at the right places within such circumstellar disks, they can remain in stable orbits within the binary star systems for eons. We are simulating the late stages of growth of terrestrial planets within binary star systems, using a new, ultrafast, symplectic integrator that we have developed for this purpose. We show that the late stages of terrestrial planet formation can indeed take place in a wide variety of binary systems and we have begun to delineate the range of parameter space for which this statement is true. Results of our initial simulations of planetary growth around each star in the alpha Centauri system and other 'wide' binary systems, as well as around both stars in very close binary systems, will be presented.

Lissauer, Jack J.↗
NASA NTRS2003

Terrestrial Planet Formation Around Close Binary Stars

Most stars reside in multiple star systems; however, virtually all models of planetary growth have assumed an isolated single star. Numerical simulations of the collapse of molecular cloud cores to form binary stars suggest that disks will form within such systems. Observations indirectly suggest disk material around one or both components within young binary star systems. If planets form at the right places within such circumstellar disks, they can remain in stable orbits within the binary star systems for eons. We are simulating the late stages of growth of terrestrial planets around close binary stars, using a new, ultrafast, symplectic integrator that we have developed for this purpose. The sum of the masses of the two stars is one solar mass, and the initial disk of planetary embryos is the same as that used for simulating the late stages of terrestrial planet growth within our Solar System and in the Alpha Centauri wide binary star system. Giant planets &are included in the simulations, as they are in most simulations of the late stages of terrestrial planet accumulation in our Solar System. When the stars travel on a circular orbit with semimajor axis of up to 0.1 AU about their mutual center of mass, the planetary embryos grow into a system of terrestrial planets that is statistically identical to those formed about single stars, but a larger semimajor axis and/or a significantly eccentric binary orbit can lead to significantly more dynamically hot terrestrial planet systems.

Lissauer, Jack J.↗
NASA NTRS2017

Initial Data Analysis Results for ATD-2 ISAS HITL Simulation

To evaluate the operational procedures and information requirements for the core functional capabilities of the ATD-2 project, such as tactical surface metering tool, APREQ-CFR procedure, and data element exchanges between ramp and tower, human-in-the-loop (HITL) simulations were performed in March, 2017. This presentation shows the initial data analysis results from the HITL simulations. With respect to the different runway configurations and metering values in tactical surface scheduler, various airport performance metrics were analyzed and compared. These metrics include gate holding time, taxi-out in time, runway throughput, queue size and wait time in queue, and TMI flight compliance. In addition to the metering value, other factors affecting the airport performance in the HITL simulation, including run duration, runway changes, and TMI constraints, are also discussed.

HITL simulation↗
NASA NTRS2002

Numerical Simulations of Vortex Generator Vanes and Jets on a Flat Plate

Numerical simulations of a single low-profile vortex generator vane, which is only a small fraction of the boundary-layer thickness, and a vortex generating jet have been performed for flows over a flat plate. The numerical simulations were computed by solving the steady-state solution to the Reynolds-averaged Navier-Stokes equations. The vortex generating vane results were evaluated by comparing the strength and trajectory of the streamwise vortex to experimental particle image velocimetry measurements. From the numerical simulations of the vane case, it was observed that the Shear-Stress Transport (SST) turbulence model resulted in a better prediction of the streamwise peak vorticity and trajectory when compared to the Spalart-Allmaras (SA) turbulence model. It is shown in this investigation that the estimation of the turbulent eddy viscosity near the vortex core, for both the vane and jet simulations, was higher for the SA model when compared to the SST model. Even though the numerical simulations of the vortex generating vane were able to predict the trajectory of the stream-wise vortex, the initial magnitude and decay of the peak streamwise vorticity were significantly under predicted. A comparison of the positive circulation associated with the streamwise vortex showed that while the numerical simulations produced a more diffused vortex, the vortex strength compared very well to the experimental observations. A grid resolution study for the vortex generating vane was also performed showing that the diffusion of the vortex was not a result of insufficient grid resolution. Comparisons were also made between a fully modeled trapezoidal vane with finite thickness to a simply modeled rectangular thin vane. The comparisons showed that the simply modeled rectangular vane produced a streamwise vortex which had a strength and trajectory very similar to the fully modeled trapezoidal vane.

Allan, Brian G.↗