Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface nonlinearity”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Nonlinear convergence in contact mechanics: Immersed boundary finite volume

In this report we present an immersed boundary finite volume (IBM) method for simulating quasistatic contact mechanics of linearly elastic domains at small strains. In IBM, all external boundaries and internal contacts of an object are represented by embedded surfaces inside a Cartesian mesh, which need not conform to the grid lines. The contact constraints consist of the non-penetrability condition and Coulomb’s friction law, which are discretized using special interpolation stencils and enforced via penalty parameters. The resulting nonlinear system depends on displacement unknowns only. To solve it, we use the Newton method but find that it diverges frequently. To understand the divergence pattern, we analyze a simplified 2-cell problem and show that the global convergence of Newton cannot be ensured for any choice of penalty parameters. We thus propose a modified Newton solver, which guarantees convergence for the 2-cell problem and is numerically verified to converge for all the challenging simulations considered herein. While both 1 st - and 2 nd -order variants of IBM, in displacement unknowns, are proposed, the modified Newton solver applies only to the 1 st -order variant.

42 ENGINEERING↗
NASA NTRS2005

Contamination and Radiation Effects on Nonlinear Crystals for Space Laser Systems

Space Lasers are vital tools for NASA s space missions and military applications. Although, lasers are highly reliable on the ground, several past space laser missions proved to be short-lived and unreliable. In this communication, we are shedding more light on the contamination and radiation issues, which are the most common causes for optical damages and laser failures in space. At first, we will present results based on the study of liquids and subsequently correlate these results to the particulates of the laser system environment. We present a model explaining how the laser beam traps contaminants against the optical surfaces and cause optical damages and the role of gravity in the process. We also report the results of the second harmonic generation efficiency for nonlinear optical crystals irradiated with high-energy beams of protons. In addition, we are proposing to employ the technique of adsorption to minimize the presence of adsorbing molecules present in the laser compartment.

Abdeldayem, Hossain A.↗
NASA NTRS1972

Analysis of an electrohydraulic aircraft control surface servo and comparison with test results

An analysis of an electrohydraulic aircraft control-surface system is made in which the system is modeled as a lumped, two-mass, spring-coupled system controlled by a servo valve. Both linear and nonlinear models are developed, and the effects of hinge-moment loading are included. Transfer functions of the system and approximate literal factors of the transfer functions for several cases are presented. The damping action of dynamic pressure feedback is analyzed. Comparisons of the model responses with results from tests made on a highly resonant rudder control-surface servo indicate the adequacy of the model. The effects of variations in hinge-moment loading are illustrated.

Edwards, J. W.↗
DOE OSTI2022

Variational encoder geostatistical analysis (VEGAS) with an application to large scale riverine bathymetry

Estimation of riverbed profiles, also known as bathymetry, plays a vital role in many applications, such as safe and efficient inland navigation, prediction of bank erosion, land subsidence, and flood risk management. The high cost and complex logistics of direct bathymetry surveys, i.e, depth imaging, have encouraged the use of indirect measurements such as surface flow velocities. However, estimating high-resolution bathymetry from indirect measurements is an inverse problem that can be computationally challenging. Here, we propose a reduced-order model (ROM) based approach that utilizes a variational autoencoder (VAE), a type of deep neural network with a narrow layer in the middle, to compress bathymetry and flow velocity information and accelerate bathymetry inverse problems from flow velocity measurements. In our application, the shallow-water equations (SWE) with appropriate boundary conditions (BCs), e.g., the discharge and/or the free surface elevation, constitute the forward problem, to predict flow velocity. Then, ROMs of the SWEs are constructed on a nonlinear manifold of low dimensionality through a variational encoder and the bathymetry inversion problem is derived on the low-dimensional latent space in a Hierarchical Bayesian setting. Further, the reformulation allows variational inference with a small number (e.g., $\mathscr{O}$ (100) of ROM runs and efficient uncertainty quantification. We have tested our inversion approach on a one-mile reach of the Savannah River, GA, USA. Once the neural network is trained (offline stage), the proposed technique can perform the inversion operation orders of magnitude faster than traditional inversion methods that are commonly based on linear projections, such as principal component analysis (PCA), or the principal component geostatistical approach (PCGA). Furthermore, tests show that the algorithm can estimate the bathymetry with good accuracy even with sparse flow velocity measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

High-order accurate finite difference discretisations on fully unstructured dual quadrilateral meshes

Here, we present a novel approach for high-order accurate numerical differentiation on unstructured meshes of quadrilateral elements. To differentiate a given function, an auxiliary function with greater smoothness properties is defined which when differentiated provides the derivatives of the original function. The method generalises traditional finite difference methods to meshes of arbitrary topology in any number of dimensions for any order of derivative and accuracy. We demonstrate the accuracy of the numerical scheme using dual quadrilateral meshes and a refinement method based on subdivision surfaces. The scheme is applied to the solution of a range of partial differential equations, including both linear and nonlinear, second and fourth order equations, and a time-dependent first order equation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A new woven composite constitutive model validated by shock wave experiments

In this paper, we present results of plate impact simulations of shock compressed woven glass fiber-reinforced plastic (GRP) performed using the Arbitrary Lagrangian–Eulerian three-dimensional finite element code. A hyperelastic large-strain-based empirical Continuum Damage Mechanics (CDM) formulation is employed to describe damage initiation and growth in the shock-compressed GRP. The model parameters calibration scheme utilizes the Velocity Interferometer System for Any Reflector normal particle velocity measurements at the free surface of the GRP target plates. The impact velocity in the experiments ranged from 8.5 to 418 m/s. The finite element model considered planar 0°/90° bidirectional plies with an individual ply thickness of 0.68 mm, stacked to reach a total laminate thickness of 6.8 mm. The anisotropic elastic strains were estimated from the experimentally determined tetragonal symmetry stiffness matrix for the GRP. The strain-based damage model captures several salient features observed in the measured free surface particle wave profiles, including the shock rise time, onset of Elastic—Elastic Cracking, and the shape of the nonlinear portion of the experimental particle velocity profiles. Furthermore, the CDM model predicts the dominant damage mode to be matrix microcracking due to shear and the associated bulk expansion (bulking) under the global compressive loading in the plate impact configuration.

42 ENGINEERING↗
DOE OSTI2022

Resistive wall tearing mode disruptions in DIII-D and ITER tokamaks

Disruptions are a serious problem in tokamaks, in which thermal and magnetic energy confinement is lost. This paper uses data from the DIII-D experiment, theory, and simulations to demonstrate that resistive wall tearing modes (RWTMs) produce the thermal quench (TQ) in a typical locked mode shot. Analysis of the linear RWTM dispersion relation shows the parameter dependence of the growth rate, particularly on the resistive wall time. Linear simulations of the locked mode equilibrium show that it is unstable with a resistive wall and stable with an ideally conducting wall. Nonlinear simulations demonstrate that the RWTM grows to sufficient amplitude to cause a complete thermal quench. The RWTM growth time is proportional to the thermal quench time. Here, the nonlinearly saturated RWTM magnetic perturbation amplitude agrees with experimental measurements. The onset condition is that the q = 2 rational surface is sufficiently close to the resistive wall. Collectively, this identifies the RWTM as the cause of the TQ. In ITER, RWTMs will produce long TQ times compared to present-day experiments. ITER disruptions may be significantly more benign than previously predicted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1977

Solar oscillations

The main observational results and related theoretical investigations concerning solar oscillations are reviewed. The normal modes of oscillation of the sun are classified according to their angular eigenvalues and the number of their radial nodes. Observations of excited normal modes are discussed, particularly in relation to five-minute oscillations, nonlinearly coupled unstable g-modes, oscillations of the sun's apparent diameter, and reported oscillations of the solar surface with a period of 2 hr 40 min. Linear stability calculations are briefly examined, and stochastic excitation of p-modes by turbulent convection is considered. The five-minute oscillations are described as the superposition of many excited nonradial p-modes.

Goldreich, P.↗
NASA NTRS1979

The effects of ground hydrology on climate sensitivity to solar constant variations

The effects of two different evaporation parameterizations on the climate sensitivity to solar constant variations are investigated by using a zonally averaged climate model. The model is based on a two-level quasi-geostrophic zonally averaged annual mean model. One of the evaporation parameterizations tested is a nonlinear formulation with the Bowen ratio determined by the predicted vertical temperature and humidity gradients near the earth's surface. The other is the linear formulation with the Bowen ratio essentially determined by the prescribed linear coefficient.

Chou, S. H.↗
NASA NTRS1992

Geometric modeling of inflatable structures for lunar base

A modular inflatable structure consisting of thin, composite membranes is presented for use in a lunar base. Results from a linear elastic analysis of the structure indicate that it is feasible in the lunar environment. Further analysis requires solving nonlinear equations and accurately specifying the geometries of the structural members. A computerized geometric modeling technique, using bicubic Bezier surfaces to generate the geometries of the inflatable structure, was conducted. Simulated results are used to create three-dimensional wire frames and solid renderings of the individual components of the inflatable structure. The component geometries are connected into modules, which are then assembled based upon the desired architecture of the structure.

Nowak, Paul S.↗
NASA NTRS1992

TranAir: A full-potential, solution-adaptive, rectangular grid code for predicting subsonic, transonic, and supersonic flows about arbitrary configurations. User's manual

The TranAir computer program calculates transonic flow about arbitrary configurations at subsonic, transonic, and supersonic freestream Mach numbers. TranAir solves the nonlinear full potential equations subject to a variety of boundary conditions modeling wakes, inlets, exhausts, porous walls, and impermeable surfaces. Regions with different total temperature and pressure can be represented. The user's manual describes how to run the TranAir program and its graphical support programs.

Johnson, F. T.↗
NASA NTRSDate not supplied

Multidisciplinary Design Optimization and Analysis of an Open Rotor Stage: Part 1

Successful design of open rotor propulsors requires effective analysis across multiple disciplines, including aerodynamics, acoustics, and structures. A viable design must not only be efficient but must also produce an acceptable level of noise and meet all static and dynamic structural requirements. For design and optimization, this is especially challenging because running high fidelity analyses is resource-intensive, and optimizing a design may require many thousands of cases to be analyzed. For this reason, the NASA team has applied design methodology that utilizes low-cost aerodynamic methods, machine learning models, and high-fidelity analyses when necessary. This includes standard two-dimensional methods such as throughflow analysis and quasi-3D blade-to-blade CFD, as well as some newly developed methods. Optimization using 3D CFD is necessary to maximize performance, and this is considered as well. All optimizations are carried out subject to structural constraints evaluated using finite element analysis. Doing this accurately requires a robust trunnion design, capable of pitching the blade between cruise and takeoff conditions while maintaining acceptable factor of safety. Hot to cold analysis must also be applied in order to correctly determine the as-manufactured shape. For acoustics, the unsteady pressures on the blade surfaces must be predicted, and this can be done either through full-annulus unsteady CFD or through a nonlinear harmonic method (NLH). NLH can run much faster, allowing some acoustic considerations to be made earlier in the design process. The design process is ongoing, and this presentation will review the current status and planned next steps. This part of the talk will focus on aerodynamic performance and be followed by a talk on structures and acoustics.

Design↗
DOE OSTI2020

Feedbacks Between Surface Deformation and Permafrost Degradation in Ice Wedge Polygons, Arctic Coastal Plain, Alaska

In the past three decades, an abrupt, pan-Arctic acceleration of ice wedge melting has transformed tundra landscapes, spurring the formation of hummock-like features known as high-centered polygons (HCPs). This rapid geomorphic transition profoundly alters regional hydrology and influences surface emissions of CO 2 and CH 4 . In Arctic Alaska, most recent instances of ice wedge degradation have arrested within 15–20 years of inception, stabilizing HCP microtopography. However, feedbacks between ground surface deformation and permafrost stability are incompletely understood, limiting our capacity to predict trajectories of landscape evolution in a still warmer future. In this work, we use field data from a site near Prudhoe Bay, Alaska, to develop a modeling-based framework for assessing the strength of positive (i.e., exacerbating) feedbacks on ice wedge degradation, focusing on the importance of heterogeneity in surface drainage and microtopographic conditions. Our simulations suggest that, when troughs are narrow, positive feedbacks on ice wedge melting (associated with thermokarst pool formation) are relatively weak. Positive feedbacks are markedly stronger beneath wide troughs, such as those that form above older, larger ice wedges. Seasonal thaw abruptly accelerates once a talik begins to form beneath wide and deep thermokarst pools. Once a talik initiates, winter severity and snowpack thickness increase in importance as predictors of thaw intensity in summer. Our results indicate that meter-scale heterogeneity in polygonal microtopography potentially exerts strong, nonlinear controls on thermokarst trajectories. These findings are useful for predicting future thermokarst dynamics and for interpreting the results from coarser-resolution land surface models operating at greater spatial and temporal scales.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Assessment of Coupling One-Dimension and Three-Dimension CFD in Conjugate Heat Transfer Modeling and Simulations for Cooling of Turbine Engine

The CHT (conjugate heat transfer) modeling and simulation work fits under a bigger area in the TTT (Transformational Tools and Technologies) project called “Cooling Tools and Technologies for Combustor and Turbine.” The existing validation data sets such as TTT GRC (Glenn Research Center) THX (Turbulent Heat FLUX) IV experiments will be used to assess the solid/fluid heat transfer development and validation. The experimental configuration was a subsonic convergent nozzle with a square exit blowing a heated exhaust over an effusion-cooled deck. The simulation solutions were obtained and compared to experimental measurements of the profiles of flow velocities and temperatures, as well as deck surface temperatures. Cooling air blowing ratios of 0, 1, and 2 were considered. The three-dimensional cooling tubes were modelled by the medial axis of the cylindrical geometries due to large length over width ratios. It represented a typical coupling of one-dimension and three-dimension simulations of conjugate heat transfer approaches. The surface temperatures on the solid surfaces and the profiles of velocities and temperatures near the surfaces compared to the experimental data were presented.

nonlinear sub-filter stresses↗
DOE OSTI2020

Simulating coupled surface–subsurface flows with ParFlow v3.5.0: capabilities, applications, and ongoing development of an open-source, massively parallel, integrated hydrologic model

Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.

58 GEOSCIENCES↗
NASA NTRS2018

Görtler Instability and Its Control via Surface Suction over an Axisymmetric Cone at Mach 6

The characteristics of Görtler instability over an axisymmetric cone with an aft concave section are studied via linear and nonlinear instability analysis and direct numerical simulations. Several options for the cone geometry have been investigated numerically, subject to a fixed forecone section and constraints on the maximum cone diameter, overall cone length, and minimum N-factor for the most amplified Görtler modes. Computations show that it is possible to design a cone with a peak N-factor of Nmax > 8 at the target Reynolds number of 12.1×10(exp 6) per meter, corresponding to the maximum quiet Reynolds number in the Boeing/AFOSR Mach-6 Quiet Tunnel at Purdue University. Direct numerical simulations show that an array of roughness elements corresponding to a peak roughness height of 0.1006 mm at the center can excite Görtler vortices that evolve into sufficiently strong streamwise streaks that may break down via high-frequency secondary instability. Thus, the selected axisymmetric configuration of interest should provide an acceptable baseline to investigate the feasibility of several aspects of laminar-flow control via boundary-layer suction. The apparatus that is being used for measurements in the Boeing/AFOSR Quiet Tunnel is described, along with some preliminary experimental results.

Gortler vortices↗
NASA NTRS1995

Thermal Stress Analysis of RCG-Tempered TUFI Tile TPS for Hypersonic Vehicles

This paper presents detailed results from linear and nonlinear finite-element thermal stress analyses of a new tile, Thermal Protection Systems (TPS) concept. A very thin coating of Reaction Cured Glass (RCG) is used to "temper" the surface of Toughened Uni-Piece Fibrous Insulation (TUFI) tiles to improve resistance to thermal shock and thermal cycling effects. The coating also serves to reduce catalytic heating and may improve waterproofing. Calculations include trajectory-based aerothermal heating environments for X-34 wing leading edge TPS designs and arc jet environments for TPS test articles. The nonlinear analyses include the high temperature plasticity of RCG to demonstrate the reuseability of the material.

Milos, Frank S.↗