Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mixed effects”

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 649 records · Page 36

Radiative heat transfer in FLiBe molten salt participating medium in a vertical heated tube under forced and mixed convection laminar flows

The contribution of radiative heat transfer (RHT) to convective heat transfer in a heated vertical pipe with laminar flow and constant wall temperature, for FLiBe molten salt, is investigated computationally. The combined effects of conduction in the fluid, forced convection, buoyancy, temperature-dependent physical properties, and thermal radiation are investigated. The P1 approximation is employed in the discretization of the Radiative Transfer Equation (RTE). The COMSOL Multiphysics software is used to generate the numerical solutions. The theoretical analysis developed here demonstrates that only for intermediate values of the optical thickness, i.e. τ D ~ O 1 , the participating media effects are expected to be important. This analysis is confirmed computationally and a value of τ D ≈ 4 is shown to lead to the highest increase in overall heat transfer behavior. Under forced convection, the Nusselt ratio N u total / N u no - rad reaches a peak value of 1.76 at τ D = 4 and z/D = 200 , and is less than 1.1 for τ D < 0 . 1 and τ D > 60 . Under aiding-flow mixed convection, N u total / N u no - rad reaches a peak of 1.34 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 4 and τ D > 36 . Under opposing-flow mixed convection, Nu total / Nu no - rad reaches a peak of 1.81 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 2 and τ D > 50 . RHT effects on the overall heat transfer are most pronounced in opposing mixed convection, where N u total / N u no - rad is 2.26 observed at z/D = 50 . The sensitivity to wall emissivity is evaluated; in forced convection, the peak N u total / N u no - rad is 1.66 at ε = 0 (reflective wall) and 2.04 at ε = 1 (absorptive wall). Finally, the sensitivity to pipe diameter is also discussed. Overall, for a vertical heated tube, radiative heat transfer effects lead to enhancement of heat transfer by as high as a factor of two, and they depend on the optical thickness of the flow, mixed convection environment ( Gr / R e 2 and direction of flow relative to the gravitational force), surface emissivity, and entrance effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Data and scripts associated with a manuscript modeling microbial regulation of priming effects

This data package is associated with the publication “Modeling Microbial Regulatory Feedback in Organic Matter Decomposition Identifies Copiotrophic Traits as Key Drivers of Positive Priming” published as a preprint on BioRXiv by Ahamed et al. (2026); https://doi.org/10.1101/2024.08.11.607483. The package contains MATLAB scripts and saved simulation outputs used to implement a cybernetic model of microbial regulation during complex organic matter (OM) decomposition governing priming effects. It includes models of (i) single microbial functional groups (copiotrophic or oligotrophic degraders) and (ii) binary consortia composed of degraders and non-degraders with contrasting or common growth traits. Simulation results were generated using Monte Carlo analyses, with randomized key model parameters across a range of environmental mixing fractions of complex and labile OM. The dataset was created to provide a transparent and reusable computational framework for systematically exploring how microbial growth traits, metabolic regulation, and community composition influence OM decomposition dynamics and priming effects. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes the variable definitions. This package includes: (1) annotated MATLAB code implementing the system of ordinary differential equations and cybernetic control laws; (2) saved output files containing data (e.g., biomass, substrates, enzyme levels, priming metrics); and (3) scripts for processing saved outputs and regenerating figures. Specifically, the data package contains three main MATLAB scripts: runPrimingModel.m, runPlotData.m, and runPlotSuppFigS1.m, along with this readme and supporting documentation. Users should begin with runPrimingModel.m, which contains the annotated code implementing the system of ordinary differential equations and cybernetic control laws. This script runs the Monte Carlo simulations of microbial OM decomposition and allows users to modify microbial trait definitions, adjust parameter distributions, or define new community configurations. Simulation outputs are automatically saved as .mat files in the folder named SavedData, which stores all pre-generated results included in this package. The second script, runPlotData.m, reads files from the SavedData folder and processes them to regenerate the figures presented in the manuscript. The third script, runPlotSuppFigS1.m, specifically generates Figure S1 in the Supplementary Material of the manuscript. The package also includes the aforementioned files in non-proprietary .txt format. If users intend to use them, they should first save the files in their respective .m or .mat formats prior to execution in MATLAB.

Biomass concentration↗

Advantages of tandem versus simultaneous operation: The case of isomerization/hydrogenation of terpinolene epoxide to Terpinen-4-ol using a Ni/TiO 2 -SiO 2 bifunctional catalyst

The conversion of terpinolene epoxide to terpinen-4-ol is a desirable chemical transformation for the production of a valuable intermediate. Here this report describes a one-pot, two-step system for the reaction that uses a bifunctional Ni/TiO 2 -SiO 2 catalyst with isomerization and hydrogenation capabilities. The TiO 2 -SiO 2 serves to isomerize the starting epoxide. The Ni component works to hydrogenate the allylic double bond. The tandem system achieves a high 64% yield of terpinen-4-ol. X-ray diffraction analysis showed that TiO 2 -SiO 2 was a mixed binary oxide, and Fourier transform infrared spectroscopy of adsorbed pyridine suggested that acid sites were mainly Lewis acid sites. Electron microscopy showed that the Ni/TiO 2 -SiO 2 catalyst consisted of spherical particles comprising a titanosilicate phase with highly dispersed zero-valent Ni or anatase on the surface. Additionally, there was strong interaction between Ni and the TiO 2 -SiO 2 mixed oxides. Control kinetic experiments demonstrated that the tandem use of the catalysts was more effective than simultaneous operation.

42 ENGINEERING↗

Probing autoionization decay lifetimes of the 4d−16ℓ core-excited states in xenon using attosecond noncollinear four-wave-mixing spectroscopy

The decay of core-excited states is a sensitive probe of autoionization dynamics and correlation effects in many-electron systems, occurring on ultrafast timescales. Xenon, with its dense manifold of autoionizing resonances that can be coupled with near-infrared light, provides a platform to investigate these processes. In this work, the autoionization decay lifetimes of 4d-16ℓ (ℓ = s, p, d, …) core-excited states in xenon atoms are probed with extreme ultraviolet (XUV) attosecond noncollinear four-wave-mixing (FWM) spectroscopy. The 4d{5/2,3/2}-16p XUV-bright states (optically dipole-allowed) exhibit decay lifetimes of ∼6 fs, which is consistent with spectator-type decay. In contrast, the 4d{5/2,3/2}-16s and 4d{5/2,3/2}-16d XUV-dark states (optically dipole forbidden) show longer decay lifetimes of ∼20 fs. Photoionization calculations confirm that all core-hole states with 4d character should decay via spectator channels in ≤6 fs, suggesting that the apparently longer dark-state decay times arise from an alternative mechanism. A few-level simulation of the FWM process shows that the inclusion of a nearby, longer-lived dark state can mimic the experimental FWM signal, suggesting population cycling with a second electronic state with non-4d character. Ab initio calculations support the presence of such multielectron excited states in the 60-70 eV range. These results demonstrate that FWM signals can encode coupled-state dynamics when probing complex systems, highlighting the importance of combining theoretical and experimental approaches to disentangle accurate core-level decay pathways and lifetimes.

Puskar, Nicolette G↗

Mechanics of vitrimer particle compression and fusion under heat press

The compression and fusion of vitrimer particles is a fundamental problem underlying the recycling of vitrimers. Incomplete particle compression can lead to voids in the fused vitrimer, thereby impacting its mechanical property. In this work, we present a two-dimensional finite element model to capture the interplay of multiple complex mechanisms including random particle packing, large deformation, inter-particle contact, thermally activated bulk stress relaxation, and interface healing. Specifically, we focus on understanding two key components of this problem: i) evolution of porosity during compression and fusion of randomly packed vitrimer particles, and ii) effective tensile modulus and strength of fused vitrimer. Using the finite element model, we first show that using smaller particles can reduce the porosity in the fused vitrimer, which only slightly increases the effective modulus, but the increase in tensile strength is more pronounced. In addition, under the same processing conditions, particles with mixed sizes can achieve better densification and higher tensile strength after fusion than uniformly sized particles with the same average radius. Regarding the effects of processing conditions, we discuss how the porosity and effective modulus of the fused vitrimer depend on the processing time, temperature and pressure, and find a general trend consistent with experimental observations, i.e. longer processing time, higher temperature or higher pressure can lead to lower porosity and higher modulus of the fused vitrimer. Finally, the theoretical insights towards the compression and fusion process as well as the simulation framework can be useful in optimizing the powder-based heat press process of vitrimer and its composites.

42 ENGINEERING↗

Strong Sensitivity of Simulated Biomass Burning Aerosol Transport and Radiative Effects Over the South Atlantic to Carbonaceous Aerosol Aging and Particle Density

Biomass burning aerosol (BBA) impacts climate through aerosol‐cloud‐radiation interactions, but models disagree on the sign and magnitude of BBA radiative effects. We quantify the sensitivity of BBA radiative effects and transport to three BBA‐relevant processes and properties: parameterized oxidative aging of organic aerosol (OA), a combined change to black carbon (BC) density and the method for calculating aerosol refractive index, and reduction in OA density. We evaluate Unified Model simulations against two aircraft campaigns from summer 2017 over the Southeast Atlantic. The model generally performs well, such that discrepancies between the observational data sets may sometimes limit the precision of the evaluation. Our newly developed aging parameterization reproduces observed OA:BC mass ratios well and allows modeled OA:BC to decrease with smoke age, but increases bias in aerosol extinction and changes the BBA radiative effect little (+0.12 W m -2 ). We calculate aerosol refractive index using either a volume‐weighted component average or the Maxwell‐Garnett (MG) mixing assumption, which represents BC as small inclusions in a host material. Compared to MG mixing, the volume‐weighted average refractive index and reduced BC density increase aerosol absorption, substantially increasing the total BBA radiative effect (+2.66 W m -2 ) and amount of BBA transported across the ocean through BC self‐lofting. Reducing OA density to better match literature values changes the total BBA radiative effect by −1.96 W m -2 . Changes to direct radiative effects exceed changes to cloud radiative effects. Our findings emphasize the sensitivity of aerosol radiative effects and transport to these processes and properties, which we suggest could be improved in climate models.

Southeast Atlantic↗

The effect on soot and its gas precursors of doping ethylene with 2,2,4,6,6-pentamethyl-heptane in the nitrogen-fuel stream of a laminar non-premixed Planar Mixing Layer Flame (PMLF)

Synthetic Aviation Turbine Fuels (SATFs) are promising for reducing soot emissions from the aviation sector and diversifying Jet Fuel (JF) sources. Accurately predicting the combustion and emissions behavior of SATFs (and other JFs) necessitates robust experimental databases to elucidate the chemistry of long-chain iso-paraffins, which can compose up to two-thirds of SATF blends and whose behavior is considered to be well-represented by that of iso-dodecane isomers. Here, this study characterizes two laminar non-premixed Planar Mixing Layer Flames (PMLFs) with mild soot loads fueled by nitrogen-diluted ethylene, pure and doped with 2,2,4,6,6-pentamethyl-heptane, respectively. The two PMLFs have the same stoichiometric mixture fraction and total hydrocarbon mole fraction in the fuel stream (X F,F =X C2H4,F +X C12H26,F = 0.260), resulting in nearly the same maximum temperature (T max ≈1800 K) and simple identification of the effects of doping. Importantly, any horizontal PMLF cross-section has a self-similar structure that can be modeled as an equivalent One-Dimensional Counterflow Flame (1D-CF) with vanishingly small strain rate (a). The cross-section at a Height Above the Burner (HAB) of 50 mm is characterized in terms of C 0 -C 18 gas species using capillary sampling followed by GC-MS analyses. Laser-Induced Emission Spectroscopy (LIES) quantifies the soot volume fraction (ƒ v ) profiles at HAB=25 and 50 mm where Elastic Laser Light Scattering (E-LLS) is performed to determine the a of the equivalent 1D-CFs and the profile of the E-LLS equivalent diameter ( d 6,3 ) of soot. The substitution of 1500 ppm of ethylene with 2,2,4,6,6-pentamethyl-heptane causes an increase of ≈1.5 in the concentrations of several polycyclic aromatic hydrocarbons and fv. Concurrently, the measured d 6,3 doubles in the oxidizer stream, yet remains the same in the fuel stream, at HAB=50 mm. Instead, at HAB= 25 mm, the iso-dodecane doping does not affect the d 6,3 profile in either stream. The experimental results partially validate the chemical reactions and soot formation kinetic model developed at Lawrence Livermore National Laboratory and provide directions to further improve its predictions.

Iso-dodecane (2,2,4,6,6-Pentamethyl-Heptane)↗

Anomalous interactions between mesons with nonzero spin and glueballs

Topologically nontrivial fluctuations control the anomalous interactions for the η and η ′ pseudoscalar mesons. We consider the anomalous interactions for mesons with higher spin, the heterochiral nonets with J P C = 1 + − and 2 − + . Under the approximation of a dilute gas of instantons, the mixing angle between nonstrange and strange mesons decreases strongly as J increases, and oscillates in sign. Anomalous interactions also open up new, rare decay channels. For glueballs, anomalous interactions indicate that the X ( 2600 ) state is primarily gluonic. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Evidence of Polaron Formation in Halide Perovskites via Carrier Effective Mass Measurements

Carrier effective mass is a central parameter in solid-state physics. It is a measure of the strength of the coupling between a carrier and excitations arising from its surrounding medium, and features prominently in transport and optical calculations. Experimental techniques employed to determine it are steady-state ones, and so are unable to detect any change in the effective mass after a strong perturbation to the system, e.g., strong optical excitation. By combining time-resolved terahertz spectroscopy and transient absorption spectroscopy, on a mixed-cation mixed-halide perovskite thin film, we observe a large and long-lived photoinduced enhancement of the carrier effective mass, and from it deduce a twofold increase of the carrier-phonon coupling constant, giving evidence of polaron formation. Our work demonstrates a new approach to track the strength and ultrafast lifetimes of photoinduced carrier-boson interactions down to picosecond timescales that can be applied to a wide range of solid-state systems.

36 MATERIALS SCIENCE↗

Influences of shock imprinting on mix in a 3D-printed porous media

Mixing of materials in porous media can cause a significant impact on fusion yield as previously demonstrated by the National Ignition Facility (NIF) MARBLE Campaign. Initially, the reactants are separated, with deuterium in the lattice struts and a tritium gas fill in the voids. Lattice parameters such as the strut thickness and relative pitch, provide a control for the mix parameters in the experiment. Los Alamos National Laboratory’s (LANL) BOSQUE project looks to better understand how the mix of the reactants and shell materials impact the fusion burn and resultant yield on various laser platforms. xRAGE’s Eulerian hydrodynamics and adaptive mesh refinement (AMR) provide the unique ability to study the impacts of multiscale features of complex lattice structures. This modeling provides the ability to measure shock front variations as the wave progress’ through a given media. Initial conditions of the lattice are essential to accurately model mix and burn measured by experiment. By varying the initial orientation and densities of these lattice regions the early time dynamics of how the shock is launched into the system is changed and advocates for the study of resulting effects. In this work, we will study the sensitivities of shock effects in varying 3D printed geometric systems and how these shocks alter the structure and mix in the lattice. We will discuss both preliminary experimental results and simulations to help plan and constrain future experiments where we will study the impact of different lattice geometries and lattice bulk densities. This work concludes with the relative impacts of lattice geometries on shock speeds at different bulk densities and the resultant mix widths due to those shock interactions. We see agreement with theory at the higher end of our bulk density study, however, as we approach lower bulk density systems the dynamics of these interactions begin to change.

2 Photon Polymerization↗

Measuring Climate and Water Risk across the Bulk Power System

As climate impacts increase and power systems transition to renewables, planners and operators need insights into climate risks to power generation and infrastructure to ensure reliable decision-making in the short and long-term. We present a standardized, consistent mechanism for utilities and system operators to evaluate the climate- and water-related risks of their current and future grid assets. Using a risk-based approach on the combined outcomes of high-fidelity climate drivers together with water and power system models, we examine the temperature and water availability impacts within the contiguous United States to power system assets at the water basin level in three different time periods and report resulting outcomes on lost capacity across different expansion scenarios and climate models. The results indicate that air temperature has the highest effect on derating. Changes in streamflow do not have a large impact on generation capacity at the national level. Electric sector buildout scenarios each have a unique regional risk profile, depending on the technology mix and total capacity, although risks from high temperatures are significant for both traditional and renewable energy generation. Stakeholders can use this approach to monitor effects of generation capacity losses and potential impacts as climate, generation mix, and infrastructure change.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Giant Isotope Effect of Thermal Conductivity in Silicon Nanowires

Isotopically purified semiconductors potentially dissipate heat better than their natural, isotopically mixed counterparts as they have higher thermal conductivity (κ). But the benefit is low for Si at room temperature, amounting to only ~10% higher κ for bulk 28 Si than for bulk natural Si ( nat Si). Here we show that in stark contrast to this bulk behavior, 28 Si (99.92% enriched) nanowires have up to 150% higher κ than nat Si nanowires with similar diameters and surface morphology. Using a first-principles phonon dispersion model, this giant isotope effect is attributed to a mutual enhancement of isotope scattering and surface scattering of phonons in nat Si nanowires, correlated via transmission of phonons to the native amorphous SiO 2 shell. The Letter discovers the strongest isotope effect of κ at room temperature among all materials reported to date and inspires potential applications of isotopically enriched semiconductors in microelectronics.

36 MATERIALS SCIENCE↗

Aqueous Phase Aldol Condensation of Formaldehyde and Acetone on Anatase TiO2(101) Surface: A Theoretical Investigation

A mechanistic understanding of catalytic reactions at solid-liquid interface is limited both experimentally and theoretically but attracts much interest. Using density functional theory calculations (DFT) and ab initio molecular dynamics (AIMD) simulations, we investigated the effect of liquid water on a-H abstraction, C-C coupling, and dehydration steps of aldol condensation of formaldehyde and acetone on an anatase TiO2(101) surface. The existence of the aqueous phase lowered the Gibbs energy of activation of dehydration step pronouncedly from 187 to 74 kJ/mol through proton transfer mechanism, making the hydrogenation pathway more favorable in the aqueous phase. In contrast, a mixed route prevails in the vapor phase. This work provides insights into the effect of the bulk water through a proton transfer mechanism on the dehydrogenation, C-C coupling, and dehydration

anatase TiO2, formaldehyde and acetone, aldol cond↗

Detecting Atmospheric Molecules of Nontransiting Temperate Terrestrial Exoplanets Using High-resolution Spectroscopy in the Mid-infrared Domain

Motivated by the development of high-dispersion spectrographs in the mid-infrared (MIR) range, we study their application to the atmospheric characterization of nearby nontransiting temperate terrestrial planets around M-type stars. We examine the detectability of CO{sub 2}, H{sub 2}O, N{sub 2}O, and O{sub 3 }features in high-resolution planetary thermal emission spectra at 12–18 μm assuming an Earth-like profile and a simplified thermal structure. The molecular line width of such planets can be comparable to or broader than the Doppler shift due to the planetary orbital motion. Given the likely difficulty in knowing the high-resolution MIR spectrum of the host star with sufficient accuracy, we propose observing the target system at two quadrature phases and extracting the differential spectra as the planetary signal. In this case, the signals can be substantially suppressed compared with the case where the host star spectrum is perfectly known, as some parts of the spectral features do not remain in the differential spectra. Despite this self-subtraction, the CO{sub 2} and H{sub 2}O features of nearby (≲5 pc) systems with mid-/late-M host stars would be feasible with a 6.5 m class cryogenic space telescope, and orbital inclination could also be constrained for some of them. For CO{sub 2} and N{sub 2}O in a 1 bar Earth-like atmosphere, this method would be sensitive when the mixing ratio is 1–10{sup 3} ppm. The detectability of molecules except O{sub 3} is not significantly improved when the spectral resolution is higher than R≳10,000, although the constraint on the orbital inclination is improved. This study provides some benchmark cases useful for assessing the value of MIR high-resolution spectroscopy in terms of characterization of potentially habitable planets.

47 OTHER INSTRUMENTATION↗

Effects of diffusion barriers on reaction wave stability in Co/Al reactive multilayers

Bimetallic, reactive multilayers are uniformly structured materials composed of alternating sputter-deposited layers that may be ignited to produce self-propagating mixing and formation reactions. These nanolaminates are most commonly used as rapid-release heat sources. The specific chemical composition at each metal/metal interface determines the rate of mass transport in a mixing and formation reaction. The inclusion of engineered diffusion barriers at each interface will not only inhibit solid-state mixing but also may impede the self-propagating reactions by introducing instabilities to wavefront morphology. Here, this work examines the effect of adding diffusion barriers on the propagation of reaction waves in Co/Al multilayers. The Co/Al system has been shown to exhibit a reaction propagation instability that is dependent on the bilayer thickness, which allows for the occurrence of unstable modes in otherwise stable designs from the inclusion of diffusion barriers. Based on the known stability criteria in the Co/Al multilayer system, the way in which the inclusion of diffusion barriers changes a multilayer's heat of reaction, thermal conductivity, and material mixing mechanisms can be determined. These factors, in aggregate, lead to changes in the wavefront velocity and stability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Improving Conversion of Captured CO 2 on Mixed-Oxide Dual-Function Materials by Low-Concentration Metal Promotion

Reactive carbon capture (RCC) is a promising nascent approach to harness the untapped resource of waste CO 2 to produce fuels and platform chemicals. However, process economics relies on the design of dual-function materials (DFMs) that exhibit maximal conversion of captured CO 2 and high selectivity toward valuable carbon products (e.g., CO). Herein, we report the use of low-concentration (ca. 0.5 wt %) metal promotors (Pt and Au) to improve the conversion of captured CO 2 on a Zn-Al mixed-oxide DFM by >30%, while also more than doubling per-cycle CO yield due to enhanced retention of captured CO 2 and H 2 activation. In situ DRIFTS of bound CO 2 on the DFMs showed distinctions in surface carbonate formation for the Pt- and Au-modified materials that were correlated with product evolution during the reactive desorption step of RCC. Particularly, Pt afforded superior H 2 activation but formed irreversibly-bound bicarbonates, compromising overall productivity, while the Na/Au/ZA DFM predominantly adsorbed CO 2 on Na sites, resulting in polydentate carbonates that are retained until higher temperatures before conversion to CO. In conclusion, these results indicate that adding a low concentration of oxidation-resistant noble metals, especially Au, is an effective strategy for improving key performance metrics of Zn-Al mixed-oxide DFMs to further advance the development of scalable RCC technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spin-singlet topological superconductivity in the attractive Rashba-Hubbard model

Fully gapped, spin-singlet superconductors with antisymmetric spin-orbit coupling in a Zeeman magnetic field provide a promising route to realize superconducting states with non-Abelian topological order and therefore fault-tolerant quantum computation. Here we use a quantum Monte Carlo dynamical cluster approximation to study the superconducting properties of a doped two-dimensional attractive Hubbard model with Rashba spin-orbit coupling in a Zeeman magnetic field. We generally find that the Rashba coupling has a beneficial effect towards s-wave superconductivity. In the presence of a finite Zeeman field, when superconductivity is suppressed by Pauli pair breaking, the Rashba coupling counteracts the spin imbalance created by the Zeeman field by mixing the spins, and thus restores superconductivity at finite temperatures. Further, we show that this favorable effect of the spin-orbit coupling is traced to a spin-flip driven enhancement of the amplitude for the propagation of a pair of electrons in time-reversed states. Moreover, by inspecting the Fermi surface of the interacting model, we show that for sufficiently large Rashba coupling and Zeeman field, the superconducting state is expected to be topologically nontrivial.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Latent variable Gaussian process models: A rank‐based analysis and an alternative approach

Abstract Gaussian process (GP) models have been extended to emulate expensive computer simulations with both qualitative/categorical and quantitative/continuous variables. Latent variable (LV) GP models, which have been recently developed to map each qualitative variable to some underlying numerical LVs, have strong physics‐based justification and have achieved promising performance. Two versions use LVs in Cartesian (LV‐Car) space and hyperspherical (LV‐sph) space, respectively. Despite their success, the effects of these different LV structures are still poorly understood. This article illuminates this issue with two contributions. First, we develop a theorem on the effect of the ranks of the qualitative factor correlation matrices of mixed‐variable GP models, from which we conclude that the LV‐sph model restricts the interactions between the input variables and thus restricts the types of response surface data with which the model can be consistent. Second, following a rank‐based perspective like in the theorem, we propose a new alternative model named LV‐mix that combines the LV‐based correlation structures from both LV‐Car and LV‐sph models to achieve better model flexibility than them. Through extensive case studies, we show that LV‐mix achieves higher average accuracy compared with the existing two.

Tao, Siyu↗