Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,621 records · Page 90

A prion-like protein regulator of seed germination undergoes hydration-dependent phase separation

Many organisms evolved strategies to survive desiccation. Plant seeds protect dehydrated embryos from various stressors and can lay dormant for millennia. Hydration is the key trigger to initiate germination, but the mechanism by which seeds sense water remains unresolved. Here we identified an uncharacterized Arabidopsis thaliana prion-like protein we named FLOE1, which phase separates upon hydration and allows the embryo to sense water stress. We demonstrate that biophysical states of FLOE1 condensates modulate its biological function in vivo in suppressing seed germination under unfavorable environments. We find intragenic, intraspecific, and interspecific natural variation in FLOE1 expression and phase separation and show that intragenic variation is associated with adaptive germination strategies in natural populations. This combination of molecular, organismal, and ecological studies uncovers FLOE1 as a tunable environmental sensor with direct implications for the design of drought-resistant crops, in the face of climate change.

59 BASIC BIOLOGICAL SCIENCES↗

Rheological and pumping behaviour of 3D printable cementitious materials with varying aggregate content

3D printing with cementitious materials is an emerging technology and provides an active area of research. Despite this, there are only limited studies that examine the pumping behaviour of high yield stress fluids such as printable concrete. We present a study on the rheological behaviour of 3D printable cementitious materials with different aggregate-to-binder ratios (a/b) and analyse the pumping and printing characteristics with the measured rheological parameters. We observed that an increase in a/b resulted in a significant increase of plastic viscosity, and a nominal increase in yield stress and storage modulus. The variation in plastic viscosity and yield stress with increasing aggregate content was analysed with the Krieger-Dougherty and Chateau-Ovarlez-Trung model, respectively. Also, the rheological properties of the paste rich lubrication layer during pumping were determined with the help of a tribometer. The yield stress and viscous constant of the lubrication layer increase with the increase in a/b. The pumping pressure predictions were carried out by a suitable theoretical model by considering appropriate flow conditions for different mixtures. The variation in critical height of printed elements with an increase in a/b was correlated with the observed rheological parameters.

36 MATERIALS SCIENCE↗

Fe Kα XANES, Fe Kβ HERFD XANES and EPMA flank method determinations of the oxidation state of Fe in garnet

The ferric to total iron ratios (Fe 3+ /ΣFe) of garnets can be paired with thermodynamic mineral activity models to quantify the oxygen fugacity of garnet-bearing rocks. However, techniques with a high analytical and spatial resolution are necessary to distinguish differences in garnet Fe 3+ /ΣFe ratios at the percent level and to accurately measure garnets that are zoned or contain inclusions. We acquired conventional Fe Kα and high-resolution energy fluorescence detection (HERFD) Fe Kβ X-ray absorption near edge structure (XANES) spectra and electron microprobe flank method analyses on a suite of 27 peridotitic and eclogitic garnets with Fe 3+ /ΣFe ratios previously determined by Mössbauer spectroscopy to evaluate the precision of each technique. We examined variations in the energy and intensity of three XANES spectral features as a function of Fe 3+ /ΣFe ratios: 1) the intensity ratio of two-post edge features (I-ratio; Fe Kα only); 2) the energy of the Fe edge at 90% normalized intensity (E 0.9 ; Fe Kα only) and 3) the pre-edge centroid energy (Fe Kα and HERFD Fe Kβ). In accordance with previous work, we find the energies of garnet pre-edge centroids are relatively insensitive to Fe 3+ /ΣFe ratios. The I-ratios of peridotitic and eclogitic garnets are offset from each other at low Fe 3+ /ΣFe ratios (≤0.13); I-ratio garnet XANES calibrations are composition-specific. The E 0.9 feature is independent of garnet major element composition in spectra that have been corrected for the effects of self-absorption. We produce two Fe Kα garnet XANES calibrations based on variations in the E0.9 feature; one calibration with all garnet reference materials included (Fe 3+ /ΣFe up to 1.0; “all garnet calibration”) and another calibration specific to garnets with low Fe 3+ /ΣFe ratios (“low ferric calibration”). Fe 3+ /ΣFe ratios calculated from the mean of up to 25 flank method measurements on eight garnet reference materials fall within 4% absolute of a one-to-one correlation with Fe 3+ /ΣFe ratios measured by Mössbauer. The standard error of the mean Fe 3+ /ΣFe ratio calculated from flank method approaches the Mössbauer-determined Fe 3+ /ΣFe ratio within estimated error (3%) after three analyses. Flank method precision is enhanced at higher beam current; however, the precision of the flank method does not approach the precision of XANES under any microprobe analytical condition tested here. Garnet reference materials detailed here are available by request to the Smithsonian Institution.

58 GEOSCIENCES↗

The effect of autoionization on the HBr + X 2 Π 3/2,1/2 state photoelectron angular distributions

A double-imaging photoelectron-photoion spectrometer and synchrotron radiation have been used to measure the HBr + X 2 Π 3/2 v + = 0, 1, 2 and the 2 Π 1/2 v + = 0 state photoelectron angular distributions, as characterized by the anisotropy parameter β , the total ion yield, and the threshold photoelectron spectrum. Particular attention has been focussed on the photon energy range between the 2 Π 3/2 and the 2 Π 1/2 spin–orbit components of the ground ionic state. This region encompasses Rydberg states, belonging to series converging onto the upper 2 Π 1/2 ionization limit, which may decay by autoionization into the 2 Π 3/2 ionization continuum. A detailed study has been performed on the effects of autoionization on the 2 Π 3/2 v + = 0 state photoelectron angular distributions. Furthermore, the observed energy dependent variations in the β -values exhibit a regular pattern that correlates with excitation into members of a very broad d-type Rydberg series. Additional rapid variations in the β -parameters, which occur over a narrow energy range, appear to coincide with sharp autoionizing Rydberg states belonging to s, p and d series. The present experimental results for the HBr + X 2 Π 3/2 v + = 0 state photoelectron anisotropy parameter are compared to previously reported theoretical predictions and to earlier studies of the Kr + 4p 5 2 P 3/2 state β -parameter. The threshold photoelectron spectrum of the X 2 Π 3/2 v + = 0 band exhibits partially resolved rotational structure. A simulation of this structure yields an ionization threshold of 11.6673 ± 0.0010 eV, which is consistent with previous measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Weighted nodal domain averages of eigenstates for quantum Monte Carlo and beyond

In this report we study the nodal properties of many-body eigenstates of stationary Schrödinger equation that affect the accuracy of real-space quantum Monte Carlo calculations. In particular, we introduce weighted nodal domain averages that provide a new probe of nodal surfaces beyond the usual expectations. Particular choices for the weight function reveal, for example, that the difference between two arbitrary fermionic eigenvalues is given by the nodal hypersurface integrals normalized by overlaps with the bosonic ground state of the given Hamiltonian. Noninteracting and fully interacting Be atom with corresponding almost exact and approximate wave functions are used to illustrate several aspects of these concepts. Variational formulations that employ different weights are proposed for prospective improvement of nodes in variational and fixed-node diffusion Monte Carlo calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Discussion on molecular dynamics (MD) simulations of the asphalt materials

The application of asphalt materials in pavement engineering has been increasingly widespread and sophisticated over the past several decades. Variations in the properties of asphalt binder during mixing, transportation, and paving can affect the performance of asphalt pavement. However, the asphalt material is a non-homogeneous and complex organic substance, consisting of various molecules with widely various molecular weights, elemental compositions, and structures. This complexity leads to difficulties for researchers to clearly and immediately understand the properties of asphalt materials and their variations. The multi-scale research approach combines macroscopic experimental data and microscopic simulation results from a practical engineering perspective. It helps to improve the understanding of asphalt materials. The molecular dynamics (MD) simulation proposes a corresponding molecular model of asphalt material based on experimental data, and the simulation algorithm is able to derive properties similar to those of real asphalt. Here, this paper provides a comprehensive review of the current studies on MD simulation of asphalt materials, including modeling, properties, and multi-scale analysis. As a key part of the computational simulation, this paper discusses the typical asphalt binder and asphalt-aggregate interface models constructed by different groups, and also presents their differences from real samples and their feasibility based on fundamental properties. After the introduction of molecular models, the extensive work made by researchers based on molecular models is categorically reviewed and discussed. The strengths and weaknesses of MD simulation methods in the study of asphalt materials are also summarized in order to provide the reader with a more comprehensive understanding of the relevant contents and to guide subsequent research.

36 MATERIALS SCIENCE↗

The second-generation Shifted Boundary Method and its numerical analysis

Recently, the Shifted Boundary Method (SBM) was proposed within the class of unfitted (or immersed, or embedded) finite element methods. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we present enhanced variational SBM formulations for the Poisson and Stokes problems with improved flexibility and robustness. These simplified variational forms allow to relax some of the assumptions required by the mathematical proofs of stability and convergence of earlier implementations. First, we show that these new SBM implementations can be proved asymptotically stable and convergent even without the rather restrictive assumption that the inner product between the normals to the true and surrogate boundaries is positive. Second, we show that it is not necessary to introduce a stabilization term involving the tangential derivatives of the solution at Dirichlet boundaries, therefore avoiding the calibration of an additional stabilization parameter. Finally, we prove enhanced L 2 -estimates without the cumbersome assumption – of earlier proofs – that the surrogate domain is convex. Instead we rely on a conventional assumption that the boundary of the true domain is smooth, which can also be replaced by requiring convexity of the true domain. The aforementioned improvements open the way to a more general and efficient implementation of the Shifted Boundary Method, particularly in complex three-dimensional geometries. We complement these theoretical developments with numerical experiments in two and three dimensions.

42 ENGINEERING↗

Linearization errors in discrete goal-oriented error estimation

This paper is concerned with goal-oriented a posteriori error estimation for nonlinear functionals in the context of nonlinear variational problems solved with continuous Galerkin finite element discretizations. A two-level, or discrete, adjoint-based approach for error estimation is considered. The traditional method to derive an error estimate in this context requires linearizing both the nonlinear variational form and the nonlinear functional of interest which introduces linearization errors into the error estimate. In this paper, we investigate these linearization errors. In particular, we develop a novel discrete goal-oriented error estimate that accounts for traditionally neglected nonlinear terms at the expense of greater computational cost. We demonstrate how this error estimate can be used to drive mesh adaptivity. Here, we show that accounting for linearization errors in the error estimate can improve its effectivity for several nonlinear model problems and quantities of interest. We also demonstrate that an adaptive strategy based on the newly proposed estimate can lead to more accurate approximations of the nonlinear functional with fewer degrees of freedom when compared to uniform refinement and traditional adjoint-based approaches.

42 ENGINEERING↗

Finite elements for Matérn-type random fields: Uncertainty in computational mechanics and design optimization

This work highlights an approach for incorporating realistic uncertainties into scientific computing workflows based on finite elements, focusing on prevalent applications in computational mechanics and design optimization. We leverage Matérn-type Gaussian random fields (GRFs) generated using the SPDE method to model aleatoric uncertainties, including environmental influences, variating material properties, and geometric ambiguities. Our focus lies on delivering practical GRF realizations that accurately capture imperfections and variations and understanding how they impact the predictions of computational models as well as the shape and topology of optimized designs. Here we describe a numerical algorithm based on solving a generalized SPDE to sample GRFs on arbitrary meshed domains. The algorithm leverages established techniques and integrates seamlessly with the open-source finite element library MFEM and associated scientific computing workflows, like those found in industrial and national laboratory settings. Our solver scales efficiently for large-scale problems and supports various domain types, including surfaces and embedded manifolds. We showcase its versatility through biomechanics and topology optimization applications, emphasizing the potential to influence these domains. The flexibility and efficiency of SPDE-based GRF generation empowers us to run large-scale optimization problems on 2D and 3D domains, including finding optimized designs on embedded surfaces, and to generate design features and topologies beyond the reach of conventional techniques. Moreover, these capabilities allow us to model and quantify geometric uncertainties on reconstructed submanifolds, such as the interpolated surfaces of cerebral aneurysms provided by postprocessing CT scans. In addition to offering benefits in these specific domains, the proposed techniques transcend specific applications and generalize to arbitrary forward and backward problems in uncertainty quantification involving finite elements.

97 MATHEMATICS AND COMPUTING↗

Micro- to nano-scale areal heterogeneity in pore structure and mineral compositions of a sub-decimeter-sized Eagle Ford Shale

Mineral and organic matter compositions & pore structures of fine-grained shale influence reservoir properties. To improve the understanding of spatial heterogeneity in core-sized samples, methods of microscale X-ray fluorescence (μXRF) mapping, (ultra-) small-angle x-ray scattering [(U)SAXS] and wide-angle X-ray scattering (WAXS) have been used to determine elemental, pore-structure variations at scales up to ~10 cm on two samples prepared at circular and rectangular orientations from a piece of Eagle Ford Shale outcrop in South Texas, USA. In addition, thin section petrography and field emission-scanning electron microscopy observations, X-ray diffraction (XRD), total organic carbon, and pyrolysis were utilized to investigate the potential spatial heterogeneity of pore types, mineral and organic matter compositions for cm-sized samples at both orientations. Overall, the siliceous-carbonate mineral contents in these two samples (8 cm×8 cm×0.08 cm and 5 cm × 8 cm × 0.08 cm, in terms of width×length×thickness) of carbonate-rich Eagle Ford Shale vary between laminations at mm scales. For the circular sample, porosity and specific surface area (SSA) variations range from 0.82 to 3.04% and 1.51 to 14.1 m 2 /g, respectively. For the rectangular sample, values for porosity and SSA vary from 0.93 to 2.50% and 3.95 to 10.8 m 2 /g. By analyzing six selected sub-samples on each of two samples with (U)SAXS and XRD techniques, nm-sized pores are mainly interparticle ones in the higher calcite regions, where the porosity is also relatively lower, while the lower calcite regions consist of both interparticle and intraparticle pore types with higher porosity. Lastly, the μXRF and (U)SAXS mappings are combined to generate porosity distribution maps to provide more insights about sample heterogeneity related to the laminations and fractures at our observational scales.

58 GEOSCIENCES↗

Uniformity and volume-representativity of consecutive submillimetre wafers of shale for SANS and USANS investigations

(Ultra) Small Angle Neutron Scattering (SANS/USANS) is a non-destructive technique suitable for characterization of porous structures from nano- to micro-pore size (1 nm to 20 μm), which has been widely used to characterize accessible and inaccessible porosity in tight unconventional hydrocarbon reservoirs (e.g., shales). Here, in this study, SANS and USANS are used to examine the degree of nano- to micropore uniformity of the rock body on sub-millimetre to centimetre scale. This study presents a test of the hypothesis of volume representativity of pore structure information measured by SANS and USANS for shale wafers (0.22–1 mm) cut from six North American unconventional shale cores of different maturities (New Albany Shale and Marcellus Shale). Experimental reproducibility of USANS is further investigated for two wafers measured using various acquisition times. All samples show that fluctuations of the relative intensity generally decrease from the USANS Q-range (pore size from around 200 nm to 20 μm) to most of the SANS Q-range (pore size from around 5 nm to 700 nm), and peaks in the SANS large-Q region (for pores smaller than about 5 nm). There is a positive correlation between the spatial fluctuation of the relative SANS intensity and the scattering power of shales. In addition, the scattering intensity and its fluctuations is larger for the less mature shale cores. Shale is known to be highly heterogenous and the hypothesis of uniformity does not apply at sub-millimetre scale even for samples that appear homogeneous at visual inspection. This study presents an approach to minimise the variation through sample selection, data interpretation methodology and instrumental reproducibility tests. The best values of the microstructural spatial uniformity for the apparently homogeneous sections of a core (i.e., the SANS and USANS intensity variation) is less than ±10%.

01 COAL, LIGNITE, AND PEAT↗

Answering old questions with new techniques: Understanding performance-limiting factors in transition metal dichalcogenide photoelectrochemical solar cells

In the late 1970s and early 1980s, several research groups reported high-efficiency and stable photoelectrochemical solar cells based on transition metal dichalcogenides (TMDS; e.g., MoS 2 , WSe 2 , and MoSe 2 ) immersed in iodide/tri-iodide electrolytes. A consensus emerged that smooth crystals were necessary for high efficiency based on significant evidence that rough crystals with exposed edge sites produced lower photocurrents and fill factors. However, anecdotal observations in the literature hinted at significant performance variation among apparently smooth crystals with the possibility of highly active “hot” edge sites. Furthermore, this mini-review article is a case study on how spatially resolved photoelectrochemical techniques developed in the 2020s are answering old questions regarding the origin of performance variation in high efficiency n-type TMD|I – ,I 3 – |Pt photoelectrochemical solar cells.

14 SOLAR ENERGY↗

End-gas autoignition fraction and flame propagation rate in laser-ignited primary reference fuel mixtures at elevated temperature and pressure

Knock in spark-ignited (SI) engines is initiated by autoignition of the unburned gasses upstream of spark-ignited, propagating, turbulent premixed flames. Knock propensity of fuel/air mixtures is typically quantified using research octane number (RON), motor octane number (MON), or methane number (MN; for gaseous fuels), which are measured using single-cylinder, variable compression ratio engines. In this study, knock propensity of SI fuels was quantified via observations of end-gas autoignition (EGAI) in unburned gasses upstream of laser-ignited, premixed flames at elevated pressures and temperatures in a rapid compression machine. Stoichiometric primary reference fuel (PRF; n-heptane/isooctane) blends of varying reactivity (50 ≤ PRF ≤ 100) were ignited using an Nd:YAG laser over a range of temperatures and pressures, all in excess of 545 K and 16.1 bar. Laser ignition produced outwardly-propagating premixed flames. High-speed pressure measurements and schlieren images indicated the presence of EGAI. The fraction of the total heat release attributed to EGAI (i.e., EGAI fraction) varied with fuel reactivity (i.e., octane number) and the time-integrated temperature of the end-gas prior to ignition. Flame propagation rates, which were measured using schlieren images, were only weakly correlated with octane number but were affected by turbulence caused by variation in piston timing. Under conditions of low turbulence, measured flame propagation rates approached one-dimensional premixed laminar flame speed computations performed at the same conditions. Experiments were simulated with a three-dimensional CONVERGE™ model using reduced chemical kinetics (121 species, 538 reactions). The simulations accurately captured the measured flame propagation rates, as well as the variation in EGAI fraction with fuel reactivity and time-integrated end-gas temperature. The simulations also revealed low-temperature heat release as well as formaldehyde and hydrogen peroxide formation in the end-gas upstream of the propagating flame, which increased the temperature and degree of chain branching in the end-gas, ultimately leading to EGAI.

42 ENGINEERING↗

Spatially-resolved soot evolution and statistics in high-pressure diesel spray flames using two-color pyrometry

Previous studies of diesel spray flames have focused on details of the sooting behavior mainly during the quasi-steady period, but few have considered the spatially-resolved transient evolution in combination with injection-to-injection variations. In this study, a 500-injection data set is utilized to investigate the temporal evolution of a spray flame during auto-ignition, the premixed burn phase, and the quasi-steady period. Spatially-resolved OH* chemiluminescence data provide ignition delay times and reaction zone locations. Two-color pyrom-etry with a vastly improved optical system is used to quantify spatially-resolved soot evolution and its statistical variations. The ambient thermodynamic con-ditions are slightly below those in modern diesel engines, resulting in longer lift-off lengths and lower overall soot production. Spatially, soot formation in the lift-off region is small, but appreciable soot forms in the jet core and jet head regions, while it oxidizes gradually on the jet periphery. Total soot mass profiles indicate that regions with larger local soot mass take longer to form, but are the first to oxidize. Probability distributions of soot mass in localized regions indicate that a few injections with high soot mass bias the average soot mass data towards higher values. Altogether, results show good agreement with previous studies employing different diagnostic techniques while providing statistical details of transient and localized soot behavior in high pressure diesel spray flames in support of the related modeling efforts.

09 BIOMASS FUELS↗

Molecular scale insight of pore morphology relation with mechanical properties of amorphous silica using ReaxFF

Porous materials are typically heterogeneous and they contain large variations of micro-/nano-pore structures, causing complicated behaviors. In continuum models, most mechanical properties of porous materials are estimated based on porosity, while the variations of micro/nano structures are ignored. That could be problematic as the microscopic heterogeneity may affect the mechanical response of porous materials. Thus, understanding micro/nano heterogeneity impact has been the focus in many scientific and engineering subjects. In the present study, we investigated the effect of nanopore structure (including pore shape and orientation) as well as porosity on mechanical properties of amorphous silica (a-SiO 2 ). The pore sizes in our simulations are comparable to the corresponding ones observed in a-SiO 2 based materials. We found that the existing of nanopores strongly influences Young’s modulus (E) and critical energy release rate (G IC ). These properties decrease with increasing porosity. Importantly, the impact of nanopores was characterized by structural parameters of porous materials. In addition to dependency on porosity, Young’s modulus also was found to vary as a function of potential energy per atom, which highly depends on nanopore shape. Furthermore, critical energy release rate was found to increase with increasing ligament length (also known as pore wall thickness). Finally, the results highlighted the importance of nanopore structures, which must be taken into account when studying fracture mechanisms in porous materials. Based on our findings, it was proposed that mechanical properties of porous materials can be controlled by nano-engineering pore structures.

36 MATERIALS SCIENCE↗

Chemical randomness, lattice distortion and the wide distributions in the atomic level properties in high entropy alloys

High entropy alloys (HEAs) consist of multiple elements present in large proportions that are randomly distributed on a crystal lattice. On the one hand, the presence of multiple elements engenders wide ranges of atomic radii, electronegativities, electron valences and magnetic moments, whereas on the other, the presence of chemical randomness creates unique nearest neighbor environments among the lattice sites. As a result, the symmetry of the energy landscape is broken essentially at each lattice site thereby resulting in highly distorted energy landscapes. At the atomistic level, the lattice distortion has been widely observed in the form of varying bond lengths. At the electronic level, a range of charge transfers result in the charge density distortion. Collectively, the distorted landscapes cause large quantitative variations of the atomic level properties; in this review, we highlight the effect of lattice distortion on point defect energetics, stacking fault energies, and dislocation mobility. Besides the well- known large HEAs phase space, the enormity of the distorted energy landscape that scales with the atomic configurations is a new consideration; understanding this coupling between composition, lattice distortion and properties’ variations thus becomes an exciting but challenging area within the field of HEAs. Furthermore, this coupling is expected to open a new door for materials design, where the materials properties could be tuned via leveraging the lattice distortion, which is essentially absent in dilute/ordered alloys.

36 MATERIALS SCIENCE↗

Myna: Connecting powder bed fusion build data to simulation tools for digital twin applications

Additive manufacturing (AM), as a digital process, can generate a detailed digital thread linking a part’s design and manufacturing to its operational performance. As AM systems advance, an increasing amount of process data is stored in manufacturing databases. In principle, this data can be utilized by simulation-based digital twin approaches, such as real-time process control and asynchronous post-processing guidance. However, few tools currently exist for systematically integrating digital thread data with computational tools. Here, in this study, we propose a software package, called Myna, for connecting data from powder bed fusion processes to simulation tools. The utility of such a platform is demonstrated using build data from the Oak Ridge National Laboratory Manufacturing Demonstration Facility “Peregrine v2023-10” public dataset to automatically configure and run 54 semi-analytical 3DThesis melt pool simulations, 78 numerical Additive FOAM melt pool simulations, and 3 ExaCA microstructure simulations. The simulated, spatially registered microstructures are then compared directly with electron backscatter diffraction characterization of the corresponding as-built part locations. The resulting simulated microstructure showed variation as a function of process parameters, particularly stripe width; however, the experimental data had little variation between the microstructure texture and grain size resulting from different processing conditions. Analysis of the discrepancies suggest that it is possible a two-phase ferritic-austenitic solidification model is needed to accurately predict grain size and texture for certain stainless steel 316L feedstock compositions under powder bed fusion conditions, providing direction for future research. As illustrated here, due to the number and complexity of the simulations involved in AM process-structure–property predictions, automated methods to connect process data and simulations will remain necessary tools for testing hypotheses and implementing digital twin applications.

Knapp, Gerald L. [Oak Ridge National Laboratory (O↗

Stochastic boundary condition effects on supersonic leading edge blowing

The present manuscript studies the robustness of leading edge blowing for active control of the upstream bow shock and overall aerothermal characteristics of a supersonic airfoil subject to variability in the injection conditions. The injection boundary conditions were considered stochastic and their randomness was propagated to the flow topology and aerothermal quantities of interest by means of non-intrusive Polynomial Chaos Expansions. The evaluations needed to build these expansions were carried out via steady two-dimensional Reynolds-Averaged Navier Stokes simulations. Two different injection port sizes were considered, with the smallest promoting a Coanda effect that gives rise to non-symmetric flow topologies. The uncertain boundary conditions lead to significant variations in the flow topology, particularly at the injection exit and the bow shock position, with variations of up to 75% of the nominal value. Here this study identifies the key operating conditions to ensure the optimal application of leading-edge injection for active control of supersonic airfoils.

97 MATHEMATICS AND COMPUTING↗