Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix 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,639 records · Page 91

A micro-scale cutting model for UD CFRP composites with thermo-mechanical coupling

Cutting a unidirectional carbon fiber-reinforced polymer (UD CFRP) structure is the basic unit for CFRP machining, which is a complex thermal-mechanically coupled process. To reveal the deformation mechanism and predict cutting force in UD CFRP micro cutting, a micro-scale fracture model for UD CFRP cutting with thermal-mechanical coupling is demonstrated in this paper, which captures the failure modes for fibers, matrix and the interface based on a micro-level RVE using a relatively simple damage based fracture method. The thermal-mechanical coupling model at the micro scale is developed on the basis of the plastic energy dissipation and frictional heating during cutting. Failure models for the fiber, matrix and interface region are applied depending on the material properties of each of these three phases. Numerical simulations based on the above model with different fiber orientations were performed to predict the deformation and forces of different components in UD CFRP. Cutting experiments with the same fiber orientations as considered in the simulations were carried out to validate the force and deformation results. The predicted force and deformation patterns match well with evidence from our experiments. In general, the cutting force is larger than the thrust force regardless of fiber orientation. The cutting force reaches a maximum as the fiber orientation approaches 90 , but thrust forces do not vary substantially across cases. When the fiber orientation is acute, the deformation of fibers is much smaller than when the cutting angle is obtuse. Surface roughness follows the same trend with cutting angle as fiber deformation.

42 ENGINEERING↗

A deep learning-based Bayesian framework for high-resolution calibration of building energy models

Calibrating building energy models (BEMs), i.e., closing discrepancy between modeling and field measurements, is of significance to support its applications in building sustainability and resilience analysis. However, as being widely used in practice, current Bayesian calibration is mostly performed in low-resolution (annual or monthly), instead of high-resolution (hourly or sub-hourly), which is crucial to support emerging BEM applications, such as building-renewable energy integration (demand response) and smart control. This is attributable to the gaps in current Bayesian calibration process, including (1) difficulty in supporting reliable high-resolution calibration with over-parameterization and multi-solution issues, (2) inadequacy of meta-model to capture temporal building dynamics in high-resolution, and (3) excessive computational burdens of covariance matrix calculation in Bayesian inference. Therefore, to close these gaps, this research proposes a novel deep learning-based Bayesian calibration framework, involving pre-calibration mechanism, Long Short-Term Memory as surrogate models, and simplified covariance matrix calculation, to calibrate BEMs in high temporal resolution (i.e., hourly) with enhanced accuracy and computational efficiency. Finally, the case study demonstrates its effectiveness to match modeling outcomes with measurements and realize CV-RMSE of < 30 % and NMBE of < 6 % in hourly resolution, as well as a significant reduction of calibration time (by > 99 %, from > 600 h to ~ 1.5 h).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Aerothermal/FEM Analysis of Hypersonic Sharp Leading Edges

Advanced hypersonic vehicles, like wave riders, will have sharp leading edges to minimize drag. These designs require accurate finite element modeling (FEM) of the thermal-structural behavior of a diboride ceramic matrix composite sharp leading edge. By coupling the FEM solver to an engineering model of the aerothermodynamic heating environment the impact of non catalytic surfaces, rarefied flow effects, and multidimensional conduction on the performance envelopes of sharp leading edges can be examined.

Kolodziej, Paul↗

Low-dimensional Representation of Error Covariance

Ensemble and reduced-rank approaches to prediction and assimilation rely on low-dimensional approximations of the estimation error covariances. Here stability properties of the forecast/analysis cycle for linear, time-independent systems are used to identify factors that cause the steady-state analysis error covariance to admit a low-dimensional representation. A useful measure of forecast/analysis cycle stability is the bound matrix, a function of the dynamics, observation operator and assimilation method. Upper and lower estimates for the steady-state analysis error covariance matrix eigenvalues are derived from the bound matrix. The estimates generalize to time-dependent systems. If much of the steady-state analysis error variance is due to a few dominant modes, the leading eigenvectors of the bound matrix approximate those of the steady-state analysis error covariance matrix. The analytical results are illustrated in two numerical examples where the Kalman filter is carried to steady state. The first example uses the dynamics of a generalized advection equation exhibiting nonmodal transient growth. Failure to observe growing modes leads to increased steady-state analysis error variances. Leading eigenvectors of the steady-state analysis error covariance matrix are well approximated by leading eigenvectors of the bound matrix. The second example uses the dynamics of a damped baroclinic wave model. The leading eigenvectors of a lowest-order approximation of the bound matrix are shown to approximate well the leading eigenvectors of the steady-state analysis error covariance matrix.

Tippett, Michael K.↗

Fatigue behavior of continuous fiber silicon-carbide-aluminum composites

Four lay-ups of continuous fiber silicon carbide (SCS2) fiber/aluminum matrix composites were tested to assess fatigue mechanisms including stiffness loss when cycled below their respective fatigue limits. The lay-ups were 0 (sub 8), 0(sub 2)/ + or - 45 (sub 2s), 0/90 (sub 2s),and 0/ + or 45/90 (subs). The data were compared with predictions from the author's previously published shakedown model which predicts fatigue-induced stiffness loss in metal matrix composites. A fifth lay-up, + or - 45 (sub 2s), was tested to compare shakedown and fatigue limits. The particular batch of silicon-carbide fibers tested in this program had a somewhat lower modulus (340 GPa) than expected and displayed poor bonding to the aluminum matrix. Good agreement was obtained between the stiffness loss model and the test data. The fatigue damage below the fatigue limit was primarily in the form of matrix cracking. The fatigue limit corresponded to the laminate shakedown for the + or - 45 (sub 2s) laminate.

Johnson, W. S.↗

Fatigue behavior of continuous-fiber silicon carbide/aluminum composites

Four lay-ups of continuous fiber silicon carbide (SCS2) fiber/aluminum matrix composites were tested to assess fatigue mechanisms including stiffness loss when cycled below their respective fatigue limits. The lay-ups were 0 (sub 8), 0(sub 2)/ + or - 45 (sub 2s), 0/90 (sub 2s), and 0/ + or - 45/90 (subs). The data were compared with predictions from the author's previously published shakedown model which predicts fatigue-induced stiffness loss in metal matrix composites. A fifth lay-up, + or - 45 (sub 2s), was tested to compare shakedown and fatigue limits. The particular batch of silicon-carbide fibers tested in this program had a somewhat lower modulus (340 GPa) than expected and displayed poor bonding to the aluminum matrix. Good agreement was obtained between the stiffness loss model and the test data. The fatigue damage below the fatigue limit was primarily in the form of matrix cracking. The fatigue limit corresponded to the laminate shakedown for the + or - 45 (sub 2s) laminate.

Johnson, W. S.↗

An Alternative Conceptual Model for the Spent Nuclear Fuel–Water Interaction in Deep Geologic Disposal Conditions

For the long-term safety assessment of direct disposal of spent nuclear fuel in deep geologic repositories, knowledge on the radionuclide release rate from the UO 2 matrix is essential. This work provides a conceptual model to explain the results of leaching experiments involving used nuclear fuel or simulant materials in confirmed reducing conditions. Key elements of this model are: direct effect of radiation from radiolytic species (including defects and excited states) in the solid and in the first water layers in contact with its surface; and excess H 2 may be produced due to processes occurring at the surface of the spent fuel and in confined water volumes, which may also play a role in keeping the spent fuel surface in a reduced state. The implication is that the fractional radionuclide release rate used in most long-term safety assessments (10 -7 year -1 ) is over estimated because it assumes that there is net UO 2 oxidation caused by radiolysis, in contrast with the alternative conceptual model presented here. Furthermore, conventional water radiolysis models and radiation chemical yields published in the literature are not directly applicable to a heterogeneous system such as the spent fuel–water interface. Suggestions are provided for future work to develop more reliable models for the long-term safety assessment of spent nuclear fuel disposal.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Analysis of local delaminations caused by angle ply matrix cracks

Two different families of graphite/epoxy laminates with similar layups but different stacking sequences, (0,theta,-theta) sub s and (-theta/theta/0) sub s were analyzed using three-dimensional finite element analysis for theta = 15 and 30 degrees. Delaminations were modeled in the -theta/theta interface, bounded by a matrix crack and the stress free edge. The total strain energy release rate, G, along the delamination front was computed using three different techniques: the virtual crack closure technique (VCCT), the equivalent domain Integral (EDI) technique, and a global energy balance technique. The opening fracture mode component of the strain energy release rate, Gl, along the delamination front was also computed for various delamination lengths using VCCT. The effect of residual thermal and moisture stresses on G was evaluated.

Salpekar, Satish A.↗

Weak Production of Strange and Charmed Ground-State Baryons in Nuclei

We present results for the quasielastic weak production of Λ and Σ hyperons induced by ν¯ scattering off nuclei in the kinematical region of interest for accelerator neutrino experiments. We employ realistic hole spectral functions and we describe the propagation of the hyperons in the nuclear medium by means of a Monte Carlo cascade. The latter strongly modifies the kinematics and the relative production rates of the hyperons, leading to a nonvanishing Σ+ cross section, to a sizable enhancement of the Λ production and to a drastic reduction of the Σ0 and Σ− distributions. We also compute the quasielastic weak Λc production cross section, paying special attention to estimate the uncertainties induced by the model dependence of the vacuum n→Λc weak matrix element. In this regard, the recent BESIII measurements of the branching ratios of Λc→Λl+νl (l=e,μ) are used to benchmark the available theoretical predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Experimental study on two-dimensional free-flying robot satellite model

The experimental study on a control method for a free flying space robotic arm was treated by means of a two-dimensional laboratory model. The main target is to develop a new control method for trajectory tracking or target capturing, considering dynamical interaction between the manipulator arm and the base vehicle in space micro-gravity environment. In order to simulate the micro-gravity environment mechanically, a laboratory model of a robot satellite supported on air bearings was developed. The model comprises a base equipped with power and air supplies and a two-link manipulator arm. This model has relatively low gravitational and frictional disturbance in planar motion. An on-line resolved motion rate control scheme with vision feedback is developed for experimenting capture operations. The scheme utilizes the Generalized Jacobian Matrix. In experiment, the acceleration environment of the model is evaluated firstly, then target capture operations are examined. The manipulator can properly chase and capture both a standing target and a moving target in spite of complex satellite/manipulator dynamical interactions. The experimental results confirm the validity of the Generalized Jacobian Matrix concept and the proposed control method.

Umetani, Yoji↗

Venus tectonic styles and crustal differentiation

Two of the most important constraints are known from Pioneer Venus data: the lack of a system of spreading rises, indicating distributed deformation rather than plate tectonics; and the high gravity/topography ratio, indicating the absence of an asthenosphere. In addition, the high depth/diameter ratios of craters on Venus indicate that Venus probably has no more crust than Earth. The problems of the character of tectonics and crustal formation and recycling are closely coupled. Venus appears to lack a recycling mechanism as effective as subduction, but may also have a low rate of crustal differentiation because of a mantle convection pattern that is more distributed, less concentrated, than Earth's. Distributed convection, coupled with the nonlinear dependence of volcanism on heat flow, would lead to much less magmatism, despite only moderately less heat flow, compared to Earth. The plausible reason for this difference in convective style is the absence of water in the upper mantle of Venus. We have applied finite element modeling to problems of the interaction of mantle convection and crust on Venus. The main emphasis has been on the tectonic evolution of Ishtar Terra, as the consequence of convergent mantle flow. The early stage evolution is primarily mechanical, with crust being piled up on the down-stream side. Then the downflow migrates away from the center. In the later stages, after more than 100 m.y., thermal effects develop due to the insulating influence of the thickened crust. An important feature of this modeling is the entrainment of some crustal material in downflows. An important general theme in both convergent and divergent flows is that of mixing vs. stratification. Models of multicomponent solid-state flow obtain that lower-density crustal material can be entrained and recycled, provided that the ration of low-density to high-density material is small enough (as in subducted slabs on Earth). The same considerations should apply in upflows; a small percent of partial melt may be carried along with its matrix and never escape to the surface. Models that assume melt automatically rising to the crust and no entrainment or other mechanism of recycling lower-density material obtain oscillatory behavior, because it takes a long time for heat to build up enough to overcome a Mg-rich low-density residuum. However, these models develop much thicker crust than consistent with estimates from crater depth/diameter ratios.

Kaula, W. M.↗

Fabrication and characterization of multifunctional Fe/TiO2 composite coatings

Highlights: • Electrodeposition of Fe/TiO{sub 2} composite coatings using methanesulfonate electrolytes. • TiO{sub 2} sol is synthesized by means of precipitation technique. • Electrodeposition kinetics obeys the improved Guglielmi's model. • Incorporation of TiO{sub 2} into Fe matrix results in an increase in microhardness. • Enhanced electrocatalytic and photocatalytic activity of Fe/TiO{sub 2} composites. - Abstract: Composite Fe/TiO{sub 2} coatings were electrodeposited using environmentally friendly colloidal methanesulfonate electrolytes. A sol of titania nano-particles was synthesized by means of precipitation technique. The electrodeposition parameters and the properties of composite coatings were compared with those typical of analogous composite electrodeposits produced from the plating bath containing suspended commercially available nano-powder Degussa P25. The electrodeposition process obeys the modified Guglielmi's model. The incorporation of TiO{sub 2} dispersed phase into iron matrix results in an appreciable increase in the microhardness due to dispersion strengthening. The microhardness of deposits obtained from electrolyte with the addition of nano-sized titania sol exceeds that which is observed for coatings produced from the bath containing aggregated Degussa P25 TiO{sub 2} particles. The introduction of dioxide titanium particles into Fe matrix yields photocatalytic properties towards the reaction of organic dye degradation and improves electrocatalytic performance towards hydrogen and oxygen evolution reactions.

36 MATERIALS SCIENCE↗

In-beam γ-ray spectroscopy of 37–42 P

The level schemes of the neutron-rich 37-42 P isotopes are investigated via in-beam gamma-ray spectroscopy following the fragmentation of a 45 Cl projectile beam at intermediate beam energies. Information on gamma-gamma coincidence relationships complemented by comparisons to shell-model calculations in the sd-pf model space were used to construct excitation level schemes for these neutron-rich nuclei. For the odd-mass 37,39 P isotopes, a level scheme is presented that appears essentially complete at low energies and exhausts the states predicted by the SDPF-MU shell-model Hamiltonian. In conclusion, simple Nilsson configurations are proposed for the low-lying excited states of 38,39,40,41 P from an analysis of the E2 transition matrix elements and moments calculated within the shell model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Optimization of micropillar sequences for fluid flow sculpting

Inertial fluid flow deformation around pillars in a microchannel is a new method for controlling fluid flow. Sequences of pillars have been shown to produce a rich phase space with a wide variety of flow transformations. Previous work has successfully demonstrated manual design of pillar sequences to achieve desired transformations of the flow cross section, with experimental validation. However, such a method is not ideal for seeking out complex sculpted shapes as the search space quickly becomes too large for efficient manual discovery. We explore fast, automated optimization methods to solve this problem. We formulate the inertial flow physics in microchannels with different micropillar configurations as a set of state transition matrix operations. These state transition matrices are constructed from experimentally validated streamtraces for a fixed channel length per pillar. This facilitates modeling the effect of a sequence of micropillars as nested matrix-matrix products, which have very efficient numerical implementations. With this new forward model, arbitrary micropillar sequences can be rapidly simulated with various inlet configurations, allowing optimization routines quick access to a large search space. We integrate this framework with the genetic algorithm and showcase its applicability by designing micropillar sequences for various useful transformations. We computationally discover micropillar sequences for complex transformations that are substantially shorter than manually designed sequences. We also determine sequences for novel transformations that were difficult to manually design. Finally, we experimentally validate these computational designs by fabricating devices and comparing predictions with the results from confocal microscopy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Learning physics-based reduced-order models from data using nonlinear manifolds

Here we present a novel method for learning reduced-order models of dynamical systems using nonlinear manifolds. First, we learn the manifold by identifying nonlinear structure in the data through a general representation learning problem. The proposed approach is driven by embeddings of low-order polynomial form. A projection onto the nonlinear manifold reveals the algebraic structure of the reduced-space system that governs the problem of interest. The matrix operators of the reduced-order model are then inferred from the data using operator inference. Numerical experiments on a number of nonlinear problems demonstrate the generalizability of the methodology and the increase in accuracy that can be obtained over reduced-order modeling methods that employ a linear subspace approximation.

97 MATHEMATICS AND COMPUTING↗

Uncertainty Models for the Hybrid Parametric Variation Method of Uncertainty Quantification; Analysis

There is some level of uncertainty in every finite element model (FEM), which flows to a level of uncertainty in predicted results. The purpose of uncertainty quantification (UQ) is to provide statistical bounds on prediction accuracy based on model uncertainty. This is distinct from model updating, which attempts to modify models to improve their accuracy. UQ does not improve the accuracy of models, but accepts that the models are inherently inaccurate and attempts to quantify the impact of that inaccuracy on predicted results. Previously, an alternate method for UQ, called the Hybrid Parametric Variation (HPV) method, was applied to Space Launch System (SLS) Hurty/Craig-Bampton (HCB) components to predict system-level statistics for launch vehicle attitude control transfer functions and core stage section loads due to buffet. The HPV method combines a parametric variation of the HCB fixed-interface (FI) modal frequencies with a nonparametric variation (NPV) method that randomly varies the HCB mass and stiffness matrices as Wishart random matrix distributions using random matrix theory (RMT). Alternatively, the most common method for modeling uncertainty in the structural dynamics community is a parametric approach, which varies physical parameters in the model. However, there are several disadvantages associated with the parametric method. Determining a reduced set of parameters that have a significant impact on the system response can be time consuming, and the selected parameter probability distributions are rarely reliably known. Therefore, in practice, the parameters are surrogates for the actual errors, and the link to parameter uncertainty is unknown. Another major drawback is that the uncertainty that can be represented is limited to the form of the nominal FEM. It is the experience of the authors that based on numerous aerospace programs, almost all FEM errors are in form rather than parameter values. This hypothesis is supported by the observation of the authors that it is almost never possible to ‘tune’ a FEM to match modal test results by only modifying model parameters. Model-form uncertainty cannot be directly represented by FEM input parameters nor included in a parametric approach. However, model-form uncertainty can be modeled using RMT, where a probability distribution is developed for the matrix ensemble of interest. The major advantage of the NPV method is that it covers errors in model form. The HPV method anchors uncertainty at the HCB component level to component modal test results by matching the HCB and test modes based on mode descriptions or other methods, and then applying differing levels of frequency variation. The specific variations depend on the confidence to which a component FEM has been validated through modal testing. The NPV method is layered on the frequency variation to match modal test self-orthogonality and cross-orthogonality (XO) results. Once the component uncertainty models are identified, they are assembled, and the uncertainty is propagated to the system level using a Monte Carlo (MC) analysis approach that generates statistics for system-level predictions This provides a UQ method that can be traced to test data, which can be updated as additional data and improved correlated models become available. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The purpose of this paper is to collect and present all of the theory for HPV that has been previously published in reports and papers and to present examples of its application. Specifically, component uncertainty models based on the dispersion of corresponding mass and stiffness matrices using proposed test/analysis correlation metrics are investigated. The first example is purely academic so that the true answers are known, and the validity of the HPV method and the corresponding uncertainty models can be determined. The second example is an application to a component that is design specific to the SLS. Based on this work and other assessments, the HPV method provides another tool to the toolset used for complex system UQ analysis. From experience gathered to date using the HPV method, additional design specific applications must be investigated to provide further confidence in the validity of the HPV method of UQ analysis.

Uncertainty quantification↗

Mechanics of interfaces in fiber reinforced SiC/RBSN ceramic matrix composites

The mechanical behavior of continuous fiber reinforced SiC/RBSN composites with strong and weak interface characteristics is evaluated. Both catastrophic and noncatastrophic failures are observed in tensile specimens. Effects of fiber/matrix interface debonding (splitting) parallel to the fibers are discussed. Micromechanical models incorporating residual stresses to calculate the critical matrix cracking strength, ultimate strength and work of pull-out are reviewed and used to predict composite response. Experimental results are compared to analytical predictions.

Chulya, Abhisak↗

Design of multivariable feedback control systems via spectral assignment using reduced-order models and reduced-order observers

The feasibility of using reduced order models and reduced order observers with eigenvalue/eigenvector assignment procedures is investigated. A review of spectral assignment synthesis procedures is presented. Then, a reduced order model which retains essential system characteristics is formulated. A constant state feedback matrix which assigns desired closed loop eigenvalues and approximates specified closed loop eigenvectors is calculated for the reduced order model. It is shown that the eigenvalue and eigenvector assignments made in the reduced order system are retained when the feedback matrix is implemented about the full order system. In addition, those modes and associated eigenvectors which are not included in the reduced order model remain unchanged in the closed loop full order system. The full state feedback design is then implemented by using a reduced order observer. It is shown that the eigenvalue and eigenvector assignments of the closed loop full order system rmain unchanged when a reduced order observer is used. The design procedure is illustrated by an actual design problem.

Mielke, R. R.↗