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 343 records · Page 19

Design analysis and computer-aided performance evaluation of shuttle orbiter electrical power system. Volume 1: Summary

Studies were conducted to develop appropriate space shuttle electrical power distribution and control (EPDC) subsystem simulation models and to apply the computer simulations to systems analysis of the EPDC. A previously developed software program (SYSTID) was adapted for this purpose. The following objectives were attained: (1) significant enhancement of the SYSTID time domain simulation software, (2) generation of functionally useful shuttle EPDC element models, and (3) illustrative simulation results in the analysis of EPDC performance, under the conditions of fault, current pulse injection due to lightning, and circuit protection sizing and reaction times.

Source record↗

Linear and nonlinear interpretation of CV-580 lightning data

Numerical models developed for the study of lightning strike data acquired by in-flight aircraft are applied to the data measured on the CV-580. The basic technique used is the three dimensional time domain finite difference solution of Maxwell's equations. Both linear and nonlinear models are used in the analysis. In the linear model, the lightning channel and the aircraft are assumed to form a linear time invariant system. A transfer function technique can then be used to study the response of the aircraft to a given lightning strike current. Conversely, the lightning current can be inferred from the measured response. In the nonlinear model, the conductivity of air in the vicinity of the aircraft is calculated and incorporated into the solution of the Maxwell's equations. The nonlinear model thus simulates corona formation and air breakdown. Results obtained from the models are in reasonable agreement with the measured data. This study provides another validation of the models and increases confidence that the models may be used to predict aircraft response to any general lightning strike.

Ng, Poh H.↗

Modeling Pilot Pulse Control

In this study, behavioral models are developed that closely reproduced pulsive control response of two pilots from the experimental pool using markedly different control techniques (styles) while conducting a tracking task. An intriguing find was that the pilots appeared to: 1) produce a continuous, internally-generated stick signal that they integrated in time; 2) integrate the actual stick position; and 3) compare the two integrations to issue and cease pulse commands. This suggests that the pilots utilized kinesthetic feedback in order to perceive and integrate stick position, supporting the hypothesis that pilots can access and employ the proprioceptive inner feedback loop proposed by Hess' pilot Structural Model. The Pulse Models used in conjunction with the pilot Structural Model closely recreated the pilot data both in the frequency and time domains during closed-loop simulation. This indicates that for the range of tasks and control styles encountered, the models captured the fundamental mechanisms governing pulsive and control processes. The pilot Pulse Models give important insight for the amount of remnant (stick output uncorrelated with the forcing function) that arises from nonlinear pilot technique, and for the remaining remnant arising from different sources unrelated to tracking control (i.e. neuromuscular tremor, reallocation of cognitive resources, etc.).

workload↗

Linearized Poststall Aerodynamic and Control Law Models of the X-31A Aircraft and Comparison with Flight Data

The X-31A aircraft has a unique configuration that uses thrust-vector vanes and aerodynamic control effectors to provide an operating envelope to a maximum 70 deg angle of attack, an inherently nonlinear portion of the flight envelope. This report presents linearized versions of the X-31A longitudinal and lateral-directional control systems, with aerodynamic models sufficient to evaluate characteristics in the poststall envelope at 30 deg, 45 deg, and 60 deg angle of attack. The models are presented with detail sufficient to allow the reader to reproduce the linear results or perform independent control studies. Comparisons between the responses of the linear models and flight data are presented in the time and frequency domains to demonstrate the strengths and weaknesses of the ability to predict high-angle-of-attack flight dynamics using linear models. The X-31A six-degree-of-freedom simulation contains a program that calculates linear perturbation models throughout the X-31A flight envelope. The models include aerodynamics and flight control system dynamics that are used for stability, controllability, and handling qualities analysis. The models presented in this report demonstrate the ability to provide reasonable linear representations in the poststall flight regime.

Stoliker, Patrick C.↗

PyJMAK: An Open-Source Python Toolkit for Modeling Solid-State Metallurgical Phase Transformations

Accurate prediction of metallurgical phase transformations is an essential basis for autonomous optimization and rapid part qualification. Several methods can be used to estimate the evolution of phase fractions such as JMAK kinetics-based models, phase-field models, thermodynamic models, and data-driven machine learning models. Thermodynamic and phase-field-based methodologies solve multiphysics equations requiring numerous calibration parameters and significant computational resources. As a result, the computation domain is limited to a point or on order of micron-meters. The data-driven models rely on large datasets from experiments and simulations. While the JMAK model only provides information about phase fraction evolution, it can predict this evolution in near real-time using thermal history and thermodynamic data without restriction on the domain. JMAK models have been popularly used by researchers to model phase transformations occuring during additive manufacturing or over arbitrary temperature profiles. Commercial proprietary software such as Abaqus and Ansys or closed-source in-house implementations offer the ability to model JMAK based kinetics to predict phase transformation. However, these software packages are not open-source or freely available for use and development in conjunction with manufacturing machines, sensors, and machine learning algorithms. In addition, the use of the model is restricted by a license token. In contrast, given temperature profiles at multiple points in the domain, this Python-based PyJMAK model can compute phase evolution in parallel due to its stand-alone modular, voxel-based structure, and it can be executed on high-performance computing resources without any license restrictions.

Prabhune, Bhagya [Oak Ridge National Laboratory (O↗

High spatial resolution thermal conductivity mapping of SiC/SiC composites

Herein, silicon carbide (SiC) ceramic matrix composites are being investigated as a new generation of fuel cladding materials due to its higher accident tolerance compared to Zircaloy. Characterization of the thermal conductivity of SiC constituents and their interfaces and interphases in SiC/SiC composites is needed as inputs to models of cladding performance. We used time-domain thermoreflectance (TDTR) to map the thermal conductivity with a spatial resolution of 2 µm. The SiC/SiC composite is comprised of Hi-Nicalon Type S fibers in a matrix made by chemical vapor infiltration (CVI matrix). The interphase material is a pyrolytic carbon/SiC multilayer. We report thermal conductivity maps of the SiC fiber and SiC matrix at temperatures of 25 °C, 90 °C, 164 °C, and 250 °C. The fiber has a uniform and isotropic thermal conductivity of 22 W m -1 K -1 at 25 °C; the thermal conductivity of fiber is approximately independent of temperature. The matrix thermal conductivity is not uniform and varies from 50 to 120 W m -1 K -1 at 25 °C; the thermal conductivity of the matrix scales approximately with the inverse of the square root of temperature, i.e., ∝ 1/T 1/2 . The thermal conductivity of the interphase at room temperature is 6 W m -1 K -1 . The calculated spatially averaged thermal conductivity of this SiC/SiC composite at 1000°C is 29 W m -1 K -1 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Design of Energy Storage Reactors for Dc-To-Dc Converters

Two methodical approaches to the design of energy-storage reactors for a group of widely used dc-to-dc converters are presented. One of these approaches is based on a steady-state time-domain analysis of piecewise-linearized circuit models of the converters, while the other approach is based on an analysis of the same circuit models, but from an energy point of view. The design procedure developed from the first approach includes a search through a stored data file of magnetic core characteristics and results in a list of usable reactor designs which meet a particular converter's requirements. Because of the complexity of this procedure, a digital computer usually is used to implement the design algorithm. The second approach, based on a study of the storage and transfer of energy in the magnetic reactors, leads to a straightforward design procedure which can be implemented with hand calculations. An equation to determine the lower-bound volume of workable cores for given converter design specifications is derived. Using this computer lower-bound volume, a comparative evaluation of various converter configurations is presented.

Chen, D. Y.↗

A method for high order linear system reduction and nonlinear system simplification

Least-squares-type algorithms for reducing the order of linear systems in the frequency domain and simplifying nonlinear systems in time domain are developed and demonstrated. The possible model structures are represented as nodes in a tree, and costs along the branches are assigned using the repeated-Gram-Schmidt orthogonalization procedure of Desrochers and Saridis (1980), permitting identification of the optimal n-term model by searching the tree to depth n, with no need for parameter identification. The efficiency and flexibility of the algorithms is shown in applications to the eighth-order linear system studied by Hsia (1972), a three-state eight-nonlinear-term aircraft-dynamics problem, and the related linear-controller problem (Garrard and Jordan, 1977).

Desrochers, A. A.↗

A wide bandwidth electrostatic field sensor for lightning research

Data obtained from UHF Radar observation of direct-lightning strikes to the NASA F-106B airplane have indicated that most of the 690 strikes acquired during direct-strike lightning tests were triggered by the aircraft. As an aid in understanding the triggered lightning process, a wide bandwidth electric field measuring system was designed for the F-106B by implementing a clamped-detection signal processing concept originated at the Air Force Cambridge Research Lab in 1953. The detection scheme combines the signals from complementary stator pairs clamped to zero volts at the exact moment when each stator pair is maximally shielded by the rotor, a process that restores the dc level lost by the charge amplifier. The new system was implemented with four shutter-type field mills located at strategic points on the airplane. The bandwidth of the new system was determined in the laboratory to be from dc to over 100 Hz, whereas past designs had upper limits of 10 Hz to 100 Hz. To obtain the undisturbed electric field vector and total aircraft charge, the airborne field mill system is calibrated by using techniques involving results from ground and flight calibrations of the F-106B, laboratory tests of a metallized model, and a finite-difference time-domain electromagnetic computer code.

Zaepfel, K. P.↗

Investigations into the triggered lightning response of the F106B thunderstorm research aircraft

An investigation has been conducted into the lightning characteristics of the NASA F106B thunderstorm research aircraft. The investigation includes analysis of measured data from the aircraft in the time and frequency domains. Linear and nonlinear computer modelling has also been performed. In addition, new computer tools have been developed, including a new enhanced nonlinear air breakdown model, and a subgrid model useful for analyzing fine details of the aircraft's geometry. Comparison of measured and calculated electromagnetic responses of the aircraft to a triggered lightning environment are presented.

Rudolph, Terence H.↗

Enhanced imaging of reflector antenna surface distortion using microwave holography

Two signal processing techniques are discussed that improve the accuracy of a microwave holographic measurement by removing unwanted signals from the aperture distribution: pattern simulation and subtraction, and time domain filtering. Pattern simulation and subtraction involves modeling unwanted scattering mechanisms and then removing them from the measured far-field data. Measurements taken on a focal point geometry and a Cassegrain geometry at 11 GHz were processed by the holographic analysis system. Pattern simulation and subtraction was applied to both geometries. Surface deformation profiles generated for the Cassegrain antenna by this system were compared to an optical measurement of the main reflector surface.

Gilmore, Sean W.↗

FDTD modeling of thin impedance sheets

Thin sheets of resistive or dielectric material are commonly encountered in radar cross section calculations. Analysis of such sheets is simplified by using sheet impedances. In this paper it is shown that sheet impedances can be modeled easily and accurately using Finite Difference Time Domain (FDTD) methods.

Luebbers, Raymond J.↗

FDTD modeling of thin impedance sheets

Thin sheets of resistive or dielectric material are commonly encountered in radar cross section calculations. Analysis of such sheets is simplified by using sheet impedances. It is shown that sheet impedances can be modeled easily and accurately using Finite Difference Time Domain (FDTD) methods. These sheets are characterized by a discontinuity in the tangential magnetic field on either side of the sheet but no discontinuity in tangential electric field. This continuity, or single valued behavior of the electric field, allows the sheet current to be expressed in terms of an impedance multiplying this electric field.

Luebbers, Raymond↗

Real-Time Parameter Estimation Using Output Error

Output-error parameter estimation, normally a post- ight batch technique, was applied to real-time dynamic modeling problems. Variations on the traditional algorithm were investigated with the goal of making the method suitable for operation in real time. Im- plementation recommendations are given that are dependent on the modeling problem of interest. Application to ight test data showed that accurate parameter estimates and un- certainties for the short-period dynamics model were available every 2 s using time domain data, or every 3 s using frequency domain data. The data compatibility problem was also solved in real time, providing corrected sensor measurements every 4 s. If uncertainty corrections for colored residuals are omitted, this rate can be increased to every 0.5 s.

Grauer, Jared A.↗

NESC GN&C TDT Workshop on 2D Image Motion Optical Transfer Functions, Pointing Performance Analysis, and Requirements

What You Will Learn: The focus is on payload imaging performance due to pointing motion. Some historical background on pointing performance analysis is given. The Optical Transfer Function (OTF) and Modulation Transfer Function (MTF) are defined. The imaging performance due to pointing motion is measured by image motion optical transfer functions (IM OTF). IM OTFs are defined for displacement, smear, and jitter motions, which are all rigorously defined. Deterministic and Statistical IM OTFs are briefly derived and graphically illustrated and compared. The IM OTFs are parameterized by pointing error metrics(PEM), which are means and covariances of displacement, smear, and jitter. Emphasis is on procedures and algorithms to evaluate the image motion optical transfer functions and pointing error metrics. Three procedures are covered, which depend on whether the pointing error data is from time-domain simulation, frequency-domain analysis, or stochastic modeling. A method to evaluate the relative contribution of disturbance sources and to identify the most significant contributors is presented. The presentation includes pertinent discussion of flexible structures and control-structure interaction. No single book can adequately cover this subject, so a book is not required for the course. A list of selected articles, reports, documents, and books is provided for reference and further study. Mathis kept to the minimum necessary to convey principles; lengthy derivations are left to the reference material. Graphics are used to illustrate concepts. As with any such learning endeavor, the knowledge gained will be retained and strengthened through actual practice.c©2019–2021 Mark E. Pittelkau— 5

NASA Engineering and Safety Center (NESC)↗

Simulation-Based Analysis and Prediction of Thrust Vector Servoelastic Coupling

A method of analysis and prediction of servoelastic coupling in launch vehicles is presented, surveying the discovery and subsequent resolution of a predicted servoelastic resonance phenomenon affecting the NASA Space Launch System launch vehicle at specific flight conditions. A physics-based linearized multibody mechanization of the governing equations is combined with first principles analysis to demonstrate that antisymmetric bending of the solid rocket motors leads to a reduction of equivalent viscous modal damping through coupling with the thrust vector control actuators. The sensitivity to parameters and the effects of the resonance phenomenon on flight control performance and stability are confirmed through extensive simulation verification in the time and frequency domain. A novel enhancement in model fidelity that accounts for Coriolis effects of fluid flow on bending within the solid rocket motor case and nozzle is shown to add sufficient damping to reduce the risk of adverse control-structure interaction.

Jeb S Orr↗

Adjoint Waveform Tomography for Crustal and Upper Mantle Structure the Middle East and Southwest Asia for Improved Waveform Simulations Using Openly Available Broadband Data

We present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. We inverted waveforms from 192 Global Centroid Moment Tensor earthquakes (MW 5.5-7.0) recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region. Spatial coverage of the available data is highly uneven due to earthquakes clustered along plate boundaries and sparse coverage of open seismic networks in the region. We considered three possible starting models: the SPiRaL global model (Simmons et al., 2021); MEC-1 (Kaviani et al., 2020); and CSEM2.0 (Noe et al., 2023). Because the SPiRaL model provides good fits to the observed waveforms measured by the time-bandwidth product of selected windows in several period bands, provides all the necessary parameters and covers the entire domain we used it for the starting model with the period band 50-100 seconds. Inversion iterations proceeded using time-frequency phase misfits in six stages and 54 total iterations reducing the minimum period to 30 seconds. Our final model, MESWA, provides improved waveform fits compared to the starting model for both the data used in the inversion and an independent validation data set of 66 events. Two metrics of waveform fit (the time-frequency phase misfit used in the optimization and normalized L2 misfit) were both reduced by nearly 60% for both data sets and MESWA provides significantly larger misfit reductions relative to the SPiRaL model than the MEC-1 or CSEM models. We also find that MESWA provides a larger time-bandwidth product of selected windows indicating that more information content of the observed waveforms is explained by MESWA than the other models. Our new model reveals tectonic features imaged by other studies and methods but in a new holistic model of shear and compressional wavespeeds (v S and v P , respectively) with anisotropy covering the crust and uppermost mantle of a larger domain. MESWA has smaller scale-length features and tends to sharpen some features relative to the SPiRaL starting model. Examples include: low crustal v S in the TurkishIranian Plateau, Zagros Mountains, Afghan Central Blocks and Sulaiman Fold Belt; low mantle vSfollowing divergent (Gulf of Aden, Red Sea) and transform (Dead Sea Fault) margins of the Arabian Plate; low and high v S in the mantle beneath the Arabian Shield and Platform, respectively. Low vS is imaged below Cenozoic volcanic centers of the Arabian Peninsula, the so-called Mecca-Madina-Nafud (MMN) Line. Positive anisotropy (v SH > v SV ) is inferred for asthenospheric depths across the region except where up/downwelling may influence fabric alignment (e.g. Afar, Red Sea, Arabian Shield). Elevated vS tracks Makran subduction under southeast Iran. MESWA resembles the SPiRaL model in its long-wavelength structure, but enhances shorter wavelengths features on the order of 200 km and smaller. The resulting model could be used for as a starting model for further improvements, say using waveforms from in-country seismic networks that are not openly available or smaller-scale studies targeting shorter period waveforms. The model also could be used for source characterization and moment tensor inversion to improve earthquake hazard studies and nuclear explosion monitoring.

58 GEOSCIENCES↗

TPSAS-NF1676L-33433-DND

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for aircraft noise reduction, and are characterized by a complex valued frequency-dependent impedance. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic scattering of geometric bodies treated with acoustic liners. This work considers using either an impedance boundary condition or admittance boundary condition to allow acoustic scattering problems to be modeled with geometries consisting of both un-lined and lined surfaces; admittance is defined to be the inverse of impedance. Three different acoustic liner models will be discussed: the Extended Helmholtz Resonator Model, the Three-Parameter Impedance Model, and a Broadband Model. In the Extended Helmholtz Resonator Model and Three-Parameter Impedance Model, the liner impedance is specified at a single frequency; the Broadband Model allows for the investigation of multiple frequencies simultaneously. The impedance and admittance boundary conditions for each model will be derived and coupled with a time-domain boundary integral equation. The scattering solution will be obtained iteratively using spatial and temporal basis functions. Stability will be demonstrated through eigenvalue analysis. Moreover, fast algorithms and high performance computing are considered to both assess and reduce the computational cost of the simulation.

Michelle E Rodio↗