Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure modeling”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Composite Analysis Solid Waste Release Data Reduction of Activity Flux from Waste Sites to the Vadose Zone for Baseline Assessment

This environmental calculation file (ECF) describes the data reduction methodology for the solid waste release datasets generated by the Composite Analysis Solid Waste Release (CASWR) model (CP-62766, Model Package Report: Composite Analysis Solid Waste Release Model [CASWR Model]) and documents the resulting reduced datasets for the updated Composite Analysis (CA) baseline assessment. This ECF is in support of the updated Hanford Site CA. The purpose of the updated Hanford Site CA is to provide an estimate of the cumulative radiological impacts from active and planned low-level radioactive waste disposals and other potentially interacting radioactive waste disposal sources that will remain following Hanford Site closure. The vadose zone modeling activities associated with the CA includes solid waste forms as contaminant sources. To support the modeling of solid waste transport from the vadose zone to the groundwater at the Hanford Site, the CASWR model was designed to generate deterministic radionuclide release rates for the Hanford Site Central Plateau solid waste disposal sites. The complete CASWR model-generated dataset for the CA Baseline Assessment consists of 2,378,831 sets of time and release rates (1,073,869 sets for the 200 East Area and 1,304,962 sets for the 200 West Area), representing a course of the simulation from year 1953 to year 12,069, 16 radiological contaminants of potential concern (COPCs), and 237 solid waste sources as documented in ECFHANFORD- 19-0112, Solid Waste Release Calculations for the Composite Analysis Baseline Assessment. For several of the STOMP vadose zone models, the number of release rate data pairs required for release at all waste sites would exceed the limitations imposed by the STOMP modelling software package used to simulate contaminant transport through the vadose zone. As a result, a data reduction of the CASWR dataset is required prior to its use as a source input for the CA vadose zone modeling activities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1986

Combustion research in the Internal Fluid Mechanics Division

The goal of this research is to bring computational fluid dynamics to a state of practical application for the aircraft engine industry. The approach is to have a strongly integrated computational and experimental program for all the disciplines associated with the gas turbine and other aeropropulsion systems by advancing the understanding of flow physics, heat transfer, and combustion processes. The computational and experimental research is integrated in the following way: the experiments that are performed provide an empirical data set so that physical models can be formulated to describe the processes that are occurring - for example, turbulence or chemical reaction. These experiments also form a data base for those who are doing code development by providing experimental data against which the codes can be verified and assesed. Models are generated as closure to some of the numerical codes, and they also provide physical insight for experiments. At the same time, codes which solve the complete Navier-Stokes equations can be used as a kind of numerical experiment from which far more extensive data can be obtained than ever could be obtained experimentally. This could provide physical insight into the complex processes that are taking place. These codes are also exercised against experimental data to assess the accuracy and applicability of models.

Mularz, Edward J.↗
NASA NTRS2001

Accuracy Assessment of Response Surface Approximations for Supersonic Turbine Design

There is a growing trend to employ CFD tools to supply the necessary information for design optimization of fluid dynamics components/systems. Such results are prone to uncertainties due to reasons including discretization. errors, incomplete convergence of computational procedures, and errors associated with physical models such as turbulence closures. Based on this type of information, gradient-based optimization algorithms often suffer from the noisy calculations, which can seriously compromise the outcome. Similar problems arise from the experimental measurements. Global optimization techniques, such as those based on the response surface (RS) concept are becoming popular in part because they can overcome some of these barriers. However, there are also fundamental issues related to such global optimization technique such as RS. For example, in high dimensional design spaces, typically only a small number of function evaluations are available due to computational and experimental costs. On the other hand, complex features of the design variables do not allow one to model the global characteristics of the design space with simple quadratic polynomials. Consequently a main challenge is to reduce the size of the region where we fit the RS, or make it more accurate in the regions where the optimum is likely to reside. Response Surface techniques using either polynomials or and Neural Network (NN) methods offer designers alternatives to conduct design optimization. The RS technique employs statistical and numerical techniques to establish the relationship between design variables and objective/constraint functions, typically using polynomials. In this study, we aim at addressing issues related to the following questions: (1) How to identify outliers associated with a given RS representation and improve the RS model via appropriate treatments? (2) How to focus on selected design data so that RS can give better performance in regions critical to design optimization? (3) How to combine NN and polynomial techniques for improving the accuracy of the RS model?

Papila, Nilay↗
NASA NTRS1980

Multiple-Time-Scale Concepts in Turbulent Transport Modeling

Progress in developing a turbulence closure employing two or more independently calculated time scales with which to characterize the rates of progress of different turbulent interactions is reported. The energy containing part of the spectrum is divided into two regions which respond at different rates and in different ways to changes in the environment. The scheme may be regarded as providing an intermediate level of approximation between the relatively simple, but fallible, single-point closures and the vastly more elaborate two-point closures which have so far been applied only to simulating homogeneous flows. The proposed approach requires only slightly more computational effort than single-scale schemes. Computational results are reported for several thin shear flows which show striking improvement in the level of agreement with experiment over that obtained with models employing only one time scale.

Hanjalic, K.↗
DOE OSTI2022

Model Package Report: Composite Analysis Solid Waste Release Model (CASWR Model)

This document describes the implementation of a solid waste form release model in GoldSim for the Hanford Site Composite Analysis (CA) Update. This Composite Analysis Solid Waste Release model (CASWR model) was designed to generate deterministic radionuclide release rates for Hanford’s Central Plateau solid waste disposal sites using single realizations of release model coefficients. Five generalized waste form types are used for the conceptual model of waste release: surplus reactor block, cement, soil-debris, grouted residual waste, and ancillary equipment. The surplus reactor block waste form consists of radionuclide leaching from surplus production graphite reactor core blocks via unspecified processes (White et al., 1984 as cited in PNNL-15965). The cement waste form represents solidified wastes whose permeability is much lower than that of the surrounding soil. The soil-debris waste form type is defined as unconsolidated waste mixed with soil material. Tanks and canyon complexes comprise the grouted residual waste form such that their solid waste will be grouted and capped with a surface barrier at the completion of their cleanup. Finally, the ancillary equipment waste form constitutes contaminant releases from ancillary and auxiliary waste form residues associated with tank farms at closure. Individual sub-models are implemented to numerically represent a respective waste form within the CASWR Model: Surplus Reactor Block Sub-model, Cement Sub-model, SoilDebris Sub-model, Grouted Residual Waste Sub-model, and Ancillary Equipment Submodel. Advection is assumed to be the primary transport process governing the release of radionuclides in the Ancillary Equipment and Soil-Debris Sub-models. Diffusion is assumed to be the primary release process in the Grouted Residual Waste and Cement Sub-models. An unspecified zero-order release process is considered in the Surplus Reactor Block Sub-model due to the lack of information regarding actual processes involved in irradiated graphite leaching. The Surplus Reactor Block, Cement, and SoilDebris Sub-models were compared against analytical solutions (PNNL-11800, Composite Analysis for Low-Level Waste Disposal in the 200 Area Plateau of the Hanford Site). The agreement between results of these analytical solutions and the corresponding waste form models verified their correct implementation in GoldSim. A 1-D modeling abstraction approach for the Grouted Residual Waste and Ancillary Equipment Sub-models was adopted from existing Performance Assessment (PA) models. Despite the simplifications made in these sub-models, they were found to be appropriate representations of the waste forms, similar to what was used in the Waste Management Area C PA model (RPP-ENV-58782, Performance Assessment of Waste Management Area C, Hanford Site, Washington, Rev. 0). A sensitivity analysis was conducted to identify the most influential parameters in each waste form sub-model. Suggestions for considering pH-dependent and redox-dependent release mechanisms are formulated through the development of a conditional constant approach in GoldSim.

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

Modelling the Liquid Waste Operation at the Savannah River Site

Modeling Successes Dounreay – Site Closure Optimization and Risk Reduction Re-structuring and streamlining a complex project plan to achieve site decommissioning in reduced time and at reduced cost, through automated task prioritization. Idaho AMWTP – Multi-stream Optimization Speeding up a complex facility with many process lines, by utilizing all of the process lines more efficiently and reducing overall plant downtime, by coordinating repairs and shutdowns. AWE –Waste Repackaging and Relocation Relocating schedule for waste material from an old storage facility to a new purpose-built one, safely, and in a timely fashion that removed the need to reassure the safety of the old facility. Sellafield – Effluent, Sludge and HAW Treatment Improving a wide range of chemical engineered processes including aqueous effluent treatment, chemically-reactive sludge management and waste vitrification processing. Strategic Petroleum Reserve – Cavern Management Linking chemistry models of petroleum prediction, through engineering cavern management to financial modeling of gasoline sales to benefit US government

Jung, Andrew W.↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
DOE OSTI2020

Data from Wan et al. (2020, JAMES, DOI: 10.1029/2019MS001982)

This data set includes the simple condensation model output and CAM4 model output used to generate figures presented in the following paper: Wan, H., Woodward, C. S., Zhang, S., Vogl, C. J., Stinis, P., Gardner, D. J., Rasch, P. J., Zeng, X., Larson, V. E., and Singh, B. (2020): Improving time-step convergence in an atmosphere model with simplified physics: the impacts of closure assumption and process coupling, Journal of Advances in Modeling Earth Systems, https://doi.org/10.1029/2019MS001982

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

CFD validation of subsonic turbulent planar shear layers

The primary objective of the present study is to assess the limitations and capabilities of RPLUS in predicting the entrainment, mixing, and burning characteristics of a high subsonic planar shear layer; these computations are performed in conjunction with an experiment being conducted at NASA Lewis Research Center. Turbulence is modeled by a two-equation k-epsilon closure modified for compressibility effects. Combustion is modeled by two different models: one in which a finite-rate laminar chemistry model is used for H2-O2 oxidation, and the other based on a composition joint pdf approach to account for the turbulence-chemistry interaction effects. Between upwind and centered differencing schemes, minimal differences are observed for the flow properties of a nonreacting shear layer case studied. In comparison with experimental data, computed growth rates are underpredicted in both reacting and nonreacting cases of a shear layer, but more severely for the reacting flows. However, the qualitative agreement between computation and experiment is reasonably good for the other observable characteristics, and indicative of potential for further improvements.

Lai, H. T.↗
NASA NTRS1978

Space ecosynthesis: An approach to the design of closed ecosystems for use in space

The use of closed ecological systems for the regeneration of wastes, air, and water is discussed. It is concluded that such systems, if they are to be used for the support of humans in space, will require extensive mechanical and physico-chemical support. The reason for this is that the buffering capacity available in small systems is inadequate, and that natural biological and physical regulatory mechanisms rapidly become inoperative. It is proposed that mathematical models of the dynamics of a closed ecological system may provide the best means of studying the initial problems of ecosystem closure. A conceptual and mathematical model of a closed ecosystem is described which treats the biological components as a farm, calculates the rates of flow of elements through the system by mass-balance techniques and control theory postulates, and can evaluate the requirements for mechanical buffering activities. It is suggested that study of the closure of ecosystems can significantly aid in the establishment of general principles of ecological systems.

Macelroy, R. D.↗
DOE OSTI2026

Wyoming CarbonSAFE Phase III: Site Characterization and Permitting Commercial-Scale Carbon Storage Complex Feasibility Study at Dry Fork Station, Wyoming

This report presents the findings of the technical and non-technical site characterization and permitting activities (“Phase III”) conducted under the Wyoming CarbonSAFE: Accelerating CCUS Commercialization and Deployment at Dry Fork Power Station (DFS) and the Wyoming Integrated Test Center project (“Wyoming CarbonSAFE”). Wyoming CarbonSAFE is part of the U.S. Department of Energy’s (DOE) National Energy Technology Laboratory (NETL) Carbon Storage Assurance Facility Enterprise (“CarbonSAFE”) Initiative. The results of Phase III demonstrate that the Wyoming CarbonSAFE project - referred to as the Northern Powder River Basin Carbon Sequestration Hub (NPRB-CSH) - meets the technical, regulatory, and commercial feasibility requirements necessary to advance toward commercial development and construction. The activities completed under this project Phase make the NPRB-CSH one of the region’s most commercially ready carbon storage sites. Completion of this Phase included the finalization of all site surface and subsurface characterization activities, completion and testing of two Class VI standard wells, 10 draft Class VI permits-to construct to address the future needs of a storage complex, finalized NEPA assessments, transportation and capture FEED studies, and a full commercialization strategy with economic modeling, operation and site closure strategies. The NPRB-CSH meets all requirements to progress to a CarbonSAFE Phase IV program or advance to full commercial operations under the development of a business partner.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Improbability of Nuclear Criticality in Compacted Criticality Control Overpacks after Room Closure by Salt Creep at Waste Isolation Pilot Plant

Based on the rationale presented, nuclear criticality is improbable after salt creep causes compaction of criticality control overpacks (CCOs) disposed at the Waste Isolation Pilot Plant, an operating repository in bedded salt for the disposal of transuranic (TRU) waste from atomic energy defense activities. For most TRU waste, the possibility of post-closure criticality is exceedingly small either because the salt neutronically isolates TRU waste canisters or because closure of a disposal room from salt creep does not sufficiently compact the low mass of fissile material. The criticality potential has been updated here because of the introduction of CCOs, which may dispose up to 380 fissile gram equivalent plutonium-239 in each container. The criticality potential is evaluated through high-fidelity geomechanical modeling of a disposal room filled with CCOs during two representative conditions: (1) large salt block fall, and (2) gradual salt compaction (without brine seepage and subsequent gas generation to permit maximum room closure). Geomechanical models of rock fall demonstrate three tiers of CCOs are not greatly disrupted. Geomechanical models of gradual room closure from salt creep predict irregular arrays of closely packed CCOs after 1000 years, when room closure has asymptotically approached maximum compaction. Criticality models of spheres and cylinders of 380 fissile gram equivalent of plutonium (as oxide) at the predicted irregular spacing demonstrate that an array of CCOs is not critical when surrounded by salt and magnesium oxide, provided the amount of hydrogenous material shipped in the CCO (usually water and plastics) is controlled or boron carbide (a neutron poison) is mixed with the fissile contents.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2010

Alfven Wave Reflection Model of Field-Aligned Currents at Mercury

An Alfven Wave Reflection (AWR) model is proposed that provides closure for strong field-aligned currents (FACs) driven by the magnetopause reconnection in the magnetospheres of planets having no significant ionospheric and surface electrical conductance. The model is based on properties of the Alfven waves, generated at high altitudes and reflected from the low-conductivity surface of the planet. When magnetospheric convection is very slow, the incident and reflected Alfven waves propagate along approximately the same path. In this case, the net field-aligned currents will be small. However, as the convection speed increases. the reflected wave is displaced relatively to the incident wave so that the incident and reflected waves no longer compensate each other. In this case, the net field-aligned current may be large despite the lack of significant ionospheric and surface conductivity. Our estimate shows that for typical solar wind conditions at Mercury, the magnitude of Region 1-type FACs in Mercury's magnetosphere may reach hundreds of kilo-Amperes. This AWR model of field-aligned currents may provide a solution to the long-standing problem of the closure of FACs in the Mercury's magnetosphere. c2009 Elsevier Inc. All rights reserved.

Lyatsky, Wladislaw↗
DOE OSTI2021

A Time-Dependent Directional Damage Theory for Brittle Rocks Considering the Kinetics of Microcrack Growth

In this work, a novel microcrack damage theory for describing the time-dependent behavior of brittle rocks is proposed. Instead of using scalar or tensorial variables to approximate the effect of microcracks, the concept of directional damage is introduced by upscaling the directional distribution density of microcracks through the noninteracting homogenization scheme. The contribution of an individual crack on global strain is first derived from the crack-opening displacements including sliding, closure, and dilation. The model is then cast in the free energy form, and the range of parameters ensuring compliance with the second law of thermodynamics is derived rigorously. Considering that the subcritical propagation of microcracks is the dominant mechanism behind the time-dependent deformation and failure of brittle rocks, the driving force for directional damage is defined as the corresponding energy release rate, and the effect of time is introduced through a damage evolution law inspired from the kinetics of subcritical crack growth. The model is applied to simulate the behavior of basalt rock under various loading conditions. In addition to predicting the stress–strain–time response, the model offers enhanced resolution in representing the anisotropic characteristics of microcrack-induced damage in brittle rocks and their time-dependent evolutions.

58 GEOSCIENCES↗
DOE OSTI2024

Revisiting a Drag Partition Model For Canopy-Like Roughness Elements

Turbulent flows over a large surface area (S) covered by n obstacles experience an overall drag due to the presence of the ground and the protruding obstacles into the flow. The drag partition between the roughness obstacles and the ground is analyzed using an analytical model proposed by Raupach (Boundary-Layer Meteorol 60:375-395, 1992) and is hereafter referred to as R92. The R92 is based on the premise that the wake behind an isolated roughness element can be described by a shelter area A and a shelter volume V. The individual sizes of A and V without any interference from other obstacles can be determined from scaling analysis for the spread of wakes. To upscale from an individual roughness element to n/S elements where wakes may interact, R92 adopted a background stress re-normalizing instead of reducing A or V with each element addition. This work demonstrates that R92’s approach results in a linear background stress reduction in A and V only when the ratio of n/S is small, due to a low probability of wake interactions. This probabilistic nature suggests that up-scaling from individual to multiple roughness elements can be re-formulated using stochastic averaging methods proposed here. The two approaches are shown to recover R92 under plausible conditions. An alternative scaling for the shelter volume is also proposed here using thermodynamic arguments of work and dissipation though the final outcome remains similar to R92. Comparisons between R92 and available data spanning more than two decades after R92 on blocks and vegetation-like roughness elements confirm the practical utility of R92. The agreement between R92 and this updated databases of experiments and simulations confirm the potential use of R92 in large-scale models provided that the relevant parameters accommodate certain features of the roughness element type (cube versus vegetation-like) and, to a lesser extent, their configuration throughout S. Last, a comparison between R92 and models based on first-order closure principles with constant mixing length suggests that R92 can outperform such models when evaluated across a wide range of roughness densities.

canopy turbulence↗
NASA NTRS1975

Computation of separated transonic turbulent flows

The two-dimensional Reynolds-averaged compressible Navier-Stokes equations are solved using MacCormack's second-order-accurate explicit finite difference method to simulate the separated transonic turbulent flow field over an airfoil. Four different algebraic eddy viscosity models are tested for viability to achieve turbulence closure for the class of flows considered. These models range from an unmodified boundary-layer mixing-length model to a relaxation model incorporating special considerations for the separation bubble region. Results of this study indicate the necessity for special attention to the separated flow region and suggest limits of applicability of algebraic turbulence models to these separated flow fields.

Deiwert, G. S.↗
NASA NTRS1981

Numerical experiments in homogeneous turbulence

The direct simulation methods developed by Orszag and Patternson (1972) for isotropic turbulence were extended to homogeneous turbulence in an incompressible fluid subjected to uniform deformation or rotation. The results of simulations for irrotational strain (plane and axisymmetric), shear, rotation, and relaxation toward isotropy following axisymmetric strain are compared with linear theory and experimental data. Emphasis is placed on the shear flow because of its importance and because of the availability of accurate and detailed experimental data. The computed results are used to assess the accuracy of two popular models used in the closure of the Reynolds-stress equations. Data from a variety of the computed fields and the details of the numerical methods used in the simulation are also presented.

Rogallo, R. S.↗
NASA NTRS1984

Direct numerical simulation of reacting flows

The objectives of this work are: (1) to extend the technique of direct numerical simulations to turbulent, chemically reacting flows, (2) to test the validity of the method by comparing computational results with laboratory data, and (3) to use the simulations to gain a better understanding of the effects of turbulence on chemical reactions. The effects of both the large scale structure and the smaller scale turbulence on the overall reaction rates are addressed. The relationship between infinite reaction rate and finite reaction rate chemistry is compared with some of the results of calculations with existing theories and laboratory data. The direct numerical simulation method involves the numerical solution of the detailed evolution of the complex turbulent velocity and concentration fields. Using very efficient numerical methods (e.g., pseudospectral methods), the fully nonlinear (possibly low pass filtered) equations of motion are solved and no closure assumptions or turbulence models are used. Statistical data are obtained by performing spatial, temporal, and/or ensemble averages over the computed flow fields.

Riley, J. J.↗