Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “orthogonal weights”

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 73 records · Page 4

Optimal control of spin-stabilized spacecraft with telescoping appendages

The control of a spin-stabilized spacecraft consisting of a rigid central hub and one or two movable offset telescoping booms (with end masses) is considered. The equations of rotational motion are linearized about either of two desired final states. A control law for the boom end mass position is sought such that a quadratic cost functional involving the weighted components of angular velocity plus the control is minimized when the final time is unspecified and involves the solution of the matrix Riccati algebraic equation. For three axis control more than one offset boom (orthogonal to each other) is required. For two-axis control with a single boom offset from a symmetrical hub, an analytic solution is obtained; when this system is used for nutation decay the time constant is one order of magnitude smaller than previously achieved using nonoptimal control logic. For the general case results are obtained numerically.

Bainum, P. M.↗

Tuning Parameters in Heuristics by Using Design of Experiments Methods

With the growing complexity of today's large scale problems, it has become more difficult to find optimal solutions by using exact mathematical methods. The need to find near-optimal solutions in an acceptable time frame requires heuristic approaches. In many cases, however, most heuristics have several parameters that need to be "tuned" before they can reach good results. The problem then turns into "finding best parameter setting" for the heuristics to solve the problems efficiently and timely. One-Factor-At-a-Time (OFAT) approach for parameter tuning neglects the interactions between parameters. Design of Experiments (DOE) tools can be instead employed to tune the parameters more effectively. In this paper, we seek the best parameter setting for a Genetic Algorithm (GA) to solve the single machine total weighted tardiness problem in which n jobs must be scheduled on a single machine without preemption, and the objective is to minimize the total weighted tardiness. Benchmark instances for the problem are available in the literature. To fine tune the GA parameters in the most efficient way, we compare multiple DOE models including 2-level (2k ) full factorial design, orthogonal array design, central composite design, D-optimal design and signal-to-noise (SIN) ratios. In each DOE method, a mathematical model is created using regression analysis, and solved to obtain the best parameter setting. After verification runs using the tuned parameter setting, the preliminary results for optimal solutions of multiple instances were found efficiently.

Arin, Arif↗

Optimal control of spin stabilized spacecraft with telescoping appendages

The control of a spin-stabilized spacecraft consisting of a rigid central hub and one or two movable offset telescoping booms (with end masses) is considered. The equations of rotational motion are linearized about either of two desired final states. A control law for the boom and mass position is sought such that a quadratic cost functional involving the weighted components of angular velocity plus the control is minimized when the final time is unspecified and involves the solution of the matrix Riccati algebraic equation. For three-axis control more than one offset boom (orthogonal to each other) is required. For two-axis control with a single boom offset from a symmetrical hub, an analytic solution is obtained; when this system is used for nutation decay the time constant is one order of magnitude smaller than previously achieved using non-optimal control logic. For the general case results are obtained numerically.

Bainum, P. M.↗

Gram-Schmidt algorithms for covariance propagation

This paper addresses the time propagation of triangular covariance factors. Attention is focused on the square-root free factorization, P = UDU/T/, where U is unit upper triangular and D is diagonal. An efficient and reliable algorithm for U-D propagation is derived which employs Gram-Schmidt orthogonalization. Partitioning the state vector to distinguish bias and colored process noise parameters increases mapping efficiency. Cost comparisons of the U-D, Schmidt square-root covariance and conventional covariance propagation methods are made using weighted arithmetic operation counts. The U-D time update is shown to be less costly than the Schmidt method; and, except in unusual circumstances, it is within 20% of the cost of conventional propagation.

Thornton, C. L.↗

Gram-Schmidt algorithms for covariance propagation

This paper addresses the time propagation of triangular covariance factors. Attention is focused on the square-root free factorization, P = UD(transpose of U), where U is unit upper triangular and D is diagonal. An efficient and reliable algorithm for U-D propagation is derived which employs Gram-Schmidt orthogonalization. Partitioning the state vector to distinguish bias and coloured process noise parameters increase mapping efficiency. Cost comparisons of the U-D, Schmidt square-root covariance and conventional covariance propagation methods are made using weighted arithmetic operation counts. The U-D time update is shown to be less costly than the Schmidt method; and, except in unusual circumstances, it is within 20% of the cost of conventional propagation.

Thornton, C. L.↗

Structured adaptive grid generation using algebraic methods

The accuracy of the numerical algorithm depends not only on the formal order of approximation but also on the distribution of grid points in the computational domain. Grid adaptation is a procedure which allows optimal grid redistribution as the solution progresses. It offers the prospect of accurate flow field simulations without the use of an excessively timely, computationally expensive, grid. Grid adaptive schemes are divided into two basic categories: differential and algebraic. The differential method is based on a variational approach where a function which contains a measure of grid smoothness, orthogonality and volume variation is minimized by using a variational principle. This approach provided a solid mathematical basis for the adaptive method, but the Euler-Lagrange equations must be solved in addition to the original governing equations. On the other hand, the algebraic method requires much less computational effort, but the grid may not be smooth. The algebraic techniques are based on devising an algorithm where the grid movement is governed by estimates of the local error in the numerical solution. This is achieved by requiring the points in the large error regions to attract other points and points in the low error region to repel other points. The development of a fast, efficient, and robust algebraic adaptive algorithm for structured flow simulation applications is presented. This development is accomplished in a three step process. The first step is to define an adaptive weighting mesh (distribution mesh) on the basis of the equidistribution law applied to the flow field solution. The second, and probably the most crucial step, is to redistribute grid points in the computational domain according to the aforementioned weighting mesh. The third and the last step is to reevaluate the flow property by an appropriate search/interpolate scheme at the new grid locations. The adaptive weighting mesh provides the information on the desired concentration of points to the grid redistribution scheme. The evaluation of the weighting mesh is accomplished by utilizing the weight function representing the solution variation and the equidistribution law. The selection of the weight function plays a key role in grid adaptation. A new weight function utilizing a properly weighted boolean sum of various flowfield characteristics is defined. The redistribution scheme is developed utilizing Non-Uniform Rational B-Splines (NURBS) representation. The application of NURBS representation results in a well distributed smooth grid by maintaining the fidelity of the geometry associated with boundary curves. Several algebraic methods are applied to smooth and/or nearly orthogonalize the grid lines. An elliptic solver is utilized to smooth the grid lines if there are grid crossings. Various computational examples of practical interest are presented to demonstrate the success of these methods.

Yang, Jiann-Cherng↗

Polymers depolymerizable by metathesis of a cleavable unit

Novel polymers are depolymerizable by metathesis of a cleavable unit. As an example, a series of linear and crosslinked polyurethanes were prepared that can be selectively depolymerized under mild conditions. Two unique polyols were synthesized bearing unsaturated units in a configuration designed to favor ring-closing metathesis to five- and six-membered cycloalkenes. These polyols were co-polymerized with toluene diisocyanate to generate linear polyurethanes and trifunctional hexamethylene- and diphenylmethane-based isocyanates to generate crosslinked polyurethanes. The polyol design is such that the ring-closing metathesis reaction cleaves the backbone of the polymer chain. Upon exposure to dilute solutions of Grubbs' catalyst under ambient conditions, the polyurethanes were rapidly depolymerized to low molecular weight, soluble products bearing vinyl and cycloalkene functionalities. These functionalities enabled further re-polymerization by traditional strategies for polymerization of double bonds. This general approach can be expanded to develop a range of chemically recyclable condensation polymers that are readily depolymerized by orthogonal metathesis chemistry.

Jones, Brad Howard↗

The Influence of Aqueous Dispersivity on Transport of Radiological Material at LANL MDA G - 20450

Aqueous dispersivity is a phenomenon by which a solute moving in groundwater mixes with water that does not contain the solute and becomes diluted. Solutes take paths of varying lengths through pore networks or experience slower velocities due to interactions with pore walls. Aqueous dispersivity leads to breakthrough curves with earlier arrival than expected in purely advective systems and persistent 'tails,' or late-arriving, dilute concentrations of the solute. Dispersivity can be longitudinal (in the direction of flow) or transverse (orthogonal to flow). LANL has disposed of radioactive material at Material Disposal Area (MDA) G since the late 1950's. Most of the disposed materials at MDA G are low-level waste (LLW), but other types of radioactive wastes are disposed or stored at the site. For example, there is a collection of remote-handled transuranic waste (RH-TRU) waste in 33 shafts, which were intended to be stored retrievably. Aqueous-phase radionuclides included in models of contaminant transport at MDA G are subject to dispersive processes in porous media of the vadose zone and saturated zone. To model dispersivity at MDA G probabilistically, we created site-specific distributions of aqueous dispersivity. First, we performed a comprehensive review of peer-reviewed publications, reports from other performance assessments of MDA G, and documentation from sites that were considered similar based on expert opinion, such as Yucca Mountain. Data were included from a wide range of sources, so that distributions developed based on the data reflect the current state of knowledge in properties of the subsurface environment. Data were grouped based on value type (e.g., determined via experiment or used in a different model), rock type, zone type (saturated or unsaturated), and site, and were weighted in the distributions according to their relevance to the site and zone type. We combined Monte Carlo simulation and bootstrapping techniques to produce approximate distributions of mean aqueous dispersivity. Distributions of dispersivity for saturated and unsaturated zone transport are differentiated by data weight and the length scales over which transport occurs. We find that underestimating or disregarding transverse dispersivity may lead to overstating simulated well capture in the saturated zone at MDA G. Neglecting transverse dispersivity is a common assumption due in part to use of one-dimensional transport models that do not consider solute spreading orthogonal to flow. Also, we find that variability in longitudinal dispersivity is associated with substantial travel time alterations in the MDA G unsaturated zone and in simulated contaminant breakthrough times at a downstream observation plane in the MDA G saturated zone. Finally, at large scales (>1000 m), a linear trend between the dispersivity coefficient and length over which transport occurs may overestimate the dispersivity coefficient. Future modeling work with application to solute spreading could include interacting advective and diffusive zones in the highly fractured basalt below MDA G. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Imaging proportional counter for X-ray astronomy

A position-sensitive proportional counter capable of imaging X-rays (0.1-3 keV) over a 10 cm x 10 cm aperture has been constructed. Positioning is obtained by sensing the signals induced by an X-ray event on the two orthogonal sets of cathode wires. Each cathode is divided into a series of cathode strips, each 0.5 cm wide. An X-ray event induces signals on several adjacent strips. Signals from each cathode strip are amplified separately and then added in an equally weighted and an unequally weighted summing amplifier. The position in each direction is obtained by dividing the output of the unequally weighted summer by that of the equally weighted summer. At 0.94 keV, the accuracy of the position sensing is 190 microns. At the same energy, the energy resolution is approximately 65% (FWHM). The proportional counter system is currently being incorporated into a sounding rocket payload having metal mirror optics, which is being constructed by the X-ray astronomy group at the California Institute of Technology.

Long, K. S.↗

An Algorithm to Generate Synthetic 3D Microstructures from 2D Exemplars

The inverse problem of constructing 3D microstructures from 2D data is an area of active research within the materials science community. This paper presents the implementation of a robust, computationally efficient algorithm, Hierarchical Algorithm for the Reconstruction of Exemplars (HARE), written in Python to reconstruct 3D features in a given microstructure from up to three orthogonal 2D exemplars using nearest neighbor matching to reproduce feature qualities, such as shape, size, and distribution. HARE’s feature sampling implements histogram reweighting to avoid both over- and undersampling. A neighborhood voting scheme allows each pixel to provisionally affect its neighbors according to its weight. Finally, the algorithm is presently configured for two-phase materials and is being extended to accommodate multiple phases. HARE is a convenient and robust base from which to generate statistically representative synthetic microstructures for use in multi-scale modeling or machine learning applications to support advanced manufacturing and materials discovery.

36 MATERIALS SCIENCE↗

Selectively Depolymerizable Polyurethanes from Unsaturated Polyols Cleavable by Olefin Metathesis

This communication describes a novel series of linear and crosslinked polyurethanes (PUs) and their selective depolymerization under mild conditions. Two unique polyols are synthesized bearing unsaturated units in a configuration designed to favor ring-closing metathesis (RCM) to five- and six-membered cycloalkenes. These polyols are co-polymerized with toluene diisocyanate to generate linear PUs and trifunctional hexamethylene- and diphenylmethane-based isocyanates to generate crosslinked PUs. The polyol design is such that the RCM reaction cleaves the backbone of the polymer chain. Upon exposure to dilute solutions of Grubbs’ catalyst under ambient conditions, the PUs are rapidly depolymerized to low molecular weight, soluble products bearing vinyl and cycloalkene functionalities. These functionalities enable further re-polymerization by traditional strategies for polymerization of double bonds. We anticipate that this general approach can be expanded to develop a range of chemically recyclable condensation polymers that are readily depolymerized by orthogonal metathesis chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Collection, Usage, and Preliminary Examination of the Apollo Sample Suite: Lessons for Artemis

Apollo Sample Collection and Usage: From 1969 to 1972 there were six Apollo missions to the surface of the Moon during which the astronauts collected 382 kg of rock and regolith (~2200 samples). The samples collected fall into these general categories: rocks (~66% by mass), rake samples (~4%), bulk regolith (~24%), and specialty regolith (deep drill cores, drive tubes, sealed bulk regolith) samples (~6%). In each category there are a variety of different subtypes available for study, e.g., among the bulk regolith samples there are also skim, trench, and (partially) shaded regolith samples each sampling unique types or depths of regolith. This variety of subsamples has enabled a multitude of different studies over the past 55 years (>3400 individual requests). We are still averaging ~50 unique requests and have allocated >500 individual Apollo samples annually for the past 10 years (2020 excepted). Looking at the 4,675 non-ANGSA (Apollo Next Generation Sample Analysis) samples allocated over the past 10 years, the proportions of allocated samples do not precisely align with the abundance (by mass) of those samples withing the collection: Rock (69.4 %); Rake (10.8 %); Bulk Regolith (15.5 %); Drive Tube (3.3 %); Core/Specialty (1.0 %). Apollo Preliminary Examination (PE): The PE process differs significantly for the various sample types enumerated above; we focus on regolith and rock samples here. During the Apollo mission era, the PE process evolved over the course of the missions; below is what was done for the Apollo 17 mission. For regolith samples, the PE process was: (1) documented bags containing regolith are opened, photographed, and described; (2) large rocks are removed and treated separately; (3) 25% to 33% of the bulk soil is scooped out, weighed, and stored in reserve; (4) the remaining sample is sieved to produce the size fractions <1, 1- 2, 2-4, and 4-10 mm, all of which are weighed. For rock samples, the process is: (1) removing rocks from the container(s) it came back from the Moon in; (2) rematching any materials that spalled off the rock to their original location; (3) numbering, weighing, and basic photographic documentation; (4) dusting with a gentle N2 gas jet; (5) Orthogonal photography; (6) detailed description of the textures and features of the rock; (7) rock modelling and measurement; (8) stereophotography; (9) determination of the orientation of the rock on the lunar surface. Drive tubes and deep drill cores were not characterized during PE beyond an initial weight and a sketch of the interior tube materials derived from 2D medical X-ray images. Catalogs: The ongoing utility of the Apollo samples is enabled by the robust cataloguing process for the samples [3-5], which allows the scientific community to accurately request samples uniquely suited to their proposed studies. A common misconception, however, is the amount of detail that goes into the initial catalog (e.g., [6]) for a collection from the preliminary examination (PE) period, versus what goes into the catalogs that come later in the life cycle of the samples from that mission (e.g., [7]). The only required data for a PE catalog is a weight, a basic photograph, and a description of the nature of the sample. Artemis PE: Over the past few years, the ongoing ANGSA project studied previously unopened Apollo 17 double drive tube samples 73001/2 [1], and a PE of the drive tubes was done. The PE took the existing core dissection process (developed during PE of Apollo cores in the 1970s, 1980s), and modernized it [7]. The main lesson from the ANGSA PE relative to future missions was that the physical work done during PE of lunar samples has not changed much over the past 5 decades. The use of “modern” technology during PE (e.g., XCT; multispectral analyses) resulted in an enhanced initial understanding of the 73001 and 73002 drive tubes, but greatly increased the time required. The lessons learned from recent astromaterial PEs (e.g., ANGSA and OREx) are important to consider when planning for Artemis, but the unique nature of the Artemis Campaign means many lessons learned from these mission will not be applicable. Given the time constraints (6 months) and likely number of samples that will be returned by Artemis (>200), the Artemis PE catalog will necessarily look much more like [4] than [6].

J Gross↗

Optimal design of structures with buckling constraints.

The paper presents an iterative, finite element method for minimum weight design of structures with respect to buckling constraints. The redesign equation is derived from the optimality criterion, as opposed to a numerical search procedure, and can handle problems that are characterized by the existence of two fundamental buckling modes at the optimal design. Application of the method is illustrated by beam and orthogonal frame design problems.

Kiusalaas, J.↗

Complementarity Between Neutrino Neutral and Charged Current Events in the Search for New Physics

At long-baseline neutrino experiments, neutral-current (NC) events accumulate in large numbers but are seldom exploited for new physics searches. We demonstrate their potential using non-standard neutrino interactions (NSI) with quarks as a case study. Charged-current (CC) analyses constrain NSI through matter effects on neutrino propagation, which probe almost exclusively the isoscalar combination of up- and down-quark couplings; the orthogonal isovector combination is suppressed by a factor of $\sim$100. Because NSI also modify NC cross sections in a flavor-dependent way, NC events become sensitive to oscillations: the far-to-near detector ratio acquires a dependence on the beam's flavor composition that probes both isoscalar and isovector couplings with comparable weight. Using existing NOvA data and DUNE projections, we derive the first bounded constraints on isovector NSI from a long-baseline experiment and show that combining CC and NC measurements resolves the individual quark couplings, breaking a degeneracy that persists in either analysis alone.

Gehrlein, Julia [Colorado State U.] (ORCID:0000000↗

Conservative projection-based data-driven model order reduction of a fluid-kinetic spectral solver

Kinetic simulations are computationally intensive due to six-dimensional phase space discretization. Many kinetic spectral solvers use the asymmetrically weighted Hermite expansion due to its conservation and fluid-kinetic coupling properties, i.e., the lower-order Hermite moments capture and describe the macroscopic fluid dynamics, and higher-order Hermite moments describe the microscopic kinetic dynamics. We leverage this structure by developing a parametric data-driven reduced-order model based on the proper orthogonal decomposition, which projects the higher-order kinetic moments while retaining the fluid moments intact. We demonstrate analytically and numerically that the method ensures local and global mass, momentum, and energy conservation. The numerical results show that the proposed method effectively replicates the high-dimensional spectral simulations at a fraction of the computational cost and memory, as validated on the weak Landau damping and two-stream instability benchmark problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Application of Taguchi methods to dual mixture ratio propulsion system optimization for SSTO vehicles

The application of advanced technologies to future launch vehicle designs would allow the introduction of a rocket-powered, single-stage-to-orbit (SSTO) launch system early in the next century. For a selected SSTO concept, a dual mixture ratio, staged combustion cycle engine that employs a number of innovative technologies was selected as the baseline propulsion system. A series of parametric trade studies are presented to optimize both a dual mixture ratio engine and a single mixture ratio engine of similar design and technology level. The effect of varying lift-off thrust-to-weight ratio, engine mode transition Mach number, mixture ratios, area ratios, and chamber pressure values on overall vehicle weight is examined. The sensitivity of the advanced SSTO vehicle to variations in each of these parameters is presented, taking into account the interaction of each of the parameters with each other. This parametric optimization and sensitivity study employs a Taguchi design method. The Taguchi method is an efficient approach for determining near-optimum design parameters using orthogonal matrices from design of experiments (DOE) theory. Using orthogonal matrices significantly reduces the number of experimental configurations to be studied. The effectiveness and limitations of the Taguchi method for propulsion/vehicle optimization studies as compared to traditional single-variable parametric trade studies is also discussed.

Stanley, Douglas O.↗

Variable Acceleration Force Calibration System (VACS)

Conventionally, force balances have been calibrated manually, using a complex system of free hanging precision weights, bell cranks, and/or other mechanical components. Conventional methods may provide sufficient accuracy in some instances, but are often quite complex and labor-intensive, requiring three to four man-weeks to complete each full calibration. To ensure accuracy, gravity-based loading is typically utilized. However, this often causes difficulty when applying loads in three simultaneous, orthogonal axes. A complex system of levers, cranks, and cables must be used, introducing increased sources of systematic error, and significantly increasing the time and labor intensity required to complete the calibration. One aspect of the VACS is a method wherein the mass utilized for calibration is held constant, and the acceleration is changed to thereby generate relatively large forces with relatively small test masses. Multiple forces can be applied to a force balance without changing the test mass, and dynamic forces can be applied by rotation or oscillating acceleration. If rotational motion is utilized, a mass is rigidly attached to a force balance, and the mass is exposed to a rotational field. A large force can be applied by utilizing a large rotational velocity. A centrifuge or rotating table can be used to create the rotational field, and fixtures can be utilized to position the force balance. The acceleration may also be linear. For example, a table that moves linearly and accelerates in a sinusoidal manner may also be utilized. The test mass does not have to move in a path that is parallel to the ground, and no re-leveling is therefore required. Balance deflection corrections may be applied passively by monitoring the orientation of the force balance with a three-axis accelerometer package. Deflections are measured during each test run, and adjustments with respect to the true applied load can be made during the post-processing stage. This paper will present the development and testing of the VASC concept.

Rhew, Ray D.↗