Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic parameter estimation”

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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

TIC 168789840: A Sextuply Eclipsing Sextuple Star System

We report the discovery of a sextuply eclipsing sextuple star system from TESS data, TIC 168789840, also known as TYC 7037-89-1, the first known sextuple system consisting of three eclipsing binaries. The target was observed in Sectors 4 and 5 during Cycle 1, with lightcurves extracted from TESS Full Frame Image data. It was also previously observed by the WASP survey and ASAS-SN. The system consists of three gravitationally bound eclipsing binaries in a hierarchical structure of an inner quadruple system with an outer binary subsystem. Follow-up observations from several different observatories were conducted as a means of determining additional parameters. The system was resolved by speckle interferometry with a 0.″42 separation between the inner quadruple and outer binary, inferring an estimated outer period of ∼2 kyr. It was determined that the fainter of the two resolved components is an 8.217 day eclipsing binary, which orbits the inner quadruple that contains two eclipsing binaries with periods of 1.570 days and 1.306 days. Markov Chain Monte Carlo (MCMC) analysis of the stellar parameters has shown that the three binaries of TIC 168789840 are “triplets,” as each binary is quite similar to the others in terms of mass, radius, and T {sub eff}. As a consequence of its rare composition, structure, and orientation, this object can provide important new insight into the formation, dynamics, and evolution of multiple star systems. Future observations could reveal if the intermediate and outer orbital planes are all aligned with the planes of the three inner eclipsing binaries.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Eco-PiNN: A Physics-informed Neural Network for Eco-toll Estimation

The eco-toll estimation problem quantifies the expected environmental cost (e.g., energy consumption, exhaust emissions) for a vehicle to travel along a path. This problem is important for societal applications such as eco-routing, which aims to find paths with the lowest exhaust emissions or energy need. The challenges of this problem are threefold: (1) the dependence of a vehicle's eco-toll on its physical parameters; (2) the lack of access to data with eco-toll information; and (3) the influence of contextual information (i.e. the connections of adjacent segments in the path) on the eco-toll of road segments. Prior work on eco-toll estimation has mostly relied on pure data-driven approaches and has high estimation errors given the limited training data. To address these limitations, we propose a novel Eco-toll estimation Physics-informed Neural Network framework (Eco-PiNN) using three novel ideas, namely, (1) a physics-informed decoder that integrates the physical laws governing vehicle dynamics into the network, (2) an attention-based contextual information encoder, and (3) a physics-informed regularization to reduce overfitting. Experiments on real-world heavy-duty truck data show that the proposed method can greatly improve the accuracy of eco-toll estimation compared with state-of-the-art methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2005

Design Optimization Tool for Synthetic Jet Actuators Using Lumped Element Modeling

The performance specifications of any actuator are quantified in terms of an exhaustive list of parameters such as bandwidth, output control authority, etc. Flow-control applications benefit from a known actuator frequency response function that relates the input voltage to the output property of interest (e.g., maximum velocity, volumetric flow rate, momentum flux, etc.). Clearly, the required performance metrics are application specific, and methods are needed to achieve the optimal design of these devices. Design and optimization studies have been conducted for piezoelectric cantilever-type flow control actuators, but the modeling issues are simpler compared to synthetic jets. Here, lumped element modeling (LEM) is combined with equivalent circuit representations to estimate the nonlinear dynamic response of a synthetic jet as a function of device dimensions, material properties, and external flow conditions. These models provide reasonable agreement between predicted and measured frequency response functions and thus are suitable for use as design tools. In this work, we have developed a Matlab-based design optimization tool for piezoelectric synthetic jet actuators based on the lumped element models mentioned above. Significant improvements were achieved by optimizing the piezoceramic diaphragm dimensions. Synthetic-jet actuators were fabricated and benchtop tested to fully document their behavior and validate a companion optimization effort. It is hoped that the tool developed from this investigation will assist in the design and deployment of these actuators.

Gallas, Quentin↗
NASA NTRS2014

Initial Evaluations of LoC Prediction Algorithms Using the NASA Vertical Motion Simulator

Flying near the edge of the safe operating envelope is an inherently unsafe proposition. Edge of the envelope here implies that small changes or disturbances in system state or system dynamics can take the system out of the safe envelope in a short time and could result in loss-of-control events. This study evaluated approaches to predicting loss-of-control safety margins as the aircraft gets closer to the edge of the safe operating envelope. The goal of the approach is to provide the pilot aural, visual, and tactile cues focused on maintaining the pilot's control action within predicted loss-of-control boundaries. Our predictive architecture combines quantitative loss-of-control boundaries, an adaptive prediction method to estimate in real-time Markov model parameters and associated stability margins, and a real-time data-based predictive control margins estimation algorithm. The combined architecture is applied to a nonlinear transport class aircraft. Evaluations of various feedback cues using both test and commercial pilots in the NASA Ames Vertical Motion-base Simulator (VMS) were conducted in the summer of 2013. The paper presents results of this evaluation focused on effectiveness of these approaches and the cues in preventing the pilots from entering a loss-of-control event.

Pilot Cue↗
NASA NTRS2013

Multiwavelength Diagnostics of the Precursor and Main Phases of an M1.8 Flare on 2011 April 22

We study the temporal, spatial and spectral evolution of the M1.8 flare, which occurred in the active region 11195 (S17E31) on 2011 April 22, and explore the underlying physical processes during the precursor phase and their relation to the main phase. The study of the source morphology using the composite images in 131Å wavelength observed by the Solar Dynamics Observatory/Atmospheric Imaging Assembly and 6-14 kiloelectronvolts [from the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI)] revealed a multi-loop system that destabilized systematically during the precursor and main phases. In contrast, hard X-ray emission (20-50 kiloelectronvolts) was absent during the precursor phase, appearing only from the onset of the impulsive phase in the form of foot-points of emitting loops. This study also revealed the heated loop-top prior to the loop emission, although no accompanying foot-point sources were observed during the precursor phase. We estimate the flare plasma parameters, namely temperature (T), emission measure (EM), power-law index (gamma) and photon turn-over energy (to), and found them to be varying in the ranges 12.4-23.4 megakelvins, 0.0003-0.6 x 10 (sup 49) per cubic centimeter, 5-9 and 14-18 kiloelectronvolts, respectively, by forward fitting RHESSI spectral observations. The energy released in the precursor phase was thermal and constituted approximately 1 percent of the total energy released during the flare. The study of morphological evolution of the filament in conjunction with synthesized T and EM maps was carried out, which reveals (a) partial filament eruption prior to the onset of the precursor emission and (b) heated dense plasma over the polarity inversion line and in the vicinity of the slowly rising filament during the precursor phase. Based on the implications from multiwavelength observations, we propose a scheme to unify the energy release during the precursor and main phase emissions in which the precursor phase emission was originated via conduction front that resulted due to the partial filament eruption. Next, the heated leftover S-shaped filament underwent slow-rise and heating due to magnetic reconnection and finally erupted to produce emission during the impulsive and gradual phases.

gamma-rays↗
NASA NTRS1981

Observed dynamics of superclusters

Radial velocity measurements of two large, rich superclusters show evidence of a gravitational slowing of the Hubble flow within the superclusters. Models assuming spherically symmetric, noncrossing mass shells yield mass density estimates for the clusters ranging from less than one to several times the critical density, and preliminary counts of galaxies and Abell clusters indicate a density enhancement within the superclusters of 10-71 times the field mass density. The combination of these two results is found to favor an open universe, assuming a standard Friedmann cosmology. It is noted that all forms of mass that are clumped on the scale of superclusters, or smaller, will be detected by these measurements, including neutrinos sufficiently massive to significantly contribute to the cosmological density parameter.

Harms, R. J.↗
NASA NTRS2005

Optimal Sampling of a Reaction Coordinate in Molecular Dynamics

Estimating how free energy changes with the state of a system is a central goal in applications of statistical mechanics to problems of chemical or biological interest. From these free energy changes it is possible, for example, to establish which states of the system are stable, what are their probabilities and how the equilibria between these states are influenced by external conditions. Free energies are also of great utility in determining kinetics of transitions between different states. A variety of methods have been developed to compute free energies of condensed phase systems. Here, I will focus on one class of methods - those that allow for calculating free energy changes along one or several generalized coordinates in the system, often called reaction coordinates or order parameters . Considering that in almost all cases of practical interest a significant computational effort is required to determine free energy changes along such coordinates it is hardly surprising that efficiencies of different methods are of great concern. In most cases, the main difficulty is associated with its shape along the reaction coordinate. If the free energy changes markedly along this coordinate Boltzmann sampling of its different values becomes highly non-uniform. This, in turn, may have considerable, detrimental effect on the performance of many methods for calculating free energies.

Pohorille, Andrew↗
DOE OSTI2019

Model-Based Localizatiion in a Shallow Ocean Environment: A Sequential Bayesian/Optimization Approach

The shallow ocean is a dynamically changing medium leading to nonstationary statistical behavior when subjected to temperature, wind, surface variations, bottom interactions, noise and extraneous disturbances as well as other conditions that render it a uniquely challenging environment — especially from a signal processing perspective. Processors must account for and adapt to such instantaneous changes in order to be effective. Thus, a processor is required to “adapt” to these environmental variations while simultaneously providing meaningful estimates that are necessary for such applications as detection, localization, inversion and enhancement. In this paper, we develop a parametrically adaptive, sequential Bayesian processor capable of jointly estimating both modal functions and environmental parameters to provide enhanced estimates for a focused optimizer capable locating a target in the noisy shallow ocean environment.

42 ENGINEERING↗
DOE OSTI2022

MEGASIM: Lifetimes and Resonances of Earth Trojan Asteroids—The Death of Primordial ETAs?

Abstract We present an analysis of lifetimes and resonances of Earth Trojan Asteroids (ETAs) in the MEGASIM data set. Trojan asteroids co-orbit the Sun with a planet, but remain bound to the Lagrange points, L4 (60° leading the planet) or L5 (60° trailing). In the circular three-body approximation, the stability of a Trojan asteroid depends on the ratio of the host planet mass and the central mass. For the inner planets, the range of stability becomes increasingly small, so perturbations from the planets have made primordial Trojans rare. To date, there have been just two ETAs (2010 TK 7 and 2020 XL 5 ), several Mars Trojans, and a Venus Trojan discovered. The estimated lifetimes of the known inner system Trojans are shorter than a million years, suggesting they are interlopers rather than members of a stable and long-lasting population. With the largest ETA n -body simulation to date, we are able to track their survival across a wide initialized parameter space. We find that the remaining fraction of ETAs over time is well fit with a stretched exponential function that, when extrapolated beyond our simulation run time, predicts zero ETAs by 2.33 Gyr. We also show correlations between ETA ejections and the periods of the Milankovitch cycles. Though Earth’s orbital dynamics dominate the instabilities of ETAs, we provide evidence that ETA ejections are linked to resonances found in the variation of the orbital elements of many if not all of the planets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Field Sensitivity Analysis of Turbulence Model Parameters for Flow Over a Wing

Reynolds-averaged-Navier-Stokes (RANS) turbulence models are a critical tool in computational-fluid-dynamics simulations of aerodynamic systems, but simulation results can be highly sensitive to RANS-model parameter choices. Sensitivity analysis can be used to quantify these impacts, and the objective of this study is to demonstrate field sensitivity analysis with respect to ten parameters in the 2003 Menter shear-stress-transport (SST) turbulence model. Here, the analysis is demonstrated for an application relevant to wind energy, namely, flow over a NACA 0015 wing at 12 degree angle of attack and a Reynolds number of 1.5 x 10 6 . We quantify sensitivity using Sobol indices and the mean-squared gradient, which are estimated using polynomial chaos and active subspace models, respectively. Our results indicate that there are substantial spatial variations in parameter sensitivities, with different sets of most-sensitive parameters near the wing, as well as in the downstream wake, consistent with the physical interpretations of the turbulence model inputs. We show that, for this particular turbulence model and flow, simultaneous dimension reduction is possible across all quantities of interest, enabling efficient exploration of model outcomes. Ultimately, this analysis provides new insights into turbulence model parameter sensitivities in incompressible flows, and also demonstrates the implementation of field sensitivity analysis for applications relevant to aerodynamics simulations.

17 WIND ENERGY↗
DOE OSTI2020

Sound Speed as a Candidate for Internal Temperature Monitoring During Solid Phase Processing of Materials

Internal temperature is an important parameter for understanding the solid phase processing (SPP) of materials and improving feedback control of these processes. Friction stir processes, shear assisted processing and extrusion (ShAPE ™ ), and ultrahigh velocity (UHV) cold spray are examples of such processes. Simple thermal models based on temperatures measured at the boundary of the deforming body cannot capture the dynamic nature of temperature generation within the deforming solid. One possibly way of understanding internal temperature is by measuring sound speed, which varies as a function of temperature for materials and directing ultrasonic wave propagation through selected material volumes of interest. In this project, longitudinal wave speed is estimated for aluminum 6061. Piezoelectric materials bonded to a specimen provided data from which longitudinal wave speed was estimated. A discussion evaluated the linear relationship of wave speed versus temperature and postulated how wave speed measurements could be used to estimate the internal temperature of selected material volumes for feedback during SPP of materials. Instrumentation, data and analysis, and a discussion are reviewed.

36 MATERIALS SCIENCE↗
NASA NTRS2001

Molecular Dynamic Simulations of Interaction of an AFM Probe with the Surface of an SCN Sample

Molecular dynamic (MD) simulations is conducted in order to estimate forces of probe-substrate interaction in the Atomic Force Microscope (AFM). First a review of available molecular dynamic techniques is given. Implementation of MD simulation is based on an object-oriented code developed at the University of Delft. Modeling of the sample material - succinonitrile (SCN) - is based on the Lennard-Jones potentials. For the polystyrene probe an atomic interaction potential is used. Due to object-oriented structure of the code modification of an atomic interaction potential is straight forward. Calculation of melting temperature is used for validation of the code and of the interaction potentials. Various fitting parameters of the probe-substrate interaction potentials are considered, as potentials fitted to certain properties and temperature ranges may not be reliable for the others. This research provides theoretical foundation for an interpretation of actual measurements of an interaction forces using AFM.

Bune, Adris↗
NASA NTRS2002

Photoevaporation of Clumps in Photodissociation Regions

We present the results of an investigation of the effects of Far Ultraviolet (FUV) radiation (6.0eV < hv < 13.6eV) from hot early type OB stars on clumps in star-forming molecular clouds. Clumps in FUV-illuminated regions (or photodissociation regions or PDRs) undergo external heating and photodissociation as they are exposed to the FUV field, resulting in a loss of cold, molecular lump mass as it is converted to warm atomic gas. The heating, if rapid, creates strong photoevaporative mass flows off the clump surfaces, and drives shocks into the clumps, compressing them to high densities. The clumps lose mass on relatively short timescales. The evolution of an individual clump is found to be sensitive to three dimensionless parameters: Nc0, the ratio of the initial column density of the clump to the column N(0) approx. 10(exp 21) cm(exp -2) of a warm FUV-heated surface region; upsilon, the ratio of the sound speed in the heated surface to that in the cold clump material: and t(FUV)t(c), the ratio of the "turn-on time" t(FUV) of the heating flux on a clump to its initial sound crossing-time t(c). The evolution also depends on whether a confining interclump medium exists, or whether the interclump region has negligible pressure, as is the case for turbulence-generated clumps. In this paper, we use spherical 1-D numerical hydrodynamic models as well as approximate analytical models to study the dependence of clump photoevaporation on the physical parameters of the clump, and to derive the dynamical evolution, mass loss rates and photoevaporative timescales of a clump for a variety of astrophysical situations. Turbulent clumps evolve so that their column densities are equal to a critical value determined by the local FUV field, and typically have short photo evaporation timescales, approx. 10(exp 4-5) years for a 1 M(solar mass) clump in a typical star-forming region (Nc0 = 10, upsilon = 10). Clumps with insufficient magnetic pressure support, and in strong FUV fields may be driven to collapse by the compressional effect of converging shock waves. We also estimate the rocket effect on photoevaporating clumps and find that it is significant only for the smallest clumps, with sizes much less than the extent of the PDR itself. Clumps that are confined by all interclump medium may either get completely photoevaporated, or may preserve a shielded core with a warm, dissociated, protective shell that, absorbs the incident FUV flux. We compare our results with observations of some well studied PDRs: the Orion Bar, M17SW NGC 2023 and the Rosette Nebula. The data are consistent with both interpretations of clump origin. turbulence and pressure confinement, with a slight indication for favouring the turbulent model for clumps over pressure-confined

Gorti, Uma↗
NASA NTRS2016

Relationships Between Vestibular Measures as Potential Predictors for Spaceflight Sensorimotor Adaptation

Introduction: During space exploration missions astronauts are exposed to a series of novel sensorimotor environments, requiring sensorimotor adaptation. Until adaptation is complete, sensorimotor decrements occur, affecting critical tasks such as piloted landing or docking. Of particularly interest are locomotion tasks such as emergency vehicle egress or extra-vehicular activity. While nearly all astronauts eventually adapt sufficiently, it appears there are substantial individual differences in how quickly and effectively this adaptation occurs. These individual differences in capacity for sensorimotor adaptation are poorly understood. Broadly, we aim to identify measures that may serve as pre-flight predictors of and individual's adaptation capacity to spaceflight-induced sensorimotor changes. As a first step, since spaceflight is thought to involve a reinterpretation of graviceptor cues (e.g. otolith cues from the vestibular system) we investigate the relationships between various measures of vestibular function in humans. Methods: In a set of 15 ground-based control subjects, we quantified individual differences in vestibular function using three measures: 1) ocular vestibular evoked myogenic potential (oVEMP), 2) computerized dynamic posturography and 3) vestibular perceptual thresholds. oVEMP responses are elicited using a mechanical stimuli approach. Computerized dynamic posturography was used to quantify Sensory Organization Tests (SOTs), including SOT5M which involved performing pitching head movements while balancing on a sway-reference support surface with eyes closed. We implemented a vestibular perceptual threshold task using the tilt capabilities of the Tilt-Translation Sled (TTS) at JSC. On each trial, the subject was passively roll-tilted left ear down or right ear down in the dark and verbally provided a forced-choice response regarding which direction they felt tilted. The motion profile was a single-cycle sinusoid of angular acceleration with a duration of 5 seconds (frequency of 0.2 Hz), which was selected as it requires sensory integration of otolith and semicircular canal cues. Stimuli direction was randomized and magnitude was determined using an adaptive sampling procedure. One hundred trials were provided and each subject's responses were fit with a psychometric curve to estimate the subject's threshold. Results: Roll tilt perceptual thresholds at 0.2 Hz ranged from 0.5 degrees to 1.82 degrees across the 15 subjects (geometric mean of 1.04 degrees), consistent with previous studies. The inter-individual variability in thresholds may be able to help explain individual differences observed in sensorimotor adaptation to spaceflight. Analysis is ongoing for the oVEMPS and computerized dynamic posturography to identify relationships between the various vestibular measures. Discussion: Predicting individual differences in sensorimotor adaptation is critical both for the development of personalized countermeasures and mission planning. Here we aim to develop a basis of vestibular tests and parameters which may serve as predictors of individual differences in sensorimotor adaptability through studying the relationship between these measures.

Clark, T. K.↗
NASA NTRS2022

Design of Hypervelocity Flow Generator (HFG) and Its Diagnostics

Ground facilities for hypersonic research are a key element for successful development of aerodynamically proven hypervelocity vehicles. Design concepts and diagnostics of a hypervelocity flow generator (HFG) were made as a test platform for hypersonic/hypervelocity spacecraft models at the NASA Langley Research Center. The HFG is a hypersonic flow field generator using optically heated gas which is blown into an 80 cu m vacuum chamber. The vacuum chamber is kept at a stable vacuum pressure with a combination of three large vacuum pumps, while the HFG is in the test mode. The HFG provides a relatively small test section with approximately a 20 cm window. This facility was designed to generate 2.45 km/sec of flow speed, and potentially generate a continuous flow with the nozzle and vacuum system. The window in test section provides a direct view of the shock wave around a model in order to measure temperature, pressure, and density profile within the shock layer. One of the key test goals under this project is to understand why the emission spectra from the standing shockwave plasma predicted by Lora-Loran codes are significantly different from the measured emission spectra from the Flight Investigation Reentry Environment (FIRE II) Flight. The correct estimate of the thermal loading on the leading edge of hypersonic vehicles greatly affects the aerodynamic design, the material selection for the vehicle, and the cooling requirement and can be obtained by the precise modeling of emission spectra from shockwave plasma. However, the estimation of thermal loading is not an easy task due to complex non-equilibrium radiative process within high temperature shock layers that still falls into a category of cold plasma. Direct flight experiments are the most desirable, but not a cost-effective approach. Analysis by computational fluid dynamics (CFD) offers many test flexibilities. However, the CFD codes must be fully tested and validated with experimental data before the codes are effectively used for practical design. Large discrepancies between experiments and codes appear in hypersonic/hypervelocity flow regimes at high altitudes of 60 km ~ 90 km. This HFG facility offers some important parameters for CFD code validation, such as collision cross-sections, relaxation times, reaction rate coefficients and transportation coefficients. The HFG test facility is based on the ejection flow of high temperature gas heated over to 3500 K through a nozzle. The tungsten gas chamber of the HFG is heated up to a desired temperature by a 60 kW optical power beam source. This system consists of an optical power source, a thermal chamber, an expansion nozzle, a test section, and an 80 cu m vacuum tank. 60 kW optical input power is obtained from the 150 kW Vortek arc lamp system (by Vortek Industries, Vancouver, Canada). This optical beam is focused to heat the gas chamber within which a flow media is heated. The maximum achievable temperature of the flow medium reached approximately 3500 K or even higher but is limited by the melting point of the chamber material used. The exhaust velocity through the nozzle was determined by the stagnation temperature and the molecular weight of the working medium at the test section. To provide design parameters, a NASA Chemical Equilibrium with Application (CEA) computer program is used for the simulation of aerothermal data. This CEA program can calculate chemical equilibrium and properties of complex mixtures using shock tube parameters. For nitrogen gas at 2666 K stagnation temperature, the maximum achievable velocity at the test section is approximately 2.45 km/sec which is within the range of the thermal velocity of 8000 K shock layer. Based on the calculation through the CEA program, the design parameters of a HFG were determined and implemented for the test section that includes an expansion nozzle. The installed test section of the experimental facility can sustain a condition of a re-entry vehicle from the Space at an altitude of 60-90 km. After installation of the HFG, the system was fully tested and its operational parameters were measured. An 80 cu m vacuum chamber of HFG was set at 1 torr level to keep a stable downstream condition. The pumping time to reach the minimum vacuum pressure (~ 1 torr) at the test chamber from atmospheric pressure was approximately an hour using Kinney (MBV-14000/MB -1600/KT-300) and Stokes (Model 1772 and 412) pumps. Such a setup condition allowed a continued stable operation of the HFG experiment with flowrates through 1-, 3-, or 5-mm diameters of nozzle throat. The flow characteristics of the HFG for various operating conditions were performed using a focused Schlieren method. At 0.8 torr chamber pressure, a barrel shock was observed at the test section. With a sphere obstacle of 12 mm diameter in the flow, a bow shock (~ 2-mm thickness) was observed by a focused Schlieren visualization method.

Hypervelocity flow↗
DOE OSTI2019

Machine learning in cardiovascular flows modeling: Predicting arterial blood pressure from non-invasive 4D flow MRI data using physics-informed neural networks

Advances in computational science offer a principled pipeline for predictive modeling of cardiovascular flows and aspire to provide a valuable tool for monitoring, diagnostics and surgical planning. Such models can be nowadays deployed on large patient-specific topologies of systemic arterial networks and return detailed predictions on flow patterns, wall shear stresses, and pulse wave propagation. However, their success heavily relies on tedious pre-processing and calibration procedures that typically induce a significant computational cost, thus hampering their clinical applicability. Here, we put forth a machine learning framework that enables the seamless synthesis of non-invasive in-vivo measurement techniques and computational flow dynamics models derived from first physical principles. We illustrate this new paradigm by showing how one-dimensional models of pulsatile flow can be used to constrain the output of deep neural networks such that their predictions satisfy the conservation of mass and momentum principles. Once trained on noisy and scattered clinical data of flow and wall displacement, these networks can return physically consistent predictions for velocity, pressure and wall displacement pulse wave propagation, all without the need to employ conventional simulators. A simple post-processing of these outputs can also provide a relatively cheap and effective way for estimating Windkessel model parameters that are required for the calibration of traditional computational models. Finally, the effectiveness of the proposed techniques is demonstrated through a series of prototype benchmarks, as well as a realistic clinical case involving in-vivo measurements near the aorta/carotid bifurcation of a healthy human subject.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2020

Molecular dynamics simulations for glass transition temperature predictions of polyhydroxyalkanoate biopolymers

Polyhydroxyalkanoates (PHAs) represent an emerging class of biosynthetic and biodegradable polyesters that exhibit considerable potential to replace petroleum-based plastics towards a sustainable future. Despite the promise, general structure–property mappings within this class of polymers remain largely unexplored. An efficient exploration of this vast chemical space calls for the development and validation of predictive methods for accurate estimation of a diverse range of properties for PHA-based polymers. Towards this aim, we present and validate in this work the results of our molecular dynamics (MD) simulation based approach aimed at predicting glass transition temperatures (T g ) of PHA-based polymers. Since generally available and widely used polymer forcefields exhibit a relatively poor performance for T g predictions, we have developed a new forcefield by modifying the polymer consistent force field (PCFF) via refining a selected set of torsion potentials of the polymer backbone using accurate density functional theory (DFT) computations. After carefully assessing the dependence of critical simulation parameters, such as, polymer chain length, number of polymer chains, supercell size, and thermal quenching rate used in the simulation, the applicability and transferability of the modified PCFF (mPCFF) is demonstrated by directly comparing the computed T g predictions of various polymers with different chemistries, polymer side chain lengths and functional groups forming the polymer side chains against the respective experimentally measured values. Furthermore, the transport properties such as self-diffusion coefficient and viscosity are computationally determined and their well-known correlation with the target properties is demonstrated. Lastly, we have employed the developed approach to predict T g values for a number of yet-to-be-synthesized PHA-based polymers with a diverse set of functional groups in the polymer side chains. The results are further rationalized by correlating the predicted T g values with the inter-chain H-bond formation tendencies of the different side chain functional groups. This work represents an important first step towards computationally guided design of PHA-based functional polymers and opens up new directions for a systematic investigation of composition- and configuration-dependent structure–property relationships in more complex binary and ternary copolymer systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Data-Driven Model for Photovoltaic Generation: Comparison with Physical Models Using a Microgrid in Puerto Rico

Photovoltaic (PV) generation is a critical component of microgrids, but its accurate modeling is challenging due to the complex and dynamic interactions between solar irradiance, temperature, and PV system installation. This paper develops a multilayer perceptron (MLP) model that inputs solar irradiance and temperature to estimate the PV generation, and it compares the proposed data-driven model’s performance to two well-known physical models: the single-diode model and the inverter model. The results demonstrate that all the models can reach high levels of accuracy. However, the MLP model outperforms the physical models on average by 4.5 to 6.6 percent in R squared scores and 220 to 290 Watts in RMSE scores, and it does not require physical system parameters. Moreover, the data-driven model can overcome the limitations of the lack of real-time PV generation data.

R pesante colón, Marcos↗