Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical modelling”

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,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Turbulence, Transition, and Numerical Method Technologies

The close collaboration in the validation experimentation effort was excellent, the complementary research efforts related tofundamental numerical model development, and the activities related to effective HPC utilization on near term architectures that are coming down the pipeline were exactly what NASA should be doing. Even though the TQR panel cannot say that the success criteria were fully met, we do feel that the research funded by TTT/RCA was successful and useful in determining the next suite of CFD validation tests and the problematic flow conditions that continue to be difficult to predict using our current suite of computational tools. Given the TQR Panel’s review criteria: 1. The deliverables meet the technical requirements; 2.The appropriate technical approaches were followed in producing the deliverables; and 3. All technical caveats and concerns have been identified and the deliverables: a.Development of more accurate physics-based methods (e.g., higher moment closure); b. Large eddy simulation (LES); c. Advanced numerical methods; d. Transition prediction and modeling; e.Validation experiments; f. Multidisciplinary analysis and design (high fidelity). TQR panel feels that the technical processes were all completed and that the review criteria and all the deliverables were met. Thus, the Milestone TACP-2016-001 has been completed.

Mujeeb R Malik↗
DOE OSTI2022

Enhancement of phase change material hysteresis model: A case study of modeling building envelope in EnergyPlus

Nowadays, buildings are expected to offer demand side services to the power grid to enhance the electrical load flexibility, which leads to the concepts of grid-interactive efficient buildings (GEBs). Phase change material (PCM)-based thermal energy storage has seen increasing attention in recent years for peak load shifting of grid-interactive efficient buildings (GEBs). Numerical models are critical tools for design and evaluation of PCM-integrated systems. Most industrial-grade PCMs are reported to melt/freeze over a temperature range instead of at a unique temperature. Such thermal hysteresis effect significantly affects the reliability of simulation results because not only the heat transfer process depends on melting and freezing temperatures, the PCM thermal properties change significantly during the phase change process as well. This study is aimed to develop a model for the PCMs used in the building envelope with the capability to accurately simulate hysteretic behaviors. Further, this model is based on a two-phase assumption and is implemented in a whole building energy performance simulation program (i.e., EnergyPlus). A comparison between numerical results and experimental data shows that during a complete phase transition, the two-phase model could achieve a good agreement with the experimental data. During a partial phase transition, the two-phase model could lead to significant improvements compared to other alternative PCM models, including the existing PCM model in EnergyPlus. Last, whole building simulations were performed to study this model's performance regarding heating/cooling loads and zone mean air temperature of a given building. The results show that the difference in hourly heating/cooling loads introduced by the models was less than 1% in design conditions, while significant changes were observed in both hourly heating/cooling loads and zone mean air temperature when the PCM envelope underwent partial phase transition processes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

On the Feasibility of Deep Geothermal Wells Using Numerical Reservoir Simulation

This study examines the geothermal energy extraction potential from the basement rock within the Denver–Julesburg Basin, focusing on the flow performance and heat extraction efficiency of different geothermal well configurations. It specifically compares U-shaped, V-shaped, inclined V-shaped, and pipe-in-pipe configurations against enhanced geothermal system setups. Through numerical modeling, we evaluated the thermal behavior of these systems under various operational scenarios and fracture conditions. The results suggest that while closed-loop systems offer moderate temperature increases, Enhanced geothermal system configurations show substantial potential for high-temperature extraction. This underscores the importance of evaluating well configurations in complex geological settings. The insights from this study aid in strategic geothermal energy planning and development, marking significant advancements in geothermal technology and setting a foundation for future explorations and optimizations.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Phase-field modeling of crack growth under coupled creep-fatigue

Crack growth under coupled creep-fatigue strongly influences service life of metallic components operating at high temperatures. However, there lack numerical models that enable direct simulations of crack growth and plasticity development under coupled creep-fatigue. Here, in this study, a unified phase-field model is developed for simulating crack growth under fatigue, creep, and coupled creep-fatigue. This model is able to reproduce the Paris laws for both brittle and ductile materials in pure fatigue mode, which shows that the degradation rate due to cyclic stress in ductile materials is generally larger than that in brittle materials. Fatigue simulations with different yield strengths demonstrate that plasticity facilitates crack growth in typical ductile materials. High-throughput simulations are performed under coupled creep-fatigue conditions, with varying hold-stress levels and hold-time. By analyzing crack growth rate under different hold-time, a creep-fatigue interaction term is obtained. Simulation results demonstrate that the creep-fatigue interaction is caused by smaller stress gradient and enhanced degradation rate ahead of the crack tips under coupled creep-fatigue, compared to the situation of pure fatigue. This work reveals the origin of different Paris-law exponents in brittle and ductile materials and demonstrates how the interplay between creep and fatigue affects crack growth.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Electric Water Heaters for Transactive Systems: Model Evaluations and Performance Quantification

Electric water heaters (EWHs) are opportune appliances for implementing demand-side control. EWH models serve as a fundamental step toward accurately estimating EWH flexibility potential and designing proper control strategies. Existing studies have adapted numerous modeling approaches in evaluating the potential of EWHs for a variety of grid applications. This paper presents an analytical study that evaluates the performance of state-of-art EWH models in terms of accuracy and computational complexity for adaptation in evaluation studies for the transactive system. The work proposes a transactive control strategy that optimally utilizes the thermal inertia of EWHs for providing grid services. Here, the performance of the control strategy and the impact of modeling accuracy is evaluated for device-level and feeder-level use cases using the IEEE 123 node test distribution system appropriately populated with EWHs. The simulation results illustrate the effectiveness of the control strategy in reducing the feeder demand during peak period by 13% and also quantities the impact of using simplified modeling approaches for determining the potential of EWHs for providing grid services.

42 ENGINEERING↗
DOE OSTI2021

Development of a Random-Forest Cloud-Regime Classification Model Based on Surface Radiation and Cloud Products

Various methods have been developed to characterize cloud type, otherwise referred to as cloud regime. These include manual sky observations, combining radiative and cloud vertical properties observed from satellite, surface-based remote sensing, and digital processing of sky imagers. While each method has inherent advantages and disadvantages, none of these cloud-typing methods actually includes measurements of surface shortwave or longwave radiative fluxes. Here, a method that relies upon detailed, surface-based radiation and cloud measurements and derived data products to train a random-forest machine-learning cloud classification model is introduced. Measurements from five years of data from the ARM Southern Great Plains site were compiled to train and independently evaluate the model classification performance. A cloud-type accuracy of approximately 80% using the random-forest classifier reveals that the model is well suited to predict climatological cloud properties. Furthermore, an analysis of the cloud-type misclassifications is performed. While physical cloud types may be misreported, the shortwave radiative signatures are similar between misclassified cloud types. From this, we assert that the cloud-regime model has the capacity to successfully differentiate clouds with comparable cloud–radiative interactions. Therefore, we conclude that the model can provide useful cloud-property information for fundamental cloud studies, inform renewable energy studies, and be a tool for numerical model evaluation and parameterization improvement, among many other applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Vadose Zone Flow and Transport Parameters Data Package for the Hanford Site Composite Analysis

This report provides a description of the basis for the development and implementation of a conceptual model for vadose zone flow and transport for the Hanford Site composite analysis (CA) groundwater pathway analysis. The parameterization for a numerical model is intimately linked to the conceptual model framework. This report describes the basis for the selection of hydraulic and transport parameters for the hydrostratigraphic units (HSUs) identified in the 200 East and 200 West Areas. The following three-step process was used to develop the hydraulic properties and transport parameters for the CA modeling: 1. Breakdown of existing database on hydraulic properties by HSUs 2. Constitutive model parameters for laboratory, core-scale hydraulic properties 3. Upscaling for macroscopic, field-scale flow, and transport parameters This document provides data and information on the following items: A brief description and technical basis for the conceptual model selected to represent vadose zone flow and transport for the CA modeling. Variation in sediment particle size distribution for representative sites in the 200 East and West Areas. A summary of effective (upscaled) moisture retention, and saturated and unsaturated hydraulic conductivity estimates for the HSUs in the 200 East and West Areas. A summary of effective transport parameters, including sediment bulk density and macrodispersivity estimates for various HSUs. The 200 East Area B Complex perched water aquifer hydraulic properties

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

Magnetically Coupled Broadband Dual Magnetic Mass/Spring Vibrational Energy Harvesting Design

Self-powering sensors and networks are a reality. The ability to extract ambient energy from the surroundings to power electronic devices has a profound impact on the realization of smart adaptable sensor networks. In this study, a magnetically coupled dual spring and magnet design has been investigated to improve the efficiency and performance bandwidth of vibration energy harvesting (VEH) sensors. Using numerical models based on traditional systems of coupled ordinary differential equations (ODE), an optimized design was developed and compared to experimental measurements. Numerical and empirical results show good agreement. Results show improvement in the bandwidth over an equivalent linear system and corresponding improvement in output power conversion efficiency. The increased bandwidth allows improved conversion sensitivity and enhanced power harvesting capabilities. This operational bandwidth coincides with the expected input spectrum for in situ applications.

42 ENGINEERING↗
NASA NTRS1981

On the amplitudes reached by baroclinically unstable disturbances

The maximum eddy energies reached by baroclinically unstable disturbances in some simplified numerical models (Gall, 1976; Simmons and Hoskins, 1978) are compared to a scale-dependent measure of the energy available to the disturbances. It is found that, except for the longest wavelengths, the ratio of the maximum eddy energies to this available energy (the eddy efficiency) tends to remain constant relative to the large variations in the maximum energies. The roles of the growth rate of the initial disturbance and the beta-effect are discussed with reference to results from a maximally simplified model. It is shown that dimensional arguments paralleling Charney (1971) lead to the maximum energy versus scale relationships found in the numerical results.

Schneider, E. K.↗
NASA NTRS1983

The construction of life prediction models for the design of Stirling engine heater components

The service life of Stirling-engine heater structures of Fe-based high-temperature alloys is predicted using a numerical model based on a linear-damage approach and published test data (engine test data for a Co-based alloy and tensile-test results for both the Co-based and the Fe-based alloys). The operating principle of the automotive Stirling engine is reviewed; the economic and technical factors affecting the choice of heater material are surveyed; the test results are summarized in tables and graphs; the engine environment and automotive duty cycle are characterized; and the modeling procedure is explained. It is found that the statistical scatter of the fatigue properties of the heater components needs to be reduced (by decreasing the porosity of the cast material or employing wrought material in fatigue-prone locations) before the accuracy of life predictions can be improved.

Petrovich, A.↗
NASA NTRS1989

The evolution of complex type B Allende inclusion - An ion microprobe trace element study

Results are presented of a detailed trace-element and isotopic analyses of the constituent phases in each of the major textural parts (mantle, core, and islands) of a Type B refractory inclusion, the USNM 5241 inclusion from Allende, first described by El Goresy et al. (1985). The REE data on 5241 were found to be largely consistent with a model in which the mantle and the core of 5241 formed sequentially out of a single melt by fractional crystallization. The numerical models of REE evolution in the 5241 melt, especially that of Eu, require that a significant mass of spinel-free island material was assimilated into the evolving melt during the last half of the solidification history of 5241. The trace element results pbtained thus strongly support the interpretation of El Goresy et al. (1985) that the spinel-free islands in the 5241 are trapped xenoliths.

Macpherson, Glenn J.↗
NASA NTRS1992

Model of an axially strained weakly guiding optical fiber modal pattern

Axial strain can be determined by monitoring the modal pattern variation of an optical fiber. The results of a numerical model developed to calculate the modal pattern variation at the end of a weakly guiding optical fiber under axial strain is presented. Whenever an optical fiber is under stress, the optical path length, the index of refraction, and the propagation constants of each fiber mode change. In consequence, the modal phase term for the fields and the fiber output pattern are also modified. For multimode fibers, very complicated patterns result. The predicted patterns are presented, and an expression for the phase variation with strain is derived.

Egalon, Claudio O.↗
NASA NTRS1998

Effects of Shear, Defocus, and Wavefront Error on the Theoretical Performance of the Composite Infrared Spectrometer for Cassini

The combined effects on performance of shear between the two arms, defocus of die detector, and difference in wavefront between the two arms of a Fourier transform spectrometer using cube corner retroreflectors were investigated. Performance was characterized by the amplitude of the fringe signals coming from a detector as the path-length difference was scanned. A closed-form expression was found for the combined effects of shear and defocus, and it was found that defocus had no effect in the absence of shear. The effect of wavefront error was modeled numerically and assumed to be independent of shear and defocus. Results were compared with measurements made on the breadboard and engineering model of the Composite Infrared Spectrometer for the Cassini mission to Saturn, and good agreement was found.

Martino, Antonio J.↗
NASA NTRS2006

The Effect of the Saharan Air Layer on the Formation of Hurricane Isabel (2003) Simulated with AIRS Data

The crucial physics of how the atmosphere really accomplishes the tropical cyclogenesis process is still poorly understood. The presence of the Saharan Air Layer (SAL), an elevated mixed layer of warm and dry air that extends from Africa to the tropical Atlantic and contains a substantial amount of mineral dust, adds more complexity to the tropical cyclogenesis process in the Atlantic basin. The impact of the SAL on tropical cyclogenesis is still uncertain. Karyampudi and Carlson (1988) conclude that a strong SAL can potentially aid tropical cyclone development while Dunion and Velden (2004) argue that the SAL generally inhibits tropical cyclogenesis and intensification. Advancing our understanding of the physical mechanisms of tropical cyclogenesis and the associated roles of the SAL strongly depends on the improvement in the observations over the data-sparse ocean areas. After the Atmospheric Infrared Sounder (AIRS), the Advanced Microwave Sounding Unit (AMSU), and the microwave Humidity Sounder of Brazil (HSB) were launched with the NASA Aqua satellite in 2002, new data products retrieved from the AIRS suite became available for studying the effect of the warm, dry air mass associated with the SAL (referred to as the thermodynamic effect). The vertical profiles of the AIRS retrieved temperature and humidity provide an unprecedented opportunity to examine the thermodynamic effect of the SAL. The observational data can be analyzed and assimilated into numerical models, in which the model thermodynamic state is allowed to relax to the observed state from AIRS data. The objective of this study is to numerically demonstrate that the thermodynamic effect of the SAL on the formation of Hurricane Isabel (2003) can be largely simulated through nudging of the AIRS data.

Wu, iguang↗
NASA NTRS2009

A Model of Thermal Conductivity for Planetary Soils: 1. Theory for Unconsolidated Soils

We present a model of heat conduction for mono-sized spherical particulate media under stagnant gases based on the kinetic theory of gases, numerical modeling of Fourier s law of heat conduction, theoretical constraints on the gas thermal conductivity at various Knudsen regimes, and laboratory measurements. Incorporating the effect of the temperature allows for the derivation of the pore-filling gas conductivity and bulk thermal conductivity of samples using additional parameters (pressure, gas composition, grain size, and porosity). The radiative and solid-to-solid conductivities are also accounted for. Our thermal model reproduces the well-established bulk thermal conductivity dependency of a sample with the grain size and pressure and also confirms laboratory measurements finding that higher porosities generally lead to lower conductivities. It predicts the existence of the plateau conductivity at high pressure, where the bulk conductivity does not depend on the grain size. The good agreement between the model predictions and published laboratory measurements under a variety of pressures, temperatures, gas compositions, and grain sizes provides additional confidence in our results. On Venus, Earth, and Titan, the pressure and temperature combinations are too high to observe a soil thermal conductivity dependency on the grain size, but each planet has a unique thermal inertia due to their different surface temperatures. On Mars, the temperature and pressure combination is ideal to observe the soil thermal conductivity dependency on the average grain size. Thermal conductivity models that do not take the temperature and the pore-filling gas composition into account may yield significant errors.

Piqueux, S.↗
DOE OSTI2020

Analysis of the performance of different culvert boundary conditions in 2D shallow flow models

Culverts allow roads to safely traverse small streams or drainage ditches, and their proper design is critical to ensure a safe and reliable transportation network. A correct modelization of these hydraulic structures becomes crucial in the assessment of flood footprints or discharge peak estimation in a risk evaluation plan. The question of how to include culverts comes up frequently when assembling a hydraulic model that requires the presence of as many singular elements as possible. In this work, three different culvert integrations with the surface domain are studied and compared in the context of a 2D shallow water (SW) model. All of them are based on the Federal Highway Administration (FHWA) formulation for the culvert discharge estimation but differ in complexity and in the interaction with the numerical model for surface flow, some of them as internal boundary conditions. Several steady and unsteady validation test cases are presented and the numerical results are compared with the predictions from HEC-RAS 1D and HY-8 software. Furthermore, the culvert area, shape and their sensitivity to the 2D computational mesh is also analyzed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Multi-dimensional modeling of mixture preparation in a direct injection engine fueled with gaseous hydrogen

With the recent advances of direct injection (DI) technology, introducing hydrogen into the combustion chamber through DI is being considered as a viable approach to circumvent backfire and pre-ignition encountered in early generations of hydrogen engines. As part of a broader vision to develop a robust numerical model to study hydrogen spark ignition (SI) combustion in internal combustion (IC) engines, the present numerical investigation focuses on mixture preparation in a hydrogen DI SI engine. This study is carried out with a single hole injector with gaseous hydrogen injected at 100 bar injection pressure. Simulations are carried out for high and low tumble configurations and validated against optical data acquired from planar laser induced fluorescence (PLIF) measurements. Varying mesh configurations are investigated for the impact on in-cylinder mixture distribution. A particular emphasis is placed on the effect of nozzle geometry and mesh orientation near the wall. Overall, the computational model is found to predict the mixture distribution in the combustion cylinder reasonably well. The results showed that the alignment of mesh with the flow direction is important to achieve good agreement between numerical analysis and optical measurement data.

08 HYDROGEN↗
NASA NTRS2004

Fire Suppression in Low Gravity Using a Cup Burner

Longer duration missions to the moon, to Mars, and on the International Space Station increase the likelihood of accidental fires. The goal of the present investigation is to: (1) understand the physical and chemical processes of fire suppression in various gravity and O2 levels simulating spacecraft, Mars, and moon missions; (2) provide rigorous testing of numerical models, which include detailed combustion-suppression chemistry and radiation sub-models; and (3) provide basic research results useful for advances in space fire safety technology, including new fire-extinguishing agents and approaches.The structure and extinguishment of enclosed, laminar, methane-air co-flow diffusion flames formed on a cup burner have been studied experimentally and numerically using various fire-extinguishing agents (CO2, N2, He, Ar, CF3H, and Fe(CO)5). The experiments involve both 1g laboratory testing and low-g testing (in drop towers and the KC-135 aircraft). The computation uses a direct numerical simulation with detailed chemistry and radiative heat-loss models. An agent was introduced into a low-speed coflowing oxidizing stream until extinguishment occurred under a fixed minimal fuel velocity, and thus, the extinguishing agent concentrations were determined. The extinguishment of cup-burner flames, which resemble real fires, occurred via a blowoff process (in which the flame base drifted downstream) rather than the global extinction phenomenon typical of counterflow diffusion flames. The computation revealed that the peak reactivity spot (the reaction kernel) formed in the flame base was responsible for attachment and blowoff of the trailing diffusion flame. Furthermore, the buoyancy-induced flame flickering in 1g and thermal and transport properties of the agents affected the flame extinguishment limits.

Takahashi, Fumiaki↗