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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Modelling galaxy cluster triaxiality in stacked cluster weak lensing analyses

Counts of galaxy clusters offer a high-precision probe of cosmology, but control of systematic errors will determine the accuracy of this measurement. Using Buzzard simulations, we quantify one such systematic, the triaxiality distribution of clusters identified with the redMaPPer optical cluster finding algorithm, which was used in the Dark Energy Survey Year-1 (DES Y1) cluster cosmology analysis. We test whether redMaPPer selection biases the clusters’ shape and orientation and find that it only biases orientation, preferentially selecting clusters with their major axes oriented along the line of sight. Modelling the richness–mass relation as log-linear, we find that the log-richness amplitude ln (A) is boosted from the lowest to highest orientation bin with a significance of 14σ, while the orientation dependence of the richness-mass slope and intrinsic scatter is minimal. We also find that the weak lensing shear-profile ratios of cluster-associated dark haloes in different orientation bins resemble a ‘bottleneck’ shape that can be quantified with a Cauchy function. We test the correlation of orientation with two other leading systematics in cluster cosmology – miscentering and projection – and find a null correlation. The resulting mass bias predicted from our templates confirms the DES Y1 finding that triaxiality is a leading source of bias in cluster cosmology. However, the richness-dependence of the bias confirms that triaxiality does not fully resolve the tension at low-richness between DES Y1 cluster cosmology and other probes. Our model can be used for quantifying the impact of triaxiality bias on cosmological constraints for upcoming weak lensing surveys of galaxy clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

OPTICAL I-BAND LINEAR POLARIMETRY OF THE MAGNETAR 4U 0142+61 WITH SUBARU

Magnetars are known to have optical and/or infrared (IR) emission, but the origin of the emission is not well understood. In order to fully study their emission properties, we have carried out for the first time optical linear polarimetry of the magnetar 4U 0142+61, which has been determined from different observations to have a complicated broadband spectrum over optical and IR wavelengths. From our I-band imaging polarimetric observation, conducted with the 8.2-m Subaru telescope, we determine the degree of linear polarization to be P = 1.0 ± 3.4%, or P ≤ 5.6% (90% confidence level). Considering models that were suggested for optical emission from magnetars, we discuss the implications of our result. The upper limit measurement indicates that, differing from radio pulsars, magnetars probably would not have strongly polarized optical emission if the emission arises from their magnetosphere as suggested.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1977

Failure detection by pilots during automatic landing - Models and experiments

A model is proposed to describe the pilot as a monitor of automatic landing systems. The failures treated are equivalent to the addition of a dynamic change in the mean of the observation process. The failure detection model of the pilot consists of two stages: a linear estimator (Kalman filter) and a decision mechanism based on sequential analysis. The filter equations are derived from a simplified version of the linearized dynamics of the airplane and the control loop. The perceptual observation noise is modified to include the effects of allocation of attention among the several instruments. The final result is a simple model consisting of a high-pass filter to produce the observation residuals and a decision function which is a pure integration of the residuals minus a bias term. The dynamics of a Boeing 707 were used to simulate the fully coupled final approach in a fixed-base simulator. Observers monitored the approaches and detected the failures; their performance was compared with the predictions of the model.

Gai, E. G.↗
NASA NTRS2000

High-Lift Flight Tunnel - Phase II Report

The High-Lift Flight Tunnel (HiLiFT) concept is a revolutionary approach to aerodynamic ground testing. This concept utilizes magnetic levitation and linear motors to propel an aerodynamic model through a tube containing a quiescent test medium. This medium (nitrogen) is cryogenic and pressurized to achieve full flight Reynolds numbers higher than any existing ground test facility world-wide for the range of 0.05 to 0.50 Mach. The results of the Phase II study provide excellent assurance that the HiLiFT concept will provide a valuable low-speed, high Reynolds number ground test facility. The design studies concluded that the HiLiFT facility is feasible to build and operate and the analytical studies revealed no insurmountable difficulties to realizing a practical high Reynolds number ground test facility. It was determined that a national HiLiFT facility, including development, would cost approximately $400M and could be operational by 2013 if fully funded. Study participants included National Aeronautics and Space Administration Langley Research Center as the Program Manager and MSE Technology Applications, Inc., (MSE) of Butte, Montana as the prime contractor and study integrator. MSE#s subcontractors included the University of Texas at Arlington for aerodynamic analyses and the Argonne National Laboratory for magnetic levitation and linear motor technology support.

Lofftus, David↗
DOE OSTI2016

Scalable implicit incompressible resistive MHD with stabilized FE and fully-coupled Newton–Krylov-AMG

Here, we discuss that the computational solution of the governing balance equations for mass, momentum, heat transfer and magnetic induction for resistive magnetohydrodynamics (MHD) systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena, as well as the significant range of time- and length-scales that the interactions of these physical mechanisms produce. This paper explores the development of a scalable, fully-implicit stabilized unstructured finite element (FE) capability for 3D incompressible resistive MHD. The discussion considers the development of a stabilized FE formulation in context of the variational multiscale (VMS) method, and describes the scalable implicit time integration and direct-to-steady-state solution capability. The nonlinear solver strategy employs Newton–Krylov methods, which are preconditioned using fully-coupled algebraic multilevel preconditioners. These preconditioners are shown to enable a robust, scalable and efficient solution approach for the large-scale sparse linear systems generated by the Newton linearization. Verification results demonstrate the expected order-of-accuracy for the stabilized FE discretization. The approach is tested on a variety of prototype problems, that include MHD duct flows, an unstable hydromagnetic Kelvin–Helmholtz shear layer, and a 3D island coalescence problem used to model magnetic reconnection. Initial results that explore the scaling of the solution methods are also presented on up to 128K processors for problems with up to 1.8B unknowns on a CrayXK7.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

ANALYTICAL THEORY FOR THE INITIAL MASS FUNCTION. III. TIME DEPENDENCE AND STAR FORMATION RATE

The present paper extends our previous theory of the stellar initial mass function (IMF) by including time dependence and by including the impact of the magnetic field. The predicted mass spectra are similar to the time-independent ones with slightly shallower slopes at large masses and peak locations shifted toward smaller masses by a factor of a few. Assuming that star-forming clumps follow Larson-type relations, we obtain core mass functions in good agreement with the observationally derived IMF, in particular, when taking into account the thermodynamics of the gas. The time-dependent theory directly yields an analytical expression for the star formation rate (SFR) at cloud scales. The SFR values agree well with the observational determinations of various Galactic molecular clouds. Furthermore, we show that the SFR does not simply depend linearly on density, as is sometimes claimed in the literature, but also depends strongly on the clump mass/size, which yields the observed scatter. We stress, however, that any SFR theory depends, explicitly or implicitly, on very uncertain assumptions like clump boundaries or the mass of the most massive stars that can form in a given clump, making the final determinations uncertain by a factor of a few. Finally, we derive a fully time dependent model for the IMF by considering a clump, or a distribution of clumps accreting at a constant rate and thus whose physical properties evolve with time. In spite of its simplicity, this model reproduces reasonably well various features observed in numerical simulations of converging flows. Based on this general theory, we present a paradigm for star formation and the IMF.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Kinetic study of strong shock waves in fully ionized plasmas

Fully kinetic simulations are used in this work to study strong shock waves in fully ionized plasmas for charge state 𝑍 ranging from 1 to 18. The ion stress tensor amplitude agrees well with classical predictions within the compression layer; however, its spatial structure is noticeably shaped by nonlocal transport effects. The electron stress tensor remains significantly smaller than that of the ions for all values of 𝑍 considered. Ion overheating in the compression layer is observed to scale approximately linearly with 𝑍. Both findings diverge from hydrodynamic models based on classical electron stress formulations. In conclusion, the electrons exhibit strong nonlocal behavior in the compression layer, with Knudsen numbers approaching unity for 𝑍 ≥ 6.

Physics - Plasma physics↗
DOE OSTI2022

Automatic Calibration of a Geomechanical Model from Sparse Data for Estimating Stress in Deep Geological Formations

Summary In this study, we demonstrate geomechanical modeling with fully automatic parameter calibration to estimate the full geomechanical stress fields of a prospective US carbon dioxide (CO2) storage site, based on sparse measurement data. The goal is to compute full stress tensor field estimates (principal stresses and orientations) that are maximally compatible with observations within the constraints of the model assumptions, thereby extending pointwise, incomplete partial stress measurement to a simulated full formation stress field, as well as a rough assessment of the associated error. We use the Perch site, located in Otsego County, Michigan, USA, as our case study. The input data consist of partial stress tensor information inferred from in-situ borehole tests, geophysical well logs, and processing of seismic data. A static earth model (SEM) of the site was developed, and geomechanical simulation functionality of the open-source MATLAB Reservoir Simulation Toolbox (MRST) was used to model the stress field. Adjoint-based nonlinear optimization was used to adjust boundary conditions and material properties to calibrate simulated results of observations. Results were interpreted through a Bayesian framework. The focus of this paper is to demonstrate how the fully automatic calibration procedure works and discuss the results obtained; it does not attempt a detailed analysis of the stress field in the context of the proposed CO2 storage initiatives. Our work is part of a larger effort to noninvasively determine in-situ stresses in deep formations considered for CO2 storage. Guided by previously published research on geomechanical model calibration, our work presents a novel calibration approach supporting a potentially large number of linear or nonlinear calibration parameters to produce results optimally agreeing with available measurements and thus extend partial pointwise estimates to full tensor fields compatible with the physics of the site.

Engineering↗
NASA NTRS1994

Effects of Tangential Edge Constraints on the Postbuckling Behavior of Flat and Curved Panels Subjected to Thermal and Mechanical Loads

A parametric study of the effects of tangential edge constraints on the postbuckling response of flat and shallow curved panels subjected to thermal and mechanical loads is presented. The mechanical loads investigated are uniform compressive edge loads and transverse lateral pressure. The temperature fields considered are associated with spatially nonuniform heating over the panels, and a linear through-the-thickness temperature gradient. The structural model is based on a higher-order transverse-shear-deformation theory of shallow shells that incorporates the effects of geometric nonlinearities, initial geometric imperfections, and tangential edge motion constraints. Results are presented for three-layer sandwich panels made from transversely isotropic materials. Simply supported panels are considered in which the tangential motion of the unloaded edges is either unrestrained, partially restrained, or fully restrained. These results focus on the effects of the tangential edge restraint on the postbuckling response. The results of this study indicate that tangentially restraining the edges of a curved panel can make the panel insensitive to initial geometric imperfections in some cases.

Lin, W.↗
NASA NTRS1997

Nonlinear Behavior of a Typical Airfoil Section with Control Surface Freeplay: A Numerical and Experimental Study

A three degree-of-freedom aeroelastic typical section with control surface freeplay is modeled theoretically as a system of piecewise linear state-space models. The system response is determined by time marching of the governing equations using a standard Runge-Kutta algorithm in conjunction with Henon's method for integrating a system of equations to a prescribed surface of phase space section. Henon's method is used to locate the "switching points" accurately and efficiently as the system moves from one linear region into another. An experimental model which closely approximates the three degree-of-freedom, typical section in two-dimensional, incompressible flow has been created to validate the theoretical model. Consideration is given to modeling realistically the structural damping present in the experimental system. The effect of the freeplay on the system response is examined numerically and experimentally. The development of the state-space model offers a low-order, computationally efficient means of modeling fully the freeplay nonlinearity and may offer advantages in future research which will investigate the effects of freeplay on the control of flutter in the typical section.

Conner, M. D.↗
DOE OSTI2016

Fourier Analysis of the Multilevel QD Method with {alpha}-Approximation For Time-Dependent Radiative Transfer Problems

Thermal radiation transfer (TRT) problems describe interaction of photon radiation with matter. They are defined by the time-dependent radiative transfer (RT) equation for the specific intensity I coupled with the energy balance (EB) equation. This class of problems is characterised by high dimensionality, multiple scales and strong nonlinearity. In this paper we analyze the multilevel quasi-diffusion (QD) method for solving TRT problems. This is a nonlinear method of moments with exact closures. The QD method enables one to reduce dimensionality of the original radiative transfer problem. It is formulated by a nonlinear system of equations consisting of (i) the high-order RT (HORT) equation for the intensity of photon radiation and (ii) low-order equations for the angular and energy moments of the intensity. The nonlinear closures are formulated by means of stable factors. The low-order problem brings the radiative transfer model to the level of multi-physics equations. The multilevel QD (MLQD) method possesses attractive features for efficient multi-physics coupling. To calculate the high-order RT solution we apply the {alpha}-approximation which assumes that the intensity changes exponentially in time over a time step. The {alpha}-approximation leads to the steady-state RT equation with a modified opacity and enables one to avoid storing the high-dimensional solution from the previous time step. The rate of change is evaluated by the solution of the time-dependent low-order QD (LOQD) equations for moments. We note that the HORT equation in the {alpha}-approximation contains an extra nonlinear term. It affects the iterative behaviour of the MLQD method. The obtained computational results indicated that the rate of convergence is slower compared to the rate when the time-dependent HORT is used to calculate the intensity. In this paper we perform an iterative stability analysis of the MLQD method in a discretized form and analyze convergence of outer (transport) iterations on a time step. We consider multigroup TRT problems in 1D slab geometry. The fully implicit scheme is used to discretize equations in time. The HORT equation is approximated with the method of step characteristics. A finite volume method is applied for spatial discretization of multigroup LOQD equations. The grey (one-group) LOQD equations are algebraically consistent with the multigroup LOQD equations. To perform Fourier analysis a model TRT problem in an infinite homogeneous medium is considered. The system of nonlinear equations of the MLQD method is linearized near an equilibrium solution at some given temperature. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

SU-E-T-19: A Comparison of the Dosimetric Effects of Brass Mesh and Superflab Boluses

Purpose: We compared the dosimetric effects of brass mesh (Whiting and Davis, Attleboro Falls, MA) and Superflab (CNMC, Nashville, TN) boluses to verify equivalence between the two. Brass mesh bolus may be a convenient alternative to traditional bolus as it better conforms to the skin surface. Methods: Measurements were taken using a 6MV 10×10cm field produced by an Infinity linear accelerator (Elekta, Stockholm, Sweden) in a solid water phantom using a parallel plate ion chamber (Model 96035, Keithley Instruments, Cleveland, OH). Measurements compared 0.5cm and 1cm of Superflab to one to six layers of brass bolus mesh. Measurements were performed at depths from 0cm (‘skin’) to 10cm. Oblique beams were not studied. Results: Four layers of brass mesh were equivalent to 0.5cm Superflab within 5% at all depths. Six layers of brass compares most favorably with 0.5cm Superflab, with doses at all depths within 3%. Six layers of brass were not fully equivalent to 1cm Superflab, although the agreement was within 5% for depths greater than 3mm. Surface dose was 12% lower than 1cm Superflab. Surface dose can be up to 34% different between Superflab and brass mesh, but is less than 5% different with 4–6 layers of brass when compared to 0.5cm Superflab. To achieve surface dose agreement better than 5% compared to 1cm Superflab, more than 6 layers would be needed. The attenuation at depth was a maximum of 0.17cm per layer. Conclusion: Between four and six layers of brass mesh can be equivalent to 0.5cm Superflab, depending on the level of agreement desired. Equivalence within 5% at all depths to 1cm Superflab was not achieved even with six layers. This data agrees with measurements taken by Utsunomiya et al. (2010). More point measurements at shallower depths should be taken prior to clinical implementation of brass bolus mesh.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Level set topology optimization of structural problems with interface cohesion

Summary This paper presents a finite element topology optimization framework for the design of two‐phase structural systems considering contact and cohesion phenomena along the interface. The geometry of the material interface is described by an explicit level set method, and the structural response is predicted by the extended finite element method. In this work, the interface condition is described by a bilinear cohesive zone model on the basis of the traction‐separation constitutive relation. The non‐penetration condition in the presence of compressive interface forces is enforced by a stabilized Lagrange multiplier method. The mechanical model assumes a linear elastic isotropic material, infinitesimal strain theory, and a quasi‐static response. The optimization problem is solved by a nonlinear programming method, and the design sensitivities are computed by the adjoint method. The performance of the presented method is evaluated by 2D and 3D numerical examples. The results obtained from topology optimization reveal distinct design characteristics for the various interface phenomena considered. In addition, 3D examples demonstrate optimal geometries that cannot be fully captured by reduced dimensionality. The optimization framework presented is limited to two‐phase structural systems where the material interface is coincident in the undeformed configuration, and to structural responses that remain valid considering small strain kinematics. Copyright © 2017 John Wiley & Sons, Ltd.

Behrou, Reza↗
DOE OSTI2026

Helicon edge heating in Proto-MPEX: Azimuthally resolved RF modeling

The Material Plasma Exposure eXperiment (MPEX) is a linear device under construction at ORNL. MPEX employs a variety of heating systems, including a helicon antenna, which is the subject of this work. Lacking a resonant absorption mechanism, as tokamaks have, modeling the RF power absorption in such antennas requires detailed consideration of cold collisional physics. Previous modeling of the heating in this antenna was done exclusively in 2D, assuming axisymmetry and a single azimuthal mode number m = 1⁠. In this study, we introduce a fully 3D model. By including the full 3D physics, with its self-consistent geometry-imposed azimuthal spectrum, we are able to reproduce the experimentally observed transition from edge heating to core heating.

Tierens, Wouter [Oak Ridge National Laboratory (OR↗
DOE OSTI2016

A parametrisation of modified gravity on nonlinear cosmological scales

Viable modifications of gravity on cosmological scales predominantly rely on screening mechanisms to recover Einstein's Theory of General Relativity in the Solar System, where it has been well tested. A parametrisation of the effects of such modifications in the spherical collapse model is presented here for the use of modelling the modified nonlinear cosmological structure. The formalism allows an embedding of the different screening mechanisms operating in scalar-tensor theories through large values of the gravitational potential or its first or second derivatives as well as of linear suppression effects or more general transitions between modified and Einstein gravity limits. Each screening or suppression mechanism is parametrised by a time, mass, and environment dependent screening scale, an effective modified gravitational coupling in the fully unscreened limit that can be matched to linear theory, the exponent of a power-law radial profile of the screened coupling, determined by derivatives, symmetries, and potentials in the scalar field equation, and an interpolation rate between the screened and unscreened limits. Along with generalised perturbative methods, the parametrisation may be used to formulate a nonlinear extension to the linear parametrised post-Friedmannian framework to enable generalised tests of gravity with the wealth of observations from the nonlinear cosmological regime.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

Cryogenic star-tracking telescope for Gravity Probe B

This paper describes the design, development and preliminary testing of the cryogenic star-tracking telescope used as an optical reference for the gyroscopes in the Gravity Probe B Relativity Gyroscope experiment. The telescope is operated at 1.8 K; it is fabricated entirely from fused quartz components held together by optical contacting; it has a physical length of 14 in., a focal length of 150 in. and an aperture of 5.6 in. Readout is by two photomultiplier chopper-detector assemblies at ambient satellite temperature. When fully operational, the telescope may be expected to have a precision approaching 0.1 milliarcsec over a linear range of 70 + or 70 milliarcsec. Its projected noise performance corresponds to an angular resolution of 1 milliarcsec in 1 Hz bandwidth. The paper includes a theoretical analysis, a description of the design and fabrication of a laboratory version of the telescope, a discussion of techniques of optical contacting, an account of vibration tests on a separate mass model of the telescope, a description of the artificial star developed for optical tests, and an account of preliminary experimental results.

Everitt, C. W. F.↗
DOE OSTI2024

Sound propagation in realistic interactive 3D scenes with parameterized sources using deep neural operators

We address the challenge of acoustic simulations in three-dimensional (3D) virtual rooms with parametric source positions, which have applications in virtual/augmented reality, game audio, and spatial computing. The wave equation can fully describe wave phenomena such as diffraction and interference. However, conventional numerical discretization methods are computationally expensive when simulating hundreds of source and receiver positions, making simulations with parametric source positions impractical. To overcome this limitation, we propose using deep operator networks to approximate linear wave-equation operators. This enables the rapid prediction of sound propagation in realistic 3D acoustic scenes with parametric source positions, achieving millisecond-scale computations. By learning a compact surrogate model, we avoid the offline calculation and storage of impulse responses for all relevant source/listener pairs. Our experiments, including various complex scene geometries, show good agreement with reference solutions, with root mean squared errors ranging from 0.02 to 0.10 Pa. Notably, our method signifies a paradigm shift as—to our knowledge—no prior machine learning approach has achieved precise predictions of complete wave fields within realistic domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The arbitrary-order virtual element method for linear elastodynamics models. Convergence, stability and dispersion-dissipation analysis.

We design the conforming virtual element method for the numerical approximation of the two dimensional elastodynamics problem. We prove stability and convergence of the semi-discrete approximation and derive optimal error estimates under $\textit{h}$-refinement in both the energy and the $L^2$ norms, and optimal error estimates under $\textit{p}$-refinement in the energy norm. The performance of the proposed virtual element method is assessed on a set of different computational meshes, including non-convex cells up to order four in the h-refinement setting. Exponential convergence is also experimentally observed under p-refinement. Finally, we present a dispersion-dissipation analysis for both the semi-discrete and fully-discrete schemes, showing that polygonal meshes behave as classical simplicial/quadrilateral grids in terms of dispersion-dissipation properties.

97 MATHEMATICS AND COMPUTING↗