Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix models”

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 883 records · Page 49

Thermodynamic and Kinetic Modeling of Co-utilization of Glucose and Xylose for 2,3-BDO Production by Zymomonas mobilis

Prior engineering of the ethanologen Zymomonas mobilis has enabled it to metabolize xylose and to produce 2,3-butanediol (2,3-BDO) as a dominant fermentation product. When co-fermenting with xylose, glucose is preferentially utilized, even though xylose metabolism generates ATP more efficiently during 2,3-BDO production on a BDO-mol basis. To gain a deeper understanding of Z. mobilis metabolism, we first estimated the kinetic parameters of the glucose facilitator protein of Z. mobilis by fitting a kinetic uptake model, which shows that the maximum transport capacity of glucose is seven times higher than that of xylose, and glucose is six times more affinitive to the transporter than xylose. With these estimated kinetic parameters, we further compared the thermodynamic driving force and enzyme protein cost of glucose and xylose metabolism. It is found that, although 20% more ATP can be yielded stoichiometrically during xylose utilization, glucose metabolism is thermodynamically more favorable with 6% greater cumulative Gibbs free energy change, more economical with 37% less enzyme cost required at the initial stage and sustains the advantage of the thermodynamic driving force and protein cost through the fermentation process until glucose is exhausted. Glucose-6-phosphate dehydrogenase (g6pdh), glyceraldehyde-3-phosphate dehydrogenase (gapdh) and phosphoglycerate mutase (pgm) are identified as thermodynamic bottlenecks in glucose utilization pathway, as well as two more enzymes of xylose isomerase and ribulose-5-phosphate epimerase in xylose metabolism. Acetolactate synthase is found as potential engineering target for optimized protein cost supporting unit metabolic flux. Pathway analysis was then extended to the core stoichiometric matrix of Z. mobilis metabolism. Growth was simulated by dynamic flux balance analysis and the model was validated showing good agreement with experimental data. Dynamic FBA simulations suggest that a high agitation is preferable to increase 2,3-BDO productivity while a moderate agitation will benefit the 2,3-BDO titer. Taken together, this work provides thermodynamic and kinetic insights of Z. mobilis metabolism on dual substrates, and guidance of bioengineering efforts to increase hydrocarbon fuel production.

23-butanediol↗

Comparison of loads from HAWC2 and OpenFAST for the IEA Wind 15 MW Reference Wind Turbine

Reference wind turbines (RWTs) that reflect the state-of-the-art of current wind energy technology are necessary in order to properly evaluate innovative methods in wind turbine design and evaluation. The International Energy Agency (IEA) Wind Technology Collaboration Platform (TCP) Task 37 has recently developed a new RWT geared towards offshore floating-foundation applications: the IEA Wind 15 MW. The model has been implemented in two aeroelastic codes, OpenFAST and HAWC2, based on an underlying common ontology. However, these toolchains result in slightly different structural parameters, and the two codes utilise different structural models. Thus, to increase the utility of the model, it is necessary to compare the aeroelastic responses. This paper compares aeroelastic loads calculated using different fidelities of the blade model in OpenFAST (ElastoDyn and BeamDyn) and HAWC2 (prismatic Timoshenko without torsion and Timoshenko with fully populated stiffness matrix), where both codes use the DTU Basic controller and the same turbulence boxes to reduce discrepancies. The aeroelastic responses to steady wind, step wind and turbulent wind (per IEC 61400-1 wind class IB) are considered. The results indicate a generally good agreement between the loads dominated by aerodynamic thrust and force, especially for the no-torsion blade models. Discrepancies were observed in other load channels, partially due to differences in the asymmetric loading of the rotor and partially due to differing closed-loop dynamics, and they will be the subject of future investigations.

17 WIND ENERGY↗

Many-Body Potential for Simulating Self-Assembly of Polymer-Grafted Nanoparticles in a Polymer Matrix

Many-body interactions between polymer-grafted nanoparticles (NPs) play a key role in promoting their assembly into low-dimensional structures within polymer melts, even when the particles are spherical and isotropically grafted. However, capturing such interactions in simulations of NP assembly is very challenging because explicit modeling of the polymer grafts and melt chains is highly computationally expensive, even using coarse-grained models. Here, we develop a many-body potential for describing the effective interactions between spherical polymer-grafted NPs in a polymer matrix through a machine-learning approach. The approach involves using permutationally invariant polynomials to fit two- and three-body interactions derived from the potential of mean force calculations. The potential developed here reduces the computational cost by several orders of magnitude, thereby, allowing us to explore assembly behavior over large length and time scales. We show that the potential not only reproduces previously known assembled phases such as 1D strings and 2D hexagonal sheets, which generally cannot be achieved using isotropic two-body potentials, but can also help discover interesting phases such as networks, clusters, and gels. We demonstrate how each of these assembly morphologies intrinsically arises from a competition between two- and three-body interactions. Our approach for deriving many-body effective potentials can be readily extended to other colloidal systems, enabling researchers to make accurate predictions of their behavior and dissect the role of individual interaction energy terms of the overall potential in the observed behavior.

36 MATERIALS SCIENCE↗

Structure of two-loop SMEFT anomalous dimensions via on-shell methods

We describe on-shell methods for computing one- and two-loop anomalous dimensions in the context of effective field theories containing higher-dimension operators. We also summarize methods for computing one-loop amplitudes, which are used as inputs to the computation of two-loop anomalous dimensions, and we explain how the structure of rational terms and judicious renormalization scheme choices can lead to additional vanishing terms in the anomalous dimension matrix at two loops. We describe the two-loop implications for the Standard Model Effective Field Theory (SMEFT). As a by-product of this analysis we verify a variety of one-loop SMEFT anomalous dimensions computed by Alonso, Jenkins, Manohar and Trott.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Distinct Mechanism of Anti-Corrosion and Swelling-Adhesion Modeling of Low-Dimensional Nylon-Fluoropolymer Composite Coatings

Compared to other polymers, composite coatings with fluoropolymers as the matrix have attracted considerable interest due to their mechanical performance, chemical stability, and low surface energy. Herein, we present an anticorrosion mechanism of fluoropolymer composite coatings achieved by dispersing varying amounts of polyamide 12 (PA-12) particles over a poly(vinylidene fluoride-co-hexafluoropropylene) (PVDF-HFP) matrix. The choice of PVDF-HFP as a coating matrix was driven by its superior processability and mechanical properties over Teflon, while PA-12 afforded swelling capacity, preventing electrolyte permeation, and inhibiting matrix stress crack propagation. The corrosion resistance of the coatings was assessed mainly by potentiodynamic polarization and impedance measurements while being immersed in a NaCl solution. Our electrochemical measurement findings showed that 0.75% w/w PA-12 in the PVDF-HFP matrix significantly enhanced corrosion protection even after 21 days of immersion. Microscopy, dielectric spectroscopy, surface analysis, and mechanical testing (American Society for Testing and Materials standards) corroborated the results. PVDF-HFP coatings containing a minimum percolation threshold of 0.75 wt % PA-12 with 1 mil thickness on mild steel exhibited a remarkable increase in resistance and film adhesion, and a notable corrosion rate decrease. Finite element analysis simulation with a bilinear traction-separation law confirmed the swelling mechanism and adhesion behavior. Finally, this study highlights the potential of nylon particle dispersion in PVDF-HFP matrices as a practical approach to developing advanced anticorrosion coatings with promising practical applications in various industries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using Methyl Bromide for Interspecies Cell–Cell Signaling and As a Reporter in a Model Soil Consortium

Synthetic microbial communities are emerging as model systems for studying consortia-scale phenotypes. To establish synthetic biology tools for studying these communities in hard-to-image environmental materials, we evaluated whether a single member of a model soil consortium (MSC) can be programmed to report on gene expression without requiring matrix disruption. For these studies, we targeted a five member MSC that includes Dyadobacter fermentans, Ensifer adhaerens, Rhodococcus sp003130705, Streptomyces sp001905665, and Variovorax beijingensis. By coupling the expression of a methyl halide transferase to a constitutive promoter, we show that Variovorax beijingensis can be programmed to synthesize methyl halides that accumulate in the soil headspace at levels that are ≥24-fold higher than all other MSC members across a range of environmentally-relevant hydration conditions. We find that methyl halide production can report on a MSC promoter that is activated by changes in water potential, and we demonstrate that a synthetic gas signal can be read out directly using gas chromatography and indirectly using a soil-derived Methylorubrum that is programmed to produce a visual output in response to methyl halides. Finally, these tools will be useful for future studies that investigate how MSC respond to dynamic hydration conditions, such as drought and flood events induced by climate change, which can alter soil water potential and induce the release of stored carbon.

59 BASIC BIOLOGICAL SCIENCES↗

Multi-energy reconstructions, central electron temperature measurements, and early detection of the birth and growth of runaway electrons using a versatile soft x-ray pinhole camera at MST

A multi-energy soft x-ray pinhole camera has been designed, built, and deployed at the Madison Symmetric Torus to aid the study of particle and thermal transport, as well as MHD stability physics. This novel imaging diagnostic technique employs a pixelated x-ray detector in which the lower energy threshold for photon detection can be adjusted independently on each pixel. The detector of choice is a PILATUS3 100 K with a 450 μm thick silicon sensor and nearly 100 000 pixels sensitive to photon energies between 1.6 and 30 keV. An ensemble of cubic spline smoothing functions has been applied to the line-integrated data for each time-frame and energy-range, obtaining a reduced standard-deviation when compared to that dominated by photon-noise. The multi-energy local emissivity profiles are obtained from a 1D matrix-based Abel-inversion procedure. Central values of Te can be obtained by modeling the slope of the continuum radiation from ratios of the inverted radial emissivity profiles over multiple energy ranges with no a priori assumptions of plasma profiles, magnetic field reconstruction constraints, high-density limitations, or need of shot-to-shot reproducibility. In tokamak plasmas, a novel application has recently been tested for early detection, 1D imaging, and study of the birth, exponential growth, and saturation of runaway electrons at energies comparable to 100 × Te,0; thus, early results are also presented.

Delgado-Aparicio, L. F. (ORCID:0000000167394380)↗

Deviations from the Porter-Thomas Distribution due to Nonstatistical 𝛾 Decay below the 150 Nd Neutron Separation Threshold

We introduce a new method for the study of fluctuations of partial transition widths based on nuclear resonance fluorescence experiments with quasimonochromatic linearly polarized photon beams below particle separation thresholds. It is based on the average branching of decays of 𝐽=1 states of an even-even nucleus to the 2$^{+}_{1}$ state in comparison to the ground state. Between 5 and 7 MeV, a constant average branching ratio for 𝛾 decays from 1 − states of 0.490(16) is observed for the nuclide 150 Nd. Assuming 𝜒 2 -distributed partial transition widths, this average branching ratio is related to a degree of freedom of 𝜈 = 1.93⁢(12), rejecting the validity of the Porter-Thomas distribution, requiring 𝜈 = 1. The observed deviation can be explained by nonstatistical effects in the 𝛾-decay behavior with contributions in the range of 9.4(10)% up to 94(10)%.

150 ≤ A ≤ 189↗

Monte Carlo control loops for cosmic shear cosmology with DES Year 1 data

Weak lensing by large-scale structure is a powerful probe of cosmology and of the dark universe. This cosmic shear technique relies on the accurate measurement of the shapes and redshifts of background galaxies and requires precise control of systematic errors. Monte Carlo control loops (MCCL) is a forward modeling method designed to tackle this problem. It relies on the ultra fast image generator (UFig) to produce simulated images tuned to match the target data statistically, followed by calibrations and tolerance loops. Here, we present the first end-to-end application of this method, on the Dark Energy Survey (DES) Year 1 wide field imaging data. We simultaneously measure the shear power spectrum $C_ℓ$ and the redshift distribution $n(z)$ of the background galaxy sample. The method includes maps of the systematic sources, point spread function (PSF), an approximate Bayesian computation (ABC) inference of the simulation model parameters, a shear calibration scheme, and a fast method to estimate the covariance matrix. We find a close statistical agreement between the simulations and the DES Y1 data using an array of diagnostics. In a nontomographic setting, we derive a set of $C_ℓ$ and $n(z)$ curves that encode the cosmic shear measurement, as well as the systematic uncertainty. Following a blinding scheme, we measure the combination of $Ω_m$, $σ_8$, and intrinsic alignment amplitude $A_{IA}$, defined as $S_8D_{IA}=σ_8(Ω_m/0.3)^{0.5}D_{IA}$, where $D_{IA}=1-0.11(A_{IA}-1)$. We find $S_8D_{IA}=0.8954_{-0.039}^{+0.054}$, where systematics are at the level of roughly 60% of the statistical errors. We discuss these results in the context of earlier cosmic shear analyses of the DES Y1 data. Our findings indicate that this method and its fast runtime offer good prospects for cosmic shear measurements with future wide-field surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Nonrenormalization and Operator Mixing via On-Shell Methods

Using on-shell methods, we present a new perturbative nonrenormalization theorem for operator mixing in massless four-dimensional quantum field theories. By examining how unitarity cuts of form factors encode anomalous dimensions, we show that longer operators are often restricted from renormalizing shorter operators at the first order where Feynman diagrams exist. The theorem applies quite generally and depends only on the field content of the operators involved. We apply our theorem to operators of dimension five through seven in the standard model effective field theory, including examples of nontrivial zeros in the anomalous-dimension matrix at one through four loops. The zeros at two and higher loops go beyond those previously explained using helicity selection rules. We also include explicit sample calculations at two loops.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A deep learning approach to real-time HIV outbreak detection using genetic data

Pathogen genomic sequence data are increasingly made available for epidemiological monitoring. A main interest is to identify and assess the potential of infectious disease outbreaks. While popular methods to analyze sequence data often involve phylogenetic tree inference, they are vulnerable to errors from recombination and impose a high computational cost, making it difficult to obtain real-time results when the number of sequences is in or above the thousands. Here, we propose an alternative strategy to outbreak detection using genomic data based on deep learning methods developed for image classification. The key idea is to use a pairwise genetic distance matrix calculated from viral sequences as an image, and develop convolutional neutral network (CNN) models to classify areas of the images that show signatures of active outbreak, leading to identification of subsets of sequences taken from an active outbreak. We showed that our method is efficient in finding HIV-1 outbreaks with R0 ≥ 2.5, and overall a specificity exceeding 98% and sensitivity better than 92%. We validated our approach using data from HIV-1 CRF01 in Europe, containing both endemic sequences and a well-known dual outbreak in intravenous drug users. Our model accurately identified known outbreak sequences in the background of slower spreading HIV. Importantly, we detected both outbreaks early on, before they were over, implying that had this method been applied in real-time as data became available, one would have been able to intervene and possibly prevent the extent of these outbreaks. This approach is scalable to processing hundreds of thousands of sequences, making it useful for current and future real-time epidemiological investigations, including public health monitoring using large databases and especially for rapid outbreak identification.

59 BASIC BIOLOGICAL SCIENCES↗

A Provably Accurate Randomized Sampling Algorithm for Logistic Regression

In statistics and machine learning, logistic regression is a widely-used supervised learning technique primarily employed for binary classification tasks. When the number of observations greatly exceeds the number of predictor variables, we present a simple, randomized sampling-based algorithm for logistic regression problem that guarantees high-quality approximations to both the estimated probabilities and the overall discrepancy of the model. Our analysis builds upon two simple structural conditions that boil down to randomized matrix multiplication, a fundamental and well-understood primitive of randomized numerical linear algebra. We analyze the properties of estimated probabilities of logistic regression when leverage scores are used to sample observations, and prove that accurate approximations can be achieved with a sample whose size is much smaller than the total number of observations. To further validate our theoretical findings, we conduct comprehensive empirical evaluations. Overall, our work sheds light on the potential of using randomized sampling approaches to efficiently approximate the estimated probabilities in logistic regression, offering a practical and computationally efficient solution for large-scale datasets.

Chowdhury, Agniva↗

Equation-Free Coarse Control of Distributed Parameter Systems via Local Neural Operators

The control of high-dimensional distributed parameter systems (DPS) remains a challenge when explicit coarse-grained equations are unavailable. Classical equation-free (EF) approaches rely on fine-scale simulators treated as black-box timesteppers. However, repeated simulations for steady-state computation, linearization, and control design are often computationally prohibitive, or the microscopic timestepper may not even be available, leaving us with data as the only resource. We propose a data-driven alternative that uses local neural operators, trained on spatiotemporal microscopic/mesoscopic data, to obtain efficient short-time solution operators. These surrogates are employed within Krylov subspace methods to compute coarse steady and unsteady-states, while also providing Jacobian information in a matrix-free manner. Krylov-Arnoldi iterations then approximate the dominant eigenspectrum, yielding reduced models that capture the open-loop slow dynamics without explicit Jacobian assembly. Both discrete-time Linear Quadratic Regulator (dLQR) and pole-placement (PP) controllers are based on this reduced system and lifted back to the full nonlinear dynamics, thereby closing the feedback loop.

93B52, 93C20, 47N70, 65J15, 65M32, 68T07, 68T20, 6↗

Climate response to off-equatorial stratospheric sulfur injections in three Earth system models – Part 1: Experimental protocols and surface changes

There is now substantial literature on climate model studies of equatorial or tropical stratospheric SO 2 injections that aim to counteract the surface warming produced by rising concentrations of greenhouse gases. Here we present the results from the first systematic intercomparison of climate responses in three Earth system models wherein the injection of SO 2 occurs at different latitudes in the lower stratosphere: CESM2-WACCM6, UKESM1.0 and GISS-E2.1-G. The first two use a modal aerosol microphysics scheme, while two versions of GISS-E2.1-G use a bulk aerosol (One-Moment Aerosol, OMA) and a two-moment (Multiconfiguration Aerosol TRacker of mIXing state, MATRIX) microphysics approach, respectively. Our aim in this work is to determine commonalities and differences between the climate model responses in terms of the distribution of the optically reflective sulfate aerosols produced from the oxidation of SO 2 and in terms of the surface response to the resulting reduction in solar radiation. A focus on understanding the contribution of characteristics of models transport alongside their microphysical and chemical schemes, and on evaluating the resulting stratospheric responses in different models, is given in the companion paper (Bednarz et al., 2023). The goal of this exercise is not to evaluate these single-point injection simulations as stand-alone proposed strategies to counteract global warming; instead we determine sources and areas of agreement and uncertainty in the simulated responses and, ultimately, the possibility of designing a comprehensive intervention strategy capable of managing multiple simultaneous climate goals through the combination of different injection locations. We find large disagreements between GISS-E2.1-G and the CESM2-WACCM6 and UKESM1.0 models regarding the magnitude of cooling per unit of aerosol optical depth (AOD) produced, which varies from 4.7 K per unit of AOD in CESM2-WACCM6 to 16.7 K in the GISS-E2.1-G version with two-moment aerosol microphysics. By normalizing the results with the global mean response in each of the models and thus assuming that the amount of SO 2 injected is a free parameter that can be managed independently, we highlight some commonalities in the overall distributions of the aerosols, in the inter-hemispheric surface temperature response and in shifts to the Intertropical Convergence Zone, as well as some areas of disagreement, such as the extent of the aerosol confinement in the equatorial region and the efficiency of the transport to polar latitudes. In conclusion, we demonstrate that it is possible to use these simulations to produce more comprehensive injection strategies in multiple climate models. However, large differences in the injection magnitudes can be expected, potentially increasing inter-model spreads in some stratospheric quantities (such as aerosol distribution) while reducing the spread in the surface response in terms of temperature and precipitation; furthermore, the selection of the injection locations may be dependent on the models' specific stratospheric transport.

54 ENVIRONMENTAL SCIENCES↗

Electrolytic gold plating, stripping, and ion transport dynamics through a solid-state iodide perovskite

The pronounced electrochemical reactivity between halide perovskites and metal electrodes can introduce mobile extrinsic metal ions which can cause device instability or enable novel functionalities. Here we systematically investigate the kinetics of gold cation (Au + ) migration in indium tin oxide (ITO)/methylammonium lead triiodide (MAPbI 3 )/Au model devices under long-term potentiostatic biasing. Scanning electron microscopy (SEM), X-ray photoelectron spectroscopy (XPS), and density functional theory (DFT) analyses reveal that Au + ions, electrochemically generated at the Au anode, traverse the perovskite layer with diffusion coefficients on the order of 10 −11 to 10 −10 cm 2 s −1 and are subsequently reduced at the cathode as Au 0 clusters, resembling metal plating behavior in electrolytic cells and solid-state batteries during charging. Furthermore, reversing the applied bias strips the plated Au 0 , revealing reversibility suitable for bipolar resistive switching devices and providing direct evidence of the electrochemical and ionic nature of Au transport within the perovskite matrix. Quantitatively determining diffusion coefficients and ion concentrations provides foundational inputs for future drift-diffusion modelling opportunities and allows us to relate our findings to implications on long term operation of devices like photovoltaic modules. These results clearly demonstrate the solid-state electrochemical nature of perovskite devices, highlight methods to be more quantitative about ion transport properties, provide and emphasize the importance of disentangling electro-, photo-, photoelectrochemical processes for understanding device performance and unlocking new functionalities.

14 SOLAR ENERGY↗

Coherence from multiorbital tunneling ionization of molecules

We present a simple and general coherence model for multiorbital tunnel ionization of molecules, which we incorporate into our previously developed density-matrix approach for sequential double ionization. Here, the influence of this coherence is investigated through simulations of single ionization and sequential double ionization of N 2 and O 2 using few-cycle near-infrared laser pulses. In the case of single ionization, our results reveal the crucial role played by this coherence in generating population inversion in N 2 + , suggesting a potential mechanism for air lasing. Regarding sequential double ionization, we observe only minor changes in the kinetic energy release spectra when the coherence is included, while noticeable differences in the angle-dependent dication yield for both N 2 and O 2 are found. Based on these findings, we recommend the inclusion of multiorbital tunnel ionization coherence in models for single ionization of general molecules, while suggesting that it can be safely neglected in the case of sequential double ionization.

74 ATOMIC AND MOLECULAR PHYSICS↗

Quantifying the impact of temporal analysis of products reactor initial state uncertainties on kinetic parameters

The temporal analysis of products (TAP) reactor, a transient kinetic tool, provides users with information as the catalyst state evolves. However, the state of the art for TAP uncertainty quantification only considers the experimental noise present in the outlet flow signal. Additional sources of uncertainty, including, but not limited to, initial surface coverages, the catalyst zone location, the inert void fraction, and the gas pulse intensity and pulse delay, are not included. For this reason, a framework for quantifying all uncertainty sources present in TAP experiments is presented and applied to a carbon monoxide oxidation case study. Herein, two methods for quantifying these sources of uncertainty are introduced. The first utilizes initial state sensitivities to approximate the parameter variances, as well as to provide insights into the structural certainty of the model. The second generates parameter confidence distributions through an ensemble-based sampling algorithm. This initial state covariance matrix can ultimately be merged with the experimental noise covariance matrix, providing a unified description of the parameter uncertainties for a single TAP experiment.

36 MATERIALS SCIENCE↗

Phase-Field Modeling of Thermally-Grown Oxide and the Induced Damage Evolution in Environmental Barrier Coatings

The advent of next-generation hydrogen-based engines necessitates materials capable of withstanding temperatures beyond the reach of current superalloys. SiC-based ceramic matrix composites, augmented with environmental barrier coatings (EBCs), present a promising materials solution. Given the active search for effective and durable EBCs, there is a pressing need for modeling tools to understand and predict damage evolution in these materials to help accelerate their development. This study introduces a phase-field model (PFM) designed to simulate the critical role of thermally grown oxides (TGO) in the degradation and failure of EBCs. The model accounts for the severe volume expansion due to oxidation, alongside phase transformations and microstructural evolution during thermal cycling, offering a comprehensive view of the damage processes. Simulation results are validated against experimental findings reported in the literature, establishing the model's potential as a significant tool for understanding and improving the resilience of EBCs in cyclic oxidative environments.

Cheng, Tianle↗