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 631 records · Page 35

Theoretical and Observational Determination of Global and Regional Radiation Budget, Forcing and Feedbacks at the Top-of-Atmosphere and Surface

Report consists of: 1. List of accomplishments 2. List of publications 3. Abstracts of published or submitted papers and 4. Subject invention disclosure. The accomplishments of the grant listed are: 1. Improved the third-order turbulence closure in cloud resolving models to remove the liquid water oscillation. 2. Used the University of California-Los Angeles (UCLA) large-eddy simulation (LES) model to provide data for radiation transfer testing. 3. Revised shortwave k-distribution models based on HITRAN 2000. 4. Developed a gamma-weighted two-stream radiative transfer model for radiation budget estimate applications. 5. Estimated the effect of spherical geometry to the earth radiation budget. 6. Estimated top-of-atmosphere irradiance over snow and sea ice surfaces. 7. Estimated the aerosol direct radiative effect at the top of the atmosphere. 8. Estimated the top-of-atmosphere reflectance of the clear-sky molecular atmosphere over ocean. 9. Developed and validated new set of Angular Distribution Models for the CERES TRMM satellite instrument (tropical) 10. Developed and validated new set of Angular Distribution Models for the CERES Terra satellite instrument (global) 11. Quantified the top-of-atmosphere direct radiative effect of aerosols over global oceans from merged CERES and MODIS observations 12 Clarified the definition of TOA flux reference level for radiation budget studies 13. Developed new algorithm for unfaltering CERES measured radiances 14. Used multiangle POLDER measurements to produce narrowband angular distribution models and examine the effect of scene identification errors on TOA albedo estimates 15. Developed and validated a novel algorithm called the Multidirectional Reflectance Matching (MRM) model for inferring TOA albedos from ice clouds using multi-angle satellite measurements. 16. Developed and validated a novel algorithm called the Multidirectional Polarized Reflectance Matching (MPRM) model for inferring particle shapes from ice clouds using multi-angle polarized satellite measurements. 17. Developed 4 advanced light scattering models including the three-dimensional (3D) uniaxial perfectly matched layer (UPML) finite-difference time-domain (FDTD) model. 18. Develop sunglint in situ measurement and study reflectance distribution in the sunglint area. 19. Lead a balloon-borne radiometer TOA albedo validation effort. 20. Developed a CERES surface UVB, UVA, and UV index product.

Loeb, Norman G.↗

Changes in Tropical Clouds and Atmospheric Circulation Associated with Rapid Adjustment Induced by Increased Atmospheric CO2 – A Multiscale Modeling Framework Study

The radiative heating increase due to increased CO2 concentration is the primary source for the rapid adjustment of atmospheric circulation and clouds. In this study, we investigate the rapid adjustment resulting from doubling of CO2 and its physical mechanism using a multiscale modeling framework (MMF). The MMF includes an advanced higher-order turbulence closure in its cloud-resolving model component and simulates realistic shallow and deep cloud climatology and boundary layer turbulence. The rapid adjustment over the tropics is characterized by 1) reduced ascent and descent strengths over the ocean, 2) increased lower tropospheric stability (LTS) over the subsidence region, 3) shoaling of planetary boundary layers over the ocean, 4) increased deep convection over lands and shift of cloud coverage from the ocean to lands, and 5) reduced sensible (SH) and latent heat (LH) fluxes over the oceanic deep convective regions. Unlike conventional general circulation models and another MMF, a reduction in the global-mean shortwave cloud radiative cooling is not simulated, due to the increase in low clouds at lower altitudes over the ocean, resulting from reduced cloud-top entrainment due to strengthened inversion. Changes in regional circulation play a key role in cloud changes and shift of cloud coverage to lands. Weaker energy transport resulting from water vapor and cloud CO2 masking effects reduces the upward motion and convective clouds in the oceanic regions. The ocean-land transports are linked to the partitioning of surface SH and LH fluxes that increases humidity over lands and enhances deep convection over the tropical lands.

Xu, Kuan-Man↗

Film Cooling Performance Prediction for Air and Supercritical CO 2

The current study explores the possibility of cooling the vanes and blades of a direct-fired sCO 2 turbine using film cooling. The operating conditions of a direct-fired sCO 2 cycle and thermophysical properties of the fluid at those conditions can alter the flow field characteristics of the coolant jet and its mixing with the mainstream. Very little information is present in the literature regarding the performance of film cooling geometries employing supercritical CO 2 . The objective of this study is to estimate the resulting film cooling effectiveness while also capturing the effects of the crossflow-to-mainstream velocity ratio on the coolant jet. Herein, a computational fluid dynamic model is used to study the coolant jet exiting a cylindrical hole located on a flat plate, with the coolant fed by an internal channel. Steady-state Reynolds-averaged Navier–Stokes equations were solved along with the (shear-stress transport) SST k–ω model to provide the turbulence closure. The operating conditions for the direct-fired sCO 2 turbine are obtained using an in-house Cooled Turbine Model. Numerical predictions revealed that the crossflow effects and jet lift-off were more pronounced in the case of sCO 2 when compared to air. Spatial distribution of flow field and cooling effectiveness are presented at different operating conditions.

30 DIRECT ENERGY CONVERSION↗

Film Cooling Performance Predictions For Air And Supercritical CO 2

The current study explores the possibility of cooling the vanes and blades of a direct-fired sCO 2 turbine using film cooling geometries. The operating conditions of a direct-fired sCO 2 cycle and thermophysical properties of the fluid at those conditions can alter the flow field characteristics of the coolant jet and its mixing with the mainstream. Very little information is present in the literature regarding the performance of film cooling geometries employing a supercritical fluid. The objective of this study is to estimate the resulting film cooling effectiveness while also capturing the effects of the crossflow-to-mainstream velocity ratio on the coolant jet. A computational fluid dynamic model is used to study the coolant jet exiting a cylindrical hole located on a flat plate, with the coolant fed by an internal channel. Steady-state RANS equations were solved along with the SST kω model to provide the turbulence closure. The operating conditions for the direct-fired sCO 2 turbine are obtained using an in-house Cooled Turbine Model. Numerical predictions revealed that the crossflow effects and jet lift-off were more pronounced in the case of sCO 2 when compared to air. Spatial distribution of flow field and cooling effectiveness will be presented at different operating conditions. Recommendations for future studies will be given wherever applicable.

42 ENGINEERING↗

Modeling Cost of Offshore Carbon Storage in Saline Reservoirs

Poster presentation to the AAPG-SEG-SPE CCUS 2024 Conference. Offshore saline reservoirs provide a significant and accessible resource for geologic carbon storage (CS). The offshore environment requires distinct approaches to site selection, operations, monitoring, and risk that affect the technoeconomic assessment of offshore CS projects. The National Energy Technology Laboratory (NETL) has developed a CS cost model for offshore saline reservoirs known as CO2_S_COM_Offshore. Based on NETL’s widely used CO2_S_COM cost model for onshore saline CS, CO2_S_COM_Offshore enables technoeconomic analysis of CS in offshore areas. This model comprehensively incorporates multiple facets of offshore CS projects, from regional evaluation and site selection to permitting, transport, operations, monitoring, site closure, and decommissioning. Developed to model cost for offshore United States (US) Exclusive Economic Zones, aspects of this model can be adapted to international projects. Presented at the SPE/AAPG/SEG Carbon Capture Utilization and Storage Conference in Houston, TX, March 11-13, 2024.

Mark-Moser, Mackenzie K.↗

Modeling MTS pyrolysis and SiC deposition kinetics using principal component analysis and neural networks

Accurate chemical kinetics modeling is crucial for improving the efficiency of chemical processing and synthesis of ceramic matrix composites. Detailed kinetic models are computationally expensive due to the large number of transported chemical species, while the simplified physics-based models, such as single-step global mechanisms, are efficient but often overlook key chemical intermediates and pathways. Recent deep learning approaches promise accurate and cost-effective models. Yet, they require additional closures for the transported nonlinear latent variables, complicating integration with existing solvers. In this work, we develop a hybrid linear—nonlinear reduced model for silicon carbide deposition from methyltrichlorosilane precursor by combining principal component analysis (PCA) and autoencoder (AE) neural network (NN) approaches. PCA is used to identify a smaller set of linear transport variables, enabling direct reuse of conventional transport solvers. NNs then reconstruct the full chemical state from these reduced variables. We demonstrate the method on a chemical vapor deposition reactor—comprising a gas-phase pyrolysis plug flow reactor and a heterogeneous surface reactor—over a wide range of temperatures, pressures, and residence times. Our PCA–AE model achieves high accuracy with only five transported scalars, achieving an eightfold cost reduction compared to detailed mechanisms, in both a priori (using data from the test set only) and a posteriori (coupled with a differential equation solver). In conclusion, notable errors arise primarily near training domain boundaries and for long residence times, indicating the need for domain shift indicators and better long-horizon predictions in future reduced chemistry model development.

autoencoder neural networks↗

Bio-Optics of the Chesapeake Bay from Measurements and Radiative Transfer Calculations

We combined detailed bio-optical measurements and radiative transfer (RT) modeling to perform an optical closure experiment for optically complex and biologically productive Chesapeake Bay waters. We used this experiment to evaluate certain assumptions commonly used when modeling bio-optical processes, and to investigate the relative importance of several optical characteristics needed to accurately model and interpret remote sensing ocean-color observations in these Case 2 waters. Direct measurements were made of the magnitude, variability, and spectral characteristics of backscattering and absorption that are critical for accurate parameterizations in satellite bio-optical algorithms and underwater RT simulations. We found that the ratio of backscattering to total scattering in the mid-mesohaline Chesapeake Bay varied considerably depending on particulate loading, distance from land, and mixing processes, and had an average value of 0.0128 at 530 nm. Incorporating information on the magnitude, variability, and spectral characteristics of particulate backscattering into the RT model, rather than using a volume scattering function commonly assumed for turbid waters, was critical to obtaining agreement between RT calculations and measured radiometric quantities. In situ measurements of absorption coefficients need to be corrected for systematic overestimation due to scattering errors, and this correction commonly employs the assumption that absorption by particulate matter at near infrared wavelengths is zero.

Tzortziou, Maria↗

Structural stability augmentation system design using BODEDIRECT: A quick and accurate approach

A methodology is presented for a modal suppression control law design using flight test data instead of mathematical models to obtain the required gain and phase information about the flexible airplane. This approach is referred to as BODEDIRECT. The purpose of the BODEDIRECT program is to provide a method of analyzing the modal phase relationships measured directly from the airplane. These measurements can be achieved with a frequency sweep at the control surface input while measuring the outputs of interest. The measured Bode-models can be used directly for analysis in the frequency domain, and for control law design. Besides providing a more accurate representation for the system inputs and outputs of interest, this method is quick and relatively inexpensive. To date, the BODEDIRECT program has been tested and verified for computational integrity. Its capabilities include calculation of series, parallel and loop closure connections between Bode-model representations. System PSD, together with gain and phase margins of stability may be calculated for successive loop closures of multi-input/multi-output systems. Current plans include extensive flight testing to obtain a Bode-model representation of a commercial aircraft for design of a structural stability augmentation system.

Goslin, T. J.↗

Machine learning for the identification of phase transitions in interacting agent-based systems: A Desai-Zwanzig example

Deriving closed-form analytical expressions for reduced-order models, and judiciously choosing the closures leading to them, has long been the strategy of choice for studying phase- and noise-induced transitions for agent-based models (ABMs). In this paper, we propose a data-driven framework that pinpoints phase transitions for an ABM—the Desai-Zwanzig model—in its mean-field limit, using a smaller number of variables than traditional closed-form models. To this end, we use the manifold learning algorithm Diffusion Maps to identify a parsimonious set of data-driven latent variables, and we show that they are in one-to-one correspondence with the expected theoretical order parameter of the ABM. We then utilize a deep learning framework to obtain a conformal reparametrization of the data-driven coordinates that facilitates, in our example, the identification of a single parameter-dependent ordinary differential equation (ODE) in these coordinates. Additionally, we identify this ODE through a residual neural network inspired by a numerical integration scheme (forward Euler). We then use the identified ODE—enabled through an odd symmetry transformation—to construct the bifurcation diagram exhibiting the phase transition.

97 MATHEMATICS AND COMPUTING↗

Composite Analysis for Low-Level Waste Disposal in the Hanford Site Central Plateau (FY 2020) (Rev. 2)

This document provides an updated Hanford Site composite analysis (CA). A Hanford Site CA was prepared and issued in 1998 (PNNL-11800) with an addendum provided in 2002 (PNNL-11800 Addendum 1). The CA was approved in 2002 (Frei, 2002) and has been maintained (DOE/RL-2000-29 and subsequent revisions) to support low-level waste disposal performance assessments (PAs) and disposal authorizations for facilities at the Hanford Site, including the following: Continued operation of the Environmental Restoration Disposal Facility (ERDF) and the 200 East and 200 West Low-Level Burial Grounds; Construction of the Integrated Disposal Facility (IDF); Forthcoming closure of tank residual waste systems such as Waste Management Area (WMA) C. The CA maintenance program resulted in a determination in 2015 (DOE/RL-2015-66) that the Hanford Site CA needed an update for the following reasons: While the initial Hanford Site CA has been maintained since 2001, the accumulation of basis changes reported in the annual summary reports over the succeeding 14 years merit evaluation in an updated analysis; The U.S. Department of Energy (DOE) Headquarters requested in a memorandum in 2015 (Gilbertson and Marcinowski, 2015) that “as soon as the relevant PAs are complete, the CA will be revised to account for all of the new information.” This updated Hanford Site CA accounts for the following new information: 1. Inclusion of a detailed Hanford Site baseline disposition that projects remedial activities through site closure. There have been significant changes through decision making in the Comprehensive Environmental Response, Compensation, and Liability Act of 1980 (CERCLA) process that were not available when the original CA was produced; 2. Inclusion of an updated inventory basis, new modeling capabilities, and new decisions reached in the associated record of decision (ROD) that was provided by issuance of DOE/EIS-0391 in fiscal year (FY) 2013. Development of a Hanford Site groundwater model from the baseline provided in a technical transfer of models for the Final Tank Closure and Waste Management Environmental Impact Statement (DOE/EIS-0391) commenced in FY 2014 and resulted in the plateau-to-river (P2R) groundwater flow model that is used in this updated Hanford Site CA; 3. Accounting for pump and treat systems, which were not evaluated in the initial CA. Pump and treat systems have had significant impact on groundwater flow system behavior, contaminant transport, and contaminant removal from Hanford Site groundwater. This process is accounted for in the use of the P2R model for this updated Hanford Site CA; 4. Inclusion of water-level data collected since the initial CA was completed. Data collected as the unconfined aquifer water levels continue to recede since the cessation of large liquid discharges in the late 1990s have led to marked improvement in understanding of the flow system for future conditions, particularly regarding northward flow potential in the critical Gable Gap area. The calibration in the latest version of the P2R model accounts for this information; 5. Inclusion of results from the revised ERDF PA completed in 2013 (WCH-5209) that account for updated inventory and expansion of the ERDF facility to about twice the size that was evaluated in the original CA; 6. Use of updated geoframeworks to provide the structural basis for numerical fate and transport models in the groundwater pathway of this updated CA. The geologic basis for groundwater models has continued to improve with additional data collection and interpretation with the creation and maintenance of the Hanford South Geoframework and the Central Plateau Vadose Zone Geoframework tools; 7. Incorporation of updated tank residual inventory estimates. Tank residual inventory estimates have improved with the incorporation of tank retrieval inventory data for those tanks that have completed retrieval. The CA inventory data package includes this updated information; 8. Incorporation of WMA C PA results. Two additional PAs for tank farm closure decisions are in preparation during the period required to prepare an updated Hanford Site CA: WMA C (FY 2016) and WMA A-AX (in preparation). The WMA C PA results are incorporated into the updated Hanford Site CA, and its grouted residuals model is used as the basis for a release model to account for the other tank farm systems modeling in this updated CA; 9. An update to the IDF PA was submitted in FY 2017 and has been reviewed and approved. The results of the IDF PA are incorporated into this updated CA; 10. Updated risk assessment scenarios. The risk assessment scenarios currently in use for Hanford Site CERCLA and Resource Conservation and Recovery Act of 1976 (RCRA) analyses differ from those evaluated in the initial Hanford Site CA. The representative person exposure scenario evaluated in this updated Hanford Site CA is consistent with recent PAs and CERCLA and RCRA analyses as the Hanford Site. This updated Hanford Site CA provides the following: A comparison of the updated Hanford Site CA all-pathway dose results with the performance measures during the compliance period, which is assumed to begin with site closure in calendar 2070 with the last scheduled disposal action and continue for 1,000 years postclosure (to calendar year 3069); A comparison of Hanford Site CA all-pathway dose results with the performance measures during the postcompliance period to address potential peaks beyond the compliance period. This is accomplished by evaluating dose in the period for 9,000 years following the compliance period (i.e., from calendar year 3070 to 12070).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A new stochastic systems approach to structural integrity

This paper develops improved stochastic models for the description of a large variety of fatigue crack growth phenomena that occur in components of considerable importance to the functionality and reliability of complex engineering structures. In essence, the models are based on the McGill-Markov and Closure-Lognormal stochastic processes. Not only do these models have the capability of predicting the statistical dispersion of crack growth rates, they also, by incorporating the concept of crack closure, have the capability of transferring stochastic crack growth properties measured under ideal laboratory conditions to situations of industrial significance, such as those occurring under adverse loading and/or environmental conditions. The primary data required in order to be in a position to estimate the pertinent parameters of these stochastic models are obtained from a statistically significant number of replicate tests. In this paper, both the theory and the experimental technique are illustrated using a Ti-6Al-4V alloy. Finally, important structural integrity, reliability, availability and maintainability concepts are developed and illustrated.

Provan, James W.↗

Experience with turbulence interaction and turbulence-chemistry models at Fluent Inc.

This viewgraph presentation discusses (1) turbulence modeling: challenges in turbulence modeling, desirable attributes of turbulence models, turbulence models in FLUENT, and examples using FLUENT; and (2) combustion modeling: turbulence-chemistry interaction and FLUENT equilibrium model. As of now, three turbulence models are provided: the conventional k-epsilon model, the renormalization group model, and the Reynolds-stress model. The renormalization group k-epsilon model has broadened the range of applicability of two-equation turbulence models. The Reynolds-stress model has proved useful for strongly anisotropic flows such as those encountered in cyclones, swirlers, and combustors. Issues remain, such as near-wall closure, with all classes of models.

Choudhury, D.↗

Aerosol-Ice Formation Closure Pilot Study (AEROICESTUDY) Field Campaign Report

The aerosol-ice formation closure pilot study (AEROICESTUDY) was a field campaign conducted at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility at the Southern Great Plains (SGP) observatory from October 7 to 28, 2019. The purpose of this campaign was to test a field observational approach for conducting an aerosol-ice formation closure study. In other words, the goal was to predict the number of ice nucleating particles (INPs) from the measurements of physical and chemical aerosol properties and compare those to measured INP number concentrations. This project was motivated by the fact that the aerosol community has widely conducted aerosol radiative closure, aerosol-cloud condensation nuclei (CCN), and CCN-droplet closure studies to test the physical models and parameterizations that cloud-resolving and climate models rely on to perform reliable simulations of the Earth system and energy budget. However, very few closure studies related to INPs have been conducted, and to our knowledge, none using robust, size-resolved, ambient aerosol composition measurements as model inputs

AEROICESTUDY↗

North Atlantic Aerosol Radiative Impacts Based on Satellite Measurements and Aerosol Intensive Properties from TARFOX and ACE-2

We estimate the impact of North Atlantic aerosols on the net short-wave flux at the tropopause by combining satellite-derived aerosol optical depth (AOD) maps with model aerosol properties determined via closure analyses in TARFOX and ACE 2. We exclude African dust, primarily by restricting latitudes to 25-60 N. The analyses use in situ aerosol composition measurements and air- and ship-borne sun-photometer measurements of AOD spectra. The aerosol model yields computed flux sensitivities (dFlux/dAOD) that agree with measurements by airborne flux radiometers in TARFOX. Its midvisible single-scattering albedo is 0.9. which is in the range obtained from in situ measurements of scattering and absorption in both TARFOX and ACE 2. Combining satellite-derived AOD maps with the aerosol model yields maps of 24-hour average net radiative flux changes. For simultaneous AVHRR, radiance measurements exceeded the sunphotometer AODs by about 0.04. However. shipboard sunphotometer and AVHRR AODs agreed Within 0.02 for data acquired during satellite overflights on two other days. We discuss attempts to demonstrate column closure within the MBL by comparing shipboard sunphotometer AODs and values calculated from simultaneous shipboard in-situ aerosol size distribution measurements. These comparisons were mostly unsuccessful, but they illustrate the difficulties inherent in this type of closure analysis. Specifically, AODs derived from near-surface in-situ size distribution measurements are extremely sensitive to the assumed hygroscopic growth model that itself requires an assumption of particle composition as a function of height and size, to the radiosonde-measured relative humidity, and to the vertical profile of particle number. We investigate further the effects of hygroscopic particle growth within the MBL by using shipboard lidar aerosol backscatter profiles together with the sunphotometer AOD.

Russell, P. B.↗

Fluidization of Group A Glass Particles: Experiments and Preliminary Validation

This work is aimed at providing reliable high-quality data from fluidization experiments. Geldart Group A glass beads are used as bed material. The fluidizing medium is air at atmospheric pressure which enters the system at 40% relative humidity. The air inlet velocity is varied over a broad range from 2 U mf to 26 U mf , where the minimum fluidization velocity, U mf , is 0.608 cm/s. There is not a significant change in the mean values of differential pressure; however, there are noticeable trends in their standard deviation values. Preliminary validation studies are performed using the open-source software Multiphase Flow with Interphase eXchanges (MFiX) developed by the National Energy Technology Laboratory (NETL). The results from Eulerian-Eulerian (MFiX-TFM, two-fluid model) and Eulerian-Lagrangian (MFiX-PIC, particle-in-cell) formulations are compared. The work presented in this report highlights the applicability of MFiX-PIC for large-scale applications, and reveals that PIC simulation offers a suitable alternative to TFM simulation when effects due to multiple components or polydispersity are significant. PIC modeling provides a readily extendible framework for such systems. However, PIC models do rely on an empirical closure for inter-particle stress. Some solution discrepancies were evident in MFiX-PIC predictions which require further analysis. As such, the study underlines the need for a systematic approach to quantify the sources of numerical uncertainty within the PIC model. However, even with known empiricism, MFiX-PIC has demonstrated considerable potential in analyzing gas-solid systems. The results obtained are encouraging and warrant further investigation to improve the predictive capability of MFiX-PIC.

36 MATERIALS SCIENCE↗

Aerosol-Ice Formation Closure Pilot Study (AEROICSTUDY) Field Campaign Report

The aerosol-ice formation closure pilot study (AEROICESTUDY) was a field campaign conducted at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility at the Southern Great Plains (SGP) observatory from October 7 to 28, 2019. The purpose of this campaign was to test a field observational approach for conducting an aerosol-ice formation closure study. In other words, the goal was to predict the number of ice nucleating particles (INPs) from the measurements of physical and chemical aerosol properties and compare those to measured INP number concentrations. This project was motivated by the fact that the aerosol community has widely conducted aerosol radiative closure, aerosol-cloud condensation nuclei (CCN), and CCN-droplet closure studies to test the physical models and parameterizations that cloud-resolving and climate models rely on to perform reliable simulations of the Earth system and energy budget. However, very few closure studies related to INPs have been conducted, and to our knowledge, none using robust, size-resolved, ambient aerosol composition measurements as model inputs.

54 ENVIRONMENTAL SCIENCES↗

Aerosol-Ice Formation Closure Pilot Study (AEROICSTUDY) Field Campaign Report

The aerosol-ice formation closure pilot study (AEROICESTUDY) was a field campaign conducted at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility at the Southern Great Plains (SGP) observatory from October 7 to 28, 2019. The purpose of this campaign was to test a field observational approach for conducting an aerosol-ice formation closure study. In other words, the goal was to predict the number of ice nucleating particles (INPs) from the measurements of physical and chemical aerosol properties and compare those to measured INP number concentrations. This project was motivated by the fact that the aerosol community has widely conducted aerosol radiative closure, aerosol-cloud condensation nuclei (CCN), and CCN-droplet closure studies to test the physical models and parameterizations that cloud-resolving and climate models rely on to perform reliable simulations of the Earth system and energy budget. However, very few closure studies related to INPs have been conducted, and to our knowledge, none using robust, size-resolved, ambient aerosol composition measurements as model inputs.

aerosol-ice formation↗

Modeling of Turbulent Swirling Flows

Aircraft engine combustors generally involve turbulent swirling flows in order to enhance fuel-air mixing and flame stabilization. It has long been recognized that eddy viscosity turbulence models are unable to appropriately model swirling flows. Therefore, it has been suggested that, for the modeling of these flows, a second order closure scheme should be considered because of its ability in the modeling of rotational and curvature effects. However, this scheme will require solution of many complicated second moment transport equations (six Reynolds stresses plus other scalar fluxes and variances), which is a difficult task for any CFD implementations. Also, this scheme will require a large amount of computer resources for a general combustor swirling flow. This report is devoted to the development of a cubic Reynolds stress-strain model for turbulent swirling flows, and was inspired by the work of Launder's group at UMIST. Using this type of model, one only needs to solve two turbulence equations, one for the turbulent kinetic energy k and the other for the dissipation rate epsilon. The cubic model developed in this report is based on a general Reynolds stress-strain relationship. Two flows have been chosen for model evaluation. One is a fully developed rotating pipe flow, and the other is a more complex flow with swirl and recirculation.

Shih, Tsan-Hsing↗