Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linearized distribution”

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,027 records · Page 57

Effect of a Nonplanar Melt-Solid Interface On Lateral Compositional Distribution During Unidirectional Solidification of a Binary Alloy with a Constant Growth Velocity V: Theory - Pt. 1

Infrared detected materials, such as Hg(1-x)Cd(x)Te, Hg(1-x)Zn(x)Te have energy gaps almost linearly proportional to their composition. Due to the wide separation of liquidus and solidus curves of their phase diagram, there are compositional segregation in both of the axial and radial directions of these crystals grown in the Bridgman system unidirectionally with constant growth rate. It is important to understand the mechanisms, which affect lateral segregation such that large radially uniform composition crystal can be produced. Following Coriel, etc's treatment, we have developed a theory to study the effect of a curved melt-solid interface shape on lateral composition distribution. The model is considered to be a cylindrical system with azimuthal symmetry and a curved melt-solid interface shape which can be expressed as a linear combination of a series of Bessell's functions. The results show that melt-solid interface shape has a dominant effect on the lateral composition distribution of these systems. For small values of beta, the solute concentration at the melt-solid interface scales linearly with interface shape with a proportional constant of the produce of beta and (1 -k), where beta = VR/D, with V as growth velocity, R as the sample radius, D as the diffusion constant and k as the distribution constant. A detailed theory will be presented. A computer code has been developed and simulations have been performed and compared with experimental results. These will be published in another paper.

Wang, Jai-Ching↗

Direct Numerical Simulation of Transitional Multicomponent-Species Gaseous and Multicomponent-Liquid Drop-Laden Mixing

A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.

tubulent multicompenent-species↗

An MILP-Based Distributed Energy Management for Coordination of Networked Microgrids

An MILP-based distributed energy management for the coordination of networked microgrids is proposed in this paper. Multiple microgrids and the utility grid are coordinated through iteratively adjusted price signals. Based on the price signals received, the microgrid controllers (MCs) and distribution management system (DMS) update their schedules separately. Then, the price signals are updated according to the generation–load mismatch and distributed to MCs and DMS for the next iteration. The iteration continues until the generation–load mismatch is small enough, i.e., the generation and load are balanced under agreed price signals. Through the proposed distributed energy management, various microgrids and the utility grid with different economic, resilient, emission and socio-economic objectives are coordinated with generation–load balance guaranteed and the microgrid customers’ privacy preserved. In particular, a piecewise linearization technique is employed to approximate the augmented Lagrange term in the alternating direction method of multipliers (ADMM) algorithm. Thus, the subproblems are transformed into mixed integer linear programming (MILP) problems and efficiently solved by open-source MILP solvers, which would accelerate the adoption and deployment of microgrids and promote clean energy. The proposed MILP-based distributed energy management is demonstrated through various case studies on a networked microgrids test system with three microgrids.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ∼1/Np with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly non-equilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Linear stability analysis of the Vlasov-Poisson equations in high density plasmas in the presence of crossed fields and density gradients

The equations for the single-particle orbits in a nonneutral high density plasma in the presence of inhomogeneous crossed fields are obtained. Using these orbits, the linearized Vlasov equation is solved as an expansion in the orbital radii in the presence of inhomogeneities and density gradients. A model distribution function is introduced whose cold-fluid limit is exactly the same as that used in many previous studies of the cold-fluid equations. This model function is used to reduce the linearized Vlasov-Poisson equations to a second-order ordinary differential equation for the linearized electrostatic potential whose eigenvalue is the perturbation frequency.

Kaup, D. J.↗

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ~1/$\sqrt{N_p}$ with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly nonequilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Study of Ion Generation by a ps-laser for External Injection into Electron Beam Ion Source

We investigated ion generation from Al, Ti, Cu, Nb and Ta target elements by a picosecond-laser in the range of 10 11 – 10 13 W/cm 2 power densities at the target surface. A ps-laser with 1.27 mJ maximum energy within a 8 ps pulse and repetition rate up to 400 Hz has been used to generate a laser-ablated plasma. Dependencies of ion current vs time, total charge of registered ions as well as ion kinetic energy distributions are characterized using a Faraday cup. Significant difference in ion current dynamics between first, second and following shots onto the same target spot was found for all five target elements. The total charge of ions registered by the Faraday cup increase linearly with increasing laser pulse energy and are almost independent on the target element and number of shots onto the same target spot for all five target elements studied. Target craters have been examined using a commercially available microscope (Nikon Eclipse LV150) and crater cross section areas were found to be in the range of 0.0076 – 0.0172 mm 2 within the first five shots onto the same target spot for different target elements. A linear dependence of crater cross section square on laser pulse energy has been observed for the Ta target. Ion energy distributions for all target elements have a plateau or a slight peak in the energy range of 10 - 100 eV and long low-intensity energy tail up to tens of keV.

43 PARTICLE ACCELERATORS↗

Time Series Dimension Reduction for Surrogate Models of Port Scanning Cyber Emulations

Surrogate model development is a key resource in the scientific modeling community for providing computational expedience when simulating complex systems without loss of great fidelity. The initial step to development of a surrogate model is identification of the primary governing components of the system. Principal component analysis (PCA) is a widely used data science technique that provides inspection of such driving factors, when the objective for modeling is to capture the greatest sources of variance inherent to a dataset. Although an efficient linear dimension reduction tool, PCA makes the fundamental assumption that the data is continuous and normally distributed. Thus, it provides ideal performance when these conditions are met. In the case for which cyber emulations provide realizations of a port scanning scenario, the data to be modeled follows a discrete time series function comprised of monotonically increasing piece-wise constant steps. The sources of variance are related to the timing and magnitude of these steps. Therefore, we consider using XPCA, an extension to PCA for continuous and discrete random variates. This report provides the documentation of the trade-offs between the PCA and XPCA linear dimension reduction algorithms, for the intended purpose to identify key components of greatest variance in our time series data. These components will ultimately provide the basis for future surrogate models of port scanning cyber emulations.

97 MATHEMATICS AND COMPUTING↗

Condition-Based Maintenance of a Circulating Water System of a Canadian Nuclear Power Plant using Machine Learning and Statistical Tools

Canada Deuterium Uranium pressurized-heavy-water reactors (PHWR) are a type of nuclear power plant that generate clean and reliable energy. The scope of this work is to automate data analysis methodologies to inform a condition-based maintenance strategy of a circulating water system (CWS) of a PHWR. The multiunit CWS provides a continuous supply of water to cool steam condensers, even during transient scenarios, thereby improving the thermal efficiency. This work aims to develop a machine learning (ML) based approach to detect anomalies in heterogeneous data of a CWS in a PHWR to help inform a predictive maintenance strategy. The heterogeneous data include textual and numeric time series data for a PHWR. Natural-language-processing (NLP)-based models are used to analyze textual data contained in work orders and operator logs and an event-timeseries correlation detection method is applied to assist anomalies diagnoses for CWS. An ML model Robust Linear Model (RLM) is also used to remove the seasonal variations in the system variable distributions based on distributions of environmental variables. A machine learning model, Density-Based Spatial Clustering of Applications with Noise (DBSCAN), trained on both original data and data without any seasonal variations will then be used to detect if an anomaly exists. Thus, by moving to an automated methodology to detect, classify, and forecast anomalies, the maintenance strategy would be based on component condition instead of a time-based schedule.

97 - MATHEMATICS AND COMPUTING↗

Wave and ion evolution downstream of quasi-perpendicular bow shocks

Distribution functions of ions heated in quasi-perpendicular bow shocks have a large perpendicular temperature anisotropy that provides free energy for the growth of Alfven ion cyclotron (AIC) waves and mirror waves. Both types of waves have been observed in the Earth's magnetosheath downstream of quasi-perpendicular shocks. We use a two-dimensional hybrid simulations to give a self-consistent description of the evolution of the wave spectra downstream of quasi-perpendicular shocks. Both mirror and AIC waves are identified in the simulated magnetosheath. They are generated at or near the shock front and convected away from it by the sheath plasma. Near the shock, the waves have a broad spectrum, but downstream of the shock, shorter-wavelength modes are heavily damped and only longer-wavelength modes persist. The characteristics of these surviving modes can be predicted with reasonable accuracy by linear kinetic theory appropriate for downstream conditions. We also follow the evolution of the ion distribution function. The shocked ions that provide the free energy for wave growth have a two-component distribution function. The halo is initially gyrophase-bunched and extremely anisotropic. Within a relatively short distance downstream of the shock (of the order of 10 ion inertial lengths), wave-particle interactions remove these features from the halo and reduce the anisotropy of the distribution to near-threshold levels for the mirror and AIC instabilities. A similar evolution has been observed for ions at the Earth's bow shock.

Mckean, M. E.↗

Characterizing the Roles of Biogeochemical Cycling and Ocean Circulation in Regulating Marine Copper Distributions

Copper (Cu) is a key micronutrient for marine phytoplankton. Its oceanic biogeochemical cycle has elicited considerable attention due to dissolved Cu exhibiting a unique linear profile with depth. Several processes have been proposed for explaining this behavior. Here, in this study, we characterize the relationships between the observed Cu, PO 4 3– , and Si on a global scale. We find that the depth profiles of Cu resemble those of Si more than of PO 4 3– in the global ocean. To understand their relationships, we couple the biogeochemical and internal circulation processes in a model, fitting optimal Cu:PO 4 3– uptake ratios and remineralization length-scales to replicate the marine Cu distributions. The modeling results suggest that Cu uptake needs in the Southern Ocean are substantially higher than those in other oceanic regions. In addition, our modeling results indicate a deep Cu remineralization in the global ocean. We offer an alternative mechanism that relies on biogeochemical cycling and internal circulation to produce the linear depth profiles of dissolved Cu. Our results suggest that diatoms are likely the major phytoplankton dominating oceanic Cu cycling.

54 ENVIRONMENTAL SCIENCES↗

Generalized Jeans' Escape of Pick-Up Ions in Quasi-Linear Relaxation

Jeans escape is a well-validated formulation of upper atmospheric escape that we have generalized to estimate plasma escape from ionospheres. It involves the computation of the parts of particle velocity space that are unbound by the gravitational potential at the exobase, followed by a calculation of the flux carried by such unbound particles as they escape from the potential well. To generalize this approach for ions, we superposed an electrostatic ambipolar potential and a centrifugal potential, for motions across and along a divergent magnetic field. We then considered how the presence of superthermal electrons, produced by precipitating auroral primary electrons, controls the ambipolar potential. We also showed that the centrifugal potential plays a small role in controlling the mass escape flux from the terrestrial ionosphere. We then applied the transverse ion velocity distribution produced when ions, picked up by supersonic (i.e., auroral) ionospheric convection, relax via quasi-linear diffusion, as estimated for cometary comas [1]. The results provide a theoretical basis for observed ion escape response to electromagnetic and kinetic energy sources. They also suggest that super-sonic but sub-Alfvenic flow, with ion pick-up, is a unique and important regime of ion-neutral coupling, in which plasma wave-particle interactions are driven by ion-neutral collisions at densities for which the collision frequency falls near or below the gyro-frequency. As another possible illustration of this process, the heliopause ribbon discovered by the IBEX mission involves interactions between the solar wind ions and the interstellar neutral gas, in a regime that may be analogous [2].

Moore, T. E.↗

X-ray linear dichroic tomography of crystallographic and topological defects

Abstract The functionality of materials is determined by their composition 1–4 and microstructure, that is, the distribution and orientation of crystalline grains, grain boundaries and the defects within them 5,6 . Until now, characterization techniques that map the distribution of grains, their orientation and the presence of defects have been limited to surface investigations, to spatial resolutions of a few hundred nanometres or to systems of thickness around 100 nm, thus requiring destructive sample preparation for measurements and preventing the study of system-representative volumes or the investigation of materials under operational conditions 7–15 . Here we present X-ray linear dichroic orientation tomography (XL-DOT), a quantitative, non-invasive technique that allows for an intragranular and intergranular characterization of extended polycrystalline and non-crystalline 16 materials in three dimensions. We present the detailed characterization of a polycrystalline sample of vanadium pentoxide (V 2 O 5 ), a key catalyst in the production of sulfuric acid 17 . We determine the nanoscale composition, microstructure and crystal orientation throughout the polycrystalline sample with 73 nm spatial resolution. We identify and characterize grains, as well as twist, tilt and twin grain boundaries. We further observe the creation and annihilation of topological defects promoted by the presence of volume crystallographic defects. The non-destructive and spectroscopic nature of our method opens the door to operando combined chemical and microstructural investigations 11,18 of functional materials, including energy, mechanical and quantum materials.

Science & Technology - Other Topics↗

Dynamic of Langmuir and Ion-Sound Waves in Type 3 Solar Radio Sources

The evolution of Langmuir and ion-sound waves in type 3 sources is investigated, incorporating linear growth, linear damping, and nonlinear electrostatic decay. Improved estimates are obtained for the wavenumber range of growing waves and the nonlinear coupling coefficient for the decay process. The resulting prediction for the electrostatic decay threshold is consistent with the observed high-field cutoff in the Langmuir field distribution. It is shown that the conditions in the solar wind do not allow a steady state to be attained; rather, bursty linear and nonlinear interactions take place, consistent with the highly inhomogeneous and impulsive waves actually observed. Nonlinear growth is found to be fast enough to saturate the growth of the parent Langmuir waves in the available interaction time. The resulting levels of product Langmuir and ion-sound waves are estimated theoretically and shown to be consistent with in situ ISEE 3 observations of type 3 events at 1 AU. Nonlinear interactions slave the growth and decay of product sound waves to that of the product Langmuir waves. The resulting probability distribution of ion-sound field strengths is predicted to have a flat tail extending to a high-field cutoff. This prediction is consistent with statistics derived here from ISEE 3 observations. Agreement is also found between the frequencies of the observed waves and predictions for the product S waves. The competing processes of nonlinear wave collapse and quasilinear relaxation are discussed, and it is concluded that neither is responsible for the saturation of Langmuir growth. When wave and beam inhomogeneities are accounted for, arguments from quasi-linear relaxation yield an upper bound on the Langmuir fields that is too high to be relevant. Nor are the criteria for direct wave collapse of the beam-driven waves met, consistent with earlier simulation results that imply that this process is not responsible for saturation of the beam instability. Indeed, even if the highest observed Langmuir fields are assumed to he part of a long-wavelength 'condensate' produced via electrostatic decay, they still fall short of the relevant requirements for wave collapse. The most stringent requirement for collapse is that collapsing wave packets not be disrupted by ambient density fluctuations in the solar wind. Fields of several mV m(exp -1) extending over several hundred km would be needed to satisfy this requirement; at 1 AU such fields are rare at best.

Robinson, P. A.↗

Transient vibration and product formation of photoexcited CS 2 measured by time-resolved x-ray scattering

Here, we have observed details of the internal motion and dissociation channels in photoexcited carbon disulfide (CS 2 ) using time-resolved x-ray scattering (TRXS). Photoexcitation of gas-phase CS 2 with a 200 nm laser pulse launches oscillatory bending and stretching motion, leading to dissociation of atomic sulfur in under a picosecond. During the first 300 fs following excitation, we observe significant changes in the vibrational frequency as well as some dissociation of the C–S bond, leading to atomic sulfur in the both 1 D and 3 P states. Beyond 1400 fs, the dissociation is consistent with primarily 3 P atomic sulfur dissociation. This channel-resolved measurement of the dissociation time is based on our analysis of the time-windowed dissociation radial velocity distribution, which is measured using the temporal Fourier transform of the TRXS data aided by a Hough transform that extracts the slopes of linear features in an image. The relative strength of the two dissociation channels reflects both their branching ratio and differences in the spread of their dissociation times. Measuring the time-resolved dissociation radial velocity distribution aids the resolution of discrepancies between models for dissociation proposed by prior photoelectron spectroscopy work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Magnetic fields in star-forming regions - Observations

We review the observational aspects of magnetic fields in dense, star-forming regions. First we discuss ways to observe the field. These include direct methods, which consist of the measurement of both linear and circular polarization of spectral line and continuum radiation; and indirect methods, consisting of the angular distribution of H2O masers on the sky and the measurement of ambipolar diffusion. Next we discuss selected observational results, focusing on detailed discussions of a small number of points rather than a generalized discussion that covers the waterfront. We discuss the Orion/BN-KL region in detail, both on the small and large scales. Next we discuss the derivation of the complete magnetic vector, including both the systematic and fluctuating component, from a large sample of Zeeman and linear polarization measurements for the L204 dark cloud. We examine the virial theorem as it applies to dark clouds in general and one dark cloud, Barnard 1, in particular. We critically discuss the numerous claims for alignment of cloud structural features with the plane-of-the-sky component of the magnetic field, and find that many of these have not been definitively established.

Heiles, Carl↗

Removing Malmquist bias from linear regressions

Malmquist bias is present in all astronomical surveys where sources are observed above an apparent brightness threshold. Those sources which can be detected at progressively larger distances are progressively more limited to the intrinsically luminous portion of the true distribution. This bias does not distort any of the measurements, but distorts the sample composition. We have developed the first treatment to correct for Malmquist bias in linear regressions of astronomical data. A demonstration of the corrected linear regression that is computed in four steps is presented.

Verter, Frances↗

Data-Conforming Data-Driven Control: Avoiding Premature Generalizations Beyond Data

Data-driven and adaptive control approaches face the problem of introducing sudden distributional shifts beyond the distribution of data encountered during learning. Therefore, they are prone to invalidating the very assumptions used in their own construction. This is due to the linearity of the underlying system, inherently assumed and formulated in most data-driven control approaches, which may falsely generalize the behavior of the system beyond the behavior experienced in the data. This article seeks to mitigate these problems by enforcing consistency of the newly designed closed-loop systems with data and slowing down any distributional shifts in the joint state-input space. This is achieved through incorporating affine regularization terms and linear matrix inequality constraints to data-driven approaches, resulting in convex semi-definite programs that can be efficiently solved by standard software packages. We discuss the optimality conditions of these programs and then conclude this article with a numerical example that further highlights the problem of premature generalization beyond data and shows the effectiveness of our proposed approaches in enhancing the safety of data-driven control methods.

97 MATHEMATICS AND COMPUTING↗