Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parametric model dependence”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Estimating Temperature Retrieval Accuracy Associated With Thermal Band Spatial Resolution Requirements for Center Pivot Irrigation Monitoring and Management

This study explores the use of synthetic thermal center pivot irrigation scenes to estimate temperature retrieval accuracy for thermal remote sensed data, such as data acquired from current and proposed Landsat-like thermal systems. Center pivot irrigation is a common practice in the western United States and in other parts of the world where water resources are scarce. Wide-area ET (evapotranspiration) estimates and reliable water management decisions depend on accurate temperature information retrieval from remotely sensed data. Spatial resolution, sensor noise, and the temperature step between a field and its surrounding area impose limits on the ability to retrieve temperature information. Spatial resolution is an interrelationship between GSD (ground sample distance) and a measure of image sharpness, such as edge response or edge slope. Edge response and edge slope are intuitive, and direct measures of spatial resolution are easier to visualize and estimate than the more common Modulation Transfer Function or Point Spread Function. For these reasons, recent data specifications, such as those for the LDCM (Landsat Data Continuity Mission), have used GSD and edge response to specify spatial resolution. For this study, we have defined a 400-800 m diameter center pivot irrigation area with a large 25 K temperature step associated with a 300 K well-watered field surrounded by an infinite 325 K dry area. In this context, we defined the benchmark problem as an easily modeled, highly common stressing case. By parametrically varying GSD (30-240 m) and edge slope, we determined the number of pixels and field area fraction that meet a given temperature accuracy estimate for 400-m, 600-m, and 800-m diameter field sizes. Results of this project will help assess the utility of proposed specifications for the LDCM and other future thermal remote sensing missions and for water resource management.

Ryan, Robert E.↗
DOE OSTI2020

A machine learning approach to emulation and biophysical parameter estimation with the Community Land Model, version 5

Abstract. Land models are essential tools for understanding and predicting terrestrial processes and climate–carbon feedbacks in the Earth system, but uncertainties in their future projections are poorly understood. Improvements in physical process realism and the representation of human influence arguably make models more comparable to reality but also increase the degrees of freedom in model configuration, leading to increased parametric uncertainty in projections. In this work we design and implement a machine learning approach to globally calibrate a subset of the parameters of the Community Land Model, version 5 (CLM5) to observations of carbon and water fluxes. We focus on parameters controlling biophysical features such as surface energy balance, hydrology, and carbon uptake. We first use parameter sensitivity simulations and a combination of objective metrics including ranked global mean sensitivity to multiple output variables and non-overlapping spatial pattern responses between parameters to narrow the parameter space and determine a subset of important CLM5 biophysical parameters for further analysis. Using a perturbed parameter ensemble, we then train a series of artificial feed-forward neural networks to emulate CLM5 output given parameter values as input. We use annual mean globally aggregated spatial variability in carbon and water fluxes as our emulation and calibration targets. Validation and out-of-sample tests are used to assess the predictive skill of the networks, and we utilize permutation feature importance and partial dependence methods to better interpret the results. The trained networks are then used to estimate global optimal parameter values with greater computational efficiency than achieved by hand tuning efforts and increased spatial scale relative to previous studies optimizing at a single site. By developing this methodology, our framework can help quantify the contribution of parameter uncertainty to overall uncertainty in land model projections.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2002

Preliminary Sizing Completed for Single- Stage-To-Orbit Launch Vehicles Powered By Rocket-Based Combined Cycle Technology

Single-stage-to-orbit (SSTO) propulsion remains an elusive goal for launch vehicles. The physics of the problem is leading developers to a search for higher propulsion performance than is available with all-rocket power. Rocket-based combined cycle (RBCC) technology provides additional propulsion performance that may enable SSTO flight. Structural efficiency is also a major driving force in enabling SSTO flight. Increases in performance with RBCC propulsion are offset with the added size of the propulsion system. Geometrical considerations must be exploited to minimize the weight. Integration of the propulsion system with the vehicle must be carefully planned such that aeroperformance is not degraded and the air-breathing performance is enhanced. Consequently, the vehicle's structural architecture becomes one with the propulsion system architecture. Geometrical considerations applied to the integrated vehicle lead to low drag and high structural and volumetric efficiency. Sizing of the SSTO launch vehicle (GTX) is itself an elusive task. The weight of the vehicle depends strongly on the propellant required to meet the mission requirements. Changes in propellant requirements result in changes in the size of the vehicle, which in turn, affect the weight of the vehicle and change the propellant requirements. An iterative approach is necessary to size the vehicle to meet the flight requirements. GTX Sizer was developed to do exactly this. The governing geometry was built into a spreadsheet model along with scaling relationships. The scaling laws attempt to maintain structural integrity as the vehicle size is changed. Key aerodynamic relationships are maintained as the vehicle size is changed. The closed weight and center of gravity are displayed graphically on a plot of the synthesized vehicle. In addition, comprehensive tabular data of the subsystem weights and centers of gravity are generated. The model has been verified for accuracy with finite element analysis. The final trajectory was rerun using OTIS (Boeing Corporation's trajectory optimization software package), and the sizing output was incorporated into a solid model of the vehicle using PRO/Engineer computer-aided design software (Parametric Technology Corporation, Waltham, MA).

Roche, Joseph M.↗
NASA NTRS1986

Thermal balance of the atmospheres of Jupiter and Uranus

Two-dimensional, radiative-convective-dynamical models of the visible atmospheres of Jupiter and Uranus are presented. Zonally-averaged temperatures and heat fluxes are calculated numerically as functions of pressure and latitude. In addition to radiative heat fluxes, the dynamical heat flux due to large-scale baroclinic eddies is included and is parametrized using a mixing length theory which gives heat fluxes similar to those of Stone. The results for Jupiter indicate that the internal heat flow is non-uniform in latitude and nearly balances the net radiative flux leaving the atmosphere. The thermal emission is found to be uniform in latitude in agreement with Pioneer and Voyager observations. Baroclinic eddies are calculated to transport only a small amount of the meridional heat flow necessary to account for the uniformity of thermal emission with latitude. The bulk of the meridional heat transfer is found to occur very deep in the stable interior of Jupiter as originally proposed by Ingersoll and Porco. The relative importance of baroclinic eddies vs. internal heat flow in the thermal balance of Uranus depends on the ratio of emitted thermal power to absorbed solar power. The thermal balance of Uranus is compared to that of Jupiter for different values of this ratio.

Friedson, A. J.↗
DOE OSTI2022

On the hot surface ignition of a wall-stagnating spray flame

The ignition of fuel sprays due to interaction with hot surfaces is an important phenomenon in the safety analysis of many engineering systems. We perform a parametric study of the hot surface ignition (HSI) of a fuel spray approaching a heated surface caused by the accidental leakage of a fuel line. Here, to this end, we employ a one-dimensional Eulerian-Eulerian formulation with a non-equilibrium evaporation model and a realistic chemical mechanism to describe $n$-dodecane fuel chemistry. We first describe and analyze the phenomenology of the unsteady processes leading to ignition using non-dimensionalized quantities. Through consideration of the temporal development of the most reactive mixture, we demonstrate that ignition occurs at a fuel-lean composition in a premixed region near the hot surface. Using non-dimensional parameters identified from the governing equations, we perform a parametric study of the time, location and local mixture composition at ignition and determine the ignition limits. We then identify the most important parametric sensitivities for physical analysis using a data-driven classification method. Our analysis demonstrates a contraction of the ignition limits with increased Stokes number and a regime of parametric insensitivity of igniting mixture composition. We also show that at high Damköhler numbers, the ignition location conforms to the parametric behavior of the thermal boundary layer, whereas at low Damköhler numbers approaching the ignition limit it reaches a near-unity value of the quenching Peclet number. We then compare the demonstrated parametric dependencies to the results of the quasi-steady asymptotic ignition literature, showing that our results are consistent with those obtained analytically within the limitations imposed by the simplified formulation of the latter.

42 ENGINEERING↗
NASA NTRS1993

Status of the CNRS-LCSR program on high pressure droplet vaporization and burning

Depending on the surrounding flow and thermodynamic conditions, a single droplet may experience several gasification regimes, ranging from the envelope flame regime to pure vaporization. In practical situations, such as rocket propulsion or diesel combustion, the size distribution of droplets is, at best, bimodal, so that the possibility exists for the simultaneous presence of various regimes. For example, very small droplets are transported by the gas phase with zero relative velocity. This picture validates then the spherical symmetry hypothesis applied to the droplet and to the diffusion flame enveloping it. On the other hand, for larger droplets, a relative velocity exists due to drag forces. The most important influence of forced convection on droplet burning is the possibility to extinguish globally the envelope flame, or to establish a flame stabilized in the wake region. The burning rates of these regimes differ strongly. The characteristic time of droplet gasification is also influenced by the surrounding pressure and temperature. A parametric investigation of single droplet burning regimes is then helpful in providing the necessary physical ideas for sub-grid models used in spray combustion numerical prediction codes. The CNRS-LCSR experimental program on droplet vaporization and burning deals with these various regimes: stagnant and convective monocomponent droplet burning convective mono and bicomponent droplet vaporization; high temperature convective mono and biocomponent droplet vaporization; burning regimes of hydrazine and hydroxyl-ammonium-nitrate based monopropellant droplets and the vaporization regimes of liquid oxygen droplets. Studies on interacting droplets and on liquid aluminum droplets will start in the near future. The influence of high pressure is a common feature of all these studies. This paper summarizes the status of the CNRS-LCSR program on the effects of high pressure on monocomponent single droplet burning and vaporization, and some recent results obtained under normal and reduced gravity conditions with suspended droplets are presented. In the work described here, parabolic flights of an aircraft is used to create a reduced gravity environment of the order of 10(exp -2) g.

Chauveau, Christian↗
DOE OSTI2025

Uncertainty quantification of a physics-informed model based on sparse identification of a Thermal Energy Distribution System

Integrated energy systems (IES)s are crucial for enhancing the economy and efficiency of power generation sources (e.g., nuclear energy) necessary to unleash American energy dominance. These systems can be integrated with thermal energy storage (TES) and intermittent renewable energies to optimize overall energy use, peak-load regulation, and demand-side responses. However, the stabilization of energy generation, transport, and utilization introduces operational complexities that exceed the challenges of managing each sub-component individually. Currently, though IESs rely on human operators for efficiency and stability, reducing human error risk and enhancing performance through automation is highly desirable. Recent advances at Idaho National Laboratory have demonstrated successful control of the Thermal Energy Distributed System (TEDS). However, the automatic control system depends on a deterministic Sparse Identification of Nonlinear Dynamics with Control (SINDyC) model, which are trained based on simulation data from physics-based simulations. Because of uncertainties in physics-based simulation, SINDyC model results in large discrepancies against experimental data and cannot be reliably used in automatic control. In this paper, we present an innovative approach to address these discrepancies by quantifying uncertainties and developing a more robust model. We first generated trajectories by using first-principles physics codes to encapsulate the experiment. Next, we trained thousands of models by randomly sampling these trajectories. We then collapsed all those models into one probabilistic SINDyC by fitting a multivariate Gaussian distribution onto the resulting coefficient’s distribution. Despite its simplicity, our approach successfully produced 95% confidence intervals that captured the experimental trajectories. It even did so with a higher probability and better U-pooling score across six of the seven relevant quantities of interest (QoIs), as compared to other classical approaches. In conclusion, ongoing research is focusing on generating new experimental trajectories to validate this approach, and on employing Bayesian calibration to refine parametric uncertainties and guide future model development efforts.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2012

Parametric 3D Atmospheric Reconstruction in Highly Variable Terrain with Recycled Monte Carlo Paths and an Adapted Bayesian Inference Engine

We describe a method for accelerating a 3D Monte Carlo forward radiative transfer model to the point where it can be used in a new kind of Bayesian retrieval framework. The remote sensing challenge is to detect and quantify a chemical effluent of a known absorbing gas produced by an industrial facility in a deep valley. The available data is a single low resolution noisy image of the scene in the near IR at an absorbing wavelength for the gas of interest. The detected sunlight has been multiply reflected by the variable terrain and/or scattered by an aerosol that is assumed partially known and partially unknown. We thus introduce a new class of remote sensing algorithms best described as "multi-pixel" techniques that call necessarily for a 3D radaitive transfer model (but demonstrated here in 2D); they can be added to conventional ones that exploit typically multi- or hyper-spectral data, sometimes with multi-angle capability, with or without information about polarization. The novel Bayesian inference methodology uses adaptively, with efficiency in mind, the fact that a Monte Carlo forward model has a known and controllable uncertainty depending on the number of sun-to-detector paths used.

Monte Carlo↗
NASA NTRS1994

An error covariance model for sea surface topography and velocity derived from TOPEX/POSEIDON altimetry

In order to facilitate the use of satellite-derived sea surface topography and velocity oceanographic models, methodology is presented for deriving the total error covariance and its geographic distribution from TOPEX/POSEIDON measurements. The model is formulated using a parametric model fit to the altimeter range observations. The topography and velocity modeled with spherical harmonic expansions whose coefficients are found through optimal adjustment to the altimeter range residuals using Bayesian statistics. All other parameters, including the orbit, geoid, surface models, and range corrections are provided as unadjusted parameters. The maximum likelihood estimates and errors are derived from the probability density function of the altimeter range residuals conditioned with a priori information. Estimates of model errors for the unadjusted parameters are obtained from the TOPEX/POSEIDON postlaunch verification results and the error covariances for the orbit and the geoid, except for the ocean tides. The error in the ocean tides is modeled, first, as the difference between two global tide models and, second, as the correction to the present tide model, the correction derived from the TOPEX/POSEIDON data. A formal error covariance propagation scheme is used to derive the total error. Our global total error estimate for the TOPEX/POSEIDON topography relative to the geoid for one 10-day period is found tio be 11 cm RMS. When the error in the geoid is removed, thereby providing an estimate of the time dependent error, the uncertainty in the topography is 3.5 cm root mean square (RMS). This level of accuracy is consistent with direct comparisons of TOPEX/POSEIDON altimeter heights with tide gauge measurements at 28 stations. In addition, the error correlation length scales are derived globally in both east-west and north-south directions, which should prove useful for data assimilation. The largest error correlation length scales are found in the tropics. Errors in the velocity field are smallest in midlatitude regions. For both variables the largest errors caused by uncertainty in the geoid. More accurate representations of the geoid await a dedicated geopotential satellite mission. Substantial improvements in the accuracy of ocean tide models are expected in the very near future from research with TOPEX/POSEIDON data.

Tsaoussi, Lucia S.↗
DOE OSTI2022

Diffstar: a fully parametric physical model for galaxy assembly history

We present Diffstar, a smooth parametric model for the in situ star formation history (SFH) of galaxies. The Diffstar model is distinct from traditional SFH models because it is parametrized directly in terms of basic features of galaxy formation physics. Diffstar includes ingredients for: the halo mass assembly history; the accretion of gas into the dark matter halo; the fraction of gas that is eventually transformed into stars, ϵ ms ; the time-scale over which this transformation occurs, τ cons ; and the possibility that some galaxies will experience a quenching event at time t q , and may subsequently experience rejuvenated star formation. We show that our model is sufficiently flexible to describe the average stellar mass histories of galaxies in both the IllustrisTNG (TNG) and UniverseMachine (UM) simulations with an accuracy of ~0.1 dex across most of cosmic time. We use Diffstar to compare TNG to UM in common physical terms, finding that: (i) star formation in UM is less efficient and burstier relative to TNG; (ii) UM galaxies have longer gas consumption time-scales, relative to TNG; (iii) rejuvenated star formation is ubiquitous in UM, whereas quenched TNG galaxies rarely experience sustained rejuvenation; and (iv) in both simulations, the distributions of ϵ ms , τ cons , and t q share a common characteristic dependence upon halo mass, and present significant correlations with halo assembly history. Finally, we conclude with a discussion of how Diffstar can be used in future applications to fit the SEDs of individual observed galaxies, as well as in forward-modelling applications that populate cosmological simulations with synthetic galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Finite-temperature electron-capture rates for neutron-rich nuclei near N = 50 and effects on core-collapse supernova simulations

In this paper the temperature dependence of stellar electron-capture (EC) rates is investigated, with a focus on nuclei near N = 50, just above Z = 28, which play an important role during the collapse phase of core-collapse supernovae (CCSN). Two new microscopic calculations of stellar EC rates are obtained from relativistic and nonrelativistic finite-temperature quasiparticle random-phase approximation approaches, for a conventional grid of temperatures and densities. In both approaches, EC rates due to Gamow-Teller transitions are included. In the relativistic calculation, contributions from first-forbidden transitions are also included and add strongly to the EC rates. The new EC rates are compared with large-scale shell-model calculations for the specific case of 86 Kr, providing insight into the finite-temperature effects on the EC rates. At relevant thermodynamic conditions for core collapse, the discrepancies between the different calculations of this paper are within about one order of magnitude. Numerical simulations of CCSN are performed with the spherically symmetric gr1d simulation code to quantify the impact of such differences on the dynamics of the collapse. These simulations also include EC rates based on two parametrized approximations. A comparison of the neutrino luminosities and enclosed mass at core bounce shows that differences between simulations with different sets of EC rates are relatively small (≈ 5 %), suggesting that the EC rates used as inputs for these simulations have become well constrained.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2004

Nucleation in Synoptically Forced Cirrostratus

Formation and evolution of cirrostratus in response to weak, uniform and constant synoptic forcing is simulated using a one-dimensional numerical model with explicit microphysics, in which the particle size distribution in each grid box is fully resolved. A series of tests of the model response to nucleation modes (homogeneous-freezing-only/heterogeneous nucleation) and heterogeneous nucleation parameters are performed. In the case studied here, nucleation is first activated in the prescribed moist layer. A continuous cloud-top nucleation zone with a depth depending on the vertical humidity gradient and one of the nucleation parameters is developed afterward. For the heterogeneous nucleation cases, intermittent nucleation zones in the mid-upper portion of the cloud form where the relative humidity is on the rise, because existent ice crystals do not uptake excess water vapor efficiently, and ice nuclei (IN) are available. Vertical resolution as fine as 1 m is required for realistic simulation of the homogeneous-freezing-only scenario, while the model resolution requirement is more relaxed in the cases where heterogeneous nucleation dominates. Bulk microphysical and optical properties are evaluated and compared. Ice particle number flux divergence, which is due to the vertical gradient of the gravity-induced particle sedimentation, is constantly and rapidly changing the local ice number concentration, even in the nucleation zone. When the depth of the nucleation zone is shallow, particle number concentration decreases rapidly as ice particles grow and sediment away from the nucleation zone. When the depth of the nucleation zone is large, a region of high ice number concentration can be sustained. The depth of nucleation zone is an important parameter to be considered in parametric treatments of ice cloud generation.

Lin, R.-F.↗
NASA NTRS2024

Exploring Standing and Reflected Slow-mode Waves in Flaring Coronal Loops: A Parametric Study Using 2.5D MHD Modeling

Abstract Recent observations of reflected propagating and standing slow-mode waves in hot flaring coronal loops have spurred our investigation into their underlying excitation and damping mechanisms. To understand these processes, we conduct 2.5D magnetohydrodynamic (MHD) simulations using an arcade active region model that includes a hot and dense loop. Our simulations allow for in-depth parametric investigations complementing and expanding our previous 3D MHD modeling results. We excite these waves using a large-amplitude, flow pulse applied at one footpoint of the loop in two distinct models as motivated by observations from the Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA). The first model (Model 1) incorporates classical compressive viscosity coefficient, while the second model (Model 2) adopts a 10-times enhanced viscosity coefficient. We obtain the following major results: (1) Our 2.5D MHD simulations reinforce previous conclusions derived from 1D and 3D MHD models that significantly enhanced viscosity is crucial for the rapid excitation of standing slow waves with damping times consistent with observations by Wang et al. (2015). (2)We uncover that nonlinearity in Model 1 delays the conversion of a reflected propagating wave into a standing wave. In contrast, Model 2 exhibits a much weak influence of nonlinearity on the excitation time of standing waves, thanks to the suppression of these effects by enhanced viscosity. (3) Our results reveal that the transverse temperature structure holds more influence on wave behavior than the density structure. In Model 1, increased loop temperature contrast significantly enhances wave trapping within the structure, mitigating the impact of temperature-dependent viscous damping. Conversely, in Model 2, the impact of temperature structure on wave behavior weakens in comparison to the effect of viscosity. (4) Model 1 displays evident nonlinear coupling to the fast and kink magnetoacoustic waves and pronounced wave leakage into the corona. Model 2 exhibits significantly weaker effects in this regard. Analyzing three observed wave events by SDO/AIA aligns with Model 2 predictions, providing further support for the substantial viscosity increase. Our 2.5D study unravels the complex interplay of wave-flow phenomena and nonlinear processes in coronal loops, extending our previous 1D modeling results to incorporate more realistic loop geometry. This provides insights into scenarios where 3D effects may be neglected, thereby enhancing our understanding of the intricate dynamics of the solar corona.

Tongjiang Wang↗
DOE OSTI2023

Parametric study of Alfvénic instabilities driven by runaway electrons during the current quench in DIII-D

Abstract To avoid or mitigate runaway electron (RE) beams in ITER, RE-driven instabilities are actively studied as a complimentary technique to massive material injection. In this work we report experimental dependencies of Alfvénic instabilities driven by REs during the current quench in DIII-D on plasma and material injection parameters. These instabilities, observed in the frequency range of 0.1–3 MHz, correlate with increased RE loss and thus may play a role in non-sustained RE beams. It was found that as the toroidal magnetic field ( B T ) decreases, the RE population becomes more energetic, the energy of instabilities increases, and no RE beam is observed when the maximum energy of REs exceeds 15 MeV (or when B T is below 1.8 T). Analysis of disruptions at plasma core temperature ( T e ) of 1 keV and 8 keV shows that the RE population is much less energetic (with the maximum energy of only about 3 MeV) when T e is high, and no instabilities are observed in this case. Besides disruptions above caused by Ar injection, cases with Ne and D 2 injections were also studied. Both Ne and D 2 injections cause no sustained RE beams, however, for different reasons. Measurements of the instability polarization indicate that it is of predominantly compressional nature at the edge, which is consistent with modeling suggesting excitation of compressional Alfvén eigenmodes. However, drive of global Alfvén eigenmodes is also possible at low frequencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2002

Transient Mixing Driven by Buoyancy Flows

Mixing driven by buoyancy-induced flows is of particular interest to microgravity processes, as the body force that governs the intensity of flow fields can be directly controlled. We consider a model experimental system to explore the dynamics of mixing which employs two miscible liquids inside a cavity separated initially by a divider. The two liquids are oriented vertically inside a rectangular cavity with constant width and height, and varying depths to span the range of a Hele-Shaw cell to a 3-[) configuration. The two miscible liquids can be sufficiently diluted and died, for example water and deuterium oxide, such that a distinct interface exists across the divider. The transient mixing characteristic of the two fluids is addressed by following the Lagrangian history of the interface for various aspect ratios in the z-plane (depth variation) as well as a range of pulling velocities of the divider. The mixing characteristics of the two fluids are quantified from measurement of the length stretch of the interface and its flow field using respectively image processing techniques and Particle Imaging Velocimetry. Scaling analysis shows that the length stretch depends on four governing parameters, namely the Grashof number (Gr), Schmidt number (Sc), aspect ratio (Ar), and Reynolds number (Re). Variation of the Schmidt number is taken into account through thermophysical property variation. Thus our problem reduces to a codimension three bifurcation in parametric space for Gr, Ar, and Re. Our experimental results show that for Gr on the order of 106 and a nominal cavity aspect ratio Ar = 0.2, the net effect of removal of the divider and the overwhelming buoyancy force causes an overturning motion which stretches and folds the interface to produce ml internal breakwave. The structure of the breakwave is similar to the ubiquitous Rayleigh-Taylor instability morphology. The breakwave is dissipated either through internal or wall collision depending on the impulsive velocity of the divider as prescribed by the Reynolds number. The decay of the collision event occurs through sloshing oscillations over a short time scale. The two fluids then become stably stratified with a diffusive band at the interface indicating local mass transport. The local bifurcation of the internal breakwave is investigated as a function of aspect ratio. Results show that for narrow cavities on the order of 2mm (Ar = 0.04) folding does not occur, the interface only stretches. As the cavity size increases folding occurs through a supercritical bifurcation. Insight into the mechanism of folding is obtained from measurement of the flow field which shows that in the neighborhood of the folding event, there exists hyperbolic points caused by multiple vortex interactions. The global length stretch of the interface as a function of time is nearly Gaussian; calculations of finite-time Liapunov exponents as well as construction of horseshoe maps indicate the likelihood of a chaotic transient.

W.M.B. Duval↗
DOE OSTI2022

Stellar population analysis of MaNGA early-type galaxies: IMF dependence and systematic effects

We study systematics associated with estimating simple stellar population (SSP) parameters – age, metallicity [M/H], α-enhancement [α/Fe], and initial mass function (IMF) shape – and associated M * /L gradients, of elliptical slow rotators (E-SRs), fast rotators (E-FRs), and S0s from stacked spectra of galaxies in the MaNGA survey. These systematics arise from (i) how one normalizes the spectra when stacking; (ii) having to subtract emission before estimating absorption line strengths; (iii) the decision to fit the whole spectrum or just a few absorption lines; (iv) SSP model differences (e.g. isochrones, enrichment, IMF). The MILES+Padova SSP models, fit to the H β , 〈Fe〉, TiO 2SDSS , and [MgFe] Lick indices in the stacks, indicate that out to the half-light radius R e : (a) ages are younger and [α/Fe] values are lower in the central regions but the opposite is true of [M/H]; (b) the IMF is more bottom-heavy in the center, but is close to Kroupa beyond about R e /2; (c) this makes M * /L about 2 × larger in the central regions than beyond R e /2. While the models of Conroy et al. return similar [M/H] and [α/Fe] profiles, the age and (hence) M * /L profiles can differ significantly even for solar abundances and a Kroupa IMF; different responses to non-solar abundances and IMF parametrization further compound these differences. There are clear (model independent) differences between E-SRs, E-FRs, and S0s: younger ages and less enhanced [α/Fe] values suggest that E-FRs and S0s are not SSPs, but relaxing this assumption is unlikely to change their inferred M * /L gradients significantly.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Hydro-chemo-mechanical phase field formulation for corrosion induced cracking in reinforced concrete

Steel corrosion induced cracking in reinforced concrete structures is often caused by volumetric straining at the steel-concrete interface. This volumetric strain is mainly controlled by the rate of generation of corrosion product, which depends on the diffusion of Cl{sup −} and dissolved O{sub 2} as well as the chemical reaction at the steel-concrete interface. While many numerical modeling techniques have been proposed to tackle this complex phenomenon, most have been limited to using staggered solution schemes, where the diffusion, chemical reaction, mechanical strain, and crack fields are solved for separately. Unfortunately, such approaches introduce numerical errors that can lead to unrealistic predictions. Therefore, a coupled solution, is needed to remove the formulation induced errors and understand this complex phenomenon. To provide a unifying solution framework for corrosion induced cracking in reinforced concrete, a phase field formulation is presented and implemented with FEM algorithms. The governing equations are derived from the system's free energy and Fick's law of diffusion. The proposed formulation was first validated in 1D by comparing the modeling results with the semi-analytical solution approximated via Laurent series to avoid singularities in the complex domain. 2D plane strain solutions to the embedded rebar undergoing corrosion induced cracking with pre-cracks are also presented and compared with the approximated solutions from linear elastic fracture mechanics theory (LEFM). The modeling of crack propagation enabled by the proposed formulation is also presented along with parametric studies to reveal the roles of the fracture toughness and permeability of the concrete on the corrosion induced cracking. A case study was also conducted on H-pile steel with concrete jackets to demonstrate the feasibility of modeling corrosion induced fracture in reinforced concrete structures with complex steel-concrete interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star

A realistic parametrization of convection and convective boundary mixing in conventional stellar evolution codes is still the subject of ongoing research. Furthermore, to improve the current situation, multidimensional hydrodynamic simulations are used to study convection in stellar interiors. Such simulations are numerically challenging, especially for flows at low Mach numbers which are typical for convection during early evolutionary stages. We explore the benefits of using a low-Mach hydrodynamic flux solver and demonstrate its usability for simulations in the astrophysical context. Simulations of convection for a realistic stellar profile are analyzed regarding the properties of convective boundary mixing. The time-implicit Seven-League Hydro (SLH) code was used to perform multidimensional simulations of convective helium shell burning based on a 25 M ⊙ star model. The results obtained with the low-Mach AUSM + -up solver were compared to results when using its non low-Mach variant AUSM B + -up. We applied well-balancing of the gravitational source term to maintain the initial hydrostatic background stratification. The computational grids have resolutions ranging from 180 × 90 2 to 810 × 540 2 cells and the nuclear energy release was boosted by factors of 3 × 10 3 , 1 × 10 4 , and 3 × 10 4 to study the dependence of the results on these parameters.

79 ASTRONOMY AND ASTROPHYSICS↗