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,675 records · Page 93

Mechanics of composite materials: Recent advances; Proceedings of the Symposium, Virginia Polytechnic Institute and State University, Blacksburg, VA, August 16-19, 1982

The present conference on the mechanics of composites discusses microstructure's influence on particulate and short fiber composites' thermoelastic and transport properties, the elastoplastic deformation of composites, constitutive equations for viscoplastic composites, the plasticity and fatigue of metal matrix composites, laminate damping mechanisms, the micromechanical modeling of Kevlar/epoxy composites' time-dependent failure, the variational characterization of waves in composites, and computational methods for eigenvalue problems in composite design. Also discussed are the elastic response of laminates, elastic coupling nonlinear effects in unsymmetrical laminates, elasticity solutions for laminate problems having stress singularities, the mechanics of bimodular composite structures, the optimization of laminated plates and shells, NDE for laminates, the role of matrix cracking in the continuum constitutive behavior of a damaged composite ply, and the energy release rates of various microcracks in short fiber composites.

Hashin, Z.↗

An Automated Meshing Framework for Progressive Damage Analysis of Fabrics Using CompDam

A meshing and model construction framework for progressive damage analysis of fabric composites is proposed. A mesoscale idealization is used such that the fiber tows are considered homogeneous and transversely isotropic continua, and the matrix is a homogeneous isotropic continuum. By using this idealization, relatively mature damage modeling techniques developed for laminates of unidirectional plies are transferrable to fabrics. Cohesive interface elements and a continuum damage mechanics (CDM) model implemented in the software CompDam are used to account for the key damage modes including cracking at the tow/matrix interface, transverse matrix cracking in the tows, and fiber breakage in the tows. Matrix degradation is accounted for by plasticity. The modeling approach is demonstrated for a 3D orthogonal woven fabric. The main novelty of the proposed model construction framework is in its flexibility resulting from a formulation that is agnostic to the fabric architecture and implementation using mesh based operations. The algorithm does not make any assumptions about the tow sizes, shape, path, or proximity to other tows, thus it is applicable to a broad range of fabric architectures. The matrix region mesh surrounding the tows is constructed from domain boundary and tow surface meshes (i.e. through a series of mesh-based operations) such that the model domain can be arbitrarily shaped.

Andrew Cole Bergan↗

Factorized three-body S-matrix restrained by the Yang–Baxter equation and quantum entanglements

This paper investigates the physical effects of the Yang–Baxter equation (YBE) to quantum entanglements through the 3-body S-matrix in entangling parameter space. The explicit form of 3-body S-matrix Ř{sub 123}(θ,φ) based on the 2-body S-matrices is given due to the factorization condition of YBE. The corresponding chain Hamiltonian has been obtained and diagonalized, also the Berry phase for 3-body system is given. It turns out that by choosing different spectral parameters the Ř(θ,φ)-matrix gives GHZ and W states respectively. The extended 1-D Kitaev toy model has been derived. Examples of the role of the model in entanglement transfer are discussed. - Highlights: • We give the relation between 3-body S-matrix and 3-qubit entanglement. • The relation between 3-qubit and 2-qubit entanglements is investigated via YBE. • 1D Kitaev toy model is derived by the Type-II solution of YBE. • The condition of YBE kills the “Zero boundary mode” in our chain model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Neutron Number Probability Distributions in a Subcritical System Using the Forward Master Equation

The subject of this paper is the determination of the equilibrium probability distribution function (PDF) of the neutron population within a subcritical multiplying system in the presence of an intrinsic random source as a result of spontaneous fission (SF). Such systems are strongly stochastic and low order moments such as the mean and variance of the neutron number are not sufficient to characterize the true state of the neutron population. We base our work on the forward Master equation in a lumped model formulation for P{sub n} (t), n = 0, 1, 2,..., the probability of there existing η neutrons at time t and obtain solutions for the steady state distributions. The key parameters in this model, neutron lifetime and fission probability, are rigorously obtained by numerically solving the k-eigenvalue problem in a one-dimensional spherical system and using standard definitions of lumped parameters in terms of production and loss rates in finite systems. Previous work has examined single-chain PDFs (a single fission and the resulting progeny) for neutrons, uncoupled neutron-photon single-chain PDFs, both using a lumped backward Master equation formulation, as well as neutron PDFs in supercritical systems. In this paper, we analyze strictly subcritical systems, for which we compare analytically derived solutions with numerical computations for different multiplicity distribution models. For subcritical systems where the neutron population remains finite at steady state, we show that the forward Master equation can be very efficiently solved numerically to obtain the equilibrium distributions. Finally, this work is motivated by potential application to passive interrogation for detection of diverted nuclear materials as well to criticality safety. The equilibrium neutron number probability distribution in subcritical systems has been numerically investigated for several models of the induced and spontaneous fission multiplicity distribution. A well known approximate but closed form solution, namely the quadratic model, was shown to be unacceptably accurate for subcritical media. The matrix inversion and the recursion methods agree well and were benchmarked with the binary fission model, but the recursion formulation proves far more computationally efficient. A qualitative change in the shape of the neutron PDF is observed for increasing source strength, as expected, and a transition from a stochastic to a deterministic system is expected for a sufficiently strong source. We are currently extending this approach to obtain the equilibrium joint distribution of neutrons within the system and those that have leaked, to explicitly represent the neutron PDF as observed by an external detector. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Evaluation and prediction of long-term environmental effects of nonmetallic materials

The properties of a number of nonmetallic materials were evaluated experimentally in simulated space environments in order to develop models for accelerated test methods useful for predicting such behavioral changes. Graphite-epoxy composites were exposed to thermal cycling. Adhesive foam tapes were subjected to a vacuum environment. Metal-matrix composites were tested for baseline data. Predictive modeling designed to include strength and aging effects on composites, polymeric films, and metals under such space conditions (including the atomic oxygen environment) is discussed. The Korel 8031-00 high strength adhesive foam tape was shown to be superior to the other two tested.

Papazian, H.↗

Diffraction Analysis of Mesh Deployable Reflector Antennas

A formulation and many representative numerical results for mesh reflector antennas are presented. The reflection coefficient matrix for the prescribed mesh configuration was determined and the local coordinate system of the mesh cells at each point on the curved reflector surface was accentuated. A novel strip aperture model was used to formulate the transmission coefficient matrix for a variety of mesh cell configurations. Numerical data are tailored to the dimensions of a conceptually designed land mobile satellite system (LMSS) which employs a large mesh deployable offset parabolic antenna. Results are shown for an offset parabolic reflector with mesh surfaces similar to the mesh surface of tracking and data relay satellite system (TDRSS).

Rahmat-Samii, Y.↗

On the formulation of a minimal uncertainty model for robust control with structured uncertainty

In the design and analysis of robust control systems for uncertain plants, representing the system transfer matrix in the form of what has come to be termed an M-delta model has become widely accepted and applied in the robust control literature. The M represents a transfer function matrix M(s) of the nominal closed loop system, and the delta represents an uncertainty matrix acting on M(s). The nominal closed loop system M(s) results from closing the feedback control system, K(s), around a nominal plant interconnection structure P(s). The uncertainty can arise from various sources, such as structured uncertainty from parameter variations or multiple unsaturated uncertainties from unmodeled dynamics and other neglected phenomena. In general, delta is a block diagonal matrix, but for real parameter variations delta is a diagonal matrix of real elements. Conceptually, the M-delta structure can always be formed for any linear interconnection of inputs, outputs, transfer functions, parameter variations, and perturbations. However, very little of the currently available literature addresses computational methods for obtaining this structure, and none of this literature addresses a general methodology for obtaining a minimal M-delta model for a wide class of uncertainty, where the term minimal refers to the dimension of the delta matrix. Since having a minimally dimensioned delta matrix would improve the efficiency of structured singular value (or multivariable stability margin) computations, a method of obtaining a minimal M-delta would be useful. Hence, a method of obtaining the interconnection system P(s) is required. A generalized procedure for obtaining a minimal P-delta structure for systems with real parameter variations is presented. Using this model, the minimal M-delta model can then be easily obtained by closing the feedback loop. The procedure involves representing the system in a cascade-form state-space realization, determining the minimal uncertainty matrix, delta, and constructing the state-space representation of P(s). Three examples are presented to illustrate the procedure.

Belcastro, Christine M.↗

Emergent Chiral Spin Liquid: Fractional Quantum Hall Effect in a Kagome Heisenberg Model

The fractional quantum Hall effect (FQHE) realized in two-dimensional electron systems under a magnetic field is one of the most remarkable discoveries in condensed matter physics. Interestingly, it has been proposed that FQHE can also emerge in time-reversal invariant spin systems, known as the chiral spin liquid (CSL) characterized by the topological order and the emerging of the fractionalized quasiparticles. A CSL can naturally lead to the exotic superconductivity originating from the condense of anyonic quasiparticles. Although CSL was highly sought after for more than twenty years, it had never been found in a spin isotropic Heisenberg model or related materials. By developing a density-matrix renormalization group based method for adiabatically inserting flux, we discover a FQHE in a spin-$ \frac{1}{2}$ isotropic kagome Heisenberg model. We identify this FQHE state as the long-sought CSL with a uniform chiral order spontaneously breaking time reversal symmetry, which is uniquely characterized by the half-integer quantized topological Chern number protected by a robust excitation gap. The CSL is found to be at the neighbor of the previously identified Z2 spin liquid, which may lead to an exotic quantum phase transition between two gapped topological spin liquids.

MAGNETIC PROPERTIES AND MATERIALS↗

Effects and detectability of quasi-single field inflation in the large-scale structure and cosmic microwave background

Quasi-single field inflation predicts a peculiar momentum dependence in the squeezed limit of the primordial bispectrum which smoothly interpolates between the local and equilateral models. This dependence is directly related to the mass of the isocurvatons in the theory which is determined by the supersymmetry. Therefore, in the event of detection of a non-zero primordial bispectrum, additional constraints on the parameter controlling the momentum-dependence in the squeezed limit becomes an important question. We explore the effects of these non-Gaussian initial conditions on large-scale structure and the cosmic microwave background, with particular attention to the galaxy power spectrum at large scales and scale-dependence corrections to galaxy bias. We determine the simultaneous constraints on the two parameters describing the QSF bispectrum that we can expect from upcoming large-scale structure and cosmic microwave background observations. We find that for relatively large values of the non-Gaussian amplitude parameters, but still well within current uncertainties, galaxy power spectrum measurements will be able to distinguish the QSF scenario from the predictions of the local model. A CMB likelihood analysis, as well as Fisher matrix analysis, shows that there is also a range of parameter values for which Planck data may be able distinguish between QSF models and the related local and equilateral shapes. Given the different observational weightings of the CMB and LSS results, degeneracies can be significantly reduced in a joint analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

The effect of interface properties on nickel base alloy composites

This program was performed to assess the extent to which mechanical behavior models can predict the properties of sapphire fiber/nickel aluminide matrix composites and help guide their development by defining improved combinations of matrix and interface coating. The program consisted of four tasks: 1) selection of the matrices and interface coating constituents using a modeling-based approach; 2) fabrication of the selected materials; 3) testing and evaluation of the materials; and 4) evaluation of the behavior models to develop recommendations. Ni-50Al and Ni-20AI-30Fe (a/o) matrices were selected which gave brittle and ductile behavior, respectively, and an interface coating of PVD YSZ was selected which provided strong bonding to the sapphire fiber. Significant fiber damage and strength loss was observed in the composites which made straightforward comparison of properties with models difficult. Nevertheless, the models selected generally provided property predictions which agreed well with results when fiber degradation was incorporated. The presence of a strong interface bond was felt to be detrimental in the NiAI MMC system where low toughness and low strength were observed.

Groves, M.↗

Multi-scale computational analysis of unidirectional carbon fiber reinforced polymer composites under various loading conditions

In this study, a multi-scale computational analysis based on representative volume element (RVE) modeling and molecular dynamics (MD) simulations is developed to investigate the microscopic failure mechanisms of unidirectional (UD) carbon fiber reinforced polymer (CFRP) composites. The average properties of the 200-nm thickness interphase region between fiber and matrix are characterized through MD simulations and an analytical gradient model. The results demonstrate that the interphase region has higher Young’s modulus and strength, compared to the bulk matrix. This stiffened interphase region influences the composite response significantly. Specifically, the traditional two-phase model with zero-thickness interface fails to capture the stress-strain behavior compared to the experimental data. However, by adding the interphase region to a modified RVE model, the accuracy of simulation results will be improved significantly. Furthermore, a coupled experimental-computational micromechanics approach is adopted to calibrate and validate the cohesive parameters of the interface. By including the cohesive interface, our modified RVE model accurately captures the failure strength of the composites. Finally, different failure mechanisms for specimens are investigated using our multi-scale computational framework. In conclusion, the results show that the failure modes of UD CFRP composites are very complex and multiple failure mechanisms co-exist depending on the loading conditions, agreeing well with our experimental analyses.

36 MATERIALS SCIENCE↗

Inhibition of soluble epoxide hydrolase attenuates hepatic fibrosis and endoplasmic reticulum stress induced by carbon tetrachloride in mice

Liver fibrosis is a pathological condition in which chronic inflammation and changes to the extracellular matrix lead to alterations in hepatic tissue architecture and functional degradation of the liver. Inhibitors of the enzyme soluble epoxide hydrolase (sEH) reduce fibrosis in the heart, pancreas and kidney in several disease models. In this study, we assess the effect of sEH inhibition on the development of fibrosis in a carbon tetrachloride (CCl{sub 4})-induced mouse model by monitoring changes in the inflammatory response, matrix remolding and endoplasmic reticulum stress. The sEH inhibitor 1-trifluoromethoxyphenyl-3-(1-propionylpiperidin-4-yl) urea (TPPU) was administered in drinking water. Collagen deposition in the liver was increased five-fold in the CCl{sub 4}-treated group, and this was returned to control levels by TPPU treatment. Hepatic expression of Col1a2 and 3a1 mRNA was increased over fifteen-fold in the CCl{sub 4}-treated group relative to the Control group, and this increase was reduced by 50% by TPPU treatment. Endoplasmic reticulum (ER) stress observed in the livers of CCl{sub 4}-treated animals was attenuated by TPPU treatment. In order to support the hypothesis that TPPU is acting to reduce the hepatic fibrosis and ER stress through its action as a sEH inhibitor we used a second sEH inhibitor, trans-4-(4-[3-(4-trifluoromethoxy-phenyl)-ureido]-cyclohexyloxy)-benzoic acid (t-TUCB), and sEH null mice. Taken together, these data indicate that the sEH may play an important role in the development of hepatic fibrosis induced by CCl{sub 4}, presumably by reducing endogenous fatty acid epoxide chemical mediators acting to reduce ER stress. - Highlights: • We administer an inhibitor of sEH in a CCl4 murine model. • sEH inhibition reduces liver collagen deposition and pro-fibrotic gene expression. • sEH inhibition induces MMP-1a activity.

60 APPLIED LIFE SCIENCES↗

3D Modeling Effect of Spherical Inclusions on the Magnetostriction of Bulk Superconductors

In this paper, the dependence of the effective magnetostriction of bulk superconductors on the elastic parameters including the volume fraction and elastic modulus ratio is studied by a three-dimensional model consisting of a spherical inclusion–superconducting matrix system. The effect of the elastic modulus and volume fraction on the magnetostriction is also obtained through the magnetostriction loop. The results indicate that the elastic modulus and volume fraction have obvious effects on the effective magnetostriction of the superconducting composite, which gives an explanation about the differences between the experimental and the theoretical results. Furthermore, it is worth pointing out that the linear field dependence of magnetostriction is unique to the Bean model by comparing the curve shapes of the magnetostriction loop with and without inclusion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nonclassical effects in a highly nonlinear generalized homogeneous Dicke model

An intensity dependent nonlinear coupling model of N two-level atoms (generalized Dicke model) interacting dispersively with a bimodal cavity field via two-photon transitions is investigated in a scenario where the rotating wave approximation is assumed. The model becomes homogeneous in the sense that the spin transition frequency is the same for all atoms and the coupling constants emerging from the collective interactions of the atomic system with the cavity field depend only on the particular radiation field mode. This allows us to represent the Dicke Hamiltonian entirely in terms of the total angular momentum J. It is assumed that, initially, the atomic system and the field are in a disentangled state where the field modes are in Glauber coherent states and the atomic system is a superposition of states |JM> (Dicke states). The model is numerically tested against simulations of normal squeezing variance of the field, squeezing factors based on the Heisenberg uncertainty principle, along with the statistical properties of the light leading to the possible production of nonclassical effects, such as degree of second-order coherence in the modes, degree of intermode correlation, as well as violation of the Cauchy-Schwartz inequality. Analytical expression of the total density operator matrix elements at t>0 shows the present nonlinear model to be strongly entangled, which is reflected in the time evolution of the linear entropy, where the superposition states are reduced to statistical mixtures. Thus, the present generalized Dicke model does not preserve the modulus of the Bloch vector. The computations, performed in the weak coupling and strong field limits, were conducted via second-order Dyson perturbative expansion of the time evolution operator matrix elements for the totality of the angular momentum states of the atomic system. - Highlights: > The present highly nonlinear generalized Dicke model involves two-photon transitions. > Quantum coherences significantly affect the dynamics of the atomic system. > At t>0 the time-dependent density operator matrix becomes strongly entangled. > Linear entropy shows that the superposition states turn into statistical mixtures. > The present Dicke model does not preserve the modulus of the Bloch vector.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effects of error covariance structure on estimation of model averaging weights and predictive performance

When conducting model averaging for assessing groundwater conceptual model uncertainty, the averaging weights are often evaluated using model selection criteria such as AIC, AICc, BIC, and KIC (Akaike Information Criterion, Corrected Akaike Information Criterion, Bayesian Information Criterion, and Kashyap Information Criterion, respectively). However, this method often leads to an unrealistic situation in which the best model receives overwhelmingly large averaging weight (close to 100%), which cannot be justified by available data and knowledge. It was found in this study that this problem was caused by using the covariance matrix, CE, of measurement errors for estimating the negative log likelihood function common to all the model selection criteria. This problem can be resolved by using the covariance matrix, Cek, of total errors (including model errors and measurement errors) to account for the correlation between the total errors. An iterative two-stage method was developed in the context of maximum likelihood inverse modeling to iteratively infer the unknown Cek from the residuals during model calibration. The inferred Cek was then used in the evaluation of model selection criteria and model averaging weights. While this method was limited to serial data using time series techniques in this study, it can be extended to spatial data using geostatistical techniques. The method was first evaluated in a synthetic study and then applied to an experimental study, in which alternative surface complexation models were developed to simulate column experiments of uranium reactive transport. It was found that the total errors of the alternative models were temporally correlated due to the model errors. The iterative two-stage method using Cek resolved the problem that the best model receives 100% model averaging weight, and the resulting model averaging weights were supported by the calibration results and physical understanding of the alternative models. Using Cek obtained from the iterative two-stage method also improved predictive performance of the individual models and model averaging in both synthetic and experimental studies.

error correlation↗

Formation of silver nanoparticles inside a soda-lime glass matrix in the presence of a high intensity Ar{sup +} laser beam

Formation and motion of the silver nanoparticles inside an ion-exchanged soda-lime glass in the presence of a focused high intensity continuous wave Ar{sup +} laser beam (intensity: 9.2 x 10{sup 4} W/cm{sup 2}) have been studied in here. One-dimensional diffusion equation has been used to model the diffusion of the silver ions into the glass matrix, and a two-dimensional reverse diffusion model has been introduced to explain the motion of the silver clusters and their migration toward the glass surface in the presence of the laser beam. The results of the mentioned models were in agreement with our measurements on thickness of the ion-exchange layer by means of optical microscopy and recorded morphology of the glass surface around the laser beam axis by using a Mirau interferometer. SEM micrographs were used to extract the size distribution of the migrated silver particles over the glass surface.

36 MATERIALS SCIENCE↗

Recent Developments in SHERPA

Some recent QCD-related developments in the SHERPA event generator are presented. In the past decades, event generators such as PYTHIA [1, 2] and HERWIG [3, 4] have been central for nearly all physics analyses at particle physics experiments at the high-energy frontier. This will also hold true at the LHC, where a large number of interesting signals for new particles or new phenomena (the Higgs boson or any other manifestation of the mechanism behind electro-weak symmetry breaking, supersymmetry, extra dimensions etc.) is hampered by a plethora of severe, sometimes overwhelming backgrounds. Nearly all of them are largely influenced by QCD. Therefore it seems fair to say that the success of the LHC in finding new physics may very well depend on a deep and detailed understanding of old physics, like QCD. Examples for this include, among others, the central-jet veto for the vector boson fusion channel for Higgs production or topologies, where gauge bosons emerge in association with many jets, a background for many search channels. In a reflection on increased needs by the experimental community, aiming at higher precision, incorporation of new physics models and so on, the work horses of old have undergone serious renovation efforts, resulting in new, improved versions of the respective codes, namely PYTHIA8 [5] and HERWIG++ [6]. In addition a completely new code, SHERPA [7], has been constructed and is in the process of maturing. The status of this code is the topic of this contribution. SHERPA's hallmark property is the inclusion of higher-order tree-level QCD contributions, leading to an improved modelling of jet production. They are introduced through a full-fledged matrix element generator, AMEGIC++ [8], which is capable of generating matrix elements and corresponding phase space mappings for processes with multi-particle final states in various models, including the Standard Model, anomalous gauge triple and quadruple couplings according to [9, 10], the Minimal Supersymmetric Standard Model with Feynman rules from [11], the ADD-model of extra dimensions [12, 13], and a model with an extra U(1) singlet coupling to the Higgs boson only [14]. The code has been thoroughly tested and validated [15]. This code, however, is limited, especially in the treatment of many ({ge} 6) external QCD particles. Therefore, in the near future, SHERPA will incorporate another, new matrix element generator, COMIX, which is based on Berends-Giele recursion relations [16] and color-dressing [17] rather than color-ordering. In Tabs. 1 and 2 some example cross sections for gg {yields} ng at fixed energies and pp {yields} b{bar b} + n jets obtained with this program are exhibited and compared to those from other programs. In addition, concerning the calculation of higher-order matrix elements and cross sections, there have been first steps towards an automation of such calculations at truly next-to leading order accuracy. They manifest themselves in the implementation of a procedure [19] to fully automatically construct and evaluate Catani-Seymour dipole subtraction terms [20] for the real part of such NLO calculations. The results from the matrix element calculations are merged with the subsequent parton shower through the formalism of [21, 22]. The results of its implementation in SHERPA [23] has recently been compared with other algorithms [24]. Although there remains some dispute about the theoretical equivalence of the different approaches, the overall results show satisfying agreement with each other, such that they can be used with confidence for data analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗