Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time domain modeling”

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 289 records · Page 16

PLATSIM: An efficient linear simulation and analysis package for large-order flexible systems

PLATSIM is a software package designed to provide efficient time and frequency domain analysis of large-order generic space platforms implemented with any linear time-invariant control system. Time domain analysis provides simulations of the overall spacecraft response levels due to either onboard or external disturbances. The time domain results can then be processed by the jitter analysis module to assess the spacecraft's pointing performance in a computationally efficient manner. The resulting jitter analysis algorithms have produced an increase in speed of several orders of magnitude over the brute force approach of sweeping minima and maxima. Frequency domain analysis produces frequency response functions for uncontrolled and controlled platform configurations. The latter represents an enabling technology for large-order flexible systems. PLATSIM uses a sparse matrix formulation for the spacecraft dynamics model which makes both the time and frequency domain operations quite efficient, particularly when a large number of modes are required to capture the true dynamics of the spacecraft. The package is written in MATLAB script language. A graphical user interface (GUI) is included in the PLATSIM software package. This GUI uses MATLAB's Handle graphics to provide a convenient way for setting simulation and analysis parameters.

Maghami, Periman↗

Precipitation Processes developed during ARM (1997), TOGA COARE(1992), GATE(1 974), SCSMEX(1998) and KWAJEX(1999): Consistent 2D and 3D Cloud Resolving Model Simulations

Real clouds and clouds systems are inherently three dimensional (3D). Because of the limitations in computer resources, however, most cloud-resolving models (CRMs) today are still two-dimensional (2D). A few 3D CRMs have been used to study the response of clouds to large-scale forcing. In these 3D simulations, the model domain was small, and the integration time was 6 hours. Only recently have 3D experiments been performed for multi-day periods for tropical cloud system with large horizontal domains at the National Center for Atmospheric Research. The results indicate that surface precipitation and latent heating profiles are very similar between the 2D and 3D simulations of these same cases. The reason for the strong similarity between the 2D and 3D CRM simulations is that the observed large-scale advective tendencies of potential temperature, water vapor mixing ratio, and horizontal momentum were used as the main forcing in both the 2D and 3D models. Interestingly, the 2D and 3D versions of the CRM used in CSU and U.K. Met Office showed significant differences in the rainfall and cloud statistics for three ARM cases. The major objectives of this project are to calculate and axamine: (1)the surface energy and water budgets, (2) the precipitation processes in the convective and stratiform regions, (3) the cloud upward and downward mass fluxes in the convective and stratiform regions; (4) cloud characteristics such as size, updraft intensity and lifetime, and (5) the entrainment and detrainment rates associated with clouds and cloud systems that developed in TOGA COARE, GATE, SCSMEX, ARM and KWAJEX. Of special note is that the analyzed (model generated) data sets are all produced by the same current version of the GCE model, i.e. consistent model physics and configurations. Trajectory analyse and inert tracer calculation will be conducted to identify the differences and similarities in the organization of convection between simulated 2D and 3D cloud systems.

Tao, W.-K.↗

Dynamic Phasor Modeling of Various Multipulse Rectifiers and a VSI Fed by 18-Pulse Asymmetrical Autotransformer Rectifier Unit for Fast Transient Analysis

To study the fast modelling of dynamic characteristics caused by power electronic switching, the dynamic phasor (DP) theory based on the time-varying Fourier decomposition and frequency shifting is firstly applied to the modelling of various multipulse rectifiers. The DP model for symmetrical 12-pulse phase-shifting reactor rectifier unit (PSR-RU) is derived with model order reduction, relations between ac and dc terminals, and Taylor series expansion. The DP model of asymmetrical 18-pulse autotransformer rectifier unit (AT-RU) is proposed based on the switching functions expressed in the DP domain. Meanwhile, the DP model of voltage source inverter (VSI) fed by asymmetrical 18-pulse AT-RU is built with the harmonic state-space (HSS) equations. Under both balanced and unbalanced conditions, the good calculation accuracy and rapid simulation speed of the developed DP models are validated by the detailed time-domain (TD) simulation.

42 ENGINEERING↗

Multiphysics Modeling of Hydrokinetic Turbine Energy Conversion System

Hydrokinetic turbines (HKTs) hold great promise as a renewable energy source, but high maintenance costs and limited energy output hinder their widespread adoption. The lack of comprehensive research on HKT drivetrain designs creates a knowledge gap in enhancing generation efficiency and cost reduction. A model of the HKT system with a focus on electric drivetrains and power converters is required to address this knowledge gap. This paper first introduces a MATLAB-averaged model integrating electrical-mechanical-thermal domains and aging behaviors within multi-time frames. Then a PLECS model, which incorporates maximum power point tracking and dq reference framed control for AC-DC-AC power converters, is enriched by a dynamic thermal model to predict fast transients accurately. These models optimize the design at the component level, resulting in improved integrated system performance. Furthermore, the validation of the models is carried out through hardware experiments for the averaged model and hardware-in-the-loop testing for the dynamic model.

ADVANCED PROPULSION SYSTEMS,HYDRO ENERGY↗

Adaptive model reduction for continuous systems via recursive rational interpolation

A method for adaptive identification of reduced-order models for continuous stable SISO and MIMO plants is presented. The method recursively finds a model whose transfer function (matrix) matches that of the plant on a set of frequencies chosen by the designer. The algorithm utilizes the Moving Discrete Fourier Transform (MDFT) to continuously monitor the frequency-domain profile of the system input and output signals. The MDFT is an efficient method of monitoring discrete points in the frequency domain of an evolving function of time. The model parameters are estimated from MDFT data using standard recursive parameter estimation techniques. The algorithm has been shown in simulations to be quite robust to additive noise in the inputs and outputs. A significant advantage of the method is that it enables a type of on-line model validation. This is accomplished by simultaneously identifying a number of models and comparing each with the plant in the frequency domain. Simulations of the method applied to an 8th-order SISO plant and a 10-state 2-input 2-output plant are presented. An example of on-line model validation applied to the SISO plant is also presented.

Lilly, John H.↗

Unsteady aerodynamic loads on pitching aerofoils represented by Gaussian body force distributions

The actuator line model (ALM) is an approach commonly used to represent lifting and dragging devices like wings and blades in large-eddy simulations (LES). The crux of the ALM is the projection of the actuator point forces onto the LES grid by means of a Gaussian regularisation kernel. The minimum width of the kernel is constrained by the grid size; however, for most practical applications like LES of wind turbines, this value is an order of magnitude larger than the optimal value that maximises accuracy. This discrepancy motivated the development of corrections for the actuator line, which, however, neglect the effect of unsteady spanwise shed vorticity. In this work we develop a model for the impact of spanwise shed vorticity on the unsteady loading of an aerofoil modelled as a Gaussian body force distribution, where the model is applicable within the regime of unsteady attached flow. The model solution is derived both in the time and frequency domain and features an explicit dependence on the Gaussian kernel width. We verify the model with ALM-LES for both pitch steps and periodic pitching. The model solution is compared with Theodorsen theory and validated with both computational fluid dynamics using body fitted grids and experiment. It is concluded that the optimal kernel width for unsteady aerodynamics is approximately 40 % of the chord. The ALM is able to predict the magnitude of the unsteady loading up to a reduced frequency of 𝑘 ≈ 0.2.

17 WIND ENERGY↗

Dynamic Modeling from Flight Data with Unknown Time Skews

A method for estimating dynamic model parameters from flight data with unknown time skews is described and demonstrated. The method combines data reconstruction, nonlinear optimization, and equation-error parameter estimation in the frequency domain to accurately estimate both dynamic model parameters and the relative time skews in the data. Data from a nonlinear F-16 aircraft simulation with realistic noise, instrumentation errors, and arbitrary time skews were used to demonstrate the approach. The approach was further evaluated using flight data from a subscale jet transport aircraft, where the measured data were known to have relative time skews. Comparison of modeling results obtained from time-skewed and time-synchronized data showed that the method accurately estimates both dynamic model parameters and relative time skew parameters from flight data with unknown time skews.

Morelli, Eugene A.↗

TEAMER: Numerical Modeling and Optimization of the iProTech Pitching Inertial Pump (PIP) Wave Energy Converter (WEC)

This project focused on developing an automated workflow to evaluate and optimize the iProTech Pitching Inertial Pump (PIP) wave energy converter (WEC) using open-source Python packages and the MATLAB/Simulink tool, WEC-Sim. The process involved parameterizing key design variables, running time-domain simulations, and performing sensitivity analyses to determine their impact on power output. The workflow, designed for the PIP device, is generalized and can be extended to optimize other WECs that can be simulated in WEC-Sim. This work establishes a foundation for future time-domain-based WEC design optimizations. Included in this submission are all figures from the final report and the model inputs required to generate them. This includes Python scripts with inputs that produce the meshes, boundary element method (BEM) models, hydrodynamic coefficients, and the WEC-Sim models used for time-domain analyses. Although data for every single run is not included to save space, all of it can be reproduced using the provided models. Detailed instructions for setting up the environment and running the codes are also included.

16 TIDAL AND WAVE POWER↗

Influence of defects on the femtosecond laser damage resistance of multilayer dielectric gratings

Multilayer dielectric (MLD) gratings with high diffraction efficiency and a high laser-induced damage (LID) threshold for pulse compressors are key to scaling the peak and average power of chirped pulse amplification lasers. However, surface defects introduced by manufacturing, storage, and handling processes can reduce the LID resistance of MLD gratings and impact the laser output. The underlying mechanisms of such defect-initiated LID remain unclear, especially in the femtosecond regime. In this Letter, we model dynamic processes in interactions of a 20-fs near-infrared (NIR) laser pulse and a MLD grating design in the presence of cylindrically symmetrical nodules and particle contaminants and cracks at the surface. Utilizing a dynamic model based on a 2D finite difference in time domain (FDTD) field solver coupled with photoionization, electron collision, and refractive index modification, we study the simulation results for the damage site distribution initiated by defects of various types and sizes and its impact on the LID threshold of the grating design.

Zhang, Simin (ORCID:0000000174948230)↗

Applied Routh approximation

The Routh approximation technique for reducing the complexity of system models was applied in the frequency domain to a 16th order, state variable model of the F100 engine and to a 43d order, transfer function model of a launch vehicle boost pump pressure regulator. The results motivate extending the frequency domain formulation of the Routh method to the time domain in order to handle the state variable formulation directly. The time domain formulation was derived and a characterization that specifies all possible Routh similarity transformations was given. The characterization was computed by solving two eigenvalue-eigenvector problems. The application of the time domain Routh technique to the state variable engine model is described, and some results are given. Additional computational problems are discussed, including an optimization procedure that can improve the approximation accuracy by taking advantage of the transformation characterization.

Merrill, W. C.↗

Real-Time Modeling of Thermal Interactions in Cryogenic Ball Bearings

Based on the integration of classical differential equations of motion of bearing elements, a real-time dynamic simulation of ball bearing performance is coupled with simultaneous modeling of thermal interactions between the bearing elements. The transient heat generations are time averaged over a significantly larger time step before they are used to compute the changes in bearing temperature fields. Such an averaging algorithm eliminates the numerical difficulties associated with simultaneous integration of mechanical thermal differential equations with vastly different time scales. Under stable operation, the step change in temperatures at each thermal time step converges to a steady-state solution. Following is a list of accomplishments in the present investigation: (1) A real-time dynamic performance of both 440C and hybrid ball bearings in a LOX environment is modeled in both transient and steady-state time domains, as the inner race accelerates to a prescribed operating speed. (2) The frictional interactions at the ball to race contacts are modeled by experimentally measured traction coefficients as a function of slide-to-roll ratios. (3) Bearing heat generation, as defined by the heat transferred to circulating LOX, is validated against experimental values determined by measured LOX temperatures under prescribed pressures and flow rates. The predicted heat generation is in good agreement with the experimental values. (4) The steady-state solutions are shown to be independent of initial conditions. This establishes convergence of the numerical integration of a differential equation of motion of the bearing elements. (5) Churning and drag losses in the circulating LOX constitutes the majority of bearing heat generation under the experimental test conditions. (6) The heat dissipated in frictional interactions in hybrid bearings is significantly lower in hybrid bearings in comparison to that simulated in all-steel bearings, while the churning and drag losses are relatively unchanged.

Gupta, P. K.↗

Results of an A109 simulation validation and handling qualities study

The results for the validation of a mathematical model of the Agusta A109 helicopter, and subsequent use of the model as the baseline for a handling qualities study of cockpit centerstick requirements, are described. The technical approach included flight test, non-realtime analysis, and realtime piloted simulation. Results of the validation illustrate a time- and frequency-domain approach to the model and simulator issues. The final A109 model correlates well with the actual aircraft with the Stability Augmentation System (SAS) engaged, but is unacceptable without the SAS because of instability and response coupling at low speeds. Results of the centerstick study support the current U.S. Army handling qualities requirements for centerstick characteristics.

Eshow, Michelle M.↗

Results of an A109 simulation validation and handling qualities study

The results for the validation of a mathematical model of the Agusta A109 helicopter, and subsequent use of the model as the baseline for a handling qualities study of cockpit centerstick requirements, are described. The technical approach included flight test, non-realtime analysis, and realtime piloted simulation. Results of the validation illustrate a time- and frequency-domain approach to the model and simulator issues. The final A109 model correlates well with the actual aircraft with the Stability Augmentation System (SAS) engaged, but is unacceptable without the SAS because of instability and response coupling at low speeds. Results of the centerstick study support the current U.S. Army handling qualities requirements for centerstick characteristics.

Eshow, Michelle M.↗

A Numerical Model of Unsteady, Subsonic Aeroelastic Behavior

A method for predicting unsteady, subsonic aeroelastic responses was developed. The technique accounts for aerodynamic nonlinearities associated with angles of attack, vortex-dominated flow, static deformations, and unsteady behavior. The fluid and the wing together are treated as a single dynamical system, and the equations of motion for the structure and flow field are integrated simultaneously and interactively in the time domain. The method employs an iterative scheme based on a predictor-corrector technique. The aerodynamic loads are computed by the general unsteady vortex-lattice method and are determined simultaneously with the motion of the wing. Because the unsteady vortex-lattice method predicts the wake as part of the solution, the history of the motion is taken into account; hysteresis is predicted. Two models are used to demonstrate the technique: a rigid wing on an elastic support experiencing plunge and pitch about the elastic axis, and an elastic wing rigidly supported at the root chord experiencing spanwise bending and twisting. The method can be readily extended to account for structural nonlinearities and/or substitute aerodynamic load models. The time domain solution coupled with the unsteady vortex-lattice method provides the capability of graphically depicting wing and wake motion.

Strganac, Thomas W.↗

Optimum power analysis of a self-reactive wave energy point absorber with mechanically-driven power take-offs

This paper studies a self-reactive ocean wave energy converter (WEC) where energy is absorbed through the relative motion between a floating buoy on the ocean surface and a submerged body. Two types of direct-drive power take off (PTO) systems are examined. One adopts a ball screw system in inverse driving mode with mechanical motion rectifier (MMR) that converts the oscillating bi-directional input heave motion into unidirectional rotation to drive an electromagnetic generator, and another one uses a ball screw system to directly drive the generator (called non-MMR). Dynamic models for both types of PTOs are established and integrated with the overall WEC system model in both the time and frequency domains. Numerical simulation is used to investigate the dynamic performance of both PTOs. The influences of the PTO inerter, which is from the rotational inertia of the flywheel and generator, and the WEC drag damping coefficient are discussed and explored. The analytical closed-form solution for the optimum condition of the non-MMR system is derived based on the criterion of maximum power absorption, and numerical simulation yields the optimum condition for the MMR system. Overall, the results under regular and irregular waves show that the PTO inerter in the MMR system can improve the power absorption for small wave periods and maintain the same performance for large wave periods. The PTO inerter in the non-MMR system shifts the peak frequency. The WEC drag damping negatively influences the power absorption in both PTO systems, especially at large wave periods.

16 TIDAL AND WAVE POWER↗

Tuning Hybrid exciton–Photon Fano Resonances in Two-Dimensional Organic–Inorganic Perovskite Thin Films

As easy-to-grow quantum wells with narrow excitonic features at room temperature, two-dimensional (2D) Ruddleson–Popper perovskites are promising for realizing novel nanophotonic devices based on exciton–photon interactions. Here, we demonstrate a distinct hybrid exciton–photon Fano resonance in (C 4 H 9 NH 3 ) 2 PbI 4 thin films prepared via spin coating. Using a classical coupled-oscillator model and finite-difference time-domain simulations, we link the Fano interference to the coupling of the exciton with the Rayleigh-like scattering of the film microstructure. Combining colloidal plasmonic cavities with the 2D perovskite films, we demonstrate tuning of the Fano resonance. In combination with silver nanoparticles, the exciton–photon Fano interference couples to the in-plane plasmonic modes with indications of Rabi splitting. By creating a nanoparticle on mirror geometry, we address the out-of-plane excitonic component, reaching an intermediate coupling regime. Furthermore, these structures suggest possible photonic targets for biomolecular self-assembly applications.

2D perovskites↗

Wide-range continuous tuning of the thermal conductivity of La0.5Sr0.5CoO3-δ films via room-temperature ion-gel gating

Abstract Solid-state control of the thermal conductivity of materials is of exceptional interest for novel devices such as thermal diodes and switches. Here, we demonstrate the ability to continuously tune the thermal conductivity of nanoscale films of La 0.5 Sr 0.5 CoO 3- δ (LSCO) by a factor of over 5, via a room-temperature electrolyte-gate-induced non-volatile topotactic phase transformation from perovskite (with δ ≈ 0.1) to an oxygen-vacancy-ordered brownmillerite phase (with δ = 0.5), accompanied by a metal-insulator transition. Combining time-domain thermoreflectance and electronic transport measurements, model analyses based on molecular dynamics and Boltzmann transport equation, and structural characterization by X-ray diffraction, we uncover and deconvolve the effects of these transitions on heat carriers, including electrons and lattice vibrations. The wide-range continuous tunability of LSCO thermal conductivity enabled by low-voltage (below 4 V) room-temperature electrolyte gating opens the door to non-volatile dynamic control of thermal transport in perovskite-based functional materials, for thermal regulation and management in device applications.

42 ENGINEERING↗

Terminus: A Versatile Simulator for Space-based Telescopes

Space-based telescopes offer unparalleled opportunities for characterizing exoplanets, solar system bodies, and stellar objects. However, observatories in low-Earth orbits (e.g., Hubble, CHaracterising ExOPlanets Satellite, Twinkle, and an ever-increasing number of cubesats) cannot always be continuously pointed at a target due to Earth obscuration. For exoplanet observations consisting of transit, or eclipse, spectroscopy, this causes gaps in the light curve, which reduces the information content and can diminish the science return of the observation. Terminus, a time-domain simulator, has been developed to model the occurrence of these gaps to predict the potential impact on future observations. The simulator is capable of radiometrically modeling exoplanet observations as well as producing light curves and spectra. Here, Terminus is baselined on the Twinkle mission, but the model can be adapted for any space-based telescope and is especially applicable to those in a low-Earth orbit. Terminus also has the capability to model observations of other targets such as asteroids or brown dwarfs.

47 OTHER INSTRUMENTATION↗