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 1,063 records · Page 59

Transport methods and interactions for space radiations

This report presents a brief history leading to the involvement of the Langley Research Center of the National Aeronautics and Space Administration (NASA) in space-radiation physics and protection. Indeed, a relatively complete summary of technical capability as of the summer of 1990 is given. The Boltzmann equations for coupled ionic and neutronic fields are presented and inversion techniques for the Boltzmann operator are discussed. Errors generated by the straight ahead approximation are derived and are shown to be negligible for most problems of space-radiation protection. A decoupling of projectile propagation from the target fields greatly simplifies the Boltzmann equations and allows an analytic solution of the target fragment transport. Analytic and numerical methods of solving the projectile transport equations are discussed. The nuclear physics underlying the coefficients in the Boltzmann equation is discussed. A coupled-channel optical model is found as a consequence of the loose binding of nuclear matter and closure of the nuclear states in high-energy reactions. Transport solutions with the developed data base are used with laboratory experiments to validate both the transport code and the data base. Numerical benchmarks and comparison with Monte Carlo calculations are also used for code validation.

Wilson, John W.↗

A Hybrid Numerical Method for Turbulent Mixing Layers

A hybrid method has been developed for simulations of compressible turbulent mixing layers. Such mixing layers dominate the flows in exhaust systems of modern day aircraft and also those of hypersonic vehicles currently under development. The method configurations in which a dominant structural feature provides an unsteady mechanism to drive the turbulent development in the mixing layer. The hybrid method uses a Reynolds-averaged Navier-Stokes (RANS) procedure to calculate wall bounded regions entering a mixing section, and a Large Eddy Simulation (LES) procedure to calculate the mixing dominated regions. A numerical technique was developed to enable the use of the hybrid RANS-LES method on stretched, non-Cartesian grids. Closure for the RANS equations was obtained using the Cebeci-Smith algebraic turbulence model in conjunction with the wall-function approach of Ota and Goldberg. The wall-function approach enabled a continuous computational grid from the RANS regions to the LES region. The LES equations were closed using the Smagorinsky subgrid scale model. The hybrid RANS-LES method is applied to a benchmark compressible mixing layer experiment. Preliminary two dimensional calculations are used to investigate the effects of axial grid density and boundary conditions. Vortex shedding from the base region of a splitter plate separating the upstream flows was observed to eventually transition to turbulence. The location of the transition, however, was much further downstream than indicated by experiments. Actual LES calculations, performed in three spatial directions, also indicated vortex shedding, but the transition to turbulence was found to occur much closer to the beginning of the mixing section. which is in agreement with experimental observations. These calculations demonstrated that LES simulations must be performed in three dimensions. Comparisons of time-averaged axial velocities and turbulence intensities indicated reasonable agreement with experimental data.

Georgiadis, Nicholas J.↗

Toward a Unified Parameterization of the Boundary Layer and Moist Convection: Lateral Mass Exchanges and Subplume-Scale Fluxes - Part 2

The dissipation parameterizations developed for higher-order closure are used to parameterize lateral entrainment and detrainment in a mass-flux model. In addition, a subplume-scale turbulence scheme is included to represent fluxes not captured in the conventional mass-flux framework. These new parameterizations are tested by simulating trade wind cumulus from the Barbados Oceanographic and Meteorological Experiment (BOMEX).

Lappan, Cara-Lyn↗

Application of Benchmark Examples to Assess the Single and Mixed-Mode Static Delamination Propagation Capabilities in ANSYS

The application of benchmark examples for the assessment of quasi-static delamination propagation capabilities is demonstrated for ANSYS. The examples are independent of the analysis software used and allow the assessment of the automated delamination propagation in commercial finite element codes based on the virtual crack closure technique (VCCT). The examples selected are based on two-dimensional finite element models of Double Cantilever Beam (DCB), End-Notched Flexure (ENF), Mixed-Mode Bending (MMB) and Single Leg Bending (SLB) specimens. First, the quasi-static benchmark examples were recreated for each specimen using the current implementation of VCCT in ANSYS . Second, the delamination was allowed to propagate under quasi-static loading from its initial location using the automated procedure implemented in the finite element software. Third, the load-displacement relationship from a propagation analysis and the benchmark results were compared, and good agreement could be achieved by selecting the appropriate input parameters. The benchmarking procedure proved valuable by highlighting the issues associated with choosing the input parameters of the particular implementation. Overall the results are encouraging, but further assessment for three-dimensional solid models is required.

Krueger, Ronald↗

Automated Symbolic Upscaling: 1. Model Generation for Extended Applicability Regimes

Abstract In porous media theory, upscaling techniques are fundamental to deriving rigorous Darcy‐scale models for flow and reactive transport in subsurface systems. Due to limitations in classical techniques, a number of ad hoc approaches have been proposed to address physical regimes in which reactive time scales are similar to, or faster than, diffusive time scales. In Part 1 of this two part series, we present a strategy for expanding the applicability of classical homogenization theory by generalizing the assumed closure form. We detail the implementation of this strategy on two reactive mass transport problems with moderately reactive physics. The strategy produces nontrivial homogenized models with emergent terms and effective parameters that couple reactive, diffusive, and advective transport. The differences in equation forms between the macroscopic and pore‐scale descriptions advise caution to further studies where the forms of macroscopic equations are assumed, as opposed to rigorously derived. Numerical validation is provided for each problem to show that the error estimates of homogenization theory are satisfied, and to justify the implemented strategy. In Part 2, the presented strategy is automated using symbolic computing to expedite its implementation.

Pietrzyk, Kyle↗

An Enhanced Virtual Crack Closure Technique for Stress Intensity Factor Calculation along Arbitrary Crack Fronts and the Application in Hydraulic Fracturing Simulation

The virtual crack closure technique (VCCT) is widely used for calculating energy release rates along crack fronts and modeling the propagation of cracks in solid materials. Although the VCCT formulation for smooth crack fronts has been sufficiently addressed in the literature, the application of VCCT to a nonsmoothed crack front with sharp corners warrants further investigation. In this study, we describe an enhanced VCCT to calculate energy release rates and stress intensity factors for cracks with arbitrary shapes in 3D domains discretized on structured grids. The formulations of the enhanced VCCT were developed and implemented into a multiphysics simulation environment capable of simulating crack propagation in the framework of linear elastic fracture mechanics. Comparisons with existing analytical/numerical solutions and other VCCT approaches were performed to verify the enhanced VCCT in terms of SIF calculation along nonsmoothed crack fronts. We then applied the enhanced VCCT to a hydraulically driven penny-shaped fracture problem to further demonstrate its capability to simulate nonsmoothed fracture propagation.

02 PETROLEUM↗

Safety Functions and Features, Events and Processes for the E-Area Performance Assessment

The DOE Technical Standard, “Disposal Authorization Statement and Tank Closure Documentation,” (DOE 2017) recommends the use of safety functions and features, events and processes (FEPs) to support development of conceptual models and identification of scenarios to be considered in a performance assessment (PA). The FEP process provides a means to describe how a PA considers and addresses the factors that could influence the performance of key barriers. Understanding the roles of barriers in terms of limiting migration helps to focus on how changes in the system could lead to a situation where those roles cannot be fulfilled and there is the potential for compromised performance. The FEPs screening and review process was used to identify FEPs that are relevant for the EArea Low-Level Waste Facility (LLWF) and specifically those FEPs that could have a detrimental impact on the effectiveness of a given safety function. For this PA, a default list of FEPs developed at the International Atomic Energy Agency (IAEA 2004) and an approach implemented for PAs at the Hanford and Idaho sites (Mehta et al. 2016, DOE-ID 2019) are used to identify processes and events that could influence the effectiveness of a given safety function for the E-Area LLWF (e.g., subsidence can impact the safety function of the cover system and lead to increased infiltration). The Hanford and Idaho PAs represent two of the most recent applications of this approach. The PA evaluates the potential impacts of changes in performance of different features of the system and demonstrates that the safety functions represent multiple and redundant barriers. Barrier analyses, assuming a safety function is not present, also test the robustness of the system in the event of the loss of one or more safety functions. Such evaluations also support a qualitative illustration of the concept of defense in depth. The safety concept for closure of the E-Area LLWF (generically referred to as “E-Area”) encompasses a variety of different features (i.e., administrative controls, natural site features, and engineered barriers) that reduce the potential impacts on human health and the environment from the residual waste that will remain after closure. These features can be represented as a collection of safety functions acting independently and as a system to provide for overall safety. In some applications, there have been attempts to assign numerical expectations to specific safety functions, but that is not the intent in this case. The concept of safety functions is used more qualitatively in two ways for this PA: 1. To illustrate the robustness of the E-Area design, operational practices and closure approach by documenting features that are and are not credited in different modeling cases. 2. To identify the roles of the different features and potential processes and events that could compromise the performance of safety features and need to be considered when developing the modeling approach. This report addresses both safety functions and FEPs for the E-Area PA.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

On the modelling of scalar and mass transport in combustor flows

Results are presented of a numerical study of swirling and nonswirling combustor flows with and without density variations. Constant-density arguments are used to justify closure assumptions invoked for the transport equations for turbulent momentum and scalar fluxes, which are written in terms of density-weighted variables. Comparisons are carried out with measurements obtained from three different axisymmetric model combustor experiments covering recirculating flow, swirling flow, and variable-density swirling flow inside the model combustors. Results show that the Reynolds stress/flux models do a credible job of predicting constant-density swirling and nonswirling combustor flows with passive scalar transport. However, their improvements over algebraic stress/flux models are marginal. The extension of the constant-density models to variable-density flow calculations shows that the models are equally valid for such flows.

Nikjooy, M.↗

Modeling and Test Data Analysis of a Tank Rapid Chill and Fill System for the Advanced Shuttle Upper Stage (ASUS) Concept

The Advanced Shuttle Upper Stage (ASUS) concept addresses safety concerns associated .with cryogenic stages by launching empty, and filling on ascent. The ASUS employs a rapid chill and fill concept. A spray bar is used to completely chill the tank before fill, allowing the vent valve to be closed during the fill process. The first tests of this concept, using a flight size (not flight weight) tank. were conducted at Marshall Space Flight Center (MSFC) during the summer of 2000. The objectives of the testing were to: 1) demonstrate that a flight size tank could be filled in roughly 5 minutes to accommodate the shuttle ascent window, and 2) demonstrate a no-vent fill of the tank. A total of 12 tests were conducted. Models of the test facility fill and vent systems, as well as the tank, were constructed. The objective of achieving tank fill in 5 minutes was met during the test series. However, liquid began to accumulate in the tank before it was chilled. Since the tank was not chilled until the end of each test, vent valve closure during fill was not possible. Even though the chill and fill process did not occur as expected, reasonable model correlation with the test data was achieved.

Flachbart, Robin↗

Charge radii in covariant density functional theory: A global view

A systematic global investigation of differential charge radii has been performed within the CDFT framework for the first time. Theoretical results obtained with conventional covariant energy density functionals and the separable pairing interaction of Tian et al. [Phys. Lett. B 676, 44 (2009)] are compared with experimental differential charge radii in the regions of the nuclear chart in which available experimental data crosses the neutron shell closures at N=28,50,82, and 126. The analysis of absolute differential radii of different isotopic chains and their relative properties indicate clearly that such properties are reasonably well described in model calculations in the cases when the mean-field approximation is justified. However, while the observed clusterization of differential charge radii of different isotopic chains is well described above the N=50 and N=126 shell closures, it is more difficult to reproduce it above the N=28 and N=82 shell closures because of possible deficiencies in the underlying single-particle structure. The impact of the latter has been evaluated for spherical shapes and it was shown that the relative energies of the single-particle states and the patterns of their occupation with increasing neutron number have an appreciable impact on the evolution of the δ$\langle$r2$\rangle$N,N' values. These factors also limit the predictive power of model calculations in the regions of high densities of the single-particle states of different origin. It is shown that the kinks in the charge radii at neutron shell closures are due to the underlying single-particle structure and due to weakening or collapse of pairing at these closures. The regions of the nuclear chart in which the correlations beyond mean field are expected to have an impact on charge radii are indicated; the analysis shows that the assignment of a calculated excited prolate minimum to the experimental ground state allows us to understand the trends of the evolution of differential charge radii with neutron number in many cases of shape coexistence even at the mean-field level. It is usually assumed that pairing is a dominant contributor to odd-even staggering (OES) in charge radii. Our analysis paints a more complicated picture. It suggests a new mechanism in which the fragmentation of the single-particle content of the ground state in odd-mass nuclei due to particle-vibration coupling provides a significant contribution to OES in charge radii.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Panel Stiffener Debonding Analysis using a Shell/3D Modeling Technique

A shear loaded, stringer reinforced composite panel is analyzed to evaluate the fidelity of computational fracture mechanics analyses of complex structures. Shear loading causes the panel to buckle. The resulting out -of-plane deformations initiate skin/stringer separation at the location of an embedded defect. The panel and surrounding load fixture were modeled with shell elements. A small section of the stringer foot, web and noodle as well as the panel skin near the delamination front were modeled with a local 3D solid model. Across the width of the stringer fo to, the mixed-mode strain energy release rates were calculated using the virtual crack closure technique. A failure index was calculated by correlating the results with a mixed-mode failure criterion of the graphite/epoxy material. The objective was to study the effect of the fidelity of the local 3D finite element model on the computed mixed-mode strain energy release rates and the failure index.

Krueger, Ronald↗

Small Scale Response and Modeling of Periodically Forced Turbulence

The response of the small scales of isotropic turbulence to periodic large scale forcing is studied using two-point closures. The frequency response of the turbulent kinetic energy and dissipation rate, and the phase shifts between production, energy and dissipation are determined as functions of Reynolds number. It is observed that the amplitude and phase of the dissipation exhibit nontrivial frequency and Reynolds number dependence that reveals a filtering effect of the energy cascade. Perturbation analysis is applied to understand this behavior which is shown to depend on distant interactions between widely separated scales of motion. Finally, the extent to which finite dimensional models (standard two-equation models and various generalizations) can reproduce the observed behavior is discussed.

Bos, Wouter↗

Assessing CLUBB PDF Closure Assumptions for a Continental Shallow-to-Deep Convective Transition Case Over Multiple Spatial Scales

Assumed-PDF (probability density function) higher-order turbulence closures (APHOCs) are now widely used for parameterizing boundary layer turbulence and shallow convection in Earth system models (ESMs). A better understanding of the resolution-dependent behavior of APHOCs is essential for improving the performance of next-generation ESMs with intended horizontal resolutions finer than 10 km. In this study, we evaluate the PDF family of Analytic double-Gaussian 1 implemented in Cloud Layers Unified By Binormals (CLUBB) over a range of spatial scales (Dx) from 2 to 100 km. A 120-km-wide large eddy simulation (LES) for a continental convection case during 2016 Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign serves as benchmark to evaluate the PDF closure using an off-line approach. We find during the shallow convection period, the CLUBB PDF closure tends to produce positive biases of cloud properties and liquid water flux near cloud base for all scales of analysis. It produces negative biases for these variables near cloud top that are more severe for Dx larger than 25 km. Results show that replacing the CLUBB-parameterized moisture and temperature skewnesses with LES-derived ones can fix most of the biases if clipping of input moments is allowed to prevent the occurrence of unrealizable solutions. Overall, the performance of the PDF closure is better for smaller Dx = 2–5 km than for larger Dx = 50–100 km; for a given grid spacing, it is better when the convective clouds become deeper in the late afternoon. Likely causes for the resolution dependence and implications for improving the PDF closure are discussed.

54 ENVIRONMENTAL SCIENCES↗

Effects of pressure gradients on turbulent premixed flames

The influence of a constant acceleration on a turbulent premixed flame is studied by direct numerical simulation. This acceleration induces a mean pressure gradient across the flame brush, leading to a modification of the turbulent flame structure due to differential buoyancy mechanisms between heavy cold fresh and light hot burnt gases. Such a pressure gradient may be encountered in practical applications in ducted flames. A favorable pressure gradient, i.e. the pressure decreases from unburnt to burnt gases, is found to decrease the flame wrinkling, the flame brush thickness, and the turbulent flame speed. A favorable pressure gradient also promotes counter-gradient turbulent transport. On the other hand, adverse pressure gradients tend to increase the flame brush thickness and turbulent flame speed, and promote classical gradient turbulent transport. The balance equation for the turbulent flux of the Favre averaged progress variable is also analyzed. The first results show that the fluctuating pressure term, cannot be neglected as generally assumed in models. Simple models assuming that a high mean pressure gradient may only be balanced by the cross-dissipation term seem too approximate. This analysis has to be continued to compare simulation data and closure schemes proposed for the transport equation. The analysis developed by Veynante et al.(1995) has been extended to imposed acceleration and mean pressure gradients. A simple model for the turbulent flux is proposed and validated from simulation data. Then, a modified criterion is derived to delineate between counter-gradient and gradient turbulent diffusion. In fact, counter-gradient diffusion may occur in most practical applications, especially for ducted flames.

Veynante, D.↗

A Dual Spring Modeling Approach for Static and Fatigue Failure Assessments of Carbon/Epoxy Composite Sub-Elements

A dual spring model is developed for the static and fatigue damage predictions of laminates interface in composite structures. Stress concentrations can be induced by the defects formed in the fabrication or service process. A conventional S-N based fatigue damage model may not be accurate to predict the fatigue life of a structure with high stress concentration. With the dual spring model, static delamination failure can be simulated using springs of cohesive type material model while fatigue delamination development can be predicted using linear springs, where the crack driving force is computed based on virtual crack closure technique (VCCT). A Paris law type fatigue growth law with its mode mixity is applied for fatigue crack growth prediction. After verified using benchmark examples, including Double Cantilever Beam (DCB), End-Notched Flexure (ENF) and Mix-Mode Bending (MMB), the proposed dual spring model is applied in the static and fatigue damage prediction of NASA/Boeing sub-elements and UTC sub-elements.

Cui, Xiaodong↗

Surprising Charge-Radius Kink in the Sc Isotopes at N = 20

For this work, charge radii of neutron deficient 40 Sc and 41 Sc nuclei were determined using collinear laser spectroscopy. With the new data, the chain of Sc charge radii extends below the neutron magic number N=20 and shows a pronounced kink, generally taken as a signature of a shell closure, but one notably absent in the neighboring Ca, K, and Ar isotopic chains. Theoretical models that explain the trend at N=20 for the Ca isotopes cannot reproduce this puzzling behavior.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Scale Separation Effects on Simulations of Plasma Turbulence

Abstract Understanding plasma turbulence requires a synthesis of experiments, observations, theory, and simulations. In the case of kinetic plasmas such as the solar wind, the lack of collisions renders the fluid closures such as viscosity meaningless and one needs to resort to higher-order fluid models or kinetic models. Typically, the computational expense in such models is managed by simulating artificial values of certain parameters such as the ratio of the Alfvén speed to the speed of light ( v A / c ) or the relative mass ratio of ions and electrons ( m i / m e ). Although, typically care is taken to use values as close as possible to realistic values within the computational constraints, these artificial values could potentially introduce unphysical effects. These unphysical effects could be significant at sub-ion scales, where kinetic effects are the most important. In this paper, we use the 10-moment fluid model in the Gkeyll framework to perform controlled numerical experiments, systematically varying the ion–electron mass ratio from a small value down to the realistic proton–electron mass ratio. We show that the unphysical mass ratio has a significant effect on the kinetic range dynamics as well as the heating of both plasma species. The dissipative process for both ions and electrons becomes more compressive in nature, although the ions remain nearly incompressible in all cases. The electrons move from being dominated by incompressive viscous-like heating/dissipation to very compressive heating/dissipation dominated by compressions/rarefactions. While the heating change is significant for the electrons, a mass ratio of m i / m e ∼ 250 captures the asymptotic behavior of electron heating.

79 ASTRONOMY AND ASTROPHYSICS↗

A linear stochastic model of the human operator

A linear stochastic model of the human operator is developed and applied to the problem of piloted control of an aircraft. The pilot and aircraft are modeled as linear time-invariant systems containing both process and measurement noise. The loop closure by the pilot is determined by formulating the problem as an optimal stochastic control problem. The solution to the optimal control problem yields not only the pilot's optimal control output which he uses to control the vehicle, but also the optimal combination of his observations of the vehicle states upon which the pilot bases his control. A method is presented so that, using experimental pilot vehicle data, the cost functional which is minimized in the optimal control problem will be numerically equal to the pilot rating that the pilot would associate with the given vehicle and task.

Durrett, J. C.↗