A preconditioned numerical solver for stiff nonlinear reaction-diffusion equations with fractional Laplacians that avoids dense matrices

The numerical solution of fractional partial differential equations poses significant computational challenges in regard to efficiency as a result of the spatial nonlocality of the fractional differential operators. The dense coefficient matrices that arise from spatial discretisation of these operators mean that even one-dimensional problems can be difficult to solve using standard methods on grids comprising thousands of nodes or more. In this work we address this issue of efficiency for one-dimensional, nonlinear space-fractional reaction–diffusion equations with fractional Laplacian operators. We apply variable-order, variable-stepsize backward differentiation formulas in a Jacobian-free Newton–Krylov framework to advance the solution in time. A key advantage of this approach is the elimination of any requirement to form the dense matrix representation of the fractional Laplacian operator. We show how a banded approximation to this matrix, which can be formed and factorised efficiently, can be used as part of an effective preconditioner that accelerates convergence of the Krylov subspace iterative solver. Our approach also captures the full contribution from the nonlinear reaction term in the preconditioner, which is crucial for problems that exhibit stiff reactions. Numerical examples are presented to illustrate the overall effectiveness of the solver.

Quantification of stability in an agility drill using linear and nonlinear measures of variability

This study implemented linear and nonlinear methods of measuring variability to determine differences in stability of two groups of skilled (n = 10) and unskilled (n = 10) participants performing 3m forward/backward shuttle agility drill. We also determined whether stability measures differed between the forward and backward segments of the drill. Finally, we sought to investigate whether local dynamic stability, measured using largest finite-time Lyapunov exponents, changed from distal to proximal lower extremity segments. Three-dimensional coordinates of five lower extremity markers data were recorded. Results revealed that the Lyapunov exponents were lower (P < 0.05) for skilled participants at all joint markers indicative of higher levels of local dynamic stability. Additionally, stability of motion did not differ between forward and backward segments of the drill (P > 0.05), signifying that almost the same control strategy was used in forward and backward directions by all participants, regardless of skill level. Furthermore, local dynamic stability increased from distal to proximal joints (P < 0.05) indicating that stability of proximal segments are prioritized by the neuromuscular control system. Finally, skilled participants displayed greater foot placement standard deviation values (P < 0.05), indicative of adaptation to task constraints. The results of this study provide new methods for sport scientists, coaches to characterize stability in agility drill performance.

Linear and weakly nonlinear stability analyses of cooperative car-following models

Stability analyses have been widely used to better understand the mechanism of traffic jam formation. In this paper, we consider the impact of cooperative systems (a.k.a. connected vehicles) on traffic dynamics and, more precisely, on flow stability. Cooperative systems are emerging technologies enabling communication between vehicles and/or with the infrastructure. In a distributed communication framework, equipped vehicles are able to send and receive information to/from other equipped vehicles. Here, the effects of cooperative traffic are modeled through a general bilateral multianticipative car-following law that improves cooperative drivers' perception of their surrounding traffic conditions within a given communication range. Linear stability analyses are performed for a broad class of car-following models. They point out different stability conditions in both multianticipative and nonmultianticipative situations. To better understand what happens in unstable conditions, information on the shock wave structure is studied in the weakly nonlinear regime by the mean of the reductive perturbation method. The shock wave equation is obtained for generic car-following models by deriving the Korteweg de Vries equations. We then derive traffic-state-dependent conditions for the sign of the solitary wave (soliton) amplitude. This analytical result is verified through simulations. Simulation results confirm the validity of the speed estimate. The variation of the soliton amplitude as a function of the communication range is provided. The performed linear and weakly nonlinear analyses help justify the potential benefits of vehicle-integrated communication systems and provide new insights supporting the future implementation of cooperative systems.

Simulated systems optimization through response surface techniques [Optimización de sistemas simulados a través de técnicas de superficie de respuesta]

The purpose of this article is to show the applicability and benefits of the techniques of design of experiments as an optimization tool for discrete simulation models. The simulated systems are computational representations of real-life systems; its characteristics include a constant evolution that follows the occurrence of discrete events along the time. In this study, a production system, designed with the business philosophy JIT (Just in Time) is used, which seeks to achieve excellence in organizations through waste reduction in all the operational aspects. The most typical tool of JIT systems is the KANBAN production control that seeks to synchronize demand with flow of materials, minimize work in process, and define production metrics. Using experimental design techniques for stochastic optimization, the impact of the operational factors on the efficiency of the KANBAN / CONWIP simulation model is analyzed. The results show the effectiveness of the integration of experimental design techniques and discrete simulation models in the calculation of the operational parameters. Furthermore, the reliability of the methodologies found was improved with a new statistical consideration.

A three-dimensional hybrid smoothed finite element method (H-SFEM) for nonlinear solid mechanics problems

This paper presents a novel three-dimensional hybrid smoothed finite element method (H-SFEM) for solid mechanics problems. In 3D H-SFEM, the strain field is assumed to be the weighted average between compatible strains from the finite element method (FEM) and smoothed strains from the node-based smoothed FEM with a parameter α equipped into H-SFEM. By adjusting α, the upper and lower bound solutions in the strain energy norm and eigenfrequencies can always be obtained. The optimized α value in 3D H-SFEM using a tetrahedron mesh possesses a close-to-exact stiffness of the continuous system, and produces ultra-accurate solutions in terms of displacement, strain energy and eigenfrequencies in the linear and nonlinear problems. The novel domain-based selective scheme is proposed leading to a combined selective H-SFEM model that is immune from volumetric locking and hence works well for nearly incompressible materials. The proposed 3D H-SFEM is an innovative and unique numerical method with its distinct features, which has great potential in the successful application for solid mechanics problems.

Study on tracheal collapsibility, compliance, and stress by considering nonlinear mechanical property of cartilage

Tracheal cartilage has been widely regarded as a linear elastic material either in experimental studies or in analytic and numerical models. However, it has been recently demonstrated that, like other fiber-oriented biological tissues, tracheal cartilage is a nonlinear material, which displays higher strength in compression than in extension. Considering the nonlinearity requires a more complex theoretical frame work and costs more to simulate. This study aims to quantify the deviation due to the simplified treatment of the tracheal cartilage as a linear material. It also evaluates the improved accuracy gained by considering the nonlinearity. Pig tracheal rings were used to exam the mechanical properties of cartilage and muscular membrane. By taking into account the asymmetric shape of tracheal cartilage, the collapse behavior of complete rings was simulated, and the compliance of airway and stress in the muscular membrane were discussed. The results obtained were compared with those assuming linear mechanical properties. The following results were found: (1) Models based on both types of material properties give a small difference in representing collapse behavior; (2) regarding compliance, the relative difference is big, ranging from 10 to 40% under negative pressure conditions; and (3) the difference in determining stress in the muscular membrane is small too: <5%. In conclusion, treating tracheal cartilage as a linear material will not cause big deviations in representing the collapse behavior, and mechanical stress in the muscular part, but it will induce a big deviation in predicting the compliance, particularly when the transmural pressure is lower than -0.5 kPa. The results obtained in this study may be useful in both understanding the collapse behavior of trachea and in evaluating the error induced by the simplification of treating the tracheal cartilage as a linear elastic material.

Nonlinear mechanical property of tracheal cartilage: A theoretical and experimental study

Background: Despite being the stiffest airway of the bronchial tree, the trachea undergoes significant deformation due to intrathoracic pressure during breathing. The mechanical properties of the trachea affect the flow in the airway and may contribute to the biological function of the lung. Method: A Fung-type strain energy density function was used to investigate the nonlinear mechanical behavior of tracheal cartilage. A bending test on pig tracheal cartilage was performed and a mathematical model for analyzing the deformation of tracheal cartilage was developed. The constants included in the strain energy density function were determined by fitting the experimental data. Result: The experimental data show that tracheal cartilage is a nonlinear material displaying higher strength in compression than in tension. When the compression forces varied from -0.02 to -0.03 N and from -0.03 to -0.04 N, the deformation ratios were 11.03±2.18% and 7.27±1.59%, respectively. Both were much smaller than the deformation ratios (20.01±4.49%) under tension forces of 0.02 to 0.01 N. The Fung-type strain energy density function can capture this nonlinear behavior very well, whilst the linear stress-strain relation cannot. It underestimates the stability of trachea by exaggerating the displacement in compression. This study may improve our understanding of the nonlinear behavior of tracheal cartilage and it may be useful for the future study on tracheal collapse behavior under physiological and pathological conditions.

Nonlinear analysis for the pre- and post-yield behaviour of a composite structure with the refined plastic hinge approach

The computational technique of the full ranges of the second-order inelastic behaviour evaluation of steel-concrete composite structure is not always sought forgivingly, and therefore it hinders the development and application of the performance-based design approach for the composite structure. To this end, this paper addresses of the advanced computational technique of the higher-order element with the refined plastic hinges to capture the all-ranges behaviour of an entire steel-concrete composite structure. Moreover, this paper presents the efficient and economical cross-section analysis to evaluate the element section capacity of the non-uniform and arbitrary composite section subjected to the axial and bending interaction. Based on the same single algorithm, it can accurately and effectively evaluate nearly continuous interaction capacity curve from decompression to pure bending technically, which is the important capacity range but highly nonlinear. Hence, this cross-section analysis provides the simple but unique algorithm for the design approach. In summary, the present nonlinear computational technique can simulate both material and geometric nonlinearities of the composite structure in the accurate, efficient and reliable fashion, including partial shear connection and gradual yielding at pre-yield stage, plasticity and strain-hardening effect due to axial and bending interaction at post-yield stage, loading redistribution, second-order P-δ and P-Δ effect, and also the stiffness and strength deterioration. And because of its reliable and accurate behavioural evaluation, the present technique can be extended for the design of the high-strength composite structure and potentially for the fibre-reinforced concrete structure.

Optimization models for expanding a railway’s theoretical capacity

Changing the topology of a railway network can greatly affect its capacity. Railway networks however can be altered in a multitude of different ways. As each way has significant immediate and long term financial ramifications, it is a difficult task to decide how and where to expand the network. In response some railway capacity expansion models (RCEM) have been developed to help capacity planning activities, and to remove physical bottlenecks in the current railway system. The exact purpose of these models is to decide given a fixed budget, where track duplications and track sub divisions should be made, in order to increase theoretical capacity most. These models are high level and strategic, and this is why increases to the theoretical capacity is concentrated upon. The optimization models have been applied to a case study to demonstrate their application and their worth. The case study evidently shows how automated approaches of this nature could be a formidable alternative to current manual planning techniques and simulation. If the exact effect of track duplications and sub-divisions can be sufficiently approximated, this approach will be very applicable.

The moderating effects of job control and selection, optimization, and compensation strategies on the age-work ability relationship

Work ability describes employees' capability to carry out their work with respect to physical and psychological job demands. This study investigated direct and interactive effects of age, job control, and the use of successful aging strategies called selection, optimization, and compensation (SOC) in predicting work ability. We assessed SOC strategies and job control by using employee self-reports, and we measured employees' work ability using supervisor ratings. Data collected from 173 health-care employees showed that job control was positively associated with work ability. Additionally, we found a three-way interaction effect of age, job control, and use of SOC strategies on work ability. Specifically, the negative relationship between age and work ability was weakest for employees with high job control and high use of SOC strategies. These results suggest that the use of successful aging strategies and enhanced control at work are conducive to maintaining the work ability of aging employees. We discuss theoretical and practical implications regarding the beneficial role of the use of SOC strategies utilized by older employees and enhanced contextual resources at work for aging employees.

The buffering effect of selection, optimization, and compensation strategy use on the relationship between problem solving demands and occupational well-being: A daily diary study

This study investigated within-person relationships between daily problem solving demands, selection, optimization, and compensation (SOC) strategy use, job satisfaction, and fatigue at work. Based on conservation of resources theory, it was hypothesized that high SOC strategy use boosts the positive relationship between problem solving demands and job satisfaction, and buffers the positive relationship between problem solving demands and fatigue. Using a daily diary study design, data were collected from 64 administrative employees who completed a general questionnaire and two daily online questionnaires over four work days. Multilevel analyses showed that problem solving demands were positively related to fatigue, but unrelated to job satisfaction. SOC strategy use was positively related to job satisfaction, but unrelated to fatigue. A buffering effect of high SOC strategy use on the demands-fatigue relationship was found, but no booster effect on the demands-satisfaction relationship. The results suggest that high SOC strategy use is a resource that protects employees from the negative effects of high problem solving demands.

Maintaining a focus on opportunities at work: The interplay between age, job complexity, and the use of selection, optimization, and compensation strategies

The concept of focus on opportunities describes how many new goals, options, and possibilities employees believe to have in their personal future at work. This study investigated the specific and shared effects of age, job complexity, and the use of successful aging strategies called selection, optimization, and compensation (SOC) in predicting focus on opportunities. Results of data collected from 133 employees of one company (mean age = 38 years, SD = 13, range 16–65 years) showed that age was negatively, and job complexity and use of SOC strategies were positively related to focus on opportunities. In addition, older employees in high-complexity jobs and older employees in low-complexity jobs with high use of SOC strategies were better able to maintain a focus on opportunities than older employees in low-complexity jobs with low use of SOC strategies.

Selection, optimization, and compensation strategies: Interactive effects on daily work engagement

The theory of selective optimization with compensation (SOC) proposes that the “orchestrated” use of three distinct action regulation strategies (selection, optimization, and compensation) leads to positive employee outcomes. Previous research examined overall scores and additive models (i.e., main effects) of SOC strategies instead of interaction models in which SOC strategies mutually enhance each other's effects. Thus, a central assumption of SOC theory remains untested. In addition, most research on SOC strategies has been cross-sectional, assuming that employees' use of SOC strategies is stable over time. We conducted a quantitative diary study across nine work days (N = 77; 514 daily entries) to investigate interactive effects of daily SOC strategies on daily work engagement. Results showed that optimization and compensation, but not selection, had positive main effects on work engagement. Moreover, a significant three-way interaction effect indicated that the relationship between selection and work engagement was positive only when both optimization and compensation were high, whereas the relationship was negative when optimization was low and compensation was high. We discuss implications for future research and practice regarding the use of SOC strategies at work.

Open-pit block sequencing optimization: A mathematical model and solution technique

This study presents a comprehensive mathematical formulation model for a short-term open-pit mine block sequencing problem, which considers nearly all relevant technical aspects in open-pit mining. The proposed model aims to obtain the optimum extraction sequences of the original-size (smallest) blocks over short time intervals and in the presence of real-life constraints, including precedence relationship, machine capacity, grade requirements, processing demands and stockpile management. A hybrid branch-and-bound and simulated annealing algorithm is developed to solve the problem. Computational experiments show that the proposed methodology is a promising way to provide quantitative recommendations for mine planning and scheduling engineers.

Integration and optimization of pressure retarded osmosis with reverse osmosis for power generation and high efficiency desalination

Salinity gradient power is proposed as a source of renewable energy when two solutions of different salinity are mixed. In particular, Pressure Retarded Osmosis (PRO) coupled with a Reverse Osmosis process (RO) has been previously suggested for power generation, using RO brine as the draw solution. However, integration of PRO with RO may have further value for increasing the extent of water recovery in a desalination process. Consequently, this study was designed to model the impact of various system parameters to better understand how to design and operate practical PRO-RO units. The impact of feed salinity and recovery rate for the RO process on the concentration of draw solution, feed pressure, and membrane area of the PRO process was evaluated. The PRO system was designed to operate at maximum power density of . Model results showed that the PRO power density generated intensified with increasing seawater salinity and RO recovery rate. For an RO process operating at 52% recovery rate and 35 g/L feed salinity, a maximum power density of 24 W/m2 was achieved using 4.5 M NaCl draw solution. When seawater salinity increased to 45 g/L and the RO recovery rate was 46%, the PRO power density increased to 28 W/m2 using 5 M NaCl draw solution. The PRO system was able to increase the recovery rate of the RO by up to 18% depending on seawater salinity and RO recovery rate. This result suggested a potential advantage of coupling PRO process with RO system to increase the recovery rate of the desalination process and reduce brine discharge.