Search NASA⌕ Search

SEARCH · Search NASA

Results for “parametric modeling”

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

Regenerative life support system research

Sections on modeling, experimental activities during the grant period, and topics under consideration for the future are contained. The sessions contain discussions of: four concurrent modeling approaches that were being integrated near the end of the period (knowledge-based modeling support infrastructure and data base management, object-oriented steady state simulations for three concepts, steady state mass-balance engineering tradeoff studies, and object-oriented time-step, quasidynamic simulations of generic concepts); interdisciplinary research activities, beginning with a discussion of RECON lab development and use, and followed with discussions of waste processing research, algae studies and subsystem modeling, low pressure growth testing of plants, subsystem modeling of plants, control of plant growth using lighting and CO2 supply as variables, search for and development of lunar soil simulants, preliminary design parameters for a lunar base life support system, and research considerations for food processing in space; and appendix materials, including a discussion of the CELSS Conference, detailed analytical equations for mass-balance modeling, plant modeling equations, and parametric data on existing life support systems for use in modeling.

Source record↗

Multi-fidelity parametric sensitivity estimation for large eddy simulation with the Spalart–Allmaras model

A computationally affordable approach to estimate parametric sensitivities of engineering relevant quantities of interest for a large eddy simulation (LES) is explored. The method is based on defining a Reynolds-averaged Navier–Stokes (RANS) problem that is constrained to reproduce the LES mean flow field. Here, the proposed method is described and assessed for a shock/boundary layer interaction problem, where the shock angle and wall temperature are considered variable or uncertain. In the current work, we show that the proposed method offers improved sensitivity predictions for certain flow features as compared to standalone RANS simulations, while using a fraction of the LES cost.

42 ENGINEERING↗

A Parametric Study of Erupting Flux Rope Rotation: Modeling the 'Cartwheel CME' on 9 April 2008

The rotation of erupting filaments in the solar corona is addressed through a parametric simulation study of unstable, rotating flux ropes in bipolar force-free initial equilibrium. The Lorentz force due to the external shear-field component and the relaxation of tension in the twisted field are the major contributors to the rotation in this model, while reconnection with the ambient field is of minor importance, due to the field's simple structure. In the low-beta corona, the rotation is not guided by the changing orientation of the vertical field component's polarity inversion line with height. The model yields strong initial rotations which saturate in the corona and differ qualitatively from the profile of rotation vs. height obtained in a recent simulation of an eruption without preexisting flux rope. Both major mechanisms writhe the flux rope axis, converting part of the initial twist helicity, and produce rotation profiles which, to a large part, are very similar within a range of shear-twist combinations. A difference lies in the tendency of twist-driven rotation to saturate at lower heights than shear-driven rotation. For parameters characteristic of the source regions of erupting filaments and coronal mass ejections, the shear field is found to be the dominant origin of rotations in the corona and to be required if the rotation reaches angles of order 90 degrees and higher; it dominates even if the twist exceeds the threshold of the helical kink instability. The contributions by shear and twist to the total rotation can be disentangled in the analysis of observations if the rotation and rise profiles are simultaneously compared with model calculations. The resulting twist estimate allows one to judge whether the helical kink instability occurred. This is demonstrated for the erupting prominence in the "Cartwheel CME" on 9 April 2008, which has shown a rotation of approximately 115 deg. up to a height of 1.5 Solar R above the photosphere. Out of a range of initial equilibria which include strongly kink-unstable (Phi = 5 pi), weakly kink-unstable (Phi = 3.5 pi), and kink-stable (Phi = 2.5 pi) configurations, only the evolution of the weakly kink-unstable flux rope matches the observations in their entirety.

initiation and propagation↗

Deep Learning Explicit Differentiable Predictive Control Laws for Buildings

We present a differentiable predictive control (DPC) methodology for learning constrained control laws for unknown nonlinear systems. DPC poses an approximate solution to multiparametric programming problems emerging from explicit nonlinear model predictive control (MPC). Contrary to approximate MPC, DPC does not require supervision by an expert controller. Instead, a system dynamics model is learned from a small dataset of recorded observations of the perturbed system's dynamics and the control law is optimized offline by interaction with the learned system model. The DPC method is based on two sequential steps, i) system identification using a constrained neural state-space model, and ii) optimization of an explicit control law parametrized by another neural network in closed-loop simulation with the identified neural state-space model. The combination of a differentiable closed-loop system and penalty methods for constraint handling of system outputs and inputs allows us to optimize the control law's parameters directly by backpropagating economic MPC loss through the learned system model. By incorporating domain knowledge and leveraging established techniques from optimal control, our method leverages deep neural networks as nonlinear function approximators for system identification and control while avoiding concomitant costs of intractably large datasets, and computationally expensive over-parametrized models. The scalability, data efficiency, and constrained optimal control capability of the proposed DPC method are demonstrated in simulation using a multi-zone building emulator.

Drgona, Jan↗

A practical method for modeling temporally-averaged ocean wave frequency-directional spectra for characterizing wave energy climates

Wave energy resource characterization investigations have been hindered because frequency-directional wave spectra are not broadly available for many coastal regions due to the insufficient spatial coverage of buoy observation networks and spectral wave model hindcast outputs. We address this problem by developing and validating a practical method for approximating temporally-averaged frequency-directional wave spectra averaged over periods of months or years from bulk wave parameters, which are available at a much greater coverage and resolution than frequency-directional wave spectra. While the temporally averaged frequency-directional wave spectrum over these periods cannot be used for analyzing a single sea state, it aggregates multiple sea states, identifies dominant wave energy systems representing wave climate, and resolves their spectral characteristics. Therefore, modelling temporally averaged frequency-directional wave spectra is of great value for planning and designing wave energy projects, e.g., resource characterization, site assessment, and conceptual design of wave energy converters. Temporally-averaged frequency-directional wave spectra and related important wave energy parameters approximated using this method are found to be more accurate than commonly used parametric wave spectrum models. This method can be applied to a wide range of wave climates given its universality and high accuracy. Also, the availability of bulk wave parameters from multi-decade high-resolution wave model hindcasts is increasing.

16 TIDAL AND WAVE POWER↗

Parametrically designed surface topography on CAD models of additively manufactured lattice structures for improved design validation

Additively Manufactured (AM) lattice structures of the laser powder bed fusion (LPBF) process generally have different surface geometries depending on the overhang angle and location. This means that the design validation stage is often challenging, considering it is based on an ideal Computer Aided Design (CAD) model that is not truly representative, since the AM part will be different in terms of dimensional accuracy and surface finish. Previous studies have relied on the design of surface textures that are independent of the overhang angle, or techniques based on directly using X-ray Computed Tomography (XCT) data. In this paper, a new technique for designing surface texture on the CAD of lattices has been investigated and correlated with areal surface roughness parameters. After extracting areal surface parameters from the XCT data of the lattice, the method allowed for the design of this surface texture on the part CAD taking in consideration if its up skin or down skin, without using XCT data as the main input. By applying this method, it is possible to obtain a lattice CAD model with a designed surface texture and geometry that is more representative of the actual AM lattice. The mean deviation between the CAD model with the designed surface and the XCT was a third of the one between the XCT and the initial CAD. The proposed method allows for designing and replicating LPBF AM surfaces on the CAD of a lattice, taking into consideration the dimensional deviation caused by AM surfaces, especially on overhangs. Assuming that an LPBF AM process is stable and produces approximately the same AM surface, this method can be used to predict the geometry of a lattice, providing a cost-effective and efficient methodology that minimizes the time needed for design validation.

42 ENGINEERING↗

Geometric optics model of Thomson scattering enhanced by parametric coupling

We formulate a theory of Thomson scattering (TS) from inhomogeneous plasmas using a ray-tracing approach. The scattered radiation is described within the eikonal approximation, including convective gain arising from scattering instabilities driven by the TS probe laser. In the limit of a weak probe and vanishing gain, the theory recovers the standard TS cross section and scattered power due to volumetric particle noise. The ray-tracing formalism is implemented in a TS fitting framework, enabling accurate analysis of inhomogeneous plasmas whose background profiles can be provided by large-scale hydrodynamic simulations or experimental diagnostics. Here, several examples of synthetic TS spectra are presented, and further extensions of the theory are discussed.

Carleton, D. E. [University of Alberta, Edmonton, ↗

A Novel Framework for Parametric Analysis of Coastal Transition Zone Modeling

Abstract Vulnerability of coastal regions to extreme events motivates an operational coupled inland‐coastal modeling strategy focusing on the coastal transition zone (CTZ), an area between the coast and upland river. To tackle this challenge, we propose a top‐down framework for investigating the contribution of different processes to the hydrodynamics of CTZs with various geometrical shapes, different physical properties, and under several forcing conditions. We further propose a novel method, called tidal vanishing point (TVP), for delineating the extent of CTZs through the upland. We demonstrate the applicability of our framework over the United States East and Gulf coasts. We categorize CTZs in the region into three classes, namely, without estuary (direct river–coast connection), triangular‐, and trapezoidal‐shaped estuary. The results show that although semidiurnal tidal constituents are dominant in most cases, diurnal tidal constituents become more prevalent in the river segment as the discharge increases. Also, decreasing the bed roughness value promotes more significant changes in the results than increasing it by the same value. Additionally, the estuary promotes tidal energy attenuation and consequently decreases the reach of tidal signals through the upland. The proposed framework is generic and extensible to any coastal region.

Chegini, Taher↗

A parametric study of dissociation and ionization models at 12 km/sec

Thermochemical nonequilibrium-solution dependence on available models for the chemical reaction rates is examined. Solutions from the Kang and Dunn (1973) reaction-rate set, the Park rate set of 1987, and the Park rate set of 1991 are compared. The blunt-nosed, axisymmetric geometry considered is a 60-deg sphere cone with nose radius of 1.07 m and cicular aft skirt. The nonequilibrium test case is 12 km/sec entry into the earth's atmosphere at 80 km altitude. The model variations are implemented into the Langley aerothermodynamics upwind relaxation algorithm code. While variations in the reaction rates have no effect on the surface pressure distribution and little effect on the convective heating, the effect on degree of ionization and radiative heating can be a factor of three.

Mitcheltree, R. A.↗

On the Numerical Formulation of Parametric Linear Fractional Transformation (LFT) Uncertainty Models for Multivariate Matrix Polynomial Problems

Robust control system analysis and design is based on an uncertainty description, called a linear fractional transformation (LFT), which separates the uncertain (or varying) part of the system from the nominal system. These models are also useful in the design of gain-scheduled control systems based on Linear Parameter Varying (LPV) methods. Low-order LFT models are difficult to form for problems involving nonlinear parameter variations. This paper presents a numerical computational method for constructing and LFT model for a given LPV model. The method is developed for multivariate polynomial problems, and uses simple matrix computations to obtain an exact low-order LFT representation of the given LPV system without the use of model reduction. Although the method is developed for multivariate polynomial problems, multivariate rational problems can also be solved using this method by reformulating the rational problem into a polynomial form.

Belcastro, Christine M.↗

Parametric Evaluation of RDE Inlet Performance using a Reduced-Order Model

This paper presents the results of a parametric study on the effects of inlet restriction, along with a number of other controlling parameters, using a previously developed reduced-order RDE inlet model. The results of this study will be used to inform subsequent high-fidelity CFD simulations and/or selection of experimental operating conditions which may be conducive to maximizing performance.

Bedick, Clinton↗

Real-Time Adaptive Drag Minimization Wind Tunnel Investigation of a Flexible Wing with Variable Camber Continuous Trailing Edge Flap System

This paper reports the results of a recently completed real-time adaptive drag minimization wind tunnel investigation of a highly flexible wing wind tunnel model equipped with the Variable Camber Continuous Trailing Flap (VCCTEF) technology at the University of Washington Aeronautical Laboratory (UWAL). The wind tunnel investigation is funded by NASA SBIR Phase II contract with Scientific Systems Company, Inc. (SSCI) and University of Washington (UW) as a subcontractor. The wind tunnel model is a sub-scale Common Research Model (CRM) wing constructed of foam core and fiberglass skin and is aeroelastically scaled to achieve a wing tip deflection of 10% of the wing semi-span which represents a typical wing tip deflection for a modern transport such as Boeing 787. The jig-shape twist of the CRM wing is optimized using a CART3D aero-structural model to achieve the minimum induced drag for the design cruise lift coefficient of 0.5. The wing is equipped with two chord wise cambered segments for each of the six span wise flap sections for a total of 12 individual flap segments that comprise the VCCTEF system. Each of the 12 flap segments is actively controlled by an electric servo-actuator. The real-time adaptive drag optimization strategy includes an on-board aerodynamic model identification, a model excitation, and a real-time drag optimization. The on-board aerodynamic model is constructed parametrically as a function of the angle of attack and flap positions to model the lift and drag coefficients of the wing. The lift coefficient models include a linear model and a second-order model. The drag coefficient models include a quadratic model and a higher-order up to 6th-order model to accurately model the drag coefficient at high angles of attack. The onboard aerodynamic model identification includes a recursive least-squares (RLS) algorithm and a batch least-squares (BLS) algorithm designed to estimate the model parameters. The model excitation method is designed to sample the input set that comprises the angle of attack and the flap positions. Three model excitation methods are developed: random excitation method, sweep method, and iterative angle-of-attack seeking method. The real-time drag optimization includes a generic algorithm developed by SSCI and several optimization methods developed by NASA which include a second-order gradient Newton-Raphson optimization method, an iterative gradient optimization method, a pseudo-inverse optimization method, an analytical optimization method, and an iterative refinement optimization method. The first wind tunnel test entry took place in September 2017. This test revealed major hardware issues and required further redesign of the flap servo mechanisms. The second test entry took place in April 2018. However, the test was not successful due to the issues with the onboard aerodynamic model identification RLS algorithm which incorrectly identified model parameters. This test also provides an experimental comparison study between the VCCTEF and a variable camber discrete trailing edge flap (VCDTEF) without the elastomer transition mechanisms. The experimental result confirms the benefit of the VCCTEF which produces lower drag by 5% than the VCDTEF. The third and final test entry took place in June 2018 after the issues with the RLS algorithm have been identified and corrected. Additional improvements were implemented. These include the BLS algorithm, the iterative angle-of-attack seeking method, the iterative gradient optimization method, and the pseudo-inverse optimization method. The test objectives were successfully demonstrated as the real-time drag optimization identifies several optimal solutions at off-design lift coefficients. The iterative gradient optimization method is found to achieve up to 4.7% drag reduction for the off-design lift coefficient of 0.7. The pseudo-inverse optimization method which does not require the drag coefficient model is found to be quite effective in reducing drag. Up to 9.4% drag reduction for the off-design lift coefficient of 0.7 is achieved with the pseudo-inverse optimization method. The wind tunnel investigation demonstrates the potential of real-time drag optimization technology. Several new capabilities are developed that could enable future adaptive wing technologies for flexible wings equipped with drag control devices such as the VCCTEF.

Drag Optimization↗

Parametric and Sensitivity Analysis of a Steam Generator Model Using Python and Machine-Learning Tools

For this study, we used Python and machine-learning tools to perform a comprehensive parametric and sensitivity analysis on a steam generator (SG) model. (The Python model was based on a previously completed MATLAB framework for the Holtec SMR-160 SG.) We investigated the influence of various input parameters (e.g., heat transfer coefficient [HTC], Nusselt number, and heat exchanger effectiveness) on the system’s output. With machine-learning tools such as the Risk Analysis Virtual Environment (RAVEN), which was developed at Idaho National Laboratory, we were then able to perform an automated analysis of the SG inputs’ effect on the HTC. The analysis results give valuable insights into the performance and optimization of SG systems. We found the inlet mass flow rate (MFR) to have the greatest impact on the HTC, followed closely by the inlet temperature, and then pressure. Shifting of the input parameters causes the location of the maximum HTC along the SG length to change incrementally. The cold leg (CL) MFR was also found to impact the HTC magnitude as well as the location of the maximum HTC. At between 0.4–0.9 of the total SG length, the input parameters experience maximum impact on the HTC, leading us to suggest that sensors be efficiently placed on the SG so as to closely and effectively monitor thermal-hydraulic properties during reactor operation. We also found that the sensitivity data calculated manually agrees with the RAVEN – based data, confirming the same range of maximum sensitivity. However, the RAVEN-based analysis showed that cold leg pressure and hot leg temperature have a greater impact on the heat transfer coefficient than the mass flow rate, implying that a manual sensitivity study taking only two samples is not accurate.

20 FOSSIL-FUELED POWER PLANTS↗

Process Cost Modeling for Multi-Disciplinary Design Optimization

For early design concepts, the conventional approach to cost is normally some kind of parametric weight-based cost model. There is now ample evidence that this approach can be misleading and inaccurate. By the nature of its development, a parametric cost model requires historical data and is valid only if the new design is analogous to those for which the model was derived. Advanced aerospace vehicles have no historical production data and are nowhere near the vehicles of the past. Using an existing weight-based cost model would only lead to errors and distortions of the true production cost. This report outlines the development of a process-based cost model in which the physical elements of the vehicle are costed according to a first-order dynamics model. This theoretical cost model, first advocated by early work at MIT, has been expanded to cover the basic structures of an advanced aerospace vehicle. Elemental costs based on the geometry of the design can be summed up to provide an overall estimation of the total production cost for a design configuration. This capability to directly link any design configuration to realistic cost estimation is a key requirement for high payoff MDO problems. Another important consideration in this report is the handling of part or product complexity. Here the concept of cost modulus is introduced to take into account variability due to different materials, sizes, shapes, precision of fabrication, and equipment requirements. The most important implication of the development of the proposed process-based cost model is that different design configurations can now be quickly related to their cost estimates in a seamless calculation process easily implemented on any spreadsheet tool. In successive sections, the report addresses the issues of cost modeling as follows. First, an introduction is presented to provide the background for the research work. Next, a quick review of cost estimation techniques is made with the intention to highlight their inappropriateness for what is really needed at the conceptual phase of the design process. The First-Order Process Velocity Cost Model (FOPV) is discussed at length in the next section. This is followed by an application of the FOPV cost model to a generic wing. For designs that have no precedence as far as acquisition costs are concerned, cost data derived from the FOPV cost model may not be accurate enough because of new requirements for shape complexity, material, equipment and precision/tolerance. The concept of Cost Modulus is introduced at this point to compensate for these new burdens on the basic processes. This is treated in section 5. The cost of a design must be conveniently linked to its CAD representation. The interfacing of CAD models and spreadsheets containing the cost equations is the subject of the next section, section 6. The last section of the report is a summary of the progress made so far, and the anticipated research work to be achieved in the future.

Bao, Han P.↗

Surface and Volume Grid Generation in Parametric Form

The algorithm for surface modeling and volume grid generation using parametric Non-Uniform Rational B-splines (NURBS) geometric representation are presented. The enhanced re-parameterization algorithm which can yield a desired physical distribution on the curve, surface, and volume is also presented. This approach bridges the gap between computer aided design surface/volume definition and surface/volume grid generation.

Yu, Tzuyi↗

An Image-Plane Approach to Gravitational Lens Modeling of Interferometric Data

Strong gravitational lensing acts as a cosmic telescope, enabling the study of the high-redshift universe. Astronomical interferometers, such as the Atacama Large Millimeter/submillimeter Array (ALMA), have provided high-resolution images of strongly lensed sources at millimeter and submillimeter wavelengths. To model the mass and light distributions of lensing and source galaxies from strongly lensed images, strong lens modeling for interferometric observations is conventionally performed in the visibility space, which is computationally expensive. In this paper, we implement an image-plane lens modeling methodology for interferometric dirty images by accounting for noise correlations. We show that the image-plane likelihood function produces accurate model values when tested on simulated ALMA observations with an ensemble of noise realizations. We also apply our technique to ALMA observations of two sources selected from the South Pole Telescope survey, comparing our results with previous visibility-based models. Our model results are consistent with previous models for both parametric and pixelated source-plane reconstructions. We implement this methodology for interferometric lens modeling in the open-source software package lenstronomy.

Zhang, Nan [Illinois U., Urbana (main)] (ORCID:000↗

A stochastic damage model for the rupture prediction of a multi-phase solid. I - Parametric studies. II - Statistical approach

A stochastic damage model for predicting the rupture of a brittle multiphase material is developed, based on the microcrack-macrocrack interaction. The model, which incorporates uncertainties in locations, orientations, and numbers of microcracks, characterizes damage by microcracking and fracture by macrocracking. A parametric study is carried out to investigate the change of the stress intensity at the macrocrack tip by the configuration of microcracks. The inherent statistical distribution of the fracture toughness arising from the intrinsic random nature of microcracks is explored using a statistical approach. For this purpose, a computer simulation model is introduced, which incorporates a statistical characterization of geometrical parameters of a random microcrack array.

Lua, Yuan J.↗