Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “State functions”

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,027 records · Page 57

Au assisted smooth ultrathin epitaxial ZnO film grown by pulsed laser deposition on sapphire(0001)

High quality ZnO has applications in photonics and electronics. Literature has shown that continuous, but rough ZnO film can be grown on sapphire substrate by pulsed laser deposition (PLD). In this work, ZnO and Au were co-deposited by PLD on sapphire(0001) substrates from a single ZnO target overlaid with an Au strip that serves as a catalyst. The reflection high-energy electron diffraction (RHEED) patterns show paired double stripes from sample of ZnO-Au where the difference of in-plane lattice parameters between ZnO and Au is 11.4 %. The two-dimensional (2D) reciprocal space maps constructed from azimuthal RHEED (ARHEED) patterns reveal parallel epitaxial relationship between ZnO and Au in the continuous ZnO-Au film on sapphire substrate. Transmission electron microscopy (TEM) cross section images and energy dispersive spectroscopy maps show a minute amount of Au nanocrystals distributed in the entire ultrathin ZnO film. Atomic force microscopy shows the root-mean-square roughness of the surface of ultrathin ZnO film with Au is in the sub-nm range. The ARHEED and TEM results show that a couple of atomic percentages of Au co-deposited with ZnO catalyzes the growth of smooth ultrathin epitaxial ZnO film.

36 MATERIALS SCIENCE↗

Influence of ethylenediaminetetraacetic acid on the long-term oxidation state distribution of plutonium

Here, spectrophotometry was used to study the effect of EDTA on plutonium oxidation state distribution as a function of time, pH, and ligand-to-metal ratio (L/M) under anoxic conditions. Novel Pu(V)-EDTA absorption bands were identified at 571, 993, 1105, and 1150 nm with molar absorption coefficients of 15 ± 1, 6 ± 1, 10 ± 1, and 10 ± 1 cm –1 M –1 , respectively. Pu(V)-EDTA spectral changes occurred at L/M < 1, indicating only Pu V O 2 (EDTA) 3- formed with logK = 3.6 ± 0.3. Time-resolved experiments showed EDTA drastically increased the Pu(V/VI) reduction rate, which we propose is driven by amine lone-pair electron donation and the oxidative decarboxylation of EDTA. Oxidation of Pu(III)-EDTA to Pu(IV)-EDTA occurred on a slower time scale (110–237 days) than previously reported (<15 min) and is hypothesized to be radiolysis driven. Pu(V/VI)-EDTA and Pu(III)-EDTA both approached Pu(IV)-EDTA stabilization over time, yet Pu(V/VI)-EDTA solubility data was ≥ 1.0 log 10 units higher than predicted by Pu(IV)-EDTA solubility models, indicating that current thermodynamic models are incomplete. Ultimately, the data show EDTA preferentially stabilizes Pu(IV) over time regardless of initial oxidation state, but Pu(V)-EDTA can persist under environmentally-relevant conditions, emphasizing the need to continue investigating redox reactions, speciation, and behavior of these complexes to support the transuranic waste disposal and surface remediation/containment efforts.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Tropical Intraseasonal Air-Sea Exchanges during the 1997 Pacific Warming

The Madden Julian Oscillations (MJO) and associated westerly wind (WW) events account for much of the tropical intraseasonal variability (TISV). The TISV has been suggested as an important stochastic forcing that may be one of the underlying causes for the observed irregularities of the El Nino-Southern Oscillation (ENSO). Recent observational studies and theories of interannual to interdecadal-scale variability suggest that ENSO may arise from different mechanisms depending on the basic states. The Pacific warming event of 1997, being associated with a period of strong MJO and WW events, serves as a natural experiment for studying the possible role of TISV in triggering an ENSO event. We have performed a combined statistical and composite analysis of surface WW events based on the assimilated surface wind and sea level pressure for the period of 1980-1993, the SSM/I wind for the period of 1988-1997, and OLR. Results indicates that extratropical forcing contribute significantly to the evolution of MJO and establishment of WW events over the Pacific warm pool. Following the major WW events, there appeared an eastward extension of equatorial warm SST anomalies from the western Pacific warm pool. Such tropical-extratropical interaction is particularly clear in the winter of 96-97 that leads to the recent warming event in 1997/98. From the above discussion, our current study on this subject is based on the hypothesis that 1) there is an enhanced air-sea interaction associated with TISV and the northerly surges from the extratropics in the initial phase of the 97/98 warming event, and 2) the relevant mechanisms are functions of the basic state of the coupled system (in terms of SST distribution and atmospheric mean circulation) that varies at the interannual and interdecadal time scale. We are analyzing the space-time structure of the northerly surges, their association with air-sea fluxes and upper ocean responses during the period of September 1996 to June 1997. The estimate of daily values of latent heat fluxes is based on NSCAT wind, SST, and ECMWF surface air temperature and SSM/I water vapor data (Chou et al. 1997). To understand the relevant mechanisms, we will analyze the origin of the northerly surges in terms of atmospheric instability associated with the extratropical circulation, and the mutual influence between the tropical heating and the extratropical circulation. In this meeting, we will report the analysis addressing the first part of the above hypothesis.

Sui, C.-H.↗

The promising role of proteomes and metabolomes in defining the single-cell landscapes of plants

The plant community has a strong track-record of RNA sequencing technology deployment, which combined with the recent advent of spatial platforms (e.g., 10x genomics), has resulted in an explosion of outstanding single cell and nuclei datasets that can be put in an in situ context within tissues (e.g., a cell atlas)1. In the genomics era, application of proteomics technologies in the plant sciences has always trailed behind that of RNA sequencing technologies, largely due to accessibility, ease-of-use and access to expertise along with depth of analysis benefits. On the other hand, the use of early analytical tools for characterizing small molecules (metabolites) from plant systems predates nucleic acid sequencing and proteomics analysis2, as the search for plant-based natural products has played a significant role in improving human health throughout history. However, the employment of proteomics and metabolomics assays for characterizing plant cell processes now remains significantly behind transcriptional approaches, even though both provide a direct functional readout of cell states and phenotypes.

Anderton, Christopher R. [BATTELLE (PACIFIC NW LAB↗

Shake loss intensities in x-ray photoelectron spectroscopy: Theory, experiment, and atomic composition accuracy for MgO and related compounds

The relative intensities of XPS core levels, scaled by their photoionization cross sections, are regularly used to determine sample atomic composition. Cross sections, however, give the intensity to all possible final states for the core ionizations, not just to the main peak. This includes all intrinsic satellite structure (shake states and, for open shell systems, the different ionic multiplets). In practice, for solids, this is usually experimentally impossible to determine accurately because such a satellite structure sits on the inelastically scattered electron background and cannot be easily separated. Therefore, usually, only the intensity of the main peak is used. This limits the ultimate possible accuracy of XPS composition determination. The purpose of the present paper is to examine the contributions that a theoretical analysis of losses of intensity can make to improve quantitation. For an MgO single crystal, we show that the correct stoichiometry of 1:1 can be recovered using the theoretical analysis of the experimental MgO peak ratio intensities. For materials with a sufficient bandgap for the XPS main peaks to be separated from the scattered background, the intensity of main peaks can often be accurately determined. Thus, if one uses theory to calculate that fraction of the total intensity lost from a main peak into all its satellite structure, the intensity of just main peaks could then be used to more accurately determine relative atom % composition. This work tests this approach using a single crystal MgO (50% Mg, 50% O) standard. Ab initio electronic structure theory of representative MgO clusters is used to determine Hartree–Fock wave functions for the ground state and final ionized states corresponding to the main Mg 2p and O1s XPS peaks of the oxide. The sudden approximation, SA, is used to determine the fractional losses from these main peaks to shake satellites, which is found to be greater for O1s than Mg2p. This results in predicted “apparent composition” for stoichiometric MgO of 55.2% Mg, 44.8% O instead of the true 50% Mg, 50% O. Equivalent theory for CaO results in a predicted apparent Ca value of 53.4%. Experimentally, using Mg2s or 2p intensity ratio to O1s, we find values between 52.2% and 56.0% Mg using two crystals and four different instrument electron pass energies. The average value of the measurements is 54.5% Mg when corrected for the presence of an adventitious carbon overlayer and slight surface hydroxide. Though this agreement with theory may be somewhat fortuitous, given the potential experimental errors, which are fully discussed, it is similar to that in our earlier study on LiF. We also present preliminary experimental data on Mg(OH) 2 and MgSO 4 , which show a similar trend of apparently higher than 50% Mg, but we have no theory values. We are not yet able to experimentally test for validation of the difference between apparent composition for MgO (55.2% Mg) and CaO (53.4% Ca), owing to significant carbonate formation at the surface of the single crystal CaO. Here, an important conclusion is that the theoretical determination of shake losses, obtained with ab initio wavefunctions and the SA, is likely to be a useful way to calibrate the accuracy and reliability of compositions obtained from XPS intensities and merits further study.

47 OTHER INSTRUMENTATION↗

The effect of copper additives on hybrid Zr-based chemical conversion coating morphology and chemical compositions

In the realm of corrosion protection, Zr-based conversion coatings offer an environmentally friendly, chromate-free alternative to conventional coating. This study uses advanced X-ray, electron microscopy and electrochemical testing techniques to better understand the impact of varying Cu²⁺ ion concentrations on the characteristics of Zr-based coatings on Fe substrates. Our findings demonstrate that within the tested conditions, higher Cu²⁺ additive concentrations at 40 ppm enhance surface characteristics, increasing stability towards anti-corrosion capability, particularly under NaCl treatment. Conversely, at a lower Cu²⁺ concentration of 20 ppm, coatings exhibited more significant dissolution of Cu clusters and increased vulnerability to chloride-induced degradation. X-ray photoelectron spectroscopy and synchrotron X-ray fluorescence (XRF) and X-ray absorption near edge structure (XANES) spectroscopy analyses revealed Cu²⁺ formation in Cu20, while Cu40 retained its metallic state (Cu⁰) with a slower reaction rate. Cu20 offers some protection but lacks durability, whereas in Cu40, lower I corr , and enhance corrosion resistance, making it ideal for protective coatings in salt-water environments. This study underscores the importance of balancing Cu²⁺ ions concentration in the coating solution to optimize performance, highlighting the role of Cu in enhancing both surface properties and long-term stability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Chemical and morphological evolution of hybrid conversion coatings in low-Earth orbit space environment

Understanding how protective coatings respond to the harsh low-Earth orbit (LEO) environment is essential for ensuring the safety, longevity, and cost-effectiveness of spacecraft. In particular, identifying environmentally friendly, non-chromate alternatives that can maintain performance under such conditions has both technological and regulatory significance. This study investigates the environmental stability of zirconium-based hybrid conversion coatings with Cu additives (Cu10 and Cu20) applied to cold-rolled steel, tested in the Materials International Space Station Experiment (MISSE) outside the International Space Station (ISS). Chemical and morphological analyses were carried out using a combination of electron microscopy and X-ray spectroscopy techniques, including scanning electron microscopy (SEM), scanning transmission electron microscopy with energy-dispersive X-ray spectroscopy (STEM-EDS), X-ray photoelectron spectroscopy (XPS), and X-ray absorption near-edge structure (XANES) spectroscopy. After exposure outside the ISS, all coatings remained structurally intact, with all exhibiting a uniform Zr-rich matrix and embedded Cu-rich clusters, while a thin Si-rich surface layer developed from interaction with space environments. Depth-resolved XPS showed a layered structure with CuO on the surface, Cu 2 O, and partial Zr(IV) reduction near Cu-rich sites, evidence of Atomic Oxygen (AO)-driven surface oxidation. These results demonstrate that Cu–Zr coatings maintain their chemical integrity and microstructure in harsh space environments, offering a non-chromate alternative for long-term aerospace protection. These insights provide valuable guidance for developing next-generation protective coatings that combine environmental sustainability with the reliability required for future aerospace and orbital applications.

36 MATERIALS SCIENCE↗

Chiral flat-band optical cavity with atomically thin mirrors

A fundamental requirement for photonic technologies is the ability to control the confinement and propagation of light. Widely used platforms include two-dimensional (2D) optical microcavities in which electromagnetic waves are confined in either metallic or distributed Bragg reflectors. Recently, transition metal dichalcogenides hosting tightly bound excitons with high optical quality have emerged as promising atomically thin mirrors. In this work, we propose and experimentally demonstrate a subwavelength 2D nanocavity using two atomically thin mirrors with degenerate resonances. Angle-resolved measurements show a flat band, which sets this system apart from conventional photonic cavities. We demonstrate how the excitonic nature of the mirrors enables the formation of chiral and tunable optical modes upon the application of an external magnetic field. Moreover, we show the electrical tunability of the confined mode. Our work demonstrates a mechanism for confining light with high-quality excitonic materials, opening perspectives for spin-photon interfaces, and chiral cavity electrodynamics.

36 MATERIALS SCIENCE↗

Sub-wavelength optical lattice in 2D materials

Recently, light-matter interaction has been vastly expanded as a control tool for inducing and enhancing many emergent nonequilibrium phenomena. However, conventional schemes for exploring such light-induced phenomena rely on uniform and diffraction-limited free-space optics, which limits the spatial resolution and the efficiency of light-matter interaction. Here, we overcome these challenges using metasurface plasmon polaritons (MPPs) to form a sub-wavelength optical lattice. Specifically, we report a “nonlocal” pump-probe scheme where MPPs are excited to induce a spatially modulated AC Stark shift for excitons in a monolayer of MoSe 2 , several microns away from the illumination spot. We identify nearly two orders of magnitude reduction for the required modulation power compared to the free-space optical illumination counterpart. Moreover, we demonstrate a broadening of the excitons’ linewidth as a robust signature of MPP-induced periodic sub-diffraction modulation. Our results will allow exploring power-efficient light-induced lattice phenomena below the diffraction limit in active chip-compatible MPP architectures.

36 MATERIALS SCIENCE↗

Pressure-induced structural phase transitions in CrSBr

There is growing interest in combining chemical complexity with external stimuli like pressure, field, and light for property control in van der Waals solids. This is because extreme conditions trigger the development of new states of matter and functionality. In this work, we bring together synchrotron-based infrared absorption, Raman scattering, and diamond anvil cell techniques with first-principles calculations of the lattice dynamics and energy landscape to reveal the series of structural phase transitions in CrSBr. By tracking how the phonons change under pressure, we uncover a remarkable chain of complex symmetry modifications, interlayer interactions, and chemical reactions. A group-subgroup analysis suggests that CrSBr undergoes an orthorhombic Pmmn → monoclinic P2/m transition at 7.6 GPa, and based upon a comparison with model oxychlorides like FeOCl and CrOCl, we propose that changes in the pendant halide groups drive the system to a P2 1 /m-like space group above 15.3 GPa. Compression above 20.2 GPa is irreversible, resulting in the formation of an entirely new compound that is metastable for months. This work opens the door to the use of pressure and possibly strain to control the properties of CrSBr.

Microscopy↗

A time-parallel multiple-shooting method for large-scale quantum optimal control

Quantum optimal control plays a crucial role in quantum computing by providing the interface between compiler and hardware. Solving the optimal control problem is particularly challenging for multi-qubit gates, due to the exponential growth in computational complexity with the system's dimensionality and the deterioration of optimization convergence. To ameliorate the computational complexity of time-integration, this paper introduces a multiple-shooting approach in which the time domain is divided into multiple windows and the intermediate states at window boundaries are treated as additional optimization variables. Further, this enables parallel computation of state evolution across time-windows, significantly accelerating objective function and gradient evaluations. Since the initial state matrix in each window is only guaranteed to be unitary upon convergence of the optimization algorithm, the conventional gate trace infidelity is replaced by a generalized infidelity that is convex for non-unitary state matrices. Continuity of the state across window boundaries is enforced by equality constraints. A quadratic penalty optimization method is used to solve the constrained optimal control problem, and an efficient adjoint technique is employed to calculate the gradients in each iteration. We demonstrate the effectiveness of the proposed method through numerical experiments on quantum Fourier transform gates in systems with 2, 3, and 4 qubits, noting a speedup of 80x for evaluating the gradient in the 4-qubit case, highlighting the method's potential for optimizing control pulses in multi-qubit quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Static and dynamic spin properties in the quantum triangular lattice antiferromagnet Ag 2 CoO 2

In Ag 2 CoO 2 , Co forms triangular lattice layers, which are separated by the metallic (Ag 2 ) block. The magnetic susceptibility and heat capacity measurements show that this material exhibits an antiferromagnetic transition at T N =17.5 K and the Weiss temperature (T Θ ) and the effective moment are -274 K and 1.62μ B , respectively, indicating that the Co ion carries spin (S) 1/2 and has a strongly frustrated state with T Θ /T N =15.7. A density functional theory calculation confirmed that the valence state of the Co ions is 2+ and the low-spin state with S=1/2 is realized at reduced on-site Coulomb interaction on Co. We performed elastic and inelastic neutron scattering experiments in a powder sample of Ag 2 CoO 2 . Although no noticeable magnetic Bragg peaks were observed below T N , distinct magnetic excitations were observed in the inelastic neutron scattering experiments. The excitations are consistent with those expected for the S=1/2 Heisenberg triangular lattice antiferromagnet. These results suggest that the ordered moment is reduced due to the quantum effect, which explains the absence of the magnetic Bragg peaks. Our results thus suggest that Ag 2 CoO 2 is a good candidate to realize a quantum Heisenberg triangular lattice antiferromagnet.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Thermodynamically consistent forward and reverse degrees of rate control in reversible reactions

The net rate of a composite reaction is the difference between the forward and reverse reaction rates, which are kinetically distinct despite sharing elementary reaction steps and therefore have different rate-controlling transition states and species. Thus, degrees of rate control defined to identify rate-controlling transition states and species for the net rate confound contributions from the forward and reverse reactions. Herein, the forward and reverse degrees of rate control are defined to independently quantify the extent that species and transition states control the forward and reverse rates in reversible reactions. These degrees of rate control are defined as the relative change in the forward and reverse reaction rates per k B T decrease in the standard-state molecular free energies of transition states and species, and they are related to experimentally measurable quantities such as steady-state fractional coverages, reaction orders, and energies and entropies of activation of the forward and reverse reaction rates. The forward/reverse degrees of rate control represent stoichiometric coefficients for species and transition states in an equilibrium between the apparent initial states and the apparent transition states of the apparent rate-controlling steps of forward and reverse overall reactions. At equilibrium, the apparent transition states for the forward reaction and reverse reaction converge, and thus the forward and reverse rate-controlling steps combine to form a single apparent rate-controlling step at equilibrium. Here, this apparent rate-controlling step is comprised of an apparent initial state, transition state, and final state, where the apparent final state of the forward reaction is the apparent initial state of the reverse reaction. The apparent rate-controlling step behaves identically to an elementary step reaction at equilibrium with a pseudo-mass-action rate function given by the transition-state-theory (TST) form rate function (Foley and Bhan, 2020) with a stoichiometric number equal to the affinity-averaged stoichiometric number, $\bar{σ}$ The ratio of the forward and reverse TST-form rate functions is identical to the overall thermodynamic equilibrium relation, consistent with the principles of microscopic reversibility and detailed balance at equilibrium.

42 ENGINEERING↗

Particle Markov Chain Monte Carlo Approach to Inference in Transient Surface Kinetics

Here, in this work, we develop a novel Bayesian approach to study the adsorption and desorption of CO onto a Pd(111) surface, a process of great importance in natural sciences. The motivation for this work comes from the recent availability of time-resolved infrared spectroscopy data and the need for model interpretability and uncertainty quantification in chemical processes. The objective is to learn the relevant parameters that characterize the process: coverage with time, rate constants, activation energies, and pre-exponential factors. Our approach consists of three main schemes: (i) a problem design and probabilistic model for the whole system, (ii) a particle Markov chain Monte Carlo sampler to learn the hidden coverages and rate constant parameters, and (iii) two Bayesian formulations to infer the activation energies and pre-exponential factors. The flexibility of the Bayesian framework allows for uncertainty quantification where possible and integration of mathematical constraints in the model to reflect the system physically. We found that our results for the activation energies and pre-exponential factor are in agreement with those reported in the experimental literature, independently, and we provide discussions on the advantages and disadvantages as well as applicability to other systems.

36 MATERIALS SCIENCE↗

Quantifying dispersity in size and shape of nanoparticles from small-angle scattering data using machine learning based CREASE

Here, we use machine learning (ML) enhanced computational reverse engineering analysis of scattering experiments (CREASE) to interpret small-angle X-ray scattering (SAXS) data obtained from a system of nanoparticles without a priori knowledge of their exact shapes (e.g. spheres or ellipsoids), sizes (0.5–50 nm) and distributions. The SAXS measurements yielded three categories of scattering profiles exhibiting 'strong', 'weak' and 'no' features. Diminishing features (e.g. broadening or disappearing peaks) in scattering profiles have always been attributed to the presence of significant dispersity in the system. Such featureless SAXS data are not suitable for traditional analysis using analytical models. If one were to fit a relevant analytical model (e.g. the lmfit analytical model for polydisperse spheres) to these 'weak' and 'no' SAXS profiles from our nanoparticle systems, one would obtain non-unique interpretations of the data. Relying on electron microscopy to identify the distributions of nanoparticle shapes and sizes is also unfeasible, especially in high-throughput synthesis and characterization loops. In such situations, to identify the distributions of particle sizes and shapes that could be present in the sample, one must rely on methods like ML-CREASE to interpret the data quickly and output all relevant interpretations about the structure present in the system. The ML-CREASE optimization loop takes the experimental scattering profile as input and outputs multiple candidate solutions whose computed scattering profiles match the SAXS profile input. The ML-CREASE method outputs distributions of relevant structural features, such as the volume fraction of the nanoparticles in the system and the mean and standard deviation of the particle size and aspect ratio, assuming a type of distribution (e.g. normal, log-normal) for size and aspect ratio. We find that, for the SAXS profiles analyzed here, accounting for the shape dispersity along with size dispersity of the nanoparticles using ML-CREASE improved the match between the computed scattering profiles and input experimental profiles.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Optomechanical Tuning of Second Harmonic Generation Anisotropy in Janus MoSSe/MoS 2 Heterostructures

Symmetry breaking in van der Waals materials enables the realization of quantum states and advanced device functionalities. Janus transition-metal dichalcogenides (TMDs) exhibit distinctive nonlinear optical properties due to their broken out-of-plane mirror symmetry. However, the dynamic control of second harmonic generation (SHG) anisotropy and resonance behavior via optical excitation remains elusive. Here, in this work, we investigate the SHG response of Janus MoSSe/MoS 2 heterostructures with 2H and 3R stackings. We can tune the SHG response by varying the incident photon wavelength from 800 to 1000 nm, which shows a resonance-dependent enhancement in intensity and a deviation from 6-fold symmetry, indicating wavelength-dependent anisotropy. The ratio between maximum and minimum intensity in the armchair directions, associated with the SHG anisotropy, reaches a value of 1.73 at the excitation wavelength of 1000 nm. Group theory analysis and first-principles calculations reveal that the observed anisotropy arises from optically induced strain. Our findings highlight the role of symmetry breaking and optical resonance contributing to the optomechanical tuning of SHG anisotropy, offering opportunities for developing Janus TMD-based photonic devices for frequency conversion, light generation, and optical switching.

Janus transition metal dichalcogenides↗

Deciphering electrocatalysts with multimodal operando approaches

Here, the electrocatalytic processes of a copper catalyst during nitrate electroreduction are unveiled by correlated operando microscopy and spectroscopy. Catalysts are vital to the modern chemical industry, yet their development has largely relied on trial-and-error approaches. Optimizing catalysts requires a fundamental understanding of their structure–chemical property behaviour under operational conditions. However, operando characterization remains challenging, especially for reactions occurring in the complex and dynamic liquid environments of electrochemical systems. Over the past decade, the rapid development of in situ and operando environmental transmission electron microscopy (ETEM) and microelectromechanical system (MEMS)-based closed-cell holders, enabling environmental studies within the vacuum environment of transmission electron microscopy (TEM), has substantively advanced understanding of heterogeneous catalysis, notably for gas-phase reactions. In situ ETEM enables the direct observation of the catalyst evolution in terms of structure, morphology and chemical state at the nano-to-atomic scale, providing insights into their correlation with catalytic performance.

36 MATERIALS SCIENCE↗

Variational Methods in Sensitivity Analysis and Optimization for Aerodynamic Applications

Variational methods (VM) sensitivity analysis, which is the continuous alternative to the discrete sensitivity analysis, is employed to derive the costate (adjoint) equations, the transversality conditions, and the functional sensitivity derivatives. In the derivation of the sensitivity equations, the variational methods use the generalized calculus of variations, in which the variable boundary is considered as the design function. The converged solution of the state equations together with the converged solution of the costate equations are integrated along the domain boundary to uniquely determine the functional sensitivity derivatives with respect to the design function. The determination of the sensitivity derivatives of the performance index or functional entails the coupled solutions of the state and costate equations. As the stable and converged numerical solution of the costate equations with their boundary conditions are a priori unknown, numerical stability analysis is performed on both the state and costate equations. Thereafter, based on the amplification factors obtained by solving the generalized eigenvalue equations, the stability behavior of the costate equations is discussed and compared with the state (Euler) equations. The stability analysis of the costate equations suggests that the converged and stable solution of the costate equation is possible only if the computational domain of the costate equations is transformed to take into account the reverse flow nature of the costate equations. The application of the variational methods to aerodynamic shape optimization problems is demonstrated for internal flow problems at supersonic Mach number range. The study shows, that while maintaining the accuracy of the functional sensitivity derivatives within the reasonable range for engineering prediction purposes, the variational methods show a substantial gain in computational efficiency, i.e., computer time and memory, when compared with the finite difference sensitivity analysis.

Ibrahim, A. H.↗