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 667 records · Page 37

Analysis of soot formation from aviation fuels in laminar counterflow flames

Combustion emissions from aviation contribute to the formation of condensation trail (contrail) that can lead to the formation of anthropogenic cirrus clouds. Ice particles that form contrails are observed to have a linear correlation with soot particle number density. Synthetic aviation fuels (SAFs) offer a promising route to mitigate the production of soot particles while also increasing energy security. Although studies have focused on combustion and spray behavior, the detailed investigation of soot formation processes for different jet fuels and their impact on models for computational fluid dynamics (CFD) applications is not well understood. Moreover, experimental measurements of soot for canonical flames using Synthetic aviation fuels (SAF) for model validation remain scarce. To address this, we use employed the Lawrence Livermore National Laboratory (LLNL) detailed soot model based on the discrete sectional method. Additionally, we develop two reduced chemical mechanisms for Jet-A and Alcohol-to-Jet (C1) that are suitable for turbulent flame simulations and couple them with the Hybrid Method of Moments (HMOM). The detailed and reduced model frameworks are validated against experimental measurements of soot volume fraction (ƒ ν ) from a counterflow burner experiment previously reported in the literature. Given the good agreement between modeling results and experimental measurements for the (1) spatial distribution of ƒ ν and (2) the non-linear variation of peak ƒ ν with strain rate, we further investigate the modeled sub-processes (nucleation, condensation, surface growth, and oxidation) using the LLNL model to analyze the assumptions in the reduced model framework. Furthermore, the results indicate a significant contribution from resonant radicals to the surface growth of soot particles, which are not accounted for in the current implementation of HMOM and could help reconcile soot predictions by the reduced model with observations.

Counterflow↗

Distributionally robust facility location problem under decision-dependent stochastic demand

While the traditional facility location problem considers exogenous demand, in some applications, locations of facilities could affect the willingness of customers to use certain types of services, e.g., carsharing, and therefore they also affect realizations of random demand. Moreover, a decision maker may not know the exact distribution of such endogenous demand and how it is affected by location choices. In this paper, we consider a distributionally robust facility location problem, in which we interpret the moments of stochastic demand as functions of facility-location decisions. We reformulate a two-stage decision-dependent distributionally robust optimization model as a monolithic formulation, and then derive exact mixed-integer linear programming reformulation as well as valid inequalities when the means and variances of demand are piecewise linear functions of location solutions. We conduct extensive computational studies, in which we compare our model with a decision-dependent deterministic model, as well as stochastic programming and distributionally robust models without the decision-dependent assumption. Here, the results show superior performance of our approach with remarkable improvement in profit and quality of service under various settings, in addition to computational speed-ups given by formulation enhancements. These results draw attention to the need of considering the impact of location decisions on customer demand within this strategic-level planning problem.

97 MATHEMATICS AND COMPUTING↗

Linear optimum trim solution for aircraft with three longitudinal control effectors

A simple and fast linear optimum method for trimming aircraft with a combination of three longitudinal aerodynamic or propulsive control effectors is developed, with application to preliminary design or in-flight scheduling. A Lagrangian formulation is used to determine the longitudinal load distribution which minimizes trim drag. Trim drag results obtained with the linear optimum method are similar to those obtained with a nonlinear optimum method, and significant trim drag reductions are found using either optimum method, in comparison with results obtained using various ad hoc methods.

Goodrich, Kenneth H.↗

Ion anisotropy instabilities in the magnetosheath

Recent observations in Earth's magnetosheath have delineated several different kinds of magnetic fluctuation spectra below the proton cyclotron frequency. This paper provides a theoretical interpretation for some of these observations describing solutions of the linear Vlasov dispersion equation for fully electromagnetic instabilities for particle distributions which model those observed in the magnetosheath. This model yields three growing modes: the proton cyclotron anisotropy the helium cyclotron anisotropy, and the mirror instabilities. The results show very good agreement with the observations of mirror-like and proton-cyclotron-like events. This agreement with observations implies that the transition between cyclotron and mirror fluctuation dominance is consistent with linear theory.

Gary, S. P.↗

Disturbance Dynamics in Transitional and Turbulent Flows

In order to expand the predictive capability of single-point turbulence closure models to account for the early-stage transition regime, a methodology for the formulation and calibration of model equations for the ensemble-averaged disturbance kinetic energy and energy dissipation rate is presented. First the decay of laminar disturbances and turbulence in mean shear-free flows is studied. In laminar flows, such disturbances are linear superpositions of modes governed by the Orr-Sommerfeld equation. In turbulent flows, disturbances are described through transport equations for representative mean quantities. The link between a description based on a deterministic evolution equation and a probability based mean transport equation is established. Because an uncertainty in initial conditions exists in the laminar as well as the turbulent regime, a probability distribution must be defined even in the laminar case. Using this probability distribution, it is shown that the exponential decay of the linear modes in the laminar regime can be related to a power law decay of both the (ensemble) mean disturbance kinetic energy and the dissipation rate. The evolution of these mean disturbance quantities is then described by transport equations similar to those for the corresponding turbulent decaying flow. Second, homogeneous shear flow, where disturbances can be described by rapid distortion theory (RDT), is studied. The relationship between RDT and linear stability theory is exploited in order to obtain a closed set of modeled equations. The linear disturbance equations are solved directly so that the numerical simulation yields a database from which the closure coefficients in the ensemble-averaged disturbance equations can be determined.

Grosch, Chester E.↗

Lower hybrid heating of ionospheric ions due to ion ring distributions in cusp

The stability of H(+) and H(++) ring distributions which have been observed downflowing into the cusp on the DE and S3-3 satellites is examined in the context of the feedback of those instabilities on plasma of ionospheric origin consisting of (H(+) and O(+). Lower hybrid waves are excited by the ring distributions in three distinct phases of wave-particle interaction: linear growth, trapping, and quasi-linear diffusion.The latter phase accounts for most particle heating. Including background O(+) and/or a He(++) ring introduces new modes not present in a pure H(+) plasma which play an important role in heating heavier ions. O(+) is heated significantly more by a He(++) ring than a H(+) ring of comparable energy density. It is suggested that lower hybrid waves generated by downflowing ion ring distributions play a role in energizing ion conics in the cusp.

Roth, I.↗

Amphiphile conformation impacts aggregate morphology and solution structure across multiple lengthscales

Although the self-assembly of amphiphiles is well-studied in aqueous solutions, much less is understood about the fundamental driving forces and structure property relationships in non-polar media. In recent work [Journal of Physical Chemistry B, 2020, 124, 10822.] the authors have studied a series of malonomide-based amphiphiles that are relevant to liquid-liquid extraction. That work demonstrated that aggregation is largely driven by local dipole-dipole interactions between molecules. Here, we build upon this observation to develop a more detailed understanding of how the balance of dipole-dipole interactions (controlled by conformation) and molecular architecture influences the morphology of the aggregates across lengthscales. Using constrained molecular dynamics about key degrees of freedom, we demonstrate that the conformation of N,N-dimethyl,N,N-dioctylhexylethoxy malonamide (DMDOHEMA) and N,N-dimethyl,N,N-dibutyltetradecyl malonamide (DMDBTDMA) has a significant impact upon self-association - where appropriate conformational sampling is essential. To quantify the aggregate morphology, several graph theoretic and persistent homology based properties are determined. The former examines the patterns of intermolecular interactions within clusters, while the latter examines the 3-dimensional spatial distribution. Based upon these analyses, we find that the morphology of aggregates, particularly at higher malonamide concentration, depends on a balance of dipole alignment and alkyl tail sterics. Dipole alignment encourages linear patterns of the intermolecular interactions within aggregates, while the alkyl tail steric interactions between the malonamide result in noticeably less linear aggregates for DMDOHEMA than DMDBTDMA. This is reflected in the spatial distribution, where more holes or voids exist between extractants within the DMDOHEMA that distribute within the solution in more of a "swiss cheese" arrangement as opposed to the more filamentous distribution of DMBDTDMA. This study links conformation and molecular structure to the morphology of amphiphile assemblies, and serves as a basis for ongoing study of multicomponent amphiphile solutions with polar and other solutes, and how these impact aggregation phenomena.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fusion techniques using distributed Kalman filtering for detecting changes in systems

A comparison is made of the performances of two detection strategies that are based on different data fusion techniques. The strategies detect changes in a linear system. One detection strategy involves combining the estimates and error covariance matrices of distributed Kalman filters, generating a residual from the used estimates, comparing this residual to a threshold, and making a decision. The other detection strategy involves a distributed decision process in which estimates from distributed Kalman filters are used to generate distributed residuals which are compared locally to a threshold. Local decisions are made and these decisions are then fused into a global decision. The performances of each of these detection schemes are compared, and it is concluded that better performance is achieved when local decisions are made and then fused into a global decision.

Belcastro, Celeste M.↗

Frequency distributions and correlations of solar X-ray flare parameters

Frequency distributions of flare parameters are determined from over 12,000 solar flares. The flare duration, the peak counting rate, the peak hard X-ray flux, the total energy in electrons, and the peak energy flux in electrons are among the parameters studied. Linear regression fits, as well as the slopes of the frequency distributions, are used to determine the correlations between these parameters. The relationship between the variations of the frequency distributions and the solar activity cycle is also investigated. Theoretical models for the frequency distribution of flare parameters are dependent on the probability of flaring and the temporal evolution of the flare energy build-up. The results of this study are consistent with stochastic flaring and exponential energy build-up. The average build-up time constant is found to be 0.5 times the mean time between flares.

Crosby, Norma B.↗

Desalination for a circular water economy

Today's water systems are enabled by ample fresh water sources, low-cost centralized treatment, and facile wastewater disposal. Climatic change, aging infrastructure, and source water contamination have exposed the vulnerabilities of this linear water paradigm. While seawater desalination enables coastal communities to augment their supply, more broadly securing water systems for municipal, industrial, and agricultural water users will require distributed desalination and fit-for-purpose reuse of nontraditional water sources. Our linear water economy must evolve into a resilient circular water economy, where water is continuously reused and “contaminants” become the feedstocks for other economically valuable processes. Additionally, technology innovation is needed to deliver autonomous, precise, resilient, intensified, modular, and electrified desalination systems that reduce the cost, improve the performance, and enhance the resilience of nontraditional water reuse systems. Meanwhile, strong federal leadership and coordination is needed to accelerate desalination research, promote information gathering efforts to direct technology development, and create an expanded role for non-profit organizations in knowledge dissemination.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Calcined Polyethyleneimine Coated Optical Fibers for Distributed pH Monitoring at High Pressures and Temperatures

In the oil and gas, CO2 sequestration, H2 subsurface storage, and geothermal energy sectors, subsurface pH measurements are critical for monitoring the geochemical conditions and estimating potential corrosion rates of wellbore systems. Real-time pH measurements in these conditions are vital for detecting and predicting corrosion deterioration of wellbore components that may jeopardize the safety and continued operation of wellbore systems. Building off of earlier TiO2 based pH sensors and the known pH sensitivity of amine-based polymers, a coating based on the secondary amine polymer polyethyleneimine (PEI) was developed. The sensor coating was treated with a high temperature (500 °C) calcination procedure in air to convert it into a more stable oxidized coating capable of withstanding hot aqueous solutions without dissolving while retaining linear pH sensitivity from pH values between 2 and 11. The sensor performance was measured using optical transmission measurements in solutions of various pHs and using optical backscatter reflectometry for distributed pH sensing demonstration in wellbore-relevant pressures (up to 1000 psi) and temperatures (80 °C). A calibration curve with strong differentiation between acidic and alkaline pH was developed for both transmission-based and distributed pH measurements, using fixed wavelength transmission and integrated linear amplitude of backscattered light for distributed measurements.

Shumski, Alexander↗

On the effect of scalar flux weighing of linearly anisotropic scattering matrices in few-group transport calculations

The majority of the codes available for homogenized group constant generation for deterministic transport calculations apply the approximation of scalar flux weighting during energy group condensation of higher-order anisotropic scattering matrices. In this paper, we discuss the effect of scalar flux weighting of the linearly anisotropic scattering matrices in the frame of S P{sub 3} and S{sub 12} calculations performed for a two-dimensional VVER-440 reactor benchmark. To compare group constants generated for two neutron energy groups, an infinite pin cell was homogenized with Serpent 2 and ERANOS ECCO. Serpent 2 applies scalar flux, while ERANOS ECCO performs current weighting of the linearly anisotropic scattering matrices during energy group condensation. For analyzing the effect of the various weighting options, three simple reactor models were built assuming different core sizes using standard rectangular assemblies with 15*15 fuel pins. Diffusion, SP{sub 3} and S{sub 12} calculations were performed for the 3 models using group constants generated with Serpent 2 and ERANOS ECCO. The effect of scalar flux weighting of linearly anisotropic scattering matrices is shown by comparing the decrease in reactivity due to the decreased reactor size, as well as the assembly power distribution to reference results obtained with Serpent 2 Monte Carlo calculations. Neglecting higher than linearly anisotropic scattering and indirect application of diffusion coefficients in higher-order transport calculations is advised if angular flux-moment spectra weighted higher-order scattering matrices cannot be generated. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Spline-based distributed system identification with application to large space antennas

A parameter and state estimation technique for distributed models is demonstrated through the solution of a problem generic to large space antenna system identification. Assuming the position of the reflective surface of the maypole (hoop/column) antenna to be approximated by the static two-dimensional, stretched-membrane partial differential equation with variable-stiffness coefficient functions, a spline-based approximation procedure is described that estimates the shape and stiffness functions from data set observations. For given stiffness functions, the Galerkin projection with linear spline-based functions is applied to project the distributed problem onto a finite-dimensional subspace wherein algebraic equations exist for determining a static shape (state) prediction. The stiffness functions are then parameterized by cubic splines and the parameters estimated by an output error technique. Numerical results are presented for data descriptive of a 100-m-diameter maypole antenna.

Banks, H. T.↗

Applied Joint-Space Torque and Stiffness Control of Tendon-Driven Fingers

Existing tendon-driven fingers have applied force control through independent tension controllers on each tendon, i.e. in the tendon-space. The coupled kinematics of the tendons, however, cause such controllers to exhibit a transient coupling in their response. This problem can be resolved by alternatively framing the controllers in the joint-space of the manipulator. This work presents a joint-space torque control law that demonstrates both a decoupled and significantly faster response than an equivalent tendon-space formulation. The law also demonstrates greater speed and robustness than comparable PI controllers. In addition, a tension distribution algorithm is presented here to allocate forces from the joints to the tendons. It allocates the tensions so that they satisfy both an upper and lower bound, and it does so without requiring linear programming or open-ended iterations. The control law and tension distribution algorithm are implemented on the robotic hand of Robonaut-2.

Abdallah, Muhammad E.↗

A Large-Particle Monte Carlo Code for Simulating Non-Linear High-Energy Processes Near Compact Objects

High-energy radiation processes in compact cosmic objects are often expected to have a strongly non-linear behavior. Such behavior is shown, for example, by electron-positron pair cascades and the time evolution of relativistic proton distributions in dense radiation fields. Three independent techniques have been developed to simulate these non-linear problems: the kinetic equation approach; the phase-space density (PSD) Monte Carlo method; and the large-particle (LP) Monte Carlo method. In this paper, we present the latest version of the LP method and compare it with the other methods. The efficiency of the method in treating geometrically complex problems is illustrated by showing results of simulations of 1D, 2D and 3D systems. The method is shown to be powerful enough to treat non-spherical geometries, including such effects as bulk motion of the background plasma, reflection of radiation from cold matter, and anisotropic distributions of radiating particles. It can therefore be applied to simulate high-energy processes in such astrophysical systems as accretion discs with coronae, relativistic jets, pulsar magnetospheres and gamma-ray bursts.

Stern, Boris E.↗

Pressure safety approach for PIP-II cryogenic distribution system and cryomodules

The Proton Improvement Plan-II (PIP-II) is a superconducting linear accelerator being built at Fermilab that will provide 800 MeV proton beam for neutrino production. The linac consists of a total of twenty-three (23) cryomodules of five (5) different types. Cooling is required at 2K, 5K and 40K. The Cryogenic Distribution System (CDS) consists of a Distribution Valve Box, ~285 m of cryogenic transfer line, modular Bayonet Cans to interface with cryomodules, and a Turnaround Can. The cryogenic system must provide protection from over-pressure by sizing pressure relief devices for all volumes and process line circuits. The cryomodule cavity circuits have dual pressure ratings, 4.1 bara when cold and 2.05 bara when warm (T>80K). Worst case relieving cases will be identified. The methods for determining heat flux will be presented. For the relieving occurring in the linac tunnel, flow must vent to outside to avoid an oxygen deficiency hazard. Also, we will present vacuum vesse l relief sizing to protect the cryogenic distribution system vacuum shells from over pressure during an internal line rupture. The project is funded by US DOE Offices of Science, High Energy Physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Pressure safety approach for PIP-II cryogenic distribution system and cryomodules

The Proton Improvement Plan-II (PIP-II) is a superconducting linear accelerator being built at Fermilab that will provide 800 MeV proton beam for neutrino production. The linac consists of a total of twenty-three (23) cryomodules of five (5) different types. Cooling is required at 2K, 5K and 40K. The Cryogenic Distribution System (CDS) consists of a Distribution Valve Box, ~285 m of cryogenic transfer line, modular Bayonet Cans to interface with cryomodules, and a Turnaround Can. The cryogenic system must provide protection from over-pressure by sizing pressure relief devices for all volumes and process line circuits. The cryomodule cavity circuits have dual pressure ratings, 4.1 bara when cold and 2.05 bara when warm (T>80K). Worst case relieving cases will be identified. The methods for determining heat flux will be presented. For the relieving occurring in the linac tunnel, flow must vent to outside to avoid an oxygen deficiency hazard. Also, we will present vacuum vessel relief sizing to protect the cryogenic distribution system vacuum shells from over pressure during an internal line rupture. The project is funded by US DOE Offices of Science, High Energy Physics.

43 PARTICLE ACCELERATORS↗

Re-Analysis of DARHT Axis 2 S4 Magnet Sweep Measurements

The Dual-Axis Radiographic Hydrodynamic Test (DARHT) facility provides flash radiography capabilities using two electron Linear Induction Accelerators (LIA’s). Understanding fundamental properties of the electron beams is essential to optimizing spot size and dose for radiography. This work describes distribution measurements of a single-kicked electron beam (16.5MeV, 1.7kA, 60ns) at a location downstream at imaging station C on DARHT Axis 2. The beam at station C is focused by the S4 solenoid located upstream. The beam distribution is measured by imaging optical transition radiation (OTR) from electrons striking a 51-µm thick titanium foil. This data and analysis were originally published by Ekdahl. The results here are a re-analysis of the asymmetric halo emittance contribution as compared to the emittance resulting from the FWHM of the distribution. This is accomplished using xtr2 fits to the data which include non-linear magnetic field effects in S4. Also examined are the calculated S4 spherical aberration contributions to the emittance using xtr.

43 PARTICLE ACCELERATORS↗