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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

DECOVALEX-2023: Task F2 Salt Final Report

The subject of Task F of DECOVALEX-2023 concerns performance assessment modelling of radioactive waste disposal in deep mined repositories. The primary objectives of Task F are to build confidence in the models, methods, and software used for performance assessment (PA) of deep geologic nuclear waste repositories, and/or to bring to the fore additional research and development needed to improve PA methodologies. In Task F2- (salt), these objectives have been accomplished through staged development and comparison of the models and methods used by participating teams in their PA frameworks. Coupled-process submodels and deterministic simulations of the entire PA model for a reference scenario for waste disposal in domal salt have been conducted. The task specification has been updated continuously since the initiation of the project to reflect the staged development of the conceptual repository model and performance metrics. Thermal, hydrological, mechanical, and chemical properties of individual components of the engineered and natural system were chosen for relevance by participating teams. The salt reference case system was characterized using data and measurements collected at relevant underground research laboratories (URLs), field sites, and simulation results from teams with specialized modelling capability. Participating teams made a wide range of model assumptions from compartmentalized networks to full 3D models of the salt formation. No single contributed model includes full-fidelity representation of all the features, events, and processes (FEPs) detailed in the task specification, but almost all features and processes are represented in at least one model. Despite differences in the modelling strategies developed by participating teams, all models indicate that salt compaction and radionuclide diffusion are key processes in the repository, and for the FEPs and model scenario considered, little of the disposed radionuclides will migrate beyond the repository seal over the 100,000 year simulations. In general, the model output quantities have the largest differences over the short term and near the waste. The models tend to be more similar further from waste and at later time. Disparities between the models are believed to be due to differing simplifications from the task specification, some of which are chosen simplifications to reduce complexity, and some are restrictions imposed by the modelling tools. A second round of this task has been accepted for DECOVALEX-2027 in conjunction with Task F1 on crystalline PA modelling. The future round includes waste package heating, improved modelling of salt creep closure, additional comparisons of coupled-process sub-models, and the impact of repository engineering design on radionuclide migration in the repository. Participants will also propose and finalize a set of uncertain inputs for the reference case simulations, propagate these uncertainties in a set of realizations, and conduct sensitivity analyses on the simulation results.

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

New results on the realizability of Reynolds stress turbulence closures

The realizability of Reynolds stress models in homogeneous turbulence is critically assessed from a theoretical standpoint. It is proven that a well known second-order closure formulated by Shih and Lumley using the strong realizability constraints of Schumann is, in fact, not a realizable model. The problem arises from the failure to properly satisfy the necessary positive second time derivative constraint when a principal Reynolds stress vanishes - a fatal flaw that becomes apparent when the non-analytic terms in their model are made single-valued as required on physical grounds. It is furthermore shown that the centrifugal acceleration generated by rotations of the principal axes of the Reynolds stress tensor can make the second derivative singular at the most extreme limits of realizable turbulence. This previously overlooked effect appears to make it impossible to identically satisfy the strong form of realizability in any version of the present generation of second-order closures. On the other hand, models properly formulated to satisfy the weak form of realizability - wherein states of one or two component turbulence are not accessible in finite time are found to be realizable. However, unlike the simpler and more commonly used second order closures, these models can be ill-behaved near the extreme limits of realizable turbulence due to the way that higher-degree nonlinearities are often unnecessarily introduced to satisfy realizability. Illustrative computations of homogeneous shear flows are presented to demonstrate these points which can have important implications for turbulence modeling.

Speziale, Charles G.↗
NASA NTRS1973

Effects of fineness and closure ratios on boattail drag of circular-arc-afterbody models with jet exhaust at Mach numbers up to 1.3

The effects of variations in fineness ratio and closure ratio on the boattail drag of circular-arc afterbody models were investigated at static conditions and at Mach numbers from 0.40 to 1.30. Jet total-pressure ratio was varied from jet off to about 6, depending on Mach number. Reynolds number based on model maximum diameter varied from about 1.2 million to 2.18 million. Angle of attack was varied from -4 deg to 8 deg. The results indicate that, at subsonic speeds before the drag rise and at a scheduled pressure ratio, the total drag (pressure plus calculated skin friction) is approximately equal for configurations with the same closure ratio. Also, for configurations with the same closure ratio, increasing fineness ratio results in an increased drag-rise Mach number; and for configurations with the same fineness ratio, increasing closure ratio results in an increased drag-rise Mach number.

Reubush, D. E.↗
NASA NTRS2007

Simulation of Boundary-Layer Cumulus and Stratocumulus Clouds using a Cloud-Resolving Model With Low- and Third-Order Turbulence Closures

The effects of subgrid-scale condensation and transport become more important as the grid spacings increase from those typically used in large-eddy simulation (LES) to those typically used in cloud-resolving models (CRMs). Incorporation of these effects can be achieved by a joint probability density function approach that utilizes higher-order moments of thermodynamic and dynamic variables. This study examines how well shallow cumulus and stratocumulus clouds are simulated by two versions of a CRM that is implemented with low-order and third-order turbulence closures (LOC and TOC) when a typical CRM horizontal resolution is used and what roles the subgrid-scale and resolved-scale processes play as the horizontal grid spacing of the CRM becomes finer. Cumulus clouds were mostly produced through subgrid-scale transport processes while stratocumulus clouds were produced through both subgrid-scale and resolved-scale processes in the TOC version of the CRM when a typical CRM grid spacing is used. The LOC version of the CRM relied upon resolved-scale circulations to produce both cumulus and stratocumulus clouds, due to small subgrid-scale transports. The mean profiles of thermodynamic variables, cloud fraction and liquid water content exhibit significant differences between the two versions of the CRM, with the TOC results agreeing better with the LES than the LOC results. The characteristics, temporal evolution and mean profiles of shallow cumulus and stratocumulus clouds are weakly dependent upon the horizontal grid spacing used in the TOC CRM. However, the ratio of the subgrid-scale to resolved-scale fluxes becomes smaller as the horizontal grid spacing decreases. The subcloud-layer fluxes are mostly due to the resolved scales when a grid spacing less than or equal to 1 km is used. The overall results of the TOC simulations suggest that a 1-km grid spacing is a good choice for CRM simulation of shallow cumulus and stratocumulus.

Xu, Kuan-Man↗
DOE OSTI2019

Updated Groundwater Flow Simulations of the Savannah River Site General Separations Area

The groundwater flow model supporting Performance Assessments (PAs) and Composite Analyses (CAs) at the Savannah River Site was significantly revised in 2016 and 2017 using new hydrostratigraphic surfaces, updated well water level calibration targets, and semi-automated model calibration with the PEST optimization code. This model is referred to as “GSA_2016”. This report documents further refinement of the GSA_2016 model in 2018 to incorporate updates to model calibration targets, closure of the H-Area Ash Basin, construction of E-Area Slit Trench operational covers, and plume information from the Mixed Waste Management Facility and Low-Level Radioactive Waste Disposal Facility. Another objective was to lower hydraulic head residuals by adding another calibration zone. The resulting model is referred to as “GSA_2018”.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1983

The mechanics of delamination in fiber-reinforced composite materials. Part 1: Stress singularities and solution structure

The fundamental mechanics of delamination in fiber composite laminates is studied. Mathematical formulation of the problem is based on laminate anisotropic elasticity theory and interlaminar fracture mechanics concepts. Stress singularities and complete solution structures associated with general composite delaminations are determined. For a fully open delamination with traction-free surfaces, oscillatory stress singularities always appear, leading to physically inadmissible field solutions. A refined model is introduced by considering a partially closed delamination with crack surfaces in finite-length contact. Stress singularities associated with a partially closed delamination having frictional crack-surface contact are determined, and are found to be diferent from the inverse square-root one of the frictionless-contact case. In the case of a delamination with very small area of crack closure, a simplified model having a square-root stress singularity is employed by taking the limit of the partially closed delamination. The possible presence of logarithmic-type stress singularity is examined; no logarithmic singularity of any kind is found in the composite delamination problem. Numerical examples of dominant stress singularities are shown for delaminations having crack-tip closure with different frictional coefficients between general (1) and (2) graphite-epoxy composites.

Wang, S. S.↗
NASA NTRS1983

The mechanics of delamination in fiber-reinforced composite materials. I - Stress singularities and solution structure

The fundamental mechanics of delamination in fiber composite laminates is studied. Mathematical formulation of the problem is based on laminate anisotropic elasticity theory and interlaminar fracture mechanics concepts. Stress singularities and complete solution structures associated with general composite delaminations are determined. For a fully open delamination with traction-free surfaces, oscillatory stress singularities always appear, leading to physically inadmissible field solutions. A refined model is introduced by considering a partially closed delamination with crack surfaces in finite-length contact. Stress singularities associated with a partially closed delamination having frictional crack-surface contact are determined, and are found to be different from the inverse square-root one of the frictionless-contact case. In the case of a delamination with very small area of crack closure, a simplified model having a square-root stress singularity is employed by taking the limit of the partially closed delamination. The possible presence of logarithmic-type stress singularity is examined; no logarithmic singularity of any kind is found in the composite delamination problem. Numerical examples of dominant stress singularities are shown for delaminations having crack-tip closure with different frictional coefficients between general (1) and (2) graphite-epoxy composites. Previously announced in STAR as N84-13221

Wang, S. S.↗
NASA NTRS1985

Simulation and modeling of homogeneous, compressed turbulence

Low Reynolds number homogeneous turbulence undergoing low Mach number isotropic and one-dimensional compression was simulated by numerically solving the Navier-Stokes equations. The numerical simulations were performed on a CYBER 205 computer using a 64 x 64 x 64 mesh. A spectral method was used for spatial differencing and the second-order Runge-Kutta method for time advancement. A variety of statistical information was extracted from the computed flow fields. These include three-dimensional energy and dissipation spectra, two-point velocity correlations, one-dimensional energy spectra, turbulent kinetic energy and its dissipation rate, integral length scales, Taylor microscales, and Kolmogorov length scale. Results from the simulated flow fields were used to test one-point closure, two-equation models. A new one-point-closure, three-equation turbulence model which accounts for the effect of compression is proposed. The new model accurately calculates four types of flows (isotropic decay, isotropic compression, one-dimensional compression, and axisymmetric expansion flows) for a wide range of strain rates.

Wu, C. T.↗
NASA NTRS1992

The remarkable ability of turbulence model equations to describe transition

This paper demonstrates how well the k-omega turbulence model describes the nonlinear growth of flow instabilities from laminar flow into the turbulent flow regime. Viscous modifications are proposed for the k-omega model that yield close agreement with measurements and with Direct Numerical Simulation results for channel and pipe flow. These modifications permit prediction of subtle sublayer details such as maximum dissipation at the surface, k approximately y(exp 2) as y approaches 0, and the sharp peak value of k near the surface. With two transition specific closure coefficients, the model equations accurately predict transition for an incompressible flat-plate boundary layer. The analysis also shows why the k-epsilon model is so difficult to use for predicting transition.

Wilcox, David C.↗
DOE OSTI2017

Spatiotemporal Variability of Turbulence Kinetic Energy Budgets in the Convective Boundary Layer over Both Simple and Complex Terrain

The assumption of sub-grid scale (SGS) horizontal homogeneity within a model grid cell, which forms the basis of SGS turbulence closures used by mesoscale models, becomes increasingly tenuous as grid spacing is reduced to a few kilometers or less, such as in many emerging high-resolution applications. Herein, we use the turbulence kinetic energy (TKE) budget equation to study the spatio-temporal variability in two types of terrain—complex (Columbia Basin Wind Energy Study [CBWES] site, north-eastern Oregon) and flat (ScaledWind Farm Technologies [SWiFT] site, west Texas) using the Weather Research and Forecasting (WRF) model. In each case six-nested domains (three domains each for mesoscale and large-eddy simulation [LES]) are used to downscale the horizontal grid spacing from 10 km to 10 m using the WRF model framework. The model output was used to calculate the values of the TKE budget terms in vertical and horizontal planes as well as the averages of grid cells contained in the four quadrants (a quarter area) of the LES domain. The budget terms calculated along the planes and the mean profile of budget terms show larger spatial variability at CBWES site than at the SWiFT site. The contribution of the horizontal derivative of the shear production term to the total production shear was found to be 45% and 15% of the total shear, at the CBWES and SWiFT sites, respectively, indicating that the horizontal derivatives applied in the budget equation should not be ignored in mesoscale model parameterizations, especially for cases with complex terrain with <10 km scale.

convective boundary layer↗
NASA NTRS1983

A nonlinear fracture mechanics approach to the growth of small cracks

An analytical model of crack closure is used to study the crack growth and closure behavior of small cracks in plates and at notches. The calculated crack opening stresses for small and large cracks, together with elastic and elastic plastic fracture mechanics analyses, are used to correlate crack growth rate data. At equivalent elastic stress intensity factor levels, calculations predict that small cracks in plates and at notches should grow faster than large cracks because the applied stress needed to open a small crack is less than that needed to open a large crack. These predictions agree with observed trends in test data. The calculations from the model also imply that many of the stress intensity factor thresholds that are developed in tests with large cracks and with load reduction schemes do not apply to the growth of small cracks. The current calculations are based upon continuum mechanics principles and, thus, some crack size and grain structure exist where the underlying fracture mechanics assumptions become invalid because of material inhomogeneity (grains, inclusions, etc.). Admittedly, much more effort is needed to develop the mechanics of a noncontinuum. Nevertheless, these results indicate the importance of crack closure in predicting the growth of small cracks from large crack data.

Newman, J. C., Jr.↗
NASA NTRS1982

Computer program for solving laminar, transitional, or turbulent compressible boundary-layer equations for two-dimensional and axisymmetric flow

A numerical algorithm and computer program are presented for solving the laminar, transitional, or turbulent two dimensional or axisymmetric compressible boundary-layer equations for perfect-gas flows. The governing equations are solved by an iterative three-point implicit finite-difference procedure. The software, program VGBLP, is a modification of the approach presented in NASA TR R-368 and NASA TM X-2458, respectively. The major modifications are: (1) replacement of the fourth-order Runge-Kutta integration technique with a finite-difference procedure for numerically solving the equations required to initiate the parabolic marching procedure; (2) introduction of the Blottner variable-grid scheme; (3) implementation of an iteration scheme allowing the coupled system of equations to be converged to a specified accuracy level; and (4) inclusion of an iteration scheme for variable-entropy calculations. These modifications to the approach presented in NASA TR R-368 and NASA TM X-2458 yield a software package with high computational efficiency and flexibility. Turbulence-closure options include either two-layer eddy-viscosity or mixing-length models. Eddy conductivity is modeled as a function of eddy viscosity through a static turbulent Prandtl number formulation. Several options are provided for specifying the static turbulent Prandtl number. The transitional boundary layer is treated through a streamwise intermittency function which modifies the turbulence-closure model. This model is based on the probability distribution of turbulent spots and ranges from zero to unity for laminar and turbulent flow, respectively. Several test cases are presented as guides for potential users of the software.

Harris, J. E.↗
DOE OSTI2012

Open-source MFIX-DEM software for gas–solids flows: Part I—Verification studies

With rapid advancements in computer hardware, it is now possible to perform large simulations of granular flows using the Discrete Element Method (DEM). As a result, solids are increasingly treated in a discrete Lagrangian fashion in the gas–solids flow community. In this paper, the open-source MFIX-DEM software is described that can be used for simulating the gas–solids flow using an Eulerian reference frame for the continuum fluid and a Lagrangian discrete framework (Discrete Element Method) for the particles. This method is referred to as the continuum discrete method (CDM) to clearly make a distinction between the ambiguity of using a Lagrangian or Eulerian reference for either continuum or discrete formulations. This freely available CDM code for gas–solids flows can accelerate the research in computational gas–solids flows and establish a baseline that can lead to better closures for the continuum modeling (or traditionally referred to as two fluid model) of gas–solids flows. In this paper, a series of verification cases is employed which tests the different aspects of the code in a systematic fashion by exploring specific physics in gas–solids flows before exercising the fully coupled solution on simple canonical problems. It is critical to have an extensively verified code as the physics is complex with highly-nonlinear coupling, and it is difficult to ascertain the accuracy of the results without rigorous verification. These series of verification tests set the stage not only for rigorous validation studies (performed in part II of this paper) but also serve as a procedure for testing any new developments that couple continuum and discrete formulations for gas–solids flows.

Computational fluid dynamics↗
DOE OSTI2012

Center for Gyrokinetic/MHD Hybrid Simulation of Energetic Particle Physics in Toroidal Plasmas (Final Report)

At Colorado University-Boulder the primary task is to extend our gyrokinetic Particle-in-Cell simulation of tokamak micro-turbulence and transport to the area of energetic particle physics. We have implemented a gyrokinetic ion/massless fluid electron hybrid model in the global δƒ-PIC code GEM, and benchmarked the code with analytic results on the thermal ion radiative damping rate of Toroidal Alfvén Eigenmodes (TAE) and with mode frequency and spatial structure from eigenmode analysis. We also performed nonlinear simulations of both a single-n mode (n is the toroidal mode number) and multiple-n modes, and in the case of single-n, benchmarked the code on the saturation amplitude vs. particle collision rate with analytical theory. Most simulations use the f method for both ions species, but we have explored the full-ƒ method for energetic particles in cases where the burst amplitude of the excited instabilities is large as to cause significant re-distribution or loss of the energetic particles. We used the hybrid model to study the stability of high-n TAEs in ITER. Our simulations show that the most unstable modes in ITER lie in the rage of 10 < n < 20. Thermal ion pressure effect and alpha particles non-perturbative effect are important in determining the mode radial location and stability threshold. The thermal ion Landau damping rate and radiative damping rate from the simulations are compared with analytical estimates. The thermal ion Landau damping is the dominant damping mechanism. Plasma elongation has a strong stabilizing effect on the alpha driven TAEs. The central alpha particle pressure threshold for the most unstable n=15 mode is about β α (0) = 0.7% for the fully shaped ITER equilibrium. We also carried nonlinear simulations of the most unstable n = 15 mode and found that the saturation amplitude for the nominal ITER discharge is too low to cause large redistribution or loss of alpha particles. To include kinetic electron effects in the hybrid model we have studied a kinetic electron closure scheme for the fluid electron model. The most important element of the closure scheme is a complete Ohm’s law for the parallel electric field E || , derived by combining the quasi-neutrality condition, the Ampere’s equation and the v || moment of the gyrokinetic equations. A discretization method for the closure scheme is studied in detail for a three-dimensional shear-less slab plasma. It is found that for long-wavelength shear Alfvén waves the kinetic closure scheme is both more accurate and robust than the previous GEM algorithm using the split-scheme, whereas for the ion-gradient-driven instability the previous algorithm is more efficient. This kinetic electron closure scheme will be implemented in GEM in the future. We have studied the beam driven Reverse Shear Alfven Eigenmodes (RSAE) observed in DIII-D discharge 142111. For this purpose a new scheme for obtaining the electric potential is implemented, i.e., by solving the gyrokinetic moment (GKM) equation, which is essentially the equation for $∂\phi/∂t$ used in GEM’s split-weight scheme, and then integrating in time. Due to charge-neutrality the ExB motions of the equilibrium densities of all species cancel each other and do not cause charge separation if there is no finite Larmor radius effect. The advantage of solving the GKM equation is that this lowest-order cancellation can be made explicit. The GKM approach is found to be more accurate and robust. GEM simulations have reproduced many features of RSAE seen in the experiment, such as frequency chirping and the chirping range. It has been reported by other simulation codes that the shearing direction of the mode structure in the poloidal plane disagrees with observation. We found that the mode structure, including the shearing in the poloidal plane, is in general sensitive to the beam distribution. Using the same beam density profile as in other codes but with a slowing-down distribution in velocity, GEM simulations reproduce the mode shearing direction seen in the experiment. We have carried out extensive nonlinear simulations of RSAE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Opening and closing of cracks at high cyclic strains

The closure behavior of cracks of different length and at different cyclic strain levels (ranging from predominantly elastic to grossly plastic strains) was studied to observe the effect of residual crack-tip plasticity on crack closure. Cracks were initiated either naturally or artificially (from electric discharge machining pits) in uniaxial test specimens of strengthened alloy steel AISI 4340 with a grain size of 0.016 mm. It was found that, at high strains, cracks closed only when the lowest stress level in the cycle was approached. The stress or the strain opening level depended upon the exact point along the crack length where the observations were made. As the plastic deformation increased, the relative crack opening level was found to decrease and approach the value of stress ratio R. The experimental results were compared with those of three analytical models of crack closure and opening, demonstrating the limitations of the currently available elastic-plastic crack growth analysis.

Iyyer, N. S.↗
DOE OSTI2018

Active learning of constitutive relation from mesoscopic dynamics for macroscopic modeling of non-Newtonian flows

We simulate complex fluids by means of an on-the-fly coupling of the bulk rheology to the underlying microstructure dynamics. In particular, a continuum model of polymeric fluids is constructed without a pre-specified constitutive relation, but instead it is actively learned from mesoscopic simulations where the dynamics of polymer chains is explicitly computed. To couple the bulk rheology of polymeric fluids and the microscale dynamics of polymer chains, the continuum approach (based on the finite volume method) provides the transient flow field as inputs for the (mesoscopic) dissipative particle dynamics (DPD), and in turn DPD returns an effective constitutive relation to close the continuum equations. In this multiscale modeling procedure, we employ an active learning strategy based on Gaussian process regression (GPR) to minimize the number of expensive DPD simulations, where adaptively selected DPD simulations are performed only as necessary. Numerical experiments are carried out for flow past a circular cylinder of a non-Newtonian fluid, modeled at the mesoscopic level by bead-spring chains. The results show that only five DPD simulations are required to achieve an effective closure of the continuum equations at Reynolds number RE = 10. Furthermore, when Re is increased to 100, only one additional DPD simulation is required for constructing an extended GPR-informed model closure. Compared to traditional message-passing multiscale approaches, applying an active learning scheme to multiscale modeling of non-Newtonian fluids can significantly increase the computational efficiency. Although the method demonstrated here obtains only a local viscosity from the polymer dynamics, it can be extended to other multiscale models of complex fluids whose macro-rheology is unknown.

Computer Science↗
NASA NTRS1993

Flow interaction experiment. Volume 1: Aerothermal modeling, phase 2

An experimental and computational study is reported for the flow of a turbulent jet discharging into a rectangular enclosure. The experimental configurations consisting of primary jets only, annular jets only, and a combination of annular and primary jets are investigated to provide a better understanding of the flow field in an annular combustor. A laser Doppler velocimeter is used to measure mean velocity and Reynolds stress components. Major features of the flow field include recirculation, primary and annular jet interaction, and high turbulence. A significant result from this study is the effect the primary jets have on the flow field. The primary jets are seen to create statistically larger recirculation zones and higher turbulence levels. In addition, a technique called marker nephelometry is used to provide mean concentration values in the model combustor. Computations are performed using three levels of turbulence closures, namely k-epsilon model, algebraic second moment (ASM), and differential second moment (DSM) closure. Two different numerical schemes are applied. One is the lower-order power-law differencing scheme (PLDS) and the other is the higher-order flux-spline differencing scheme (FSDS). A comparison is made of the performance of these schemes. The numerical results are compared with experimental data. For the cases considered in this study, the FSDS is more accurate than the PLDS. For a prescribed accuracy, the flux-spline scheme requires a far fewer number of grid points. Thus, it has the potential for providing a numerical error-free solution, especially for three-dimensional flows, without requiring an excessively fine grid. Although qualitatively good comparison with data was obtained, the deficiencies regarding the modeled dissipation rate (epsilon) equation, pressure-strain correlation model, and the inlet epsilon profile and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗