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 793 records · Page 44

Navier-Stokes calculations on multi-element airfoils using a chimera-based solver

A study of Navier-Stokes calculations of flows about multielement airfoils using a chimera grid approach is presented. The chimera approach utilizes structured, overlapped grids which allow great flexibility of grid arrangement and simplifies grid generation. Calculations are made for two-, three-, and four-element airfoils, and modeling of the effect of gap distance between elements is demonstrated for a two element case. Solutions are obtained using the thin-layer form of the Reynolds averaged Navier-Stokes equations with turbulence closure provided by the Baldwin-Lomax algebraic model or the Baldwin-Barth one equation model. The Baldwin-Barth turbulence model is shown to provide better agreement with experimental data and to dramatically improve convergence rates for some cases. Recently developed, improved farfield boundary conditions are incorporated into the solver for greater efficiency. Computed results show good comparison with experimental data which include aerodynamic forces, surface pressures, and boundary layer velocity profiles.

Jasper, Donald W.↗

Synchrotron self-Compton in a radiative-adiabatic fireball scenario: modelling the multiwavelength observations in some Fermi /LAT bursts

Energetic GeV photons expected from the closest and the most energetic Gamma-ray bursts (GRBs) provide a unique opportunity to study the very-high-energy emission as well as the possible correlations with lower energy bands in realistic GRB afterglow models. In the standard GRB afterglow model, the relativistic homogeneous shock is usually considered to be fully adiabatic, however, it could be partially radiative. Based on the external forward-shock scenario in both stellar wind and constant-density medium, we present a radiative-adiabatic analytical model of the synchrotron self-Compton (SSC) and synchrotron processes considering an electron energy distribution with a power-law index of $1\lt p\lt 2$ and $2\le p$. We show that the SSC scenario plays a relevant role in the radiative parameter $\epsilon$, leading to a prolonged evolution during the slow cooling regime. In a particular case, we derive the Fermi/LAT light curves together with the photons with energies $\ge 100$ MeV in a sample of nine bursts from the second Fermi/LAT GRB catalogue that exhibited temporal and spectral indices with $\gtrsim 1.5$ and $\approx 2$, respectively. These events can hardly be described with closure relations of the standard synchrotron afterglow model, and also exhibit energetic photons above the synchrotron limit. We have modelled the multiwavelength observations of our sample to constrain the microphysical parameters, the circumburst density, the bulk Lorentz factor, and the mechanism responsible for explaining the energetic GeV photons.

acceleration of particles↗

Factors Governing Cloud Growth and Entrainment Rates in Shallow Cumulus and Cumulus Congestus During GoAmazon2014/5

Shallow cumulus and cumulus congestus clouds play an important role in the large-scale tropical circulation by mixing heat and moisture vertically and preconditioning the environment for deeper convection. Different representations of these shallow clouds account for much of the spread in General Circulation Model (GCMs) climate sensitivity, potentially because of how entrainment is represented in GCM parameterizations. This study uses observations from the Department of Energy's Atmospheric Radiation Measurement (ARM) mobile facility deployed at Manacapuru, Brazil, during the Green Ocean Amazon (GoAmazon2014/5) Campaign. Environmental thermodynamic profiles and observations of cloud top height (CTH) are used to constrain an entraining plume model to estimate bulk entrainment rates (ERs). Estimates of CTH are obtained from a combination of vertically pointing W-band ARM cloud radar and 1,290 MHz Radar Wind Profiler observations. A combination of radiosonde, microwave radiometer profiler, and microwave radiometer observations provides new best estimates of the environmental thermodynamic state. We quantify uncertainty in ERs considering uncertainties in estimated CTH, environmental thermodynamic properties, and assumed initial parcel characteristics. We find ERs ranging from 0.16 to 2.8 km -1 with an average of 0.58 ± 0.10 km -1 over a selected population of 469 shallow cumulus and cumulus congestus clouds. Using the retrieved estimates of ER, we evaluate several entrainment closures that are currently used in atmospheric models or have been proposed based on theory or large eddy simulation. Finally, entrainment rates in cumulus clouds are weakly correlated with low-level buoyancy, cloud depth, and cloud size.

54 ENVIRONMENTAL SCIENCES↗

Reynolds stress closures: Status and prospects

The basic pattern of Reynolds closures and the reason they look as an attractive type of model for practical shear flow calculations is summarized. The relationship between the organized structures and Reynolds stress closures is discussed. An outline given of fundamental developments that are being introduced to extend the modest reliability of Reynolds stress closures.

Launder, B. E.↗

Modeling deep convection in the Greenland Sea

The development of deep convective events in the high-latitude ocean is studied using a three-dimensional, coupled ice-ocean model. Oceanic mixing is described according to the level 2.5 turbulence closure scheme in which convection occurs in a continuous way, i.e., convective adjustment is not invoked. The model is forced by strong winds and surface cooling. Strong upwelling at the multilyear ice edge and consequent entrainment of warm Atlantic waters into the mixed layer is produced by winds parallel to the ice edge. Concomitant cooling drives deep convection and produces chimneylike structures. Inclusion of a barotropic mean flow over topography to the model provides important preconditioning and selects the location of deep convection. The most efficient preconditioning occurs at locations where the flow ascends a slope. In a stratified environment similar to the Greenland Sea with a 12 m/s wind the model simulations show that localized deep convection takes place after about 10 days to depths of 1000 m.

Hakkinen, S.↗

Precision mass measurement of lightweight self-conjugate nucleus 80 Zr

Protons and neutrons in the atomic nucleus move in shells analogous to the electronic shell structures of atoms. The nuclear shell structure varies due to changes of the nuclear mean field with the number of neutrons N and protons Z, and these variations can be probed by measuring mass differences between nuclei. The N=Z=40 self-conjugate nucleus 80 Zr is of particular interest as its proton and neutron shell structures are expected to be very similar, and its ground state is highly deformed. Here, we provide evidence for the existence of a deformed double shell closure in 80Zr through high precision Penning trap mass measurements of 80-83 Zr. Our mass values show that 80Zr is significantly lighter, and thus more strongly bound than predicted. This can be attributed to the deformed shell closure at N=Z=40 and the large Wigner energy. A statistical Bayesian model mixing analysis employing several global nuclear mass models demonstrates difficulties with reproducing the observed mass anomaly using current theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Integral boundary-layer models for turbulent separated flows

Two-dimensional, integral boundary-layer methods of the entrainment and dissipation type, with closure relationships which account for turbulent separated flows through the velocity and turbulence models, are developed. Kinematic shape factor relations used in the two methods are evaluated by comparison with recently obtained experimental data for separated flows. Prescribed-displacement-thickness computations are made with both methods and compared to experiment for a series of separated flows, including low speed and transonic shock-induced separation. Effects of upstream history on the turbulent shear stress development are shown to be important for the accurate prediction of the separated flows considered and the experimental trends are reasonably modeled by a rate equation derived from the turbulent kinetic energy equation.

Thomas, J. L.↗

Uncertainties in estimating global potential yields and their impacts for long-term modeling

Estimating realistic potential yields by crop type and region is challenging; such yields depend on both biophysical characteristics (e.g., soil characteristics, climate, etc.), and the crop management practices available in any site or region (e.g., mechanization, irrigation, crop cultivars). A broad body of literature has assessed potential yields for selected crops and regions, using several strategies. In this study we first analyze future potential yields of major crop types globally by two different estimation methods, one of which is based on historical observed yields (“Empirical”), while the other is based on biophysical conditions (“Simulated”). Potential yields by major crop and region are quite different between the two methods; in particular, Simulated potential yields are typically 200% higher than Empirical potential yields in tropical regions for major crops. Applying both of these potential yields in yield gap closure scenarios in a global agro-economic model, GCAM, the two estimates of future potential yields lead to very different outcomes for the agricultural sector globally. In the Simulated potential yield closure scenario, Africa, Asia, and South America see comparatively favorable outcomes for agricultural sustainability over time: low land use change emissions, low crop prices, and high levels of self-sufficiency. In contrast, the Empirical potential yield scenario is characterized by a heavy reliance on production and exports in temperate regions that currently practice industrial agriculture. At the global level, this scenario has comparatively high crop commodity prices, and more land allocated to crop production (and associated land use change emissions) than either the baseline or Simulated potential yield scenarios. This study highlights the importance of the choice of methods of estimating potential yields for agro-economic modeling.

60 APPLIED LIFE SCIENCES↗

Delamination Modeling of Composites for Improved Crash Analysis

Finite element crash modeling of composite structures is limited by the inability of current commercial crash codes to accurately model delamination growth. Efforts are made to implement and assess delamination modeling techniques using a current finite element crash code, MSC/DYTRAN. Three methods are evaluated, including a straightforward method based on monitoring forces in elements or constraints representing an interface; a cohesive fracture model proposed in the literature; and the virtual crack closure technique commonly used in fracture mechanics. Results are compared with dynamic double cantilever beam test data from the literature. Examples show that it is possible to accurately model delamination propagation in this case. However, the computational demands required for accurate solution are great and reliable property data may not be available to support general crash modeling efforts. Additional examples are modeled including an impact-loaded beam, damage initiation in laminated crushing specimens, and a scaled aircraft subfloor structures in which composite sandwich structures are used as energy-absorbing elements. These examples illustrate some of the difficulties in modeling delamination as part of a finite element crash analysis.

Fleming, David C.↗

Delamination Modeling of Composites for Improved Crash Analysis

Finite element crash modeling of composite structures is limited by the inability of current commercial crash codes to accurately model delamination growth. Efforts are made to implement and assess delamination modeling techniques using a current finite element crash code, MSC/DYTRAN. Three methods are evaluated, including a straightforward method based on monitoring forces in elements or constraints representing an interface; a cohesive fracture model proposed in the literature; and the virtual crack closure technique commonly used in fracture mechanics. Results are compared with dynamic double cantilever beam test data from the literature. Examples show that it is possible to accurately model delamination propagation in this case. However, the computational demands required for accurate solution are great and reliable property data may not be available to support general crash modeling efforts. Additional examples are modeled including an impact-loaded beam, damage initiation in laminated crushing specimens, and a scaled aircraft subfloor structures in which composite sandwich structures are used as energy-absorbing elements. These examples illustrate some of the difficulties in modeling delamination as part of a finite element crash analysis.

Fleming, David C.↗

Ion-Exchange Modeling of Crystalline Silicotitanate for Cesium Removal - 20283

The Tank Closure Cesium Removal (TCCR) system is a Savannah River Site (SRS) demonstration 'at-tank' process designed to remove {sup 137}Cs from the high-level aqueous tank waste so that the decontaminated solution can be disposed as low-level waste. Cesium is removed by ion exchange (IX) columns using engineered IONSIV{sup TM} R9120-B form of the Crystalline Silicotitanate (CST) media. The TCCR system is deployed at Tanks 10 and Tank 11 in the SRS H Tank Farm. Water is added to the salt-cake in Tank 10 H to dissolve it. The dissolved salt solution waste is pumped out of Tank 10H (feed tank), through filters and IX columns. The decontaminated salt solution is transferred to Tank 11 (receipt tank), and on to Tank 50H for final disposal in the Saltstone Production Facility. The current TCCR can accommodate lead-lag (two-column) or lead-lag-guard (three-column) configurations to optimize media utilization and achieve the target decontamination. To assist the TCCR operations, a parametric study was conducted to evaluate the impact of different parameters (e.g., column configurations (single column, two columns or three columns in series), waste characteristics, operating temperature, process flow rate, CST average particle size) on the IX column performance including CST bed utilization. The initial results indicate that the IX column performance is improved at slower process flow rate, at lower operating temperature, and with smaller CST average particle size. Multi-column configurations are recommended, because the single-column configuration does not utilize CST bed effectively. This paper demonstrates the versatility of the ion exchange modeling to evaluate the effects of CST characteristics and operational parameters on IX column performances. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Aerosol–Ice Formation Closure: A Southern Great Plains Field Campaign

Prediction of ice formation in clouds presents one of the grand challenges in the atmospheric sciences. Immersion freezing initiated by ice-nucleating particles (INPs) is the dominant pathway of primary ice crystal formation in mixed-phase clouds, where supercooled water droplets and ice crystals coexist, with important implications for the hydrological cycle and climate. However, derivation of INP number concentrations from an ambient aerosol population in cloud-resolving and climate models remains highly uncertain. We conducted an aerosol-ice formation closure pilot study using a field-observational approach to evaluate the predictive capability of immersion freezing INPs. The closure study relies on co-located measurements of the ambient size-resolved and single-particle composition and INP number concentrations. The acquired particle data serve as input in several immersion freezing parameterizations, that are employed in cloud-resolving and climate models, for prediction of INP number concentrations. We discuss in detail one closure case study in which a front passed through the measurement site, resulting in a change of ambient particle and INP populations. We achieved closure in some circumstances within uncertainties, but we emphasize the need for freezing parameterization of potentially missing INP types and evaluation of the choice of parameterization to be employed. Overall, this closure pilot study aims to assess the level of parameter details and measurement strategies needed to achieve aerosol-ice formation closure. The closure approach is designed to accurately guide immersion freezing schemes in models, and ultimately identify the leading causes for climate model bias in INP predictions.

ice formation↗

Closure measurements of naturally initiating small cracks

The initiation and growth of cracks in smooth 2024-T3 aluminum specimens have been investigated using acetate replicas. Crack opening displacements were measured as a function of load at several positions along the crack as it grew across the thickness of the specimen and became a through-the-thickness crack of a few millimeters in length. Tests run for R-ratios of 0.5, 0.0 and -1.0 showed that closure loads were well predicted by Newman's (1984) model. It is noted that small cracks grow slightly faster than would be predicted by long-crack data, and it is suggested that effects other than plasticity-induced closure must be taken into account.

Sharpe, W. N., Jr.↗

On reduced modelling of the modulational dynamics in magnetohydrodynamics

This paper explores structure formation in two-dimensional magnetohydrodynamic (MHD) turbulence as a modulational instability (MI) of turbulent fluctuations. We focus on the early stages of structure formation and consider simple backgrounds that allow for a tractable model of the MI while retaining the full chain of modulational harmonics. This approach allows us to systematically examine the validity of popular closures such as the quasilinear approximation and other low-order truncations. We find that, although such simple closures can provide quantitatively accurate approximations of the MI growth rates in some regimes, they can fail to capture the modulational dynamics in adjacent regimes even qualitatively, falsely predicting MI when the system is actually stable. We find that this discrepancy is due to the excitation of propagating spectral waves (PSWs) which can ballistically transport energy along the modulational spectrum, unimpeded until dissipative scales, thereby breaking the feedback loops that would otherwise sustain MIs. The PSWs can be self-maintained as global modes with real frequencies and drain energy from the primary structure at a constant rate until the primary structure is depleted. To describe these waves within a reduced model, we propose an approximate spectral closure that captures them and MIs on the same footing. We also find that introducing corrections to ideal MHD, conservative or dissipative, can suppress PSWs and reinstate the accuracy of the quasilinear approximation. In this sense, ideal MHD is a ‘singular’ system that is particularly sensitive to the accuracy of the closure within mean-field models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A unified view of convective transports by stratocumulus clouds, shallow cumulus clouds, and deep convection

A bulk planetary boundary layer (PBL) model was developed with a simple internal vertical structure and a simple second-order closure, designed for use as a PBL parameterization in a large-scale model. The model allows the mean fields to vary with height within the PBL, and so must address the vertical profiles of the turbulent fluxes, going beyond the usual mixed-layer assumption that the fluxes of conservative variables are linear with height. This is accomplished using the same convective mass flux approach that has also been used in cumulus parameterizations. The purpose is to show that such a mass flux model can include, in a single framework, the compensating subsidence concept, downgradient mixing, and well-mixed layers.

Randall, David A.↗

A physically based fractional cloudiness parameterization

A bulk PBL model with a simple internal vertical structure and a simple second-order closure has been developed for use as a PBL parameterization in a large-scale model. A method to determine the value of sigma, the fractional area covered by rising motion, is given. Large-eddy simulations of Moeng (1984, 1986) are used to validate various aspects of the parameterization.

Randall, David A.↗

Aircraft Observations of Turbulence in Cloudy and Cloud-Free Boundary Layers Over the Western North Atlantic Ocean From ACTIVATE and Implications for the Earth System Model Evaluation and Development

This study examines boundary layer turbulence derived from high temporal resolution meteorological measurements from 40 research flights over the western North Atlantic Ocean during the 2020 deployments of ACTIVATE. Frequency distributions of various turbulent quantities reveal stronger turbulence during the winter deployment than in summer and for cloud‐topped than in cloud‐free boundary layers during the summer deployment. Maximum turbulence kinetic energy (TKE) is most often within cloud from observations in winter and summer, whereas it is mostly below cloud in both seasons by a global model turbulence parameterization. Further, bivariate frequency distributions are consistent with the bivariate Gaussian probability distribution functions assumed for the closure of higher‐order turbulence/shallow convection parameterizations used by some global models. Turbulence simulated by the Community Atmosphere Model version 6 and the Energy Exascale Earth System Model Atmosphere Model version 2 using such parameterizations is not as strong as observed, with more TKE going into vertical wind perturbations rather than into zonal wind perturbations as observed, suggesting that the treatment of turbulence in Earth system models still needs to be further improved.

54 ENVIRONMENTAL SCIENCES↗

Better practices for inferring ecosystem water use strategy from eddy covariance data

Eddy covariance data are critical for inferring ecosystem water use strategies. Yet, such inferences are sensitive to a range of assumptions applied across studies, hindering our understanding of water use strategies within and across eddy covariance sites. A recent analysis across 151 FLUXNET2015 and AmeriFlux-FLUXNET datasets found that poor model performance was the key driver of non-robust inferences of ecosystem water use strategies. Here, we leverage this previous analysis to (i) identify the specific assumptions that improve inference model performance across most sites, (ii) explain the mechanisms behind the performance improvements, and (iii) check whether better performance improves water use inference. We find that the common practice of fitting a model to canopy conductance (G c ) derived from the evapotranspiration (ET) observations, rather than to observed ET itself, artificially amplifies data errors and degrades the model performance. Next, accounting for vegetation dynamics by applying a growing season filter or incorporating satellite LAI data improves performance, but the former practice may remove soil water stress periods. Lastly, using the leaf-to-air vapor pressure deficit (VPD l ) derived from ET observations as a model input may artificially inflate performance. Based on these results, we recommend selecting observed ET (rather than derived G c ) as the response variable, carefully accounting for vegetation dynamics, and avoiding derived VPD l as a model input; these best practices improve model performance by c. 20% and robustness by c. 80% across all eddy covariance sites. Nevertheless, the performance improvements do not always correspond to more robust inference of water use strategies, as model parameter selection and surface energy budget closure corrections still strongly influence the ecosystem water use parameter estimation in a site-specific manner.

AmeriFlux↗