Engineering Papers⌕ Search

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.

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

At least 1,801 records · Page 100
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Computer modeling of pulsed CO2 lasers for lidar applications

The experimental results will enable a comparison of the numerical code output with experimental data. This will ensure verification of the validity of the code. The measurements were made on a modified commercial CO2 laser. Results are listed as following. (1) The pulse shape and energy dependence on gas pressure were measured. (2) The intrapulse frequency chirp due to plasma and laser induced medium perturbation effects were determined. A simple numerical model showed quantitative agreement with these measurements. The pulse to pulse frequency stability was also determined. (3) The dependence was measured of the laser transverse mode stability on cavity length. A simple analysis of this dependence in terms of changes to the equivalent fresnel number and the cavity magnification was performed. (4) An analysis was made of the discharge pulse shape which enabled the low efficiency of the laser to be explained in terms of poor coupling of the electrical energy into the vibrational levels. And (5) the existing laser resonator code was changed to allow it to run on the Cray XMP under the new operating system.

Spiers, Gary D.↗
NASA NTRS2013

Assessing the Performance Limits of Internal Coronagraphs Through End-to-End Modeling

As part of the NASA ROSES Technology Demonstrations for Exoplanet Missions (TDEM) program, we conducted a numerical modeling study of three internal coronagraphs (PIAA, vector vortex, hybrid bandlimited) to understand their behaviors in realistically-aberrated systems with wavefront control (deformable mirrors). This investigation consisted of two milestones: (1) develop wavefront propagation codes appropriate for each coronagraph that are accurate to 1% or better (compared to a reference algorithm) but are also time and memory efficient, and (2) use these codes to determine the wavefront control limits of each architecture. We discuss here how the milestones were met and identify some of the behaviors particular to each coronagraph. The codes developed in this study are being made available for community use. We discuss here results for the HBLC and VVC systems, with PIAA having been discussed in a previous proceeding.

Inner working angle (IWA)↗
DOE OSTI2024

Coupled Experimental/Computational Investigation of the Dynamics of Interacting Magnetized Plasmas

The interaction, or interpenetration, of magnetized plasmas of different density and/or pressure occurs in a wide variety of natural and man-made systems. Such systems include extragalactic jets propagating into the intergalactic medium, solar coronal mass ejections into background solar wind, compact toroid (CT) fueling of magnetic fusion plasmas, and jets of capsule shell impurities into DT fusion fuel, which can lead to enhanced impurity mix in inertial fusion implosions. These plasmas may take the form of jets, with open, helical magnetic structure, or plasma “bubbles” with closed magnetic fields (B-fields), such as spheromaks or CT’s. Such structures, both open and closed B-field cases, can transport heat, particles and magnetic flux or magnetic helicity into background plasma regions. For example, the origin of extragalactic magnetic fields may be due, at least in part, to transport by astrophysical jets. The goal of this proposed work was to elucidate the detailed plasma and magnetic field dynamics of high-density plasma jets (open B-field) and bubbles (closed B-field) propagating into lower density background magnetized plasma through controlled laboratory experiments and closely coupled nonlinear MHD modeling. These experiments were conducted in the HelCat (Helicon-Cathode) linear plasma device at the University of New Mexico (UNM). Plasma jets and bubbles were launched via an existing compact coaxial plasma gun, mounted on the HelCat device. This gun produced plasmas tens of cm in scale and lasting tens of microseconds, thereby allowing detailed multipoint, space- and time-resolved measurements to be made routinely. The experiments were directly modeled using the extended magnetohydrodynamic (XMHD) PERSEUS code, developed at Cornell University [23,24]. Both experimental and numerical modeling work are ongoing. The main results to date are reported here. Additional supplemental funding for one year (8/1/2019 – 7/31/2020) supported numerical investigation of photoionization processes important in many low temperature plasmas, including the HelCat device. Initial results of this modeling work is also reported.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

Precipitation Estimates for Hydroelectricity

Hydroelectric plants require precise and timely estimates of rain, snow and other hydrometeors for operations. However, it is far from being a trivial task to measure and predict precipitation. This paper presents the linkages between precipitation science and hydroelectricity, and in doing so it provides insight into current research directions that are relevant for this renewable energy. Methods described include radars, disdrometers, satellites and numerical models. Two recent advances that have the potential of being highly beneficial for hydropower operations are featured: the Global Precipitation Measuring (GPM) mission, which represents an important leap forward in precipitation observations from space, and high performance computing (HPC) and grid technology, that allows building ensembles of numerical weather and climate models.

Tapiador, Francisco J.↗
DOE OSTI2019

The Relationship Between U.S. East Coast Sea Level and the Atlantic Meridional Overturning Circulation: A Review

Abstract Scientific and societal interest in the relationship between the Atlantic Meridional Overturning Circulation (AMOC) and U.S. East Coast sea level has intensified over the past decade, largely due to (1) projected, and potentially ongoing, enhancement of sea level rise associated with AMOC weakening and (2) the potential for observations of U.S. East Coast sea level to inform reconstructions of North Atlantic circulation and climate. These implications have inspired a wealth of model‐ and observation‐based analyses. Here, we review this research, finding consistent support in numerical models for an antiphase relationship between AMOC strength and dynamic sea level. However, simulations exhibit substantial along‐coast and intermodel differences in the amplitude of AMOC‐associated dynamic sea level variability. Observational analyses focusing on shorter (generally less than decadal) timescales show robust relationships between some components of the North Atlantic large‐scale circulation and coastal sea level variability, but the causal relationships between different observational metrics, AMOC, and sea level are often unclear. We highlight the importance of existing and future research seeking to understand relationships between AMOC and its component currents, the role of ageostrophic processes near the coast, and the interplay of local and remote forcing. Such research will help reconcile the results of different numerical simulations with each other and with observations, inform the physical origins of covariability, and reveal the sensitivity of scaling relationships to forcing, timescale, and model representation. This information will, in turn, provide a more complete characterization of uncertainty in relevant relationships, leading to more robust reconstructions and projections.

Little, Christopher M.↗
DOE OSTI2019

Analytical Model of High-Velocity Impact of a Deformable Projectile Against Textile-Based Composites

Textile-based composites are used in manufacturing of ballistic shields due to their favorable combination of mechanical properties which enables weight saving and efficiency while retaining or even improving ballistic performance. Predictive models are an important tool in reducing experimental tests which are still an important step in ballistic shield design. Among those, analytical models provide less accurate results but require far less computational cost than numerical models and can be used for preliminary analyses. Many authors in the literature have developed analytical models for high-velocity impact on textile-based composites based on wave propagation theory and an energy equilibrium approach. The main limitation of these models is the assumption of rigid projectile which leads to the underestimation of the ballistic limit for soft core projectiles, which are a very common type of threat for this kind of protection. For this reason, an enhanced version of these models is reported in this paper, which releases the assumption of rigid projectile and adds erosion and deformation of the projectile to the formulation. Ballistic tests were performed on panels manufactured from layers of Kevlar 29 plain weave fabric embedded in an epoxy matrix. An insight into the different kinds of energy absorption considered by the model shows that projectile deformation plays an important role in the dissipation of the kinetic energy of the projectile in case of intermediate impact velocity. For this reason, the present model delivers more accurate results than the models present in the literature. In particular, the estimation of the ballistic limit is drastically improved which is a key step in ballistic shields optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Bidirectional reflection of semitransparent polytetrafluoroethylene (PTFE) sheets on a silver film

In our recent work, record-high solar reflectance (exceeding 99%) has been demonstrated by placing a layer of highly scattering, semitransparent polytetrafluoroethylene (PTFE) sheet on top of a silver film. This can significantly improve the performance of daytime radiative cooling system. Although a PTFE slab of 10-mm thick acts as an excellent diffuse reflector, it is unclear how this dual-layer structure made of a thin PTFE on a mirror-like reflector will alter the bidirectional reflection. Here, we report measurements of the bidirectional reflectance distribution function (BRDF) of the dual-layer structure using a laser scatterometer at the wavelength of 635 nm. Measurements of the BRDF and the bidirectional transmittance distribution function (BTDF) of free-standing PTFE sheets are also performed. The measured BRDF for a PTFE sheet as thin as 236 μm on Ag is essentially diffuse except at large incidence angles. A Monte Carlo ray-tracing method incorporating both surface scattering and volumetric scattering is developed to theoretically model the BRDF in order to elucidate the underlying mechanisms. Finally, the numerical model captures the main features and can be used to study the bidirectional properties for PTFE sheets thinner than what is experimentally available.

42 ENGINEERING↗
DOE OSTI2024

A Variable Eddington Factor Model for Thermal Radiative Transfer with Closure Based on Data-Driven Shape Function

Here, a new variable Eddington factor (VEF) model is presented for nonlinear problems of thermal radiative transfer (TRT). The VEF model is data-driven and acts on known (a-priori) radiation-diffusion solutions for material temperatures in the TRT problem. A linear auxiliary problem is constructed for the radiative transfer equation (RTE) whose emission source and opacities are evaluated at these known material temperatures. The solution to this RTE approximates the specific intensity distribution in phase-space and time. It is applied as a shape function to define the Eddington tensor for the presented VEF model. The shape function computed via the auxiliary RTE problem will capture some degree of transport effects within the TRT problem. The VEF moment equations closed with this approximate Eddington tensor will thus carry with them these captured transport effects. In this study, the temperature data comes from multigroup P 1 , P 1/3 , and flux-limited diffusion radiative transfer models. The proposed VEF model can be interpreted as a transport-corrected diffusion reduced-order model. Numerical results are presented on the Fleck-Cummings test problem which models a supersonic wavefront of radiation. The VEF model is shown to improve accuracy by 1–2 orders of magnitude compared to the considered radiation-diffusion model solutions to the TRT problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Multiscale Electrothermal Design of a Modular Multilevel Converter for Grid-Tied Applications

As a key feature of modular multilevel converters (MMCs), a large number of semiconductor devices are employed in the converter and distributed over a stack of submodules. Each submodule has a strict temperature limit that imposes constraints on the operating range of the converter. Different loading conditions/losses in the submodules lead to unavoidable temperature variations inside the MMC, which consequently affect the system-level performance and reliability. This paper is focused on the electrothermal analysis and design of a medium-voltage silicon carbide (SiC)-based MMC system, from submodule power semiconductors to the overall MMC system integration of multiple submodules. Loss calculations are performed to estimate the cooling requirements and aid thermal design at different levels. The performance of forced-air cooling approach is analyzed within a numerical modeling framework. Maximum temperature of the SiC power modules is predicted using numerical tools. Thermal design of the MMC cabinet with various arrangement of air inlet(s) and outlet(s) is investigated and compared from a cooling performance perspective.

CFD↗
DOE OSTI2020

Advanced High-Performance Computational Modeling of the Seismic Response of High-Hazard and/or Nuclear Facilities and Critical Infrastructure at the NNSS

New methods for predicting the amplitude and variability of ground shaking from earthquakes (and explosions) are needed for seismic hazard analysis for buildings, nuclear power plants, and critical infrastructure at the NNSS. We are comparing existing 1-D and new 3-D geophysical methods for estimating the shear-wave velocity structure in the upper 30 meters of the ground surface (Vs30), which plays a major role in ground motion amplification and seismic response of buildings. We evaluate the performance of these methodologies at the U1a Complex at the NNSS and develop simple 1-D and high-resolution 3-D Vs30 models. We then emplace these high-resolution models into a background seismic velocity model. We will collaborate with Lawrence Livermore National Laboratory (LLNL) to conduct numerical modeling of the ground shaking at the NNSS using their high-performance computing technology and state-of-the-art ground motion simulation methodology. The primary work that was completed in FY 2019 was to acquire the seismic systems and familiarize staff at the NNSS with their use. We also worked on developing a collection plan with the Device Assembly Facility (DAF) at the NNSS, but due to time constraints and other ongoing projects at the DAF, we had to use U1a as a backup. We were able to coordinate the seismic survey, and we will complete the collection of seismic data in FY 2020. Additionally, during FY 2019, we completed the geologic framework model (GFM) for the U1a Complex and modeled the Yucca fault. LLNL worked with us through FY 2019 to prepare the data files for their modeling software and tested the software for reliability. In FY 2020 we will develop the end-to-end capability so that any facility could easily be modeled and the expected shaking from a local earthquake understood. The work in FY 2020 will include building fault models from the GFM and finalizing the velocity model analysis. The final simulations will be run for multiple rupture models, and final assessments will demonstrate the seismic hazard at the U1a Complex.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1985

Features of a reattaching turbulent shear layer in divergent channel flow

Experimental data have been obtained in an incompressible turbulent flow over a rearward-facing step in a diverging channel flow. Mean velocities, Reynolds stresses, and triple products that were measured by a laser Doppler velocimeter are presented for two cases of tunnel wall divergence. Eddy viscosities, production, convection, turbulent diffusion, and dissipation (balance of kinetic energy equation) terms are extracted from the data. These data are compared with various eddy-viscosity turbulence models. Numerical calculations incorporating the k-epsilon and algebraic-stress turbulence models are compared with the data. When determining quantities of engineering interest, the modified algebraic-stress model (ASM) is a significant improvement over the unmodified ASM and the unmodified k-epsilon model; however, like the others, it dramatically overpredicts the experimentally determined dissipation rate.

Driver, D. M.↗
NASA NTRS1985

Finite element algorithms for compressible flow computation on a supercomputer

Methods of applying computational fluid dynamics techniques to engineering problems are discussed. The linking factors between theoretical and applied research are the shape and weight functions which produce the Galerkin finite element schemes. The flow is unsteady, viscous, heat conducting and compressible, and steady-state flows are the asymptotic limit of unsteady flows. All flows are described by the time-averaged Navier-Stokes equations (NSE) with mass and energy conservation. Local curvilinear intrinsic coordinates are applied to discretize the NSE in arbitrary geometric domains. A modified weighted residuals approach defines the discrete analogs of the physical systems modeled and a general interpolants method is used to derive families of numerical models, both implicit and explicit, finite difference and finite elemnt, from a single point of departure. Time-dependent element approximations are achieved with a progressive assembly of generalized elements method. Hyperbolic steady-state Euler algorithms and quasi-parabolic are spatial marching algorithms for solving the discretized equations. Application of the techniques is illustrated through calculation of the compressible flow around a hypersonic flight vehicle at 100,000 ft at zero angle of attack at speeds of Mach 4-10.

Spradley, L. W.↗
NASA NTRS1994

Explosion of comet Shoemaker-Levy 9 on entry into the Jovian atmosphere

We use the astrophysical hydrocode ZEUS to compute high-resolution models of the disruption and deceleration of cometary fragments striking Jupiter. We find that simple analytic and semianalytic models work well for kilometer-size impactors. We show that previous numerical models that placed the explosion much deeper in the atmosphere failed to fully resolve important gasdynamical instabilities. These instabilities tear the comet apart, greatly increase its effective cross section, and bring it to an abrupt halt. A 1 km diameter fragment loses over 90% of its kinetic energy within a single scale height at an atmospheric pressure of order 10 bars. For all practical purposes, it explodes.

Mac Low, Mordecai-Mark↗
NASA NTRS1996

Current Results and Proposed Activities in Microgravity Fluid Dynamics

The Institute for Problems in Mechanics' Laboratory work in mathematical and physical modelling of fluid mechanics develops models, methods, and software for analysis of fluid flow, instability analysis, direct numerical modelling and semi-empirical models of turbulence, as well as experimental research and verification of these models and their applications in technological fluid dynamics, microgravity fluid mechanics, geophysics, and a number of engineering problems. This paper presents an overview of the results in microgravity fluid dynamics research during the last two years. Nonlinear problems of weakly compressible and compressible fluid flows are discussed.

Polezhaev, V. I.↗
NASA NTRS2005

Verification of a Multiphysics Toolkit against the Magnetized Target Fusion Concept

In the spring of 2004 the Advanced Concepts team at MSFC embarked on an ambitious project to develop a suite of modeling routines that would interact with one another. The tools would each numerically model a portion of any advanced propulsion system. The tools were divided by physics categories, hence the name multiphysics toolset. Currently most of the anticipated modeling tools have been created and integrated. Results are given in this paper for both a quarter nozzle with chemically reacting flow and the interaction of two plasma jets representative of a Magnetized Target Fusion device. The results have not been calibrated against real data as of yet, but this paper demonstrates the current capability of the multiphysics tool and planned future enhancements

Thomas, Scott↗
NASA NTRS2020

Lagrangian Modeling and Simulation of the Free Surface-Affected Dynamics of Underwater Vehicles

The paper deals with the development of a novelsimulator for underwater vehicles that takes into account theinteraction between the submerged body and the free watersurface. By using a Lagrangian approach, the potential flowtheory is combined with the 6-DoF equations of motions, inorder to derive a mathematical model of the system dynamics tobe solved in time-domain. A numerical model is then specificallydeveloped and an extensive simulation campaign was carried out,leading to the Submerged Bodies Simulator (SubBoS). The resultshighlight an extension of the prediction validity with respect tothe state-of-art modeling, typically based on stringent hypotheseson body motion (e.g. small displacements, pure surge motion).Index Terms—Lagrangian dynamics, modeling, simulation,potential theory, naval architecture.

Quadrelli, Marco B.↗
DOE OSTI2016

Simulating Regional Groundwater Flow in the Vicinity of the Nevada National Security Site, Nye County, Nevada - 16104

The Death Valley Regional Groundwater Flow System (DVRFS) occupies an area of about 70,000 km{sup 2} in south-central Nevada and southeastern California. Interest in this system stems from the need to (1) understand the groundwater flow paths and travel times associated with potential migration of radioactive material from the Nevada National Security Site (NNSS) where underground nuclear testing was conducted from 1956 to 1992; (2) characterize the groundwater system in the vicinity of the proposed high-level radioactive waste repository at Yucca Mountain, Nev.; and (3) address a variety of potential effects from NNSS activities on users down gradient from the NNSS and Yucca Mountain. A single groundwater flow model of the DVRFS has been developed that integrates data and results of two previous regional-scale models and supports Department of Energy (DOE) programs at the NNSS. The numerical model was constrained by organizing hydrogeologic data and interpretations in three-dimensional (3D) geographic information systems and constructing a digital three-dimensional (3D) hydrogeologic framework model (HFM) to represent 27 hydrogeologic units (HGUs) and major structures in the DVRFS region. The 3D hydrogeologic data sets were discretized to 1,500 meters (m) grid cell resolution for input arrays required for the model. The DVRFS was simulated using the USGS 3D groundwater-flow modeling code MODFLOW-2000. Two examples showing how the DVRFS regional model is used by the DOE as part of waste management activities at the NNSS: (1) use of the regional model to provide the boundary conditions for site-scale models, and (2) use of the regional model to assess changes in hydraulic gradient and flow paths in the regional flow system resulting from changes in pumping in various areas of interest. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Multi-dimensional rheology-based two-phase model for sediment transport and applications to sheet flow and pipeline scour

Sediment transport is fundamentally a two-phase phenomenon involving fluid and sediments; however, many existing numerical models are one-phase approaches, which are unable to capture the complex fluid-particle and inter-particle interactions. In the last decade, two-phase models have gained traction; however, there are still many limitations in these models. For example, several existing two-phase models are confined to one-dimensional problems; in addition, the existing two-dimensional models simulate only the region outside the sand bed. This paper develops a new three-dimensional two-phase model for simulating sediment transport in the sheet flow condition, incorporating recently published rheological characteristics of sediments. The enduring-contact, inertial, and fluid viscosity effects are considered in determining sediment pressure and stresses, enabling the model to be applicable to a wide range of particle Reynolds number. A k − ε turbulence model is adopted to compute the Reynolds stresses. In addition, a novel numerical scheme is proposed, thus avoiding numerical instability caused by high sediment concentration and allowing the sediment dynamics to be computed both within and outside the sand bed. The present model is applied to two classical problems, namely, sheet flow and scour under a pipeline with favorable results. For sheet flow, the computed velocity is consistent with measured data reported in the literature. For pipeline scour, the computed scour rate beneath the pipeline agrees with previous experimental observations. However, the present model is unable to capture vortex shedding; consequently, the sediment deposition behind the pipeline is overestimated. Sensitivity analyses reveal that model parameters associated with turbulence have strong influence on the computed results.

42 ENGINEERING↗

Refine your search to explore more results.