Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLOW MODELS”

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,387 records · Page 77

Dynamic Inland Propagation of Thinning Due to Ice Loss at the Margins of the Greenland Ice Sheet

Mass-balance analysis of the Greenland ice sheet based on surface elevation changes observed by the European Remote-sensing Satellite (ERS) (1992-2002) and Ice, Cloud and land Elevation Satellite (ICESat) (2003-07) indicates that the strongly increased mass loss at lower elevations (<2000 m) of the ice sheet, as observed during 2003-07, appears to induce interior ice thinning at higher elevations. In this paper, we perform a perturbation experiment with a three-dimensional anisotropic ice-flow model (AIF model) to investigate this upstream propagation. Observed thinning rates in the regions below 2000m elevation are used as perturbation inputs. The model runs with perturbation for 10 years show that the extensive mass loss at the ice-sheet margins does in fact cause interior thinning on short timescales (i.e. decadal). The modeled pattern of thinning over the ice sheet agrees with the observations, which implies that the strong mass loss since the early 2000s at low elevations has had a dynamic impact on the entire ice sheet. The modeling results also suggest that even if the large mass loss at the margins stopped, the interior ice sheet would continue thinning for 300 years and would take thousands of years for full dynamic recovery.

ice dynamics↗

High-Energy X-Ray Imaging of the Pulsar Wind Nebula MSH 15-52: Constraints on Particle Acceleration and Transport

We present the first images of the pulsar wind nebula (PWN) MSH 15−52 in the hard X-ray band (8 keV), as measured with the Nuclear Spectroscopic Telescope Array (NuSTAR). Overall, the morphology of the PWN as measured by NuSTAR in the 3-7 keV band is similar to that seen in Chandra high-resolution imaging. However, the spatial extent decreases with energy, which we attribute to synchrotron energy losses as the particles move away from the shock. The hard-band maps show a relative deficit of counts in the northern region toward the RCW 89 thermal remnant, with significant asymmetry. We find that the integrated PWN spectra measured with NuSTAR and Chandra suggest that there is a spectral break at 6 keV, which may be explained by a break in the synchrotron emitting electron distribution at approximately 200 TeV and/or imperfect cross calibration. We also measure spatially resolved spectra, showing that the spectrum of the PWN softens away from the central pulsar B1509−58, and that there exists a roughly sinusoidal variation of spectral hardness in the azimuthal direction. We discuss the results using particle flow models. We find non-monotonic structure in the variation with distance of spectral hardness within 50 of the pulsar moving in the jet direction, which may imply particle and magnetic-field compression by magnetic hoop stress as previously suggested for this source. We also present two-dimensional maps of spectral parameters and find an interesting shell-like structure in the N(sub H) map. We discuss possible origins of the shell-like structure and their implications.

individual objects: G320.4-1.2↗

Aerodynamic Conceptual Design of Boundary Layer Ingestion Propulsor Systems: A Quasi-2D Through Flow Analysis Method and Multi-Fidelity Propulsor Design Framework

A propulsor design framework for maximizing the benefits of boundary layer ingestion is presented. The performance of BLI is strongly affected by all propulsor components, including the boundary layer characteristics (displacement thickness and form factor) of the ingested boundary layer at the inlet, the radical loading characteristics if the fan and exit guide vane (EGV), the area contraction from inlet to nozzle, and flow expansion at the exhaust cone. A strategy and it's associated multi-fidelity design framework are proposed for an efficient conceptual design of the BLI proplusor which inherently differs from the conventional engine. In the frame work, a quasi-2D through flow model served as the underlying fidelity model is introduced to incorporate the radial effect of the boundary layer entering the propulsor. Multi-fidelity design work is conducted to maximize the predefined performance metrics. On top of this efficient quasi-2D model. computational fluid dynamics based 3-D propulsor models are implemented to refine and validate the design. A BLI propulsion system integrated with fuselage is designed to showcase the framework. The performance of the resulting BLI propulsor system is evaluated via body-force model in unstructured Reynolds-Averaged Navier-Stoke (RANS) CFD and improvement is presented.

Lee, B. J.↗

NASA SpaceCube Edge TPU SmallSat Card for Autonomous Operations and Onboard Science-Data Analysis

Using state-of-the-art artificial intelligence (AI)frameworks onboard spacecraft is challenging because common spacecraft processors cannot provide comparable performance to datacenters with server-grade CPUs and GPUs available for terrestrial applications and advanced deep-learning networks. This limitation makes small, lo w-p o we r AI microchip architectures, such as the Google Coral Edge Tensor Processing Unit (TPU), attractive for space missions where the application-specific design enables both high-performance and power-efficient computing for AI applications. To address these challenging considerations for space deployment, this research introduces the design and capabilities of a CubeSat-sized Edge TPU-based co-processor card, known as the SpaceCube Low-power Ed g e Artificial Intelligence Resilient Node (SC-LEARN). This design conforms to NASA’s CubeSat Card Specification (CS2) for integration into next-generation SmallSat and CubeSat systems. This paper describes the overarching architecture and design of the SC-LEARN, as well as, the supporting test card designed for rapid prototyping and evaluation. The SC-LEARN was developed with three operational modes: (1) a high-performance parallel-processing mode,(2)a fault-tolerant mode for onboard resilience, and (3) a power-saving mode with cold spares. Importantly, this research also elaborates on both training and quantization of Tensor Flow models for the SC-LEARN for use onboard with representative, open-source datasets. Lastly, we describe future research plans, including radiation-beam testing and flight demonstration.

Advanced avionics↗

South Repeater Building (Swmu 121) Pilot Study Work Plan, Kennedy Space Center, Florida

This document is a Pilot Study Work Plan for a groundwater remediation system at the South Repeater Building (N6-1118) located at the John F. Kennedy Space Center (KSC), Florida. Aqueous fire-fighting foam was applied near the building in the late 1990s during a wildfire. Per- and polyfluoroalkyl substances (PFAS) are components of the foam. Sampling completed at and near the South Repeater Building identified PFAS in groundwater, surface water and soil. KSC is planning to design and install a groundwater remediation system to address the potential for off-Center migration of PFAS compounds in groundwater. This Pilot Study Work Plan outlines subsurface hydrogeologic data that will be collected to use in the construction of a groundwater flow model that will guide the design of a full-scale remediation system. This document was prepared by AECOM Technical Services, Inc., for NASA under Indefinite Delivery Indefinite Quantity Contract 80KSC019D0010, Task Order 80KSC019F0096.

PFAS↗

Occurrence of Large Geomagnetically Induced Currents Within the EPRI SUNBURST Monitoring Network

Space weather, a natural hazard, can adversely impact human technological assets. High-voltage electric power transmission grids constitute one of the most critical technological systems vulnerable to space weather driven geomagnetically induced currents (GICs). One of the major challenges pertaining to the study of GICs over the continental United States has been the availability of GIC measurements, which are critical for validation of geoelectric field and power flow models, for example. In this study, we analyze GIC measurements collected at 17 Electrical Power Research Institute (EPRI) SUNBURST transformer locations across the United States for which a GIC value of 10 A or greater was recorded. This data set includes 52 individual geomagnetic storms with Kp index 6 and above during the period from 2010 to 2021. The analysis confirms that there is a good correlation between the number of geomagnetic storms per year and the number of recorded GIC events. Our results also show that about 76% of the top 17 GIC events are associated with the storm main phase, while only 24% are attributed to storm sudden commencements. In addition, it is shown, for the first time, that mid-latitude positive bays can cause large GICs over the continental United States. Finally, this study shows that the largest measured GIC event in the data set was associated with a localized intense dB/dt structure, which could be attributed to substorm activity.

Chigomezyo M. Ngwira↗

PV High Wind User Test Facility

The creation of a PV high wind test facility at the NREL Flatirons Campus that will enable private and public research efforts aimed at storm hardening PV systems. The facility will be instrumented to monitor effects of high wind conditions on modules and other array components. These in-field tests will coordinate with flow models, wind tunnel testing, validation of PV aeroelastic design codes, and post-testing module and system component analysis.

14 SOLAR ENERGY↗

The Functor system: a new on-the-fly take on Material Properties based on C++ functions

In the context of solving multiphysics problems, the discretization of the partial differential equations (PDE) at hand often takes the spotlight. However, for most engineering users and even application developers, the discretization of the equations has already been performed. Instead, they are tasked with implementing specific closure relations and material properties. MOOSE has long enabled this using the Materials system. This system relied on the pre-computation of all properties before they are used in the PDE or in postprocessing. In this talk we will introduce the Functor system, which was deployed in MOOSE in 2021, then present a few applications of functors in flow modeling simulations by the NEAMS program. Functors first offer great flexibility in their evaluation. Rather than storing various arrays for material properties, they are evaluated on the fly at the location and state, e.g. current or old value, requested. Unlike regular material properties, several operations such as the time derivative, the divergence and the curl can be requested from a functor. Similar to material properties, functors can be made to depend on arbitrary combinations of variables, functions, postprocessors and other properties. However, unlike material properties, any of these can be substituted for a functor material property. Thanks to this, objects no longer need to be duplicated based on the types of their parameters.

97 - MATHEMATICS AND COMPUTING↗

Analytical Model for Steady Flow through a Finite Channel with One Porous Wall with Arbitrary Variable Suction or Injection

This paper presents an exact solution of two-dimensional laminar flow through a finite length channel with one porous wall. It improves upon previous solutions by (1) satisfying the no-slip boundary condition at the channel dead end, (2) adding a turbulent term to the porous wall boundary condition, (3) allowing for arbitrary variable suction or injection across the porous wall, and (4) model validation against new cryogenic liquid hydrogen and oxygen experimental data. Of particular interest in the current work is the modeling of cryogenic propellant flow through a porous liquid acquisition device (LAD) screen and channel inside a propellant tank. First, a detailed review of the literature is presented for previously attempted solutions to channel flow with one porous wall. Next, the governing equations, boundary conditions, and model assumptions are used to derive the analytical flow solution and present general model results for pressure and velocity fields within the channel. Then, the model solution is compared with horizontal LAD channel flow data in liquid oxygen as well as vertical LAD channel flow data in an inverted outflow configuration in liquid hydrogen. Model results are used to update the static cryogenic bubble point pressure model with a dynamic bubble point term which factors in enhanced convection and cooling at the screen during propellant outflow. Convective heat transfer at the LAD screen during outflow is also quantified by comparing model and data. The new analytical flow solution with the dynamic bubble point model is shown to compare well with available cryogenic experimental data

Navier Stokes Equations↗

Improved melt model for power flow

Accelerators that drive z-pinch experiments transport current densities in excess of 1 MA/cm 2 in order to melt or ionize the target and implode it on axis. These high current densities stress the transmission lines upstream from the target, where rapid electrode heating causes plasma formation, melt, and possibly vaporization. These plasmas negatively impact accelerator efficiency by diverting some portion of the current away from the target, referred to as “current loss”. Simulations that are able to reproduce this behavior may be applied to improving the efficiency of existing accelerators and to designing systems operating at ever higher current densities. The relativistic particle-in-cell code CHICAGO ® is the primary code for modeling power flow on Sandia National Laboratories’ Z accelerator. We report here on new algorithms that incorporate vaporization and melt into the standard power-flow simulation framework. Taking a hybrid approach, the CHICAGO® kinetic/multi-fluid treatment has been expanded to include vaporization while the quasi-neutral equation-of-motion has been updated for melt at high current-densities. For vaporization, a new one-dimensional substrate model provides a more accurate calculation of electrode thermal, mass, and magnetic field diffusion as well as a means of emitting absorbed contaminants and vaporized metal ions. A quasi-fluid model has been implemented expressly to mimic the motion of imploding liners for accurate inductance histories. For melt, a multi-ion Hall-MHD option has been implemented and benchmarked against Alegra MHD. This new model is described with sufficient detail to reproduce these algorithms in any hybrid kinetic code. Physics results from the new code are also presented. A CHICAGO ® Hall-MHD simulation of a radial transmission line demonstrates that Hall physics, not included in Alegra, has no significant impact on the diffusion of electrode material. When surface contaminant desorption is mocked in as a hydrogen surface plasma, both the surface and bulk-material plasmas largely compress under the influence of the j × B force. Similar results are seen in Alegra, which also shows magnetic and material diffusion scaling with peak current. Test vaporization simulations using MagLIF and a power-flow experimental geometry show Fe + ions diffuse only a few hundred µm from the electrodes, so present models of Z power flow remain valid.

43 PARTICLE ACCELERATORS↗

Unified Model Deformation and Flow Transition Measurements

The number of optical techniques that may potentially be used during a given wind tunnel test is continually growing. These include parameter sensitive paints that are sensitive to temperature or pressure, several different types of off-body and on-body flow visualization techniques, optical angle-of-attack (AoA), optical measurement of model deformation, optical techniques for determining density or velocity, and spectroscopic techniques for determining various flow field parameters. Often in the past the various optical techniques were developed independently of each other, with little or no consideration for other techniques that might also be used during a given test. Recently two optical techniques have been increasingly requested for production measurements in NASA wind tunnels. These are the video photogrammetric (or videogrammetric) technique for measuring model deformation known as the video model deformation (VMD) technique, and the parameter sensitive paints for making global pressure and temperature measurements. Considerations for, and initial attempts at, simultaneous measurements with the pressure sensitive paint (PSP) and the videogrammetric techniques have been implemented. Temperature sensitive paint (TSP) has been found to be useful for boundary-layer transition detection since turbulent boundary layers convect heat at higher rates than laminar boundary layers of comparable thickness. Transition is marked by a characteristic surface temperature change wherever there is a difference between model and flow temperatures. Recently, additional capabilities have been implemented in the target-tracking videogrammetric measurement system. These capabilities have permitted practical simultaneous measurements using parameter sensitive paint and video model deformation measurements that led to the first successful unified test with TSP for transition detection in a large production wind tunnel.

Burner, Alpheus W.↗

CFD modeling of turbulent flows around the SSME main injector assembly using porosity formulation

Hot gas turbulent flow distribution around the main injector assembly of the Space Shuttle Main Engine (SSME) and liquid oxygen (LOX) flow distributions through the LOX posts have great effect on the combustion phenomenon inside the main combustion chamber. An advanced computational fluid dynamics (CFD) analysis will help to provide more accurate and efficient characterization of this type of flow field. In order to design a CFD model to be an effective engineering analysis tool with good computational turn-around time and still maintain good accuracy in describing the flow features, the concept of porosity is employed to describe the effects of blockage and drag force due to the presence of the LOX posts in the turbulent flow field around the main injector assembly of the SSME. A validated non-isotropic porosity model is developed and incorporated into an existing Navier-Stokes flow solver (FDNS). Volume and surface porosity parameters, which are based on the configurations of local LOX post clustering, are introduced into the governing equations.

Cheng, Gary C.↗

Evaluation of Wall-Modeled Les for Flow Over A Multi-Element Airfoil

We perform wall-modeled large-eddy simulations (WMLES) of turbulent flow over a 30P30N multielement airfoil at various angles of attack to assess its ability to accurately predict lift using an unstructured node-based finite-volume methodology. The flow conditions are based on the experiments reported in Klausmeyer & Lin [NASA /TM-112858]. The oncoming flow Mach number is 0.2, and the stowed chord-based Reynolds number (𝑅𝑒𝑐) is 9 million. The angle of attack (𝛼) is varied between 8◦ and 23◦ with the expected stall angle of around 21◦. The simulated geometry is periodic in the span, thus any three-dimensional effects present in the experiment are ignored. We first analyze the predictions with baseline isotropic and anisotropic grids containing over 16 grid points per boundary-layer thickness. The near-wall streamwise/nominal spacing ratio is 0.4 at the leading/trailing edges of the slat, main element and flap, and 1 elsewhere, and the first wall-normal/nominal spacing ratio is 0.375 over most of the airfoil for the anisotropic grid, while these are unity for the isotropic grid. It is found that the predicted lift and wall pressures, near-wall velocity profiles, and the location of onset of resolved turbulence are sensitive to the near-wall grid anisotropy. Further sensitivities were explored at lower angles of attack by independently varying the streamwise and wall-normal grid anisotropies. The effects of both off-wall and near-wall grid refinement were also quantified. For the grid resolutions used here, our results suggest that the first wall-normal spacing, which coincides with the WMLES exchange location, has a dominant effect on the predictions for this flow configuration.

Computational Fluid Dynamics↗

Rates of Sea‐Level Rise Are Highly Sensitive to Ice Viscosity Parameters in Model Benchmarks

Glacier flow plays a major role in current and future rates of globally averaged sea-level rise. The viscosity of glacial ice, controlling the rate of flow, decreases as stress increases and is highly sensitive to the value of the stress exponent, $n$, in the constitutive equation for viscous flow. Glaciologists and climate modelers almost exclusively assume $n=3$ when modeling ice flow and projecting sea-level rise through forward modeling. However, recent work suggests that $n\approx 4$ better fits observations, prompting the question: How sensitive are projections of sea-level rise to the value of $n$? We use an established community ice flow model and standard benchmark experiments designed as an idealized representation of Pine Island Glacier, West Antarctica. While initializing an $n=3$ model to match observations of an $n=4$ ice sheet is possible, we find that incorrectly assuming $n=3$ when in fact $n=4$ dramatically underestimates rates of sea-level rise. The scale of this error grows nonlinearly with the magnitude of the climate forcing, acting to increase projection uncertainties. Additionally, we find that models often account for this stress-dependent rheology mismatch during model initialization in a way that masks this rheological effect in the short term while leaving model outputs vulnerable to larger biases in longer-term projections. Initializations to observations of Pine Island Glacier display similar rheology-mismatch fingerprints to our idealized example.

climate sensitivity↗

Developments in Computational Modeling of Turbulent Flows

In this paper, we will discuss some recent model developments in two turbulence closure schemes. One is the Reynolds-stress algebraic equation model and the other is the Reynolds-stress transport equation model. Various model constraints required by the rapid distortion theory, the invariant theory and the realizability principle, etc. will be described in the model development.

Shih, Tsan-Hsing↗

Kinetic and Transport Modeling for Entry Flows in Hydrogen-Helium Atmospheres

Several different thermochemical models for H2/He mixtures are used to simulate a 1-D shock in a viscous CFD code at conditions relevant for hypersonic entry into the ice and gas giants. These models range from a 6-species, 1-temperature (1-T) model, to more complex 17 and 25-species state-to-state (StS) models. Corresponding radiance profiles for each of these models are then computed using a radiation code, and results are compared to experimental data from the NASA Ames Electric Arc Shock Tube (EAST) facility. Overall, the 6-species model with a quasi-steady state (QSS) solver for atomic H and modified ionization rates, along with the 17 and 25-species StS models, are found to reproduce the magnitude of the radiance seen in the EAST experiments reasonably accurately. Two additional changes to the modeling of species diffusion fluxes are also investigated: the exclusion of a ∇W term, as well as the inclusion of Soret/ Dufour effects. It is determined that these changes to the transport modeling have a negligible impact on the flowfield and therefore the radiance predictions for the H2/He shock conditions considered in this study.

Outer planets↗

Numerical Prediction of Non-Reacting and Reacting Flow in a Model Gas Turbine Combustor

The three-dimensional, viscous, turbulent, reacting and non-reacting flow characteristics of a model gas turbine combustor operating on air/methane are simulated via an unstructured and massively parallel Reynolds-Averaged Navier-Stokes (RANS) code. This serves to demonstrate the capabilities of the code for design and analysis of real combustor engines. The effects of some design features of combustors are examined. In addition, the computed results are validated against experimental data.

Davoudzadeh, Farhad↗

Numerical simulation of supersonic and hypersonic turbulent compression corner flows using relaxation models

Relaxation turbulence eddy viscosity models are incorporated to solve the complete Navier-Stokes equations for supersonic and hypersonic flows over a two-dimensional compression corner. The system of equations is solved by a time-split, second-order accurate numerical scheme. Details of relaxation process are studied, and several relaxation models are tested and compared. Good improvement in the prediction of upstream pressure propagation is obtained for a Mach number of 2.96, and Reynolds number of 10,000,000, with wedge angle of 25 deg. However, the application of the relaxation models to hypersonic flow at Mach number 8.66 and Reynolds number of 22,000,000, with highly cooled wall, shows unfavorable effects on heat transfer and skin friction in the reattachment and recompression regions.

Hung, C. M.↗