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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Computer Simulations of Bottle Brushes: From Melts to Soft Networks

We use a combination of Molecular dynamics simulations and analytical calculations, and study dens bottle-brush systems in a melt and network State. Analysis of our simulation results shows that bottle-brush macromolecules in melt behave as ideal chains with effective Kuhn length b K . Simulations show that the bottle-brush-induced bending rigidity is due to an entropy decrease caused by redistribution of the side chains upon backbone bending. The Kuhn length of the bottle:brushes increases with increasing the side-chain degree of polymerization n sc as b K proportional to n sc 0.46 . Moreover, this model of bottle brush macromolecules is extended to describe mechanical properties of bottle brush networks in linear and nonlinear deformation regimes. In the linear deformation regime, the network shear modulus scales with the degree of polymerization of the side chains as G 0 proportional to (n sc + 1) -1 as long as the ratio of the Kuhn length, b K , to the size of the fully extended bottle-brush backbone between cross-links, R-max, is smaller than unity, b K /R max << 1. Bottle-brush networks With b K /R max proportional to 1 demonstrate behavior similar to that of networks Of semiflexible chains with G 0 proportional to n sc -0.5 . Finally, in the nonlinear network deformation regime, the deformation-dependent shear modulus is a universal function of the first strain invariant I 1 and bottle-brush backbone deformation ratio beta describing stretching ability of the bottle-brush backbone between cross-links.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Spin wave Hamiltonian and anomalous scattering in NiPS 3

Here, we report a comprehensive spin wave analysis of the semiconducting honeycomb van der Waal antiferromagnet NiPS 3 . Using single-crystal inelastic neutron scattering, we map out the full Brillouin zone and fit the observed modes to a spin wave model with rigorously defined uncertainty. We find that the third-nearest-neighbor exchange J 3 dominates the Hamiltonian, a feature which we fully account for by ab initio density functional theory calculations. We also quantify the degree to which the threefold rotation symmetry is broken and account for the Q = 0 excitations observed in other measurements, yielding a spin exchange model which is consistent across multiple experimental probes. We also identify a strongly reduced static ordered moment and reduced low-energy intensity relative to the linear spin wave calculations, signaling unexplained features in the magnetism which requires going beyond the linear spin wave approximation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Multi-Segment Decentralized Control Strategies for Renewables-Rich Microgrids in Extreme Conditions

Microgrids provide a promising approach to accommodating various distributed energy resources (DERs), while requiring significant communication infrastructures that may be affected by extreme conditions such as natural disasters and cyber-attacks. In this paper, a fully decentralized control strategy without the need for communication is proposed for islanded microgrids with high renewables penetration. First, special multi-segment power/frequency characteristic curves are designed, so that different DERs can be automatically coordinated in a prioritized manner such as renewables first to maintain power balance while DER frequencies are regulated at their reference values. Second, piecewise linear served load versus frequency models are designed to prioritize loads according to their significance, so that only noncritical loads will be curtailed as needed while critical loads are supplied without any interruptions during power deficiency. The proposed strategy can effectively deal with various normal and extreme system conditions including 100% renewables penetration, loads and renewables variations, power deficiencies requiring load curtailments, disconnection of existing DERs, connection of new DERs as well as network sectionalization and reconfiguration. As a result, the proposed control strategy is validated in the real-time digital simulator (RTDS) model of the IIT Campus Microgrid to demonstrate its effectiveness in enhancing the resiliency of renewables-rich microgrids in extreme conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1995

Aeroelasticity and structural optimization of composite helicopter rotor blades with swept tips

This report describes the development of an aeroelastic analysis capability for composite helicopter rotor blades with straight and swept tips, and its application to the simulation of helicopter vibration reduction through structural optimization. A new aeroelastic model is developed in this study which is suitable for composite rotor blades with swept tips in hover and in forward flight. The hingeless blade is modeled by beam type finite elements. A single finite element is used to model the swept tip. Arbitrary cross-sectional shape, generally anisotropic material behavior, transverse shears and out-of-plane warping are included in the blade model. The nonlinear equations of motion, derived using Hamilton's principle, are based on a moderate deflection theory. Composite blade cross-sectbnal properties are calculated by a separate linear, two-dimensional cross section analysis. The aerodynamic loads are obtained from quasi-steady, incompressible aerodynamics, based on an implicit formulation. The trim and steady state blade aeroelastic response are solved in a fully coupled manner. In forward flight, where the blade equations of motion are periodic, the coupled trim-aeroelastic response solution is obtained from the harmonic balance method. Subsequently, the periodic system is linearized about the steady state response, and its stability is determined from Floquet theory.

Yuan, K. A.↗
NASA NTRS2020

The HadGEM3-GA7.1 radiative kernel: the importance of a well-resolved stratosphere

We present top-of-atmosphere and surface radiative kernels based on the atmospheric component (GA7.1) of the HadGEM3 general circulation model developed by the UK Met Office. We show that the utility of radiative kernels for forcing adjustments in idealised CO2 perturbation experiments is greatest where there is sufficiently high resolution in the stratosphere in both the target climate model and the radiative kernel. This is because stratospheric cooling to a CO2 perturbation continues to increase with height, and low-resolution or low-top kernels or climate model output are unable to fully resolve the full stratospheric temperature adjustment. In the sixth phase of the Coupled Model Intercomparison Project (CMIP6), standard atmospheric model data are available up to 1 hPa on 19 pressure levels, which is a substantial advantage compared to CMIP5. We show in the IPSL-CM6A-LR model where a full set of climate diagnostics are available that the HadGEM3-GA7.1 kernel exhibits linear behaviour and the residual error term is small, as well as from a survey of kernels available in the literature that in general low-top radiative kernels underestimate the stratospheric temperature response. The HadGEM3-GA7.1 radiative kernels are available at https://doi.org/10.5281/zenodo.3594673 (Smith, 2019).

radiative kernels↗
DOE OSTI2022

A fully Lagrangian, non-parametric bias model for dark matter halos

Here, we present a non-parametric Lagrangian biasing model and fit the ratio of the halo and mass densities at the field level using the mass-weighted halo field in the AbacusSummit simulations at z=0.5. Unlike the perturbative halo bias model that has been widely used in interpreting the observed large-scale structure traced by galaxies, we find a non-negative halo-to-mass ratio that increases monotonically with the linear overdensity δ 1 in the initial Lagrangian space. The bias expansion, however, does not guarantee non-negativity of the halo counts, and may lead to rising halo number counts at negative overdensities. The shape of the halo-to-mass ratio is unlikely to be described by a polynomial function of δ 1 and other quantities. Especially for massive halos with 6×10 12 h -1 M ⊙ , the halo-to-mass ratio starts soaring up at δ 1 >0, substantially different from the predictions of the bias expansion. We show that for the halo masses we consider (M>3×10 11 h -1 M ⊙ ) a non-parametric halo-to-mass ratio as a function of δ 1 and its local derivative ∇^2δ 1 can recover the halo power spectra to sub-percent level accuracy for wavenumbers k=0.01-0.1 h Mpc -1 given a proper smoothing scale to filter the initial density field, even though we do not fit the power spectrum directly. However, there is mild dependence of the recovery of the halo power spectrum on the smoothing scale and other input parameters. At k<0.01 h Mpc -1 and for massive halos with M>6x10 12 h -1 M ⊙ , our non-parametric model leads to a few percent overestimation of the halo power spectrum, indicating the need for larger or multiple smoothing scales. The halo-to-mass ratios obtained qualitatively agree with intuitions from extended Press-Schechter theory. We compare our framework to the bias expansion and discuss possible extensions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

A composite neural network that learns from multi-fidelity data: Application to function approximation and inverse PDE problems

Presently the training of neural networks relies on data of comparable accuracy but in real applications only a very small set of high-fidelity data is available while inexpensive lower fidelity data may be plentiful. We propose a new composite neural network (NN) that can be trained based on multi-fidelity data. It is comprised of three NNs, with the first NN trained using the low-fidelity data and coupled to two high-fidelity NNs, one with activation functions and another one without, in order to discover and exploit nonlinear and linear correlations, respectively, between the low-fidelity and the high-fidelity data. We first demonstrate the accuracy of the new multi-fidelity NN for approximating some standard benchmark functions but also a 20-dimensional function that is not easy to approximate with other methods, e.g. Gaussian process regression. Subsequently, we extend the recently developed physics-informed neural networks (PINNs) to be trained with multi-fidelity data sets (MPINNs). MPINNs contain four fully-connected neural networks, where the first one approximates the low-fidelity data, while the second and third construct the correlation between the low- and high-fidelity data and produce the multi-fidelity approximation, which is then used in the last NN that encodes the partial differential equations (PDEs). Specifically, by decomposing the correlation into a linear and nonlinear part, the present model is capable of learning both the linear and complex nonlinear correlations between the low- and high-fidelity data adaptively. By training the MPINNs, we can: (1) obtain the correlation between the low- and high-fidelity data, (2) infer the quantities of interest based on a few scattered data, and (3) identify the unknown parameters in the PDEs. In particular, we employ the MPINNs to learn the hydraulic conductivity field for unsaturated flows as well as the reactive models for reactive transport. The results demonstrate that MPINNs can achieve relatively high accuracy based on a very small set of high-fidelity data. Despite the relatively low dimension and limited number of fidelities (two-fidelity levels) for the benchmark problems in the present study, the proposed model can be readily extended to very high-dimensional regression and classification problems involving multi-fidelity data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

A fast and robust computational modeling approach for density and shape predictions in powder metallurgy hot isostatic pressing

Powder metallurgy hot isostatic pressing (PM-HIP) is an advanced manufacturing process that produces near-net-shape parts with high material utilization and uniform microstructures. PM-HIP is frequently used for producing small-scale parts with complicated geometries and is potentially economical for producing large-scale parts. However, excessive post-HIP shape distortions can reduce its effectiveness and economic advantage, especially for larger parts. A PM-HIP computational model can predict and help mitigate these distortions. However, due to complex deformation mechanisms and thermo-mechanical coupling present in PM-HIP processes, these non-linear computational models sometimes become numerically unstable. The numerical instabilities in these models can lead to very slow convergence or no convergence at all, which often translates to slow and unreliable models. These limitations are more pronounced in large models with complicated geometries. Hence, in this work, an alternative modeling approach is presented that improves numerical stability and computational performance. The presented approach achieves these improvements through approximating the fully coupled thermo-mechanical PM-HIP model as a decoupled model and adding inertial damping to the model’s mechanical part. In conclusion, a comparison with the fully coupled model indicated a slight dip in prediction accuracy (<5% error) but significant improvements in numerical stability (>20 times larger time step size) and computational performance (5-10 times speed-up with less computational resource usage) when using the presented approach.

Hot isostatic pressing↗
NASA NTRS1978

Aircraft model prototypes which have specified handling-quality time histories

Several techniques for obtaining linear constant-coefficient airplane models from specified handling-quality time histories are discussed. The pseudodata method solves the basic problem, yields specified eigenvalues, and accommodates state-variable transfer-function zero suppression. The algebraic equations to be solved are bilinear, at worst. The disadvantages are reduced generality and no assurance that the resulting model will be airplane like in detail. The method is fully illustrated for a fourth-order stability-axis small motion model with three lateral handling quality time histories specified. The FORTRAN program which obtains and verifies the model is included and fully documented.

Johnson, S. H.↗
NASA NTRS1990

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 2: User's guide

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 2 is the User's Guide, and describes the program's general features, the input and output, the procedure for setting up initial conditions, the computer resource requirements, the diagnostic messages that may be generated, the job control language used to run the program, and several test cases.

Towne, Charles E.↗
NASA NTRS1990

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 3: Programmer's reference

A new computer code was developed to solve the 2-D or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating-direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 3 is the Programmer's Reference, and describes the program structure, the FORTRAN variables stored in common blocks, and the details of each subprogram.

Towne, Charles E.↗
NASA NTRS2001

Finite Element Methods for real-time Haptic Feedback of Soft-Tissue Models in Virtual Reality Simulators

We have applied the linear elastic finite element method to compute haptic force feedback and domain deformations of soft tissue models for use in virtual reality simulators. Our results show that, for virtual object models of high-resolution 3D data (>10,000 nodes), haptic real time computations (>500 Hz) are not currently possible using traditional methods. Current research efforts are focused in the following areas: 1) efficient implementation of fully adaptive multi-resolution methods and 2) multi-resolution methods with specialized basis functions to capture the singularity at the haptic interface (point loading). To achieve real time computations, we propose parallel processing of a Jacobi preconditioned conjugate gradient method applied to a reduced system of equations resulting from surface domain decomposition. This can effectively be achieved using reconfigurable computing systems such as field programmable gate arrays (FPGA), thereby providing a flexible solution that allows for new FPGA implementations as improved algorithms become available. The resulting soft tissue simulation system would meet NASA Virtual Glovebox requirements and, at the same time, provide a generalized simulation engine for any immersive environment application, such as biomedical/surgical procedures or interactive scientific applications.

Frank, Andreas O.↗
NASA NTRS2011

Estuarine Sediment Deposition during Wetland Restoration: A GIS and Remote Sensing Modeling Approach

Restoration of the industrial salt flats in the San Francisco Bay, California is an ongoing wetland rehabilitation project. Remote sensing maps of suspended sediment concentration, and other GIS predictor variables were used to model sediment deposition within these recently restored ponds. Suspended sediment concentrations were calibrated to reflectance values from Landsat TM 5 and ASTER using three statistical techniques -- linear regression, multivariate regression, and an Artificial Neural Network (ANN), to map suspended sediment concentrations. Multivariate and ANN regressions using ASTER proved to be the most accurate methods, yielding r2 values of 0.88 and 0.87, respectively. Predictor variables such as sediment grain size and tidal frequency were used in the Marsh Sedimentation (MARSED) model for predicting deposition rates for three years. MARSED results for a fully restored pond show a root mean square deviation (RMSD) of 66.8 mm (<1) between modeled and field observations. This model was further applied to a pond breached in November 2010 and indicated that the recently breached pond will reach equilibrium levels after 60 months of tidal inundation.

Newcomer, Michelle↗
DOE OSTI2014

Experimental characterization and modeling of non-linear coupling of the LHCD power on Tore Supra

To achieve steady state operation on future tokamaks, in particular on ITER, the unique capability of a LHCD system to efficiently drive off-axis non-inductive current is needed. In this context, it is of prime importance to study and master the coupling of LH wave to the core plasma at high power density (tens of MW/m{sup 2}). In some specific conditions, deleterious effects on the LHCD coupling are sometimes observed on Tore Supra. At high power the waves may modify the edge parameters that change the wave coupling properties in a non-linear manner. In this way, dedicated LHCD experiments have been performed using the LHCD system of Tore Supra, composed of two different conceptual designs of launcher: the Fully Active Multijunction (FAM) and the new Passive Active Multijunction (PAM) antennas. A nonlinear interaction between the electron density and the electric field has been characterized in a thin plasma layer in front of the two LHCD antennas. The resulting dependence of the power reflection coefficient with the LHCD power, leading occasionally to trips in the output power, is not predicted by the standard linear theory of the LH wave coupling. Therefore, it is important to investigate and understand the possible origin of such non-linear effects in order to avoid their possible deleterious consequences. The PICCOLO-2D code, which self-consistently treats the wave propagation in the antenna vicinity and its interaction with the local edge plasma density, is used to simulate Tore Supra discharges. The simulation reproduces very well the occurrence of a non-linear behavior in the coupling observed in the LHCD experiments. The important differences and trends between the FAM and the PAM antennas, especially a larger increase in RC for the FAM, are also reproduced by the PICCOLO-2D simulation. The working hypothesis of the contribution of the ponderomotive effect in the non-linear observations of LHCD coupling is therefore validated through this comprehensive modeling for the first time on the FAM and PAM antennas on Tore Supra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

The development and verification of a highly accurate collision prediction model for automated noncoplanar plan delivery

Purpose: Significant dosimetric benefits had been previously demonstrated in highly noncoplanar treatment plans. In this study, the authors developed and verified an individualized collision model for the purpose of delivering highly noncoplanar radiotherapy and tested the feasibility of total delivery automation with Varian TrueBeam developer mode. Methods: A hand-held 3D scanner was used to capture the surfaces of an anthropomorphic phantom and a human subject, which were positioned with a computer-aided design model of a TrueBeam machine to create a detailed virtual geometrical collision model. The collision model included gantry, collimator, and couch motion degrees of freedom. The accuracy of the 3D scanner was validated by scanning a rigid cubical phantom with known dimensions. The collision model was then validated by generating 300 linear accelerator orientations corresponding to 300 gantry-to-couch and gantry-to-phantom distances, and comparing the corresponding distance measurements to their corresponding models. The linear accelerator orientations reflected uniformly sampled noncoplanar beam angles to the head, lung, and prostate. The distance discrepancies between measurements on the physical and virtual systems were used to estimate treatment-site-specific safety buffer distances with 0.1%, 0.01%, and 0.001% probability of collision between the gantry and couch or phantom. Plans containing 20 noncoplanar beams to the brain, lung, and prostate optimized via an in-house noncoplanar radiotherapy platform were converted into XML script for automated delivery and the entire delivery was recorded and timed to demonstrate the feasibility of automated delivery. Results: The 3D scanner measured the dimension of the 14 cm cubic phantom within 0.5 mm. The maximal absolute discrepancy between machine and model measurements for gantry-to-couch and gantry-to-phantom was 0.95 and 2.97 cm, respectively. The reduced accuracy of gantry-to-phantom measurements was attributed to phantom setup errors due to the slightly deformable and flexible phantom extremities. The estimated site-specific safety buffer distance with 0.001% probability of collision for (gantry-to-couch, gantry-to-phantom) was (1.23 cm, 3.35 cm), (1.01 cm, 3.99 cm), and (2.19 cm, 5.73 cm) for treatment to the head, lung, and prostate, respectively. Automated delivery to all three treatment sites was completed in 15 min and collision free using a digital Linac. Conclusions: An individualized collision prediction model for the purpose of noncoplanar beam delivery was developed and verified. With the model, the study has demonstrated the feasibility of predicting deliverable beams for an individual patient and then guiding fully automated noncoplanar treatment delivery. This work motivates development of clinical workflows and quality assurance procedures to allow more extensive use and automation of noncoplanar beam geometries.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Supernova constraints on multi-coupled dark energy

The persisting consistency of ever more accurate observational data with the predictions of the standard ΛCDM cosmological model puts severe constraints on possible alternative scenarios, but still does not shed any light on the fundamental nature of the cosmic dark sector. As large deviations from a ΛCDM cosmology are ruled out by data, the path to detect possible features of alternative models goes necessarily through the definition of cosmological scenarios that leave almost unaffected the background and — to a lesser extent — the linear perturbations evolution of the universe. In this context, the Multi-coupled DE (McDE) model was proposed by Baldi [9] as a particular realization of an interacting Dark Energy field characterized by an effective screening mechanism capable of suppressing the effects of the coupling at the background and linear perturbation level. In the present paper, for the first time, we challenge the McDE scenario through a direct comparison with real data, in particular with the luminosity distance of Type Ia supernovae. By studying the existence and stability conditions of the critical points of the associated background dynamical system, we select only the cosmologically consistent solutions, and confront their background expansion history with data. Confirming previous qualitative results, the McDE scenario appears to be fully consistent with the adopted sample of Type Ia supernovae, even for coupling values corresponding to an associated scalar fifth-force about four orders of magnitude stronger than standard gravity. Our analysis demonstrates the effectiveness of the McDE background screening, and shows some new non-trivial asymptotic solutions for the future evolution of the universe. Clearly, linear perturbation data and, even more, nonlinear structure formation properties are expected to put much tighter constraints on the allowed coupling range. Nonetheless, our results show how the background expansion history might be highly insensitive to the fundamental nature and to the internal complexity of the dark sector.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Electron temperature gradient instability and transport analysis in NSTX and NSTX-U plasmas

Extensive linear and nonlinear simulations to study electron temperature gradient (ETG) stability and thermal transport in National Spherical Torus Experiment (NSTX) and NSTX-U plasmas were performed using the fully electromagnetic gyrokinetic code CGYRO. Linear simulations were performed to determine ETG thresholds in different discharges, showing that ETG modes in spherical tokamaks can present different scalings compared to conventional aspect-ratio tokamaks. Nonlinear gyrokinetic simulations were conducted for selected cases to calculate electron thermal transport and compare to experimental values. Results are also compared with those of ETG modes in the multi-mode model and the Trapped-Gyro-Landau-Fluid reduced model codes, to better understand their applicability in spherical tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Hybrid Fixed-Point Fixed-Stress Splitting Method for Linear Poroelasticity

Efficient and accurate poroelasticity models are critical in modeling geophysical problems such as oil exploration, gas-hydrate detection, and hydrogeology. We propose an efficient operator splitting method for Biot’s model of linear poroelasticity based on fixed-point iteration and constrained stress. In this method, we eliminate the constraint on time step via combining our fixed-point approach with a physics-based restraint between iterations. Three different cases are considered to demonstrate the stability and consistency of the method for constant and variable parameters. The results are validated against the results from the fully coupled approach. In case I, a single iteration is used for continuous coefficients. The relative error decreases with an increase in time. In case II, material coefficients are assumed to be linear. In the single iteration approach, the relative error grows significantly to 40% before rapidly decaying to zero. This is an artifact of the approximate solutions approaching the asymptotic solution. The error in the multiple iterations oscillates within 10 − 6 before decaying to the asymptotic solution. Nine iterations per time step are enough to achieve the relative error close to 10 − 7 . In the last case, the hybrid method with multiple iterations requires approximately 16 iterations to make the relative error 5 × 10 − 6 .

Delgado, Paul M.↗