Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear systems”

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 595 records · Page 33

Modal dynamics and flutter analysis of floating offshore vertical axis wind turbines

Floating platforms provide vertical axis wind turbines (VAWTs) significant advantages over offshore horizontal axis wind turbines, which has led to increased research interest in offshore VAWT technology. Here, this paper examines dynamic stability of floating VAWTs and presents structural dynamic and aero-elastic flutter analysis of utility-scale 5 MW floating VAWTs by coupling a rotor finite element model with a linearized model for the floating system. The coupled floating VAWT model is utilized to study the effects of tower height, number of blades, blade tapering scheme, and impacts of the floating platform on modal properties and flutter. The study provides important new results and guidance to floating VAWT designers: (1) tower height increase leads to lower tower frequencies more prone to resonance, (2) by moving from land-based foundations to a floating TLP foundation, tower-mode frequency is increased significantly (by 49%–68%) which mitigates resonance and increases the flutter RPM (154% increase), (3) increasing the number of blades reduces the tower mode frequencies for both land-based and floating configuration, and (4) blade chord tapering can have a significant impact on increasing the tower mode flutter RPM. Further, analysis of low-frequency rigid body modes shows a safe, higher flutter RPM relative to the operating RPM.

17 WIND ENERGY↗

Model-based and Model-free Designs for an Extended Continuous-time LQR with Exogenous Inputs

We present an extended linear quadratic regulator (LQR) design for continuous-time linear time-invariant (LTI) systems in the presence of exogenous inputs. We first propose a model-based solution with cost minimization guarantees for states and inputs using dynamic programming (DP). The control law consists of a combination of the optimal state feedback and an additional optimal term dependent on the exogenous inputs. The control gains for the two components are obtained by solving a set of matrix differential equations. We provide these solutions for both finite horizons and steady-state cases. In the second part of the paper, we formulate a reinforcement learning (RL) based algorithm which does not need any model information except the input matrix, and can compute an approximate steady-state LQR gain using measurements of the states, the control inputs, and the exogenous inputs. Both model-based and data-driven optimal control algorithms are tested with a numerical example under different exogenous inputs showcasing the effectiveness of the designs.

Mukherjee, Sayak↗

Status and future plans for C 3 R&D

C 3 is an opportunity to realize an e + e - collider for the study of the Higgs boson at √s = 250 GeV, with a well defined upgrade path to 550 GeV while staying on the same short facility footprint. C 3 is based on a fundamentally new approach to normal conducting linear accelerators that achieves both high gradient and high efficiency at relatively low cost. Given the advanced state of linear collider designs, the key system that requires technical maturation for C 3 is the main linac. This paper presents the staged approach towards a facility to demonstrate C 3 technology with both Direct (source and main linac) and Parallel (beam delivery, damping ring, ancillary component) R&D. The primary goal of the C 3 Demonstration R&D Plan is to reduce technical and cost risk by building and operating the key components of C 3 at an adequate scale. This R&D plan starts with the engineering design, and demonstration of one cryomodule and will culminate in the construction of a 3 cryomodule linac with pre-production prototypes. This R&D program would also demonstrate the linac rf fundamentals including achievable gradient and gradient stability over a full electron bunch train and breakdown rates. It will also investigate beam dynamics including energy spread, wakefields, and emittance growth. This work will be critical to confirm the suitability of the C 3 beam parameters for the physics reach and detector performance in preparation for a Conceptual Design Report (CDR), as well as for follow-on technology development and industrialization. The C 3 Demonstration R&D Plan will open up significant new scientific and technical opportunities based on development of high-gradient and high-efficiency accelerator technology. It will push this technology to operate both at the GeV scale and mature the technology to be reliable and provide high-brightness electron beams. The timeline for progressing with C 3 technology development will be governed by practical limitations on both the technical progress and resource availability. It consists of four stages: Stage 0) Ongoing fundamental R&D on structure prototypes, damping and vibrations. Stage 1) Advancing the engineering maturity of the design and developing start-to-end simulations including space-charge and wakefield effects. This stage will include testing of strucutres operating at cryogenic temperatures. Beam tests would be performed with high beam current to test full beam loading. Stage 2) Production and testing of the first cryomodule at cryogenic temperatures. This would provide sufficient experimental data to compile a CDR and it is anticipated for Stage 2 to last 3 years and to culminate with the transport of photo-electrons through the first cryomodule. Stage 3) Updates to the engineering design of the cryomodules, production of the second and third cryomodule and their installation. Lower charge and lower emittance beams will be used to investigate emittance growth. The successful full demonstration of the 3 cryomodules to deliver up to a 3 GeV beam and achieve the C 3 five gradient will allow a comprehensive and robust evaluation of the technical design of C 3 as well as mitigate technical, schedule, and cost risks required to proceed with a Technical Design Report (TDR).

radiation hardened magnets↗

Simulating quantum-classical interfaces via the Lindblad master equation

In hybrid quantum systems, the interface between quantum and classical domains is essential for the generation, control, and measurement of quantum states. Quantum-classical interfaces (QCIs) are ubiquitous in devices such as optical modulators, quantum sensors, and signal processors, where classical signals influence quantum dynamics. In this paper, we employ the Lindblad master equation to simulate the evolution of a quantum system interacting with a classical control system. Our model captures both linear and nonlinear interactions by incorporating first- and second-order susceptibilities, and it quantifies the influence of externally applied control parameters on decoherence and state evolution. As an illustrative example, we analyze an optical modulator and demonstrate how variations in material response and drive conditions affect photon statistics, coherence, and phase-space distributions. In conclusion, the findings offer a path to an all-encompassing model for understanding and optimizing QCIs, with wide-ranging implications for the performance, design, and robustness of next-generation quantum devices.

Quantum engineering↗

Super-swelling behavior of stacked lipid bilayer systems

Bilayer systems comprising lipid mixtures are the most well-studied model of biological membranes. While the plasma membrane of the cell is a single bilayer, many intra- and extra-cellular biomembranes comprise stacks of bilayers. Most bilayer stacks in nature are periodic, maintaining a precise water layer separation between bilayers. That equilibrium water separation is governed by multiple inter-bilayer forces and is highly responsive. Biomembranes re-configure inter-bilayer spacing in response to temperature, composition, or mass transport cues. In synthetic bilayer systems for applications in cosmetics or topical treatments, control of the hydration level is a critical design handle. Herein we investigate a binary lipid system that leverages key inter-bilayer forces leading to unprecedented levels of aqueous swelling while maintaining a coherent multilamellar form. We found that combining cationic lipids with bicontinuous cubic phase-forming lipids (lipids with positive Gaussian modulus), results in the stabilization of multilamellar phases against repulsive steric forces that typically lead to bilayer delamination at high degrees of swelling. Using ultra-small-angle X-ray scattering alongside confocal laser scanning microscopy, we characterized various super-swelled states of 1,2-dioleoyl-3-trimethylammonium-propane (DOTAP) and glycerol monooleate (GMO) lipids, as well as other analogous systems, at varied concentration and molar ratios. Through these experiments we established swelling profiles of various binary lipid systems that were near-linear with decreasing lipid volume fraction, showing maximum swelling with periodicity well above 200 nanometers.

42 ENGINEERING↗

Idealized Large-Eddy Simulations of Stratocumulus Advecting over Cold Water. Part I: Boundary Layer Decoupling

Abstract We explore the decoupling physics of a stratocumulus-topped boundary layer (STBL) moving over cooler water, a situation mimicking warm-air advection (WADV). We simulate an initially well-mixed STBL over a doubly periodic domain with the sea surface temperature decreasing linearly over time using the System for Atmospheric Modeling large-eddy model. Due to the surface cooling, the STBL becomes increasingly stably stratified, manifested as a near-surface temperature inversion topped by a well-mixed cloud-containing layer. Unlike the stably stratified STBL in cold-air advection (CADV) that is characterized by cumulus coupling, the stratocumulus deck in the WADV is unambiguously decoupled from the sea surface, manifested as weakly negative buoyancy flux throughout the subcloud layer. Without the influxes of buoyancy from the surface, the convective circulation in the well-mixed cloud-containing layer is driven by cloud-top radiative cooling. In such a regime, the downdrafts propel the circulation, in contrast to that in CADV regime for which the cumulus updrafts play a more determinant role. Such a contrast in convection regime explains the difference in many aspects of the STBLs including the entrainment rate, cloud homogeneity, vertical exchanges of heat and moisture, and lifetime of the stratocumulus deck, with the last being subject to a more thorough investigation in Part II. Finally, we investigate under what conditions a secondary stratus near the surface (or fog) can form in the WADV. We found that weaker subsidence favors the formation of fog whereas a more rapid surface cooling rate does not. Significant Statement The low-lying blanketlike clouds, called stratocumulus (Sc), reflect much incoming sunlight, substantially modulating Earth’s temperature. While much is known about how the Sc evolves when it moves over warmer water, few studies examine the opposite situation of Sc moving over colder water. We used a high-resolution numerical model to simulate such a case. When moving over cold water, the Sc becomes unambiguously decoupled from the water surface, distinctive from its warm counterpart in which the Sc interacts with the water surface via intermittent cauliflower-like clouds called cumulus clouds. Such decoupling influences many aspects of the Sc–sea surface system, which combine to alter the ability of the Sc to reflect sunlight, thereby influencing the climate. This work laid the foundation for future work that quantifies the contribution of such a decoupled Sc regime to Earth’s radiative budget and climate change.

54 ENVIRONMENTAL SCIENCES↗

Forced Component Estimation Statistical Method Intercomparison Project (ForceSMIP)

Anthropogenic climate change is unfolding rapidly, yet its regional manifestation can be obscured by internal variability. A primary goal of climate science is to identify the externally forced climate response from among the noise of internal variability. Separating the forced response from internal variability can be addressed in climate models by using a large ensemble to average over different possible realizations of internal variability. However, with only one realization of the real world, it is a major challenge to isolate the forced response directly in observations. In the Forced Component Estimation Statistical Method Intercomparison Project (ForceSMIP), contributors used existing and newly developed statistical and machine learning methods to estimate the forced response over 1950–2022 within individual realizations of the climate system. Participants used neural networks, linear inverse models, fingerprinting methods, and low-frequency component analysis, among other approaches. These methods were trained using large ensembles from multiple climate models and then applied to observations. Here, we evaluate method performance within large ensembles and investigate the estimates of the forced response in observations. Our results show that many different types of methods are skillful for estimating the forced response in climate models, though the relative skill of individual methods varies depending on the variable and evaluation metric. Methods with comparable skill in models can give a wide range of estimates of the forced response pattern in observations, illustrating the epistemic uncertainty in forced response estimates. ForceSMIP gives new insights into the forced response in observations, its uncertainty, and methods for its estimation.

Climate attribution↗

Signal Preconditioning to Minimize Impulse Response Contribution

A study was performed to identify a method to minimize the effect of a linear time-invariant (LTI) system impulse response on an input. Three methods were studied: Wiener filter, the N4SID algorithm and transfer function estimation, the latter two using functions from MATLAB’s System Identification Toolbox. Although all three methods were able estimate an unknown forward impulse response given an input/output time series pair, only the Wiener filter was able to estimate a system inverse which satisfactorily solved the problem using a cosine similarity measure.

42 ENGINEERING↗

Extending Analysis of Asay Windows to High Areal Mass Range [Slides]

The accretion layer addition to Asay window analysis results in new equations of motion that depend on the location of collisions within the layer. Asay window equations of motion can be represented in terms of experiment observables. Balance of linear momentum for variable mass systems provides a framework for developing equations of motion for Asay window diagnostics under a variety of operating assumptions. The growth in the accretion layer thickness during an experiment substantially influences mass and density measurements when high accumulation of areal mass is present.

97 MATHEMATICS AND COMPUTING↗

Multi-Fidelity Stochastic Economic Dispatch for Operating Low Carbon Power Grids

Grid operators can address the inherently stochastic nature of renewables by solving a two-stage stochastic programming model that minimizes the cost of dispatch decisions while accounting for the complex grid dynamics. It is common to use a sample average approximation to estimate the expectation of the second stage costs in this model. However, the large sample count needed for numerical accuracy makes effective modeling large-scale electric grids computationally intractable. We introduce a control variate multi-fidelity estimator for the second-stage recourse that enables high quality dispatch decisions in real-time with a reduced computational burden. We obtain a hierarchy of model fidelities by linearizing the AC power flow system representation to DC power flow, and by relaxing transmission and voltage network constraints. We evaluate the performance of our proposed method on a synthetic grid with 73 buses against a deterministic baseline with persistence forecast and a high-fidelity reference. Our analysis shows a computational speed-up of 7.62x with a minimal loss in accuracy. The multi-fidelity method is well suited to fidelity combinations that use a simplified network topology in their lower fidelity model and is an attractive option for applications where accurate grid modeling needed on a limited computational budget.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Extending the Asay window measurement model to high areal mass ranges [Slides]

The accretion layer addition to Asay window analysis results in new equations of motion that depend on the location of collisions within the layer. Asay window equations of motion can be represented in terms of experiment observables. Balance of linear momentum for variable mass systems provides a framework for developing equations of motion for Asay window diagnostics under a variety of operating assumptions. The growth in the accretion layer thickness during an experiment substantially influences mass and density measurements when high accumulation of areal mass is present.

97 MATHEMATICS AND COMPUTING↗

C$^3$ Demonstration Research and Development Plan

C$^3$ is an opportunity to realize an e$^+$e$^-$ collider for the study of the Higgs boson at $\sqrt{s} = 250$ GeV, with a well defined upgrade path to 550 GeV while staying on the same short facility footprint. C$^3$ is based on a fundamentally new approach to normal conducting linear accelerators that achieves both high gradient and high efficiency at relatively low cost. Given the advanced state of linear collider designs, the key system that requires technical maturation for C$^3$ is the main linac. This white paper presents the staged approach towards a facility to demonstrate C$^3$ technology with both Direct (source and main linac) and Parallel (beam delivery, damping ring, ancillary component) R&D. The white paper also includes discussion on the approach for technology industrialization, related HEP R&D activities that are enabled by C$^3$ R&D, infrastructure requirements and siting options.

43 PARTICLE ACCELERATORS↗

Marine Alternative Fuel Pricing, Supply, and Demand

The International Maritime Organization low sulfur fuel rules and greenhouse gas-reduction strategy may create a significant perturbation to global refinery operations, impacting volumes and prices of marine fuel globally and within the United States and incentivizing the transition to alternative fuel use. The goals of this project are to (1) enhance understanding of how very low-sulfur fuel-primarily fuel oil and diesel-and low carbon fuel requirements and promising biofuel processes will affect the marine fuel supply chain; and (2) explore how these perturbations interact with and impact indicators such as pricing, number of trips, and demand behavior, along with the potential to meet low sulfur and carbon fuel demand with biofuel supply chains, (3) merge innovative thinking in the area of marine fuels within research centers of DOE (NREL) and the U.S. Department of Transportation (Volpe). These goals will be accomplished through the combination of detailed refinery, marine fuel burn, system dynamics, and geospatially explicit linear programming models.

biofuel↗

Projection-based multifidelity linear regression for data-scarce applications

Surrogate modeling for systems with high-dimensional quantities of interest remains challenging, particularly when training data are costly to acquire. This work develops multifidelity methods for multiple-input multiple-output linear regression targeting data-limited applications with high-dimensional outputs. Multifidelity methods integrate many inexpensive low-fidelity model evaluations with limited, costly high-fidelity evaluations. We introduce two projection-based multifidelity linear regression approaches with linear and nonlinear features that leverage principal component basis vectors for dimensionality reduction and combine multifidelity data through: (i) a direct data augmentation using low-fidelity data, and (ii) a data augmentation incorporating explicit linear corrections between low-fidelity and high-fidelity data. The data augmentation approaches combine high-fidelity and low-fidelity data into a unified training set and train the linear regression model through weighted least squares with fidelity-specific weights. We introduce a proximity-based weighting scheme with automatic weight selection strategy through cross-validation. Here, the proposed multifidelity linear regression methods are demonstrated on approximating the surface pressure field of a hypersonic vehicle in flight and the temperature field on an aircraft disc braking system. In an ultra low-data regime of no more than twelve high-fidelity samples, multifidelity linear regression achieves approximately 2% – 12% improvement in median accuracy and a higher R 2 score relative to single-fidelity methods at comparable computational cost.

data augmentation↗

Gate-tunable heavy fermion quantum criticality in a moiré Kondo lattice

We propose a realization of Kondo lattice physics in moiré superlattices at the interface between a WX 2 homobilayer and MoX 2 monolayer (X = S,Se). Under appropriate gating conditions, the interface WX 2 layer forms a triangular lattice of local moments that couple to itinerant electrons in the other WX 2 layer via a gate-tunable Kondo exchange interaction. Using a parton mean-field approach, we identify a range of twist angles which support a gate-tuned quantum phase transition between a heavy Fermi liquid with large anomalous Hall conductance and a chiral spin liquid coexisting with a light Fermi liquid and describe experimental signatures to distinguish among competing theoretical scenarios.

36 MATERIALS SCIENCE↗

Missing data in multi-omics integration: Recent advances through artificial intelligence

Biological systems function through complex interactions between various ‘omics (biomolecules), and a more complete understanding of these systems is only possible through an integrated, multi-omic perspective. This has presented the need for the development of integration approaches that are able to capture the complex, often non-linear, interactions that define these biological systems and are adapted to the challenges of combining the heterogenous data across ‘omic views. A principal challenge to multi-omic integration is missing data because all biomolecules are not measured in all samples. Due to either cost, instrument sensitivity, or other experimental factors, data for a biological sample may be missing for one or more ‘omic techologies. Recent methodological developments in artificial intelligence and statistical learning have greatly facilitated the analyses of multi-omics data, however many of these techniques assume access to completely observed data. A subset of these methods incorporate mechanisms for handling partially observed samples, and these methods are the focus of this review. We describe recently developed approaches, noting their primary use cases and highlighting each method's approach to handling missing data. We additionally provide an overview of the more traditional missing data workflows and their limitations; and we discuss potential avenues for further developments as well as how the missing data issue and its current solutions may generalize beyond the multi-omics context.

97 MATHEMATICS AND COMPUTING↗

Quasi-linear theory of forced magnetic reconnection for the transition from the linear to the Rutherford regime

Abstract Using the in-viscid two-field reduced magneto-hydrodynamic (MHD) model, a new analytical theory is developed to unify the Hahm–Kulsrud–Taylor (HKT) linear solution and the Rutherford quasi-linear regime. Adopting a quasi-linear approach, we obtain a closed system of equations for plasma response in the Taylor problem. An integral form of an analytical solution is obtained for the forced magnetic reconnection, uniformly valid throughout the entire regimes from the HKT linear solution to the Rutherford quasi-linear solution. In particular, the quasi-linear effect can be described by a single coefficient K s ∝ S 8 / 5 ψ c 2 , where S = τ R τ A and ψ c are the Lundquist number and amplitude of external magnetic perturbation, respectively. The HKT linear solution for the response can be recovered when the index K s → 0. On the other hand, the quasi-linear effect plays a key role in the island growth when K s ~ 1. Our new analytical solution has also been compared with reduced MHD simulations with agreement.

Physics↗

Start-up Optimization Considering Integrated Power Electronic Systems

The startup of systems composed of power electronic systems is becoming more complex. As the systems grow, standard pre-configuration sequencing methods will struggle. This work proposes a linear programming optimization method to couple power electronic systems in a series of startup commands. This formulation is demonstrated in a real-time implementation on hardware.

Starke, Michael↗