SEARCH · Search NASA

Results for “Initialization 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Solution of partial differential equations on vector and parallel computers

The present status of numerical methods for partial differential equations on vector and parallel computers was reviewed. The relevant aspects of these computers are discussed and a brief review of their development is included, with particular attention paid to those characteristics that influence algorithm selection. Both direct and iterative methods are given for elliptic equations as well as explicit and implicit methods for initial boundary value problems. The intent is to point out attractive methods as well as areas where this class of computer architecture cannot be fully utilized because of either hardware restrictions or the lack of adequate algorithms. Application areas utilizing these computers are briefly discussed.

Ortega, J. M.↗

NASA NTRS1985

Solution of partial differential equations on vector and parallel computers

Ortega, J. M.↗

NASA NTRS2017

Laser Pulse Bidirectional Reflectance from CALIPSO Mission

In this Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) study, we present a simple way of determining laser pulse bidirectional reflectance over snow/ice surface using the Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) 532 nanometer polarization channels' measurements. The saturated laser pulse returns from snow and ice surfaces are recovered based on surface tail information. The method overview and initial assessment of the method performance will be presented. The retrieved snow surface bidirectional reflectance is compared with reflectance from both CALIOP cloud cover regions and Moderate Resolution Imaging Spectroradiometer (Earth Observing System (EOS)) (MODIS) Bi-directional Reflectance Distribution Function (BRDF) / Albedo model parameters. The comparisons show that the snow surface bidirectional reflectance over Antarctica for saturation region are generally reliable with a mean value of about 0.90 plus or minus 0.10, while the mean surface reflectance from cloud cover region is about 0.84 plus or minus 0.13 and the calculated MODIS reflectance at 555 nanometers from the BRDF / Albedo model with near nadir illumination and viewing angles is about 0.96 plus or minus 0.04. The comparisons here demonstrate that the snow surface reflectance underneath the cloud with cloud optical depth of about 1 is significantly lower than that for a clear sky condition.

Lu, Xiaomei↗

NASA NTRS1997

Efficient Gradient-Based Shape Optimization Methodology Using Inviscid/Viscous CFD

The formerly developed preconditioned-biconjugate-gradient (PBCG) solvers for the analysis and the sensitivity equations had resulted in very large error reductions per iteration; quadratic convergence was achieved whenever the solution entered the domain of attraction to the root. Its memory requirement was also lower as compared to a direct inversion solver. However, this memory requirement was high enough to preclude the realistic, high grid-density design of a practical 3D geometry. This limitation served as the impetus to the first-year activity (March 9, 1995 to March 8, 1996). Therefore, the major activity for this period was the development of the low-memory methodology for the discrete-sensitivity-based shape optimization. This was accomplished by solving all the resulting sets of equations using an alternating-direction-implicit (ADI) approach. The results indicated that shape optimization problems which required large numbers of grid points could be resolved with a gradient-based approach. Therefore, to better utilize the computational resources, it was recommended that a number of coarse grid cases, using the PBCG method, should initially be conducted to better define the optimization problem and the design space, and obtain an improved initial shape. Subsequently, a fine grid shape optimization, which necessitates using the ADI method, should be conducted to accurately obtain the final optimized shape. The other activity during this period was the interaction with the members of the Aerodynamic and Aeroacoustic Methods Branch of Langley Research Center during one stage of their investigation to develop an adjoint-variable sensitivity method using the viscous flow equations. This method had algorithmic similarities to the variational sensitivity methods and the control-theory approach. However, unlike the prior studies, it was considered for the three-dimensional, viscous flow equations. The major accomplishment in the second period of this project (March 9, 1996 to March 8, 1997) was the extension of the shape optimization methodology for the Thin-Layer Navier-Stokes equations. Both the Euler-based and the TLNS-based analyses compared with the analyses obtained using the CFL3D code. The sensitivities, again from both levels of the flow equations, also compared very well with the finite-differenced sensitivities. A fairly large set of shape optimization cases were conducted to study a number of issues previously not well understood. The testbed for these cases was the shaping of an arrow wing in Mach 2.4 flow. All the final shapes, obtained either from a coarse-grid-based or a fine-grid-based optimization, using either a Euler-based or a TLNS-based analysis, were all re-analyzed using a fine-grid, TLNS solution for their function evaluations. This allowed for a more fair comparison of their relative merits. From the aerodynamic performance standpoint, the fine-grid TLNS-based optimization produced the best shape, and the fine-grid Euler-based optimization produced the lowest cruise efficiency.

Baysal, Oktay↗

NASA NTRS1980

Evolution of solar magnetic fields - A new approach to MHD initial-boundary value problems by the method of nearcharacteristics

A method of analysis for the MHD initial-boundary problem is presented in which the model's formulation is based on the method of nearcharacteristics developed by Werner (1968) and modified by Shin and Kot (1978). With this method, the physical causality relationship can be traced from the perturbation to the response as in the method of characteristics, while achieving the advantage of a considerable reduction in mathematical procedures. The method offers the advantage of examining not only the evolution of nonforce free fields, but also the changes of physical conditions in the atmosphere accompanying the evolution of magnetic fields. The physical validity of the method is demonstrated with examples, and their significance in interpreting observations is discussed.

Nakagawa, Y.↗

NASA NTRS1981

Elasto visco-plastic flow with special attention to boundary conditions

A simple but nontrivial steady-state creeping elasto visco-plastic (Maxwell fluid) radial flow problem is analyzed, with special attention given to the effects of the boundary conditions. Solutions are obtained through integration of a governing equation on stress using the Runge-Kutta method for initial value problems and finite differences for boundary value problems. A more general approach through the finite element method, an approach that solves for the velocity field rather than the stress field and that is applicable to a wide range of problems, is presented and tested using the radial flow example. It is found that steady-state flows of elasto visco-plastic materials are strongly influenced by the state of stress of material as it enters the region of interest. The importance of this boundary or initial condition in analyses involving materials coming into control volumes from unusual stress environments is emphasized.

Shimazaki, Y.↗

NASA NTRS1986

Anisotropic constitutive model for nickel base single crystal alloys: Development and finite element implementation

A tool for the mechanical analysis of nickel base single crystal superalloys, specifically Rene N4, used in gas turbine engine components is developed. This is achieved by a rate dependent anisotropic constitutive model implemented in a nonlinear three dimensional finite element code. The constitutive model is developed from metallurigical concepts utilizing a crystallographic approach. A non Schmid's law formulation is used to model the tension/compression asymmetry and orientation dependence in octahedral slip. Schmid's law is a good approximation to the inelastic response of the material in cube slip. The constitutive equations model the tensile behavior, creep response, and strain rate sensitivity of these alloys. Methods for deriving the material constants from standard tests are presented. The finite element implementation utilizes an initial strain method and twenty noded isoparametric solid elements. The ability to model piecewise linear load histories is included in the finite element code. The constitutive equations are accurately and economically integrated using a second order Adams-Moulton predictor-corrector method with a dynamic time incrementing procedure. Computed results from the finite element code are compared with experimental data for tensile, creep and cyclic tests at 760 deg C. The strain rate sensitivity and stress relaxation capabilities of the model are evaluated.

Dame, L. T.↗

NASA NTRS1990

Swept wing ice accretion modeling

An effort to develop a three-dimensional modeling method was initiated. This first step towards creation of a complete aircraft icing simulation code builds on previously developed methods for calculating three-dimensional flow fields and particle trajectories combined with a two-dimensional ice accretion calculation along coordinate locations corresponding to streamlines. This work is a demonstration of the types of calculations necessary to predict a three-dimensional ice accretion. Results of calculations using the 3-D method for a MS-317 swept wing geometry are projected onto a 2-D plane normal to the wing leading edge and compared to 2-D results for the same geometry. It is anticipated that many modifications will be made to this approach, however, this effort will lay the groundwork for future modeling efforts. Results indicate that the flow field over the surface and the particle trajectories differed for the two calculations. This led to lower collection efficiencies, convective heat transfer coefficients, freezing fractions, and ultimately ice accumulation for the 3-D calculation.

Potapczuk, Mark G.↗

NASA NTRS1990

Swept wing ice accretion modeling

An effort to develop a three-dimensional ice accretion modeling method is initiated. This first step toward creation of a complete aircraft icing simulation code builds on previously developed methods for calculating three-dimensional flowfields and particle trajectories combined with a two-dimensional ice accretion calculation along coordinate locations corresponding to streamlines. This work is intended as a demonstration of the types of calculations necessary to predict a three-dimensional ice accretion. Results of calculations using the 3D method for a MS-317 swept wing geometry are projected onto a 2D plane normal to the wing leading edge and compared to 2D results for the same geometry. These results indicate that the flowfield over the surface and the particle trajectories differed for the two calculations. This led to lower collection efficiencies, convective heat transfer coefficients, freezing fractions, and ultimately ice accumulation for the 3D calculation.

Potapczuk, M. G.↗

NASA NTRS1991

Numerical simulation of ice growth on a MS-317 swept wing geometry

An effort to develop a 3-D ice accretion modeling method was initiated. This first step towards creation of a complete aircraft icing simulation code builds on previously developed methods for calculating 3-D flow fields and particle trajectories combined with a 2-D ice accretion calculation along coordinate locations corresponding to streamlines. The types of calculations necessary to predict 3-D ice accretion is demonstrated. Results of calculations using 3-D method for a MS-317 swept wing geometry are projected onto a 2-D plane parallel to the free stream direction and compared to experimental results for the same geometry. It is anticipated that many modifications will be made to this approach, however this effort will lay the groundwork for future modeling efforts. Results indicate that rime ice shapes indicate a difficulty in accurately calculating the ice shape in the runback region.

Potapczuk, M. G.↗

NASA NTRS1991

Numerical simulation of ice growth on a MS-317 swept wing geometry

An effort to develop a 3-D ice accretion modeling method was initiated. This first step towards creation of a complete aircraft icing simulation code builds on previously developed methods for calculating 3-D flow fields and particle trajectories combined with a 2-D ice accretion calculation along coordinate locations corresponding to streamlines. The types of calculations necessary to predict three-dimensional ice accretion is demonstrated. Results of calculations using 3-D method for a MS-317 swept wing geometry are projected onto a 2-D plane parallel to the free stream direction and compared to experimental results for the same geometry. It is anticipated that many modifications will be made to this approach, however, this effort will lay the groundwork for future modeling efforts. Results indicate that rime ice shapes indicate a difficulty in accurately calculating the ice shape in the runback region.

Potapczuk, M. G.↗

NASA NTRS1975

Relaxation methods in fluid mechanics

The present work considers the iterative solution of a coupled set of difference equations and examines methods that carry successive approximates to a state that is invariant with further iteration and independent of the initial guess. Methods are studied with regard to their efficiency and economy of computer resources. The basic principles of classical relaxation are set forth, with attention confined to linear elliptic equations. This discussion involves the evaluation of the spectral radius that is the magnitude of the eigenvalue with largest modulus. The subject of relaxation is then related to the study of ordinary differential equations and hyperbolic partial differential equations. Problems that occur when linearly dependent eigenvectors appear in the relaxation matrix are discussed, leading to multiply connected eigenvalues in the Jordan canonical form. Finally, a brief survey of relaxation methods used in aerodynamics is given.

Lomax, H.↗

NASA NTRS1994

Nucleation of Crystals From Solution in Microgravity (USML-1 Glovebox (GBX) Investigation)

A new method for initiating nucleation from solutions in microgravity which avoids nucleation on container walls and other surfaces is described. This method consists of injecting a small quantity of highly concentrated, heated solution into the interior of a lightly supersaturated, cooler host gowth solution. It was tested successfully on USML-I, producing a large number of LAP crystals whose longest dimension averaged 1 mm.

Kroes, Roger L.↗

NASA NTRS1974

Application of ultrasonic signature analysis for fatigue detection in complex structures

Ultrasonic signature analysis shows promise of being a singularly well-suited method for detecting fatigue in structures as complex as aircraft. The method employs instrumentation centered about a Fourier analyzer system, which features analog-to-digital conversion, digital data processing, and digital display of cross-correlation functions and cross-spectra. These features are essential to the analysis of ultrasonic signatures according to the procedure described here. In order to establish the feasibility of the method, the initial experiments were confined to simple plates with simulated and fatigue-induced defects respectively. In the first test the signature proved sensitive to the size of a small hole drilled into the plate. In the second test, performed on a series of fatigue-loaded plates, the signature proved capable of indicating both the initial appearance and subsequent growth of a fatigue crack. In view of these encouraging results it is concluded that the method has reached a sufficiently advanced stage of development to warrant application to small-scale structures or even actual aircraft.

Zuckerwar, A. J.↗

NASA NTRS2002

Control of Initialized Fractional-Order Systems

Due to the importance of historical effects in fractional-order systems, this paper presents a general fractional-order control theory that includes the time-varying initialization response. Previous studies have not properly accounted for these historical effects. The initialization response, along with the forced response, for fractional-order systems is determined. Stability properties of fractional-order systems are presented in the complex Airplane, which is a transformation of the s-plane. Time responses are discussed with respect to pole positions in the complex Airplane and frequency response behavior is included. A fractional-order vector space representation, which is a generalization of the state space concept, is presented including the initialization response. Control methods for vector representations of initialized fractional-order systems are shown. Nyquist, root-locus, and other input-output control methods are adapted to the control of fractional-order systems. Finally, the fractional-order differintegral is generalized to continuous order-distributions that have the possibility of including a continuum of fractional orders in a system element.

Hartly, Tom T.↗

NASA NTRS2002

Control of Initialized Fractional-Order Systems

Due to the importance of historical effects in fractional-order systems, this paper presents a general fractional-order control theory that includes the time-varying initialization response. Previous studies have not properly accounted for these historical effects. The initialization response, along with the forced response, for fractional-order systems is determined. Stability properties of fractional-order systems are presented in the complex w-plane, which is a transformation of the s-plane. Time responses are discussed with respect to pole positions in the complex w-plane and frequency response behavior is included. A fractional-order vector space representation, which is a generalization of the state space concept, is presented including the initialization response. Control methods for vector representations of initialized fractional-order systems are shown. Nyquist, root-locus, and other input-output control methods are adapted to the control of fractional-order systems. Finally, the fractional-order differintegral is generalized to continuous order-distributions that have the possibility of including a continuum of fractional orders in a system element.

Hartley, Tom T.↗

NASA NTRS2022

Microgreens Root Zone/Shoot Zone Partitioned Planting Box

To enable sustainable food production in future human exploration missions, plant growth is being studied by the Space Crop Production Team at KSC. Microgreens are good candidates for food supplements and contain specific nutrients that are lacking in the prepackaged diet, including vitamin C and vitamin K. Because they are densely sown, typical growth methods do not allow the ability to distinguish between the levels of evaporation from the rooting substrate and transpiration from the leaves. With larger plants, the root and shoot zone can be separated to distinguish these fluxes and accurate transpiration measurements of plant canopies are feasible. Furthermore, separation of the root and shoot zone may also be beneficial when harvesting microgreens in microgravity as it may reduce microbial contamination of the edible biomass by the roots, which have high microbial loads. Using a root and shoot separator box when harvesting may help with microbial contamination, but these tiny plants are challenging to handle in microgravity, so harvest management remains an open question. The innovation proposed here is a microgreen root/shoot partitioned planting box, which offers a solution to these challenges with accurate gas exchange measurements and a safe microgreen harvest in low gravity environments. Being able to measure transpiration of a microgreen canopy will be important for modeling plant growth in reduced gravity environments, so the first objective was to develop a planting unit with a seal between the root zone and the shoot zone. With an unsuitable harvesting technique, freshly harvested microgreens may add debris to the cabin, so the second objective was to test different harvesting techniques and management approaches associated with this innovation. These two objectives were pursued in parallel since many goals were the same: develop a planting unit that 1) separates the shoots from the roots, 2) allows acceptable germination rate, and 3) allows for seedlings to emerge and develop. What differed was the need to have a seal, which was only applicable for our gas exchange goal, and the need to have an embedded harvesting mechanism and bagging method which only applied to our harvesting goal. Testing of the various harvesting mechanisms and bagging methods was performed during a series of parabolic flights. All parabolic flight procedures took place inside a rented secondary containment chamber (e.g., glovebox) that was developed by the University of Louisville specifically for experiments involving fluids and other materials that may become airborne during reduced gravity flight. Three different harvesting methods and two different bagging collection methods were tested for microgreens. A third bagging method was initially tested but found to be unsuitable. Human factors were also taken into consideration, to identify which harvest and bagging collection methods would be easiest to use with favorable results in microgravity. Three parabolic flight tests were performed in total, one in November 2021 and two in December 2021.

Gioia Massa↗

NASA NTRS2020

Riley County Water Resources Project Summary - Comparing Runoff Curve Calculation Methods to Inform Local Resiliency Initiatives in Riley County, Kansas

Riley County, Kansas, has observed increased levels of flooding, potentially due to changes in land use/land cover (LULC) and seasonal vegetation variation. This study contrasts two methods of generating runoff curve numbers (CN) from 2006-2020. (1) The traditional Soil Conservation Service CN calculation method uses a look-up table and tracked LULC to determine runoff changes. These tables allow for land cover-specific CN and account for various farming techniques but lack flexibility in calculations for various seasons or plant health. (2) A dynamic method employs normalized difference vegetation index (NDVI) compiled over the rainy season each year to calculate CN using seasonal vegetation. This method allows for a more precise analysis of runoff variability within and between rainy seasons because it can be updated with greater temporal detail and captures higher spatial resolutions by using NDVI as a proxy for LULC. This study further uses inputs from the United States Geological Survey (USGS) National Land Cover Database (NLCD), the United States Department of Agriculture (USDA) Cropland Data Layer, and Landsat imagery to create more precise LULC raster datasets including both urban cover and crop-specific land use and curve number maps of the area. Results can guide decision makers in the City of Manhattan, Riley County Department of Planning and Development, Riley County Conservation District, the Kansas Forest Service, and the Kansas Department of Health and Environment toward informed decisions on resiliency strategies to address future flooding.

DEVELOP Project Summary↗