Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter continuation method”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Measuring Two Key Parameters of H3 Color Centers in Diamond

A method of measuring two key parameters of H3 color centers in diamond has been created as part of a continuing effort to develop tunable, continuous-wave, visible lasers that would utilize diamond as the lasing medium. (An H3 color center in a diamond crystal lattice comprises two nitrogen atoms substituted for two carbon atoms bonded to a third carbon atom. H3 color centers can be induced artificially; they also occur naturally. If present in sufficient density, they impart a yellow hue.) The method may also be applicable to the corresponding parameters of other candidate lasing media. One of the parameters is the number density of color centers, which is needed for designing an efficient laser. The other parameter is an optical-absorption cross section, which, as explained below, is needed for determining the number density. The present method represents an improvement over prior methods in which optical-absorption measurements have been used to determine absorption cross sections or number densities. Heretofore, in order to determine a number density from such measurements, it has been necessary to know the applicable absorption cross section; alternatively, to determine the absorption cross section from such measurements, it has been necessary to know the number density. If, as in this case, both the number density and the absorption cross section are initially unknown, then it is impossible to determine either parameter in the absence of additional information.

Roberts, W. Thomas↗
DOE OSTI2025

Process intensification approach to enhancing heat and mass transfer: Radio frequency (RF) assisted drying of paper and board

Conventional multi-cylinder drying of paper and board typically relies on both conductive drying from steam-heated cylinders and convective drying, where heated air flows over the paper web surface. Conduction primarily contributes to heat transfer, while convection is the main driver of mass transfer. However, conventional drying systems are heavily dependent on steam, usually powered by fossil fuels, and are often energy-inefficient with high levels of waste. Additionally, these surface-driven processes result in a lower percentage of energy absorption compared to the energy supplied, leading to significant energy losses. To improve this long-standing process, an experimental system was developed to investigate a process intensification approach involving the integration of Radio Frequency (RF) heating, a volumetric electromagnetic technology, alongside traditional conduction and convection drying methods. This study also emphasizes the use of sensors to continuously monitor key parameters such as moisture content, supply system temperatures, sample temperatures, air flows, and drying rates in real-time. The effect of RF as an auxiliary energy source in localized environments at varying moisture levels was explored to optimize industrial drying systems, quantify potential improvements, and provide insights for future studies. Experimental results from trials combining RF with convection and with the base case alternating conduction-convection drying processes are presented. It is shown that RF is a viable process intensification approach for paper drying improving the drying rate and energy intensity at higher moisture contents. These findings offer valuable insights for process intensification and contribute to the process development, modeling, and simulation of advanced paper drying techniques.

42 ENGINEERING↗
DOE OSTI2024

Colloquium: Eigenvector continuation and projection-based emulators

Eigenvector continuation is a computational method for parametric eigenvalue problems that uses subspace projection with a basis derived from eigenvector snapshots from different parameter sets. It is part of a broader class of subspace-projection techniques called reduced-basis methods. In this Colloquium, the development, theory, and applications of eigenvector continuation and projection-based emulators are presented. In conclusion, the basic concepts are introduced, the underlying theory and convergence properties are discussed, and recent applications for quantum systems and future prospects are presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2012

Automated and Manual Rocket Crater Measurement Software

An update has been performed to software designed to do very rapid automated measurements of craters created in sandy substrates by rocket exhaust on liftoff. The previous software was optimized for pristine lab geometry and lighting conditions. This software has been enhanced to include a section for manual measurements of crater parameters; namely, crater depth, crater full width at half max, and estimated crater volume. The tools provide a very rapid method to measure these manual parameters to ease the burden of analyzing large data sets. This software allows for rapid quantization of the rocket crater parameters where automated methods may not work. The progress of spreadsheet data is continuously saved so that data is never lost, and data can be copied to clipboards and pasted to other software for analysis. The volume estimation of a crater is based on the central max depth axis line, and the polygonal shape of the crater is integrated around that axis.

Metzger, Philip↗
NASA NTRS1994

The Influence of Sound Cues on the Maintenance of Temporal Organization in the Sprague-Dawley Rat

Temporal organization is a fundamental property of living matter. From single cells to complex animals including man, most physiological systems undergo daily periodic changes in concert with environmental cues (e.g., light, temperature etc.). It is known that pulsed Environmental synchronizers, zeitgebers, (e.g. light) can modify rhythm parameters. Rhythm stability is a necessary requirement for most animal experiments. The extent to which sound can influence the circadian system of laboratory rats is poorly understood. This has implications to animal habitats in the novel environments of the Space-Laboratory or Space Station. A series of three white noise (88+/-0.82 db) zeitgeber experiments were conducted (n=6/experiment).The sound cue was introduced in the circadian free-running phase (DD-NQ) and in one additional case sound was added to the usual photoperiod (12L:12D) to determine masking effects. Circadian rhythm parameters of drinking frequency, feeding frequency, and gross locomotor activity were continuously monitored. Data analysis for these studies included macroscopic and microscopic methods. Raster plots to visually detect entrainment versus free-running period, were plotted for each animal, for all three parameters, during all sound perturbation tested. These data were processed through a series of detrending (robust locally weighted regression analyses) and complex demodulation analyses. In summary, these findings show that periodic "white" noise "influences" the rats circadian system but does not "entrain" the feeding, drinking or locomotor activity rhythms.

Winget, C. M.↗
NASA NTRS1986

Focusing anisotropy of solar cosmic rays

Analysis of 50-MeV particle data from ICE and 1-GeV data from eight high-latitude neutron monitors of the February 16, 1984 solar cosmic ray event indicates that an exponential with scattering mean free path of greater than 2 AU and with index in the 1.3-1.7 range provides a better description of particle anisotropies than the first-order form, confirming the prediction from cosmic ray transport theories which incorporate the effect of adiabatic focusing into the usual pitch angle scattering formalism. The time profile of the event is largely a reflection of coronal transport processes, with good fits provided by a coronal diffusion coefficient of about 10 to the 18th sq cm/s and with velocity-dependent escape. The simple method of predicting the interplanetary scattering parameters has application where particle injection is nearly continuous, but where previously accelerated particles have filled the heliosphere with a large omnidirectional background.

Bieber, J. W.↗
NASA NTRS2017

SpaRibs Geometry Parameterization for Wings with Multiple Sections using Single Design

The SpaRibs topology of an aircraft wing has a significant effect on its structural behavior and stability as well as the flutter performance. The development of additive manufacturing techniques like Electron Beam Free Form Fabrication (EBF3) has made it feasible to manufacture aircraft wings with curvilinear spars, ribs (SpaRibs) and stiffeners. In this article a new global-local optimization framework for wing with multiple sections using curvilinear SpaRibs is described. A single design space is used to parameterize the SpaRibs geometry. This method has been implemented using MSC-PATRAN to create a broad range of SpaRibs topologies using limited number of parameters. It ensures C0 and C1 continuities in SpaRibs geometry at the junction of two wing sections with airfoil thickness gradient discontinuity as well as mesh continuity between all structural components. This method is advantageous in complex multi-disciplinary optimization due to its potential to reduce the number of design variables. For the global-local optimization the local panels are generated by an algorithm which is totally based on a set algebra on the connectivity matrix data. The great advantage of this method is that it is completely independent of the coordinates of the nodes of the finite element model. It is also independent of the order in which the elements are distributed in the FEM. The code is verified by optimizing of the CRM Baseline model at trim condition at Mach number equal to 0.85 for five different angle of attack (-2deg, 0deg,2deg,4deg and 6deg). The final weight of the wing is 19,090.61 lb. This value is comparable to that obtained by Qiang et al. 6 (19,269 lb).

wing sections↗
NASA NTRS1994

Continuous monitoring of bacterial attachment

A major concern with the Space Station Freedom (SSF) water supply system is the control of longterm microbial contamination and biofilm development in the water storage and distribution systems. These biofilms have the potential for harboring pathogens as well as microbial strains containing resistance factors that could negatively influence crew health. The proposed means for disinfecting the water system on SSF (iodine) may encourage the selection of resistant strains. In fact, biofilm bacteria were observed in water lines from the Space Shuttle Columbia (OV-102); therefore, an alternative remediation method is required to disinfect spacecraft water lines. A thorough understanding of colonization events and the physiological parameters that will influence bacteria adhesion is required. The limiting factor for development of this technology is the ability to continuously monitor adhesion events and the effects of biocides on sessile bacteria. Methods were developed to allow bacterial adhesion and subsequent biocidal treatment to be monitored continuously. This technique couples automated image analysis with a continuous flow of a bacterial suspension through an optical flow cell. A strain of Pseudomonas cepacia isolated from the water supply of the Space Shuttle Discovery (OV-103) during STS-39 was grown in a nitrogen-limited continuous culture. This culture was challenged continuously with iodine during growth, and the adhesion characteristics of this strain was measure with regard to flow rate. Various biocides (ozone, hypochlorite, and iodine) were added to the flow stream to evaluate how well each chemical removed the bacteria. After biocide treatment, a fresh bacterial suspension was introduced into the flow cell, and the attachment rate was evaluated on the previously treated surface. This secondary fouling was again treated with biocide to determine the efficacy of multiple batch chemical treatments in removing biofilm.

Koeing, D. W.↗
NASA NTRS2013

Closed-Loop Process Control for Electron Beam Freeform Fabrication and Deposition Processes

A closed-loop control method for an electron beam freeform fabrication (EBF(sup 3)) process includes detecting a feature of interest during the process using a sensor(s), continuously evaluating the feature of interest to determine, in real time, a change occurring therein, and automatically modifying control parameters to control the EBF(sup 3) process. An apparatus provides closed-loop control method of the process, and includes an electron gun for generating an electron beam, a wire feeder for feeding a wire toward a substrate, wherein the wire is melted and progressively deposited in layers onto the substrate, a sensor(s), and a host machine. The sensor(s) measure the feature of interest during the process, and the host machine continuously evaluates the feature of interest to determine, in real time, a change occurring therein. The host machine automatically modifies control parameters to the EBF(sup 3) apparatus to control the EBF(sup 3) process in a closed-loop manner.

Taminger, Karen M.↗
NASA NTRS1979

Improved prediction of laminar leading edge separation

Research was conducted to provide a definite criterion for the prediction of the bubble burst on airfoils typical of those used for fighter wings. The approach taken was to correlate existing airfoil bubble burst data using various parameters at the laminar separation point. The method due to Weber was modified to provide a continuous analytic solution for the velocity distribution around the airfoil leading edge. Coupling the modified Weber method with the Stratford laminar separation prediction method leads to a universal chart giving the conditions at separation as a function of stagnation location and leading edge radius. Application of the combined method to available two-dimensional airfoil data resulted in an empirical criterion presenting the limiting local velocity gradient at separation as a function of the boundary layer momentum thickness at separation for bubble burst. The correlation leads as well to the qualitative explanation of two types of laminar stall: thin airfoil and leading edge. The validity of the correlation is demonstrated by predicting the lift coefficient and angle of attack for stall on airfoils with leading edge or trailing edge flaps.

Herring, R. N.↗
NASA NTRS2008

Non-Linear System Identification for Aeroelastic Systems with Application to Experimental Data

Representation and identification of a non-linear aeroelastic pitch-plunge system as a model of the NARMAX class is considered. A non-linear difference equation describing this aircraft model is derived theoretically and shown to be of the NARMAX form. Identification methods for NARMAX models are applied to aeroelastic dynamics and its properties demonstrated via continuous-time simulations of experimental conditions. Simulation results show that (i) the outputs of the NARMAX model match closely those generated using continuous-time methods and (ii) NARMAX identification methods applied to aeroelastic dynamics provide accurate discrete-time parameter estimates. Application of NARMAX identification to experimental pitch-plunge dynamics data gives a high percent fit for cross-validated data.

Kukreja, Sunil L.↗
NASA NTRS2008

Nonlinear System Identification for Aeroelastic Systems with Application to Experimental Data

Representation and identification of a nonlinear aeroelastic pitch-plunge system as a model of the Nonlinear AutoRegressive, Moving Average eXogenous (NARMAX) class is considered. A nonlinear difference equation describing this aircraft model is derived theoretically and shown to be of the NARMAX form. Identification methods for NARMAX models are applied to aeroelastic dynamics and its properties demonstrated via continuous-time simulations of experimental conditions. Simulation results show that (1) the outputs of the NARMAX model closely match those generated using continuous-time methods, and (2) NARMAX identification methods applied to aeroelastic dynamics provide accurate discrete-time parameter estimates. Application of NARMAX identification to experimental pitch-plunge dynamics data gives a high percent fit for cross-validated data.

Kukreja, Sunil L.↗
DOE OSTI2024

Multi-parametric analysis for mixed integer linear programming: An application to transmission upgrade and congestion management

Upgrading the capacity of existing transmission lines is essential for meeting the growing energy demands, facilitating the integration of renewable energy, and ensuring the security of the transmission system. This study focuses on the selection of lines whose capacities and by how much should be expanded from the perspective of the Independent System Operators (ISOs) to minimize the total system cost. We employ advanced multi-parametric programming and an enhanced branch-and-bound algorithm to address complex mixed-integer linear programming (MILP) problems, considering multi-period time constraints and physical limitations of generators and transmission lines. To characterize the various decisions in transmission expansion, we model the increased capacity of existing lines as parameters within a specified range. This study first relaxes the binary variables to continuous variables and applies the Lagrange method and Karush-Kuhn-Tucker (KKT) conditions to obtain optimal solutions and identify critical regions associated with active and inactive constraints. Moreover, we extend the traditional branch-and-bound (B&B) method by determining the problem’s upper and lower bounds at each node of the B&B decision tree, helping to manage computational challenges in large-scale MILP problems. Here, we compare the difference between the upper and lower bounds to obtain an approximate optimal solution within the decision-makers’ tolerable error range. In addition, the first derivative of the objective function on the parameters of each line is used to inform the selection of lines for easing congestion and maximizing social welfare. Finally, the capacity upgrades are selected by weighing the reductions in system costs against the expense of upgrading line capacities. The findings are supported by numerical simulations and provide transmission-line planners with decision-making guidance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

A provably stable numerical method for the anisotropic diffusion equation in confined magnetic fields

We present a novel numerical method for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box which is accurate and provably stable. We derive energy estimates of the solution of the continuous initial boundary value problem. A discrete formulation is presented using operator splitting in time with the summation by parts finite difference approximation of spatial derivatives for the perpendicular diffusion operator. Weak penalty procedures are derived for implementing both boundary conditions and parallel diffusion operator obtained by field line tracing. We prove that the fully-discrete approximation is unconditionally stable. Discrete energy estimates are shown to match the continuous energy estimate given the correct choice of penalty parameters. A nonlinear penalty parameter is shown to provide an effective method for tuning the parallel diffusion penalty and significantly minimises rounding errors. Several numerical experiments, using manufactured solutions, the “NIMROD benchmark” problem and a single island problem, are presented to verify numerical accuracy, convergence, and asymptotic preserving properties of the method. Finally, we present a magnetic field with chaotic regions and islands and show the contours of the anisotropic diffusion equation reproduce key features in the field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2015

Three-Dimensional Piecewise-Continuous Class-Shape Transformation of Wings

Class-Shape Transformation (CST) is a popular method for creating analytical representations of the surface coordinates of various components of aerospace vehicles. A wide variety of two- and three-dimensional shapes can be represented analytically using only a modest number of parameters, and the surface representation is smooth and continuous to as fine a degree as desired. This paper expands upon the original two-dimensional representation of airfoils to develop a generalized three-dimensional CST parametrization scheme that is suitable for a wider range of aircraft wings than previous formulations, including wings with significant non-planar shapes such as blended winglets and box wings. The method uses individual functions for the spanwise variation of airfoil shape, chord, thickness, twist, and reference axis coordinates to build up the complete wing shape. An alternative formulation parameterizes the slopes of the reference axis coordinates in order to relate the spanwise variation to the tangents of the sweep and dihedral angles. Also discussed are methods for fitting existing wing surface coordinates, including the use of piecewise equations to handle discontinuities, and mathematical formulations of geometric continuity constraints. A subsonic transport wing model is used as an example problem to illustrate the application of the methodology and to quantify the effects of piecewise representation and curvature constraints.

Olson, Erik D.↗
NASA NTRS2023

KMT-2021-BLG-1150Lb: Microlensing Planet Detected Through a Densely Covered Planetary-Caustic Signal

Aims. Recently, there have been reports of various types of degeneracies in the interpretation of planetary signals induced by planetary caustics. In this work we check whether such degeneracies persist in the case of well-covered signals by analyzing the lensing event KMT-2021-BLG-1150, the light curve of which exhibits a densely and continuously covered short-term anomaly. Methods. In order to identify degenerate solutions, we thoroughly investigated the parameter space by conducting dense grid searches for the lensing parameters. We then checked the severity of the degeneracy among the identified solutions. Results. We identify a pair of planetary solutions resulting from the well-known inner-outer degeneracy, and find that interpreting the anomaly is not subject to any degeneracy other than the inner-outer degeneracy. The measured parameters of the planet separation (normalized to the Einstein radius) and mass ratio between the lens components are ( s , q ) in ∼ (1.297, 1.10 × 10 −3 ) for the inner solution and ( s , q ) out ∼ (1.242, 1.15 × 10 −3 ) for the outer solution. According to a Bayesian estimation, the lens is a planetary system consisting of a planet with a mass M p = 0.88 +0.38 −0.36 M J and its host with a mass M h = 0.73 +0.32 −0.30 M ⊙ lying toward the Galactic center at a distance D L = 3.8 +1.3 −1.2 kpc. By conducting analyses using mock data sets prepared to mimic those obtained with data gaps and under various observational cadences, we find that gaps in data can result in various degenerate solutions, while the observational cadence does not pose a serious degeneracy problem as long as the anomaly feature can be delineated.

Gravitational microlensing↗
NASA NTRS1988

Domain decomposition methods for the parallel computation of reacting flows

Domain decomposition is a natural route to parallel computing for partial differential equation solvers. Subdomains of which the original domain of definition is comprised are assigned to independent processors at the price of periodic coordination between processors to compute global parameters and maintain the requisite degree of continuity of the solution at the subdomain interfaces. In the domain-decomposed solution of steady multidimensional systems of PDEs by finite difference methods using a pseudo-transient version of Newton iteration, the only portion of the computation which generally stands in the way of efficient parallelization is the solution of the large, sparse linear systems arising at each Newton step. For some Jacobian matrices drawn from an actual two-dimensional reacting flow problem, comparisons are made between relaxation-based linear solvers and also preconditioned iterative methods of Conjugate Gradient and Chebyshev type, focusing attention on both iteration count and global inner product count. The generalized minimum residual method with block-ILU preconditioning is judged the best serial method among those considered, and parallel numerical experiments on the Encore Multimax demonstrate for it approximately 10-fold speedup on 16 processors.

Keyes, David E.↗
NASA NTRS1992

A piecewise continuous Timoshenko beam model for the dynamic analysis of tapered beam-like structures

Distributed parameter modeling offers a viable alternative to the finite element approach for modeling large flexible space structures. The introduction of the transfer matrix method into the continuum modeling process provides a very useful tool to facilitate the distributed parameter model applied to some more complex configurations. A uniform Timoshenko beam model for the estimation of the dynamic properties of beam-like structures has given comparable results. But many aeronautical and aerospace structures are comprised of non-uniform sections or sectional properties, such as aircraft wings and satellite antennas. This paper proposes a piecewise continuous Timoshenko beam model which is used for the dynamic analysis of tapered beam-like structures. A tapered beam is divided into several segments of uniform beam elements. Instead of arbitrarily assumed shape functions used in finite element analysis, the closed-form solution of the Timoshenko beam equation is used. Application of the transfer matrix method relates all the elements as a whole. By corresponding boundary conditions and compatible conditions a characteristic equation for the global tapered beam has been developed, from which natural frequencies can be derived. A computer simulation is shown in this paper, and compared with the results obtained from the finite element analysis. While piecewise continuous Timoshenko beam model decreases the number of elements significantly; comparable results to the finite element method are obtained.

Shen, Ji Yao↗