Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic 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 1,513 records · Page 84

The influence of unsteady aerodynamics on hingeless rotor ground resonance

Calculations of the model frequency and damping for a hingeless rotor on a gimballed support in hover are compared with measured results for two configurations (differing in blade flap stiffness). Good correlation is obtaned when an inflow dynamics model is used to account for the influence of the unsteady aerodynamics. The effect of the unsteady aerodynamics is significant for this rotor system. The inflow dynamics model introduces additional states corresponding to perturbations of the wake-induced velocity at the rotor disk. The calculations confirm the experimental observation that the inflow mode introduced by these additional states is measurable for one configuration but not for the other.

Johnson, W.↗

Status on the Component Models Developed in the Modelica Framework: High-Temperature Steam Electrolysis Plant & Gas Turbine Power Plant

This report has been prepared as part of an effort to design and build a modeling and simulation (M&S) framework to assess the economic viability of a nuclear-renewable hybrid energy system (N-R HES). In order to facilitate dynamic M&S of such an integrated system, research groups in multiple national laboratories have been developing various subsystems as dynamic physics-based components using the Modelica programming language. In fiscal year (FY) 2015, Idaho National Laboratory (INL) performed a dynamic analysis of two region-specific N-R HES configurations, including the gas-to-liquid (natural gas to Fischer-Tropsch synthetic fuel) and brackish water reverse osmosis desalination plants as industrial processes. In FY 2016, INL has developed two additional subsystems in the Modelica framework: a high-temperature steam electrolysis (HTSE) plant and a gas turbine power plant (GTPP). HTSE has been proposed as a high priority industrial process to be integrated with a light water reactor (LWR) in an N-R HES. This integrated energy system would be capable of dynamically apportioning thermal and electrical energy (1) to provide responsive generation to the power grid and (2) to produce alternative industrial products (i.e., hydrogen and oxygen) without generating any greenhouse gases. A dynamic performance analysis of the LWR/HTSE integration case was carried out to evaluate the technical feasibility (load-following capability) and safety of such a system operating under highly variable conditions requiring flexible output. To support the dynamic analysis, the detailed dynamic model and control design of the HTSE process, which employs solid oxide electrolysis cells, have been developed to predict the process behavior over a large range of operating conditions. As first-generation N-R HES technology will be based on LWRs, which provide thermal energy at a relatively low temperature, complementary temperature-boosting technology was suggested for integration with the HTSE process that requires higher temperature input. Simulation results involving several case studies show that the suggested control scheme could maintain the controlled variables (including the steam utilization factor, cathode stream inlet composition, and temperatures of the process streams at various locations) within desired limits under various plant operating conditions. The results also indicate that the proposed HTSE plant could provide operational flexibility to participate in energy management at the utility scale by dynamically optimizing the use of excess plant capacity within an N-R HES. A natural-gas fired GTPP has been proposed as a secondary energy supply to be included in an N-R HES. This auxiliary generator could be used to cover rapid dynamics in grid demand that cannot be met by the remainder of the N-R HES. To evaluate the operability and controllability of the proposed process during transients between load (demand) levels, the dynamic model and control design were developed. Special attention was given to the design of feedback controllers to regulate the power frequency, and exhaust gas and turbine firing temperatures. Several case studies were performed to investigate the system responses to the major disturbance (power load demand) in such a control system. The simulation results show that the performance of the proposed control strategies was satisfactory under each test when the GTPP experienced high rapid variations in the load.

08 HYDROGEN↗

Extended dynamic mode decomposition for model reduction in fluid dynamics simulations

High computational cost and storage/memory requirements of fluid dynamics simulations constrain their usefulness as a predictive tool. Reduced-order models (ROMs) provide a viable solution to this challenge by extracting the key underlying dynamics of a complex system directly from data. We investigate the efficacy and robustness of an extended dynamic mode decomposition (xDMD) algorithm in constructing ROMs of three-dimensional cardiovascular computations. Focusing on the ROMs' accuracy in representation and interpolation, we relate these metrics to the truncation rank of singular value decomposition, which underpins xDMD and other approaches to ROM construction. Our key innovation is to relate the truncation rank to the singular values of the original flow problem. This result establishes a priori guidelines for the xDMD deployment and its likely success as a means of data compression and reconstruction of the system's dynamics from dominant spatiotemporal structures present in the data.

Mechanics↗

Transient and thermo-economic analysis of MED-MVC desalination system

An exergo-economic model is used to assess the performance of a multi-effect desalination plant integrated to a mechanical vapor compressor unit (MED-MVC) with a water production capacity of 1500 $m^3$/day. The results show that the second law efficiency ($\eta_{II}$) 2.8%. The MVC and evaporator units are responsible for about 39 and 52% of the total exergy destruction, respectively. The total water price (TWP) is 1.70 $\$/m^3$ when calculated using a simple conventional economic model and 1.63 $\$/m^3$ when calculated using an exergy-based cost model. Increasing the number of effects from 1 to 6 results in a 39% reduction in the specific power consumption (SPC), a 70% increase in $\eta_{II}$ and a 24% decrease in TWP. A dynamic model is developed to investigate the effect of fluctuations of compressor work ($\dot{W}_c$) and inlet seawater temperature ($T_{sw}$) on the plant behavior and performance. The dynamic model results show that the disturbance in $\dot{W}$ has a significant effect on the plant transient behavior and may cause the plant to cease operation while a disturbance in $T_{sw}$ has only a moderate impact. Increasing $T_{sw}$ above a certain value of the steady-state condition without proper control on the plant response could lead to evaporator dry out. In term of performance, a reduction in $\dot{W}_c$ causes a decrease in the plant production capacity and SPC, while it increases the plant performance ratio (PR). On the other hand, a reduction in the inlet $T_{sw}$ causes a reduction in the plant production capacity and PR and an increase in SPC for the same compressor work. Furthermore, a comparison between a MED-MVC system and a MED integrated to a thermal vapor compressor system (MED-TVC) reveals that the latter system is rather sensitive to the reduction in $T_{sw}$ due to the presence of the condenser unit in the MED-TVC. The response of the MED-MVC system is slower than the MED-TVC which is due to the high thermal capacity of the preheaters for the feed in the MED-MVC.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Quasi‐dynamic breathing model of the lung incorporating viscoelasticity of the lung tissue

Abstract We advanced a novel model to calculate viscoelastic lung compliance and airflow resistance in presence of mucus, accounting for the quasi‐linear viscoelastic stress–strain response of the parenchyma (alveoli) tissue. We adapted a continuum‐based numerical modeling approach for the lung, integrating the fluid mechanics of the airflow within individual generations of the bronchi and alveoli. The model accounts for elasticity of the deformable bronchioles, resistance to airflow due to the presence of mucus within the bronchioles, and subsequent mucus flow. Simulated quasi‐dynamic inhalation and expiration cycles were used to characterize the net compliance and resistance of the lung, considering the rheology of the mucus and viscoelastic properties of the parenchyma tissue. The structure and material properties of the lung were identified to have an important contribution to the lung compliance and airflow resistance. The secondary objective of this work was to assess whether a higher frequency and smaller volume of harmonic air flow rate compared to a normal ventilator breathing cycle enhanced mucus outflow. Results predict, lower mucus viscosity and higher excitation frequency of breathing are favorable for the flow of mucus up the bronchi tree, towards the trachea.

60 APPLIED LIFE SCIENCES↗

Reduced evolutionary rates in HIV-1 reveal extensive latency periods among replicating lineages

Background HIV-1 can persist for the duration of a patient's life due in part to its ability to hide from the immune system, and from antiretroviral drugs, in long-lived latent reservoirs. Latent forms of HIV-1 may also be disproportionally involved in transmission. Thus, it is important to detect and quantify latency in the HIV-1 life cycle. Results We developed a novel molecular clock–based phylogenetic tool to investigate the prevalence of HIV-1 lineages that have experienced latency. The method removes alternative sources that may affect evolutionary rates, such as hypermutation, recombination, and selection, to reveal the contribution of generation-time effects caused by latency. Our method was able to recover latent lineages with high specificity and sensitivity, and low false discovery rates, even on relatively short branches on simulated phylogenies. Applying the tool to HIV-1 sequences from 26 patients, we show that the majority of phylogenetic lineages have been affected by generation-time effects in every patient type, whether untreated, elite controller, or under effective or failing treatment. Furthermore, we discovered extensive effects of latency in sequence data ( gag, pol , and env ) from reservoirs as well as in the replicating plasma population. To better understand our phylogenetic findings, we developed a dynamic model of virus-host interactions to investigate the proportion of lineages in the actively replicating population that have ever been latent. Assuming neutral evolution, our dynamic modeling showed that under most parameter conditions, it is possible for a few activated latent viruses to propagate so that in time, most HIV-1 lineages will have been latent at some time in their past. Conclusions These results suggest that cycling in and out of latency plays a major role in the evolution of HIV-1. Thus, no aspect of HIV-1 evolution can be fully understood without considering latency - including treatment, drug resistance, immune evasion, transmission, and pathogenesis.

Virology↗

Dynamic hysteretic sensing model of bending-mode Galfenol transducer

A dynamic hysteretic sensing model has been developed to predict the dynamic responses of the magnetic induction, the stress, and the output voltage for a bending-mode Galfenol unimorph transducer subjected simultaneously to acceleration and bias magnetic field. This model is obtained by coupling the hysteretic Armstrong model and the structural dynamic model of the Galfenol unimorph beam. The structural dynamic model of the beam is founded based on the Euler-Bernouli beam theory, the nonlinear constitutive equations, and the Faraday law of electromagnetic induction. Comparisons between the calculated and measured results show the model can describe dynamic nonlinear voltage characteristics of the device, and can predict hysteretic behaviors between the magnetic induction and the stress. Moreover, the model can effectively analyze the effects of the bias magnetic field, the acceleration amplitude, and frequency on the root mean square voltage of the device.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mode specific dynamics of the H{sub 2} + CH{sub 3} → H + CH{sub 4} reaction studied using quasi-classical trajectory and eight-dimensional quantum dynamics methods

An eight-dimensional quantum dynamical model is proposed and applied to the title reaction. The reaction probabilities and integral cross sections have been determined for both the ground and excited vibrational states of the two reactants. The results indicate that the H{sub 2} stretching and CH{sub 3} umbrella modes, along with the translational energy, strongly promote the reactivity, while the CH{sub 3} symmetric stretching mode has a negligible effect. The observed mode specificity is confirmed by full-dimensional quasi-classical trajectory calculations. The mode specificity can be interpreted by the recently proposed sudden vector projection model, which attributes the enhancement effects of the reactant modes to their strong couplings with the reaction coordinate at the transition state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coupling fast fluid dynamics and multizone airflow models in Modelica Buildings library to simulate the dynamics of HVAC systems

Multizone models are widely used in building airflow and energy performance simulations due to their fast computing speed. However, multizone models assume that the air in a room is well mixed, consequently limiting their application. In specific rooms where this assumption fails, the use of computational fluid dynamics (CFD) models may be an alternative option. Previous research has mainly focused on coupling CFD models and multizone models to study airflow in large spaces. While significant, most of these analyses did not consider the coupled simulation of the building airflow with the building’s Heating, Ventilation, and Air-Conditioning (HVAC) systems. This paper tries to fill the gap by integrating the models for HVAC systems with coupled multizone and CFD simulations for airflows, using the Modelica simulation platform. To improve the computational efficiency, we incorporated a simplified CFD model named fast fluid dynamics (FFD). We first introduce the data synchronization strategy and implementation in Modelica. Then, we verify the implementation using two case studies involving an isothermal and a non-isothermal flow by comparing model simulations to experiment data. Afterward, we study another three cases that are deemed more realistic. Furthermore, this is done by attaching a variable air volume (VAV) terminal box and a VAV system to previous flows to assess the capability of the models in studying the dynamic control of HVAC systems. Finally, we discuss further research needs on the coupled simulation using the models.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Pseudo-viscous modeling of transport in dense granular flows for thermal energy storage applications

Dense, granular flows were examined to effectively capture and model bulk viscous properties in thin packed beds. A modified Couette cell with particle image velocimetry was used to experimentally determine pseudo-viscosity properties of four particulate media with varying morphologies: (1) iron oxide-coated SiO 2 particles, (2) CARBOBEAD CP30-60 particles, (3) CARBOBEAD CP40-100 particles, and (4) Al 2 O 3 beads. The pseudo-viscosity functions were fitted using a power law to correlate the measured shear stress as a function of measured shear rate. The pseudo-viscous functions were used as inputs to computation fluid dynamics models for a single-phase viscous fluid to predict granular flow profiles. Steady-state free surface velocity profiles at angular velocities <7 rad/s predicted by the model were in good agreement with the experimental particle image velocimetry measurements, resulting in Pearson correlation coefficients of 0.97 for iron-oxide coated SiO 2 particles and 0.95 for CP30-60 particles. As a result, this alternative approach to measuring pseudo-viscous properties under shearing and modeling bulk transport behavior of granular flow using computation fluid dynamics model offered significant reduction in computational load compared to discrete element methods.

14 SOLAR ENERGY↗

SPAR analysis of LDEF vibration characteristics

This paper presents a structural dynamic modeling of the Long Duration Exposure Facility (LDEF), which is a Space Shuttle payload of passive scientific experiments contained in trays mounted on a large cylindrically shaped structure. Special detailed finite element modeling, using the SPAR system of computer programs was required to obtain good agreement between analytical and test vibrations modes. Experimental trays contributed significantly to overall LDEF stiffnesses, and these contributions were realistically represented for each tray by the stiffness matrix of an equivalent orthotropic panel in the overall LDEF SPAR model. Orthotropic stiffnesses for this panel were obtained from finely detailed statically loaded tray SPAR models in which stiffness coupling was accounted for along with partial relative sliding allowed by the tray clamping attachments. Joint boundary conditions were also significant in the structural dynamic modeling of LDEF, and static data proved valuable in assessing modeling of local end fittings.

Edighoffer, H. H.↗

Carbon cycle instability as a cause of the late Pleistocene ice age oscillations - Modeling the asymmetric response

A dynamical model of the Pleistocene ice ages is presented, which incorporates many of the qualitative ideas advanced recently regarding the possible role of ocean circulation, chemistry, temperature, and productivity in regulating long-term atmospheric carbon dioxide variations. This model involves one additional term (and free parameter) beyond that included in a previous model (Saltzman and Sutera, 1987), providing the capacity for an asymmetric response. It is shown that many of the main features exhibited by the delta(O-18)-derived ice record and the Vostok core/delta(C-13)-derived carbon dioxide record in the late Pleistocene can be deduced as a free oscillatory solution of the model.

Saltzman, Barry↗

LDEF environment modeling updates

An updated gas dynamics model for gas interactions around the LDEF is presented that includes improved scattering algorithms. The primary improvement is more accurate predictions of surface fluxes in the wake region. The code used is the Integrated Spacecraft Environments Model (ISEM). Additionally, initial results of a detailed ISEM prediction model of the Solar Array Passive LDEF Experiment (SAMPLE), A0171, is presented. This model includes details of the A0171 geometry and outgassing characteristics of the many surfaces on the experiment. The detailed model includes the multiple scattering that exists between the ambient atmosphere, LDEF outgassing, and atomic oxygen erosion products. Predictions are made for gas densities, surface fluxes and deposition at three different time periods of the LDEF mission.

Gordon, Tim↗

Kalman Orbit Optimized Loop Tracking

Under certain conditions of low signal power and/or high noise, there is insufficient signal to noise ratio (SNR) to close tracking loops with individual signals on orbiting Global Navigation Satellite System (GNSS) receivers. In addition, the processing power available from flight computers is not great enough to implement a conventional ultra-tight coupling tracking loop. This work provides a method to track GNSS signals at very low SNR without the penalty of requiring very high processor throughput to calculate the loop parameters. The Kalman Orbit-Optimized Loop (KOOL) tracking approach constitutes a filter with a dynamic model and using the aggregate of information from all tracked GNSS signals to close the tracking loop for each signal. For applications where there is not a good dynamic model, such as very low orbits where atmospheric drag models may not be adequate to achieve the required accuracy, aiding from an IMU (inertial measurement unit) or other sensor will be added. The KOOL approach is based on research JPL has done to allow signal recovery from weak and scintillating signals observed during the use of GPS signals for limb sounding of the Earth s atmosphere. That approach uses the onboard PVT (position, velocity, time) solution to generate predictions for the range, range rate, and acceleration of the low-SNR signal. The low- SNR signal data are captured by a directed open loop. KOOL builds on the previous open loop tracking by including feedback and observable generation from the weak-signal channels so that the MSR receiver will continue to track and provide PVT, range, and Doppler data, even when all channels have low SNR.

Young, Lawrence E.↗

Multiscale Cloud System Modeling

The central theme of this paper is to describe how cloud system resolving models (CRMs) of grid spacing approximately 1 km have been applied to various important problems in atmospheric science across a wide range of spatial and temporal scales and how these applications relate to other modeling approaches. A long-standing problem concerns the representation of organized precipitating convective cloud systems in weather and climate models. Since CRMs resolve the mesoscale to large scales of motion (i.e., 10 km to global) they explicitly address the cloud system problem. By explicitly representing organized convection, CRMs bypass restrictive assumptions associated with convective parameterization such as the scale gap between cumulus and large-scale motion. Dynamical models provide insight into the physical mechanisms involved with scale interaction and convective organization. Multiscale CRMs simulate convective cloud systems in computational domains up to global and have been applied in place of contemporary convective parameterizations in global models. Multiscale CRMs pose a new challenge for model validation, which is met in an integrated approach involving CRMs, operational prediction systems, observational measurements, and dynamical models in a new international project: the Year of Tropical Convection, which has an emphasis on organized tropical convection and its global effects.

Tao, Wei-Kuo↗

Unique Challenges for Modeling Defect Dynamics in Concentrated Solid-Solution Alloys

Recently developed concentrated solid solution alloys (CSAs) are shown to have improved irradiation performance that strongly depends on the number of alloying elements, alloying species and their concentrations. In contrast to conventional dilute alloys, CSAs are comprised of multiple principal elements randomly situated in a simple crystalline lattice. As a result, the intrinsic disorder has a profound influence on the energy dissipation pathways and defect evolution when these CSAs are subjected to ion irradiation. Extraordinary irradiation resistance, including suppression of void formation by two orders magnitude at an elevated temperature, has been achieved with increasing compositional complexity in CSAs. Yet, the intrinsic disorder poses great challenges to theoretical modelling at the electronic and atomic level. Based on recent computer simulation results from a novel set of Ni-based face-centered cubic CSAs, we demonstrate in this review that how the disorder can be tackled in theoretical modeling and what is the impact of disorder on the defect dynamics. Finally, the unique challenges of describing the energy dissipation and defect evolution in CSAs are discussed.

36 MATERIALS SCIENCE↗

Active learning of constitutive relation from mesoscopic dynamics for macroscopic modeling of non-Newtonian flows

We simulate complex fluids by means of an on-the-fly coupling of the bulk rheology to the underlying microstructure dynamics. In particular, a continuum model of polymeric fluids is constructed without a pre-specified constitutive relation, but instead it is actively learned from mesoscopic simulations where the dynamics of polymer chains is explicitly computed. To couple the bulk rheology of polymeric fluids and the microscale dynamics of polymer chains, the continuum approach (based on the finite volume method) provides the transient flow field as inputs for the (mesoscopic) dissipative particle dynamics (DPD), and in turn DPD returns an effective constitutive relation to close the continuum equations. In this multiscale modeling procedure, we employ an active learning strategy based on Gaussian process regression (GPR) to minimize the number of expensive DPD simulations, where adaptively selected DPD simulations are performed only as necessary. Numerical experiments are carried out for flow past a circular cylinder of a non-Newtonian fluid, modeled at the mesoscopic level by bead-spring chains. The results show that only five DPD simulations are required to achieve an effective closure of the continuum equations at Reynolds number RE = 10. Furthermore, when Re is increased to 100, only one additional DPD simulation is required for constructing an extended GPR-informed model closure. Compared to traditional message-passing multiscale approaches, applying an active learning scheme to multiscale modeling of non-Newtonian fluids can significantly increase the computational efficiency. Although the method demonstrated here obtains only a local viscosity from the polymer dynamics, it can be extended to other multiscale models of complex fluids whose macro-rheology is unknown.

Computer Science↗

Dynamics of the abrupt change in Pacific Plate motion around 50 million years ago

Changes in Pacific Plate motion combined near equally with hotspot drift to generate the prominent bend in the Hawaiian-Emperor seamount chain some 50 million years ago, according to kinematic plate reconstruction and global dynamic models. A drastic change in plate tectonics and mantle convection occurred around 50 Ma as exemplified by the prominent Hawaiian-Emperor Bend. Both an abrupt Pacific Plate motion change and a change in mantle plume dynamics have been proposed to account for the Hawaiian-Emperor Bend, but debates surround the relative contribution of the two mechanisms. Here we build kinematic plate reconstructions and high-resolution global dynamic models to quantify the amount of Pacific Plate motion change. We find Izanagi Plate subduction, followed by demise of the Izanagi-Pacific Ridge and Izu-Bonin-Mariana subduction initiation alone, is incapable of causing a sudden change in plate motion, challenging the conventional hypothesis on the mechanisms of Pacific Plate motion change. Instead, Palaeocene slab pull from Kronotsky intraoceanic subduction in the northern Pacific exerts a northward pull on the Pacific Plate, while its Eocene demise leads to a sudden 30-35 degrees change in plate motion, accounting for about half of the Hawaiian-Emperor Bend. We suggest the Pacific Plate motion change and hotspot drift due to plume dynamics could have contributed nearly equally to the formation of the Hawaiian-Emperor Bend. Such a scenario is consistent with available constraints from global plate circuits, palaeomagnetic data and geodynamic models.

Hu, Jiashun↗