Search NASA⌕ Search

SEARCH · Search NASA

Results for “fully linear model”

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 523 records · Page 29

Pure quark and gluon observables in collinear drop

We construct a class of pure quark and gluon observables by using the collinear drop grooming technique. The construction is based on linear combinations of multiple cumulative distributions of the jet mass in collinear drop, whose specific weights are fully predicted perturbatively. This yields observables which obtain their values purely from quarks (or purely from gluons) in a wide region of phase space. We demonstrate this by showing that these observables are effective in two phase space regions, one dominated by perturbative resummation and one dominated by nonperturbative effects. The nonperturbative effects are included using shape functions which only appear as a common factor in the linear combinations constructed. We test this construction using a numerical analysis with next-to-leading logarithmic resummation and various shape function models, as well as analyzing these observables with PYTHIA and VINCIA. Choices for the collinear drop parameters are optimized for experimental use.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Research In Nonlinear Flight Control for Tiltrotor Aircraft Operating in the Terminal Area

The research during the first year of the effort focused on the implementation of the recently developed combination of neural net work adaptive control and feedback linearization. At the core of this research is the comprehensive simulation code Generic Tiltrotor Simulator (GTRS) of the XV-15 tilt rotor aircraft. For this research the GTRS code has been ported to a Fortran environment for use on PC. The emphasis of the research is on terminal area approach procedures, including conversion from aircraft to helicopter configuration. This report focuses on the longitudinal control which is the more challenging case for augmentation. Therefore, an attitude command attitude hold (ACAH) control augmentation is considered which is typically used for the pitch channel during approach procedures. To evaluate the performance of the neural network adaptive control architecture it was necessary to develop a set of low order pilot models capable of performing such tasks as, follow desired altitude profiles, follow desired speed profiles, operate on both sides of powercurve, convert, including flaps as well as mastangle changes, operate with different stability and control augmentation system (SCAS) modes. The pilot models are divided in two sets, one for the backside of the powercurve and one for the frontside. These two sets are linearly blended with speed. The mastangle is also scheduled with speed. Different aspects of the proposed architecture for the neural network (NNW) augmented model inversion were also demonstrated. The demonstration involved implementation of a NNW architecture using linearized models from GTRS, including rotor states, to represent the XV-15 at various operating points. The dynamics used for the model inversion were based on the XV-15 operating at 30 Kts, with residualized rotor dynamics, and not including cross coupling between translational and rotational states. The neural network demonstrated ACAH control under various circumstances. Future efforts will include the implementation into the Fortran environment of GTRS, including pilot modeling and NNW augmentation for the lateral channels. These efforts should lead to the development of architectures that will provide for fully automated approach, using similar strategies.

Calise, A. J.↗

A linear shock cell model for jets of arbitrary exit geometry

The shock cell structures of single supersonic non-ideally expanded jets with arbitrary exit geometry are studied. Both vortex sheets and realistic mean profiles are considered for the jet shear layer. The boundary element method is used to predict the shock spacing and screech tones in a vortex sheet model of a single jet. This formulation enables the calculations to be performed only on the vortex sheet. This permits the efficient and convenient study of complicated jet geometries. Results are given for circular, elliptic and rectangular jets and the results are compared with analysis and experiment. The agreement between the predictions and measurements is very good but depends on the assumptions made to predict the geometry of the fully expanded jet. A finite diffference technique is used to examine the effect of finite mixing layer thickness for a single jet. The finite thickness of the mixing layer is found to decrease the shock spacing by approximately 20 percent over the length of the jet potential core.

Morris, P. J.↗

Pairwise interaction extended point-particle model for a random array of monodisperse spheres

This study introduces a new point-particle force model that attempts to account for the hydrodynamic influence of the neighbouring particles in an Eulerian–Lagrangian simulation. In previous point-particle models the force on a particle depends only on Reynolds number and mean volume fraction. Thus, as long as the mean local volume fraction is the same, the force on different particles will be estimated to be the same, even though the precise arrangement of neighbours can be vastly different. From direct numerical simulation (DNS) it has been observed that in a random arrangement of spheres that were distributed with uniform probability, the particle-to-particle variation in force can be as large as the mean drag. Since the Reynolds number and mean volume fraction of all the particles within the array are the same, the standard models fail to account for the significant particle-to-particle force variation within the random array. Here, we develop a model which can compute the drag and lateral forces on each particle by accounting for the precise location of a few surrounding neighbours. A pairwise interaction is assumed where the perturbation flow induced by each neighbour is considered separately, then the effects of all neighbours are linearly superposed to obtain the total perturbation. Faxén correction is used to quantify the force perturbation due to the presence of the neighbours. The single neighbour perturbations are mapped in the vicinity of a reference sphere and stored as libraries. We test the pairwise interaction extended point-particle (PIEP) model for random arrays at two different volume fractions of $\unicode[STIX]{x1D719}=0.1$ and 0.21 and Reynolds numbers in the range $16.5\leqslant Re\leqslant 170$ . The PIEP model predictions are compared against drag and lift forces obtained from the fully resolved DNS simulations performed using the immersed boundary method. Although not perfect, we observe the PIEP model prediction to correlate much better with the DNS results than the classical mean drag model prediction.

Mechanics↗

Extended Lagrangian Born–Oppenheimer molecular dynamics for orbital-free density-functional theory and polarizable charge equilibration models

We report extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) is formulated for orbital-free Hohenberg–Kohn density-functional theory and for charge equilibration and polarizable force-field models that can be derived from the same orbital-free framework. The purpose is to introduce the most recent features of orbital-based XL-BOMD to molecular dynamics simulations based on charge equilibration and polarizable force-field models. These features include a metric tensor generalization of the extended harmonic potential, preconditioners, and the ability to use only a single Coulomb summation to determine the fully equilibrated charges and the interatomic forces in each time step for the shadow Born–Oppenheimer potential energy surface. The orbital-free formulation has a charge-dependent, short-range energy term that is separate from long-range Coulomb interactions. This enables local parameterizations of the short-range energy term, while the long-range electrostatic interactions can be treated separately. The theory is illustrated for molecular dynamics simulations of an atomistic system described by a charge equilibration model with periodic boundary conditions. The system of linear equations that determines the equilibrated charges and the forces is diagonal, and only a single Ewald summation is needed in each time step. The simulations exhibit the same features in accuracy, convergence, and stability as are expected from orbital-based XL-BOMD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An mRNA vaccine encoding the Ebola virus glycoprotein induces high neutralizing antibody titers and provides strong protection against lethal infections in mouse models

Ebola virus (EBOV) is the causative agent of Ebola disease (EBOD), a viral hemorrhagic fever with a notably high case fatality rate. Current treatments for EBOD are limited to monoclonal antibodies or two licensed viral vector vaccines, a recombinant vesicular stomatitis virus (rVSV)-vectored vaccine or an adenovirus and modified vaccinia Ankara regimen. However, comparisons of protection, efficacy, and durability with alternative nucleotide platforms remain understudied. Here, we evaluated the immunogenicity of an mRNA vaccine expressing the EBOV glycoprotein (GP) in parallel with rVSV- and DNA-based vaccine platforms. The mRNA EBOV-GP vaccine, formulated in lipid nanoparticles, elicited significantly higher levels of total IgG and neutralizing antibody titers compared to the rVSV-EBOV-GP vaccine. Linear antibody epitope analysis indicated a preference for targeting the mucin-like domain in EBOV-GP1 following rVSV-based vaccination, while the mRNA platform distinctly targeted the internal fusion loop of EBOV-GP2. After characterizing the immunogenicity of the mRNA vaccine, two models of EBOD were used to demonstrate its protective efficacy: a surrogate rVSV-based challenge model of EBOD using type-I interferon deficient C57BL/6 mice and infection of BALB/c mice with authentic mouse-adapted EBOV. In both studies, the EBOV mRNA vaccine fully protected the mouse cohorts against morbidity and mortality. Additionally, the EBOV mRNA vaccine produced greater neutralizing antibody titers compared to the DNA EBOV-GP vaccine. These results suggest that an mRNA vaccine expressing EBOV-GP can induce robust, functional humoral responses that are protective against EBOD, warranting further development as an alternative to, or as part of a vaccine strategy including, viral vectored vaccines.

DNA vaccines↗

A novel conditional generative model for efficient ensemble forecasts of state variables in large-scale geological carbon storage

Integrating monitoring data to efficiently update reservoir pressure and CO 2 plume distribution forecasts presents a significant challenge in geological carbon storage (GCS) applications. Inverse modeling techniques are commonly used to fuse observational data and refine reservoir model parameters, thereby improving state variable forecasts. However, these techniques often rely on linear or Gaussian assumptions, which can limit their effectiveness in accurately predicting state variables. Moreover, simulating large-scale three-dimensional (3D) GCS problems is computationally expensive, making iterative runs in inverse problems prohibitive. To address these challenges, we propose a conditional generative model utilizing the score-based diffusion method for real-time 3D pressure and saturation field distribution predictions. Our approach involves solving the score function with a mini-batch-based Monte Carlo estimator to generate labeled data. This data is subsequently employed to train a fully connected neural network, enabling it to learn the conditional sample generator within a supervised learning framework. This method enables the rapid generation of a large ensemble of predictions, facilitating comprehensive uncertainty quantification of state variables. Here we applied our method to forecast the dynamic 3D distributions of pressure and saturation fields over a 30-year injection period. The statistical assessment with low root mean square error (RMSE) values demonstrates that our method can accurately predict the spatiotemporal distributions of both pressure and saturation fields. Moreover, the developed conditional generative model shows high computational efficiency by generating 100 ensemble forecasts of 3D state variables in less than 10 min. The consistency between ensemble averages and ground truth values further illustrates the model’s capability to capture state variable dynamics during the CO 2 plume injection process. Notably, the ground truth values fall within the ensemble forecasts, indicating that our uncertainty quantification effectively captures variability and potential noise in the observations. Thus, the developed conditional generative model proves to be a more efficient, accurate, and practical tool for GCS applications, facilitating timely risk analysis and informed decision-making.

58 GEOSCIENCES↗

Dislocation dynamics in hexagonal close-packed crystals

Extensions of the dislocation dynamics methodology necessary to enable accurate simulations of crystal plasticity in hexagonal close-packed (HCP) metals are presented. They concern the introduction of dislocation motion in HCP crystals through linear and non-linear mobility laws, as well as the treatment of composite dislocation physics. Formation, stability and dissociation of and other dislocations with large Burgers vectors defined as composite dislocations are examined and a new topological operation is proposed to enable their dissociation. Furthermore, the results of our simulations suggest that composite dislocations are omnipresent and may play important roles both in specific dislocation mechanisms and in bulk crystal plasticity in HCP materials. While fully microscopic, our bulk DD simulations provide wealth of data that can be used to develop and parameterize constitutive models of crystal plasticity at the mesoscale.

36 MATERIALS SCIENCE↗

Morphological Expression of the Coherence and Relative Phase of Optical Inputs to the Photoelectrodeposition of Nanopatterned Se–Te Films

Highly anisotropic and ordered nanoscale lamellar morphologies can be spontaneously generated over macroscopic areas, without the use of a photomask or any templating agents, via the photoelectrodeposition of Se–Te alloy films. To form such structures, the light source can be a single, linearly polarized light source that need not necessarily be highly coherent. In this work, the variation in the morphologies produced by this deposition process was evaluated in response to differences in the coherence and relative phase between multiple simultaneous linearly polarized illumination inputs. Specifically, the morphologies of photoelectrodeposits were evaluated when two tandem same-wavelength sources with discrete linear polarizations, both either mutually incoherent or mutually coherent (with defined phase differences), were used. Additionally, morphologies were simulated via computer modeling of the interfacial light scattering and absorption during the photoelectrochemical growth process. The morphologies that were generated using two coherent, in-phase sources were equivalent to those generated using only a single source. In contrast, the use of two coherent, out-of-phase sources produced a range of morphological patterns. For small out-of-phase addition of orthogonal polarization components, lamellar-type patterns were observed. When fully out-of-phase orthogonal sources (circular polarization) were used, an isotropic, mesh-type pattern was instead generated, similar to that observed when unpolarized illumination was utilized. In intermediate cases, anisotropic lamellar-type patterns were superimposed on the isotropic mesh-type patterns, and the relative height between the two structures scaled with the amount of out-of-phase addition of the orthogonal polarization components. Similar results were obtained when two incoherent sources were utilized. In every case, the long axis of the lamellar-type morphology component aligned parallel to the intensity-weighted average polarization orientation. The observations consistently agreed with computer simulations, indicating that the observed morphologies were fully determined by the nature of the illumination utilized during the growth process. Here, the collective data thus indicated that the photoelectrodeposition process exhibits sensitivity toward the coherency, relative phase, and polarization orientations of all optical inputs and that this sensitivity is physically expressed in the morphology of the deposit.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

A dual-input nonlinear system analysis of autonomic modulation of heart rate

Linear analyses of fluctuations in heart rate and other hemodynamic variables have been used to elucidate cardiovascular regulatory mechanisms. The role of nonlinear contributions to fluctuations in hemodynamic variables has not been fully explored. This paper presents a nonlinear system analysis of the effect of fluctuations in instantaneous lung volume (ILV) and arterial blood pressure (ABP) on heart rate (HR) fluctuations. To successfully employ a nonlinear analysis based on the Laguerre expansion technique (LET), we introduce an efficient procedure for broadening the spectral content of the ILV and ABP inputs to the model by adding white noise. Results from computer simulations demonstrate the effectiveness of broadening the spectral band of input signals to obtain consistent and stable kernel estimates with the use of the LET. Without broadening the band of the ILV and ABP inputs, the LET did not provide stable kernel estimates. Moreover, we extend the LET to the case of multiple inputs in order to accommodate the analysis of the combined effect of ILV and ABP effect on heart rate. Analyzes of data based on the second-order Volterra-Wiener model reveal an important contribution of the second-order kernels to the description of the effect of lung volume and arterial blood pressure on heart rate. Furthermore, physiological effects of the autonomic blocking agents propranolol and atropine on changes in the first- and second-order kernels are also discussed.

NASA Discipline Regulatory Physiology↗

The Relaxation Matrix for Symmetric Tops with Inversion Symmetry. I. Effects of Line Coupling on Self-Broadened v (sub 1) and Pure Rotational Bands of NH3

The Robert-Bonamy formalism has been commonly used to calculate half-widths and shifts of spectral lines for decades. This formalism is based on several approximations. Among them, two have not been fully addressed: the isolated line approximation and the neglect of coupling between the translational and internal motions. Recently, we have shown that the isolated line approximation is not necessary in developing semi-classical line shape theories. Based on this progress, we have been able to develop a new formalism that enables not only to reduce uncertainties on calculated half-widths and shifts, but also to model line mixing effects on spectra starting from the knowledge of the intermolecular potential. In our previous studies, the new formalism had been applied to linear and asymmetric-top molecules. In the present study, the method has been extended to symmetric-top molecules with inversion symmetry. As expected, the inversion splitting induces a complete failure of the isolated line approximation. We have calculated the complex relaxation matrices of selfbroadened NH3. The half-widths and shifts in the ν1 and the pure rotational bands are reported in the present paper. When compared with measurements, the calculated half-widths match the experimental data very well, since the inapplicable isolated line approximation has been removed. With respect to the shifts, only qualitative results are obtained and discussed. Calculated off-diagonal elements of the relaxation matrix and a comparison with the observed line mixing effects are reported in the companion paper (Paper II).

Trajectory models↗

Augmented Method to Improve Thermal Data for the Figure Drift Thermal Distortion Predictions of the JWST OTIS Cryogenic Vacuum Test

The JWST Optical Telescope Element (OTE) assembly is the largest optically stable infrared-optimized telescope currently being manufactured and assembled, and is scheduled for launch in 2018. The JWST OTE, including the 18 segment primary mirror, secondary mirror, and the Aft Optics Subsystem (AOS) are designed to be passively cooled and operate near 45K. These optical elements are supported by a complex composite backplane structure. As a part of the structural distortion model validation efforts, a series of tests are planned during the cryogenic vacuum test of the fully integrated flight hardware at NASA JSC Chamber A. The successful ends to the thermal-distortion phases are heavily dependent on the accurate temperature knowledge of the OTE structural members. However, the current temperature sensor allocations during the cryo-vac test may not have sufficient fidelity to provide accurate knowledge of the temperature distributions within the composite structure. A method based on an inverse distance relationship among the sensors and thermal model nodes was developed to improve the thermal data provided for the nanometer scale WaveFront Error (WFE) predictions. The Linear Distance Weighted Interpolation (LDWI) method was developed to augment the thermal model predictions based on the sparse sensor information. This paper will encompass the development of the LDWI method using the test data from the earlier pathfinder cryo-vac tests, and the results of the notional and as tested WFE predictions from the structural finite element model cases to characterize the accuracies of this LDWI method.

Thermal Modeling↗

Evaluation of NorESM-OC (versions 1 and 1.2), the ocean carbon-cycle stand-alone configuration of the Norwegian Earth System Model (NorESM1)

Idealised and hindcast simulations performed with the stand-alone ocean carbon-cycle configuration of the Norwegian Earth System Model (NorESM-OC) are described and evaluated. We present simulation results of three different model configurations (two different model versions at different grid resolutions) using two different atmospheric forcing data sets. Model version NorESM-OC1 corresponds to the version that is included in the NorESM-ME1 fully coupled model, which participated in CMIP5. The main update between NorESM-OC1 and NorESM-OC1.2 is the addition of two new options for the treatment of sinking particles. We find that using a constant sinking speed, which has been the standard in NorESM's ocean carbon cycle module HAMOCC (HAMburg Ocean Carbon Cycle model), does not transport enough particulate organic carbon (POC) into the deep ocean below approximately 2000 m depth. The two newly implemented parameterisations, a particle aggregation scheme with prognostic sinking speed, and a simpler scheme that uses a linear increase in the sinking speed with depth, provide better agreement with observed POC fluxes. Additionally, reduced deep ocean biases of oxygen and remineralised phosphate indicate a better performance of the new parameterisations. For model version 1.2, a re-tuning of the ecosystem parameterisation has been performed, which (i) reduces previously too high primary production at high latitudes, (ii) consequently improves model results for surface nutrients, and (iii) reduces alkalinity and dissolved inorganic carbon biases at low latitudes. We use hindcast simulations with prescribed observed and constant (pre-industrial) atmospheric CO 2 concentrations to derive the past and contemporary ocean carbon sink. As a result, for the period 1990–1999 we find an average ocean carbon uptake ranging from 2.01 to 2.58 Pg C yr -1 depending on model version, grid resolution, and atmospheric forcing data set.

54 ENVIRONMENTAL SCIENCES↗

Structural behavior of scientific balloons - Finite element simulation and verification

An off-the-shelf nonlinear finite element code was used to analyze fully inflated scientific balloons. The thin balloon film was modeled by shell bending elements. Numerical difficulties caused by insignificant bending stiffness terms were overcome by introducing some artificial bending stiffness. This approximation is justified by the fact that in thin shells with nonzero Gaussian curvature the membrane solution component is essentially independent of the bending solution component. Perturbation of the coveraged solution by increasing the bending stiffness by a full decade verified this assertion. This analytical approach was experimentally verified. As a result of this verification process it was discovered that the generally accepted linearly visco-elastic model for polyethylene film is inappropriate for a significant planar (as opposed to uniaxial) stress state. A linear elastic model presents a good approximation for planar stress states.

Schur, Willi W.↗

Particle simulations of mode conversion between slow mode and fast mode in lower hybrid range of frequencies

The propagation and mode conversion of lower hybrid waves in an inhomogeneous plasma are investigated by using the nonlinear δf algorithm in a two-dimensional particle-in-cell simulation code based on the gyrokinetic electron and fully kinetic ion (GeFi) scheme [Lin et al., Plasma Phys. Controlled Fusion 47, 657 (2005)]. The characteristics of the simulated waves, such as wavelength, frequency, phase, and group velocities, agree well with the linear theoretical analysis. It is shown that a significant reflection component emerges in the conversion process between the slow mode and the fast mode when the scale length of the density variation is comparable to the local wavelength. The dependences of the reflection coefficient on the scale length of the density variation are compared with the results based on the linear full wave model for cold plasmas. It is indicated that the mode conversion for the waves with a frequency of 2.45 GHz (ω ∼ 3ω{sub LH}, where ω{sub LH} represents the lower hybrid resonance) and within Tokamak relevant amplitudes can be well described in the linear scheme. As the frequency decreases, the modification due to the nonlinear term becomes important. For the low-frequency waves (ω ∼ 1.3ω{sub LH}), the generations of the high harmonic modes and sidebands through nonlinear mode-mode coupling provide new power channels and thus could reduce the reflection significantly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling jet and outflow feedback during star cluster formation

Powerful jets and outflows are launched from the protostellar disks around newborn stars. These outflows carry enough mass and momentum to transform the structure of their parent molecular cloud and to potentially control star formation itself. Despite their importance, we have not been able to fully quantify the impact of jets and outflows during the formation of a star cluster. The main problem lies in limited computing power. We would have to resolve the magnetic jet-launching mechanism close to the protostar and at the same time follow the evolution of a parsec-size cloud for a million years. Current computer power and codes fall orders of magnitude short of achieving this. In order to overcome this problem, we implement a subgrid-scale (SGS) model for launching jets and outflows, which demonstrably converges and reproduces the mass, linear and angular momentum transfer, and the speed of real jets, with ∼1000 times lower resolution than would be required without the SGS model. We apply the new SGS model to turbulent, magnetized star cluster formation and show that jets and outflows (1) eject about one-fourth of their parent molecular clump in high-speed jets, quickly reaching distances of more than a parsec, (2) reduce the star formation rate by about a factor of two, and (3) lead to the formation of ∼1.5 times as many stars compared to the no-outflow case. Most importantly, we find that jets and outflows reduce the average star mass by a factor of ∼ three and may thus be essential for understanding the characteristic mass of the stellar initial mass function.

79 ASTRONOMY AND ASTROPHYSICS↗

A Lagrangian dynamic subgrid-scale model turbulence

A new formulation of the dynamic subgrid-scale model is tested in which the error associated with the Germano identity is minimized over flow pathlines rather than over directions of statistical homogeneity. This procedure allows the application of the dynamic model with averaging to flows in complex geometries that do not possess homogeneous directions. The characteristic Lagrangian time scale over which the averaging is performed is chosen such that the model is purely dissipative, guaranteeing numerical stability when coupled with the Smagorinsky model. The formulation is tested successfully in forced and decaying isotropic turbulence and in fully developed and transitional channel flow. In homogeneous flows, the results are similar to those of the volume-averaged dynamic model, while in channel flow, the predictions are superior to those of the plane-averaged dynamic model. The relationship between the averaged terms in the model and vortical structures (worms) that appear in the LES is investigated. Computational overhead is kept small (about 10 percent above the CPU requirements of the volume or plane-averaged dynamic model) by using an approximate scheme to advance the Lagrangian tracking through first-order Euler time integration and linear interpolation in space.

Meneveau, C.↗

EXPERIMENTAL STUDIES OF NONLINEAR INTEGRABLE OPTICS

State-of-the-art accelerators at energy and intensity frontiers require increasingly bright and powerful particle beams. In conventional linear lattices, intense beams suffer from collective instabilities, resulting in beam losses and maximum beam intensity limits. This thesis presents experimental studies of a novel lattice design concept, the nonlinear integrable optics (NIO), aimed at enhancing beam stability limits with little to no beam dynamics degradation. Single-particle beam dynamics measurements of two NIO devices, the quasi-integrable octupole system and the fully integrable Danilov-Nagaitsev system, were carried out at the purpose-built Integrable Optics Test Accelerator (IOTA) at Fermilab. Their simulation, hardware design, and commissioning process are presented. Extensive model and analysis algorithm development and benchmarking is described. Electron beam data from two scientific runs is analyzed, yielding frequency and phase space dynamics consisten t with m odels. These results demonstrate viability and advantages of the NIO design, providing the groundwork for proton studies in the strong space-charge regime and future integrable accelerators.

Kuklev, Nikita↗