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 937 records · Page 52

Optical Design and Validation of an Infrared Transmissive Spectrum Splitting Concentrator Photovoltaic Module

A new modular, hybrid solar power system is designed to generate both electrical and thermal energy by utilizing the full solar spectrum. The key element, an infrared-transparent concentrator photovoltaic (CPV) module, acts as a spectrum splitter, dividing solar radiation into two parts. The ultraviolet and visible light (“in-band”) are converted to electricity with high efficiency in CPV cells, while the infrared light (“out-of-band”) is transmitted directly to a thermal receiver, where thermal power may be converted to electricity by a suitable heat engine or used directly for industrial process heat applications whenever needed. Here, we describe the optical design, modeling, fabrication, and performance validation of this novel spectrum splitting CPV module. A transfer matrix style approach, cumulative transmission model, is built to study the reflection, absorption, and transmission in each layer of the CPV module. To optimize the optical performance, different materials for module superstrate/substrate, encapsulant, cell substrate, and cooling fluids are compared in order to enhance the transmission of out-of-band light through the CPV module by minimizing absorption. Six antireflection coatings along with front and backside electrical contact grids are designed to maximize transmittance of in-band light to the cell and out-of-band light to the thermal receiver. Here, the final design, currently being prototyped, predicts out-of-band light transmission to the thermal receiver of 74.1% (for the passively cooled version) and 65.3% (for the actively cooled version). When epitaxial liftoff technology is applied, the transmission will change to 80.8% (passively cooled) and 71.9% (actively cooled). Experimental prototypes show good agreement with modeled optical performance.

14 SOLAR ENERGY↗

Effects of fiber and interfacial layer architectures on the thermoplastic response of metal matrix composites

Examined here is the effect of fiber and interfacial layer morphologies on thermal fields in metal matrix composites (MMCs). A micromechanics model based on an arbitrarily layered concentric cylinder configuration is used to calculate thermal stress fields in MMCs subjected to spatially uniform temperature changes. The fiber is modelled as a layered material with isotropic or orthotropic elastic layers, whereas the surrounding matrix, including interfacial layers, is treated as a strain-hardening, elastoplastic, von Mises solid with temperature-dependent parameters. The solution to the boundary-value problem of an arbitrarily layered concentric cylinder under the prescribed thermal loading is obtained using the local/global stiffness matrix formulation originally developed for stress analysis of multilayered elastic media. Examples are provided that illustrate how the morphology of the SCS6 silicon carbide fiber and the use of multiple compliant layers at the fiber/matrix interface affect the evolution of residual stresses in SiC/Ti composites during fabrication cool-down.

Pindera, Marek-Jerzy↗

Near-Threshold Rovibrational Excitation of HF by Electron Impact

The pronounced threshold peaks found in cross sections for electron-impact excitation of the vibrational levels of HF have been studied using an R-matrix method, in which electron exchange, correlation effects and polarization are included. Vibration is modelled inside the R-matrix sphere using the energy-modified adiabatic phase matrix method (EMAP). Rotation of the target molecule is treated outside the R-matrix sphere, using the frame transformation method. The external rovibronic close-coupling equations are solved by analytic methods. Analysis of the multichannel threshold structures indicates that the peaks found for rovibrational excitation of nu = 1 can be interpreted as virtual states while those for nu = 2 correspond to nuclear-excited Feshbach resonances. These two different results correspond simply to a change of sign of a small parameter that is determined by the threshold analysis. Whereas the overall shapes of the excitation peaks agree with earlier theoretical results, we have found detailed rotational fine structures for scattering energies on the order of the rotational constant above and below the opening of the vibrational channels. The partial and total cross sections for rovibrational state-to-state transitions in the energy region of the peaks are analysed. Characteristic cusp structures and excitation peaks are found at the rovibrational thresholds for transitions with small angular momentum transfer (Delta)j =0, 1 and 2. Differential cross sections are similar to recent measurements. These experimental data do not show the relatively sharp rise at the vibrational thresholds indicated by theory. Theory and experiment are in good agreement for the shape and magnitude of the angular dependence of cross sections for rovibrational state-to-state transitions, especially at larger energies.

Thummel, Helmar T.↗

Ultralight Door Design, Development and Demonstration Project (Final Scientific Report)

The project objectives associated with the Ultralight Door Project were to develop a lightweight driver's side door, 42.5% lighter than a current model year high-volume baseline door, at an incremental cost not to exceed $ \$ $ 5 per pound saved, while maintaining equivalent functionality, safety, corrosion, durability and sound transmission performance of the baseline vehicle. The project approach included the selection of a door architecture, development of a key performance target matrix and incremental cost model. Several unique door structures were developed and evaluated using computer-aided engineering (CAE) and cost modeling tools. The aluminum-intensive door structure was the only design capable of meeting both the performance and cost criteria objectives. The door module, glass, interior trim, latch and electrical subsystems were functionally integrated with the aluminum door structure, enabling additional mass reduction beyond material substitution. The lightweight door design was prototyped, integrated with baseline vehicles and tested to validate equivalent functionality and performance. The project results demonstrate a mass reduction of 15.2 kg (40%) at an incremental cost of $ \$ $2.81 per pound saved, while maintaining functionality and performance of the baseline door. The environmental benefit associated with the Ultralight Door enables a 6g CO 2 eq. per kilometer environmental benefit, for a 4-door vehicle powered by a downsized fuel injected gasoline combustion engine. The project results were presented and published in the 2018 SAE International Conference proceedings and in the Journal of Life Cycle Assessment.

36 MATERIALS SCIENCE↗

Comparison of Coarse-Grained Approaches in Predicting Polymer Nanocomposite Phase Behavior

Because of the considerable parameter space, efficient theoretical and simulation methods are required to predict the morphology and guide experiments in polymer nanocomposites (PNCs). Unfortunately, theoretical and simulation methods are restricted in their ability to accurately map to experiments based on necessary approximations and numerical limitations. In this study, we provide direct comparisons of two recently developed coarse-grained approaches for modeling polymer nanocomposites (PNCs): polymer nanocomposite field theory (PNC-FT) and dynamic mean-field theory (DMFT). These methods are uniquely suited to efficiently capture mesoscale phase behavior of PNCs in comparison to other theoretical and simulation frameworks. We demonstrate the ability of both methods to capture macrophase separation and describe the thermodynamics of PNCs. We systematically test how the nanoparticle morphology in PNCs is affected by a uniform probability distribution of grafting sites, common in field-based methods, versus random discrete grafting sites on the nanoparticle surface. We also analyze the accuracy of the mean-field approximation in capturing the phase behavior of PNCs. Moreover, the DMFT method introduces the ability to describe nonequilibrium phase behavior while the PNC-FT method is strictly an equilibrium method. With the DMFT method we are able to show the evolution of nonequilibrium states toward their equilibrium state and a qualitative assessment of the dynamics in these systems. These simulations are compared to experiments consisting of polystyrene grafted gold nanorods in a poly(methyl methacrylate) matrix to ensure the model gives results that qualitatively agree with the experiments. This study reveals that nanoparticles in a relatively high matrix molecular weight are trapped in a nonequilibrium state and demonstrates the utility of the DMFT framework in capturing nonequilibrium phase behavior of PNCs. In conclusion, both the PNC-FT and DMFT framework are promising methods to describe the thermodynamic and nonequilibrium phase behavior of PNCs.

36 MATERIALS SCIENCE↗

Hanle-Zeeman Scattering Matrix for Magnetic Dipole Transitions

The polarization of the light that is scattered by the coronal ions is influenced by the anisotropic illumination from the photosphere and the magnetic field structuring in the solar corona. The properties of the coronal magnetic fields can be well studied by understanding the polarization properties of coronal forbidden emission lines that arise from magnetic dipole ( M 1) transitions in the highly ionized atoms that are present in the corona. We present the classical scattering theory of the forbidden lines for a more general case of arbitrary-strength magnetic fields. We derive the scattering matrix for M 1 transitions using the classical magnetic dipole model of Casini and Lin and applying the scattering matrix approach of Stenflo. We consider a two-level atom model and neglect collisional effects. The scattering matrix so derived is used to study the Stokes profiles formed in coronal conditions in those regions where the radiative excitations dominate collisional excitations. To this end, we take into account the integration over a cone of an unpolarized radiation from the solar disk incident on the scattering atoms. Furthermore, we also integrate along the line of sight to calculate the emerging polarized line profiles. We consider radial and dipole magnetic field configurations and spherically symmetric density distributions. For our studies we adopt the atomic parameters corresponding to the [Fe xiii] 10747 Å coronal forbidden line. We also discuss the nature of the scattering matrix for M 1 transitions and compare it with that for the electric dipole ( E 1) transitions.

79 ASTRONOMY AND ASTROPHYSICS↗

Correlation functions from tensor network influence functionals: The case of the spin-boson model

We investigate the application of matrix product state (MPS) representations of the influence functionals (IFs) for the calculation of real-time equilibrium correlation functions in open quantum systems. Focusing specifically on the unbiased spin-boson model, we explore the use of IF-MPSs for complex time propagation, as well as IF-MPSs for constructing correlation functions in the steady state. We examine three different IF approaches: one based on the Kadanoff–Baym contour targeting correlation functions at all times, one based on a complex contour targeting the correlation function at a single time, and a steady state formulation, which avoids imaginary or complex times, while providing access to correlation functions at all times. We show that within the IF language, the steady state formulation provides a powerful approach to evaluate equilibrium correlation functions.

Chemistry↗

Sphere and disk partition functions in Liouville and in matrix integrals

We compute the sphere and disk partition functions in semiclassical Liouville and analogous quantities in double-scaled matrix integrals. The quantity sphere/disk2 is unambiguous and we find a precise numerical match between the Liouville answer and the matrix integral answer. An application is to show that the sphere partition function in JT gravity is infinite.

2D Gravity↗

A study of pilot modeling in multi-controller tasks

A modeling approach, which utilizes a matrix of transfer functions to describe the human pilot in multiple input, multiple output control situations, is studied. The approach used was to extend a well established scalar Wiener-Hopf minimization technique to the matrix case and then study, via a series of experiments, the data requirements when only finite record lengths are available. One of these experiments was a two-controller roll tracking experiment designed to force the pilot to use rudder in order to coordinate and reduce the effects of aileron yaw. One model was computed for the case where the signals used to generate the spectral matrix are error and bank angle while another model was computed for the case where error and yaw angle are the inputs. Several anomalies were observed to be present in the experimental data. These are defined by the descriptive terms roll up, break up, and roll down. Due to these algorithm induced anomalies, the frequency band over which reliable estimates of power spectra can be achieved is considerably less than predicted by the sampling theorem.

Whitbeck, R. F.↗

A two-phase model that unifies and extends the classical models of membrane transport

Two models describe solvent transport through swollen, nonporous membranes. The pore-flow model, based on fluid mechanics, works for porous membranes, whereas the solution-diffusion model invokes molecular diffusion to treat nonporous membranes. Both approaches make valid arguments for swollen polymer membranes, but they disagree in their predictions of intramembrane pressure and concentration profiles. Using a fluid-solid model that treats the solvent and membrane matrix as separate phases, we show both classical models to be valid, to represent complementary approaches to the same phenomenon, and to make identical predictions. The fluid-solid model clarifies recent reverse osmosis measurements; provides a predictive and mechanistic basis for empirical high-pressure limiting flux phenomena, in quantitative agreement with classic measurements; and gives a framework to treat nonporous but mechanically heterogeneous membrane materials.

Science & Technology - Other Topics↗

Gogny-force-derived effective shell-model Hamiltonian

The density-dependent finite-range Gogny force has been used to derive the effective Hamiltonian for the shell-model calculations of nuclei. The density dependence simulates an equivalent three-body force, while the finite range gives a Gaussian distribution of the interaction in the momentum space and hence leads to an automatic smooth decoupling between low-momentum and high-momentum components of the interaction, which is important for finite-space shell-model calculations. Two-body interaction matrix elements, single-particle energies, and the core energy of the shell model can be determined by the unified Gogny force. The analytical form of the Gogny force is advantageous to treat cross-shell cases, while it is difficult to determine the cross-shell matrix elements and single-particle energies using an empirical Hamiltonian by fitting experimental data with a large number of matrix elements. In this paper, we have applied the Gogny-force effective shell-model Hamiltonian to the p- and sd-shell nuclei. The results show good agreements with experimental data and other calculations using empirical Hamiltonians. The experimentally known neutron drip line of oxygen isotopes and the ground states of typical nuclei 10B and 18N can be reproduced, in which the role of three-body force is non-negligible. Here, the Gogny-force-derived effective Hamiltonian has also been applied to the cross-shell calculations of the sd–pf shell.

20 ≤ A ≤ 38↗

50 years of CP violation — What have we learned?

Early after the discovery of CP violation, the explanation of how the Standard Model of particle physics could allow CP violation was quickly given, but it took many years for the original observation to be unequivocally explained on that basis. It was also proposed that this observation opened up the possibility that we could now explain the fact that the universe is made of matter. Remarkably, 50 years later we have no evidence in particle physics that there is any CP violation except that of the Kobayashi Maskawa mechanism of the standard model. Yet we fail completely to explain the baryon asymmetry of the Universe through that mechanism. After reviewing the main points in the history I describe the present experimental attempts to find CP violation beyond the standard model, and explain the theoretical attempts to explain the matter in the Universe.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Tensile Creep Behavior of SiC(f)/SiC Ceramic Matrix Minicomposites

Single fiber-tow minicomposites represent the major load-bearing element of woven and laminate ceramic matrix composites (CMCs). To understand the effects of fiber type, fiber content, and matrix cracking on tensile creep in SiCf/SiC CMCs, single-tow SiCf/SiC minicomposites with different fiber types and contents were investigated. The minicomposites studied contained either Hi-Nicalon™ or Hi-Nicalon™ Type S SiC fibers with a boron nitride (BN) interphase and a chemical-vapor-infiltrated-silicon-carbide (CVI-SiC) matrix. Tensile creep was performed at 1200 °C in air. A bottom-up creep modeling approach was applied where creep parameters of the fibers and matrix were obtained separately at 1200 °C. Next, a theoretical model based on the rule of mixtures was derived to model the fiber and matrix creep-time-dependent stress redistribution. Fiber and matrix creep parameters, load transfer model results, and numerical modeling were used to construct a creep strain model to predict creep damage evolution of minicomposites with different fiber types and contents.

SiC/SiC Minicomposites Creep↗

Kinetic equations for baryogenesis via sterile neutrino oscillation

We investigate baryogenesis in the νMSM (neutrino Minimal Standard Model), which is the SM extended by three right-handed neutrinos with masses below the electroweak scale. The baryon asymmetry of the universe can be generated by the mechanism via flavor oscillation of right-handed (sterile) neutrinos which are responsible to masses of active neutrinos confirmed by various experiments. We present the kinetic equations for the matrix of densities of leptons which describe the generation of asymmetries. Especially, the momentum dependence of the matrix of densities is taken into account. By solving these equations numerically, it is found that the momentum distribution is significantly distorted from the equilibrium one, since the production for the modes with lower momenta k << T (T is the temperature of the universe) is enhanced, while suppressed for higher modes. As a result, the most important mode for the yields of sterile neutrinos as well as the baryon asymmetry is k ≅ 2T, which is smaller than (k) inferred from the thermal average. The comparison with the previous works is also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

ICAN/PART: Particulate composite analyzer, user's manual and verification studies

A methodology for predicting the equivalent properties and constituent microstresses for particulate matrix composites, based on the micromechanics approach, is developed. These equations are integrated into a computer code developed to predict the equivalent properties and microstresses of fiber reinforced polymer matrix composites to form a new computer code, ICAN/PART. Details of the flowchart, input and output for ICAN/PART are described, along with examples of the input and output. Only the differences between ICAN/PART and the original ICAN code are described in detail, and the user is assumed to be familiar with the structure and usage of the original ICAN code. Detailed verification studies, utilizing dim dimensional finite element and boundary element analyses, are conducted in order to verify that the micromechanics methodology accurately models the mechanics of particulate matrix composites. ne equivalent properties computed by ICAN/PART fall within bounds established by the finite element and boundary element results. Furthermore, constituent microstresses computed by ICAN/PART agree in average sense with results computed using the finite element method. The verification studies indicate that the micromechanics programmed into ICAN/PART do indeed accurately model the mechanics of particulate matrix composites.

Goldberg, Robert K.↗

Steady state RANS simulations of temperature fluctuations in single phase turbulent mixing

Single phase turbulent mixing in nuclear power plant circuits where a strong temperature gradient is present is known to precipitate pipe failure due to thermal fatigue. Experiments in a square mixing channel offer the opportunity to study the phenomenon under simple and easily reproducible boundary conditions. Measurements of this kind have been performed extensively at the Paul Scherrer Inst. in Switzerland with a high density of instrumentation in the Generic Mixing Experiment (GEMIX). As a fundamental mixing phenomena study closely related to the thermal fatigue problem, the experimental results from GEMIX are valuable for the validation of CFD codes striving to accurately simulate both the temperature and velocity fields in single phase turbulent mixing. In the experiments two iso-kinetic streams meet at a shallow angle of 3 degrees and mix in a straight channel of square cross-section under various degrees of density, temperature, and viscosity stratification over a range of Reynolds numbers ranging from 5*10{sup 3} to 1*10{sup 5}. Conductivity measurements, using wire-mesh and wall sensors, as well as optical measurements, using particle image velocimetry, were conducted with high temporal and spatial resolutions (up to 2.5 kHz and 1 mm in the case of the wire mesh sensor) in the mixing zone, downstream of a splitter plate. The present paper communicates the results of RANS modeling of selected GEMIX tests. Steady-state CFD calculations using a RANS turbulence model represent an inexpensive method for analyzing large and complex components in commercial nuclear reactors, such as the downcomer and reactor pressure vessel heads. Crucial to real world applicability, however, is the ability to model turbulent heat fluctuations in the flow; the Turbulent Heat Flux Transport model developed by ANSYS CFX is capable, by implementation of a transport equation for turbulent heat fluxes, of readily modeling these values. Furthermore, the closure of the turbulent heat flux transport equation evokes a transport equation for the variance of the enthalpy. It is therefore possible to compare the modeled fluctuations of the liquid temperature directly with the scalar fluctuations recorded experimentally with the wire-mesh. Combined with a working Turbulent Heat Flux Transport model, complex mixing problems in large geometries could be better understood. We aim for the validation of Reynolds Stress based RANS simulations extended by the Turbulent Heat Flux Transport model by modeling the GEMIX experiments in detail. Numerical modeling has been performed using both BSL and SSG Reynolds Stress Models in a test matrix comprising experimental trials at the GEMIX facility. We expand on the turbulent mixing RANS CFD results of (Manera 2009) in a few ways. In the GEMIX facility we introduce density stratification in the flow while removing the characteristic large scale vorticity encountered in T-junctions and therefore find better conditions to check the diffusive conditions in the model. Furthermore, we study the performance of the model in a very different, simpler scalar fluctuation spectrum. The paper discusses the performance of the model regarding the dissipation of the turbulent kinetic energy and dissipation of the enthalpy variance. A novel element is the analyses of cases with density stratification. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Local-global analysis of crack growth in continuously reinfoced ceramic matrix composites

This paper describes the development of a mathematical model for predicting the strength and micromechanical failure characteristics of continuously reinforced ceramic matrix composites. The local-global analysis models the vicinity of a propagating crack tip as a local heterogeneous region (LHR) consisting of spring-like representation of the matrix, fibers and interfaces. Parametric studies are conducted to investigate the effects of LHR size, component properties, and interface conditions on the strength and sequence of the failure processes in the unidirectional composite system.

Ballarini, Roberto↗