Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ordering transformation”

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,441 records · Page 80

The H(sub infinity) optimal controller design and reduction for the inertial hold mode of the attitude control system of the XTE spacecraft

The Inertial Hold Mode (IHM) is one mode of the attitude control system of the X-ray Timing Explorer spacecraft that is disturbed by both parametric uncertainties and external torque disturbance. The IHM model is modified into a typical H-infinity mixed-sensitivity problem through choosing suitable weighting functions W(sub 1)(s) and W(sub 3)(s). The controller is designed by the H-infinity optimization technique with the transformation of shifting the imaginary axis. It can stabilize the plant with uncertainties from the natural frequencies of the flexible body. The gain margin and phase margin of the system are 24.03 db and 55.04 deg, respectively. The step response attenuates to zero within 150 seconds. These show that the controller satisfies the specified requirements. Since the order of the controller appears high, it is reduced to fourth order one. The results show that the stability and the performance of the system with the reduced controller are retained perfectly.

Xu, Zhong Ling↗

The Ocean Colour Climate Change Initiative: I. A Methodology for Assessing Atmospheric Correction Processors Based on In-Situ Measurements

The Ocean Colour Climate Change Initiative intends to provide a long-term time series of ocean colour data and investigate the detectable climate impact. A reliable and stable atmospheric correction procedure is the basis for ocean colour products of the necessary high quality. In order to guarantee an objective selection from a set of four atmospheric correction processors, the common validation strategy of comparisons between in-situ and satellite derived water leaving reflectance spectra, is extended by a ranking system. In principle, the statistical parameters such as root mean square error, bias, etc. and measures of goodness of fit, are transformed into relative scores, which evaluate the relationship of quality dependent on the algorithms under study. The sensitivity of these scores to the selected database has been assessed by a bootstrapping exercise, which allows identification of the uncertainty in the scoring results. Although the presented methodology is intended to be used in an algorithm selection process, this paper focusses on the scope of the methodology rather than the properties of the individual processors.

Climate Change↗

A new concept of stability in orbit propagation, useful for quantifying numerical errors

We present the concept of topological stability in the numerical propagation of orbits, and show how it results in a useful new method for measuring the global numerical error of an orbit propagation. The concept applies to any problem in orbital dynamics. Moreover, it can be extended to any three-dimensional system of di erential equations of second order. In order to assess the topological stability of a given integration a special metric is introduced, which can be used to estimate the numerical errors robustly. The method is particularly well suited for dealing with strongly perturbed and chaotic systems. The construction is based on the constraint imposed by the Hopf map that supports the Kustaanheimo-Stiefel transformation. Generic concepts of stability are translated to KS space.

Pelaez, Jesus↗

Sparse Cholesky factorization for solving nonlinear PDEs via Gaussian processes

In recent years, there has been widespread adoption of machine learning-based approaches to automate the solving of partial differential equations (PDEs). Among these approaches, Gaussian processes (GPs) and kernel methods have garnered considerable interest due to their flexibility, robust theoretical guarantees, and close ties to traditional methods. They can transform the solving of general nonlinear PDEs into solving quadratic optimization problems with nonlinear, PDE-induced constraints. However, the complexity bottleneck lies in computing with dense kernel matrices obtained from pointwise evaluations of the covariance kernel, and its partial derivatives, a result of the PDE constraint and for which fast algorithms are scarce. The primary goal of this paper is to provide a near-linear complexity algorithm for working with such kernel matrices. We present a sparse Cholesky factorization algorithm for these matrices based on the near-sparsity of the Cholesky factor under a novel ordering of pointwise and derivative measurements. The near-sparsity is rigorously justified by directly connecting the factor to GP regression and exponential decay of basis functions in numerical homogenization. We then employ the Vecchia approximation of GPs, which is optimal in the Kullback-Leibler divergence, to compute the approximate factor. This enables us to compute ϵ-approximate inverse Cholesky factors of the kernel matrices with complexity O(N log d (N/ϵ)) in space and O(N log 2d (N/ϵ)) in time. We integrate sparse Cholesky factorizations into optimization algorithms to obtain fast solvers of the nonlinear PDE. We numerically illustrate our algorithm’s near-linear space/time complexity for a broad class of nonlinear PDEs such as the nonlinear elliptic, Burgers, and Monge-Ampère equations. In summary, we provide a fast, scalable, and accurate method for solving general PDEs with GPs and kernel methods.

97 MATHEMATICS AND COMPUTING↗

Ion diffusion retarded by diverging chemical susceptibility

Abstract For first-order phase transitions, the second derivatives of Gibbs free energy (specific heat and compressibility) diverge at the transition point, resulting in an effect known as super-elasticity along the pressure axis, or super-thermicity along the temperature axis. Here we report a chemical analogy of these singularity effects along the atomic doping axis, where the second derivative of Gibbs free energy (chemical susceptibility) diverges at the transition point, leading to an anomalously high energy barrier for dopant diffusion in co-existing phases, an effect we coin as super-susceptibility. The effect is realized in hydrogen diffusion in vanadium dioxide (VO 2 ) with a metal-insulator transition (MIT). We show that hydrogen faces three times higher energy barrier and over one order of magnitude lower diffusivity when it diffuses across a metal-insulator domain wall in VO 2 . The additional energy barrier is attributed to a volumetric energy penalty that the diffusers need to pay for the reduction of latent heat. The super-susceptibility and resultant retarded atomic diffusion are expected to exist universally in all phase transformations where the transformation temperature is coupled to chemical composition, and inspires new ways to engineer dopant diffusion in phase-coexisting material systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Parametric tuning of quantum phase transitions in ultracold reactions

Advances in atomic physics have led to the possibility of a coherent transformation between ultracold atoms and molecules including between completely bosonic condensates. Such transformations are enabled by the magneto-association of atoms at a Feshbach resonance which results in a passage through a quantum critical point. In this study, we show that the presence of generic interaction between the constituent atoms and molecules can fundamentally alter the nature of the critical point, change the yield of the reaction and the order of the consequent phase transition. We find that the correlations introduced by this interaction induce nontrivial many-body physics such as coherent oscillations between atoms and molecules, and a selective formation of squeezed molecular quantum states and quantum cat states. We provide analytical and numerical descriptions of these effects, along with scaling laws for the reaction yield in non-adiabatic regimes.

36 MATERIALS SCIENCE↗

Low-dimensional physics of clay particle size distribution and layer ordering

Abstract Clays are known for their small particle sizes and complex layer stacking. We show here that the limited dimension of clay particles arises from the lack of long-range order in low-dimensional systems. Because of its weak interlayer interaction, a clay mineral can be treated as two separate low-dimensional systems: a 2D system for individual phyllosilicate layers and a quasi-1D system for layer stacking. The layer stacking or ordering in an interstratified clay can be described by a 1D Ising model while the limited extension of individual phyllosilicate layers can be related to a 2D Berezinskii–Kosterlitz–Thouless transition. This treatment allows for a systematic prediction of clay particle size distributions and layer stacking as controlled by the physical and chemical conditions for mineral growth and transformation. Clay minerals provide a useful model system for studying a transition from a 1D to 3D system in crystal growth and for a nanoscale structural manipulation of a general type of layered materials.

42 ENGINEERING↗

Proactive Data Containers for Scientific Storage (Final Report)

Emerging HPC systems are expected to be deployed with an unprecedented level of complexity, due to a deep system memory/storage hierarchy and heterogeneity of the storage hardware. This hierarchy is expected to range from CPU cache through several levels of volatile memory to nonvolatile memory, traditional hard disks, and tape. Simple and efficient methods of data management and movement through this hierarchy is critical for scientific applications using exascale systems. Existing storage system and I/O (SSIO) technologies face severe challenges in dealing with these requirements. POSIX and MPI I/O standards that are the basis for existing I/O libraries and parallel file systems present fundamental challenges in the areas of scalable metadata operations, semantics-based data movement performance tuning, asynchronous operation, and support for scalable consistency of distributed operations. Moving toward new paradigms for SSIO in the extreme-scale era, we have proposed to investigate novel object-based data abstractions and storage mechanisms that take advantage of the deep storage hierarchy and enable proactive automated performance tuning. In order to achieve these overarching goals, we initiated an effort to develop a fundamental new data abstraction, called Proactive Data Containers (PDC). A PDC is a container within a locus of storage (memory, NVRAM, disk, etc.) that stores science data in an object-centric manner. Managing data as objects enables powerful optimization opportunities for data movement and transformations. The R&D focus of this project are: 1) formulation of object-oriented PDCs and their mapping in different levels of the exascale storage hierarchy; 2) efficient strategies for moving data in deep storage hierarchies using PDCs; 3) techniques for transforming and reorganizing data based on application requirements; and 4) novel analysis paradigms for enabling data transformations and user-defined analysis on data in PDCs. Toward achieving these overarching goals, we designed an object-centric application programing interface (API) for HPC, scalable metadata management for object-centric storage systems, and data movement optimizations such as Data Elevator for moving data between two levels of storage devices and TAPIOCA for efficient aggregation of data on compute nodes. We then implemented several components of the PDC system. They include metadata management, data placement services, remote procedure calls, data aggregation, etc. We have put them together into the overall PDC framework.

97 MATHEMATICS AND COMPUTING↗

A critique of some recent second-order turbulence closure models for compressible boundary layers

Computations based on two recently developed second-order turbulence closure models are compared with a series of boundary-layer experiments and with predictions of these experiments using an algebraic mixing length model. One of the models employs an eddy viscosity, whereas the other evaluates components of the Reynolds stress tensor. For flat plates, the computations are compared with the van Driest skin-friction transformation to assess the handling of compressibility. For boundary layers in pressure gradients, four experiments at Mach 4 and one at Mach 6.7 are used as the bases for comparison. In general, both models represent mean velocities and skin friction reasonably well, but represent the turbulence shear stress less accurately.

Rubesin, M. W.↗

Sulfur isotope measurements of submicrometer sulfate aerosol particles over the Pacific Ocean

Stable isotopes were used to analyze the submicron-size sulfate aerosol particles in the atmosphere over the Pacific Ocean, together with the air-mass back trajectories, in order to test the hypothesis of Charlson et al. (1987) who suggested that, over the remote ocean areas, the primary source of atmospheric nonseasalt (NSS) sulfate is marine emissions of dimethylsulfide (DMS). The observed results of isotopic fractionation between the seawater sulfate and NSS sulfate fractions was found to be consistent with the isotopic fractionation predicted for the transformation of the seawater sulfate to the atmospheric NSS sulfate via a DMS path way, supporting the hypothesis of Charlson et al.

Calhoun, Julie A.↗

Connectivity strategies for higher-order neural networks applied to pattern recognition

Different strategies for non-fully connected HONNs (higher-order neural networks) are discussed, showing that by using such strategies an input field of 128 x 128 pixels can be attained while still achieving in-plane rotation and translation-invariant recognition. These techniques allow HONNs to be used with the larger input scenes required for practical pattern-recognition applications. The number of interconnections that must be stored has been reduced by a factor of approximately 200,000 in a T/C case and about 2000 in a Space Shuttle/F-18 case by using regional connectivity. Third-order networks have been simulated using several connection strategies. The method found to work best is regional connectivity. The main advantages of this strategy are the following: (1) it considers features of various scales within the image and thus gets a better sample of what the image looks like; (2) it is invariant to shape-preserving geometric transformations, such as translation and rotation; (3) the connections are predetermined so that no extra computations are necessary during run time; and (4) it does not require any extra storage for recording which connections were formed.

Spirkovska, Lilly↗

An Upgrade of the Aeroheating Software "MINIVER"

Many software packages assist engineers with performing flight vehicle analysis, but some of these packages have gone many years without updates or significant improvements to their workflows. One such software, known as MINIVER, is a powerful yet lightweight tool that is used for aeroheating analyses. However, it is an aging program that has not seen major improvements within the past decade. As part of a collaborative effort with Florida Institute of Technology, MINIVER has received a major user interface overhaul, a change in program language, and will be continually receiving updates to improve its capabilities. The user interface update includes a migration from a command-line interface to that of a graphical user interface supported in the Windows operating system. The organizational structure of the preprocessor has been transformed to clearly defined categories to provide ease of data entry. Helpful tools have been incorporated, including the ability to copy sections of cases as well as a generalized importer which aids in bulk data entry. A visual trajectory editor has been included, as well as a CAD Editor which allows the user to input simplified geometries in order to generate MINIVER cases in bulk. To demonstrate its continued effectiveness, a case involving the JAXA OREX flight vehicle will be included, providing comparisons to captured flight data as well as other computational solutions. The most recent upgrade effort incorporated the use of the CAD Editor, and current efforts are investigating methods to link MINIVER projects with SINDA/Fluint and Thermal Desktop.

Louderback, Pierce M.↗

A Upgrade of the Aeroheating Software "MINIVER"

Many software packages assist engineers with performing flight vehicle analysis, but some of these packages have gone many years without updates or significant improvements to their workflows. One such software package, known as MINIVER, is a powerful yet lightweight tool used for aeroheating analyses. However, it is an aging program that has not seen major improvements within the past decade. As part of a collaborative effort with the Florida Institute of Technology, MINIVER has received a major user interface overhaul, a change in program language, and will be continually receiving updates to improve its capabilities. The user interface update includes a migration from a command-line interface to that of a graphical user interface supported in the Windows operating system. The organizational structure of the pre-processor has been transformed to clearly defined categories to provide ease of data entry. Helpful tools have been incorporated, including the ability to copy sections of cases as well as a generalized importer which aids in bulk data entry. A visual trajectory editor has been included, as well as a CAD Editor which allows the user to input simplified geometries in order to generate MINIVER cases in bulk. To demonstrate its continued effectiveness, a case involving the JAXA OREX flight vehicle will be included, providing comparisons to captured flight data as well as other computational solutions. The most recent upgrade effort incorporated the use of the CAD Editor, and current efforts are investigating methods to link MINIVER projects with SINDA/Fluint and Thermal Desktop.

Louderback, Pierce↗

Low temperature phase transformations in the metallic phases of iron and stony-iron meteorites

The nickel content and the structure of kamacite and decomposed taenite (clear taenite 1, CT-1; the cloudy zone, CZ; and clear taenite 2, CT-2) in the metallic phases of meteorites were determined using X-ray microanalysis techniques in the AEM. The kamacite near the CT-1 interface was found to contain about 4 wt pct Ni. The CT-1 structure contains 51.4-45.6 wt pct Ni; it is ordered FeNi with the L1(0) superstructure. The CZ structure consists of two phases: a globular phase (ordered FeNi containing 50.9 wt pct Ni) and a surrounding honeycomb martensitic phase containing 11.7 wt pct Ni. The CT-2 was found in all of the iron meteorite groups studied and in the pallasites, but not in the mesosiderites. Based on the preliminary evidence, this region is believed to be ordered Fe3Ni. Possible mechanisms for the decomposition of taenite are discussed.

Reuter, K. B.↗

The Delplot kinetic method applied to systems with adsorbates: Hydrodeoxygenation of benzofuran on a bimetallic CoPd phosphide catalyst supported on KUSY

The Delplot method provides a means of analyzing conversion and selectivity data to derive a reaction sequence and has been applied so far to gas-phase or liquid-phase species. This study analyses the applicability of the method for catalytic reactions and considers for the first time adsorbed intermediates. The method is applied to a contact time study of the hydrodeoxygenation (HDO) of benzofuran at 0.5 MPa and 350 °C over a catalyst consisting of bimetallic CoPd phosphide supported on a potassium ion-exchanged ultra-stable Y (KUSY) zeolite. Both simple and detailed reaction networks were derived. The simple networks considered only gas-phase species and were modeled by a first-order sequence of steps, whereas the detailed network took into account adsorbed intermediates and was simulated using a rake mechanism. The networks were compared using F-statistics, which accounted for the differences in the number of fitting parameters. The detailed network gave a better fit to the experimental data because it represented a more realistic description of the transformation. A suggested sequence from the Delplot method was consistent with the simple network but not with the detailed network. The lack of applicability of the Delplot analysis to the network with adsorbates was linked to overly large equilibrium constants, a determination supported by Delplot fitting for a model sequence. Further, this study indicated the inapplicability of the Delplot method in determining reaction sequences that involve adsorbed species with large equilibrium constants for formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Self-organization in collisionless, high- β turbulence

The magnetohydrodynamic (MHD) equations, as a collisional fluid model that remains in local thermodynamic equilibrium (LTE), have long been used to describe turbulence in myriad space and astrophysical plasmas. Yet, the vast majority of these plasmas, from the solar wind to the intracluster medium (ICM) of galaxy clusters, are only weakly collisional at best, meaning that significant deviations from LTE are not only possible but common. Recent studies have demonstrated that the kinetic physics inherent to this weakly collisional regime can fundamentally transform the evolution of such plasmas across a wide range of scales. Here, we explore the consequences of pressure anisotropy and Larmor-scale instabilities for collisionless, β $\gg$ 1, turbulence, focusing on the role of a self-organizational effect known as ‘magneto-immutability’. We describe this self-organization analytically through a high-β, reduced ordering of the Chew–Goldberger–Low-MHD (CGL-MHD) equations, finding that it is a robust inertial-range effect that dynamically suppresses magnetic-field-strength fluctuations, anisotropic-pressure stresses and dissipation due to heat fluxes. As a result, the turbulent cascade of Alfvénic fluctuations continues below the putative viscous scale to form a robust, nearly conservative, MHD-like inertial range. These findings are confirmed numerically via Landau-fluid CGL-MHD turbulence simulations that employ a collisional closure to mimic the effects of microinstabilities. We find that microinstabilities occupy a small (~5%) volume-filling fraction of the plasma, even when the pressure anisotropy is driven strongly towards its instability thresholds. We discuss these results in the context of recent predictions for ion-vs-electron heating in low-luminosity accretion flows and observations implying suppressed viscosity in ICM turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhanced beam-beam modeling to include longitudinal variation during weak-strong simulation

Beam-beam interactions pose substantial challenges in the design and operation of circular colliders, significantly affecting their performance. In particular, the weak-strong simulation approach is pivotal for investigating single-particle dynamics during the collider design phase. This paper evaluates the limitations of existing models in weak-strong simulations, noting that while they accurately account for energy changes due to slingshot effects, they fail to incorporate longitudinal coordinate changes ( z variation). To address this gap, we introduce two novel transformations that enhance Hirata’s original framework by including both z variation and slingshot effect-induced energy changes. Through rigorous mathematical analysis and extensive weak-strong simulation studies, we validate the efficacy of these enhancements in achieving a more precise simulation of beam-beam interactions. Our results reveal that although z variation constitutes a higher-order effect and does not substantially affect the emittance growth rate within the specific design parameters of the Electron-Ion Collider, the refined model offers improved accuracy, particularly in scenarios involving the interaction between beam-beam effects and other random diffusion processes, as well as in simulations incorporating realistic lattice models. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Model of ramp compression of diamond from ab initio simulations

Ramp compression experiments characterize high-pressure states of matter at temperatures well below those present in shock compression. However, because temperature is typically not directly measured during ramp compression, it is uncertain how much heating occurs under these shock-free conditions. Here, we performed a series of ab initio simulations on carbon in order to match the density-stress measurements of Smith et al. [Smith et al., Nature (London) 511, 330 (2014)]. We considered isotropically as well as uniaxially compressed solid carbon in the diamond and BC8 phases, with and without defects, as well as liquid carbon. Our idealized model ascribes heating during ramp compression to an initially uniaxially compressed cell transforming isochorically into an isotropically (hydrostatic equivalent) compressed state having lower internal energy, hence higher temperature so as to conserve energy. Multiple such heating events can occur during a single ramp experiment, leading to higher temperatures than with isentropic compression. Comparison with experiments shows that heating alone does not explain the equation of state measurements on diamond, instead implying that a significant uniaxial stress component remains present at high compression. The temperature predictions of our ramp compression model remain to be verified with future laboratory measurements.

58 GEOSCIENCES↗