Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “closure 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,027 records · Page 57

Sensitivity Analysis of Wind and Turbulence Predictions With Mesoscale‐Coupled Large Eddy Simulations Using Ensemble Machine Learning

Abstract Coupling between mesoscale models and large‐eddy simulation (LES) models is increasingly used to more realistically represent the wide range of scales of atmospheric motions affecting boundary layer winds and turbulence that need to be simulated accurately for applications such as wind energy. However, such mesoscale‐to‐microscale coupled modeling frameworks are potentially affected by a large number of uncertain closure parameters. Here, we investigate the sensitivity associated with six closure parameters related to a 1.5‐order subgrid‐scale turbulence closure for an ensemble of mesoscale‐coupled LES. The simulations are performed using the Weather Research and Forecasting model nested from horizontal resolutions of greater than a kilometer down to tens of meters. Closure parameters are varied to generate perturbed parameter ensembles for two case studies of highly sheared, convective boundary layers observed in the Columbia Basin of Oregon and Washington during the Second Wind Forecast Improvement Project. Machine learning algorithms are used to explore the sensitivity of LES predictions, considering the effects of the perturbed physical parameters alongside categorical factors such as the case study identity, measurement location, and LES resolution. For the conditions we examine, a single parameter, the eddy viscosity coefficient, is the dominant source of parametric sensitivity and its importance is comparable to the categorical factors for several of the simulation response variables we examine.

54 ENVIRONMENTAL SCIENCES↗

2016 KIVA-hpFE Development: A Robust and Accurate Engine Modeling Software

Los Alamos National Laboratory and its collaborators are facilitating engine modeling by improving accuracy and robustness of the modeling, and improving the robustness of software. We also continue to improve the physical modeling methods. We are developing and implementing new mathematical algorithms, those that represent the physics within an engine. We provide software that others may use directly or that they may alter with various models e.g., sophisticated chemical kinetics, different turbulent closure methods or other fuel injection and spray systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Overview of NASA's Integrated Design and Engineering Analysis (IDEA)Environment

Historically, the design of subsonic and supersonic aircraft has been divided into separate technical disciplines (such as propulsion, aerodynamics and structures) each of which performs their design and analysis in relative isolation from others. This is possible in most cases either because the amount of interdisciplinary coupling is minimal or because the interactions can be treated as linear. The design of hypersonic airbreathing vehicles, like NASA s X-43, is quite the opposite. Such systems are dominated by strong non-linear interactions between disciplines. The design of these systems demands that a multi-disciplinary approach be taken. Furthermore, increased analytical fidelity at the conceptual design phase is highly desirable as many of the non-linearities are not captured by lower fidelity tools. Only when these systems are designed from a true multi-disciplinary perspective can the real performance benefits be achieved and complete vehicle systems be fielded. Toward this end, the Vehicle Analysis Branch at NASA Langley Research Center has been developing the Integrated Design & Engineering Analysis (IDEA) Environment. IDEA is a collaborative environment for parametrically modeling conceptual and preliminary launch vehicle configurations using the Adaptive Modeling Language (AML) as the underlying framework. The environment integrates geometry, configuration, propulsion, aerodynamics, aerothermodynamics, trajectory, closure and structural analysis into a generative, parametric, unified computational model where data is shared seamlessly between the different disciplines. Plans are also in place to incorporate life cycle analysis tools into the environment which will estimate vehicle operability, reliability and cost. IDEA is currently being funded by NASA s Hypersonics Project, a part of the Fundamental Aeronautics Program within the Aeronautics Research Mission Directorate. The environment is currently focused around a two-stage-to-orbit configuration with a turbine based combined cycle (TBCC) first stage and reusable rocket second stage. This paper provides an overview of the development of the IDEA environment, a description of the current status and detail of future plans.

Robinson, Jeffrey S.↗

Present-day kinematics of the Rivera plate and implications for tectonics in southwestern Mexico

A model for the present-day motion of the Rivera plate relative to the North America, Cocos, and Pacific plates is derived using new data from the Pacific-Rivera rise and Rivera transform fault, together with new estimates of Pacific-Rivera motions. The results are combined with the closure-consistent NUVEL-1 global plate motion model of DeMets et al. (1990) to examine present-day deformation in southwestern Mexico. The analysis addresses several questions raised in previous studies of the Rivera plate. Namely, do plate motion data from the northern East Pacific rise require a distinct Rivera plate? Do plate kinematic data require the subduction of the Rivera plate along the seismically quiescent Acapulco trench? If so, what does the predicted subduction rate imply about the earthquake recurrence interval in the Jalisco region of southwestern Mexico?

Demets, Charles↗

Refinement of Modeling Techniques for the Structural Evaluation of Hanford Single-Shell Nuclear Waste Storage Tanks

Abstract: A total of 149 tanks out of 177 at the Hanford Site (in Washington State, USA) belong to the first generation of underground nuclear waste storage tanks known as single shell tanks (SSTs). These tanks were constructed between 1943 and 1964 and are well beyond their design life. All the SSTs had been removed from active service by November 1980 and have been later interim stabilized by removing the pumpable liquids. The remaining waste in the tanks is in the form of salt cake and sludge awaiting r permanent disposal.. The evaluation of the structural integrity of these tanks is of utmost importance not only for the continued safe storage of the waste until waste retrieval and closure, but also to assure safe retrieval and closure operations. This article discusses the structural analysis approach, modeling challenges and issues encountered during the ongoing analysis of record (AOR) for evaluating the structural integrity of the SSTs. There are several geometrical and material nonlinearities and uncertainties to be dealt with while performing the modern finite element analysis of these tanks. Several studies were conducted to refine the models in order to minimize modeling artifacts introduced by soil arching, boundary effects, concrete cracking, and concrete-soil interface behavior. The analysis takes into account the temperature history of the tanks and allowable mechanical operating loads of these tanks for proper estimation of creep strains and thermal degradation of material properties. The loads imposed in the AOR models also include anticipated loads that these tanks may see during waste retrieval and closure. Due to uncertainty in a number of inputs to the models, sensitivity studies were conducted to address questions related to the boundary conditions to realistically or conservatively represent the influence of surrounding tanks in a tank farm, the influence of backfill excavation slope, the extent of backfill and the total extent of undisturbed soil surrounding the backfill. The article also discusses the criteria and design standards used for evaluating the structural integrity of these underground concrete tanks. Because of the non-availability of complete data on the thermal and operating history for many of the individual tanks, some of the data was assumed or interpolated. However, the models developed for the analysis of record represent the bounding scenarios and include the worst and extreme loading cases that the tanks were subjected to or anticipated. The modeling refinement techniques followed in the AOR resulted in conservative estimates for force and moment demands at various sections in the concrete tanks. The SSTs are classified into 4 types as per their configuration and capacity. This article discusses the modeling aspects related to two types of SSTs that have been analyzed until now. The TOLA results combined with seismic demands from seismic analysis for the analysis of record indicate that the tanks analyzed are structurally stable as per the evaluation criteria established. These results are presented in a separate article. The modeling techniques, methodology and evaluation criteria developed for evaluating the structural integrity of SSTs at Hanford are in general applicable to any similar tanks or underground concrete storage structures.

Single Shell Tank, Structural analysis, Structural↗

Refinement of Modeling Techniques for the Structural Evaluation of Hanford Single-Shell Nuclear Waste Storage Tanks

ABSTRACT Several tanks at the Hanford Site (in Washington State, USA) belong to the first generation of underground nuclear waste storage tanks known as single shell tanks (SSTs). These tanks were constructed between 1943 and 1964 and are well beyond their design life. This article discusses the structural analysis approach and modeling challenges encountered during the ongoing analysis of record (AOR) for evaluating the structural integrity of the SSTs. There are several geometrical and material nonlinearities and uncertainties to be dealt with while performing the modern finite element analysis of these tanks. The analysis takes into account the temperature history of the tanks and allowable mechanical operating loads of these tanks for proper estimation of creep strains and thermal degradation of material properties. The loads prescribed in the AOR models also include anticipated loads that these tanks may see during waste retrieval and closure. Due to uncertainty in a number of inputs to the models, sensitivity studies were conducted to address questions related to the boundary conditions to realistically or conservatively represent the influence of surrounding tanks in a tank farm, the influence of backfill excavation slope, the extent of backfill and the total extent of undisturbed soil surrounding the backfill. Because of the limited availability of data on the thermal and operating history for many of the individual tanks, some of the data was assumed or interpolated. However, the models developed for the analysis of record represent the bounding scenarios and include the loading conditions that the tanks were subjected to or anticipated. The modeling refinement techniques followed in the AOR resulted in conservative estimates for force and moment demands at various sections in the concrete tanks. This article discusses the modeling aspects related to Type-II and Type-III SSTs. The modeling techniques, methodology and evaluation criteria developed for evaluating the structural integrity of SSTs at Hanford are in general applicable to any similar tanks or underground concrete storage structures.

Waste Managament, Single Shell Tanks, Structural E↗

Proposed Semi-Analytic Benchmark for Coupled Neutronics/Thermal-Hydraulics Feedback in Multiphysics Simulations

The past decade has seen considerable advancements in multiscale and multiphysics simulations that account for nonlinear effects of feedback between different physical processes within a single calculation. One area of particular interest has been multiphysics simulations that include both neutronics and thermal-hydraulics effects to accurately model the thermal feedback effects present in nuclear reactors. During this time, many groups have successfully modeled thermal feedback effects by coupling a variety of different types of neutron transport and thermal-hydraulic solvers together. However, formal verification and validation of these coupled code systems has proven challenging due to a lack of benchmark problems with well-characterized reference solutions. One approach is to confirm gross energy balance using bulk coolant properties (outlet temperature/enthalpy) for code verification purposes. Alternatively, code-to-code comparisons of results for a predefined reference/challenge problem have been used to establish reasonability of local results (neutron flux, heating, enthalpy rise) for various implementations, for example the VERA Core Physics Progression Problems. While community solutions are valuable for assessing the relative performance and fidelity of different implementations, they are not guaranteed to provide a measure of the absolute accuracy of a given implementation. In fact, these code-to-code comparisons can give a false sense of confidence by masking the effects of errors or approximations that are common to all of the solution methods tested. This is a particular concern when many methods derive from a common approach, such as Picard iteration. To this end, an analytical benchmark would be very valuable for verifying individual coupled physics implementations, as well as for providing a formal assessment of algorithm accuracy and sensitivity to model approximations such as mesh size and number of convergence iterations. However, creation of an analytical benchmark is difficult due to the complex nature of the eigenvalue form of the neutron transport equation (which does not permit a simple analytical solution for non-trivial cases) as well as the variety of different empirical closure relationships used in thermal-hydraulic modeling. In this paper we propose a benchmark for coupled neutronics/thermal-hydraulics simulations. The proposed benchmark overcomes the challenges associated with a multiphysics benchmark problem by using a fixed-source (rather than eigenvalue) formulation of the neutron transport equation and limiting fluid properties to the single-phase regime where the density and temperature of the coolant are related to enthalpy by a known relationship, while still preserving significant nuclear/thermal feedback effects within the simplified problem. The new benchmark is easy to simulate using existing coupled physics implementations, and allows for a direct comparison with a reference solution, as provided in this summary.

42 ENGINEERING↗

Defining Modeling and Analysis Complexity Needs for Remediation - 18425

Many subsurface remedy decisions rely on model predictions or other types of data analyses to support site characterization and remedy selection, design, implementation, performance assessment, and closure. Thus, determining the appropriate level of complexity for modeling and data analysis is important to facilitate subsurface remedy decisions. Multiple factors are important in selecting the appropriate level of complexity in modeling and data analysis. For instance, the type and number of inputs required typically increases with increased complexity of analysis, and the ability to obtain or verify those inputs needs to be considered. The approach needs to establish the relationship between parameters or processes and contaminant transport and the detail of the description needed to meet the analysis objectives. Assessing the modeling complexity associated with the appropriate level of detail can be facilitated by considering complexity of the 1) modeling approach, 2) description of process, and 3) description of heterogeneity. As demonstrated by case studies, using a framework to guide determination of modeling and analysis complexity facilitates selection of the appropriate approach for the site situation. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Performance Assessment Program for the Savannah River Site Liquid Waste Facilities - 13610

The Liquid Waste facilities at the U.S. Department of Energy's (DOE) Savannah River Site (SRS) are operated by Liquid Waste Operations contractor Savannah River Remediation LLC (SRR). A separate Performance Assessment (PA) is prepared to support disposal operations at the Saltstone Disposal Facility and closure evaluations for the two liquid waste tank farm facilities at SRS, F-Tank Farm and H-Tank Farm. A PA provides the technical basis and results to be used in subsequent documents to demonstrate compliance with the pertinent requirements identified in operations and closure regulatory guidance. The Saltstone Disposal Facility is subject to a State of South Carolina industrial solid waste landfill permit and the tank farms are subject to a state industrial waste water permit. The three Liquid Waste facilities are also subject to a Federal Facility Agreement approved by the State, DOE and the Environmental Protection Agency (EPA). Due to the regulatory structure, a PA is a key technical document reviewed by the DOE, the State of South Carolina and the EPA. As the waste material disposed of in the Saltstone Disposal Facility and the residual material in the closed tank farms is also subject to reclassification prior to closure via a waste determination pursuant to Section 3116 of the Ronald W. Reagan National Defense Authorization Act of Fiscal Year 2005, the U.S. Nuclear Regulatory Commission (NRC) is also a reviewing agency for the PAs. Pursuant to the Act, the NRC also has a continuing role to monitor disposal actions to assess compliance with stated performance objectives. The Liquid Waste PA program at SRS represents a continual process over the life of the disposal and closure operations. When the need for a PA or PA revision is identified, the first step is to develop a conceptual model to best represent the facility conditions. The conceptual model will include physical dimensions of the closed system, both the engineered and natural system, and modeling input parameters associated with the modeled features, both initial values (at the time of facility closure) and degradation rates/values. During the development of the PA, evaluations are conducted to reflect not only the results associated with the best available information at the time but also to evaluate potential uncertainties and sensitivities associated with the modeled system. While the PA will reflect the modeled system results from the best available information, it will also identify areas for future work to reduce overall PA uncertainties moving forward. DOE requires a PA Maintenance Program such that work continues to reduce model uncertainties, thus bolstering confidence in PA results that support regulatory decisions. This maintenance work may include new Research and Development activities or modeling as informed by previous PA results and other new information that becomes available. As new information becomes available, it is evaluated against previous PAs and appropriate actions are taken to ensure continued confidence in the regulatory decisions. Therefore, the PA program is a continual process that is not just the development of a PA but seeks to incorporate new information to reduce overall model uncertainty and provide continuing confidence in regulatory decisions. (author)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The Water Balance Representation in Urban‐PLUMBER Land Surface Models

Abstract Urban Land Surface Models (ULSMs) simulate energy and water exchanges between the urban surface and atmosphere. However, earlier systematic ULSM comparison projects assessed the energy balance but ignored the water balance, which is coupled to the energy balance. Here, we analyze the water balance representation in 19 ULSMs participating in the Urban‐PLUMBER project using results for 20 sites spread across a range of climates and urban form characteristics. As observations for most water fluxes are unavailable, we examine the water balance closure, flux timing, and magnitude with a score derived from seven indicators expecting better scoring models to capture the latent heat flux more accurately. We find that the water budget is only closed in 57% of the model‐site combinations assuming closure when annual total incoming fluxes (precipitation and irrigation) fluxes are within 3% of the outgoing (all other) fluxes. Results show the timing is better captured than magnitude. No ULSM has passed all water balance indicators for any site. Models passing more indicators do not capture the latent heat flux more accurately refuting our hypothesis. While output reporting inconsistencies may have negatively affected model performance, our results indicate models could be improved by explicitly verifying water balance closure and revising runoff parameterizations. By expanding ULSM evaluation to the water balance and related to latent heat flux performance, we demonstrate the benefits of evaluating processes with direct feedback mechanisms to the processes of interest.

Jongen, H. J.↗

Geomechanical Simulation of Bayou Choctaw Strategic Petroleum Reserve - Model Calibration.

A finite element numerical analysis model has been constructed that consists of a realistic mesh capturing the geometries of Bayou Choctaw (BC) Strategic Petroleum Reserve (SPR) site and multi - mechanism deformation ( M - D ) salt constitutive model using the daily data of actual wellhead pressure and oil - brine interface. The salt creep rate is not uniform in the salt dome, and the creep test data for BC salt is limited. Therefore, the model calibration is necessary to simulate the geomechanical behavior of the salt dome. The cavern volumetric closures of SPR caverns calculated from CAVEMAN are used for the field baseline measurement. The structure factor, A 2 , and transient strain limit factor, K 0 , in the M - D constitutive model are used for the calibration. The A 2 value obtained experimentally from the BC salt and K 0 value of Waste Isolation Pilot Plant (WIPP) salt are used for the baseline values. T o adjust the magnitude of A 2 and K 0 , multiplication factors A2F and K0F are defined, respectively. The A2F and K0F values of the salt dome and salt drawdown skins surrounding each SPR cavern have been determined through a number of back fitting analyses. The cavern volumetric closures calculated from this model correspond to the predictions from CAVEMAN for six SPR caverns. Therefore, this model is able to predict past and future geomechanical behaviors of the salt dome, caverns, caprock , and interbed layers. The geological concerns issued in the BC site will be explained from this model in a follow - up report .

58 GEOSCIENCES↗

Development and Documentation of kL-Based Linear, Nonlinear, and Full Reynolds Stress Turbulence Models

The development and implementation of kL-based Reynolds Average Navier-Stokes (RANS) turbulence models are reported herein. The kL is based on Abdol-Hamid's closure and Menter's modi cation to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's kL equation is well suited for term-by-term modeling and displays useful features compared to other scale formulation. One of the important di erences is the inclusion of higher order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady ows in URANS mode. The present report documents the formulation of three model levels of turbulence models as implemented in the CFD code FUN3D. Methodology and calibration examples are shown in detail. The levels are the linear k-kL and the two-equation algebraic Reynolds stress model (ARSM) as well as the full Reynolds Stress Model (RSM). Attached, separated and corner ow cases are documented and compared with experimental, theoretical and other turbulence model data. The results show generally very good comparisons with canonical and experimental data. The results from this formulation are similar or better than results using the SST two- equation turbulence model. ARSM shows great promise with similar level of computational resources as general two equation turbulence models.

Abdol-Hamid, Khaled S.↗

Geomechanical Model Calibration Using Field Measurements for a Petroleum Reserve

A finite element numerical analysis model has been constructed that consists of a mesh that effectively captures the geometries of Bayou Choctaw (BC) Strategic Petroleum Reserve (SPR) site and multimechanism deformation (M-D) salt constitutive model using the daily data of actual wellhead pressure and oil–brine interface location. The salt creep rate is not uniform in the salt dome, and the creep test data for BC salt are limited. Therefore, the model calibration is necessary to simulate the geomechanical behavior of the salt dome. The cavern volumetric closures of SPR caverns calculated from CAVEMAN are used as the field baseline measurement. The structure factor, A 2 , and transient strain limit factor, K o , in the M-D constitutive model are used for the calibration. The value of A 2 , obtained experimentally from BC salt, and the value of K o , obtained from Waste Isolation Pilot Plant salt, are used for the baseline values. To adjust the magnitude of A 2 and K0, multiplication factors A 2 F and K o F are defined, respectively. The A 2 F and K0F values of the salt dome and salt drawdown skins surrounding each SPR cavern have been determined through a number of back analyses. The cavern volumetric closures calculated from this model correspond to the predictions from CAVEMAN for six SPR caverns. Therefore, this model is able to predict behaviors of the salt dome, caverns, caprock, and interbed layers. In conclusion, the geotechnical concerns associated with the BC site from this analysis will be explained in a follow-up paper.

02 PETROLEUM↗

Final Report to the U.S. Department of Energy for studies of Evaluation of Turbulence Parameterizations for Cloud-Resolving Models

The intermediately-prognostic higher-order turbulence closure (IPHOC) introduces a joint double-Gaussian distribution of liquid water potential temperature (θ l ), total water mixing ratio (q t ), and vertical velocity (w ) to represent any skewed turbulence circulations .The distribution is inferred from the first-, second-, and third-order moments of the variables given above, and is used to diagnose cloud fraction and grid-mean liquid water mixing ratio, as well as the buoyancy and fourth-order terms in the equations describing the evolution of the second- and third-order moments. Only three third-order moments (those of θ l , q t , and w ) are predicted in the IPHOC.

99 GENERAL AND MISCELLANEOUS↗

A minimal model of self-sustaining turbulence

In this work, we examine the turbulence maintained in a Restricted Nonlinear (RNL) model of plane Couette flow. This model is a computationally efficient approximation of the second order statistical state dynamics obtained by partitioning the flow into a streamwise averaged mean flow and perturbations about that mean, a closure referred to herein as the RNL{sub ∞} model. The RNL model investigated here employs a single member of the infinite ensemble that comprises the covariance of the RNL{sub ∞} dynamics. The RNL system has previously been shown to support self-sustaining turbulence with a mean flow and structural features that are consistent with direct numerical simulations (DNS). Regardless of the number of streamwise Fourier components used in the simulation, the RNL system’s self-sustaining turbulent state is supported by a small number of streamwise varying modes. Remarkably, further truncation of the RNL system’s support to as few as one streamwise varying mode can suffice to sustain the turbulent state. The close correspondence between RNL simulations and DNS that has been previously demonstrated along with the results presented here suggest that the fundamental mechanisms underlying wall-turbulence can be analyzed using these highly simplified RNL systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A boundary-layer model for Mars - Comparison with Viking lander and entry data

A 1D boundary-layer model of Mars based on a momentum equation that describes friction, pressure gradient, and Coriolis forces is presented. Frictional forces and convective heating are computed using the level-2 turbulence closure theory of Mellor and Yamada (1974). The model takes into account the radiative effects of CO2 gas and suspended dust particles. Both radiation and convection depend on surface temperatures which are computed from a surface heat budget. Model predictions are compared with available observations from Viking landers. It is concluded that, in general, the model reproduces the basic features of the temperature data. The agreement is particularly good at entry time for the V L-2 site, where the model and observations are within several degrees at all levels for which data are available.

Haberle, Robert M.↗

A passive model for the evolution of subgrid-scale instabilities in turbulent flow regimes

In this work, a zero-dimensional model for the development of turbulent mixing layers driven by Rayleigh–Taylor (RT) or Richtmyer–Meshkov (RM) instabilities is presented, for application in Reynolds-Averaged Navier–Stokes (RANS) simulations. The model approximates the behavior of the turbulent mixing layer at a passively advected tracer particle that sits on the material interface of interest in the simulation, and is used to provide initial conditions for the RANS model once the mixing layer width has developed to a size that may be resolved on the computational mesh. The model relies on approximations to the Besnard–Harlow–Rauenzahn (BHR) RANS closure equations and is thus consistent the turbulent development of the mixing layer, unlike previous models of this type which have focused on the early-time laminar development of RT or RM instabilities. These instabilities will often grow through an initially laminar state before developing into turbulence, and criteria for switching from an existing laminar model to the proposed turbulent model are considered. The proposed model reproduces the mixing layer development of RT and RM instabilities seen in mixed width resolving BHR RANS simulations to a reasonable degree of accuracy, and relaxes the need for the RANS simulation to resolve the mixed width as soon as the mixing layer becomes turbulent, which may be impractical in cases where the mixing layer grows through a wide range of scales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Advancing representations of turbulence in Earth System Models [Slides]

We have created a suite of coastally focused Large Eddy Simulation (LES), including a first ever investigation of the influence of the ocean bottom on turbulence. We have added important forcing not previously considered in LES models. We have evaluated a new more physically realistic and advanced closure against our LES database of simulations. Test cases will constrain constants in the new model and help us transition the mixing model from the deep ocean to the coastal ocean for the first time.

58 GEOSCIENCES↗