SEARCH · Engineering Papers

Results for “numerical model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,711 records · Page 95

All sources NASA NTRS DOE OSTI

NASA NTRS1976

Martian atmospheric lee waves

Areas of extensive mountain lee wave phenomenon in the northern midlatitudes during late winter are observed in Mariner 9 television pictures of Mars. The lee wave pattern is due to the discrete spectrum of waves generated by a crater ridge in an airstream. The cloud patterns resulting from the waves generated by the flow across a mountain or crater are dependent on the velocity profile of the airstream and the vertical stability of the atmosphere. A two-layer velocity model of the airstream is used in calculations based on the theory of mountain lee waves. Regardless of the simple model used for the airstream, the calculations appear to be capable of reproducing the observed wave patterns. The many observations of lee waves suggest that the near surface wind speeds as implied by the lower stream velocities do not greatly exceed, and most of the time are less than, 40 m/sec in the northern midlatitudes during late winter except possibly on a very localized scale. Results are in agreement with circulation models, especially the Leovy and Mintz (1969) two-layer numerical model.

Pirraglia, J. A.↗

NASA NTRS1984

The Use of Available Potential Energy to Evaluate Two GLAS FGGE Data Sets

Meteorological data derived from satellites provided an essential component of the FGGE data sets, especially in data sparse regions of the globe. The impact of satellite derived observations on objective analyses and numerical model predictions are investigated. The FGGE data sets and resulting numerical forecasts and the ability to define the available potential energy of the atmosphere, and its relationship to synoptic scale features such as extratropical cyclones are studied.

Horn, L. H.↗

NASA NTRS1988

Mesoscale modeling of the severe thunderstorm environment

The abilities of limited-area mesoscale models to provide accurate predictions of the environment of midlatitude severe thunderstorms and the possible feedback effects of the storms upon their environment are reviewed. Mesoaplha-scale models, mesobeta models, and terrain-induced mesoscale systems are discussed. The importance of the initial state and model numerics and physics is examined. It is found that mesoscale models must be run locally if they are to be used for short-range forecasting.

Koch, Steven E.↗

NASA NTRS1993

On the penetration of the 660 km phase change by mantle downflows

We present a simple analytic model of the interaction of cold convective downwelling currents with an endothermic phase change. The model describes the ponding and lateral spreading of downflows along the phase transition interface. A simple comparison of the vertical forces on the ponding material provides a necessary condition for a downflow to penetrate the phase boundary. This condition is fundamentally dependent on the geometry of the downflow. For planar downwellings, the model predicts a minimum ponding time before the structure can penetrate the phase boundary. For columnar (axisymmetric) downflows, there is no minimum time of spreading required before penetration can proceed. The model thus provides an explanation for the observation that in numerical models of three-dimensional convection with an endothermic phase change, cylindrical downflows penetrate the phase interface while planar ones do not. Since descending slabs in the Earth's mantle display a wide spectrum of geometries between planar and cylindrical (given various trench curvatures, as well as intersections of two or more subduction zones), this phenomenon may explain, in part, why some slabs appear to extend into the lower mantle while others are deflected at the 660 km discontinuity.

Bercovici, David↗

NASA NTRS2007

Space Plasma Ion Processing of the Lunar Soil: Modeling of Radiation-Damaged Rim Widths on Lunar Grains

Chemically and microstructurally complex altered rims around grains in the finest size fraction (<20 micron) of the lunar regolith are the result of multi-stage processes involving both solar ion radiation damage and nanoscale deposition of impact or sputter-derived vapors. The formation of the rims is an important part of the space weathering process, and is closely linked to key changes in optical reflectance and other bulk properties of the lunar surface. Recent application of field-emission scanning transmission electron microscope techniques, including energy dispersive X-ray spectral imaging, is making it easier to unravel the "nano-stratigraphy" of grain rims, and to delineate the portions of rims that represent Radiation-Amorphized (RA) host grain from overlying amorphous material that represents vapor/sputter deposits. For the portion of rims formed by host grain amorphization (henceforth called RA rims), we have been investigating the feasibility of using Monte Carlo-type ion-atom collision models, combined with experimental ion irradiation data, to derive predictive numerical models linking the width of RA rims to the grain s integrated solar ion radiation exposure time.

Chamberlin, S.↗

DOE OSTI2015

The shallow water equations as a hybrid flow model for the numerical and experimental analysis of hydro power stations

The numerical analysis of hydro power stations is an important method of the hydraulic design and is used for the development and optimisation of hydro power stations in addition to the experiments with the physical submodel of a full model in the hydraulic laboratory. For the numerical analysis, 2D and 3D models are appropriate and commonly used.The 2D models refer mainly to the shallow water equations (SWE), since for this flow model a large experience on a wide field of applications for the flow analysis of numerous problems in hydraulic engineering already exists. Often, the flow model is verified by in situ measurements. In order to consider 3D flow phenomena close to singularities like weirs, hydro power stations etc. the development of a hybrid fluid model is advantageous to improve the quality and significance of the global model. Here, an extended hybrid flow model based on the principle of the SWE is presented. The hybrid flow model directly links the numerical model with the experimental data, which may originate from physical full models, physical submodels and in-situ measurements. Hence a wide field of application of the hybrid model emerges including the improvement of numerical models and the strong coupling of numerical and experimental analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2006

A Thermodynamically Consistent Damage Model for Advanced Composites

A continuum damage model for the prediction of damage onset and structural collapse of structures manufactured in fiber-reinforced plastic laminates is proposed. The principal damage mechanisms occurring in the longitudinal and transverse directions of a ply are represented by a damage tensor that is fixed in space. Crack closure under load reversal effects are taken into account using damage variables established as a function of the sign of the components of the stress tensor. Damage activation functions based on the LaRC04 failure criteria are used to predict the different damage mechanisms occurring at the ply level. The constitutive damage model is implemented in a finite element code. The objectivity of the numerical model is assured by regularizing the dissipated energy at a material point using Bazant's Crack Band Model. To verify the accuracy of the approach, analyses of coupon specimens were performed, and the numerical predictions were compared with experimental data.

Maimi, Pere↗

NASA NTRSDate not supplied

Development of an Extended Reality (XR) Tool for Earth Science Visualization

In this presentation, we will discuss our work in adapting the NASA open source XR software, the Mixed Reality Exploration Toolkit (MRET), to an earth science domain. MRET is a NASA open source XR software for rapidly building extended reality (XR) environments for NASA domain problems, e.g., pulling in CAD models of thermal vac chamber and Roman Space Telescope to do fit checks. Primarily used for hardware integration & test, we have been adapting and extending MRET for science problems. Traditionally, scientists view and analyze the result of calculated or measured observables with static 1-D, 2-D or 3-D plots. It can be difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. Additionally, numerical models, such as the NASA GEOS climate model, are almost exclusively formulated and analyzed on Eulerian grids with points fixed in space and time. However, atmospheric phenomena such as convective clouds, hurricanes and wildfire smoke plumes move with the 3-D flow field, and it is often difficult and unnatural to understand these phenomena in an Eulerian reference frame as opposed to the Lagrangian reference frame in which nature operates. As part of an Earth Science Technology Office (ESTO) proposal, we have been adapting MRET to be a scientific exploration and analysis XR tool with integrated Lagrangian Dynamics (LD) for the Goddard Earth Observing System (GEOS) numerical weather prediction model. We believe this will help scientists identify, track, and understand the evolution of Earth Science phenomena. This presentation will discuss current results in our work in developing this XR tool for Earth Science.

Thomas G Grubb↗

NASA NTRSDate not supplied

Development of an Extended Reality (XR) Tool for Earth Science Visualization

Thomas Grubb↗

NASA NTRS2001

Geodynamo Modeling of Core-Mantle Interactions

Angular momentum exchange between the Earth's mantle and core influences the Earth's rotation on time scales of decades and longer, in particular in the length of day (LOD) which have been measured with progressively increasing accuracy for the last two centuries. There are four possible coupling mechanisms for transferring the axial angular momentum across the core-mantle boundary (CMB): viscous, magnetic, topography, and gravitational torques. Here we use our scalable, modularized, fully dynamic geodynamo model for the core to assess the importance of these torques. This numerical model, as an extension of the Kuang-Bloxham model that has successfully simulated the generation of the Earth's magnetic field, is used to obtain numerical results in various physical conditions in terms of specific parameterization consistent with the dynamical processes in the fluid outer core. The results show that depending on the electrical conductivity of the lower mantle and the amplitude of the boundary topography at CMB, both magnetic and topographic couplings can contribute significantly to the angular momentum exchange. This implies that the core-mantle interactions are far more complex than has been assumed and that there is unlikely a single dominant coupling mechanism for the observed decadal LOD variation.

Kuang, Wei-Jia↗

DOE OSTI2024

Simulation of Multiphase Flow and Poromechanical Effects Around Injection Wells in CO 2 Storage Sites

In geological CO 2 storage operations, wellbore deformations and leakage pathways formations can occur around injection and abandoned wells subjected to high rates and long-term CO 2 injection. To guide engineering design and prevent CO 2 leakage risks, a full understanding of the underlying physics and robust numerical models is necessary to evaluate the response of underground formations in the near wellbore region and in the reservoir. In this study, a multi-scale and multi-physics open-source simulator (GEOS) is used to simulate multiphase flow and poromechanical deformations over time in three dimensions. The governing equations for mechanical deformations of the rock body and multiphase compositional fluid flow within the rock matrix are solved with a fully coupled finite element and finite volume approach. The Drucker–Prager model with friction hardening is applied to simulate elastoplastic deformation and a multiphase fluid model with power-law correlations for relative permeability is used to model the migration of CO 2 plume, which are coupled with numerical implicit scheme. Simulation results are verified against multiple analytical solutions for multiphase flow and wellbore problems, thus demonstrating the accuracy of this advanced simulator. In two engineering applications, here we highlight the impact of elastoplastic deformation and coupled modeling for assessing induced displacements and stress perturbations, which are more pronounced in the near wellbore regions. This work focuses on short-term processes in the vicinity of injection wells where stress evolutions, rock deformations and multiphase compositional flow and transport are simulated jointly to ensure wellbore stability and prevent damage. This fully coupled geomechanical model can simulate multiphase flow and any associated poromechanical effects within the CO 2 storage site and in the surrounding formations. Such a large-scale, long-term, multi-physics simulation model is useful in many ways: it can guide operational decisions for CO 2 injection, assess the containment potential and risks of a site, and analyze the wellbore stability and integrity during and after CO 2 injection.

58 GEOSCIENCES↗

DOE OSTI2020

Null Space Monte Carlo Evaluation of the Plateau to River Model

The Plateau to River Groundwater Model (P2R Model) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by CH2M HILL Plateau Remediation Company at the Hanford Site in Washington State. Figure 1-1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model provides a computational framework to simulate the F&T of contaminants in groundwater associated with the 200-PO-1, 200-UP-1, 200-BP-5, and 200-ZP-1 Groundwater Operable Units (OUs) in the Hanford Site Central Plateau. In addition, the model includes adjacent areas and facilities (e.g., the State Approved Land Disposal Site). Intended and anticipated uses of the model include calculating water levels, hydraulic gradients, and groundwater flows throughout the model domain (encompassing the 200 West and 200 East Areas) for use in subsequent F&T calculations for contaminants of concern and developing scale-appropriate, telescopic-mesh refinement models for detailed evaluation of areas within the model domain where required. The overall objective of the modeling effort is to provide a basis for making informed remedial action decisions based on descriptions of current and expected future contaminant concentrations in groundwater at decision points within the OU boundaries. The objective for the model development phase is to create a common modeling platform that can be used for investigations of the Central Plateau groundwater OUs and areas downgradient toward the Columbia River. The P2R Model calibration to historical data observed at the Hanford Site is documented in CP-57037, Model Package Report for the Plateau to River Model Version 8.3. The purpose of this environmental calculation is to describe a null space Monte Carlo (NSMC) evaluation was conducted with the historic calibration of the P2R model. Use of numerical groundwater models is always accompanied with uncertainty in the results produced by a model because models are approximations of reality. Thus, by definition, lack the detail to fully represent observed behavior. Use of numerical techniques, such as a NSMC analysis, can help in identifying and quantifying the potential uncertainties associated with a numerical model such as the P2R Model. Use of the NSMC approach results in 100 groundwater flow models that are variants of the calibrated P2R Model. These variant models can be used to evaluate uncertainty in model predictions made by the calibrated P2R Model for other analyses. A secondary purpose of the environmental calculation is to establish these variant models for use with other applications.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1987

Navier-Stokes calculations and turbulence modeling in the trailing edge region of a circulation control airfoil

The accurate prediction of turbulent flows over curved surfaces in general and over the trailing edge region of circulation control airfoils in particular requires the coupled efforts of turbulence modelers, numerical analysts and experimentalists. The purpose of the research program in this area is described. Then, the influence on turbulence modeling of the flow characteristics over a typical circulation control wing is discussed. Next, the scope of this effort to study turbulence in the trailing edge region of a circulation control airfoil is presented. This is followed by a brief overview of the computation scheme, including the grid, governing equations, numerical method, boundary conditions and turbulence models applied to date. Then, examples of applications of two algebraic eddy viscosity models to the trailing edge region of a circulation control airfoil is presented. The results from the calculations is summarized, and conclusions drawn based on examples. Finally, the future directions of the program is outlined.

Viegas, John R.↗

DOE OSTI2021

Computational general relativistic force-free electrodynamics: I. Multi-coordinate implementation and testing

General relativistic force-free electrodynamics is one possible plasma-limit employed to analyze energetic outflows in which strong magnetic fields are dominant over all inertial phenomena. The amazing images of black hole (BH) shadows from the Galactic Center and the M87 galaxy provide a first direct glimpse into the physics of accretion flows in the most extreme environments of the universe. The efficient extraction of energy in the form of collimated outflows or jets from a rotating BH is directly linked to the topology of the surrounding magnetic field. We aim at providing a tool to numerically model the dynamics of such fields in magnetospheres around compact objects, such as BHs and neutron stars. To do so, we probe their role in the formation of high energy phenomena such as magnetar flares and the highly variable teraelectronvolt emission of some active galactic nuclei. In this work, we present numerical strategies capable of modeling fully dynamical force-free magnetospheres of compact astrophysical objects. Here, we provide implementation details and extensive testing of our implementation of general relativistic force-free electrodynamics in Cartesian and spherical coordinates using the infrastructure of the E INSTEIN T OOLKIT . The employed hyperbolic/parabolic cleaning of numerical errors with full general relativistic compatibility allows for fast advection of numerical errors in dynamical spacetimes. Such fast advection of divergence errors significantly improves the stability of the general relativistic force-free electrodynamics modeling of BH magnetospheres.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

An Eulerian multimaterial framework for simulating high-explosive aquarium tests

Aquarium tests of cylindrical high-explosive charges provide optical data of the detonation front velocity and shape, propagation of the shock wave in the surrounding water, and expansion rates of the detonation products behind the front. Data from aquarium experiments is often used for calibration of reactive burn models based on phenomenological equations of state (EOS) and reaction rate laws. This paper presents a multimaterial numerical modeling framework to solve the 2D axisymmetric reactive Euler equations for high-explosive aquarium tests, in particular for ammonium nitrate - fuel oil (ANFO) explosives. An extension of the Ghost Fluid Method (GFM) is used to handle the dynamic material interfaces for the ANFO explosion products, the charge-confining material (polymethyl methacrylate PMMA), and the surrounding water. This study analyzes the sensitivity of calculations (both computational efficiency and numerical accuracy) to different algorithms for the material interface models including the original GFM versus Riemann solver-based strategies. A novel method for defining the left and right states in the interfacial Riemann problem eliminates the need for sorting or nodal interpolation during the projection along the material interface. Numerical tests indicate that populating the interface node values using the Riemann solution mitigate the overheating error observed in steady-state calculations. Solution convergence and computational efficiency are explored as a function of the spatial and temporal order of the schemes. Results from the computational model with analytical equations of state and fitted reaction rate parameters show very good quantitative agreement with experimentally observed detonation front velocity, reaction products expansion, and shock wave propagation in the surrounding water for a cylindrical ANFO charge encased in PMMA. Finally, the proposed modeling framework, in conjunction with experimental tests, provides a reliable tool to assess equations of state and reaction rate expressions for reactive burn models of confined high explosives.

42 ENGINEERING↗

DOE OSTI2023

PyApprox: A software package for sensitivity analysis, Bayesian inference, optimal experimental design, and multi-fidelity uncertainty quantification and surrogate modeling

PyApprox is a Python-based one-stop-shop for probabilistic analysis of numerical models such as those used in the earth, environmental and engineering sciences. Easy to use and extendable tools are provided for constructing surrogates, sensitivity analysis, Bayesian inference, experimental design, and forward uncertainty quantification. The algorithms implemented represent a wide range of methods for model analysis developed over the past two decades, including recent advances in multi-fidelity approaches that use multiple model discretizations and/or simplified physics to significantly reduce the computational cost of various types of analyses. An extensive set of Benchmarks from the literature is also provided to facilitate the easy comparison of new or existing algorithms for a wide range of model analyses. Here, this paper introduces PyApprox and its various features, and presents results demonstrating the utility of PyApprox on a benchmark problem modeling the advection of a tracer in groundwater.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Electron-beam-ion-source (EBIS) modeling progress at FAR-TECH, Inc

FAR-TECH, Inc. has been developing a numerical modeling tool for Electron-Beam-Ion-Sources (EBISs). The tool consists of two codes. One is the Particle-Beam-Gun-Simulation (PBGUNS) code to simulate a steady state electron beam and the other is the EBIS-Particle-In-Cell (EBIS-PIC) code to simulate ion charge breeding with the electron beam. PBGUNS, a 2D (r,z) electron gun and ion source simulation code, has been extended for efficient modeling of EBISs and the work was presented previously. EBIS-PIC is a space charge self-consistent PIC code and is written to simulate charge breeding in an axisymmetric 2D (r,z) device allowing for full three-dimensional ion dynamics. This 2D code has been successfully benchmarked with Test-EBIS measurements at Brookhaven National Laboratory. For long timescale (< tens of ms) ion charge breeding, the 2D EBIS-PIC simulations take a long computational time making the simulation less practical. Most of the EBIS charge breeding, however, may be modeled in 1D (r) as the axial dependence of the ion dynamics may be ignored in the trap. Where 1D approximations are valid, simulations of charge breeding in an EBIS over long time scales become possible, using EBIS-PIC together with PBGUNS. Initial 1D results are presented. The significance of the magnetic field to ion dynamics, ion cooling effects due to collisions with neutral gas, and the role of Coulomb collisions are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Reviews and syntheses: Four decades of modeling methane cycling in terrestrial ecosystems

Abstract. Over the past 4 decades, a number of numerical models have been developed to quantify the magnitude, investigate the spatial and temporal variations, and understand the underlying mechanisms and environmental controls of methane (CH4) fluxes within terrestrial ecosystems. These CH4 models are also used for integrating multi-scale CH4 data, such as laboratory-based incubation and molecular analysis, field observational experiments, remote sensing, and aircraft-based measurements across a variety of terrestrial ecosystems. Here we summarize 40 terrestrial CH4 models to characterize their strengths and weaknesses and to suggest a roadmap for future model improvement and application. Our key findings are that (1) the focus of CH4 models has shifted from theoretical to site- and regional-level applications over the past 4 decades, (2) large discrepancies exist among models in terms of representing CH4 processes and their environmental controls, and (3) significant data–model and model–model mismatches are partially attributed to different representations of landscape characterization and inundation dynamics. Three areas for future improvements and applications of terrestrial CH4 models are that (1) CH4 models should more explicitly represent the mechanisms underlying land–atmosphere CH4 exchange, with an emphasis on improving and validating individual CH4 processes over depth and horizontal space, (2) models should be developed that are capable of simulating CH4 emissions across highly heterogeneous spatial and temporal scales, particularly hot moments and hotspots, and (3) efforts should be invested to develop model benchmarking frameworks that can easily be used for model improvement, evaluation, and integration with data from molecular to global scales. These improvements in CH4 models would be beneficial for the Earth system models and further simulation of climate–carbon cycle feedbacks.

54 ENVIRONMENTAL SCIENCES↗