Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation decomposition”

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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Subharmonic resonance in a mixing layer

The subharmonic resonance phenomenon in a spatially evolving mixing layer is studied using direct simulations of the 2-D Navier-Stokes equations. The computational domain extends to + or - infinity in the cross-stream direction with U(1) = 1.25 and U(2) = 0.25 imposed at + infinity and - infinity respectively. The domain is finite in the streamwise direction with inflow and outflow boundary conditions imposed at x/delta sub omega = 0 and 100, respectively. A hyperbolic-tangent mean velocity profile is assumed at the inlet and the Reynolds number based on the inlet vorticity thickness and velocity difference is Re = 600. It is observed that the phase angle between the fundamental and its subharmonic plays a key role in the spatial development of these modes. Contour plots of vorticity show that varying the phase will have a dramatic effect on the dynamics of the vortices. Pairing or shredding is observed depending on the phase. Fourier decomposition of the time traces show that the fundamental grows, saturates and decays with the downstream distance. The subharmonic has a similar behavior. However, the level at which the modes will saturate is affected by the phase. At 0 deg phase, it was found that as the fundamental saturates, the growth rate of the subharmonic is enhanced. At 90 deg phase, it was found that as the fundamental saturates, the growth rate of the subharmonic is inhibited. In the later case, the growth rate of the subharmonic recovers after saturation of the fundamental. These results are in qualitative agreement with experimental data.

Mansour, N. N.↗
DOE OSTI2020

Some observations concerning "laminarization" in heated vertical tubes

For low "turbulent" Reynolds numbers in strongly-heated vertical gas flow, fundamental results of existing direct numerical simulation (DNS) databases have been examined to deduce differences between cases which "revert" to turbulent and those that yield integral parameters which correspond to laminar flow (henceforth laminarizing or laminarized). The objectives of this paper are (1) to examine the streamwise evolution of near-wall flow structures and (2) to determine which, if any, proposed laminarization parameters could be used to discriminate between turbulent and laminarizing flows resulting. Views of streamlines in r-θ cross sections showed evidence of a ring of irregular-shaped streamwise vortices near the wall at all locations. The transverse advection of these vortices is hypothesized to provide a path of least resistance to the transport of streamwise momentum in the wall-normal direction. It is demonstrated that a steady streamwise laminar vortex can provide apparent turbulent momentum transport (as a consequence of the Reynolds decomposition) such that a reasonable mean velocity profile is predicted near the wall. For the present cases, the values of the non-dimensional radius (y c,w + = Re τ,w ) and of the modified wall Reynolds number did discriminate whether turbulent or laminarizing flow resulted.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Sensitivity analysis applied to SiC failure probability in TRISO modeled with BISON

Here, a sensitivity analysis on the failure probability of the Silicon Carbide (SiC) layer in tristructural isotropic (TRISO) nuclear fuel during transient conditions predicted by the BISON fuel performance code is performed. The principal goal of the analysis is to understand the most important parameters dictating SiC failure behavior in BISON during Reactivity Initiated Accidents (RIAs). SiC brittle fracture probability is modeled using Weibull statistics. A total of seven inputs related to SiC failure has been selected for the analysis, including the Weibull statistics parameters, elastic moduli for SiC and Pyrolitic Carbon (PyC) and SiC stress-free temperature. A 1D TRISO BISON model has been established for various reactivity insertions performed at the Nuclear Safety Research Reactor (NSRR). The principal advantage associated with the 1D TRISO model developed in this work is its computational efficiency. The Sobol variance decomposition method is used, and the sensitivity indices are presented for eight different values of the energy deposition. The results show that the two most important parameters impacting the predicted SiC failure probability are the Weibull modulus and the characteristic stress, and a co-variance amongst these parameters is obtained for low reactivity insertions. An additional new finding of this work is that the relative importance of Weibull parameters depends on the energy deposition, and thus reactivity, regime. For low energy depositions, two parameters are of influence on SiC failure probability, while for high energy depositions only one parameter impacts failure probability results. Moreover, optimization of the Weibull modulus and characteristic stress is performed by minimizing the RMSE between BISON failure probability predictions and experimental failure fractions for each energy deposition. This work also demonstrates the validity of the NSRR tests BISON simulations and of the respective sensitivity analysis results as conservative, yet indicative for slower transients characterized by lower deposited energies. Such verification is achieved through the partial extension of the analysis to a group Control Rod Withdrawal reproduced from a previous study. Another novel result of this analysis is that no single set of Weibull parameters can reproduce all reactivity insertion experimental failure results, which is related to the intrinsic nature of SiC failure properties, and a new range for the parameters is proposed to produce a failure probability envelope that encompasses the experimental fractions. Additionaly, this work proposes a new approach for future failure analysis with BISON, consisting in the use of Weibull parameters ranges, rather than fixed sets, along with failure envelopes generation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Modeling Heat Transfer and Reaction Kinetics of Biomass in Pyrolysis Feeding Systems

Biomass feedstocks can be used to produce low carbon-intensity chemicals and liquid fuels by displacing fossil-fuel sources. Pyrolysis is the thermal decomposition of biomass in the absence of oxygen and can be used to produce crude biofuels. A typical feeding mechanism for pyrolysis reactors is by lock hopper followed by a horizontal auger feeder. Particle agglomeration and plugging of the auger screw are regularly occurring phenomena that is detrimental to the process since it results in the need for periodic inspection and cleaning. The underlying mechanisms leading to agglomeration and plugging are unclear, but a working hypothesis is that heat from the reactor raises the temperature of the biomass in the feeder to the point at which preliminary decomposition reactions occur and produce "sticky" products. Using principles of heat and mass transfer, we derived one-dimensional differential equations for heat flow from the pyrolysis reactor through the auger feeder. Temperature profiles were solved from the system of one-dimensional differential equations using the shooting method. Subsequently, we used the temperature profiles as a one-way coupling to a reaction model for the biomass in the feeder. Kinetic models were adopted from the literature and used to predict the formation of intermediates, products, char, extractives, and metaplastics as well as phase changes along the auger feeder. Species were identified that correlated with softening of the biomass, causing it to become sticky and plug the feeder. We will discuss the details of this work, our findings, and the direction of future work that includes fully coupled transport and reaction-kinetics simulations.

BIOMASS FUELS↗
NASA NTRS2018

A Singular Value Decomposition Framework for Retrievals with Vertical Distribution Information from Greenhouse Gas Column Absorption Spectroscopy Measurements

We review the singular value decomposition (SVD) framework and use it for quantifying and discerning vertical information in greenhouse gas retrievals from column integrated absorption measurements. While the commonly used traditional Bayesian optimal estimation (OE) assumes a prior distribution in order to regularize the inversion problem, the SVD approach identifies principal components that can be retrieved from the measurement without explicitly specifying a prior mean and prior covariance matrix. We review the SVD method, explicitly recognize the use of an uninformative prior and show it to incur no bias from the choice of the prior. We also make the connection between the SVD method and the pseudo-inverse, which makes it more intuitive and easy to understand. We illustrate the use of the SVD method on an integrated path differential absorption CO2 lidar measurement model and verify our derivations and bias-free properties versus optimal estimation using numerical simulations. In contrast, traditional OE retrievals exhibit bias when the prior mean used in the retrieval differs from the true mean. Hence, the SVD method is particularly useful for situations in which knowledge of the prior mean and prior covariance of the true state (e.g., greenhouse gas profiles) is inadequate.

Ramanathan, Anand K.↗
NASA NTRS2015

Thermochemical Ablation Analysis of the Orion Heatshield

The Orion Multi-Purpose Crew Vehicle will one day carry astronauts to the Moon and beyond, and Orion's heatshield is a critical component in ensuring their safe return to Earth. The Orion heatshield is the structural component responsible for absorbing the intense heating environment caused by re-entry to Earth's atmosphere. The heatshield is primarily composed of Avcoat, an ablative material that is consumed during the re-entry process. Ablation is primarily characterized by two processes: pyrolysis and recession. The decomposition of in-depth virgin material is known as pyrolysis. Recession occurs when the exposed surface of the heatshield reacts with the surrounding flow. The Orion heatshield design was changed from an individually filled Avcoat honeycomb to a molded block Avcoat design. The molded block Avcoat heatshield relies on an adhesive bond to keep it attached to the capsule. In some locations on the heatshield, the integrity of the adhesive bond cannot be verified. For these locations, a mechanical retention device was proposed. Avcoat ablation was modelled in CHAR and the in-depth virgin material temperatures were used in a Thermal Desktop model of the mechanical retention device. The retention device was analyzed and shown to cause a large increase in the maximum bondline temperature. In order to study the impact of individual ablation modelling parameters on the heatshield sizing process, a Monte Carlo simulation of the sizing process was proposed. The simulation will give the sensitivity of the ablation model to each of its input parameters. As part of the Monte Carlo simulation, statistical uncertainties on material properties were required for Avcoat. Several properties were difficult to acquire uncertainties for: the pyrolysis gas enthalpy, non-dimensional mass loss rate (B´c), and Arrhenius equation parameters. Variability in the elemental composition of Avcoat was used as the basis for determining the statistical uncertainty in pyrolysis gas enthalpy and B´c. A MATLAB program was developed to allow for faster, more accurate and automated computation of Arrhenius reaction parameters. These parameters are required for a material model to be used in the CHAR ablation analysis program. This MATLAB program, along with thermogravimetric analysis (TGA) data, was used to generate uncertainties on the Arrhenius parameters for Avcoat. In addition, the TGA fitting program was developed to provide Arrhenius parameters for the ablation model of the gap filler material, RTV silicone.

Sixel, William↗
NASA NTRS2011

Analysis of the Radiative Transfer Equation with Highly Asymmetric Phase Function

This paper considers a scalar radiative transfer problem with high scattering anisotropy, Two computational methods are presented based on decomposition of the diffuse light field into a regular and anisotropic part. The first algorithm (DOMAS) singles out the anisotropic radiance in the forward scattering peak using the Small-Angle Modification of RTE. The second algorithm (DOM2+) separates the single scattering radiance as an anisotropic part, which largely defines the fine detail of the total radiance in the backscattering directions. In both cases, the anisotropic part is represented analytically. With anisotropy subtraction, the regular part of the signal. which requires a numerical solution, is essentially smoothed as a function of angles. Further, the transport equation is obtained for the regular part that contains an additional source function from the anisotropic part of the signal. This equation is solved with the discrete ordinates method. A conducted numerical analysis of this work showed that algorithm DOMAS has a strong advantage as compared to the standard discrete ordinates method for simulation of the radiance transmission, and DOM2 + is the best of the three for the reflection computations. Both algorithms offer at least a factor of three acceleration of convergence of the azimuthal series for highly anisotropic phase functions.

Korkin, Sergey V.↗
DOE OSTI2022

The eROSITA Final Equatorial-Depth Survey (eFEDS)

Here, we investigate the physical properties – such as the stellar mass (M * ), star-formation rate, infrared (IR) luminosity (L IR ), X-ray luminosity (L X ), and hydrogen column density (N H ) – of mid-IR (MIR) galaxies and active galactic nuclei (AGN) at z < 4 in the 140 deg 2 field observed by eROSITA on SRG using the Performance-and-Verification-Phase program named the eROSITA Final Equatorial Depth Survey (eFEDS). By cross-matching the WISE 22 μm (W4)-detected sample and the eFEDS X-ray point-source catalog, we find that 692 extragalactic objects are detected by eROSITA. We have compiled a multiwavelength dataset extending from X-ray to far-IR wavelengths. We have also performed (i) an X-ray spectral analysis, (ii) spectral-energy-distribution fitting using X-CIGALE, (iii) 2D image-decomposition analysis using Subaru Hyper Suprime-Cam images, and (iv) optical spectral fitting with QSFit to investigate the AGN and host-galaxy properties. For 7088 WISE 22 μm objects that are undetected by eROSITA, we have performed an X-ray stacking analysis to examine the typical physical properties of these X-ray faint and probably obscured objects. We find that (i) 82% of the eFEDS–W4 sources are classified as X-ray AGN with log L X > 42 erg s –1 ; (ii) 67 and 24% of the objects have log(L IR /L ⊙ ) > 12 and 13, respectively; (iii) the relationship between L X and the 6 μm luminosity is consistent with that reported in previous works; and (iv) the relationship between the Eddington ratio and N H for the eFEDS–W4 sample and a comparison with a model prediction from a galaxy-merger simulation indicates that approximately 5.0% of the eFEDS–W4 sources in our sample are likely to be in an AGN-feedback phase, in which strong radiation pressure from the AGN blows out the surrounding material from the nuclear region. Thanks to the wide area coverage of eFEDS, we have been able to constrain the ranges of the physical properties of the WISE 22 μm-selected sample of AGNs at z < 4, providing a benchmark for forthcoming studies on a complete census of MIR galaxies selected from the full-depth eROSITA all-sky survey.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A solid electrolyte interphase to protect the sulfurized polyacrylonitrile (SPAN) composite for Li–S batteries: computational approach addressing the electrolyte/SPAN interfacial reactivity

This study addresses the reactivity of multiple solvents and lithium bis(fluorosulfonyl)imide (LiFSI) at the interface with sulfurized polyacrylonitrile (SPAN) in multiple stages of lithiation via ab initio molecular dynamics simulations. Here, both ether 1,3-dioxolane (DOL) and dimethyl carbonate (DMC) proved stable on the lithiated SPAN surface regardless of the lithium content, meaning that neither of these species likely contributes to growing a solid electrolyte interphase (SEI) coating to protect the SPAN composite. Conversely, cyclic carbonates, ethylene carbonate (EC) and fluoroethylene carbonate (FEC) are shown to be very active. The EC reduction occurs only on a highly lithiated surface with a 3.0 Li/S molar ratio, the equivalent of the SPAN composite in an over-discharge regime with voltages close to 0.0 V vs. Li/Li + . The FEC reduction starts with a 2.0 Li/S molar ratio and above, suggesting that FEC could act as a useful additive in electrolyte formulations with EC. Both EC and FEC follow multiple reduction mechanisms to produce complex reduction products and LiF in the FEC case. We provide a mechanistic description for each detected decomposition path. The LiFSI salt also proves reactive against the lithiated SPAN surface. The FSI – defluorination is the dominant reduction path. However, the SO 2 NSO 2 F – and SO 2 NSO 2 2– species proved stable against S–N cleavage. This behavior makes the LiFSI salt a potential candidate for SPAN-based Li–S batteries because it produces LiF without releasing SO 2 .

25 ENERGY STORAGE↗
DOE OSTI2021

Anisotropic strength behavior of single-crystal TATB

High-rate strength behavior plays an important role in the shock initiation of high explosives, with plastic deformation serving to localize heat into hot spots and as a mechanochemical means to enhance reactivity. Recent simulations predict that detonation-like shocks produce highly reactive nanoscale shear bands in the layered crystalline explosive TATB (1,3,5-triamino-2,4,6-trinitrobenzene), but the thresholds leading to this response are poorly understood. We utilize molecular dynamics to simulate the high-rate compressive stress–strain response of TATB, with a focus on understanding flow behavior. Here, the dependence of strength on pressure and loading axis (crystal orientation) is explored. The deformation mechanisms fall broadly into two categories, with compression along crystal layers activating a buckling/twinning mode and compression normal to the layers producing nanoscale shear bands. Despite the complexity of the underlying mechanisms, the crystal exhibits relatively straightforward stress–strain curves. Most of the crystal orientations studied show rapid strain softening following the onset of yielding, which settles to a steady flow state. Trajectories are analyzed using five metrics for local states and structural order, but most of these metrics yield similar distributions for these deformation mechanisms. On the other hand, a recently proposed measure of intramolecular strain energy is found to most cleanly distinguish between these mechanisms, while also providing a plausible connection with mechanochemically accelerated decomposition kinetics. Localization of intramolecular strain energy is found to depend strongly on crystal orientation and pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Causal Inference Model Based on Random Forests to Identify the Effect of Soil Moisture on Precipitation

Soil moisture influences precipitation mainly through its impact on land–atmosphere interactions. Understanding and correctly modeling soil moisture–precipitation (SM–P) coupling is crucial for improving weather forecasting and subseasonal to seasonal climate predictions, especially when predicting the persistence and magnitude of drought. However, the sign and spatial structure of SM–P feedback are still being debated in the climate research community, mainly due to the difficulty in establishing causal relationships and the high degree of nonlinearity in land–atmosphere processes. To this end, we developed a causal inference model based on the Granger causality analysis and a nonlinear machine learning model. This model includes three steps: nonlinear anomaly decomposition, nonlinear Granger causality analysis, and evaluation of the quality of SM–P feedback, which eliminates the nonlinear response of interannual and seasonal variability and the memory effects of climatic factors and isolates the causal relationship of local SM–P feedback. We applied this model by using National Climate Assessment–Land Data Assimilation System (NCA-LDAS) datasets over the United States. Here, the results highlight the importance of nonlinear atmosphere responses in land–atmosphere interactions. In addition, the strong feedback over the southwestern United States and the Great Plains both highlight the impacts of topographic factors rather than only the sensitivity of evapotranspiration to soil moisture. Furthermore, the SM–P index defined by our framework is used to benchmark Earth system models (ESMs), which provides a new metric for efficiently identifying potential model biases in modeling local land–atmosphere interactions and may help the development of ESMs in improving simulations of water cycle variability and extremes.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1999

Automated Development of Accurate Algorithms and Efficient Codes for Computational Aeroacoustics

The simulation of sound generation and propagation in three space dimensions with realistic aircraft components is a very large time dependent computation with fine details. Simulations in open domains with embedded objects require accurate and robust algorithms for propagation, for artificial inflow and outflow boundaries, and for the definition of geometrically complex objects. The development, implementation, and validation of methods for solving these demanding problems is being done to support the NASA pillar goals for reducing aircraft noise levels. Our goal is to provide algorithms which are sufficiently accurate and efficient to produce usable results rapidly enough to allow design engineers to study the effects on sound levels of design changes in propulsion systems, and in the integration of propulsion systems with airframes. There is a lack of design tools for these purposes at this time. Our technical approach to this problem combines the development of new, algorithms with the use of Mathematica and Unix utilities to automate the algorithm development, code implementation, and validation. We use explicit methods to ensure effective implementation by domain decomposition for SPMD parallel computing. There are several orders of magnitude difference in the computational efficiencies of the algorithms which we have considered. We currently have new artificial inflow and outflow boundary conditions that are stable, accurate, and unobtrusive, with implementations that match the accuracy and efficiency of the propagation methods. The artificial numerical boundary treatments have been proven to have solutions which converge to the full open domain problems, so that the error from the boundary treatments can be driven as low as is required. The purpose of this paper is to briefly present a method for developing highly accurate algorithms for computational aeroacoustics, the use of computer automation in this process, and a brief survey of the algorithms that have resulted from this work. A review of computational aeroacoustics has recently been given by Lele.

Goodrich, John W.↗
DOE OSTI2021

Polarizable Water Potential Derived from a Model Electron Density

A new empirical potential for efficient, large scale molecular dynamics simulation of water is presented. The HIPPO (Hydrogen-like Intermolecular Polarizable POtential) force field is based upon the model electron density of a hydrogen-like atom. This framework is used to derive and parametrize individual terms describing charge penetration damped permanent electrostatics, damped polarization, charge transfer, anisotropic Pauli repulsion, and damped dispersion interactions. Initial parameter values were fit to Symmetry Adapted Perturbation Theory (SAPT) energy components for ten water dimer configurations, as well as the radial and angular dependence of the canonical dimer. The SAPT-based parameters were then systematically refined to extend the treatment to water bulk phases. The final HIPPO water model provides a balanced representation of a wide variety of properties of gas phase clusters, liquid water, and ice polymorphs, across a range of temperatures and pressures. This water potential yields a rationalization of water structure, dynamics, and thermodynamics explicitly correlated with an ab initio energy decomposition, while providing a level of accuracy comparable or superior to previous polarizable atomic multipole force fields. The HIPPO water model serves as a cornerstone around which similarly detailed physics-based models can be developed for additional molecular species.

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

Spinodal electronic phase separation during insulator-metal transitions

Electronic phase transitions such as insulator-metal transitions are common in strongly correlated systems. Here, using a combination of thermodynamic linear-stability analysis and phase-field simulations and employing VO 2 as a prototypical example, we predict that an insulator-metal transition driven by photoexcitation may involve an intermediate, modulated charge density state with a temperature-dependent characteristic wavelength. It is shown that such an intermediate two-phase electronic state is formed through a spinodal mechanism and that its formation can be generic for insulator-metal transitions driven by fast stimuli. This transient electronic phase separation is expected to stimulate future experimental and computational efforts.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Icarus Application to Dragonfly Heatshield

Icarus, an in-house material response solver developed at NASA Ames, is applied to perform 2-D bondline-temperature analysis on the heatshield surface of the latest Dragonfly entry system design for a mission to Titan. Two axisymmetric domains are used: the near-shoulder region and the entire heatshield. The heatshield consists of a PICA TPS bonded to a layered system of Aluminum honeycomb wrapped by an M55J carbon-fiber facesheet. The backshell, which is also included for more accurate analysis, consists of similar materials layered. Icarus simulations using both domains indicate that the in-plane thermal conductivity property of the M55J facesheet is observed to play a dominant role in bondline temperature. Similarly, the maximum bondline temperature is found around the PICA-tile-interface region instead of the near-shoulder region or the stagnation point based on the current trajectory using the orthotropic properties of M55J. Therefore, these findings indicate the significance of 2-D or higher-dimensional material modeling analysis in fully understanding how bondline temperature behaves across the heatshield and obtaining a basis for TPS design locations. Verification at different radial locations ahead of the shoulder, where thermal conduction is close to 1-D, between Icarus (2-D), Icarus (1-D), and FIAT is performed and good agreement is observed. With the entire heatshield domain, the total mass loss of the aeroshell materials is approximately 3% of their initial mass potentially due to thermal decomposition of phenolic resins inside PICA under the current high-heating environment.

TPS material modeling↗
NASA NTRS2002

The Delta x B = 0 Constraint Versus Minimization of Numerical Errors in MHD Simulations

The MHD equations are a system of non-strictly hyperbolic conservation laws. The non-convexity of the inviscid flux vector resulted in corresponding Jacobian matrices with undesirable properties. It has previously been shown by Powell et al. (1995) that an 'almost' equivalent MHD system in non-conservative form can be derived. This non-conservative system has a better conditioned eigensystem. Aside from Powell et al., the MHD equations can be derived from basic principles in either conservative or non-conservative form. The Delta x B = 0 constraint of the MHD equations is only an initial condition constraint, it is very different from the incompressible Navier-Stokes equations in which the divergence condition is needed to close the system (i.e., to have the same number of equations and the same number of unknown). In the MHD formulations, if Delta x B = 0 initially, all one needs is to construct appropriate numerical schemes that preserve this constraint at later time evolutions. In other words, one does not need the Delta x B condition to close the MHD system. We formulate our new scheme together with the Cargo & Gallice (1997) form of the MHD approximate Riemann solver in curvilinear grids for both versions of the MHD equations. A novel feature of our new method is that the well-conditioned eigen-decomposition of the non-conservative MHD equations is used to solve the conservative equations. This new feature of the method provides well-conditioned eigenvectors for the conservative formulation, so that correct wave speeds for discontinuities are assured. The justification for using the non-conservative eigen-decomposition to solve the conservative equations is that our scheme has a better control of the numerical error associated with the divergence of the magnetic condition. Consequently, computing both forms of the equations with the same eigen-decomposition is almost equivalent. It will be shown that this approach, using the non-conservative eigensystem when solving the conservative equations, also works well in the context of standard shock-capturing schemes.

Yee, H. C.↗
NASA NTRS1995

The coupling of fluids, dynamics, and controls on advanced architecture computers

This grant provided for the demonstration of coupled controls, body dynamics, and fluids computations in a workstation cluster environment; and an investigation of the impact of peer-peer communication on flow solver performance and robustness. The findings of these investigations were documented in the conference articles.The attached publication, 'Towards Distributed Fluids/Controls Simulations', documents the solution and scaling of the coupled Navier-Stokes, Euler rigid-body dynamics, and state feedback control equations for a two-dimensional canard-wing. The poor scaling shown was due to serialized grid connectivity computation and Ethernet bandwidth limits. The scaling of a peer-to-peer communication flow code on an IBM SP-2 was also shown. The scaling of the code on the switched fabric-linked nodes was good, with a 2.4 percent loss due to communication of intergrid boundary point information. The code performance on 30 worker nodes was 1.7 (mu)s/point/iteration, or a factor of three over a Cray C-90 head. The attached paper, 'Nonlinear Fluid Computations in a Distributed Environment', documents the effect of several computational rate enhancing methods on convergence. For the cases shown, the highest throughput was achieved using boundary updates at each step, with the manager process performing communication tasks only. Constrained domain decomposition of the implicit fluid equations did not degrade the convergence rate or final solution. The scaling of a coupled body/fluid dynamics problem on an Ethernet-linked cluster was also shown.

Atwood, Christopher↗
NASA NTRS2006

Aeroelastic Flight Data Analysis with the Hilbert-Huang Algorithm

This report investigates the utility of the Hilbert Huang transform for the analysis of aeroelastic flight data. It is well known that the classical Hilbert transform can be used for time-frequency analysis of functions or signals. Unfortunately, the Hilbert transform can only be effectively applied to an extremely small class of signals, namely those that are characterized by a single frequency component at any instant in time. The recently-developed Hilbert Huang algorithm addresses the limitations of the classical Hilbert transform through a process known as empirical mode decomposition. Using this approach, the data is filtered into a series of intrinsic mode functions, each of which admits a well-behaved Hilbert transform. In this manner, the Hilbert Huang algorithm affords time-frequency analysis of a large class of signals. This powerful tool has been applied in the analysis of scientific data, structural system identification, mechanical system fault detection, and even image processing. The purpose of this report is to demonstrate the potential applications of the Hilbert Huang algorithm for the analysis of aeroelastic systems, with improvements such as localized online processing. Applications for correlations between system input and output, and amongst output sensors, are discussed to characterize the time-varying amplitude and frequency correlations present in the various components of multiple data channels. Online stability analyses and modal identification are also presented. Examples are given using aeroelastic test data from the F-18 Active Aeroelastic Wing airplane, an Aerostructures Test Wing, and pitch plunge simulation.

Brenner, Martin J.↗