Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonlinear model reduction”

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 235 records · Page 13

A method for obtaining reduced-order control laws for high-order systems using optimization techniques

A method of synthesizing reduced-order optimal feedback control laws for a high-order system is developed. A nonlinear programming algorithm is employed to search for the control law design variables that minimize a performance index defined by a weighted sum of mean-square steady-state responses and control inputs. An analogy with the linear quadractic Gaussian solution is utilized to select a set of design variables and their initial values. To improve the stability margins of the system, an input-noise adjustment procedure is used in the design algorithm. The method is applied to the synthesis of an active flutter-suppression control law for a wind tunnel model of an aeroelastic wing. The reduced-order controller is compared with the corresponding full-order controller and found to provide nearly optimal performance. The performance of the present method appeared to be superior to that of two other control law order-reduction methods. It is concluded that by using the present algorithm, nearly optimal low-order control laws with good stability margins can be synthesized.

Mukhopadhyay, V.↗

New Uses for Sensitivity Analysis: How Different Movement Tasks Effect Limb Model Parameter Sensitivity

Original results for a newly developed eight-order nonlinear limb antagonistic muscle model of elbow flexion and extension are presented. A wider variety of sensitivity analysis techniques are used and a systematic protocol is established that shows how the different methods can be used efficiently to complement one another for maximum insight into model sensitivity. It is explicitly shown how the sensitivity of output behaviors to model parameters is a function of the controller input sequence, i.e., of the movement task. When the task is changed (for instance, from an input sequence that results in the usual fast movement task to a slower movement that may also involve external loading, etc.) the set of parameters with high sensitivity will in general also change. Such task-specific use of sensitivity analysis techniques identifies the set of parameters most important for a given task, and even suggests task-specific model reduction possibilities.

Winters, J. M.↗

A Computational and Experimental Study of Coflow Laminar Methane/Air Diffusion Flames: Effects of Fuel Dilution, Inlet Velocity, and Gravity

The influences of fuel dilution, inlet velocity, and gravity on the shape and structure of laminar coflow CH4-air diffusion flames were investigated computationally and experimentally. A series of nitrogen-diluted flames measured in the Structure and Liftoff in Combustion Experiment (SLICE) on board the International Space Station was assessed numerically under microgravity (mu g) and normal gravity (1g) conditions with CH4 mole fraction ranging from 0.4 to 1.0 and average inlet velocity ranging from 23 to 90 cm/s. Computationally, the MC-Smooth vorticity-velocity formulation was employed to describe the reactive gaseous mixture, and soot evolution was modeled by sectional aerosol equations. The governing equations and boundary conditions were discretized on a two-dimensional computational domain by finite differences, and the resulting set of fully coupled, strongly nonlinear equations was solved simultaneously at all points using a damped, modified Newton's method. Experimentally, flame shape and soot temperature were determined by flame emission images recorded by a digital color camera. Very good agreement between computation and measurement was obtained, and the conclusions were as follows. (1) Buoyant and nonbuoyant luminous flame lengths are proportional to the mass flow rate of the fuel mixture; computed and measured nonbuoyant flames are noticeably longer than their 1g counterparts; the effect of fuel dilution on flame shape (i.e., flame length and flame radius) is negligible when the flame shape is normalized by the methane flow rate. (2) Buoyancy-induced reduction of the flame radius through radially inward convection near the flame front is demonstrated. (3) Buoyant and nonbuoyant flame structure is mainly controlled by the fuel mass flow rate, and the effects from fuel dilution and inlet velocity are secondary.

flames↗

Global burned area increasingly explained by climate change

Fire behaviour is changing in many regions worldwide. However, nonlinear interactions between fire weather, fuel, land use, management and ignitions have impeded formal attribution of global burned area changes. Here, in this work, we demonstrate that climate change increasingly explains regional burned area patterns, using an ensemble of global fire models. The simulations show that climate change increased global burned area by 15.8% (95% confidence interval (CI) [13.1–18.7]) for 2003–2019 and increased the probability of experiencing months with above-average global burned area by 22% (95% CI [18–26]). In contrast, other human forcings contributed to lowering burned area by 19.1% (95% CI [21.9–15.8]) over the same period. Moreover, the contribution of climate change to burned area increased by 0.22% (95% CI [0.22–0.24]) per year globally, with the largest increase in central Australia. Our results highlight the importance of immediate, drastic and sustained GHG emission reductions along with landscape and fire management strategies to stabilize fire impacts on lives, livelihoods and ecosystems.

54 ENVIRONMENTAL SCIENCES↗

Atlantic Tropical Cyclone Monitoring with AMSU-A: Estimation of Maximum Sustained Wind Speeds

The first Advanced Microwave Sounding Unit temperature sounder (AMSU-A) was launched on the NOAA-15 satellite on 13 May 1998. The AMSU-A's higher spatial and radiometric resolutions provide more useful information on the strength of the middle- and upper-tropospheric warm cores associated with tropical cyclones than have previous microwave temperature sounders. The gradient wind relationship suggests that the temperature gradient near the core of tropical cyclones increases nonlinearly with wind speed. The gradient wind equation is recast to include AMSU-A-derived variables, Stepwise regression is used to determine which of these variables is most closely related to maximum sustained winds (V(sub max)). The satellite variables investigated include the radially averaged gradients at two spatial resolutions of AMSU-A channels 1-10 T(sub b) data (delta(sub r)T(sub B)), the squares of these gradients, a channel-15-based scattering index (SI(sub 89)), and area-averaged T(sub B). Calculations of T(sub B) and delta(sub r)T(sub B) from mesoscale model simulations of Andrew reveal the effects of the AMSU spatial sampling on the cyclone warm core presentation. Stepwise regression of 66 AMSU-A terms against National Hurricane Center V(sub max) estimates from the 1998 and 1999 Atlantic hurricane season confirms the existence of a nonlinear relationship between wind speed and radially averaged temperature gradients near the cyclone warm core. Of six regression terms, four are dominated by temperature information, and two are interpreted as correcting for hydrometeor contamination. Jackknifed regressions were performed to estimate the algorithm performance on independent data. For the 82 cases that had in situ measurements of V(sub max), the average error standard deviation was 4.7 m/s. For 108 cases without in situ wind data, the average error standard deviation was 7.5 m/s Operational considerations, including the detection of weak cyclones and false alarm reduction, are also discussed.

Spencer, Roy W.↗

Atlantic Tropical Cyclone Monitoring with AMSU-A: Estimation of Maximum Sustained Wind Speeds

The first Advanced Microwave Sounding Unit temperature sounder (AMSU-A) was launched on the NOAA-15 satellite on 13 May 1998. The AMSU-A's higher spatial and radiometric resolutions provide more useful information on the strength of the middle and upper tropospheric warm cores associated with tropical cyclones than have previous microwave temperature sounders. The gradient wind relationship suggests that the temperature gradient near the core of tropical cyclones increases nonlinearly with wind speed. We recast the gradient wind equation to include AMSU-A derived variables. Stepwise regression is used to determine which of these variables is most closely related to maximum sustained winds (V(sub max)). The satellite variables investigated include the radially averaged gradients at two spatial resolutions of AMSU-A channels 1 through 10 T(sub b) data (delta(sub r)T(sub b)), the squares of these gradients, a channel 15 based scattering index (SI-89), and area averaged T(sub b). Calculations of Tb and delta(sub r)T(sub b) from mesoscale model simulations of Andrew reveal the effects of the AMSU spatial sampling on the cyclone warm core presentation. Stepwise regression of 66 AMSU-A terms against National Hurricane Center (NHC) V(sub max) estimates from the 1998 and 1999 Atlantic hurricane season confirms the existence of a nonlinear relationship between wind speed and radially averaged temperature gradients near the cyclone warm core. Of six regression terms, four are dominated by temperature information, and two are interpreted as correcting for hydrometeor contamination. Jackknifed regressions were performed to estimate the algorithm performance on independent data. For the 82 cases that had in situ measurements of V(sub max), the average error standard deviation was 4.7 m/s. For 108 cases without in situ wind data, the average error standard deviation was 7.5 m/s. Operational considerations, including the detection of weak cyclones and false alarm reduction are also discussed.

Spencer, Roy↗

Nonlinear Buckling Analysis of Tow-Steered Composite Cylinders with Cutouts

The buckling and postbuckling behavior of two composite tow-steered shells with cutouts of different sizes is assessed using nonlinear finite element (FE) analysis and compared to experimental measurements. The cylindrical shells are manufactured using an automated fiber placement system, where the shells’ fiber orientation angles vary continuously around the shell circumference from ±10 degrees on the crown and keel to ±45 degrees on the sides. One shell features thickness variations due to tow overlaps that result from application of all 24 tows during each pass of the fiber placement system. The second shell uses the system’s tow drop/add capability to achieve a more uniform wall thickness without overlaps. Unreinforced cutouts of two different sizes—the first smaller cutout representing a passenger door on a commercial aircraft and the second larger cutout a cargo door—were machined into each of the two cylinders resulting in a total of four test cases. These cylinders were tested in axial compression and buckled elastically in previous work and are now analyzed using nonlinear FE models to compare bifurcation buckling loads as well as the load-displacement response in the prebuckling and postbuckling regimes. For all four shells analyzed, the prebuckling stiffness, buckling load, and deformation mode sequence throughout the loading-unloading cycle is accurately reproduced by the models. In particular, the shells first buckle locally around the cutouts in a stable (super-critical) manner with only a slight decrease in axial stiffness, which occurs due to the favorable load redistribution facilitated by tow steering. The shells then buckle globally in an unstable (sub-critical) manner with diamond-shaped buckles forming to the left and right of the cutouts. The buckling load of all shells with cutouts is at least 82% of the buckling load of the pristine shells without cutouts. Overall, the ability to sustain local buckling phenomena, and the relatively small reductions in global buckling load compared to pristine shells without cutouts, demonstrates the great potential of using tow steering to mitigate the adverse effects of cutouts in axially-compressed shell structures.

composites↗

Nonlinear Buckling Analysis of Tow-Steered Composite Cylinders with Cutouts

The buckling and postbuckling behavior of two composite tow-steered shells with cutouts of different sizes is assessed using nonlinear finite element (FE) analysis and compared to experimental measurements. The cylindrical shells are manufactured using an automated fiber placement system, where the shells’ fiber orientation angles vary continuously around the shell circumference from 10 degrees on the crown and keel to 45 degrees on the sides. One shell features thickness variations due to tow overlaps that result from application of all 24 tows during each pass of the fiber placement system. The second shell uses the system’s tow drop/add capability to achieve a more uniform wall thickness without overlaps. Unreinforced cutouts of two different sizes—the first smaller cutout representing a passenger door on a commercial aircraft and the second larger cutout a cargo door—were machined into each of the two cylinders resulting in a total of four test cases. These cylinders were tested in axial compression and buckled elastically in previous work and are now analyzed using nonlinear FE models to compare bifurcation buckling loads as well as the load-displacement response in the prebuckling and postbuckling regimes. For all four shells analyzed, the prebuckling stiffness, buckling load, and deformation mode sequence throughout the loading-unloading cycle is accurately reproduced by the models. In particular, the shells first buckle locally around the cutouts in a stable (super-critical) manner with only a slight decrease in axial stiffness, which occurs due to the favorable load redistribution facilitated by tow steering. The shells then buckle globally in an unstable (sub-critical) manner with diamond-shaped buckles forming to the left and right of the cutouts. The buckling load of all shells with cutouts is at least 82% of the buckling load of the pristine shells without cutouts. Overall, the ability to sustain local buckling phenomena, and the relatively small reductions in global buckling load compared to pristine shells without cutouts, demonstrates the great potential of using tow steering to mitigate the adverse effects of cutouts in axially-compressed shell structures.

Buckling↗

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

Sustained advances in the mathematics of modeling and simulation have resulted in the capability today for routine simulation of a number of large scale complex DOE-relevant systems. As remarkable as this capability for solving the so-called forward problem is, it is typically only the first step-an inner loop within an outer loop that explores the simulation model's parameter space and decision space to characterize uncertainty in the model's predictions, learn unknown model parameters from data, design the most informative experiments, determine optimal control strategies, and create optimal designs. Broadly, what unifies all of these outer loop problems is that they are, in one form or another, optimization problems over parameter/control/design space that are constrained by complex uncertain models. To fully realize the power of scientific simulation as a basis for scientific discovery, technological innovation, and rational decision-making, it is imperative to move beyond simulation to tackle the outer loop of optimization for learning from data, experimental design, and control with complex uncertain models. When the models under consideration are large-scale and complex, and when the optimization variable and uncertain parameter spaces are high (or infinite) dimensional, this constitutes a grand challenge of the highest order, and is intractable with conventional methods. To overcome these challenges, the AEOLUS Center was established to develop a unified mathematical, computational, and statistical framework for (1) Learning predictive models from complex data via Bayesian inference and optimization, and (2) Optimizing experiments, processes, and designs using the resulting uncertain models. These problems are intractable with conventional methods, for several reasons: (1) The simulation problems that govern the inner loops of the optimization problems are expensive to execute (due to severe nonlinearity, heterogeneity, multiphysics/multiscale coupling); (2) The optimization variable and uncertain parameter spaces are high dimensional, often stemming from discretizations of infinite dimensional fields such as initial conditions, sources, or material properties. We argue that the key to overcoming these challenges is to develop new mathematical, computational, and statistical methods that exploit the structure of the Bayesian inference and optimization problems mediated by their underlying complex uncertain models. This structure includes the regularity, sparsity, geometry, low intrinsic dimensionality, and multifidelity nature of the maps from uncertain parameter/optimization variable spaces to the specific objectives targeted: Bayesian inference, optimal experimental design, and optimal control design. Black box methods developed as generic tools are incapable of exploiting this structure. To be successful, we must create, integrate, and cross-fertilize ideas across multiple areas of applied math--including approximation theory, Bayesian inference, data science, experimental design, information theory, machine learning, model reduction, optimal control theory, parallel algorithms, PDE-constrained optimization, randomized algorithms, stochastic optimization, and uncertainty quantification--all while exploiting the structure of the problems at hand. With this goal in mind, we have marshaled a team of leading authorities in these areas. While the methods we develop will be broadly applicable across a wide spectrum of DOE problems in which experiments inform models and the systems those models describe must be optimized under uncertainty, we have chosen a specific area, advanced manufacturing and materials, to drive our work. AMM is characterized by complex models across multiple scales, and is a rich source of challenging problems in inference, experimental design, and optimal control, requiring multifaceted and integrated advances in applied mathematics. As such, AMM serves as an excellent vehicle to motivate and demonstrate the advances in applied mathematics developed by our center.

97 MATHEMATICS AND COMPUTING↗

Impact of Color Space and Color Resolution on Vehicle Recognition Models

In this study, we analyze both linear and nonlinear color mappings by training on versions of a curated dataset collected in a controlled campus environment. We experiment with color space and color resolution to assess model performance in vehicle recognition tasks. Color encodings can be designed in principle to highlight certain vehicle characteristics or compensate for lighting differences when assessing potential matches to previously encountered objects. The dataset used in this work includes imagery gathered under diverse environmental conditions, including daytime and nighttime lighting. Experimental results inform expectations for possible improvements with automatic color space selection through feature learning. Moreover, we find there is only a gradual decrease in model performance with degraded color resolution, which suggests the need for simplified data collection and processing. By focusing on the most critical features, we could see improved model generalization and robustness, as the model becomes less prone to overfitting to noise or irrelevant details in the data. Such a reduction in resolution will lower computational complexity, leading to quicker training and inference times.

47 OTHER INSTRUMENTATION↗

Acoustic Analysis of Sensor Ports

Acquiring data in rocket engines that is representative of the actual dynamic environment can often be difficult due to a multitude of influences. One source of contamination that is often not considered entirely is the response associated with the acoustic cavity created by a sensor offset. It is common to offset a sensor due to various reasons such as for mounting accessibility, thermal isolation, shock reduction, or prevention of debris impingement. While estimating the natural frequency of the acoustic cavity is straightforward, limited analysis has been described on the determination of the overall frequency response. The sensor port design approach usually attempts to ensure the port is short enough such that the acoustic response is negligible near the frequency-of-interest, but this requires knowledge of the frequency response and simple rules-of-thumb are not always guaranteed. Data correction and/or data interpretation is also often desired for an unsatisfactory response. The limited response analysis in the literature only offers approximations or neglects important contributions. A new approach is devised theoretically and computationally that captures the true acoustic response of a sensor port. This paper summarizes the acoustics background, the port response theoretical development, and provides comparisons of a port acoustic response using an analytical model and computational acoustics. The effects of nonlinear acoustics are also examined. Additionally, the paper summarizes the design of a specialized filter using the predicted sensor port response that can be applied to data for correction.

Casiano, M. J.↗

Predicting nonequilibrium Green’s function dynamics and photoemission spectra via nonlinear integral operator learning

Understanding the dynamics of nonequilibrium quantum many-body systems is an important research topic in a wide range of fields across condensed matter physics, quantum optics, and high-energy physics. However, numerical studies of large-scale nonequilibrium phenomena in realistic materials face serious challenges due to intrinsic high-dimensionality of quantum many-body problems and the absence of time-invariance. The nonequilibrium properties of many-body systems can be described by the dynamics of the correlator, or the Green's function of the system, whose time evolution is given by a high-dimensional system of integro-differential equations, known as the Kadanoff–Baym equations (KBEs). The time-convolution term in KBEs, which needs to be recalculated at each time step, makes it difficult to perform long-time numerical simulation. In this paper, we develop an operator-learning framework based on recurrent neural networks (RNNs) to address this challenge. We utilize RNNs to learn the nonlinear mapping between Green's functions and convolution integrals in KBEs. By using the learned operators as a surrogate model in the KBE solver, we obtain a general machine-learning scheme for predicting the dynamics of nonequilibrium Green's functions. Besides significant savings per each time step, the new methodology reduces the temporal computational complexity from $O(N_t^3)$ to $O(N_t)$ where N t is the number of steps taken in a simulation, thereby making it possible to study large many-body problems which are currently infeasible with conventional KBE solvers. Through various numerical examples, we demonstrate the effectiveness of the operator-learning based approach in providing accurate predictions of physical observables such as the reduced density matrix and time-resolved photoemission spectra. Moreover, our framework exhibits clear numerical convergence and can be easily parallelized, thereby facilitating many possible further developments and applications.

97 MATHEMATICS AND COMPUTING↗

Scalable edge clustering of dynamic graphs via weighted line graphs

Timestamped relational datasets consisting of records (or connections) between pairs of entities are ubiquitous in network science. For applications like peer-to-peer communication, email, various social network interactions, and computer network security, it is useful to organize these records into groups based on how and when they are occurring. Weighted line graphs offer a natural way to model how records are related in such datasets but for large real-world graph topologies, building and utilizing the line graph is prohibitively expensive. Here, we present the framework to cluster the edges of a dynamic graph via the associated line graph that contains two major contributions. The first is a method to work with the line graph implicitly and the second is a distributed scale implementation of an agglomerative hierarchical graph clustering algorithm. We outline a novel hierarchical dynamic graph edge clustering approach that efficiently breaks massive relational datasets into small sets of edges containing events at various timescales. This is in stark contrast to traditional graph clustering algorithms that prioritize highly connected (clique-like) community structures. Our approach relies on constructing a sufficient subgraph of a weighted line graph and applying a hierarchical agglomerative clustering. This approach is related to scalable techniques from spatial clustering, nonlinear-dimension reduction, topological data analysis, and draws particular inspiration from HDBSCAN. As an edge clustering, this method yields an overlapping node clustering. Our algorithm is parallelizable and we demonstrate efficient clustering of a billion-scale, real-world dynamic graph into small edge sets that correlate in topology and time. The entire clustering process for a graph with tens of billions of edges takes just a few minutes of run time on 256 nodes of a distributed compute environment. We argue how the output of the edge clustering is useful for a multitude of data visualization and powerful machine learning tasks, both involving the original massive dynamic graph data and metadata associated with the nodes and edges. Finally, we describe how this approach can be extended to dynamic hypergraphs and dynamic graphs/hypergraphs with unstructured data living on vertices and edges.

Data Analysis↗

Initial Uncertainty Analysis of Carbon Tetrachloride Contamination and Remediation in the Ringold A and Lower Mud Units at the Central Plateau

The long-term effectiveness of groundwater cleanup at the Hanford Site Central Plateau depends on predictive models that can capture key uncertainties in contaminant fate and transport. Carbon tetrachloride (CCl 4 ), a persistent and toxic compound, presents particular challenges due to variability in degradation rates, uncertainty in initial plume distribution, and subsurface heterogeneity. These uncertainties directly influence plume persistence, migration pathways, and remedy performance, and thus must be systematically evaluated to support long-term remediation planning. To address these gaps, a large-scale Monte Carlo analysis was conducted using the Plateau to River (P2R) model framework. The modeling approach parameterized three primary uncertainty factors: (1) degradation rate, (2) initial plume distribution, and (3) hydraulic conductivity. Degradation was represented as a first-order process, with half-lives ranging from 70 to 700 years. Initial plume distributions were created using a geostatistical simulation method (sgsim), which generates many equally plausible versions of how contaminants might be distributed underground. From this, 100 different scenarios were mapped onto the P2R grid. Variability in hydraulic conductivity was represented in a similar way, with 100 scenarios each for the Ringold Lower Mud and Ringold A units (layers 6 and 7), based on fitted exponential variograms and conditioned to well data. In total, more than 1000 realizations were simulated to assess plume behavior under uncertainty. Results demonstrate that degradation kinetics exert the strongest control over plume persistence: Shorter half-lives produced rapid mass reduction, while longer half-lives yielded persistent plumes with limited attenuation. A nonlinear response was observed, with steep mass reductions at half-lives greater than 200 years and near-linear declines beyond this threshold, reflecting interactions between degradation and pumping. The initial plume distribution strongly influenced early transport patterns, with broader sources generating larger plume footprints, although pump-and-treat operations constrained plume migration to managed areas. By comparison, hydraulic conductivity variability in the Ringold units had only a secondary influence, modifying spreading behavior without altering the dominant migration pathways governed by source configuration and hydraulic controls. Overall, the analysis highlights that uncertainty in degradation rate and initial plume configuration are the primary drivers of variability in plume predictions, while conductivity heterogeneity plays a limited role. These findings underscore the need for improved site-specific data on degradation processes and source characterization to enhance the reliability of long-term performance assessments and to better inform remedial decision-making at the Central Plateau.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Nonlinear dynamics of global atmospheric and Earth-system processes

General Circulation Model (GCM) studies of the atmospheric response to change boundary conditions are discussed. Results are reported on an extensive series of numerical studies based on the National Center for Atmospheric Research (NCAR) Community Climate Model (CCM) general circulation model. In these studies the authors determined the response to systematic changes in atmospheric CO2 ranging from 100 to 1000 ppm; to changes in the prescribed sea surface temperature (SST) in the Gulf of Mexico, such as occurred during the deglaciation phase of the last ice age; to changes in soil moisture over North America; and to changes in sea ice extent in the Southern Hemisphere. Study results show that the response of surface temperature and other variables is nearly logarithmic, with lower levels of CO2 implying greater sensitivity of the atmospheric state to changes in CO2. It was found that the surface temperature of the Gulf of Mexico exerts considerable control over the storm track and behavior of storm systems over the North Atlantic through its influence on evaporation and the source of latent heat. It was found that reductions in soil moisture can play a significant role in amplifying and maintaining North American drought, particularly when a negative soil moisture anomaly prevails late in the spring.

Saltzman, Barry↗

Landslide Hazard is Projected to Increase Across High Mountain Asia

High Mountain Asia has long been known as a hotspot for landslide risk, and studies have suggested that landslide hazard is likely to increase in this region over the coming decades. Extreme precipitation may become more frequent, with a nonlinear response relative to increasing global temperatures. However, these changes are geographically varied. This article maps probable changes to landslide hazard, as shown by a landslide hazard indicator (LHI) derived from downscaled precipitation and temperature. In order to capture the nonlinear response of slopes to extreme precipitation, a simple machine-learning model was trained on a database of landslides across High Mountain Asia to develop a regional LHI. This model was applied to statistically downscaled data from the 30 members of the Seamless System for Prediction and Earth System Research large ensembles to produce a range of possible outcomes under the Shared Socioeconomic Pathways 2-4.5 and 5-8.5. The LHI reveals that landslide hazard will increase in most parts of High Mountain Asia. Absolute increases will be highest in already hazardous areas such as the Central Himalaya, but relative change is greatest on the Tibetan Plateau. Even in regions where landslide hazard declines by year 2100, it will increase prior to the mid-century mark. However, the seasonal cycle of landslide occurrence will not change greatly across High Mountain Asia. Although substantial uncertainty remains in these projections, the overall direction of change seems reliable. These findings highlight the importance of continued analysis to inform disaster risk reduction strategies for stakeholders across High Mountain Asia.

Thomas A Stanley↗

Human influence on joint changes in temperature, rainfall and continental aridity

Despite the pervasive impact of drought on human and natural systems, the large-scale mechanisms conducive to regional drying remain poorly understood. Here we use a multivariate approach to identify two distinct externally forced fingerprints from multiple ensembles of Earth system model simulations. The leading fingerprint, F(M1)(x), is characterized by global warming, intensified wet–dry patterns and progressive large-scale continental aridification, largely driven by multidecadal increases in greenhouse gas (GHG) emissions. The second fingerprint, F(M2)(x), captures a pronounced interhemispheric temperature contrast, associated meridional shifts in the intertropical convergence zone and correlated anomalies in precipitation and aridity over California, the Sahel and India. F(M2)(x) exhibits nonlinear temporal behaviour: the intertropical convergence zone moves southwards before 1975 in response to increases in hemispherically asymmetric sulfate aerosol emissions, and it shifts northwards after 1975 due to reduced sulfur dioxide emissions and the GHG-induced warming of Northern Hemisphere landmasses. Both fingerprints are statistically identifiable in observations of joint changes in temperature, rainfall and aridity during 1950–2014. We show that the reliable simulation of these changes requires combined forcing by GHGs, direct and indirect effects of aerosols, and large volcanic eruptions. Our results suggest that GHG-induced aridification may be modulated regionally by future reductions in sulfate aerosol emissions.

Human influence↗

Optimal guidance law development for an advanced launch system

The objective of this research effort was to develop a real-time guidance approach for launch vehicles ascent to orbit injection. Various analytical approaches combined with a variety of model order and model complexity reduction have been investigated. Singular perturbation methods were first attempted and found to be unsatisfactory. The second approach based on regular perturbation analysis was subsequently investigated. It also fails because the aerodynamic effects (ignored in the zero order solution) are too large to be treated as perturbations. Therefore, the study demonstrates that perturbation methods alone (both regular and singular perturbations) are inadequate for use in developing a guidance algorithm for the atmospheric flight phase of a launch vehicle. During a second phase of the research effort, a hybrid analytic/numerical approach was developed and evaluated. The approach combines the numerical methods of collocation and the analytical method of regular perturbations. The concept of choosing intelligent interpolating functions is also introduced. Regular perturbation analysis allows the use of a crude representation for the collocation solution, and intelligent interpolating functions further reduce the number of elements without sacrificing the approximation accuracy. As a result, the combined method forms a powerful tool for solving real-time optimal control problems. Details of the approach are illustrated in a fourth order nonlinear example. The hybrid approach is then applied to the launch vehicle problem. The collocation solution is derived from a bilinear tangent steering law, and results in a guidance solution for the entire flight regime that includes both atmospheric and exoatmospheric flight phases.

Calise, Anthony J.↗