Engineering PapersSearch

SEARCH · Engineering Papers

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 73 records · Page 4

Nonlinear behavior of solar gravity modes driven by He-3 in the core. II - Numerical simulations

The nonlinear behavior of gravity-mode oscillations driven by He-3-destroying reactions in the solar core has been examined by numerically integrating equations describing a very simplified model. The results of a previous bifurcation analysis, which suggest that such oscillations are unlikely to attain amplitudes sufficient to trigger core convection, are verified. These results are extended to models whose nuclear reaction rates and thermal stratification represent the core somewhat more accurately. Nonlinear processes give rise to a preference for the oscillations to develop as standing waves rather than traveling waves, thus breaking the degeneracy between these two types of motion which exists in linearized theory. Study of the large-amplitude behavior of the oscillations is hindered by a tendency for the model to become thermally unstable.

Merryfield, William J.

Decadal Polar Motion of the Earth Excited by the Convective Outer Core From Geodynamo Simulations

Long time geodetic observation records show that the orientation of the Earth's rotation axis with respect to the terrestrial reference frame, or polar motion, changes on a broad range of timescales. Apart from external torques from the luni-solar tides, these changes are excited by interactions among different components of the Earth system. The convective fluid outer core has long been conjectured a likely contributor to the observed polar motion on timescales upward of decades, such as the 30 year Markowitz wobble. We investigated the electromagnetic coupling scenario across the core-mantle boundary via numerical geodynamo simulation for different geodynamo parameters (Rayleigh numbers and magnetic Rossby numbers). Our simulated polar motion varies strongly with the dynamo parameters, while its excitation on decadal timescales appear to converge asymptotically within the adopted range of numerical Rossby numbers. Three strongest asymptotic modes emerge from numerical results, with periods around 30, 40, and 60 years for the prograde excitation and around 24, 30, and 60 years for the retrograde excitation. Their amplitudes are all larger than 5 times 10 (sup minus 8), or approximately 10 milliseconds of arc.The results suggest that the electromagnetic core-mantle coupling could explain a substantial portion, if not all, of the observed decadal polar motion. In particular, the predicted 60 year polar motion deserves special attention for future observations and studies.

Kuang, W.

Historical and Future Black Carbon Deposition on the Three Ice Caps: Ice Core Measurements and Model Simulations from 1850 to 2100

Ice core measurements in conjunction with climate model simulations are of tremendous value when examining anthropogenic and natural aerosol loads and their role in past and future climates. Refractory black carbon (BC) records from the Arctic, the Antarctic, and the Himalayas are analyzed using three transient climate simulations performed with the Goddard Institute for Space Studies ModelE. Simulations differ in aerosol schemes (bulk aerosols vs. aerosol microphysics) and ocean couplings (fully coupled vs. prescribed ocean). Regional analyses for past (1850-2005) and future (2005-2100) carbonaceous aerosol simulations focus on the Antarctic, Greenland, and the Himalayas. Measurements from locations in the Antarctic show clean conditions with no detectable trend over the past 150 years. Historical atmospheric deposition of BC and sulfur in Greenland shows strong trends and is primarily influenced by emissions from early twentieth century agricultural and domestic practices. Models fail to reproduce observations of a sharp eightfold BC increase in Greenland at the beginning of the twentieth century that could be due to the only threefold increase in the North American emission inventory. BC deposition in Greenland is about 10 times greater than in Antarctica and 10 times less than in Tibet. The Himalayas show the most complicated transport patterns, due to the complex terrain and dynamical regimes of this region. Projections of future climate based on the four CMIP5 Representative Concentration Pathways indicate further dramatic advances of pollution to the Tibetan Plateau along with decreasing BC deposition fluxes in Greenland and the Antarctic.

Antarctic Regions

Propellant flow rate through simulated liquid-core nuclear rocket fuel bed.

Experimental investigation of the validity of Zuber and Finlay's (1965) gas flow-rate formula for a two-phase flow in a rotating cylinder under high centrifugal acceleration. This formula was originally derived from tests in a 1-g environment in pipes. In the light of the investigation results obtained, the formula is valid also for a high-g environment in the rotating chamber tested.

Mcguirk, J. P.

Relevence of 3D Simulations and Sandwich Core Topology for the Modeling of Honeycomb Core Sandwich Panels Undergoing Mode I Delamination

A recently developed cohesive zone traction-separation law, which includes the effects of fiber bridging in a novel way, is extended from 2D to 3D. The proposed cohesive model is applied to low fidelity (i.e. homogenized core representation) and high fidelity (i.e. directly accounting for the core topology) finite element models of a composite panel comprised of carbon fiber reinforced plastic facesheets and a honeycomb sandwich core. This enables the investigation of 2D to 3D parameter transferability, width-dependent effects such as thumbnail-shaped crack growth, and the verification of plane strain / plane stress assumptions. A pronounced curvature of the initial interface-related crack front is observed, while the bridging-related crack front is straight. Furthermore, it is found that the cohesive parameters can easily be transferred from 2D to 3D under plane stress assumptions, but not under plane strain assumptions. The numerical predictions are compared to experimental load-displacement and R-curves.

Composites

Kernel fusion in atomistic spin dynamics simulations on Nvidia GPUs using tensor core

In atomistic spin dynamics simulations, the time cost of constructing the space- and time-displaced pair correlation function in real space increases quadratically as the number of spins N, leading to significant computational effort. The GEMM subroutine can be adopted to accelerate the calculation of the dynamical spin-spin correlation function, but the computational cost of simulating large spin systems (>40000 spins) on CPUs remains expensive. In this work, we perform the simulation on the graphics processing unit (GPU), a hardware solution widely used as an accelerator for scientific computing and deep learning. Here we show that GPUs can accelerate the simulation up to 25-fold compared to multi-core CPUs when using the GEMM subroutine on both. To hide memory latency, we fuse the element-wise operation into the GEMM kernel using CUTLASS that can improve the performance by 26% ~ 33% compared to implementation based on cuBLAS. Furthermore, we perform the on-the-fly calculation in the epilogue of the GEMM subroutine to avoid saving intermediate results on global memory, which makes the large-scale atomistic spin dynamics simulation feasible and affordable.

97 MATHEMATICS AND COMPUTING

A GPU‐Based Ocean Dynamical Core for Routine Mesoscale‐Resolving Climate Simulations

Abstract We describe an ocean hydrostatic dynamical core implemented in Oceananigans optimized for Graphical Processing Unit (GPU) architectures. On 64 A100 GPUs, equivalent to 16 computational nodes in current state‐of‐the‐art supercomputers, our dynamical core can simulate a decade of near‐global ocean dynamics per wall‐clock day at an 8‐km horizontal resolution; a resolution adequate to resolve the ocean's mesoscale eddy field. Such efficiency, achieved with relatively modest hardware resources, suggests that climate simulations on GPUs can incorporate fully eddy‐resolving ocean models. This removes a major source of systematic bias in current IPCC coupled model projections, the parameterization of ocean eddies, and represents a major advance in climate modeling. We discuss the computational strategies, focusing on GPU‐specific optimization and numerical implementation details that enable such high performance.

Silvestri, Simone [Massachusetts Institute of Tech

Advancements in NEAMS Tool Capabilities for Multiphysics Simulation of Fast Reactor Core Bowing and Identification of Validation Test Data

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves Multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Last year, thermal fluids and mechanics were coupled on a multi-assembly benchmark problem based on ABR-1000 design. This year, the reactor physics code Griffin was assessed for readiness of core bowing calculations. Preliminary integration of Griffin’s ring-heterogeneous model with thermal fluids and thermal mechanics solvers was performed. Specifically, thermal-mechanics and reactor physics were coupled for single- and multi-assembly problems, and reactor physics and subchannel methods were coupled for a single assembly model. Finally, the workflow of all three physics was preliminarily demonstrated on a single assembly model. Caveats and future development needed have been identified. To supplement the multiphysics demonstration, verification and assessment efforts of thermos-mechanical capabilities for modeling thermo-mechanical core bowing behavior were continued by analyzing IAEA Verification Problem 5 which includes radiation swelling and creep. Additionally, a small core reactor physics benchmark defined by Japan Atomic Energy Agency (JAEA) was performed to assess neutronics models for estimating reactivity feedback. Finally, Fast Flux Test Facility (FFTF) validation test data for core bowing phenomena has been identified and summarized, with a recommended path forward for validation once this capability is mature.

22 GENERAL STUDIES OF NUCLEAR REACTORS

Surface Time-Variable Gravity Signals and Possible Sources Including Core Mass Flow

Over two decades of geodetic satellite-laser-ranging (SLR) data show that the variation of the Earth's oblateness parameter J2 has a clear seasonal signal of amplitude of about 3e-10 and a secular decrease of about -2.8e-11/year, superimposed on some interesting interannual fluctuations. Physically, any change in mass distribution or/inside the Earth will be reflected in the time-variable gravity signal obtained outside the Earth, according to Newton s gravitational law. Therefore, such signal contains contributions from all geophysical sources that redistribute mass, on all temporal and spatial scales, including those from the core. Besides Earth rotation and geomagnetic field variations, the time-variable gravity also contains information linking Earth surface observations with internal core dynamical processes. The time scales of the gravity signal are critical in helping differentiate different contributions. The atmosphere and hydrosphere are responsible for the seasonal and much of the interannual and intraseasoanl fluctuations, while the secular trend is due mainly to the post-glacial rebound but possibly core mass flow. To estimate the latter effect, we use our MoSST (Modular, Scalable, Self-consistent, Three-dimensional) core dynamics model to forward simulate the core flow, and density variation due to the core convection. Our results suggest that, when upward continued to the surface, the J2 component of the core mass redistribution can reach an overall amplitude of e-11/year, approaching the SLR detectability and significant in geophysical terms. We also find a general westward drift of the mass flow, with a speed comparable to that of the geomagnetic westward drift.

Chao, Benjamin F.

Concept for coring from a low-mass rover

Future Mars missions, such as the Mars Sample Return (MSR) mission, may benefit from core sample acquisition from a low-mass rover where the rover cannot be assumed to be stationary during a coring operation. Manipulation from Mars rovers is currently done under the assumption that the rover acts as a stationary, stable platform for the arm. An MSR mission scenario with a low-mass rover has been developed and the technology needs have been investigated. Models for alternative types of coring tools and tool-environment interaction have been developed and input along with wheel-soil interaction models into the Stanford Simulation & Active Interfaces (SAI) simulation environment to enable simulation of coring operations from a rover. Coring tests using commercial coring tools indicate that the quality of the core is a critical criterion in the system design. Current results of the models, simulation, and coring tests are provided.

rovers

ECAR-7932 Rev 0 Large Eddy Simulation of MARVEL Reactor Core Subchannel to Evaluate Model Uncertainty of Reynolds-Averaged Navier-Stokes Equation Based Computational Fluid Dynamics Analysis

In the previous work (ECAR-7210), the peak cladding temperature of the MARVEL microreactor has been evaluated by steady-state Reynolds-Averaged Navier-Stokes (RANS) based computational fluid dynamics (CFD) simulations. Although numerical uncertainties of RANS-based CFD simulations has been assessed in ECAR-7210, the model uncertainty of RANS turbulence models must be investigated to resolve the issues related to inaccurate prediction of turbulent heat flux and flow pulsation in a tight lattice rod bundle using the steady-state RANS simulations. Consequently, this ECAR conducted a high-fidelity CFD analysis utilizing Large Eddy Simulation (LES) to generate reference data and investigated the model uncertainty of RANS-based CFD simulations.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN

Automatic Tip Vortex Core Profiling for Numerical Flow Simulations of Rotorcraft in Hover

An automated approach is presented that extracts visual and quantitative data from vortex cores produced by Navier-Stokes simulations of rotorcraft in hover mode. This approach extracts contiguous rotor tip vortex-core trajectories, cross-flow velocity profiles, and vortex-core diameter variation with wake age (azimuth angle). This automated approach is faster and more accurate than a conventional manual approach. Moreover, this new approach allows for an efficient way to quantitatively compare vortex-core profiles from different flow simulations, e.g., grid resolution studies, and validate computed results with experimental data

Kao, David L.