Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic closures”

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 199 records · Page 11

Large Eddy Simulation Of Gravitational Effects In Transitional And Turbulent Gas-Jet Diffusion Flames

The influence of gravity on the spatial and the compositional structures of transitional and turbulent hydrocarbon diffusion flames are studies via large eddy simulation (LES) and direct numerical simulation (DNS) of round and planar jets. The subgrid-scale (SGS) closures in LES are based on the filtered mass density function (FMDF) methodology. The FMDF represents the joint probability density function (PDF) of the SGS scalars, and is obtained by solving its transport equation. The fundamental advantage of LES/FMDF is that it accounts for the effects of chemical reaction and buoyancy exactly. The methodology is employed for capturing some of the fundamental influences of gravity in equilibrium flames via realistic chemical kinetic schemes. Some preliminary investigation of the gravity effects in non-equilibrium flames is also conducted, but with idealized chemical kinetics models.

Jaberi, Farhad A.↗

Effects of expansion ratio on the calculation of parallel-walled backward-facing step flows - Comparison of four models of turbulence

This paper presents a numerical study of the effects of expansion ratio on two-dimensional separating and reattaching flows in plane backward-facing step geometries with parallel walls. Closure of the Reynolds equations was achieved by four different turbulence models: k-epsilon, 'modified' k-epsilon, algebraic stress, and 'modified' algebraic stress models. The k-epsilon model relates the Reynolds stresses to the mean rate of strain through the definition of an isotropic turbulent viscosity. The more advanced algebraic stress model calculates the stresses from implicit algebraic relationships containing the stresses themselves, the mean rate of strain, and the turbulent kinetic energy and its dissipation rate. 'Modified' versions of the models employ a new dissipation rate equation whose production term was made more sensitive to streamwise curvature effects. A new nonequilibrium wall function treatment proposed by Chieng and Launder (1980) was also incorporated into each model.

Sindir, M. M.↗

Characterizing Turbulence at a Forest Edge: Comparing Sub-Filter Scale Turbulence Models in Simulations of Flow over a Canopy

In wildfires, atmospheric turbulence plays a major role in the transfer of turbulent kinetic energy. Understanding how turbulence feeds back into a dynamical system is important, down to the varying small scales of fuel structures (i.e. pine needles, grass). Large eddy simulations (LES) are a common way of numerically representing turbulence. The Smagorinsky model (1963) serves as one of the most studied sub-grid scale representations in LES. In this investigation, the Smagorinsky model was implemented in HIGRAD/FIRETEC, LANL’s coupled fire-atmosphere model. This study was motivated by the need to quantitatively investigate the vorticity budget equation in HIGRAD/FIRETEC. The Smagorinsky turbulent kinetic energy (TKE) was compared to FIRETEC’s 1.5-order TKE eddy-viscosity subgrid-scale model, known as the Linn turbulence model. This was done in simulations of flow over flat terrain with a homogeneous, cuboidal canopy in the center of the domain. Examinations of the modeled vertical TKE profile and turbulent statistics at the leading edge, and throughout the canopy, show that the Smagorinsky model provides comparable results to that of the original closure model posed in FIRETEC.

58 GEOSCIENCES↗

Optimized Carbon Fiber Intermediate Development to Enable High-Volume Manufacturing of Lightweight Automotive Composites

The ongoing pursuit of improved fuel economy and reduced greenhouse gas emissions has resulted in sustained interest for lightweight materials technologies. In this context, carbon fiber composites have captured the imagination of automotive engineers due to the potential to achieve substantial mass reduction when compared to traditional steel construction. That stated, the use of carbon fiber composites in automotive has been limited for the most part to premium supercars and other derivative platforms. In these cases, manufacturing costs are less of an obstacle to implementation, and the performance benefits of carbon fiber have enabled production of structures offering more than 50% weight savings. In practice, translating these low-volume demonstrations onto high-volume vehicle platforms has remained challenging. This can be attributed to several factors, with the absence of suitable high throughput production methodologies being a key impediment. To date, structural, crash critical components have relied upon manufacturing techniques born out of the aerospace industry. This has created a disconnect between automotive production systems that are accustomed to manufacturing multiple parts per minute and the aerospace technologies that have cycle times in the order of hours. Consequently, the focus of this project is the development of manufacturing process technology for carbon fiber composites that can support a mainstream vehicle program at an assumed throughput of 100,000 vehicles per year. In practice, this translates to a part-to-part cycle time of less than 3 minutes. Project participation included contributions from a broad range of academic, industrial and national lab partners. The primary scope of work, being the development of new carbon fiber epoxy compounds that are stable at room temperature and suited to high throughput automated processing. For project management, the work streams were divided into six key areas, with the lead organization in parentheses. • Carbon fiber/epoxy materials formulation development and scale up (Dow) • Simulation of discontinuous near isotropic meso-structure intermediates (Purdue) • Simulation of mechanical performance of compression molded components (Purdue) • Meso-Scale morphological analysis and correlation with structural performance (UTK) • Paint and adhesion durability analysis (MSU) • Demonstrator part design, prototype production, and validation testing (Ford). The primary goal at the commencement of the project was development of a chopped carbon fiber sheet molding compound (SMC) that offered a three times improvement in tensile modulus over a comparable glass-based SMC. In addition to meeting mechanical performance targets, the resin kinetics were modified to achieve a processing cycle time of less than 3 minutes. Other critical-to-quality (CTQ) specifications were also stipulated to account for a broad range of materials and processing characteristics. To achieve the above, staff scientists at Dow Chemical created an extensive series of new epoxy blends for testing and validation. Throughout this development, a key challenge was attaining material performance goals without impacting processing behavior and paintability of finished components. The latter required a new internal mold release system being developed by Dow that was designed to complement the kinetics of the rapid cure epoxy. As a complement to work studies at the industrial partners, the teams from academia executed a series of analytical and experimental studies to investigate potential factors influencing CF-SMC performance. Unit cell models were developed to capture the meso-scale representations of the fiber matrix architecture. Results of this analysis and subsequent morphological investigations led to the design of a novel composite derivative comprising carbon fiber platelets embedded in an epoxy matrix; the platelet size and aspect ratio playing significant role in final composite properties. This approach was a departure from previous research in CF-SMC development whereby bulk filamentization or disassembly of the carbon fiber rovings had been considered the most effective means of achieving both fiber wet through and wet-out. As the course of the academia studies progressed, the aspect ratio of the fiber constituents was further optimized before finalizing material attributes and processing conditions. For the purposes of technology validation, the Ford team led a work stream devoted to the design, fabrication and testing of demonstration components. The carbon fiber SMC material has the potential to displace numerous stampings and castings on an automotive structure but ultimately vehicle closure applications were selected to showcase the abilities of the CF-SMC to achieve both mass reduction and business case for large complex structures. Using target properties established by the Dow staff scientists, the complete closure system for a full-size sedan decklid and a mid-size wagon liftgate were engineered. Prototypes for both applications were fabricated using production representative processing methods to allow for physical testing and performance validation of the CF-SMC structures. Following completion of a testing program that concluded with a FMVSS301 55 mph offset rear crash, the CF-SMC formulation was declared by the Ford team to have met all engineering requirements. To summarize, the joint development activities during this project led to significant technical breakthroughs and achievement of all milestones. The result was the development of a novel, tack-free carbon fiber molding compound that is suited to automated processing. This combined room temperature stability, fast cure kinetics, and internal mold release system facilitates cycle times that are conducive to high-volume production. The VORAFUSE M6400 successfully passed technology validation at Ford and is now eligible for consideration on future production commercial vehicle programs.

36 MATERIALS SCIENCE↗

Solutions of the chemical kinetic equations for initially inhomogeneous mixtures.

Following the recent discussions by O'Brien (1971) and Donaldson and Hilst (1972) of the effects of inhomogeneous mixing and turbulent diffusion on simple chemical reaction rates, the present report provides a more extensive analysis of when inhomogeneous mixing has a significant effect on chemical reaction rates. The analysis is then extended to the development of an approximate chemical sub-model which provides much improved predictions of chemical reaction rates over a wide range of inhomogeneities and pathological distributions of the concentrations of the reacting chemical species. In particular, the development of an approximate representation of the third-order correlations of the joint concentration fluctuations permits closure of the chemical sub-model at the level of the second-order moments of these fluctuations and the mean concentrations.

Hilst, G. R.↗

Comments on the present state and future directions of PDF methods

The one point probability density function (PDF) method is examined in light of its use in actual engineering problems. The PDF method, although relatively complicated, appears to be the only format available to handle the nonlinear stochastic difficulties caused by typical reaction kinetics. Turbulence modeling, if it is to play a central role in combustion modeling, has to be integrated with the chemistry in a way which produces accurate numerical solutions to combustion problems. It is questionable whether the development of turbulent models in isolation from the peculiar statistics of reactant concentrations is a fruitful line of development as far as propulsion is concerned. There are three issues for which additional viewgraphs are prepared: the one point pdf method; the amplitude mapping closure; and a hybrid strategy for replacing a full two point pdf treatment of reacting flows by a single point pdf and correlation functions. An appeal is made for the establishment of an adequate data base for compressible flow with reactions for Mach numbers of unity or higher.

Obrien, E. E.↗

Numerical Mixing Suppresses Submesoscale Baroclinic Instabilities Over Sloping Bathymetry

Abstract The impacts of spurious numerical salinity mixing on the larger‐scale flow and tracer fields are characterized using idealized simulations. The idealized model is motivated by realistic simulations of the Texas‐Louisiana shelf and features oscillatory near‐inertial wind forcing. can exceed the physical mixing from the turbulence closure in frontal zones and within the mixed layer. This suggests that simulated mixing processes in frontal zones are driven largely by . Near‐inertial alongshore wind stress amplitude is varied to identify a base case that maximizes the ratio of to in simulations with no prescribed horizontal mixing. We then test the sensitivity of the base case with three tracer advection schemes (MPDATA, U3HC4, and HSIMT) and conduct ensemble runs with perturbed bathymetry. Instability growth is evaluated using the volume‐integrated eddy kinetic energy and available potential energy . While all schemes have similar total mixing, the HSIMT simulations have over double the volume‐integrated and 20% less relative to other schemes, which suppresses the release of and reduces the by roughly 25%. This results in reduced isohaline variability and steeper isopycnals, evidence that enhanced suppresses instability growth. Differences in and between the MPDATA and U3HC4 simulations are marginal. However, the U3HC4 simulations have 25% more . Experiments with variable horizontal viscosity and diffusivity coefficients show that small amounts of prescribed horizontal mixing improve the representation of the ocean state for all advection schemes by reducing the and increasing the .

54 ENVIRONMENTAL SCIENCES↗

Toward a Unified Parameterization of the Boundary Layer and Moist Convection: Simulations of Clear and Cloudy Convection - Part 3

A model that employs a new form of mass-flux closure (described in Part I of this paper) is applied to a variety of clear and cloudy planetary boundary layers (PBLs) including dry convection from the Wangara Experiment, trade wind cumulus from the Barbados Oceanographic and Meteorological Experiment (BOMEX), and marine stratocumulus from the Atlantic Stratocumulus Experiment (ASTEX). For Wangara, the simulated variances and fluxes match that expected from similarity arguments, while the mean state is a little less mixed than the observations. In the BOMEX simulation. the shape and magnitude of the fluxes and the turbulence kinetic energy budget agree with LES results and observations. However, the liquid water mixing ratio is too large. This is attributed to an underprediction of the skewness. In agreement with observations from the ASTEX experiment, many of the model-simulated fields distinctly reflect a regime in transition between the trade wind cumulus and the classic stratocumulus-topped boundary layers. In general, the simulated entrainment rate tends to be a little underpredicted in regimes where there is little cloud-top radiative cooling (Wangara and BOMEX), while it is overpredicted in regimes where this process is more critical (e.g., ASTEX). Prior work suggests that this may be related to the manner in which the pressure terms are parameterized in the model. Overall, the model is able to capture some key physical features of these PBL regimes, and appears to have the potential to represent both cloud and boundary layer processes. Thus, this approach is a first step toward unifying these processes in large-scale models.

Lappen, Cara-Lyn↗

Shock noise features using the SCIPVIS code

The numerical computational supersonic jet plume model SCIPVIS was analyzed in order to determine its effectiveness for predicting the complicated flow structure associated with supersonic shock containing jet plumes. The analysis focused on assessing the applicability of using this code for predicting the broadband shock noise radiation. The analysis was carried out for several exit Mach number cases (including 1.00, 1.41, and 1.99) covering both underexpanded and overexpanded flow conditions. The three primary turbulence closure models incorporated in the SCIPVIS code were utilized for comparisons to measured plume static pressure distributions. Among other results, the prediction of plume static pressures for all three nozzle exit Mach numbers was found to perform exceptionally well when compared with experimental data for both underexpanded and overexpanded cases. Computations performed using this code for calculating the average shock wavelength for the 1.99 exit Mach number case were shown to provide an excellent prediction of the peak frequency for broadband shock noise. The spatial variation of the turbulent kinetic energy was found to be in qualitative agreement with wedge hot-film measurements.

Seiner, J. M.↗

Review and Summary of Corrosion Behavior for Aluminum-Clad Spent Nuclear Fuel in Dry Storage

A description of the corrosion behavior of aluminum alloys used for cladding on aluminum-based, aluminum-clad nuclear fuel used in research reactors under potential dry storage conditions has been compiled. An evaluation is made of potential additional corrosion with water postulated to not be removed during drying. The relative humidity (RH) produced by typical dryness criteria for industrial spent nuclear fuel drying is expected to be low, particularly at high temperatures (17% RH at 20°C and lower at higher temperatures). Existing corrosion data suggests that negligible vapor-phase corrosion of aluminum is expected at such low humidity, even for nominally bare aluminum surfaces. That is, the predicted relative humidity from free water (e.g., 17% RH) is below the “critical” relative humidity (~40% RH at room temperature), below which no significant corrosion is observed for bare aluminum. Furthermore, reported room-temperature corrosion rates are very low even under saturated (100% RH) water vapor. Existing results showing significant vapor-phase corrosion corresponded specifically to conditions of high temper ed with high relative humidity. A relatively large reservoir of (chemically bound) water exists in the aluminum (oxy)hydroxide films on the SNF cladding surface. The estimated total could saturate the gas even at relatively high temperature if fully released; however, its release as molecular water is expected to only be plausible if the temperature during storage exceeds both the drying temperature and the threshold for thermal decomposition of the trihydroxides (~220°C). Therefore, the combination of high temperature and high relative humidity that could drive significant corrosion is considered implausible during sealed dry storage following an appropriate drying process, to include >220°C drying for canisters that may approach or exceed this temperature during storage. Vapor-phase testing of aluminum samples with an adherent (oxy)hydroxide layer, prepared in liquid water to resemble those on actual aluminum-clad spent nuclear fuel (ASNF), did not observe evidence of additional corrosion even at combined high-temperature (up to 180°C) and high-humidity (up to 100% RH) conditions. Instead, small net mass losses were observed for most specimens and attributed to dehydration of the samples, which had been air-dried only prior to testing. ASNF being moved to dry storage is expected to have an existing (oxy)hydroxide film, suggesting that the cladding will be less prone to additional corrosion than bare aluminum metal. If significant corrosion did occur during sealed storage, the overall extent and impact is expected to be small. Corrosion post-closure of a canister would not alter the total amount of hydrogen in the canister, so it would have no effect on the maximum H 2 release already assessed in existing bounding calculations. In addition, the total amount of water in a 25-µm dense bayerite film would consume a only ~8 µm additional aluminum metal on average, if fully consumed by oxidizing aluminum metal to Al 2 O 3 . These conclusions are consistent with ASNF-in-canister simulations to date, which included a reaction pathway for corrosion using kinetics from previous literature and believed to be conservative and predicted very low rates of corrosion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

On the shear-current effect: toward understanding why theories and simulations have mutually and separately conflicted

ABSTRACT The shear-current effect (SCE) of mean-field dynamo theory refers to the combination of a shear flow and a turbulent coefficient β21 with a favourable negative sign for exponential mean-field growth, rather than positive for diffusion. There have been long-standing disagreements among theoretical calculations and comparisons of theory with numerical experiments as to the sign of kinetic ($\beta ^u_{21}$) and magnetic ($\beta ^b_{21}$) contributions. To resolve these discrepancies, we combine an analytical approach with simulations, and show that unlike $\beta ^b_{21}$, the kinetic SCE $\beta ^u_{21}$ has a strong dependence on the kinetic energy spectral index and can transit from positive to negative values at $\mathcal {O}(10)$ Reynolds numbers if the spectrum is not too steep. Conversely, $\beta ^b_{21}$ is always negative regardless of the spectral index and Reynolds numbers. For very steep energy spectra, the positive $\beta ^u_{21}$ can dominate even at energy equipartition urms ≃ brms, resulting in a positive total β21 even though $\beta ^b_{21}\lt 0$. Our findings bridge the gap between the seemingly contradictory results from the second-order-correlation approximation versus the spectral-τ closure, for which opposite signs for $\beta ^u_{21}$ have been reported, with the same sign for $\beta ^b_{21}\lt 0$. The results also offer an explanation for the simulations that find $\beta ^u_{21}\gt 0$ and an inconclusive overall sign of β21 for $\mathcal {O}(10)$ Reynolds numbers. The transient behaviour of $\beta ^u_{21}$ is demonstrated using the kinematic test-field method. We compute dynamo growth rates for cases with or without rotation, and discuss opportunities for further work.

79 ASTRONOMY AND ASTROPHYSICS↗

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗

A numerical study of a marine subtropical stratus cloud layer and its stability

A model for numerical simulation of stratus cloud layers is constructed by combining a second-order closure, turbulent transfer model with a thermal radiative transfer model. The turbulent transfer model allows water vapor saturation. The combined turbulence-radiation model is applied to both a horizontally uniform one-dimensional case and a horizontally nonuniform two-dimensional case. In the latter, the dynamics of mesoscale circulations are also incorporated. Results of the two-dimensional simulation show that the layer cloud instability occurs where the sea surface temperature is high and the large-scale subsidence is weak. The simulated instability is analyzed in view of an instability criterion, the eddy kinetic energy budget, and evaporative cooling near the cloud top.

Moeng, C.-H.↗

Reduced Stress Tensor and Dissipation and the Transport of Lamb Vector

We develop a methodology to ensure that the stress tensor, regardless of its number of independent components, can be reduced to an exactly equivalent one which has the same number of independent components as the surface force. It is applicable to the momentum balance if the shear viscosity is constant. A direct application of this method to the energy balance also leads to a reduction of the dissipation rate of kinetic energy. Following this procedure, significant saving in analysis and computation may be achieved. For turbulent flows, this strategy immediately implies that a given Reynolds stress model can always be replaced by a reduced one before putting it into computation. Furthermore, we show how the modeling of Reynolds stress tensor can be reduced to that of the mean turbulent Lamb vector alone, which is much simpler. As a first step of this alternative modeling development, we derive the governing equations for the Lamb vector and its square. These equations form a basis of new second-order closure schemes and, we believe, should be favorably compared to that of traditional Reynolds stress transport equation.

Wu, Jie-Zhi↗

Product distribution, kinetics, and aerosol formation from the OH oxidation of dimethyl sulfide under different RO 2 regimes

Abstract. The atmospheric oxidation of dimethyl sulfide (DMS) represents a major natural source of atmospheric sulfate aerosols. However, there remain large uncertainties in our understanding of the underlying chemistry that governs the product distribution and sulfate yield from DMS oxidation. Here, chamber experiments were conducted to simulate gas-phase OH-initiated oxidation of DMS under a range of reaction conditions. Most importantly, the bimolecular lifetime (τbi) of the peroxy radical CH3SCH2OO was varied over several orders of magnitude, enabling the examination of the role of peroxy radical isomerization reactions on product formation. An array of analytical instruments was used to measure nearly all sulfur-containing species in the reaction mixture, and results were compared with a near-explicit chemical mechanism. When relative humidity was low, “sulfur closure” was achieved under both high-NO (τbi<0.1 s) and low-NO (τbi>10 s) conditions, though product distributions were substantially different in the two cases. Under high-NO conditions, approximately half the product sulfur was in the particle phase, as methane sulfonic acid (MSA) and sulfate, with most of the remainder as SO2 (which in the atmosphere would eventually oxidize to sulfate or be lost to deposition). Under low-NO conditions, hydroperoxymethyl thioformate (HPMTF, HOOCH2SCHO), formed from CH3SCH2OO isomerization, dominates the sulfur budget over the course of the experiment, suppressing or delaying the formation of SO2 and particulate matter. The isomerization rate constant of CH3SCH2OO at 295 K is found to be 0.13±0.03 s−1, in broad agreement with other recent laboratory measurements. The rate constants for the OH oxidation of key first-generation oxidation products (HPMTF and methyl thioformate, MTF) were also determined (kOH+HPMTF=2.1×10-11 cm3 molec.−1 s−1, kOH+MTF=1.35×10-11 cm3 molec.−1 s−1). Product measurements agree reasonably well with mechanistic predictions in terms of total sulfur distribution and concentrations of most individual species, though the mechanism overpredicts sulfate and underpredicts MSA under high-NO conditions. Lastly, results from high-relative-humidity conditions suggest efficient heterogenous loss of at least some gas-phase products.

Ye, Qing (ORCID:0000000337978988)↗

Toward closure between predicted and observed particle viscosity over a wide range of temperatures and relative humidity

Abstract. Atmospheric aerosols can exist in amorphous semi-solid or glassy phase states whose viscosity varies with atmospheric temperature and relative humidity. The temperature and humidity dependence of viscosity has been hypothesized to be predictable from the combination of a water–organic binary mixing rule of the glass transition temperature, a glass-transition-temperature-scaled viscosity fragility parameterization, and a water uptake parameterization. This work presents a closure study between predicted and observed viscosity for sucrose and citric acid. Viscosity and glass transition temperature as a function of water content are compiled from literature data and used to constrain the fragility parameterization. New measurements characterizing viscosity of sub-100 nm particles using the dimer relaxation method are presented. These measurements extend the available data of temperature- and humidity-dependent viscosity to −28 ∘C. Predicted relationships agree well with observations at room temperature and with measured isopleths of constant viscosity at ∼107 Pa s at temperatures warmer than −28 ∘C. Discrepancies at colder temperatures are observed for sucrose particles. Simulations with the kinetic multi-layer model of gas–particle interactions suggest that the observed deviations at colder temperature for sucrose can be attributed to kinetic limitations associated with water uptake at the timescales of the dimer relaxation experiments. Using the available information, updated equilibrium phase-state diagrams (-80∘C

54 ENVIRONMENTAL SCIENCES↗

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.↗

Data-Driven Closures and Assimilation for Stiff Multiscale Random Dynamics

Here, we introduce a data-driven and physics-informed framework for propagating uncertainty in stiff, multiscale random ordinary differential equations (RODEs) driven by correlated (colored) noise. Unlike systems subjected to Gaussian white noise, a deterministic equation for the joint probability density function (PDF) of RODE state variables does not exist in closed form. Moreover, such an equation would require as many phase-space variables as there are states in the RODE system. To alleviate this curse of dimensionality, we instead derive exact, albeit unclosed, reduced-order PDF (RoPDF) equations for low-dimensional observables/quantities of interest. The unclosed terms take the form of state-dependent conditional expectations, which are directly estimated from data at sparse observation times. However, for systems exhibiting stiff, multiscale dynamics, data sparsity introduces regression discrepancies that compound during RoPDF evolution. This is overcome by introducing a kinetic-like defect term to the RoPDF equation, which is learned by assimilating in sparse, low-fidelity RoPDF estimates. Two assimilation methods are considered, namely nudging and deep neural networks, which are successfully tested against Monte Carlo simulations.

97 MATHEMATICS AND COMPUTING↗